Rayleigh wave propagation in curved waveguides

A JWKB asymptotic expansion describing inplane elastic waves is used to approximate a Rayleigh-like wave guided within a curved elastic waveguide whose curvature is small and changes slowly over a wavelength. The two lowest eigenmodes in a curved guide, taken together, constitute the Rayleigh-like wave. It is shown that this wave lies in the shadows of four, closely spaced, virtual caustics, two caustics per constituent eigenmode. If the curvature becomes too large one or more of the caustics ceases to be virtual and enters the guide after which a Rayleigh-like wave cannot propagate. The overall disturbance is shown to have an amplitude that is modulated because the wavenumbers of the constituent eigenmodes differ by a small amount. Moreover, the disturbance is shown to propagate with a wavenumber that, to leading order, has a linear dependence on the curvature causing the phase to be modulated, as well. Passing from a thin guide to a very thick one suppresses the amplitude modulation, making the phase modulation evident. Propagation into an environment of increasing curvature, for both thin and thick, shallowly curved guides is studied so that the modulations may be observed.     

Theoretical and numerical investigation of HF elastic wave propagation in two-dimensional periodic beam lattices

The elastic wave propagation phenomena in two-dimensional periodic beam lattices are studied by using the Bloch wave transform. The numerical modeling is applied to the hexagonal and the rectangular beam lattices, in which, both the in-plane （with respect to the lattice plane） and out-of-plane waves are considered. The dispersion relations are obtained by calculating the Bloch eigenfrequencies and eigenmodes. The frequency bandgaps are observed and the influence of the elastic and geometric properties of the primitive cell on the bandgaps is studied. By analyzing the phase and the group velocities of the Bloch wave modes, the anisotropic behaviors and the dispersive characteristics of the hexagonal beam lattice with respect to the wave prop- agation are highlighted in high frequency domains. One im- portant result presented herein is the comparison between the first Bloch wave modes to the membrane and bend- ing/transverse shear wave modes of the classical equivalent homogenized orthotropic plate model of the hexagonal beam lattice. It is shown that, in low frequency ranges, the homog- enized plate model can correctly represent both the in-plane and out-of-plane dynamic behaviors of the beam lattice, its frequency validity domain can be precisely evaluated thanks to the Bloch modal analysis. As another important and original result, we have highlighted the existence of the retro- propagating Bloch wave modes with a negative group veloc- ity, and of the corresponding ＂retro-propagating＂ frequency bands.     

Fast Fourier homogenization for elastic wave propagation in complex media

In the context of acoustic or elastic wave propagation, the non-periodic asymptotic homogenization method allows one to determine a smooth effective medium and equations associated with the wave propagation in a given complex elastic or acoustic medium down to a given minimum wavelength. By smoothing all discontinuities and fine scales of the original medium, the homogenization technique considerably reduces meshing difficulties as well as the numerical cost associated with the wave equation solver, while producing the same waveform as for the original medium (up to the desired accuracy). Nevertheless, finding the effective medium requires one to solve the so-called “cell problem”, which corresponds to an elasto-static equation with a finite set of distinct loadings. For general elastic or acoustic media, the cell problem is a large problem that has to be solved on the whole domain and its resolution implies the use of a finite element solver and a mesh of the fine scale medium. Even if solving the cell problem is simpler than solving the wave equation in the original medium (because it is time and source independent, based on simple tetrahedral meshes and embarrassingly parallel) it is still a challenge. In this work, we present an alternative method to the finite element approach for solving the cell problem. It is based on a well-known method designed by H. Moulinec and P. Suquet in 1998 in structural mechanics. This iterative technique relies on Green functions of a simple reference medium and extensively uses Fast Fourier Transforms. It is easy to implement, very efficient and relies on a simple regular gridding of the medium. Through examples we show that the method gives excellent results, even, under some conditions, for discontinuous media.     

Quasi-P-waves at a corrugated interface between two dissimilar monoclinic elastic half-spaces

In this paper we have derived reflection and transmission coefficients of qP-waves at a corrugated interface between two different elastic half-spaces of monoclinic type. Using Rayleigh’s method, the expressions for reflection and transmission coefficients are derived in closed form for a specific interface and for the first order approximation of the corrugation. Numerical computations are performed for a specific model and the results obtained have been shown graphically. The variation of the modulus of reflection and transmission coefficients with the angle of incidence, frequency and corrugation of the interface are shown separately. These coefficients are found to be strongly influenced by the angle of incidence, frequency, elastic parameters and amplitude of the corrugation of the interface. It is found that (i) the modulus of reflection and transmission coefficients at the plane interface are independent of corrugation of the interface and that of frequency of the incident wave, (ii) the reflection and transmission coefficients of regularly reflected and transmitted waves are found to be greater than that of irregularly reflected and transmitted waves, (iii) the coefficients of irregularly reflected and transmitted waves are found to increase and decrease with increase of corrugation and frequency parameters respectively. The results of Singh and Khurana [Singh, S.J., Khurana, S., 2001. Reflection and transmission of P- and SV-waves at the interface two between monoclinic elastic half-spaces. Proc. Natl. Acad. Sci. India 71(A) (IV), 305–319] have been reduced from the present problem.     

Long wavelength inner-resonance cut-off frequencies in elastic composite materials

We revisit an ancient paper (Auriault and Bonnet, 1985) which points out the existence of cut-off frequencies for long acoustic wavelength in high-contrast elastic composite materials, i.e. when the wavelength is large with respect to the characteristic heterogeneity length. The separation of scales enables the use of the method of multiple scale expansions for periodic structures, a powerful upscaling technique from the heterogeneity scale to the wavelength scale. However, the results remain valid for non-periodic composite materials which show a Representative Elementary Volume (REV). The paper extends the previous investigations to three-component composite materials made of hard inclusions, coated with a soft material, both of arbitrary geometry, and embedded in a connected stiff material. The equivalent macroscopic models are rigorously established as well as their domains of validity. Provided that the stiffness contrast within the soft and the connected stiff materials is of the order of the squared separation of scales parameter, it is demonstrated (i) that the propagation of long wave may coincide with the resonance frequencies of the hard inclusions/soft material system and (ii) that the macroscopic model presents a series of cut-off frequencies given by an eigenvalue problem for the resonating domain in the cell. These results are illustrated in the case of stratified composites and the possible microstructures of heterogeneous media in which the inner dynamics phenomena may occur are discussed.     

Inelastic deformation and energy dissipation in ceramics: A mechanism-based constitutive model

A mechanism-based constitutive model is presented for the inelastic deformation and fracture of ceramics. The model comprises four essential features: (i) micro-crack extension rates based on stress-intensity calculations and a crack growth law, (ii) the effect of the crack density on the stiffness, inclusive of crack closure, (iii) plasticity at high confining pressures, and (iv) initial flaws that scale with the grain size. Predictions of stress/strain responses for a range of stress states demonstrate that the model captures the transition from deformation by micro-cracking at low triaxiality to plastic slip at high triaxialities. Moreover, natural outcomes of the model include dilation (or bulking) upon micro-cracking, as well as the increase in the shear strength of the damaged ceramic with increasing triaxiality. Cavity expansion calculations are used to extract some key physics relevant to penetration. Three domains have been identified: (i) quasi-static, where the ceramic fails due to the outward propagation of a compression damage front, (ii) intermediate velocity, where an outward propagating compression damage front is accompanied by an inward propagating tensile (or spallation) front caused by the reflection of the elastic wave from the outer surface and (iii) high velocity, wherein plastic deformation initiates at the inner surface of the shell followed by spalling within a tensile damage front when the elastic wave reflects from the outer surface. Consistent with experimental observations, the cavity pressure is sensitive to the grain size under quasi-static conditions but relatively insensitive under dynamic loadings.     

Nonlocal shell model for elastic wave propagation in single- and double-walled carbon nanotubes

This paper investigates the transverse and torsional wave in single- and double-walled carbon nanotubes (SWCNTs and DWCNTs), focusing on the effect of carbon nanotube microstructure on wave dispersion. The SWCNTs and DWCNTs are modeled as nonlocal single and double elastic cylindrical shells. Molecular dynamics (MD) simulations indicate that the wave dispersion predicted by the nonlocal elastic cylindrical shell theory shows good agreement with that of the MD simulations in a wide frequency range up to the terahertz region. The nonlocal elastic shell theory provides a better prediction of the dispersion relationships than the classical shell theory when the wavenumber is large enough for the carbon nanotube microstructure to have a significant influence on the wave dispersion. The nonlocal shell models are required when the wavelengths are approximately less than 2.36×10⁻⁹ and 0.95×10⁻⁹ m for transverse wave in armchair (15,15) SWCNT and torsional wave in armchair (10,10) SWCNT, respectively. Moreover, an MD-based estimation of the scale coefficient e₀ for the nonlocal elastic cylindrical shell model is suggested. Due to the small-scale effects of SWCNTs and the interlayer van der Waals interaction of DWCNTs, the phase difference of the transverse wave in the inner and outer tube can be observed in MD simulations in wave propagation at high frequency. However, the van der Waals interaction has little effect on the phase difference of transverse wave.     

Effect of Thermal Modulation on the Onset of Convection in a Porous Medium Layer Saturated by a Nanofluid

The effect of time-periodic temperature modulation at the onset of convection in a Boussinesq porous medium saturated by a nanofluid is studied analytically. The model used for the nanofluid incorporates the effects of Brownian motion. Three types of boundary temperature modulations are considered namely, symmetric, asymmetric, and only the lower wall temperature is modulated while the upper wall is held at constant temperature. The perturbation method is applied for computing the critical Rayleigh and wave numbers for small amplitude temperature modulation. The shift in the critical Rayleigh number is calculated as a function of frequency of modulation, concentration Rayleigh number, porosity, Lewis number, and thermal capacity ratio. It has been shown that it is possible to advance or delay the onset of convection by time-periodic modulation of the wall temperature. The nanofluid is found to have more stabilizing effect when compared to regular fluid. Low frequency is destabilizing, while high frequency is always stabilizing for symmetric modulation. Asymmetric modulation and only lower wall temperature modulation is stabilizing for all frequencies when concentration Rayleigh number is greater than one.     

Transverse wave at a plane interface between isotropic elastic and unsaturated porous elastic solid half-spaces

Based on the poroelasticity theory, this article investigates the reflection and transmission characteristics of an incident plane transverse wave at a plane interface between an isotropic elastic half-space and an unsaturated poroelastic solid half-space. For this purpose, the effect of the saturation degree and frequency on the properties of the four bulk waves in unsaturated porous medium, i.e., three longitudinal waves and one transverse wave, are discussed at first. Two general cases of mode conversion are considered: (i) The initial transverse wave is incident from an unsaturated poroelastic half-space to the interface, and (ii) the initial transverse wave is incident from an elastic solid half-space to the interface. The expressions for the partition of energy at the interface during transmission and reflection process of waves are presented in explicit forms. At last, numerical computations are performed for these two cases and the results obtained are depicted, respectively. The variation of the amplitude ratios and energy ratios with the saturation degree and incident angle is illustrated in detail. It is also verified that, at the interface, the sum of energy ratios is approximately equal to unity as expected.     

Inertial waves in a rotating annulus with inclined inner cylinder: comparing the spectrum of wave attractor frequency bands and the eigenspectrum in the limit of zero inclination

We investigate theoretically inertial waves inside a liquid confined between two co-rotating coaxial cylinders of finite length. We consider the case of small viscosity and high angular velocity (i.e., small Ekman numbers), a parameter range of interest for many geophysical applications. In this case, inertial waves propagating in the container show multiple reflections at the walls before the waves can be damped by weak diffusion. We allow for the inner cylinder wall to be parallel or inclined with respect to the annulus’ vector of rotation (truncated cone). For the limit of zero viscosity, the wave propagation is governed by a boundary value problem that is composed of a linear second-order hyperbolic partial differential equation and the impermeability boundary conditions. For the special case of vertical cylinder walls (no inclination of the inner cylinder), this boundary value problem is separable, the corresponding eigenmodes can analytically be found and they are regular. However, when the inner cylinder wall is inclined, the hyperbolicity of the governing equation leads to internal shear layers (corresponding to singularities for the inviscid case). The geometrical structure of the shear layers can be explained by inertial waves, trapped on limit cycles denoted as wave attractors. The shape of the limit cycles depends on the wave frequency. In fact, the spectrum of regular modes, existing for the case of vertical cylinder walls, vanishes almost completely when the inner wall is inclined. Instead of a spectrum of discrete frequencies and regular eigenmodes, a spectrum of wave attractor frequency bands and singular eigenmodes exist. The question addressed here is whether the spectrum of wave attractor intervals collapses to the discrete frequency spectrum when the inclination angle of the inner cylinder goes to zero. To answer this question, the attractor frequency intervals are evaluated numerically for a series of decreasing cylinder inclination angles and are compared with the analytically found eigenspectrum for the case of zero inclination. Goal is to better understand the asymptotic behavior of the problem for decreasing inclination angles. This understanding helps to interpret results from laboratory experiments with geometries that differ from the perfect annulus with parallel cylinder walls.     

Effect of Thermal/Gravity Modulation on the Onset of Convection in a Maxwell Fluid Saturated Porous Layer

The effect of thermal/gravity modulation on the onset of convection in a Maxwell fluid saturated porous layer is investigated by a linear stability analysis. Modified Darcy–Maxwell model is used to describe the fluid motion. The regular perturbation method based on the small amplitude of modulation is employed to compute the critical Rayleigh number and the corresponding wavenumber. The stability of the system characterized by a correction Rayleigh number is calculated as a function of the viscoelastic parameter, Darcy–Prandtl number, normalized porosity, and the frequency of modulation. It is found that the low frequency symmetric thermal modulation is destabilizing while moderate and high frequency symmetric modulation is always stabilizing. The asymmetric modulation and lower wall temperature modulations are, in general, stabilizing while the system becomes unstable for large values of Darcy–Prandtl number and for small frequencies. It is shown that in general the gravity modulation produces a stabilizing effect on the onset of convection for moderate and high frequency. The small frequency gravity modulation is found to have destabilizing effect on the stability of the system.     

Plane waves in linear elastic materials with voids

The behavior of plane harmonic waves in a linear elastic material with voids is analyzed. There are two dilational waves in this theory, one is predominantly the dilational wave of classical linear elasticity and the other is predominantly a wave carrying a change in the void volume fraction. Both waves are found to attenuate in their direction of propagation, to be dispersive and dissipative. At large frequencies the predominantly elastic wave propagates with the classical elastic dilational wave speed, but at low frequencies it propagates at a speed less than the classical speed. It makes a smooth but relatively distinct transition between these wave speeds in a relatively narrow range of frequency, the same range of frequency in which the specific loss has a relatively sharp peak. Dispersion curves and graphs of specific loss are given for four particular, but hypothetical, materials, corresponding to four cases of the solution.     

A physical perspective of the length scales in gradient elasticity through the prism of wave dispersion

The goal of this study is to understand the physical meaning and evaluate the intrinsic length scale parameters, featured in the theories of gradient elasticity, by deploying the analytical treatment and experimental measurements of the dispersion of elastic waves. The developments are focused on examining the propagation of longitudinal waves in an aluminum rod with periodically varying cross-section. First, the analytical solution for the dispersion relationship, based on the periodic cell analysis of a bi-layered laminate and Bloch theorem, is compared to two competing models of gradient elasticity. It is shown that the customary gradient elastic model with two length-scale parameters is able to capture the dispersion accurately up to the beginning of the first band gap. On the other hand, the gradient elastic model with an additional length scale (affiliated with the fourth-order time derivative in the field equation) is shown to capture not only the first dispersion branch before the band gap, but also the band gap itself and the preponderance of the second branch. Closed form relations between the microstructure parameters and the intrinsic length scales are obtained for both gradient elasticity models. By way of the asymptotic treatment in the limit of a weak contrast between the laminae, a clear physical meaning and scaling of the length-scale parameters was established in terms of: (i) the microstructure (given by the size of the unit cell and the contrast between the laminae), and (ii) thus induced dispersion relationship (characterized by the location and the width of the band gap). The analysis is verified through an experimental observation of wave dispersion, and wave attenuation within the band gap. A comparison between the analytical treatment, the gradient elastic model with three intrinsic length scales, and experimental measurements demonstrates a good agreement over the range of frequencies considered.     

Wave propagation analysis of porous functionally graded piezoelectric nanoplates with a visco-Pasternak foundation

This study investigates the size-dependent wave propagation behaviors under the thermoelectric loads of porous functionally graded piezoelectric(FGP) nanoplates deposited in a viscoelastic foundation. It is assumed that(i) the material parameters of the nanoplates obey a power-law variation in thickness and(ii) the uniform porosity exists in the nanoplates. The combined effects of viscoelasticity and shear deformation are considered by using the Kelvin-Voigt viscoelastic model and the refined hi...     

A Gradient-Based Optimization Method for the Design of Layered Phononic Band-Gap Materials

Phononic materials with specific band-gap characteristics at desired frequency ranges are in great demand for vibration and noise isolation, elastic wave filters, and acoustic devices.The attenuation coefficient curve depicts both the frequency range of band gap and the attenuation of elastic wave, where the frequency ranges corresponding to the none-zero attenuation coefficients are band gaps. Therefore, the band-gap characteristics can be achieved through maximizing the attenuation coefficient at the corresponding frequency or within the corresponding frequency range. Because the attenuation coefficient curve is not smooth in the frequency domain,the gradient-based optimization methods cannot be directly used in the design optimization of phononic band-gap materials to achieve the maximum attenuation within the desired frequency range. To overcome this difficulty, the objective of maximizing the attenuation coefficient is transformed into maximizing its Cosine, and in this way, the objective function is smoothed and becomes differentiable. Based on this objective function, a novel gradient-based optimization approach is proposed to open the band gap at a prescribed frequency range and to further maximize the attenuation efficiency of the elastic wave at a specific frequency or within a prescribed frequency range. Numerical results demonstrate the effectiveness of the proposed gradient-based optimization method for enhancing the wave attenuation properties.     

Determination of relaxation and retardation spectrum by inverse functional filtering

The article is devoted for the determination of the relaxation and retardation spectrum (RRS) from monotonic time- and frequency-domain material functions by the inverse functional filters executing discrete convolution algorithms for geometrically spaced data. It is shown that the problem of RRS determination from a wide variety of material functions leads to the three inverse filtering tasks on a logarithmic time or frequency scale with the three specific frequency responses concerning: (i) the time-domain functions, (ii) the real parts and (iii) the imaginary parts of the frequency-domain functions, and three algorithms (having the versions with even and odd number of coefficients) are to be applied to: (i) time-domain compliance and modulus functions, (ii) their derivatives, and (iii) frequency-domain functions. It is demonstrated that ill-posedness of an inverse filter manifests as large sampling-rate-dependent noise amplification coefficients. A novel regularization strategy allowing to ensure the desired noise immunity is proposed based on choosing sampling rate for geometrically spaced data. The performance of the algorithms is investigated. Optimal sampling rates are disclosed for specific material functions. The frequency range of 2–3 decades is established to be optimal for the recovery of a single RRS point estimate ensuring maximum accuracy with reasonable noise immunity. Practical algorithms are proposed for recovering RRS from the real and imaginary parts of frequency-domain functions. Some known non-parametric methods are compared with the suggested functional filters.     

Modulation of an internal gravity wave packet in a stratified shear flow

Modulations of an internal gravity wave packet in a stratified shear flow are discussed in the weakly nonlinear and weakly dispersive context. It is shown that the modulations are described by a variable coefficient nonlinear Schrödinger equation when the modulations are confined to the direction of wave propagation. Transverse modulations couple the nonlinear Schrödinger equation to the mean flow equations. For long waves, it is shown that the modulation equations may be somewhat simplified. An Appendix describes the equations governing long wave resonance.