Abstract A comparison is made of solutions of the continuity equations for the motion of electrons and ions in a strong electric field using the methods of Euler, Runge-Kutta, Lax-Wendroff, characteristics and the flux-corrected transport (FCT) algorithm “Phoenical LPE Shasta” developed by Boris and Book( J. Comput. Phys. 20 (1976) , 397), with their flux limiter and with the new flux limiting algorithms developed by Zalesak ( J. Comput. Phys. 31 (1979) , 335). Results with and without ionization, and for uniform and non-uniform electric fields, are compared. Because ionization by electrons is strongly dependent on the electric field, considerable care needs to be taken to choose an optimum numerical scheme for non-uniform fields. For example, the calculation of the properties of corona discharges, or the cathode-fall of a glow discharge, requires fine spatial resolution. It is found that the Lax-Wendroff method and the method of characteristics give acceptable results; however, the Phoenical LPE Shasta algorithm with the flux limiting algorithm of Zalesak gives the best results, with the added advantage of suppressing spurious oscillations.
[1]
S. Zalesak.
Fully multidimensional flux-corrected transport algorithms for fluids
,
1979
.
[2]
David Potter.
Computational physics
,
1973
.
[3]
J. Boris,et al.
Flux-corrected transport. III. Minimal-error FCT algorithms
,
1976
.
[4]
E. Turkel,et al.
On the Practical Use of High-Order Methods for Hyperbolic Systems*
,
1980
.
[5]
L. Kline.
Calculations of discharge initiation in overvolted parallel‐plane gaps
,
1974
.
[6]
R. Wooldridge,et al.
Advanced Engineering Mathematics
,
1967,
The Mathematical Gazette.
[7]
Gilbert Strang,et al.
Accurate partial difference methods
,
1964
.
[8]
R. D. Richtmyer,et al.
Difference methods for initial-value problems
,
1959
.
[9]
A. J. Davies,et al.
THE THEORY OF IONIZATION GROWTH IN GASES UNDER PULSED AND STATIC FIELDS
,
1973
.