Isothermal NAPL infiltration problem where LNAPL contaminates the unsaturated and the saturated groundwater zone. More...

#include <infiltrationproblem.hh>

Inheritance diagram for Opm::InfiltrationProblem< TypeTag >:

Public Member Functions
	InfiltrationProblem (Simulator &simulator)

void	finishInit ()
	Called by the Opm::Simulator in order to initialize the problem. More...

Problem parameters
bool	shouldWriteRestartFile () const
	Returns true if a restart file should be written to disk. More...

std::string	name () const
	The problem name. More...

void	endTimeStep ()
	Called by the simulator after each time integration. More...

template<class Context >
Scalar	temperature (const Context &, unsigned, unsigned) const

template<class Context >
const DimMatrix &	intrinsicPermeability (const Context &context, unsigned spaceIdx, unsigned timeIdx) const

template<class Context >
Scalar	porosity (const Context &, unsigned, unsigned) const

template<class Context >
const MaterialLawParams &	materialLawParams (const Context &, unsigned, unsigned) 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 &context, unsigned spaceIdx, unsigned timeIdx) const
	Evaluate the initial value 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...

Static Public Member Functions
static void	registerParameters ()

Detailed Description

template<class TypeTag>
class Opm::InfiltrationProblem< TypeTag >

Isothermal NAPL infiltration problem where LNAPL contaminates the unsaturated and the saturated groundwater zone.

The 2D domain of this test problem is 500 m long and 10 m deep, where the lower part represents a slightly inclined groundwater table, and the upper part is the vadose zone. A LNAPL (Non-Aqueous Phase Liquid which is lighter than water) infiltrates (modelled with a Neumann boundary condition) into the vadose zone. Upon reaching the water table, it spreads (since lighter than water) and migrates on top of the water table in the direction of the slope. On its way through the vadose zone, it leaves a trace of residually trapped immobile NAPL, which can in the following dissolve and evaporate slowly, and eventually be transported by advection and diffusion.

Left and right boundaries are constant hydraulic head boundaries (Dirichlet), Top and bottom are Neumann boundaries, all no-flow except for the small infiltration zone in the upper left part.

Constructor & Destructor Documentation

◆ InfiltrationProblem()

template<class TypeTag >

Opm::InfiltrationProblem< TypeTag >::InfiltrationProblem ( Simulator & simulator )

inline

Member Function Documentation

◆ boundary()

template<class TypeTag >

template<class Context >

void Opm::InfiltrationProblem< TypeTag >::boundary	(	BoundaryRateVector &	values,
		const Context &	context,
		unsigned	spaceIdx,
		unsigned	timeIdx
	)		const

inline

Evaluate the boundary conditions for a boundary segment.

Parameters

values	Stores the fluxes over the boundary segment.
context	The object representing the execution context from which this method is called.
spaceIdx	The local index of the spatial entity which represents the boundary segment.
timeIdx	The index used for the time discretization

◆ endTimeStep()

template<class TypeTag >

void Opm::InfiltrationProblem< 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)

◆ finishInit()

template<class TypeTag >

void Opm::InfiltrationProblem< 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()

◆ initial()

template<class TypeTag >

template<class Context >

void Opm::InfiltrationProblem< TypeTag >::initial	(	PrimaryVariables &	values,
		const Context &	context,
		unsigned	spaceIdx,
		unsigned	timeIdx
	)		const

inline

Evaluate the initial value for a control volume.

Parameters

values	Stores the primary variables.
context	The object representing the execution context from which this method is called.
spaceIdx	The local index of the spatial entity which represents the boundary segment.
timeIdx	The index used for the time discretization

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

◆ intrinsicPermeability()

template<class TypeTag >

template<class Context >

const DimMatrix & Opm::InfiltrationProblem< TypeTag >::intrinsicPermeability	(	const Context &	context,
		unsigned	spaceIdx,
		unsigned	timeIdx
	)		const

inline

◆ materialLawParams()

template<class TypeTag >

template<class Context >

const MaterialLawParams & Opm::InfiltrationProblem< TypeTag >::materialLawParams	(	const Context &	,
		unsigned	,
		unsigned
	)		const

inline

Referenced by Opm::InfiltrationProblem< TypeTag >::initial().

◆ name()

template<class TypeTag >

std::string Opm::InfiltrationProblem< 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::InfiltrationProblem< TypeTag >::porosity	(	const Context &	,
		unsigned	,
		unsigned
	)		const

inline

◆ registerParameters()

template<class TypeTag >

static void Opm::InfiltrationProblem< TypeTag >::registerParameters ( )

inlinestatic

◆ shouldWriteRestartFile()

template<class TypeTag >

bool Opm::InfiltrationProblem< TypeTag >::shouldWriteRestartFile ( ) const

inline

Returns true if a restart file should be written to disk.

The default behavior is to write one restart file every 10 time steps. This method should be overwritten by the implementation if the default behavior is deemed insufficient.

This problem writes a restart file after every time step.

◆ source()

template<class TypeTag >

template<class Context >

void Opm::InfiltrationProblem< TypeTag >::source	(	RateVector &	rate,
		const Context &	,
		unsigned	,
		unsigned
	)		const

inline

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

Parameters

rate	Stores the values of the volumetric creation/anihilition rates of the conserved quantities.
context	The object representing the execution context from which this method is called.
spaceIdx	The local index of the spatial entity which represents the boundary segment.
timeIdx	The 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::InfiltrationProblem< TypeTag >::temperature	(	const Context &	,
		unsigned	,
		unsigned
	)		const

inline

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

infiltrationproblem.hh

Public Member Functions

Static Public Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ InfiltrationProblem()

Member Function Documentation

◆ boundary()

◆ endTimeStep()

◆ finishInit()

◆ initial()

◆ intrinsicPermeability()

◆ materialLawParams()

◆ name()

◆ porosity()

◆ registerParameters()

◆ shouldWriteRestartFile()

◆ source()

◆ temperature()