D3DEnergetic / FIDASIM

A Neutral Beam and Fast-ion Diagnostic Modeling Suite
http://d3denergetic.github.io/FIDASIM/
Other
29 stars 19 forks source link

Aladdin cross sections for E < 2 keV #148

Closed lstagner closed 6 years ago

lstagner commented 6 years ago

This implements the changes to the atomic cross sections discussed at the last meeting. The change is to the Hydrogen Charge exchange from n=1->m=2 cross section below 2 keV. The lower energy limit of the new cross sections fit is 19 eV. Since we usually provide tables down to 1 eV the cross section was linearly extrapolated in log-space. The cross sections are also no longer normalized to the total Janev cross sections. The new cross sections are shown below. Solid lines are the old cross sections, dashed lines are the new cross sections, and the fine-dashed line is linear extrapolation.

image

Because of the way we fill out the CX table this change has some rippling effects. Shown below are the other cross sections that are affected by this change. image

image

I also used this as an opportunity to upgrade the table calculations. You can now specify which tables you want to calculate and you can now specify more than one impurity charge so you can get tables for multiple impurity species.

@geigerb @shaunhaskey

lstagner commented 6 years ago

I added back the normalization to the Janev total cross sections to avoid the differences at high energy. With the normalization, the only difference to the tables is for the cx from the n=1->2 (new cross sections are shown below (dashed lines)) At a later date, we should do a survey of the cross sections we use compared to those suggested by the Aladdin database.

I will do a quick benchmark to see the effects of this change and I will merge this shortly thereafter. image

lstagner commented 6 years ago

From the quick benchmark that I did this seems to have little effect. I will go ahead and merge this now.