Dear Etienne Ahusborde,

Sorry for the late reply.
I encountered the same problem in the beginning. There is a known bug in Dumux. 
The porosity in the 1p problem is evaluated only with the current solution, so 
the derivate of the storage term is wrong and that can explain why the 
Mandel-Cryer effect is not observed.

You may use the branch "temp/poromech-storagederivs-hacky-fix”, where the 
porosity error is corrected.

Best regards,
Yue Wang

On 17. Oct 2022, at 10:41, Etienne Ahusborde 
<[email protected]<mailto:[email protected]>> wrote:

Dear Dumux community,

With my colleagues, we are trying to perform an hydro-mechanical coupling model 
with Dumux.

We have first investigated the option to use the model "el2P" supplied in the 
test section of the multidomain coupling in version 3.4. In order to validate 
the coupling, we have been interesting in reproducing Mandel 's problem 
benchmark by taking into account of what has been carried out in Martin Beck 's 
publications and source code from 
https://git.iws.uni-stuttgart.de/dumux-pub/beck2019a .

Unfortunately, we are not capable to obtain good results on the pressure field 
where the Mandel-Cryer effect should be captured. As depicted on the images 
that we join to this mail, we are obtaining good results for horizontal and 
vertical displacements but not for the pressure which decreases directly at 
early times.

We have no idea of where comes from the issue, this is why we are kindly asking 
if someone has already been facing this situation.

Thank you for your help



<comp_pressure.pdf><Hdisplacement.pdf><Vdisplacement.pdf>_______________________________________________
DuMux mailing list
[email protected]<mailto:[email protected]>
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux

_______________________________________________
DuMux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux

Reply via email to