porenetwork/2p/model.hh File Reference

A two-phase-flow, isothermal pore-network model using the fully implicit scheme. More...

#include <dumux/common/properties.hh>
#include <dumux/flux/porenetwork/advection.hh>
#include <dumux/porenetwork/properties.hh>
#include <dumux/material/fluidmatrixinteractions/porenetwork/throat/transmissibility1p.hh>
#include <dumux/material/fluidmatrixinteractions/porenetwork/throat/transmissibility2p.hh>
#include <dumux/material/fluidmatrixinteractions/porenetwork/pore/2p/multishapelocalrules.hh>
#include <dumux/porousmediumflow/immiscible/localresidual.hh>
#include <dumux/porousmediumflow/nonisothermal/model.hh>
#include <dumux/porousmediumflow/nonisothermal/indices.hh>
#include <dumux/porousmediumflow/nonisothermal/iofields.hh>
#include <dumux/porousmediumflow/2p/model.hh>
#include "fluxvariablescache.hh"
#include "gridfluxvariablescache.hh"
#include "iofields.hh"
#include "volumevariables.hh"
#include "spatialparams.hh"

Go to the source code of this file.

Description

A two-phase-flow, isothermal pore-network model using the fully implicit scheme.

A mass balance equation is formulated for each pore body \(i\) and each phase \(\alpha\):

\[ V_i \frac{\partial (\varrho_\alpha S_\alpha)_i}{\partial t} + \sum_j (\varrho_\alpha Q_\alpha)_{ij} = (V q_\alpha)_i ~. \]

\(V_i\) is the pore body volume, and the advective mass flow \((\varrho_\alpha Q_\alpha)_{ij}\) through throat \(ij\) can be based on the fluid phase density \(\varrho\) either of the upstream pore body \(i\) or \(j\) (upwinding) or on the respective averaged value. \(q_\alpha\) is a mass sink or source term defined on pore body \(i\).

Per default, the volume flow rate \(Q_{\alpha,ij}\) follows a linear Hagen-Poiseuille-type law (PoreNetworkModel::CreepingFlow) which is only valid for \(Re < 1\):

\[ Q_{\alpha,ij} = g_{\alpha, ij} (p_{\alpha, i} - p_{\alpha, j} + \Psi_\alpha) ~. \]

\(g_{\alpha,ij}\) is a suitable throat conductance value that takes into account the presence/saturation of the individual phases while \(p_{\alpha,i}\) and \(p_{\alpha,j}\) are averaged pore body phase pressures.

The (optional) influence of gravity is given by

\[ \Psi_\alpha = \varrho_\alpha \mathbf{g} (\mathbf{x_i} - \mathbf{x_j}) ~, \]

where \(\mathbf{x_i} - \mathbf{x_j}\) is the distance vector between the centers of pore bodies \(i\) and \(j\) and \(\mathbf{g}\) is the gravitational acceleration.

The primary variables are the wetting phase pressure and the the nonwetting phase saturation ( \(p_w\) and \(S_n\)) or the nonwetting phase pressure and the the wetting phase saturation ( \(p_n\) and \(S_w\)), depending on the chose formulation (see TwoPModel).

Classes
struct	Dumux::Properties::TTag::PNMTwoP
struct	Dumux::Properties::TTag::PNMTwoPNI
	The type tags for the corresponding non-isothermal problems. More...
struct	Dumux::Properties::VolumeVariables< TypeTag, TTag::PNMTwoP >
	Set the volume variables property. More...
struct	Dumux::Properties::FluxVariablesCache< TypeTag, TTag::PNMTwoP >
	The flux variables cache. More...
struct	Dumux::Properties::GridFluxVariablesCache< TypeTag, TTag::PNMTwoP >
	The grid flux variables cache vector class. More...
struct	Dumux::Properties::SpatialParams< TypeTag, TTag::PNMTwoP >
struct	Dumux::Properties::AdvectionType< TypeTag, TTag::PNMTwoP >
	The advection type. More...
struct	Dumux::Properties::EnergyLocalResidual< TypeTag, TTag::PNMTwoP >
struct	Dumux::Properties::IOFields< TypeTag, TTag::PNMTwoP >
struct	Dumux::Properties::ModelTraits< TypeTag, TTag::PNMTwoPNI >
	model traits of the non-isothermal model. More...
struct	Dumux::Properties::VolumeVariables< TypeTag, TTag::PNMTwoPNI >
	Set the volume variables property. More...
struct	Dumux::Properties::IOFields< TypeTag, TTag::PNMTwoPNI >
	Set the vtk output fields specific to the non-isothermal two-phase model. More...
struct	Dumux::Properties::ThermalConductivityModel< TypeTag, TTag::PNMTwoPNI >
	Use the average for effective conductivities. More...

Namespaces
namespace	Dumux
namespace	Dumux::Properties
namespace	Dumux::Properties::TTag
	Type tag for numeric models.

Include dependency graph for porenetwork/2p/model.hh: