🪨 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

Configure inputs and click Compute to view results.

📘 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