Steepest-descent algorithm for simulating plasma-wave caustics via metaplectic geometrical optics

Sean M Donnelly,Nicolas A Lopez,I Y Dodin,Sean M. Donnelly,Nicolas A. Lopez,I. Y. Dodin
DOI: https://doi.org/10.1103/PhysRevE.104.025304
IF: 2.707
2021-08-14
Physical Review E
Abstract:The design and optimization of radiofrequency-wave systems for fusion applications is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at wave cutoffs and caustics. To accurately model the wave behavior in these regions, more advanced and computationally expensive "full-wave" simulations are typically used, but this is not strictly necessary. A new generalized formulation called metaplectic geometrical optics (MGO) has been proposed that reinstates GO near caustics. The MGO framework yields an integral representation of the wavefield that must be evaluated numerically in general. We present an algorithm for computing these integrals using Gauss-Freud quadrature along the steepest-descent contours. Benchmarking is performed on the standard Airy problem, for which the exact solution is known analytically. The numerical MGO solution provided by the new algorithm agrees remarkably well with the exact solution and significantly improves on previously derived analytical approximations of the MGO integral. https://doi.org/10.1103/PhysRevE.104.025304 ©2021 American Physical Society
physics, fluids & plasmas, mathematical
What problem does this paper attempt to address?