[FLASH-USERS] Negative 3T Internal Energies

[赵旭]

Hi Adam,

I want to follow this issue a bit.

My FLASH version is 4.8, and I use the native SESAME reader in FLASH.

(not sure if it is means setting eos_targSubType   = "sesame" in .par)

I also use PE (7173) as target in laserSlab, but same I get similar negative pressure error.

I have tried set initial Te/Tion/Trad = 2000K, but after a while the error still.

Another thing is when setting higher temperature, the density / Te map become wired (see attached figure).

I dont know how to explain/ fix these.


Best,
Xu  

----- 原始邮件 -----
发件人: "Adam Reyes" <adam.reyes at rochester.edu>
收件人: "Jonas Kuhlke" <jonas.kuhlke at uni-rostock.de>
抄送: "flash-users" <flash-users at flash.rochester.edu>
发送时间: 星期一, 2024年 5 月 06日 下午 4:36:02
主题: Re: [FLASH-USERS] Negative 3T Internal Energies

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 | 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 
[ mailto:flash-users at flash.rochester.edu | flash-users at flash.rochester.edu ] 

For list info, including unsubscribe: 
[ https://flash.rochester.edu/mailman/listinfo/flash-users | https://flash.rochester.edu/mailman/listinfo/flash-users ] 


_______________________________________________
flash-users mailing list
flash-users at flash.rochester.edu

For list info, including unsubscribe:
https://flash.rochester.edu/mailman/listinfo/flash-users
-- 
Zhao Xu
Laboratory for Laser Plasmas (MoE)
Shanghai Jiao Tong University
800 Dongchuan Rd, Shanghai 200240
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Te.png
Type: image/png
Size: 243881 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20240924/8fc031d7/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: dens.png
Type: image/png
Size: 404016 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20240924/8fc031d7/attachment-0003.png>