💧 Droplet Evaporation & d²-Law

Calculate single droplet lifetime, Godsave/Spalding evaporation constant K, convective Ranz-Marshall correction, and mass history.

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

📊 Solver Telemetry ● ACTIVE

👁️ Consultations 42

⚡ Calculs faits 35

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

📅 Mise en service Jun 2026

⏱️ Latence < 1 ms

💧 Droplet Shrinkage, Thermal Halo & Vapor Plume Dynamics

Real-time d²-Law Godsave/Spalding evaporation simulation

📝 Configuration & Presets

💧 Water Mist (150°C Air) 🔥 Diesel Combustion (850°C) 🥛 Spray Drying (200°C) ❄️ Cryogenic LN₂ Droplet

📐 Initial Droplet & Kinematics

Initial Diameter d₀ [µm]:

Relative Velocity U_rel [m/s]:

🌡️ Temperatures & Phase Change

Ambient Gas Temp T_∞ [°C]:

Droplet Surface Temp T_s [°C]: Typically saturation or wet-bulb temperature

Latent Heat of Vaporization L_v [kJ/kg]:

Liquid Density ρ_L [kg/m³]:

💨 Gas Transport Properties

Gas Thermal Conductivity k_g [W/(m·K)]:

Gas Specific Heat C_p,g [J/(kg·K)]:

Gas Density ρ_g [kg/m³]:

Gas Viscosity μ_g [Pa·s]:

Key Formulations:
• d²-Law: d²(t) = d₀² − K · t
• Lifetime: τ_life = d₀² / K
• Spalding Number: B_T = C_p,g (T_∞ − T_s) / L_v
• Constant: K = [8 k_g / (ρ_L C_p,g)] ln(1 + B_T) (1 + 0.3 Re^1/2 Pr^1/3)

📊 Evaporation Results

📊 Output Summary
💾 Fortran Source

Total Droplet Lifetime (τ_life) 130.24 ms 0.1302 seconds

Evaporation Constant (K) 0.019 mm²/s 1.920e-8 m²/s

Spalding Transfer Number (B_T) 0.0427 Ranz-Marshall Factor = 1.818x

Initial Mass Evap. Rate 1.505e-9 kg/s Stop Distance = 17.55 mm

📈 Normalized (d/d₀)² vs Time (d²-Law Verification)

📉 Droplet Diameter d(t) vs Time

=================================================================
 THERMOFLUIDCALC — DROPLET EVAPORATION & D2-LAW REPORT
=================================================================
Case Title                 : Water Mist Evaporation in Hot Air Stream
Initial Diameter (d0)      : 50.00 um (5.000e-5 m)
Ambient Temperature (T_inf): 150.00 C (423.15 K)
Surface Temperature (Ts)   : 55.00 C (328.15 K)
Relative Velocity (U_rel)  : 5.00 m/s
-----------------------------------------------------------------
Prandtl Number (Pr)        : 0.7454
Initial Reynolds (Re0)     : 9.043e+0
Spalding Number (BT)       : 0.0427
Convective Factor (F_conv) : 1.8180
-----------------------------------------------------------------
Evaporation Constant (K)   : 0.0192 mm2/s (1.9196e-8 m2/s)
Total Lifetime (tau_life)  : 130.238 ms (0.1302 s)
Initial Mass Evaporation   : 1.5046e-9 kg/s (1.505 ug/s)
Aerodynamic Stop Distance  : 17.55 mm
=================================================================

📘 Calculation Methodology & Engineering Theory

The Classical $d^2$-Law

Diffusion-controlled droplet evaporation under quasi-steady conditions follows Godsave and Spalding's classic linear diameter-squared relation:

                    d²(t) = d₀² − K · t,    τlife = d₀² / K
                

Where the evaporation constant $K$ depends on gas thermal conductivity and the thermodynamic driving force $\ln(1 + B_T)$.

Convective Enhancement (Ranz-Marshall)

When relative motion exists between droplet and surrounding gas ($U_{rel} > 0$), forced convection thins the boundary layer:

                    Nu = 2.0 + 0.6 · Red1/2 · Pr1/3
                

This accelerates heat transfer and mass evaporation rate by a convective multiplier $F_{conv}$.

Key Engineering Assumptions

Spherically symmetric liquid core throughout evaporation lifetime.
Droplet temperature remains near equilibrium wet-bulb/saturation value $T_s$.
Gas phase quasi-steady state assumption ($\tau_{gas} \ll \tau_{droplet}$).
Ideal gas mixture behavior in ambient gas film.

ThermoFluidCalc Engineering Report