Effects of adiabatic exponent on Richtmyer-Meshkov instability

We present a systematical numerical study of the effects of adiabatic exponent gamma on Richtmyer-Meshkov instability (RMI) driven by cylindrical shock waves, based on the gamma model for the multi-component problems and numerical simulation with high-order and high-resolution method for compressible Euler equations. The results show that the RMI of different gamma across the interface exhibits different evolution features with the case of single gamma. Moreover, the large gamma can hold back the development of nonlinear structures, such as spikes and bubbles.

Microstructural evolution and formation of nanocrystalline intermetallic compound during surface mechanical attrition treatment of cobalt

Nanocrystalline intermetallic Co3Fe7 was produced on the surface of cobalt via surface mechanical attrition (SMA). Deformationinduced diffusion entailed the formation of a series of solid solutions. Phase transitions occurred depending on the atomic fraction of Fe in the surface solid solutions: from hexagonal close-packed (<4% Fe) to face-centered cubic (fcc) (4-11% Fe), and from fcc to body-centered cubic (>11% Fe). Nanoscale compositional probing suggested significantly higher Fe contents at grain boundaries and triple junctions than grain interiors. Short-circuit diffusion along grain boundaries and triple junctions dominate in the nanocrystalline intermetallic compound. Stacking faults contribute significantly to diffusion. Diffusion enhancement due to high-rate deformation in SMA was analyzed by regarding dislocations as solute-pumping channels, and the creation of excess vacancies. Non-equilibrium, atomic level alloying can then be ascribed to deformation-induced intermixing of constituent species. The formation mechanism of nanocrystalline intermetallic grains on the SMA surface can be thought of as a consequence of numerous nucleation events and limited growth. (C) 2007 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Dynamic equations for curved submerged floating tunnel

In virtue of reference Cartesian coordinates, geometrical relations of spatial curved structure are presented in orthogonal curvilinear coordinates. Dynamic equations for helical girder are derived by Hamilton principle. These equations indicate that four generalized displacements are coupled with each other. When spatial structure degenerates into planar curvilinear structure, two generalized displacements in two perpendicular planes are coupled with each other. Dynamic equations for arbitrary curvilinear structure may be obtained by the method used in this paper.

Energy principle of ferroelectric ceramics and single domain mechanical model

Many physical experiments have shown that the domain switching in a ferroelectric material is a complicated evolution process of the domain wall with the variation of stress and electric field. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of ferroelectric ceramic and used to study the nonlinear constitutive behavior of ferroelectric body in this paper. The principle of stationary total energy is put forward in which the basic unknown quantities are the displacement u (i) , electric displacement D (i) and volume fraction rho (I) of the domain switching for the variant I. Mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion. On the basis of the domain switching criterion, a set of linear algebraic equations for the volume fraction rho (I) of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. Then a single domain mechanical model is proposed in this paper. The poled ferroelectric specimen is considered as a transversely isotropic single domain. By using the partial experimental results, the hardening relation between the driving force of domain switching and the volume fraction of domain switching can be calibrated. Then the electromechanical response can be calculated on the basis of the calibrated hardening relation. The results involve the electric butterfly shaped curves of axial strain versus axial electric field, the hysteresis loops of electric displacement versus electric filed and the evolution process of the domain switching in the ferroelectric specimens under uniaxial coupled stress and electric field loading. The present theoretic prediction agrees reasonably with the experimental results given by Lynch.

Numerical analysis of spatial evolution of the small signal gain in a chemical oxygen-iodine laser operating without primary buffer gas

A chemical oxygen iodine laser (COIL) that operates without primary buffer gas has become a new way of facilitating the compact integration of laser systems. To clarify the properties of spatial gain distribution, three-dimensional (3-D) computational fluid dynamics (CFD) technology was used to study the mixing and reactive flow in a COIL nozzle with an interleaving jet configuration in the supersonic section. The results show that the molecular iodine fraction in the secondary flow has a notable effect on the spatial distribution of the small signal gain. The rich iodine condition produces some negative gain regions along the jet trajectory, while the lean iodine condition slows down the development of the gain in the streamwise direction. It is also found that the new configuration of an interleaving jet helps form a reasonable gain field under appropriate operation conditions. (c) 2007 Elsevier Ltd. All rights reserved.

Collective evolution characteristics and computer simulation of short fatigue cracks

Fatigue testing was conducted using a kind of triangular isostress specimen to obtain the short-fatigue-crack behaviour of a weld low-carbon steel. The experimental results show that short cracks continuously initiate at slip bands within ferrite grain domains and the crack number per unit area gradually increases with increasing number of fatigue cycles. The dispersed short cracks possess an orientation preference, which is associated with the crystalline orientation of the relevant slip system. Based on the observed collective characteristics, computer modelling was carried out to simulate the evolution process of initiation, propagation and coalescence of short cracks. The simulation provides progressive displays which imitate the appearance of experimental observations. The results of simulation indicate that the crack path possesses a stable value of fractal dimension whereas the critical value of percolation covers a wide datum band, suggesting that the collective evolution process of short cracks is sensitive to the pattern of crack site distribution.

Vortex dynamics in the studies of looping in tropical cyclone tracks

This paper proposes criteria for predicting the tendency of looping in tropical cyclone tracks using the approach of vortex dynamics. We model the asymmetric structure of a cyclone by a system of vortex patches. The evolution of such system of vortices is simulated by the method of contour dynamics. A new set of exact analytic formulas for contour dynamics calculations is derived, which is shown to be more computationally effective. Based on point-vortex models, we derive analytic formulas for the criteria of looping in a cyclone track. From numerical experiments, the simulated trajectories obtained from the point-vortex system and vortex patch system agree quite well. Hence, the looping criteria obtained from the point-vortex system can be applied by forecasters to stay alert for tendency of looping in a cyclone track. To demonstrate the applicability of the proposed criteria, the trajectory of Typhoon Yancy (9012), whose field data are available from ''TCM-90'', is simulated. The case study shows that the asymmetric structure similar to the pattern of a beta gyre is responsible for its recurvature when Yancy landed Fujian Province, China on 20 August 1990.

The mixed-type integral equations for transient elastodynamic plane crack problems

In the present paper, by use of the boundary integral equation method and the techniques of Green fundamental solution and singularity analysis, the dynamic infinite plane crack problem is investigated. For the first time, the problem is reduced to solving a system of mixed-typed integral equations in Laplace transform domain. The equations consist of ordinary boundary integral equations along the outer boundary and Cauchy singular integral equations along the crack line. The equations obtained are strictly proved to be equivalent with the dual integral equations obtained by Sih in the special case of dynamic Griffith crack problem. The mixed-type integral equations can be solved by combining the numerical method of singular integral equation with the ordinary boundary element method. Further use the numerical method for Laplace transform, several typical examples are calculated and their dynamic stress intensity factors are obtained. The results show that the method proposed is successful and can be used to solve more complicated problems.

A stochastic jump and deterministic dynamics model of impact failure evolution with rate effect

Motivated by the observation of the rate effect on material failure, a model of nonlinear and nonlocal evolution is developed, that includes both stochastic and dynamic effects. In phase space a transitional region prevails, which distinguishes the failure behavior from a globally stable one to that of catastrophic. Several probability functions are found to characterize the distinctive features of evolution due to different degrees of nucleation, growth and coalescence rates. The results may provide a better understanding of material failure.

Experimental and theoretical study on numerical density evolution of short fatigue cracks

Fatigue testing was performed using a kind of triangular shaped specimen to obtain the characteristics of numerical density evolution for short cracks at the primary stage of fatigue damage. The material concerned is a structural alloy steel. The experimental results show that the numerical density of short cracks reaches the maximum value when crack length is slightly less than the average grain diameter, indicating grain boundary is the main barrier for short crack extension. Based on the experimental observations and related theory, the expressions for growth velocity and nucleation rate of short cracks have been proposed. With the solution to phase space conservation equation, the theoretical results of numerical density evolution for short cracks were obtained, which were in agreement with our experimental measurements.

Numerical simulation of evolution-induced catastrophe

A numerical simulation of damage evolution in a two-dimensional system of micocracks is presented. It reveals that the failure is induced by a cascade of coalescences of microcracks, and the fracture surface appears fractal. A model of evolution-induced catastrophe is introduced. The fractal dimension is found to be a function of evolution rule only. This result could qualitatively explain the correlation of fractal dimension and fracture toughness discovered in experiments.

Characteristics and microstructure in the evolution of shear localization in ti-6al-4v alloy

A new interrupting method was proposed and the split Hopkinson torsional bar (SHTB) was modified in order to eliminate the effect of loading reverberation on post-mortem observations. This makes the comparative study of macro- and microscopic observations on tested materials and relevant transient measurement of tau - gamma curve possible. The experimental results of the evolution of shear localization in in Ti-6Al-4V alloy studied with the modified SHTB are reported in the paper. The collapse of shear stress seems to be closely related to the appearance of a certain critical coalescence of microcracks. The voids may form within the localized shear zone at a quite early stage. Finally, void coalescence results in elongated cavities and their extension leads to fracture along the shear band.

Statistical description of pattern evolution in damage-fracture

A model of dynamical process and stochastic jump has been put forward to study the pattern evolution in damage-fracture. According to the final states of evolution processes, the evolution modes can be classified as globally stable modes (GS modes) and evolution induced catastrophic modes (ElC modes); the latter are responsible for fracture. A statistical description is introduced to clarify the pattern evolution in this paper. It is indicated that the appearance of fracture in disordered materials should be depicted by probability distribution function.

Evolution induced catastrophe

Fracture due to coalescence of microcracks seems to be catalogued in a new model of evolution induced catastrophe (EIC). The key underlying mechanism of the EIC is its automatically enlarging interaction of microcracks. This leads to an explosively evolving catastrophe. Most importantly, the EIC presents a fractal dimension spectrum which appears to be dependent on the interaction.

Upwind compact scheme with dispersion control for aerodynamic equations

A compact upwind scheme with dispersion control is developed using a dissipation analogy of the dispersion term. The term is important in reducing the unphysical fluctuations in numerical solutions. The scheme depends on three free parameters that may be used to regulate the size of dissipation as well as the size and direction of dispersion. A coefficient to coordinate the dispersion is given. The scheme has high accuracy, the method is simple, and the amount of computation is small. It also has a good capability of capturing shock waves. Numerical experiments are carried out with two-dimensional shock wave reflections and the results are very satisfactory.