Regarding the issue of linear interpolation in spatial databases

PyLith
1

Dear everyone

I want to calculate the shear stress variation on another fault (Fault 1) through the coseismic slip of one fault (Fault 2 and Fault 3) on the intersecting fault
At the very beginning, I set up friction on fault1 and assigned simple SimpleDB spatial databases to fault2 and fault3. However, strangely enough, the nearest interpolation program could run, but the linear interpolation program got stuck at this step.（I think the data in the database is relatively evenly distributed.）

json (2).zip (17.2 KB)

7bca7f9c-7304-4ec3-8cc9-c1ccd8f095ef
7bca7f9c-7304-4ec3-8cc9-c1ccd8f095ef1145×713 400 KB

Share1
Share1845×451 90.6 KB

My Settings

[pylithapp.timedependent]
interfaces = [fault1,fault2,fault3]

[pylithapp.timedependent.interfaces]
fault1 = pylith.faults.FaultCohesiveDyn
fault2 = pylith.faults.FaultCohesiveKin
fault3 = pylith.faults.FaultCohesiveKin

[pylithapp.timedependent.interfaces.fault2]
label = fault2
edge = fault_slab_edge2
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.dimension = 2

[pylithapp.timedependent.interfaces.fault2.eq_srcs.rupture.slip_function]
slip.iohandler.filename =Slip_distribution2.spatialdb
slip.query_type = linear
slip.label = Final slip
slip_time = spatialdata.spatialdb.UniformDB
slip_time.label = Slip time
slip_time.values = [slip-time]
slip_time.data = [0.0*year]

[pylithapp.timedependent.interfaces.fault3]

label = fault3
edge = fault_slab_edge3
id=102
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.dimension = 2

[pylithapp.timedependent.interfaces.fault3.eq_srcs.rupture.slip_function]
slip.iohandler.filename =Slip_distribution3.spatialdb
slip.query_type = linear
slip.label = Final slip
slip_time = spatialdata.spatialdb.UniformDB
slip_time.label = Slip time
slip_time.values = [slip-time]
slip_time.data = [0.0*year]

[pylithapp.timedependent.interfaces.fault1]
label = fault1
edge = fault_slab_edge1
id=101
quadrature.cell = pylith.feassemble.FIATSimplex
quadrature.cell.dimension = 2

friction = pylith.friction.StaticFriction
friction.label = Static friction
friction.db_properties = spatialdata.spatialdb.UniformDB
friction.db_properties.label = Static friction
friction.db_properties.values = [friction-coefficient, cohesion]
friction. Db_propertiproperties.data = [0.4, 50.0*MPa]

Is there any way to solve this problem? In addition, I want to know whether setting burial edges at the intersection will affect the shear stress of fault1. The result I obtained by using nearest neighbor interpolation is several times larger than that of other studies, especially in the intersection part.

Share3
Share3569×482 18.7 KB

2

In your spatial databases for slip, you have points along a line, but you set data-dim = 2. For points on a line, the data-dim is 1 (0=point, 1=line, 2=area, 3=volume).

The buried edges of a fault should be along any of its edges that are not on external boundaries of the domain. Figure 19 in the current PyLith manual shows an example of buried edges. At the intersection of faults 1, 2, and 3, faults 2 and 3 should have a buried edge along the intersection with fault 1 (faults 2 and 3 terminate at fault 1), but fault 1 should only have a buried edge along its bottom and lateral edges.