An error in the fluid compressibility computation in RichardsMechanics
As @joergbuchwald pointed out in his comment about MR. !3166 (merged), auto const beta_LR = 1 / rho ...
should be
auto const beta_LR = 1 / rho_LR
. It is correct. Therefore, the following two lines in RichardsMechanicsFEM-impl.h have to be corrected:
https://gitlab.opengeosys.org/ogs/ogs/-/blob/master/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h#L524-524 https://gitlab.opengeosys.org/ogs/ogs/-/blob/master/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h#L935-935