Opened 22 months ago

Last modified 21 months ago

#358 new defect - wrong answer

PDR sims affected by induced processes

Reported by: Gary J. Ferland Owned by: Gary J. Ferland
Priority: blocker Milestone: C17_branch
Component: radiative transfer Version: trunk
Keywords: Cc:


tests reveal unexplained behavior in PDR sims with no ionizing radiation. Modify pdf_leiden_f4 to stop in zone 2. In its default state we get this hydrogen ionization for these two zones:

1.00E+00	5.18E-08
1.00E+00	3.51E-05 

the ionization goes up by ~3 dex.

Add the command "no induced processes" and we get

1.00E+00	5.18E-08
1.00E+00	5.19E-08

the answer does not change and is equal to the first zone, with induced processes. It is counterintuitive that induced processes would make the ionization higher in the second zone, which should be more shielded.

Attachments (1)

optimize-species-read.patch (3.7 KB) - added by mchatzikos 22 months ago.
patch for conditional loading of atomic data

Download all attachments as: .zip

Change History (3)

Changed 22 months ago by mchatzikos

Attachment: optimize-species-read.patch added

patch for conditional loading of atomic data

comment:1 Changed 22 months ago by mchatzikos

The patch that triggered the discovery of the unexplained behavior is attached (optimize-species-read.patch). It does not load atomic data for species that have been disabled.

comment:2 Changed 21 months ago by Gary J. Ferland

tried removing NO INDUCED and replacing with CASE B NO PHOTOIONIZATION. The latter disables photoionization of 2s. The large change in ionization still occurred, showing that photoionization of 2s is not responsible.

So we know: turning induced processes off and on disables / enables the anomalous ionization

it is not due to 2s metastables.

Note: See TracTickets for help on using tickets.