[CIG-SHORT] Single event dynamic earthquake rupture modeling in PyLith

[Brad Aagaard]

Mamum,

The traction perturbation to nucleate rupture is described in the SCEC 
benchmark description. These are posted at 
http://scecdata.usc.edu/cvws/benchmark_descriptions.html with the 
details for TPV102 at
http://scecdata.usc.edu/cvws/download/SCEC_validation_ageing_law.pdf.

The traction perturbations for nucleation in these benchmarks are mostly 
ad hoc in terms of their spatial and temporal evolution. Some of the 
first ones used square distributions with discontinuities. Obviously, 
this can cause ambiguity and numerical issues in implementing them. 
TPV102 uses a smoother distribution as do most of the later ones.

Regards,
Brad


On 05/25/2016 04:29 PM, Mamun Miah wrote:
> Dear Brad,
>
> I have been able to simulate dynamic rupture of a fault from fluid
> injection after following your links. I am using TOUGH2 for fluid
> simulation. However, I have a question in regard to the computation of
> spatial and temporal variation of tractions on the fault surface. In the
> Python script 'create_traction.py', you used a function F =
> mask*numpy.exp(r**2/(r**2-R**2)) and G =
> mask1*numpy.exp((t-T)**2/(t*(t-2*T))) + mask2*1.0. I am a bit unclear
> about how the combination of these two functions is replicating the
> stress transfer in the form of wave propagation to simulate dynamic
> rupture on the fault. I would highly appreciate if you kindly elaborate
> on this.
>
> Thanks,
> Mamun
>
> On Wed, Nov 4, 2015 at 11:42 AM, Brad Aagaard <baagaard at usgs.gov
> <mailto:baagaard at usgs.gov>> wrote:
>
>     The SCEC Dynamic Rupture Benchmark TPV102 is such an example. The
>     problem description is posted at
>     http://scecdata.usc.edu/cvws/cgi-bin/cvws.cgi
>
>     The PyLith input files (no mesh) are in the GitHug repo
>     https://github.com/geodynamics/pylith_benchmarks
>     in dynamic/scecdynrup/tpv102.
>
>     Note that this examples requires CUBIT and a computer cluster to
>     run. I am not sure the input files have been completely updated to
>     v2.1.0.
>
>     Regards,
>     Brad
>
>
>
>     On 11/4/15 11:36 AM, Mamun Miah wrote:
>
>         Hello Developers, I hope I am not too much trouble.  I really
>         appreciate
>         all the support that you guys are giving me.
>
>         I would like to model the dynamic rupture of a single event in
>         PyLith
>         governed by rate- and state-dependent friction law. Is there any
>         tutorial on a single dynamic rupture? I can't seem to find it in the
>         manual. Thanks!
>
>         Mamun
>         --
>         Thanks & Best Regards,
>         Mamun Miah
>         Graduate Research Assistant
>         Department of Civil Engineering
>         The University of Mississippi
>         Oxford, MS 38655
>         Phone: (662)-801-1904 <tel:%28662%29-801-1904>
>         Website: http://www.linkedin.com/pub/mamun-miah-eit/21/53/339
>
>
>         _______________________________________________
>         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
>
>
>     _______________________________________________
>     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
>
>
>
>
> --
> Thanks & Best Regards,
> Mamun Miah
> Graduate Research Assistant
> Department of Civil Engineering
> The University of Mississippi
> Oxford, MS 38655
> Phone: (662)-801-1904
> Website: http://www.linkedin.com/pub/mamun-miah-eit/21/53/339
>
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>