Stability of travelling wave solutions for coupled surface and grain boundary motion

We investigate the spectral stability of the travelling wave solution for the coupled motion of a free surface and grain boundary that arises in materials science. In this problem a grain boundary, which separates two materials that are identical except for their crystalline orientation, evolves according to mean curvature. At a triple junction, this boundary meets the free surfaces of the two crystals, which move according to surface diffusion. The model is known to possess a unique travelling wave solution. We study the linearization about the wave, which necessarily includes a free boundary at the location of the triple junction. This makes the analysis more complex than that of standard travelling waves, and we discuss how existing theory applies in this context. Furthermore, we compute numerically the associated point spectrum by restricting the problem to a finite computational domain with appropriate physical boundary conditions. Numerical results strongly suggest that the two-dimensional wave is stable with respect to both two- and three-dimensional perturbations.     

An atomic scale characterization of coupled grain boundary motion in silicon bicrystals

The mechanical response of symmetric tilt grain boundaries (GBs) in silicon bicrystals under shear loading are characterized using molecular dynamics simulations. It is seen that under shear, high-angle GBs namely Σ5 and Σ13 having a rotation axis [0 0 1] demonstrate coupled GB motion, such that the displacement of grains parallel to the GB interface is accompanied by normal GB motion. An atomic-scale characterization revealed that concerted rotations of silicon tetrahedra within the GB are the primary mechanisms leading to the coupled GB motion. Interestingly, so far, this phenomenon has only been examined in detail for metallic systems. A distinguishing feature of the coupled GB motion observed for the silicon symmetric tilt bicrystals as compared to metallic bicrystals is the fact that in the absence of shear, spontaneous coupled motion is not observed at high temperatures.     

Grain boundary interdiffusion in the case of concentration-dependent grain boundary diffusion coefficient

The grain boundary diffusion in a binary system which exhibits a grain boundary phase transition is considered in the framework of Fisher's model. The kinetic law of the growth of the grain boundary phase and the distribution of the diffusant near the grain boundary are calculated. The method of determining of the concentration dependence of the grain boundary diffusion coefficient from the experimentally measured penetration profiles of the diffusant along the grain boundaries is suggested. The experimental results on Zn diffusion in Fe(Si) bicrystals, Ni diffusion in Cu bicrystals and grain boundary grooving in Al in the presence of liquid In are discussed in light of the suggested model.     

Theory and simulation of coupled grain boundary migration and grain rotation in low angle grain boundaries

Abstract

Microstructure evolution due to coupled grain boundary migration and grain rotation in low angle grain boundaries is studied through a combination of molecular dynamics and phase field modeling. We have performed two dimensional molecular dynamics simulations on a bicrystal with a circular grain embedded in a larger grain. Both size and orientation of the embedded grain are observed to evolve with time. The shrinking embedded grain is observed to have two regimes: constant dislocation density on the grain boundary followed by constant rate of increase in dislocation density. Based on these observations from the molecular dynamics simulations, a theoretical formulation of the kinetics of coupled grain rotation is developed. The grain rotation rate is derived for the two regimes of constant dislocation density and constant rate of change of dislocation density on the grain boundary during evolution. The theoretical calculation of the grain rotation rate shows strong dependence on the grain size and compares very well with the molecular dynamics simulations. A multi-order parameter based phase field model with coupled grain rotation is developed using the theoretical formulation to model polycrystalline microstructure evolution.     

The stability of triple junctions during the grain boundary diffusion process

The grain boundary diffusion in a system with triple junctions is considered in such a geometry, in which the flows of diffusing atoms meet at the triple line. The solutions of the diffusion equation is given in the frameworks of Fisher's model and under the assumption of quasi-stationary distribution of the diffusing atoms along the grain boundaries. The change of the mechanical equilibrium at the triple junction due to the increase of the concentration of solute atoms is considered. It is shown that under some circumstances the triple junction looses its stability with respect to migration in the direction to the diffusion source. The stability diagrams in the segregation-diffusivity parameter space are plotted.     

A study of the determination of grain boundary diffusivity and energy through the thermally grown oxide ridges on a Fe-22Cr alloy surface

Abstract

The grain boundaries (GBs) present in polycrystalline materials are important with respect to materials behaviour and properties. During the transient stage of oxidation, the higher GB diffusivity results in heterogeneous oxidation structures in the form of oxide ridges that emerge along the alloy GBs. In an attempt to delve into the more fundamental aspects of the GBs, such as GB energy, the size of the oxide ridges was quantitatively measured by atomic force microscopy on the post oxidation surface of a Fe-22 wt % Cr alloy after an oxidation exposure at 800 °C in dry air. The GB diffusivity was calculated utilising the ridge size data and the relationship between the GB diffusivity and the GB characteristics was determined. Furthermore, the GB energy was calculated from the GB diffusivity data, also to make comparison with the data available in the literature. The absolute value of the calculated GB energy was quite close to the values reported in the literature. However, compared to the extremely low temperature (0 K) data-set from the literature, the data-set obtained from this study showed much less spread. The smaller variation range may be attributed to the higher temperature condition (1073 K) in this study.     

Anisotropic rheology during grain boundary diffusion creep and its relation to grain rotation,grain boundary sliding and superplasticity

The response of periodic microstructures to deformation can be analysed rigorously and this provides guidance in understanding more complex microstructures. When deforming by diffusion creep accompanied by sliding, irregular hexagons are shown to be anisotropic in their rheology. Analytic solutions are derived in which grain rotation is a key aspect of the deformation. If grain boundaries cannot support shear stress, the polycrystal viscosity is extremely anisotropic. There are two orthogonal directions of zero strength: sliding and rotation cooperate to allow strain parallel to these directions to be accomplished without any dissolution or plating. When a linear velocity/shear stress relationship is introduced for grain boundaries, the anisotropy is less extreme, but two weak directions still exist along which polycrystal strength is controlled only by the grain boundary “viscosity”. Irregular hexagons are characterised by four parameters. A particular subset of hexagons defined by two parameters, which includes regular hexagons as well as some elongate shapes, shows singular behaviour. Grain shapes that are close to that of the subset may exhibit large grain rotation rates and have no well-defined rheology unless there is a finite grain boundary viscosity. This new analysis explains why microstructures based on irregular but near equiaxed grains show high rotation rates during diffusion creep and it provides a framework for understanding strength anisotropy during diffusion creep.     

Grain boundary misorientation changes during grain growth in pure aluminium

A new method is described for data-logging large amounts of grain boundary misorientation information from channelling patterns in the scanning electron microscope (SEM). The method relies on producing specimens where the grain size is larger than the specimen thickness and where the grain boundary planes are perpendicular to the specimen plane (the so-called columnar structure). Results for grain growth in pure aluminium at 460 and 500°C are presented. There is an increase in the proportion of low angle boundaries at the expense of high angle boundaries during growth times of up to a few hours. The reasons are thought to be partly connected with lower low angle boundary mobility compared with high angle boundaries. However, the growth kinetics appear to be normal over the entire growth time range.     

Significant decrease in interfacial energy of grain boundary through serrated grain boundary transition

The mechanism of serrated grain boundary formation and its effect on liquation behaviour have been studied in a wrought nickel-based superalloy – Alloy 263. It was newly discovered that grain boundaries are considerably serrated in the absence of γ?′-phase or M₂₃C₆ at the grain boundaries. An electron energy-loss spectroscopy study suggests that serration is triggered by the discontinuous segregation of C and Cr atoms at grain boundaries for the purpose of relieving the excessive elastic strain energy. The grain boundaries serrate to have specific segments approaching one {111} low-index plane at a boundary so that the interfacial free energy of the grain boundary can be decreased, which may be responsible for the driving force of the serration. The serrated grain boundaries effectively suppress grain coarsening and are highly resistant to liquation due to their lower wettability resulting from a lower interfacial energy of the grain boundary.     

Analysis of kinetics regimes in grain boundary self-diffusion†

The evolution of self-diffusion concentration depth profiles into a bulk polycrystalline material from a thin-film tracer source was studied using the lattice Monte Carlo method. The onset of the Type-B kinetics regime from the Type-A kinetics regime was investigated for a cubic grain model, which can be expected to capture the three-dimensional nature of the problem better than the usual grain boundary slab model. The solutions of Whipple–Le Claire (derived for the constant source condition), Suzuoka (derived for the instantaneous source condition), Bokstein and colleagues (derived for the constant source condition) and Levine and MacCallum (derived for the constant source condition) are tested. All these solutions describe the tail region of the profiles well for this model (provided that a simple empirical factor in the case of the first three solutions is introduced). The intermediate Type-AB kinetics regime proposed by Divinski and Larikov between the Type-A and Type-B kinetics regimes was also investigated but could not be located from analyses of best-fit exponents. Forcing the fitting nonetheless to the tail regions shows that their analysis needs some further refinement before application to this model.     

Dependence of grain boundary character distribution on the initial grain size of 304 austenitic stainless steel

Abstract

In order to study the dependence of the grain boundary character distributions (GBCD) on the grain size, annealing treatment was carried out on 304 austenitic stainless steel with different initial grain sizes. The evolution of the GBCD was analysed by electron backscatter diffraction. The experimental results showed that abnormal grain growth (AGG) occurred when grain size was small. With a smaller initial grain size, the number density of abnormally large grains and the fraction of low-Σ CSL boundaries increased but the size of abnormally large grains decreased and the random boundaries presented a continuous network. With a larger initial grain size, the fraction of low-Σ CSL boundaries also increased as well as the size of abnormally large grains but the number density of abnormally large grains decreased and the connectivity of random boundary network was disrupted by low-Σ CSL boundaries, especially Σ3ⁿ (n = 1, 2, 3) boundaries. However, with a very large initial grain size, normal grain growth (NGG) occurred, which had no effect on the fraction of low-Σ CSL boundaries and the connectivity of random boundary network.     

Coercivity enhancement of sintered Nd-Fe-B magnets by grain boundary diffusion with Pr80-xAlxCu20 alloys

A grain boundary diffusion (GBD) process with Pr_80-xAl_xCu₂₀ (x = 0, 10, 15, 20) low melting point alloys was applied to commercial 42M sintered Nd-Fe-B magnets. The best coercivity enhancement of a diffused magnet was for the Pr₆₅Al₁₅Cu₂₀ GBD magnet, from 16.38 kOe to 22.38 kOe. Microstructural investigations indicated that increase in the Al content in the diffusion source can form a continuous grain boundary (GB) phase, optimizing the microstructure to enhance the coercivity. The coercivity enhancement is mainly due to the formation of a continuous GB phase to separate the main phase grains. Exchange decoupling between the adjacent main phase grains is enhanced after the GBD process. Meanwhile, the introduction of Al can effectively promote the infiltration of Pr into the magnet, which increases the diffusion rate of rare-earth elements within a certain range. This work provides a feasible method to enhance coercivity and reduce the use of rare-earth resources by partial replacement of rare-earth elements with non-rare-earth elements in the diffusion source.     

The Harrison diffusion kinetics regimes in solute grain boundary diffusion

Knowledge of the limits of the principal Harrison kinetics regimes (Types A, B and C) for grain boundary diffusion is very important for the correct analysis of depth profiles in a tracer diffusion experiment. These regimes for self‐diffusion have been extensively studied in the past by making use of the phenomenological lattice Monte Carlo (LMC) method with the result that the limits are now well established. However, the relationship of these self‐diffusion limits to the corresponding ones for solute diffusion in the presence of solute segregation to the grain boundaries remains unclear. In the present study, the influence of solute segregation on the limits was investigated with the LMC method for the well‐known parallel grain boundary slab model by showing the equivalence of two diffusion models. It is shown which diffusion parameters are useful for identifying the limits of the Harrison kinetics regimes for solute grain boundary diffusion. It is also shown how the measured segregation factor from the diffusion experiment in the Harrison Type‐B kinetics regime may differ from the global segregation factor.     

Grain Boundary Motion during Ag and Cu Grain Boundary Diffusion in Cu Polycrystals

Grain boundary (GB) motion in high-purity Cu material (5N8 and 5N Cu) is investigated using the results of radiotracer GB diffusion measurements with tracers exhibiting fundamental differences in the solute-matrix atom interactions. The results on GB solute diffusion of Ag (revealing a miscibility gap in the Ag-Cu phase diagram) and Au (forming intermetallic compounds with Cu) in Cu and on Cu self-diffusion are analyzed.The initial parts of the Ag and Cu penetration profiles turned out to be substantially curved. The profile curvature is explained via the effect of GB motion during ^110m Ag and ⁶⁴Cu GB penetration. The activation enthalpies of GB motion in these two independent measurements occurred to be very close, 95 and 103 kJ/mol, respectively. Moreover, these values turn out to be close, but still somewhat larger than the activation enthalpy of Cu GB self-diffusion in Cu material of the same very high purity, Q ^Cu _gb = 72 kJ/mol. Although tracer diffusion measurements of Au GB diffusion in Cu yielded only limited information on GB motion, the absolute values of GB velocities are consistent with those calculated from the Ag and Cu GB diffusion data.     

The study of grain orientation distributions and grain boundary misorientation distributions in 304 and 316L stainless steels

Orientation distribution functions in two recrystallized austenitic stainless steels (AISI types 304 and 316L) with known grain boundary misorientation distributions have been studied. Previously obtained data on grain boundary spectra in these steels have been re-examined and analyzed from the point of view of texture analysis.The results obtained have shown that there is no unambiguous relatonship between grain boundary misorientation distribution and grain orientation distribution (ODF) determined by the X-ray analysis in the materials under study. This ambiguity is due to the following reason. In the grain boundary misorientation statistics only nearest-neighbor grains are taken into account, but in the orientation distribution function orientations are averaged over the entire volume of the specimen independent as to whether the grains are adjacent or not. Two main results were established for the steels under study: (i) Textures of the two steels differ, though their grain boundary misorientation distributions are similar; and (ii) misorientations of the majority of grain boundaries can be described as rotations about the axes close to 110.     

Simultaneous enhancement in coercivity and remanence of Nd2Fe14B permanent magnet by grain boundary diffusion process using NdHx

We have successfully developed a Dy-free grain boundary diffusion process with neodymium hydride (NdH_x) alloy to the permanent magnet Nd₂Fe₁₄B powders using hydrogenation – disproportionation – desorption – recombination (HDDR) method. All the diffusion treatments were performed at 700–800 °C for various annealing time under the high vacuum with rotating diffusion method that effectively control the abnormal grain growth. The coercivities of Dy-treated Nd₂Fe₁₄B powders were varied from 9.5 kOe to 13.2 kOe but the remanence was decreased to 8.1 kG (10% reduction) depending on dysprosium hydride (DyH_x) content and diffusion treated time. However, the coercivity and remanence of Dy-free diffusion treated powder have been increased to 12.2 kOe (28.5% enhancement) and 11.1 kG (22% enhancement) at the optimal diffusion treatment (800 °C for 3 h), respectively. This unique simultaneous enhancement is to isolate the magnetic coupling between Nd₂Fe₁₄B grains by creating non-magnetic Nd grain boundaries and enhance the alignment of the Nd₂Fe₁4B hard magnetic phase, fabricated by optimal diffusion conditions.     

Modelling of grain boundary impurity segregation during high temperature plastic deformation

A modified theoretical model is proposed to predict the grain boundary segregation of impurity atoms during high temperature plastic deformation. The model is based on the supersaturated vacancy-impurity complex created by plastic deformation and involves quasi-thermodynamics and kinetics. Model predictions are made for phosphorus grain boundary segregation during plastic deformation in ferrite steel. The results reveal that phosphorus segregates at grain boundaries during plastic deformation. At a given te...     

Co-energy of surface and grain boundary in Ag film

A model is proposed for calculating the co-energy of surface and grain boundary (GB) by the modified analytical embedded atom method (MAEAM). As example, the energy densities E_d of the two adjacent grains are obtained when their (0 0 1) twist GB meets the free surface (h₁ k₁ 0)/(h₂ k₂ 0) of Ag film. The period along the boundary line on the surface is found and the energy density is calculated for the situations either with or without periodicity. The results show that, the energy value achieved via this model can be stable even for most grains with boundary line smaller than 100 nm. Among the grains with (h k 0) surface, (1 0 0) grains should be favored and grow fastest when they meet (1 1 0) grains.     

高温塑性变形引起的P非平衡晶界偏聚

阐述了高温塑性变形引起的非平衡晶界偏聚的准热力学和动力学,并使用该模型预测了低合金结构钢中高温塑性变形导致的P在奥氏体晶界的非平衡偏聚.研究发现：当变形为20%,应变速率为1×10^-3 s^-1时,在800 ℃左右会出现一个P的晶界偏聚浓度峰值;在1000 ℃变形为20％时,晶界偏聚浓度随着应变速率的增加而增加.预测结果与现有的实验结果基本一致. 关键词： 非平衡偏聚 晶界 塑性变形 磷