Macroscopic and microscopic examination of the relationship between crack velocity and path and Rayleigh surface wave speed in single crystal silicon

The speed of Rayleigh surface waves, denoted C_R, is the accepted upper limit for Mode I crack velocity in monolithic solids. In the current contribution, we discuss several critical issues associated with the velocity of Rayleigh surface waves and crack velocity in single crystal (SC) brittle solids, and the global and local influence of C_R on crack path selection in particular.Recent cleavage experiments in SC silicon showed that crack velocity at certain cleavage planes and crystallographic orientations cannot exceed a small fraction of C_R, and thereafter the crack deflects to other cleavage planes. Indeed, C_R defined by the continuum mechanics ignores atomistic phenomena occurring during rapid crack propagation, and therefore is limited in predicting the crack velocity. Examination of these anomalies shows that this limitation lies in microstructural lattice arrangement and in anisotropic phonon radiation during rapid crack propagation. Globally, C_R has no influence on the crack deflection phenomenon. However, the misfit in C_R between the original plane of propagation and the deflected plane generates local instabilities along the deflection zone.     

Non-linear modulation of SH waves in a two-layered plate and formation of surface SH waves

Non-linear modulation of shear horizontal (SH) waves in a two-layered elastic plate of uniform thickness is considered. Both layers are assumed to be homogeneous, isotropic and incompressible elastic and having different mechanical properties. The problem is investigated by a perturbation method and in the analysis it is assumed that between the linear shear velocities of the top layer, c₁, and the bottom layer, c₂, the inequality c₁<c₂ is valid. In the layered structure then an SH wave exists if the wave velocity c of the wave satisfies either the condition c₁<c?c₂ or the one c₁<c₂?c. Here the problem is examined under the former condition and it is shown that the non-linear modulation of SH waves is governed by a non-linear Schrödinger equation. In this case the formation of surface SH (Love) waves is also revealed if the top layer is thinner when compared with the bottom layer. Then the stability condition is discussed and the existence of bright (envelope) and dark solitons are manifested.     

Effect of rigid boundary on propagation of torsional surface waves in porous elastic layer

The paper presents the effect of a rigid boundary on the propagation of torsional surface waves in a porous elastic layer over a porous elastic half-space using the mechanics of the medium derived by Cowin and Nunziato (Cowin, S. C. and Nunziato, J. W. Linear elastic materials with voids. Journal of Elasticity, 13(2), 125–147 (1983)). The velocity equation is derived, and the results are discussed. It is observed that there may be two torsional surface wave fronts in the medium whereas three wave fronts of torsional surface waves in the absence of the rigid boundary plane given by Dey et al. (Dey, S., Gupta, S., Gupta, A. K., Kar, S. K., and De, P. K. Propagation of torsional surface waves in an elastic layer with void pores over an elastic half-space with void pores. Tamkang Journal of Science and Engineering, 6(4), 241–249 (2003)). The results also reveal that in the porous layer, the Love wave is also available along with the torsional surface waves. It is remarkable that the phase speed of the Love wave in a porous layer with a rigid surface is different from that in a porous layer with a free surface. The torsional waves are observed to be dispersive in nature, and the velocity decreases as the oscillation frequency increases.     

A new framework for studying the stability of genus-1 and genus-2 KP patterns

The Kadomtsev-Petviashvili equation - or KP equation - is a model equation for waves that are weakly two-dimensional in a horizontal plane, and models water waves in shallow water with weak three-dimensionality. It has a vast array of interesting genus—k pattern solutions which can be obtained explicitly in terms of Riemann theta functions. However the linear or nonlinear stability of these patterns has not been studied. In this paper, we present a new formulation of the KP model as a Hamiltonian system on a multi-symplectic structure. While it is well-known that the KP model is Hamiltonian - as an evolution equation in time - multi-symplecticity assigns a distinct symplectic operator for each spatial direction as well, and is independent of the integrability of the equation. The multi-symplectic framework is then used to formulate the linear stability problem for genus-1 and genus-2 patterns of the KP equation; generalizations to genus—k with k > 2 are also discussed.     

Body wave propagation in rotating elastic media

We investigate the propagation of elastic waves through an elastic medium submitted to an angular rotation Ω. Wave propagation is shown to be directly related to the Kibel number Ki=ω/Ω, where ω is the wave frequency. Two dispersive waves W₁ and W₂ are obtained which tend to the classical dilatational and shear waves, respectively, when Ki tends to infinity. Wave W₁ shows a cutoff frequency ω_c=Ω below which it does not propagate. The case of small angular rotation Ω is also studied. The corrections to be introduced to dilatational and shear waves are then shown to be of order O(Ki⁻¹).     

Hopf bifurcation and heteroclinic orbit in a 3D autonomous chaotic system

In this paper we revisit a 3D autonomous chaotic system, which contains both the modified Lorenz system and the conjugate Chen system, presented in [Huang and Yang, Chaos Solitons Fractals 39:567–578, 2009]. First by citing two examples to show the errors and limitations for the local stability of the equilibrium point S ₊ obtained in this literature, we formulate a complete determining criterion for the local stability of S ₊ of this system. Although the local bifurcation problem of this system, mainly for Hopf bifurcation, etc., has been studied, the invoking of incorrect proposition leads to an incorrect result for Hopf bifurcation. We then renew the study of the Hopf bifurcation of this system by utilizing the Project Method. The global bifurcation problem, relatively speaking, should be more difficult than the local bifurcation problem for a given system. However, the global bifurcation problem of this system, to the best of our knowledge, has not been investigated yet in the literatures. So next we consider the global bifurcation problem for this system, mainly for the existence of homoclinic and heteroclinic orbits. Our results, one of which shows the existence of two heteroclinic orbits, not only correct and further supplement the ones obtained in the literature, but also give something new to theoretically help fully understand the occurrence of chaos.     

Global nonlinear stability for a triply diffusive convection in a porous layer

A triply convective-diffusive fluid mixture saturating a porous horizontal layer in the Darcy–Oberbeck–Boussinesq scheme is studied. The nonlinear stability analysis of the conduction solution is performed when the layer is heated from below and salted from above by one salt and below by another salt. Denoting by P _i , (i = 1, 2), the salts Prandtl numbers, it is shown that in the cases {P ₁ = 1; P ₂ = 1; P ₁ = P ₂} do not exist subcritical instabilities and the thermal Rayleigh critical number of global stability in a simple closed form is given. The methodology used and the results obtained appear to be new in the existing literature and useful for the applications.     

Global Solution to the Generalized Riemann Problem in the Presence of a Boundary and Contact Discontinuities

This work is a continuation of our previous work. In the present paper we study the global structure stability of the Riemann solution $u=U(\frac{x}{t})$ containing only contact discontinuities for general n×n quasilinear hyperbolic systems of conservation laws in the presence of a boundary. We prove the existence and uniqueness of a global piecewise C ¹ solution containing only contact discontinuities to a class of the generalized Riemann problems for general n×n quasilinear hyperbolic systems of conservation laws in a half space. Our result indicates that this kind of Riemann solution $u=U(\frac{x}{t})$ mentioned above for general n×n quasilinear hyperbolic systems of conservation laws in the presence of a boundary possesses a global nonlinear structure stability. Some applications to quasilinear hyperbolic systems of conservation laws occurring in physics and other disciplines, particularly to the system describing the motion of the relativistic string in Minkowski space R ^1?+?n , are also given.     

Instability of channel flow of a shear-thinning White–Metzner fluid

We consider the inertialess planar channel flow of a White–Metzner (WM) fluid having a power-law viscosity with exponent n. The case n = 1 corresponds to an upper-convected Maxwell (UCM) fluid. We explore the linear stability of such a flow to perturbations of wavelength k⁻¹. We find numerically that if n < n_c ≈ 0.3 there is an instability to disturbances having wavelength comparable with the channel width. For n close to n_c, this is the only unstable disturbance. For even smaller n, several unstable modes appear, and very short waves become unstable and have the largest growth rate. If n exceeds n_c, all disturbances are linearly stable. We consider asymptotically both the long-wave limit which is stable for all n, and the short-wave limit for which waves grow or decay at a finite rate independent of k for each n.The mechanism of this elastic shear-thinning instability is discussed.     

Numerical solution of the regularized long wave equation using nonpolynomial splines

In this paper, we employ nonpolynomial spline (NPS) basis functions to obtain approximate solutions of the regularized long wave (RLW) equation. By considering suitable relevant parameters, it is shown that the local truncation error behaves O(k ²+h ²) with respect to the time and space discretization. Numerical stability of the method is investigated by using a linearized stability analysis. To illustrate the applicability and efficiency of the aforementioned basis, we compare obtained numerical results with other existing recent methods. Motion of single solitary wave and double and triple solitary waves, wave undulation, generation of solitary waves using the Maxwellian initial condition and conservation properties of mass, energy, and momentum of numerical solutions of the equation are dealt with.     

Stability of Periodic Solutions¶of Conservation Laws with Viscosity:¶Analysis of the Evans Function

Nonclassical conservation laws with viscosity arising in multiphase fluid and solid mechanics exhibit a rich variety of traveling-wave phenomena, including homoclinic (pulse-type) and periodic solutions along with the standard heteroclinic (shock, or front-type) solutions. Here, we investigate stability of periodic traveling waves within the abstract Evans-function framework established by R. A. Gardner. Our main result is to derive a useful stability index analogous to that developed by Gardner and Zumbrun in the traveling-front or -pulse context, giving necessary conditions for stability with respect to initial perturbations that are periodic on the same period T as the traveling wave; moreover, we show that the periodic-stability index has an interpretation analogous to that of the traveling-front or -pulse index in terms of well-posedness of an associated Riemann problem for an inviscid medium, now to be interpreted as allowing a wider class of measure-valued solutionsor, alternatively, in terms of existence and nonsingularity of a local “mass map” from perturbation mass to potential time-asymptotic T-periodic states. A closely related calculation yields also a complementary long-wave stability criterion necessary for stability with respect to periodic perturbations of arbitrarily large period NT, N → ∞. We augment these analytical results with numerical investigations analogous to those carried out by Brin in the traveling-front or -pulse case, approximating the spectrum of the linearized operator about the wave.The stability index and long-wave stability criterion are explicitly evaluable in the same planar, Hamiltonian cases as is the index of Gardner and Zumbrun, and together yield rigorous results of instability similar to those obtained previously for pulse-type solutions; this is established through a novel dichotomy asserting that the two criteria are in certain cases logically exclusive. In particular, we obtain results bearing on the nature and mechanism for formation of highly oscillatory Turing-like patterns observed numerically by Frid and Liu and ?ani? and Peters in models of multiphase flow. Specifically, for the van der Waals model considered by Frid and Liu, we show instability of all periodic waves such that the period increases with amplitude in the one-parameter family of nearby periodic orbits, and in particular of large- and small-amplitude waves; for the standard, double-well potential, this yields instability of all periodic waves.Likewise, for a quadratic-flux model like that considered by ?ani? and Peters, we show instability of large-amplitude waves of the type lying near observed patterns, and of all small-amplitude waves; our numerical results give evidence that intermediate-amplitude waves are unstable as well. These results give support for an alternative mechanism for pattern formation conjectured by Azevedo, Marchesin, Plohr, and Zumbrun, not involving periodic waves.     

Global non-linear stability in double diffusive convection via hidden symmetries

Via a new strategy, the global non-linear stability in double diffusive convection in porous layers (heated from below and salted from above), is investigated. The condition, necessary and sufficient, either for the absence of subcritical instabilities or for the global non-linear stability is obtained by (1) discovering symmetries hidden in the Darcy–Boussinesq equations; (2) introducing “symmetrizable double or triply diffusive convections”.     

线性黏弹性球面波的特征线分析

基于ZWT黏弹性本构方程建立了体现高应变率效应的黏弹性球面波的控制方程组,包含5个偏微分方程,解5个未知量v、σr、σθ、εr和εθ。采用特征线法,问题转化为解3族特征线上的5个常微分方程,物理上图像清晰,数学上易于求解。特征线数值分析显示,黏弹性球面波的衰减和弥散效应超过线弹性球面波。球面扩散引起的环向拉应力是导致介质拉伸破坏的主因。进一步还针对强间断黏弹性球面波得出其衰减特性的解析表达式,表明这种更强的衰减特性是几何扩散效应和本构黏性效应两者共同作用的后果。     

Propagation of Love waves in non-homogeneous substratum over initially stressed heterogeneous half-space

The paper studies the propagation of Love waves in a non-homogeneous substratum over an initially stressed heterogeneous half-space. The dispersion equation of phase velocity is derived. The velocities of Love waves are calculated numerically as a function of kH and presented in a number of graphs, where k is the wave number, and H is the thickness of the layer. The case of Gibson’s half-space is also considered. It is observed that the speed of Love waves is finite in the vicinity of the surface of the half-space and vanishes as the depth increases for a particular wave number. It is also observed that an increase in compressive initial stresses causes decreases of Love waves velocity for the same frequency, and the tensile initial stress of small magnitude in the half-space causes increase of the velocity.     

Transient heat and mass transfer in a natural circulation loop

Comprehensive work has been performed by theoretical and numerical methods in order to study the steady state, transient and stability characteristics of a double diffusive natural circulation loop. It was found that the behavior of the flow in the system depends on the initial conditions and on the location of the state in the seven-parameter space of the thermal and saline Rayleigh numbers,Ra _T ,Ra _S , the modified Prandtl and Schmidt numbers,Pr, Sc, the dimensionless heat and mass transfer coefficients,H _T ,H _S , and the “aspect ratio” (between the height and width) of the loop, γ. Numerical results are presented here, showing the flow in each of the five regions formed in the stability chart. The steady state solutions include convection (constant velocity flow), conduction (no-flow) and periodic with constant amplitude and frequency. Two main new results were obtained: long term periodic oscillations where the amplitude is not symmetric around the conduction solution, and an overshoot of the velocity in transients before reaching the stable convection solutions. In the monotonic instability region of the conduction solution, convection solutions (constant velocity flow) develop, and in the global stability region the flow decays to the conduction solution (no flow), regardless of the initial conditions.     

Phase portraits and bifurcations of the non-linear oscillator: ẍ + (α + γx2 + βx + δx3 = 0

The non-linear oscillator ? + αx + γx²x + βx + δx³ = 0 is studied using the methods of differentiable dynamics to obtain qualitative behaviour. The case x, β<0; γ, δ> 0 is considered in some detail; it has physical relevance as a simple model in certain now-induced structural vibration problems in which the structural non-linearities act to maintain overall stability. The presence of local and global bifurcations is detected and their physical significance discussed.     

The effects of the rock bridge ligament angle and the confinement on crack coalescence in rock bridges: An experimental study and discrete element method

The present article investigates the influences of the rock bridge ligament angle, β, and the confinement on crack coalescence patterns by conducting laboratory and numerical tests on rock-like specimens. Laboratory tests show that no coalescence in the rock bridge occurred for low β. With an increase of β, tensile-shear coalescence and tensile coalescences subsequently occurred. In addition, the increase in the confinement first promoted shear coalescence and then restrained crack coalescence for low β, whereas the tensile coalescence was restrained by the increase in confinement for high β. The numerical results corroborate the laboratory tests in the coalescence patterns. In addition, the numerical study shows that tensile and shear cracks subsequently initiated near crack tips because of the concentrated tensile and shear stresses, respectively. Regarding the influence of β on crack coalescence, tensile or shear stress failed to concentrate in rock bridges for low β. Therefore, the cracks failed to coalesce, whereas with the increase in β, tensile and shear stress concentrations occurred in the bridge and led to either tensile shear or tensile coalescence. Regarding the influence of confinement on crack coalescence, the increase in confinement restrained the tensile stress concentrations and further hindered tensile crack coalescence in rock bridges for high values of β.     

Linear theory of gravity waves on a Voigt viscoelastic medium

Linear surface gravity waves on a semi-infinite incompressible Voigt medium are studied in this paper. Three dimensionless parameters, the dimensionless viscoelastic parameter ϑ, the dimensionless wave number and the dimensionless surface tension are introduced. A dimensionless characteristic equation describing the waves is derived. This is a sixth order complex algebraic equation which is solved to give the complex dispersion relation. Based on the numerical solution, two critical values of ϑ, ϑ _A =0.607 and ϑ _B =2.380, which represent the appearance of the cutoff region and the disappearance of the strong dispersion region, are found. The effects of ϑ on the characteristic equation and the properties of the waves are discussed. The project supported by the National Natural Science Foundation of China (59709006)     

Large-scale edge waves generated by a moving atmospheric pressure

Long waves generated by a moving atmospheric pressure distribution, associated with a storm, in coastal region are investigated numerically. For simplicity the moving atmospheric pressure is assumed to be moving only in the alongshore direction and the beach slope is assumed to be a constant in the on-offshore direction. By solving the linear shallow water equations we obtain numerical solutions for a wide range of physical parameters, including storm size (2a), storm speed (U), and beach slope (α). Based on the numerical results, it is determined that edge wave packets are generated if the storm speed is equal to or greater than the critical velocity, U_cr, which is defined as the phase speed of the fundamental edge wave mode whose wavelength is scaled by the width of the storm size. The length and the location of the positively moving edge wave packet is roughly Ut/2 ≤ y ≤ Ut, where y is in the alongshore direction and t is the time. Once the edge wave packet is generated, the wavelength is the same as that of the fundamental edge wave mode corresponding to the storm speed and is independent of the storm size, which can, however, affect the wave amplitude. When the storm speed is less than the critical velocity, the primary surface signature is a depression directly correlated to the atmospheric pressure distribution.     

On cauchy problems for the RLW equation in two space dimensions

A cnoidal wave solution of the two dimensional RLW equation of are obtained by elliptic integral method, and the some estimations the uniqueness and the stability of the periodic solution with both x, y to the Cauchy problem are proved by the priori estimations. Biography: Huang Zheng-hong (1956-)