<html><head><meta http-equiv="Content-Type" content="text/html charset=windows-1252"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">Hello Everyone,<div><br></div><div>I have been trying to run an MHD simulation of a magnetized, laser-generated plasma. The original, non magnetized simulations, were based off of the LaserSlab example and I was able to run that fine. I then created an MHD simulation based on the FieldLoop example, but incorporating all of the necessary modules for the 3T and laser packages. I have a running simulation and the initial outputs look fine, but the simulation only runs 10-15 time steps before crashing.</div><div><br></div><div>When I run in the full MHD mode, whether I use eos_tab or eos_gam, the error that eventually crashes it is:</div><div><br></div><div> DRIVER_ABORT: [hy_uhd_ragelike] Negative 3T internal energy, CHECK LOG</div><div><br></div><div>however, before every step I get the overflow error from multiple processors:</div><div>…</div><div><div> Overflow! Mesh PE          24 , gettingFrom           1 , recvCount          60 , max particles       20000</div><div> Overflow! Mesh PE          25 , gettingFrom           1 , recvCount          12 , max particles       20000</div><div> Overflow! Mesh PE          24 , gettingFrom           2 , recvCount          69 , max particles       20000</div></div><div>…</div><div>but this does not seem to crash the run, it just seems odd (computationally unhealthy?)</div><div><br></div><div>I have tried changing the resolution of the simulation and set my cfl=0.2, but the negative internal energy error still comes up. I increased ed_maxRayCount to 40000 to see if I could get rid of the Overflow! error and all that changes is “… max particle   40000”. </div><div><br></div><div>I have tried a different polystyrene table (both tables from the original LaserSlab example), and it still crashes with a ’Negative 3T internal energy' error. </div><div><br></div><div>I changed eosModeInit from “dens_temp_gather” to “dens_ie_gather” and now the internal energy error is gone, but a different EOS error pops up ending in:</div><div><br></div><div>DRIVER_ABORT: [Eos] Error: too many iteration in Newton-Raphson</div><div><br></div><div>I have increased eos_maxNewton=2000 and set eos_tolerance=1.0e-07, but this doesn’t seem to fix anything.</div><div><br></div><div>I tried the ForcedHydroLimit=.true. scenario and I still end up with the same error behavior.</div><div><br></div><div><br></div><div>—> With the same geometric setup and EOS tables, the previous simulation using the +uhd3t solver seems to run okay for the most part. However, I have not gotten any version of the simulation to run with the +usm3t solver even when I use the ForcedHydroLimit. </div><div><br></div><div><br></div><div>Any and all suggestions are much appreciated. :)</div><div><br></div><div>Cheers,</div><div>Mario Manuel</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br><div apple-content-edited="true">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;"><div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">______________________________________</div><div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">                                                                                                                               <span class="Apple-converted-space"> </span></div><div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">Dr. Mario Manuel<br>Einstein Postdoctoral Fellow, Smithsonian Astrophysical Observatory<br><a href="mailto:mmanuel@umich.edu">mmanuel@umich.edu</a><br>617-899-8381<br><br>University of Michigan<br>Space Research Building<br>2455 Hayward Street<br>Ann Arbor, MI 48109<br>Office: 734-764-4585<br></div><div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">______________________________________</div></div>
</div>
<br></div></body></html>