[CIG-SHORT] Dynamic rupture on realistic mesh

[Huihui Weng]

Hi Brad,

I have rerun this case by single computer core. In this situation, the progress will abort with error: RuntimeError: WARNING! Fault opening with nonzero traction., v_fault: 7990, opening: 1.0026e-07, normal traction: -0.00370101. Please see the attached file for detailed log. 

The most important problem I concern about is why the normal and strike slips always accompany with the dip slip in the curved fault, even only assigning the shear stress at the dip direction. How can I inhibit these unnecessary slip? Could I do it by adjust the parameters of PETSc? Thanks.

Best,
Huihui


> On 9 Nov 2017, at 1:58 AM, Brad Aagaard <baagaard at usgs.gov> wrote:
> 
> Huiui,
> 
> My guess is that the "freeze" is occurring when PyLith has encountered an error and at least one process has aborted, but the MPI implementation is not detecting it. When this has happened to me, I will see there errors sent to stdout or stderr, which may be on a local compute node if you are on a cluster (this is MPI/machine dependent).
> 
> I will look deeper into your description when I get a chance.
> 
> Regards,
> Brad
> 
> 
> On 11/8/17 4:48 AM, Huihui Weng wrote:
>> Hi Brad,
>> 
>> When I try to simulate dynamic rupture with realistic mesh, I have a
>> similar problem with the one Hongfeng asked several years ago. But in my
>> case, the normal stress is large enough to hold the fault (attached
>> pictures).
>> 
>> For test, I let the dip shear stress at one node reaching the fault
>> strength, but I find it also has slips at both normal and strike
>> directions (as the attached pictures show, information at left bottom),
>> i.e., dip slip ~ 0.05m, strike slip ~ 0.005m, normal slip ~ 8e-5m (fault
>> open with normal stress large enough!). The simulating progress of this
>> case will freeze at some time step without any errors, and the freezing
>> time step depends on the parameter of zero_tolerance (e.g., 1.0e-7).
>> Once the progress freezes, I need to kill it by hand by Linux commands.
>> 
>> Since only the dip shear stress and normal stress are not zeros in this
>> case, how could it lead to strike and normal slip? It does seem it is
>> the abnormal normal slip that “freezes" the simulating progress. How can
>> eliminate the slip at the normal direction for such realistic fault? Thanks.
>> 
>> Best,
>> Huihui
>> 
> 
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org <mailto:CIG-SHORT at geodynamics.org>
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20171109/a5c27fe1/attachment-0002.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: log
Type: application/octet-stream
Size: 228334 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20171109/a5c27fe1/attachment-0001.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20171109/a5c27fe1/attachment-0003.html>