🪨 Porous Media Convection (Darcy-Forchheimer)

Compute convective heat transfer in saturated porous media: porous Nusselt number, Darcy-Forchheimer non-linear pressure drop, filtration pumping power, and porous Rayleigh Ra_K.

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

Convection

📊 Solver Telemetry ● ACTIVE

👁️ Views 37

⚡ Solves 30

💾 Downloads 470 📦 Fortran Code 10.8 KB

📅 Released Jun 2026

⏱️ Latency < 1 ms

💡 Hardware & Sizing Partner: Need to size a control valve or check ASME flange bolt torque for this line?

FLOWKS IEC 60534 Sizer ↗ • API 6D Valve Catalog →

🪨 Saturated Porous Matrix & Interstitial Seepage Flow

Real-time visual simulation of fluid filtration through packed grains, Darcy pressure drop & thermal dispersion

📝 Configuration & Presets

♨️ Geothermal Aquifer 🧪 Catalytic Packed Bed 🧽 Metal Foam Heat Sink 🧱 Wall Fibrous Insulation

📐 Bed Geometry & Porous Structure

Bed Flow Length H [m]:

Bed Cross Width W [m]:

Permeability K [m²]: 1 Darcy ≈ 9.87×10⁻¹³ m²

Porosity φ (0 to 1):

⚡ Flow Velocity & Temperatures

Superficial Darcy Velocity U_D [m/s]:

Bed Depth L [m]:

Hot Boundary Temp T_h [°C]:

Cold Boundary Temp T_c [°C]:

🧪 Solid & Fluid Thermal Properties

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

Fluid Conductivity k_f [W/(m·K)]:

Darcy-Forchheimer Convection Formulation:
• Effective Conductivity: k_e = φ k_f + (1 − φ) k_s
• Pressure Drop: ΔP = [ (μ/K) U_D + ρ (C_F/√K) U_D² ] · H
• Porous Nusselt: Nu = √(1 + 0.318 Pe_K) | Pe_K = U_D H / α_e
• Porous Rayleigh: Ra_K = (ρ g β ΔT K H) / (μ α_e)

📊 Porous Convection Results

📊 Output Summary
💾 Fortran Source

Total Heat Transfer Rate (Q̇)

Q̇ = 163.4 W

Effective HTC: 5.0 W/(m²·K) | Nu = 8.47

k_e = 0.71 W/m·K

Forchheimer Pressure Drop (ΔP) 623.1 Pa Hydraulic resistance

Pumping Power (P_pump) 12.46 W Mechanical energy required

Porous Rayleigh (Ra_K) 1.99e+0 Peclet Pe_K = 222.2

Effective Conductivity (k_e) 0.713 W/(m·K) Volumetric mixture

📈 Porous Nusselt Nu vs Darcy Velocity U_D (m/s)

📉 Forchheimer Pressure Drop ΔP vs Velocity U_D (m/s)

=================================================================
 THERMOFLUIDCALC — POROUS MEDIA CONVECTION (DARCY-FORCHHEIMER) REPORT
=================================================================
Case Title                 : Catalytic Packed Bed Gas Phase Chemical Reactor
Porous Bed Dimensions      : Height H = 1.20 m, Width W = 0.50 m, Depth L = 0.50 m
Porous Matrix Properties   : Permeability K = 5.000e-9 m2, Porosity phi = 0.42
Effective Conductivity     : ke_eff = 0.713 W/(m.K) (ks = 1.20, kf = 0.04)
Filtration Conditions      : Darcy Velocity UD = 0.0800 m/s, DeltaT = 130.0 C
-----------------------------------------------------------------
POROUS NUSSELT NUMBER (Nu) : 8.466
Effective HTC (h_porous)   : 5.03 W/(m2.K)
TOTAL HEAT TRANSFER RATE Q : 163.43 W (0.163 kW)
Forchheimer Pressure Drop  : 623.12 Pa (0.623 kPa)
Filtration Pumping Power   : 12.462 W
=================================================================

📘 Calculation Methodology & Porous Convection Standards

Darcy-Forchheimer Model

Combines linear viscous drag (Darcy law) and non-linear quadratic form drag (Forchheimer inertial term):

                    ∇P = − (μ/K) UD − ρ (CF/√K) UD²
                

Effective Thermal Conductivity

Represents the parallel/series mixture volume average of solid matrix and interstitial saturated fluid ($k_e = \phi k_f + (1-\phi) k_s$).

Key Engineering Assumptions

Homogeneous, isotropic porous matrix.
Local thermal equilibrium between fluid and solid phases ($T_f = T_s$).
Laminar and transitional filtration velocities ($Re_K = \frac{U_D \sqrt{K}}{\nu} < 10$).

ThermoFluidCalc Engineering Report