Open access peer-reviewed chapter - ONLINE FIRST

# Using the Principles of Nonequilibrium Thermodynamics for the Analysis of Phase Transformations in Iron-Carbon Alloys

By Bobyr Sergiy Volodimyrovych

Submitted: May 18th 2018Reviewed: December 19th 2018Published: May 30th 2019

DOI: 10.5772/intechopen.83657

## Abstract

Using the principles of nonequilibrium thermodynamics, a technique has been developed for calculating diffusion flows during phase transformations in iron-carbon alloys. Expressions for the calculation of cross coefficients, driving forces, and flows in Onsager equations for the model thermodynamic system are given; examples of the use of the developed technique are given for the processes of graphitization and the formation of carbides in chromium steel during tempering. The nonequilibrium thermodynamics analysis of the eutectoid transformation is executed into carbon steel. Onsager’s equations of motion are built for the model thermodynamics system describing eutectoid transformation. The basic kinetic parameters of process are growth rate of perlite and between inter-plates distance for the stationary process of eutectoid transformation. We founded dependencies of basic kinetic parameters of process from the size of supercooling. A nonequilibrium thermodynamic model of the austenite nondiffusion transformation in iron and alloys based on it is developed, taking into account internal stresses in the system. Onsager motion equations are found for a model thermodynamic system describing a nondiffusion transformation and kinetic equations for changing deformations and growth rates of the α-phase. A scheme of austenitic nondiffusion transformations is constructed, including normal and martensitic transformations, as limiting cases.

### Keywords

• nonequilibrium thermodynamics
• the iron-based alloys
• transformation of austenite
• diffusion
• equations of motion
• nondiffusion transformation

## 1. Introduction

The study of phase transformations is one of the most important problems in the physics of metals [1, 2, 3]. Phase transformations are divided into diffusion and nondiffusion [1]. If the kinetics of phase transformation in steels and cast irons is determined by the diffusion of carbon, this allows them to be attributed to conversions controlled by diffusion [1, 2, 3, 4]. Such transformations in iron-carbon alloys include pearlitic transformation of austenite, and transformations occurring during tempering, graphitization of undoped cementite, separation of carbides in alloyed steels, and others [4, 5, 6].

When the rate of transformation of austenite is determined by the rate at which the interface separates, differing only in its crystalline structure, the transformation is called nondiffusion [1]. Kinetically, the normal polymorphic and martensitic transformations of austenite are distinguished. When the temperature of the normal transformation decreases, its velocity first increases and then decreases. The kinetics of the martensitic transformation is characterized by a very high rate of growth of individual crystals and the maximum space velocity at the initial moment of transformation under isothermal conditions.

In addition to martensite, at least two other structural components are known, which are formed with a shear (“martensitic”) morphology of crystal formation—ferrite side-plates and acicular ferrite. They can also be attributed, with some simplifying assumptions, to the products of the nondiffusion transformation of austenite. In addition, in some alloys martensitic and normal transformations occur at the same temperature [1]. The consistent theory of nondiffusion transformations should explain this phenomenon. Thus, the theoretical description of the processes of phase transformations in iron-carbon alloys is a complex and urgent task of modern metal physics.

Nonequilibrium thermodynamics provides the necessary apparatus for analyzing the processes of phase transformations in iron-carbon alloys [7, 8, 9]. In the general case, the thermodynamic equations of motion have the form [7]:

Ji=к=1NLikXki=1..N,E1

where Ji are flows, Xk are the thermodynamic forces, Lik = Lki are the Onsager kinetic coefficients [9], and i, k are the charge numbers (transfer substrates).

The main driving forces of phase transformations in nonequilibrium thermodynamics are gradients of the chemical potentials of their components [6, 7, 8, 9]. When discontinuous systems are considered, the finite differences of chemical potentials (−Δμi,) as the transition from a metastable state to a stable state are used as thermodynamic forces [10, 11]. Equations of nonequilibrium thermodynamics were first used in the physics of metals to describe the process of graphitization of nonalloyed iron-carbon alloys [6, 11].

As is known, unalloyed cementite in iron-carbon alloys at normal pressure is a metastable phase, its activity in phases with it in equilibrium exceeds the solubility of graphite, a stable phase [11]. Therefore, at a sufficiently high temperature, graphitization of such alloys takes place, that is, phase transition from metastable to stable equilibrium. Despite the seeming simplicity of this process, its theoretical description is a complex task.

If two values are used as charges of the graphitization process-carbon and iron concentrations, then, according to (1), the equations of motion take the form:

J1=L11Х1+L12Х2E2
J2=L21Х1+L22Х2,E3

where J1 is the carbon flow characterizing the rate of the graphitization process, J2 is the flow of iron, and X1 = (−ΔμFe) and X2 = (−ΔμC) are the thermodynamic forces of iron and carbon. The potential drop has a “+” sign as it increases, and the flow is directed toward a decrease in the potential, so the expressions for the forces contain the sign “.”

The main question that must be solved when using the Onsager Eqs. (1)(3) is the values of the cross coefficients.

In [5], for the first time on the basis of a special variational procedure, an expression for the cross coefficients in the Onsager equations was proposed in the form:

L21=L12=L11×L22,E4

and the sign–before the root is chosen on the basis that the observed flux of iron with respect to the flow of carbon had a negative sign.

As shown in [6, 11], in the complex process with two flows, an increase in the potential of one of the charges is observed, that is, one process is “leading,” and the other is “driven.” The “driven” process in itself, i.e., in isolation from the “lead,” is not possible, since thermodynamically not beneficial. In the system of Eqs. (2) and (3), the thermodynamic force (−ΔμFe) is negative and inhibits the process as a whole, the diffusion of iron is a forced process, and the leading one is the diffusion of carbon.

Thus, the graphitization process must be accompanied by a very intensive transfer of a solid solution (mainly iron), which makes it possible for the phase with a low-density graphite to grow in it. The authors of [6, 11] assumed that the factor contributing to graphitization is the pressure that arises in the austenite matrix under the action of graphite inclusions that expand it. However, in [12], considering the mechanism of graphitization of cast irons during thermocyclic treatment, K.P. Bunin with AA. Baranov came to the conclusion that the absolute value of the contact pressures is an order of magnitude less than the necessary for the dislocation creep mechanism under the influence of contact pressure. Since graphite films in pores cannot possess super strong properties, the evacuation of matrix atoms is apparently carried out by another mechanism.

In [5], using nonequilibrium thermodynamic methods, it was shown that under the conditions of the system’s striving for dynamic equilibrium, the concentration of vacancies in graphite inclusion becomes less than the vacancy concentration at the γ-phase-graphite boundary. This can occur as a result of approaching the γ-phase boundary—graphite of austenitic vacancies. In this case, the thermodynamic force (−Δμ0v) prevents the graphitization, and the reduced graphitization force (−Δμ*С) decreases to zero and can even take a negative value.

Thus, the goal of this paper is to show how the methods of nonequilibrium thermodynamics can be used fruitfully to solve the theoretical problems of metal physics, namely, the analysis of phase transformations. Let us further consider the application of the principles of nonequilibrium thermodynamics to the analysis of specific cases of phase transformations in iron-carbon alloys.

## 2. Formation of carbides in chrome steel during tempering

Consider the process of separation of carbides in a low-carbon steel system of iron-carbon-chromium with 0.15% carbon and about 5% chromium at 600°C. In this model system, there are two phases—the doped α-phase (F) and carbides (K), in which carbon, iron, chromium, and vacancies flows (Figure 1). As charges, we will use four quantities—the concentrations of carbon, iron, chromium, and vacancies. The flow of vacancies in the carbide phase will be assumed to be equal to the flow of vacancies in the ferrite.

In the absence of a change in the volume of the system, for flows in the doped α phase, condition [13] is fulfilled:

JFe+JCr+Jv=0,E5

so one of the threads (in our case—Jv) is a dependent quantity. According to (1), the thermodynamic equations for flows in the carbide phase take the form:

JFe=L11ΔμFeL12ΔμcL13ΔμCrE6
JС=L21ΔμFeL22ΔμcL23ΔμCrE7
JCr=L31ΔμFeL32ΔμcL33ΔμCr,E8

where JFe, JС, and JCr are the flows of iron, carbon, and chromium, respectively.

Based on the general principles of nonequilibrium thermodynamics, we can find the values of the thermodynamic forces −∆μFe, ∆μCr, and −∆μС, as well as the values of the kinetic coefficients L12, L13, and L23, as it was previously performed in [5] for a system with two flows. In the conditions of complete equilibrium, ∆μFe = 0, ∆μС = 0, and ∆μCr = 0. However, for a linear thermodynamic system, there is also the possibility of dynamic equilibrium, in which all the flows are 0, but some thermodynamic forces in the system are not equal to zero (there are their variations) [5, 7].

Let us consider this possibility for a triple thermodynamic system. From Eqs. (6)(8), it follows that near equilibrium, in the presence of variations of thermodynamic forces, the following conditions must be fulfilled:

JFe=0L11δμFe+L12δμС+L13δμCr=0,E9
JС=0L21δμFe+L22δμС+L13δμCr=0,E10
JСr=0L31δμFe+L32δμC+L33δμСr=0,E11

where the index δμ denotes the coordinated variations of the thermodynamic forces that ensure the dynamic equilibrium of the system. It follows from the system of Eqs. (9)(11) that the expressions for the flows of iron, chromium, and carbon are connected: the cross rates L12, L13, and L23 in expressions for the flows must have values such that the determinant of the matrix A composed of the coefficients of this system was equal to 0. In this case, the values of the flows of iron and chromium can significantly increase due to cross-kinetic coefficients in comparison with the independent diffusion of these elements [7, 16].

Let us find the expressions for the cross coefficients, which make it possible to obtain a nontrivial solution of the system of Eqs. (9)(11). From the first Eq. (9), we establish a connection between the variations of forces:

δμFe=L12/L11δμСL13/L11δμCr.E12

Substituting (12) into Eqs. (10) and (11), we find

JС=L22L122/L11δμС+L23L12L13/L11δμCr=0,E13
JСr=L32L13L12/L11δμС+L33L132/L11δμСr=0,E14

For independent variations δμС and δμСr, the linear system of Eqs. (13) and (14) is compatible if the coefficients of δμС and δμСr are equal to 0, from which we immediately find the relation between Onsager’s kinetic coefficients:

Lik=Lki=±Lii×Lkk,i,k=13E15

and the sign before the root is selected based on the sign (direction) of the flows under consideration (see Figure 1). The considered procedure of variation allows us to find cross rates in the Onsager equations after the direct kinetic coefficients are calculated. In this case, the established connection (15) is satisfied for systems not very far from equilibrium and for the real system is approximate.

## 3. Calculation of thermodynamic forces and kinetic coefficients

Let us find the values of the thermodynamic forces and kinetic coefficients for the steel of the Fe-C-Cr system with 0.15% C at 600°C. We will assume that in a solid α-solution, there is chromium with a concentration of СCr = 0.05 and a carbon with a concentration of СC = 0.007, an iron concentration of СFe = 0.943. In cementite-type carbide, chromium is found with a mass fraction of ∼20% (with a concentration of СCr = 0.2) and carbon with a carbon concentration of 0.25, an iron concentration in the carbide C’Fe = 0.55.

It is known from the experimental data that carbon is removed very rapidly (approximately 1 minute) from the α-solution of alloyed steel at a temperature of 550–650°C and, consequently, the formation of carbide inclusions is primarily due to carbon diffusion [14].

The thermodynamic force for carbon can be calculated from the formula [11]:

ΔμС=RTlnаСКаСα,E16

where аαС is the thermodynamic activity of carbon in α-solution, аКС is the thermodynamic activity of carbon in cementite, R is the universal gas constant, and T is the temperature of the alloy.

The change in the thermodynamic activity of carbon in the alloy upon doping with component i can be found by the method of [15, 16] from the equation:

lnаС/аС0=βiNi,E17

where βi is the coefficient of the element’s influence on the thermodynamic activity of carbon in the alloy, Ni is the content of the element in the alloy in atomic fractions, and аС0 is the thermodynamic activity of carbon for the alloy in the standard state.

We will assume that for our steel in the standard state аαС0 =аКС0 = аС0, i.e., unalloyed cementite in steel with 0.15%, C is stable and in equilibrium with the solid solution at a tempering temperature of 600°C [13]. Using this condition and Eqs. (16) and (17), we find:

lnаКС/аαС=βКCrNКCrβαCrNαCrE18

The value of βi is calculated through the interfacial distribution coefficient of the alloying element Ki = Ni (K)/Ni (α) and the atomic fraction of carbon in the alloy Nc [15, 16]:

βi=Ki1+NcКKiNcαKi1Nc+NcКKiNcα.E19

With a slight error for low-alloyed alloys, we can take Nc (K) = 0.25, Nc (α) ≈ 0.001—the carbon content in the undoped phases of steel, taken from the Fe-C state diagram.

Using the coefficient of chromium distribution between the α-phase and the carbide KCr, equal to 4, we find the equations for calculating the coefficients of influence βCr:

βСr = −3.246/(3,0Nc + 0.246),E20

whence βαСr = −12.16 and βКСr = −3.26.

Then from expressions (16)(18), one can find the values.

lnаКС/аαС=0.6085+0.652=0.0425andμС=308.47Joule.E21

The work done in the diffusion of carbon from the α-phase to cementite is positive. For the diffusion of iron, it is not possible to calculate the difference of thermodynamic potentials, since the coefficient of iron activity in carbide is unknown. However, from the experimental data and the thermodynamics of the process, it is known that diffusion of carbon is the leading one, the diffusion of chromium accompanies the diffusion of carbon, and the diffusion of iron is forced, since it is directed toward increasing the concentration of iron.

With this in mind, we find the values of the kinetic coefficients Lii in the Onsager equations.

As is known [8, 13], the kinetic coefficients Lii are related to the diffusion coefficients Di by the relation:

Lii=СiDi/RT,E22

where C1 is the concentration of iron in the alloy (0.943), C2 is the concentration of carbon in the alloy (0.007), and C3 is the concentration of chromium in the alloy (0.05).

Dependences of the diffusion coefficients of chromium and carbon in doped chromium ferrite on the temperature have the form [14, 17]:

DСα=0,177exp88230RTсм2/сек,E23
DFeα=2,9104exp251000RTсм2/сек,E24
DСrα=3,05exp358000RTсм2/сек.E25

At a temperature of 600° C:

D1 = DαFe ≈ 3.03·10−19 cm2/s, D2 = DαC ≈ 1.02 10−6 cm2/s, and D3 = DαCr ≈ 1.38 10−21 cm2/s.

Using relations (23)(25) and (15), we find the values of the kinetic coefficients for our system:

L11 = 0.394 × 10−22, L22 = 0.984 × 10−13, L12 = −1.97 × 10−17, L33 = 0.95 × 10−26, L13 = −0.611 × 10−24, and

L23 = 0.306 × 10−19. Consequently, the system of Eqs. (6)(8) takes the form:

JFe=0.394×1022ΔμFe1.97×1017ΔμС0.611×1024Δμсr,E26
JС=1.97×1017ΔμFe+0.984×1013ΔμС+0.306×1019ΔμCr,E27
JСr=0.611×1024ΔμFe+0.306×1019ΔμС+0.95×1026ΔμCr.E28

It follows from Eqs. (26)(28) that the values of iron and chromium fluxes increase substantially due to the cross-ratios L12 and L32 of a significant thermodynamic force (−ΔμC). The value of the carbon flux having a positive sign is determined mainly by the intrinsic coefficient L22. The thermodynamic forces of iron and chromium make an insignificant contribution to the fluxes, because of the small value of the kinetic coefficients and their influence can be neglected. Then, as direct calculations show:

JFe = − 6.08 × 10−15, JC = 3.04 × 10−11, JCr = 0.94 × 10−17, and

JV = 6.07 × 10−15 cm2/s.

It was established in [18] that during the tempering period, a certain amount of nanoparticles of special chromium carbide with a size of ∼100 nm can be formed in the steel, which were detected experimentally.

## 4. The nonequilibrium thermodynamics analysis of the eutectoid transformation

In [19], a generalization of the equations characterizing the growth of the pearlite colony is proposed, based on the application of nonequilibrium thermodynamic methods.

To this end, Eq. (19) from [20], which characterizes the growth rate of a perlite colony, is represented as:

dXdt=DхСф'Сц'/Сф'Сф+Сф'Сц'/СцСц'/=Dх/φ,E29

where Dx is the carbon diffusion coefficient in austenite along the x axis at a given temperature T, Δ is the thickness of a layer of austenite with different concentration of carbon, C’Φ and С’y is the carbon concentration in the austenite near the ferrite and cementite plates, respectively, at a temperature T (Figure 2), Сy is the carbon content in cementite (∼6.67%), is the carbon content in the ferrite at a given temperature T, and −Δφ is the thermodynamic force of perlite lateral growth. It is determined by the carbon concentrations in ferrite and cementite and has a dimensionless value.

The second equation characterizing our system—the heat balance Eqs. (23)(25) from [20]–is written in the form:

Cγ/dt=αТqγ//dtE30

where α is the heat transfer coefficient, C is the specific heat, ΔТ is the temperature difference between the sample (T) and the cooling medium, q is the specific amount of heat released during the formation of perlite, and γ is the density of steel.

If two quantities are used as charges for the eutectoid transformation of austenite-the temperature of the sample T and the thickness of the plates of perlite X, then, according to (4), the Onsager motion equations must have symmetric forms (2) and (3),

where J1 = − dХ/dt is the flow of the pearlitic layer (with increasing absolute value of the thermodynamic growth force of perlite, the flow increases in absolute value), and J2 = − Cγ/dt is the heat flow in the sample (with a drop in sample temperature, the flow is positive), Х1 = (−∆φ), Х2 = (−∆Т/Т) is the thermodynamic forces of perlite growth and temperature [14].

In order for Eq. (29) to correspond to Eq. (2), it must contain an additional term L12 (−∆Т/Т); with the value of the coefficient L11 = (/∆):

J1=Dх/φ+L12Т/ТE31

where L12 is a cross ratio whose value is not yet known. Thus, we introduce (we assume) an additional dependence of the growth rate of the perlite layer not only on the carbon concentrations in the phases but also on the temperature.

Substituting expression (34) into the energy balance Eq. (33), we find the expression for the heat flow J2:

J2=qγ/Dх/φ+L12Т/ТαТ∆Т/Т=qγDх/2∆φ+qγL12/αТ∆Т/ТE32

Relating Eqs. (3) and (32) to each other, we obtain:

L21=qγDх/2E33

Using for the kinetic coefficients, the Onsager reciprocity relations Lik = Lki [9], we find that

L12=L21=qγDх/2,E34

whereas

L22=q2γ2Dх/3αТ.E35

The system of Eqs. (31) and (32) takes the form:

J1=Dх/φ+qγDх/2Т/ТE36
J2=qγDх/2φ+q2γ2Dх/3αТ∆Т/Т.E37

In accordance with (36), the perlite growth rate is affected not only by the concentration thermodynamic force, but also by the temperature difference between the sample and the environment. Let us further consider the phase transformation of austenite under special conditions of steady growth of the pearlite colony, when it can be assumed that ΔT ≈ сonst, /dt ≈ 0. In this case, Eq. (37) takes the following form:

J2=qγDх/2φ+q2γ2Dх/3αТ∆Т/Т=0.E38

For small ΔT, we can write approximately, as was done in [20]:

φ=к∆Т/Т,E39

where k is the proportionality coefficient.

By analogy with the previously obtained solutions [21], we introduce the following notation:

Δ0=кqγDx/αТ;E40
whereΔl=/kis the characteristic parameter of the system.E41

Eq. (38) can now be represented in the form:

3Δ02Δ02Δl=0E42

For Δl = 0, as expected, the solution of Eq. (31)Δ = Δ0. We obtain a well-known solution for the pearlite transformation of austenite [20]. In the real domain, there is one solution of Eq. (44). For small Δl (<0.5), the root Xk is in the region close to 1 (Xk → D0), with increasing Δl (in units of D0), the root value increases. For large values of Δl, the root of Xk is approximately equal to

ХkΔl/Δ03.E43

The between interplate distance of perlite for a stationary growth process is found from the formula:

S0=2Хк×Δ0=2кХkΔlDx/αТ.E44

Using Eqs. (36) and (43), (44), we also find an improved expression for the perlite growth rate for an isothermal transformation

dXdt=kDxS0ΔTT1+2ΔlS0E45

The formula (45) is a more precise expression for determining the growth rate of perlite in the eutectoid transformation, than the expression obtained earlier by the authors of [20].

We use the well-known dependence of the diffusion coefficient on temperature [17]:

D = A exp.(−Q/RT),E46

где Q is the activation energy, (Q ≈ 134 кJ/mol), and R is a constant (R = 8314 J/(mol⋅К)).

After substituting the known values of the steel parameters and taking into account that to 2.0, we find the calculated dependence of the perlite growth rate on the supercooling value of the alloy (Figure 3). In this figure, the dependence of the perlite growth rate on the supercooling value, calculated according to Zener’s formula (1) [22, 23], is given for comparison.

According to the constructed model, the perlite growth rate in the direction of the X axis has a maximum value at supercooling ΔТ = 140.0°С. The perlite growth rate calculated according to Zener’s formula has a theoretical maximum value at overcooling ΔТ = 96.0°С. Consequently, the theoretical expressions (31) and (32) make it possible to describe with greater accuracy the maximum and the course of the experimental curve for the perlite formation rate presented in [3, 24] for high-purity eutectoid steel.

The expression for perlite growth rate obtained in this section has a significant value at supercooling of 300–400°С, thereby determining the possibility of perlite formation in this temperature range. Indeed, the formation of perlite in carbon steels in the temperature range 375–325°С was revealed in [24].

The calculated dependence of the between interplate distance of perlite by formula (46) on the magnitude of the supercooling of steel is shown in Figure 4. The same figure shows the experimental points from [24].

A fairly good agreement of the calculated dependence with the results of the latest experiments is observed, which indicates the adequacy of the proposed model.

## 5. Application of the positions of nonequilibrium thermodynamics to the analysis of the nondiffusion transformation of austenite

Martensite is the basis of hardened steel, so studying the mechanism and kinetics of its transformation is still of extreme interest for the theory and practice of heat treatment.

In the works of G.V. Kurdyumov and coworkers, the martensitic transformation is considered as the usual phase transformation in a one-component system, further complicated by the influence of a strong interatomic interaction, which leads to the development of significant stresses in the martensite crystal and matrix [25].

In accordance with the alternative mechanism, the martensitic transformation takes place by means of an instantaneous shift of atomic planes that does not require thermal activation and is not associated with thermodynamic transformation stimuli [1], [26]. In this case, the stress initiating the transformation is believed to be the stresses arising from the sharp cooling of the sample (quenching) [26].

Considering the martensitic transformation as a thermally activated process, B.Ya. Lyubov used the equations of normal transformation obtained on the basis of the positions of nonequilibrium thermodynamics to describe his kinetics [3].

Changes in a complex or composite system under constant external conditions can be described as the process of increasing entropy. The rate of increase of entropy σ can be represented as the sum of the flux products and the corresponding forces for all transfer substrates in an amount of N [7, 8, 9, 10]:

σ=dSdtirrev=к=1NJkXkk=1..N,E47

In the general case, the flows can be represented in the form (1).

The irreversible change in the entropy dSirrev is equal to the sum of entropy changes in the system and the environment:

dSirrev=dS+dSeE48

Under isothermal conditions, when the released heat is absorbed by the environment and the temperature remains constant:

dSe = − dQ/T, dQ = dU + PdV
dSirrev=dSdUa+PdVT1=TdSdUPdVT1.E49

Since dU + PdV–TdS = dG, and if we take into account the low compressibility of bodies in the condensed state and relatively small pressures, then.

dSdtirrev=Т1dGdtТ1dFdt,E50

where F is the free energy of the system.

The change in free energy in a system with a variable number of particles and internal stresses can be represented in the form [3], p. 142:

dF=dFε+dFn=σikik+φldnl,E51

where dFε is the change in free energy in the system related to internal stresses, dFn is the change in the free energy in the system, determined by the variable number of particles of type l, σik is the stress tensor, εik is the strain tensor of the system, φl is the chemical potential of the lth element of the system, and nl is the number of particles of the lst element of the system per unit volume, l = 1, N.

We now introduce some simplifying assumptions. First, for the nondiffusion transformation of austenite, only one kind of particles, the α-phase of iron nα, will be taken into account. Approximately, this is also true for alloys of iron with close elements (nickel, chromium, cobalt). Of course, φ is some effective (averaged) chemical potential of the atoms of the alloy.

Secondly, we assume that the deformation of the system is a triaxial compression-expansion, and in the expression for dFε, only the diagonal components of stress and strain tensors are taken into account:

σik=εik=0,ik.σii=σ,εiiεE52

The change in internal energy can then be represented as:

dF=3σdε+φdnαE53

but the change in entropy:

dSdtirrev=Т13σdεαdt+φdεαdt.E54

Thus, in our system, in addition to the particle flow from the γ phase to the α-phase of J1 = dnα/dt, we will also take into account the change in the strain of the sample with time J2 = dε/dt. These flows are related to the driving forces by the Eq. (1).

If, as charges of the process of nondiffusion transformation of austenite, the two quantities are the concentration of α-phase particles and the strain value, then, according to (1), the equations of motion take the form:

J1=L11Х1+L12Х2E55
J2=L21Х1+L22Х2,E56

where X1 = Δφ is the thermodynamic force for iron, the change in the chemical potential at the transition of particles from the γ-phase to the α-phase, and X2 = Δσ is the change in the internal stress during the transition from the γ-phase to the α-phase.

The system of Eqs. (55) and (56) describes the contribution of stresses and deformations to the nondiffusion transformation of austenite. However, we do not yet know the coefficients of the equations in it. We now find expressions for the coefficients of the system of Eqs. (55) and (56). The coefficient L11 characterizes the normal transformation:

J1=dnαdt=L11φγφαE57

In the normal kinetics of the phase transformation, the formation of the center (particle) of the α-phase occurs through separate (independent) acts of detachment of particles from the γ-phase and the attachment of atoms to the ferrite center. If we consider the process of formation of an α-phase close to the process of self-diffusion of iron in the γ-phase, then the coefficient L11 has the form [13]:

L11=DγRTE58

where is the self-diffusion coefficient of iron in the γ-phase (or the effective coefficient of self-diffusion in the γ-phase of the alloy), T is the transformation temperature, and R is the gas constant [27].

The self-diffusion coefficient of iron is taken in the usual notation [17]:

Dγ=D0eUkT=4.58·104ехр252000/RTE59

where D0 is a multiplier and U is the activation energy of diffusion.

The coefficient L22 characterizes the direct relationship:

J2=dεαdt=L22σγσα.E60

Let σγ = 0. Let us take into account that for triaxial compression stretching [28]:

where σα is the stress in the α phase and σγ is the stress in the γ phase.

Let σγ = 0. Let us take into account that for triaxial compression-stretching [28]:

σα=E312μΔVV=E12μεα,E61

where E is the modulus of elasticity of steel (∼2.17·105 МPa) and μ is the Poisson ratio (∼ 0.26).

Then, expression (60) can be transformed as follows:

dεαdt=L22σα=L22E12μεα=vLε,E62

where the following values are entered:

v is the propagation velocity of the microdeformation in sample (∼1000 m/с) [3] and L is the characteristic distance over which the microdeformation of the shear is propagated (the size of the martensitic strips or plates). At the initial stage of the formation of the shear structure, it has a magnitude of the order of the diameter of the austenite grain (∼ 100 μm), and then decreases with decreasing temperature [1].

From Eq. (62), we find that the coefficient L22 is equal to:

L22=v12μ.E63

The cross-coefficients L12 = L21 for a nonequilibrium thermodynamic system are found with sufficient accuracy by the formulas proposed in [5]:

L12=L11L22=DγRTv12μE64

Thus, we obtained simple differential equations for a nonequilibrium thermodynamic system describing the nondiffusion transformation of austenite taking into account the influence of internal stresses.

Let us write the equations of motion of our system in the form:

dnαdt=L11φ+L12σγL12σα.E65
dεαdt=L21φ+L22σγL22σα.E66

We first transform Eq. (66) taking into account expression (62). We have:

dεαdt+α/L=L21φ+L22σγ,E67

where εα is the magnitude of deformations of the α-phase. The differential Eq. (66) with constant coefficients (temperature) has a solution:

εα=L21Δφ+L22σγ12μL22Е1evLt.E68

This kinetic equation describes the change in the magnitude of the deformation of the α-phase in time. At t = 0, εα = 0. When the time is counted, a fast (∼ 10−6 s) process of transition to deformation occurs:

εα=εγ+Δεα=εγ+L21Δφ12μL22Е.E69

Eq. (69) shows that the residual deformation of the α-phase after the transient process consists of the deformation of the austenite εγ and the additional deformation Δεα. This additional deformation determines the change in the volume of the sample as γ → α-transformation:

ΔVγαVγ=3nαΔεα.E70

Then, substituting expression (53) into Eq. (51.1), we find:

dnαdt=L11φ+L12σγL12L21Δφ+L22σγL221evLt=L11φ+L12σγL11Δφ+L12σγ1e=vLt=L11Δφ+L12σγevLt.E71

It can be concluded from expression (71) that the growth rate of α-phase particles depends on the stresses in the γ-phase. The greater the value of tensile stresses in the γ phase, the higher the growth rate of ferrite particles. The rate of growth of the α-phase particles at a constant temperature very rapidly (exponentially) decreases in time, determining the incompleteness of the transformation.

Integration of Eq. (71) with time-independent coefficients L11 and L12 allows us to obtain the kinetic equation for nα:

nα=L11Δφ+L12σγLv1evLt.E72

In accordance with Eq. (72), the amount of α-phase formed depends not only on the thermodynamic force Δφ, but also on the magnitude of the stresses in the γ-phase.

## 6. Scheme of the nondiffusion transformation of austenite based on the constructed model.

Before discussing the equations obtained, we introduce some more useful relations characterizing the γ → α transformation. With the γ → α transformation, the effective atomic volume of the iron lattice changes in the sample under consideration, characterized by ΔVγ → α and the relative volume change ΔVγαVγ. According to the data of [3]:

ΔVγα=0.2681.62104Т,sm3/molE73

We will assume that with the formation of the α-phase, the relative change in volume is determined by the additional deformation: ΔVγαVγ= 3Δεα, and the compressive stress arising in the α-phase has the value.

σα=E12μΔεα.E74

When the alloy sample is cooled by ΔT, a deformation occurs in its surface layer: εγ ∼ αΔТ and the tensile stress σγ corresponding to this deformation:

σγ=E12μαΔТ.E75

Comparing the values of thermodynamic forces among themselves, it is possible to classify the types of nondiffusion transformation according to the kinetic criterion. As shown in [1], p. 208, for small deviations of the system from equilibrium, the growth of crystals is more likely, controlled by self-diffusion, at large–cooperative growth. The same phase transition in a single-component system under different external conditions can take place with an independent (or slightly dependent) temperature growth rate (martensitic kinetics) and with a rate that exponentially depends on the temperature at an activation energy close to the activation energy of self-diffusion (normal kinetics). The parameter characterizing the deviation of the system from equilibrium is the supercooling of the alloy ΔT = Ac3T, where Ac3 is the temperature of the end α → γ of the conversion upon heating, and T is the transformation temperature. The transformation scheme for the constructed model is shown in Figure 5.

Ac1 is the temperature of the beginning of α → γ transformation when the alloy is heated and Mni is the temperature of the onset of the formation of isothermal martensite upon supercooling of the alloy. Mn is the temperature of the onset of athermal martensite formation upon supercooling of the alloy. Mk is the temperature of the end of martensite formation upon supercooling of the alloy.

Thus, for small

ΔT:L11φ>L12σα>L12σγ,E76

then the growth of α-phase crystals is determined by self-diffusion by the normal mechanism. However, as follows from Eq. (72), in this case too, the contribution of deformations (and stresses) to the conversion kinetics is very significant. In order that the condition (76) is satisfied, it is necessary that the stress level in the γ- and α-phases be small; for the α-phase, this is possible only in the case of relaxation of internal stresses in the alloy at high temperature by the mechanism of recrystallization.

With increasing supercooling of the alloy, the thermodynamic stimulus and the rate of normal transformation increase.

For a larger

ΔТ:L11φL12σα>L12σγE77

The existing thermal stresses in the γ phase (75) contribute to the formation of the α-phase by the shear mechanism, and the stresses arising in the α-phase compensate thermal stresses in the γ-phase. With a certain amount of α-phase, the stress equals σα = σγ arises and the further formation of the α-phase occurs according to the normal mechanism with the relaxation of the arising stresses by recrystallization. Consequently, the condition (77) corresponds to the transformation of the γ-phase by a mixed mechanism, and also to the formation of a ferrite side plates (Widmanstätten), followed by the release of the α-phase by the normal mechanism [1].

With a certain supercooling of ΔTi, stress compensation occurs only when the γ-phase is completely transformed into ferrite by a shearing mechanism. In this case:

L11∆φL12σα,σα=σγ

The temperature corresponding to this supercooling is the starting point for the formation of the isothermal martensite Mni (Figure 5). Below the point Mni, the formation of the α-phase occurs by a shearing mechanism. However, the normal component of the process still has a significant value, affecting the morphology of the resulting precipitates. When supercooling a greater ΔTi, L11∆φ < L12σα, σα < σγ

At temperatures below Mni, isothermal martensite or acicular ferrite is formed with a “reticular” or acicular morphology of precipitates. Finally, for large ΔT (below Мna):

L11φL12σγ,E78

Inequality (78) determines the condition for the formation of “athermal” martensite, when the normal component does not affect the formation of the shear structure. The main effect on the rate of the γ → α transformation, in accordance with expression (71), is due to thermal stresses in the γ phase. Thus, the constructed model of the nondiffusion austenite transformations allows us to consider the normal and martensitic transformations, as limiting cases.

## 7. Conclusions

Based on the possibility of dynamic equilibrium, expressions are found for calculating the cross-kinetic coefficients of a thermodynamic system consisting of two and three components. The values ​​of the thermodynamic force for diffusion of carbon, kinetic coefficients and flows of a thermodynamic system describing the kinetics of carbide precipitation during the tempering of chromium steel are calculated. It has been established that the values ​​of iron and chromium fluxes increase substantially due to the cross ratios and the significant magnitude of the thermodynamic force (−ΔμC).

Analysis of the eutectoid transformation of austenite using the relations of nonequilibrium thermodynamics allowed us to generalize the equations of motion of the system obtained earlier by the authors of [20] and to find more accurate theoretical expressions for the perlite growth rate and its between interplate distance on the magnitude of the supercooling of steel. According to the constructed model, the perlite growth rate in the direction of the X axis has a maximum value at supercooling ΔТ = 140.0°С. The perlite growth rate calculated according to Zener’s formula has a theoretical maximum value at overcooling ΔТ = 96.0°С. Consequently, the theoretical expressions (31) and (32) make it possible to describe with greater accuracy the maximum and the course of the experimental curve for the perlite formation for high-purity eutectoid steel.

The application of nonequilibrium thermodynamics to the analysis of the nondiffusion transformation of austenite made it possible to obtain a system of equations for the thermodynamic system and to generalize the results obtained earlier by B.Ya. Lyubov the equations for a normal transformation. The theoretical expression for the growth rate of the α-phase, obtained in this paper, takes into account the influence of stresses on the process of austenite transformation. It is shown that the rate of growth of α-phase particles at a constant temperature very rapidly (exponentially) decreases in time, determining the incompleteness of the transformation. According to the constructed model, a scheme of nondiffusion austenite transformations was developed, including normal and martensitic transformations, as limiting cases.

Thus, the use of the principles of nonequilibrium thermodynamics makes it possible to obtain completely new results in the analysis of phase transformations in iron-carbon alloys.

chapter PDF

## More

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## How to cite and reference

### Cite this chapter Copy to clipboard

Bobyr Sergiy Volodimyrovych (May 30th 2019). Using the Principles of Nonequilibrium Thermodynamics for the Analysis of Phase Transformations in Iron-Carbon Alloys [Online First], IntechOpen, DOI: 10.5772/intechopen.83657. Available from: