test problem for the sequential one-phase model. More...
#include <test/porousmediumflow/1p/sequential/test_1pproblem.hh>
test problem for the sequential one-phase model.
Public Member Functions | |
TestProblemOneP (TimeManager &timeManager, Grid &grid) | |
void | boundaryTypesAtPos (BoundaryTypes &bcTypes, const GlobalPosition &globalPos) const |
Specifies which kind of boundary condition should be used for which equation on a given boundary segment. More... | |
void | dirichlet (PrimaryVariables &values, const Intersection &intersection) const |
Evaluate the boundary conditions for a dirichlet control volume. More... | |
void | dirichletAtPos (PrimaryVariables &values, const GlobalPosition &globalPos) const |
Evaluate the boundary conditions for a dirichlet control volume. More... | |
void | neumannAtPos (PrimaryVariables &values, const GlobalPosition &globalPos) const |
Evaluate the boundary conditions for a neumann boundary segment. More... | |
void | sourceAtPos (PrimaryVariables &values, const GlobalPosition &globalPos) const |
Evaluate the source term for all phases within a given sub-control-volume. More... | |
void | initial (PrimaryVariables &values, const Element &element) const |
Evaluate the initial value for a control volume. More... | |
void | initialAtPos (PrimaryVariables &values, const GlobalPosition &globalPos) const |
Evaluate the initial value for a control volume. More... | |
void | init () |
Called by the TimeManager in order to initialize the problem. More... | |
void | preTimeStep () |
Called by TimeManager just before the time integration. More... | |
void | timeIntegration () |
Called by TimeManager in order to do a time integration on the model. More... | |
void | postTimeStep () |
Called by TimeManager whenever a solution for a timestep has been computed and the simulation time has been updated. More... | |
void | advanceTimeLevel () |
Called by the time manager after everything which can be done about the current time step is finished and the model should be prepared to do the next time integration. More... | |
Scalar | maxTimeStepSize () const |
Returns the user specified maximum time step size. More... | |
Scalar | timeStepSize () const |
Returns the current time step size [seconds]. More... | |
void | setTimeStepSize (Scalar dt) |
Sets the current time step size [seconds]. More... | |
Scalar | nextTimeStepSize (Scalar dt) |
Called by TimeManager whenever a solution for a timestep has been computed and the simulation time has been updated. More... | |
void | setOutputTimeInterval (const Scalar timeInterval) |
Sets a time interval for Output. More... | |
void | setOutputInterval (int interval) |
Sets the interval for Output. More... | |
bool | shouldWriteOutput () const |
Returns true if the current solution should be written to disk (i.e. as a VTK file) More... | |
void | writeOutput (bool verbose=true) |
Write the fields current solution into an VTK output file. More... | |
void | episodeEnd () |
Called when the end of an simulation episode is reached. More... | |
void | setName (const std::string &newName) |
Set the problem name. More... | |
const GridView & | gridView () const |
The GridView which used by the problem. More... | |
const VertexMapper & | vertexMapper () const |
Returns the mapper for vertices to indices. More... | |
const ElementMapper & | elementMapper () const |
Returns the mapper for elements to indices. More... | |
const GlobalPosition & | bBoxMin () const |
The coordinate of the corner of the GridView's bounding box with the smallest values. More... | |
const GlobalPosition & | bBoxMax () const |
The coordinate of the corner of the GridView's bounding box with the largest values. More... | |
TimeManager & | timeManager () |
Returns TimeManager object used by the simulation. More... | |
const TimeManager & | timeManager () const |
Returns TimeManager object used by the simulation. More... | |
Variables & | variables () |
Returns variables object. More... | |
const Variables & | variables () const |
Returns variables object. More... | |
Model & | model () |
Returns numerical model used for the problem. More... | |
const Model & | model () const |
Returns numerical model used for the problem. More... | |
Problem parameters | |
Scalar | temperature (const Element &element) const |
Returns the temperature within the domain. More... | |
Scalar | referencePressure (const Element &element) const |
Returns the reference pressure for evaluation of constitutive relations. More... | |
const GravityVector & | gravity () const |
Returns the acceleration due to gravity. More... | |
SpatialParams & | spatialParams () |
Returns the spatial parameters object. More... | |
const SpatialParams & | spatialParams () const |
Returns the spatial parameters object. More... | |
Problem parameters | |
std::string | name () const |
The problem name. More... | |
bool | shouldWriteRestartFile () const |
void | addOutputVtkFields () |
Scalar | temperatureAtPos (const GlobalPosition &globalPos) const |
Returns the temperature within the domain. More... | |
Scalar | referencePressureAtPos (const GlobalPosition &globalPos) const |
Returns the reference pressure for evaluation of constitutive relations. More... | |
void | source (PrimaryVariables &values, const Element &element) const |
source term [kg/(m^3 s)] More... | |
void | boundaryTypes (BoundaryTypes &bcType, const Intersection &intersection) const |
Returns the type of boundary condition. More... | |
void | dirichletAtPos (PrimaryVariables &values, const GlobalPosition &globalPos) const |
return dirichlet condition (pressure, [Pa]) More... | |
void | neumann (PrimaryVariables &values, const Intersection &intersection) const |
return neumann condition (flux, [kg/(m^2 s)]) More... | |
Restart mechanism | |
void | serialize () |
This method writes the complete state of the problem to the harddisk. More... | |
void | restart (double tRestart) |
This method restores the complete state of the problem from disk. More... | |
VtkMultiWriter & | resultWriter () |
VtkMultiWriter & | resultWriter () const |
|
inline |
|
inline |
|
inlineinherited |
Called by the time manager after everything which can be done about the current time step is finished and the model should be prepared to do the next time integration.
|
inlineinherited |
The coordinate of the corner of the GridView's bounding box with the largest values.
|
inlineinherited |
The coordinate of the corner of the GridView's bounding box with the smallest values.
|
inline |
Returns the type of boundary condition.
BC can be dirichlet (pressure) or neumann (flux).
|
inlineinherited |
Specifies which kind of boundary condition should be used for which equation on a given boundary segment.
bcTypes | The boundary types for the conservation equations |
globalPos | The position of the center of the boundary intersection |
|
inlineinherited |
Evaluate the boundary conditions for a dirichlet control volume.
values | The dirichlet values for the primary variables |
intersection | The boundary intersection |
For this method, the values parameter stores primary variables.
|
inlineinherited |
Evaluate the boundary conditions for a dirichlet control volume.
values | The dirichlet values for the primary variables |
globalPos | The position of the center of the boundary intersection |
For this method, the values parameter stores primary variables.
|
inline |
return dirichlet condition (pressure, [Pa])
|
inlineinherited |
Returns the mapper for elements to indices.
|
inlineinherited |
Called when the end of an simulation episode is reached.
|
inlineinherited |
Returns the acceleration due to gravity.
If the EnableGravity
property is true, this means \(\boldsymbol{g} = ( 0,\dots,\ -9.81)^T \), else \(\boldsymbol{g} = ( 0,\dots, 0)^T \)
|
inlineinherited |
The GridView which used by the problem.
|
inlineinherited |
Called by the TimeManager in order to initialize the problem.
|
inlineinherited |
Evaluate the initial value for a control volume.
values | The initial values for the primary variables |
element | The element |
For this method, the values parameter stores primary variables.
|
inlineinherited |
Evaluate the initial value for a control volume.
values | The dirichlet values for the primary variables |
globalPos | The position of the center of the finite volume for which the initial values ought to be set (in global coordinates) |
For this method, the values parameter stores primary variables.
|
inlineinherited |
Returns the user specified maximum time step size.
Overload in problem for custom needs.
|
inlineinherited |
Returns numerical model used for the problem.
|
inlineinherited |
Returns numerical model used for the problem.
|
inline |
The problem name.
This is used as a prefix for files generated by the simulation.
|
inline |
return neumann condition (flux, [kg/(m^2 s)])
|
inlineinherited |
Evaluate the boundary conditions for a neumann boundary segment.
values | The neumann values for the conservation equations [kg / (m^2 *s )] |
globalPos | The position of the center of the boundary intersection |
For this method, the values parameter stores the mass flux in normal direction of each phase. Negative values mean influx.
|
inlineinherited |
Called by TimeManager whenever a solution for a timestep has been computed and the simulation time has been updated.
|
inlineinherited |
Called by TimeManager whenever a solution for a timestep has been computed and the simulation time has been updated.
This is used to do some janitorial tasks like writing the current solution to disk.
|
inlineinherited |
Called by TimeManager just before the time integration.
|
inlineinherited |
Returns the reference pressure for evaluation of constitutive relations.
element | The element |
|
inline |
Returns the reference pressure for evaluation of constitutive relations.
|
inlineinherited |
This method restores the complete state of the problem from disk.
It is the inverse of the serialize() method.
|
inlineprotectedinherited |
|
inlineprotectedinherited |
|
inlineinherited |
|
inlineinherited |
Set the problem name.
This function sets the simulation name, which should be called before the application problem is declared! If not, the default name "sim" will be used.
|
inlineinherited |
Sets the interval for Output.
The default is 1 -> Output every time step
|
inlineinherited |
Sets a time interval for Output.
The default is 0.0 -> Output determined by output number interval (setOutputInterval(int)
)
|
inlineinherited |
Sets the current time step size [seconds].
|
inlineinherited |
Returns true if the current solution should be written to disk (i.e. as a VTK file)
The default behaviour is to write out every the solution for very time step. This file is intented to be overwritten by the implementation.
|
inline |
|
inline |
source term [kg/(m^3 s)]
|
inlineinherited |
Evaluate the source term for all phases within a given sub-control-volume.
values | The source and sink values for the conservation equations |
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 rate mass generated or annihilate per volume unit. Positive values mean that mass is created, negative ones mean that it vanishes.
|
inlineinherited |
Returns the spatial parameters object.
|
inlineinherited |
Returns the spatial parameters object.
|
inlineinherited |
Returns the temperature within the domain.
element | The element |
|
inline |
Returns the temperature within the domain.
This problem assumes a temperature of 10 degrees Celsius.
|
inlineinherited |
Called by TimeManager in order to do a time integration on the model.
|
inlineinherited |
Returns TimeManager object used by the simulation.
|
inlineinherited |
Returns TimeManager object used by the simulation.
|
inlineinherited |
Returns the current time step size [seconds].
|
inlineinherited |
Returns variables object.
|
inlineinherited |
Returns variables object.
|
inlineinherited |
Returns the mapper for vertices to indices.
|
inlineinherited |
Write the fields current solution into an VTK output file.