Opened 5 years ago

Last modified 5 years ago

#254 new defect - convergence

user sim does not conserve energy

Reported by: Gary J. Ferland Owned by: nobody
Priority: major Milestone:
Component: etc Version: trunk
Keywords: Cc:

Description (last modified by Gary J. Ferland)

user alexei.baskin on yahoo reported following on 2013 Jan 13:


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 – 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 (, 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

Attachments (1) (91.3 KB) - added by Gary J. Ferland 5 years ago.
slightly modified version of original post

Download all attachments as: .zip

Change History (3)

Changed 5 years ago by Gary J. Ferland

Attachment: added

slightly modified version of original post

comment:1 Changed 5 years ago by Gary J. Ferland

Description: modified (diff)

comment:2 Changed 5 years ago by Gary J. Ferland

further note on c13

Hi Gary,

A quick update regarding trying my sims on v13 RC1. I ran a sim with more (reasonable) 
AGN-like parameters i.e., log U=3 and log n=8. The code still fails to converge 
pressure on a second iteration. The pressure "oscillates" by a factor of ~4 between 1 
and 4, while density and temperature remain approximately constant (see for example 
last zone in the second iteration in the attached files). For this sim, the line 
pressure contribution to the total pressure jumps from ~30-50% to 90% in adjacent 
zones (see the end of line pressure file).

BTW, running the attached sim on v10 also fails, but this time because the calculation 
stops due to zone thickness becoming too thin.

Alexei Baskin
Note: See TracTickets for help on using tickets.