<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">
Hello all,
<div><br>
</div>
<div>Thank you for all the recent mails on the development of Aspect, it's really nice to be kept in the loop so. </div>
<div><br>
</div>
<div>Recently I reencountered a problem with the strain rate values in the top corners of our viscoplastic models as was earlier stated by one of our students (see messages below). It greatly complicates some of our models, but we hope we've found the cause. </div>
<div><br>
</div>
<div>The model is after Kaus 2010: it is a plastic medium with an inclusion at the bottom in extension (or compression). We enforce extension by using "Prescribed velocity boundary indicators: 0:function, 1:function" and then prescribing a certain horizontal
velocity and zero vertical velocity. The bottom boundary is free slip, the top boundary not specified.</div>
<div><br>
</div>
<div>What this should produce are two conjugate shears as in the first figure. What it actually produces is shown in the second figure; a broad shear zone reaching into the top corners plus the desired shears. When we prescribe extension velocities on the bottom
boundary as well, we get something close to what we want, but still with some strain rate anomaly in the upper corners (figure 3, different scaling). </div>
<div><br>
</div>
<div>Now, the first 2 figures are produced with another code, one that can prescribe boundary conditions per component of velocity: in figure 1 the vertical velocity component is free slip and only the horizontal component is prescribed. In figure 2 both are
prescribed (zero vertical velocity), as is the case with the Aspect model. We believe that because some material enters through the top near the corners with a vertical velocity component, but no vertical component is allowed along the vertical boundary, strain
rate is increased locally, triggering an additional shear band. </div>
<div><br>
</div>
<div>So, my question is, am I correct that right now you have to specify all velocity components when prescribing velocity and if so, would you be willing to change this?</div>
<div><br>
</div>
<div>Thank you for your help and hope to see you at AGU,</div>
<div><br>
</div>
<div>Anne</div>
<div><br>
<div><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><span class="Apple-style-span" style="border-collapse: separate; font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; ">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">
<span class="Apple-style-span" style="border-collapse: separate; font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; ">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">
<div><b style="color: rgb(51, 0, 153); "><span lang="NL" style="font-size: 7.5pt; "><font face="Verdana">Anne Glerum, MSc </font></span></b><a name="x_SafeHtmlFilter_13a4ab32afe5707b_OLE_LINK2" style="color: rgb(153, 0, 0); "><span lang="NL" style="font-size: 8pt; "><font face="Verdana">|</font></span></a><font face="Verdana" style="color: rgb(51, 0, 153); "><span><span lang="NL" style="font-size: 8pt; "> PhD
candidate</span></span><span><span lang="NL" style="font-size: 8pt; "> </span></span></font><a name="x_SafeHtmlFilter_13a4ab32afe5707b_OLE_LINK2" style="color: rgb(153, 0, 0); "><span lang="NL" style="font-size: 8pt; "><font face="Verdana">|</font></span></a><font face="Verdana" style="color: rgb(51, 0, 153); "><span><span lang="NL" style="font-size: 8pt; "></span></span></font><font face="Verdana" style="color: rgb(51, 0, 153); "><span><span lang="NL" style="font-size: 8pt; "></span></span><span><span lang="NL" style="font-size: 8pt; "></span></span><span><span lang="NL" style="font-size: 8pt; "> Department</span></span><span><span lang="NL" style="font-size: 8pt; "> of
Earth Sciences</span></span><span><span lang="NL" style="font-size: 8pt; "> </span></span></font><a name="x_SafeHtmlFilter_13a4ab32afe5707b_OLE_LINK2" style="color: rgb(153, 0, 0); "><span lang="NL" style="font-size: 8pt; "><font face="Verdana">|</font></span></a><font face="Verdana" style="color: rgb(51, 0, 153); "><span><span lang="NL" style="font-size: 8pt; "></span></span></font><font face="Verdana" style="color: rgb(51, 0, 153); "><span><span lang="NL" style="font-size: 8pt; "></span></span><span><span lang="NL" style="font-size: 8pt; "> </span></span><span><span lang="NL" style="font-size: 8pt; ">Utrecht
University</span></span><span><span lang="NL" style="font-size: 8pt; "></span></span></font><a name="x_SafeHtmlFilter_13a4ab32afe5707b_OLE_LINK2" style="color: rgb(153, 0, 0); "><span lang="NL" style="font-size: 8pt; "><font face="Verdana">|</font></span></a><font face="Verdana" style="color: rgb(51, 0, 153); "><span><span lang="NL" style="font-size: 8pt; "></span></span></font><font face="Verdana" style="color: rgb(51, 0, 153); "><span><span lang="NL" style="font-size: 8pt; "> </span></span><span><span lang="NL" style="font-size: 8pt; "> Budapestlaan
4</span></span><span><span lang="NL" style="font-size: 8pt; "></span></span><span><span lang="NL" style="font-size: 8pt; ">, 3584 CD Utrecht</span></span><span><span lang="NL" style="font-size: 8pt; "> </span></span></font><a name="x_SafeHtmlFilter_13a4ab32afe5707b_OLE_LINK2" style="color: rgb(153, 0, 0); "><span lang="NL" style="font-size: 8pt; "><font face="Verdana">|</font></span></a><font face="Verdana" style="color: rgb(51, 0, 153); "><span><span lang="NL" style="font-size: 8pt; "></span></span></font><font face="Verdana" style="color: rgb(51, 0, 153); "><span><span lang="NL" style="font-size: 8pt; "></span></span><span><span lang="NL" style="font-size: 8pt; "> </span></span><span><span lang="NL" style="font-size: 8pt; ">Room
Z.204</span></span><span><span lang="NL" style="font-size: 8pt; "> </span></span><span><span lang="NL" style="font-size: 8pt; "></span></span></font><a name="x_SafeHtmlFilter_13a4ab32afe5707b_OLE_LINK2"><span lang="NL" style="font-size: 8pt; "><font face="Verdana"><font class="Apple-style-span" color="#990000">|<span class="Apple-converted-space"> </span></font><font class="Apple-style-span" color="#3d00ff">A.C.Glerum@uu.nl</font><font class="Apple-style-span" color="#990000"> </font></font></span></a><a name="x_SafeHtmlFilter_13a4ab32afe5707b_OLE_LINK2" style="color: rgb(153, 0, 0); "><span lang="NL" style="font-size: 8pt; "><font face="Verdana">|</font></span></a></div>
<div style="color: rgb(0, 0, 0); "><br>
</div>
</div>
</span><img height="528" width="1188" apple-width="yes" apple-height="yes" id="76bea6d3-1d7f-4c3f-b7cb-8493ecd8a57d" src="cid:7B26DFB5-1AE7-4471-82C3-83E51BD04467@geo.uu.nl"></div>
</span><br class="Apple-interchange-newline">
</span><br class="Apple-interchange-newline">
</div>
<div><img height="549" width="1119" apple-width="yes" apple-height="yes" id="15f9dafd-1344-47ba-a08c-a78b14c7fce1" src="cid:5531CE59-9443-4D3C-8B2B-9380BA1B6623@geo.uu.nl"></div>
<div><img height="425" width="857" apple-width="yes" apple-height="yes" id="1c4c8f1e-565b-4713-9c40-0a6b9718d709" src="cid:54A2DE02-B0EE-448E-8AEA-AE1C079E0688@geo.uu.nl"></div>
<br>
<div>
<div>On Jul 1, 2013, at 1:36 AM, A.B.M. Graas wrote:</div>
<br class="Apple-interchange-newline">
<blockquote type="cite">
<div dir="ltr">Hi there,<br>
<div><br>
</div>
<div style="">No worries. Actually, I found out later that the problem is not caused by the material model but rather by strain rate values in the corners. This is more or less the same problem and should be easy to reproduce. An demonstration input file is
here: <a href="http://pastebin.com/5LCvM9fT">http://pastebin.com/5LCvM9fT</a>. Just tested it with the latest version. It might be necessary to alter the visual range of your viewer to find the strain rate values, since they are extremely small. However, in
other experiments they can become significantly large and even invalidate the set-up.</div>
<div style=""><br>
</div>
<div style="">Thanks for looking at this and let me know when you need any additional information,</div>
<div style=""><br>
</div>
<div style="">Adriaan</div>
<div style=""><br>
</div>
<div style=""><br>
</div>
<div style=""><br>
</div>
</div>
<div class="gmail_extra"><br>
<br>
<div class="gmail_quote">2013/6/30 Timo Heister <span dir="ltr"><<a href="mailto:heister@math.tamu.edu" target="_blank">heister@math.tamu.edu</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear Adriaan,<br>
<br>
sorry for net getting back to your email earlier.<br>
<div class="im"><br>
>>> - In many of my models I model plasticity by using a Material Model<br>
>>> implementing the new Interface. Often, though not always, there are<br>
>>> unexpected<br>
>>> irregularities (viscosity jumps) in the corners. I found that the<br>
>>> evaluate()<br>
>>> method of the Material Model is never called with a position Point in the<br>
>>> corners before the Stokes equation is solved. However, afterwards the<br>
>>> evaluate() method is called by the Viscosity postprocessor.<br>
>><br>
>><br>
>> You mean the point is never in a corner when assembling the linear<br>
>> systems, but it may be in a corner when using a postprocessor? That is<br>
>> correct -- we assemble linear systems using Gauss quadrature formulas where<br>
>> all quadrature points are in the interior of a cell. But we do output data<br>
>> using evaluation points that are the vertices of cells.<br>
>><br>
>> Can you explain how this leads to problems in your models?<br>
><br>
> If I understand this correctly the problem is limited to the visualization?<br>
> (here is an example VTU screenshot:<br>
> <a href="http://img196.imageshack.us/img196/5954/corners.png" target="_blank">http://img196.imageshack.us/img196/5954/corners.png</a>)<br>
<br>
</div>
Can you explain what the material model is and what you expect to see<br>
here? This image highlights some bug for sure (be it in your material<br>
model, the solver, or in the output).<br>
<div class="im"><br>
> What I intend to do is flipping on and off the "Nonlinear iteration"<br>
> parameter in the input file. However, even when "Nonlinear iteration" is<br>
> false, and "Nonlinear solver scheme" is set (f.e. to "iterated IMPES")<br>
> nonlinear iterations are executed. This is the case for box.prm and the<br>
> parameter file I send you.<br>
<br>
</div>
"Nonlinear iteration" was a left over parameter that did not have any<br>
effect. Therefore, I removed this parameter. You should use "Nonlinear<br>
solver scheme".<br>
<br>
Best,<br>
Timo<br>
<br>
--<br>
Timo Heister<br>
<a href="http://www.math.tamu.edu/~heister/" target="_blank">http://www.math.tamu.edu/~heister/</a><br>
<div class="HOEnZb">
<div class="h5">_______________________________________________<br>
Aspect-devel mailing list<br>
<a href="mailto:Aspect-devel@geodynamics.org">Aspect-devel@geodynamics.org</a><br>
<a href="http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel" target="_blank">http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel</a><br>
</div>
</div>
</blockquote>
</div>
<br>
</div>
_______________________________________________<br>
Aspect-devel mailing list<br>
<a href="mailto:Aspect-devel@geodynamics.org">Aspect-devel@geodynamics.org</a><br>
http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel</blockquote>
</div>
<br>
</div>
</body>
</html>