<div dir="ltr">Dear Flash-users,<div><br></div><div>I have been running into a few errors while trying to run a basic StirTurb problem. What I essentially did was that, from the original StirTurb problem package, I only made a few noticeable changes to flash.par and added "REQUIRES physics/Hydro/HydroMain/<wbr>unsplit/Hydro_Unsplit" to the Config file. Then I setup the problem using the command</div><div><br></div><div>./setup —auto -3d StirTurb -nxb=64 -nyb=64 -nzb=64 +ug<br></div><div><div><br></div><div>and it successfully compiled.</div><div><br></div><div>Running it however it proving to be a problem for me. When I set gamma to 1.4, I get a successful run but when I set it to 1.000 (Hoping to achieve isothermal conditions), I get an error that says:</div><div><br></div><div>
<p class="gmail-m_9098461374110142968gmail-p1"><b><span class="gmail-m_9098461374110142968gmail-s1">DRIVER_ABORT: </span>[hy_uhd_<wbr>eigenParameters-A]: Zero or imaginary sound speed has obtained! Please try other (more diffusive) slope limiter, flux, order, cfl, etc.</b></p><p class="gmail-m_9098461374110142968gmail-p1"><br></p><p class="gmail-m_9098461374110142968gmail-p1">Could you please help me decipher where I am going wrong?</p><p class="gmail-m_9098461374110142968gmail-p1">Thank in advance.</p><p class="gmail-m_9098461374110142968gmail-p1">Kezman</p><p class="gmail-m_9098461374110142968gmail-p1"><b><br></b></p><p class="gmail-m_9098461374110142968gmail-p1"><b>Provided below is the flash.par information I am using;</b></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff"># runtime parameters, nxb = nyb = nzb = 64</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">256^3 grid, gamma = 1.000</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">xmax = 1.</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">xmin = 0.</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">ymax = 1.</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">ymin = 0.</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">zmax = 1.</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">zmin = 0.</font></span></p><p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><span class="gmail-s1"></span><br></font></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">basenm = "driventurb_3d_"</font></span></p><p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)">restart = .false.<br><span class="gmail-s1"></span></font></p><p class="gmail-p1"><font color="#0000ff" style="background-color:rgb(255,255,255)"><br></font></p><p class="gmail-p1"><font color="#0000ff" style="background-color:rgb(255,255,255)"># file numbers - if you restart you have to change the checkpointFileNumber</font></p><p class="gmail-p2"><font color="#0000ff"><span class="gmail-s1" style="background-color:rgb(255,255,255)"></span></font></p><p class="gmail-p2"><font color="#0000ff"><span class="gmail-s1" style="background-color:rgb(255,255,255)"></span></font></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">checkpointFileNumber = 0</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">plotFileNumber = 0</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">particleFileNumber = 0</font></span></p><p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><span class="gmail-s1"></span><br></font></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff"># set the time between dumps</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">checkpointFileIntervalTime = 0.01</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">plotFileIntervalTime = 0.</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">particleFileIntervalTime = 0.25</font></span></p><p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><br></font></p><p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"># set the number of steps between dumps<br><span class="gmail-s1"></span></font></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">checkpointFileIntervalStep = 0</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">plotFileIntervalStep = 0</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">particleFileIntervalStep = 0</font></span></p><p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><br></font></p><p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)">useParticles= .false.<br><span class="gmail-s1"></span></font></p><p class="gmail-p2"><font color="#0000ff"><span class="gmail-s1" style="background-color:rgb(255,255,255)"></span></font></p><p class="gmail-p2"><font color="#0000ff"><span class="gmail-s1" style="background-color:rgb(255,255,255)"></span></font></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">pt_numX = 5</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">pt_numY = 5</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">pt_numZ = 5</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">pt_maxPerProc=1000</font></span></p><p class="gmail-m_9098461374110142968gmail-p1">
</p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">plot_var_1 = "dens"</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">plot_var_2 = "pres"</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">plot_var_3 = "temp"</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">plot_var_4 = "velx"</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">plot_var_5 = "vely"</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">plot_var_6 = "velz"</font></span></p><p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><span class="gmail-s1"></span><br></font></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">plot_grid_var_1 = "mvrt"</font></span></p><p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><span class="gmail-s1"></span><br></font></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">gamma = 1.000</font></span></p><p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><span class="gmail-s1"></span><br></font></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">cfl = 0.8</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">nend = 10000</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">tmax = 0.05 # was originally 20.0</font></span></p><p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><span class="gmail-s1"></span><br></font></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">xl_boundary_type = "periodic"</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">xr_boundary_type = "periodic"</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">yl_boundary_type = "periodic"</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">yr_boundary_type = "periodic"</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">zl_boundary_type = "periodic"</font></span></p><p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><span class="gmail-s1"></span>
</font></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">zr_boundary_type = "periodic"</font></span></p></div>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">st_stirmax = 18.8478</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">st_stirmin = 6.2832</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">st_energy = 0.1</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">st_decay = 0.5</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">st_freq = 1</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">eintSwitch = 1.</font></span></p>
<p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><span class="gmail-s1"></span><br></font></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">dtinit = 1.e-5</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">dtmin = 1.e-12 # This parameter must be << minimum timestep</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff"> # in order to avoid numerical instability</font></span></p>
<p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><span class="gmail-s1"></span><br></font></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">smallt = 1.e-15</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">smalle = 1.e-10</font></span></p>
<p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><span class="gmail-s1"></span><br></font></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff"># AMR refinement parameters</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">#lrefine_max = 6</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">#refine_var_1 = "dens"</font></span></p>
<p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><span class="gmail-s1"></span><br></font></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">#These parameters below are only necessary for the Uniform Grid</font></span></p>
<p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><span class="gmail-s1"></span><br></font></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">iProcs = 4 #num procs in i direction</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">jProcs = 4 #num procs in j direction</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">kProcs = 4</font></span></p>
<p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><span class="gmail-s1"></span><br></font></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff"># When using UG, iProcs, jProcs and kProcs must be specified.</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff"># These are the processors along each of the dimensions</font></span></p>
<p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><span class="gmail-s1"></span><br></font></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">#FIXEDBLOCKSIZE mode ::</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff"># When using fixed blocksize, iGridSize etc are redundant in</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff"># runtime parameters. These quantities are calculated as</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff"># iGridSize = NXB*iprocs</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff"># jGridSize = NYB*jprocs</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff"># kGridSize = NZB*kprocs</font></span></p>
<p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><span class="gmail-s1"></span><br></font></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">#NONFIXEDBLOCKSIZE mode ::</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff"># iGridSize etc must be specified. They constitute the global</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff"># number of grid points in the physical domain without taking</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff"># the guard cell into account. The local blocksize is calculated</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff"># as iGridSize/iprocs etc.</font></span></p><p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><br><span class="gmail-s1"></span></font></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">## -------------------------------------------------------------##</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">## SWITCHES SPECIFIC TO THE UNSPLIT HYDRO SOLVER ##</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff"># I. INTERPOLATION SCHEME:</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">order = 3 # Interpolation order (first/second/third/fifth order)</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">slopeLimiter = "vanLeer" # Slope limiters (minmod, mc, vanLeer, hybrid, limited)</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">LimitedSlopeBeta= 1. # Slope parameter for the "limited" slope by Toro</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">charLimiting = .true. # Characteristic limiting vs. Primitive limiting</font></span></p><p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><span class="gmail-s1"></span><br></font></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">use_avisc = .false. # use artificial viscosity (originally for PPM)</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">#cvisc = 0.1 # coefficient for artificial viscosity</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">use_flattening = .false. # use flattening (dissipative) (originally for PPM)</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">use_steepening = .false. # use contact steepening (originally for PPM)</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">use_upwindTVD = .false. # use upwind biased TVD slope for PPM (need nguard=6)</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">use_hybridOrder = .true.</font></span></p><p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><span class="gmail-s1"></span><br></font></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">
</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff"># Magnetic and Electric Fields:</font></span></p>
<div>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">E_modificaton = .true.</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">energyFix = .true.</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">ForceHydroLimit = .false.</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">prolMethod = "injection_prol"</font></span></p>
<p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><span class="gmail-s1"></span><br></font></p>
<p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><span class="gmail-s1"></span><br></font></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff"># II. RIEMANN SOLVERS:</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">RiemannSolver = "Hybrid" # Roe, HLL, HLLC, LLF, Marquina</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">entropy = .true. # Entropy fix for the Roe solver</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">EOSforRiemann = .true. # Call EOS in Riemann flux calculations</font></span></p>
<p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)"><span class="gmail-s1"></span><br></font></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff"># III. STRONG SHOCK HANDELING SCHEME:</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">shockDetect = .true. # Shock Detect for numerical stability</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">## -------------------------------------------------------------##</font></span></p>
<p class="gmail-p2"><font color="#0000ff" style="background-color:rgb(255,255,255)">## ---------------------------------------------------------------##<br><span class="gmail-s1"></span></font></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">## SWITCHES SPECIFIC TO THE SUPER-TIME-STEPPING (STS) ALGORITHM ##</font></span></p>
<p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">## NOTE: For details on using STS runtime parameters, please ##</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">## refer to user's guide (Driver chapter). ##</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">useSTS = .false.</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">nstepTotalSTS = 5</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">nuSTS = 0.2</font></span></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff">## ---------------------------------------------------------------##</font></span></p><p class="gmail-p1"><br></p><p class="gmail-p1"><br></p><p class="gmail-p1"><span class="gmail-s1" style="background-color:rgb(255,255,255)"><font color="#0000ff"><br></font></span></p><p class="gmail-p2"><span class="gmail-s1"></span><br></p><p class="gmail-p1"><span class="gmail-s1">
</span></p><p class="gmail-p2"><span class="gmail-s1"></span><br></p></div><div><br></div><div><br></div></div></div>