Local Thermal Equilibrium
By making the local thermal equilibrium hypothesis for solid and fluid phases (Equation 4-40), The Heat Transfer in Porous Media Interface solves for the following version of the heat equation (Ref. 17), reformulated using a common temperature, T:
(4-41)
(4-42)
The different quantities appearing here are:
ρf is the fluid density.
Cp,f is the fluid heat capacity at constant pressure.
Cp)eff is the effective volumetric heat capacity at constant pressure, defined by
εp is the porosity.
θs is the solid matrix volume fraction.
ρs is the solid matrix density.
Cp,s is the solid matrix heat capacity at constant pressure.
keff is the effective thermal conductivity (a scalar or a tensor if the thermal conductivity is anisotropic).
q is the conductive heat flux.
u is the velocity field, either an analytic expression or computed from a Fluid Flow interface. It should be interpreted as the Darcy velocity, that is, the volume flow rate per unit cross sectional area. The average linear velocity (the velocity within the pores) can be calculated as uf = up, where εp is the fluid’s volume fraction, or equivalently the porosity.
Q is the heat source (or sink). Add one or several heat sources as separate physics features.
For a steady-state problem the temperature does not change with time, and the terms with time derivatives of Equation 4-41 disappear.
The effective thermal conductivity of the solid-fluid system, keff, is related to the conductivity of the solid, ks, and to the conductivity of the fluid, kf, and depends in a complex way on the geometry of the medium. In Ref. 13, Ref. 14, and Ref. 18, the following models are proposed:
This volume average model provides an upper bound for the effective thermal conductivity.
This reciprocal average model provides a lower bound for the effective thermal conductivity.
A good estimate is given by the weighted geometric mean of kf and ks, as long as kf and ks are not too different from each other
Note that when kf and ks are equal, all these models give the same effective thermal conductivity. In addition, the effective conductivity is equal to kf when the porosity is 1, and to ks when the porosity is 0.