Roll waves structure in two-layer Hele–Shaw flows

The mathematical model of inhomogeneous fluid motion in a Hele–Shaw cell is proposed. Based on this model the equations for describing two-layer flows and development of roll waves at the interface are derived. Conditions of roll waves existence are formulated in terms of Whitham criterion. Numerical calculations of the interface position are provided. It is shown that small perturbations of the interface in the inlet section of the channel lead to the roll waves for certain parameters of the flow. Two-parametric class of exact solutions corresponding to the roll waves regime is obtained. Diagrams of critical depths of roll waves development are constructed.     

激光热弹激励流-固泄漏Lamb波的有限元模拟

根据变分原理,得到热弹体运动方程和热传导方程相对应的有限元方程. 通过数值积分方法求解有限元方程,得到脉冲激光线源在水/铝、空气/铝这两种流-固界面上热弹激励的泄漏Lamb波瞬态波形. 计算结果表明,泄漏Lamb波不但存在于液-固界面,而且存在于气-固界面;和Lamb波相反,泄漏Lamb波的S_0模态是反对称的,而A_0模态是对称的;但由于这两种流-固界面的性质差异导致泄漏Lamb波的波形和幅度不同.      

Generation of free-surface gravity waves by an unsteady Stokeslet

The interaction of unsteady Stokeslets with the free surface of an initially quiescent incompressible fluid of infinite depth is investigated analytically for two- and three-dimensional cases. The disturbed flows are generated by an unsteady singular force moving perpendicularly downwards away from the surface. The analysis is based on the assumption that the motion satisfies the linearized unsteady Navier–Stokes equations with linear kinematic and dynamic boundary conditions. Firstly, the asymptotic representation for the transient free-surface waves due to an instantaneous Stokeslet is derived for a large time with a fixed distance-to-time ratio. As is well known, the corresponding inviscid waves predicted by the potential theory do not decay to zero as the time goes to infinity. In the present study, the transient waves predicted by the viscous theory eventually vanish due to the presence of viscosity, which is consistent with reality from the physical point of view. Secondly, the asymptotic solutions are obtained for the unsteady free-surface waves due to a harmonically oscillating Stokeslet. It is found that the unsteady waves can be decomposed into steady-state and transient responses. The steady state can be attained as time approaches infinity. It is shown that the viscosity of the fluid plays an important role in the evolution of the singularity-induced waves.     

Interfacial viscoelastic SH waves

Shear horizontal waves, in the form of transient perturbations, are considered at the interface between two different viscoelastic solids. The admissibility of these interfacial waves is studied via the asymptotic expansion of the Laplace transform of the viscoelastic kernel. The compatibility condition is reduced to a set of algebraic systems which can be solved iteratively to the desired order in the asymptotic expansion. Two classes of solutions are found which correspond to transient waves decaying away from the interface and attenuated along the propagation direction. Numerical examples are given to illustrate the results.     

Full field analysis of an anti-plane interface crack subjected to dynamic body forces

In this study, the transient full field response of an interface crack between two different media subjected to dynamic body force at one material is investigated. For time t < 0, the bimaterial medium is stress free and at rest. At t = 0, a concentrated anti-plane dynamic point loading is applied at the medium as shown in Fig. 1. The total wave field is due to the effect of this point loading and the scattering of the incident waves by the interface crack. An alternative methodology that is different from the conventional superposition method is used to construct the reflected, refracted and diffracted wave fields. A useful fundamental solution is proposed in this study and the full field solution is determined by superposition of the fundamental solution in the Laplace transform domain. The proposed fundamental problem is the problem of applying an exponentially distributed traction (in the Laplace transform domain) on the interfacial crack faces. The Cagniard–de Hoop method of Laplace inversion is used to obtain the transient solution in time domain. Exact transient closed form solutions for stresses and stress intensity factors are obtained. Numerical results for the time history of stresses and stress intensity factors during the transient process are discussed in detail.     

非均质流固耦合介质轴对称动力问题时域解

将地基视为流固两相介质并考虑其非均质成层特性，推导了多层地基动力问题时域解．文中首先建立了一组解耦的两相介质动力控制方程；而后利用Ｌａｐｌａｃｅ－Ｈａｎｋｅｌ变换推导了单层地基象空间初参数解答；再利用初参数法及传递矩阵技术导出了任意多层地基瞬态解的一般解析算式．本文获得的解答可方便地退化为现有理想弹性介质的解答     

Computation of interface wave motions by reciprocity considerations

The current theoretical study deals with computation of Stoneley waves along a solid–solid interface and Scholte waves (also called Scholte-Gogoladze) along a solid–liquid interface by reciprocity considerations. Closed-form solutions of the wave motions generated by a time-harmonic line load applied in two bonded elastic half-spaces of different material properties are derived in a simple manner. In order to perform direct applications of reciprocity theorems, we introduce in this article new expressions for the displacements of free interface waves. Reciprocity relations between an actual state, interface wave motion generated by a time-harmonic line load, and a virtual state, an appropriately chosen free wave traveling along the interface, are derived. Scattered amplitudes of Stoneley waves and Scholte waves due to the load are thus computed. To show application of the obtained results, scattering of Stoneley wave by a delamination at the interface is then studied.     

Generalized solutions of transient thermal shock problem with bounded boundaries

In this work, the generalized thermoelastic solutions with bounded boundaries for the transient shock problem are proposed by an asymptotic method. The governing equations are taken in the context of the generalized thermoelasticity with one relaxation time (L–S theory). The general solutions for any set of boundary conditions are obtained in the physical domain by the Laplace transform techniques. The corresponding asymptotic solutions for a thin plate with finite thickness, subjected to different sudden temperature rises in its two boundaries, are obtained by means of the limit theorem of Laplace transform. In the context of these asymptotic solutions, two specific problems with different boundary conditions have been conducted. The distributions of displacement, temperature and stresses, as well as the propagations, intersections and reflections of two elastic waves, named as thermoelastic wave and thermal wave separately, are obtained and plotted. These results are agreed with the results obtained in the existing literatures.     

Non-linear temporal dynamics of two-mode interactions of magnetized flow

This paper considers, in the frame work of the model of two superposed layers of viscous-potential incompressible magnetic fluids, the problem on formation of resonant waves of two modes on the interface between fluids that arisen as a result of second-harmonic resonance. The fluids moving with uniform velocities parallel to their interface are stressed by a tangential magnetic field. The analysis includes the linear, as well as the non-linear effects where the analytical solutions are constructed using the method of multiple scales, in both space and time, and hence the solvability conditions correspond to the uniform (convergent) solutions are obtained. The solvability conditions are then exploited to derive a more general system of non-linear partial differential equations with complex coefficients governing the amplitudes of the resonant waves. These equations are examined for solutions corresponding to sinusoidal wavetrains consequently different kinds of instabilities are demonstrated. The stability criterion in each case is derived and discussed both analytically and graphically.     

结构中微裂纹与超声波的混频非线性作用数值仿真研究

针对结构中微裂纹检测难题,本文对结构中微裂纹与超声波的混频非线性作用进行了数值仿真研究。基于经典非线性理论,得到了两列超声纵波相互作用产生混频效应的理论条件。通过有限元仿真,研究了两列纵波与微裂纹相互作用产生混频的条件,并分析了界面处静应力、摩擦系数和裂纹方向对混频效应的影响。研究发现,超声波与微裂纹相互作用产生混频非线性效应的发生条件仍符合经典非线性理论下的混频产生条件。裂纹界面处施加的静应力对差频横波幅值有明显影响;当施加静应力与无裂纹模型得到的最大应力值接近时,混频非线性效应最强;裂纹界面的摩擦系数对超声波的混频非线性效应影响较小;透射差频横波传播方向与经典非线性理论预测的理论差频分量方向基本一致,且几乎不受裂纹方向变化的影响,而反射差频横波的传播方向随裂纹方向的改变而有所不同。本文研究工作为微裂纹检出及方向识别做了有益探索。     

Static and Dynamic Green’s Functions in Peridynamics

We derive the static and dynamic Green’s functions for one-, two- and three-dimensional infinite domains within the formalism of peridynamics, making use of Fourier transforms and Laplace transforms. Noting that the one-dimensional and three-dimensional cases have been previously studied by other researchers, in this paper, we develop a method to obtain convergent solutions from the divergent integrals, so that the Green’s functions can be uniformly expressed as conventional solutions plus Dirac functions, and convergent nonlocal integrals. Thus, the Green’s functions for the two-dimensional domain are newly obtained, and those for the one and three dimensions are expressed in forms different from the previous expressions in the literature. We also prove that the peridynamic Green’s functions always degenerate into the corresponding classical counterparts of linear elasticity as the nonlocal length tends to zero. The static solutions for a single point load and the dynamic solutions for a time-dependent point load are analyzed. It is analytically shown that for static loading, the nonlocal effect is limited to the neighborhood of the loading point, and the displacement field far away from the loading point approaches the classical solution. For dynamic loading, due to peridynamic nonlinear dispersion relations, the propagation of waves given by the peridynamic solutions is dispersive. The Green’s functions may be used to solve other more complicated problems, and applied to systems that have long-range interactions between material points.     

A dynamical study of the evolution of pressure waves propagating through a semi-infinite region of homogeneous gas combustion subject to a time-harmonic signal at the boundary

In this paper, we study the evolution of pressure waves propagating in a region of gas combustion subject to a time-harmonic signal at the boundary. The problem is modeled by a non-linear, hyperbolic partial differential equation. Steady-state behavior is investigated using the perturbation method to ensure that enough time has passed for transient effects to have dissipated. The zeroth-order and first two approximations are obtained. Furthermore, the behavior of the following quantities is investigated, with particular attention paid to the low and high-frequency limits: the location of the peak of the first-order approximation, dispersion relations and phase speeds. Additionally, a maximum value of the perturbation parameter is determined ensuring boundedness of the solution. Next, approximate solutions are obtained in the low and high frequency limits and a comparison is made with the corresponding perturbation solution. Finally, the solution obtained from the perturbation method is compared with the long-time solution obtained by a non-standard finite-difference scheme.     

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.     

Localized transition waves in bistable-bond lattices

Discrete two-dimensional square- and triangular-cell lattices consisting of point particles connected by bistable bonds are considered. The bonds follow a trimeric piecewise linear force-elongation diagram. Initially, Hooke's law is valid as the first branch of the diagram; then, when the elongation reaches the critical value, the tensile force drops to the other. The latter branch can be parallel with the former (mathematically this case is simpler) or have a different inclination. For a prestressed lattice the dynamic transition is found analytically as a wave localized between two neighboring lines of the lattice particles. The transition wave itself and dissipation waves carrying energy away from the transition front are described. The conditions are determined which allow the transition wave to exist. The transition wave speed as a function of the prestress is found. It is also found that, for the case of the transition leading to an increased tangent modulus of the bond, there exists nondivergent tail waves exponentially localized in a vicinity of the transition line behind the transition front. The previously obtained solutions for crack dynamics in lattices appear now as a partial case corresponding to the second branch having zero resistance. At the same time, the lattice-with-a-moving-crack fundamental solutions are essentially used here in obtaining those for the localized transition waves in the bistable-bond lattices. Steady-state dynamic regimes in infinite elastic and viscoelastic lattices are studied analytically, while numerical simulations are used for the related transient regimes in the square-cell lattice. The numerical simulations confirm the existence of the single-line transition waves and reveal multiple-line waves. The analytical results are compared to the ones obtained for a continuous elastic model and for a related version of one-dimensional Frenkel-Kontorova model.     

Theoretical transient analysis and wave propagation of piezoelectric bi-materials

The transient response of piezoelectric bi-materials subjected to a dynamic anti-plane concentrated force or electric charge with perfectly bonded interface is examined in the present study. The problem is solved by using the Laplace transform method and the inverse Laplace transform is evaluated by means of Cagniard’s method. Exact transient full-field solutions of the contribution for each wave are expressed in explicit closed forms. The transient behavior of field quantities is examined in detail by numerical calculations. The existence condition of a propagating surface wave along the interface is discussed in detail. A surface wave can be guided by the interface of two semi-infinite materials in contact if one, at least, of these two materials is piezoelectric. The propagation velocity of the surface wave is explicitly expressed and is found to be less than the lower shear wave velocity of the two materials. The existence of the surface wave for piezoelectric–piezoelectric bi-materials is restricted to the situation that the shear waves of the two piezoelectric materials are very close. The possibility for the existence of the surface wave for piezoelectric–elastic bi-materials is much greater than that of the piezoelectric–piezoelectric bi-materials.     

Nonlinear Plane Waves in Finite Deformable Infinite Mooney Elastic Materials

For infinite perfectly elastic Mooney materials, nonlinear plane waves are examined in both two and three dimensions. In two dimensions, longitudinal and shear plane waves are examined, while in three dimensions, longitudinal and torsional plane waves are considered. These exact dynamic deformations, applying to the incompressible perfectly elastic Mooney material, can be viewed as extensions of the corresponding static deformations first derived by Adkins [1] and Klingbeil and Shield [2]. Furthermore, the Mooney strain-energy function is the most general material admitting nontrivial dynamic deformations of this type. For two dimensions the determination of plane wave solutions reduces to elementary mathematical analysis, while in three dimensions an integral of the governing system of highly nonlinear ordinary differential equations is determined. In the latter case, solutions corresponding to particular parameter values are shown graphically. This revised version was published online in June 2006 with corrections to the Cover Date.     

Fully three-dimensional direct numerical simulation of a plunging breakerSimulation numérique directe tridimensionnelle du déferlement plongeant

The scope of this paper is to show the results obtained for simulating three-dimensional breaking waves by solving the Navier–Stokes equations in air and water. The interface tracking is achieved by a Lax–Wendroff TVD scheme (Total Variation Diminishing), which is able to handle interface reconnections. We first present the equations and the numerical methods used in this work. We then proceed to the study of a three-dimensional plunging breaking wave, using initial conditions corresponding to unstable periodic sinusoidal waves of large amplitudes. We compare the results obtained for two simulations, a longshore depth perturbation has been introduced in the solution of the flow equations in order to see the transition from a two-dimensional velocity field to a fully three-dimensional one after plunging. Breaking processes including overturning, splash-up and breaking induced vortex-like motion beneath the surface are presented and discussed. To cite this article: P. Lubin et al., C. R. Mecanique 331 (2003).     

A note on sinusoidal motion of a viscoelastic non-Newtonian fluid

In this note, the exact solutions of velocity field and associated shear stress corresponding to the flow of second-grade fluid in a cylindrical pipe, subject to a sinusoidal shear stress, are determined by means of Laplace and finite Hankel transform. These solutions are written as sum of steady-state and transient solutions, and they satisfy governing equations and all imposed initial and boundary conditions. The corresponding solutions for the Newtonian fluid, performing the same motion, can be obtained from our general solutions. At the end of this note, the effects of different parameters are presented and discussed by showing flow profiles graphically.     

Analysis of transient wave scattering from an inclusion with an unilateral smoothly contact interface by BEM

Summary A 2D time-domain Boundary Element Method (BEM) is applied to solve the problem of transient scattering of plane waves by an inclusion with a unilateral smooth contact interface. The incident wave is assumed strong enough so that localized separations take place along the interface. The present problem is indeed a nonlinear boundary value problem since the mixed boundary conditions involve unknown intervals (separation and contact regions). In order to determine the unknown intervals, an iterative technique is developed. As an example, we consider the scattering of plane waves by the cross section of a circular cylinder embedded in an infinite solid. Numerical results for the near field solutions are presented. The distortion of the response waves and the variation of the interface states are discussed. The financial support by the China National Natural Science Foundation under Grant No. 19872001 and No. 59878004 is gratefully acknowledged. The second author is also grateful to the support of the National Science Fund for Distinguished Young Scholars under Grant No. 10025211.     

A study of transient elastic wave propagation in a bimaterial modeling the thorax

The transient response of an elastic bimaterial, made out of a “hard” medium and a “soft” medium, welded at a plane interface, have been investigated by using an integral transform technique that permits isolation of the pressure and shear waves contributions to the wave-field. The method, often referred to as the generalized ray/Cagniard-de Hoop method (GR/CdH), is briefly presented. The wave motion is generated alternatively by a buried point source of strain rate and by a point force perpendicular to the free surface of the bimaterial. New simplified solutions are derived for points located on the axis perpendicular to the interfaces and passing through the source. Owing to the formalism, an approximation to the strain energy is shown to be readily obtained. The numerical schemes for the implementation of the exact three-dimensional GR/CdH are presented. Numerical examples are concerned with the propagation of an impact wave in the thorax modeled as a bimaterial (thoracic wall–lung). The effects of the weak coupling between the thoracic wall and the lung are investigated. The distributions of transient strain energies, respectively carried by the pressure and the shear waves in the media representing the lung, are plotted.