❄️ Freeze Drying Sublimation (Pikal Lyophilization)

Simulate pharmaceutical primary freeze-drying sublimation cycle times, moving ice front temperature, dry cake vapor resistance (Rp), and vial heat flux.

⚡ Fortran 90 Engine Double Precision (IEEE 754) ✓ ISO / ASME Validated
📊 Solver Telemetry ● ACTIVE
👁️ Consultations 43
⚡ Calculs faits 33
💾 Téléchargements 397 📦 Code Fortran 6.4 KB
📅 Mise en service Jun 2026
⏱️ Latence < 1 ms
⚡ Outils & Rapports :
💾 Télécharger Fortran 90

❄️ Pharmaceutical Vial Sublimation Front Progression & Vapor Flow

Real-time visual simulation of ice sublimation interface descending through porous cake under chamber vacuum

📝 Configuration & Presets

💉 Biopharma mAb (6R Vial) 🧪 Vaccine mRNA LNP (-33°C) 💊 Antibiotic Cake (10R) 🔬 Diagnostic Enzyme
🧪 Vial Formulation & Geometry
2R: ~2.0 cm², 6R: ~3.8 cm², 10R: ~5.5 cm²
❄️ Lyophilizer Operating Conditions
100 mTorr = 0.133 mbar = 13.3 Pa
Pikal Lyophilization Formulations:
• Sublimation Mass Flux: ṁsub = (Psat(Tice) − Pch) / Rp [g/(cm²·h)]
• Shelf Heat Flux: q̇ = Kv · (Tshelf − Tbottom) = ṁsub · ΔHsub
• Primary Drying Time: tp = ∫ [ ρice (1 − cs) / ṁsub ] dLdry
• Collapse Constraint: Tice < Tcollapse

📊 Primary Drying Results

📊 Output Summary
💾 Fortran Source

Primary Drying Cycle Duration
tprimary = 2.5 hours
Ice Front Temperature: -18.0 °C (⚠️ WARNING: Micro-Collapse Risk!)
L₀ = 9.1 mm
Ice Sublimation Vapor Pressure 935 mTorr Chamber: 150 mTorr
Average Sublimation Flux 0.336 g/(cm²·h) 3.4 kg/(m²·h)
Shelf-to-Vial Heat Flux 504.0 W/m² Kv = 28 W/(m²·K)
Safety Margin to Collapse 0.0 °C Limit: -18.0 °C

📈 Sublimation Front Depth L_dry (mm) vs Drying Time (h)

📉 Primary Drying Duration (hours) vs Shelf Temp (°C)

=================================================================
 THERMOFLUIDCALC — PHARMACEUTICAL LYOPHILIZATION REPORT
=================================================================
Case Title                 : Small Molecule Antibiotic High-Concentration Lyophile
Vial Geometry              : 5.00 mL fill in 5.50 cm2 base (Cake Depth = 9.09 mm)
Formulation Solids Content : 15.0 wt% (Collapse limit = -18.0 °C)
Chamber Vacuum / Shelf T   : 150.0 mTorr / 0.0 °C
-----------------------------------------------------------------
PRIMARY DRYING DURATION    : 2.48 hours
Ice Interface Temperature  : -18.00 °C
Cake Structural Integrity  : ⚠️ WARNING: Micro-Collapse Risk!
Ice Saturation Vapor P     : 935.2 mTorr
Average Sublimation Flux   : 0.3365 g/(cm2.h)
Shelf-to-Vial Heat Flux    : 504.00 W/m2
=================================================================

📘 Calculation Methodology & Pikal Lyophilization Standards

Pikal Coupled Heat & Mass Balance

During primary drying, latent heat of ice sublimation ($\Delta H_{sub} = 2835\,\text{J/g}$) supplied by shelf conduction must equal the vapor sublimation mass transfer:

Kv(Tshelf − Tbottom) = ΔHsub · [ (Psat(Tice) − Pch) / Rp ]

Micro-Collapse Prevention

If the ice sublimation interface exceeds the glass transition temperature $T_g'$ or collapse temperature $T_{collapse}$, the porous cake collapses, resulting in product degradation.

Key Engineering Assumptions

  • Planar 1D moving sublimation front progressing from top to bottom.
  • Vapor transport governed by dry cake Knudsen & viscous resistance $R_p$.
  • Equilibrium Clausius-Clapeyron ice vapor pressure.