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

Permasalahan Stefanmerupakanmodelmatematikayangmenggambarkanprosesperu- bahan fase,sepertipelelehandanpembekuan,denganantarmukabergerakyangposi- sinyaberubahterhadapwaktu.Keberadaanantarmukabergerakmenyebabkansolusi masalah inisulitditentukansecaralangsung.Padapenelitianinidigunakanformulasi entalpi sehinggaposisiantarmukatidakperludilacaksecaraeksplisit.Modeltersebutdi- diskritisasi menggunakanmetodevolumehinggadandiimplementasikanpadasimulasi permasalahan Stefansatudanduadimensidenganberbagaikondisibatas. Hasil validasiskemanumerikpadakasuspencairansatudimensimenunjukkanhasil yang konsistendengansolusianalitik.Simulasipembekuansatudimensimenunjukkan bahwajeniskondisibatasmemengaruhikarakteristikpergerakanantarmuka,dimana syarat batasDirichletmenghasilkanperubahanyanglebihcepatpadaawalwaktudiban- dingkan syaratbatasNeumann.PerluasanskemapadapermasalahanStefanduadimensi diterapkan padasimulasipembekuanareatanahhomogendisekitarbantaransungai akibat fluktuasisuhuudaratahunandinegaraempatmusim.Hasilsimulasimenunjukkan bahwaareadekatpermukaanmengalamipembekuanpadarentangNovember–Februari, sementara lapisanyanglebihdalamrelatifstabilpadasuhurata-ratatahunansekitar 8.8?C