<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>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>
</body>
</html>