[CIG-SHORT] Fault opening error
Birendra jha
bjha7333 at yahoo.com
Tue Aug 6 08:18:59 PDT 2013
Hi Charles
The curved segment of the fault smoothly bends approx. 17 deg clockwise from the segment that is parallel to the y_pos boundary. Is that angle too much?
I sent all the files to run the case in my first email on 8/2/13. If you have 15 minutes, I'd really appreciate if you could run the case and visualize it in Paraview. I think it runs as long as there is no slip which can be achieved with a high friction coefficient.
Thanks and regards
Birendra
--------------------------------------------
On Tue, 8/6/13, Charles Williams <C.Williams at gns.cri.nz> wrote:
Subject: Re: [CIG-SHORT] Fault opening error
To: cig-short at geodynamics.org
Date: Tuesday, August 6, 2013, 5:08 AM
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. 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
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.
-----Inline Attachment Follows-----
_______________________________________________
CIG-SHORT mailing list
CIG-SHORT at geodynamics.org
http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short
More information about the CIG-SHORT
mailing list