An iterative method for including Doppler shift in global wave solvers using FEM decomposition

A method for calculating the wave field for spatial dispersive media is proposed suitable for FEM. The method is based on operator splitting by separating the induced current and wave field calculations, and solving the system by means of iterations. In order to take into account several coexisting waves with different poloidal mode numbers when calculating the induced current the wave field is decomposed into wavelets, for which the current is calculated assuming the plasma to be weakly non-uniform.