Complex geometric optics theory applied to an inhomogeneous plasma: A numerical model

By introducing a complex dispersion relation, the scope of applicability of the Geometric Optics approximation for the propagation of electromagnetic wave in a nondegenerate inhomogeneous plasma is extended to include absorption of the propagating waves at harmonics of the ion cyclotron frequency. Complex ray-tracing equations are obtained that are Hamiltonian in character. The equations describing the evolution of the wave amplitude and its phase along the ray paths are introduced as well. As complement to the theoretical work, the energy density conservation law is obtained. Since the resulting equations are complex, they are analytically intractable, so a numerical solution seems to be the only reasonable approach that might be taken at this stage. For this purpose, a ray-tracing code incorporating ion magnetization and focusing effects is developed. To accomplish this, a very accurate asymptotic expression for the dispersion relation is derived and employed to solve numerically the ray-tracing and wave-amplitude equations