Because of the two simulations (or models) that we want to couple, Richards and Transport, and the fact that they can have different time steps, it is necessary to store their solutions and be able to provide intermediate results by some evaluation policy (i.e. linear interpolation, previous solution, or last solution).
dune-modelling provide a solution, nevertheless it has two drawbacks:
- Simple containers always do a hard copy of the solutions. (easy to change)
SolutionStoragerequire a complete interface made with dune modelling with the definition of
In order to continue using the same interface, I would like to require that a container of several solutions to mirror or extend a
GridFunction such that is easy to be reused in other contexts. For example, that way will be easy to store the solution in arbitrary times with the
VTKWriter (see discussion on #102). On the other hand, It will select the right solution when a solvers asks the local operators to set the time.
How to test the implementation?
- An indirect form to test this is to check that mass is being conserved on the solute solution