[CIG-SHORT] pylith install error

Charles Williams willic3 at gmail.com
Sun May 25 19:42:49 PDT 2014 

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

More information about the CIG-SHORT mailing list