Definition of a problem, where CO2 is injected in a reservoir. More...
#include <test/porousmediumflow/co2/implicit/problem.hh>
Definition of a problem, where CO2 is injected in a reservoir.
The domain is sized 200m times 100m and consists of four layers, a permeable reservoir layer at the bottom, a barrier rock layer with reduced permeability, another reservoir layer and at the top a barrier rock layer with a very low permeablility.
CO2 is injected at the permeable bottom layer from the left side. The domain is initially filled with brine.
The grid is unstructered and permeability and porosity for the elements are read in from the grid file. The grid file also contains so-called boundary IDs which can be used assigned during the grid creation in order to differentiate between different parts of the boundary. These boundary ids can be imported into the problem where the boundary conditions can then be assigned accordingly.
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. The default setting for useMoles is false.
To run the simulation execute the following line in shell (works with the box and cell centered spatial discretization method): ./test_ccco2
or ./test_boxco2
Public Types | |
using | SpatialParams = GetPropType< TypeTag, Properties::SpatialParams > |
Export spatial parameter type. More... | |
Public Member Functions | |
template<class SpatialParams > | |
HeterogeneousProblem (std::shared_ptr< const GridGeometry > gridGeometry, std::shared_ptr< SpatialParams > spatialParams) | |
template<class VTKWriter > | |
void | addFieldsToWriter (VTKWriter &vtk) |
Appends all quantities of interest which can be derived from the solution of the current time step to the VTK writer. More... | |
void | setName (const std::string &newName) |
Set the problem name. More... | |
Problem parameters | |
const std::string & | name () const |
The problem name. More... | |
Scalar | temperatureAtPos (const GlobalPosition &globalPos) const |
Returns the temperature within the domain. More... | |
Boundary conditions | |
BoundaryTypes | boundaryTypes (const Element &element, const SubControlVolume &scv) const |
Specifies which kind of boundary condition should be used for which equation on a given boundary segment. More... | |
BoundaryTypes | boundaryTypes (const Element &element, const SubControlVolumeFace &scvf) const |
Specifies which kind of boundary condition should be used for which equation on a given boundary segment. More... | |
PrimaryVariables | dirichletAtPos (const GlobalPosition &globalPos) const |
Evaluates the boundary conditions for a Dirichlet boundary segment. More... | |
NumEqVector | neumann (const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const ElementFluxVariablesCache &elemFluxVarsCache, const SubControlVolumeFace &scvf) const |
Evaluates the boundary conditions for a Neumann boundary segment. More... | |
Volume terms | |
PrimaryVariables | initialAtPos (const GlobalPosition &globalPos) const |
Evaluates the initial values at a position. More... | |
Physical parameters for porous media problems | |
Scalar | temperature () const |
Returns the temperature \mathrm{[K]} for an isothermal problem. More... | |
const GravityVector & | gravityAtPos (const GlobalPosition &pos) const |
Returns the acceleration due to gravity \mathrm{[m/s^2]}. More... | |
const GravityVector & | gravity () const |
Returns the acceleration due to gravity \mathrm{[m/s^2]}. More... | |
SpatialParams & | spatialParams () |
Returns the spatial parameters object. More... | |
const SpatialParams & | spatialParams () const |
Returns the spatial parameters object. More... | |
GravityVector | gravity_ |
The gravity acceleration vector. More... | |
std::shared_ptr< SpatialParams > | spatialParams_ |
Boundary conditions and sources defining the problem | |
BoundaryTypes | boundaryTypesAtPos (const GlobalPosition &globalPos) const |
Specifies which kind of boundary condition should be used for which equation on a given boundary segment. More... | |
PrimaryVariables | dirichlet (const Element &element, const SubControlVolumeFace &scvf) const |
Evaluate the boundary conditions for a dirichlet control volume face. More... | |
PrimaryVariables | dirichlet (const Element &element, const SubControlVolume &scv) const |
Evaluate the boundary conditions for a dirichlet control volume. More... | |
NumEqVector | neumannAtPos (const GlobalPosition &globalPos) const |
Evaluate the boundary conditions for a neumann boundary segment. More... | |
NumEqVector | source (const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const SubControlVolume &scv) const |
Evaluate the source term for all phases within a given sub-control-volume. More... | |
NumEqVector | sourceAtPos (const GlobalPosition &globalPos) const |
Evaluate the source term for all phases within a given sub-control-volume. More... | |
void | addPointSources (std::vector< PointSource > &pointSources) const |
Applies a vector of point sources. The point sources are possibly solution dependent. More... | |
void | pointSource (PointSource &source, const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const SubControlVolume &scv) const |
Evaluate the point sources (added by addPointSources) for all phases within a given sub-control-volume. More... | |
void | pointSourceAtPos (PointSource &pointSource, const GlobalPosition &globalPos) const |
Evaluate the point sources (added by addPointSources) for all phases within a given sub-control-volume. More... | |
template<class MatrixBlock > | |
void | addSourceDerivatives (MatrixBlock &block, const Element &element, const FVElementGeometry &fvGeometry, const VolumeVariables &volVars, const SubControlVolume &scv) const |
Add source term derivative to the Jacobian. More... | |
NumEqVector | scvPointSources (const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const SubControlVolume &scv) const |
Adds contribution of point sources for a specific sub control volume to the values. Caution: Only overload this method in the implementation if you know what you are doing. More... | |
void | computePointSourceMap () |
Compute the point source map, i.e. which scvs have point source contributions. More... | |
const PointSourceMap & | pointSourceMap () const |
Get the point source map. It stores the point sources per scv. More... | |
void | applyInitialSolution (SolutionVector &sol) const |
Applies the initial solution for all degrees of freedom of the grid. More... | |
template<class Entity > | |
PrimaryVariables | initial (const Entity &entity) const |
Evaluate the initial value for an element (for cell-centered models) or vertex (for box / vertex-centered models) More... | |
template<class ElementSolution > | |
Scalar | extrusionFactor (const Element &element, const SubControlVolume &scv, const ElementSolution &elemSol) const |
Return how much the domain is extruded at a given sub-control volume. More... | |
Scalar | extrusionFactorAtPos (const GlobalPosition &globalPos) const |
Return how much the domain is extruded at a given position. More... | |
const GridGeometry & | fvGridGeometry () const |
The finite volume grid geometry. More... | |
const GridGeometry & | gridGeometry () const |
The finite volume grid geometry. More... | |
const std::string & | paramGroup () const |
The parameter group in which to retrieve runtime parameters. More... | |
static constexpr bool | enableInternalDirichletConstraints () |
If internal Dirichlet contraints are enabled Enables / disables internal (non-boundary) Dirichlet constraints. If this is overloaded to return true, the assembler calls problem.hasInternalDirichletConstraint(element, scv). This means you have to implement the following member function. More... | |
Implementation & | asImp_ () |
Returns the implementation of the problem (i.e. static polymorphism) More... | |
const Implementation & | asImp_ () const |
Returns the implementation of the problem (i.e. static polymorphism) More... | |
|
inherited |
Export spatial parameter type.
|
inline |
|
inline |
Appends all quantities of interest which can be derived from the solution of the current time step to the VTK writer.
|
inlineinherited |
Applies a vector of point sources. The point sources are possibly solution dependent.
pointSources | A vector of PointSource s that contain source values for all phases and space positions. |
For this method, the values method of the point source has to return the absolute rate values in units [ \textnormal{unit of conserved quantity} / s ] . Positive values mean that the conserved quantity is created, negative ones mean that it vanishes. E.g. for the mass balance that would be a mass rate in [ kg / s ] .
|
inlineinherited |
Add source term derivative to the Jacobian.
|
inlineinherited |
Applies the initial solution for all degrees of freedom of the grid.
sol | the initial solution vector |
|
inlineprotectedinherited |
Returns the implementation of the problem (i.e. static polymorphism)
|
inlineprotectedinherited |
Returns the implementation of the problem (i.e. static polymorphism)
|
inline |
Specifies which kind of boundary condition should be used for which equation on a given boundary segment.
element | The finite element |
scv | The sub-control volume |
|
inline |
Specifies which kind of boundary condition should be used for which equation on a given boundary segment.
element | The finite element |
scvf | The sub-control volume face |
|
inlineinherited |
Specifies which kind of boundary condition should be used for which equation on a given boundary segment.
globalPos | The position of the finite volume in global coordinates |
As a default, i.e. if the user's problem does not overload any boundaryTypes method set Dirichlet boundary conditions everywhere for all primary variables
|
inlineinherited |
Compute the point source map, i.e. which scvs have point source contributions.
|
inlineinherited |
Evaluate the boundary conditions for a dirichlet control volume.
element | The finite element |
scv | the sub control volume |
The method returns the boundary types information.
|
inlineinherited |
Evaluate the boundary conditions for a dirichlet control volume face.
element | The finite element |
scvf | the sub control volume face |
The method returns the boundary types information.
|
inline |
Evaluates the boundary conditions for a Dirichlet boundary segment.
globalPos | The global position |
|
inlinestaticconstexprinherited |
If internal Dirichlet contraints are enabled Enables / disables internal (non-boundary) Dirichlet constraints. If this is overloaded to return true, the assembler calls problem.hasInternalDirichletConstraint(element, scv). This means you have to implement the following member function.
bool hasInternalDirichletConstraint(const Element& element, const SubControlVolume& scv) const;
which returns a bool signifying whether the dof associated with the element/scv pair is contraint. If true is returned for a dof, the assembler calls problem.internalDiririchlet(element, scv). This means you have to additionally implement the following member function
PrimaryVariables internalDiririchlet(const Element& element, const SubControlVolume& scv) const;
which returns the enforced Dirichlet values the dof associated with the element/scv pair.
|
inlineinherited |
Return how much the domain is extruded at a given sub-control volume.
This means the factor by which a lower-dimensional (1D or 2D) entity needs to be expanded to get a full dimensional cell. The default is 1.0 which means that 1D problems are actually thought as pipes with a cross section of 1 m^2 and 2D problems are assumed to extend 1 m to the back.
|
inlineinherited |
Return how much the domain is extruded at a given position.
This means the factor by which a lower-dimensional (1D or 2D) entity needs to be expanded to get a full dimensional cell. The default is 1.0 which means that 1D problems are actually thought as pipes with a cross section of 1 m^2 and 2D problems are assumed to extend 1 m to the back.
|
inlineinherited |
The finite volume grid geometry.
|
inlineinherited |
Returns the acceleration due to gravity \mathrm{[m/s^2]}.
This method is used for problems where the gravitational acceleration does not depend on the spatial position. The default behaviour is that if the ProblemEnableGravity
property is true, \boldsymbol{g} = ( 0,\dots,\ -9.81)^T holds, else \boldsymbol{g} = ( 0,\dots, 0)^T .
|
inlineinherited |
Returns the acceleration due to gravity \mathrm{[m/s^2]}.
This is discretization independent interface. By default it just calls gravity().
|
inlineinherited |
The finite volume grid geometry.
|
inlineinherited |
Evaluate the initial value for an element (for cell-centered models) or vertex (for box / vertex-centered models)
entity | The dof entity (element or vertex) |
|
inline |
Evaluates the initial values at a position.
globalPos | The global position |
|
inline |
The problem name.
This is used as a prefix for files generated by the simulation.
|
inline |
Evaluates the boundary conditions for a Neumann boundary segment.
This is the method for the case where the Neumann condition is potentially solution dependent and requires some quantities that are specific to the fully-implicit method.
element | The finite element |
fvGeometry | The finite-volume geometry |
elemVolVars | All volume variables for the element |
elemFluxVarsCache | Flux variables caches for all faces in stencil |
scvf | The sub-control volume face |
For this method, the values parameter stores the flux in normal direction of each phase. Negative values mean influx. E.g. for the mass balance that would the mass flux in [ kg / (m^2 \cdot s)] .
|
inlineinherited |
Evaluate the boundary conditions for a neumann boundary segment.
globalPos | The position of the boundary face's integration point in global coordinates |
Negative values mean influx. E.g. for the mass balance that would be the mass flux in [ kg / (m^2 \cdot s)] .
As a default, i.e. if the user's problem does not overload any neumann method return no-flow Neumann boundary conditions at all Neumann boundaries
|
inlineinherited |
The parameter group in which to retrieve runtime parameters.
|
inlineinherited |
Evaluate the point sources (added by addPointSources) for all phases within a given sub-control-volume.
This is the method for the case where the point source is solution dependent
source | A single point source |
element | The finite element |
fvGeometry | The finite-volume geometry |
elemVolVars | All volume variables for the element |
scv | The sub control volume |
For this method, the values() method of the point sources returns the absolute conserved quantity rate generated or annihilate in units [ \textnormal{unit of conserved quantity} / s ] . Positive values mean that the conserved quantity is created, negative ones mean that it vanishes. E.g. for the mass balance that would be a mass rate in [ kg / s ] .
|
inlineinherited |
Evaluate the point sources (added by addPointSources) for all phases within a given sub-control-volume.
This is the method for the case where the point source is space dependent
pointSource | A single point source |
globalPos | The point source position in global coordinates |
For this method, the values() method of the point sources returns the absolute conserved quantity rate generated or annihilate in units [ \textnormal{unit of conserved quantity} / s ] . Positive values mean that the conserved quantity is created, negative ones mean that it vanishes. E.g. for the mass balance that would be a mass rate in [ kg / s ] .
|
inlineinherited |
Get the point source map. It stores the point sources per scv.
|
inlineinherited |
Adds contribution of point sources for a specific sub control volume to the values. Caution: Only overload this method in the implementation if you know what you are doing.
|
inlineinherited |
Set the problem name.
This static method sets the simulation name, which should be called before the application problem is declared! If not, the default name "sim" will be used.
newName | The problem's name |
|
inlineinherited |
Evaluate the source term for all phases within a given sub-control-volume.
This is the method for the case where the source term is potentially solution dependent and requires some quantities that are specific to the fully-implicit method.
element | The finite element |
fvGeometry | The finite-volume geometry |
elemVolVars | All volume variables for the element |
scv | The sub control volume |
For this method, the return parameter stores the conserved quantity rate generated or annihilate per volume unit. Positive values mean that the conserved quantity is created, negative ones mean that it vanishes. E.g. for the mass balance that would be a mass rate in [ kg / (m^3 \cdot s)] .
|
inlineinherited |
Evaluate the source term for all phases within a given sub-control-volume.
globalPos | The position of the center of the finite volume for which the source term ought to be specified in global coordinates |
For this method, the values parameter stores the conserved quantity rate generated or annihilate per volume unit. Positive values mean that the conserved quantity is created, negative ones mean that it vanishes. E.g. for the mass balance that would be a mass rate in [ kg / (m^3 \cdot s)] .
As a default, i.e. if the user's problem does not overload any source method return 0.0 (no source terms)
|
inlineinherited |
Returns the spatial parameters object.
|
inlineinherited |
Returns the spatial parameters object.
|
inlineinherited |
Returns the temperature \mathrm{[K]} for an isothermal problem.
This is not specific to the discretization. By default it just throws an exception so it must be overloaded by the problem if no energy equation is used.
|
inline |
Returns the temperature within the domain.
globalPos | The global position |
This problem assumes a geothermal gradient with a surface temperature of 10 degrees Celsius.
|
protectedinherited |
The gravity acceleration vector.
|
protectedinherited |