[CIG-LONG] Low viscosity material around dike

Taichi SATO taichix at aori.u-tokyo.ac.jp
Mon May 3 00:38:46 PDT 2010


Hi  Walter,

Thank you very much!

This input file works but the outputs are something strange.  I found 
that the viscosity field does not change. This situation does not change 
even if I changed parameters in Non-newtonian section.

Sincerely,

Taichi

(2010/05/03 13:32), Walter Landry wrote:
> 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
>    
>
>
> _______________________________________________
> CIG-LONG mailing list
> CIG-LONG at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long
>    

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


More information about the CIG-LONG mailing list