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>
Go to the source code of this file.
A single-phase, multi-component model considering mineralization processes.
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 ( \varrho_f X^\kappa \phi )} {\partial t} - \text{div} \left\{ \varrho_f X^\kappa \frac{k_{r}}{\mu} \mathbf{K} (\text{grad}\, p - \varrho_{f} \mathbf{g}) \right\} - \text{div} \left\{{\bf D_{pm}^\kappa} \varrho_{f} \text{grad}\, X^\kappa \right\} - q_\kappa = 0 \qquad \kappa \in \{w, a,\cdots \} \]
The solid or mineral phases are assumed to consist of a single component. Their mass balance consist only of a storage and a source term:
\[ \frac{\partial ( \varrho_\lambda \phi_\lambda )} {\partial t} = q_\lambda \]
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}}\).
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)\) Additionally, the permeability is updated depending on the current porosity.
Namespaces | |
namespace | Dumux |
Adaption of the non-isothermal two-phase two-component flow model to problems with CO2. | |
namespace | Dumux::Properties |
namespace | Dumux::Properties::TTag |
Type tag for numeric models. | |