RegularizedVanGenuchten.hpp File Reference

Implementation of the regularized van Genuchten's capillary pressure / relative permeability <-> saturation relation. More...

#include "VanGenuchten.hpp"
#include "RegularizedVanGenuchtenParams.hpp"
#include <opm/material/common/Spline.hpp>
#include <algorithm>
Include dependency graph for RegularizedVanGenuchten.hpp:
This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes

class  Opm::RegularizedVanGenuchten< TraitsT, ParamsT >
 Implementation of the regularized van Genuchten's capillary pressure / relative permeability <-> saturation relation. More...
 

Namespaces

 Opm
 

Detailed Description

Implementation of the regularized van Genuchten's capillary pressure / relative permeability <-> saturation relation.

This class bundles the "raw" curves as static members and doesn't concern itself converting absolute to effective saturations and vice versa.

In order to avoid very steep gradients the marginal values are "regularized". This means that in stead of following the curve of the material law in these regions, some linear approximation is used. Doing this is not worse than following the material law. E.g. for very low wetting phase values the material laws predict infinite values for $p_c$ which is completely unphysical. In case of very high wetting phase saturations the difference between regularized and "pure" material law is not big.

Regularizing has the additional benefit of being numerically friendly: Newton's method does not like infinite gradients.

The implementation is accomplished as follows:

  • check whether we are in the range of regularization
    • yes: use the regularization
    • no: forward to the standard material law.

An example of the regularization of the capillary pressure curve is shown below:

regularizedVanGenuchten.png
See also
VanGenuchten