ccitonline.com

Simulasi Konduksi Panas Satu Dimensi Menggunakan CFDSOF

Awareness of I

Dalam dunia teknik, metode numerik digunakan untuk membantu menjelaskan berbagai fenomena fisika yang sulit diamati secara langsung. Salah satu fenomena yang sering dianalisis adalah perpindahan panas secara konduksi. Konduksi panas merupakan proses perpindahan energi termal melalui suatu material akibat adanya perbedaan temperatur tanpa disertai perpindahan zat atau materialnya.

Melalui pembelajaran ini, saya memahami bahwa proses perpindahan panas tidak hanya berkaitan dengan penggunaan persamaan matematika, tetapi juga berkaitan dengan cara memahami keteraturan fenomena di sekitar kita. Ketika terdapat dua daerah dengan suhu yang berbeda, energi panas akan berpindah dari daerah bersuhu lebih tinggi menuju daerah yang bersuhu lebih rendah. Fenomena tersebut tampak sederhana, namun sebenarnya dapat dijelaskan secara ilmiah menggunakan pendekatan matematis dan simulasi numerik.

Pada tahap ini saya menyadari bahwa memahami konduksi panas memerlukan proses yang bertahap. Tidak cukup hanya mengetahui bentuk persamaannya, tetapi juga perlu memahami arti fisik dari setiap parameter yang digunakan. Dengan demikian, hasil simulasi tidak hanya dipandang sebagai kumpulan angka ataupun visualisasi, melainkan sebagai representasi dari kondisi fisik yang benar-benar terjadi.

Intention

Tujuan dari simulasi ini adalah untuk membangun model matematika perpindahan panas satu dimensi serta menerapkannya menggunakan perangkat lunak CFDSOF. Model tersebut digunakan untuk melihat bagaimana distribusi temperatur terjadi pada suatu medium ketika diberikan kondisi batas tertentu.

Selain itu, simulasi ini bertujuan untuk membantu memahami keterkaitan antara teori perpindahan panas dengan visualisasi hasil simulasi yang dihasilkan perangkat lunak. Distribusi temperatur dapat diamati melalui grafik maupun kontur temperatur sehingga proses perpindahan panas menjadi lebih mudah dipahami secara visual.

Melalui simulasi ini juga dapat dipahami bahwa kondisi batas sangat memengaruhi hasil akhir perhitungan. Nilai temperatur pada dinding, karakteristik material, serta bentuk domain simulasi akan menentukan pola penyebaran panas yang terbentuk. Oleh karena itu, setiap langkah dalam proses simulasi perlu dilakukan dengan cermat agar hasil yang diperoleh dapat mendekati kondisi sebenarnya.

Intentional Thinking

Dasar teori utama dalam perpindahan panas secara konduksi adalah Hukum Fourier. Hukum ini menyatakan bahwa laju perpindahan panas dipengaruhi oleh konduktivitas termal material serta gradien temperatur terhadap posisi.

Secara matematis, Hukum Fourier dapat dituliskan sebagai berikut:$$q = -k \frac{\partial T}{\partial x}$$q=−k∂x∂T​

Keterangan:

• $q$q = laju perpindahan panas per satuan luas $(W/m^2)$(W/m2)
• $k$k = konduktivitas termal material $(W/m.K)$(W/m.K)
• $\frac{\partial T}{\partial x}$∂x∂T​ = gradien temperatur terhadap arah $x$x

Tanda negatif pada persamaan menunjukkan bahwa panas selalu bergerak dari temperatur yang lebih tinggi menuju temperatur yang lebih rendah. Dengan kata lain, arah aliran panas berlawanan dengan arah kenaikan temperatur.

Selanjutnya, prinsip konservasi energi digunakan untuk menurunkan persamaan konduksi panas. Misalnya terdapat elemen kecil pada batang dengan panjang $\Delta x$Δx, luas penampang $A$A, massa jenis $\rho$ρ, dan kalor jenis $c_p$cp​. Energi panas yang masuk dikurangi energi yang keluar akan sama dengan perubahan energi yang tersimpan di dalam elemen tersebut.

Secara sederhana dapat dituliskan:$$\text{Energi masuk} – \text{Energi keluar} = \text{Perubahan energi dalam}$$Energi masuk−Energi keluar=Perubahan energi dalam

Perubahan energi dalam elemen dapat dinyatakan sebagai:$$\rho c_p A \Delta x \frac{\partial T}{\partial t}$$ρcp​AΔx∂t∂T​

Dengan memasukkan Hukum Fourier ke dalam persamaan konservasi energi, diperoleh persamaan konduksi panas satu dimensi:$$\rho c_p \frac{\partial T}{\partial t} = k \frac{\partial^2 T}{\partial x^2}$$ρcp​∂t∂T​=k∂x2∂2T​

Jika nilai konduktivitas termal dianggap konstan, maka persamaan tersebut dapat disederhanakan menjadi:$$\frac{\partial T}{\partial t} = \alpha \frac{\partial^2 T}{\partial x^2}$$∂t∂T​=α∂x2∂2T​

dengan:$$\alpha = \frac{k}{\rho c_p}$$α=ρcp​k​

Besaran $\alpha$α disebut sebagai difusivitas termal, yaitu parameter yang menunjukkan seberapa cepat panas dapat menyebar di dalam suatu material. Semakin besar nilai difusivitas termal, maka semakin cepat pula perubahan temperatur terjadi di dalam material tersebut.

Idealization

Agar proses pemodelan lebih sederhana dan mudah dianalisis, beberapa asumsi digunakan dalam simulasi ini. Penyederhanaan tersebut dilakukan bukan untuk menghilangkan fenomena sebenarnya, tetapi untuk mempermudah proses analisis numerik.

Pada simulasi ini, perpindahan panas diasumsikan hanya terjadi dalam satu arah, yaitu arah $x$x. Material dianggap homogen sehingga sifat termalnya sama di seluruh bagian material. Selain itu, nilai konduktivitas termal diasumsikan tetap selama proses berlangsung. Material juga dianggap tidak mengalami perubahan bentuk selama perpindahan panas terjadi.

Sebelum simulasi dilakukan, kondisi awal dan kondisi batas perlu ditentukan terlebih dahulu. Kondisi awal menunjukkan distribusi temperatur pada saat awal pengamatan, sedangkan kondisi batas menggambarkan keadaan temperatur ataupun fluks panas pada bagian tepi domain.

Kondisi awal dapat dituliskan sebagai:$$T(x,0) = T_i$$T(x,0)=Ti​

Persamaan tersebut menunjukkan bahwa temperatur awal di sepanjang batang bernilai $T_i$Ti​.

Sedangkan kondisi batas temperatur tetap dapat dituliskan sebagai:$$T(0,t) = T_1$$T(0,t)=T1​ $$T(L,t) = T_2$$T(L,t)=T2​

Artinya, ujung kiri batang memiliki temperatur $T_1$T1​, sedangkan ujung kanan memiliki temperatur $T_2$T2​.

Selain temperatur tetap, terdapat pula kondisi batas adiabatik. Pada kondisi ini tidak terjadi perpindahan panas melalui dinding. Secara matematis dapat dituliskan sebagai:$$\frac{\partial T}{\partial x} = 0$$∂x∂T​=0

Dengan adanya kondisi awal dan kondisi batas tersebut, model matematika kemudian dapat diselesaikan secara numerik menggunakan perangkat lunak CFDSOF.

Instruction Set

Model matematika utama yang digunakan dalam simulasi konduksi panas satu dimensi adalah:$$\frac{\partial T}{\partial t} = \alpha \frac{\partial^2 T}{\partial x^2}$$∂t∂T​=α∂x2∂2T​

Persamaan tersebut menjadi dasar utama dalam proses simulasi menggunakan CFDSOF. Langkah pertama yang dilakukan adalah membuat domain simulasi. Domain ini merepresentasikan area tempat proses perpindahan panas berlangsung. Setelah domain dibuat, ukuran panjang, lebar, dan jumlah cell ditentukan sesuai kebutuhan simulasi.

Jumlah cell berpengaruh terhadap tingkat ketelitian hasil simulasi. Semakin banyak jumlah cell yang digunakan, maka distribusi temperatur yang dihasilkan akan semakin detail. Akan tetapi, penggunaan cell yang terlalu banyak juga menyebabkan waktu komputasi menjadi lebih lama. Oleh sebab itu, jumlah cell perlu disesuaikan dengan kebutuhan analisis agar hasil yang diperoleh tetap efisien.

Setelah domain selesai dibuat, setiap sisi domain diberikan nama seperti W1, W2, W3, W4, dan W5. Masing-masing dinding memiliki fungsi yang berbeda dalam simulasi. Beberapa dinding digunakan sebagai batas adiabatik, sebagian lainnya berfungsi sebagai dinding konduksi, dan ada pula yang bertindak sebagai sumber panas dalam sistem simulasi tersebut.