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$.     

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.     

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$.     

The two-dimensional free-space Green’s function and its derivatives in a tangent-normal system

The two-dimensional free-space Green’s function, G⁽²⁾$G^{\left(2\right)}$, and its derivatives, are used extensively in the formulation of scattering and diffraction problems through its presence in single- and double-layer potentials, and their use in integral equations. The vast majority of the results from elementary classical mathematical physics for G⁽²⁾$G^{\left(2\right)}$ is based on Cartesian coordinate-space, either directly as a Hankel function in coordinate-space or through a transform, such as the Weyl transform, also based on Cartesian coordinate-space. However, if the geometry of the problem is not Cartesian, for example in scattering from a rough surface, there are difficulties in using a transform representation for G⁽²⁾$G^{\left(2\right)}$ which depends on Cartesian geometry, as the standard Weyl transform does. Here we formulate transform-space representations using a tangent-normal coordinate system. The result for G⁽²⁾$G^{\left(2\right)}$ is a new Weyl-type tangent-normal transform representation from which the results for the vector derivatives of the single-layer potential, the double-layer potential, and the vector derivatives of the double-layer potential follow quite simply. The latter three results can be expressed in terms of two new spectral functions in tangent-normal space, _S1$S_1$ and _S2$S_2$. The overall results are new representations for G⁽²⁾$G^{\left(2\right)}$ and its derivatives which may be useful in integral equation formulations of scattering problems for non-Cartesian geometries.     

Wake response of a stationary finite-sized particle in a turbulent channel flow

Here we consider the effect of a finite-sized stationary particle in a channel flow of modest turbulence at _Reτ=178.12${\mathit{Re}}_{\tau }=178.12$. The size of particle is varied such that the particle Reynolds number ranges from about 40 to 450. The location of the particle is chosen to be either in the buffer layer (_yp+=17.81)$(y_{p+}=17.81)$ or at the channel center. Fully resolved direct numerical simulations of the turbulent channel flow around the particles is performed. Here the ambient turbulence intensity relative to the mean velocity seen by the particle is large (I=23.16%)$(I=23.16\%)$ in the buffer region, while it is substantially lower (I=4.09%)$(I=4.09\%)$ at the channel center. We present results on turbulence modulation due to the particle in terms of wake dynamics and vortex shedding.     

Lamb wave mode decomposition for structural health monitoring

Lamb waves propagate over large distances in plate-like thin structures and they have received great attention in the structural health monitoring (SHM) field as an efficient means to inspect a large area of a structure by using only a small number of sensors. The times-of-flight of the Lamb wave modes are useful for detecting damage generated in a structure. However, due to the dispersive and multi-mode nature of Lamb waves, it is very challenging to decompose Lamb wave modes into symmetric and anti-symmetric modes for potential applications to structural health monitoring. Thus, we propose an efficient Lamb wave mode decomposition method based on two fundamental rules: the group velocity ratio rule and the mode amplitude ratio rule. The group velocity ratio rule means that the ratio of the group velocities of _A0$A_0$ and _S0$S_0$ modes must be constant. The mode amplitude ratio rule means that the ratio of the magnitudes of _A0$A_0$ and _S0$S_0$ modes in a measured response signal must be always greater than one once the center frequency of the input signal is determined, such that the magnitude of the _A0$A_0$ mode in the excited signal is larger than that of the _S0$S_0$ mode, and vice versa. The proposed method is verified through experiments conducted for a plate specimen.     

Supersonic responses induced by point load moving steadily on an anisotropic half-plane

Supersonic responses of an anisotropic half-plane solid induced by a point load moving steadily on the half-plane boundary are investigated. Analytic expressions for the responses of the displacements and stresses for field points either inside or on the surface of the half-plane solid are given for general anisotropic materials. For the special cases of monoclinic materials with symmetry plane at $x_3=0$ and orthotropic materials, the supersonic as well as subsonic responses of the displacements and stresses are further expressed explicitly in terms of elastic stiffnesses. Responses for the case of isotropic materials known in the literature are recoverable from present results.     

Influence of grading on the undrained behavior of granular materials

This paper aims at investigating the influence of grading on the undrained behavior of granular materials. Series of undrained triaxial tests were carried out with two different materials, glass balls and Hostun sand. For each material, samples with different gradings and similar relative densities were prepared. Experimental results show that the undrained shear strength decreases when the coefficient of uniformity C_u=d₆₀/d₁₀$C_{\mathrm{u}}=d_{60}/d_{10}$ increases from 1.1 to 20. The conditions of instability for the selected granular materials were also analyzed, based on the sign of the second-order work during undrained triaxial loading. The results demonstrate a significant influence of the grading: increasing the coefficient of uniformity heightens the potential of static liquefaction and the materials become more unstable. Furthermore, numerical tests using the three-dimensional discrete element method (DEM) were conducted on assemblies of spheres. The DEM inter-particle parameters were derived from the experimental test results on glass balls. The DEM simulations showed similar behaviors compared to experimental results and confirmed the influence of the grain size distribution on the stress–strain relationship and instability phenomena.     

Extraction of bulk wave characteristics from a pulsed ultrasonic polar scan

The pulsed ultrasonic polar scan (P-UPS) technique provides a map with characteristic contours reflecting the critical bulk wave angles, which directly relate to the in-plane elastic properties of the insonified material spot. Besides these contours, additional information is captured in a P-UPS experiment which has particular importance for the inspection of anisotropic materials. By extensive post-processing on a single P-UPS experiment, we successfully extracted the local directional (i) critical bulk wave angles, (ii) phase velocity profiles, (iii) slowness curves, (iv) walk-off angles and (v) energy velocity profiles, for the three different polarization states of bulk waves. The procedure is demonstrated and discussed for an isotropic aluminum sample as well as an autoclave manufactured orthotropic _[0]8${\left[0\right]}_8$ carbon/epoxy laminate. In addition, a numerical and experimental investigation of the role of the immersion liquid as a lens for the P-UPS result is performed. This lensing effect permits to zoom in/out on the P-UPS fingerprint, hence broadening the class of materials which can be inspected by the P-UPS technique.     

Study of the rotational diffusivity coefficient of fibres in planar contracting flows with varying turbulence levels

The Fokker–Planck equation is solved by describing the evolution of a 3D fibre orientation state along a planar contraction. A constant value of the effective rotational diffusion coefficient was determined for four different turbulent flow cases in planar contractions, reported experimentally in the literature. Two hypotheses for the non-dimensional rotational diffusivity are presented, each based on two different turbulent time scales, i.e. the Kolmogorov time scales and the time scale associated with large energy bearing eddies. These hypotheses are dependent on either the Reynolds number, based on the Taylor micro-scale, and/or a non-dimensional fibre length. The hypothesis, based on the assumption of long fibres, _Lf/η?25$L_{\text{f}}/\eta ?25$, compared to the Kolmogorov scale and in the limit of large _Reλ${\mathit{Re}}_{\lambda }$ seems to capture the basic trends presented in the literature. This hypothesis has also the feature of predicting effects of varying fibre length within certain limits. Accordingly, by modeling the variation of turbulent quantities along the contraction in a CFD analysis, local values of rotational diffusivity can be evaluated with the mentioned hypothesis, based on either Kolmogorov time scale or Eulerian integral time scale.     

A new secular equation for slip waves along the interface of two dissimilar anisotropic elastic half-spaces in sliding contact

We consider wave propagation along the interface of two dissimilar anisotropic elastic half-spaces that are in sliding contact. A new secular equation is obtained that covers all special cases in one equation. One special case is when a Rayleigh wave (called the R$R$-wave) can propagate in both half-spaces with the same wave speed. Another special case is when a slip wave (called the S$S$-wave) can propagate in each of the half-spaces with the same wave speed. If a Rayleigh wave and a slip wave can propagate in one of the half-spaces it is called the RS-wave. In this case an interfacial slip wave exists in which the other half-space is at rest unless an RS-wave can also propagate in the other half-space. The results for general anisotropic elastic materials are applied to orthotropic materials.     

A penny-shaped crack in transversely isotropic magneto-electro-thermo-elastic medium subjected to uniform symmetric heat flux

The problem of a penny-shaped crack subjected to symmetric uniform heat flux in an infinite transversely isotropic magneto-electro-thermo-elastic medium is investigated. The exact solution in the full space is in terms of line integrals and the exact solution in the crack plane also is obtained. Although we start our derivations with magneto-electro-thermo-elastic, the solution presented in this paper is also applicable for linear transversely isotropic thermopiezoelectric, thermomagnetoelastic,thermoelastic materials (see Appendix E). The solution in the crack plane, which shows a great agreement with the solution for a transversely isotropic medium obtained by Tsai (1983), indicates that _σx,_σy,_Dx,_Dy,_Bx${\sigma }_x\text{,}{\sigma }_y\text{,}{D}_x\text{,}{D}_y\text{,}{B}_x$, and _By$B_y$ along the crack rim are of the same singularity of the normal stress or its equivalent quantities. To illustrate how the applied symmetric heat fluxes affect the whole fields, a numerical example is also given.     

Molecular modeling of cracks at interfaces in nanoceramic composites

Toughness in Ceramic Matrix Composites (CMCs) is achieved if crack deflection can occur at the fiber/matrix interface, preventing crack penetration into the fiber and enabling energy-dissipating fiber pullout. To investigate toughening in nanoscale CMCs, direct atomistic models are used to study how matrix cracks behave as a function of the degree of interfacial bonding/sliding, as controlled by the density of C interstitial atoms, at the interface between carbon nanotubes (CNTs) and a diamond matrix. Under all interface conditions studied, incident matrix cracks do not penetrate into the nanotube. Under increased loading, weaker interfaces fail in shear while stronger interfaces do not fail and, instead, the CNT fails once the stress on the CNT reaches its tensile strength. An analytic shear lag model captures all of the micromechanical details as a function of loading and material parameters. Interface deflection versus fiber penetration is found to depend on the relative bond strengths of the interface and the CNT, with CNT failure occurring well below the prediction of the toughness-based continuum He–Hutchinson model. The shear lag model, in contrast, predicts the CNT failure point and shows that the nanoscale embrittlement transition occurs at an interface shear strength scaling as _τs~_εf,CNTσCNT${\tau }_s~{\epsilon }_{f,CNT}{\sigma }_{CNT}$ rather than _τs~_σCNT${\tau }_s~{\sigma }_{CNT}$ typically prevailing for micron scale composites, where _εf,CNT${\epsilon }_{f,CNT}$ and _σCNT${\sigma }_{CNT}$ are the CNT failure strain and stress, respectively. Interface bonding also lowers the effective fracture strength in SWCNTs, due to formation of defects, but does not play a role in DWCNTs having interwall coupling, which are weaker than SWCNTs but less prone to damage in the outerwall.     

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.     

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.     

Active manipulation of a particle-laden jet

The interaction of a particle-laden jet (_ReUe=6600)$({\mathit{Re}}_{U_{\text{e}}}=6600)$ with a single synthetic jet or a continuous control jet located upstream of the main jet exit (i.e., within the main jet nozzle) was examined experimentally using PIV and PTV. A reduction technique was used to calculate 3-D flow fields from multiple 2-D measurement planes to study the complex 3-D interactions. The synthetic jet was shown to influence the particles both directly and indirectly through the manipulation of the carrier fluid’s drag force on the particles. The synthetic jet impulse directly vectors the particles away from the synthetic jet, while the formation of large vortical structures indirectly affects the particles, spreading throughout the measurement domain. By comparison, a continuous control jet only vectors the particles away from itself. The lowest Stokes number particles respond similarly to the carrier fluid, while higher Stokes number particles are less responsive to the control and only follow the strong vortical structures (i.e., higher circulation), which suggests that the preferential concentration concept depends on both the Stokes number as well as the strength of the coherent structures.