[CIG-LONG] Gale: NonNewtonian viscosity
Walter Landry
walter at geodynamics.org
Fri Nov 6 13:24:41 PST 2009
<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
More information about the CIG-LONG
mailing list