[CIG-LONG] Gale: NonNewtonian viscosity

[Guillaume.Duclaux at csiro.au]

Hi Walter,

please find the input file attached. Min and maxViscosity are commented. I ran this example on the win32 binary (v1.4.1)

Cheers

Gilly 
 


Dr Guillaume Duclaux
Structural Geologist/Modeller
CSIRO Exploration and Mining
Visiting address: 6G13, ARRC, 26 Dick Perry Avenue, Kensington, Western Australia, 6151
Postal address: PO Box 1130, Bentley, Western Australia, 6102, Australia
Ph: + 61 8 6436 8728  Fax: + 61 8 6436 8555  Web: www.em.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, 7 November 2009 5:25 AM
To: cig-long at geodynamics.org
Subject: Re: [CIG-LONG] Gale: NonNewtonian viscosity

<Guillaume.Duclaux at csiro.au> wrote:
> 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.

For some reason I never committed the fix for this.  Done (rev 15938).

> 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. 

Are you sure?  I am looking at the code, and it seems correct.  Do you
have an input file?

Cheers,
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: NonNewtonianTest.xml
Type: text/xml
Size: 28710 bytes
Desc: NonNewtonianTest.xml
Url : http://geodynamics.org/pipermail/cig-long/attachments/20091107/8904d819/attachment-0001.bin