Traveling waves in trimer granular lattice I: Bifurcation structure of traveling waves in the unit-cell model

Present paper is the first one in the series devoted to the dynamics of traveling waves emerging in the uncompressed, tri-atomic granular crystals. This work is primarily concerned with the dynamics of one-dimensional periodic granular trimer (tri-atomic) chains in the state of acoustic vacuum. Each unit cell consists of three spherical particles of different masses subject to periodic boundary conditions. Hertzian interaction law governs the mutual interaction of these particles. Under the assumption of zero pre-compression, this interaction is modeled as purely nonlinear, which means the absence of linear force component. The dynamics of such chains is governed by the two system parameters that scale the mass ratios between the particles of the unit cell. Such a system supports two different classes of periodic solutions namely the traveling and standing waves. The primary objective of the present study is the numerical analysis of the bifurcation structure of these solutions with emphasis on the dynamics of traveling waves. In fact, understanding of the bifurcation structure of the traveling wave solutions emerging in the unit-cell granular trimer is rather important and can shed light on the more complex nonlinear wave phenomena emerging in semi-infinite trimer chains. (c) 2016 Elsevier B.V. All rights reserved.

The Structure and Evolution of a Two-Dimensional H2/O2/Ar Cellular Detonation

This paper reports on two-dimensional numerical simulation of cellular detonation wave in a / / mixture with low initial pressure using a detailed chemical reaction model and high order WENO scheme. Before the final equilibrium structure is produced, a fairly regular but still non-equilibrium mode is observed during the early stage of structure formation process. The numerically tracked detonation cells show that the cell size always adapts to the channel height such that the cell ratio is fairly independent of the grid sizes and initial and boundary conditions. During the structural evolution in a detonation cell, even as the simulated detonation wave characteristics suggest the presence of an ordinary detonation, the evolving instantaneous detonation state indicates a mainly underdriven state. As a considerable region of the gas mixture in a cell is observed to be ignited by the incident wave and transverse wave, it is further suggested that these two said waves play an essential role in the detonation propagation.

Extended self-similar scaling law of multi-scale eddy structure in wall turbulence

The longitudinal fluctuating velocity of a turbulent boundary layer was measured in a water channel at a moderate Reynolds number. The extended self-similar scaling law of structure function proposed by Benzi was verified. The longitudinal fluctuating velocity, in the turbulent boundary layer was decomposed into many multi-scale eddy structures by wavelet transform. The extended self-similar scaling law of structure function for each scale eddy velocity was investigated. The conclusions are I) The statistical properties of turbulence could be self-similar not only at high Reynolds number, but also at moderate and low Reynolds number, and they could be characterized by the same set of scaling exponents xi (1)(n) = n/3 and xi (2)(n) = n/3 of the fully developed regime. 2) The range of scales where the extended self-similarity valid is much larger than the inertial range and extends far deep into the dissipation range,vith the same set of scaling exponents. 3) The extended selfsimilarity is applicable not only for homogeneous turbulence, but also for shear turbulence such as turbulent boundary layers.

Shock-Induced Near-Wall Two-Phase Flow Structure Over a Micron-Sized Particles Bed

The present paper studies numerical modelling of near-wall two-phase flows induced by a normal shock wave moving at a constant speed, over a micronsized particles bed. In this two-fluid model, the possibility of particle trajectory intersection is considered and a full Lagrangian formulation of the dispersed phase is introduced. The finiteness of the Reynolds and Mach numbers of the flow around a particle as well as the fineness of the particle sizes are taken into account in describing the interactions between the carrier- and dispersed- phases. For the small mass-loading ratio case, the numerical simulation of flow structure of the two phases is implemented and the profiles of the particle number density are obtained under the constant-flux condition on the wall. The effects of the shock Mach number and the particle size and material density on particle entrainment motion are discussed in detail.The obtained results indicate that interphase non-equilibrium in the velocity and temperature is a common feature for this type of flows and a local particle accumulation zone may form near the envelope of the particle trajectory family.

Step Structure in Cold-Rolled Deformed Nanocrystalline Nickel

The microstructure characteristic of the cold-rolled deformed nanocrystalline nickel metal is studied by transmission electron microscopy. The results show that there are step structures nearby the grain boundary (GB), and the contrast of stress field in front of the step corresponds to the step in the shape. It is indicated that the interaction between twins and dislocations is not a necessary condition to realizing the deformation. In the later stage of the deformation when the grain size becomes about 100nm, the deformation can depend upon the moving of the boundary of the stack faults (SFs) which result from the partial dislocations emitted from GBs. However, when the size of SFs grows up, the local internal stress which is in front of the step gradually becomes higher. When this stress reaches a critical value which stops the gliding of the partial dislocations, the SFs will stop to grow up and leave a step structure behind.

DNS of Compressible Turbulent Boundary Layer Over a Blunt Wedge

Direct numerical simulation of spatially evolving compressible boundary layer over a blunt wedge is performed in this paper. The free-stream Mach number is 6 and the disturbance source produced by wall blowing and suction is located downstream of the sound-speed point. Statistics are studied and compared with the results in incompressible flat-plate boundary layer. The mean pressure gradient effects on the vortex structure are studied.

Finite Boundary Effects in Problem of a Crack Perpendicular to and Terminating at a Bimaterial Interface

In this paper, the problem of a crack perpendicular to and terminating at an interface in bimaterial structure with finite boundaries is investigated. The dislocation simulation method and boundary collocation approach are used to derive and solve the basic equations. Two kinds of loading form are considered when the crack lies in a softer or a stiffer material, one is an ideal loading and the other one fits to the practical experiment loading. Complete solutions of the stress field including the T stress are obtained as well as the stress intensity factors. Influences of T stress on the stress field ahead of the crack tip are studied. Finite boundary effects on the stress intensity factors are emphasized. Comparisons with the problem presented by Chen et al. (Int. J. Solids and Structure, 2003, 40, 2731-2755) are discussed also.

Dns Of Compressible Turbulent Boundary Layer Around A Sharp Cone

Direct numerical simulation of the turbulent boundary layer over a sharp cone with 20 degrees cone angle (or 10 degrees half-cone angle) is performed by using the mixed seventh-order up-wind biased finite difference scheme and sixth-order central difference scheme. The free stream Mach number is 0.7 and free stream unit Reynolds number is 250000/inch. The characteristics of transition and turbulence of the sharp cone boundary layer are compared with those of the flat plate boundary layer. Statistics of fully developed turbulent flow agree well with the experimental and theoretical data for the turbulent flat-plate boundary layer flow. The near wall streak-like structure is shown and the average space between streaks (normalized by the local wall unit) keeps approximately invariable at different streamwise locations. The turbulent energy equation in the cylindrical coordinate is given and turbulent energy budget is studied. The computed results show that the effect of circumferential curvature on turbulence characteristics is not obvious.

Supersonic dusty-gas flows with Knudsen effect in interphase momentum exchange

On the basis of the two-continuum model of dilute gas-solid suspensions, the dynamic behavior of inertial particles in supersonic dusty-gas flows past a blunt body is studied for moderate Reynolds numbers, when the Knudsen effect in the interphase momentum exchange is significant. The limits of the inertial particle deposition regime in the space of governing parameters are found numerically under the assumption of the slip and free-molecule flow regimes around particles. As a model problem, the flow structure is obtained for a supersonic dusty-gas point-source flow colliding with a hypersonic flow of pure gas. The calculations performed using the full Lagrangian approach for the near-symmetry-axis region and the free-molecular flow regime around the particles reveal a multi-layer structure of the dispersed-phase density with a sharp accumulation of the particles in some thin regions between the bow and termination shock waves.

Two-way coupling model for shock-induced laminar boundary-layer flows of a dusty gas

The present paper describes a numerical two-way coupling model for shock-induced laminar boundary-layer flows of a dust-laden gas and studies the transverse migration of fine particles under the action of Saffman lift force. The governing equations are formulated in the dilute two-phase continuum framework with consideration of the finiteness of the particle Reynolds and Knudsen numbers. The full Lagrangian method is explored for calculating the dispersed-phase flow fields (including the number density of particles) in the regions of intersecting particle trajectories. The computation results show a significant reaction of the particles on the two-phase boundary-layer structure when the mass loading ratio of particles takes finite values.

An improvement on the theoretical model for the large-scale coherent structures in the outer region of a turbulent boundary layer

According to the experimental results, there exist large-scale coherent structures in the outer region of a turbulent boundary layer, which have been studied by many authors.As experimental results, Antonia (1990) showed the phase- aver aged streamlines and isovorticity lines of the large-scale coherent structures in a turbulent boundary layer for different Reynolds numbers. Based on the hydrodynamic stability theory, the 2-D theoretical model for the large-scale structures was proposed by Luo and Zhou, in which the eddy viscosity was defined as a complex function of the position in the normal direction. The theoretical results showed in ref. were in agreement with those in ref. However, there were two problems in the results. One is that in the experimental results, there were divergent focuses between two saddle points in the streamlines, but in the theoretical results, there were centers. The other is that the stretched parts of the isovorticity lines appear at the location of centers in the theoretical results, while in the experimental results they located somewhere between the focuses and saddle points. The reason is that the computations were based on a 2-D model.

Mechanical modeling of grain boundary sliding of polycrystals

Using a variational method, a general three-dimensional solution to the problem of a sliding spherical inclusion embedded in an infinite anisotropic medium is presented in this paper. The inclusion itself is also a general anisotropic elastic medium. The interface is treated as a thin interface layer with interphase anisotropic properties. The displacements in the matrix and the inclusion are expressed as polynomial series of the cartesian coordinate components. Using the virtual work principle, a set of linear algebraic equations about unknown coefficients are obtained. Then the general sliding spherical inclusion problem is accurately solved. Based on this solution, a self-consistent method for sliding polycrystals is proposed. Combining this with a two-dimensional model of an aggregate polycrystal, a systematic analysis of the mechanical behaviour of sliding polycrystals is given in detail. Numerical results are given to show the significant effect of grain boundary sliding on the overall mechanical properties of aggregate polycrystals.

Asymmetric structure of an isolated flux arch

This paper extends two-dimensional model of symmetric magnetostatic flux arches confined in stratified atmospheres (Zhang and Hu, 1992, 1993) to asymmetric models. Numerical results show that the flux structure is influenced greatly by the boundary condition of magnetic field, the force-free factor, the atmospheric pressure distribution and the position of footpoints (especially the width ratio of outlet to entrance, which differs from symmetric case).

Equilibrium of isolated magnetostatic flux structure

Two-dimensional magnetostatic models of flux structure confined in stratified atmosphere are discussed in the present paper. The magnetic field in the flux structure is assumed to be force-free at the first step. Numerical solutions for this nonlinear free boundary problem are obtained by finite element method. Results show clearly the relation between the inside fields and outside pressure, especially the influence of atmospheric pressure distribution on the flux structure.

Edge and Size Effects in Micro Composite Structure

Because of the load transfer effect of interface layer, the stress distribution inside the composite structure of film/substrate can be very different from the Timoshenko's model. In this paper, we give the derivation and analysis of such load transfer effect of shear-lag (S-L) model. The micro-structure size (boundary conditions) effect together with interface load transfer effect becomes more and more important as the microstructure size including the three dimensions of thickness, width and length shrinks. The microstructure size is also responsible for the so-called edge-induced stress. The edge effect and difference of S-L model and Timoshenko model are also demonstrated.