[FLASH-USERS] Wrong sum of fraction in ionization unit

[Suoqing Ji]

Hi,

I’m simulating low density plasma (~1.e-4 cm^-3 in number density) with non-equilibrium ionization. The simulation lasts for ~10 timesteps and then aborts with the following epsilon check error:

---------------------------------------
 Element I.D.:           10  elsum =  0.99976923704261000     
 error in routine neimn: wrong sum of population fraction
 DRIVER_ABORT: error in routine neimn: wrong sum of pop. fraction
  
 Element I.D.:           10  elsum =  0.99976923704261000     
 error in routine neimn: wrong sum of population fraction
 DRIVER_ABORT: error in routine neimn: wrong sum of pop. fraction
---------------------------------------

I keep all the parameters in ionization unit as default ones except dneimin = 1.e-6.

Is there a way to fix or circumvent this issue? Say, will it be appropriate to renormalize the fraction in someway?

Thanks in advance for any suggestions!

Best wishes,
--
Suoqing JI
Ph.D Candidate
Department of Physics
University of California, Santa Barbara
http://web.physics.ucsb.edu/~suoqing

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20150803/6e0681e5/attachment.htm>