<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
</style><style type="text/css" id="owaParaStyle"></style>
</head>
<body dir="ltr" fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">
<div><br>
</div>
Tony,
<div><br>
</div>
<div>It looks like your Helmholtz tolerance is too tight.</div>
<div><br>
</div>
<div>If your problem is scaled in convective time units with characteristic length scale 1, timescale 1, and visc ~ 1/Re,</div>
<div>then typical Helmholtz tolerances are around 1.e-7.</div>
<div><br>
</div>
<div>You might also try turning on projection ---</div>
<div><br>
</div>
<div>set p93=20 (and lmxprev=20 in SIZE), p94=5 and p95=5 in the rea file.</div>
<div><br>
</div>
<div>No guarantees at this point because there seems to be something flaky in the projection routines</div>
<div>these days under certain runtime combinations (e.g., I think splitting + heat transfer + </div>
<div>velocity projection does not work)... It's on my todo list to clean these up.</div>
<div><br>
</div>
<div>Note that your helmholtz solves are pretty well converged... residual is 1.7e-9, etc.</div>
<div>Projection should easily put you over the top</div>
<div><br>
</div>
<div>Paul</div>
<div><br>
</div>
<div>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF632791" style="direction: ltr;"><font face="Tahoma" size="2" color="#000000"><b>From:</b> nek5000-users-bounces@lists.mcs.anl.gov [nek5000-users-bounces@lists.mcs.anl.gov] on behalf of nek5000-users@lists.mcs.anl.gov [nek5000-users@lists.mcs.anl.gov]<br>
<b>Sent:</b> Thursday, January 29, 2015 9:50 AM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> [Nek5000-users] Error in Hmholtz solving<br>
</font><br>
</div>
<div></div>
<div>
<div id="divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#FFFFFF; font-family:Calibri,Arial,Helvetica,sans-serif">
<p>Dear Neks,</p>
<p><br>
</p>
<p>I'm getting the error message of '**ERROR**: Failed in HMHOLTZ: ' when running my simulation. I've checked the subroutine 'cggo' in hmholtz.f where the error message is printed. I understand that when the tolerance is not satisfied, the code would output
this error message. But what I'm wondering is why the maximum iteration number is only 100. In the subroutine 'cggo', niter is equal to min(maxit,maxcg) where maxit=200 and maxcg=900. I can't see other places where maxit is restricted to only 100. So why the
maximum iteration number is 100 when solving the Hmholtz equation?</p>
<p><br>
</p>
<p>Also, is it right that if I increase the maximum number of iterations to certain point where the tolerance is satisfied, I won't get the error message?</p>
<p><br>
</p>
<p>Hope anyone can help me on this. Thank you very much in advance.</p>
<p><br>
</p>
<p>Best regards,<br>
</p>
<div id="Signature">
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
Tony<br>
<br>
<br>
Example of the error message:<br>
<br>
Step 52, t= 2.5600000E+01, DT= 5.0000000E-02, C= 2.689 5.6747E+02 1.0237E+01 <br>
Solving for fluid T T T <br>
52 PRES gmres: 29 9.7057E-07 1.0000E-06 1.5219E-03 1.0926E+00 2.1480E+00 F<br>
52 Hmholtz VELX: 99 8.5954E-10 6.0898E-03 1.0000E-09 <br>
52 100 **ERROR**: Failed in HMHOLTZ: VELY 1.7088E-09 2.5049E-02 1.0000E-09 <br>
52 100 **ERROR**: Failed in HMHOLTZ: VELZ 1.4665E-09 2.5165E-02 1.0000E-09 <br>
L1/L2 DIV(V) : 2.7573E-18 9.1661E-04 <br>
L1/L2 QTL : 0.0000E+00 0.0000E+00 <br>
L1/L2 DIV(V)-QTL: 2.7573E-18 9.1661E-04 <br>
52 2.5600E+01 1.0094E+01 Fluid done
<br>
Step 53, t= 2.5650000E+01, DT= 5.0000000E-02, C= 2.697 5.7766E+02 1.0186E+01 <br>
Solving for fluid T T T <br>
53 PRES gmres: 28 9.6076E-07 1.0000E-06 1.1322E-03 1.0558E+00 2.0542E+00 F<br>
53 Hmholtz VELX: 98 9.4891E-10 5.8783E-03 1.0000E-09 <br>
53 100 **ERROR**: Failed in HMHOLTZ: VELY 1.6450E-09 2.4941E-02 1.0000E-09 <br>
53 100 **ERROR**: Failed in HMHOLTZ: VELZ 1.4160E-09 2.5049E-02 1.0000E-09 <br>
L1/L2 DIV(V) : 2.4695E-18 1.0845E-03 <br>
L1/L2 QTL : 0.0000E+00 0.0000E+00 <br>
L1/L2 DIV(V)-QTL: 2.4695E-18 1.0845E-03 <br>
53 2.5650E+01 9.9794E+00 Fluid done <br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>