Modeling acoustic waves in locally enhanced meshes with a staggered-grid finite difference approach

Finite difference is a well-suited technique for modeling acoustic wave propagation in heterogeneous media as well as for imaging and inversion. Typically, the method aims at solving a set of partial differential equations for the unknown pressure field by using a regularly spaced grid. Although finite differences can be fast and cheap to implement, the accuracy of the solution is always restricted by the computational resources. This is a fundamental key point to treat when dealing with large-scale problems. In this work, we present and test a method that uses a non-uniform distribution of grid points to improve on accuracy or to reduce the required computational resources. The applied grid is generated through a coordinate transformation. Differential geometry and generalized coordinates are used to handle and analyze the effect of using a non-uniform grid. Results obtained with the presented method show that the applied transformation as well as the number of points-per-wavelength influences the stability and dispersion in the solution. We exploit this observation to locally improve the accuracy of our simulations. The work presented in this paper allows us to conclude that differential geometry for finite differences can be used to reduce dispersion and hence improve the accuracy when modeling acoustic wave propagation in heterogeneous media. In addition, it can be used to avoid oversampling through the optimization of the number of grid nodes required to have an accurate solution or just honor to the boundaries.     

Floquet–Bloch decomposition for the computation of dispersion of two-dimensional periodic,damped mechanical systems

Floquet–Bloch theorem is widely applied for computing the dispersion properties of periodic structures, and for estimating their wave modes and group velocities. The theorem allows reducing computational costs through modeling of a representative cell, while providing a rigorous and well-posed spectral problem representing wave dispersion in undamped media. Most studies employ the Floquet–Bloch approach for the analysis of undamped systems, or for systems with simple damping models such as viscous or proportional damping. In this paper, an alternative formulation is proposed whereby wave heading and frequency are used to scan the k-space and estimate the dispersion properties. The considered approach lends itself to the analysis of periodic structures with complex damping configurations, resulting for example from active control schemes, the presence of damping materials, or the use of shunted piezoelectric patches. Examples on waveguides with various levels of damping illustrate the performance and the characteristics of the proposed approach, and provide insights into the properties of the obtained eigensolutions.     

Resonant waves in elastic structured media: Dynamic homogenisation versus Green’s functions

We address an important issue of dynamic homogenisation in vector elasticity for a doubly periodic mass-spring elastic lattice. The notion of logarithmically growing resonant waves is used in the analysis of star-shaped wave forms induced by an oscillating point force. We note that the dispersion surfaces for Floquet–Bloch waves in the elastic lattice may contain critical points of the saddle type. Based on the local quadratic approximations of a dispersion surface, where the radian frequency is considered as a function of wave vector components, we deduce properties of a transient asymptotic solution associated with the contribution of the point source to the wave form. The notion of local Green’s functions is used to describe localised wave forms corresponding to the resonant frequency. The special feature of the problem is that, at the same resonant frequency, the Taylor quadratic approximations for different groups of the critical points on the dispersion surfaces (and hence different Floquet–Bloch vectors) are different. Thus, it is shown that for the vector case of micro-structured elastic medium there is no uniformly defined dynamic homogenisation procedure for a given resonant frequency. Instead, the continuous approximation of the wave field can be obtained through the asymptotic analysis of the lattice Green’s functions, presented in this paper.     

应力波在变截面体中的传播特性

对应力波在变截面体中的传播特性进行了理论研究和数值分析。以杆中一维纵波波动理论和谐波分析法为基础，研究截面变化所导致的应力波的波形弥散和波幅变化。推导了与截面变化相关的应力波演化因子，并对由于截面变化所造成的几何弥散等二维效应进行了分析，同时计算了变截面体的几何特征参数和截面变化等因素影响应力波演化规律的特点。     

A second-order radiation boundary condition for the shallow water wave equations on two-dimensional unstructured finite element grids

A second-order radiation boundary condition (RBC) is derived for 2D shallow water problems posed in ‘wave equation’ form and is implemented within the Galerkin finite element framework. The RBC is derived by matching the dispersion relation for the interior wave equation with an approximate solution to the exterior problem for outgoing waves. The matching is correct to second order, accounting for curvature of the wave front and the geometry. Implementation is achieved by using the RBC as an evolution equation for the normal gradient on the boundary, coupled through the natural boundary integral of the Galerkin interior problem. The formulation is easily implemented on non-straight, unstructured meshes of simple elements. Test cases show fidelity to solutions obtained on extended meshes and improvement relative to simpler first-order RBCs.     

On Boussinesq's paradigm in nonlinear wave propagation

Boussinesq's original derivation of his celebrated equation for surface waves on a fluid layer opened up new horizons that were to yield the concept of the soliton. The present contribution concerns the set of Boussinesq-like equations under the general title of ‘Boussinesq's paradigm’. These are true bi-directional wave equations occurring in many physical instances and sharing analogous properties. The emphasis is placed: (i) on generalized Boussinesq systems that involve higher-order linear dispersion through either additional space derivatives or additional wave operators (so-called double-dispersion equations); and (ii) on the ‘mechanics’ of the most representative localized nonlinear wave solutions. Dissipative cases and two-dimensional generalizations are also considered. To cite this article: C.I. Christov et al., C. R. Mecanique 335 (2007).     

Three-dimensional Gouyon waves

Three-dimensional steady-state periodic waves in deep water with weak (of order of the wave amplitude) vorticity are considered. The solution describing the wave properties is constructed by the perturbation theory method in modified Lagrangian coordinates. The wave structure and dispersion properties are found correct to the square of the wave amplitude.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 125–130, July–August, 1996.     

Acoustic Waves of Various Geometry in Multi-Fraction Bubbly Liquids

The propagation of acoustic waves of various geometry in mixtures of a liquid and a disperse phase consisting of small bubbles which differ from one another by both the radii and the thermophysical properties is investigated. A systemof differential equations of motion of the mixture is written and the dispersion relation is derived. The dispersion curves are constructed and damping of the pressure pulses is compared for the plane, cylindrical, and spherical waves in the bubbly liquids considered. The theory is compared with the experimental data.     

Prediction on dispersion in elastoplastic unsaturated granular media

This study focuses on the propagation of the plane wave in the elastoplastic unsaturated granular media, and the wave equations and dispersion equations are derived for the media under the framework of Cosserat theory. Due to symmetry, five different wave modes are considered and predicted for the elastoplastic unsaturated granular media based on the Cosserat theory, including two longitudinal waves, one rotational longitudinal wave and two coupled transverse–rotational transverse waves. The correspondence is discussed between these Cosserat wave modes and the classical wave modes. Based on the dispersion equations, the dispersion behaviors are obtained for the five Cosserat wave modes. The results indicated that the different stress-strain stages,including the elastic, hardening and softening stages, have obvious effect on the dispersion behaviors of the Cosserat wave modes.     

Propagation of longitudinal elastic waves in composites with a random set of spherical inclusions (effective field approach)

The work is devoted to the problem of plane monochromatic longitudinal wave propagation through a homogeneous elastic medium with a random set of spherical inclusions. The effective field method and quasicrystalline approximation are used for the calculation of the phase velocity and attenuation factor of the mean (coherent) wave field in the composite. The hypotheses of the method reduce the diffraction problem for many inclusions to a diffraction problem for one inclusion and, finally, allow for the derivation of the dispersion equation for the wave vector of the mean wave field in the composite. This dispersion equation serves for all frequencies of the incident field, properties and volume concentrations of inclusions. The long and short wave asymptotics of the solution of the dispersion equation are found in closed analytical forms. Numerical solutions of this equation are constructed in a wide region of frequencies of the incident field that covers long, middle, and short wave regions of propagating waves. The phase velocities and attenuation factors of the mean wave field are calculated for various elastic properties, density, and volume concentrations of the inclusions. Comparisons of the predictions of the method with some experimental data are presented; possible errors of the method are indicated and discussed.     

非圆截面杆中的非线性扭转波及其行波解

研究了非圆截面杆中非线性扭转波的传播特性.由于非圆截面杆的扭转运动会伴随有横截面的翘曲,这种翘曲运动将引起扭转波的弥散.如果同时考虑有限扭转变形和翘曲弥散的共同作用,将会得到非线性扭转波的方程.在相平面上,对非线性扭转波动方程进行定性分析,结果表明,在一定条件下方程存在同宿轨道或异宿轨道,分别相应于方程的孤波解或冲击波解.本文利用Jacobi椭圆函数展开法,对该非线性方程进行求解,得到了非线性波动方程的三类准确周期解及相应的孤波解和冲击波解,讨论了这些解存在的必要条件.这些条件与定性分析的结果相一致.     

Wave dispersion in gradient elastic solids and structures: A unified treatment

Analytical wave propagation studies in gradient elastic solids and structures are presented. These solids and structures involve an infinite space, a simple axial bar, a Bernoulli–Euler flexural beam and a Kirchhoff flexural plate. In all cases wave dispersion is observed as a result of introducing microstructural effects into the classical elastic material behavior through a simple gradient elasticity theory involving both micro-elastic and micro-inertia characteristics. It is observed that the micro-elastic characteristics are not enough for resulting in realistic dispersion curves and that the micro-inertia characteristics are needed in addition for that purpose for all the cases of solids and structures considered here. It is further observed that there exist similarities between the shear and rotary inertia corrections in the governing equations of motion for bars, beams and plates and the additions of micro-elastic (gradient elastic) and micro-inertia terms in the classical elastic material behavior in order to have wave dispersion in the above structures.     

Analysis of wave propagation in piezoelectric coupled cylinder affected by transverse shear and rotary inertia

In this paper, we examine the wave propagation in a piezoelectric coupled cylindrical shell affected by the shear effect and rotary inertia. A complete mathematical analysis of wave propagation solution in this piezoelectric coupled cylindrical shell is provided. The dispersion characteristics are derived through the solving an eigenvalue problem. Results are validated by the classical solution of a metallic cylinder. Besides providing and discussing the dispersion curves for different wave modes, we also examine the piezoelectric effects on the dispersion curves. Further to the above investigation, comparison of dispersion solutions from different shell theories is also conducted. This work may serve as a benchmark for wave propagation in piezoelectric coupled cylindrical shells.     

Effective medium method in the problem of axial elastic shear wave propagation through fiber composites

The effective medium method (EMM) is applied to the solution of the problem of monochromatic elastic shear wave propagation through matrix composite materials reinforced with cylindrical unidirected fibers. The dispersion equations for the wave numbers of the mean wave field in such composites are derived using two different versions of the EMM. Asymptotic solutions of these equations in the long and short wave regions are found in closed analytical forms. Numerical solutions of the dispersion equations are constructed in a wide region of frequencies of the incident field that covers long, middle and short wave regions of the mean wave field. Velocities and attenuation factors of the mean wave fields in the composites obtained by different versions of the EMM are compared for various volume concentrations and properties of the inclusions. The main discrepancies in the predictions of different versions of the EMM are indicated, analyzed and discussed.     

The wave propagation and dynamic response of rectangular functionally graded material plates with completed clamped supports under impulse load

In this paper, the wave propagation and dynamic response of the rectangular FGM plates with completed clamped supports under impulse load are analyzed. The effective material properties of functionally graded materials (FGMs) for the plate are assumed to vary continuously through the plate thickness and be distributed according to a volume fraction power law along the plate thickness. Considering the effects of transverse shear deformation and rotary inertia, the governing equations of the wave propagation in the functionally graded plate are derived by using the Hamilton’s principle. A complete discussion of dispersion of the FGM plates is given. Using the dispersion relation and integral transforms, exact integral solutions for the FGM plates under impulse load are obtained. The influence of volume fraction distributions on wave propagation and dynamic response of the FGM plates is analyzed.     

Wave propagation in carbon nanotubes embedded in an elastic matrix

This paper reports the results of an investigation into the characteristics of wave propagation in carbon nanotubes embedded in an elastic matrix, based on an exact shell model. Each of the concentric tubes of multi-walled carbon nanotubes is considered as an individual elastic shell and coupled together through the van der Waals forces between two adjacent tubes. The matrix surrounding carbon nanotubes is described as a spring element defined by the Winkler model. The effects of rotatory inertia and elastic matrix on the wave velocity, the critical frequency, and the amplitude ratio between two adjacent tubes are described and discussed through numerical examples. The results obtained show that wave propagation in carbon nanotubes may appear in a critical frequency at which the wave velocity changes suddenly; the elastic matrix surrounding carbon nanotubes debases the critical frequency and the wave velocity, and changes the vibration modes between two adjacent tubes; the rotatory inertia based on an exact shell model obviously influences the wave velocity at some wave modes. Finally, a comparison of dispersion solutions from different shell models is given. The present work may serve as a useful reference for the application and the design of nano-electronic and nano-drive devices, nano-oscillators, and nano-sensors, in which carbon nanotubes act as basic elements.     

Waves on the interface between a nematic liquid crystal and an ideal isotropic fluid

The problem of surface wave propagation over the interface between a nematic liquid crystal and an ideal isotropic fluid is considered. For the nematic liquid crystal the Frank-Oseen model with an isotropic viscous stress tensor is used. Anisotropic surface tension is described by the Rapini model. In this formulation, for the problem of harmonic small-amplitude surface wave propagation, in the case of infinite depths of both phases, an analytical solution is obtained. The dispersion relation is derived and its properties are investigated.     

Rayleigh-type wave propagation through liquid layer over corrugated substrate

The propagation of the Rayleigh-type wave in a fluid layer overlying a corrugated substrate is studied. The corrugated substrate is considered as a fluid saturated poroelastic substrate and a quadratically heterogeneous isotropic elastic substrate in Case I and Case II, respectively. Closed form expressions of dispersion relation for Case I and Case II are obtained. The influence of corrugation, porosity, and heterogeneity on the phase velocity of Rayleigh-type wave, for both cases, is highlighted and demonstrated through numerical computation and graphical discussion. Neglecting corrugation at the common interface, expressions of phase velocity of the Rayleigh-type wave for both cases are derived in a closed form as a special case of the problem. Comparison between the presence and the absence of both heterogeneity and poroelasticity in the substrate of the composite structure is a key in the present study.     

Axi-symmetric wave propagation in a cylinder coated with a piezoelectric layer

This paper presents the wave propagation in a cylinder coated with a thin piezoelectric layer. The piezoelectric coupling effects are fully modeled in the mechanics model for this piezoelectric coupled cylindrical shell with bending resistance. The decoupled torsional wave velocity and the dispersion curves for the two- mode shell model are obtained theoretically. The cut-off frequency and phase velocities at limit wave number are also derived. The numerical simulations are conducted to present the results of wave propagation in this cylindrical shell and as well as to compare the results by the current bending theory and the membrane shell theory. From the comparisons, the results display obvious deference of wave propagations in terms of dispersion characteristics by different shell theories when thicker piezoelectric layer are used and when higher wave number is considered. The results of this paper can serve as a reference for future study on wave propagation in coupled structures as well as in the design of smart structures incorporating piezoelectric materials.     

Elastic waves in piezoceramic cylinders of sector cross-section

The propagation of elastic waves in piezoceramic cylindrical waveguides of circular cross-sections with sector cut is investigated on the basis of the linear theory of electroelasticity. Dispersion functions are obtained from boundary conditions in an analytical form of functional determinants for each value of the generalized wave number. A selected set of numerical results including real, imaginary and complex branches of full dispersion spectrums with various symmetry of wave movements is presented to describe the essential characteristics of the waves. Leading effects of spectrums transformation by change of waveguide’s angular measure are enlightened, and wave asymptotic behavior is analyzed. The variation of the cross-section is considered as a mechanism to control the dispersion characteristics of waveguides.