Two-reaction test case (start here)

graphThis test case corresponds to an Ar plasma consisting of electrons, atomic ions, and neutrals. The charged particles are supposed to be generated by direct electron impact ionization and lost by 3-body recombination:

  • e + Ar = e + e + Ar+, k1
  • e + Ar+ + Ar = Ar + Ar, k2

The BOLSIG+ solver is evoked to calculate k1, the rate coefficient for direct ionization. The recombination rate coefficient assumed to be constant k2 = 10-25 cm6/s, and the reduced field strength, E/N = 50 Td, is constant. This system has a simple asymptotic solution (see the Figure).

The input text file kinet.inp can be written using any editor. This input file is converted to a Fortran module using the PREPROCESSOR utility provided in the download package. An example of the Fortran code corresponding to above-mentioned kinet.inp file is zdplaskin_m.F90. The generated Fortran code is intended to be used as a “black box”.

An example of a user code test_2reac.F90 shows how the “black box” Fortran code could be incorporated in a master code.

Requred cross sections for electron processes with species indicated in section BOLSIG have to be stored in bolsigdb.dat data file. In the present example, these electron cross sections were retrieved from ELECTRON SCATTERING DATABASE, SIGLO database.

This test case must be compiled, using a Fortran 90/95 compiler, along with DVODE solver (dvode_f90_m.F90), ZDPLASKIN module (zdplaskin_m.F90) and BOLSIG+ library (bolsig.xxx) provided with the ZDPlasKin package. Screen output is in out.txt.

You can download all these files in a zip archive.