digilib@itb.ac.id +62 812 2508 8800

Dalam tesis ini, solusi dari masalah Stefan, baik satu dimensi maupun dua dimensi akan dipelajari. Untuk memperoleh solusi numerik akan digunakan formulasi entalpi yang didiskritisasi dengan metode Godunov. Dengan menggunakan formulasi entalpi, persamaan pengatur dari sistem adalah sama, apapun fasenya, liquid atau solid. Untuk masalah Stefan satu dimensi akan dibahas dua kasus. Kasus pertama, efek-efek perubahan densitas selama perubahan fase diabaikan. Kasus ini akan diselesaikan dengan metode entalpi biasa. Kasus kedua, efek perubahan densitas seperti ekspansi volume dan pergerakan material akan diperhitungkan. Untuk kasus kedua ini akan diselesaikan dengan metode entalpi yang telah dimodifikasi. Hasil-hasil numerik untuk kedua kasus ini akan dibandingkan dengan solusi eksaknya. Metode serupa akan diterapkan ke dalam masalah Stefan dua dimensi. Hasil-hasil simulasi dua dimensi untuk cryosurgery (penggunaan temperatur sangat dingin untuk memusnahkan sel-sel kanker) akan ditunjukkan. Dalam simulasi tersebut diperlihatkan profil temperatur dan posisi interface yang merupakan informasi penting untuk mengontrol proses cryosurgery.