[CIG-SHORT] Dynamic faulting under gravity

[luigi vadacca]

Hi Birendra,

you can try to:

1 - refine the mesh on the fault;
2 - increase the number of the friction_ksp_max_it

Moreover, if you switch the choesion value from 0 to ...Pa then you
can help the fault to stay locked at time step=0. It could help the
nonlinear solver to converge.


Luigi




2012/11/19 Birendra jha <bjha7333 at yahoo.com>:
> Dear all
>
> I have a relatively simple case of a slanted fault plane (along y axis) inside a layered box. I am using FaultCohesiveDyn with static friction model (please see the attached file). +x, -x boundaries are fixed in x,y directions and bottom is fixed in z direction. y faces are free.
> I have gravity on and I use initial_stress.spatialdb with a density of 2400 kg/m3.
>
> The nonlinear solver is not converging. Can you please help? What am I missing? There are four thin layers in the grid sandwiched between thick layers above and below. Is that the problem?
>
> I am using the Pylith 1.8.0 windows binary.
>
> Thanks and regards
> Birendra
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>



-- 
Luigi Vadacca
PhD student
Istituto Nazionale di Geofisica e Vulcanologia
via di Vigna Murata, 605 - 00143 Rome, Italy
phone: +39 06 51860751
skype:  luigigeo25
email: luigi.vadacca at ingv.it