Two-way Fourier Split Step Algorithm over Variable Terrain with Narrow and Wide Angle Propagators

Ozgun, Ozlem
Apaydin, Gökhan
Kuzuoğlu, Mustafa
Sevgi, Levent
Helmholtz's wave equation can be approximated by means of two differential equations, corresponding to forward and backward propagating waves each of which is in parabolic wave equation (PWE) form. The standard PWE is very suitable for marching-type numerical solutions. The one-way Fourier split-step parabolic equation algorithm (SSPE) is highly effective in modeling electromagnetic (EM) wave propagation above the Earth's irregular surface through inhomogeneous atmosphere. The two drawbacks of the standard PWE are: (i) It handles only the forward-propagating waves, and cannot account for the backscattered ones. The forward waves are usually adequate for typical longrange propagation scenarios. However, the backward waves become significant in the presence of obstacles that redirect the incoming wave. Hence, this necessitates the accurate estimation of the multipath effects to model the tropospheric wave propagation over terrain, (ii) It is a narrow-angle approximation, which consequently restricts the accuracy to propagation angles up to 10°-15° from the paraxial direction. To handle propagation angles beyond these values, wide-angle propagators have been introduced. Recently, a two-way SSPE algorithm was implemented to incorporate the backwardpropagating waves into the standard one-way SSPE, through a recursive forwardbackward scheme to model the tropospheric electromagnetic propagation over a staircase-approximated terrain. This algorithm has employed the standard narrowangle propagators in its implementation. The primary goal of this paper is to present the improved version of the algorithm based on wide-angle propagators, and to demonstrate the results of the comparison tests performed in some canonical scenarios, together with more complex scenarios involving variable terrains.


Two methods for the numerical integration of the time-dependent Schrodinger equation with given initial conditions (initial wave packet) are presented. The first method (method A) is based on the Schrodinger representation of the quantum-dynamical system while the second one (method B) is based upon the intermediate representation. In both cases the quantum dynamical equation is transformed into a system of Hamilton-Jacobi type equations of motion as occurring in multi particle classical dynamics, i.e. stan...
Improved analytical approximation to arbitrary l-state solutions of the Schrodinger equation for the hyperbolical potential
IKHDAİR, SAMEER; Sever, Ramazan (2009-04-01)
A new approximation scheme to the centrifugal term is proposed to obtain the l not equal 0 bound-state solutions of the Schrodinger equation for an exponential-type potential in the framework of the hypergeometric method. The corresponding normalized wave functions are also found in terms of the Jacobi polynomials. To show the accuracy of the new proposed approximation scheme, we calculate the energy eigenvalues numerically for arbitrary quantum numbers n and l with two different values of the potential par...
Two-level finite element method with a stabilizing subgrid for the incompressible MHD equations
Aydın Bayram, Selma; Nesliturk, A. I.; Tezer, Münevver (2010-01-20)
We consider the Galerkin finite element method (FEM) for the incompressible magnetohydrodynamic (MHD) equations in two dimension. The domain is discretized into a set of regular triangular elements and the finite-dimensional spaces employed consist of piecewise continuous linear interpolants enriched with the residual-free bubble functions. To find the bubble part of the solution, a two-level FEM with a stabilizing subgrid of a single node is described and its application to the MHD equations is displayed. ...
Two approximation schemes to the bound states of the Dirac-Hulthen problem
IKHDAİR, SAMEER; Sever, Ramazan (IOP Publishing, 2011-09-02)
The bound-state (energy spectrum and two-spinor wavefunctions) solutions of the Dirac equation with the Hulthen potential for all angular momenta based on the spin and pseudospin symmetry are obtained. The parametric generalization of the Nikiforov-Uvarov method is used in the calculations. The orbital dependence (spin-orbit-and pseudospin-orbit-dependent coupling too singular 1/r(2)) of the Dirac equation are included to the solution by introducing a more accurate approximation scheme to deal with the cent...
A Karhunen-Loeve-based approach to numerical simulation of transition in Rayleigh-Benard convection
Tarman, HI (2003-06-01)
A Karhunen-Loeve ( K - L) basis is generated empirically, using a database obtained by numerical integration of Boussinesq equations representing Rayleigh - Benard convection in a weakly turbulent state in a periodic convective box with free upper and lower surfaces. This basis is then used to reduce the governing partial differential equation (PDE) into a truncated system of amplitude equations under Galerkin projection. In the generation and implementation of the basis, the symmetries of the PDE and the g...
Citation Formats
O. Ozgun, G. Apaydin, M. Kuzuoğlu, and L. Sevgi, “Two-way Fourier Split Step Algorithm over Variable Terrain with Narrow and Wide Angle Propagators,” 2010, Accessed: 00, 2020. [Online]. Available: