* allow setting a fixed value anywhere in the domain
* do not use boundaryType but a separate function to check for fixed cell
* set 1.0 on the diagonal in the system matrix and 0.0 for the off-diagonal entries
  to improve numerical accuracy

Bernd Flemisch authored 6 years ago and Kilian Weishaupt committed 6 years ago

Implement decision made in !1137:
- 1p gas: `gas`
- 1p liquid: `liq`
- 2p gas/liquid: `gas` and `liq`
- 2p liquid/liquid: `aq` and `napl`
- 3p: `gas`, `aq` and `napl`
- solid: `s`

Add `IsAqueous` to components to detect if a component is aqueous or not.
Adapt reference solutions.
Adjust `primaryVariableName` only where necessary to make the tests pass.
Will be treated in general in !1212.

* only local element needs to be bound

[test][navierstokesnc] Fix test test_msfreeflow. Make binaryDiffusionCoefficient have the same argument as in fluidsystem base class.

Kilian Weishaupt authored 6 years ago and Thomas Fetzer committed 6 years ago

* remove now obsolete NavierStokesNCIndices

* index of the pressure has to equal the one of the main component
* advection law does not need outflow information
* fix some reference solutions

Melanie Lipp authored 6 years ago and Dennis Gläser committed 6 years ago

With contributions of Farid Mohammadi <farid.mohammadi@iws.uni-stuttgart.de>.

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 6 years ago and Timo Koch committed 6 years ago

Contributions by Timo Koch, Kilian Weishaupt.