GitLab

Katharina Heck authored Jun 28, 2019 and Ned Coltman committed Mar 26, 2020

johansen

* Use l in krw/krn instead of hard-coded l=0.5
* Correct the non-wetting relperm to comply with the Mualem model
* Correct the non-wetting relperm derivative

The parameter is used in the relative permeability curve in the Mualem model.
Mualem (1976) suggested that 0.5 is a good value for many soils. Sometimes a
value different from 0.5 is provided, which is why this patch implements it as additional variable.

This will lead to a more unified interface for effective laws.
Any dependence on spatial parameters ca be introduced via
the volume variables and seems misplaced here. Also, effective
diffusion laws do not have such a blown-up interface anyway.

Also replace all calls to isLiquid by equivalent calls to isGas

Bernd Flemisch authored Jun 19, 2018 and Timo Koch committed Jun 20, 2018

Add a corresponding python script `bin/util/create_cmakelists.py`.
Move the folder `bin/moduleutil` to `bin/util`.
Run the script.

Remove buggy and obsolete PhilToPhobLaw -> use wettingPhase() in spatial params to set the wetting phase index

This speeds up the simulations significantly as reconstruction of
sn only has to be done at material interfaces. Having an equality
comparison operator in material parameters allow their identification.
This commit introduces a comparison operator for the van-genuchten parameters.

Fluid systems exported phase indices for wetting and nonwetting components. The fluid system
is not the right place to decide about wettability. This introduces wettingPhase() in the spatial
params which returns the index of the phase in the fluidsystem we consider the wetting phase.
This way wettability can change temporally and spatially.
Indices throughout the 2p/2p2c/2pnc nc and mpnc where renamed accordingly.

Dennis Gläser authored Apr 05, 2018 and Timo Koch committed Apr 15, 2018

Contributions by Timo Koch, Kilian Weishaupt.