Quantum Monte Carlo study of the CO interaction with a dimer model surface for Cr(110)

The chemisorption of CO on a Cr( 110) surface is investigated using the quantum Monte Carlo method in the diffusion Monte Carlo (DMC) variant and a model Cr2CO cluster. The present results are consistent with the earlier ab initio HF study with this model that showed the tilted/ near-parallel orientation as energetically favoured over the perpendicular arrangement. The DMC energy difference between the two orientations is larger (1.9 eV) than that computed in the previous study. The distribution and reorganization of electrons during CO adsorption on the model surface are analysed using the topological electron localization function method that yields electron populations, charge transfer and clear insight on the chemical bonding that occurs with CO adsorption and dissociation on the model surface.

Atomistic simulations of crack nucleation and intergranular fracture in bulk nanocrystalline nickel

We report large scale molecular dynamics simulations of dynamic cyclic uniaxial tensile deformation of pure, fully dense nanocrystalline Ni, to reveal the crack initiation, and consequently intergranular fracture is the result of coalescence of nanovoids by breaking atomic bonds at grain boundaries and triple junctions. The results indicate that the brittle fracture behavior accounts for the transition from plastic deformation governed by dislocation to one that is grain-boundary dominant when the grain size reduces to the nanoscale. The grain-boundary mediated plasticity is also manifested by the new grain formation and growth induced by stress-assisted grain-boundary diffusion observed in this work. This work illustrates that grain-boundary decohesion is one of the fundamental deformation mechanisms in nanocrystalline Ni.

Diffusion dominated process for the crystal growth of a binary alloy

The pure diffusion process has been often used to study the crystal growth of a binary alloy in the microgravity environment. In the present paper, a geometric parameter, the ratio of the maximum deviation distance of curved solidification and melting interfaces from the plane to the radius of the crystal rod, was adopted as a small parameter, and the analytical solution was obtained based on the perturbation theory. The radial segregation of a diffusion dominated process was obtained for cases of arbitrary Peclet number in a region of finite extension with both a curved solidification interface and a curved melting interface. Two types of boundary conditions at the melting interface were analyzed. Some special cases such as infinite extension in the longitudinal direction and special range of Peclet number were reduced from the general solution and discussed in detail.

A new dynamic model for study of dislocation pattern formation

Based on the principle given in nonlinear diffusion-reaction dynamics, a new dynamic model for dislocation patterning is proposed by introducing a relaxation time to the relation between dislocation density and dislocation flux. The so-called chemical potential like quantities, which appear in the model can be derived from variation principle for free energy functional of dislocated media, where the free energy density function is expressed in terms of not only the dislocation density itself but also their spatial gradients. The Linear stability analysis on the governing equations of a simple dislocation density shows that there exists an intrinsic wave number leading to bifurcation of space structure of dislocation density. At the same time, the numerical results also demonstrate the coexistence and transition between different dislocation patterns.

An anisotropic elastic-plastic constitutive model for single and polycrystalline metals .1. Theoretical developments

An anisotropic elastic-plastic constitutive model for single and polycrystalline metals is proposed. The anisotropic hardening of single crystals, at first, is discussed with the viewpoint of yield surface and a new formulation of it is proposed. Then, a model for the anisotropic hardening of polycrystals is suggested by increasing the number of slip systems and incorporating the interaction of all slip systems. The interaction of grains through grain boundaries is shown to be similar to, and incorporated into, the interaction of slip systems in grains. The numerical predictions and their comparisons with experiments will follow in Part II of this paper.

Molecular dynamics simulation of interaction of a dislocation array from a crack tip with grain boundaries

The interaction of a dislocation array emitted from a crack tip under mode II loading with asymmetric tilt grain boundaries (GBs) is analysed by the molecular dynamics method. The GBs can generally be described by planar and linear matching zones and unmatching zones. All GBs are observed to emit dislocations. The GBs migrated easily due to their planar and linear matching structure and asymmetrical type. The diffusion induced by stress concentration is found to promote the GB migration. The transmissions of dislocations are either along the matched plane or along another plane depending on tilt angle theta. Alternate processes of stress concentration and stress relaxation take place ahead of the pileup. The stress concentration can be released either by transmission of dislocations, by atom diffusion along GBs, or by migration of GBs by formation of twinning bands. The simulated results also unequivocally demonstrate two processes, i.e. asymmetrical GBs evolving into symmetrical ones and unmatching zones evolving into matching ones during the loading process.

Numerical simulation of hypersonic viscous flow around space shuttle with method of diffusion analogy

Hypersonic viscous flow around a space shuttle with M(infinity) = 7, Re = 148000 and angle of attack alpha = 5-degrees is simulated numerically with the special Jacobian matrix splitting technique and simplified diffusion analogy method. With the simplified diffusion analogy method the efficiency of computation and resolution of the shock can be improved.

A perturbational h4 exponential finite-difference scheme for the convective diffusion equation

A perturbational h4 compact exponential finite difference scheme with diagonally dominant coefficient matrix and upwind effect is developed for the convective diffusion equation. Perturbations of second order are exerted on the convective coefficients and source term of an h2 exponential finite difference scheme proposed in this paper based on a transformation to eliminate the upwind effect of the convective diffusion equation. Four numerical examples including one- to three-dimensional model equations of fluid flow and a problem of natural convective heat transfer are given to illustrate the excellent behavior of the present exponential schemes, the h4 accuracy of the perturbational scheme is verified using double precision arithmetic.

H4 exponential finite-difference scheme for convective diffusion-equations

Perturbations are applied to the convective coefficients and source term of a convection-diffusion equation so that second-order corrections may be applied to a second-order exponential scheme. The basic Structure of the equations in the resulting fourth-order scheme is identical to that for the second order. Furthermore, the calculations are quite simple as the second-order corrections may be obtained in a single pass using a second-order scheme. For one to three dimensions, the fourth-order exponential scheme is unconditionally stable. As examples, the method is applied to Burgers' and other fluid mechanics problems. Compared with schemes normally used, the accuracies are found to be good and the method is applicable to regions with large gradients.

Diffusion analogy and shock capturing for solving aerodynamic equations

Improving the resolution of the shock is one of the most important subjects in computational aerodynamics. In this paper the behaviour of the solutions near the shock is discussed and the reason of the oscillation production is investigated heuristically. According to the differential approximation of the difference scheme the so-called diffusion analogy equation and the diffusion analogy coefficient are defined. Four methods for improving the resolution of the shock are presented using the concept of diffusion analogy.

Diffusion sliding rate in the binary gaseous mixture

radiation incident upon a test cell filled with gaseous SF6 has