Parameters< Number, LType > Class Template Reference

Detailed Description

template<typename Number, typename LType>
class Parameters< Number, LType >

Program solving water flow and contaminant transport in porous medium. One phase unsaturated water flow is driven by Richards equation, variable is hydraulic head and uses van Genuchten-Mualem model for water retention curve and hydraulic permeability. Contaminant part has two components C0, C1 representing concentrations. They are transported via advection and diffusion (no dispersion).

Whole system is split into two parts: water flow and contaminant transport. Water flow does not depend on contaminants and can be solved separately. The contaminant system is split by operator splitting technique into two parts: advection and diffusion. Second order (in time) symmetric Strang splitting is implemented. Firstly we take half step of convection, then full diffusion step follower by second convection half step. Operator splitting is iterative and we stop iterating when the correction size is small.

All three parts use finite volume elements on a regular 2D axiparallel grid. Time stepping scheme is implicit Euler.

Examples

recipe-operator-splitting.cc.

---

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

• doc/Recipes/recipe-operator-splitting.cc