version 3.9-dev
porousmediumflow/nonisothermal/model.hh File Reference

The implicit non-isothermal model. More...

Go to the source code of this file.

Description

This model implements a generic energy balance for single and multi-phase transport problems. Currently the non-isothermal model can be used on top of the 1p2c, 2p, 2p2c and 3p3c models. Comparison to simple analytical solutions for pure convective and conductive problems are found in the 1p2c test. Also refer to this test for details on how to activate the non-isothermal model.

For the energy balance, local thermal equilibrium is assumed. This results in one energy conservation equation for the porous solid matrix and the fluids,

\begin{align*} \frac{\partial (\sum_\alpha \phi \varrho_\alpha u_\alpha S_\alpha )}{\partial t} & + \frac{\partial \left((\left( 1 - \phi \right) \varrho_s c_s T\right)}{\partial t} - \sum_\alpha \nabla \cdot \left\{ \varrho_\alpha h_\alpha \frac{k_{r\alpha}}{\mu_\alpha} \mathbf{K} \left( \nabla p_\alpha - \varrho_\alpha \mathbf{g} \right) \right\} \\ & - \nabla \cdot \left(\lambda_{pm} \nabla T \right) - q^h = 0, \end{align*}

where:

  • \( \phi \) is the porosity of the porous medium,
  • \( S_\alpha \) represents the saturation of phase \( \alpha \),
  • \( \rho_\alpha \) is the mass density of phase \( \alpha \),
  • \( h_\alpha \) is the specific enthalpy of phase \( \alpha \),
  • \( u_\alpha \) is the specific internal energy of phase \( \alpha \),
  • \( \lambda_{pm}\) is the effective heat conductivity in the porous medium,
  • \( T \) is the temperature,
  • \( \rho_s \) is the mass density of the solid phase,
  • \( c_s \) is the heat capacity of the solid,
  • \( k_{r\alpha} \) is the relative permeability of phase \( \alpha \),
  • \( \mu_\alpha \) is the dynamic viscosity of phase \( \alpha \),
  • \( \mathbf{K} \) is the intrinsic permeability tensor,
  • \( p_\alpha \) is the pressure of phase \( \alpha \),
  • \( \mathbf{g} \) is the gravitational acceleration vector,
  • \( q^h \) is a source or sink term.

Classes

struct  Dumux::PorousMediumFlowNIModelTraits< IsothermalT, TDM >
 Specifies a number properties of non-isothermal porous medium flow models based on the specifics of a given isothermal model. More...
 

Namespaces

namespace  Dumux
 
Include dependency graph for porousmediumflow/nonisothermal/model.hh: