<div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Hi all,</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">I've been doing MHD simulations involving shocks in 2D (cylindrical symmetry) and have been occasionally running into cases where the runs aborts with:</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><div class="gmail_default"> dtCheck=  0.000000000000000E+000</div><div class="gmail_default"> DRIVER_ABORT: [Hydro]: Computed dt is not positive! Aborting!</div><div class="gmail_default"><br></div><div class="gmail_default">I believe that this is occurring for cells where the plasma beta, i.e. Pthermal/Pmagnetic, is very small.  These cells are not near any boundaries, so that doesn't seem to be the problem.  I do see that the dt_hydro is dropping quickly just before, and drops below the dtmin that I set.  So it's likely that the proximate cause is the violation of the courant condition in the sense that the timestep is too long given the signal speed.  What I really need to know then is how to prevent the very low density with high magnetic field condition from occurring which leads to the very high alfven speed that leads to this problem.</div><div class="gmail_default"><br></div><div class="gmail_default">Note also that the simulations that I'm doing include cooling, which decreases the thermal pressure.  Some parameters that I've set, partially based on information in the Flash User's Guide are:</div><div class="gmail_default"><div class="gmail_default">killdivb = .true.</div><div class="gmail_default">flux_correct = .true.</div><div class="gmail_default">E_modification = .true.</div><div class="gmail_default">energyFix = .true.</div><div class="gmail_default">shockLowerCFL = .true.</div><div class="gmail_default">RiemannSolver = "HYBRID"</div><div class="gmail_default">prolMethod = "BALSARA_PROL"</div><div class="gmail_default"><br></div><div class="gmail_default">I'm wondering if anyone else has run into problems like this with MHD runs and if anyone has suggestions for getting around this.  There are many possible parameter settings and it could be that I just need to find the best ones for my problem.</div><div class="gmail_default"><br></div><div class="gmail_default">Regards,</div><div class="gmail_default">Jon</div></div></div><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr">________________________________________________________<br>Jonathan D. Slavin                 Harvard-Smithsonian CfA<br><a href="mailto:jslavin@cfa.harvard.edu" target="_blank">jslavin@cfa.harvard.edu</a>       60 Garden Street, MS 83<br>phone: (617) 496-7981       Cambridge, MA 02138-1516<br>cell: (781) 363-0035             USA<br>________________________________________________________<br><br></div></div></div></div>
</div>