GEOS-DEV / GEOS

GEOS Simulation Framework
GNU Lesser General Public License v2.1
203 stars 80 forks source link

Stress artifacts in parallel mechanics solves #3119

Open ryar9534 opened 1 month ago

ryar9534 commented 1 month ago

Describe the bug I am seeing some odd stress behavior in my poromechanical initialization of HI24L. In particular, there were some oscillations in stress on the boundaries between MPI ranks. On further testing I have narrowed it down to the mechanics solver, which seems to exhibit the same behavior. Here is an xml for a mechanics only sim on the underburden of HI24L: HI24L_mech_debug.zip

I ran with 1 rank (left) and 4 ranks (right), and slightly different results are obtained:

Screenshot 2024-05-10 093737

And the differences seem to line up with the MPI rank boundaries:

image

I have yet to see this effect on any other mesh, so perhaps its HI specific.