[CIG-SHORT] pylith install error
Birendra jha
bjha7333 at yahoo.com
Tue May 27 17:28:26 PDT 2014
Hi Charles
Few things I want to share with the community:
1. Please correct me on this, but I found that it is not necessary to set zero load at the top boundary, which is at a depth of 500 m. I can represent stresses at depth without including all the overburden in my domain. So I can keep the original initialstress, topboundary, and rightboundary spatialdb files which had overburden at 500 m.
2. It was also not necessary to increase my critical slip parameter from 5 mm to 5 m. However, I understand its importance once the fault starts to slip.
3. The only crucial part in getting the nonlinear solver converged to the zero solution was to decrease the ksp tolerances to the order of 1e-20.
Thanks
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