 |
3.1-git
DUNE for Multi-{Phase, Component, Scale, Physics, ...} flow and transport in porous media
|
| |
3.1-git
DUNE for Multi-{Phase, Component, Scale, Physics, ...} flow and transport in porous media
|
Definition of a problem, for the 1pnc problem. More...
#include <test/porousmediumflow/1pnc/implicit/1p2c/nonisothermal/conduction/problem.hh>
Definition of a problem, for the 1pnc problem.
Nitrogen is dissolved in the water phase and is transported with the water flow from the left side to the right.
The model domain is 1m times 1m with a discretization length of 0.05m and homogeneous soil properties ( \( \mathrm{K=10e-10, \Phi=0.4, \tau=0.28}\)). Initially, the domain is filled with pure water.
At the left side, a Dirichlet condition defines a nitrogen mole fraction of 0.3mol/mol. The water phase flows from the left side to the right due to the applied pressure gradient of 1e5Pa/m. The nitrogen is transported with the water flow and leaves the domain at the right boundary, where again Dirichlet boundary conditions are applied. Here, the nitrogen mole fraction is set to 0.0mol/mol.
This problem uses the 1pnc model.
To run the simulation execute the following line in shell: ./test_box1pnc -parameterFile ./test_box1pnc.input or ./test_cc1pnc -parameterFile ./test_cc1pnc.input
Public Types | |
| using | SpatialParams = GetPropType< TypeTag, Properties::SpatialParams > |
| Export spatial parameter type. More... | |
Public Member Functions | |
| OnePTwoCNIConductionProblem (std::shared_ptr< const GridGeometry > gridGeometry) | |
| const std::vector< Scalar > & | getExactTemperature () |
| Get the analytical temperature. More... | |
| void | updateExactTemperature (const SolutionVector &curSol, Scalar time) |
| Udpate the analytical temperature. More... | |
| const std::string & | name () const |
| The problem name. More... | |
| void | setName (const std::string &newName) |
| Set the problem name. More... | |
Problem parameters | |
| Scalar | temperature () const |
| Returns the temperature within the domain [K]. More... | |
Boundary conditions | |
| 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 | dirichletAtPos (const GlobalPosition &globalPos) const |
| Evaluates the boundary conditions for a Dirichlet boundary segment. More... | |
| NumEqVector | neumannAtPos (const GlobalPosition &globalPos) const |
| Evaluates the boundary conditions for a Neumann boundary segment. More... | |
Volume terms | |
| NumEqVector | sourceAtPos (const GlobalPosition &globalPos) const |
| Evaluates the source term for all phases within a given sub-control volume. More... | |
| PrimaryVariables | initialAtPos (const GlobalPosition &globalPos) const |
| Evaluates the initial value for a control volume. More... | |
Physical parameters for porous media problems | |
| Scalar | temperatureAtPos (const GlobalPosition &globalPos) const |
| Returns the temperature \(\mathrm{[K]}\) at a given global position. 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 | 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 | 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 | neumann (const Element &element, const FVElementGeometry &fvGeometry, const ElementVolumeVariables &elemVolVars, const ElementFluxVariablesCache &elemFluxVarsCache, const SubControlVolumeFace &scvf) 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... | |
| 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 |
|
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)
|
inlineinherited |
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 |
|
inlineinherited |
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 |
|
inline |
Specifies which kind of boundary condition should be used for which equation on a given boundary segment.
| globalPos | The position for which the bc type should be evaluated |
|
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 position for which the bc type should be evaluated |
|
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.
|
inline |
Get the analytical temperature.
|
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 value for a control volume.
| globalPos | The position for which the initial condition should be evaluated |
For this method, the values parameter stores primary variables.
|
inlineinherited |
The problem name.
This is used as a prefix for files generated by the simulation. It could be either overwritten by the problem files, or simply declared over the setName() function in the application file.
|
inlineinherited |
Evaluate the boundary conditions for a neumann boundary segment.
This is the method for the case where the Neumann condition is potentially solution dependent
| 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 |
Negative values mean influx. E.g. for the mass balance that would be the mass flux in \( [ kg / (m^2 \cdot s)] \).
|
inline |
Evaluates the boundary conditions for a Neumann boundary segment.
|
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)] \).
|
inline |
Evaluates the source term for all phases within a given sub-control volume.
For this method, the priVars parameter stores the rate mass of a component is generated or annihilated per volume unit. Positive values mean that mass is created, negative ones mean that it vanishes.
The units must be according to either using mole or mass fractions (mole/(m^3*s) or kg/(m^3*s)).
|
inlineinherited |
Returns the spatial parameters object.
|
inlineinherited |
Returns the spatial parameters object.
|
inline |
Returns the temperature within the domain [K].
This problem assumes a temperature of 20 degrees Celsius.
|
inlineinherited |
Returns the temperature \(\mathrm{[K]}\) at a given global position.
This is not specific to the discretization. By default it just calls temperature().
| globalPos | The position in global coordinates where the temperature should be specified. |
|
inline |
Udpate the analytical temperature.
|
protectedinherited |
The gravity acceleration vector.
|
protectedinherited |
1.9.3