Index of /soft/plasma

Evgenii B. Rudnyi

 Name                    Last modified       Size  

 Parent Directory        25-Mar-2009 15:34      -  
 al.dat                  02-Dec-1993 01:17     8k  
 c2f4.dat                13-Jun-1994 23:44    17k  
 getmatr1.for            16-Sep-1994 17:37     3k  
 gettsiv1.for            11-Sep-1994 18:51     1k  
 kf.dat                  19-Apr-1996 14:23    23k  
 mo.dat                  15-Sep-1994 17:34    24k  
 mo.for                  31-Dec-2004 13:52    12k  
 mo.old.for              02-Feb-1995 11:37    11k  
 read.mo.txt             30-Dec-2004 13:00     4k  
 testRead.for            31-Dec-2004 11:46     1k  
 w.dat                   19-Nov-1994 13:04    24k

http://Evgenii.Rudnyi.Ru

Evgenii B. Rudnyi, 2004
Use at your own risk.


Influence of Negative Ions on Electron Concentration in Air Plasma
==================================================================

The archive contains a sample code that I have written in 1993-1998 to
compute equilibrium composition for my papers
(http://Evgenii.Rudnyi.Ru/publications.html#plasma).

Requirements: VCS subroutine from [1]. Insert COMMON /MY/ in its body.

mo.old.for is the original code that was written for Lahey compiler. 

I have fixed the code and now mo.for can be compiled by g77 3.4 on my
Mac. I failed to convert

			DATA (E(I), I = 1, Natom) /'e', 'O', 'H', 'N', 'C', 'Mo', 'K'/

to g77 and I have commented this statement out. Oh, I have forgotten
Fortran. But this is not important. This was to make an internal VCS
output nicer.

Much more puzzling is that I needed to change format from I2 to F2.0
to read DOUBLE PRECISION correctly. testRead.for shows what happens
when you do not do it. Weird, isn't it?
 
Basically, my program is a simple loop to run VCS several times to
produce a text file with results for plotting. File read.mo.txt
explains the program behavior.

The program reads the Gibbs energies from the file. In order to make
the data file, I have used gettsiv1.for to convert the output of CHET,
and then getmatrix1.for to insert a formula matrix. This is just for
your information. All data files are already in the final form.

al.dat contains the Gibbs energies for [2]
c2f4.dat contains the Gibbs energies for [3]
mo.dat and w.dat contain the Gibbs energies for [4]
kf.dat contains the Gibbs energies for [5]

mo.for is tuned to work with mo.dat. You may need to modify it to work 
with other files (file name, PARAMETER, DATA, and probably something else).

Feel free to contact me if you have any question. My e-mail address 
can be found at http://Evgenii.Rudnyi.Ru/. 

[1] Smith W.R., Missen R.W. Chemical Reaction Equilibrium Analysis:
Theory and Algorithms. 1982.

[2] N. N. Pilyugin, E. B. Rudnyi. Calculation of the equilibrium
composition of air plasma with product of ablation of aluminum for
conditions of ballistic experiments. Tepl. Vys. Temp., 1994, v. 32, N
4, p. 511-517 (in Russian). High Temp. Russia, 1994, v. 32, N 4, p.
476-482 (English transl.).

[3] N. N. Pilyugin, E. B. Rudnyi. Calculation of the equilibrium
composition of air plasma with admixture of C2F4 for conditions of
ballistic experiments. Tepl. Vys. Temp., 1995, v. 33, N 4, p. 532-538
(in Russian). High Temp. Russia, 1995, v. 33, N 4 (English transl.).

[4] N. N. Pilyugin, E. B. Rudnyi. Influence of ablation on the
electron concentration in air plasma. Tepl. Vys. Temp., 1996, V. 34, N
3, p. 355-364 (in Russian). High Temp. Russia, 1996, v. 34, N 3
(English transl.).

[5] N. N. Pilyugin, V. V. Kuz'menko, E. B. Rudnyi. Combined Effect of
Impurities of Water Vapors, Potassium and Flourine on the Electron
Number Density of Air Plasma. Tepl. Vys. Temp., 1998, v. 36, N 2, p.
201-205 (in Russian). High Temp., v. 36, N 2, p. 183-188 (English
transl.).