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.

Evolution induced catastrophe in a nonlinear dynamical model of material failure

In order to study the failure of disordered materials, the ensemble evolution of a nonlinear chain model was examined by using a stochastic slice sampling method. The following results were obtained. (1) Sample-specific behavior, i.e. evolutions are different from sample to sample in some cases under the same macroscopic conditions, is observed for various load-sharing rules except in the globally mean field theory. The evolution according to the cluster load-sharing rule, which reflects the interaction between broken clusters, cannot be predicted by a simple criterion from the initial damage pattern and even then is most complicated. (2) A binary failure probability, its transitional region, where globally stable (GS) modes and evolution-induced catastrophic (EIC) modes coexist, and the corresponding scaling laws are fundamental to the failure. There is a sensitive zone in the vicinity of the boundary between the GS and EIC regions in phase space, where a slight stochastic increment in damage can trigger a radical transition from GS to EIC. (3) The distribution of strength is obtained from the binary failure probability. This, like sample-specificity, originates from a trans-scale sensitivity linking meso-scopic and macroscopic phenomena. (4) Strong fluctuations in stress distribution different from that of GS modes may be assumed as a precursor of evolution-induced catastrophe (EIC).