[FLASH-USERS] non convergence in Riemann subroutine

[Tomasz Plewa]

Latif -

The gammas close to 1 are a possible cause of this failure.

What equation of state are you using?

Tomek
--
M.A. Latife wrote:
> 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.
>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: tplewa.vcf
Type: text/x-vcard
Size: 296 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20081103/6b6fa0cc/attachment.vcf>