There are two files labeled C1 & C2 as follows:
C1coronaIntroduction.doc
and C2coronaNitrogen.xls
The Introduction.doc file which is this file you have just opened describes:
1 Ionization Fractions, denoted by a’s
2. Effective charge of ions, denoted by Zeff and
3. Specific Heat Ratio, denoted by g.
Uses a Corona Model (such as discussed by H Griem in Plasma Spectroscopy 6-11) as a better approximation than the Saha Equations.
To run the programme:
Click on the programme file C2coronaNitrogen.xls.
The Active Sheet opens.
Type in the starting temperature and the maximum temperature to specify the range you want. To illustrate, change the T Max value to a small one such as 1,000,000.
Place the cursor in any unused cell. Operate ‘Ctrl’+’A’.
The code will calculate the values of a’s within the range you have specified, in this case, 10,000 to 1,000,000K.
Results are displayed in columns, T in K, T in eV, log T
Then follow the columns with the a’s 0-7 for Nitrogen,
Followed by a column for Zeff, a repeat column for TeV
And a separated column for the specific heat ratio g.
The graphs for Zeff and g are also displayed.
The gas has been set for Nitrogen.
To make
changes to the code:
Click on Tools
Click on Macro
Click on Macros
Click on’Step Into’
The code will be displayed.
Some notes on the code:
Dimension has been set as 54 to accommodate gases up to Xenon.
The Ionization Energies (in eV, each referenced to Ground State) of Nitrogen are entered as I(1)……I(7). M is set as 7.
If you wish to change to another gas, say Argon; you will need to set M=18
And key in the values of I(1) to I(18) using the correct values for argon.
Also go to the statements “For L= 0 to 6” and change to “For L=0-17”
And go to statements “For L=1 to 7” and change to “For L=1 to 18”
And go to statement “For L=0 to 7” and change to “For L=0 to 18”
After the changes, exit from code by clicking on red square “X” in upper right hand corner. Click OK to message “Command will stop the debugger”
The Active Sheet will appear and you are ready to compute with the revised code.
Notes on the Graphs
The graphs are formatted on the Active Sheet. The format may need to be adjusted when a new gas is computed. The formatting uses regular EXCEL techniques.