[CIG-LONG] CIG-LONG Digest, Vol 45, Issue 12
Sephiroth AvALoN
sephiroth42 at gmail.com
Sat Aug 21 12:24:53 PDT 2010
Hi,
For the non-newtonian rheology, after checking at the unit of A in gale, it
seems that the unit is not in Pa^(-1) . s^(-1/n) as indicated in gale manual
but in Pa^(-n) . s(-1) (Same as ranalli, 1997).
Indeed when using translated values from Frederiksen et al., 2001, we obtain
A before scaling = 1.63e-26, then we have to scale directly with the scaling
factor for gale (~1e25). Ascaled becomes 0.163 and works just fine !
[A translated from MPa . s^(1/n) to Pa^(-n) . s(-1) by multiplying *^(-n)
----> 208 MPa . s^(1/n) becomes 208e6^(-3.1) = 1.63e-26]
Cheers
Nicolas Riel
On Wed, Aug 18, 2010 at 9:00 PM, <cig-long-request at geodynamics.org> wrote:
> Send CIG-LONG mailing list submissions to
> cig-long at geodynamics.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long
> or, via email, send a message with subject or body 'help' to
> cig-long-request at geodynamics.org
>
> You can reach the person managing the list at
> cig-long-owner at geodynamics.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of CIG-LONG digest..."
>
>
> Today's Topics:
>
> 1. Wrong "Piece Extend" coordinates on field*.pvtu files
> (Patrice Rey)
> 2. Re: Wrong "Piece Extend" coordinates on field*.pvtu files
> (Guillaume.Duclaux at csiro.au)
> 3. Non-Newtonian rheology scaling ! (Nicolas RIEL)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Wed, 18 Aug 2010 14:05:06 +1000
> From: Patrice Rey <patrice.rey at me.com>
> Subject: [CIG-LONG] Wrong "Piece Extend" coordinates on field*.pvtu
> files
> To: cig-long at geodynamics.org
> Message-ID: <B424876C-0F6B-444F-8AEE-8E3A681E5BAF at me.com>
> Content-Type: text/plain; charset=us-ascii
>
> G'day all
>
> I have noticed that, when running Gale-1.5 on a cluster, the field*.pvtu
> files are coming out with the wrong "Piece Extend" coordinates. The vtu
> files called in .pvtu files have the same coordinates.
>
> ...
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.0.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.1.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.2.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.3.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.4.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.5.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.6.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.7.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.8.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.9.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.10.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.11.00000.vts"/>
> ....
>
> This makes loading in Paraview a bit strenuous as one has to upload all
> individual *.vtu files, instead of calling the *.pvtu.
>
> The solution I am using at the moment is to edit the .pvts file, ok when
> running on a few CPUs, not so ok when running on many.
>
> Interestingly, this does not happen when I run the same input script on the
> 8 cores of my MacPro.
>
> I was wondering if other users have experienced the same problem?
>
> Cheers
>
> Patrice
>
>
>
> Associate Professor Patrice F. Rey
> EarthByte Research Group
> School of Geosciences
> The University of Sydney
> +61 2 9351 2067
> patrice.rey at mac.com
> patrice.rey at sydney.edu.au
>
>
>
>
>
>
>
>
> ------------------------------
>
> Message: 2
> Date: Wed, 18 Aug 2010 15:16:19 +0800
> From: <Guillaume.Duclaux at csiro.au>
> Subject: Re: [CIG-LONG] Wrong "Piece Extend" coordinates on
> field*.pvtu files
> To: <patrice.rey at me.com>
> Cc: cig-long at geodynamics.org
> Message-ID: <00A83AFA-5C48-481F-9663-900FB0D3C5EC at csiro.au>
> Content-Type: text/plain; charset="us-ascii"
>
> Patrice,
>
> I had a similar issue with Underworld few months ago... and similarly to
> what you describe this occurred on one machine only..
> I reported this to Luke Hodkinson in Monash but he didn't observe this
> 'feature' with his install, all was working fine!
> My guess is that it has something to do with the version of VTK on the
> machine.
> Workaround with Underworld: use another server or switch to the xdmf
> outputs for visualisation...
>
>
>
> On 18/08/2010, at 12:05 PM, Patrice Rey wrote:
>
> G'day all
>
> I have noticed that, when running Gale-1.5 on a cluster, the field*.pvtu
> files are coming out with the wrong "Piece Extend" coordinates. The vtu
> files called in .pvtu files have the same coordinates.
>
> ...
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.0.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.1.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.2.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.3.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.4.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.5.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.6.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.7.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.8.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.9.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.10.00000.vts"/>
> <Piece Extent="0 28 0 21 0 0"
> Source="fields.11.00000.vts"/>
> ....
>
> This makes loading in Paraview a bit strenuous as one has to upload all
> individual *.vtu files, instead of calling the *.pvtu.
>
> The solution I am using at the moment is to edit the .pvts file, ok when
> running on a few CPUs, not so ok when running on many.
>
> Interestingly, this does not happen when I run the same input script on the
> 8 cores of my MacPro.
>
> I was wondering if other users have experienced the same problem?
>
> Cheers
>
> Patrice
>
>
>
> Associate Professor Patrice F. Rey
> EarthByte Research Group
> School of Geosciences
> The University of Sydney
> +61 2 9351 2067
> patrice.rey at mac.com<mailto:patrice.rey at mac.com>
> patrice.rey at sydney.edu.au
>
>
>
>
>
>
> _______________________________________________
> CIG-LONG mailing list
> CIG-LONG at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long
>
>
> Dr Guillaume Duclaux
>
> Mineral Down Under Flagship & AuScope Grid
> CSIRO Earth Science and Resource Engineering
>
> Phone: +61 8 6436 8728 | Fax: +61 8 6436 8559 | Mobile: +61 422 289 732
>
> guillaume.duclaux at csiro.au<mailto:guillaume.duclaux at csiro.au> |
> www.csiro.au<http://www.csiro.au/> |
> Address: Australian Resources Research Centre, 26 Dick Perry Avenue,
> Kensington WA 6151
>
> PLEASE NOTE
> The information contained in this email may be confidential or privileged.
> Any unauthorised use or disclosure is prohibited. If you have received this
> email in error, please delete it immediately and notify the sender by return
> email. Thank you. To the extent permitted by law, CSIRO does not represent,
> warrant and/or guarantee that the integrity of this communication has been
> maintained or that the communication is free of errors, virus, interception
> or interference.
>
> Please consider the environment before printing this email.
>
>
>
>
>
> ------------------------------
>
> Message: 3
> Date: Wed, 18 Aug 2010 17:20:19 +0200
> From: Nicolas RIEL <nicolas.riel at free.fr>
> Subject: [CIG-LONG] Non-Newtonian rheology scaling !
> To: cig-long at geodynamics.org
> Message-ID:
> <AANLkTimF7AA1jjOfbcz+QwFTXTpJp-uunMLPp+Jh2Fjn at mail.gmail.com<AANLkTimF7AA1jjOfbcz%2BQwFTXTpJp-uunMLPp%2BJh2Fjn at mail.gmail.com>
> >
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hi,
>
> I'm encountering some problems when scaling the "A" (creeping parameter)
> for
> the non-newtonian rheology.
> In gale user manual (2.2.8.2) it is said that the unit of A parameter is *
> s(-1/n).Pa(-1)
> *For other authors the unit of this parameter is most of the time in *
> Pa(-n).s(-1)* or *Pa.s(1/n).*
>
> I will take an example in order to understand :
>
>
> A = 262800* Pa.s(1/n)*
> n = 4*
>
> *If I want to translate this value for gale, first I have to change the
> units for *s(-1/n).Pa(-1) *by multplying with ^(-1) ?
>
> Then A unscaled becomes for gale 262800^(-1) = 3.8051e-6
>
> Based on that and following gale user manual A scaled for gale is 3.8051e-6
> * (10e25)^(1/n) = 6.7574
> However in the dyke example the scaling of A is around 1e-12.
>
> I'm wondering if there a mistake for the given scaling equation. Maybe this
> equation should be Anew = Aold(10e25)^(-1/n) ?
>
> When I use this, I obtain Anew = 2.1368e-12 and this works
>
> Do you have any suggestions about this ?
>
> Nicolas Riel
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://geodynamics.org/pipermail/cig-long/attachments/20100818/a12b49f2/attachment.html
>
> ------------------------------
>
> _______________________________________________
> CIG-LONG mailing list
> CIG-LONG at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long
>
>
> End of CIG-LONG Digest, Vol 45, Issue 12
> ****************************************
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/cig-long/attachments/20100821/9a6c4d79/attachment.htm
More information about the CIG-LONG
mailing list