<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Hi Walter,<br>
<br>
Thank you very much!<br>
<br>
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.<br>
<br>
Sincerely,<br>
<br>
Taichi <br>
<br>
(2010/05/03 13:32), Walter Landry wrote:
<blockquote
cite="mid:20100502.213224.513173052708537441.walter@geodynamics.org"
type="cite">
<pre wrap="">Hi Taichi,
Taichi SATO <a class="moz-txt-link-rfc2396E" href="mailto:taichix@aori.u-tokyo.ac.jp"><taichix@aori.u-tokyo.ac.jp></a> wrote:
</pre>
<blockquote type="cite">
<pre wrap="">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
</pre>
</blockquote>
<pre wrap="">
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.
</pre>
<blockquote type="cite">
<pre wrap="">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> -->
</pre>
</blockquote>
<pre wrap="">
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
<a class="moz-txt-link-abbreviated" href="mailto:walter@geodynamics.org">walter@geodynamics.org</a>
</pre>
<pre wrap="">
<fieldset class="mimeAttachmentHeader"></fieldset>
_______________________________________________
CIG-LONG mailing list
<a class="moz-txt-link-abbreviated" href="mailto:CIG-LONG@geodynamics.org">CIG-LONG@geodynamics.org</a>
<a class="moz-txt-link-freetext" href="http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long">http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long</a>
</pre>
</blockquote>
<br>
</body>
</html>