[FLASH-USERS] 3T MHD with multi species EOS tables - negative internal energy

[長島圭祐]

Hello Mario,

Sorry for asking about the post of 2 years ago but I'm working on 3T MHD 
simulation of laser generated , magnetized plasma with multi species 
just like your case and I got the same error message "Negative 3T 
internal energy" when I run my code.

Following is the error message written in the log file. Do you still 
remember how you fixed this error? Is it the problem of runtime 
parameters or bad writing in Simulation initBlock file?

   1  [ 11-13-2017  21:18:05.266 ] [hy_uhd_ragelike] Negative internal 
energies in 3T update
   2  [ 11-13-2017  21:18:05.266 ]   Falling back to using internal 
energy ratios instead of pressure ratios
   3  [ 11-13-2017  21:18:05.266 ]
   4  [ 11-13-2017  21:18:15.866 ] [hy_uhd_unsplitUpdateMultiTemp] ERROR 
DETECTED (adv_err_chk)
   5  [ 11-13-2017  21:18:15.866 ]   Zero or negative specific internal 
energy detected
   6  [ 11-13-2017  21:18:15.866 ]   DENS_OLD = 1.3388944513161E-04
   7  [ 11-13-2017  21:18:15.866 ]   DENS_NEW = 1.0000000000000E-10
   8  [ 11-13-2017  21:18:15.866 ]   DUTOT_ADV = -1.8113834390984E+10
   9  [ 11-13-2017  21:18:15.866 ]   DUTOT_NEW = -5.7850199808167E+09
  10  [ 11-13-2017  21:18:15.866 ]   UTOT_OLD = 5.7850199808167E+09
  11  [ 11-13-2017  21:18:15.866 ]   UTOT_NEW = 1.0000000000000E-20
  12  [ 11-13-2017  21:18:15.866 ]     UELE_OLD = 5.7839162928446E+09
  13  [ 11-13-2017  21:18:15.867 ]     UION_OLD = 1.1026662495695E+06
  14  [ 11-13-2017  21:18:15.867 ]     URAD_OLD = 1.0217225176679E+03
  15  [ 11-13-2017  21:18:15.867 ]     UELE_NEW = 3.0683578349976E+08
  16  [ 11-13-2017  21:18:15.867 ]     UION_NEW = 5.2776254928130E+04
  17  [ 11-13-2017  21:18:15.867 ]     URAD_NEW = 5.0566403937702E+02
  18  [ 11-13-2017  21:18:15.867 ]     uele_adv  = -1.2326462540022E+10
  19  [ 11-13-2017  21:18:15.867 ]     uion_adv  = -2.3499163683085E+06
  20  [ 11-13-2017  21:18:15.867 ]     urad_adv  = -1.9537770163699E+03
  21  [ 11-13-2017  21:18:15.867 ]     duele_adv = -1.8110378832867E+10
  22  [ 11-13-2017  21:18:15.867 ]     duion_adv = -3.4525826178780E+06
  23  [ 11-13-2017  21:18:15.867 ]     durad_adv = -2.9754995340377E+03
  24  [ 11-13-2017  21:18:15.867 ]   TELE_VAR = 3.1148090218563E+05
  25  [ 11-13-2017  21:18:15.867 ]   TION_VAR = 1.4810060769597E+03
  26  [ 11-13-2017  21:18:15.867 ]   TRAD_VAR = 1.9169934413387E+04
  27  [ 11-13-2017  21:18:15.867 ]   CELL X =  4.1796875000000E-03
  28  [ 11-13-2017  21:18:15.867 ]   CELL Y =  1.1976562500000E-02
  29  [ 11-13-2017  21:18:15.867 ]   CELL Z =  0.0000000000000E+00
  30  [ 11-13-2017  21:18:15.867 ]   You might be able to get things to 
keep running by trying
  31  [ 11-13-2017  21:18:15.867 ]   more diffusive hydro options 
(reduce CFL, order, etc...)

I appreciate any kinds of advice and suggestion.

Thanks,

-- 
長島圭祐
Keisuke Nagashima
Advanced Space Propulsion Laboratory
Department of Energy Science and Engineering
School of Engineering
Kyushu University