Advances in Materials Science and Engineering
Volume 2019 (2019), Article ID 5985617, 12 pages
https://doi.org/10.1155/2019/5985617

Research Article

Modeling of the Growth Kinetics of Boride Layers in Powder-Pack Borided ASTM A36 Steel Based on Two Different Approaches

M. Ortiz-Domínguez,¹ O. A. Gómez-Vargas,² G. Ares de Parga,³ G. Torres-Santiago,² R. Velázquez-Mancilla,² V. A. Castellanos-Escamilla,² J. Mendoza-Camargo,² and R. Trujillo-Sánchez²

¹Universidad Autónoma del Estado de Hidalgo, Escuela Superior de Ciudad Sahagún-Ingeniería Mecánica, Carretera Cd. Sahagún-Otumba s/n, CP. 43990, Hidalgo, Mexico
²Instituto Tecnológico de Tlalnepantla-ITTLA. Av., Instituto Tecnológico, S/N. Col. La Comunidad, Tlalnepantla de Baz, CP. 54070, Estado de Mexico, Mexico
³Instituto Politécnico Nacional-ESFM, U.P. Adolfo López Mateos, Zacatenco, CP. 07738. Cd, De Mexico, Mexico

Correspondence should be addressed to M. Ortiz-Domínguez; martin_ortiz@uaeh.edu.mx

Received 11 October 2018; Accepted 7 July 2019; Published 7 October 2019

Academic Editor: Giorgio Pia

Copyright © 2019 M. Ortiz-Domínguez et al. 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.

An indispensable tool to choose the suitable process parameters for obtaining boride layer of an adequate thickness is the modeling of the boriding kinetics. In this work, two mathematical approaches were used in order to determine the value of activation energy in the Fe₂B layers on ASTM A36 steel during the iron powder-pack boriding in the temperature range of 1123–1273 K for treatment times between 2 and 8 h. The first approach was based on the mass balance equation at the interface (Fe₂B/substrate) and the solution of Fick’s second law under steady state (without time dependent). The second approach was based on the same mathematical principles as the first approach for one-dimensional analysis under non-steady-state condition. The measurements of the thickness (Fe₂B), for different temperatures of boriding, were used for calculations. As a result, the boron activation energy for the ASTM A36 steel was estimated as 161 kJ·mol⁻¹. This value of energy was compared between both models and with other literature data. The Fe₂B layers grown on ASTM A36 steel were characterized by use of the following experimental techniques: X-ray diffraction (XRD), scanning electron microscopy (SEM), and energy dispersive X-ray Spectroscopy (EDS). Finally, the experimental value of Fe₂B layer’s thickness obtained at 1123 K with an exposure time of 2.5 h was compared with the predicted thicknesses by using these two approaches. A good concordance was achieved between the experimental data and the simulated results.

1. Introduction

Nowadays and due to the increasing technological development, it is necessary to have metallic materials with specific features that must be maintained in critical service conditions: for example, the metal dies are used in the different hot and cold working metallurgical processes, which given the working conditions require high toughness and high surface hardness. The thermochemical treatments applied to steel are those in which the composition of the surface of the workpiece is altered by the addition of carbon, nitrogen, sulphur, boron, aluminium, zinc, chromium, or other elements. The most common treatments in the industry are carburization, nitriding, carbonitriding, and boriding. Boriding is a thermochemical treatment controlled by the diffusion of boron atoms, which modifies the properties of the material generating hard surfaces. Likewise, the boriding process has a positive effect on the tribological applications: abrasive, adhesive, fatigue, and corrosion wear in acid and alkaline media. The process involves heating a ferrous or nonferrous alloy for a temperature range of 700°C to 1000°C with a treatment time of 1 to 12 hours for powder-pack methods. When boron diffuses in a substrate, a monolayer (Fe₂B) or a double layer (FeB–Fe₂B) can be formed, depending on the chemical composition of the substrate and the chemical potential of boron. These phases consist of orthorhombic and tetragonal lattices (body-centred), respectively. One basic advantage of boride layers is that they can reach high hardness values (between 1800 and 2000 HV), kept at high temperatures [1–4]. Abrasion and adhesion wear are forms of wear by contact between a particle and solid material, being the characteristic result of almost all types of mechanical stress. Borided steels are resistant to abrasion because of their extreme hardness on the surface; this characteristic makes them suitable to be applied in pneumatic conveying systems; dies for stamping; components of plastic processing machines, such as extrusion screws; bearings for oil extraction pumps; ball valves; plungers for use in manufacturing glass; and components in textile machinery. In practice, there are many techniques of surface modification by boriding, such as powder-pack boriding [5], paste boriding [6], gaseous boriding [7], plasma boriding [8], plasma paste boriding [9], and laser boriding [10]. However, the most frequently used method in industry is the powder-pack boriding, which demands a low investment cost of equipment and an easy handling. From a kinetic point of view, several approaches [3, 5–7, 11–33] were developed in the objective of optimizing the thicknesses of borided layers in order to meet the functional requirements during industrial use of borided steels. Some of these models that estimate the thickness of the monolayer (Fe₂B) or a double layer (FeB–Fe₂B) are based on the solution of Fick’s second law without time dependent ( ⟶ steady state) [3, 6, 7, 16–18, 20–22, 24–26, 30] and some others on the solution of Fick’s second law with time dependent ( ⟶ non-steady state) [5, 24, 27, 28, 31, 33].

The ASTM A36 steel has a good machinability with an acceptable wear resistance. ASTM A36 steel has the following applications:

(i)It is used in bolted, riveted, or welded construction of bridges, buildings, and oil rigs(ii)It is used in forming tanks, bins, bearing plates, fixtures, rings, templates, jigs, sprockets, cams, gears, base plates, forgings, ornamental works, stakes, brackets, automotive and agricultural equipment, frames, and machinery parts(iii)It is used for various parts obtained by flame cutting such as in parking garages, walkways, boat landing ramps, and trenches

In the present work, two different models were suggested for simulating the growth kinetics of monolayer (Fe₂B) on ASTM A36 steel in the range of boriding temperature 1123–1273 K. The parabolic growth constants of Fe₂B were determined. The boron diffusion coefficient in the boride layers was estimated from these two approaches based on the conditions of the boriding process in the Fe-B system. X-ray diffraction (XRD), scanning electron microscopy (SEM), and energy dispersive X-ray spectroscopy (EDS) were conducted on material borided to characterize the presence of the Fe₂B layer and the distribution of heavy elements in the surface of the ASTM A36 steel. Finally, the experimental value of Fe₂B layer’s thickness obtained at 1123 K with an exposure time of 2.5 h was compared with the predicted result from these two approaches.

2. Mathematical Approaches

One basic parameter that characterizes the Fe₂B layers is the thickness, since the properties of the coating depend on it, such as resistance to wear, fatigue, hardness, and dynamic loads, as well as to a large extent determining the grip with the substrate. Having an expression that allows estimating the layer thickness during the boriding process facilitates the appropriate selection of the technological parameters in order to guarantee the desired properties. The layer thickness exhibits a time dependent such that

2.1. Derivation of the Parabolic Growth Law

In diffusion in solid, parabolic kinetics occurs when the mass gain on a sample is proportional to the square root of time. In general, parabolic kinetics indicates that diffusion of reactants (such as boron) through a growing layer is rate-determining. If the diffusion of B atoms is rate-determining, the layer rate is proportional to the flux through the substrate:

El flux, , can be written as

where is the boron concentration profile in mol/m³ and is the velocity of Fe₂B layer in m/s, and giving units of mol/m²·s. The velocity of a particle is proportional to the force, F, on the particle:

where is the mobility of the boron atoms. Writing the chemical potential as , this force is written as

for a Fe₂B layer with thickness x. If combining equations (4) and (5), then equation (3) can be written as

From the relationship,

where is the Boltzmann’s constant, we can write

In an ideal system, the concentration, , is equivalent to activity, . Substituting equations (8) into (6), we get

As shown in equation (2),

so that a combination of equations (2) and (9) gives

If we assume that the potential is fixed at each boundary of the Fe₂B layer, we can replace in equation (11) with the slope (). We then introduce the parabolic growth constant , and set

Combining equations (11) and (12) gives

Equation (13) can be rewritten as

Upon integration of equation (14),

We arrive at the parabolic growth law (with ):

where represents the Fe₂B layer thickness.

2.2. First Approach: Steady-State Diffusion Model

The first approach is based on the diffusion model proposed by Flores-Rentería et al. [3], where a mathematical model has been applied based on the mass balance equation at the (Fe₂B/substrate) interface by assuming a linear boron concentration profile through the Fe₂B layer ( ⟶ steady state). This approach was used to simulate the kinetics of formation of Fe₂B layer on ASTM A36 steel with the presence of boride incubation time. Steady state means that there will not be any change in the composition profile with time. and denote the upper and lower boron concentrations in the Fe₂B phase [3, 34, 35]. represents the boron solubility in the matrix and can be neglected [3, 34, 35]. The assumptions made during the mathematical formulation of the diffusion model are given in the reference work [3].

The mass balance equation [3], describing the evolution of displacement of growing interface with respect to the time, is given by

The linear boron concentration profile [3], through the Fe₂B layer, is given by the solution of Fick’s second law without time dependent ( ⟶ steady state ( (is called the Laplace operator or Laplacian in one dimension)) and is rewritten as

By substituting equations (18) and (16) into equation (17),

2.3. Second Approach: Non-Steady-State Diffusion Model in One Dimension

The second approach [5] was applied to analyze the kinetics of formation of the monolayer (Fe₂B) generated at the surface of pack-borided AISI 1045 steel. In the present work, this mathematical model was adopted for studying the boriding kinetics of ASTM A36 steel. Likewise, the mass balance equation [5], describing the evolution of displacement of growing interface with respect to the time, is given by

The linear boron concentration profile [5], through the Fe₂B layer, is given by the solution of Fick’s second law with time dependent ( ⟶ non-steady state) and is deduced as follows:

By substituting equations (21) into (20), the following equation is obtained:

Substituting the expression of the parabolic growth law obtained from equation (16) () into equation (22), we have

The diffusion coefficient () can be estimated numerically from equation (23) by the Newton–Raphson method. An illustrative representation of the parabolic growth law of Fe₂B layer thickness () is represented in Figure 1. In addition, is the effective growth time of the Fe₂B layer and t is the boriding time [3, 5].

Figure 1: A schematic diagram of the parabolic growth law of Fe₂B layer thickness against boriding time [3, 5].

3. Materials and Methods

3.1. Powder-Pack Boriding Process

ASTM A36 steel was used for investigation. It had a nominal chemical composition of 0.25–0.29% C, 0.20–0.28% Si, 0.85–1.35% Mn, 0.15–0.20% Cu, 0.035–0.040% P, and 0.050% S. The steel samples were sectioned into small cubes with the following dimensions: 10 mm × 10 mm × 10 mm. Prior to the boriding process, the steel samples were grinded with SiC abrasive paper up to grit 2500 and cleaned using a multistage ultrasonic bath with n-heptane and ethanol for 20 min. The mean hardness of the substrate was 170 HV. The ASTM A36 steel samples were immersed in a closed cylindrical case made of AISI 316L steel as shown in Figure 2, using Ekabor 2 as a boron-rich mixture.

Figure 2: Schematic view of the stainless steel AISI 316L container for the pack-powder boriding treatment (1: lid; 2: powder boriding medium (B₄C + KBF₄ + SiC); 3: sample; 4: container) (millimeter scale) [29, 32].

The thermochemical process was carried out in a conventional furnace model Nabertherm N 250/85 HA (this type of furnace is gas-tight; it is equipped with direct heating depending on the temperature. It is excellent for maintaining an atmosphere defined by an inert gas), maintaining a pure argon atmosphere, to eliminate the oxidation of the boron released in the chemical reaction of the boriding medium [1]. The thermochemical treatment was carried out at boriding temperatures of 1123, 1173, 1223, and 1273 K for a variable time (2, 4, 6, and 8 h). The treatment temperatures were selected according to the Fe-B phase diagram.

3.2. Microscopical Observations of Boride Layers

The hardened samples were sectioned and prepared metallographically (the samples were polished using a diamond suspension with a particle size of 6 μm, finishing with a particle size of 3 μm), using a GX51 Olympus equipment. Likewise, the borided samples were analyzed through scanning electron microscope. The equipment used was the Quanta 3D FEG-FEI JSM7800-JOEL. Figure 3 shows the cross sections of Fe₂B layers formed on the surfaces of ASTM A36 steel at different exposure times (2, 4, 6, and 8 h) and for 1173 K of boriding temperature. The mechanical properties of borided alloys depend on the composition and structure of the boride layer. The images obtained from the scanning electron microscope (Figure 3) present a sawtooth morphology; this characteristic is typical of ARMCO pure iron, and low and medium carbon steels [24, 36].

(a)

(b)

(c)

(d)

Figure 3: SEM micrographs of the cross sections of ASTM A36 steel samples borided at 1123 K during different exposure times: (a) 2 h, (b) 4 h, (c) 6 h, and (d) 8 h.

When the alloying elements and/or the carbon content of the steel increases, the layer thickness tends to favor the formation of iron borides with flat growth fronts. The alloying elements have obvious effects on the formation of the layer thickness, restricted diffusion of the boron atoms, thus forming a diffusion barrier. Because of the diffusion of boron atoms, there is a segregation of the alloying elements from the surface to the (Fe₂B/substrate) interface. Some alloying elements tend to form compounds with the boron atoms and others cannot interact with them. Boron has some weird and wonderful chemistry. The alloying elements cannot form compounds and tend to concentrate at the tips of boride columns, decreasing the layer thickness [36].

Figure 3 shows that the boride layer thickness increases with respect to the boriding time for a predetermined temperature. For simulating the growth kinetics of the Fe₂B layer grown on ASTM A36 steel, an average measurement of the boride layer thickness was made, where the longest tips of the sawtooth morphology were taken into account (see Figure 4), and the software used was MSQ Materials Analysis. Fifty measurements were collected from the boride surface to the longest tips of boride columns of the borided ASTM A36 steel, as plotted in Figure 4 [3, 5]. The identification of phases formed on the surface of the borided sample was conducted through the X-ray diffraction technique (XRD). The equipment used for the study was an INEL EQUINOX 2000 X-ray diffractometer, using CoKα radiation of 0.179 nm wavelength, operated at 30 mA and 20 kV.

Figure 4: A schematic diagram illustrating the procedure for estimation of boride layer thickness in ASTM A36 steel [3, 5].

Likewise, Match version 3.3 was the software for phase identification. In addition, the elemental distribution of the transition elements within the cross section of boride layer was determined by using scanning electron microscopy-energy dispersive X-ray spectroscopy (SEM-EDS). The equipment used for the study was a Quanta™ 3D FEG scanning electron microscope, with an accelerating voltage of 200 V–30 kV and a magnification of 30 X–1280 kX in “quad” mode.

4. Results and Discussion

4.1. SEM Observations and EDS Analysis

The metallography of coating/substrate formed in ASTM A36 borided steel at different exposure times (2, 4, 6, and 8 h) and for 1123 K of boriding temperature is shown in Figure 5.

(a)

(b)

Figure 5: SEM micrograph of the cross sections of the borided ASTM A36 steel micrograph image of microstructure of the ASTM A36 boride layer obtained at 1123 K with exposure time of 8 h (a) and (b) EDS analysis on cross section of the borided sample.

The EDS analysis obtained by SEM is shown in Figures 5(a) and 5(b). The results show in Figure 5(a) that the manganese (Mn) negatively affects the boride layer thickness and morphology. Likewise, as can be seen in Figure 5(b), carbon (C) and silicon (Si) do not dissolve and diffusive through the boride layer. During the boriding process, carbon is transmitted from the boride surface to the matrix and forms, together with boron, borocementite Fe₃(B, C) (or, more appropriately, Fe₃(B0.67C0.33)) [21, 32]. In addition to carbon, silicon is also not soluble in the boride layer. This element is expelled from the surface by boron atoms to the growth interface (Fe₂B/substrate), forming iron-silicoborates (FeSi0.4B0.6 and Fe₅SiB₂) [37].

4.2. X-Ray Diffraction Analysis

Figure 6 shows the diffractogram recorded on the surface of the borided ASTM A36 steel at a temperature of 1123 K for a treatment time of 8 h. The X-ray diffraction patters (see Figure 6) show the intensity of X-rays scattered at different angles by a borided sample, where the presence of the Fe₂B phase is confirmed.

Figure 6: XRD patterns obtained at the surface of the borided ASTM A36. The boriding temperature of 1123 K for 8 h of treatment.

4.3. First Approach: Estimation of Boron Activation Energy with Steady-State Model

To study the growth kinetic of Fe₂B on ASTM A36 steel, a simple diffusion model inspired from the mass balance equation at the (Fe₂B/substrate) interface and the solution of Fick’s second law under steady-state condition for one-dimensional analysis was applied [3, 6, 7, 16–18, 20–22, 24–26, 30]. The Fe₂B layer thickness () obeys the parabolic growth law given by equation (16), with the presence of boride incubation time ( after transferring the sample to the furnace) associated with the formation of the Fe₂B layer. Figure 7 displays the time dependent of the squared value of Fe₂B layer thickness for different temperatures, the slopes of each of the straight lines supply the values of the parabolic growth constants ().

Figure 7: Averaged squares of the boride layers () vs. boriding time () at increasing temperatures.

Table 1 summarizes the experimental values of parabolic growth constants at the (Fe₂B/substrate) interface in the temperature range of 1123–1273 K with the associated boride incubation time. Also, the plots presented in Figure 7 demonstrate a diffusion-controlled process (Brackman et al. [21]).

Table 1: Experimental values of parabolic growth constants at the (Fe₂B/substrate) interface in the range of 1123–1273 K and estimated values of boron diffusion coefficients in the Fe₂B layers by using the first approach.

Temperature (K)
Type of layer
Growth constants (m²s⁻¹)
equation (19) (m²s⁻¹)

1123
Fe₂B
2.810 × 10⁻¹
4.28 × 10⁻¹¹

1173
6.120 × 10⁻¹
9.32 × 10⁻¹¹

1223
12.80 × 10⁻¹
1.95 × 10⁻¹⁰

1273
20.89 × 10⁻¹
3.18 × 10⁻¹⁰

The value of activation energy (is the energy required for a reaction to proceed) for boron diffusion in Fe₂B ( = 161.00061 kJ·mol⁻¹) and preexponential factor ( = 1.361592 × 10⁻³ m²/s) can be calculated after a linear fitting of these data according to Arrhenius relationship from the slope and interception of the straight line, respectively, shown in coordinate system (see equation (24)): vs , it is displayed in Figure 8:

with R = 8.314 Jmol⁻¹K⁻¹ and T the absolute temperature in Kelvin.

Figure 8: Arrhenius relationship for the boron diffusion coefficient () from the first approach in the Fe₂B layer.

4.4. Second Approach: Estimation of Boron Activation Energy with Non-Steady-State Diffusion Model

The values of the growth constants () at each temperature reported in Table 1 are also used in the second approach and in Table 2, the values of the boron diffusion coefficients () for Fe₂B are gathered, which were calculated by the Newton–Raphson method using equation (23).

Table 2: Estimated values of boron diffusion coefficients in the Fe₂B layer using the second approach (non-steady-state diffusion model).

Temperature (K)
Type of layer
equation (23) (m²s⁻¹)

1123
Fe₂B
4.28076 × 10⁻¹¹

1173
9.3056 × 10⁻¹¹

1223
1.9467 × 10⁻¹⁰

1273
3.1822 × 10⁻¹⁰

In Figure 9, the vs is plotted, the values of activation energy and pre exponential factor for boron diffusion in Fe₂B from the second approach were  = 160. 9922 kJ·mol⁻¹ and  = 1.35914 × 10⁻³ m²/s, respectively. Equation (25) was deducted from a fitting of data according to the Arrhenius relationship with the coefficient of determination close to unity.

with R = 8.314 Jmol⁻¹K⁻¹ and T the absolute temperature in Kelvin.

Figure 9: Arrhenius relationship for the boron diffusion coefficient () from the second approach in the Fe₂B layer (non-steady-state diffusion model).

4.5. The Two Diffusion Models

This section describes the differences between the two diffusion models that have been used to compute the growth kinetics of boride layers. It is noticed that the estimated values of boron activation energy for ASTM A36 steel from the first approach (see equation (24)) and the second approach (see equation (25)) is approximately the same value for both diffusion models. Likewise, in the estimated values of pre exponential factor from the first approach (D₀ = 1.361592 × 10⁻³ m²/s) and the second approach (D₀ = 1.35914 × 10⁻³ m²/s), there is a small variation. To find out how this similarity is possible in the diffusion coefficients obtained by two different models, we first focus our attention on equation (23). The error function is a monotonically increasing function of . Its Maclaurin series (for small ) is given by [38, 39]:

According to the numerical value of the , equation (26) can be rewritten as

Similarly for the real exponential function : ℝ ⟶ ℝ can be characterized in a variety of equivalent ways. Most commonly, it is defined by the following power series [38, 39]:

Thus, equation (28) can be modified in the following manner:

By substituting Equations (27) and (29) into equation (23), we have

where

The result obtained by equation (31) is the same as that obtained in equation (19) estimated by a steady-state diffusion model. The result from equation (19) would appear to imply that the non-steady-state diffusion model is superior to the steady-state diffusion model and so should always be used. However, in many interesting cases, the models are equivalent as in this case.

Table 3 presents the values of activation energies for boron diffusion in some different alloys and Armco iron along with the estimated value of boron activation energy in ASTM A36 steel. From Table 3, it is seen that the estimated value of activation energy in ASTM A36 steel is consistent with the literature data [3, 5–7, 27, 31, 33, 40–43].

Table 3: Reported values of activation energies for boron diffusion in ASTM A36 steel with other borided Armco iron and steels.

Material
Boriding method
Boron activation energy (kJ·mol⁻¹)
Temperature range (°C)
Method of calculation
References

Armco iron
Gaseous
73.08 (FeB)
120.65 (Fe₂B)
800–1000
Diffusion model
[7]

Armco iron
Powder
157.5
850–1000
Diffusion model
[27]

AISI M2 steel
Paste
257.5 (FeB)
210.0 (Fe₂B)
920–1000
Diffusion model
[6]

AISI 1018 steel
Electrochemical
172.75 ± 8.6
850–1000
Parabolic growth law
[40]

AISI 1026 steel
Powder
178.4
850–1000
Diffusion model
[3]

AISI 1045 steel
Powder
180.0
850–1000
Diffusion model
[5]

AISI 8620
Plasma paste boriding
124.7–138.5
700–800
Parabolic growth law
[41]

AISI 4340 steel
Salt bath
234.0
800–1000
Parabolic growth law
[42]

AISI D2 steel
Salt bath
170.0
800–1000
Parabolic growth law
[42]

AISI 1018 steel
Powder
159.3 (Fe₂B)
850–1000
Diffusion model
[31]

AISI D2 steel
Powder
201.5 (Fe₂B)
850–1000
Diffusion model
[33]

AISI P20 steel
Powder
200 (FeB + Fe₂B)
800–950
Parabolic growth law
[43]

ASTM A36 steel
Powder
161.00 (Fe₂B)
850–1000
Diffusion models
This work

4.6. Experimental Validation of the Kinetic Models

The suggested models were validated by comparing the experimental value of Fe₂B layers’ thickness with the theoretical result. By substituting equation (19) into (16), we have

Figure 10 shows the SEM micrographs of the cross sections of borided ASTM A36 steels at 1123 K for 2.5 h. Table 4, presents the values of the theoretical result of Fe₂B layers’ thickness with the experimental data. A good concordance was the observed value between the experimental result and the simulated value for the given boriding condition.

Figure 10: Scanning electron microscopy (SEM) image of the microstructure of boride layer (Fe₂B) formed on the ASTM A36 steel during powder pack at 1123 K with an exposure time of 2.5 h.

Table 4: Estimated value of ASTM A36 boride layer thickness () obtained at the 1123 K of boriding temperature with an exposure time of 2.5 h.

Temperature (K)
Type of layer
Boride layer thicknesses (μm) estimated by equation (32)
Experimental boride layer thickness (μm) for an exposure time of 2.5 h

1123
Fe₂B
Exposure time of 2.5 h
 = 46.081 μm
 = 42.066 ± 6.9947 μm

5. Conclusions

In this work, the ASTM A36 steel was pack-borided in the temperature range of 1123–1273 K for a variable exposure time ranging between 2 and 8 h. The kinetics data on treated ASTM A36 steel by the powder-pack boriding were used to estimate the value of activation energy for boron diffusion in the boride layer (Fe₂B) by means of two different mathematical approaches. In the first approach, the mass balance equation was formulated by assuming a linear boron concentration profile in the boride layer () for an upper boron content in Fe₂B of 60 × 10³ mol·m⁻³. The second approach was based on the same mathematical principles as the first approach for one-dimensional analysis under non-steady-state condition by assuming the solution of Fick’s second law of diffusion when the bulk concentration is greater than the surface concentration (), was successfully applied to the boriding kinetics of ASTM A36 steel by considering the effect of boride incubation time. As a main result, the estimated value of activation energy for boron diffusion in the Fe₂B layer by the two approaches was estimated as 161 kJ·mol⁻¹. In a future prospect, these diffusion approaches are in general not identical, both are equivalent models, and this fact can be used as a tool to select optimum values of layers’ thicknesses for practical utilization of any borided steels to produce boride layers with sufficient thicknesses that meet the requirements during service life.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The work described in this paper was supported by a grant of National Council of Science and Technology (CONACyT) and PRODEP México. Also, the authors want to thank Dr. Jorge Zuno Silva, who is a research professor and principal at Escuela Superior de Ciudad Sahagún-UAEH. Likewise, the authors would like to thank Dr. José Solís Romero for his helpful advice on various technical issues examined in this manuscript and Dr. Adrian Leyland for his advice and comments.

References

M. Graf von Matuschka, Boronizing, Carl Hanser Verlag, Munich, Germany, 1st edition, 1980.

J. R. Davis, Surface Hardening of Steels Understanding the Basics, ASM International, Novelty, OH, USA, 2002.

M. A. Flores-Rentería, M. Ortiz-Domínguez, M. Keddam et al., “A simple kinetic model for the growth of Fe₂B layers on AISI 1026 steel during the powder-pack boriding,” High Temperature Materials and Processes, vol. 34, no. 1, pp. 1–11, 2015.

R. Kara, F. Çolak, and Y. Kayali, “Investigation of wear and adhesion behaviors of borided steels,” Transactions of the Indian Institute of Metals, vol. 69, no. 6, pp. 1169–1177, 2016.

J. Zuno-Silva, M. Ortiz-Domínguez, M. Keddam et al., “Boriding kinetics of Fe₂B layers formed on AISI 1045 steel,” Journal of Mining and Metallurgy, Section B: Metallurgy, vol. 50, no. 2, pp. 101–107, 2014.

I. Campos, G. Ramírez, U. Figueroa, J. Martínez, and O. Morales, “Evaluation of boron mobility on the phases FeB, Fe₂B and diffusion zone in AISI 1045 and M2 steels,” Applied Surface Science, vol. 253, no. 7, pp. 3469–3475, 2007.

M. Keddam, M. Kulka, N. Makuch, A. Pertek, and L. Małdziński, “A kinetic model for estimating the boron activation energies in the FeB and Fe₂B layers during the gas-boriding of Armco iron: effect of boride incubation times,” Applied Surface Science, vol. 298, pp. 155–163, 2014.

E. R. Cabeo, G. Laudien, S. Biemer, K. T. Rie, and S. Hoppe, “Plasma-assisted boriding of industrial components in a pulsed D. C. glow discharge,” Surface and Coatings Technology, vol. 116–119, pp. 229–233, 1999.

M. Keddam, R. Chegroune, M. Kulka et al., “Characterization, tribological and mechanical properties of plasma paste borided AISI 316 steel,” Transactions of the Indian Institute of Metals, vol. 71, no. 1, pp. 79–90, 2018.

M. Kulka, N. Makuch, and A. Pertek, “Microstructure and properties of laser-borided 41Cr4 steel,” Optics & Laser Technology, vol. 45, pp. 308–318, 2013.

B. Mebarek, M. Keddam, and M. Aboshighiba, “LS-SVM approach for modeling the growthe kinetics of FeB and Fe₂B layers formed on Armco iron,” Ingénierie des Systèmes d’Information, vol. 23, no. 5, pp. 29–41, 2018.

I. Campos, M. Islas, E. González, P. Ponce, and G. Ramírez, “Use of fuzzy logic for modeling the growth of Fe₂B boride layers during boronizing,” Surface and Coatings Technology, vol. 201, no. 6, pp. 2717–2723, 2006.

K. Rayane and O. Allaoui, “Application of Artificial Neural Network for prediction of boride layer depth obtained on XC38 steel in molten salts,” Defect and Diffusion Forum, vol. 365, pp. 194–199, 2015.

I. Campos, M. Islas, G. Ramírez, C. VillaVelázquez, and C. Mota, “Growth kinetics of borided layers: artificial neural network and least square approaches,” Applied Surface Science, vol. 253, no. 14, pp. 6226–6231, 2007.

R. D. Ramdan, T. Takaki, K. Yashiro, and Y. Tomita, “The effects of structure orientation on the growth of Fe,” Materials Transactions, vol. 51, no. 1, pp. 62–67, 2010.

M. A. Doñu, N. López-Perrusquia, D. Sánchez-Huerta et al., “Growth kinetics of boride coatings formed at the surface AISI M2 during dehydrated paste pack boriding,” Thin Solid Films, vol. 596, pp. 147–154, 2015.

M. Keddam and M. Kulka, “A kinetic model for the boriding kinetics of AISI D2 steel during the diffusion annealing process,” Protection of Metals and Physical Chemistry of Surfaces, vol. 54, no. 2, pp. 282–290, 2018.

M. Keddam and M. Kulka, “Simulation of the growth kinetics of FeB and Fe₂B layers on AISI D2 steel by the integral method,” Physics of Metals and Metallography, vol. 119, no. 9, pp. 860–869, 2018.

V. I. Dybkov, L. V. Goncharuk, V. G. Khoruzha, K. A. Meleshevich, A. V. Samelyuk, and V. R. Sidorko, “Diffusional growth kinetics of boride layers on iron-chromium alloys,” Solid State Phenomena, vol. 138, pp. 181–188, 2008.

M. Kulka, N. Makuch, A. Pertek, and L. Maldzinski, “Simulation of the growth kinetics of boride layers formed on Fe during gas boriding in H₂-BCl₃ atmosphere,” Journal of Solid State Chemistry, vol. 199, pp. 196–203, 2013.

C. M. Brakman, A. W. J. Gommers, and E. J. Mittemeijer, “Boriding Fe and Fe-C, Fe-Cr, and Fe-Ni alloys; Boride-layer growth kinetics,” Journal of Materials Research, vol. 4, no. 6, pp. 1354–1370, 1989.

L. G. Yu, X. J. Chen, K. A. Khor, and G. Sundararajan, “FeB/Fe₂B phase transformation during SPS pack-boriding: boride layer growth kinetics,” Acta Materialia, vol. 53, no. 8, pp. 2361–2368, 2005.

I. Campos-Silva, M. Ortiz-Domínguez, N. López-Perrusquia, A. Meneses-Amador, R. Escobar-Galindo, and J. Martínez–Trinidad, “Characterization of AISI 4140 borided steels,” Applied Surface Science, vol. 256, no. 8, pp. 2372–2379, 2010.

M. Ortiz-Domínguez, “Contribución de la Modelación Matemática en el Tratamiento Termoquímico de Borurización,” SEPI-ESIME from the Instituto Politécnico Nacional, Mexico, 2013, Ph.D. thesis.

I. Campos-Silva, D. Bravo-Bárcenas, A. Meneses-Amador et al., “Growth kinetics and mechanical properties of boride layers formed at surface of the ASTM F–75 biomedical alloy,” Surface and Coatings Technology, vol. 237, pp. 402–414, 2013.

I. Campos-Silva, M. Ortiz-Domínguez, O. Bravo-Bárcenas et al., “Formation and kinetics of FeB/Fe₂B layers and diffusion zone at the surface of AISI 316 borided steels,” Surface and Coatings Technology, vol. 205, pp. 403–412, 2010.

M. Elias-Espinosa, M. Ortiz-Domínguez, M. Keddam et al., “Growth kinetics of the Fe₂B layers and adhesion on Armco iron substrate,” Journal of Materials Engineering and Performance, vol. 23, pp. 2943–2952, 2014.

M. Ortiz-Domínguez, M. Keddam, M. Elias-Espinosa et al., “Investigation of boriding kinetics of AISI D2 steel,” Surface Engineering, vol. 30, pp. 490–497, 2014.

M. Ortiz-Domínguez, M. A. Flores-Rentería, M. Keddam et al., “Simulation of growth kinetics of Fe₂B layers formed on gray cast iron during the powder-pack boriding,” Materiali in Tehnologije/Materials and Technology, vol. 48, no. 6, pp. 905–916, 2014.

I. Campos-Silva, M. Ortiz-Domínguez, C. Tapia-Quintero, G. Rodriguez-Castro, M. Y. Jimenez-Reyes, and E. Chavez-Gutierrez, “Kinetics and boron diffusion in the FeB/Fe₂B layers formed at the surface of borided high-alloy steel,” Journal of Materials Engineering and Performance, vol. 21, no. 8, pp. 1714–1723, 2012.

M. Ortiz-Domínguez, J. Zuno-Silva, M. Keddam, O. Damián-Mejía, and M. Elias-Espinosa, “Diffusion model and characterisation of Fe₂B layers on AISI 1018 steel,” International Journal of Surface Science and Engineering, vol. 9, no. 4, pp. 281–297, 2015.

E. E. Vera-Cárdenas, R. Lewis, A. I. Martínez-Pérez et al., “Characterization and wear performance of boride phases over tool steel substrates,” Advances in Mechanical Engineering, vol. 8, no. 2, pp. 1–10, 2016.

M. Ortiz-Domínguez, M. Elias-Espinosa, M. Keddam et al., “Growth kinetics and mechanical properties of Fe₂B layers formed on AISI D2 steel,” Indian Journal of Engineering and Materials Science, vol. 22, no. 2, pp. 231–243, 2015.

T. B. Massalski, Binary Alloy Phase Diagrams, ASM International, Materials Park, OH, USA, 1990.

H. Okamoto, “B-Fe (Boron-Iron),” Journal of Phase Equilibria and Diffusion, vol. 25, pp. 297-298, 2004.

G. Palombarini and M. Carbucicchio, “Growth of boride coatings on iron,” Journal of Materials Science Letter, vol. 6, pp. 415-416, 1987.

I. S. Dukarevich, M. V. Mozharov, and A. S. Shigarev, “Redistribution of elements in boride coatings,” Metal Sciencie and Heat Treatment, vol. 15, no. 2, pp. 64–66, 1973.

G. Ferraro, The Rise and Development of the Theory of Series up to the Early 1820s, Springer, Denmark, 2008.

W. Rubin, Real and Complex Analysis, McGraw-Hill, USA, 1987.

G. Kartal, O. Eryilmaz, G. Krumdick, A. Erdemir, and S. Timur, “Kinetics of electrochemical boriding of low carbon steel,” Applied Surface Science, vol. 257, no. 15, pp. 6928–6934, 2011.

G. Ibrahim, T. Sukru, B. Cuma, Y. Yilmaz, U. Sukru, and K. Yusuf, “Investigation of diffusion kinetics of plasma paste borided AISI 8620 steel using a mixture of B₂O₃ paste and B₄C/SiC,” Sádhaná, vol. 38, no. 3, pp. 513–526, 2013.

S. Sen, U. Sen, and C. Bindal, “An approach to kinetic study of borided steels,” Surface and Coatings Technology, vol. 191, no. 2-3, pp. 274–285, 2005.

I. Uslu, H. Comert, M. Ipek, O. Ozdemir, and C. Bindal, “Evaluation of borides formed on AISI P20 steel,” Materials & Design, vol. 28, no. 1, pp. 55–67, 2007.