[CIG-SHORT] pylith install error
Birendra jha
bjha7333 at yahoo.com
Mon May 26 12:03:56 PDT 2014
Hi Charles
Thanks for the input. I also increased the slip-weakening-parameter, i.e. the critical slip from 5 mm to 5 m and now it converges immediately to zero solution. I had used the critical slip of 5 mm for a much finer mesh of the same problem and didn't raise it for this coarse mesh.
I also decreased the ksp_rtol from 1e-10 to 1e-15, snes_rtol from 1e-10 to 1e-12, and snes_atol from 1e-9 to 1e-10.
Thanks for the help.
Birendra
--------------------------------------------
On Mon, 5/26/14, Charles Williams <willic3 at gmail.com> wrote:
Subject: Re: [CIG-SHORT] pylith install error
To: cig-short at geodynamics.org
Date: Monday, May 26, 2014, 8:12 AM
Hi Birendra,
Your gravitational stresses,
initial stresses, and traction BC are balancing now (you can
see this by turning off the fault). As far as why the
fault solution isn’t converging, I’m not sure. I
haven’t been through your parameters in detail.
Actually, I just looked again at your parameters, and
noticed you have the normalizer length scale set to 1 m,
which is much smaller than any of your element sizes. Try
increasing that to 100 m.
Cheers,
Charles
On
26/05/2014, at 1:36 pm, Birendra jha <bjha7333 at yahoo.com>
wrote:
> Hi Charles
>
> So if I set the
initialstress = 0 at z = -500 m and initialstress =
rho*g*2000 at z = -2500 m, and modify the topboundary and
rightboundary spatialdb accordingly (files attached), I
should get zero displacement solution? I don't get that.
>
> Thanks and
regards
> Birendra
>
--------------------------------------------
> On Mon, 5/26/14, Charles Williams <willic3 at gmail.com>
wrote:
>
> Subject:
Re: [CIG-SHORT] pylith install error
>
To: cig-short at geodynamics.org
> Date: Monday, May 26, 2014, 5:25 AM
>
> Hi Birendra,
>
> The gravitational
stresses are
> independent of any
vertical coordinate shift, which means
>
stresses at the surface should still be zero. If you
want
> to represent stresses at depth,
you will need to include the
> overlying
material.
>
>
Cheers,
> Charles
>
>
> On
> 26/05/2014, at 10:27 am, Birendra jha
<bjha7333 at yahoo.com>
> wrote:
>
>> Hi Charles
>>
>> Yes, I have
>
prescribed those three types of stresses, and I would
like
> to get a zero displacement
solution by balancing them out.
>
>> The top depth of the model is 500
m.
> Therefore, I have a non-zero
traction at the top boundary to
> balance
the initial stress, which is also non-zero at the
> top i.e. at the 500 m depth. On the right,
I am trying to
> represent a horizontal
compression which is 0.7 times the
>
vertical stress (overburden), i.e. a normal faulting
> regime.
>>
>> Does this
> make
sense?
>>
>>
> Thanks for looking into it
>> Birendra
>>
>
--------------------------------------------
>> On Mon, 5/26/14, Charles Williams
<willic3 at gmail.com>
> wrote:
>>
>> Subject:
> Re:
[CIG-SHORT] pylith install error
>>
> To: cig-short at geodynamics.org
>> Date: Monday, May 26, 2014, 3:43 AM
>>
>> Hi
Birendra,
>>
>>
I’m having
> difficulty
>> figuring out what you are
> trying to represent here. It
>> appears
> that you
are using:
>>
>> 1. Gravitational stresses.
>> 2. Initial stresses (presumably
equal to
> the
>>
gravitational stresses).
>> 3.
Normal
>> tractions
> on the top and right boundaries.
>>
>> Is this
correct?
>>
>>
Actually, looking a bit more,
>> these
stresses don’t seem to balance.
> You
have nonzero
>> normal stresses on
> your top boundary (they should be
zero),
>> as well as at the top of
your right
> side. Also, you have
>> nonzero initial
>
stresses at the top of your mesh, which
>>
> won’t balance
with gravitational stresses. It appears
>> you have an inconsistent stress
state.
>>
>>
Cheers,
>> Charles
>>
>>
>> On
>>
22/05/2014, at 3:24 pm, Birendra jha
>
<bjha7333 at yahoo.com>
>> wrote:
>>
>>> Hello Pylith
>>
> developers and
users,
>>>
>>> Can you please help me with this
issue
> of
>>
non-convergence of the nonlinear
>
solver? I attached the
>>
simulations
> files in my earlier email.
Why is pylith not
>> able to converge
to zero displacement
> solution at t=0
when
>> initial stresses,
> boundary tractions, and the body forces
>>
> are supposed to
balance out? I am using pylith-1.9.0
>> installed from the binary.
>>>
>>>
Thanks and
> best regards
>>> Birendra
>>>
>>
>
--------------------------------------------
>>> On Wed, 5/21/14, Birendra jha
<bjha7333 at yahoo.com>
>> wrote:
>>>
>>> Subject:
>> Re:
> [CIG-SHORT]
pylith install error
>>>
>
>> To: cig-short at geodynamics.org
>>> Date: Wednesday, May 21, 2014,
8:07
> PM
>>>
>>> Hi
>
Charles, Matt,
>>>
>>> I tried running the
>> attached problem in pylith-1.9.0
and
> the
>>>
nonlinear solver does not
> converge at
t=0.
>> This problem has
>>> a dynamic fault but
>> there should be no slip at t=0
because
>>>
>> the initial
>
stresses, boundary tractions, and gravity
>> should
>>>
result in
> zero displacement at
>> t=0. I don't
>
understand why I
>>>
>> don't get it. Can you please take
a
> look? This runs in
>> few
>>>
seconds.
>>>
>>> Thanks and regards
>>>
>>
Birendra
>>>
>>>
>>>
>>
>
--------------------------------------------
>>> On Wed, 5/21/14, Charles
Williams
> <willic3 at gmail.com>
>>> wrote:
>>>
>>>
Subject: Re: [CIG-SHORT]
> pylith
install
>> error
>>> To: cig-short at geodynamics.org
>>> Date: Wednesday, May 21,
> 2014, 1:56
>> PM
>>>
>
>>> Something is
>> strange.
>>> Do you not have these
>> symbols in libdl?
>>>
>
>>> Also, why are you
> using pylith-1.8.0?
>>>
>>>
Charles
>>>
>>>
>>>
On 21/05/2014, at 8:00
> pm, Birendra
jha
>> <bjha7333 at yahoo.com>
>>> wrote:
>>>
>>>>
Hi
>>>>
>>>> I am
>
trying to
>>>
>>
> install
pylith-1.8.0 and I am getting the following error
>>> during make:
>>>>
>>
>>>>
>>> /bin/bash
>> ../../libtool
>
--tag=CXX --mode=link
>>> g++ -g -O2
> -DMPICH_IGNORE_CXX_SEEK
>> -ldl
>>>
-L/home/bjha/pylith-1.8.0/lib
>
>> -o DllPlugInTester
>>> DllPlugInTester.o
>> CommandLineParser.o
>>>
>>
>
../../src/cppunit/libcppunit.la -lm
>>>>
>>
>>> mkdir .libs
>>>>
>> g++
-g
> -O2
>>>
-DMPICH_IGNORE_CXX_SEEK
> -o
>> .libs/DllPlugInTester
>>>
>>
>
DllPlugInTester.o CommandLineParser.o -ldl
>>>
-L/home/bjha/pylith-1.8.0/lib
>>>
../../src/cppunit/.libs/libcppunit.so
>> -lm -Wl,--rpath
>>>
>>
>
-Wl,/home/bjha/pylith-1.8.0/lib
>>>>
>
>>>
../../src/cppunit/.libs/libcppunit.so:
>> undefined reference
>>> to `dlsym'
>>>>
>>>
>>
../../src/cppunit/.libs/libcppunit.so:
>
undefined
>> reference
>>> to `dlopen'
>>>>
>>>
>>
../../src/cppunit/.libs/libcppunit.so:
>
undefined
>> reference
>>> to `dlclose'
>>>> collect2: ld returned
>>> 1 exit status
>>>>
>>
make[5]:
> ***
>>> [DllPlugInTester]
> Error
>> 1
>>>>
>
make[5]:
>>>
>>
> Leaving
directory
>>>
>>
>
`/home/bjha/build/pylith-1.8.0/cppunit-build/src/DllPlugInTester'
>>>> make[4]: ***
[all-recursive]
> Error
>> 1
>>>>
> make[4]: Leaving directory
>>>
>>
>
`/home/bjha/build/pylith-1.8.0/cppunit-build/src'
>>>> make[3]: ***
[all-recursive]
> Error
>> 1
>>>>
> make[3]: Leaving directory
>>>
>>
>
`/home/bjha/build/pylith-1.8.0/cppunit-build'
>>>> make[2]: *** [cppunit] Error
2
>>>> make[2]: Leaving
directory
>>>
`/home/bjha/build/pylith-1.8.0'
>>>>
>>> make[1]: ***
>> [installed_cppunit] Error 2
>>>>
>>
make[1]:
> Leaving directory
>>>
>> `/home/bjha/build/pylith-1.8.0'
>>>>
>>
>>> make: ***
> [installed_pylith] Error
>> 2
>>>>
>>>
>>>>
>>>> I
>
have
>> attached
>>> the log file.
>>>>
>>>>
>>> Thanks
>>>>
>>>
>>
>
Birendra<config.log>_______________________________________________
>>>> CIG-SHORT mailing list
>>>> CIG-SHORT at geodynamics.org
>>>> http://lists.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
>>>
>>>
>>
>
_______________________________________________
>>> 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
>>
>> 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
>>
>>
>
_______________________________________________
>> 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
>
> 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
>
>
_______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short<initialstress.spatialdb><rightboundarytract.spatialdb><topboundarytract.spatialdb>_______________________________________________
>
CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.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
_______________________________________________
CIG-SHORT mailing list
CIG-SHORT at geodynamics.org
http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
More information about the CIG-SHORT
mailing list