A single-phase, multi-component model considering mineralization processes. More...

#include <dumux/porousmediumflow/1pnc/model.hh>
#include <dumux/porousmediumflow/1pnc/indices.hh>
#include <dumux/porousmediumflow/1pnc/volumevariables.hh>
#include <dumux/material/solidstates/compositionalsolidstate.hh>
#include <dumux/porousmediumflow/mineralization/model.hh>
#include <dumux/porousmediumflow/mineralization/localresidual.hh>
#include <dumux/porousmediumflow/mineralization/volumevariables.hh>
#include <dumux/porousmediumflow/mineralization/iofields.hh>
#include <dumux/porousmediumflow/nonisothermal/indices.hh>
#include <dumux/porousmediumflow/nonisothermal/iofields.hh>
#include <dumux/material/fluidmatrixinteractions/1p/thermalconductivityaverage.hh>

Description

This model implements one-phase n-component flow of a compressible fluid composed of the n components \(\kappa \) in combination with mineral precipitation and dissolution of the solid phases. The standard multi-phase Darcy approach is used 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 components, one gets one transport equation for each component,

\[ \frac{\partial ( \phi \varrho_f X^\kappa )} {\partial t} - \nabla \cdot \left\{ \varrho_f X^\kappa \frac{k_{r}}{\mu} \mathbf{K} (\nabla p - \varrho_{f} \mathbf{g}) \right\} - \nabla \cdot \left\{{\bf D_{pm}^\kappa} \varrho_{f} \nabla X^\kappa \right\} - q_\kappa = 0 \qquad \kappa \in \{w, a,\cdots \}, \]

where:

\( \phi \) is the porosity,
\( \varrho_f \) is the mass density of the fluid,
\( X^\kappa \) is the mass fraction of component \( \kappa \) in the fluid,
\( k_{r} \) is the relative permeability,
\( \mu \) represents the dynamic viscosity,
\( \mathbf{K} \) is the intrinsic permeability tensor,
\( p \) is the pressure,
\( \mathbf{g} \) is the gravitational acceleration vector,
\( {\bf D_{pm}^\kappa} \) is the effective diffusivity in the porous medium,
\( q_\kappa \) is a source or sink term.

The solid or mineral phases are assumed to consist of a single component. Their mass balance consists of only a storage and a source term,

\[ \frac{\partial ( \varrho_\lambda \phi_\lambda )} {\partial t} = q_\lambda, \]

where:

\( \varrho_\lambda\) represents the mass density of the solid phase,
\( \phi_\lambda \) is the porosity of the solid,
\( q_\lambda \) is a source or sink term representing.

The primary variables are the pressure \(p\) and the mole fractions of the dissolved components \(x^k\). The primary variable of the solid phases is the volume fraction \(\phi_\lambda = \frac{V_\lambda}{V_{total}}\),

where:

\( V_\lambda \) is the volume of phase \( \lambda \),
\( V_{\text{total}} \) is the total volume of the system.

The source an sink terms link the mass balances of the n-transported component to the solid phases. The porosity \(\phi\) is updated according to the reduction of the initial (or solid-phase-free porous medium) porosity \(\phi_0\) by the accumulated volume fractions of the solid phases, \( \phi = \phi_0 - \sum (\phi_\lambda),\)

where:

\( \phi \) represents the remaining porosity in the system,
\( \phi_0 \) is the initial porosity,
\( \phi_\lambda \) denotes the volume fraction of phase \( \lambda \). Additionally, the permeability is updated depending on the current porosity.

Classes
struct	Dumux::Properties::TTag::OnePNCMin

struct	Dumux::Properties::TTag::OnePNCMinNI

struct	Dumux::Properties::VolumeVariables< TypeTag, TTag::OnePNCMin >
	use the mineralization volume variables together with the 1pnc vol vars More...

struct	Dumux::Properties::LocalResidual< TypeTag, TTag::OnePNCMin >

struct	Dumux::Properties::ModelTraits< TypeTag, TTag::OnePNCMin >
	Use non-mineralization model traits with 1pnc traits. More...

struct	Dumux::Properties::SolidState< TypeTag, TTag::OnePNCMin >
	The two-phase model uses the immiscible fluid state. More...

struct	Dumux::Properties::IOFields< TypeTag, TTag::OnePNCMin >
	Use the mineralization vtk output fields. More...

struct	Dumux::Properties::IOFields< TypeTag, TTag::OnePNCMinNI >
	Non-isothermal vtk output. More...

struct	Dumux::Properties::ModelTraits< TypeTag, TTag::OnePNCMinNI >
	The non-isothermal model traits. More...

struct	Dumux::Properties::VolumeVariables< TypeTag, TTag::OnePNCMinNI >
	use the mineralization volume variables together with the 1pnc vol vars More...

struct	Dumux::Properties::ThermalConductivityModel< TypeTag, TTag::OnePNCMinNI >
	Use the average for effective conductivities. 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/1pncmin/model.hh:

Description

Classes

Namespaces