<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p>Hi Flash Users,</p>
<p><br>
</p>
<p>I'm currently trying simulated Laser shock compression of polyethylene targets using LaserSlab. For the EOS tables I have converted SESAME tables into the IONMIX4 format using
<span>opacplot2</span> (<a href="https://github.com/flash-center/opacplot2" class="OWAAutoLink" id="LPlnk916438">https://github.com/flash-center/opacplot2</a>) and this works great for other target materials like polystyrene, but not for PE for some reason.
The error message given is 'DRIVER_ABORT: [hy_uhd_unsplitUpdateMultiTemp] Negative 3T internal energy, CHECK LOG'.
<br>
<br>
I've noticed from mail archive, that I'm not the only one with this problem and I tried the proposed solutions there. Most of the time this seems to be related to some numerical instability, but I'm not sure if this is the case here, because the simulation
crashes in the first time step (I think). Also PS works and PE doesn't with the same parameters (except material specific ones).<br>
<br>
The log file says:</p>
<p></p>
<div>[ 05-06-2024 16:13:53.033 ] WARNING after gc filling: min. unk(ENER_VAR)=-371477695.1070602 PE=0 block=2 type=2<br>
[ 05-06-2024 16:13:53.033 ] WARNING after gc filling: min. unk(EINT_VAR)=-371477695.1070602 PE=0 block=2 type=2<br>
[ 05-06-2024 16:13:53.033 ] WARNING after gc filling: min. unk(ENER_VAR)=-371477695.1070602 PE=0 block=4 type=2<br>
[ 05-06-2024 16:13:53.034 ] WARNING after gc filling: min. unk(EINT_VAR)=-371477695.1070602 PE=0 block=4 type=2</div>
<div><br>
</div>
<div>This makes me think something is wrong with the spatial grid, boundary conditions or my IONMIX EOS for PE. Has someone ran into similar issues or knows a way to verify my EOS tables?</div>
<div><br>
</div>
<div>Also this is my setup command:</div>
<div>
<div>./setup -auto LaserSlab -1d -nxb=16 +hdf5typeio \ species=cham,targ +mtmmmt +laser +uhd3t \ -parfile=PE_1D.par<br>
<br>
Kind regards,<br>
Jonas<br>
</div>
<br>
</div>
<br>
<p></p>
</div>
</body>
</html>