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

Analysis of Kic and its temperature dependence of metals by a simplified dislocation model

The relative Kic values of metals are calculated with a simplified dislocation model. It is found that the ratio of KIIc to KIc and the temperature dependence of fracture toughness of some metals estimated with this model are consistent with the experimental results.

Numerical-Value Perturbation Method with Application of Solving Convective-Diffusion Integral Equation

The convective--diffusion equation is of primary importance in such fields as fluid dynamics and heat transfer hi the numerical methods solving the convective-diffusion equation, the finite volume method can use conveniently diversified grids (structured and unstructured grids) and is suitable for very complex geometry The disadvantage of FV methods compared to the finite difference method is that FV-methods of order higher than second are more difficult to develop in three-dimensional cases. The second-order central scheme (2cs) offers a good compromise among accuracy, simplicity and efficiency, however, it will produce oscillatory solutions when the grid Reynolds numbers are large and then very fine grids are required to obtain accurate solution. The simplest first-order upwind (IUW) scheme satisfies the convective boundedness criteria, however. Its numerical diffusion is large. The power-law scheme, QMCK and second-order upwind (2UW) schemes are also often used in some commercial codes. Their numerical accurate are roughly consistent with that of ZCS. Therefore, it is meaningful to offer higher-accurate three point FV scheme. In this paper, the numerical-value perturbational method suggested by Zhi Gao is used to develop an upwind and mixed FV scheme using any higher-order interpolation and second-order integration approximations, which is called perturbational finite volume (PFV) scheme. The PFV scheme uses the least nodes similar to the standard three-point schemes, namely, the number of the nodes needed equals to unity plus the face-number of the control volume. For instanc6, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D problems, 2~Dand 3-D flow model equations. Comparing with other standard three-point schemes, The PFV scheme has much smaller numerical diffusion than the first-order upwind (IUW) scheme, its numerical accuracy are also higher than the second-order central scheme (2CS), the power-law scheme (PLS), the QUICK scheme and the second-order upwind(ZUW) scheme.

Formation of dendritic nanostructures in Pyrex glass anodically bonded to silicon coated with an aluminum thin film

Anodic bonding with thin films of metal or alloy as an intermediate layer, finds increasing applications in micro/nanoelectromechanical systems. At the bonding temperature of 350 degrees C, voltage of 400 V, and 30 min duration, the anodic bonding is completed between Pyrex glass and crystalline silicon coated with an aluminum thin film with a thickness comprised between 50 and 230 nm. Sodium-depleted layers and dendritic nanostructures were observed in Pyrex 7740 glass adjacent to the bonding interface. The sodium depletion width does not increase remarkably with the thickness of aluminum film. The dendritic nanostructures result from aluminum diffusion into the Pyrex glass. This experimental research is expected to enhance the understanding of how the depletion layer and dendritic nanostructures affect the quality of anodic bonding. (C) 2007 Elsevier B.V. All rights reserved.

Twinning partial multiplication at grain boundary in nanocrystalline fcc metals

Most deformation twins in nanocrystalline face-centered cubic fcc metals have been observed to form from grain boundaries. The growth of such twins requires the emission of Shockley partials from the grain boundary on successive slip planes. However, it is statistically improbable for a partial to exist on every slip plane. Here we propose a dislocation reaction and cross-slip mechanism on the grain boundary that would supply a partial on every successive slip plane for twin growth.This mechanism can also produce a twin with macrostrain smaller than that caused by a conventional twin.