Computation of the Kirkendall velocity and displacement fields in a one-dimensional binary diffusion couple with a moving interface

Author(s):  
W.J Boettinger ◽  
J.E Guyer ◽  
C.E Campbell ◽  
G.B McFadden

The moving interface problem in a one-dimensional binary α/β diffusion couple is studied using sharp and diffuse interface (Cahn–Hilliard) approaches. With both methods, we calculate the solute field and the Kirkendall marker velocity and displacement fields. In the sharp interface treatment, the velocity field is generally discontinuous at the interphase boundary, but can be integrated to obtain a displacement field that is continuous everywhere. The diffuse interface approach avoids this discontinuity, simplifies the integration and yet gives the same qualitative behaviour. Special features observed experimentally and reported in the literature are also studied with the two methods: (i) multiple Kirkendall planes, where markers placed on the initial compositional discontinuity of the diffusion couple bifurcate into two locations, and (ii) a Kirkendall plane that coincides with the interphase interface. These situations occur with special values of the interdiffusion coefficients and starting couple compositions. The details of the deformation in these special situations are given using both methods and are discussed in terms of the stress-free strain rate associated with the Kirkendall effect.

2015 ◽  
Vol 2015 ◽  
pp. 1-11 ◽  
Author(s):  
Marek Danielewski ◽  
Henryk Leszczyński

The problem of Kirkendall’s trajectories in finite, three- and one-dimensional ternary diffusion couples is studied. By means of the parabolic transformation method, we calculate the solute field, the Kirkendall marker velocity, and displacement fields. The velocity field is generally continuous and can be integrated to obtain a displacement field that is continuous everywhere. Special features observed experimentally and reported in the literature are also studied: (i) multiple Kirkendall’s planes where markers placed on an initial compositional discontinuity of the diffusion couple evolve into two locations as a result of the initial distribution, (ii) multiple Kirkendall’s planes where markers placed on an initial compositional discontinuity of the diffusion couple move into two locations due to composition dependent mobilities, and (iii) a Kirkendall plane that coincides with the interphase interface. The details of the deformation (material trajectories) in these special situations are given using both methods and are discussed in terms of the stress-free strain rate associated with the Kirkendall effect. Our nonlinear transform generalizes the diagonalization method by Krishtal, Mokrov, Akimov, and Zakharov, whose transform of diffusivities was linear.


Author(s):  
Erasmo Carrera ◽  
Enrico Zappino

One-dimensional models are widely used in mechanical design. Classical models, Euler-Bernoulli or Timoshenko, ensure a low computational cost but are limited by their assumptions, many refined models were proposed to overcome these limitations and extend one-dimensional models at the analysis of complex geometries or advanced materials. In this work a new approach is proposed to couple different kinematic models. A new finite element is introduced in order to connect one-dimensional elements with different displacement fields. The model is derived in the frameworks of the Carrera Unified Formulation (CUF), therefore the formulation can be written in terms of fundamental nuclei. The results show that the use variable kinematic models allows the computational costs to be reduced without reduce the accuracy, moreover, refined-one dimensional models can be used in the analysis of complex structures.


2013 ◽  
Vol 5 (04) ◽  
pp. 548-568 ◽  
Author(s):  
Tao Lin ◽  
Yanping Lin ◽  
Xu Zhang

AbstractThis article extends the finite element method of lines to a parabolic initial boundary value problem whose diffusion coefficient is discontinuous across an interface that changes with respect to time. The method presented here uses immersed finite element (IFE) functions for the discretization in spatial variables that can be carried out over a fixed mesh (such as a Cartesian mesh if desired), and this feature makes it possible to reduce the parabolic equation to a system of ordinary differential equations (ODE) through the usual semi-discretization procedure. Therefore, with a suitable choice of the ODE solver, this method can reliably and efficiently solve a parabolic moving interface problem over a fixed structured (Cartesian) mesh. Numerical examples are presented to demonstrate features of this new method.


Author(s):  
T. Philippe ◽  
H. Henry ◽  
M. Plapp

At equilibrium, the shape of a strongly anisotropic crystal exhibits corners when for some orientations the surface stiffness is negative. In the sharp-interface problem, the surface free energy is traditionally augmented with a curvature-dependent term in order to round the corners and regularize the dynamic equations that describe the motion of such interfaces. In this paper, we adopt a diffuse interface description and present a phase-field model for strongly anisotropic crystals that is regularized using an approximation of the Willmore energy. The Allen–Cahn equation is employed to model kinetically controlled crystal growth. Using the method of matched asymptotic expansions, it is shown that the model converges to the sharp-interface theory proposed by Herring. Then, the stress tensor is used to derive the force acting on the diffuse interface and to examine the properties of a corner at equilibrium. Finally, the coarsening dynamics of the faceting instability during growth is investigated. Phase-field simulations reveal the existence of a parabolic regime, with the mean facet length evolving in t , with t the time, as predicted by the sharp-interface theory. A specific coarsening mechanism is observed: a hill disappears as the two neighbouring valleys merge.


2012 ◽  
Vol 323-325 ◽  
pp. 491-496 ◽  
Author(s):  
Soma Prasad ◽  
Aloke Paul

Diffusion couple technique is used to study interdiffusion in Nb-Mo, Nb-Ti and Nb-Zr systems. Interdiffusion coefficients at different temperatures and compositions are determined using the relation developed by Wagner. The change in activation energy for interdiffusion with composition is determined. Further, impurity diffusion coefficient of the species are determined and compared with the available data in literature.


10.30544/308 ◽  
2017 ◽  
Vol 23 (3) ◽  
pp. 197-211 ◽  
Author(s):  
Yuanrong Liu ◽  
Weimin Chen ◽  
Jing Zhong ◽  
Ming Chen ◽  
Lijun Zhang

The previously developed numerical inverse method was applied to determine the composition-dependent interdiffusion coefficients in single-phase finite diffusion couples. The numerical inverse method was first validated in a fictitious binary finite diffusion couple by pre-assuming four standard sets of interdiffusion coefficients. After that, the numerical inverse method was then adopted in a ternary Al-Cu-Ni finite diffusion couple. Based on the measured composition profiles, the ternary interdiffusion coefficients along the entire diffusion path of the target ternary diffusion couple were obtained by using the numerical inverse approach. The comprehensive comparisons between the computations and the experiments indicate that the numerical inverse method is also applicable to high-throughput determination of the composition-dependent interdiffusion coefficients in finite diffusion couples.


2021 ◽  
Vol 27 ◽  
pp. 14
Author(s):  
Andreas Rätz ◽  
Matthias Röger

Standard diffuse approximations of the Willmore flow often lead to intersecting phase boundaries that in many cases do not correspond to the intended sharp interface evolution. Here we introduce a new two-variable diffuse approximation that includes a rather simple but efficient penalization of the deviation from a quasi-one dimensional structure of the phase fields. We justify the approximation property by a Gamma convergence result for the energies and a matched asymptotic expansion for the flow. Ground states of the energy are shown to be one-dimensional, in contrast to the presence of saddle solutions for the usual diffuse approximation. Finally we present numerical simulations that illustrate the approximation property and apply our new approach to problems where the usual approach leads to an undesired behavior.


Sign in / Sign up

Export Citation Format

Share Document