🪨 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

📊 Solver Telemetry ● ACTIVE

👁️ Consultations 37

⚡ Calculs faits 30

💾 Téléchargements 358 📦 Code Fortran 10.8 KB

📅 Mise en service Jun 2026

⏱️ Latence < 1 ms

🪨 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̇ = 0.1 W

Effective HTC: 0.0 W/(m²·K) | Nu = 2.23

k_e = 0.03 W/m·K

Forchheimer Pressure Drop (ΔP) 4.4 Pa Hydraulic resistance

Pumping Power (P_pump) 0.00 W Mechanical energy required

Porous Rayleigh (Ra_K) 1.03e+1 Peclet Pe_K = 12.5

Effective Conductivity (k_e) 0.026 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                 : Building Wall Fibrous Insulation Air Infiltration
Porous Bed Dimensions      : Height H = 2.50 m, Width W = 1.00 m, Depth L = 0.15 m
Porous Matrix Properties   : Permeability K = 1.000e-9 m2, Porosity phi = 0.95
Effective Conductivity     : ke_eff = 0.026 W/(m.K) (ks = 0.04, kf = 0.03)
Filtration Conditions      : Darcy Velocity UD = 0.0001 m/s, DeltaT = 32.0 C
-----------------------------------------------------------------
POROUS NUSSELT NUMBER (Nu) : 2.230
Effective HTC (h_porous)   : 0.02 W/(m2.K)
TOTAL HEAT TRANSFER RATE Q : 0.11 W (0.000 kW)
Forchheimer Pressure Drop  : 4.38 Pa (0.004 kPa)
Filtration Pumping Power   : 0.000 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