[CIG-SHORT] fault meshing, SCEC mesh file

[Brad Aagaard]

Birendra,

If you want to use LaGriT, I strongly recommend that you contact Carl 
Gable (gable at lanl.gov) at LANL. He is one of the developers and can 
provide much better help than we can. I haven't used LaGriT in a long time.

Regards,
Brad


On 08/26/2015 12:46 PM, Birendra jha wrote:
> Thank you Brad and Charles. Right now, I am trying Lagrit instead of Trelis, given its flexibility in honoring interior surfaces (faults) of different sizes. I am using the mesh_tet4_1000m.lagrit file from pylith. I understand the basic commands after reading lagrit.lanl.gov/docs/commands.
>
> Can you please tell me how I can import a file in Lagrit to create a surface? Here is what I do:
>
> 1. Import fault.facet file in Trelis with make_elements on so as to get the triangles
> 2. Create a block from the surface
> 3. Export as Abaqus inp file, fault.inp
> 4. In the lagrit file, use
> cmo/create/mo_fault
> surface/surf_fault/intrface/sheet/mo_fault/< fault.inp
>
> I get the error Sheet CMO has the wrong dimensionality: mo_fault
>
> Thanks
> Birendra
>
> --------------------------------------------
> On Wed, 8/26/15, Charles Williams <willic3 at gmail.com> wrote:
>
>   Subject: Re: [CIG-SHORT] fault meshing, SCEC mesh file
>   To: cig-short at geodynamics.org
>   Date: Wednesday, August 26, 2015, 3:23 PM
>
>   There are two methods you
>   can use.  The first one involves creating an interior
>   volume.  You first cut the domain volume with the interior
>   volume (chop command), then cut the interior volume with the
>   fault surface.  The other method is to just cut the domain
>   volume with an extended fault surface, but only use a subset
>   of the fault nodes to define the fault nodeset.  I hope
>   that helps.
>
>   Cheers,
>   Charles
>
>
>   > On 27/08/2015, at 4:11 am, Brad Aagaard
>   <baagaard at usgs.gov>
>   wrote:
>   >
>   > On
>   08/26/2015 04:52 AM, Birendra jha wrote:
>   >> Hi
>   >>
>   >> I want to mesh a domain with a fault
>   that does not cut it. Similarly
>   >> to
>   Fig. 9 in Aagaard et al, JGR 2013. How do I honor the
>   fault
>   >> surface in Trelis without
>   using webcut? Are imprint and merge
>   >>
>   commands sufficient?
>   >>
>   >> I see that Fig. 10 in the above paper
>   does this. So I am looking for
>   >> the
>   mesh file for SCEC benchmark TPV13. I am interested in a 3D
>   tet
>   >> mesh. If there is a mesh file
>   for that, it will help me a lot.
>   >
>   > Figure 10 is just a diagram for the
>   simulation. It does not show the actual geometry used with
>   CUBIT. We are unaware of any methods for creating interior
>   surfaces in CUBIT that will be honored by the mesh. If
>   someone does know how to do this, please let us know!
>   >
>   > In the geometry.jou
>   file (https://github.com/geodynamics/pylith_benchmarks/blob/master/dynamic/scecdynrup/tpv210/geometry.jou)
>   we create an extended fault surface that spans across the
>   volume and use webcut to create the fault surface. We do
>   scribe lines on that surface so that we can create a nodeset
>   confined to the portion of the surface corresponding to the
>   fault and edges of cells line up with discontinuities in the
>   simulation parameters.
>   >
>   > WARNING: The journal files may or may not
>   work exactly as is for your version of CUBIT/Trelis. The ids
>   sometimes need to be updated as a result of how entities are
>   split by CUBIT. Using Idless journal files does remedy this
>   problem. We have tried to migrate the examples bundled with
>   PyLith to be independent of the CUBIT/Trelis version, but
>   have not migrated all of the benchmark journal files.
>   >
>   > Regards,
>   > Brad
>   >
>   >
>   _______________________________________________
>   > CIG-SHORT mailing list
>   > CIG-SHORT at geodynamics.org
>   > http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>
>   Charles A. Williams
>   Scientist
>   GNS Science
>   1 Fairway Drive, Avalon
>   PO Box
>   30368
>   Lower Hutt  5040
>   New
>   Zealand
>   ph (office): 0064-4570-4566
>   fax (office): 0064-4570-4600
>   C.Williams at gns.cri.nz
>
>   _______________________________________________
>   CIG-SHORT mailing list
>   CIG-SHORT at geodynamics.org
>   http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>