Perturbation Formula for Phase Velocity of Rayleigh Waves in Prestressed Anisotropic Media

Herein we consider Rayleigh waves propagating along the free surface of a macroscopically homogeneous, anisotropic, prestressed half-space. We adopt the formulation of linear elasticity with initial stress and assume that the deviation of the prestressed anisotropic medium from a comparative ‘unperturbed’, unstressed and isotropic state, as formally caused by the initial stress and by the anisotropic part of the incremental elasticity tensor, be small. No assumption, however, is made on the material anisotropy of the incremental elasticity tensor. With the help of the Stroh formalism, we derive a first-order perturbation formula for the shift of phase velocity of Rayleigh waves from its comparative isotropic value. Our perturbation formula does not agree totally with that which was derived some years ago by Delsanto and Clark, and we provide another argument as further support for our version of the formula. According to our first-order formula, the anisotropy-induced velocity shifts of Rayleigh waves, taken in totality of all propagation directions on the free surface, carry information only on 13 elastic constants of the anisotropic part of the incremental elasticity tensor. The remaining eight elastic constants are those which would become zero if were monoclinic with the two-fold symmetry axis normal to the free surface of the material half-space in question.     

Perturbation Formulas for Polarization Ratio and Phase Shift of Rayleigh Waves in Prestressed Anisotropic Media

This paper completes an earlier study (Tanuma and Man, Journal of Elasticity, 85, 21–37, 2006) where we derive a first-order perturbation formula for the phase velocity of Rayleigh waves that propagate along the free surface of a macroscopically homogeneous, anisotropic, prestressed half-space. We adopt the formulation of linear elasticity with initial stress and assume that the deviation of the prestressed anisotropic medium from a suitably-chosen, comparative, unstressed and isotropic state be small. No assumption, however, is made on the material anisotropy of the incremental elasticity tensor. With the help of the Stroh formalism, here we derive first-order perturbation formulas for the changes in polarization ratio and phase shift of Rayleigh waves from their respective comparative isotropic value. Examples are given, which show that the perturbation formulas for phase velocity and polarization ratio can serve as a starting point for investigations on the possible advantages of using Rayleigh-wave polarization, as compared with using wave speed, for acoustoelastic measurement of stress.      

Propagation of Waves on the Free Surface of a Two-Phase Mixture

The problem of the propagation of waves on the free surface of layer of a two-phase mixture is considered. An analytic solution in the form of damped steady waves is found in the linear approximation. The dispersion relation, an expression for the decay rate, and the shape of the free surface are determined. The effect of the dispersed phase on the wave velocity is found.     

Focusing and Scattering of a Nonstationary Wave by the Free Surface of an Anisotropic Elastic Medium

The interaction of a plane discontinuous wave with the free curvilinear surface of an anisotropic elastic medium is studied. The ray-path method is used to describe the wavefront. The nonlinear Snell's equations are solved by combining Newton's method and the parameter continuation algorithm. The cases of focusing and scattering of wave rays are studied     

Thermal Capacity Effect on Transient Free Convection Adjacent to a Vertical Surface in a Porous Medium

In this paper we analyse how the presence of the thermal capacity of a vertical flat plate of finite thickness, which is embedded in a porous medium affects the transient free convection boundary-layer flow. At the time t = 0, the plate is suddenly loaded internally with a constant heat flux rate q, so that a transient boundary-layer flow is initiated adjacent to the plate. Initially, the transient effects due to the imposition of the uniform heat flux rate at the plate are confined to a thin fluid region near to the surface and are described by a small time solution. These effects continue to penetrate outwards and eventually evolve into a new steady state flow. Analytical solutions have been derived for these transient (small time) and steady state (large time) flow regimes, which are then matched by a numerical solution of the full boundary-layer equations. It has been found that the non-dimensional fluid temperature (or fluid velocity) profiles are reduced when the thermal capacity effects, described by a parameter Q ^*, are reduced. For small values of Q ^*, the approach of these profiles to their steady state values is monotonic. However, for large values of Q ^*, the temperature profiles are observed to locally exceed (pass through a maximum value) the final steady state values at certain distances from the plate. In general, the maxima in the temperature profiles increase in size as Q ^* increases and the time taken to approach the steady state solutions increases significantly.     

Dispersion of Small Amplitude Waves in a Pre-Stressed, Compressible Elastic Plate

The dispersion of harmonic waves, propagating along a principal direction in a pre-stressed, compressible elastic plate, is investigated in respect of the most general isotropic strain-energy function. Different cases, dependent on the choice of material parameters and pre-stress, are analysed. A complete long and short wave asymptotic analysis is carried out, with the approximations obtained giving phase speed (and frequency) as explicit functions of wave and mode number. Various wave fronts, both associated with the short wave limit of harmonics and arising through the combination of harmonics in a narrow wave speed region, are discussed. It is mentioned that the case of high compressibility is of particular interest. In contrast with the classical (un-strained) case, the longitudinal body wave speed may be less than the corresponding shear wave speed. In consequence, the short wave limit of all harmonics may be the appropriate longitudinal wave speed; contrasting with the classical case for which this limit is necessarily associated with a shear wave front. A further possible short wave limit is also shown to exist for which the associated wave normal has a component in the direction normal to the plate. Particularly novel numerical results are presented when the longitudinal and shear wave speeds are equal. The analysis is illustrated by numerical calculations for various strain-energy functions.     

Nonlinear Capillary-Gravity Waves on the Charged Surface of an Ideal Fluid

A second-order asymptotic expression for the profile of a capillary-gravity wave traveling over the charged surface of an ideal incompressible fluid is calculated analytically. Two types of steady-state profiles of nonlinear periodic capillary-gravity waves are found. For a certain fixed dimensionless surface charge the shape of the tops of the nonlinear waves changes: from blunt to pointed for short waves and from pointed to blunt for long waves.     

Tollmien-Schlichting Waves in a Hypersonic Boundary Layer Flow in the Neighborhood of a Cooled Surface

A nonlinear time-dependent model of the development of longwave perturbations in a hypersonic boundary layer flow in the neighborhood of a cooled surface is constructed. The pressure in the flow is assumed to be induced the combined variation of the thicknesses of the near-wall and main parts of the boundary layer. Numerical and analytic solutions are obtained in the linear approximation. It is shown that if the main part of the boundary layer is subsonic as a whole, its action reduces the perturbation damping upstream and the perturbation growth downstream, while a supersonic, as a whole, main part of the boundary layer creates the opposite effects. An analysis of the solutions obtained makes it possible to conclude that the asymptotic model proposed can describe the three-dimensional instability of the Tollmien-Schlichting waves.     

Equivalence Theorems on the Propagation of Small-Amplitude Waves in Prestressed Linearly Elastic Materials with Internal Constraints

By extending the procedure of linearization for constrained elastic materials in the papers by Marlow and Chadwick et al., we set up a linearized theory of constrained materials with initial stress (not necessarily based on a nonlinear theory). The conditions of propagation are characterized for small-displacement waves that may be either of discontinuity type of any given order or, in the homogeneous case, plane progressive. We see that, just as in the unconstrained case, the laws of propagation of discontinuity waves are the same as those of progressive waves. Waves are classified as mixed, kinematic, or ghost. Then we prove that the analogues of Truesdell"s two equivalence theorems on wave propagation in finite elasticity hold for each type of wave. This revised version was published online in August 2006 with corrections to the Cover Date.     

Symmetrical Representation of Stresses in the Stroh Formalism and its Application to a Dislocation and a Dislocation Dipole in an Anissotropic Elastic Medium

Symmetrical stress representation in the Stroh formalism for anisotropic elastic bodies is introduced and the range of its applicability is analysed. By making use of this stress representation new formulae for influence functions giving stresses in an infinite anisotropic medium subjected to a straight dislocation and a straight dislocation dipole are derived. The advantage of the new formulae is that they explicitly show the symmetrical structure of these influence functions not referred to previously. Relations of these influence functions to influence functions giving stresses and Airy stress function due to a straight wedge disclination, whose explicit expressions are also introduced, are derived. Application of these results in computation of stresses by the hypersingular and regularized Somigliana stress identities is discussed. This revised version was published online in August 2006 with corrections to the Cover Date.     

Localization of wave motions in a fluid-filled cylinder

The properties of harmonic surface waves in a fluid-filled cylinder made of a compliant material are studied. The wave motions are described by a complete system of dynamic equations of elasticity and the equation of motion of a perfect compressible fluid. An asymptotic analysis of the dispersion equation for large wave numbers and a qualitative analysis of the dispersion spectrum show that there are two surface waves in this waveguide system. The first normal wave forms a Stoneley wave on the inside surface with increase in the wave number. The second normal wave forms a Rayleigh wave on the outside surface. The phase velocities of all the other waves tend to the velocity of the shear wave in the cylinder material. The dispersion, kinematic, and energy characteristics of surface waves are analyzed. It is established how the wave localization processes differ in hard and compliant materials of the cylinder __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 4, pp. 72–86, April 2008.     

Investigation of the Dynamic Behavior of a Griffith Crack in a Piezoelectric Material Strip Subjected to the Harmonic Elastic Anti-plane Shear Waves by Use of the Non-local Theory

In this paper, the dynamic behavior of a Griffith crack in a piezoelectric material strip subjected to the harmonic anti-plane shear waves is investigated by use of the non-local theory for impermeable crack surface conditions. To overcome the mathematical difficulties, a one-dimensional non-local kernel is used instead of a two-dimensional one for the anti-plane dynamic problem to obtain the stress and the electric displacement near at the crack tip. By means of the Fourier transform, the problem can be solved with the help of two pairs of dual integral equations. These equations are solved using the Schmidt method. Contrary to the classical solution, it is found that no stress and electric displacement singularity is present near the crack tip. The non-local dynamic elastic solutions yield a finite hoop stress near the crack tip, thus allowing for a fracture criterion based on the maximum dynamic stress hypothesis. The finite hoop stress at the crack tip depends on the crack length, the thickness of the strip, the circular frequency of incident wave and the lattice parameter.     

Effect of electrohydrodynamic interaction on the stability of a flat plate laminar boundary layer

The stability of a unipolarly charged electrohydrodynamic boundary layer on a flat dielectric plate along which an electric current flows between electrodes located on the plate is investigated within the framework of the linear theory. The solution of the steady-state problem is obtained on the basis of methods developed earlier for conditions typical of aerodynamical experiments and various electric currents and electrode voltages. The effect of the interaction between perturbations of the electric and hydrodynamic flow parameters on the flow stability is estimated within the framework of the locally homogeneous approximation. This effect turns out to be insignificant under the conditions considered. It is shown that steady-state electrohydrodynamic action on the main flow makes it possible to obtain “accelerating” velocity profiles with increased absolute values of the second derivative in the transverse direction. This ensures a significant increase in the critical Reynolds numbers of loss of stability and a narrowing of the growing perturbation wavenumber range.     

Dynamics of a prestressed stiff layer on an elastic half space: filtering and band gap characteristics of periodic structural models derived from long-wave asymptotics

Anti-plane and plane-strain, time-harmonic, small-amplitude vibrations of an elastic layer on an elastic half space are considered, superimposed upon a state of finite, uniform stress and strain. A (compressible) elastic material is considered, orthotropic with orthotropy axes aligned parallel and orthogonal both to the layer and the prestress principal directions. A non-uniform mass density is assumed in the layer. A formal long-wave asymptotic solution is derived under the assumptions of high contrast between the stiffnesses of the layer and the half space and between certain prestress components and the current elastic shear modulus.It is shown that (i) the layer asymptotically behaves as a beam subject to transversal and axial vibrations; (ii) the response of the half space can be found in a closed-form, under the assumption of plane wave motion (which becomes consistent when the density of the layer is uniform), otherwise it is represented by a hypersingular integral equation; (iii) if the nonlocality introduced by the hypersingular integral equation is restricted to an influence area of finite extent, the integral can be analytically approximated, so that a Winkler-type spring model representing the half space is rigorously derived. For uniform density of the layer, the constants defining the spring model are given as functions of the prestress and anisotropy parameters of the half space; and, finally, (iv) the asymptotic solution provides new analytical expressions for incremental displacement of the layer, which, compared to the exact numerical solution (also included), are shown to perform quite well, even for values of parameters much beyond the limits imposed by the asymptotic analysis.The asymptotic analysis allows us to explore, for the first time, dynamic properties of a periodic layer bonded to an elastic half space and subject to a uniform prestress state. We find that the system exhibits band gaps (ranges of forbidden frequencies) and that the prestress can be used as a parameter tuning the filtering properties of the structure, an effect which may have important consequences in the design of resonant devices.     

Oscillation of a Punch Moving on the Free Surface of an Elastic Half Space

The contact problem concerning oscillation of a circular rigid punch, moving uniformly at sub-Rayleigh speed along the surface of an elastic half space, is investigated using a three-dimensional formulation. Slow motion of the punch is considered, which implies that the characteristic time for the external loading is much larger than the time interval necessary for shear waves to propagate across the punch. An asymptotic solution for the vertical oscillation of the punch is given. It is shown that the vertical displacement of the punch can approximately be described by the equation of dynamics for a system of one degree of freedom with viscous friction. The dependence of the coefficients for effective viscosity and stiffness, occurring in this equation, on the speed of the punch and Poisson ratio of the half space, is investigated. The solution for the non-stationary problem concerning a suddenly applied moving point load is also obtained, correcting and extending the result known so far. Mathematics Subject Classifications (2000) 74H10, 74J05, 74M15.     

Separation Impact on a Plate Floating on the Surface of an Ideal Incompressible Fluid in a Bounded Tank

The paper studies the planar problem of separation impact on a plate floating on the surface of an ideal incompressible fluid in a bounded tank. The problem is solved using an asymptotic method under the assumption that the immovable rigid walls of the tank are at a large distance from the plate. It is concluded that the tank walls of arbitrary shape have an ambiguous effect on the fluid particle separation zone formed on the plate surface is revealed. Examples of solutions are given.     

Nonlinear Stability Analysis of Thin Viscoelastic Film Flowing Down on the Inner Surface of a Rotating Vertical Cylinder

This paper investigates the stability of thin viscoelastic liquid film flowing down on the inner surface of a rotating vertical cylinder by means of the long wave perturbation. After proving the insufficiency of the linear model in characterization of certain flow behaviors, a generalized nonlinear kinematic model is then derived to represent the physical system. This model is solved through the following procedure. In the first step, the normal mode method is used to characterize the linear behaviors. The amplitude growth rates and the threshold conditions are characterized subsequently and summarized as the by-products of the linear solutions. In the second step, a nonlinear film flow model is solved by using the method of multiple scales to characterize flow behaviors at various states of sub-critical stability, sub-critical instability, supercritical stability, and supercritical explosion. The modeling results indicate that with the increase in the rotation speed Ω and the radius of cylinder R, the film flow system will be more stable.     

Determination of the molecular weight distribution of linear polymers by inversion of a blending law on complex viscosities

The method presented in this paper allows to calculate the molecular weight distribution (MWD) of linear homopolymer melts from the complex shear modulus data measured in a wide frequency domain. An empirical blending law on complex viscosities is first developed; as a consequence, the variations of the storage and loss modulus as a function of MWD are presented. This simulation demonstrates also the role of the shape of the MWD itself, and shows that one should not postulate a priori the shape of the MWD. An efficient numerical approach based on a Tikhonov regularization method with constraint is used to solve this ill-posed problem; the MWD is hence derived without any assumption on its shape. This method is first applied on simulated data to prove its numerical efficiency. Then the inversion method is applied on complex moduli data of various commercial polymers (polypropylene, polyethylene and polystyrene) and on an artificial mixture of polystyrene that have been presented in the literature. For amorphous polymers, the coupling of the terminal relaxation domains with the transition region at higher frequency leads to errors in the low molecular weight tail: one way to solve this problem is to cut off the experimental data at the high frequencies. This general method needs only a few physical parameters, namely the scaling law for the Newtonian viscosity η₀=f(M _w ) and the plateau modulus G _N ⁰, and leads to reasonable results with respect to the simplicity of the viscoelastic model used. Received: 27 October 1997 Accepted: 24 February 1998     

Loss Reduction by Means of Two-Dimensional Roughness Elements on the Suction Surface of a Linear Turbine Rotor Cascade

This paper reports the results of experimental investigations carried out to reduce pressure losses by means of two-dimensional roughness elements (in the form of stainless steel tubes of different diameters). The roughness elements are fixed at various axial stations on the suction surface of 120° turning, 175 mm chord, impulse turbine rotor blades. Flow measurements are carried out at the exit of the cascade at five axial stations, using a five-hole probe operating in non-nulling mode. In addition, the blade surface static pressure distribution is measured. The data from the five-hole probe measurements are used to calculate pressure, velocity and flow angle distributions at the cascade exit and these results are used to calculate mass averaged values and integral parameters such as wake half-width, loss coefficient, etc. The static pressure distribution is altered very little except near the roughness element. The lift coefficient remains almost constant for all configurations and the drag coefficient is reduced for some configurations. The non-dimensional total pressure defects in the wake for all configurations followed Gaussian distribution. A two-dimensional roughness element of 0.6 mm diameter placed at 0.65 chord on the suction surface showed an appreciable reduction in pressure losses.