Siapa pun yang pernah menguji suhu wajan besi cor atau baja tahan karat dengan memercikkan beberapa tetes air ke permukaannya pasti pernah menyaksikan salah satu paradoks paling menakjubkan dalam fisika termal. Ketika suhu permukaan wajan hanya sedikit di atas titik didih air — sekitar $110^\circ\text{C}$ hingga $130^\circ\text{C}$ — tetesan air yang jatuh langsung menyebar, mendesis keras, dan menguap habis dalam waktu kurang dari satu detik.

Namun, jika wajan yang sama terus dipanaskan hingga melampaui $200^\circ\text{C}$ atau $300^\circ\text{C}$, perilaku air berbalik secara drastis. Alih-alih langsung lenyap menjadi uap, tetesan air justru menggumpal menjadi bulatan berkilau bagaikan bola merkuri. Tetesan itu tidak lagi mendesis, sama sekali tidak membasahi logam, melainkan meluncur dan memantul ke segala arah di atas wajan dengan gesekan hampir nol, serta bertahan hidup selama beberapa menit — hingga seratus kali lebih lama dibandingkan pada suhu $120^\circ\text{C}$!

Fenomena yang melawan intuisi ini dikenal sebagai efek Leidenfrost (Leidenfrost effect). Di atas ambang suhu kritis yang disebut suhu Leidenfrost ($T_L$), bagian bawah tetesan air menguap begitu cepat saat mendekati substrat panas sehingga membentuk bantalan uap (lapisan uap / vapor cushion) kontinu setebal $10$ hingga $100\,\mu\text{m}$. Lapisan uap mikroskopis ini bekerja ganda: secara hidrodinamis sebagai hovercraft yang menopang berat tetesan melalui tekanan pelumasan viskos tanpa kontak padat-cair, dan secara termal sebagai isolator ulung karena konduktivitas termal uap air ($\kappa_v \approx 0.03\,\text{W/(m}\cdot\text{K)}$) lebih dari dua puluh kali lebih rendah daripada air cair ($\kappa_l \approx 0.68\,\text{W/(m}\cdot\text{K)}$).

Dalam artikel ini, kita menurunkan fisika efek Leidenfrost secara sistematis dari prinsip pertama: mulai dari risalah bersejarah Johann Gottlob Leidenfrost (1756) dan kurva pendidihan Shiro Nukiyama (1934), melewati batas hidrodinamis Fluks Panas Kritis (CHF) Novak Zuber dan Fluks Panas Minimum (MHF) Paul Berenson, hingga penurunan matematis hukum penskalaan pelumasan Biance, Clanet, dan Quéré (2003), vorteks termokapiler internal Marangoni (“roda Leidenfrost”), osilasi poligon bintang Rayleigh-Lamb, ratchet berpenggerak mandiri Heiner Linke (2006), serta aplikasi kritis pada pendinginan darurat kecelakaan kehilangan pendingin (LOCA) di reaktor nuklir.

1. “Bola Kristal” di Atas Sendok Besi Membara: Penemuan Leidenfrost (1756)

Meskipun dokter dan kimiawan Belanda Herman Boerhaave sempat mencatat perilaku alkohol di atas besi membara pada tahun 1732, studi eksperimental kuantitatif pertama diterbitkan pada tahun 1756 oleh dokter dan filsuf alam Jerman Johann Gottlob Leidenfrost (1715–1794) di Universitas Duisburg.

Dalam monograf bahasa Latinnya De Aquae Communis Nonnullis Qualitatibus Tractatus (Risalah tentang Beberapa Sifat Air Biasa), Leidenfrost mendeskripsikan eksperimennya menggunakan sendok besi yang dipanaskan di atas bara api hingga berpijar merah:

“Ke dalam sendok yang berpijar ini, teteskan satu tetes air biasa yang murni… Air tersebut tidak menempel pada besi dan tidak pula menyebar; sebaliknya, tepat pada saat bersentuhan, ia mengerut menjadi bulatan sferis yang terang dan transparan layaknya sebuah bola kristal (globulus crystallinus). Ia berputar cepat dan menari dengan ringan di dasar sendok tanpa mengeluarkan suara maupun uap yang tampak. Jauh dari langsung menghilang, ukurannya menyusut perlahan selama tiga puluh atau empat puluh detik — bahkan satu menit penuh — hingga akhirnya sebutir partikel amat kecil lenyap dengan bunyi letupan halus.”

Dengan mengukur waktu hidup tetesan menggunakan bandul seiring mendinginnya sendok, Leidenfrost membuktikan bahwa laju penguapan turun drastis di atas suhu kritis akibat adanya lapisan tipis uap (vaporis strato). Pada tahun 1842, fisikawan Prancis Pierre Hippolyte Boutigny menamai kondisi ini “keadaan sferoidal” dan bahkan mendemonstrasikan pembekuan air menjadi es di dalam cawan platina yang membara menggunakan sulfur dioksida cair dalam keadaan Leidenfrost.

2. Kurva Pendidihan Shiro Nukiyama (1934) dan Empat Rezim Pendidihan Kolam

Untuk memahami secara kuantitatif mengapa tetesan bertahan bermenit-menit pada $250^\circ\text{C}$ tetapi lenyap kurang dari satu detik pada $125^\circ\text{C}$, kita meninjau hubungan antara fluks panas permukaan $q''$ $[\text{W/m}^2]$ dan kelebihan suhu dinding (wall superheat) $\Delta T_e \equiv T_w - T_{\text{sat}}$ (dengan $T_{\text{sat}} = 100^\circ\text{C}$ untuk air pada $1\,\text{atm}$), yang pertama kali dipetakan secara utuh pada tahun 1934 oleh Shiro Nukiyama di Universitas Kekaisaran Tōhoku.

  flowchart TD
    subgraph BoilingCurve["Kurva Pendidihan Shiro Nukiyama (q'' vs. Superheat Dinding ΔTe = Tw - Tsat)"]
        R1["Rezim I: Konveksi Alami<br/>ΔTe < 5°C (Tw < 105°C)<br/>Tanpa nukleasi gelembung; konveksi satu fasa<br/>q'' ∝ ΔTe^(5/4)"]
        R2["Rezim II: Pendidihan Inti (Nucleate Boiling)<br/>5°C ≤ ΔTe ≤ 30°C (Tw = 105–130°C)<br/>Nukleasi gelembung hebat di rongga mikro & penguapan lapisan mikro<br/>q'' ∝ ΔTe^3 (Efisiensi pendinginan maksimum)"]
        CHF["Puncak: Fluks Panas Kritis (CHF / DNB)<br/>ΔTe ≈ 30°C, q''_max ≈ 1.1–1.26 MW/m²<br/>Batas Hidrodinamis Zuber (runtuhnya kolom uap)"]
        R3["Rezim III: Pendidihan Transisi<br/>30°C < ΔTe < 100–120°C<br/>Selimut uap parsial tidak stabil & pembasahan intermiten<br/>Fluks panas TURUN saat suhu naik! (∂q''/∂ΔTe < 0)"]
        LP["Minimum: Titik Leidenfrost (MHF)<br/>ΔTe,L ≈ 93–120°C (Tw ≈ 193–220°C)<br/>q''_min ≈ 18–25 kW/m² (50 kali lebih rendah dari CHF!)"]
        R4["Rezim IV: Pendidihan Film Stabil<br/>ΔTe ≥ ΔTe,L (Tw > 200°C)<br/>Bantalan uap isolator kontinu (κ_v ≈ 0.03 W/m·K)<br/>Radiasi termal (∝ Tw^4) dominan pada Tw > 500°C"]
    end

    R1 -->|"Awal Pendidihan Inti (ONB)"| R2
    R2 -->|"Penggabungan gelembung & krisis arus lawan"| CHF
    CHF -->|"Jalur suhu terkontrol"| R3
    R3 -->|"Stabilisasi penuh lapisan uap"| LP
    LP -->|"Rezim pendidihan film"| R4
    CHF -.->|"Lompatan Burnout pada fluks panas terkontrol (ΔTw > +800°C)"| R4

1. Rezim I: Konveksi Alami ($\Delta T_e < 5^\circ\text{C}$): Akibat tekanan Laplace $\Delta P = 2\gamma / r_c$, embrio uap di dalam celah mikro berjari-jari $r_c$ memerlukan $\Delta T_e > \Delta T_{\text{ONB}} \approx \frac{2\gamma T_{\text{sat}}}{\rho_v L_{lv} r_c}$ untuk tumbuh. Di bawah nilai ini, perpindahan panas berlangsung melalui konveksi alami satu fasa ($q'' \propto \Delta T_e^{5/4}$).
2. Rezim II: Pendidihan Inti ($5^\circ\text{C} \le \Delta T_e \le 30^\circ\text{C}$) dan Fluks Panas Kritis (CHF) Zuber: Penguapan lapisan mikro di bawah gelembung dan agitasi mikro-konvektif meningkatkan fluks panas sebanding dengan pangkat tiga superheat (korelasi Rohsenow: $q'' \propto \Delta T_e^3$). Pada $\Delta T_e \approx 30^\circ\text{C}$ ($T_w \approx 130^\circ\text{C}$), gelembung bergabung membentuk kolom uap kontinu. Pada tahun 1958, Novak Zuber menggabungkan ketidakstabilan Rayleigh-Taylor ($\lambda_{RT} = 2\pi\sqrt{3\gamma/[g(\rho_l - \rho_v)]}$) dengan ketidakstabilan Kelvin-Helmholtz untuk membuktikan rumus Fluks Panas Kritis (CHF): $$q''_{\max} = \frac{\pi}{24} \rho_v L_{lv} \left[ \frac{\gamma g (\rho_l - \rho_v)}{\rho_v^2} \right]^{1/4} \approx 0.131 \rho_v L_{lv} \left[ \frac{\gamma g (\rho_l - \rho_v)}{\rho_v^2} \right]^{1/4} \approx 1.1\text{–}1.26\,\text{MW/m}^2$$ Pada pemanasan listrik kawat platina, mencapai $q''_{\max}$ memicu lonjakan suhu seketika dari $130^\circ\text{C}$ ke atas $1100^\circ\text{C}$ (krisis pendidihan / burnout).
3. Rezim III: Pendidihan Transisi ($30^\circ\text{C} < \Delta T_e < 100\text{–}120^\circ\text{C}$): Lapisan uap parsial yang tidak stabil menyelimuti dinding; semakin tinggi $T_w$, semakin besar fraksi uap kering isolator sehingga fluks panas justru menurun saat suhu naik ($\partial q'' / \partial \Delta T_e < 0$).
4. Rezim IV: Pendidihan Film Stabil dan Titik Leidenfrost (MHF Berenson): Pada titik Leidenfrost ($T_L \approx 193^\circ\text{C}\text{–}220^\circ\text{C}$) kontak padat-cair hilang sepenuhnya. Pada tahun 1961, Paul J. Berenson menurunkan rumus Fluks Panas Minimum (MHF): $$q''_{\min} = 0.09 \rho_v L_{lv} \left[ \frac{\gamma g (\rho_l - \rho_v)}{(\rho_l + \rho_v)^2} \right]^{1/4} \approx 18\text{–}25\,\text{kW/m}^2$$ Nilai ini 50 kali lebih kecil dibandingkan puncak pendidihan inti!

3. Statika Kapiler dan Hukum Penskalaan Hidrodinamis Bantalan Uap (Biance, Clanet, & Quéré 2003)

3.1 Panjang Kapiler $\ell_c$ dan Morfologi Tetesan Levitasi

Dengan sudut kontak semu $\theta_c = 180^\circ$, keseimbangan antara tekanan kapiler Laplace ($\sim \gamma / R$) dan tekanan hidrostatik ($\sim \rho_l g R$) mendefinisikan panjang kapiler $\ell_c$:

Berdasarkan bilangan Bond $\text{Bo} = (R/\ell_c)^2$, terdapat tiga rezim bentuk:

1. Tetesan Kuasi-Sferis ($R < \ell_c \approx 2.5\,\text{mm}$): Berbentuk bola berjari-jari $R$ dengan bidang kontak datar kecil di bagian bawah berjari-jari $r_c = \sqrt{\frac{2}{3}}\frac{R^2}{\ell_c}$.
2. Genangan / Panekuk Leidenfrost ($\ell_c < R < 4\ell_c$): Gravitasi memipihkan tetesan menjadi cakram dengan ketebalan universal konstan $h_p = 2\ell_c \sin(180^\circ/2) = 2\ell_c \approx 5.0\,\text{mm}$, memiliki kantong uap cekung di tengah (dimple) yang dikelilingi leher cincin sempit (neck).
3. Cerobong Uap Rayleigh-Taylor ($R > R_{\max} \approx 3.83\ell_c \approx 9.6\,\text{mm}$): Untuk $R \gtrsim 1\,\text{cm}$, antarmuka uap-air mengalami ketidakstabilan Rayleigh-Taylor sehingga gelembung uap besar meletus secara periodik menembus bagian tengah genangan (vapor chimneys).

3.2 Penurunan Ketebalan Lapisan Uap $e$ dan Waktu Hidup $\tau$

Dalam makalah klasik tahun 2003 (Physics of Fluids), Anne-Laure Biance, Christophe Clanet, dan David Quéré menggabungkan konduksi Fourier dengan pelumasan viskos untuk genangan berjari-jari $R > \ell_c$:

1. Konduksi Panas Fourier Melintasi Ketebalan $e$: $$Q = \kappa_v \frac{\Delta T}{e} \pi R^2 = \dot{m} L_{lv} \quad \Longrightarrow \quad \dot{m} = \frac{\kappa_v \Delta T}{L_{lv} e} \pi R^2 \qquad \cdots (1)$$
2. Teori Pelumasan Viskos (Aliran Radial Poiseuille): Persamaan pelumasan $\mu_v \frac{\partial^2 u_r}{\partial z^2} = \frac{dP}{dr}$ menghasilkan profil kecepatan parabolik $u_r(r,z) = \frac{1}{2\mu_v}\left(-\frac{dP}{dr}\right)z(e-z)$ dan laju aliran volumetrik $Q_v(r) = \frac{\pi r e^3}{6\mu_v}\left(-\frac{dP}{dr}\right) = \frac{\dot{m}}{\rho_v}\frac{r^2}{R^2}$. Integrasi tekanan memberikan tekanan berlebih rata-rata bantalan uap: $$\Delta P_{\text{cushion}} = \frac{3 \mu_v \dot{m}}{2 \pi \rho_v e^3} \qquad \cdots (2)$$
3. Keseimbangan Hidrostatik dan Hukum Penskalaan $e \propto R^{1/2}$: Menyamakan $\Delta P_{\text{cushion}}$ dengan tekanan hidrostatik genangan $\Delta P_{\text{weight}} = 2\rho_l g \ell_c$ dan mensubstitusikan (1) ke (2) menghasilkan hukum Biance-Clanet-Quéré: $$e = \left(\frac{3}{4}\right)^{1/4} \left( \frac{\kappa_v \mu_v \Delta T}{\rho_v \rho_l g L_{lv} \ell_c} \right)^{1/4} R^{1/2} \approx 10\text{–}100\,\mu\text{m} \qquad \cdots (3)$$ Untuk $R = 5\,\text{mm}$ dan $\Delta T = 200\,\text{K}$, diperoleh $e \approx 54\,\mu\text{m}$ (sedangkan untuk tetesan kecil $R < \ell_c$ berlaku $e_{\text{small}} \propto \Delta T^{1/4} R^{5/4}$). Mengintegrasikan $-\frac{dM}{dt} = \dot{m} \propto R^{3/2}$ menghasilkan $R(t) = R_0(1 - t/\tau)^2$ dan total waktu hidup $\tau \propto \Delta T^{-3/4} R_0^{3/2}$ (atau $\tau \propto R_0^{5/4}$ hingga $R_0^2$).

4. Vorteks Termokapiler Marangoni, Bintang Rayleigh-Lamb, dan Ratchet Leidenfrost

  flowchart TD
    subgraph ActiveDynamics["Dinamika Aktif Non-Kesetimbangan pada Tetesan Leidenfrost"]
        Grad["Gradien Termal Vertikal Internal<br/>Dasar pada Tsat = 100°C vs. Puncak Dingin Evaporatif Tapex ≈ 95–98°C"]

        subgraph Marangoni["1. Konveksi Termokapiler Internal Marangoni"]
            M1["Gradien tegangan permukaan ∂γ/∂T ≈ -0.18 mN/(m·K)<br/>Menarik antarmuka dari dasar hangat menuju puncak dingin"]
            M2["Bilangan Marangoni tinggi (Ma ~ 10^4) memicu vorteks toroidal<br/>→ Pematahan simetri spontan menjadi satu 'Roda Leidenfrost'<br/>→ Penggerak mandiri di atas permukaan datar sempurna!"]
        end

        subgraph StarModes["2. Osilasi Bintang Rayleigh-Lamb"]
            S1["Fluktuasi tekanan uap di leher cincin (f_vapor ≈ 2 f_n)"]
            S2["Resonansi parametrik subharmonik Faraday<br/>ω_n² = n(n² - 1)γ / (ρ_l R³)<br/>Membentuk bintang berdenyut dengan n = 2, 3, 4, 5 ... 12 sudut!"]
        end

        subgraph Ratchet["3. Ratchet Leidenfrost Berpenggerak Mandiri (Linke et al. 2006)"]
            L1["Substrat kuningan bergigi gergaji asimetris (periode ~1.5 mm, kedalaman ~0.3 mm)"]
            L2["Rektifikasi geometris aliran uap viskos Poiseuille<br/>Tegangan geser viskos τ_visc ~ μ_v u_v / e mendorong tetesan MENANJAK pada 10–30 cm/s!"]
        end
    end

    Grad --> M1 --> M2
    Grad --> S1 --> S2
    Grad --> L1 --> L2

4.1 Konveksi Termokapiler Marangoni dan “Roda Leidenfrost”

Perbedaan suhu antara dasar ($100^\circ\text{C}$) dan puncak ($95\text{–}98^\circ\text{C}$) menghasilkan, melalui $\frac{\partial \gamma}{\partial T} \approx -0.18\,\text{mN/(m}\cdot\text{K)}$, tegangan geser Marangoni $\tau_s = \frac{\partial \gamma}{\partial T}\nabla_s T$ ($\text{Ma} = \frac{|\partial \gamma/\partial T|\Delta T_{\text{int}} R}{\mu_l \alpha_l} \sim 2 \times 10^4 \gg \text{Ma}_c$) yang menghilangkan efek cincin kopi (coffee-ring effect). Selain itu, Bouillant, Quéré et al. (2018, Nature Physics) menemukan bahwa untuk $R \lesssim 1.5\,\text{mm}$ vorteks toroidal secara spontan mematahkan simetrinya menjadi satu sel berputar (“roda Leidenfrost”) yang memiringkan bantalan uap dan menggelindingkan tetesan secara mandiri di atas bidang horizontal.

4.2 Osilasi Poligon Bintang Rayleigh-Lamb

Genangan besar ($R \sim 5\text{–}15\,\text{mm}$) membentuk gelombang berdiri azimut berbentuk poligon bintang yang diatur oleh relasi dispersi Rayleigh-Lamb:

yang dieksitasi melalui resonansi parametrik subharmonik Faraday ($f_n = \frac{\omega_n}{2\pi} \approx \frac{1}{2} f_{\text{vapor}}$) oleh fluktuasi tekanan aliran uap di leher.

4.3 Ratchet Berpenggerak Mandiri Heiner Linke (2006)

Pada tahun 2006, Heiner Linke et al. (Phys. Rev. Lett.) mendemonstrasikan bahwa di atas pelat panas dengan alur gigi gergaji asimetris (periode $\sim 1.5\,\text{mm}$, kedalaman $\sim 0.3\,\text{mm}$), tetesan Leidenfrost bergerak mandiri dengan kecepatan $10\text{–}30\,\text{cm/s}$ menaiki anak tangga curam: meniskus dasar menutup celah vertikal dan mengarahkan aliran uap Poiseuille ke belakang sepanjang lereng landai, menghasilkan gaya dorong geser viskos $F_{\text{ratchet}} \sim \frac{\mu_v u_{\text{vapor}}}{e}\pi R^2 \sim 30\,\mu\text{N}$.

5. Aplikasi Industri: Keselamatan Reaktor Nuklir (LOCA), Quenching Baja, dan Nitrogen Cair

1. Pendinginan Darurat Reaktor Nuklir Saat Kecelakaan Kehilangan Pendingin (LOCA): Jika fluks panas pada kelongsong bahan bakar Zircaloy melampaui CHF, sistem melompat ke rezim pendidihan film Leidenfrost (Rezim IV) dan suhu kelongsong melonjak di atas $1000^\circ\text{C}$ dalam hitungan detik. Saat sistem pendingin darurat (ECCS) menyemprotkan air, bantalan uap menolak cairan ($T_w \gg T_L$) dan garis depan pembasahan (quench front) hanya dapat maju melalui konduksi aksial di dalam logam kelongsong (model Yamanouchi: $u_q \propto (T_L - T_{\text{sat}})/(T_{\text{init}} - T_L)$).
2. Quenching Metalurgi Baja: Saat mencelupkan baja austenit $850^\circ\text{C}$, selimut uap Leidenfrost awal memperlambat pendinginan hingga suhu turun di bawah $T_L$, di mana pendidihan inti membentuk martensit. Runtuhnya lapisan uap yang tidak merata memicu distorsi dan retak quenching.
3. Jari Basah di Timbal Cair dan Nitrogen Cair di Kulit: Jari basah yang dicelupkan sepersekian detik ke dalam timbal cair $400^\circ\text{C}$ (sangat berbahaya, jangan pernah dicoba!) terlindungi oleh bantalan uap air instan; serupa dengan itu, nitrogen cair ($T_{\text{sat}} = 77\,\text{K} = -196^\circ\text{C}$) merasakan kulit manusia ($307\,\text{K}$, $\Delta T_e = 230\,\text{K}$) sebagai wajan membara dan langsung melayang di atas bantalan gas nitrogen.

Referensi

1. Leidenfrost, J. G. (1756). De Aquae Communis Nonnullis Qualitatibus Tractatus. Duisburg.
2. Nukiyama, S. (1934). “The Maximum and Minimum Values of the Heat Transmitted from Metal to Boiling Water under Atmospheric Pressure”. J. Jpn. Soc. Mech. Eng., 37(206), 367–374.
3. Zuber, N. (1958). “On the stability of boiling heat transfer”. Trans. ASME, 80(3), 711–720.
4. Berenson, P. J. (1961). “Film-boiling heat transfer from a horizontal surface”. J. Heat Transf., 83(3), 351–356.
5. Biance, A.-L., Clanet, C., & Quéré, D. (2003). “Leidenfrost drops”. Physics of Fluids, 15(6), 1632–1637.
6. Linke, H., et al. (2006). “Self-Propelled Leidenfrost Droplets”. Physical Review Letters, 96(15), 154502.
7. Quéré, D. (2013). “Leidenfrost Dynamics”. Annual Review of Fluid Mechanics, 45, 197–215.
8. Bouillant, A., et al. (2018). “Leidenfrost wheels”. Nature Physics, 14(12), 1188–1192.