<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p>Thank you, this solved the problem!</p>
<p>Now all scalars behave identically as they should and no crash.</p>
<p><br>
</p>
<p>How should I understand the reason for the previous crash?<br>
</p>
<br>
<br>
<div style="color: rgb(0, 0, 0);">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" 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> Sunday, May 14, 2017 7:30 PM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> Re: [Nek5000-users] Crash when implementing additional scalars</font>
<div> </div>
</div>
<div>
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">
<div><br>
</div>
<div>Dear Johan,</div>
<div><br>
</div>
Try setting param 94=0 in your .rea file.
<div><br>
</div>
<div>Best,<br>
<div><br>
</div>
<div>Paul</div>
<div><br>
<div><br>
</div>
<div><br>
<div style="font-family:Times New Roman; color:#000000; font-size:16px">
<hr tabindex="-1">
<div id="divRpF208107" 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> Sunday, May 14, 2017 12:26 PM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> [Nek5000-users] Crash when implementing additional scalars<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 get the following crash when implementing 3 scalar fields (temperature + 2 additional scalar fields):</p>
<p><br>
</p>
<p>**************<br>
</p>
<br>
<p>Step 4, t= 4.0000000E-03, DT= 1.0000000E-03, C= 0.331 7.1728E+01 2.2114E+01<br>
Solving for heat<br>
Solving for fluid F T T<br>
4 2.51148E-02 6.07580E-04 4.13357E+01 1 TEMP PROJECT<br>
4 Hmholtz TEMP: 4 3.0551E-08 8.0302E-04 1.0000E-07<br>
4 2.51148E-02 1.60663E-08 1.56319E+06 2 PS 1 PROJECT<br>
4 Hmholtz PS 1: 0 3.0509E-08 3.0509E-08 1.0000E-07<br>
4 2.51148E-02 3.38605E-02 7.41712E-01 3 PS 2 PROJECT<br>
4 100 **ERROR**: Failed in HMHOLTZ: PS 2 1.6090E+10 5.4264E-02 1.0000E-07<br>
4proj_ortho: 3 4 PS 2 Detect rank deficiency: NaN 1.0000E+00<br>
4proj_ortho: 2 4 PS 2 Detect rank deficiency: NaN 1.0000E+00<br>
4proj_ortho: 1 4 PS 2 Detect rank deficiency: NaN 1.0000E+00<br>
4 4.0000E-03 1.2605E+01 Heat done<br>
Reorthogonalize Basis:<br>
4 2.21548E-02 2.27058E-02 9.75731E-01 1 velx PROJECT<br>
4 Hmholtz velx: 5 7.8903E-08 4.8255E-02 1.0000E-07<br>
Reorthogonalize Basis:<br>
4 2.21455E-02 2.26942E-02 9.75823E-01 1 vely PROJECT<br>
4 Hmholtz vely: 5 7.8987E-08 4.8248E-02 1.0000E-07<br>
Reorthogonalize Basis:<br>
4 2.54684E-01 2.57237E-01 9.90076E-01 1 velz PROJECT<br>
4 Hmholtz velz: 6 1.7695E-08 5.4216E-01 1.0000E-07<br>
4 U-PRES gmres: 13 9.8513E-07 1.0000E-06 3.0328E-03 1.1984E+00 4.7670E+00<br>
4 4.0000E-03 8.4136E+00 Fluid done<br>
4 4 qfilt: 6.2753E-08 1.4966E-07 4.9326E-10<br>
5.7232E-14 5.7288E-14 2.1974E+09</p>
<p><br>
</p>
<p>*****************<br>
<br>
</p>
<p>The strange thing is that this does not happen when I run with only temperature.</p>
<p>I develop my simulation in small steps, but at the step when I also simulate the scalars (that are treated identically as the temperature), I get this crash.<br>
</p>
<p>I wonder if there is somthing that I missed when increasing the number of scalars?</p>
<p><br>
</p>
<p>Below I describe how I set up the simulation.</p>
<p><br>
</p>
<p>Best,</p>
<p><br>
</p>
<p>Johan<br>
</p>
<p><br>
</p>
<p>*****************************</p>
<p><br>
</p>
<p>*** Only temperature case ***</p>
<p><br>
</p>
<p></p>
<p>I generate a mesh from six .box files using genbox. The corresponding .rea files were modified slightly, merged and extended to 3D using n2to3.</p>
<p>The 3D .rea files were merged using nekmerge2.</p>
<br>
<p></p>
<p><br>
</p>
<p>If I start from .box files with only temperature I get a .rea file with p023 = 0 (NPSCAL).</p>
<p>I set the Peclet number to -233, and the field settings to</p>
<p><br>
T T F F F F F F F F F IFNAV & IFADVC (convection in P.S. fields)<br>
F F T T T T T T T T T T IFTMSH (IF mesh for this field is T mesh)</p>
<p><br>
</p>
<p>end everything works excellently.</p>
<p><br>
</p>
<p></p>
<p>*** 3 scalar case ***<br>
</p>
<p><br>
</p>
<p>If modify the previously mentioned .box files by adding 2 to the number of fields, and adding the same BC:s to them.</p>
<p>I get a .rea file with p023 = 2 (NPSCAL).<br>
</p>
<p><br>
</p>
<p>I set the Peclet number for the new scalars to -233, so that they are the same as for the temperature.</p>
<p>I change the field settings to:</p>
<p><br>
T T T T F F F F F F F IFNAV & IFADVC (convection in P.S. fields)<br>
F F F F T T T T T T T T IFTMSH (IF mesh for this field is T mesh)</p>
<p><br>
</p>
<p>Now I get the crash described above.</p>
<p><br>
</p>
<p>*****************************<br>
</p>
<p><br>
</p>
<p>*** The .usr file***<br>
</p>
<p><br>
</p>
<p>The temperature is read out in userf as:<br>
</p>
<p><br>
</p>
<p> ffz = -t(ix,iy,iz,e,1)/1.32</p>
<p><br>
</p>
<p>All the scalars have the same boundary and initial conditions (they enter useric and userf only as "temp").</p>
<p><br>
</p>
<p>*** The SIZE file ***</p>
<p><br>
</p>
<p>I use the maximum number of passive scalars (ldimt) = 4 for both cases.<br>
<br>
</p>
<p></p>
<p><br>
</p>
<p><br>
</p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>