- 23 Nov, 2016 40 commits
-
-
-
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.
-
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!
-
we reintroduce this property tag because for mpfa methods there can be saved a lot of computational time.
-
-
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).
-
The stencil method is not necessary in the interaction volumes. Furthermore, we store the global scvf indices and return a const reference instead of copying every time.
-
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).
-
Instead of specializing the fill function we specialize the whole class.
-
-
-
-
-
-
-
-
-
-
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.
-
-
-
-
-
-
Some methods will be required by the non-caching, as well as the caching specialization. These methods are now outsourced into a common helper class.
-
for multiphasic problems this led to a wrong size of the neumann flux vector.
-
-
-
-
-
-
-
-
-
Treatment of interior boundaries will be specialized in the multidimension module, where a coupling with a different on the element facets is realized.
-
-
-
-
-
-
