Abstract

We suggest a numerical algorithm for complex ray tracing. Such an algorithm is intended for the computation of a wave field in the framework of complex geometrical optics. The main advantage of the complex method is the possibility to take into account diffraction effects by use of only ordinary differential equations of geometrical optics, thus reducing the calculation time. The efficiency of the suggested algorithm is illustrated by several numerical examples that allow comparison with known analytic solutions: the field of a plane wave behind a caustic in a linear layer, uniform field asymptotics on a caustic in a linear layer, and a Gaussian beam field in a homogeneous medium. It is pointed out that the approach under consideration can be readily applied to a great variety of real wave problems that have an analytical solution: nonplane waves, nonplane-stratified media, and the like. In particular, a numerical solution for Gaussian beam propagation through inhomogeneities of Gaussian form is presented.

© 2001 Optical Society of America

Diffraction of a Gaussian beam in a three-dimensional smoothly inhomogeneous medium: an eikonal-based complex geometrical-optics approach

Pawel Berczynski, Konstantin Yu. Bliokh, Yuri A. Kravtsov, and Andrzej Stateczny
J. Opt. Soc. Am. A 23(6) 1442-1451 (2006)

Gaussian beam and pulsed-beam dynamics: complex-source and complex-spectrum formulations within and beyond paraxial asymptotics

Ehud Heyman and Leopold B. Felsen
J. Opt. Soc. Am. A 18(7) 1588-1611 (2001)

Using rays better. I. Theory for smoothly varying media

G. W. Forbes and M. A. Alonso
J. Opt. Soc. Am. A 18(5) 1132-1145 (2001)

Dynamical approach to inhomogeneous rays

Sam L. Robinson
J. Opt. Soc. Am. A 13(10) 2021-2026 (1996)

Ray-based diffraction calculations using stable aggregates of flexible elements

Miguel A. Alonso
J. Opt. Soc. Am. A 30(6) 1223-1235 (2013)