Adaption of the fully implicit model to the one-phase n-component flow model. More...
#include <dumux/common/properties.hh>
#include <dumux/material/fluidmatrixinteractions/1p/thermalconductivityaverage.hh>
#include <dumux/material/fluidstates/compositional.hh>
#include <dumux/porousmediumflow/properties.hh>
#include <dumux/porousmediumflow/1p/model.hh>
#include <dumux/porousmediumflow/compositional/localresidual.hh>
#include <dumux/porousmediumflow/nonisothermal/model.hh>
#include <dumux/porousmediumflow/nonisothermal/indices.hh>
#include <dumux/porousmediumflow/nonisothermal/iofields.hh>
#include <dumux/porousmediumflow/nonequilibrium/model.hh>
#include <dumux/porousmediumflow/nonequilibrium/volumevariables.hh>
#include <dumux/material/fluidmatrixinteractions/diffusivitymillingtonquirk.hh>
#include "indices.hh"
#include "volumevariables.hh"
#include "iofields.hh"
Go to the source code of this file.
Adaption of the fully implicit model to the one-phase n-component flow model.
This model implements a one-phase flow of a compressible fluid, that consists of n components, using a standard Darcy approach as the equation for the conservation of momentum. For details on Darcy's law see dumux/flux/darcyslaw.hh. Gravity can be enabled or disabled via the property system. By inserting Darcy's law into the continuity equation, one gets
\[ \phi\frac{\partial \varrho}{\partial t} - \text{div} \left\{ \varrho \frac{\textbf K}{\mu} \left(\textbf{grad}\, p - \varrho {\textbf g} \right) \right\} = q \;, \]
The transport of the components \(\kappa \in \{ w, a, ... \}\) is described by the following equation:
\[ \phi \frac{ \partial \varrho X^\kappa}{\partial t} - \text{div} \left\lbrace \varrho X^\kappa \frac{{\textbf K}}{\mu} \left( \textbf{grad}\, p - \varrho {\textbf g} \right) + \varrho D^\kappa_\text{pm} \textbf{grad} X^\kappa \right\rbrace = q. \]
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 mole fraction of dissolved components \(x^\kappa\).
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. | |