GE REDUCED PROBLEM
The reduced problem is from
Lawler & Kortshagen, J. Phys D: Appl Phys 32, pp 3188-3198 (1999),
"Self-consistent Monte Carlo simulations of the positive column of gas
discharges."
The reduced problem was run overnight on a 600 MHz Pentium III machine.
It ran for 40 microsecs. Equilibrium was reached at about 3 microsecs
as can seen from the Average Density(t) and Ez(t) diagnostics below. On the
600 MHz Pentium III machine this corresponds to under 1 hour of run time
to get equilibrium. Below are some figures, They are very close to the
expected output. The simulation starts off with initial Bessel function
electron and ion density profiles with axial density = 1.64e9/cc.
The following is a link to the input file:
Input File.
INPUTS
Gas is argon with model cross sections (one lumped excited state)
R=1 cm
Gas temperature = 273K
Gas pressure = 2.83 mTorr (2.83e-3 Torr)
Gas density = 1e14/cc (important!)
Axial current = 15.2 mA
EXPECTED OUTPUT (from Lawler and Kortshagen paper)
Ez=1.0 V/cm (from Table 1)
Electron average energy on axis = 14.7 V (from Table 1)
Electron density on axis = 1.64e9/cc (from Table 1)
Wall potential = 37 V (by inspection of Figure 1)
XPDC1 OUTPUT
Ez=0.65 V/cm
Electron average energy on axis = 12 V
Electron density on axis = 1.75e9/cc
Wall potential = 27 V
