Next Article in Journal
Moving Pearl Vortices in Thin-Film Superconductors
Previous Article in Journal
From SrTiO3 to Cuprates and Back to SrTiO3: A Way Along Alex Müller’s Scientific Career
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Modeling of Synthesis Process of Boron Carbide Based Nanocomposites

by
Levan Chkhartishvili
1,2,*,
Levan Antashvili
1,
Lasha Dalakishvili
1,
Roin Chedia
3,
Otar Tsagareishvili
1 and
Archil Mikeladze
1
1
Ferdinand Tavadze Metallurgy and Materials Science Institute, 10 E. Mindeli Str., Tbilisi 0186, Georgia
2
Department of Engineering Physics, Georgian Technical University, 77 M. Kostava Ave., Tbilisi 0160, Georgia
3
Peter Melikishvili Institute of Physical and Organic Chemistry, 31 A. Politkovskaya Str., Tbilisi 0186, Georgia
*
Author to whom correspondence should be addressed.
Submission received: 6 January 2021 / Revised: 14 January 2021 / Accepted: 15 January 2021 / Published: 18 January 2021

Abstract

:
Nanocomposites based on boron carbide B4C are hard materials with wide field of applications in modern technologies. A system of first-order ordinary differential equations that simulates the process of chemical synthesis of nanopowders of B4C–TiB2 compositions containing titanium diboride (TiB2) as an additional phase is suggested and resolved numerically for a typical ratio of reaction constants. Reagents and products concentrations are found as time-functions. In this way, the optimal route of production technology of boron carbide-based nanomaterials can be identified.

1. Introduction

Because of complex properties useful for diverse industrial applications, hard nanocomposites based on boron carbide (with the approximate chemical formula B4C) serve for an important class of advanced materials. In previous works [1,2,3,4,5,6], we developed an effective method of their chemical synthesizing in nanopowdered form from corresponding liquid charges. As is known (e.g., see [7]), optimization of the parameters of chemical technologies is possible through their mathematical modeling. Here, we develop a set of ordinary differential equations describing evolution of concentrations of reactants (reagents, intermediate substances, and products) of chemical reactions, in which the boron carbide/titanium diboride B4C–TiB2 nanocomposites are formed.
In the original technology developed, for boron (B), sources include amorphous boron (a-B), boron oxide (B2O3), or boric acid (H3BO3). A carbon source is some liquid organic compound CxHyOz containing carbon (C), hydrogen (H), and oxygen (O). As for the titanium (Ti) source, usually a titanium dioxide (TiO2) powder is used for this purpose. The reagents are obtained in the process of a multi-stage heat treatment of the aqueous suspension of these precursors.
The B-containing reagent is always a powder of amorphous boron. Thus, the boron source practically does not change when a-B serves as the precursor material. Organic solvent pyrolysis reduces the most part of carbon into carbon black with turbostratic graphite structure. Besides, carbon is emitted with carbon dioxide (CO2). At temperatures higher than ~700 °C, C reacts with CO2 forming carbon monoxide (CO): CO2 + C → 2 CO. Note that the production of a carbide-based composite requires free carbon in excess. Under corresponding conditions, C-containing reagents are not CO2 and CO, but C and CO. As for carbon, during the heat pretreatment, TiO2 is retained in powdered form almost without changes.

2. Formally Simple Chemical Reactions Representing Technology

According to the standard mathematical modeling method [8] of the technological process, the corresponding physical-chemical transformation has to be imagined as a set of (formally) simple, irreversible, homogeneous, and isothermal chemical reactions conducted in a closed volume. Transformation of reagents B, C, CO, and TiO2 into products B4C and TiB2 seems to occur through the following reactions:
C + TiO2 → TiO + CO↑,
3 CO + TiO → TiC + 2 CO2↑,
2 B + TiC → TiB2 + C,
4 B + C → B4C.
Hence, in addition to boron carbide (B4C) and titanium diboride (TiB2) the product is carbon dioxide (CO2) as well. Free carbon (C) being a reagent in Reactions (1) and (4) is an intermediate substance in Reaction (3). Similarly, carbon monoxide (CO) is a reagent in Reaction (2) and an intermediate substance in Reaction (1). Finally, titanium oxide (TiO)—according to Reactions (1) and (2)—andtitanium carbide (TiC)—according to Reactions (2) and (3)—serve only as intermediate substances.
By definition, the order of a chemical reaction is the total number of particles—“molecules”—participating in an elementary collision act. Therefore, for Reactions (1)–(4), the orders are 2, 4, 3, and 5, respectively. As practice shows, the order of a simple chemical reaction cannot exceed 3 and, therefore, Reactions (2) and (4) seem to be complex. However, the establishment of real, i.e., effective, orders for these reactions requires more detailed consideration of their chemical mechanisms.
First, note that the icosahedral cluster B12 is the main structural motif for any condensed phase of elemental boron (as well as boron-rich compounds), including amorphous boron [9]. It remains stable only in the material’s bulk and, upon (in our case, high-temperature) ablation from the powder particle surface, takes a quasiplanar form with broken bonds and, accordingly, greatly increased chemical activity [10]. For this reason, it is natural to think that, not four individual B atoms, but a B12 quasiplanar cluster enter into the reaction. Clusters formation is also characteristic of carbon ablation process. One of the highest mass spectrometric peaks corresponds to C3 clusters [11]. Therefore, the mechanism of Reaction (4) can be represented as the interaction of only two particles: clusters B12 and C3.
Similarly, the mechanism of Reaction (3), in which elemental boron also acts as a reagent, seems to be the interaction of B12 cluster with the (TiC)6 fragment of the titanium carbide structure.
At high temperatures, in a medium containing both carbon and oxygen, various oxycarbon structures are formed. Among the small oxycarbonic clusters [12], only (CO)3 is metastable and, for this reason, hyperactive chemically. Therefore, the mechanism of Reaction (2) can be imaged as the interaction of a metastable cluster (CO)3 with a fragment of the titanium oxide structure in the form of a diatomic TiO molecule.
In its original form, the order of Reaction (1) is already 2. However, it is most likely carried out not by the interaction between C atom and TiO2 molecule [13], but some small carbon cluster, for example the same C3 and the corresponding fragment (TiO2)3 of the titanium dioxide structure.
Thus, the production of B4C–TiB2 nanocomposite powder is (formally) represented by the following combination of simple chemical reactions:
C3 + (TiO2)3 → 3 TiO + (CO)3↑,
(CO)3 + TiO → TiC + 2 CO2↑,
B12 + (TiC)6 → 6 TiB2 + 2 C3,
B12 + C3 → 3 B4C.
Note that all of them are of secondorder, since reagents are represented by pairs of clusters.
Because all reactants (reagents, intermediates, and products) are in dispersive (nanopowdered or even free cluster) or gaseous forms, their mixing is assumed to be good enough to consider distribution of reactants concentrations in the reaction space, i.e., chamber volume, as almost homogeneous.
Now, we have to discuss another key assumption of the model used: chemical closeness of the reaction chamber or lack of exchange of matter with environment. Most reactants are in the solid phase and, consequently, remain in the chamber during whole technology process. There are only two reactants, CO and CO2, in the gas phase. In Reaction (2), gaseous CO produced from the precursor mixture serves for a reagent together with intermediate CO produced in Reaction (1). Thus, at the properly chosen charge composition, no presence of free (not reacted) CO is expected in the system. CO2, the final gaseous product, is carried away from the chamber by an inert gas (e.g., argon) flow. As mentioned above, in principle, carbon dioxide can react with carbon. However, at the properly chosen charge composition, no presence of free (not reacted) C is expected. It means that CO2 left in the system would behave as an inert gas and then dispute the exchange of CO2 gas with environment, the reaction chamber effectively remaining closed.

3. Reactions Rates

From the acting masses law, the rates V 1 , V 2 , V 3 , and V 4 of second-order chemical Reactions (5)–(8) are square functions of corresponding reagents or intermediates concentrations (number of formula unit “molecules” per unit volume) [ B ] , [ C ] , [ CO ] , [ TiO 2 ] , [ TiO ] , and [ TiC ] :
V 1 = k 1 [ C ] [ TiO 2 ] ,
V 2 = k 2 [ CO ] [ TiO ] ,
V 3 = k 3 [ B ] [ TiC ]
V 4 = k 4 [ B ] [ C ] .
By definition, the reaction constants k 1 , k 2 , k 3 , and k 4 are positive and not dependent on the reagent concentrations. However, their temperature-dependence (Arrhenius law) is quite significant, and, when the thermal treatment is not isothermal, reaction constants are also complex time-functions. In the technology under consideration, at the heatingstage, reagents are formed from precursor charge, while the basic chemical Reactions (5)–(8) proceed too slowly. As for the coolingstage, it starts after all the reagents and products are completely consumed and formed, respectively. Therefore, the assumption that the technology is almost isothermal and the corresponding reaction rates are time-constants is acceptable. Such assumption significantly simplifies the system of chemical kinetics equations and, consequently, its solution.
The rate of change of a reagent concentration is opposite in sign to the reaction rate, i.e., negative. On the contrary, the rate of change of a product concentration should be positive. In both cases, the reaction rate magnitude is multiplied by the stoichiometric coefficients of corresponding reagents. Thus, the partial rates of changing of reactants concentrations with time t are:
d [ C ] ( 1 ) d t = V 1 = k 1 [ C ] [ TiO 2 ] , d [ TiO 2 ] ( 1 ) d t = V 1 = k 1 [ C ] [ TiO 2 ] , d [ TiO ] ( 1 ) d t = + V 1 = + k 1 [ C ] [ TiO 2 ] , d [ CO ] ( 1 ) d t = + V 1 = + k 1 [ C ] [ TiO 2 ] , d [ CO ] ( 2 ) d t = 3 V 2 = 3 k 2 [ CO ] [ TiO ] , d [ TiO ] ( 2 ) d t = V 2 = k 2 [ CO ] [ TiO ] , d [ TiC ] ( 2 ) d t = + V 2 = + k 2 [ CO ] [ TiO ] , d [ CO 2 ] ( 2 ) d t = + 2 V 2 = + 2 k 2 [ CO ] [ TiO ] , d [ B ] ( 3 ) d t = 2 V 3 = 2 k 3 [ B ] [ TiC ] , d [ TiC ] ( 3 ) d t = V 3 = k 3 [ B ] [ TiC ] , d [ TiB 2 ] ( 3 ) d t = + V 3 = + k 3 [ B ] [ TiC ] , d [ C ] ( 3 ) d t = + V 3 = + k 3 [ B ] [ TiC ] , d [ B ] ( 4 ) d t = 4 V 4 = 4 k 4 [ B ] [ C ] , d [ C ] ( 4 ) d t = V 4 = k 4 [ B ] [ C ] , d [ B 4 C ] ( 4 ) d t = + V 4 = + k 4 [ B ] [ C ] .
Here [ B 4 C ] , [ TiB 2 ] , and [ CO 2 ] are the product concentrations. Summarizing these partial rates, we obtain the total rates of changing in concentrations of all the reactants:
d [ B ] d t = d [ B ] ( 3 ) d t + d [ B ] ( 4 ) d t = 2 k 3 [ B ] [ TiC ] 4 k 4 [ B ] [ C ] , d [ C ] d t = d [ C ] ( 1 ) d t + d [ C ] ( 3 ) d t + d [ C ] ( 4 ) d t = k 1 [ C ] [ TiO 2 ] + k 3 [ B ] [ TiC ] k 4 [ B ] [ C ] , d [ TiO 2 ] d t = d [ TiO 2 ] ( 1 ) d t = k 1 [ C ] [ TiO 2 ] , d [ CO ] d t = d [ CO ] ( 1 ) d t + d [ CO ] ( 2 ) d t = + k 1 [ C ] [ TiO 2 ] 3 k 2 [ CO ] [ TiO ] , d [ TiO ] d t = d [ TiO ] ( 1 ) d t + d [ TiO ] ( 2 ) d t = k 1 [ C ] [ TiO 2 ] k 2 [ CO ] [ TiO ] , d [ TiC ] d t = d [ TiC ] ( 2 ) d t + d [ TiC ] ( 3 ) d t = k 2 [ CO ] [ TiO ] k 3 [ B ] [ TiC ] , d [ B 4 C ] d t = d [ B 4 C ] ( 4 ) d t = + k 4 [ B ] [ C ] , d [ TiB 2 ] d t = d [ TiB 2 ] ( 3 ) d t = + k 3 [ B ] [ TiC ] , d [ CO 2 ] d t = d [ CO 2 ] ( 2 ) d t = + 2 k 2 [ CO ] [ TiO ] .
The numbers of “molecules” of intermediate substances reacted and produced per unit time, respectively, are equal each to other, so that the total rates of change in TiO and TiC concentrations are zero:
d [ TiO ] d t = d [ TiO ] d t = 0 .
Consequently, the following conditions should be fulfilled:
k 1 [ C ] [ TiO 2 ] = k 2 [ CO ] [ TiO ] = k 3 [ B ] [ TiC ] .

4. System of Chemical Kinetics Equations

Based on the above consideration, we obtain the system of chemical kinetic equations:
d [ B ] d t = 2 k 1 [ C ] [ TiO 2 ] 4 k 4 [ B ] [ C ] ,
d [ C ] d t = k 4 [ B ] [ C ] ,
d [ TiO 2 ] d t = k 1 [ C ] [ TiO 2 ] ,
d [ CO ] d t = 2 k 1 [ CO ] [ TiO 2 ] ,
d [ B 4 C ] d t = + k 4 [ B ] [ C ] ,
d [ TiB 2 ] d t = + k 1 [ C ] [ TiO 2 ] ,
d [ CO 2 ] d t = + 2 k 1 [ C ] [ TiO 2 ] .
If we denote [ B ] N 1 ( t ) , [ C ] N 2 ( t ) , [ TiO 2 ] N 3 ( t ) , [ CO ] N 4 ( t ) , [ B 4 C ] N 5 ( t ) , [ TiB 2 ] N 6 ( t ) , and [ CO 2 ] N 7 ( t ) , this system takes the new form:
d N 1 ( t ) d t = 4 k 4 N 1 ( t ) N 2 ( t ) 2 k 1 N 2 ( t ) N 3 ( t ) ,
d N 2 ( t ) d t = k 4 N 1 ( t ) N 2 ( t ) ,
d N 3 ( t ) d t = k 1 N 2 ( t ) N 3 ( t ) ,
d N 4 ( t ) d t = 2 k 1 N 3 ( t ) N 4 ( t ) ,
d N 5 ( t ) d t = + k 4 N 1 ( t ) N 2 ( t ) ,
d N 6 ( t ) d t = + k 1 N 2 ( t ) N 3 ( t ) ,
d N 7 ( t ) d t = + 2 k 1 N 2 ( t ) N 3 ( t ) .
From these equations, it is clear that total rates of change in concentrations of reagents ( N 1 ( t ) , N 2 ( t ) , N 3 ( t ) , and N 4 ( t ) ) and products ( N 5 ( t ) , N 6 ( t ) , and N 7 ( t ) ) are definitely negative and positive, respectively. Thus, they are monotonously decreasing or increasing functions of time. Accordingly, at the relevant choice of charge composition at the final ( t = ) and initial ( t = 0 ) time moments, respectively, these concentrations are expected to equal zero:
N 1 ( ) = N 2 ( ) = N 3 ( ) = N 4 ( ) = N 5 ( 0 ) = N 6 ( 0 ) = N 7 ( 0 ) = 0 .

5. Initial Conditions

To solve the system of first-order ordinary differential equations describing the reagent and product concentration evolution with time, it is necessary to provide their concentrations at the initial ( t = 0 ) and final ( t = ) time moments, respectively: N 1 ( 0 ) , N 2 ( 0 ) , N 3 ( 0 ) , N 4 ( 0 ) , N 5 ( ) , N 6 ( ) , and N 7 ( ) . Note that some of these values can be related to each other.
Let us introduce the stoichiometric coefficients 1 x and x for the product (B4C)1−x(TiB2)x components, where x < 1 . Thus,
N 5 ( ) N 6 ( ) = 1 x 1 .
Since boron carbide- (not titanium diboride-)based composites are considered, the inequality 1 x > x or x < 1 / 2 must hold. Analogously, if 1 y and y , where y < 1 , are the stoichiometric coefficients of the reagent components: C1−y(CO)y, i.e., COy, the following relation must hold:
N 2 ( 0 ) N 4 ( 0 ) = 1 y 1 .
If N denotes the final concentration of the product’s basic component—boron carbide (B4C)—then:
N 5 ( ) = N .
That of the product’s other component—titanium diboride (TiB2)—is:
N 6 ( ) = x N 1 x .
To find the rest of initial conditions, we need mass balance relations for all reagent and product constituents’ chemical elements—B, C, Ti, and O. As [ B ] 0 initial and [ B ] final concentrations of boron are equal:
N 1 ( 0 ) = [ B ] 0 = [ B ] = 4 N 5 ( ) + 2 N 5 ( ) = 2 ( 2 x ) N 1 x , N 1 ( 0 ) = 2 ( 2 x ) N 1 x .
As [ C ] 0 initial and [ C ] final concentrations of carbon are equal:
N 2 ( 0 ) + N 4 ( 0 ) = [ C ] 0 = [ C ] = N + N 7 ( ) , N 2 ( 0 ) + N 4 ( 0 ) N 7 ( ) = N .
As [ Ti ] 0 initial and [ Ti ] final concentrations of titanium are equal:
N 3 ( 0 ) = [ Ti ] 0 = [ Ti ] = N 6 ( ) , N 3 ( 0 ) = N 6 ( ) = x N 1 x .
As [ O ] 0 initial and [ O ] final concentrations of oxygen are equal:
2 x N 1 x + N 4 ( 0 ) = 2 N 3 ( 0 ) + N 4 ( 0 ) = [ O ] 0 = [ O ] = 2 N 7 ( ) , 2 N 7 ( ) N 4 ( 0 ) = 2 x N 1 x .
From theabove-obtained system of three equations for three unknown quantities N 2 ( 0 ) , N 4 ( 0 ) , and N 7 ( ) , we find their values:
N 2 ( 0 ) = 2 y N ( 1 x ) ( 1 + y ) , N 4 ( 0 ) = 2 ( 1 y ) N ( 1 x ) ( 1 + y ) , N 7 ( ) = ( 1 + x ( 1 x ) y ) N ( 1 x ) ( 1 + y ) .
Finally, summarizing the results, we formulate the initial conditions to the system of Equations (20)–(26):
N 1 ( 0 ) = 2 ( 2 x ) N 1 x ,
N 2 ( 0 ) = 2 y N ( 1 x ) ( 1 + y ) ,
N 3 ( 0 ) = x N 1 x ,
N 4 ( 0 ) = 2 ( 1 y ) N ( 1 x ) ( 1 + y ) ,
N 5 ( ) = N ,
N 6 ( ) = x N 1 x ,
N 7 ( ) = ( 1 + x ( 1 x ) y ) N ( 1 x ) ( 1 + y ) .

6. Converting Problem into Dimensionless Form

It is expedient to analyze the problem in dimensionless form. Let us denote:
κ 4 k 4 k 1 ,
τ k 1 k 4 N t ,
n i ( τ ) N i ( t ) N .
In dimensionless kinetic-parameter κ , time-argument τ , and unknown concentration-functions n i ( τ ) , the system of equations describing chemical kinetics takes the following form:
d n 1 ( τ ) d τ = 2 κ n 1 ( τ ) n 2 ( τ ) 4 κ n 2 ( τ ) n 3 ( τ ) ,
d n 2 ( τ ) d τ = κ 2 n 1 ( τ ) n 2 ( τ ) ,
d n 3 ( τ ) d τ = 2 κ n 2 ( τ ) n 3 ( τ ) ,
d n 4 ( τ ) d τ = 4 κ n 2 ( τ ) n 3 ( τ ) ,
d n 5 ( τ ) d τ = + κ 2 n 1 ( τ ) n 2 ( τ ) ,
d n 6 ( τ ) d τ = + 2 κ n 2 ( τ ) n 3 ( τ ) ,
d n 7 ( τ ) d τ = + 4 κ n 2 ( τ ) n 3 ( τ ) .
The initial conditions for these dimensionless equations are:
n 01 n 1 ( 0 ) = 2 ( 2 x ) 1 x ,
n 02 n 2 ( 0 ) = 2 y ( 1 x ) ( 1 + y ) ,
n 03 n 3 ( 0 ) = x 1 x ,
n 04 n 4 ( 0 ) = 2 ( 1 y ) ( 1 x ) ( 1 + y ) ,
n 5 n 5 ( ) = 1 ,
n 6 n 6 ( ) = x 1 x ,
n 7 n 7 ( ) = ( 1 + x ( 1 x ) y ) ( 1 x ) ( 1 + y ) .
These conditions are formulated such that they guarantee the final concentrations of reagents n 1 n 1 ( ) , n 2 n 2 ( ) , n 3 n 3 ( ) , n 4 n 4 ( ) and the initial concentrations of products n 05 n 5 ( 0 ) , n 06 n 6 ( 0 ) , n 07 n 7 ( 0 ) equalzero in the frames of numerical calculation errors:
n 1 = n 2 = n 3 = n 4 = n 05 = n 06 = n 07 0 .
Note that the subsystem of Equations (38), (39) and (40) can be resolved independently as they contain only three unknown functions n 1 ( τ ) , n 2 ( τ ) , and n 3 ( τ ) . When these functions are found, the rest of the unknown functions are easily determinable from the rest of equations of full system. Excluding functions n 1 ( τ ) and n 2 ( τ ) , we obtain the equation determining the function n 3 ( τ ) :
d n 3 ( τ ) d τ + κ n 3 2 ( τ ) 1 κ + ( n 01 2 n 03 κ n 03 1 κ 1 κ ) n 3 1 + κ ( τ ) κ = 0 .
Its solution is:
τ ( n 3 )   = n 03 n 3 ( τ ) d n κ n 2 1 κ + ( n 01 2 n 03 κ n 03 1 κ 1 κ ) n 1 + κ κ .
In particular, for the parameter κ = 1 / 2 , we can explicitly evaluate this integral:
τ ( n 3 ) = 1 2 n 03 n 3 ( τ ) d n n 2 + ( n 01 2 n 03 2 n 03 ) n 3 / 2 = 2 2 n 03 n 01 4 n 03 ( n 03 n 3 1 + 2 n 03 n 01 4 n 03 ln 1 + 2 n 03 n 01 4 n 03 n 03 n 3 + 2 n 03 n 01 4 n 03 ) .
However, to obtain the explicit function n 3 ( τ ) , we need to invert the latter expression for τ ( n 3 ) ,which is impossible.
Instead, below, we solve the system of ordinary differential Equations (38)–(44) with initial Conditions (45)–(51) numerically at κ = 1 / 2 and; then, to check the correctness of the analytical Solution (55), we calculate that solution at discrete values of n 3 obtained numerically. After that, we build an inverse plot, the discrete set of values calculated by Equation (55) vs. these discrete values of n 3 , together with the numerical solution for n 3 . If Solution (55) is correct, these two plots must coincide in the frames of numerical calculation errors.

7. Numerical Solution Code

To solve Equations (38)–(44) numerically, we first need to fix the dimensionless parameter κ contained in this set of equations. Formula (35) determines it by the ratio of reaction constants k 1 and k 4 .
On the one hand, k 1 is the constant of Reaction (1), in which carbon C reacts with titanium dioxide TiO2. The chain of Reactions (1)–(3) leads to the formation of one of the final products, namely titanium diboride TiB2. A TiB2 molecule is formed per C atom. On the other hand, 4 k 4 is the constant of Reaction (4), in which carbon C reacts with boron B. That reaction directly leads to the formation of other final product, namelyboron carbide B4C. Again, a B4C molecule is formed per C atom.
Thus, obtaining B4C–TiB2 compositions means that above reaction rates meet the following condition: k 1 > 4 k 4 or 0 < κ < 1 (note that the integral in above-obtained Solution (54) is divergent at κ = 1 ). For test-calculations, let us take this parameter at the interval midpoint: κ = 1 / 2 .
We applied Python to solve Equations (38)–(44) numerically. More specifically, a Scipy Library was used to obtain the results and Matplotlibto draw the corresponding graphs. The methods were Odeint and Cumtrapz.
First, we created the function that describes the system and returns the 1 × 3 matrix of derivatives of reduced concentrations n 1 , n 2 , and   n 3 of reagents. Taking κ as the fixed number of 0.5, we used the relation
y = 1 x 1 + x
for loop to change the x parameter from 0 to 1. The y parameter that depends on x changed accordingly.
For each value of parameter x , we have different initial conditions and transmit our function to the Odeint solver as an argument that gives us the solutions written as three column matrices with rows of the same number as dimensionless time-steps. Each column corresponds to the concentration of reagents for the different time values. From these three concentrations, we find the change of the fourth reagent concentration n 4 change in time with the help of Cumtrapz algorithm. We use the same scheme (Cumtrapz) to obtain the results for reduced concentrations n 5 , n 6 , and n 7 of products. After that, we draw the corresponding graphs and then choose the ones that satisfy initial (final) conditions.
In addition, we compared the results to the analytical solution that we defined as another function, Tau, that depends on n 3 and the initial concentrations of the first and third reagents: n 01 and n 03 . It was impossible for us to determine analytically the solutions of n 3 from the equation that contained variable Tau. Thus, to check our results, we did the inverse operation and put the values of its arguments n 3 , n 01 , and n 03 into analytically obtained function Tau and compared it with the numerical solution of the same function. The results show that they matched perfectly.

8. Test Solutions

Table 1 summarizes the initial conditions used in numerical solutions.
The corresponding curves are presented in Figure 1.
As expected, the concentrations of reagents and products are monotonously increasing and decreasing functions of time, respectively. The initial concentrations of all products are zero. As for the final concentrations of reagents, in most cases, they vanish as well. However, at stoichiometric ratio x close to 1 ( x 0.8 ), their saturation levels exceed 0. The explanation seems to be related with impossibility at too high x to meet mass balance relations for all chemical elements participating in the reactions, in particular carbon.

9. Possibility of Method Extending to Other Nanocomposites

To not be restricted only by boron carbide/titanium diboride nanocomposites, let us discuss the generalization possibilities of the developed method of mathematical modeling of the chemical reaction kinetics taking place during production of a nanopowder material. There seem to be two different directions of such analysis: separate production of boron carbide as basic component of this class of composites and simultaneous production of components of other biphasic nanocomposites.
A comprehensive review on boron carbide synthesis and consolidation methods given in [14] can serve for obtaining data explainable/interpretable within the model proposed here for the chemical synthesis of boron carbide-based nanocomposites. The main methods of boron carbide production arecarbo- or magnesiothermic reductions, synthesis from elements or polymer precursors, and vapor phase reactions. However, commercially, boron carbide is produced only by the first two: carbothermic reduction of boric acid in electric furnace or its magnesiothermy in the presence of carbon. Our modeling method is applicable for carbothermic reduction, synthesis from elements, and synthesis from polymer precursors. Thus, it is expedient to characterize them briefly.
The yield of carbothermic reduction is relatively low due to losses of boron due to boron oxides formation. Despite this, carbothermic reduction is usually adopted as a commercial method because of the simple equipment and cheap raw material used. However, it is not suitable for the production of boron carbide enriched in boron isotopes because such a raw material is too expensive. The overall reaction, 4 H3BO3 + 7 C → B4C + 6 CO + 6 H2O, of carbothermic reduction of boric acid (H3BO3 (and boron trioxide B2O3) proceeds in three stages: 4 H3BO3 → 2 B2O3 + 6 H2O↑, B2O3 + 3 CO↑ → 2 B + 3 CO2↑, and 4 B + C → B4C. On heating, H3BO3 converts into B2O3 (releasing water (H2O vapor), reduction of which with CO is thermodynamically feasible above 1400 °C. However, the temperature in the reaction area is usually maintained at higher temperatures, ~2000 °C, to enhance the reaction rate as the boron carbide formation is a highly endothermic process (needing 16,800 kJ/mol).
According to the thermodynamic studies of gas phase carbothermic reduction of boron anhydride, B2O2 and BO are formed from B2O3 in reactions B2O3 + C → B2O2 + CO and B2O3 + C → 2 BO + CO and then reduced to free boron or target product—boron carbide. In particular, B4C nanopowders can be prepared by reduction of B2O3 vapor by carbon black at 1350 °C.
Synthesis from elements is mostly considered as a non-commercial route of boron carbide production due to the high cost of elemental boron. It is preferred only for obtaining chemically very pure or enriched in B-isotopes boron carbides. In this technology, B and C powders are mixed to form an almost homogeneous mixture, pellets of which are then reacted at temperatures of >1500 °C in vacuum or inert atmosphere. There are some reported attempts of the preparation of nanostructure boron carbide in this way.
Synthesis of boron carbide from polymer precursors is an alternative to the above-described high temperature techniques. Some B-containing organic compounds such as carborane(C2BnHn+2), triphenylborane, polyvinyl pentaborane, borazine, etc. on pyrolysis in the temperature range 1000–1500 °C in vacuum or inert atmosphere yield B4C powder. In addition, a free-flowing boron carbide powder can be obtained from boric acid and sugar mixture dissolved in ethylene glycol, dried in air at ~180 °C, and then heated in hydrogen at ~700 °C.
Reaction kinetics and, consequently, quality of B4C product are strongly affected by the feed charge porosity/dispersity, chemical composition, type of carbon used, heating/cooling rates, and treatment temperature and duration in the reactor core. A comparative summary of details of various experimental data summarized in [14] can serve for refinement reaction constants needed for mathematical modeling of the chemical technologies for obtaining boron carbide powders.
As for the biphasic nanocomposites, below, we consider the case of functional compositions of strongly correlated electron oxides, materials with important potential applications in spintronics and energy devices.
Recently, a low-temperature (and thus also low-cost) and scalable sol-based synthesis route has been reported [15] for prototypical bimagnetic LaFeO3–CoFe2O4 (LFO–CFO) nanocomposites. In the corresponding reaction mechanism, both phases are formed simultaneously and the CFO heat of formation assists the LFO phase crystallization. Besides, a symbiotic phenomenon is observed: LFO matrix restricts the growth of CFO particles. In this technology, at first, two separate sols are prepared by dissolving at room temperature the precursors of LFO (La(NO3)3·6H2O and Fe(NO3)3·9H2O) and CFO (Co(NO3)2·6H2O and Fe(NO3)3·9H2O) in distilled water added with glycine acting as the chelating agent. These two solutions are mixed, and then the composite sol is heated to 80 °C. Next, the temperature is increased and maintained at 150 °C until the formation of a gel. Finally, the temperature is increased up to 250 °C, when a self-combustion reaction occurs and yields a powder product.
Multiferroic Pb(Zr,Ti)O3–CoFe2O4 (PZT–CFO) nanocomposites have been prepared [16] in a similar synthesis technique that maximizes the surface contact area and, hence, interfacial coupling between the ferroelectric (PZT) and ferrimagnetic (CFO) phases. First, sols of PZT and CFO are prepared separately. PZT precursor solution is prepared by dissolving at room temperature lead acetate (Pb(CH3CO2)2·3H2O), zirconium propoxide (Zr(CH2CH2CH3O) and titanium isopropoxide (Ti(OCH(CH3)2) in 2-methoxyethanol and acetic acid, and then zirconium propoxide and titanium isopropoxide are added. Similarly, CFO precursor solution is prepared by dissolving at room temperature cobalt acetate tetrahydrate ((CH3COO)2Co·4H2O) and iron nitrate nonahydrate (Fe(NO3)3·9H2O) in 2-methoxyethanol and acetic acid, and then ethanolamine and ethanolamine are added. The mixed precursor solution is first stirred at room temperature and then at 80 °C. Next, the temperature is increased and maintained at ~150 °C until the gel formation. Then, the temperature is increased to ~300 °C, when a self-combustion reaction occurs yielding a powder.
In such biphasic composite, PZT component can be replaced [17] by strontium hexaferrite (SrFe12O19) (SFO) component, also prepared by the sol–gel combustion. First, Fe(NO3)3·9H2O and Sr(NO3)2·6H2O precursors are dissolved in deionized water. Then, a citric acid aqueous dispersion is added at room temperature. NH3 is added drop-wise. Next, sol is heated to 80 °C to form a gel. After the entire solution is converted into a gel, the temperature is increased to 300 °C inducing self-combustion and obtaining the dry powders. This sol–gel synthesis procedure can be directly adapted to simultaneously prepare the exchange-coupled hard/soft biphasic nanocomposite SrFe12O19–CoFe2O4 (SFO–CFO).
Thus, sol–gel self-combustion synthesis is an efficient and reproducible approach to synthesize biphasic nanocomposites. Optimizing of the synthesis conditions decreases the thermal treatment temperature, allows controlling the nanocrystallites size, and, in this way, improves the product quality. These nanocomposites reveal synthesis conditions and charge composition effect on their magnetic properties.
Despite the radical differences in chemical compositions, the obtained technologies of two biphasic nanocomposite powders, such as boron carbide/titanium diboride and strongly correlated electron oxides, show similarities to each other: both materials are obtained as a result of multi-stage heat-treatments of sol charges—precursors solved in water and liquid organic compounds. Therefore, if we write down the kinetic equations for the chemical reactions occurring during the sol–gel synthesis of strongly correlated electron oxides (LFO–CFO, PZT–CFO, SFO–CFO, etc.) [15,16,17], then the mathematical modeling of this technology can be carried out using the method developed by us for the synthesis of boron carbide/titanium diboride nanocomposite powders.

10. Concluding Remarks

In summary, based on the analysis of mechanisms of chemical reactions taking place during synthesis from liquid charge of the nanopowder composites of boron carbide and titanium diboride, a system of ordinary differential equations simulating kinetics of this technology is suggested. B, C, TiO2, and CO are considered as reagents, while the expected products are B4C, TiB2, and CO2.
This set of equations is resolved numerically for a typical ratio of reaction constants at different stoichiometric ratios x of components B4C and TiB2 in target composites (B4C)1−x(TiB2)x. In this way, reagent and product concentrations are found as time-functions. As expected, the concentrations of reagents and products are monotonously increasing and decreasing functions of time, respectively. The correctness of the specially elaborated computing code is proved by perfect matching of an analytical solution with the corresponding numerical one.
This method of mathematical modeling will lead to the optimization of parameters of production technology of hard nanocomposites based on boron carbide. Besides, it can be extended to separate production of boron carbide as basic component of this class of composites, as well as synthesizing other biphasic nanocomposites.

Author Contributions

Conceptualization, L.C. and R.C.; methodology, L.C.; software, L.A. and L.D.; validation, L.A. and L.D.; formal analysis, L.C. and R.C.; investigation, O.T. and A.M.; resources, O.T.; data curation, A.M.; writing—original draft preparation, L.C., L.A. and L.D.; writing—review and editing, L.C.; visualization, L.A. and L.D.; supervision, L.C.; project administration, L.C. and O.T.; funding acquisition, L.C., O.T. and A.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Shota Rustaveli National Science Foundation of Georgia, grant number AR–18–1045: “Obtaining of boron carbide-based nanostructured heterophase ceramic materials and products with improved performance characteristics”.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chkhartishvili, L.; Mikeladze, A.; Tsagareishvili, O.; Gachechiladze, A.; Oakley, A.; Margiev, B. Boron-Containing Nanocrystalline Ceramic and Metal-Ceramic Materials. In Handbook of Nanomaterials for Industrial Applications; Hussain, C.M., Ed.; Elsevier: Amsterdam, The Netherlands, 2018; pp. 13–35. [Google Scholar]
  2. Tsagareishvili, O.; Mikeladze, A.; Chedia, R.; Chkhartishvili, L. Method of Obtaining Boron Carbide Based Hard Nanocomposite Materials. National Center for Intellectual Property of Georgia “Geo Patent”—Patent # GE P 2018 6709 B, 11 October 2018. (In Georgian). [Google Scholar]
  3. Mikeladze, A.; Tsagareishvili, O.; Chkhartishvili, L.; Chedia, R.; Darchiashvili, M. Production of titanium-containing metal-ceramic composites based on boron carbide in the nanocrystalline state. Adv. Appl. Ceram. 2019, 118, 196–208. [Google Scholar] [CrossRef]
  4. Mikeladze, A.; Tsagareishvili, O.; Chkhartishvili, L.; Chedia, R. Obtaining of some boron-containing and related nanocrystalline systems from solutions and suspensions. In Proceedings of the Book of International Symposium on Boron, Nevsehir, Turkey, 17–19 April 2019; pp. 181–191. [Google Scholar]
  5. Barbakadze, N.; Sarajishvili, K.; Chedia, R.; Chkhartishvili, L.; Tsagareishvili, O.; Mikeladze, A.; Darchiashvili, M.; Ugrekhelidze, V. Obtaining of ultrafine powders of some boron carbide based nanocomposites using liquid precursors. Nanotechnol. Percept. 2019, 15, 1–14. [Google Scholar]
  6. Chkhartishvili, L.; Mikeladze, A.; Chedia, R.; Tsagareishvili, O.; Barbakadze, N.; Sarajishvili, K.; Darchiashvili, M.; Ugrekhelidze, V.; Korkia, T. Synthesizing fine-grained powders of complex compositions B4C–TiB2–WC–Co. Solid State Sci. 2020, 108, 1–8. [Google Scholar] [CrossRef]
  7. Bogomazova, N.V. Modeling and Optimization of Chemical-Technological Processes in Industry; Belarus State Technology University: Minsk, Belarus, 2018; 83p. (In Russian) [Google Scholar]
  8. Stromber, A.G.; Semchenko, D.P. Physical Chemistry; Higher Schoo: Moscow, Russia, 1988; 496p. (In Russian) [Google Scholar]
  9. Tsagareishvili, G.V.; Antadze, M.E.; Tavadze, F.N. Obtaining and Structure of Boron; Metsniereba: Tbilisi, Georgia, 1991; 144p. (In Russian) [Google Scholar]
  10. Chkhartishvili, L. Relative stability of planar clusters B11, B12, and B13 in neutral- and charged-states. Charact. Appl. Nanomater. 2018, 2, 1–7. [Google Scholar] [CrossRef] [Green Version]
  11. Lifshitz, C. Carbon clusters. Int. J. Mass Spectrom. 2000, 200, 423–442. [Google Scholar] [CrossRef]
  12. Schröder, D.; Schwarz, H.; Dua, S.; Blanksby, S.J.; Bowie, J.H. Mass spectrometric studies of the oxocarbons CnOn (n = 3–6). Int. J. Mass Spectrom. 1999, 188, 17–25. [Google Scholar] [CrossRef]
  13. Ramana, M.V.; Phillips, D.H. A computational study of the TiO2 molecule. J. Chem. Phys. 1988, 88, 2637–2640. [Google Scholar] [CrossRef]
  14. Suri, A.K.; Subramanian, C.; Sonber, J.K.; Murthy, T.S.R.C. Synthesis and consolidation of boron carbide: A review. Int. Mater. Rev. 2010, 55, 4–40. [Google Scholar] [CrossRef]
  15. Sayed, F.; Kotnana, G.; Muscas, G.; Locardi, F.; Comite, A.; Varvaro, G.; Peddis, D.; Barucca, G.; Mathieu, R.; Sarkar, T. Symbiotic, low-temperature, and scalable synthesis of bi-magnetic complex oxide nanocomposites. Nanoscale Adv. 2020, 2, 851–859. [Google Scholar] [CrossRef] [Green Version]
  16. Kotnana, G.; Sayed, F.; Joshi, D.; Barucca, G.; Peddis, D.; Mathieu, R.; Sarkar, T. Novel mixed precursor approach to prepare multiferroic nanocomposites with enhanced interfacial coupling. J. Magn. Magn. Mater. 2020, 511, 166792. [Google Scholar] [CrossRef]
  17. Maltoni, P.; Sarkar, T.; Varvaro, G.; Barucca, G.; Ivanov, S.; Peddis, D.; Mathieu, R. Towards bi-magnetic nanocomposites as permanent magnets through the optimization of the synthesis and magnetic properties of SrFe12O19 nanocrystallites. J. Phys. D Appl. Phys. 2020. [Google Scholar] [CrossRef]
Figure 1. Calculated time-dependences of reagents and products concentrations at different stoichiometric ratio x of synthesized (B4C)1−x(TiB2)x composite.
Figure 1. Calculated time-dependences of reagents and products concentrations at different stoichiometric ratio x of synthesized (B4C)1−x(TiB2)x composite.
Condensedmatter 06 00003 g001
Table 1. Initial and final reduced concentrations of reagents n 0 i and products n i , respectively, at different stoichiometric ratio x of synthesized (B4C)1−x(TiB2)x composite.
Table 1. Initial and final reduced concentrations of reagents n 0 i and products n i , respectively, at different stoichiometric ratio x of synthesized (B4C)1−x(TiB2)x composite.
x 0.10.20.30.40.50.60.70.80.9
n 01 4.224.504.865.336.007.008.6712.0022.00
n 02 1.001.001.001.001.001.001.001.001.00
n 03 0.110.250.430.671.001.502.334.009.00
n 04 0.220.500.861.332.003.004.678.0018.00
n 5 1.001.001.001.001.001.001.001.001.00
n 6 0.110.250.430.671.001.502.334.009.00
n 7 0.220.500.861.332.003.004.678.0018.00
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chkhartishvili, L.; Antashvili, L.; Dalakishvili, L.; Chedia, R.; Tsagareishvili, O.; Mikeladze, A. On Modeling of Synthesis Process of Boron Carbide Based Nanocomposites. Condens. Matter 2021, 6, 3. https://0-doi-org.brum.beds.ac.uk/10.3390/condmat6010003

AMA Style

Chkhartishvili L, Antashvili L, Dalakishvili L, Chedia R, Tsagareishvili O, Mikeladze A. On Modeling of Synthesis Process of Boron Carbide Based Nanocomposites. Condensed Matter. 2021; 6(1):3. https://0-doi-org.brum.beds.ac.uk/10.3390/condmat6010003

Chicago/Turabian Style

Chkhartishvili, Levan, Levan Antashvili, Lasha Dalakishvili, Roin Chedia, Otar Tsagareishvili, and Archil Mikeladze. 2021. "On Modeling of Synthesis Process of Boron Carbide Based Nanocomposites" Condensed Matter 6, no. 1: 3. https://0-doi-org.brum.beds.ac.uk/10.3390/condmat6010003

Article Metrics

Back to TopTop