Due to the extremely large size of realistic three-dimensional ocean problems, there is an understandable desire to incorporate marching methods into the solution algorithm for the inherently elliptic (frequency-domain) problem. Recognizing that typical ocean propagation problems are essentially scattering problems in terms of a transition region and transversely inhomogeneous asymptotic half-spaces, wave field splitting, invariant imbedding, and phase space methods reformulate the problem in terms of an operator scattering matrix characteristic of the transition region. The subsequent equations for the reflection and transmission operators are first order in range, nonlinear (Riccati-like), and, in general, nonlocal. The system is well-posed, but stiff. The reflected and transmitted wave fields can be computed in a very efficient manner, while the wave field in the transition region (if desired) can be computed by essentially a layer-stripping algorithm. In principle, the transition region can be divided into subregions, allowing for parallel computations and subsequent recombination. Numerical examples will be discussed. [Work supported by NSF, AFOSR, ONR, ASEE, JEWC.]
Read full abstract