[CIG-LONG] Low viscosity material around dike

[Walter Landry]

Hi Taichi,

Taichi SATO <taichix at aori.u-tokyo.ac.jp> wrote:
> I also tried to run non-newtoinan viscosity model but it does not
> work. but it does not work. Please see attached file.
> the error is,
> Error in NonNewtonian: T, n, and A must all be non-zero:
>     T=0
>     n=3
>     A=1e+10

It turns out that if you do not evolve the temperature, Gale does not
apply the initial conditions for the temperature.  So you need to
uncomment the energy equation components.  If you want to fix the
temperature, then you need to set the BCs.

> I set the parameters for non-newtonian and mantle viscous as below
> <struct name="nonnewViscosity"> -->
> <param name="Type">NonNewtonian</param> -->
> <param
> name="StrainRateInvariantField">StrainRateInvariantField</param> -->
> <param name="TemperatureField">TemperatureField</param> -->
> <param name="n">3.0</param> --> <!--power law exponent-->
> <param name="T_0">1473</param> --> <!--T_0=Q/(n*R)=20045.37 -->
> <param name="A">1.0e10</param> --> <!--power law constant-->
> <param name="refStrainRate">1.0e22</param> -->
> <param name="minViscosity">1.0e19</param> -->
> <param name="maxViscosity">1.0e25</param> -->
> </struct> -->

You need to scale minViscosity and maxViscosity.  For your scaling I
believe it is 1e-3 and 1e3.

I have attached a new input file.  There were some odd comment
characters that I took out while trying to fix the file.  Let me know
if that works for you.

Cheers,
Walter Landry
walter at geodynamics.org
-------------- next part --------------
A non-text attachment was scrubbed...
Name: magma-tecto_new1_nonnewtoian_taichi.xml
Type: text/xml
Size: 35809 bytes
Desc: not available
Url : http://geodynamics.org/pipermail/cig-long/attachments/20100502/ab4e5d1d/attachment-0001.bin