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
Dumux@listserv.uni-stuttgart.de
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
