A hybrid fast Hankel transform algorithm for electromagnetic modeling

A hybrid fast Hankel transform algorithm has been developed that uses several complementary features of two existing algorithms: Anderson’s digital filtering or fast Hankel transform (FHT) algorithm and Chave’s quadrature and continued fraction algorithm. A hybrid FHT subprogram (called HYBFHT) written in standard Fortran-77 provides a simple user interface to call either subalgorithm. The hybrid approach is an attempt to combine the best features of the two subalgorithms in order to minimize the user’s coding requirements and to provide fast execution and good accuracy for a large class of electromagnetic problems involving various related Hankel transform sets with multiple arguments. Special cases of Hankel transforms of double‐order and double‐argument are discussed, where use of HYBFHT is shown to be advantageous for oscillatory kernel functions.

[1]  A. Oppenheim,et al.  An algorithm for the numerical evaluation of the Hankel transform , 1978, Proceedings of the IEEE.

[2]  J. Kauahikaua.  ELECTROMAGNETIC FIELDS ABOUT A HORIZONTAL ELECTRIC WIRE SOURCE OF ARBITRARY LENGTH , 1978 .

[3]  Gerald W. Hohmann,et al.  Electromagnetic modeling of three-dimensional bodies in layered earths using integral equations , 1983 .

[4]  Alan D. Chave,et al.  Numerical integration of related Hankel transforms by quadrature and continued fraction expansion , 1983 .

[5]  P. Cornille,et al.  Computation of Hankel Transforms , 1972 .

[6]  H. Johansen,et al.  AN INTERACTIVE COMPUTER/GRAPHIC-DISPLAY-TERMINAL SYSTEM FOR INTERPRETATION OF RESISTIVITY SOUNDINGS* , 1975 .

[7]  G. Keller,et al.  Frequency and transient soundings , 1983 .

[8]  A. Siegman.  Quasi fast Hankel transform. , 1977, Optics letters.

[9]  Eric W. Hansen,et al.  Fast Hankel transform algorithm , 1985, IEEE Trans. Acoust. Speech Signal Process..

[10]  S. Candel,et al.  An algorithm for the Fourier-Bessel transform , 1981 .

[11]  Walter L. Anderson,et al.  Computer program; numerical integration of related Hankel transforms of orders O and 1 by adaptive digital filtering , 1979 .

[12]  Stanley H. Ward,et al.  Inductive Sounding of a Layered Earth with a Horizontal Magnetic Dipole , 1970 .

[13]  H. Frank Morrison,et al.  ELECTROMAGNETIC FIELDS ABOUT A LOOP SOURCE OF CURRENT , 1970 .

[14]  D. Mook,et al.  An algorithm for the numerical evaluation of the Hankel and Abel transforms , 1983 .

[15]  Robert Piessens,et al.  Approximation for bessel-functions and their application in the computation of hankel-transforms , 1982 .

[16]  Alan D. Chave,et al.  Controlled electromagnetic sources for measuring electrical conductivity beneath the oceans: 1. Forward problem and model study , 1982 .

[17]  D. Ghosh,et al.  COMPUTATION OF TYPE CURVES FOR ELECTROMAGNETIC DEPTH SOUNDING WITH A HORIZONTAL TRANSMITTING COIL BY MEANS OF A DIGITAL LINEAR FILTER , 1972 .

[18]  D. Ghosh.  THE APPLICATION OF LINEAR FILTER THEORY TO THE ' DIRECT INTERPRETATION OF GEOELECTRICAL RESISTIVITY SOUNDING MEASUREMENTS * , 1971 .