[FLASH-USERS] [External Sender] Negative Internal Energies in the 2D, AMR, Z-Pinch problem

[Paul Campbell]

Thanks for the advice, I've resubmitted a job with the extra parameters you've suggested.  So far the log files are still complaining about Negative internal energy but I will let it run over night and see how far it gets.  I should have included in my first email that I'm using version 4.8 (although for your suggestion I don't think the difference in version matters).  For reference this is what the warning from the log files look like:

"
[ 11-04-2024  13:44:58.918 ] [Heatexchange] Warning: Negative internal energy detected
[ 11-04-2024  13:44:58.920 ] [Heatexchange] Correcting for this...
[ 11-04-2024  13:44:58.922 ] [Heatexchange] OLD EION/TION:  7.292804E+10 5.244171E+03
[ 11-04-2024  13:44:58.923 ] [Heatexchange] OLD EELE/TELE:  5.715019E+08 7.379651E+03
[ 11-04-2024  13:44:58.925 ] [Heatexchange] BAD EION/TION:  7.721564E+10 5.546998E+03
[ 11-04-2024  13:44:58.926 ] [Heatexchange] BAD EELE/TELE: -3.716098E+09 5.546998E+03
[ 11-04-2024  13:44:58.927 ] [Heatexchange] CORRECTED EION:  7.721564E+10
[ 11-04-2024  13:44:58.929 ] [Heatexchange] CORRECTED EELE:  2.116656E+08
[ 11-04-2024  13:44:58.930 ] [Heatexchange] FINAL EION:  7.329861E+10
[ 11-04-2024  13:44:58.931 ] [Heatexchange] FINAL EELE:  2.009281E+08
[ 11-04-2024  13:44:58.932 ]
 [ 11-04-2024  13:44:58.933 ] [Heatexchange] Note: Further warnings suppressed
[ 11-04-2024  16:00:05.980 ] Species  1: Lookup for T > Tmax(table): 35 times,first for (3.1052449E+08,1.0000000E+16),highest T was 2.3375415E+10

...... (lots of look ups for the 3 species here) .......


[ 11-04-2024  16:27:26.386 ] [hy_uhd_ragelike] Negative internal energies in 3T update
[ 11-04-2024  16:27:26.387 ]   Falling back to using internal energy ratios instead of pressure ratios
"
Paul Campbell


From: Sutcliffe, Graeme <sutcliffe2 at llnl.gov>
Sent: Tuesday, November 5, 2024 1:56 PM
To: Paul Campbell <Paul.Campbell at zap.energy>
Cc: flash-users at flash.rochester.edu
Subject: Re: [External Sender] [FLASH-USERS] Negative Internal Energies in the 2D, AMR, Z-Pinch problem

Hi Paul,

You can try the following par file options if you haven't already, which helped me to get some laserslab runs working back in FLASH4.7 (adjust accordingly if you're not doing a 3T sim):
dr_usePOSdefComputeDt = .true.
dr_posdefdtfactor = 0.1
dr_numPosdefVars = 3
dr_posdefVar_1 = "eele"
dr_posdefVar_2 = "eion"
dr_posdefVar_3 = "erad"

The dr_posdefdtfactor of 0.1 might be a bit aggressive (i.e. too small), I haven't played with it.

Graeme
________________________________
From: flash-users <flash-users-bounces at flash.rochester.edu<mailto:flash-users-bounces at flash.rochester.edu>> on behalf of Paul Campbell <Paul.Campbell at zap.energy<mailto:Paul.Campbell at zap.energy>>
Sent: Tuesday, November 5, 2024 1:37 PM
To: flash-users at flash.rochester.edu<mailto:flash-users at flash.rochester.edu> <flash-users at flash.rochester.edu<mailto:flash-users at flash.rochester.edu>>
Subject: [FLASH-USERS] Negative Internal Energies in the 2D, AMR, Z-Pinch problem


Hello,



I am attempting to get the 2D Z-pinch problem working with AMR but have been unable to get the sim to run without crashing .  From the logs, I appear to be running into many negative internal energy warnings.  It looks like both the electron and ion temperatures/energies need to be continuously corrected.  To try and address this issue, I have tried refining on density, electron, and ion temperatures but the warnings persist.  In terms of the AMR settings, I am using a 32x32 grid with 4x4 cells, a lrefine_max=4, and I have left the refinement cutoffs and filters unchanged from default.  My questions are:



  1.  Is this an issue with the AMR settings, would a finer grid fix this problem?
  2.  Should I be refining on different parameters?
  3.  Is this an issue with a non-AMR setting?



Paul Campbell






________________________________

This email has been scanned for spam and viruses by Proofpoint Essentials. Click here<https://us5.proofpointessentials.com/app/report_spam.php?mod_id=11&mod_option=logitem&report=1&type=easyspam&k=k1&payload=53616c7465645f5fc5d6dc464f575bf4ef6bc95108f1f0ef2d9880e40aeae0c25b815b2655a594578ace579db3464caf5dada9100109655a34a927175d17dbbd15270519ea6b4378b2174a5d71d94beb73288d62ad5883b3763b544b71ced5ec7beac30299ad987ac4165a4cf13951b55989d3ea9d728e15983b1e3c4ecf5ff0fb2be42c66752dd61d0c29892653871ae11a65ef57a0a3ee8faa74976014d659> to report this email as spam.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20241106/ce0fd307/attachment-0001.htm>