Mathematical Problems in Engineering

Volume 2015 (2015), Article ID 650452, 11 pages

http://dx.doi.org/10.1155/2015/650452

## Computation of Trajectories and Displacement Fields in a Three-Dimensional Ternary Diffusion Couple: Parabolic Transform Method

^{1}AGH University of Science and Technology, Mickiewicza 30, 30-059 Kraków, Poland^{2}Institute of Mathematics, University of Gdańsk, Wita Stwosza 57, 80-952 Gdańsk, Poland

Received 27 February 2015; Accepted 18 May 2015

Academic Editor: Hang Xu

Copyright © 2015 Marek Danielewski and Henryk Leszczyński. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

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.

#### 1. Introduction

The Kirkendall effect in solids involves diffusion and distortion [1]. The distortion of a solid material is the change in position with time of a set of inert markers (real or imagined) fixed in the material. The Kirkendall effect means the changes in the observed marker positions as they move in reaction to the locally nonbalanced diffusion fluxes. During one-dimensional diffusion the material deforms unidirectionally. The markers that are initially fixed in a plane perpendicular to the diffusion direction are likely to move all together, defining marker planes. As shown by Stephenson [1], the deformation that accompanies diffusion is described by a stress-free strain rate.

Essential for the assessment of the Kirkendall effect is the definition of the reference frame (RF) for the components diffusion fluxes (intrinsic diffusion). Darken postulated defining it as a lattice RF [2]. Inert (or imagined) ideal markers do not take part in the diffusion and their velocity equals the local lattice velocity. It implies that intrinsic diffusion fluxes can be defined with respect to this set of markers. Although the Kirkendall effect has been known for 60 years [2, 3], a number of interesting theoretical and experimental effects continue to be revealed which warrant further considerations, including (i) multiple Kirkendall’s planes, (ii) stability of individual Kirkendall’s planes and (iii) trajectory of Kirkendall’s plane in three-dimensional diffusion couples.

Diffusion couples are usually used for experimental and theoretical studies. For a ternary alloy maintained at a constant temperature, the initial state consists of two end members joined with an initial planar discontinuity at a position we take to be . The end members are usually uniform and of different compositions that may also be of different phases. When the couple is held at a fixed temperature, diffusion occurs over time in the -direction. Experimentally, markers are typically placed only on the initial discontinuity of the diffusion couple; however, they can be placed throughout the sample to measure the entire field of deformation field [4, 5]. The rate of change of the position of the entire set of markers in the sample determines the Kirkendall velocity as a field defined throughout the sample.

In the analysis of diffusion couples, there are several quantities of interest, mainly associated with the initial planar discontinuity. The Matano plane is usually defined by proper equal-area constructions established on a concentration profile [6] and remains at for all . Kirkendall’s plane(s) problem from the origin of coordinates at represents the motion of inert markers tracking the deformation; that is, a Kirkendall plane is a material surface. The Matano plane is not a material surface. The Kirkendall plane(s), the material trajectories, and the Matano interface are generally associated with different locations in the sample once diffusion takes place.

This paper is stimulated by a variety of theoretical and experimental efforts. The Kirkendall markers, that is, those markers placed initially in the plane , are occasionally found on two different planes after heat treatment. The possibility of such multiple Kirkendall’s planes was first described theoretically for single-phase diffusion couples by Cornet and Calais [7], using similarity solutions for the composition and Kirkendall’s velocity fields in a binary diffusion couple. A necessary factor for multiple Kirkendall’s planes is a strong variation of intrinsic diffusion coefficients with composition; for example, multiple Kirkendall’s planes can occur if the difference in diffusivities changes sign with composition across the diffusion couple. Stability of the Kirkendall planes was studied by van Dal et al. [5, 8, 9]. They relate stability or instability to the sign, negative or positive, respectively, of the spatial divergence of the Kirkendall velocity at the position of the Kirkendall plane. The criterion established by Cornet and Calais [7] was verified experimentally for solid solutions of NiPd and FePd with single Kirkendall’s planes [5]. Double Kirkendall’s planes were found during interdiffusion in -NiAl [10]. The term “stability” refers to whether markers placed on the plane of the initial discontinuity of a diffusion couple remain sharply concentrated on a plane or whether they become dispersed in the diffusion direction. Computation of the trajectories in solid solutions has been performed by Höglund and Ågren [11].

The parabolic transform method will be presented and used in this paper for a finite-length diffusion couple. The stress generation due to diffusion in constrained geometries will not be discussed in the paper for many reasons. It was shown that in the case, analyzed here (all diffusion fluxes are linearly correlated), it is possible to obtain a symmetrical matrix of transport coefficients for the interdiffusion fluxes (even for a general case of a multicomponent system) and the expression for the entropy production. Thus, the drift generation is a stress-free deformation of the solid in one-dimensional diffusion couples. In complex geometries stresses play an important role and definitely are nonnegligible. It is possible to develop this approach, through more adequate formulations of the free energy functional [12, 13].

The purpose of this paper is to describe the peculiarities of the deformation in three- and one-dimensional binary diffusion couples under various conditions. We compute the trajectories (position dependent on time) of a full set of markers, focusing on cases with multiple Kirkendall’s planes and sole Kirkendall’s planes that coincide with the moving initial contact interface. Numerical experiments are compared with Belova et al. [14].

#### 2. Model

##### 2.1. Diffusion and Material Velocity

To permit applications to nonideal solid solutions, we present expressions for the standard Darken approach in terms of gradients of potentials of the components in the three-dimensional space . The mass conservation laws have the formwhere and are the component concentrations and their fluxes measured in the laboratory RF. The most convenient RF in the case of a finite (bounded) diffusion couple considered in this work is, for example, the boundary of the diffusion couple.

The central Darken postulate is defining the component diffusion fluxes via a drift velocity field that is measured in the laboratory reference frame. It is also called* convection velocity* [1],* Darken velocity* [15],* lattice velocity* [2], or* material velocity* [16, 17]. It characterizes the solid solution motion and equals the velocity of a set of markers fixed in the material. If the material is crystalline, the markers move with the lattice, and the process occurs due to the existence of vacancy sources and sinks, for example, one- and two-dimensional defects. The diffusion fluxes are called intrinsic (or lattice) fluxes and are related to those measured in the laboratory frame by the formulaThe overall molar volume and mole fractions are defined bywhere and is the th component partial molar volume. From (3), we getIn the multicomponent solid solutions the following relations hold:For isothermal processes, when all partial molar volumes are constant and stress effects are negligible, the relations (5) reduce toVariable partial molar volumes in binary systems were treated by Sauer and Freise [18] and are not considered here.

##### 2.2. Kirkendall Velocity and Displacement

Multiplying the conservation law for each species, (1), by matching partial molar volume, results in the volume continuity formulaThe intrinsic fluxes, governed by the Nernst-Planck formula, are assumed to be proportional to the gradients of their chemical potentials (diffusion potentials in general [19]):where the intrinsic (or lattice) mobilities may depend on the concentration. From relations (2), (3), (7), and (8), the Darken velocity field ( Darken’s drift) is given byIf the velocity field is not one-dimensional, additional equations might be required to determine the flow [16]. The mass conservation law can be written aswhere the drift velocity fulfills (9). In the following, we will consider a ternary system () and constant partial molar volumes such that for . Consequently, using the relations from (3), (9) and (10) becomeThe Darken velocity field characterizes a deformation of the material during diffusion. If one considers diffusion in a crystalline material, the markers are positions in an imperfect lattice in which production and/or annihilation of vacancies occur at rates that are sufficient to maintain an equilibrium chemical potential of vacancies everywhere; that is, the local thermodynamical equilibrium is sustained in the material [20].

The same equations apply to diffusion in fluids and amorphous solids, where the drift velocity describes the convective motion of the substance as indicated by the motion of small inert particles in a fluid as a result of different diffusion mobilities. The effective diffusion coefficient in gases, in the manner of Darken [2], describing the net convective motion, was developed already in 1899 in the kinetic theory of gases [20].

##### 2.3. Diffusion Potentials and Diffusion Coefficients

The term* diffusion potential* was introduced by Larché and Cahn [21]. They derived a potential as a sum of the elastic energy and chemical potential. In general, the diffusion potential is a sum all potentials that account for the storage of energy due to the addition of a diffusant [19]. The regular solutions can be treated by analogous formalism. In this work we focus on an ideal solid solution:where is the reference potential; is the gas constant; is the temperature. From relations (3), (8), and (12), the intrinsic diffusion fluxes are defined byEquations (11) becomeFrom (14) in the closed and bounded system the drift velocity is given by [22]Upon denoting , , and , from (15) and (16) we havewith . Values of belong to the interval . The initial conditions , , and are imposed for . Boundary Neumann conditions are given byon the lateral boundary , , where is the outward normal vector and is the operator of normal derivative; see [23].

#### 3. Reduction to Parabolic Equations in

In three-component models with three distinct diffusion coefficients the third variable can be expressed in terms of the other variables , . DenoteAssume that , which means that is the smallest diffusion coefficient (it is convenient to assume that ). Thus, the ternary system (17) reduces to the following:with , on the lateral boundary , and and for . We reduce this to a parabolic system by virtue of a change of dependent variables. Takeand denotethenSimilarly, one calculatesfor andfor where , and are the second-order derivatives of in . Consider the parabolic operatorwhere “lower order terms” will be specified later. The second-order terms vanish, provided that is an eigenvalueand is an eigenvector. The characteristic polynomialhas positive and distinct roots, because ; thusand , whereThe eigenvalues are equal toTheir eigenvectors fulfill the algebraic equationHence, we have the hyperbolic PDEIts characteristicshave solutions . Thus, we get two families of curves . for and . for . Introduce new dependent variablesand the inverse functionswherever they exist. We find their gradients from the systemwhich can be inverted aswhere and for . The ultimate parabolic PDE system looks as follows:for , where . The right-hand side of (39) consists of the “lower order terms” from (26).

#### 4. Trajectories in

An arbitrary smooth trajectory is governed by the ordinary differential problemwhere depends on and ; confer [23]. Recall that . It is obvious that trajectories of material points are governed by ordinary differential equations. The only question is as follows: how to derive these ODEs? We draw the reader’s attention to the microscopic setting based on stochastic diffusion equations in [24].

##### 4.1. Material Trajectories

In Darken’s method the function is equal to the Darken (material) velocity :Suppose that the marker at is placed at a surface for . We want to trace its trajectory for which constitutes planes for . Since the characteristics are governed by (40), fulfills the hyperbolic initial-value problembecauseIn particular, Darken’s trajectory fulfills the equationand Kirkendall’s plane satisfies the PDEwhere ; see [23].

The Kirkendall plane given by (45) with characteristics coinciding with Darken’s trajectories (44) fulfills the following elegant property.

Proposition 1. *For any regular region with the boundary defined by where is the solution of the PDE (45), or equivalently being constituted by the characteristics obeying (44) starting from the initial plane , the total mass of each component plus its overall inflow through the boundary is constant:**Analogous identities are valid for and .*

*Proof. *Denote by the left-hand side of (46). Then, we haveWe differentiate as follows:Since , we have . This proves the assertion.

This property proves that in ternary and higher solid solutions the Kirkendall plane allows finding the Matano plane that is identical for all components.

##### 4.2. Displacement Field

Since the velocity field is described in terms of laboratory coordinates, the displacement field is a function of the laboratory coordinates and time as well. The relation between the reference and laboratory coordinates is defined by a mapping which gives a reference (initial) position of the marker that is currently at position ; the inverse mapping is providing the current position of the marker initially placed at position . The displacement field satisfies the relationThat is, is the difference between the current and initial positions of the marker which is now at position in the laboratory frame. The displacement and velocity fields satisfy the equation

##### 4.3. Component Trajectories

Consider any regular subdomain in with a smooth boundary . If we are interested in the evolution of starting from then we obtain a family of subdomains described by the equations . The mass balance leads to the integral identityThe integration by parts leads to the formulawhere is the unit normal vector, proportional to , and is the surface element. A similar relation holds for and . We can get the evolution of a region which contains some amount of :

*(a) Balance of Mass in 1D.* Consider the one-dimensional case with the central position separating substances at . Thus, . Suppose that the separating position moves according to the mass-balance lawwhich means that the total mass of is preserved on the left of the separating position . If we differentiate the identity (54) with respect to we obtainTaking into account the governing equations, we gethence,In this case, we arrive at the equation of the separating position trajectoryor equivalentlyThe respective plane is a solution of the PDE

*(b) Tracer Trajectory.* Let us return to the trajectory of which starts from a spike-shaped droplet. Suppose that we measure the middle position of this marker; hence, we have the equationAlike (54), we arrive at the differential equationThe ideal marker means a fictitious observer (indicator) within solid moving with the local drift velocity [2] and not altering the diffusive mass flow. In practical realizations nonideal markers affect diffusion for the reason that they are often metal wires or submicron inclusions of a stable oxide embedded in the diffusion zone. The method presented by Belova et al. [14] and analyzed in the present paper has an advantage of a negligible impact on the diffusive mass transport (when concentrations of marker atoms are low). An obvious drawback of both methods (Belova’s and ours) is diffusion of the marker itself. This process can be to some extent delayed by an introduction of marker atoms whose diffusivity in the material is low; see Figure 1. Observe that the trajectory of is very close to the ideal marker, since the maximum of is attained at a point such that . We draw a useful conclusion: . We can also derive similar equations to (62) for and . Having in mind a severe singularity at , we can analyse the front movement of and as follows. Suppose that is initially concentrated on the left . We study an evolution defined bywith some small (to avoid the singularity). Hence, we get the differential equationSince is decreasing, we have and moves to the right. In an analogous way initially placed at moves to the left according to the following: