Hi Edscott,
the momentum balance for the geomechanics is formulated in an incremental way. That meens it only considers changes with respect to the initial state. Nonetheless, the initial stress state matters for post-processing routines such a the evaluation of the total principal stresses and the failure criteria. I have another suspicion for your troubles: The standard test in the 2.12-release uses the el2p-amgbackend, which seems to cause problems for some cases. For that reason, I switched to a direct solver (SuperLU) for my examples. Maybe this is reason why your problem now works, too? One more thing: I managed to implement Mandel's problem this week for the fully coupled version of the code. It works fine and reproduces the analytical solution quite well. I was trying to implement it also for the decoupled version, but I ran into some problems. Anyway, I will commit the fully coupled version to the Beck2019a module on Monday (I am out of office today), so you can look into it. Cheers, Martin ________________________________ From: Dumux <[email protected]> on behalf of Ed Scott Wilson Garcia <[email protected]> Sent: Thursday, February 28, 2019 11:50:39 PM To: DuMuX User Forum Cc: Manuel Coronado Gallardo Subject: [DuMuX] el2p issue almost clear One further question, the force that presses down on the grid from the top, is that not a boundary condition? I see from the code that it is used to set the initial total stress field (isotropic, lithostatic), (i.e. stress[] = brineDensity_ * porosity * gravity * (depthBOR_ - globalPos[dim-1]) + (1 - porosity) * rockDensity * gravity * (depthBOR_ - globalPos[dim-1]);). If I do not use the force to set a Neumann boundary condition for uy, and instead use the values Beck2019 uses (all Neumann values set to zero), the initialization of the El2P_TestProblem problem runs just fine with data from the Sangnimnuan 2018 paper (injection episode turned off). If I am not wrong about this, then my whole trouble was with that misleading fact. Code from the fullycoupled examples has been most useful, and the decoupled is very interesting. And I thank Bernd and Beck very much. Best regards, Edscott _______________________________________________ Dumux mailing list [email protected] https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
_______________________________________________ Dumux mailing list [email protected] https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
