Inhomogeneous plane waves in compressible viscous fluids

The propagation of inhomogeneous plane waves in a compressible viscous fluid is considered. The frequency and the slowness vector are both allowed to be complex. There are seen to be two types of solutions: (a) two transverse waves, which involve no density or pressure fluctuations, (b) a longitudinal wave, which involves no fluctuations in vorticity. For each type, a propagation condition is obtained giving the (complex) squared length of the slowness vector as a function of frequency. Each depends also on the viscosities. It is seen how to recover the incompressible case as the limit in which the inviscid acoustic wave speed tends to infinity. Each wave is shown to be linearly stable for real frequencies. These waves are attenuated in space and time but nevertheless it is possible to define constant weighted mean values (over a cycle of the propagating part of the wave) of the energy density, energy flux and dissipation. The energy-dissipation equation and the propagation conditions are used to derive relationships between these constant weighted means, some of which are generalizations to compressible fluids of previously known results for incompressible fluids. Explicit expressions in terms of frequency are given for the weighted means.     

饱和黏弹性多孔介质中的平面波及能量耗散

研究了流体饱和不可压黏弹性多孔介质中的非均匀平面波及其能量流和能量耗散规律. 在流 相和固相物质微观不可压、固相骨架宏观服从积分型本构关系和小变形的假定下，利用 Helmholtz分解，得到了饱和黏弹性多孔介质中非均匀平面波的一般解以及纵波、横波相速 度和衰减率等的解析表达式，分析了平面波传播矢量和衰减矢量之间的关系. 数值结果表明 孔隙流体与固相骨架间的相互作用以及固相骨架的黏性对波的相速度、衰减率等有着显著的 影响. 同时，得到了饱和黏弹性多孔介质的能量方程，给出了能量流矢量和能量耗散率. 对 非均匀平面纵波和横波，推导了平均能量流矢量和平均能量耗散率的解析表达式.     

Superposition of finite-amplitude waves propagating in a principal plane of a deformed Mooney–Rivlin material

Here we consider finite-amplitude wave motions in Mooney–Rivlin elastic materials which are first subjected to a static homogeneous deformation (prestrain). We assume that the time-dependent displacement superimposed on the prestrain is along a principal axis of the prestrain and depends on two spatial variables in the principal plane orthogonal to this axis. Thus all waves considered here are linearly polarized along this axis. After retrieving known results for a single homogeneous plane wave propagating in a principal plane, a superposition of an arbitrary number of sinusoidal homogeneous plane waves is shown to be a solution of the equations of motion. Also, inhomogeneous plane wave solutions with complex wave vector in a principal plane and complex frequency are obtained. Moreover, appropriate superpositions of such inhomogeneous waves are also shown to be solutions. In each case, expressions are obtained for the energy density and energy flux associated with the wave motion.     

Inhomogeneous waves in anisotropic dissipative bodies

Wave propagation in anisotropic dissipative bodies is considered through the form of inhomogeneous waves. The dissipativity of the body is characterized by a restriction of thermodynamic character on the viscoelastic tensor. As a consequence, the divergence of the usual (time-averaged) energy flux is proved to be negative. This in turn is shown to imply that the amplitude of the wave decays in the direction of the energy flux and the angle, subtended by the imaginary part of the wave vector and the energy flux, is acute. Moreover, the symmetry of the viscoelastic tensor and the positive definiteness of its real part imply that also the imaginary part of the wave vector subtends an acute angle with the energy flux. Such properties are shown to hold in both solids and fluids. Because of anisotropy, the angle subtended by the real and imaginary parts of the wave vector need not be acute, which means that the amplitude may increase in the direction of phase propagation.     

Singularities on wave fronts of slow waves in anisotropic fluid-saturated porous media

Energy focusing is found on the wave fronts of slow waves, which is a new propagation characteristic for slow waves in fluid-saturated porous materials. The material parameters, as well as the propagation directions, are chosen as the control parameters. Combined with the two axial variables, the influence of the anisotropy of the solid skeleton and pore fluid parameters on the propagation characteristic of slow waves in anisotropic fluid-saturated porous materials is discussed. The correspondence between the focusing on the wave fronts and the contours of zero Gaussian curvature on the slowness surface is explored. The development of the focusing patterns is investigated and the distinct trends in the energy flux focusing structures are revealed. This is helpful in further understanding the roles of the pore fluid in the damage of the fluid-saturated porous media.     

Nonlinear elastic effects in graphite/epoxy: An analytical and numerical prediction of energy flux deviation

Manipulating acoustic wave propagation through a material have several interdisciplinary applications. Here we predict shift in energy flux deviation for acoustic waves propagating in unidirectional graphite/epoxy due to applied normal and shear stresses using both an analytical model, using acoustoelastic continuum theory, and a finite element discrete numerical model. The acoustoelastic theory predicts that the quasi-transverse (QT) wave exhibits larger shifts in energy flux deviation compared to quasi-longitudinal (QL) or the pure transverse (PT) due to an applied shear stress for fiber orientation angle ranging from 0° to 60°. Due to an applied shear stress the QT wave exhibits a shift in energy flux deviation at 0° fiber orientation angle as compared to normal stress case where the flux deviation and its load induced shift are both zero. A finite element model (FEM) is developed where equations of motion include the effect of nonlinear elastic coefficients. Element equations were integrated in time using Newmark’s method to determine the shift in energy flux deviations in graphite/epoxy for different loading cases. The energy flux shift of QT waves predicted by FEM for fiber orientation angles from 0° to 60° for applied shear stress case is in excellent agreement with acoustoelastic theory. Because energy shift magnitudes are not small, it is possible to experimentally measure these shifts and calibrate shifts with respect to load type (normal/shear) and magnitude.     

Computational studies of inertia-gravity waves radiated from upper tropospheric jets

The generation and physical characteristics of inertia-gravity waves radiated from an unstable forced jet at the tropopause are investigated through high-resolution numerical simulations of the three-dimensional Navier–Stokes anelastic equations. Such waves are induced by Kelvin–Helmholtz instabilities on the flanks of the inhomogeneously stratified jet. From the evolution of the averaged momentum flux above the jet, it is found that gravity waves are continuously radiated after the shear-stratified flow reaches a quasi-equilibrium state. The time–vertical coordinate cross-sections of potential temperature show phase patterns indicating upward energy propagation. The sign of the momentum flux above and below the jet further confirms this, indicating that the group velocity of the generated waves is pointing away from the jet core region. Space–time spectral analysis at the upper flank level of the jet shows a broad spectral band, with different phase speeds. The spectra obtained in the stratosphere above the jet show a shift toward lower frequencies and larger spatial scales compared to the spectra found in the jet region. The three-dimensional character of the generated waves is confirmed by analysis of the co-spectra of the spanwise and vertical velocities. Imposing the background rotation modifies the polarization relation between the horizontal wind components. This out-of-phase relation is evidenced by the hodograph of the horizontal wind vector, further confirming the upward energy propagation. The background rotation also causes the co-spectra of the waves high above the jet core to be asymmetric in the spanwise modes, with contributions from modes with negative wavenumbers dominating the co-spectra. Dedicated to the memory of our colleague Dr. Binson Joseph     

Wave propagation and reflection-transmission in a stratified viscoelastic solid

The paper investigates time-harmonic wave propagation in continuously stratified solids and provides the results of a reflection-transmission process generated by a layer sandwiched between homogeneous half-spaces. The layer is continuously stratified and allows for jump discontinuities at a finite number of planes. The dissipative effects are accounted for through the classical Boltzmann law of viscoelasticity. By using displacement and traction as convenient vector variables, the governing equations are considered in a vector Volterra integral equation and the solution is determined by means of a matricant. Next the matricant is applied to determine the reflection and transmission coefficients of a layer, with a generic piecewise continuous profile of the material properties. The reflection-transmission process produced by an obliquely incident wave, is considered for horizontally-polarized waves. The low-frequency approximation is derived for the reflection and transmission coefficients. Next, the high-frequency approximation is investigated by a WKB-like procedure which involves a complex valued frequency-dependent shear modulus. The displacement solution is obtained for the forward- and the backward-propagating waves in the layer along with the reflection and transmission coefficients.     

Analytical Solution for Isothermal Flow in a Shock Tube Containing Rigid Granular Material

Analytical solution of shock wave propagation in pure gas in a shock tube is usually addressed in gas dynamics. However, such a solution for granular media is complex due to the inclusion of parameters relating to particles configuration within the medium, which affect the balance equations. In this article, an analytical solution for isothermal shock wave propagation in an isotropic homogenous rigid granular material is presented, and a closed-form solution is obtained for the case of weak shock waves. Fluid mass and momentum equations are first written in wave and (mathematical) non-conservation forms. Afterwards by redefining the sound speed of the gas flowing inside the pores, an analytical solution is obtained using the classical method of characteristics, followed by Taylor’s series expansion based on the assumption of weak flow which finally led to explicit functions for velocity, density and pressure. The solution enables plotting gas velocity, density and pressure variations in the porous medium, which is of high interest in the design of granular shock isolators.     

On energy balance and the structure of radiated waves in kinetics of crystalline defects

Traveling waves, with well-known closed form expressions, in the context of the defects kinetics in crystals are excavated further with respect to their inherent structure of oscillatory components. These are associated with, so called, Frenkel–Kontorova model with a piecewise quadratic substrate potential, corresponding to the symmetric as well as asymmetric energy wells of the substrate, displacive phase transitions in bistable chains, and brittle fracture in triangular lattice strips under mode III conditions. The paper demonstrates that the power expended theorem holds so that the sum of rate of working and the rate of total energy flux into a control strip moving steadily with the defect equals the rate of energy sinking into the defect, in the sense of N.F. Mott. In the conservative case of the Frenkel–Kontorova model with asymmetric energy wells, this leads to an alternative expression for the mobility in terms of the energy flux through radiated lattice waves. An application of the same to the case of martensitic phase boundary and a crack, propagating uniformly in bistable chains and triangular lattice strips, respectively, is also provided and the energy release is expressed in terms of the radiated energy flux directly. The equivalence between the well-known expressions and their alternative is established via an elementary identity, which is stated and proved in the paper as the zero lemma. An intimate connection between the three distinct types of defects is, thus, revealed in the framework of energy balance, via a structural similarity between the corresponding variants of the ‘zero’ lemma containing the information about radiated energy flux. An extension to the dissipative models, in the presence of linear viscous damping, is detailed and analog of the zero lemma is proved. The analysis is relevant to the dynamics of dislocations, brittle cracks, and martensitic phase boundaries, besides possible applications to analogous physical contexts which are marked by macroscopic energy release through emission of waves and possibly linear viscous damping.     

Finite-amplitude Love waves in a pre-stressed compressible elastic half-space with a double surface layer

The dispersive behavior of finite-amplitude time-harmonic Love waves propagating in a pre-stressed compressible elastic half-space overlaid with two compressible elastic surface layers of finite thickness is investigated. The half-space and layers are made of different pre-stressed compressible neo-Hookean materials. The dispersion relation which relates wave speed and wavenumber is obtained in explicit form. Results for the energy density and energy flux of the waves are also presented. The special case where the interfaces between the layers and the half-space are principal planes of the left Cauchy–Green deformation tensor is also investigated. Numerical results are presented showing the variation of the Love wave speed with the pre-stress and the propagation angle.     

Nonreflecting boundary conditions based on nonlinear multidimensional characteristics

Nonlinear characteristic boundary conditions based on nonlinear multidimensional characteristics are proposed for 2‐ and 3‐D compressible Navier–Stokes equations with/without scalar transport equations. This approach is consistent with the flow physics and transport properties. Based on the theory of characteristics, which is a rigorous mathematical technique, multidimensional flows can be decomposed into acoustic, entropy, and vorticity waves. Nonreflecting boundary conditions are derived by setting corresponding characteristic variables of incoming waves to zero and by partially damping the source terms of the incoming acoustic waves. In order to obtain the resulting optimal damping coefficient, analysis is performed for problems of pure acoustic plane wave propagation and arbitrary flows. The proposed boundary conditions are tested on two benchmark problems: cylindrical acoustic wave propagation and the wake flow behind a cylinder with strong periodic vortex convected out of the computational domain. This new approach substantially minimizes the spurious wave reflections of pressure, density, temperature, and velocity as well as vorticity from the artificial boundaries, where strong multidimensional flow effects exist. The numerical simulations yield accurate results, confirm the optimal damping coefficient obtained from analysis, and verify that the method substantially improves the 1‐D characteristics‐based nonreflecting boundary conditions for complex multidimensional flows. Copyright © 2009 John Wiley & Sons, Ltd.     

Analysis of the wave solution of the elastokinetic equations of a Cosserat continuum for the case of bulk plane waves

A study is made of waves in a Cosserat continuum, whose strain state is characterized by independent displacement and rotation vectors. The propagation of longitudinal and transverse bulk waves is considered. Wave solutions are sought in the form of wave trains specified by a Fourier spectrum of arbitrary shape. It is shown that if the solution is sought in the form of three components of the displacement vector and three components of the rotation vector which depend on time and the longitudinal coordinate, the initial system is split into two systems, one of which describes longitudinal waves, and the other transverse waves. For waves of both types, dispersion relations and analytical solutions in displacement are obtained. The dispersion characteristics of the solutions obtained differ from the dispersion characteristics of the corresponding classical elastic solutions. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 2, pp. 196–203, March–April, 2008.     

Spherically imploding detonation waves initiated by two-step divergent detonation

The detonation chamber developed by K. Terao and H. G. Wagner in Göttingen has been improved, in order to concentrate the combustion energy more effectively to a focus, so that imploding detonation waves are initiated by two-step divergent detonation waves in a hemispherical space having an effective diameter of 500 mm. The imploding detonation waves are investigated by measuring their propagation velocity using ion probes and pressure variations at different positions in the space by a piezoelectric pressure transducer, while the temperature in the implosion center is measured by a laser light scattering method. The results suggest that the peak pressure at the detonation front increases with the propagation to the center more rapidly than that in the Göttingen apparatus, while the propagation velocity is almost the same. A temperature from 10⁷ K to 10⁸ K is also observed in the focus of the imploding detonation waves.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

A general approximate‐state Riemann solver for hyperbolic systems of conservation laws with source terms

An approximate‐state Riemann solver for the solution of hyperbolic systems of conservation laws with source terms is proposed. The formulation is developed under the assumption that the solution is made of rarefaction waves. The solution is determined using the Riemann invariants expressed as functions of the components of the flux vector. This allows the flux vector to be computed directly at the interfaces between the computational cells. The contribution of the source term is taken into account in the governing equations for the Riemann invariants. An application to the water hammer equations and the shallow water equations shows that an appropriate expression of the pressure force at the interface allows the balance with the source terms to be preserved, thus ensuring consistency with the equations to be solved as well as a correct computation of steady‐state flow configurations. Owing to the particular structure of the variable and flux vectors, the expressions of the fluxes are shown to coincide partly with those given by the HLL/HLLC solver. Computational examples show that the approximate‐state solver yields more accurate solutions than the HLL solver in the presence of discontinuous solutions and arbitrary geometries. Copyright © 2006 John Wiley & Sons, Ltd.     

Energy and dissipation of inhomogeneous plane waves in thermoelasticity

Inhomogeneous small-amplitude plane waves of (complex) frequency ω are propagated through a linear dissipative material. For thermoelasticity we derive an energy-dissipation equation that contains all the quadratic dependence on the field quantities, see Eq. (10). In addition, we derive a new energy-dissipation equation (Eq. (22)) involving the total energy density which contains terms linear in the field quantities as well as the usual quadratic terms. The terms quadratic in the small quantities in the energy density, energy flux and dissipation give rise to inhomogeneous plane waves of frequency 2ω and to (attenuated) constant terms. Usually these quadratic quantities are time-averaged and only the attenuated constant terms remain. We derive a new result in thermoelasticity for these terms, see Eq. (54). The present innovation is to retain the terms of frequency 2ω, since they are comparable in magnitude to the attenuated constant terms, and a new result, see Eq. (44), is derived for a general energy-dissipation equation that connects the amplitudes of the terms of the energy density, energy flux and dissipation that have frequency 2ω. Furthermore, for dissipative waves or inhomogeneous conservative waves the (complex) group velocity is related to these amplitudes rather than to the attenuated constant terms as it is for homogeneous waves in conservative materials.     

Finite amplitude spherically symmetric wave propagation in a prestressed hyperelastic shell

Spherically symmetric finite amplitude wave propagation in a prestressed compressible hyperelastic spherical shell is considered. The prestress results from quasi-static application of internal pressure and a numerical solution for this elastostatic problem is obtained first. Dynamic change of the internal pressure results in the propagation of a spherically symmetric wave. A Godunov type finite difference scheme is proposed for the solution of the wave propagation problem and numerical results, which are valid until the first reflection, are presented for a particular isotropic strain energy function and for the special cases of sudden removal and sudden increase of the internal pressure.     

多个椭圆柱波浪力的一种解析解

波浪在大尺寸结构表面产生不可忽略的散射波, 该散射波在多柱体体系中继续传播, 并在同体系中的其他柱体上产生高次散射波. 本文基于椭圆坐标系和绕射波理论首先推导了波浪作用下椭圆单柱体产生的散射波压力公式, 随后考虑该散射波在多柱体系中的传播, 将其视为第二次入射波, 推导出柱体上第二次散射波压力公式, 同理可以推导出高次散射波压力公式, 最后得到椭圆多柱体波浪力解析解, 并用数值解验证了本文解析方法的正确性. 本文以双柱体和四柱体体系为例, 分析了不同参数(波数、净距、波浪入射角度等)下, 高次散射波对柱体上波浪作用的影响. 结果表明: 波数较大的情况下, 高次散射波引起柱体上的波浪力不能忽略; 结构间距较大的情况下, 虽然高次波的作用有减小的趋势但仍然明显; 高次散射波来自多个柱体对入射波的散射, 柱体数目的增加后, 高次波的影响会增加, 结构所受的高次波作用因参数变化而起的波动会变剧烈; 高次波对上游柱体波浪力的贡献较对下游柱体的贡献大.      

Complex surface rays associated with inhomogeneous skimming and Rayleigh waves

The Rayleigh wave, that propagates at the free surface of semi-infinite anisotropic medium, is composed of three inhomogeneous partial waves, each propagating along the surface with a different attenuation along the depth. Since this wave does not exhibit an attenuation on the surface, let us call it the homogeneous Rayleigh wave. The associated slowness corresponds to the real solution of the Rayleigh dispersion equation. Besides this classical solution, an infinite number of complex solutions of the Rayleigh dispersion equation exits. For such particular Rayleigh waves, the slowness vector, i.e. the identical component on the surface of the slowness of each partial waves, is taken to be complex. Thus, these Rayleigh waves are attenuated on the surface and as shown here, their attenuation is normal to the ray direction (or the energy velocity direction). Similarly to the infinite inhomogeneous plane waves which can be associated with complex rays, we call these waves, inhomogeneous Rayleigh waves. We use the inhomogeneous skimming waves, which are inhomogeneous plane waves, and the inhomogeneous Rayleigh waves to explain differently the usual diffraction phenomena on the free surface which cannot be explained by the real ray theory. For example, the arrival time of the wave packet observed beyond the cusp is in perfect accordance with the arrival time of some specific inhomogeneous Rayleigh waves. We show that these results are in agreement with the computation of the Green function. They apply to the theory of surface waves in linear elastodynamics with intrinsic anisotropy as well as to the theory of surface waves in linearised (incremental) elastodynamics with strain-induced anisotropy (also known as small-amplitude waves superimposed on the large static homogeneous deformation of a non-linear solid).