An unconditionally stable time-domain discretization on cartesian meshes for the simulation of nonuniform magnetized cold plasma

In this paper an unconditionally stable, spatially and temporally implicit time-domain discretization for nonuniform magnetized cold plasma is developed. The discrete dispersion relation is free of spurious solutions and approximates the continuous dispersion relation for well-resolved wavelengths and frequencies (kΔ ≪ π, ωΔt ≪ π). For a specific choice of parameters, the discrete dispersion relation approximates the continuous dispersion relation for all wavelengths and frequencies up to the Nyquist limit. A few examples, amongst them one involving mode conversion, illustrate the new method.