Structural Evolution and Property Changes in Nd_{60}Al_{10}Fe_{20}Co_{10} Bulk Metallic Glass During Crystallization

The structural evolution and property changes in Nd60Al10Fe20Co10 bulk metallic glass (BMG) upon crystallization are investigated by the ultrasonic method, x-ray diffraction, density measurement, and differential scanning calorimetry. The elastic constants and Debye temperature of the BMG are obtained as a function of annealing temperature. Anomalous changes in ultrasonic velocities, elastic constants, and density are observed between 600–750 K, corresponding to the formation of metastable phases as an intermediate product in the crystallization process. The changes in acoustic velocities, elastic constants, density, and Debye temperature of the BMG relative to its fully crystallized state are much smaller, compared with those of other known BMGs, the differences being attributed to the microstructural feature of the BMG.

Spallation analysis with a closed trans-scale formulation damage evolution

A closed, trans-scale formulation of damage evolution based on the statistical microdamage mechanics is summarized in this paper. The dynamic function of damage bridges the mesoscopic and macroscopic evolution of damage. The spallation in an aluminium plate is studied with this formulation. It is found that the damage evolution is governed by several dimensionless parameters, i.e., imposed Deborah numbers De* and De, Mach number M and damage number S. In particular, the most critical mode of the macroscopic damage evolution, i.e., the damage localization, is deter-mined by Deborah number De+. Deborah number De* reflects the coupling and competition between the macroscopic loading and the microdamage growth. Therefore, our results reveal the multi-scale nature of spallation. In fact, the damage localization results from the nonlinearity of the microdamage growth. In addition, the dependence of the damage rate on imposed Deborah numbers De* and De, Mach number M and damage number S is discussed.

Numerical Investigation of Richtmyer-Meshkov Instability Driven by Cylindrical Shocks

In this paper, a numerical method with high order accuracy and high resolution was developed to simulate the Richtmyer-Meshkov(RM) instability driven by cylindrical shock waves. Compressible Euler equations in cylindrical coordinate were adopted for the cylindrical geometry and a third order accurate group control scheme was adopted to discretize the equations. Moreover, an adaptive grid technique was developed to refine the grid near the moving interface to improve the resolution of numerical solutions. The results of simulation exhibited the evolution process of RM instability, and the effect of Atwood number was studied. The larger the absolute value of Atwood number, the larger the perturbation amplitude. The nonlinear effect manifests more evidently in cylindrical geometry. The shock reflected from the pole center accelerates the interface for the second time, considerably complicating the interface evolution process, and such phenomena of reshock and secondary shock were studied.

Chaotic dynamic analysis of viscoelastic plates

The dynamic buckling of viscoelastic plates with large deflection is investigated in this paper by using chaotic and fractal theory. The material behavior is given in terms of the Boltzmann superposition principle. in order to obtain accurate computation results, the nonlinear integro-differential dynamic equation is changed into an autonomic four-dimensional dynamical system. The numerical time integrations of equations are performed by using the fourth-order Runge-Kutta method. And the Lyapunov exponent spectrum, the fractal dimension of strange attractors and the time evolution of deflection are obtained. The influence of geometry nonlinearity and viscoelastic parameter on the dynamic buckling of viscoelastic plates is discussed.

Numerical solution of the incompressible Navier-Stokes equations with an upwind compact difference scheme

A new finite difference method for the discretization of the incompressible Navier-Stokes equations is presented. The scheme is constructed on a staggered-mesh grid system. The convection terms are discretized with a fifth-order-accurate upwind compact difference approximation, the viscous terms are discretized with a sixth-order symmetrical compact difference approximation, the continuity equation and the pressure gradient in the momentum equations are discretized with a fourth-order difference approximation on a cell-centered mesh. Time advancement uses a three-stage Runge-Kutta method. The Poisson equation for computing the pressure is solved with preconditioning. Accuracy analysis shows that the new method has high resolving efficiency. Validation of the method by computation of Taylor's vortex array is presented.

Nonequilibrium microstructures and their evolution in a Fe-Cr-W-Ni-C laser clad coating

The laser-solidified microstructural and compositional characterization and phase evolution during tempering at 963 K were investigated using an analytical transmission electron microscope with energy dispersive X-ray analysis. The cladded alloy, a powder mixture of Fe, Cr, W, Ni, and C with a weight ratio of 10:5:1:1:1, was processed with a 3 kW continuous wave CO2 laser. The processing parameters were 16 mm/s beam scanning speed, 3 mm beam diameter. 2 kW laser power, and 0.3 g/s feed rate. The coating was metallurgically bonded to the substrate, with a maximum thickness of 730 mu m, a microhardness of about 860 Hv and a volumetric dilution ratio of about 6%. Microanalyses revealed that the cladded coating possessed the hypoeutectic microstructure comprising the primary dendritic gamma-austenite and interdendritic eutectic consisted of gamma-austenite and M7C3 carbide. The gamma-austenite was a non-equilibrium phase with extended solid solution of alloying elements and a great deal of defect structures, i.e. a high density of dislocations, twins, and stacking faults existed in gamma phase. During high temperature aging, in situ carbide transformation occurred of M7C3 to M23C6 and M6C. The precipitation of M23C6, MC and M2C carbides from austenite was also observed.

Microstructure and Evolution of Mechanically-Induced Ultrafine Grain in Surface Layer of Al-Alloy Subjected to USSP

Experiments were conducted to investigate the ultrafine-grained (UFG) microstructures in the surface layer of an aluminum alloy 7075 heavily worked by ultrasonic shot peening. Conventional and high-resolution electron microscopy was performed at various depths of the deformed layer. Results showed that UFG structures were introdued into the surface layer of 62 μm thick. With increasing strain, the various microstructural features, e.g., the dislocation emission source, elongated microbands, dislocation cells, dislocation cell blocks, equiaxed submicro-, and nano-crystal grains etc., were successively produced. The grain subdivision into the subgrains was found to be the main mechanism responsible for grain refinement. The simultaneous evolution of high boundary misorientations was ascribed to the subgrain boundary rotation for accommodating further strains. Formed microstructures were highly nonequilibratory.  2002 Acta Materialia Inc. Published by Elsevier Science Ltd. All rights reserved.

Characteristics and Sub-Characteristics of Two-Dimensional Parabolized Stability Equations

特征分析表明：对原始扰动量的抛物化稳定性方程组(PSE),它在亚超音速区分别具有椭圆和抛物特性,给出PSE特征对马赫数的依赖关系,阐明PSE仅把信息对流-扩散传播特性抛物化,而保留了信息对流-扰动传播特性,因此PSE应称为扩散抛物化稳定性方程(DPSE)。

Generation of large-scale vortex dislocations in a three-dimensional wake-type flow

Numerical study of three-dimensional evolution of wake-type flow and vortex dislocations is performed by using a compact finite diffenence-Fourier spectral method to solve 3-D incompressible Navier-Stokes equations. A local spanwise nonuniformity in momentum defect is imposed on the incoming wake-type flow. The present numerical results have shown that the flow instability leads to three-dimensional vortex streets, whose frequency, phase as well as the strength vary with the span caused by the local nonuniformity. The vortex dislocations are generated in the nonuniform region and the large-scale chain-like vortex linkage structures in the dislocations are shown. The generation and the characteristics of the vortex dislocations are described in detail.

Perturbational Finite Volume Method For The Solution of 2-D Navier-Stokes Equations on Unstructured and Structured Colocated Meshes

Based on the first-order upwind and second-order central type of finite volume( UFV and CFV) scheme, upwind and central type of perturbation finite volume ( UPFV and CPFV) schemes of the Navier-Stokes equations were developed. In PFV method, the mass fluxes of across the cell faces of the control volume (CV) were expanded into power series of the grid spacing and the coefficients of the power series were determined by means of the conservation equation itself. The UPFV and CPFV scheme respectively uses the same nodes and expressions as those of the normal first-order upwind and second-order central scheme, which is apt to programming. The results of numerical experiments about the flow in a lid-driven cavity and the problem of transport of a scalar quantity in a known velocity field show that compared to the first-order UFV and second-order CFV schemes, upwind PFV scheme is higher accuracy and resolution, especially better robustness. The numerical computation to flow in a lid-driven cavity shows that the under-relaxation factor can be arbitrarily selected ranging from 0.3 to 0. 8 and convergence perform excellent with Reynolds number variation from 102 to 104.

Damage evolution, localization and failure of solids subjected to impact loading

In order to reveal the underlying mesoscopic mechanism governing the experimentally observed failure in solids subjected to impact loading, this paper presents a model of statistical microdamage evolution to macroscopic failure, in particular to spallation. Based on statistical microdamage mechanics and experimental measurement of nucleation and growth of microcracks in an Al alloy subjected to plate impact loading, the evolution law of damage and the dynamical function of damage are obtained. Then, a lower bound to damage localization can be derived. It is found that the damage evolution beyond the threshold of damage localization is extremely fast. So, damage localization can serve as a precursor to failure. This is supported by experimental observations. On the other hand, the prediction of failure becomes more accurate, when the dynamic function of damage is fitted with longer experimental observations. We also looked at the failure in creep with the same idea. Still, damage localization is a nice precursor to failure in creep rupture.

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.

The evolution of shear bands of saturated soil

The development of the shear bands of saturated soil in coupling-rate- and pore-pressure-dependent simple shear has been discussed, using a simple model and a matching technique at the moving boundary of a shear band. Tt is shown that the development of shear bands are dominated by the coupling-rate and pore-pressure effect of the material. The strength of the soil acts as a destabilizer, whilst pore pressure diffusion makes the band expand. The theory is discussed and some computational solutions have been presented.

Morphological Evolution During Synthesis of New AlPO4 Center dot H2O Crystal

A new crystal of aluminophosphate, AIPO(4)(.)H(2)O, is synthesized from two-batch aqueous solution under hydrothermal conditions. Three types of the crystal habits, i.e. the tetragonal double pyramid, the tetragonal prism and the plate-type tetragonal prism, are found from batch-A solution. Two types of the crystal habits, i.e. the hexagonal pyramid and the strip-type tetragonal prism, are found from batch-B solution. The change of crystal morphology is originated from the fluctuation of the synthesis conditions, such as the supersaturation, the temperature and the impurity content. It causes change of the step energies, the defect density and the step roughness, and further, change of the growth rates. Since the crystal morphology is sensitive to the mass transport mechanism, the crystal habits could be changed under the microgravity.

Contingent sensitivity to configuration in a nonlinear evolution system

A trans-scopic sensitivity of macroscopic failure to slight differentiation in the meso-scopic structure of a system with nonlinear evolution is reported. A periodical chain following a non-local load-sharing evolution was applied as a propotype in failure study. The results demonstrate that there is a transition region composed of globally stable (GS) and evolution induced catastrophic (EIC) modes. That is different from a critical threshold as predicted by percolation and renormalization group theories. Moreover, the EIC mode shows a distinctive sample specific behaviour. For instance, some neighbouring initial states may evolve into completely different final states, though different initial states can evolve into the same final states. As an example, a marginal configuration of EIC mode, a quasi-Fibonacci skeleton, is constructed.