[CIG-SHORT] Dike Intrusion

[Surendra Nadh Somala]

To add to what Brad said, db_initial_state for friction laws is used to
specify 'State Variable', which denotes the average age of contact
asperities in case of rate-and-state friction.  For other friction laws
like slip-weakening and time-weakening artificial state variables were
introduced.  The former has cumulative-slip and previous-slip while the
latter has elapsed-time.
Static friction, on the other hand, does not have any state variable.

Initial traction for faults is set using 'db_initial' instead of
'db_initial_state' which is used in your example.  In addition to
'db_initial', you can also set 'db_rate' or 'db_change' as detailed in
table 6.11 of Pylith manual version 2.2.0

-SS

___________
*Surendra Nadh Somala*
Assistant Professor
Department of Civil Engineering (Block-E, 208)
Indian Institute of Technology (IIT) Hyderabad
Kandi, Sangareddy
Telangana, India - 502285
Phone : +91 (0)40 2301-8457
http://www.iith.ac.in/~surendra/

On Tue, May 23, 2017 at 11:03 PM, Brad Aagaard <baagaard at usgs.gov> wrote:

> 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
>>
>>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20170523/01372a03/attachment.html>