Numerical analysis of eikonal equation

D. S. Kulyabov,A. V. Korolkova,T. R. Velieva,M. N. Gevorkyan

DOI: https://doi.org/10.1117/12.2525142

2019-06-23

Abstract:The Maxwell equations have a fairly simple form. However, finding solutions of Maxwell's equations is an extremely difficult task. Therefore, various simplifying approaches are often used in optics. One such simplifying approach is to use the approximation of geometric optics. The approximation of geometric optics is constructed with the assumption that the wavelengths are small (short-wavelength approximation). The basis of geometric optics is the eikonal equation. The eikonal equation can be obtained from the wave equation (Helmholtz equation). Thus, the eikonal equation relates the wave and geometric optics. In fact, the eikonal equation is a quasi-classical approximation (the Wentzel-Kramers-Brillouin method) of wave optics. This paper shows the application of geometric methods of electrodynamics to the calculation of optical devices, such as Maxwell and Luneburg lenses. The eikonal equation, which was transformed to the ODE system by the method of characteristics, is considered. The resulting system is written for the case of Maxwell and Luneburg lenses.

Computational Physics,Numerical Analysis

What problem does this paper attempt to address?

This paper aims to solve how to use the characteristics method in geometric optics to solve the eikonal equation, so as to calculate and simulate the propagation paths of rays and wavefronts in optical devices (such as Maxwell lenses and Luneberg lenses). Specifically, the paper converts the ray equation into a system of ordinary differential equations (ODE system) and uses numerical methods (such as the classic Runge - Kutta method) to solve these equations, thereby achieving accurate modeling and visualization of the beam trajectories and wavefronts inside the lenses.

The main contributions of the paper are as follows:
1. **Mathematical transformation**: It describes in detail how to convert the ray equation into a system of ordinary differential equations, and this process involves the application of the characteristics method.
2. **Numerical simulation**: It uses the Julia programming language to realize the numerical simulation of the beam trajectories inside Maxwell lenses and Luneberg lenses, and provides specific numerical methods and initial condition settings.
3. **Result visualization**: It shows the images of the beam trajectories and wavefronts emitted by point light sources under different lenses, which verifies the effectiveness of the proposed method.

Through the above methods, the paper not only solves the numerical solution problem of the ray equation in specific optical devices, but also provides an effective computational tool and reference for researchers in related fields.