The operation of various types of space plasma thrusters relies on heating the plasma with electromagnetic waves, typically in the MHz or GHz range. This is the case of electrodeless plasma thrusters.
To understand and characterize the propagation and absorption of electromagnetic power by the magnetoplasma, a commonly-used approach is to model the plasma as a linear, anisotropic, gyrotropic, dissipative medium based on the cold-plasma assumption.
Simulation codes ranging from finite differences (e.g. Yee scheme) to finite elements (e.g. Nédélec edge elements) are used to find the response of the wavefields to a given geometry, plasma density, and applied magnetic field profiles.
This presentation reviews recent efforts at Universidad Carlos III de Madrid to simulate electrodeless plasma thrusters with finite element schemes, and the occurrence of spurious wave solutions whenever the propagation regime includes so-called resonant cones. To characterize the spurious solutions, we set up a verification simulation case based on the Fresnel problem for a planar wave propagating in vacuum and into a plasma, and compare against the analytical solution. We find perfect agreement with our finite element solution, except for parametric regimes where the plasma has resonant cones.
Analysis of the analytical and numerical dispersion relations helps identify the sources of error, and potential countermeasures are proposed.