Wave scattering from a rough surface: Solution by an iterative method

Two iterative solutions of the Helmholtz equation for a scalar field in above a rough surface that admits the Dirichlet boundary condition are derived. The bases for the two iterative methods are two different boundary integral equations that represent the solution. The first integral equation is classified as a Fredholm integral equation of the first kind. The second is classified as a Fredholm integral equation of the second kind. This classification suggests that it is easier to find stable solution methods to the second equation. In both methods, the boundary integral was separated into a major part which is easy to calculate and a local residual part. The major part is a convolution and thus can be calculated using FFT in complexity O(N log N), where N is the number of surface points in which the surface height and its first derivatives together with the incoming wave and its normal derivative are all known. The residual element of the equations can be approximated efficiently only for surfaces where their amplitude is less than the wavelength of the incoming wave. The iterative schemes were tested numerically against a reference solution in order to examine the applicability range, the error estimation and the stability of the schemes. All tests supported the superiority of the second method. In particular the error estimation and stability tests indicated good performance for surfaces with slope up to 1. Yet, being an equation in the scattered field alone, makes the first method useful as a benchmark solution in its domain of applicability.     

Scattering of elastic waves by an arbitrary small imperfection in the surface of a half-space

T , the first of two articles, is concerned with the scattering of elastic waves by arbitrary surface-breaking or near surface defects in an isotropic half-plane. We present an analytical solution, by the method of matched asymptotic expansions, when the parameter , which is the ratio of a typical length scale of the imperfection to the incident radiation's wavelength, is small. The problem is formulated for a general class of small defects, including cracks, surface bumps and inclusions, and for arbitrary incident waves. As a straightforward example of the asymptotic scheme we specialize the defect to a two-dimensional circular void or protrusion, which breaks the free surface, and assume Rayleigh wave excitation ; this inner problem is exactly solvable by conformal mapping methods. The displacement field is found uniformly to leading order in , and the Rayleigh waves which are scattered by the crack are explicitly determined. In the second article we use the method given here to tackle the important problem of an inclined edge-crack. In that work we show that the scattered field can be found to any asymptotic order in a straightforward manner, and in particular the Rayleigh wave coefficients are given to O(²).     

Focusing of an ultrasonic beam by a curved interface

We calculate the wavefields near the caustics, and their cusps, formed when a well collimated, ultrasonic beam scatters from a concave fluid-solid interface. The radius of curvature of the interface is assumed to be sufficiently large and the angle of incidence sufficiently small that only reflection and transmission of the beam need be considered. The incident beam is modeled as a two-dimensional wavefield whose initial profile is rectangular. The aperture is assumed to be sufficiently large (in wavelengths) that a well collimated beam is radiated, and the interface is assumed to lie close enough to the aperture that it is struck by a wavefield that has not yet evolved into a cylindrical wave. The scattered wavefields are represented as multiple integrals and are evaluated using a combination of asymptotic and numerical analysis. Special attention is given to the sometimes competing effects of the shadow boundary of each scattered beam with its corresponding caustic, and cusp if one is formed.     

A lattice Boltzmann method for viscous free surface waves in two dimensions

We propose a new method based on the combination of the lattice Boltzmann equation (LBE) and the kinematic boundary condition (KBC) method to simulate viscous free surface wave in two dimensions. In our method, the flow field is modeled by LBE, whereas the free surface is explicitly tracked by the local height function, which is calculated by the KBC method. The free surface boundary condition (FSBC) for LBE is revised from previous researches. Interpolation‐supplemented lattice Boltzmann (ISLB) method is introduced, which enables our approach to be applied on arbitrary, nonuniform mesh grids. Five cases are simulated respectively to validate the LBE–KBC method: the stationary flow and the solitary waves simulated by the revised‐FSBC are more accurate than the one obtained by the former‐FSBC; numerical results of standing waves show that our method is compatible to the existing two‐dimensional finite‐volume scheme; cases of small amplitude Stokes wave and waves traveling over a submerged bar show good agreement on wave celerity, wavelength, wave amplitude and wave period between numerical results and corresponding analytical solutions and/or experiment data.Copyright © 2012 John Wiley & Sons, Ltd.     

Anti-plane waves near an interface between two piezoelectric half-spaces

We show the existence of certain new waves that can propagate near an interface between two half-spaces of different piezoelectric ceramics, where the interface is modeled by a membrane with the surface/interface elasticity [1]. The current configuration can be reduced into a number of well-known results as special cases, such as Love wave, Bleustein and Gulyaev wave. Together with our previous work for the imperfect interface [2], a full range of consideration of the interface affecting the anti-plane waves is now completed.     

Scattering of electromagnetic waves by an acoustic disturbance in the atmosphere

Summary The variations in density associated with an acoustic wave are shown to influence the propagation of an electromagnetic wave. The variation in electric permittivity of the atmosphere caused by an acoustic wave is expressed by the power density of the wave. Under the assumption that the acoustic waves are spherical, the scattered field at large distances is given in terms of definite integrals which are evaluated by means of the method of Stationary Phase. The limiting case of plane acoustic waves is discussed, and two numerical examples are given.     

Acoustomechanical constitutive theory for soft materials

Acoustic wave propagation from surrounding medium into a soft material can generate acoustic radiation stress due to acoustic momentum transfer inside the medium and material, as well as at the interface between the two. To analyze acoustic-induced deformation of soft materials, we establish an acoustomechanical constitutive theory by com-bining the acoustic radiation stress theory and the nonlinear elasticity theory for soft materials. The acoustic radiation stress tensor is formulated by time averaging the momen-tum equation of particle motion, which is then introduced into the nonlinear elasticity constitutive relation to construct the acoustomechanical constitutive theory for soft materials. Considering a specified case of soft material sheet subjected to two counter-propagating acoustic waves, we demonstrate the nonlinear large deformation of the soft material and ana-lyze the interaction between acoustic waves and material deformation under the conditions of total reflection, acoustic transparency, and acoustic mismatch.     

Numerical simulation of flow‐induced noise using LES/SAS and Lighthill's acoustic analogy

We present a finite element (FE) formulation of Lighthill's acoustic analogy for the hybrid computation of noise generated by turbulent flows. In the present approach, the flow field is computed using large eddy simulation and scale adaptive simulation turbulence models. The acoustic propagation is obtained by solving the variational formulation of Lighthill's acoustic analogy with the FE method. In order to preserve the acoustic energy, we compute the inhomogeneous part of Lighthill's wave equation by applying the FE formulation on the fine flow grid. The resulting acoustic nodal loads are then conservatively interpolated to the coarser acoustic grid. Subsequently, the radiated acoustic field can be solved in both time and frequency domains. In the latter case, an enhanced perfectly matched layer technique is employed, allowing one to truncate the computational domain in the acoustic near field, without compromising the numerical solution. Our hybrid approach is validated by comparing the numerical results of the acoustic field induced by a corotating vortex pair with the corresponding analytical solution. To demonstrate the applicability of our scheme, we present full 3D numerical results for the computed acoustic field generated by the turbulent flow around square cylinder geometries. The sound pressure levels obtained compare well with measured values. Copyright © 2009 John Wiley & Sons, Ltd.     

The effect of an insoluble surfactant on the skin friction of a bubble

In this paper, we develop a novel moving mesh method suitable for solving axisymmetric free-boundary problems, including the Marangoni effect induced by surfactant or temperature variation. This method employs a body-fitted grid system where the gas–liquid interface is one line of the grid system. We model the surfactant equation of state with a non-linear Langmuir law, and, for simplicity, we limit ourselves to the situation of an insoluble surfactant. We solve complicated dynamic boundary conditions accurately on the gas–liquid interface in the framework of finite-volume methods. Our method is used to study the effect of a surfactant on the skin friction of a bubble in a uniaxial flow. For the limiting case where the surface diffusivity is zero, the effect of a tangential stress generated by the surface tension gradient, allows us to explain a new phenomenon in high concentration regimes: larger surface tension, but also larger deformation. Furthermore, this condition leads to the formation of boundary layers and flow separation at high Reynolds numbers. The influence of these complex flow patterns is examined.     

An integrated smoothed particle hydrodynamics model for complex interfacial flows with large density ratios

In this paper, an integrated smoothed particle hydrodynamics (SPH) model for complex interfacial flows with large density ratios is developed. The discrete continuity equation and acceleration equation are obtained by considering the time derivative of the volume of particle and Eckart's continuum Lagrangian equation. A continuum surface force model is used to meet the fact that surface force may not be distributed uniformly on each side of the interface. An improved boundary condition is imposed to model wall free-slip and no-slip condition for interfacial flows with large density ratios. Particle shifting algorithm (PSA) is added for interfacial flows by imposing the normal correction near the interface, called as Interface-PSA. Then four representative numerical examples, including droplet deformation, Rayleigh-Taylor instability, dam breaking, and bubble rising, are presented and compared well with reference data. It is demonstrated that inherent interfacial flow physics can be well captured, including surface tension and the dynamic evolution of the complex interfaces.     

Scattered fracture of porous materials with brittle skeleton

A model of damage accumulation in a porous medium with a brittle skeleton saturated with a compressible fluid is formulated in the isothermal approximation. The model takes account of the skeleton elastic energy transformation into the surface energy of microcracks. In the case of arbitrary deformations of an anisotropic material, constitutive equations are obtained in a general form that is necessary and sufficient for the objectivity and thermodynamic consistency principles to be satisfied. We also formulate the kinetics equation ensuring that the scattered fracture dissipation is nonnegative for any loading history. For small deviations from the initial state, we propose an elastic potential which permits describing the principal characteristics of the behavior of a saturated porous medium with a brittle skeleton. We study the acoustic properties of the material under study and find their relationship with the strength criterion depending on the accumulated damage and the material current deformation. We consider the problem of scattered fracture of a saturated porous material in a neighborhood of a spherical cavity. We show that the cavity failure occurs if the Hadamard condition is violated.     

流固耦合地震波动问题的显式谱元模拟方法

流固耦合地震波动问题主要研究由流体和固体构成的复杂系统中地震波传播特性及其规律. 传统模拟方法中一般以声波方程、弹性波方程的数值解分别描述理想流体和弹性固体中的波动, 并实时地处理两种不同性质介质之间的相互耦合作用, 数值格式复杂且限制数值模拟精度与计算效率. 本文采用谱元法结合多次透射公式人工边界条件实现了一种流固耦合地震波动问题的高阶显式数值计算方法. 该方法利用了流固耦合问题统一计算框架,可将饱和多孔介质的Biot波动方程分别退化为理想流体的声波方程和弹性固体的弹性波方程. 通过P波垂直入射的水平成层理想流体-饱和多孔介质-弹性固体场地模型、P波斜入射的不规则层状界面以及任意形状界面的理想流体-饱和多孔介质-弹性固体场地模型等三个算例, 与传递函数法解析解以及集中质量有限元法计算结果进行对比分析, 证明了本文方法的正确性与有效性. 数值模拟结果表明, 本文方法相较传统有限元法可以少得多的节点数量获得更高的数值精度, 并且在较宽的频率范围内都能可靠地模拟出流固耦合系统的动力响应, 充分体现出本文方法兼顾高精度、计算效率和复杂场地建模灵活的特点.      

Surface waves in a piezoelectric–semiconductor composite structure

The present work deals with the propagation of interfacial surface waves in a composite consisting of homogeneous, transversely isotropic, piezoelectric halfspace underlying a thin layer of non-piezoelectric semiconductor material. The mathematical model of the problem is depicted by partial differential equations of motion for the structure and boundary conditions to be satisfied at the interface and free surface of the composite. After obtaining formal wave solution of the model the secular equation that governs the propagation of surface waves in the considered composite structure has been derived in compact form. The numerical solution of secular equation is being carried out for the composites Si–CdSe, Ge–CdSe and Ge–PZT by employing functional iteration method along with irreducible Cardano method using MATLAB programming. The computer simulated results in respect of dispersion curves, attenuation coefficient and specific loss factor of energy dissipation are presented graphically for Si–CdSe composite to illustrate the analytical developments. We have extended our analysis to Ge–CdSe and Ge–PZT composites also. However, to avoid clustering of profiles and also to have clear understanding of the variations, the computer simulated values of phase velocity and attenuation coefficient are presented in tabular form for all three considered composite structures. This work may be useful for designing and construction of surface acoustic wave (SAW) devices and electronics industry.     

Bandgap calculation for mixed in-plane waves in 2D phononic crystals based on Dirichlet-to-Neumann map

In this paper, a method based on the Dirichlet-to-Neumann map is developed for bandgap calculation of mixed in-plane waves propagating in 2D phononic crystals with square and triangular lattices. The method expresses the scattered fields in a unit cell as the cylindrical wave expansions and imposes the Bloch condition on the boundary of the unit cell. The Dirichlet-to-Neumann (DtN) map is applied to obtain a linear eigenvalue equation, from which the Bloch wave vectors along the irreducible Brillouin zone are calculated for a given frequency. Compared with other methods, the present method is memory-saving and time-saving. It can yield accurate results with fast convergence for various material combinations including those with large acous- tic mismatch without extra computational cost. The method is also efficient for mixed fluid-solid systems because it considers the different wave modes in the fluid and solid as well as the proper fluid-solid interface condition.     

Sound scattering and its cancellation by an elastic spherical shell in free space and near a free surface

Sound scattering by an elastic spherical shell is analysed using linear acoustics and linear structural dynamics. It is suggested to utilize the shell’s structural dynamics to reduce or even eliminate the scattered sound field, thus making it practically acoustically invisible. This can be achieved using a prescribed external pressure distribution acting on the shell’s wall. Exact analytical solutions are found for that external pressure distribution, eliminating the scattered wave when the sphere is in free space or near a free surface and is subject to an incoming planar monochromatic sound wave. The latter is assumed to propagate in a direction perpendicular to the free surface (if it exists). The case of a few pressure-actuators acting on the shell’s wall is also modelled and an optimal solution which reduces the sound scattering by these actuators is found. An aluminium shell of 1 m radius and 5 mm thickness, situated in fresh water is analysed for sound frequencies of up to 10 kHz. The scattered wave fields are presented as well as the external pressure distributions that eliminate these scattered sound field, i.e. achieving acoustic cloaking. Significant reduction in the scattered wave energy and the target strength of more than 10 dB are also realized using a few pressure-actuators as long as the distance between the actuators is no more than three times the incident wave length for the investigated cases.     

水下声波传播过程的时域间断Galerkin有限元方法

本文构建了声压波动方程的改进时域间断Galerkin有限元方法.传统时域连续有限元方法在计算高梯度、强间断特征水中声波传播问题时往往会出现虚假数值振荡现象,这些数值振荡会影响正常波动的计算精度.为了解决这一问题,本文通过引入人工阻尼的方式构建了改进的时域间断Galerkin有限元方法,并针对具有高梯度、强间断特征的多障...     

The low frequency scatter of SH waves by a rough interface in an elastic plate

The study of the reflection and transmission of low frequency SH waves incident upon a rough interface in an elastic plate is undertaken by employing a theory of acoustic wave scattering from rough surfaces originally due to Biot and subsequently generalised to the case of elastic media. In this theory the interface is replaced by a distribution of voids/asperities whose individual size is small compared to the excitation wavelength. We plot the absolute values of the reflection and transmission coefficients versus frequency when a single symmetric SH plate mode is used as the input excitation. The different types of inclusions are used to simulate the rough surface are the hollow, fluid filled and aluminum spheres. Lastly, the loss of energy due to scattering is also estimated for the different inclusion distributions considered.     

Bleustein–Gulyaev waves in 6 mm piezoelectric materials loaded with a viscous liquid layer of finite thickness

In this paper, we analyze the propagation of Bleustein–Gulyaev waves in an unbounded piezoelectric half-space loaded with a viscous liquid layer of finite thickness within the linear elastic theories. Exact solutions of the phase velocity equations are obtained in the cases of both electrically open circuit and short circuit by solving the equilibrium equations of piezoelectric materials and the diffusion equation of viscous liquid. A PZT-5H/Glycerin system is selected to perform the numerical calculation. The results show that the mass density and the viscous coefficient have different effects on the propagation attenuation and phase velocity under different electrical boundary conditions. In particular, the penetration depth of the waves is of the same order as the wavelength in the case of electrically short circuit. These effects can be used to manipulate the behavior of the waves and have implications in the application of acoustic wave devices.     

界面脱胶圆夹杂对SH波散射的远场解

采用Green函数方法和复变函数法研究了SH波对界面脱胶圆夹杂的散射问题,并给出 了远场解答. 首先,沿双质材料界面将整个空间分成上下两部分, 在下半空间,给出了在水 平表面上任意一点承受时间谐和的出平面线源载荷作用时的位移函数,取该位移函数作为 Green函数. 其次,在下半空间,利用相关文献给出的Green函数,在上下空间连接时在双质材料界面处满足 连续性条件,构造出半圆形脱胶裂纹,进而求出应力和位移的表达式,建立积分方 程组, 给出了散射波远场位移模式和散射截面的解答, 分析了在不同参数变化时SH 波散射的远场特性. 结果表明,脱胶结构的存在对位移和散射截面有较大的放大作用.     

The Sub-Harmonic Bifurcation of Stokes Waves

Steady periodic water waves on the free surface of an infinitely deep irrotational flow under gravity without surface tension (Stokes waves) can be described in terms of solutions of a quasi-linear equation which involves the Hilbert transform and which is the Euler-Lagrange equation of a simple functional. The unknowns are a 2π-periodic function w which gives the wave profile and the Froude number, a dimensionless parameter reflecting the wavelength when the wave speed is fixed (and vice versa). Although this equation is exact, it is quadratic (with no higher order terms) and the global structure of its solution set can be studied using elements of the theory of real analytic varieties and variational techniques. In this paper it is shown that there bifurcates from the first eigenvalue of the linearised problem a uniquely defined arc-wise connected set of solutions with prescribed minimal period which, although it is not necessarily maximal as a connected set of solutions and may possibly self-intersect, has a local real analytic parametrisation and contains a wave of greatest height in its closure (suitably defined). Moreover it contains infinitely many points which are either turning points or points where solutions with the prescribed minimal period bifurcate. (The numerical evidence is that only the former occurs, and this remains an open question.) It is also shown that there are infinitely many values of the Froude number at which Stokes waves, having a minimal wavelength that is an arbitrarily large integer multiple of the basic wavelength, bifurcate from the primary branch. These are the sub-harmonic bifurcations in the paper's title. (In 1925 Levi-Civita speculated that the minimal wavelength of a Stokes wave propagating with speed c did not exceed 2πc ²/g. This is disproved by our result on sub-harmonic bifurcation, since it shows that there are Stokes waves with bounded propagation speeds but arbitrarily large minimal wavelengths.) Although the work of Benjamin & Feir} and others [9, 10] has shown Stokes waves on deep water to be unstable, they retain a central place in theoretical hydrodynamics. The mathematical tools used to study them here are real analytic-function theory, spectral theory of periodic linear pseudo-differential operators and Morse theory, all combined with the deep influence of a paper by Plotnikov [36]. Accepted: December 6, 1999