Flow and Heat Transfer of Basalt Melt in the Feeder of the Smelter Furnace

Vyacheslav Kremnev, Boris Basok, Boris Davydenko, Andriy Timoshchenko, Aleksandra Timoshchenko

Institute of Engineering Thermophysics of the National Academy of Science of Ukraine, Kiev, Ukraine.

DOI: 10.4236/jamp.2019.711174                

Abstract

The flow and heat transfer of the basalt melt in the boundary layer on a flat plate is considered. The conditions of formation of the layer and the intensity of heat transfer are determined. A self-similar analysis using the symmetry method was used. A system of ordinary differential equations in self-similar form is obtained. The fluid flow and heat transfer of molten basalt at a laminar steady-state flow in the feeder furnaces are numerically researched. The term “protective layer” on the interface “basalt melt-lining” is introduced. The dependences for the calculation of dimensionless shear stresses and the Nusselt number on the lining surface are obtained. The conditions of rational organization of the technological process of basalt melt feeding in the furnace feeder are formulated.

Keywords

Basalt Fiber, Smelter Furnace, Energy Efficiency, Heat Technological Processes

1. Introduction

It is known that the technology of production of high-performance thermal insulation based on staple basalt fiber has the potential for modernization [1]. In particular, due to the rational organization of hydrodynamic and heat transfer processes in the elements of technological equipment. According to current concepts, basalt melt is a multiphase heterogeneous system containing liquid, crystalline and gaseous phases. The structure and action of this system are largely determined by the temperature and pressure at which it is located [2] [3]. The peculiarity of the basalt melt is its ability to restore in time its complex crystal structure from the liquid state when found in a certain temperature range—between the lower and upper limits of crystallization. In relation to the basalts of Ukrainian deposits [4], the critical temperature range—between the lower and upper crystallization limits—is between 900˚C and 1270˚C. The occurrence in the smelter furnace of conditions under which the basalt melt acquires a critical temperature leads to the process of secondary crystallization of the melt. At the point of contact between the melt and the furnace lining, secondary crystallization contributes to the formation of a protective layer that impedes the destruction of the latter—the advantage—guaranteeing the long-term operation of the lining. On the other hand, the occurrence of conditions under which it is possible to recrystallize the basalt melt at the place of forming coarse fibers is a problem that leads to a decrease in the quality of fibers, a violation of the technological process, a decrease in plant productivity and, as a consequence, the increased production price of the final product. Studying the flow and heat transfer of basalt melt on the flat surface gives qualitative and quantitative characteristics of the process, and opens up the possibility of developing technical proposals for process control.

The theme of fluid flow and heat transfer on a flat plate—both in the approximation of the boundary layer, analytically, and in its full formulation, numerically—is widely discussed in numerous publications and monographs [5] [6]. The peculiarity of the formulation under consideration is that the classical phenomenological approach is applied to a deliberately heterogeneous system, with a significant dependence of viscosity on temperature—basalt melt. As will be shown below, the viscosity of the basalt melt in the temperature range under consideration varies by 14 orders of magnitude. The purpose of the study was to quantify the processes of flow and heat transfer of the basalt melt and obtain the characteristics of resistance, heat transfer intensity, heat flux values, and sizes of the zone of secondary crystallization of the melt. The task was solved analytically and numerically.

2. Study of the Flow and Heat Transfer of Basalt Melt on a Flat Plate

2.1. Physical Model

The influence of basalt melt viscosity on the flow and heat transfer in the laminar boundary layer on a flat surface under the condition of the formation of a sedentary liquid layer near the wall is considered. The viscosity of the basalt melt is a nonlinear function that depends significantly on the temperature [7] (Figure 1).

2.2. Mathematical Model

The stationary problem of flow and heat transfer in the boundary layer of a liquid on a flat surface is considered. For this problem a system of equations of hydrodynamics and heat transfer in a boundary layer on a flat plate in a two-dimensional approximation [5] is considered

Figure 1. Dependence of the coefficient of dynamic viscosity on temperature for basalt from the Ivanova Dolina deposit.

$\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0$, (1)

$\rho \left(u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}\right)=\frac{\partial }{\partial y}\left(\mu \frac{\partial u}{\partial y}\right)$, (2)

$\rho c_p\left(u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}\right)=\frac{\partial }{\partial y}\left(\lambda \frac{\partial T}{\partial y}\right)$. (3)

Equations (1)-(3) are analyzed under the following boundary conditions

$y=0:u=v=0;T=T_w$ (4)

$y\to \infty :u=U_{\infty };v=0;T=T_{\infty }$ (5)

To solve the system of differential equations with partial derivatives (1)-(3), we convert them to ordinary differential equations. For this purpose, we use symmetry analysis (analysis of Lie groups) [8]. The algorithm for symmetry analysis for heat and mass transfer and hydrodynamics problems is given in [9]. The transformation of the system of Equations (1)-(3) into a system of ordinary differential equations is considered in the work [10]. The following system of ordinary differential equations in self-similar form was obtained

$\frac{1}{2}\cdot f\left(\eta \right)\cdot f^{″}\left(\eta \right)+M^{\prime }\left[\Theta \left(\eta \right)\right]\cdot {\Theta }^{\prime }\left(\eta \right)\cdot f^{″}\left(\eta \right)+M\left[\Theta \left(\eta \right)\right]\cdot f^{\left(3\right)}\left(\eta \right)=0$, (6)

$\frac{1}{2}\cdot P{r}_{\infty }\cdot f\left(\eta \right)\cdot {\Theta }^{\prime }\left(\eta \right)+{\Theta }^{″}\left(\eta \right)=0$. (7)

The self-similar variable was defined by the expression

$\eta =y\sqrt{\frac{U_{\infty }\rho }{{\mu }_{\infty }x}}$. (8)

The velocity u in self-similar variables are represented in the form

$u=U_{\infty }{f}^{\prime }\left(\eta \right)$. (9)

The velocity v in self-similar variables are represented in the form

$v=\frac{1}{2}\sqrt{\frac{U_{\infty }\cdot {\mu }_{\infty }}{\rho \cdot x}}\left(f^{\prime }\left(\eta \right)\cdot \eta -f\left(\eta \right)\right)$. (10)

Self-similar function for temperature is introduced as

$T=T_w+\left(T_{\infty }-T_w\right)\cdot \Theta \left(\eta \right)$. (11)

The dependence of the dynamic viscosity coefficient of the liquid on the temperature is described by the equation

$\mu \left(T\right)=3.83\times {10}^{-3}\cdot {\text{e}}^{\frac{2.41\times {10}^6}{{\left(\sqrt[4]{T}\right)}^7}}$. (12)

In the self-similar form equation (12) is presented as

$\mu \left(\Theta \left(\eta \right)\right)={\mu }_{\infty }\cdot M\left[\Theta \left(\eta \right)\right]$, (13)

where ${\mu }_{\infty }$ is liquid viscosity outside the boundary layer.

The system of Equations (6), (7) was solved under the following boundary conditions

$\eta =0:f\left(0\right)=0;f^{\prime }\left(0\right)=0;\Theta \left(0\right)=0;$ (14)

$\eta \to \infty :f^{\prime }\left(\infty \right)=1,\Theta \left(\infty \right)=1.$ (15)

2.3. Results

Based on the system of Equations (6), (7) with boundary conditions (14), (15), numerical simulation of fluid motion and heat transfer in the boundary layer over a flat surface in a wide range of parameters μ and Θ is performed. Numerical solutions were obtained using the 4th order Runge-Kutta method [11].

For the numerical solution of system (6), (7), we replace the variables and perform its transformation

$f\left(\eta \right)=y_1,f^{\prime }\left(\eta \right)=\frac{\text{d}{y}_1}{\text{d}\eta }=y_2,f^{″}\left(\eta \right)=\frac{\text{d}{y}_2}{\text{d}\eta }=y_3,\Theta \left(\eta \right)=t_1,{\Theta }^{\prime }\left(\eta \right)=\frac{\text{d}{t}_1}{\text{d}\eta }=t_2.$ (16)

The system takes the form

$\frac{\text{d}{y}_1}{\text{d}\eta }=y_2$, (17)

$\frac{\text{d}{y}_2}{\text{d}\eta }=y_3$, (18)

$\frac{\text{d}{y}_3}{\text{d}\eta }=-\frac{1}{2M\left[t_1\right]}\cdot y_1\cdot y_3-\frac{M^{\prime }\left[t_1\right]}{M\left[t_1\right]}\cdot t_2\cdot y_3$, (19)

$\frac{\text{d}{t}_1}{\text{d}\eta }=t_2$, (20)

$\frac{\text{d}{t}_2}{\text{d}\eta }=-\frac{1}{2}\cdot P{r}_{\infty }\cdot y_1\cdot t_2$. (21)

The system was closed by the following boundary conditions

$\eta =0,f\left(0\right)=y_1\left(0\right)=0,f^{\prime }\left(0\right)=y_2\left(0\right)=0,\Theta \left(0\right)=t_1\left(0\right)=0;$ (22)

the quantities $f^{″}\left(0\right)=y_3\left(0\right)$ and ${\Theta }^{\prime }\left(0\right)=t_2\left(0\right)$ were determined by the shooting method [11], with the aim of satisfying conditions (15).

The grid spacing along η was taken equal to $h=1.0e-06$.

The results of the calculations are presented in the graphical form. The increased viscosity of the medium and its dependence on the temperature on the surface of the plate significantly affects the velocity distribution in the boundary layer. In Figure 2, the comparison of the classical solution of Blasius for a self-similar velocity function $f^{\prime }\left(\eta \right)$ with the solution of the system (6), (7) for different temperatures on the wall surface is displayed. The deformation of the velocity profiles within the temperature boundary layer at decrease of T_w is clearly visible. Curve 1 is the classic Blasius solution obtained for isothermal flow on the plate surface. Curves 2-4 are the solutions obtained for the nonisothermal flow of the melt on the plate at different surface temperatures—1300˚C, 1100˚C, 900˚C, respectively. The melt temperature and the Prandtl number in the free flow are T_∞ = 1450˚C; Pr_∞ = 5180. Accepted conventions are also relevant to Figure 3.

To calculate the shear stress on the wall, the dependence of the second derivative of the self-similar velocity function on the wall $f^{″}{\left(\eta \right)}_{\eta =0}$ on the wall temperature T_w is obtained (Figure 3).

The dimensionless shear stress on the wall is determined by the expression

$\frac{1}{2}{c}_f=\frac{f^{″}{\left(\eta \right)}_{\eta =0}}{\sqrt{R{e}_x}}$. (23)

Tabulated values $f^{″}{\left(\eta \right)}_{\eta =0}$ is presented in Table 1.

Figure 2. Distribution of velocity profiles in the boundary layer at: 1—Blasius solution; 2—T_w = 1300˚C; 3—T_w = 1100˚C; 4—T_w = 900˚C; T_∞ = 1450˚C; Pr_∞ = 5180.

Figure 3. Dependence of the second derivative of the self-similar function on the wall temperature T_w (˚C): T_∞ = 1450˚C; Pr_∞ = 5180.

Table 1. Dependence of the second derivative of the self-similar function on the wall temperature.

Figure 4 shows the dependence of the self-similar temperature function on the self-similar coordinate. The classical solution obtained by Levek is compared with the solution of the system (6), (7) for different temperatures on the wall surface. As the temperature on the wall decreases, the temperature distribution in the temperature boundary layer approaches the linear one, which is typical for non-movable layers.

As already noted, near a flat surface with a temperature T_w < T_∞ basalt melt forms a sedentary layer that protects the surface from destruction. We formalize the concept of “sedentary layer”. We will term the melt layer near the surface, in which the velocity of the medium does not exceed 1% of the velocity in the undisturbed melt flow a protective layer. Using the dependence of the self-similar velocity function on the self-similar coordinate Figure 2, we calculate the value of the protective layer formed on a flat surface for different temperatures T_w. The calculated area is selected taking into account the actual geometric dimensions of the melting plants. The results of calculations are presented in graphical form in Figure 5. It should be noted that in the temperature range T_w = 900˚C - 1100˚C the thickness of the protective layer exceeds the size of the feeder in the melting plant. Based on this, it was concluded that it was not advisable to maintain such temperatures. On the other hand, in the temperature range T_w = 1300˚C - 1400˚C the size of the protective layer of the same order with the geometric dimensions of the feeder in the melting plant and is sufficient to fulfill its technological purpose. The specified temperature range is advisable to choose as the determining, when calculating the value of the heat flux density on the surface of the lining.

The dependence of the first derivative of the self-similar temperature function at η = 0 on the temperature on the wall surface, in order to calculate the intensity of heat exchange on the plate is shown in Figure 6. Traditionally, the heat transfer intensity on a flat surface is described by the dimensionless Nusselt number. In relation to the problem under consideration, the Nusselt number is determined by the expression

$N{u}_x={\left(\frac{\text{d}\Theta }{\text{d}\eta }\right)}_{\eta =0}\cdot {\left(R{e}_x\right)}^{1/2}$, (24)

where ${\left(\frac{\text{d}\Theta }{\text{d}\eta }\right)}_{\eta =0}$ is obtained from the solution of the system (6), (7)

${\left(\frac{\text{d}\Theta }{\text{d}\eta }\right)}_{\eta =0}=0.218\cdot {\text{e}}^{2.269\times {10}^{-3}{T}_w}$. (25)

The average heat flux density on the wall q_ave is determined by integrating an expression for the local heat flux density

$q_x=N{u}_x\cdot \left(\frac{{\lambda }_{\infty }}{x}\right)\cdot \left(T_{\infty }-T_w\right)$, (26)

$q_{ave}=\frac{1}{L}{\displaystyle {\int }_0^L{q}_x\text{d}x}$. (27)

Figure 7 shows the dependence of the average heat flux density q_ave on the temperature on the wall T_w.

As it can be noted, the q_ave dependence on the T_w temperature is nonlinear.

Figure 4. Distribution of temperature in the boundary layer at: 1—Leveque solution; 2—T_w = 1300˚C; 3—T_w = 1100˚C; 4—T_w = 900˚C; T_∞ = 1450˚C; Pr_∞ = 5180.

Figure 5. Dependence of the thickness of the protective layer δ (m) on the x (m) coordinate (for different temperatures on the wall T_w (˚C)): T_∞ = 1450˚C; Pr_∞ = 5180.

Figure 6. Dependence of the first derivative of the self-similar function on the temperature on the wall T_w (˚C): T_∞ = 1450˚C; Pr_∞ = 5180.

Figure 7. Dependence of the average heat flux density q_ave (W/m²) on the temperature on the wall T_w (˚C): T_∞ = 1450˚C; Pr_∞ = 5180.

The resulting nonlinearity is explained by the change in the heat transfer mechanism in the layer. Heat transfer changes from predominantly convective (molar), at a temperature on the surface of the wall close to the temperature of the undisturbed flow, to predominantly conductive (molecular) transfer, at surface temperatures below 1100˚C.

Its maximum value of 813 W/m² q_ave reaches at a temperature on the surface T_w = 1007˚C.

Of practical interest are the values of q_ave corresponding to the temperature range on the wall surface T_w = 1300˚C - 1400˚C. They constitute, respectively, q_ave = 532 - 218 W/m².

3. Conclusions

The concept of “protective layer” in relation to the sedentary area of the melt on the interface “basalt melt-lining” is formulated. The thermophysical conditions determining the formation and size of the protective layer on the interface “basalt melt-lining” are identified. The temperature range on the wall surface is 1300˚C - 1400˚C. The corresponding range of heat flux densities is 532 - 218 W/m².

The dependences for the calculation of dimensionless shear stresses and the Nusselt number on the interface “basalt melt-lining” are obtained.

Acknowledgements

We express our words of appreciation to the directorate and to the plant’s employees LLC CE “Chernivtsi Plant of Heat-Insulating Materials” (Chernivtsi, Ukraine) on the basis of which all experimental-industrial research has been conducted.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Kremnev, V., Basok, B., Timoshchenko, A. and Tymchyshyn, S. (2018) Energy Saving in Construction by Wide Application of High-Quality Insulation Based on Basalt Fibers. Journal of Modern Physics, 9, 1724-1734. https://doi.org/10.4236/jmp.2018.99108
[2]	Dzhigiris, D.D. and Makhova, M.F. (2002) Osnovy proizvodstva bazal’tovykh volokon i izdeliy [Basics of Production of Basalt Fibers and Articles]. Teploenergetik [Thermal Power Engineer], Moscow, 416 p. (In Rus.)
[3]	Ablesimov, N.Y. and Zemtsov, A.N. (2010) Relaksatsionnyye effekty v neravnovesnykh kondensirovannykh sistemakh. Bazal’ty: Ot izverzheniya do volokna. RAN Dal’nevostochnoye otdeleniye. Institut tektoniki i geofiziki im. YU.A. Kosygina, Moskva, 400 p. (In Rus.)
[4]	Trefilov, V.I., Makhova, M.F. and Dzhigiris, D.D. (1992) Mineral’no-syr’yevaya baza gornykh porod Ukrainy dlya proizvodstva volokon. Analiticheskiy obzor [Mineral and Raw Material Base of Rocks of Ukraine for the Production of Fibers. Analytical Review]. Promyshlennost’ stroitel’nykh materialov. Seriya 6. Promyshlennost’ polimernykh, myagkikh krovel’nykh i teploizolyatsionnykh stroitel’nykh materialov. Vypusk 2, 88 p. (In Rus.)
[5]	Shlihting, G. (1969) Teoriya pogranichnogo sloya [The Theory Boundary Layer]. Nauka, Moscow, 387 p. (In Rus.)
[6]	Loytsyanskiy, L. (2003) Mekhanika zhidkosti i gaza: Ucheb. dlya vuzov [Mechanics of Fluid and Gas: Textbook for Universities]. 7-ye izd., ispr. Drofa, Moskwa, 840 p., (Klassiki otechestvennoy nauki). (In Rus.)
[7]	Polevoj, P.P. (1982) Teploobmennye processy pri poluchenii volokna razduvom mineral’nyh nitej vysokotemperaturnymi gazovymi potokami [Heat Transfer Processes in the Production of Fiber Blown Mineral Fibers by High-Temperature Gas Flows]. Kiev: Diss. kand. Nauk, 193 p. (In Rus.)
[8]	Avramenko, A.A., Basok, B.I. and Kuznetsov, A.V. (2003) Gruppovye metody v teplofizike [Group Methods in Thermal Physics]. Kiev. Naukova dumka, 484 p. (In Rus.)
[9]	Avramenko, A.A., Blinov, D.G. and Shevchuk, I.V. (2011) Self-Similar Analysis of Fluid Flow and Heat-Mass Transfer of Nanofluids in Boundary Layer. Physics of Fluids, 23, Article ID: 082002. https://doi.org/10.1063/1.3623432
[10]	Timoshchenko, A., Dmitrenko, N. and Kovetska, M. (2018) Analysis of Flow and Heat Transfer in the Boundary Layer of Fluid with Essential Dependence of Viscosity from Temperature. Thermophysics and Thermal Power Engineering, 40, 13-18. https://doi.org/10.31472/ihe.4.2018.02
[11]	Amosov, A.A., Dubinskiy, Y.U.A. and Kopchenova, N.V. (1994) Vychislitel’nyye metody dlya inzhenerov: Ucheb. posobiye [Computational Methods for Engineers: Textbook.]. Vyssh. shk., Moskva, 544 p. (In Rus.)