A two-phase one-component flow model using the fully implicit scheme. More...

#include <dune/common/fvector.hh>
#include <dumux/common/properties.hh>
#include <dumux/material/fluidmatrixinteractions/2p/thermalconductivity/somerton.hh>
#include <dumux/material/fluidstates/compositional.hh>
#include <dumux/porousmediumflow/properties.hh>
#include <dumux/porousmediumflow/2p/formulation.hh>
#include <dumux/porousmediumflow/compositional/switchableprimaryvariables.hh>
#include <dumux/porousmediumflow/nonisothermal/model.hh>
#include <dumux/porousmediumflow/nonisothermal/indices.hh>
#include <dumux/porousmediumflow/nonisothermal/iofields.hh>
#include "darcyslaw.hh"
#include "iofields.hh"
#include "localresidual.hh"
#include "indices.hh"
#include "volumevariables.hh"

Description

Note: The 2p1c model requires the use of the non-isothermal extension found in dumux/porousmediumflow/nonisothermal.

This model is designed for simulating two fluid phases with water as the only component. It is particularly suitable for the simulation of steam injection in saturated conditions.

The model implements the flow of two phases and one component, i.e. a pure liquid (e.g. water) and its vapor (e.g. steam), \(\alpha \in \{ w, n \}\) using a standard multi-phase Darcy approach as the equation for the conservation of momentum. For details on Darcy's law see dumux/flux/darcyslaw.hh.

By inserting Darcy's law into the equations for the conservation of the phase mass, one gets

\[ \frac{\partial\ \sum_\alpha (\phi \rho_\alpha S_\alpha)}{\partial t} \\-\sum \limits_ \alpha \nabla \cdot \left \{\rho_\alpha \frac{k_{r\alpha}}{\mu_\alpha} \mathbf{K} (\nabla p_\alpha - \rho_\alpha \mathbf{g}) \right \} -q^w =0, \]

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 \),
\( 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^w \) is a source or sink term.

By using constitutive relations for the capillary pressure \(p_c = p_n - p_w\) and relative permeability \(k_{r\alpha}\) and taking advantage of the fact that \(S_w + S_n = 1\), the number of unknowns can be reduced to two. The model features a primary variable switch. If only one phase is present, \(p_n\) and \(T\) are the primary variables. In the presence of two phases, \(p_n\) and \(S_w\) become the new primary variables.

Classes
struct	Dumux::TwoPOneCNIModelTraits< f >
	Specifies a number properties of models considering two phases with water as a single component. More...

struct	Dumux::TwoPOneCVolumeVariablesTraits< PV, FSY, FST, SSY, SST, PT, MT >
	Traits class for the two-phase model. More...

struct	Dumux::Properties::TTag::TwoPOneCNI

struct	Dumux::Properties::FluidState< TypeTag, TTag::TwoPOneCNI >
	The fluid state which is used by the volume variables to store the thermodynamic state. This should be chosen appropriately for the model ((non-)isothermal, equilibrium, ...). This can be done in the problem. More...

struct	Dumux::Properties::Formulation< TypeTag, TTag::TwoPOneCNI >
	Set the default formulation to pw-sn. More...

struct	Dumux::Properties::UseBlockingOfSpuriousFlow< TypeTag, TTag::TwoPOneCNI >
	Do not block spurious flows by default. More...

struct	Dumux::Properties::LocalResidual< TypeTag, TTag::TwoPOneCNI >
	The specific local residual (i.e. balance equations). More...

struct	Dumux::Properties::AdvectionType< TypeTag, TTag::TwoPOneCNI >
	Use a modified version of Darcy's law which allows for blocking of spurious flows. More...

struct	Dumux::Properties::VolumeVariables< TypeTag, TTag::TwoPOneCNI >
	Set the volume variables property. More...

struct	Dumux::Properties::PrimaryVariables< TypeTag, TTag::TwoPOneCNI >
	The primary variables vector for the 2p1cni model. More...

struct	Dumux::Properties::ThermalConductivityModel< TypeTag, TTag::TwoPOneCNI >
	Somerton is used as default model to compute the effective thermal heat conductivity. More...

struct	Dumux::Properties::ModelTraits< TypeTag, TTag::TwoPOneCNI >
	Set the non-isothermal model traits. More...

struct	Dumux::Properties::IOFields< TypeTag, TTag::TwoPOneCNI >
	The non-isothermal vtk output fields. 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/2p1c/model.hh:

Description

Classes

Namespaces