Hankel transform based on QWE in the electromagnetic numerical computation
Hu Junhua1, Chang Yanjun1, Zuo Qiankun1, Sun Weibin2, Wu Guanghai3
1. Hubei Subsurface Multi-scale Imaging Key Laboratory, Institute of Geophysics and Geomatics, China University of Geosciences (Wuhan), Wuhan, Hubei 430074, China;
2. Integrated Geophysical and Geochemical Service, BGP, Zhuozhou, Hebei 072751, China;
3. The Second Institute, State Oceanic Administration, Hangzhou, Zhejiang 310000, China
Abstract:A high-performance numerical filtering method is usually needed to calculate Hankel integral in electromagnetic numerical simulations in order to obtain good calculation accuracy and computation speed. We propose in this paper a direct numerical integration method, quadrature-with-extrapolation (QWE) based on Shanks transform, and use it to compute one dimensional transient electromagnetic response of coincident-loop configuration and electromagnetic tensor Green's function of an electric dipole source. For the former, the QWE method can obtain more accurate imaginary components of vertical magnetic field at high frequencies and lead to more precise transient response in early and mediate times correspondingly in comparison with the fast Hankel transform (FHT) method, while its efficiency does not differ with FHT. For the latter, the QWE method with 5 nodes and 100 intervals has only 0.1% error, and saves half the time of the FHT computation under the same precision. These results show the QWE algorithm has the potential advantage and application value in the numerical computation of electromagnetic simulation.
Hohmann G W. Three-dimensional EM modeling.Geophysical Surveys, 1983, 6(1-2): 27-53.
[2]
Wannamaker P E, Hohmann G W, Sanfilipo W A. Electromagnetic modeling of three-dimensional bodiesin layered earths using integral equations. Geophy-sics, 1984, 49(1): 60-74.
[3]
Anderson W L. Computation of Green's tensor inte-grals for three-dimensional electromagnetic problemsusing fast Hankel transforms. Geophysics, 1984, 49(10): 1754-1759.
[4]
Guptasarma D, Singh B. New digital linear filters for Hankel J0 and J1 transforms. Geophysical Prospecting, 1997, 45(5): 745-762.
[5]
Longman I M. Note on a Method for Computing Infinite Integrals of Oscillatory Functions. Cambridge Univ Press, 1956.
[6]
Chave A D. Numerical integration of related Hankel transforms by quadrature and continued fraction expansion. Geophysics, 1983, 48(12): 1671-1686.
[7]
Michalski K A. Extrapolation methods for Sommerfeld integral tails. Antennas and Propagation, IEEE Transactions on, 1998, 46(10): 1405-1418.
[8]
Lucas S K. Evaluating infinite integrals involving products of Bessel functions of arbitrary order. Journal of Computational and Applied Mathematics, 1995, 64(3): 269-282.
[9]
Key K. Is the fast Hankel transform faster than quadrature? Geophysics, 2012, 77(3): F21-F30.
[10]
Singh S, Singh R. On the use of Shank's transform to accelerate the summation of slowly converging series. Microwave Theory and Techniques, IEEE Transactions on, 1991, 39(3): 608-610.