❄️ Stefan Moving Boundary Solidification Front

Solve two-phase moving boundary Stefan phase change problems: solid front position s(t), total freezing time, solid Stefan number (Stes), and instantaneous cooling heat flux.

⚡ Fortran 90 Engine Double Precision (IEEE 754) ✓ ISO / ASME Validated

📊 Solver Telemetry ● ACTIVE

👁️ Consultations 43

⚡ Calculs faits 35

💾 Téléchargements 401 📦 Code Fortran 3.4 KB

📅 Mise en service Jun 2026

⏱️ Latence < 1 ms

❄️ Two-Phase Moving Solidification Front Interface s(t)

Real-time visual simulation: Frozen solid crystal layer advancing dynamically into liquid domain

📝 Configuration & Presets

🧊 Ice Rink Freezing Plate 🔋 Paraffin Wax PCM Battery 🥩 Industrial Food Freezing 🏭 Continuous Steel Casting

📐 Geometry & Phase Change Temperatures

Target Thickness L [mm]:

Phase Change Temp T_m [°C]:

Cold Wall Temp T₀ [°C]:

Initial Liquid Temp T_i [°C]:

🧪 Latent Heat & Solid Properties

Latent Heat L_f [kJ/kg]:

Solid Conductivity k_s [W/(m·K)]:

Solid Density ρ_s [kg/m³]:

Solid Specific Heat c_p,s [J/(kg·K)]:

Stefan Problem Formulation:
• Solid Stefan Number: Ste_s = c_p,s (T_m − T₀) / L_f
• Front Constant: λ · e^λ² · erf(λ) = Ste_s / √π
• Front Position: s(t) = 2 λ √(α_s · t) [mm]
• Freezing Time: t_freeze = L² / (4 λ² α_s) [hours]

📊 Stefan Solidification Results

📊 Output Summary
💾 Fortran Source

Total Freezing / Solidification Time

t = 1.26 hours (4523 s)

Front Growth Constant: λ = 0.3430 | 1-hr Thickness: 35.7 mm

Ste_s = 0.255

Solid Stefan Number (Ste_s) 0.255 Sensible / Latent heat

Solid Thermal Diffusivity (α_s) 7.52e-7 m²/s k_s / (ρ_s · c_p,s)

1-Hour Solidified Layer 35.7 mm 2 λ √(α_s · 3600)

Liquid Superheat (Ste_l) 0.244 c_p,l (T_i − T_m) / L_f

📈 Solidification Front Position s(t) [mm] vs Time (hours)

📉 Wall Heat Flux q″ [kW/m²] vs Time (hours)

=================================================================
 THERMOFLUIDCALC — STEFAN MOVING BOUNDARY PHASE CHANGE REPORT
=================================================================
Case Title                 : Industrial Cryogenic Tunnel Food Solidification
Target Thickness (L)       : 40.0 mm
Phase Change Temperatures  : Tm = -1.5 C, Wall T0 = -35.0 C, Liquid Ti = 15.0 C
Latent Heat & Solid Props  : Lf = 250.0 kJ/kg, ks = 1.40 W/m.K, rhos = 980.0 kg/m3 (alpha_s = 7.519e-7 m2/s)
-----------------------------------------------------------------
SOLID STEFAN NUMBER (Ste_s): 0.2546
LIQUID STEFAN NUMBER (Ste_l: 0.2442
FRONT GROWTH CONSTANT (lam): 0.3430
TOTAL SOLIDIFICATION TIME  : 1.256 hours (4522.8 seconds)
Layer Thickness after 1 hr : 35.69 mm
=================================================================

📘 Calculation Methodology & Stefan Problem Standards

Neumann Transcendental Solution

Solves the exact non-linear energy balance at the moving phase change interface $x = s(t)$:

                    λ · eλ² · erf(λ) = Stes / √π
                

Parabolic Growth Law

The solid front advances with the square root of time: $s(t) = 2\lambda \sqrt{\alpha_s t}$. Freezing rate decelerates as solid thermal resistance grows.

Key Engineering Assumptions

1D semi-infinite or planar slab geometry with sharp phase interface.
Constant solid thermophysical properties ($\rho_s, k_s, c_{p,s}$).
Applicable to ice makers, PCM thermal storage, casting, and food freezing.

ThermoFluidCalc Engineering Report