[CIG-SHORT] Slip-weakening friction model for cohesion?

[Brad Aagaard]

Ekaterina,

Are you using the PyLith binary distribution, the PyLith installer, or 
did you clone PyLith from GitHub?

For either of the first two, please give the version number. If cloning 
from GitHub, what command did you use to clone?

Thanks,
Brad


On 09/21/2017 03:32 PM, Ekaterina Bolotskaya wrote:
> Dear developers,
> 
> I'm currently trying to compile the friction template.
> I'm receiving the following error when running the configure file:
> 
> ./configure: line 16672: CIT_HEADER_MPI: command not found
> ./configure: line 16681: syntax error near unexpected token `3.6.1'
> ./configure: line 16681: `CIT_PATH_PETSC(3.6.1)'
> 
> I have everything mentioned in README installed.
> Could you tell me what the problem might be?
> 
> Thanks!
> 
> Best regards,
> Ekaterina Bolotskaya
> 
> PhD in Geophysics,
> Earth, Atmospheric and Planetary Science Department,
> Massachusetts Institute of Technology
> E-mail. bolee at mit.edu
> Mob. +1 (857) 284-2805
>           +7 (963) 995-36-33
> 
> ________________________________________
> From: CIG-SHORT [cig-short-bounces at geodynamics.org] on behalf of Brad Aagaard [baagaard at usgs.gov]
> Sent: Monday, August 28, 2017 12:00
> To: cig-short at geodynamics.org
> Subject: Re: [CIG-SHORT] Slip-weakening friction model for cohesion?
> 
> On 08/28/2017 04:14 AM, Ekaterina Bolotskaya wrote:
>> Dear Pylith developers,
>>
>> I'm currently running a 2D simulation which is a 2m (height) by 3m
>> (width) elastic block.
>> The block is initially compressed both vertically (2MPa) and
>> horizontally (2MPa) and then sheared with a rate of 0.045 MPa/year for
>> 100 years. Vertical displacement is fixed for all the boundaries: both
>> horizontal and vertical.
>>
>> The block has a horizontal fault in the middle which goes almost all the
>> way to the edges of the block (2.98m - fault length).
>> I want to have an ~1m long section with a smaller friction coefficient
>> (0.6) in the middle of the fault which would slip first and I want the
>> slip to then propagate to the two side regions that initially had higher
>> friction coefficient (0.8) but that should have same 0.6 coefficient
>> once the slip started happening. I'm implementing this using the
>> slip-weakening friction database.
>>
>> I'm now interested in implementing a similar thing with cohesion. I want
>> it to have zero value for the middle section of the fault and some
>> finite value for the two locked (side) sections and then change to zero
>> (for the side sections) in a linear or stepwise fashion as the slip
>> becomes larger than some threshold value (similar to what the
>> coefficient of friction does in slip-weakening model). Is there a way to
>> implement it with Pylith?
> 
> The physics you describe is different from the slip-weakening friction
> model implemented in PyLith. As Surendra mentioned, we discuss how to
> add friction models to PyLith in the "Extending PyLith" chapter of the
> manual, which discusses an example in templates/friction.
> 
> Regards,
> Brad
> _______________________________________________
> 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
>