Dear Pylith developers,
I’m modeling dynamic rupture propagation in a 2D elastic block along a planar 1D fault. I overstress the block and observe the rupture propagating from the middle weaker section of the fault along the two stronger side section.
I’m using slip-weakening, rate and state and some custom friction models with spatially varying parameters: e.g. with slip weakening both fric coefficients = 0.6 for the middle section and mu_s = 0.65, mu_d = 0.6 for the side sections.
I think I’m facing an issue with the built in slip weakening model. It looks like the fault heals back to its original strength once my slip rate drops to zero. I have friction.force_healing = False
I originally thought it’s a mesh resolution issue, but the effect gets even worse with refining the mesh.
Please attached see the plot of shear stress, slip rate and slip along the fault.
Both simulations are the same, right side - twice finer mesh. Both are the same time in the simulation. It’s the first time step where the coarser fault seemed to heal, while the finer fault has been healing for multiple time steps by now (asymmetrical, non-elliptical slip, larger shear stress in places, where it should be weak).
Please let me know what you think.