[FLASH-USERS] Negative 3T Internal Energies

[Reyes, Adam]

Hi Jonas,

The most common problem is having tables with negative values of energy or pressure, which some sesame tables have. Opacplot makes no attempts to sanitize the table data, so if they are in there it will propagate into the generated ionmix files. You can try the native SESAME reader in FLASH which will try to remove any negative entries it finds, or initialize your simulation at a higher temperature to try and avoid those negative values that are usually around solid density/temperature.
*********************************************
Adam Reyes


Code Group Leader, Flash Center for Computational Science  
Research Scientist, Dept. of Physics and Astronomy
University of Rochester
River Campus: Bausch and Lomb Hall, 369  
500 Wilson Blvd. PO Box 270171, Rochester, NY 14627
Email adam.reyes at rochester.edu
Web https://flash.rochester.edu
 (he / him / his)


*********************************************



> On May 6, 2024, at 4:52 PM, Jonas Kuhlke <jonas.kuhlke at uni-rostock.de> wrote:
> 
> Hi Flash Users,
> 
> 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 opacplot2(https://github.com/flash-center/opacplot2) 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'.
> 
> 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).
> 
> The log file says:
> [ 05-06-2024  16:13:53.033 ] WARNING after gc filling: min. unk(ENER_VAR)=-371477695.1070602 PE=0 block=2 type=2
>  [ 05-06-2024  16:13:53.033 ] WARNING after gc filling: min. unk(EINT_VAR)=-371477695.1070602 PE=0 block=2 type=2
>  [ 05-06-2024  16:13:53.033 ] WARNING after gc filling: min. unk(ENER_VAR)=-371477695.1070602 PE=0 block=4 type=2
>  [ 05-06-2024  16:13:53.034 ] WARNING after gc filling: min. unk(EINT_VAR)=-371477695.1070602 PE=0 block=4 type=2
> 
> 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?
> 
> Also this is my setup command:
> ./setup -auto LaserSlab -1d -nxb=16 +hdf5typeio \ species=cham,targ +mtmmmt +laser +uhd3t \ -parfile=PE_1D.par
> 
> Kind regards,
> Jonas
> 
> 
> <matr_07180_polyethylene.cn4><matr_07180_polyethylene.ses><PE_lasslab.log><PE_1D.par><matr_07590_polystyrene.cn4>_______________________________________________
> flash-users mailing list
> flash-users at flash.rochester.edu <mailto:flash-users at flash.rochester.edu>
> 
> For list info, including unsubscribe:
> https://flash.rochester.edu/mailman/listinfo/flash-users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20240506/dd47f2b5/attachment-0002.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: FLASH.jpg
Type: image/jpeg
Size: 23876 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20240506/dd47f2b5/attachment-0002.jpg>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: FLASH-pride-sml.png
Type: image/png
Size: 12732 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20240506/dd47f2b5/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: smime.p7s
Type: application/pkcs7-signature
Size: 1391 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20240506/dd47f2b5/attachment-0002.p7s>