[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