Improve Richards model - Add RichardsNC model

See merge request !360

[test][richardsnc] Add mass-conservation-output. Use heavy water tracer. Flow rate in kg/s instead of mol/s

The end point capillary pressure is the capillary pressure at sw=1.
For Brooks-Corey laws it is also called entry pressure.

The components are only in the water phase. The number of components
is statically extracted from the fluidsystem. Simply changing the
fluid system and boundary and initial conditions is sufficient to setup
a new model with n tracers.

* a main component and
* a conservative tracer component

Fix/2p2c elem sol

See merge request !384

mixeddimension-facetcoupling

See merge request !380

Dennis Gläser authored 8 years ago and Timo Koch committed 8 years ago

by passing the implementation of the actual method's interaction volume type
to the o-method interaction volume we are able to overload functionalities in
derived classes and still use different interaction volume types on the boundary and
inside the domain. This was not possible through the property system as it was unclear
which interaction volume is the actual implementation.

Dennis Gläser authored 8 years ago and Timo Koch committed 8 years ago

We now call the Implementation from the o-method interaction volume because we need
to incorporate some terms in the problems with facet coupling.
Todo: Problems using interior interaction volumes that inherit from the o-method
specialization AND use the o-method on the boundaries do not work at the moment. We
would have to tell the o-method the actual implementation somehow introducing a new
template parameter for the implementation.