[CIG-SHORT] Dike Intrusion

[Charles Williams]

Hi Nick,

I don’t know if this is possible or not, but I believe the error you are getting is a missing letter in this line:

friction.db_initial_state = spatialdata.spatialdb.UniforDB

That should be UniformDB, I think.

Cheers,
Charles



> On 23/05/2017, at 3:54 PM, Nick Wogan <nicholaswogan at gmail.com> wrote:
> 
> Hi,
> 
> I was looking at some Pylith emails back in 2012 titled "Kinematic faults vs dynamic faults and fault opening.” In these emails was talk of including the capability to impose a stress boundary condition on a fault interface to simulate a dike intrusion.
> 
> I would like to make a simulation with an imposed pressure on a fault interface, which opens the fault in a dike like fashion. I tried this using the following settings.
> 
> [pylithapp.timedependent.interfaces]
> fault0 = pylith.faults.FaultCohesiveDyn
> 
> [pylithapp.timedependent.interfaces.fault0]
> label = fault0
> edge = fault0_edge
> id = 100
> zero_tolerance = 1.0e-5
> friction = pylith.friction.StaticFriction
> friction.label = Static friction
> quadrature.cell = pylith.feassemble.FIATSimplex
> quadrature.cell.dimension = 2
> friction.db_properties = spatialdata.spatialdb.UniformDB
> friction.db_properties.label = Static friction
> friction.db_properties.values = [friction-coefficient,cohesion]
> friction.db_properties.data = [0.0,0.0*Pa]
> friction.db_initial_state = spatialdata.spatialdb.UniforDB
> friction.db_initial_state.label = fault friction
> friction.db_initial_state.values = [traction-normal]
> friction.db_initial_state.data = [100000.0*Pa]
> 
> I get the following error:
> Assertion failed: (numDBStateVars > 0), function initialize, file ../../../pylith-2.2.0/libsrc/pylith/friction/FrictionModel.cc, line 187.
> 
> Is this simulation possible with polity? If so, what am I doing wrong?
> 
> Thanks,
> Nick
> 
> 
> 
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short

Charles Williams I Geodynamic Modeler
GNS Science I Te Pῡ Ao
1 Fairway Drive, Avalon 5010, PO Box 30368, Lower Hutt 5040, New Zealand
Ph 0064-4-570-4566 I Mob 0064-22-350-7326 I Fax 0064-4-570-4600
http://www.gns.cri.nz/ <http://www.gns.cri.nz/> I Email: C.Williams at gns.cri.nz <mailto:your.email at gns.cri.nz>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20170523/3208b35b/attachment.html>