[CIG-LONG] Low viscosity material around dike

[Guillaume.Duclaux at csiro.au]

In the literature A is usually expressed in MPa^(-n).s^(-1), consequently it need to be scaled:

Agale= A^(1/n)*timeScalingFactor^(1/n)

I've been using values from Ranalli 1997, and for a quartz, A=1e-3 MPa^(-n).s^(-1), n=2, Ea=167 kJ.mol^(-1). No 1e-16 in my recent scripts, though.

Attached is a basic scaling spreadsheet I prepared for scaling the input parameters for the NonNewtonian viscosity implementation in Gale.
How it works: Fill the cells in orange (Edot, n, scalingFactor, A, Ea) and get T_0 and A_Gale (in red) scaled for the input file...
Use it at your own risks! I don't guarantee the correctness of the script. 
If you still want to use it and find an error, please let me know ;)

Cheers

Gilly

________________________________________________
 
Dr Guillaume Duclaux
CSIRO Earth Science and Resource Engineering
Visiting address: ARRC, 26 Dick Perry Av., Kensington WA 6151
Postal address: PO Box 1130, Bentley WA 6102, Australia
Ph: + 61 8 6436 8728    Fax: + 61 8 6436 8555    Web: www.csiro.au
 

-----Original Message-----
From: cig-long-bounces at geodynamics.org [mailto:cig-long-bounces at geodynamics.org] On Behalf Of Walter Landry
Sent: Saturday, 8 May 2010 7:18 AM
To: cig-long at geodynamics.org
Subject: Re: [CIG-LONG] Low viscosity material around dike

<Guillaume.Duclaux at csiro.au> wrote:
> Folks,
> 
> I'm currently out of office and I didn't bring my scaling scripts
> with me... The NewNewtonian rheology in Gale implies some scaling
> (in order to obtain an effective viscosity of the order of one), you
> can't directly use values from the literature.

The scaling for A is fairly straighforward.  The units for A are
1/viscosity, so you scale A inversely to what you scale viscosity.
What I do not know (and what I think Taichi is asking) is where did
the original numbers come from.  What is the original paper that you
got the numbers from?  Presumably people will want to use their own
numbers for their own runs, but it would be nice to know where the
example numbers come from.

Thanks,
Walter Landry
walter at geodynamics.org
_______________________________________________
CIG-LONG mailing list
CIG-LONG at geodynamics.org
http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long
-------------- next part --------------
A non-text attachment was scrubbed...
Name: scaleViscosity_Gale.xls
Type: application/vnd.ms-excel
Size: 37376 bytes
Desc: scaleViscosity_Gale.xls
Url : http://geodynamics.org/pipermail/cig-long/attachments/20100509/62c11292/attachment-0001.xls