Abstract
The discrete Bessel representation in cylindrical (or spherical) geometry is presented that is analogous to the discrete Fourier scheme for a Cartesian coordinate. In other words, the integral Fourier–Bessel or Hankel transform is adequately discretized, and the underlying quadrature algorithm achieves Gaussian accuracy. The relevant discrete Bessel transform can be straightforwardly extended to the azimuthally dependent case. Efficient pseudospectral schemes are thus made possible in the resolution of the scalar Helmholtz and paraxial-wave equations in cylindrical coordinates. The superiority of the discrete Bessel scheme is illustrated by the calculation of the propagation constants for a multimode optical fiber.
© 1997 Optical Society of America
Full Article | PDF ArticleMore Like This
Charles Cerjan
J. Opt. Soc. Am. A 24(6) 1609-1616 (2007)
L. Niggl, T. Lanzl, and Max Maier
J. Opt. Soc. Am. A 14(1) 27-33 (1997)
Anthony A. Tovar and Glen H. Clark
J. Opt. Soc. Am. A 14(12) 3333-3340 (1997)