1pnc

Single-phase, multi-component flow. More...

Description

Adaption of the fully implicit scheme to the one-phase n-component pore network model.

A mass balance equation is formulated for each pore body $i$ and each component $\kappa$:

$$V_i\frac{\partial ({\varrho }_i{X}^{\kappa })}{\partial t}+\sum _j(\varrho X^{\kappa }Q{)}_{ij}=(V{q}^{\kappa }{)}_i.$$

$V_i$ is the pore body volume, and the advective mass flow $(\varrho Q{X}^{\kappa }{)}_{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_i$ is a mass sink or source term defined on pore body $i$.

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

$$Q_{ij}=g_{ij}(p_i-p_j+\mathrm{\Psi }).$$

$g_{ij}$ is a suitable throat conductance value (see e.g. PoreNetwork::TransmissibilityPatzekSilin) while $p_i$ and $p_j$ are averaged pore body pressures.

The (optional) influence of gravity is given by

$$\mathrm{\Psi }=\varrho \mathbf{g}({\mathbf{x}}_{\mathbf{i}}-{\mathbf{x}}_{\mathbf{j}}),$$

where ${\mathbf{x}}_{\mathbf{i}}-{\mathbf{x}}_{\mathbf{j}}$ is the distance vector between the centers of pore bodies $i$ and $j$ and $\mathbf{g}$ is the gravitational acceleration.

The model is able to use either mole or mass fractions. The property useMoles can be set to either true or false in the problem file. Make sure that the according units are used in the problem setup. useMoles is set to true by default.

The primary variables are the pressure $p$ and the mass or mole fraction of dissolved components $X^{\kappa }$ or $x^{\kappa }$.

Files
file	porenetwork/1pnc/iofields.hh
file	porenetwork/1pnc/model.hh
	Adaption of the fully implicit scheme to the one-phase n-component pore network model.
file	porenetwork/1pnc/volumevariables.hh

Classes
class	Dumux::PoreNetwork::OnePNCIOFields< FluidSystem >
	Adds output fields specific to the PNM 1pnc model. More...
class	Dumux::PoreNetwork::OnePNCVolumeVariables< Traits >
	Contains the quantities which are are constant within a finite volume in the one-phase, n-component model. More...