[meta] Reconstruct physical flux from DG solution
Note: This is a meta-task. It bundles several tasks together and is only closed once all these tasks are finished.
The DG solution is discontinuous by nature, but handles fluxes arising from solution jumps at interfaces are handled appropriately in the discretization scheme. This is not the case if fluxes are evaluated from the current solution without considering an appropriate reconstruction based on interface jumps. The reconstruction is also necessary to compute solute transport, since it requires a continuous physical flux as input.
For flux reconstruction we use an implementation from the dune-modelling module, which must be incorporated into the requirements and our testing image.
[ ] #66 (closed): Add dune-modelling as dependencyDelayed usage of
dune-modellingmodule until further development
- !52 (closed) Compute reconstructed fluxes in the solver routine
- #62 (closed): Replace current flux with reconstructed fluxes in the VTK output