A fault crossing two materials in 3-D

PyLith
, ,
1

Dear PyLith developing team,

Thanks for your continuous effort.

I am building a toy model of a fault embedded in a box cutting two layers of materials. The setup looks like this:

IMG_0002
IMG_0002652×605 38.6 KB

I took care to mark the fault edges and exclude the curve intersecting two layers (the red dashed line in the figure). But PyLith doesn’t seem to like it and failed to initialize the fault. Attached are the meshing script and cfg files.

debug.zip (6.6 MB)
Would it be possible to mesh this kind of geometry in 3-D? What should I do with my GMSH script?

Thanks in advance!
Yifan

2

What fault edges are you marking (can you label them in your diagram). I do not understand what you mean by “exclude the cuve intersecting the two layers”? The fault doesn’t care about the boundary between the two layers. In Gmsh, I would embed the fault in each of the volumes as demonstrated in the crustal-strikeslip-3d example .