porousmediumflow/richardsextended/model.hh File Reference

This model implements a variant of the extended Richards' equation for quasi-twophase flow (see e.g. Vanderborght et al. 2017). More...

#include <dune/common/fvector.hh>
#include <dumux/common/properties.hh>
#include <dumux/porousmediumflow/immiscible/localresidual.hh>
#include <dumux/porousmediumflow/compositional/switchableprimaryvariables.hh>
#include <dumux/material/fluidmatrixinteractions/diffusivitymillingtonquirk.hh>
#include <dumux/material/fluidmatrixinteractions/2p/thermalconductivity/somerton.hh>
#include <dumux/material/components/simpleh2o.hh>
#include <dumux/material/fluidsystems/h2oair.hh>
#include <dumux/material/fluidstates/immiscible.hh>
#include <dumux/porousmediumflow/properties.hh>
#include <dumux/porousmediumflow/nonisothermal/model.hh>
#include <dumux/porousmediumflow/nonisothermal/indices.hh>
#include <dumux/porousmediumflow/nonisothermal/iofields.hh>
#include <dumux/porousmediumflow/richards/balanceequationopts.hh>
#include <dumux/porousmediumflow/richards/model.hh>
#include <dumux/porousmediumflow/richards/velocityoutput.hh>
#include "indices.hh"
#include "volumevariables.hh"
#include "iofields.hh"
#include "localresidual.hh"

Go to the source code of this file.

Description

The extended Richards' equation

$$\frac{\partial (\varphi S_w{\varrho }_w)}{\partial t}+\frac{\partial (\varphi (1-S_w){\varrho }_n{X}_n^w)}{\partial t}-\mathrm{\nabla }\cdot \{{\varrho }_w\frac{k_{rw}}{{\mu }_w}\mathbf{K}(\text{nabla}{p}_w-{\varrho }_w\mathbf{\text{g}})+{\mathbf{D}}_{\mathbf{n}\mathbf{,}\mathbf{p}\mathbf{m}}^{\mathbf{w}}{\varrho }_n\mathrm{\nabla }{X}_n^w\}=q_w,$$

where:

• $\varphi$ is the porosity of the porous medium,
• $S_w$ represents the saturation of the wetting-phase,
• ${\varrho }_w$ is the mass density of the wetting phase,
• ${\varrho }_n$ is the mass density of the non-wetting phase,
• $k_{rw}$ is the relative permeability of the wetting phase,
• ${\mu }_w$ is the dynamic viscosity of the wetting phase,
• $\mathbf{K}$ is the intrinsic permeability tensor,
• $p_w$ is the pressure of the wetting phase,
• $\mathbf{g}$ is the gravitational acceleration vector,
• ${\mathbf{D}}_{\mathbf{n}\mathbf{,}\mathbf{p}\mathbf{m}}^{\mathbf{w}}$ is the effective diffusivity of water in the non-wetting phase,
• $X_n^w$ is the mass fraction of water in the non-wetting phase,
• $q_w$ is a source or sink term in the wetting phase,

additionally models water vapor diffusion in the gas phase. The model is derived based on the two-phase flow equations based on the assumption that the gas phase does not move but and remains at constant pressure.

Classes
struct	Dumux::ExtendedRichardsModelTraits
	Specifies a number properties of the extended Richards model. More...
struct	Dumux::ExtendedRichardsVolumeVariablesTraits< PV, FSY, FST, SSY, SST, PT, MT, DT, EDM >
	Traits class for the Richards model. More...
struct	Dumux::Properties::TTag::ExtendedRichards
struct	Dumux::Properties::TTag::ExtendedRichardsNI
struct	Dumux::Properties::LocalResidual< TypeTag, TTag::ExtendedRichards >
	The local residual operator. More...
struct	Dumux::Properties::IOFields< TypeTag, TTag::ExtendedRichards >
	Set the vtk output fields specific to this model. More...
struct	Dumux::Properties::ModelTraits< TypeTag, TTag::ExtendedRichards >
	The model traits. More...
struct	Dumux::Properties::VolumeVariables< TypeTag, TTag::ExtendedRichards >
	Set the volume variables property. More...
struct	Dumux::Properties::EffectiveDiffusivityModel< TypeTag, TTag::ExtendedRichards >
	Use the model after Millington (1961) for the effective diffusivity. More...
struct	Dumux::Properties::PrimaryVariables< TypeTag, TTag::ExtendedRichards >
	The primary variables vector for the richards model. More...
struct	Dumux::Properties::ModelTraits< TypeTag, TTag::ExtendedRichardsNI >
	set non-isothermal model traits More...
struct	Dumux::Properties::IOFields< TypeTag, TTag::ExtendedRichardsNI >
	Set the vtk output fields specific to th non-isothermal model. More...
struct	Dumux::Properties::VolumeVariables< TypeTag, TTag::ExtendedRichardsNI >
	Set the volume variables property. More...
struct	Dumux::Properties::ThermalConductivityModel< TypeTag, TTag::ExtendedRichardsNI >
	Somerton is used as default model to compute the effective thermal heat conductivity. More...

Namespaces
namespace	Dumux
namespace	Dumux::Properties
	The energy balance equation for a porous solid.
namespace	Dumux::Properties::TTag
	Type tag for numeric models.

Include dependency graph for porousmediumflow/richardsextended/model.hh: