[CIG-LONG] Gale: NonNewtonian viscosity

[Guillaume.Duclaux at csiro.au]

Hi all,

Few month ago I pointed out a small error in the encoding of the NonNewtonian viscosity to Walter. I just paste below a part of Walter's answer to my message.

> the line

viscosity=exp(T_0/(2*n*T))/(2*A);

has an extraneous factor of 2, so it should be:

viscosity=exp(T_0/(n*T))/(2*A);

Also, to correct this small feature, one can multiply T_0 with a factor of 2 in the input file.

I found today that the minViscosity and maxViscosity properties are inverted. (Never used these properties before...)
Indeed, minViscosity applies an upper threshold (max viscosity) to the effective Viscosity, and maxViscosity a lower one (min Viscosity). I guess it should the other way around.

Hope that helps.

Cheers

Gilly

________________________________________________

Dr Guillaume Duclaux
Structural Geologist /  Modeller
CSIRO Exploration and Mining
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<http://www.csiro.au/>


-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/cig-long/attachments/20091106/331a5d82/attachment.htm