[CIG-SHORT] question for dynamic simulation via PyLith

[Matthew Knepley]

On Thu, Jan 26, 2017 at 10:54 AM, Brad Aagaard <baagaard at usgs.gov> wrote:

> On 1/24/17 7:13 AM, Yasemin Korkusuz wrote:
>
>> Hello,
>>
>> I am writing you from Turkey. I would like to use PyLith code in my PhD
>> thesis study in order to obtain "a _dynamic earthquake fracture
>> simulation_ in the Sea of Marmara concerning a _non-planar fault
>> geometry_ and _different (heterogenous) stress orientations_ for each
>> segment".
>>
>> /     1. Could you please make some suggestions about how should I start?/
>>
>
> Read through the manual and find examples that are most similar to what
> you want to do. Change them piece by piece so that they get closer to your
> problem. You may also want to look at the dyanmic spontaneous rupture
> benchmarks (although they are not all updated for the latest version of
> PyLith) which are in the pylith_benchmarks repository (
> https://github.com/geodynamics/pylith_benchmarks).
>
> I would like to use medium information, and initial stress amounts as
>> inputs and obtain stress amounts with slip and slip rates after the
>> propagation of the rupture (as outputs).
>> /
>> /
>> /      2. Thus, which PyLith example is the best for my purpose?/
>>
>
> See above.
>
> Also, as I am not living in US, I can not use CUBIT for the mesh
>> generation and LaGrit has many library requirements which did not work
>> on my computer. Actually, I used GMSH for another code last year.
>>
>
> You may be able to use Trelis from csimsoft. This is for academic and
> other use.
>
>       3.  /May I use GMSH for triangular mesh for PyLith?/
>>
>
It does not work now, so for your initial runs, please do what Brad
suggests. However, we could possibly get this
going. The problem would be making sure boundaries are properly marked.

  Thanks,

     Matt


> /     4.   If so, what do you suggest about conversion of the mesh
>> output file ?/
>>
>
> You would need to convert the mesh to the ASCII format PyLith uses in the
> trivial examples or the preferred method of the ExodusII form of NetCDF (
> https://cubit.sandia.gov/public/13.2/help_manual/WebHelp/
> finite_element_model/exodus/exodus2_file_specification.htm).
>
> /     5.  Also, I am not sure whether all types of meshes (triangular,
>> square,cubic) can be used in all examples of PyLith?/
>>
>
> This is discussed in the manual.
>
>
> Regards,
> Brad
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> 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/20170126/de06f230/attachment-0001.html>