Nonlocal Boundary Conditions for Finite-Difference Parabolic Equation Solvers


  1. Yevick, D.
  2. Thomson, D.J.
Corporate Authors
Defence Research Establishment Atlantic, Dartmouth NS (CAN)
Theoretical/numerical models of underwater sound propagation usually incorporate a downgoing radiation condition on the transverse component of the acoustic field. Most one-way or parabolic equation (PE) solvers approximate this radiation condition by appending an absorbing layer to the computational mesh and setting the field to zero at the base of this layer. Papadakis et al.[ J. Acoust. Soc. Am. 92, 2030-2038 (1992) replaced such approximate treatments with a non-local boundary condition (NLBC) that exactly transforms the semi-infinite PE problem to an equivalent one in a bounded domain. Papadakis' approach requires the evaluation of a spectral (wavenumber) integral whose integrand is inversely proportional to the impedance of the subbottom medium. In this paper, an alternate procedure is analyzed for obtaining NLBCs directly from the z-space Crank-Nicolson solvers for both the Tappert and Claerbout PEs. Formulas for the field Psi at range r+increment r are derived in terms of the known field at the previously calculated range values from 0 to r by expanding the appropriate "vertical wave number" operator for the downgoing field in powers of the translation operator R=exp(increment r partial derivative r). The effectiveness of these NLBCs is examined for several numerical examples relevant to one-way underwater sound propagation.
Boundary conditions;Finite difference theory
Date of publication
06 Apr 1999
Number of Pages
Hardcopy;Document Image stored on Optical Disk

Permanent link

Document 1 of 1

Date modified: