Sound wave propagation through incompressible flows

This paper derives a variable coefficient, multidimensional Burgers equation which models the propagation of a nonlinear sound wave through an incompressible background flow. The equation is derived from the compressible Euler equations by the combination of a weakly nonlinear acoustics expansion for the sound wave and an incompressible expansion for the background flow. The main effect of the incompressible flow on the sound wave is the advection of the sound wave by the transverse velocity component of the flow.     

有限变形弹性杆中的非线性波及周期解

利用Ham ilton变分原理,导出了计及有限变形和横向Possion效应的弹性杆中非线性纵向波动方程.利用Jacob i椭圆正弦函数展开和第三类Jacob i椭圆函数展开法,对该方程和截断的非线性方程进行求解,得到了非线性波动方程的两类准确周期解及相应的孤波解和冲击波解,讨论了这些解存在的必要条件.     

梁中非线性弯曲波传播特性的研究

研究了梁中的非线性弯曲波的传播特性,同时考虑了梁的大挠度引起的几何非线性效应和 梁的转动惯性导致的弥散效应,利用Hamilton变分法建立了梁中非线性弯曲波的波动方程. 对该方程进行了定性分析,在不同的条件下,该方程在相平面上存在同宿轨道或异宿轨道, 分别对应于方程的孤波解或冲击波解. 利用Jacobi椭圆函数展开法,对该非线性方程进行 求解,得到了非线性波动方程的准确周期解及相对应的孤波解和冲击波解,讨论了这些解存 在的必要条件,这与定性分析的结果完全相同. 利用约化摄动法从非线性弯曲波动方程中导 出了非线性Schr\"{o}dinger方程,从理论上证明了考虑梁的大挠度和转动惯性时梁中存在 包络孤立波.     

On the Interphase Heat Transfer Effect on Nonlinear Wave Propagation in a Gas-Liquid Mixture

The effect of interphase heat transfer on shock wave propagation is investigated. A multiwave nonlinear equation which in the limiting case of the absence of heat transfer decomposes into two classic generalizations of the Boussinesq equations is derived. Quasi-isothermal and quasi-adiabatic propagation regimes for which the heat transfer is fairly intense are considered. For both regimes, nonlinear equations describing the wave propagation are obtained. The equation describing the first regime is investigated in detail. Exact analytic solutions of this equation are given and used to study the shock wave structures and the solitary wave behavior. Formulas for the dependence of the heat transfer rate on the equilibrium-mixture parameters are obtained.     

Anti-plane shear waves in a fibre-reinforced composite with a non-linear imperfect interface

The propagation of non-linear elastic anti-plane shear waves in a unidirectional fibre-reinforced composite material is studied. A model of structural non-linearity is considered, for which the non-linear behaviour of the composite solid is caused by imperfect bonding at the “fibre–matrix” interface. A macroscopic wave equation accounting for the effects of non-linearity and dispersion is derived using the higher-order asymptotic homogenisation method. Explicit analytical solutions for stationary non-linear strain waves are obtained. This type of non-linearity has a crucial influence on the wave propagation mode: for soft non-linearity, localised shock (kink) waves are developed, while for hard non-linearity localised bell-shaped waves appear. Numerical results are presented and the areas of practical applicability of linear and non-linear, long- and short-wave approaches are discussed.     

Rayleigh waves in an incompressible elastic half-space overlaid with a water layer under the effect of gravity

This paper is concerned with the propagation of Rayleigh waves in an incompressible isotropic elastic half-space overlaid with a layer of non-viscous incompressible water under the effect of gravity. The authors have derived the exact secular equation of the wave which did not appear in the literature. Based on it the existence of Rayleigh waves is considered. It is shown that a Rayleigh wave can be possible or not, and when a Rayleigh wave exists it is not necessary unique. From the exact secular equation the authors arrive immediately at the first-order approximate secular equation derived by Bromwich [Proc. Lond. Math. Soc. 30:98–120, 1898]. When the layer is assumed to be thin, a fourth-order approximate secular equation is derived and of which the first-order approximate secular equation obtained by Bromwich is a special case. Some approximate formulas for the velocity of Rayleigh waves are established. In particular, when the layer being thin and the effect of gravity being small, a second-order approximate formula for the velocity is created which recovers the first-order approximate formula obtained by Bromwich [Proc. Lond. Math. Soc. 30:98–120, 1898]. For the case of thin layer, a second-order approximate formula for the velocity is provided and an approximation, called global approximation, for it is derived by using the best approximate second-order polynomials of the third- and fourth-powers.     

NONLINEAR WAVES AND PERIODIC SOLUTION IN FINITE DEFORMATION ELASTIC ROD

A nonlinear wave equation of elastic rod taking account of finite deformation, transverse inertia and shearing strain is derived by means of the Hamilton principle in this paper. Nonlinear wave equation and truncated nonlinear wave equation are solved by the Jacobi elliptic sine function expansion and the third kind of Jacobi elliptic function expansion method. The exact periodic solutions of these nonlinear equations are obtained, including the shock wave solution and the solitary wave solution. The necessary condition of exact periodic solutions, shock solution and solitary solution existence is discussed.     

Stretched Cartesian grids for solution of the incompressible Navier–Stokes equations

Two Cartesian grid stretching functions are investigated for solving the unsteady incompressible Navier–Stokes equations using the pressure–velocity formulation. The first function is developed for the Fourier method and is a generalization of earlier work. This function concentrates more points at the centre of the computational box while allowing the box to remain finite. The second stretching function is for the second‐order central finite difference scheme, which uses a staggered grid in the computational domain. This function is derived to allow a direct discretization of the Laplacian operator in the pressure equation while preserving the consistent behaviour exhibited by the uniform grid scheme. Both functions are analysed for their effects on the matrix of the discretized pressure equation. It is shown that while the second function does not spoil the matrix diagonal dominance, the first one can. Limits to stretching of the first method are derived for the cases of mappings in one and two directions. A limit is also derived for the second function in order to prevent a strong distortion of a sine wave. The performances of the two types of stretching are examined in simulations of periodic co‐flowing jets and a time developing boundary layer. Copyright © 2000 John Wiley & Sons, Ltd.     

Plane strain dynamics of elastic solids with intrinsic boundary elasticity, with application to surface wave propagation

In this paper, in a development of the static theory derived by Steigmann and Ogden (Proc. Roy. Soc. London A 453 (1997) 853), we establish the equations of motion for a non-linearly elastic body in plane strain with an elastic surface coating on part or all of its boundary. The equations of (linearized) incremental motions superposed on a finite static deformation are then obtained and applied to the problem of (time-harmonic) surface wave propagation on a pre-stressed incompressible isotropic elastic half-space with a thin coating on its plane boundary. The secular equation for (dispersive) wave speeds is then obtained in respect of a general form of incompressible isotropic elastic strain-energy function for the bulk material and a general energy function for the coating material. Specialization of the form of strain-energy function enables the secular equation to be cast as a quartic equation and we therefore focus on this for illustrative purposes. An explicit form for the secular equation is thereby obtained. This involves a number of material parameters, including residual stress and moment in the properties of the coating. It is shown how this equation relates to previous work on waves in a half-space with an overlying thin layer set in the classical theory of isotropic elasticity and, in particular, the significant effect of omission of the rotatory inertia term, even at small wave numbers, is emphasized. Corresponding results for a membrane-type coating, for which the bending moment, inertia and residual moment terms are absent, are also obtained. Asymptotic formulas for the wave speed at large wave number (high frequency) are derived and it is shown how these results influence the character of the wave speed throughout the range of wave number values. A bifurcation criterion is obtained from the secular equation by setting the wave speed to zero, thereby generalizing the bifurcation results of Steigmann and Ogden (Proc. Roy. Soc. London A 453 (1997) 853) to the situation in which residual stress and moment are present in the coating. Numerical results which show the dependence of the wave speed on the various material parameters and the finite deformation are then described graphically. In particular, features which differ from those arising in the classical theory are highlighted.     

Geometrical nonlinear waves in finite deformation elastic rods

IntroductionSomenewphenomenaofnonlinearwavesinthesolidmediumsuchasshockwave ,solitarywaveetc.arepaidmoreattentiontoincreasinglybyresearchersbecausetheytakeonalotofimportantproperties.ItistheoreticallyanalyzedinRefs.[1 -6]thattheformationmechanismsofshockwaveandsolitarywaveintheelasticthinrodsaswellastheirpropagationproperties.TheexistenceofsolitarywaveintheelasticmediumsuchasarodandaplatehasbeenverifiedinRef.[7]byexperiments.Shockwaveandsolitarywavearesteadilypropagatingtraveling_wavesgenerat…     

Evolution of weak noise and regular waves on dissipative shock fronts described by the Burgers model

The interaction of weak noise and regular signals with a shock wave having a finite width is studied in the framework of the Burgers equation model. The temporal realization of the random process located behind the front approaches it at supersonic speed. In the process of moving to the front, the intensity of noise decreases and the correlation time increases. In the central region of the shock front, noise reveals non-trivial behaviour. For large acoustic Reynolds numbers the average intensity can increase and reach a maximum value at a definite distance. The behaviour of statistical characteristics is studied using linearized Burgers equation with variable coefficients reducible to an autonomous equation. This model allows one to take into account not only the finite width of the front, but the attenuation and diverse character of initial profiles and spectra as well. Analytical solutions of this equation are derived. Interaction of regular signals of complex shape with the front is studied by numerical methods. Some illustrative examples of ongoing processes are given. Among possible applications, the controlling the spectra of signals, in particular, noise suppression by irradiating it with shocks or sawtooth waves can be mentioned.     

弹性介质中几何和运动非线性问题的类孤波解

将弹性介质 的几何和运动非线性方程简化成具有电磁场中的Ｂｏｒｎ－Ｉｎｆｅｌｄ方程的形式，并证明了该方程的类孤波解的存在．     

位移浅水内孤立波

研究两层浅水系统中的内孤立波,该系统由两层常密度不可压缩无黏性水组成。利用Lagrange坐标和Hamilton原理,推导了两层浅水系统的位移浅水内波方程,并进一步导出了两层浅水系统的位移内孤立波解。数值实验表明,位移内孤立波与经典的KdV内孤立波吻合很好,说明Lagrange坐标和Hamilton方法适用于内波分析,可以为构造内波分析的保辛方法提供一种途径。     

Nonlinear nonequilibrium kinetic model of the Boltzmann equation for a gas with power-law interaction

A model kinetic equation approximating the Boltzmann equation on a wide range of the intensities of nonequilibrium states of gases is derived to describe rarefied gas flows. The kinetic model is based on a distribution function dependent on the absolute velocity of gas particles. Themodel kinetic equation possesses a high computational efficiency and the problem of shock wave structure is solved on its basis. The calculated and experimental data for argon are compared.     

Behaviour of the acceleration and shock waves in materials with internal state variables

The explicit expressions for the change in the amplitudes of one-dimensional acceleration and shock waves propagating through arbitrary homogeneous materials described by the strain and internal state variables/parameters/are derived. The existence of a critical amplitude β for the acceleration wave and a critical strain gradient λ for the shock wave is established. For an infinitesimal shock wave the general form of the solution of the governing differential equation is furnished. The differential equations for the amplitudes of these two kind of waves are applied to an elastic-viscoplastic material.     

Modelling a shock wave source for seismic well survey

The shock wave seismic source has specific advantages for reservoir survey, particularly in the technique of vertical seismic profiling. This paper is concerned with the modelling of the characteristics of a shock wave device under conditions typical of its use as a surface seismic source. Such modelling necessitates numerical solution of the time-dependent conservation equations for axi-symmetric geometry under conditions where the fluid in the surface borehole may be air or water. A version of the piecewise parabolic method has been used to take account of the two-phase behaviour based on the formal properties of the well-known Tait equation of state for water. Results are presented which enable predictions to be made for use in the field of the pressure and velocity signatures of the shock wave source in terms of axial and radial profiles. Such information is significant in the assessment of the degrees of compressional- and shear-wave energies delivered by the source in vertical seismic profiling surveys of reservoirs. Received 6 November 1997 / Accepted 3 March 1998     

Growth of discontinuities in fluids with internal relaxation

The growth equation for weak discontinuities headed by wave fronts of arbitrary shape in a relaxing gas flow is derived along the orthogonal trajectories of the wave fronts. An explicit criteria for the growth and decay of weak discontinuities along their orthogonal trajectories is given. It is concluded that the internal relaxation processes in the flow as well as the wave front curvature both have a stabilizing effect on the tendency of the wave surface to grow into a shock in the sense that they cause the shock formation time to increase.     

Shock structure in bubbly liquids: comparison of direct numerical simulations and model equations

Shock wave structure in a bubbly mixture composed of a cluster of gas bubbles in a quiescent liquid with initial void fractions around 10% inside a 3D rectangular domain excited by a sudden increase in the pressure at one boundary is investigated using the front tracking/finite volume method. The effects of bubble/bubble interactions and bubble deformations are, therefore, investigated for further modeling. The liquid is taken to be incompressible while the bubbles are assumed to be compressible. The gas pressure inside the bubbles is taken uniform and is assumed to vary isothermally. Results obtained for the pressure distribution at different locations along the direction of propagation show the characteristics of one-dimensional unsteady shock propagation evolving towards steady-state. The steady-state shock structures obtained by the present direct numerical simulations, which show a transition from A-type to C-type steady-state shock structures, are compared with those obtained by the classical Rayleigh–Plesset equation and by a modified Rayleigh–Plesset equation accounting for bubble/bubble interactions in the mean-field theory.      

Incompressible SPH simulation of wave breaking and overtopping with turbulence modelling

In this paper a truly incompressible version of the smoothed particle hydrodynamics (SPH) method is presented to investigate the surface wave overtopping. SPH is a pure Lagrangian approach which can handle large deformations of the free surface with high accuracy. The governing equations are solved based on the SPH particle interaction models and the incompressible algorithm of pressure projection is implemented by enforcing the constant particle density. The two‐equation k–ε model is an effective way of dealing with the turbulence and vortices during wave breaking and overtopping and it is coupled with the incompressible SPH numerical scheme. The SPH model is employed to reproduce the experiment and computations of wave overtopping of a sloping sea wall. The computations are validated against the experimental and numerical data found in the literatures and good agreement is observed. Besides, the convergence behaviour of the numerical scheme and the effects of particle spacing refinement and turbulence modelling on the simulation results are also investigated in further detail. The sensitivity of the computed wave breaking and overtopping on these issues is discussed and clarified. Copyright © 2005 John Wiley & Sons, Ltd.     

Interaction of isotropic turbulence with a normal shock wave

The evolution of incompressible and compressible isotropic 2-d turbulent fields interacting with a normal shock wave up to Mach numbers of 2.4 was investigated by means of direct numerical simulation using an ENO scheme. A comparison of statistics with linear analysis results is presented. Vorticity amplification in the DNS agrees well with the linear theory. Energy spectra are enhanced more in the small scales than in the large scales for incoming incompressible turbulence. The amplification rate for initially compressible turbulence is comparatively small.