<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:DengXian;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:"\@DengXian";
panose-1:2 1 6 0 3 1 1 1 1 1;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.msonormal0, li.msonormal0, div.msonormal0
{mso-style-name:msonormal;
mso-margin-top-alt:auto;
margin-right:0in;
mso-margin-bottom-alt:auto;
margin-left:0in;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
span.EmailStyle18
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">Hi, Michel:<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">Thank you!<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">I used the tabulated EOS provided by FLASH center in my simulation. If the problem comes from EOS, should I change the EOS?<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">Best regards,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">Wenyi<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> flash-users-bounces@flash.uchicago.edu [mailto:flash-users-bounces@flash.uchicago.edu]
<b>On Behalf Of </b>Michel de Peretti<br>
<b>Sent:</b> Wednesday, April 17, 2019 8:00 PM<br>
<b>To:</b> flash-users@flash.uchicago.edu<br>
<b>Subject:</b> Re: [FLASH-USERS] flash-users Digest, Vol 138, Issue 7<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">I am not fully familiar with Flash but very often this kind of problem comes from your equation of state....try to check it.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Hoping that helps you<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Michel<o:p></o:p></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">Le mer. 17 avr. 2019 à 19:00, <<a href="mailto:flash-users-request@flash.uchicago.edu">flash-users-request@flash.uchicago.edu</a>> a écrit :<o:p></o:p></p>
</div>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-right:0in">
<p class="MsoNormal">Send flash-users mailing list submissions to<br>
<a href="mailto:flash-users@flash.uchicago.edu" target="_blank">flash-users@flash.uchicago.edu</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://flash.uchicago.edu/mailman/listinfo/flash-users" target="_blank">
http://flash.uchicago.edu/mailman/listinfo/flash-users</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:flash-users-request@flash.uchicago.edu" target="_blank">flash-users-request@flash.uchicago.edu</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:flash-users-owner@flash.uchicago.edu" target="_blank">flash-users-owner@flash.uchicago.edu</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of flash-users digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. The problem of Negative Internal Energy in 3T MHD simulation<br>
with laser generated plasmas (Huo Wenyi)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Wed, 17 Apr 2019 07:43:27 +0000<br>
From: Huo Wenyi <<a href="mailto:Wenyi.Huo@eli-beams.eu" target="_blank">Wenyi.Huo@eli-beams.eu</a>><br>
Subject: [FLASH-USERS] The problem of Negative Internal Energy in 3T<br>
MHD simulation with laser generated plasmas<br>
To: "<a href="mailto:flash-users@flash.uchicago.edu" target="_blank">flash-users@flash.uchicago.edu</a>" <<a href="mailto:flash-users@flash.uchicago.edu" target="_blank">flash-users@flash.uchicago.edu</a>><br>
Message-ID: <<a href="mailto:4b8b49f0356446ab9b4326bb39c14550@eli-beams.eu" target="_blank">4b8b49f0356446ab9b4326bb39c14550@eli-beams.eu</a>><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
Dear all:<br>
<br>
I am using FLASH to simulate the magnetohydrodynamics in laser generated plasmas. However, I always encounter a problem of negative internal energy both in FLASH 4.5 and FLASH4.6. Does anyone encounter such a similar problem? The error information in the
Log file is as follows:<br>
<br>
<br>
[ 04-16-2019 14:34:35.488 ] [hy_uhd_ragelike] ERROR DETECTED<br>
[ 04-16-2019 14:34:35.488 ] Zero or negative advected specific internal energy detected<br>
[ 04-16-2019 14:34:35.488 ] DENS_OLD = 1.0000904290291E-09<br>
[ 04-16-2019 14:34:35.488 ] DENS_NEW = 1.0000983091349E-09<br>
[ 04-16-2019 14:34:35.488 ] UTOT_NEW = -3.7386802822434E-01<br>
[ 04-16-2019 14:34:35.488 ] UEXTRA = -1.7491864582091E+01<br>
[ 04-16-2019 14:34:35.488 ] UELE_NEW = -3.7135596091566E-01<br>
[ 04-16-2019 14:34:35.488 ] UION_NEW = -2.5120492887951E-03<br>
[ 04-16-2019 14:34:35.488 ] URAD_NEW = -1.8019885718467E-08<br>
[ 04-16-2019 14:34:35.488 ] EELE_VAR = 1.7001438966846E+10<br>
[ 04-16-2019 14:34:35.488 ] EION_VAR = 1.1498940630502E+08<br>
[ 04-16-2019 14:34:35.488 ] ERAD_VAR = 8.2498749422145E+02<br>
[ 04-16-2019 14:34:35.488 ] EELE = -3.7131945682109E+08<br>
[ 04-16-2019 14:34:35.488 ] EION = -2.5118023556785E+06<br>
[ 04-16-2019 14:34:35.488 ] ERAD = -1.8018114373230E+01<br>
[ 04-16-2019 14:34:35.488 ] uele_adv = 1.7002978536045E+01<br>
[ 04-16-2019 14:34:35.488 ] uion_adv = 1.1501719275908E-01<br>
[ 04-16-2019 14:34:35.488 ] urad_adv = 8.2506210304961E-07<br>
[ 04-16-2019 14:34:35.488 ] PELE_VAR = 1.2145270532940E+00<br>
[ 04-16-2019 14:34:35.488 ] PION_VAR = 7.6660475248193E-02<br>
[ 04-16-2019 14:34:35.488 ] PRAD_VAR = 2.7502069901319E-07<br>
[ 04-16-2019 14:34:35.488 ] PELE_ADV = 1.2145343502140E+00<br>
[ 04-16-2019 14:34:35.488 ] PION_ADV = 7.6672066071296E-02<br>
[ 04-16-2019 14:34:35.488 ] PRAD_ADV = 2.7502070101654E-07<br>
[ 04-16-2019 14:34:35.488 ] TELE_VAR = 6.0996867323016E+02<br>
[ 04-16-2019 14:34:35.488 ] TION_VAR = 3.6907447766733E+00<br>
[ 04-16-2019 14:34:35.489 ] TRAD_VAR = 1.0219015193890E+02<br>
[ 04-16-2019 14:34:35.489 ] CELL X = 3.3222656250000E-01<br>
[ 04-16-2019 14:34:35.489 ] CELL Y = 6.9890625000000E-01<br>
[ 04-16-2019 14:34:35.489 ] CELL Z = 0.0000000000000E+00<br>
[ 04-16-2019 14:34:35.489 ] You might be able to get things to keep running by trying<br>
[ 04-16-2019 14:34:35.489 ] more diffusive hydro options (reduce CFL, order, etc...)<br>
[ 04-16-2019 22:00:17.496 ] [hy_uhd_ragelike] Negative internal energies in 3T update<br>
[ 04-16-2019 22:00:17.496 ] Falling back to using internal energy ratios instead of pressure ratios<br>
<br>
<br>
<br>
<br>
I have tried many different RiemannSolvers. It does not help. I guess that this problem is related to the grid size, but I am not sure.<br>
<br>
Can anyone provide a solution to this problem? Thank you very much!<br>
<br>
Best regards,<br>
<br>
Wenyi Huo<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://flash.uchicago.edu/pipermail/flash-users/attachments/20190417/a66bdae6/attachment-0001.html" target="_blank">
http://flash.uchicago.edu/pipermail/flash-users/attachments/20190417/a66bdae6/attachment-0001.html</a>
<br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
flash-users mailing list<br>
<a href="mailto:flash-users@flash.uchicago.edu" target="_blank">flash-users@flash.uchicago.edu</a><br>
<a href="http://flash.uchicago.edu/mailman/listinfo/flash-users" target="_blank">http://flash.uchicago.edu/mailman/listinfo/flash-users</a><br>
<br>
<br>
End of flash-users Digest, Vol 138, Issue 7<br>
*******************************************<o:p></o:p></p>
</blockquote>
</div>
</div>
</body>
</html>