PETOOL: MATLAB-based one-way and two-way split-step parabolic equation tool for radiowave propagation over variable terrain

Abstract A MATLAB-based one-way and two-way split-step parabolic equation software tool (PETOOL) has been developed with a user-friendly graphical user interface (GUI) for the analysis and visualization of radio-wave propagation over variable terrain and through homogeneous and inhomogeneous atmosphere. The tool has a unique feature over existing one-way parabolic equation (PE)-based codes, because it utilizes the two-way split-step parabolic equation (SSPE) approach with wide-angle propagator, which is a recursive forward–backward algorithm to incorporate both forward and backward waves into the solution in the presence of variable terrain. First, the formulation of the classical one-way SSPE and the relatively-novel two-way SSPE is presented, with particular emphasis on their capabilities and the limitations. Next, the structure and the GUI capabilities of the PETOOL software tool are discussed in detail. The calibration of PETOOL is performed and demonstrated via analytical comparisons and/or representative canonical tests performed against the Geometric Optic (GO) + Uniform Theory of Diffraction (UTD). The tool can be used for research and/or educational purposes to investigate the effects of a variety of user-defined terrain and range-dependent refractivity profiles in electromagnetic wave propagation. Program summary Program title: PETOOL (Parabolic Equation Toolbox) Catalogue identifier: AEJS_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEJS_v1_0.html Program obtainable from: CPC Program Library, Queenʼs University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 143 349 No. of bytes in distributed program, including test data, etc.: 23 280 251 Distribution format: tar.gz Programming language: MATLAB (MathWorks Inc.) 2010a. Partial Differential Toolbox and Curve Fitting Toolbox required Computer: PC Operating system: Windows XP and Vista Classification: 10 Nature of problem: Simulation of radio-wave propagation over variable terrain on the Earthʼs surface, and through homogeneous and inhomogeneous atmosphere. Solution method: The program implements one-way and two-way Split-Step Parabolic Equation (SSPE) algorithm, with wide-angle propagator. The SSPE is, in general, an initial-value problem starting from a reference range (typically from an antenna), and marching out in range by obtaining the field along the vertical direction at each range step, through the use of step-by-step Fourier transformations. The two-way algorithm incorporates the backward-propagating waves into the standard one-way SSPE by utilizing an iterative forward–backward scheme for modeling multipath effects over a staircase-approximated terrain. Unusual features: This is the first software package implementing a recursive forward–backward SSPE algorithm to account for the multipath effects during radio-wave propagation, and enabling the user to easily analyze and visualize the results of the two-way propagation with GUI capabilities. Running time: Problem dependent. Typically, it is about 1.5 ms (for conducting ground) and 4 ms (for lossy ground) per range step for a vertical field profile of vector length 1500, on Intel Core 2 Duo 1.6 GHz with 2 GB RAM under Windows Vista.

[1]  P. Liu,et al.  Numerical modelling of wave propagation using parabolic approximation with a boundary-fitted co-ordinate system , 1989 .

[2]  R. H. Hardin.  Application of the split-step Fourier method to the numerical solution of nonlinear and variable coefficient wave equations , 1973 .

[3]  Parabolic equation techniques for scattering , 2000 .

[4]  M. Feit,et al.  Light propagation in graded-index optical fibers. , 1978, Applied optics.

[5]  J. Keller,et al.  Geometrical theory of diffraction. , 1962, Journal of the Optical Society of America.

[6]  James R. Wait.  Transmission and reflection of electromagnetic waves in the presence of stratified media , 1958 .

[7]  G. D. Dockery.  Modeling electromagnetic wave propagation in the troposphere using the parabolic equation , 1988 .

[8]  C. Mias,et al.  Fast Computation of the Nonlocal Boundary Condition in Finite Difference Parabolic Equation Radiowave Propagation Simulations , 2008, IEEE Transactions on Antennas and Propagation.

[9]  Sherman W. Marcus.  A hybrid (finite difference-surface Green's function) method for computing transmission losses in an inhomogeneous atmosphere over irregular terrain , 1992 .

[10]  M. Skolnik,et al.  Introduction to Radar Systems , 2021, Advances in Adaptive Radar Detection and Range Estimation.

[11]  Richard B. Evans,et al.  A two‐way parabolic equation for acoustic backscattering in the ocean , 1992 .

[12]  Levent Sevgi,et al.  Complex electromagnetic problems and numerical simulation approaches , 2001 .

[13]  Levent Sevgi,et al.  A MATLAB-based two-dimensional parabolic equation radiowave propagation package , 2005 .

[14]  L. Felsen,et al.  Wave propagation inside a two-dimensional perfectly conducting parallel-plate waveguide: hybrid ray-mode techniques and their visualizations , 2004, IEEE Antennas and Propagation Magazine.

[15]  K. H. Craig,et al.  Propagation modelling in the troposphere: parabolic equation method , 1988 .

[16]  Michael D. Collins,et al.  Two-way parabolic equation techniques for diffraction and scattering problems☆ , 2000 .

[17]  J. Kuttler,et al.  Differences between the narrow-angle and wide-angle propagators in the split-step Fourier solution of the parabolic wave equation , 1999 .

[18]  G. Daniel Dockery,et al.  Theoretical description of the parabolic approximation/Fourier split-step method of representing electromagnetic propagation in the troposphere , 1991 .

[19]  M. Levy.  Parabolic Equation Methods for Electromagnetic Wave Propagation , 2000 .

[20]  Ramakrishna Janaswamy,et al.  A curvilinear coordinate-based split-step parabolic equation method for propagation predictions over terrain , 1998 .

[21]  D. E. Kerr,et al.  Propagation of Short Radio Waves , 1989 .

[22]  N. R. Chapman,et al.  A wide‐angle split‐step algorithm for the parabolic equation , 1983 .

[23]  Mustafa Kuzuoglu,et al.  Two-way fourier split step algorithm over variable terrain with narrow and wide angle propagators , 2010, 2010 IEEE Antennas and Propagation Society International Symposium.

[24]  L. Sevgi,et al.  The Split-Step-Fourier and Finite-Element-Based Parabolic-Equation Propagation-Prediction Tools: Canonical Tests, Systematic Comparisons, and Calibration , 2010, IEEE Antennas and Propagation Magazine.

[25]  M. D. Collins.  A split‐step Padé solution for the parabolic equation method , 1993 .

[26]  D. J. Donohue,et al.  Propagation modeling over terrain using the parabolic wave equation , 2000 .

[27]  Pyotr Yakovlevich Ufimtsev.  Theory of edge diffraction in electromagnetics , 2003 .

[28]  Mustafa Kuzuoglu,et al.  A Novel Two-Way Finite-Element Parabolic Equation Groundwave Propagation Tool: Tests With Canonical Structures and Calibration , 2011, IEEE Transactions on Geoscience and Remote Sensing.

[29]  Jon F. Claerbout,et al.  Fundamentals of Geophysical Data Processing: With Applications to Petroleum Prospecting , 1985 .

[30]  A. Barrios.  A terrain parabolic equation model for propagation in the troposphere , 1994 .

[31]  L. Sevgi,et al.  Groundwave Modeling and Simulation Strategies and Path Loss Prediction Virtual Tools , 2007, IEEE Transactions on Antennas and Propagation.

[32]  M. F. Levy.  Diffraction studies in urban environment with wide-angle parabolic equation method , 1992 .

[33]  Peter D. Holm,et al.  Wide-Angle Shift-Map PE for a Piecewise Linear Terrain—A Finite-Difference Approach , 2007 .

[34]  Young-Ho Sohn,et al.  Light Propagation in Multimode GRIN(graded-index) Fibers with Intrusion Sensing Capability , 2002 .

[35]  Dehua Huang.  Finite element solution to the parabolic wave equation , 1988 .

[36]  M. F. Levy.  Parabolic equation modelling of propagation over irregular terrain , 1990 .

[37]  H. A. Dell,et al.  Introduction to Radar Systems. Merrill I. Skolnik. McGraw-Hill Book Co., London and New York. 1962. 648 pp. Illustrated. £5 12s. 6d. , 1963, The Journal of the Royal Aeronautical Society.

[38]  L. Sevgi,et al.  FEM-Based Surface Wave Multimixed-Path Propagator and Path Loss Predictions , 2009, IEEE Antennas and Wireless Propagation Letters.

[39]  R. Kouyoumjian,et al.  A uniform geometrical theory of diffraction for an edge in a perfectly conducting surface , 1974 .

[40]  L. Sevgi,et al.  Numerical Investigations of and Path Loss Predictions for Surface Wave Propagation Over Sea Paths Including Hilly Island Transitions , 2010, IEEE Transactions on Antennas and Propagation.

[41]  K. G. Budden.  The Wave-guide mode theory of wave propagation , 1963 .

[42]  O. Ozgun,et al.  Recursive Two-Way Parabolic Equation Approach for Modeling Terrain Effects in Tropospheric Propagation , 2009, IEEE Transactions on Antennas and Propagation.

[43]  Mustafa Kuzuoglu,et al.  Two-way split-step parabolic equation algorithm for tropospheric propagation: Tests and comparisons , 2010, 2010 10th Mediterranean Microwave Symposium.