Assalamu’alaikum warahmatullahi wabarakatuh
Perkenalkan saya Naysilla Salsabillah dengan NPM 2306155400 dari kelas Metode Numerik 03. Terima kasih telah meluangkan waktu untuk mengunjungi blog saya. Selamat membaca, semoga membawa keberkahan dan bermanfaat bagi semua.
Pemahaman metode numerik dalam analisis konduksi panas dua dimensi (2D) menjadi krusial dalam rekayasa termal, terutama ketika solusi analitis sulit diperoleh akibat kondisi batas yang kompleks. Salah satu pendekatan yang digunakan adalah DAI5 Framework, yang memungkinkan penyelesaian masalah numerik dengan struktur sistematis berbasis intention, initial condition, idealization, dan instruction set.
Konduksi panas dalam dua dimensi (2D) menggambarkan bagaimana energi termal berpindah melalui suatu material dalam dua arah, yaitu arah x dan arah y. Fenomena ini dikendalikan oleh persamaan diferensial parsial yang diturunkan dari Hukum Fourier, yang menyatakan bahwa laju perpindahan panas dalam suatu medium sebanding dengan gradien suhu dalam medium tersebut.
Secara matematis, konduksi panas dalam kondisi tunak (steady-state) dan tanpa adanya sumber panas internal dinyatakan sebagai berikut:
Persamaan ini menunjukkan bahwa perubahan suhu dalam arah x dan y harus memenuhi kondisi kesetimbangan termal, di mana panas yang masuk ke suatu elemen kecil dalam material sama dengan panas yang keluar, sehingga tidak ada akumulasi energi dalam sistem.
Jika terdapat sumber panas internal q (misalnya elemen pemanas atau reaksi eksotermik dalam material), maka persamaan di atas dimodifikasi menjadi:
dengan k adalah konduktivitas termal material (W/m-K), yang menentukan seberapa cepat panas dapat menyebar dalam material tersebut.
Karena persamaan diferensial parsial ini sulit diselesaikan secara analitis untuk kasus dengan kondisi batas kompleks, metode numerik seperti Finite Difference Method (FDM) digunakan untuk mendiskritisasi persamaan dalam bentuk persamaan aljabar yang dapat diselesaikan secara iteratif. Dalam pendekatan metode beda hingga, turunan kedua suhu terhadap x dan y dapat didekati dengan ekspresi:
sehingga persamaan konduksi panas diskritisasi menjadi:
Persamaan ini menunjukkan bahwa suhu di titik tengah (T i,j) adalah rata-rata dari suhu di sekitarnya, yang menggambarkan bagaimana panas menyebar dalam material secara dua dimensi.
Dalam studi ini, persamaan konduksi panas 2D diterapkan pada balok aluminium berukuran 1m × 1m dengan kondisi batas tertentu. Simulasi yang dilakukan menggunakan CFDSOF memberikan distribusi suhu yang sesuai dengan teori, di mana panas mengalir dari daerah bersuhu tinggi ke daerah bersuhu lebih rendah. Dengan pemodelan ini, analisis konduksi panas dapat dilakukan secara lebih akurat untuk memahami bagaimana energi termal berpindah dalam suatu material berdasarkan fungsi matematis yang mendasarinya.
dengan koefisien konduktivitas termal aluminium sebesar k=235k = 235 W/m-K. Persamaan ini menunjukkan bahwa distribusi suhu dalam material bergantung pada perubahan suhu dalam dua arah (x dan y), serta kondisi batas yang diterapkan pada sistem.
SOAL SIMULASI HEAT CONDUCTION 2D PADA CFDSOF
Pada simulasi ini, sebuah balok aluminium dengan dimensi 1 meter × 1 meter digunakan sebagai objek penelitian. Sumber panas diberikan di bagian atas balok dengan suhu 500°C, sementara bagian bawah memiliki suhu tetap sebesar 350°C. Selain itu, bagian samping kiri dan kanan balok masing-masing memiliki suhu tetap sebesar 425°C.
Material yang digunakan dalam simulasi adalah aluminium, dengan konduktivitas termal sebesar 235 W/m-K, yang menunjukkan kemampuan material dalam menghantarkan panas dengan efisiensi tinggi. Simulasi dilakukan dalam kondisi steady-state, di mana suhu di setiap titik dalam material tetap konstan seiring waktu. Selain itu, asumsi idealisasi yang digunakan adalah tidak adanya kesalahan relatif dalam nilai suhu maupun dimensi balok, sehingga hasil perhitungan numerik dapat memberikan gambaran yang lebih akurat tentang distribusi suhu dalam material.
Simulasi ini diawali dengan pembuatan geometri balok dalam perangkat lunak CFDSOF dengan ukuran 1 meter × 1 meter. Setelah itu, properti material ditetapkan sesuai dengan karakteristik aluminium, termasuk konduktivitas termalnya sebesar 235 W/m-K.
Setelah mendefinisikan geometri dan material, kondisi batas diterapkan sesuai skenario yang telah ditentukan. Suhu sebesar 500°C diberikan di bagian atas balok sebagai sumber panas, sementara bagian bawah balok dipertahankan pada suhu 350°C. Bagian samping kiri dan kanan memiliki suhu tetap sebesar 425°C. Dengan kondisi ini, panas akan menyebar dari bagian atas menuju bagian bawah dan samping, sesuai dengan hukum konduksi panas.
Langkah selanjutnya adalah melakukan meshing pada geometri balok untuk mendiskritisasi domain menjadi elemen-elemen kecil. Proses ini bertujuan untuk meningkatkan akurasi solusi numerik. Setelah itu, solver steady-state digunakan untuk menghitung distribusi suhu dalam balok, dengan pendekatan numerik yang sesuai. Terakhir, hasil simulasi divisualisasikan dalam bentuk kontur suhu atau grafik distribusi panas, yang akan memberikan gambaran bagaimana panas menyebar dalam sistem.
Hasil simulasi dalam CFDSOF menunjukkan kontur suhu yang mengilustrasikan bagaimana panas menyebar dari sumber panas utama (500°C) ke seluruh balok. Bagian tengah balok akan memiliki suhu di antara 350°C hingga 500°C, bergantung pada gradien termal yang terbentuk. Hasil ini menunjukkan bahwa panas menyebar secara merata dalam material, mengingat aluminium memiliki konduktivitas termal tinggi, yang memungkinkan transfer energi yang efisien.
Distribusi suhu yang dihasilkan dari simulasi ini juga menunjukkan pola konduksi yang stabil, sesuai dengan prinsip dasar hukum Fourier. Dengan mengetahui pola ini, para insinyur dapat memahami bagaimana panas bergerak dalam suatu sistem dan mengoptimalkan desain termal untuk berbagai aplikasi industri.
Melalui pendekatan metode numerik dengan DAI5 Framework, analisis konduksi panas dalam sistem dua dimensi dapat dilakukan secara lebih sistematis dan akurat. Simulasi yang diimplementasikan dalam CFDSOF menunjukkan bagaimana panas menyebar dari sumber panas utama menuju area sekitarnya, dengan gradien suhu yang bergantung pada konduktivitas termal aluminium serta kondisi batas yang diterapkan.
Pemanfaatan metode numerik seperti Finite Difference Method (FDM) memungkinkan penyelesaian persamaan diferensial parsial secara lebih fleksibel, sehingga dapat diterapkan dalam berbagai kasus teknik termal lainnya. Dengan pemahaman yang lebih dalam terhadap metode numerik dalam heat conduction 2D dapat mengoptimalkan desain termal untuk berbagai aplikasi industri, mulai dari sistem pendinginan, material insulasi termal, hingga efisiensi energi dalam bangunan hijau.
