Dear Etienne, unfortunately there is currently no example of such a simulation setup with two-phase flow. But there is a simple example for one-phase flow + transport example: https://git.iws.uni-stuttgart.de/dumux-repositories/dumux/-/tree/master/examples/1ptracer
Consider two things: (1) The "multi-domain” framework with coupling managers is only really needed for assembling coupled multi-domain Jacobians. In the sequential case you describe this should not be necessary. Basically the first solver dumps everything the second solver needs to know into some containers (std::vectors for example) and the second problem/spatialparams take these values. In the main fail you just solve two problems after each other and exchange data in-between. This is for example done in the 1ptracer example (only that there the tracer quantities don’t couple back but that’s the same process). (2) There is unfortunately no very consistent way of accessing the time discretisation information. This means if for example for your reactive transport problem you need both the saturation at the old time step and at the new time step, you need to manually overload the storage term evaluation in the base local residual class and make sure that the correct quantities for the respective time level are obtained from e.g. the spatial params (if you decide to store saturation from the other problem there). The spatial params would store both old and new saturations from the two phase flow model in this example. Hope this is good as a starting point Timo > On 8. Mar 2022, at 18:13, Etienne Ahusborde <etienne.ahusbo...@univ-pau.fr> > wrote: > > Hello DuMuX users, > > With my colleague we work on the simulation of thermo-hydro-chemical > processes in porous media. > > We developed a fully implicit scheme and now we would like to compare it with > a sequential scheme. > > For this we would like to solve a non isothermal compositionnal two phase > flow and then a reactive transport problem. > > For instance, the first subproblem would compute saturations, pressure, > temperature while the chemistry will be explicited. Then, these quantities > would be given to the second subproblem that would calculate mole fractions > of all chemical species and concentrations of minerals. The > dissolution/precipitation of minerals could modify the new porosity that > would be given to the first subproblem > > We looked how DUMUX could do this. It seems that the different coupling modes > currently available (boundary, embedded mixed dimension, conforming mixed > dimension facet) are not suitable but maybe we are wrong. > > We think that a strategy similar to the one used to coupled flow and > geomechanics (/geomechanics/poroelastic/couplingmanager.hh) could be a good > point of departure? > > Could you please confirm us what in your opinion could be the best way to > proceed using the current DuMuX capacities? > > Thanks in advance > > Regards > > Etienne > > > _______________________________________________ > DuMux mailing list > DuMux@listserv.uni-stuttgart.de > https://listserv.uni-stuttgart.de/mailman/listinfo/dumux _______________________________________________ DuMux mailing list DuMux@listserv.uni-stuttgart.de https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
