[CIG-SHORT] Fault opening error

Mon Aug 5 16:38:18 PDT 2013

Sorry about that, Birendra.  I didn't actually look at the test before sending it.  It seems to me this one may not be what you need.  The examples/2d/subduction folder may be more what you need.

How curved is the fault in your problem?

Cheers,
Charles

On 6/08/2013, at 11:09 AM, Birendra jha wrote:

> Thank you Charles. Two comments/questions:
> 1. There is nothing about solver settings in the pylithapp.cfg, no petsc settings either.
> 2. Density is specified as 2500.0*m/2.
> 
> regards
> Birendra
> 
> --------------------------------------------
> On Tue, 8/6/13, Charles Williams <willic3 at gmail.com> wrote:
> 
> Subject: Re: [CIG-SHORT] Fault opening error
> To: cig-short at geodynamics.org
> Date: Tuesday, August 6, 2013, 2:53 AM
> 
> Hi Birendra,
> It's possible the tests don't get
> distributed with the binaries.  I'm attaching a
> tarball with the directory Brad mentioned.
> Cheers,Charles
> 
> On 6/08/2013, at 5:48 AM, Birendra jha
> wrote:
> Hi Brad
> 
> I did not find the tests/2d/nonplanar folder in the pylith
> directory. Did you mean examples/2d/subduction?
> 
> Also, can you please send me the pylithapp.cfg with modified
> tolerances?
> 
> Thanks a lot
> Birendra
> 
> --------------------------------------------
> On Mon, 8/5/13, Birendra jha <bjha7333 at yahoo.com>
> wrote:
> 
>  Subject: Re: [CIG-SHORT] Fault opening error
>  To: "Brad Aagaard" <baagaard at usgs.gov>
>  Date: Monday, August 5, 2013, 8:16 PM
> 
>  Brad
> 
>  Is it possible for you to send me the pylithapp.cfg, if
>  that's the one you modified?
> 
>  I changed my ksp_rtol and snes_rtol to 1e-20 and I am
> still
>  getting the error.
> 
>  Thank you
>  Birendra
> 
>  --------------------------------------------
>  On Mon, 8/5/13, Brad Aagaard <baagaard at usgs.gov>
>  wrote:
> 
>   Subject: Re: [CIG-SHORT] Fault opening error
>   To: "Birendra jha" <bjha7333 at yahoo.com>
>   Cc: cig-short at geodynamics.org
>   Date: Monday, August 5, 2013, 7:44 PM
> 
>   Birendra,
> 
>   I ran a few time steps with the tighter tolerances
> and
>   didn't get an 
>   error. I am out of the office for another week and
> don't
>   have time to 
>   run anything until I return.
> 
>   There is a curved fault dynamic simulation test in
>   tests/2d/nonplanar.
> 
>   Regards,
>   Brad
> 
> 
>   On 8/5/13 7:11 AM, Birendra jha wrote:
> Hi Brad
> 
> I did that. It did not
> solve the fault opening
>   problem.
> 
> Do you have some time
> to run this case? It should
>  take
>   1-2 minutes. The main problem is that when slip is
>  non-zero
>   on the curved fault surface, slip vector on one of
> these
>   nodes is assigned some non zero value in the normal
>   direction during numerical iteration. Since fault is
> under
>   compression everywhere, it gives an error from
>   FaultCohesiveDyn::integrateResidual. Is it possible
> that
>   after the iterations have converged, there is no
> slip
>   component in the normal direction? If so, then this
> error
>   check in inetgrateResidual is probably too strict?
> 
> Do you have a simple
> case where spontaneous fault
>  slip
>   is modeled on a curved fault surface?
> 
> Thanks and best
> regards
> Birendra
> 
> 
> --------------------------------------------
> On Sun, 8/4/13, Brad
> Aagaard <baagaard at usgs.gov>
>   wrote:
> 
>    Subject: Re: [CIG-SHORT] Fault
> opening
>   error
>    To: cig-short at geodynamics.org
>    Date:
> Sunday, August 4, 2013, 8:23 PM
> 
>    Birendra,
> 
>    Your
> solver tolerances need
>   adjustment. The FaultCohesiveDyn
>    object in a
>    quasi-static simulation
> requires use
>   of absolute tolerances
>    in order to
>    properly detect stick/slip
> with the
>   iterative solver (see
>    the friction
>    session in the CDM2013
> tutorial for
>   details). You should
>    lower the
>    ksp_rtol and snes_rtol
> tolerances to
>   1.0e-20 to force use of
>    the
>    absolute tolerances.
> 
>    Brad
> 
>    On
> 8/2/13 2:35 PM, Birendra jha
>   wrote:
>    > Hi everyone,
>    >
>    >
> I am trying to run a model with a
>   fault
>    (FaultCohesiveDyn)
> in it. The fault
>   surface is not planar
>    but is curved. I
> am getting fault
>   opening error at some
>    node, which stops
> the simulation. Why
>   does this happen? What
>    should I do to
> avoid getting this
>   error?
>    >
>    >
> Please see the attached files
>   which should run.
>    >
>    >
> Here is what I have tried so
>   far:
>    > 1. different
> fault dips
>    >
> 2. different fault traces (plane
>   to curved)
>    > 3. different
> fault friction
>   coefficients
>    >
>    >
> My objective is to get some
>   reverse slip on this curved
>    fault surface. The
> slip does not have
>   to be on the curved
>    part, it can be on
> the part that
>   strikes parallel to the
>    y_pos boundary.
>    >
>    >
> Thank you very much for any
>   help.
>    >
>    >
> Birendra
>    >
>    >
>    >
>    >
> 
>  _______________________________________________
>    > CIG-SHORT
> mailing list
>    >
> CIG-SHORT at geodynamics.org
>    >
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>    >
> 
>    _______________________________________________
>    CIG-SHORT mailing list
>    CIG-SHORT at geodynamics.org
>    http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short
> 
> 
> 
> 
> 
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short
> 
> 
> Charles A. WilliamsScientist
> GNS Science
> 1 Fairway Drive, Avalon
> PO Box 30368
> Lower Hutt  5040
> New Zealand
> ph (office): 0064-4570-4566
> fax (office): 0064-4570-4600C.Williams at gns.cri.nz
> 
> 
> 
> -----Inline Attachment Follows-----
> 
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short

Charles A. Williams
Scientist
GNS Science
1 Fairway Drive, Avalon
PO Box 30368
Lower Hutt  5040
New Zealand
ph (office): 0064-4570-4566
fax (office): 0064-4570-4600
C.Williams at gns.cri.nz

Notice: This email and any attachments are confidential.
If received in error please destroy and immediately notify us.
Do not copy or disclose the contents.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://geodynamics.org/pipermail/cig-short/attachments/20130806/662f8ecc/attachment.html>