[FLASH-USERS] non convergence in Riemann subroutine

M.A. Latife latife at astro.rug.nl
Mon Nov 3 11:34:31 EST 2008


Hi
Dear All,
I am running 3-D simulations including hydro +  particles + cosmology + 
eos + gravity + cooling. My simulations crash after some time (1.0E15) 
and give the following error. it seems me that there is something going 
wrong in Riemann.F90 with pressure. What are possible causes of this 
problem? How can i avoid this problem?
Any idea
Kind Regards
Latif
Nonconvergence in subroutine rieman!
 
 Processor, block =            7          20
 Zone index       =            5
 Transverse index =            6           5
 Zone center      =   1.170820312500000E+025
 Sweep direction  = x
 Iterations tried =           27
 Pressure error   =    5.45416055415388    
 
 pL       =   4.117027057197573E-022  pR       =  3.938433206566095E-022
 uL       =   -12152456.2730107       uR       =  -10892043.1131041    
 cL       =   1.041411422316767E-026  cR       =  1.028454571441834E-026
 gamma_eL =    1.00765479619375       gamma_eR =   1.00766156258855    
 gamma_cL =    1.00765479619375       gamma_cR =   1.00766156258855    
 
 Iteration history:
 
   n                    p*
   1    0.100000000000E-29
   2    0.178214204180E-25
   3    0.191075724471E-27
   4    0.377139479396E-27
   5    0.178214203785E-25
   6    0.559316960640E-27
   7    0.737728667286E-27
   8    0.178214203773E-25
   9    0.912490167868E-27
  10    0.108371234987E-26
  11    0.178214204060E-25
  12    0.125150165427E-26
  13    0.141596029646E-26
  14    0.178214203879E-25
  15    0.157718647365E-26
  16    0.173527456037E-26
  17    0.178214204639E-25
  18    0.189031529261E-26
  19    0.204239594204E-26
  20    0.178214204634E-25
  21    0.219160047921E-26
  22    0.233800972837E-26
  23    0.178214204400E-25
  24    0.248170151347E-26
  25    0.262275079601E-26
  26    0.178214204915E-25
  27    0.276122980549E-26
 
 Terminating execution.




More information about the flash-users mailing list