[CIG-SHORT] pylith install error
Charles Williams
willic3 at gmail.com
Tue May 27 21:14:11 PDT 2014
Hi Birendra,
As far as point #1, it appears you are replacing the overburden with an equivalent traction load. That should be OK for infinitesimal strain problems, I think.
Did changing your normalizer length scale influence your convergence?
Cheers,
Charles
On 28/05/2014, at 12:28 pm, Birendra jha <bjha7333 at yahoo.com> wrote:
> 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
> _______________________________________________
> 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
More information about the CIG-SHORT
mailing list