Hi everyone,
I am trying to simulate rainfallinduced landslide using the Richards mechanics process and successful doing so in the 2D model and now I am trying to replicate all the procedures that I have performed on the 2D model to the 3D using highperformance computing (PETSc). But I have encountered a problem with constitutive material properties’ definition. AFIK, there are two ways to define material properties in OpenGeoSys.

The first one is by defining material ID in the constitutive equation definition stage at the initial part of the project file. But this does not work with MFront if two or more material layers present, even on the 2D model and nonparallel model. Actually, this is a known issue raised in GitLab several months ago but it is not solved yet. (Issue: https://gitlab.opengeosys.org/ogs/ogs//issues/3063)

The second method is by defining constitutive parameter by using group_id_properties (which I have used in the 2D model). This works pretty well in 2D and non_parallel models until you are using the same constitutive equation for all layers but have a problem with the parallel model. Actually, it does not give any error while running just the simulation keeps showing running but nothing happens. This problem happens not only while using MFront but also with the constitutive equation native to the OpenGeoSys.
Does anyone has encountered this problem before or know the way around it? I know parallel FEM scheme with PETSc cannot be used by all processes using linear and quadratic function spaces simultaneously. I am using a linear approximation for both displacement and pressure in this case.
Way to reproduce the problem
 Compile OGS enabling MFront and PETSc.
 Run OGS with
Trial_1.prj
to reproduce problems related to point 1 andTrial_2.prj
to reproduce problems related to point 2.
Input files:
3D_Trial_Model.zip (147.8 KB)
Thanks