user sim does not conserve energy
user alexei.baskin on yahoo reported following on 2013 Jan 13:
Hello,
I have encountered a peculiar pressure convergence failure in
Cloudy v10.00. For several zones, the code “locks” onto a large
pressure oscillation, and fails to converge. An example is
described below, using the attached files:
debug-8-3-8.out – the code output
debug-8-3-8.hist – pressure history
debug-8-3-8.temp – temperature history
debug-8-3-8.line.press – line pressure
For 2 iteration, zone 407 (i.e., the last zone calculated by the
code) the current pressure jumps repeatedly between ~4.6e1 and
6.5e-1 (see debug-8-3-8.hist, line 6628 onward). The current
pressure jumps, while the evaluated density and temperature stay
approximately constant, ~6.7e9 and 3e5, respectively (for
temperature, see debug-8-3-8.temp, line 28619 onward).
I think that the problem is in the evaluated line pressure. Looking
at the line pressure output (debug-8-3-8.line.press, starting from
line 1174), it appears there are two distinct line pressure cases
for adjacent zones. The first case (probably the correct one) has
line pressure which is ~5% from the total pressure, and Cloudy
lists the contributing lines. The second case (probably the
erroneous one) has line pressure that is ~99% from the total
pressure, but Cloudy does not list any contributor.
Note that the occurrence of convergence failure somehow depends on
the particular stopping criterion, which is found beyond the
problematic zone (i.e., deeper into the slab). For example, if the
attached model is ran with “stop temperature 1e5K” as the stopping
criterion, the pressure does converge.
Is it a problem in the code, or am I doing something wrong?
Alexei Baskin
This was with C10. With C13 at r7213 we do not conserve energy:
badprt: all entries with greater than 2.00% of incident continuum follow.
badprt: Intensities are relative to total energy in incident continuum.
Sum of energy in recombination lines is 7.87e+09 relative to total incident radiation is 1.43e-02
all strong information lines
line wl ener/total
all strong cooling lines
line wl ener/total
Sum of coolants (abs) = 1.89e+12 (rel)= 3.44e+00
HFBc 0 0.850 c
Hind 0 0.440 c
HFFc 0 0.325 c
MeFB 0 0.023 c
Clin 912.0A 0.130 c
CION 0 0.388 c
Clin 506.0A 0.029 c
Fe2c 0 0.791 c
H-FB 0 0.291 c
H2+ 0 0.145 c
all strong heating lines
line wl ener/total
Sum of heat (abs) = 2.32e+12 (rel)= 4.22e+00
BFH1 0 1.912 h
BFHe 0 0.195 h
TotM 0 0.307 h
CT H 0 0.039 h
Hlin 912.0A 1.502 h
H- H 0 0.248 h