Molecular dynamics simulation on dislocation emission process

A correlative reference model for a computer simulation of molecular dynamics is proposed in this paper. Based on this model, a flexible displacement boundary scheme is naturally introduced and the dislocations emitted from a crack tip are presumed to continuously pass through the border of an inner discrete atomic region to pile up at an outer continuum region. The simulations for a Mo crystal show that the interaction between a crack and emitted dislocations results in the decrease in local stress intensity factor gradually.

扩散处理热浸镀铝钢高温抗氧化行为的研究

Analysis of two-dimensional finite solids with microcracks

A general method is presented for solving the plane elasticity problem of finite plates with multiple microcracks. The method directly accounts for the interactions between different microcracks and the effect of outer boundary of a finite plate. Analysis is based on a superposition scheme and series expansions of the complex potentials. By using the traction-free conditions on each crack surface and resultant forces relations along outer boundary, a set of governing equations is formulated. The governing equations are solved numerically on the basis of a boundary collocation procedure. The effective Young's moduli for randomly oriented cracks and parallel cracks are evaluated for rectangular plates with microcracks. The numerical results are compared with those from various micromechanics models and experimental data. These results show that the present method provides a direct and efficient approach to deal with finite solids containing multiple microcracks.

Effect of diffusion on coating microstructure and oxidation resistance of aluminizing steel

The effect of diffuse treatment on coating microstructure and oxidation resistance at high-temperature of hot-dip aluminum were studied by means of TEM, SEM and XRD. The results show that, the diffusion temperature has significant effect on structure of coatings and its oxidation resistance. After diffusion at 750 degreesC, the coating consists of thick outer surface layer (Fe2Al5+ FeAl2), thin internal layer (FeAl + stripe FeAl2), and its oxidation resistance is poor. After diffusion at 950 degreesC, the outer surface layer is composed of single FeAl2 phase, the internal layer is composed of FeAl phase, and its oxidation resistance declines due to the occurrence of early stage internal oxidation cracks in the coating. After diffusion at 850 degreesC, the outer surface layer becomes thinner and consists of FeAl2 Fe2Al5(small amount), the internal layer becomes thicker and consists of FeAl+spherical FeAl2, and the spheroidized FeAl2 phase in the internal layer and its existing in FeAl phase steadily improve the oxidation resistance of the coating.

Some advances in study of crack identification problems

In the present paper, the crack identification problems are investigated. This kind of problems belong to the scope of inverse problems and are usually ill-posed on their solutions. The paper includes two parts: (1) Based on the dynamic BIEM and the optimization method and using the measured dynamic information on outer boundary, the identification of crack in a finite domain is investigated and a method for choosing the high sensitive frequency region is proposed successfully to improve the precision. (2) Based on 3-D static BIEM and hypersingular integral equation theory, the penny crack identification in a finite body is reduced to an optimization problem. The investigation gives us some initial understanding on the 3-D inverse problems.

Numerical Simulation and Experimental Investigation about Internal and External Flows

In this paper, TASCflow3D is used to solve inner and outer 3D viscous incompressible turbulent flow (R-e = 5.6 X 10(6)) around axisymmetric body with duct. The governing equation is a RANS equation with standard k-epsilon turbulence model. The discrete method used is a finite volume method based on the finite element approach. In this method, the description of geometry is very flexible and at the same time important conservative properties are retained. The multi-block and algebraic multi-grid techniques are used for the convergence acceleration. Agreement between experimental results and calculation is good. It indicates that this novel approach can be used to simulate complex flow such as the interaction between rotor and stator or propulsion systems containing tip clearance and cavitation.

Growth Mechanism And Corrosion Behavior Of Ceramic Coatings On Aluminum Produced By Autocontrol Ac Pulse Peo

Ceramic coatings are produced on aluminum alloy by autocontrol AC pulse Plasma Electrolytic Oxidation (PEO) with stabilized average current. Transient signal gathering system is used to study the current, voltage, and the transient wave during the PEO process. SEM, OM, XRD and EDS are used to study the coatings evolution of morphologies, composition and structure. TEM is used to study the micro profile of the outer looser layer and inner compact layer. Polarization test is used to study the corrosion property of PEO coatings in NaCl solution. According to the test results, AC pulse PEO process can be divided into four stages with different aspects of discharge phenomena, voltage and current. The growth mechanism of AC PEO coating is characterized as anodic reaction and discharge sintering effect. PEO coating can increase the corrosion resistance of aluminum alloy by one order or two; however, too long process time is not necessarily needed to increase the corrosion resistance. In condition of this paper, PEO coating at 60 min is the most protective coating for aluminum alloy substrate. (C) 2008 Elsevier B.V. All fights reserved.

Vertically forced surface wave in weakly viscous fluids bounded in a circular cylindrical vessel

Two-time scale perturbation expansions were developed in weakly viscous fluids to investigate surface wave motions by linearizing the Navier-Stokes equation in a circular cylindrical vessel which is subject to a vertical oscillation. The fluid field was divided into an outer potential flow region and an inner boundary layer region. A linear amplitude equation of slowly varying complex amplitude, which incorporates a damping term and external excitation, was derived for the weakly viscid fluids. The condition for the appearance of stable surface waves was obtained and the critical curve was determined. In addition, an analytical expression for the damping coefficient was determined and the relationship between damping and other related parameters (such as viscosity, forced amplitude, forced frequency and the depth of fluid, etc.) was presented. Finally, the influence both of the surface tension and the weak viscosity on the mode formation was described by comparing theoretical and experimental results. The results show that when the forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when the forcing frequency is high, the surface tension of the fluid is prominent.

Numerical simulation of Richtmyer-Meshkov instability

The compressible Navier-Stokes equations discretized with a fourth order accurate compact finite difference scheme with group velocity control are used to simulate the Richtmyer-Meshkov (R-M) instability problem produced by cylindrical shock-cylindrical material interface with shock Mach number Ms = 1.2 and density ratio 1:20 (interior density/outer density). Effect of shock refraction, reflection, interaction of the reflected shock with the material interface, and effect of initial perturbation modes on R-M instability are investigated numerically. It is noted that the shock refraction is a main physical mechanism of the initial phase changing of the material surface. The multiple interactions of the reflected shock from the origin with the interface and the R-M instability near the material interface are the reason for formation of the spike-bubble structures. Different viscosities lead to different spike-bubble structure characteristics. The vortex pairing phenomenon is found in the initial double mode simulation. The mode interaction is the main factor of small structures production near the interface.

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.

Numerical study of bifurcation solutions of spherical Taylor-Couette flow

The steady bifurcation flows in a spherical gap (gap ratio sigma=0.18) with rotating inner and stationary outer spheres are simulated numerically for Re(c1)less than or equal to Re less than or equal to 1 500 by solving steady axisymmetric incompressible Navier-Stokes equations using a finite difference method. The simulation shows that there exist two steady stable flows with 1 or 2 vortices per hemisphere for 775 less than or equal to Re less than or equal to 1 220 and three steady stable flows with 0, 1, or 2 vortices for 1 220

Stratification by solid particle suspensions

Turbulence was generated by an oscillating grid above a bed of sediment of spherical glass beads. As expected, part of the sediment was lifted up by the grid action and a suspension layer of depth D formed above the grid. This depth was found remaining independent of grid action but varying with the sediment layer depth when the grid action was kept constant. Volume concentration measurements show the existence of only weak concentration gradients over the layer depth with a rapid fall off in concentration at the outer edge. The theoretical analysis based on a concentration flux model is in good qualitative agreement with observations.

Hierarchal structure of the simplified Navier-Stokes equations (SNSE) and its mechanical connotation and application

The hierarchial structure and mathematical property of the simplified Navier-Stokesequations (SNSE) are studied for viscous flow over a sphere and a jet of compressible flu-id. All kinds of the hierarchial SNSE can be divided into three types according to theirmathematical property and also into five groups according to their physical content. Amultilayers structure model for viscous shear flow with a main stream direction is pre-sented. For the example of viscous incompressible flow over a flat plate there existthree layers for both the separated flow and the attached flow; the character of thetransition from the three layers of attached flow to those of separated flow is elucidated.A concept of transition layer being situated between the viscous layer and inviscidlayer is introduced. The transition layer features the interaction between viscous flow andinviscid flow. The inner-outer-layers-matched SNSE proposed by the present author inthe past is developed into the layers matched (LsM)-SNSE.

Simplified Navier-Stokes equations (SNSE)

Ten kinds of the simplified Navier-Stokes equations (SNSE) are reviewed and also used to calculate the Jeffery-Hamel flow as well as to analyze briefly the seven kinds of flows to which the exact solutions of the complete Navier-Stokes equations (CNSE) have been found. Analysis shows that the actual differences among the solutions of the different SNSE can go beyond the range of the order of magnitude of Re-1/2 and result even in different flow patterns, therefore, how to choose the viscous terms included in the SNSE is worthy of notice where Re=S∞u∞ L/μ∞ is the Reynolds numbers. For the aforesaid eight kinds of flows, the solutions to the inner-outer-layer-matched SNSE and to the thin-layer-2-order SNSE agree completely with the exact solutions to CNSE. But the solutions to all the other SNSE are not completely consistent with the exact solutions to CNSE and not a few of them are actually the solutions of the classical boundary layer theory. The innerouter-layer-matched SNSE contains the shear stress causing angular displacement of the inormal axis with respect to the streamwise axis and the normal stress causing expansion-contraction in the direction of the normal axis and the viscous terms being of the order of magnitude of the normal stress; and it can also reasonably treat the inertial terms as well as the relation between the viscous and inertial terms. Therefore, it seems promising in respects of both mechanics and mathematics.

Numerical Simulation and Experimental Investigation about Internal and External Flows

In this paper, TASCflow3D is used to solve inner and outer 3D viscous incompressible turbulent flow (R-e = 5.6 X 10(6)) around axisymmetric body with duct. The governing equation is a RANS equation with standard k-epsilon turbulence model. The discrete method used is a finite volume method based on the finite element approach. In this method, the description of geometry is very flexible and at the same time important conservative properties are retained. The multi-block and algebraic multi-grid techniques are used for the convergence acceleration. Agreement between experimental results and calculation is good. It indicates that this novel approach can be used to simulate complex flow such as the interaction between rotor and stator or propulsion systems containing tip clearance and cavitation.