Rayleigh waves with impedance boundary conditions in anisotropic solids

The paper is concerned with the propagation of Rayleigh waves in an elastic half-space with impedance boundary conditions. The half-space is assumed to be orthotropic and monoclinic with the symmetry plane _x3=0$x_3=0$. The main aim of the paper is to derive explicit secular equations of the wave. For the orthotropic case, the secular equation is obtained by employing the traditional approach. It is an irrational equation. From this equation, a new version of the secular equation for isotropic materials is derived. For the monoclinic case, the method of polarization vector is used for deriving the secular equation and it is an algebraic equation of eighth-order. When the impedance parameters vanish, this equation coincides with the secular equation of Rayleigh waves with traction-free boundary conditions.     

Explicit expression of polarization vector for surface waves,slip waves,Stoneley waves and interfacial slip waves in anisotropic elastic materials

We present explicit expression of the polarization vector for surface waves and slip waves in an anisotropic elastic half-space, and Stoneley waves and interfacial slip waves in two dissimilar anisotropic elastic half-spaces. An unexpected result is that, in the case of interfacial slip waves, the polarization vector for the material in the half-space _x2≥0$x_2\ge 0$ does not depend explicitly on the material property in the half-space _x2≤0$x_2\le 0$. It depends on the material property in the half-space _x2≤0$x_2\le 0$ implicitly through the interfacial slip wave speed υ$\upsilon$. The same is true for the polarization vector for the material in the half-space _x2≤0$x_2\le 0$.     

Squirt-Flow in Fluid-Saturated Porous Media: Propagation of Rayleigh Waves

Propagation of attenuated waves is studied in a squirt-flow model of porous solid permeated by two different pore regimes saturated with same viscous fluid. Presence of soft compliant microcracks embedded in the grains of stiff porous rock defines the double-porosity formation. Microcracks and pores respond differently to the compressional effect of a propagating wave, which induces the squirt-flow from microcracks to pores. Elastodynamics of constituent particles in porous aggregate is represented through a single-porosity formulation, which involves the frequency-dependent complex moduli. This formulation is deduced as a special case of double-porosity formation allowing the wave-induced flow of pore-fluid. This squirt-flow model of porous solid supports the attenuated propagation of two compressional waves and one shear wave. Superposition of these body waves, subject to stress-free surface, defines the propagation of Rayleigh wave. This wave is governed by a complex irrational dispersion equation, which is solved numerically after rationalising into an algebraic equation. For existence of Rayleigh wave, a complex solution of the dispersion equation should represent a leaky wave, which decays for propagation along any direction in the semi-infinite medium. A numerical example is solved to analyse the effects of squirt-flow on phase velocity, attenuation and polarisation of the Rayleigh waves, for different combinations of parameters. Numerical results suggest the existence of an additional (second) Rayleigh wave in the squirt-flow model of dissipative porous solids.     

The discrete Lamb problem: Elastic lattice waves in a block medium

We study the propagation of transient waves under the action of a vertical step point load on the surface of a half-space filled by a block medium. The block medium is modeled by a square lattice of masses connected by springs in the directions of the axes x,y$x,y$, and in the diagonal directions. The problem is solved by two methods. Analytically, we obtain asymptotic solutions in the vicinity of the Rayleigh wave at large time intervals. Numerically, we obtain a solution for any finite time interval. We compare these solutions with each other and with the solution to the Lamb problem for an elastic medium.     

Inversion for weak triclinic anisotropy from acoustic axes

Acoustic axes are directions in anisotropic elastic media, in which phase velocities of two or three plane waves (P$P$, S1$S1$ or S2$S2$ waves) coincide. Acoustic axes are important, because they can cause singularities in the field of polarization vectors and anomalies in the shape of the slowness surface. The maximum number of acoustic axes in triclinic anisotropy is 16, and their directions depend on anisotropy parameters in a complicate way. Under weak anisotropy approximation this dependence simplifies and the directions of acoustic axes can be used for the inversion for anisotropy parameters. The maximum acoustic axes under weak anisotropy is 16, the minimum number of acoustic axes is zero. In the inversion, we can retrieve 13 combinations of anisotropy parameters provided we use directions of 7 acoustic axes at least. Under weak anisotropy approximation, the directions of acoustic axes are insensitive to strength of anisotropy; hence we cannot invert for absolute values of weak anisotropy parameters, but only for their relative values. Numerical tests have shown that the inversion is applicable only to very weak anisotropy with strength of less than 5%, provided that the acoustic axes used in the inversion are determined with an accuracy of 0.1^°$0.1^{°}$ or better. In this case the inversion yields an average error for elastic parameters of less than 10%. In order to invert for the total set of 21 anisotropy parameters it is necessary to combine the measurements of the directions of the acoustic axes with measurements of other attributes of elastic waves in anisotropic media.     

Nonlinear interfacial waves in a circular cylindrical container subjected to a vertical excitation

Singular perturbation theory of two-time scale expansions was developed in inviscid fluids to investigate the motion of single interface standing wave in a two-layer liquid-filled circular cylindrical vessel, which is subjected to a vertical periodical oscillation. It is assumed that the fluid in the circular cylindrical vessel is inviscid, incompressible and the motion is irrotational, a nonlinear amplitude equation including cubic nonlinear and vertically forced terms, was derived by the method of expansion of two-time scales without taking the influence of surface tension into account. By numerical computation, it is shown that different patterns of interface standing wave can be excited for different driving frequency and amplitude. We found that the interface wave mode become more and more complex as increasing of upper to lower layer density ratio γ$\gamma$. The traits of the standing interface wave were proved theoretically. In addition, the dispersion relation and nonlinear amplitude equation obtained in this article can reduce to the known results for a single fluid when γ=0,_h2→_h1$\gamma =0,h_2\to h_1$.     

川藏公路地质环境与整治改建方案的思考

川藏公路由于地质环境复杂、建设标准低、后遗病害多,抗灾能力差,泥石流、滑坡、山崩、雪害、水毁等自然灾害频繁发生,公路阻车断道严重。国家投入巨资进行整治改建,并取得了明显的效果,但由于自然环境特殊、影响因素复杂,许多特大型、大型工程地质病害问题还没有可行、可靠的解决方案。本文通过分析川藏公路沿线的地质环境和灾害特点,总结历年整治改建和经验的教训,提出川藏公路建设的途径、可能达到的目标和应采用的原则。     

A mixture theory analysis for the surface-wave propagation in an unsaturated porous medium

The mixture theory is employed to the analysis of surface-wave propagation in a porous medium saturated by two compressible and viscous fluids (liquid and gas). A linear isothermal dynamic model is implemented which takes into account the interaction between the pore fluids and the solid phase of the porous material through viscous dissipation. In such unsaturated cases, the dispersion equations of Rayleigh and Love waves are derived respectively. Two situations for the Love waves are discussed in detail: (a) an elastic layer lying over an unsaturated porous half-space and (b) an unsaturated porous layer lying over an elastic half-space. The wave analysis indicates that, to the three compressional waves discovered in the unsaturated porous medium, there also correspond three Rayleigh wave modes (R1, R2, and R3 waves) propagating along its free surface. The numerical results demonstrate a significant dependence of wave velocities and attenuation coefficients of the Rayleigh and Love waves on the saturation degree, excitation frequency and intrinsic permeability. The cut-off frequency of the high order mode of Love waves is also found to be dependent on the saturation degree.     

An Eulerian–Eulerian model for gravity currents driven by inertial particles

In this work we introduce an Eulerian–Eulerian formulation for gravity currents driven by inertial particles. The model is based on the equilibrium Eulerian approach and on an asymptotic expansion of the two-phase flow equations. The final model consists of conservation equations for the continuum phase (carrier fluid), an algebraic equation for the disperse phase (particles) velocity that accounts for settling and inertial effects, and a transport equation for the disperse phase volume fraction. We present highly resolved two-dimensional (2D) simulations of the flow for a Reynolds number of Re=3450$\mathit{Re}=3450$ (this particular choice corresponds to a value of Grashof number of Gr=Re²/8=1.5×10⁶$\mathit{Gr}={\mathit{Re}}^2/8=1.5\times {10}^6$) in order to address the effect of particle inertia on flow features. The simulations capture physical aspects of two-phase flows, such as particle preferential concentration and particle migration down turbulence gradients (turbophoresis), which modify substantially the structure and dynamics of the flow. We observe the migration of particles from the core of Kelvin–Helmholtz vortices shed from the front of the current as well as their accumulation in the current head. This redistribution of particles in the current affects the propagation speed of the front, bottom shear stress distribution, deposition rate and sedimentation. This knowledge is helpful for the interpretation of the geologic record.     

Body waves in fractured porous media saturated by two immiscible newtonian fluids

A study of body waves in fractured porous media saturated by two fluids is presented. We show the existence of four compressional and one rotational waves. The first and third compressional waves are analogous to the fast and slow compressional waves in Biot's theory. The second compressional wave arises because of fractures, whereas the fourth compressional wave is associated with the pressure difference between the fluid phases in the porous blocks. The effects of fractures on the phase velocity and attenuation coefficient of body waves are numerically investigated for a fractured sandstone saturated by air and water phases. All compressional waves except the first compressional wave are diffusive-type waves, i.e., highly attenuated and do not exist at low frequencies.Now at Izmir Institute of Technology, Faculty of Engineering, Gaziosmanpasa Bulvari, No.16, Cankaya, Izmir, Turkey.     

Temporal behavior of laser induced elastic plate resonances

This paper investigates the dependence on Poisson’s ratio of local plate resonances in low attenuating materials. In our experiments, these resonances are generated by a pulse laser source and detected with a heterodyne interferometer measuring surface displacement normal to the plate. The laser impact induces a set of resonances that are dominated by Zero Group Velocity (ZGV) Lamb modes. For some Poisson’s ratio, thickness-shear resonances are also detected. These experiments confirm that the temporal decay of ZGV modes follows a t^−0.5$t^{-0.5}$ law and show that the temporal decay of the thickness resonances is much faster. Similar decays are obtained by numerical simulations achieved with a finite difference code. A simple model is proposed to describe the thickness resonances. It predicts that a thickness mode decays as t^−1.5$t^{-1.5}$ for large times and that the resonance amplitude is proportional to D^−1.5$D^{-1.5}$ where D$D$ is the curvature of the dispersion curve ω(k)$\omega \left(k\right)$ at k=0$k=0$. This curvature depends on the order of the mode and on the Poisson’s ratio, and it explains why some thickness resonances are well detected while others are not.     

Spatiotemporal breather in diffraction decreasing media

The similarity transformation between the (3+1$3+1$)-dimensional nonlinear Schrödinger equation with different distributed transverse diffraction and the standard nonlinear Schrödinger equation is found, and a spatiotemporal breather solution is given based on this transformation. The control for the evolutional behaviors of a spatiotemporal breather is discussed. Our results manifest that the relation between the maximum accumulated time _Tm$T_m$ and the accumulated time, _T0$T_0$, with the maximum amplitude, is the basis to realize the control and manipulation of propagation behaviors of breathers, such as fast and slow excitations, sustainment and restraint. These results are potentially useful for future experiments in the optical communications and Bose–Einstein condensations.     

Kinetic modeling of multiple scattering of elastic waves in heterogeneous anisotropic media

In this paper we develop a multiple scattering model for elastic waves in random anisotropic media. It relies on a kinetic approach of wave propagation phenomena pertaining to the situation whereby the wavelength is comparable to the correlation length of the weak random inhomogeneities—the so-called weak coupling limit. The waves are described in terms of their associated energy densities in the phase space position  ×$\times$  wave vector. They satisfy radiative transfer equations in this scaling, characterized by collision operators depending on the correlation structure of the heterogeneities. The derivation is based on a multi-scale asymptotic analysis using spatio-temporal Wigner transforms and their interpretation in terms of semiclassical operators, along the same lines as Bal (2005). The model accounts for all possible polarizations of waves in anisotropic elastic media and their interactions, as well as for the degeneracy directions of propagation when two phase speeds possibly coincide. Thus it embodies isotropic elasticity which was considered in several previous publications. Some particular anisotropic cases of engineering interest are derived in detail.     

Numerical simulation of 3D ductile cracks formation using recent improved Lode-dependent plasticity and damage models combined with remeshing

Damage to fracture transition has become a popular topic in the ductile fracture scientific community. Indeed, the transition from a damage continuous approach to a discontinuous fracture is not straightforward both from mechanical and numerical points of view. In the present study, a new improved Lode dependent phenomenological coupled damage model is used to investigate the ductile fracture in different mechanical tests. The remeshing and elements erosion techniques are employed to propagate the ductile cracks in 3D models using Forge® finite element code. This code is based on a mixed velocity–pressure formulation using the MINI element P1+/P1$P1+/P1$. In addition, the plasticity behavior is modeled by a Lode-dependent plasticity criterion. Applications to different mechanical tests at different loading configurations, using identified damage model parameters, show good agreement in terms of fracture prediction between experimental and numerical results.     

Wave propagation in transversely isotropic porous piezoelectric materials

Wave propagation in porous piezoelectric material (PPM), having crystal symmetry 6 mm, is studied analytically. Christoffel equation is derived for the propagation of plane harmonic waves in such a medium. The roots of this equation give four complex wave velocities which can propagate in such materials. The phase velocities of propagation and the attenuation quality factors of all these waves are described in terms of complex wave velocities. Phase velocities and attenuation of the waves in PPM depend on the phase direction. Numerical results are computed for the PPM BaTiO₃. The variation of phase velocity and attenuation quality factor with phase direction, porosity and the wave frequency is studied. The effects of anisotropy and piezoelectric coupling are also studied. The phase velocities of two quasi dilatational waves and one quasi shear waves get affected due to piezoelectric coupling while that of type 2 quasi shear wave remain unaffected. The phase velocities of all the four waves show non-dispersive behavior after certain critical high frequency. The phase velocity of all waves decreases with porosity while attenuation of respective waves increases with porosity of the medium. The characteristic curves, including slowness curves, velocity curves, and the attenuation curves, are also studied in this paper.     

Complex wavenumber Fourier analysis of the B-spline based finite element method

We present the results of one-dimensional complex wavenumber Fourier analysis of the B-spline variant of Finite Element Method (FEM). Generally, numerical results of elastic wave propagation in solids obtained by FEM are polluted by dispersion and attenuation. It was shown for the higher-order B-spline based FEM, that the optical modes did not occur in the case of infinite domains, unlike the higher-order Lagrangian and Hermitian finite elements, and also the dispersion errors are smaller. The paper’s main focus is on the wave propagation through B-spline multi-patch/segment discretization with the C⁰$C^0$ connection of B-spline segments and, chiefly, to the determining of dispersion and attenuation dependences. The numerical approach employed leads to substantial minimization of dispersion errors. Furthermore, the errors decrease in line with the increasing order of the B-spline elements/segments, with the local refinement, and also by the particular choice of the positions of control points through the optimizing procedure.     

A non-homogeneous constitutive model for human blood: Part II. Asymptotic solution for large Péclet numbers

In tube flow of healthy human blood the formed elements typically migrate away from vessel walls, leaving a plasma-rich, cell-depleted region there. In larger tubes (corresponding in size to arteries, for example) and at physiologically realistic flow rates, very thin wall boundary layers may develop which, nonetheless, have an impact upon the bulk flow properties. In this paper the non-homogeneous blood model of Moyers-Gonzalez et al. [M. Moyers-Gonzalez, R.G. Owens and J. Fang, A non-homogeneous constitutive model for human blood. Part I. Model derivation and steady flow, submitted for publication] is used in combination with a novel matched asymptotic method, to study the boundary layer behaviour of the steady tube flow of blood at high Péclet numbers Pe$Pe$ and in vessels of diameters corresponding to those of small arteries. A boundary layer thickness of O(Pe^−1/2)$O(P{\text{e}}^{-1/2})$ is predicted. In the absence of stress diffusion (the homogeneous case, with Pe=∞$Pe=\infty$) no cell migration takes place and the size and number density of red cell aggregates along the axis of symmetry remains constant at all flow rates. In the non-homogeneous case, however, even at very high values of Pe$Pe$, particles migrate, introducing a thin apparent slip layer next to the wall and affecting the aggregate distribution throughout the flow, even on the axis of symmetry.     

Chaining of weakly interacting particles suspended in viscoelastic fluids

An asymptotic theory based on multipole expansions is presented for multiparticle interactions in unbounded, weakly viscoelastic, creeping flows. The theory accounts for non-Newtonian sphere–sphere interactions that are of order O(De(a/R)²)$\text{O}(De{(a/R)}^2)$, where De is the Deborah number, a the sphere radius and R is the sphere–sphere separation. Analytic expressions are derived for the non-Newtonian correction to the multisphere mobility matrix for non-neutrally buoyant sedimenting spheres, and for neutrally buoyant spheres suspended in a shear flow. It is shown that these expressions give rise to particle chaining in simulations of interacting spherical particles.