[CIG-SHORT] Fault opening error

[Birendra jha]

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