GitLab

Correct handling of gravitational flux contributions requires the
mpfa framework to be completely re-designed.

this commit furthermore introduces that the problem is given to the global containers
for volume variables and flux variables in the constructor.

TODO: implement ficks and fouriers law, as well as the interaction volume data handles
accordingly for problems involving diffusion/heat conduction

the lambdas are used to obtain the tensors related to a law, i.e. an
expression for a flux (advective, diffusive, conductive). We now introduce
a lambda factory that passes lambdas for the different physical processes.

The flux variables cache now inherits from caches defined in the different laws.
These also define a respective filler class that is called when updating the flux
variables caches. That enables us to define different caches for different laws.

Dennis Gläser authored Dec 14, 2016 and Timo Koch committed Dec 16, 2016

In this commit we introduce an mpfa implementation for fick's law
as well as a new implementation of the flux variables cache filler.
The filler implementation now inherits from process-related filler
classes which handle the setup of the advective and diffusive quantities
in the caches separately.

The outsideScvIdx in the faces is now different for two
scvfs on the same intersection, thus, boundary conditions will be
specified on a finer resolution. This commit furthermore adjusts the binding
of the elementvolumevariables and the fluxvariablescache to this setting.

Dennis Gläser authored Oct 22, 2016 and Timo Koch committed Nov 23, 2016

Up to now the types for the inner and boundary interaction volume seeds
were always identical. For the mpfa-l and other methods, it would be more efficient
if these types could differ, as e.g. the sizes of the containers are known at compile time
for the l method. This commit introduces the capability of the types to differ.

Dennis Gläser authored Oct 21, 2016 and Timo Koch committed Nov 23, 2016

when UseTpfaBoundaries is deactivated, the scaling of the neumann flux
has to reoccur using the deflected vol vars. otherwise the flux calculation
in the fluxVariables goes wrong. We have to keep in mind, that using the global
flux var cache for the mpfa also means that the neuman boundaries are assumed to be
constant over time!

Dennis Gläser authored Oct 19, 2016 and Timo Koch committed Nov 23, 2016

Updating the flux var cache is very expensive for mpfa schemes.
We thus introduce a new update() method in the elementfluxvarcaches.
Instead of binding again, the localjacobian calls this method. In the
tpfa case this forwards directly to bind when parameters are solution dependent
as we know we only store the darcy transmissibilities as the others are by default
solution dependent and are calculated on the fly. For mpfa, the fluxvarcachefiller
decides what to update (in case of compositional flow only the diffusive transmissibilities
will be updated when the permeability is not solution dependent).

Dennis Gläser authored Oct 19, 2016 and Timo Koch committed Nov 23, 2016

That means, that the transmissibilities are the same for each phase and we
do not need to store them for each phase. The interface, however, stays the same
in case one desires to implement a specialization of that (or mixed BC).

Dennis Gläser authored Oct 19, 2016 and Timo Koch committed Nov 23, 2016

Instead of specializing the fill function we specialize the whole class.

Dennis Gläser authored Oct 18, 2016 and Timo Koch committed Nov 23, 2016

the flux var cache filler will be used by the non-caching specialization as well,
which is why we put this in a new header.

Dennis Gläser authored Oct 12, 2016 and Timo Koch committed Nov 23, 2016

The treatment of mixed BC is questionable for cc methods. Not allowing
for it enables us to not have to store n interaction volume seeds at the boundary.

Dennis Gläser authored Aug 17, 2016 and Timo Koch committed Nov 23, 2016

The routine to fill the actual caches can be reused for the specialization
of the caches being disabled

Dennis Gläser authored Aug 15, 2016 and Timo Koch committed Nov 23, 2016

in mpfa schemes, k*gradh is recovered on neumann boundaries by dividing the neumann flux
to be applied by the ratio density/viscosity. On neumann boundaries we
thus have to make sure that this ratio is multiplied onto the calculated
fluxes in the flux variables.