Opm::FingerProblem< TypeTag > Class Template Reference

Two-phase problem featuring some gravity-driven saturation fingers. More...

#include <fingerproblem.hh>

Inheritance diagram for Opm::FingerProblem< TypeTag >:
Inheritance graph

Public Types

using RestrictProlongOperator = CopyRestrictProlong< Grid, MaterialLawParamsContainer >
 

Public Member Functions

 FingerProblem (Simulator &simulator)
 
Soil parameters
template<class Context >
Scalar temperature (const Context &, unsigned, unsigned) const
 
template<class Context >
const DimMatrix & intrinsicPermeability (const Context &, unsigned, unsigned) const
 
template<class Context >
Scalar porosity (const Context &, unsigned, unsigned) const
 
template<class Context >
MaterialLawParams & materialLawParams (const Context &context, unsigned spaceIdx, unsigned timeIdx)
 
template<class Context >
const MaterialLawParams & materialLawParams (const Context &context, unsigned spaceIdx, unsigned timeIdx) const
 
Boundary conditions
template<class Context >
void boundary (BoundaryRateVector &values, const Context &context, unsigned spaceIdx, unsigned timeIdx) const
 Evaluate the boundary conditions for a boundary segment. More...
 
Volumetric terms
template<class Context >
void initial (PrimaryVariables &values, const Context &, unsigned, unsigned) const
 Evaluate the initial value for a control volume. More...
 
template<class Context >
void constraints (Constraints &constraints, const Context &context, unsigned spaceIdx, unsigned timeIdx) const
 Evaluate the constraints for a control volume. More...
 
template<class Context >
void source (RateVector &rate, const Context &, unsigned, unsigned) const
 Evaluate the source term for all phases within a given sub-control-volume. More...
 

Auxiliary methods

RestrictProlongOperator restrictProlongOperator ()
 return restriction and prolongation operator More...
 
std::string name () const
 The problem name. More...
 
void finishInit ()
 Called by the Opm::Simulator in order to initialize the problem. More...
 
void endTimeStep ()
 Called by the simulator after each time integration. More...
 
static void registerParameters ()
 

Detailed Description

template<class TypeTag>
class Opm::FingerProblem< TypeTag >

Two-phase problem featuring some gravity-driven saturation fingers.

The domain of this problem is sized 10cm times 1m and is initially dry. Water is then injected at three locations on the top of the domain which leads to gravity fingering. The boundary conditions used are no-flow for the left and right and top of the domain and free-flow at the bottom. This problem uses the Parker-Lenhard hystersis model which might lead to non-monotonic saturation in the fingers if the right material parameters is chosen and the spatial discretization is fine enough.

Member Typedef Documentation

◆ RestrictProlongOperator

template<class TypeTag >
using Opm::FingerProblem< TypeTag >::RestrictProlongOperator = CopyRestrictProlong< Grid, MaterialLawParamsContainer >

Constructor & Destructor Documentation

◆ FingerProblem()

template<class TypeTag >
Opm::FingerProblem< TypeTag >::FingerProblem ( Simulator simulator)
inline

Member Function Documentation

◆ boundary()

template<class TypeTag >
template<class Context >
void Opm::FingerProblem< TypeTag >::boundary ( BoundaryRateVector &  values,
const Context &  context,
unsigned  spaceIdx,
unsigned  timeIdx 
) const
inline

Evaluate the boundary conditions for a boundary segment.

Parameters
valuesStores the fluxes over the boundary segment.
contextThe object representing the execution context from which this method is called.
spaceIdxThe local index of the spatial entity which represents the boundary segment.
timeIdxThe index used for the time discretization

◆ constraints()

template<class TypeTag >
template<class Context >
void Opm::FingerProblem< TypeTag >::constraints ( Constraints &  constraints,
const Context &  context,
unsigned  spaceIdx,
unsigned  timeIdx 
) const
inline

Evaluate the constraints for a control volume.

Parameters
constraintsStores the values of the primary variables at a given spatial and temporal location.
contextThe object representing the execution context from which this method is called.
spaceIdxThe local index of the spatial entity which represents the boundary segment.
timeIdxThe index used for the time discretization

References Opm::FingerProblem< TypeTag >::constraints().

Referenced by Opm::FingerProblem< TypeTag >::constraints().

◆ endTimeStep()

template<class TypeTag >
void Opm::FingerProblem< TypeTag >::endTimeStep ( )
inline

Called by the simulator after each time integration.

This method is intended to do some post processing of the solution. (e.g., some additional output)

References Opm::FingerProblem< TypeTag >::materialLawParams().

◆ finishInit()

template<class TypeTag >
void Opm::FingerProblem< TypeTag >::finishInit ( )
inline

Called by the Opm::Simulator in order to initialize the problem.

If you overload this method don't forget to call ParentType::finishInit()

References Opm::Parameters::reset().

◆ initial()

template<class TypeTag >
template<class Context >
void Opm::FingerProblem< TypeTag >::initial ( PrimaryVariables &  values,
const Context &  ,
unsigned  ,
unsigned   
) const
inline

Evaluate the initial value for a control volume.

Parameters
valuesStores the primary variables.
contextThe object representing the execution context from which this method is called.
spaceIdxThe local index of the spatial entity which represents the boundary segment.
timeIdxThe index used for the time discretization

◆ intrinsicPermeability()

template<class TypeTag >
template<class Context >
const DimMatrix & Opm::FingerProblem< TypeTag >::intrinsicPermeability ( const Context &  ,
unsigned  ,
unsigned   
) const
inline

◆ materialLawParams() [1/2]

template<class TypeTag >
template<class Context >
MaterialLawParams & Opm::FingerProblem< TypeTag >::materialLawParams ( const Context &  context,
unsigned  spaceIdx,
unsigned  timeIdx 
)
inline

◆ materialLawParams() [2/2]

template<class TypeTag >
template<class Context >
const MaterialLawParams & Opm::FingerProblem< TypeTag >::materialLawParams ( const Context &  context,
unsigned  spaceIdx,
unsigned  timeIdx 
) const
inline

◆ name()

template<class TypeTag >
std::string Opm::FingerProblem< TypeTag >::name ( ) const
inline

The problem name.

This is used as a prefix for files generated by the simulation. It is highly recommend to overwrite this method in the concrete problem which is simulated.

◆ porosity()

template<class TypeTag >
template<class Context >
Scalar Opm::FingerProblem< TypeTag >::porosity ( const Context &  ,
unsigned  ,
unsigned   
) const
inline

◆ registerParameters()

template<class TypeTag >
static void Opm::FingerProblem< TypeTag >::registerParameters ( )
inlinestatic

References Opm::Properties::dim.

◆ restrictProlongOperator()

template<class TypeTag >
RestrictProlongOperator Opm::FingerProblem< TypeTag >::restrictProlongOperator ( )
inline

return restriction and prolongation operator

Note
This method has to be overloaded by the implementation.

◆ source()

template<class TypeTag >
template<class Context >
void Opm::FingerProblem< TypeTag >::source ( RateVector &  rate,
const Context &  ,
unsigned  ,
unsigned   
) const
inline

Evaluate the source term for all phases within a given sub-control-volume.

Parameters
rateStores the values of the volumetric creation/anihilition rates of the conserved quantities.
contextThe object representing the execution context from which this method is called.
spaceIdxThe local index of the spatial entity which represents the boundary segment.
timeIdxThe index used for the time discretization

For this problem, the source term of all components is 0 everywhere.

◆ temperature()

template<class TypeTag >
template<class Context >
Scalar Opm::FingerProblem< TypeTag >::temperature ( const Context &  ,
unsigned  ,
unsigned   
) const
inline


The documentation for this class was generated from the following file: