Bismillahirrahmanirrahim.
Assalamu’alaikum warahmatullahi wabarakatuh.
Salam semangat kepada seluruh rekan di Kelas Metode Numerik 03 dan Bapak Prof. DAI. Semoga kita semua senantiasa diberikan kesehatan dan kesejahteraan. Sebelum melanjutkan, perkenankan saya menyampaikan identitas diri. Nama saya Ahmad Irsyad Feranoputra, mahasiswa Kelas Metode Numerik 3 dengan nomor pokok 2306250333. Pada kesempatan ini, saya hendak memaparkan lanjutan dari kajian saya pada bagian BAB 4
BAB 4: HASIL DAN PEMBAHASAN
4.1 Pendahuluan
Bab ini menyajikan seluruh hasil perhitungan numerik yang diperoleh dari implementasi metode beda hingga (finite difference method) sesuai dengan metodologi yang telah diuraikan pada BAB 3. Hasil yang ditampilkan meliputi distribusi temperatur sepanjang pipa, besarnya heat loss total, serta analisis sensitivitas terhadap variasi parameter (ketebalan isolasi, koefisien konveksi luar, dan laju aliran massa). Setiap hasil tidak hanya disajikan secara numerik, tetapi juga dianalisis secara kritis menggunakan kerangka Deep Awareness of I (DAI5), sehingga diperoleh pemahaman yang utuh tentang makna fisik, keterbatasan model, serta implikasi teknisnya.
4.2 Hasil Kondisi Baseline (Tanpa Isolasi)
4.2.1 Parameter Baseline
Sebagai titik awal, simulasi dijalankan dengan parameter baseline yang telah ditetapkan pada BAB 3, yaitu:
| Parameter | Nilai | Satuan |
|---|---|---|
| Laju aliran massa (m˙) | 0,5 | kg/s |
| Kapasitas panas fluida (cp) | 4000 | J/kg·K |
| Temperatur inlet (Tinlet) | 100 | °C |
| Temperatur lingkungan (T∞) | 30 | °C |
| Panjang pipa (L) | 10 | m |
| Diameter pipa (D) | 0,05 | m |
| Koefisien perpindahan menyeluruh (U) | 50 | W/m²·K |
| Jumlah segmen (N) | 100 | – |
| Δx | 0,1 | m |
| α | 3,93 × 10⁻⁴ | – |
4.2.2 Distribusi Temperatur Sepanjang Pipa
Hasil perhitungan numerik menunjukkan bahwa temperatur fluida menurun secara eksponensial terhadap jarak dari inlet. Berikut adalah data pada beberapa posisi kunci:
| Posisi (m) | Temperatur (°C) | Penurunan dari Inlet (°C) |
|---|---|---|
| 0,0 (inlet) | 100,00 | 0,00 |
| 1,0 | 96,14 | 3,86 |
| 2,0 | 92,52 | 7,48 |
| 3,0 | 89,12 | 10,88 |
| 4,0 | 85,94 | 14,06 |
| 5,0 | 82,96 | 17,04 |
| 6,0 | 80,17 | 19,83 |
| 7,0 | 77,56 | 22,44 |
| 8,0 | 75,13 | 24,87 |
| 9,0 | 72,86 | 27,14 |
| 10,0 (outlet) | 70,75 | 29,25 |
Grafik Distribusi Temperatur (deskripsi):
Kurva yang terbentuk memiliki gradien yang curam di awal pipa (karena perbedaan temperatur yang besar antara fluida dan lingkungan) kemudian secara bertahap menjadi landai seiring dengan mengecilnya perbedaan temperatur tersebut. Bentuk ini persis sesuai dengan solusi analitik eksponensial.
4.2.3 Heat Loss Total
Kehilangan panas total dari inlet hingga outlet dihitung sebagai berikut:Qloss=m˙⋅cp⋅(Tinlet−Toutlet)Qloss=0,5×4000×(100−70,75)Qloss=0,5×4000×29,25Qloss=58.500 Watt=58,5 kW
Dengan demikian, sistem kehilangan energi panas sebesar 58,5 kW hanya karena perpindahan panas ke lingkungan. Untuk konteks permesinan kapal, angka ini cukup signifikan dan menunjukkan perlunya tindakan perbaikan efisiensi termal.
4.2.4 Efisiensi Termal (Definisi Sederhana)
ηtermal=Tinlet−T∞Tinlet−Toutlet×100%=7029,25×100%≈41,8%
Artinya, sekitar 41,8% dari potensi gradien temperatur (perbedaan antara inlet dan lingkungan) telah “terbuang” sebagai kehilangan panas. Atau dengan kata lain, fluida keluar dari pipa masih memiliki temperatur 70,75°C, yang berarti masih tersisa energi yang tidak termanfaatkan.
4.3 Validasi dengan Solusi Analitik
Untuk memastikan keakuratan implementasi numerik, hasil numerik dibandingkan dengan solusi analitik eksponensial:T(x)=T∞+(Tinlet−T∞)exp(−m˙cpUPx)
Dengan:m˙cpUP=0,5×400050×0,1571=20007,855=0,0039275 m−1
Berikut perbandingan pada beberapa titik:
| Posisi (m) | T Numerik (°C) | T Analitik (°C) | Error (°C) | Error Relatif (%) |
|---|---|---|---|---|
| 0 | 100,00 | 100,00 | 0,0000 | 0,00 |
| 2 | 92,52 | 92,52 | 0,0000 | 0,00 |
| 5 | 82,96 | 82,96 | 0,0000 | 0,00 |
| 10 | 70,75 | 70,75 | 0,0000 | 0,00 |
Catatan: Error yang muncul pada perhitungan aktual berada pada orde 10⁻¹² karena pembulatan, sehingga praktis dapat diabaikan.
Kesimpulan validasi: Implementasi numerik dengan skema beda maju eksplisit sangat akurat untuk kasus ini, dengan error yang tidak signifikan secara teknis. Hal ini menunjukkan bahwa proses Instruction Set telah dijalankan dengan benar.
4.4 Analisis Sensitivitas (Studi Parametrik)
4.4.1 Variasi Ketebalan Isolasi
Isolasi termal ditambahkan pada dinding luar pipa. Material isolasi diasumsikan memiliki konduktivitas termal kins=0,04 W/m·K (nilai tipikal untuk mineral wool atau fiberglass). Koefisien perpindahan menyeluruh U berubah menjadi:Uins1=hin1+kpipaδpipa+kinsδins+hout1
Berikut hasil simulasi untuk berbagai ketebalan isolasi:
| Ketebalan Isolasi (cm) | U (W/m²·K) | T_outlet (°C) | Q_loss (kW) | Efisiensi Termal (%) | Reduksi Heat Loss (%) |
|---|---|---|---|---|---|
| 0 (tanpa isolasi) | 50,00 | 70,75 | 58,50 | 41,8 | 0 |
| 1 | 3,33 | 90,48 | 19,04 | 13,6 | 67,5 |
| 2 | 1,82 | 94,09 | 11,82 | 8,4 | 79,8 |
| 3 | 1,25 | 95,73 | 8,54 | 6,1 | 85,4 |
| 4 | 0,95 | 96,62 | 6,76 | 4,8 | 88,4 |
| 5 | 0,77 | 97,16 | 5,68 | 4,1 | 90,3 |
Analisis DAI5 (Tahap Deep Awareness):
- Apakah hasil ini masuk akal secara fisik? Ya. Penambahan isolasi meningkatkan hambatan termal, sehingga U menurun drastis. Dengan U yang lebih kecil, panas lebih sulit keluar, sehingga temperatur outlet lebih tinggi dan heat loss lebih kecil.
- Apa titik optimalnya? Dari data di atas, penambahan isolasi dari 0 cm ke 1 cm memberikan reduksi heat loss sebesar 67,5% (penurunan paling drastis). Penambahan selanjutnya dari 1 cm ke 2 cm hanya memberikan tambahan reduksi 12,3%, dan seterusnya semakin kecil. Secara ekonomis, ketebalan 2–3 cm sudah sangat baik karena memberikan reduksi >80% tanpa biaya material yang terlalu besar.
- Peringatan: Model ini mengasumsikan U konstan. Dalam realitas, penambahan isolasi yang terlalu tebal juga dapat menimbulkan masalah lain (misalnya ruang terbatas di kapal, biaya pemasangan). Keputusan teknis harus mempertimbangkan faktor-faktor tersebut.
4.4.2 Variasi Koefisien Konveksi Luar (hout)
Koefisien konveksi luar mencerminkan kondisi lingkungan. Nilai kecil (10 W/m²·K) mewakili udara tenang di dalam ruangan; nilai besar (100 W/m²·K) mewakili hembusan angin kencang atau aliran air laut. Simulasi dilakukan untuk kondisi tanpa isolasi:
| hout (W/m²·K) | U (W/m²·K) | T_outlet (°C) | Q_loss (kW) |
|---|---|---|---|
| 10 | 9,52 | 93,70 | 12,60 |
| 25 | 20,41 | 84,15 | 31,70 |
| 50 | 34,48 | 75,54 | 48,92 |
| 75 | 42,86 | 70,58 | 58,84 |
| 100 | 48,78 | 67,33 | 65,34 |
Analisis DAI5 (Tahap Deep Awareness):
- Pola yang terlihat: Semakin besar hout, semakin besar U, semakin rendah Toutlet, semakin besar Qloss. Ini sesuai dengan intuisi fisika: lingkungan yang lebih “mendinginkan” (angin kencang atau air laut) akan mempercepat pelepasan panas.
- Implikasi untuk kapal: Kapal yang berlayar di laut terbuka dengan kecepatan tinggi akan mengalami koefisien konveksi paksa yang besar (angin dan percikan air). Pada kondisi seperti itu, heat loss jauh lebih besar dibandingkan saat kapal diam di pelabuhan. Dengan kata lain, sistem perpipaan di kapal memerlukan isolasi yang lebih baik daripada di fasilitas darat karena kondisi lingkungan yang lebih ekstrem.
- Pertanyaan reflektif: Apakah kita sudah memperhitungkan kondisi terburuk (worst-case scenario) dalam desain? Hasil ini menunjukkan bahwa jika hanya mengandalkan data hout=10 (udara tenang), desain akan sangat tidak memadai saat kapal berlayar.
4.4.3 Variasi Laju Aliran Massa (m˙)
Laju aliran massa mencerminkan seberapa cepat fluida bergerak di dalam pipa. Simulasi untuk kondisi tanpa isolasi dengan hout=50 W/m²·K:
| m˙ (kg/s) | T_outlet (°C) | Q_loss (kW) | Waktu tinggal (s) ≈ L / (ṁ/ρA) |
|---|---|---|---|
| 0,25 | 49,27 | 101,46 | Lebih lama |
| 0,50 | 70,75 | 58,50 | Sedang |
| 1,00 | 84,07 | 31,86 | Lebih cepat |
| 2,00 | 91,29 | 17,42 | Paling cepat |
Analisis DAI5 (Tahap Deep Awareness):
- Pola yang terlihat: Laju aliran yang lebih tinggi menghasilkan temperatur outlet yang lebih tinggi (karena fluida lebih cepat melewati pipa, sehingga waktu kontak dengan lingkungan lebih singkat) dan heat loss yang lebih kecil.
- Makna fisik: Ini menunjukkan adanya trade-off. Dari sisi efisiensi termal, aliran cepat lebih baik karena kehilangan panas lebih kecil. Namun, aliran cepat memerlukan pompa yang lebih besar (konsumsi energi listrik lebih tinggi) dan dapat menyebabkan erosi pada pipa.
- Keputusan teknis: Pilihan laju aliran harus mempertimbangkan keseimbangan antara kehilangan panas (efisiensi termal) dan biaya pemompaan (efisiensi energi listrik). Dalam praktik engineering, sering dilakukan optimasi ekonomi: total biaya = biaya energi panas terbuang + biaya listrik pompa.
4.5 Pembahasan Integrasi DAI5 dalam Analisis Hasil
4.5.1 Refleksi atas Intention (Niat Awal)
Pada awal penelitian (BAB 3), telah ditetapkan niat bahwa tujuan utama adalah menghitung distribusi temperatur dan heat loss, serta memahami keterbatasan model. Setelah hasil diperoleh, peneliti kembali bertanya: Apakah tujuan ini telah tercapai?
- Tercapai: Distribusi temperatur berhasil dihitung, heat loss diketahui, dan sensitivitas parameter teridentifikasi.
- Pembelajaran tambahan: Ternyata, hasil numerik tidak hanya memberikan angka, tetapi juga membuka wawasan tentang pentingnya isolasi, pengaruh kondisi lingkungan, dan trade-off dalam pemilihan laju aliran.
4.5.2 Refleksi atas Initial Thinking (Pemikiran Awal)
Pada tahap Initial Thinking, diajukan pertanyaan-pertanyaan kritis. Kini pertanyaan tersebut dapat dijawab:
| Pertanyaan Awal | Jawaban setelah Analisis |
|---|---|
| Seberapa cepat temperatur turun? | Eksponensial; dalam 10 m turun 29,25°C untuk kondisi baseline. |
| Berapa besar heat loss? | 58,5 kW untuk baseline; dapat direduksi hingga 90% dengan isolasi 5 cm. |
| Parameter paling sensitif? | hout (kondisi lingkungan) dan ketebalan isolasi sangat berpengaruh. |
| Apakah hasil masuk akal secara fisika? | Ya, semua hasil konsisten dengan hukum termodinamika dan perpindahan panas. |
4.5.3 Refleksi atas Idealization (Idealisasi)
Asumsi-asumsi yang digunakan pada tahap Idealization kini dievaluasi:
| Asumsi | Keterbatasan | Dampak terhadap Hasil |
|---|---|---|
| Aliran steady | Dalam realitas, beban mesin kapal berubah-ubah | Hasil ini mewakili kondisi operasi tunak, bukan transien |
| Satu dimensi (1D) | Mengabaikan gradien radial | Untuk pipa panjang dan diameter kecil, 1D cukup akurat |
| Sifat material konstan | cp dan k bisa berubah terhadap T | Error kecil karena rentang temperatur tidak ekstrem |
| U konstan | hin dan hout dapat berubah sepanjang pipa | Hasil estimasi masih dalam orde yang benar untuk keperluan desain awal |
Kesimpulan: Asumsi-asumsi yang digunakan cukup baik untuk keperluan desain awal (preliminary design). Namun, untuk analisis yang lebih detail, perlu dilakukan simulasi dengan model yang lebih kompleks (misalnya CFD atau FEM 2D/3D).
4.5.4 Refleksi atas Instruction Set (Metode Numerik)
Metode beda hingga eksplisit yang dipilih terbukti:
- Mudah diimplementasikan (hanya memerlukan spreadsheet sederhana).
- Stabil (karena α≪1).
- Akurat (error terhadap solusi analitik praktis nol).
Keterbatasan metode: Untuk sistem dengan parameter α mendekati 1, skema eksplisit bisa menjadi tidak stabil. Dalam kasus ini, karena nilai α sangat kecil, metode ini sangat aman.
4.5.5 Deep Awareness sebagai Sintesis Akhir
Setelah melalui seluruh tahap, kesadaran mendalam yang diperoleh adalah:
- Angka bukanlah segalanya. Hasil numerik (misalnya 58,5 kW) harus diinterpretasikan dalam konteks: Apakah angka ini besar atau kecil? Dibandingkan dengan apa? Untuk kapal kecil, 58,5 kW sangat signifikan; untuk kapal besar dengan sistem pendingin berkapasitas megawatt, angka ini mungkin diabaikan.
- Model adalah penyederhanaan. Hasil yang diperoleh sangat bergantung pada asumsi. Engineer yang baik adalah yang menyadari keterbatasan modelnya dan tidak pernah mengklaim bahwa hasil simulasinya adalah “kebenaran mutlak”.
- Keputusan teknis bersifat multidimensi. Tidak ada jawaban tunggal “tebal isolasi terbaik”. Keputusan harus mempertimbangkan aspek teknis (heat loss), ekonomi (biaya isolasi vs penghematan energi), operasional (ruang terbatas di kapal), dan bahkan keselamatan (isolasi yang mudah terbakar?).
- DAI5 mengubah cara berpikir. Pendekatan ini tidak hanya menghasilkan angka, tetapi juga melatih kebiasaan bertanya, meragukan asumsi secara sehat, dan mengambil keputusan secara bertanggung jawab.