[CIG-SHORT] Dike Intrusion

[Brad Aagaard]

Nick,

Your db_initial_state database is trying to supply initial state 
variables for a friction model that does not have state variables.

See examples/3d/hex8/step20 for how to specify a normal traction 
perturbation for modeling a dike intrusion.

Brad


On 05/23/2017 07:23 AM, Nick Wogan wrote:
> Hi Charles,
> 
> Thanks for the response. I fixed the typo and I still get the same 
> error. Is there another possible source for this sort of error?
> 
> Assertion failed: (numDBStateVars > 0), function initialize, file 
> ../../../pylith-2.2.0/libsrc/pylith/friction/FrictionModel.cc, line 187.
> 
> Thanks,
> Nick
> 
> 
> 
> On Tue, May 23, 2017 at 1:13 AM, Charles Williams <willic3 at gmail.com 
> <mailto:willic3 at gmail.com>> wrote:
> 
>     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
>>     <mailto: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
>>     <http://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 <mailto:CIG-SHORT at geodynamics.org>
>>     http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>>     <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/**I**Email: **C.Williams at gns.cri.nz*
>     <mailto:your.email at gns.cri.nz>
> 
> 
>     _______________________________________________
>     CIG-SHORT mailing list
>     CIG-SHORT at geodynamics.org <mailto:CIG-SHORT at geodynamics.org>
>     http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>     <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
>