CC BY
111
DOI: 10.17277/amt.2018.03.pp.057-062
Ensuring a Temperature Conditions in a Reactor for Functionalization of Carbon Nanotubes
D.V. Tarov*, T.P. Dyachkova, E.N. Tugolukov, I.N. Shubin, V.P. Tarov
Tambov State Technical University; 106, Sovetskaya ul., Tambov, 392000, Russia
* Corresponding author. Tel.: +7 915 660 26 53. E-mail: postmaster@kma.tstu.ru
Abstract
The problem of providing a temperature conditions in a stirrer-equipped capacitive device is considered for the functionalization of carbon nanotubes using titanium stearate in a nanodispersed liquid medium, with a supply of carbon dioxide. The analysis of mixing methods and set-up of the temperature conditions in the reactor was carried out. Based on the differential equations of heat conductivity, a mathematical model of the temperature field of a flow moving in the mode of ideal displacement through a channel of constant cross section formed by the surfaces of the half-tube and the body was developed. The calculation of the temperature field over the current time interval includes multiple solutions of heat conductivity problems with subsequent consideration of all the components of the local heat balance. Based on the results of the calculations, the dependences of changes in the temperature of the nanodispersed medium in the capacitive reactor on the beginning time and the relative half-tube length were plotted.
Keywords
Capacitive reactor; casing; functionalization; carbon nanotubes; temperature field; mathematical model; assumptions; nanodispersed system.
© D.V. Tarov, T.P. Dyachkova, E.N. Tugolukov, I.N. Shubin, V.P. Tarov, 2018
Introduction
Any complex technology passes through a series of stages constituting a kind of life cycle - from the origin of the idea to the immediate application. Forming the reactor shape makes it possible to determine the ways of solving problems, for which the technology of functionalization is developed [1].
The reactor is understood herein as a combination of indicative basic properties and relationships between its elements, which determine its capabilities and mechanisms of its implementation. The reactor is a set of elements and means adapted and technically suitable for protection against negative environmental influences, which are in relationships with each other and form a certain integrity and unity.
In the present work, we will characterize a batch functionalization reactor, the components of which are as follows:
- components reflecting the arrangement of the reactor elements and the nanodispersed (working) medium in space;
- components reflecting the state of the nanodispersed medium, zones of environmental impact, reactor loading/unloading, etc.
Features of functionalization of carbon nanotubes with titanium stearate
The initial (stock) solution of triethanolamine stearate and triethanolamine titanate has a slightly alkaline reaction and is stable at this pH value [2, 3]. This system also acts as a surfactant for dispersion of carbon nanotubes (CNTs). When the pH is lowered, which is achieved by saturation of the nanodispersed medium (solution mixture) with carbon dioxide, the system loses stability, and titanium stearate precipitates on the surface of hydrophilic CNTs. The loss of stability of the system is probably due to protonation of the amine group bound to the chelate complex together with the titanium atom. In this case, the stability of the complex decreases, and the titanium ion binds to the stearate anion and carboxyl groups on the CNT surface.
The triethanolamine is transferred to its carbonate and bicarbonate. Theoretically, all the titanium from the initial solution should be deposited over the CNT surface.
In the reactor considered herein, the nanodispersed mixture containing CNTs was saturated with carbon dioxide at the interface of the contacting phases, without gas bubbling. The amount of the liquid (fluid) adsorbed by the gas depends, first of all, on the properties and size of the contact surface [4, 7].
The paper [5] showed the emergence of destructive changes during the CNT covalent functionalization, which can negatively affect the quality of CNTs and the mechanical and electrophysical properties of composites modified with them. In this case, it is possible to ensure conditions for chemical treatment, under which CNTs with a given value of the degree of functionalization retain the structure of graphene layers.
According to the given scheme, CNTs functionalized with titanium stearate with the required parameters can be obtained by strictly observing the temperature conditions, including the heating of the reaction mass and maintenance of its temperature during the process being implemented [6].
Maintenance of the required temperature conditions of the working medium in the reactor depends on the way it is stirred, namely:
- hydrodynamic mixing by means of a pump and bypass contour providing the feeding of the nanodispersed medium from the bottom of the reactor to its upper zone and subsequent film flow of the medium with adsorption of carbon dioxide;
- stirring the nanodispersed medium by means of a stirrer located inside the reactor.
In the first case, the temperature conditions of the nanodispersed medium in the functionalization reactor are maintained by means of a typical external heat exchanger placed in the bypass loop. From outside, the reactor shell is covered with a thermal insulation layer.
In the second case, the temperature conditions of the nanodispersed medium are maintained by heating the vessel through the heat exchange casing. The latter is also covered externally by the above-mentioned layer. Here, using a capacitive cylindrical device of a vertical design with an elliptical bottom, a mixing device and outer casing seems the most rational option.
Construction of a mathematical model of the reactor temperature conditions
The purpose of mathematical modeling of heat transfer in the reactor casing is to determine the operating modes of the device that ensure the
implementation of the liquid-phase functionalization of CNTs under industrial conditions.
The problem of optimization of thermal processes in capacitive reactors can be formulated as follows. It is necessary to find the surface of the casing (Fc), the form (vt) and the initial temperature (t0cool) of the coolant in the casing, and the duration of reactor operation cycles requiring supply or removal of a large amount of heat (t), under which the coolant operating costs are minimal when using the capacitive reactor.
When compiling a mathematical model, differential equations of heat conductivity are used. The choice of this approach is based on the following assumptions [8-10]:
- the most reliable data on thermal processes taking place in capacitive reactors can be obtained based on the mathematical modeling of temperature fields of a d medium, a coolant and structural elements that affect the course of thermal processes;
- the calculation of local heat fluxes (flows) through the structural elements of the capacitive reactor makes it possible to recalculate the temperatures of the product and the coolant on the basis of local thermal balances under different conditions of their motion (these calculations are feasible if the temperature fields inside the above-mentioned reactor elements are known);
- the spatial non-stationary temperature field is described by the differential Fourier-Kirchhoff equation, which describes the field based on the fundamental laws of heat transfer in space, taking into account all the thermal effects that occur in the operation of the capacitive reactor;
- the partition of the thermal process into spatially-temporal discrete regions eliminates the need for setting up a non-linear heat conduction problem. Hereafter the isolated discrete region will be called local.
Let us consider the derivation of an equation describing the temperature situation in the casing, which provides the operating temperature in the nanodispersed medium when supplying carbon dioxide from above. The coolant moves along the channel formed by the outer wall of the reactor vessel and the inner wall of the casing half-tubes. In deriving the equation, the following assumptions are used:
- the temperature of the coolant along the halftube channel section is constant;
- the coolant represents an incompressible fluid;
- the thermophysical characteristics of the coolant do not depend on the temperature;
- the cross-sectional area of the annular channel is constant.
The assumption of the independence of the coolant thermophysical characteristics on the temperature is due to the insignificant size of the half-tube channel. Let us consider the non-stationary temperature conditions. The coolant flushes the channel walls at two different temperatures varying in time and along the length of the elementary region. The temperature field of the capacitive reactor coated with the half-tube casing presents a combination of the following temperature fields:
- temperature field of the working medium in the reactor;
- temperature field of the coolant flow in the halftube casing channels;
- the temperature field of the thermally insulated sections of the reactor vessel not covered with the halftube casing channels;
- the temperature field of the reactor vessel sections in contact with the coolant flow in the halftube casing;
- the temperature field of the heat-insulated halftube wall.
When developing a mathematical model, the following assumptions are made:
- the working medium in the device is considered in the ideal mixing mode, its temperature varies only in time;
- the coolant in the half-tube casing channel is considered in the turbulent conditions of ideal displacement, its temperature varies along the channel length and in time;
- the above-mentioned temperature fields are considered as sets of local temperature fields for discretized space-time regions, within each of which the thermal physical characteristics of substances and structural materials, as well as the kinetic characteristics of thermal processes, are considered constant, corresponding to the average temperatures of the components in the considered regions.
For each local region, the following assumptions can be made:
- the temperature field of the working medium in the reactor is determined from the heat balance conditions taking into account all operating heat sources;
- the temperature field of the coolant flow in the half-tube casing channels is simulated by the solution of the corresponding Cauchy problem;
- the temperature field of the heat-insulated sections of the device vessel not covered by the halftube casing channels is simulated by the problem of stationary thermal conductivity for a hollow two-layered cylinder;
- the temperature field of the device vessel sections in contact with the coolant flow in the halftube casing is simulated by the problem of stationary thermal conductivity for a hollow single-layered cylinder;
- the temperature field of the heat-insulated halftube wall is simulated by the problem of stationary thermal conductivity for a hollow two-layer cylinder;
- the longitudinal heat flow in the device vessel wall is simulated by the solution of the corresponding boundary value problem.
The temperature field of a flow moving in the conditions of ideal displacement through a channel of constant cross section formed by two surfaces with different temperatures is described by the following
equation
dt (x )
d x
'-+Kt (x)=S,
K=
a1 p +a2 P2 S =a1 P1 tF1 + a2 P2 tF2
Gc Gc
(1)
, (2)
where x is the spatial coordinate in the flow direction, t(x) is the current fluid temperature, G is the mass flow
rate of the liquid, c is the heat capacity of the liquid,
Pt is the perimeter of the i-th channel wall, at is the convective heat transfer coefficient from the i-th wall channel to the liquid, tFi is the temperature of the i-th channel wall, i = 1, 2 is the channel surface index.
Under the initial condition of t(0) = t0 , the solution of equation (1) has the following form:
t (x)=V+(t0-V )exp(- Kx),
y = a1 P1 tF1 + a 2 P2 tF 2
a1 P1 +a2 P2
(3)
(4)
The average temperature of the liquid in the section of length Ax is
A x
t =— f t(x)dx = V +(t° V)(1-exp(-Kx)). (5) A x * K Ax
K Ax
The problem of stationary thermal conductivity for a hollow two-layered unbounded cylinder is as follows:
d 2t, (r)+1 dt, (r) d rf r d r
dt1(R0) d r
i = 1,2, R, _1 < r, < R,; (6)
-a1 (t1 (R0 )-tc1 )=0;
(7)
A,
dt2 R )
d r2
+ a2 (R)-tc2)=0;
(8)
ti (Ri )=12 (Ri ); X! ^=X2
d — • (9) dri d r2
Here, ti(r,) is the temperature field of the i-th layer of the cylinder, oC, as a function of the radial coordinate ri, m; Ri is the outer radius of the i-th layer of the cylinder, m; Xi is the thermal conductivity of the material of the i-th layer of the cylinder, W/(m-K); a1, a2 are the coefficients of heat transfer from the cylinder surfaces to the internal and external media, res pectively, W/(m2-K); tci, tc2 are the medium temperatures from the inner and outer sides, respectively, oC.
The solution of the problem (6) - (9) has the form:
ti (ri )=A + Bi ln(r )-
4 X;
where
B, =
1 f d ^
tc2 - tci +V
ln
R,
V r0 J
+ klnf r2 i + x,
A1 = tci +
q, Ro
X 2 V r1
Ro 1
2 X1 a1
- B
1 1
-+-
V r2 a 2 r0 a1 J
X
(10)
-; (11)
ln(R0 )--^ ; (12)
a1 R0 J
b =xlb + R12 (q2- q1 );
X 2 2 X 2
A2 = tc2 +
q2 R2 f R2 1 ^ f
-+-
2 X 2 a2
- B2
ln(R2 )-
V=-qM(2-R12)- R- (q2 -q1)ln( R2 +q2 R2
4 X2
r1
2 X 2
R
X 2 a2 R2 J
q2 R2
(13)
; (14)
1 j
2 a 2
2 a2 R2
(q2 - q )+-X- (( - R02 )+
4 Xi 2 a>
(15)
The problem of stationary thermal conductivity for a hollow single-layered unbounded cylinder is as follows:
d 2t(r ) 1 dt (r )
dr2 r dr
= 0, R1 < r < R2; (16)
X -a1 (t (R1 )-tc1 ) = 0; dR
XddRr + a2 (t(R2)-tc2 ) = 0.
dR
(17)
(18)
Here t(r) is the temperature field of the cylinder, oC, as a function of the radial coordinate r, m; R1, R2 are the inner and outer radii of the cylinder, respectively, m; X is the thermal conductivity of the cylinder material, W/(m-K); a1, a2 are the coefficients of heat transfer from the cylinder surfaces to the internal and external media, respectively, W/(m2-K); tc1, tc2 are the medium temperatures from the inner and outer sides, respectively, oC.
The solution of the problem (16) - (18) has the form:
t(r)=Aln(r)+B ; (19)
f
B=td + A
A=-
X
Va1 R1
- ln(R 1 )
A
tc2 - tc1
X
a2 R2
- + -
a1 R1
-+ln
f R V
(20)
(21)
V R1
V 1J
The longitudinal heat flow in the vessel wall is determined on the basis of mathematical modeling of the longitudinal temperature distribution in the corresponding structural element.
d 2 T (x ) d x2
- k 2 T (x)=0;
T (0 )=t (0 )-
dT (Lx ) = 0; d x
(22)
(23)
(24)
k =a1 +a 2 ; m =
(a111 + a 212).
T (x)=t (x)-m (25) k
X h
Here, x is the spatial longitudinal coordinate, m; T(x) is the temperature field, oC; t1, t2 are the ambient temperatures, 0C; X is the coefficient of thermal conductivity of the structural material W/(m-K); a1, a2 are the coefficients of convective heat transfer from the external surfaces of the structural element to the environment, W/(m2-K); Lx is the length of the section considered in the heat flow direction, m; h is the thickness of the structural element, m.
The solution of problem (22) - (24) has the form:
where
T (x) = C1 ch (kx)+C2sh (kx),
C1=t (0 )-; k
C2 =-C1 th (kLx ).
(26)
(27)
(28)
X
2
X
2
2
t, °C
90
t, ° c
119 118,5 118 117,5 117 116,5 116 115,5 115 114,5 114
т, min
0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
т, min
Fig. 1. Temperature changes in the reactor
Thus, the calculation of the temperature field of the device for the current time interval includes a multiple solution of the above-mentioned heat conduction problems with subsequent consideration of all the components of the local heat balance. An iterative algorithm involving multiple uses of the solutions of the problems considered is due to the necessity of matching the temperature fields of the interacting components of the system through the appropriate boundary conditions.
In the modeling and design of the capacitive reactor, two main tasks can be distinguished, which include the calculation of the non-stationary temperature field:
- determination of the non-stationary temperature field in the working area according to the geometric dimensions of the reactor, the coolant flow and the initial temperature distribution in all the reactor elements;
- determination of the conditions and design parameters ensuring the functioning of the capacitive reactor. This task is solved most often.
The adequacy of the mathematical model of the capacitive reactor was verified by comparing the calculated data with the results of experiments carried out under laboratory conditions. The deviation of the experimental values from the calculated values does not exceed 10 %. Taking into account the results of the calculation of the non-stationary temperature field in the half-tube casing (number of approaches-4, number of turns in the approach-2), a rector for preparing a nanodispersed suspension was designed to functionalize CNTs having a weight of 300 kg, with the AMT-300 oil used as coolant.
Fig. 2. Change in temperature by the relative length
of the half-tube (0-1) as a function of time (0-50 min)
Fig. 1 demonstrates the calculated temperature changes in the reactor resulting from the calculation (abscissa axis - time, min, ordinate axis - temperature, °С).
The graph shown in Fig. 2 indicates a high uniformity in the heating of the nanodispersed medium and a relatively fast beginning of the conditions (the abscissa axis is the relative length of the half-tube approach, and the ordinate axis is the temperature, °C).
At the initial moment of time, the temperature difference was 5.5 °C for the entire half-tube length, whereat after 50 min, at the time of the beginning of the conditions, it was only 1.8 °C due to the heating of the structure.
Conclusions
When using a capacitive device, a large number of operations are carried out for processing liquid products: heating, cooling, aging of liquids; distillation of volatile fractions; chemical transformations; dissolution of granular and bulk materials; mixing of liquids; and partial loading of components. In real production processes, various combinations of these operations take place. As a rule, all of them occur in non-stationary temperature conditions. The design of the half-tube casing makes it possible to control the temperature conditions in the reactor by disconnecting the section (approach) both at full and partial loading of the reactor volume.
References
1. Dyachkova T.P., Tkachev A.G. [Methods of functionalization and modification of carbon nanotubes]. Moscow, Spectrum, 2013, 152 p. (Rus)
119,5
B0
70
60
50
40
30
20
0
2. Khokhlov P.A., Melezhik A.V., Tkachev AG. [Functionalization of carbon nanotubes by organotitanates]. Vestnik VGU, 2013, issue 1, p. 52-55. (Rus)
3. Tarov D.V., Gurova T.V., Shubin I.N. [Instrumental registration of the functionalization of carbon nanotubes by titanium stearate]. Vestnik TGTU, 2015, vol. 21, issue 2, pp.360-366. (Rus)
4. Kondratyuk I.V., Shubin I.N., Tarov D.V. [Technology for the production of hydrophilic carbon nanotubes]. Fundamental and applied scientific research: a collection of articles of the international scientific and practical conference. In two parts of Part 1. Ufa OMEGA SAINS, 2016, pp. 3-5. (Rus)
5. Dyachkova T.P., Rukhov A.V., Tugolukov E.N., Usol'tseva N.V., Khan Yu. A., Chapaksov N.A. Studies on the changes in the structure of carbon nanotubes grapheme layers after functionalization using Raman spectroscopy-based method. Liquid Crystals and Their Application, 2017, issue 4, pp. 83-89.
6. Dyachkova T.P., Tugolukov E.N., Anosova I.V., Tkachev A.G. Kinetic Study of Functionalization of Carbon nanomaterials. Nanotech France 2015: International Nanotechnology Conference. Pôle
Universitaire Léonard de Vinci, La Défense, Paris, France, 2015, p. 57.
7. Tarov D.V., Solodovnikova S.V., Shubin I.N. A line for the preparation of multifunctional oleophilic multi-walled carbon nanotubes. Increasing the efficiency of resource use in the production of agricultural products - new technologies and technology of a new generation for crop production and livestock: a collection of scientific papers. Tambov: R.V. Pershin's Publishing House, 2015, pp. 227-230 (Rus)
8. Shubin I.N. et al Technological foundations for the functionalization of multi-walled carbon nanotubes with titanium stearate for subsequent production of composites. Vestnik TGTU, 2016, vol. 3, issue 22, pp. 501-509. (Rus)
9. Tugolukov E.N. Solving the problems of thermal conductivity by the method of finite integral transformations in the automated design of technological equipment of the chemical industry. Tambov, TSTU Pulishing House. 2006. 116 p. (in Rus)
10.Tugolukov E.N. Mathematical Modeling of Process Equipment for Multi-Assortment Chemical Production. Moscow, Mashinostroenie-1, 2004, 100 p. (Rus)
