Generalization of the rotated parabolic equation to variable slopes

The rotated parabolic equation [J. Acoust. Soc. Am. 87, 1035-1037 (1990)] is generalized to problems involving ocean-sediment interfaces of variable slope. The approach is based on approximating a variable slope in terms of a series of constant slope regions. The original rotated parabolic equation algorithm is used to march the field through each region. An interpolation-extrapolation approach is used to generate a starting field at the beginning of each region beyond the one containing the source. For the elastic case, a series of operators is applied to rotate the dependent variable vector along with the coordinate system. The variable rotated parabolic equation should provide accurate solutions to a large class of range-dependent seismo-acoustics problems. For the fluid case, the accuracy of the approach is confirmed through comparisons with reference solutions. For the elastic case, variable rotated parabolic equation solutions are compared with energy-conserving and mapping solutions.     

Parabolic equation solution of seismo-acoustics problems involving variations in bathymetry and sediment thickness

Recent improvements in the parabolic equation method are combined to extend this approach to a larger class of seismo-acoustics problems. The variable rotated parabolic equation [J. Acoust. Soc. Am. 120, 3534-3538 (2006)] handles a sloping fluid-solid interface at the ocean bottom. The single-scattering solution [J. Acoust. Soc. Am. 121, 808-813 (2007)] handles range dependence within elastic sediment layers. When these methods are implemented together, the parabolic equation method can be applied to problems involving variations in bathymetry and the thickness of sediment layers. The accuracy of the approach is demonstrated by comparing with finite-element solutions. The approach is applied to a complex scenario in a realistic environment.     

A mapping approach for handling sloping interfaces

A mapping approach for handling sloping interfaces in parabolic equation solutions is developed and tested. At each range, the medium is rigidly translated vertically so that a sloping interface becomes horizontal. To simplify the approach, the slope is assumed to be small and the extra terms that arise in the wave equation under the mapping are neglected. The effects of these terms can be approximately accounted for by applying a leading-order correction to the phase. The mapping introduces variations in topography, which are relatively easy to handle for the case of a pressure-release boundary condition. The accuracy of the approach is demonstrated for problems involving fluid sediments. The approach should also be accurate for problems involving elastic sediments and should be useful for solving three-dimensional problems involving variable topography.     

An energy-conserving spectral solution

An energy-conserving spectral solution is derived and tested. A range-dependent medium is approximated by a sequence of range-independent regions. In each region, the acoustic field is represented in terms of the horizontal wave-number spectrum. A condition corresponding to energy conservation is derived for the vertical interfaces between regions. The accuracy of the approach is demonstrated for problems involving sloping ocean bottoms. The energy-conserving spectral solution is less efficient than the energy-conserving parabolic equation solution, but might be suitable for generalization to problems involving elastic bottoms.     

Single-scattering parabolic equation solutions for elastic media propagation, including Rayleigh waves

The parabolic equation method with a single-scattering correction allows for accurate modeling of range-dependent environments in elastic layered media. For problems with large contrasts, accuracy and efficiency are gained by subdividing vertical interfaces into a series of two or more single-scattering problems. This approach generates several computational parameters, such as the number of interface slices, an iteration convergence parameter τ, and the number of iterations n for convergence. Using a narrow-angle approximation, the choices of n=1 and τ=2 give accurate solutions. Analogous results from the narrow-angle approximation extend to environments with larger variations when slices are used as needed at vertical interfaces. The approach is applied to a generic ocean waveguide that includes the generation of a Rayleigh interface wave. Results are presented in both frequency and time domains.     

A two-way parabolic equation that accounts for multiple scattering

A two-way parabolic equation that accounts for multiple scattering is derived and tested. A range-dependent medium is divided into a sequence of range-independent regions. The field is decomposed into outgoing and incoming fields in each region. The conditions between vertical interfaces are implemented using rational approximations for the square root of an operator. Rational approximations are also used to relate fields between neighboring interfaces. An iteration scheme is used to solve for the outgoing and incoming fields at the vertical interfaces. The approach is useful for solving problems involving scattering from waveguide features and compact objects.     

Three-dimensional parabolic equation model for seismo-acoustic propagation: Theoretical development and preliminary numerical implementation

A three-dimensional(3D) parabolic equation(PE) model for sound propagation in a seismo-acoustic waveguide is developed in Cartesian coordinates, with x, y, and z representing the marching direction, the longitudinal direction, and the depth direction, respectively. Two sets of 3D PEs for horizontally homogenous media are derived by rewriting the 3D elastic motion equations and simultaneously choosing proper dependent variables. The numerical scheme is for now restricted to the y-independent bathymetry. Accuracy of the numerical scheme is validated, and its azimuthal limitation is analyzed. In addition, effects of horizontal refraction in a wedge-shaped waveguide and another waveguide with a polyline bottom are illustrated. Great efforts should be made in future to provide this model with the ability to handle arbitrarily irregular fluid-elastic interfaces.     

A single-scattering correction for large contrasts in elastic layers

The single-scattering solution is implemented in a formulation that makes it possible to accurately handle solid-solid interfaces with the parabolic equation method. Problems involving large contrasts across sloping stratigraphy can be handled by subdividing a vertical interface into a series of two or more scattering problems. The approach can handle complex layering and is applicable to a large class of seismic problems. The solution of the scattering problem is based on an iteration formula, which has improved convergence in the new formulation, and the transverse operator of the parabolic wave equation, which is implemented efficiently in terms of banded matrices. Accurate solutions can often be obtained by using only one iteration.     

Rayleigh wave at the boundary of an inhomogeneous nonlinear viscoelastic half-space

In the approximation of weak nonlinearity and weak viscosity of the medium, we obtain an equation describing the spectral density of the particle horizontal velocity for a Rayleigh wave propagating along the boundary of a half-space. The coefficients of nonlinear interaction between the wave harmonics are found on the assumption that the third-order elastic moduli arbitrarily depend on the depth. We find expressions for the complex correction to the wave frequency due to small relaxation corrections to the elastic moduli and small variations in the medium density, which arbitrarily depend on the depth as well. The imaginary part of this correction to the frequency determines the decay of the linear Rayleigh wave due to small relaxation corrections to the elastic moduli arbitrarily dependent on the depth. Using numerical simulation (with allowance for the interaction of 500 harmonics), we study distortions of an initially harmonic Rayleigh wave for a particular dependence of variations in the nonlinear moduli on the depth. An integral equation is derived for the nonlinear elastic moduli as functions of the depth. It is shown that for independent spatio-temporal distributions of the viscous moduli, functions determining the dependence of the viscosity on the depth are described by an analogous integral equation. __________ Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 50, No. 3, pp. 212–226, March 2007.     

海洋声学中三维抛物方程非均匀网格模型

在海洋声学中,三维抛物方程模型可以有效考虑三维空间的声传播效应。然而,采用三维抛物方程模型分析三维空间内的声传播问题时,计算时间较长,并且需要消耗较大的计算机内存,因此给远距离声场的快速精确计算带来了很大困难。为此,将非均匀网格Galerkin离散化方法用于三维直角坐标系下的水声抛物方程模型中,深度算子和水平算子Galerkin离散方式由均匀网格变为非均匀网格。仿真结果表明,三维直角坐标系下非均匀网格离散的抛物方程模型,在保持计算精度、提高计算速度的同时,可以实现远距离声场的快速预报。另外,针对远距离局部海底地形与距离有关的三维声传播问题,给出了声场快速计算方法;在海底保持水平的区域,采用经典Kraken模型,重构抛物方程算法的初始场,随后依次递推求解地形与距离有关海底下的三维声场。采用改进模型,证明了远距离楔形波导声强增强效应。      

A single-scattering correction for the seismo-acoustic parabolic equation

An efficient single-scattering correction that does not require iterations is derived and tested for the seismo-acoustic parabolic equation. The approach is applicable to problems involving gradual range dependence in a waveguide with fluid and solid layers, including the key case of a sloping fluid-solid interface. The single-scattering correction is asymptotically equivalent to a special case of a single-scattering correction for problems that only have solid layers [Ku?sel et al., J. Acoust. Soc. Am. 121, 808-813 (2007)]. The single-scattering correction has a simple interpretation (conservation of interface conditions in an average sense) that facilitated its generalization to problems involving fluid layers. Promising results are obtained for problems in which the ocean bottom interface has a small slope.     

Researches on the distribution law of vector sound field in elastic wedge bottom

The method based on elastic parabolic equation method for calculating the sound vector field has been studied.The vector field in water and corresponding seismic wave field had been calculated for infra-sound in oceanic environment with elastic wedge bottom.The effects on sound field distribution for different frequency and depth of sound source had been researched,result shows that there is sound energy leakage into the bottom,the position where leakage occurred can be determined by the ratio of the ocean depth to the wavelength,as compared with normal mode theory.     

A wide angle and high Mach number parabolic equation

Various parabolic equations for advected acoustic waves have been derived based on the assumptions of small Mach number and narrow propagation angles, which are of limited validity in atmospheric acoustics. A parabolic equation solution that does not require these assumptions is derived in the weak shear limit, which is appropriate for frequencies of about 0.1 Hz and above for atmospheric acoustics. When the variables are scaled appropriately in this limit, terms involving derivatives of the sound speed, density, and wind speed are small but can have significant cumulative effects. To obtain a solution that is valid at large distances from the source, it is necessary to account for linear terms in the first derivatives of these quantities [A. D. Pierce, J. Acoust. Soc. Am. 87, 2292-2299 (1990)]. This approach is used to obtain a scalar wave equation for advected waves. Since this equation contains two depth operators that do not commute with each other, it does not readily factor into outgoing and incoming solutions. An approximate factorization is obtained that is correct to first order in the commutator of the depth operators.     

Multiscale finite-volume method for parabolic problems arising from compressible multiphase flow in porous media

The multiscale finite-volume (MSFV) method was originally developed for the solution of heterogeneous elliptic problems with reduced computational cost. Recently, some extensions of this method for parabolic problems have been proposed. These extensions proved effective for many cases, however, they are neither general nor completely satisfactory. For instance, they are not suitable for correctly capturing the transient behavior described by the parabolic pressure equation. In this paper, we present a general multiscale finite-volume method for parabolic problems arising, for example, from compressible multiphase flow in porous media. Opposed to previous methods, here, the basis and correction functions are solutions of full parabolic governing equations in localized domains. At the same time, to enhance the computational efficiency of the scheme, the basis functions are kept pressure independent and do not have to be recalculated as pressure evolves. This general approach requires no additional assumptions and its good efficiency and high accuracy is demonstrated for various challenging test cases. Finally, to improve the quality of the results and also to extend the scheme for highly anisotropic heterogeneous problems, it is combined with the iterative MSFV (i-MSFV) method for parabolic problems. As one iterates, the i-MSFV solutions of compressible multiphase problems (parabolic problems) converge to the corresponding fine-scale reference solutions in the same way as demonstrated recently for incompressible cases (elliptic problems). Therefore, the proposed MSFV method can also be regarded as an efficient linear solver for parabolic problems and studies of its efficiency are presented for many test cases.     

Disorder driven roughening transitions of elastic manifolds and periodic elastic media

The simultaneous effect of both disorder and crystal-lattice pinning on the equilibrium behavior of oriented elastic objects is studied using scaling arguments and a functional renormalization group technique. Our analysis applies to elastic manifolds, e.g., interfaces, as well as to periodic elastic media, e.g., charge-density waves or flux-line lattices. The competition between both pinning mechanisms leads to a continuous, disorder driven roughening transition between a flat state where the mean relative displacement saturates on large scales and a rough state with diverging relative displacement. The transition can be approached by changing the impurity concentration or, indirectly, by tuning the temperature since the pinning strengths of the random and crystal potential have in general a different temperature dependence. For D dimensional elastic manifolds interacting with either random-field or random-bond disorder a transition exists for 2<D<4, and the critical exponents are obtained to lowest order in . At the transition, the manifolds show a superuniversal logarithmic roughness. Dipolar interactions render lattice effects relevant also in the physical case of D=2. For periodic elastic media, a roughening transition exists only if the ratio p of the periodicities of the medium and the crystal lattice exceeds the critical value . For p<p _c the medium is always flat. Critical exponents are calculated in a double expansion in and and fulfill the scaling relations of random field models. Received 28 August 1998     

Guided elastic waves in a pre-stressed compressible interlayer

The propagation of guided elastic waves in a pre-stressed elastic compressible layer embedded in a different compressible material is examined. The waves propagate parallel to the planar layer interfaces as a superposed dynamic stress state on the statically pre-stressed layer and host material. The underlying stress condition in the two materials is characterized by equibiaxial in-plane deformations with common principal axes of strain, one of the axes being perpendicular to the layering. Both materials have arbitrary strain energy functions. The dispersion equation is derived in explicit form. Analysis of the dispersion equation reveals the propagation characteristics and their dependence on frequency, material parameters and stress parameters. Combinations of these parameters are also defined for which guided waves cannot propagate.     

Dyadic Green's functions for the perfectly conducting parabolic cylinder

Solutions of the scalar wave equation for parabolic cylinder coordinate system are discussed here. Dyadic Green's functions of the magnetic type for free space and for a perfectly conducting parabolic cylinder are developed. These functions are of fundamental importance for the solution of electromagnetic problems developed in the parabolic cylinder coordinate system, particularly those in the presence of perfectly conducting parabolic cylinders such as that of a parabolic cylinder reflector.     

Elasticity and dislocations in quasicrystals with 18-fold symmetry

Elasticity problems of quasicrystals with 18-fold rotational symmetry are studied. Constitutive equations and governing equations are obtained. For static elasticity problems, the displacement vectors in two phason fields are expressed in terms of two pairs of associated harmonic functions or two analytic functions. For dynamic problems, the displacement vectors can be represented in terms of an auxiliary function satisfying a fourth-order partial differential equation. A general solution of phasons is given by the solution of two diffusion equations. Phason elastic fields induced by a dislocation in a quasicrystal with 18-fold symmetry are determined and exhibit an inverse singularity.     

Wave field localization in a prestressed functionally graded layer

Characteristic features of wave field formation caused by a surface source of harmonic vibration in a prestressed functionally graded layer are investigated. It is assumed that the elastic moduli and the density of the material vary with depth according to arbitrary laws. The initial material of the medium is represented by a model hyperelastic material with third-order elastic moduli. The boundary-value problem for a set of Lamè equations is reduced to a set of Cauchy problems with initial conditions, which is solved by the Runge–Kutta–Merson method modified to fit the specific problem under study. Considering shear vibrations of a functionally graded layer as an example, effects of the type of its inhomogeneity, variations in its properties, and nature of its initial stressed state on the displacement distribution in depth are investigated. Special attention is paid to characteristic features of displacement localization in a layer with an interface-type inclusion near critical frequencies. A direct relation between the inhomogeneous layer structure and the type of displacement localization in depth is demonstrated. It is found that the role of initial stresses and variations in material parameters considerably increases in the vicinities of critical frequencies.     

Range-dependent waveguide scattering model calibrated for bottom reverberation in a continental shelf environment

An analytic model is developed for scattering from random inhomogeneities in range-dependent ocean waveguides using the Rayleigh-Born approximation to Green's theorem. The expected scattered intensity depends on statistical moments of fractional changes in compressibility and density, which scatter as monopoles and dipoles, respectively, and the coherence volume of the inhomogeneities. The model is calibrated for ocean bottom scattering using data acquired by instantaneous wide-area ocean acoustic waveguide remote sensing (OAWRS) and geophysical surveys of the ONR Geoclutter Program. The scattering strength of the seafloor on the New Jersey shelf, a typical continental shelf environment, is found to depend on wave number k, medium coherence volume V(c), and seabed depth penetration factor F(p) following a 10 log(10)(F(p)V(c)k(4)) dependence. A computationally efficient numerical approach is developed to rapidly compute bottom reverberation over wide areas using the parabolic equation by exploiting correlation between monopole and dipole scattering terms and introducing seafloor depth penetration factors. An approach is also developed for distinguishing moving clutter from statistically stationary background reverberation by tracking temporal and spatial fluctuations in OAWRS intensity images.