Dear Dumux community, First, we wish you an happy new year.
Since our last questions about flow and poromechanical coupling on Mandel problem, thanks to the help of Yue Wang, we have been able to perform simulations by using the coupling manager for fully implicit approach and we greatly thank you for that.
Then we tried to reproduce in version 3.5 the test case « 4.2 CO2 injection heterogeneous » proposed by Beck in Accuracy of fully coupled and sequential approaches for modeling hydro- and geomechanical processes, Computational Geosciences (2020) 24:1707–1723.
We used the implicit model el2p from multidomain example. We looked in details the source of the code in version 2.12 but when we tried to consider the same boundary conditions in version 3.5, we obtain some convergence issues. The only solution to reach solver convergence is to impose Dirichlet conditions on all boundaries for all displacements but results are obviously different.
Has somebody already tried to perform this case in a recent version of DUMUX?
Moreover, we would like to compare the results using implicit and sequential approach by using a fixed stress scheme where the two sub-problems (flow and poromechanics) are decoupled and solved sequentially. Unfortunately we are facing some difficulties on how to set up communications between the two sub-problems.
Like the flow problem is solved with a cell-centered scheme, we need to extrapolate the pressure in order to express it from cell to vertex. This is necessary to pass it as argument to the poromechanics problem which is formulated with the box method (vertex).
The same problem occur in the other direction, from poromechanical problem to flow problem where an interpolation is necessary to pass from vertex data to cell data.
Are there some existing tools in DUMUX to make some conversions in a simple way ?
Thank you for your answer Best regards
_______________________________________________ DuMux mailing list [email protected] https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
