<div dir="ltr">Dear Klaus Weide,<div><br></div><div>Thank you for your response</div><div><br></div><div>I have <span style="color:rgb(33,33,33)">already tried</span><span style="color:rgb(33,33,33)"> </span><span style="color:rgb(33,33,33)">CFL=1.0 and u</span><span style="color:rgb(33,33,33)">sing eintSwich.</span></div><div><span style="color:rgb(33,33,33)">But I still got the error message "Negative 3T </span><span style="color:rgb(33,33,33)">internal energy".</span><br></div><div><span style="color:rgb(33,33,33)"><br></span></div><div>Sincerely yours,</div></div><br><div class="gmail_quote"><div dir="ltr">On Wed, Feb 28, 2018 at 5:55 AM Klaus Weide <<a href="mailto:klaus@flash.uchicago.edu">klaus@flash.uchicago.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">On Wed, 28 Feb 2018, Kazuki Matsuo wrote:<br>
><br>
> I'm working on 3T MHD simulation of laser-generated magnetized plasma.<br>
><br>
> My code works well with a 5T magnetic field in 2D cylindrical<br>
> geometry. However,<br>
> if I increase the magnetic field to 100 T, I got the error message<br>
> "Negative 3T internal energy".<br>
><br>
> Do you know how to deal with such a strong magnetic field?<br>
<br>
Have you already tried the following?<br>
<br>
- Enforce a smaller timestep, for example by setting CFL=1.0 or by<br>
setting dtmax.<br>
<br>
- Use eintSwich, perhaps with a value of 1.e-3 or 1.e-4 (requires<br>
use_auxEintEqn=.TRUE.)<br>
<br>
<br>
These are rather generic suggestions, hopefully others can give more<br>
specific advice relevant to strong magnetic fields.<br>
<br>
Klaus<br>
<br>
</blockquote></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><p dir="ltr">------</p>
<p dir="ltr">Kazuki Matsuo</p><p dir="ltr" style="font-family:"helvetica neue",helvetica,arial,sans-serif">Institute of Laser Engineering Osaka University</p><p dir="ltr" style="font-family:"helvetica neue",helvetica,arial,sans-serif">2-6 Yamada-oka, Suita, Osaka, 565-0871, Japan</p><p dir="ltr" style="font-family:"helvetica neue",helvetica,arial,sans-serif"><span style="color:rgb(0,0,0);font-family:helvetica;font-size:12px;line-height:normal">Webpage:   </span><a href="http://lf-lab.net/" target="_blank" style="font-family:helvetica;font-size:12px;line-height:normal">http://lf-lab.net</a><br></p>
<p dir="ltr">Facebook:  <a href="https://www.facebook.com/kazuki.matsuo.102"> https://www.facebook.com/kazuki.matsuo.102 </a></p><p dir="ltr">E-mail:   <a href="mailto:matsuo-k@ile.osaka-u.ac.jp?src=inbox">matsuo-k@ile.osaka-u.ac.jp</a></p>
</div></div>