Next Article in Journal
Environmental Assessment of Municipal Solid Waste by Two-Stage Plasma Gasification
Previous Article in Journal
Adaptive Damping Control Strategy of Wind Integrated Power System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Theoretical and Methodological Approaches to the Study of Heat Transfer in Coal Dust Combustion

Institute of Engineering and Technology, South Ural State University, Chelyabinsk 454080, Russia
*
Author to whom correspondence should be addressed.
Energies 2019, 12(1), 136; https://doi.org/10.3390/en12010136
Submission received: 7 November 2018 / Revised: 15 December 2018 / Accepted: 26 December 2018 / Published: 1 January 2019
(This article belongs to the Section I: Energy Fundamentals and Conversion)

Abstract

:
The existing theories of heat transfer in combustion chambers of boiler units fail to take into account a number of important factors that affect the reliability of results, and the methodological approaches to optimizing combustion processes can be revised in view of the spatial and temporal parameters of flame. Hence, the research aimed to improve the fundamental theoretical and methodological principles of studying heat transfer in coal dust combustion in the combustion chambers of industrial steam generators. The authors proposed to extend the theory of heat transfer with the mathematical description of particle size distribution of coal dust. In addition, the authors used the developed mathematical model of coal dust combustion based on a continuous curve of the particle size distribution in the ensemble. The mathematical model is consistent with the aeromechanical and thermal characteristics of flame. This work introduced a concept of flame continuum as a continuous medium, where the processes of combustion and heat transfer are studied. To achieve the research aim, in this paper, the methods of combustion chamber zoning, the equations of stationary and non-stationary heat conduction, radiation, and convective heat transfer, were used. These methods were tested on a number of high-temperature units.

1. Introduction

Fuel combustion and heat transfer in the sections of a boiler unit (BU) are optimized and harmful emissions are minimized at the input of the BU gas path and in the combustion preparation system [1,2,3,4]. In turn, the BU arrangement and the fuel type must be adjusted to the system for preparing and supplying combustion components, the design of fuel-burning appliances (FBA), and the design of the intense combustion zone (ICZ). Thus, the system for preparing and supplying the combustion components and the FBA design provide the conditions for stable ignition and effective fuel burn-out in the ICZ, which is related to the flame temperature and heat transfer in the ICZ that is described in the paper [5]. There is a similar description in the paper [6]. The composition of solid fuels was studied in [7], the fluidized bed fuel delivery system was described in [8].
Mathematical models, on the basis of which researchers conduct computer simulations, are shown in the papers [8,9]. In the paper [8], computer fluid dynamics (CFD) was considered in a specialized program based on a mathematical model developed by the authors of the paper [8]. In the paper [9], the simulation of the jet outflow process consisting of a mixture of air and coal dust was considered. The process of fuel preparation is described in some detail, namely it’s grinding, mixing, and feeding into the furnace through low-emission burners. In this part, the physical modeling is well analyzed. In addition, in the same paper, the authors consider mathematical modeling, on the basis of which the computer modeling process is built.
The system for preparing and supplying the combustion components and the FBA design provide the conditions for stable ignition and effective fuel burn-out in the ICZ, which is related to the flame temperature and heat exchange in the ICZ that is described in the paper [5]. There is a similar description in the paper [6]. The composition of solid fuels was studied in [7], the fluidized bed fuel delivery system was described in [8], and the influence of the design of burners on the formation of nitrogen oxides (NOx) and underburning of fuel was shown in [9].
Studies made by authors such as L. Rosendahl support the theory of weak methodological assurance of heat exchange processes in high-temperature installations; this researcher, together with M. Mando, published an article with his own mathematical model [9]. In the article [8], T. Asotani and his co-authors also consider a model of the behavior of coal particles mixed with air, and the proposed mathematical model correlates well with [9], as well as with some propositions from the authors of [5]. The mathematical models proposed in the review do not contradict each other, as well as the fundamental laws of physics, and heat transfer in particular.
The widely applied analysis of the particle size distribution of fuel dust is based on calculating the total residue according to the Rosin-Rammler formula Rδ = exp(−n), where b and n are constant coefficients which characterize the comminution fineness and the homogeneity of the grain-size distribution of fuel dust, respectively. Application of the total residue characteristic Rδ facilitates reducing all dust fractions to 1.0, taking into account the losses in the analysis, which is consistent with the mass conservation law. This analysis is also consistent with the accepted standard method for sizing disperse solid particles [10,11] and is widely applied.
But, when moving in the flame continuum during heating, volatilization, and burning of volatiles and coke, the fuel dust flows in as a mixture according to its particle size, small-sized particles of δmin being the first and large-sized particles of δmax being the last. Therefore, to analyze the thermal processes in the flame continuum, it is necessary to relate the analysis of the total residue to the local particle size value in the ensemble of fuel particles F(δi).
The generally applied method is based on the analysis of particle size distribution F(δ1/δ2) in the range of δ-value between δ1 and δ2, which corresponds to the discrete process of sizing the fuel dust. To obtain a continuous dependence of F(δi), it is necessary to choose a representative point for the interval (δ1;δ2). The method for particle size analysis recommended in [10,11], where F(δ + Δδ) is divided by the variation interval in the particle size distribution from δ to Δδ, that is, by Δδ, is not adequate to the research task, since, according to the definition, F[δ/(δ + Δδ)] = R(δ) − R(δ + Δδ), and
F [ δ / ( δ + Δ δ ) ] Δ δ = R ( δ ) R ( δ + Δ δ ) Δ δ = exp ( b δ n ) exp ( b ( δ + Δ δ ) n ) Δ δ .
For the value of the total residue R(δ) that monotonically decreases with the increasing δ, Formula (1) yields a relative decrease in the total residue in the interval [δ, δ + Δδ], which is conventionally determined in fractions of the total fuel mass. This means that Formula (1) yields the variation in the total residue with the change Δδ. This is the rate of change of the total residue R(δ) in the interval of variation δ with a change of this interval, which is almost unrelated to the change in the residue of coal dust particles in this interval.
Our developed methods were tested on a number of high-temperature units, including boiler units with natural circulation, and direct-flow boilers, when burning fuels of different thermalphysic characteristics. As a result, we proved that the parameters of the heterogeneous flame depend on the mathematical description of the particle size distribution of coal dust; the authors determined the characteristics of the combustion process that previously had not been taken into account in the theory of heat transfer; the authors provided theoretical evidence for new methodological approaches to the theory of heat transfer in the furnace of a boiler unit. In addition, the authors theoretically confirmed the possibility of considering a flame as a continuous medium and showed the prospects of determining the adiabatic temperature distribution along the length and height of the flame continuum in the intense combustion zone in the combustion chamber of a boiler unit. The conceptual issues of gas-solid heat transfer with respect to the conditions of the combustion chamber of the boiler unit have been solved. Hot flue gases and screen tubes are considered as part of a single continuum environment.
High-temperature installations should include boiler units with natural circulation, direct-flow boilers. It should be noted, the analysis of the above scientific works showed that when burning fuels with different thermal characteristics in high-temperature installations, there is currently no single methodological basis for studying heat transfer processes during exothermic combustion reactions in boiler units. Researchers consider the boiler unit as a set of individual elements that are not related to each other. Thus, there is a scientific basis for the development of new theoretical and methodological approaches in the study of heat transfer processes in high-temperature installations. Therefore the aim of our investigation is to improve the fundamental theoretical and methodological principles of studying heat transfer in coal dust combustion in the combustion chambers of industrial steam generators.
According to known mathematical dependencies, for example [12], the confidence values of the argument of a monotonically decreasing function can be determined in various ways. At the suggestion of the authors, the function should be presented as continuous, that in the sieve analysis neither in theory nor in practice, the researchers did not apply. The continuous function allows operations such as differentiation and integration.

2. Problem Statement

It is rational to use the arithmetic mean value of two intervals F(δ1/δ2) and F(δ2/δ3) to determine the representative value of the residue of coal dust particles in the interval of size change δ. The authors should bear in mind that there is a relation ysmygmyarm between the arithmetic mean of any value yarm, the geometric mean value ygm, and the harmonic mean value ysm [12,13]. Therefore
F i = 0.5 [ F ( δ 1 / δ 2 ) + F ( δ 2 / δ 3 ) ] = 0.5 [ exp ( b δ 1 n ) exp ( b δ 3 n ) ] .
The representative points Fi thus obtained make it possible to plot a continuous curve of the particle size distribution of the fuel particle ensemble, which allows application of differential and integral calculus to mathematically describe the behavior of particle ensemble in a furnace. The reliability of this approximation, that is, the adequacy of replacing the step function F[δ/(δ + Δδ)] with the continuous curve Fi, increases due to the increase in the number of measurements taken with screens for sizing fuel dust. It is also noteworthy that a sufficiently large number of particles correspond to a single point in the particle size distribution; since 1% of the weighed sample mass, equal to 25 g and 200 μm, corresponds to 330 particles of this size. It is also noteworthy that a sufficiently large number of particles correspond to a single point in the particle size distribution; since 1% of the weighed sample mass, equal to 25 g and 200 μm, corresponds to 330 particles of this size, that is confirmed by experimental studies in the laboratory of Heat Power Engineering of the scientific and technical Center under the Ministry of energy of the Kyrgyz Republic.
Figure 1 shows a screen for screening coal dust. Figure 2 shows a rotating coal (biofuel) vibrating screen.

3. Development of Fundamental Mathematics

The qualitative analysis in the framework of the special probability distribution theory [14,15] shows that the particle size distribution of solid fuel dust can be described by the equation of hypergeometric distribution, which leads to a binomial distribution if the number of measurements is increased. As shown above, when transformed to the continuous probability distribution, it corresponds to the Gauss curve with the density of normal distribution of probability to deviate from the mean value
ϕ ( u ) = exp ( 0.5 u 2 ) / 2 π ,
where the standard normal deviation u = (δξ)/(ξ·σ) is introduced as the argument; ξ corresponds to the particle size δ marking the maximum on the curve Fi = f(δ); σ2 is the dispersion, σ is the deviation standard. Normalization, in this case, means placing the ensemble of fuel particles according to their number near the maximum ξ and reducing all values to the dimensionless form; standardization is reducing the ensemble characteristics to the form of dependence on the deviation standard σ. The error function φ(u), which in this case should be referred to as the deviation function, is tabulated and given in the reference books [14,15]. Figure 3 and Figure 4 show the deviation function for the conditions of the research problem with help of different graphs. It is noteworthy that the standard deviation as a generalizing argument uniquely determines function (3) and the values u = 1 and u = −1 correspond to inflection points on the curve plotted according to the dependence (3). σ2 = (δξ)2/(ξ·u)2. Furthermore, the curve φ(u) plotted in accordance with the particle size distribution makes it possible to easily determine the dispersion σ2 = (δξ)2/(ξ·u)2.
The authors calculated the correlation coefficient ry1,y2 that characterizes the tightness of the close linear connection between y1 and y2 to determine how reliable the equation of the Gaussian curve (3), plotted according to the proposed technique ycalcy2, was for the description of the practical particle size distribution of fuel dust ypry1. Thus, it is rational to subject the dependence (3) to the logarithmic linearization
y 2 = ln ϕ ( u ) = ( 0.92 + 0.5 u 2 )
In the papers [14,15], they try to solve similar mathematical problems, as in the sieving of particles. Having determined the correlation coefficient [16,17] −1.0 ≤ r1,2 ≤ 1.0 between the practical results of fuel dust sizing and the Gaussian curve in a conventional way [18], we quantified the reliability of description (3). This technique can also be applied at the F(δi) determination stage. In [19] a solution was proposed to solve similar problems using computer simulation, and in [20,21] the authors approached the exact definition of the desired quantities using specific methods, but they obtained quite realistic and practically applicable results.
The ensemble of fuel particles with characteristics (1)–(3) moves with the air flow at the velocity wl that varies along the flame length l due to the flame expansion and the change in its temperature and dimensions through fuel heating and burn-out. The flame expansion can be registered by the drop in the average kinetic energy El/E0 = 0.56d0/0.15l of the circular jet with the exit diameter of the fuel-burning device d0 [22,23]. Then the local velocity of the flame continuum without taking into account the change in its temperature is the following wl2 = 3.73w02d0/l, whence it follows that wl = 1.93w0(l/d0)−0.5. Similar dependencies are proposed for flat and swirling flames.
The temperature factor in the velocity variation can be taken into account as a first approximation basing on the following relations. Since El/E0 = ρw20w02, the temperature dependence can be determined as El,T/E0 = (ρ0w0wT)/(ρ0w0w0) = wT/w0 when the mass ρw is released in the flux density ρw2, which is constant according to the continuity equation ρw = ρ0w0 = const, and hence
w l , T = 1.93 w 0 ( T / T 0 ) ( l / d 0 ) 0.5 .
Thus, within the scope of the stated problem, the velocity of the flame continuum increases in direct proportion to the increase in its absolute temperature and decreases in inverse proportion to the square root of the relative length in calibers l/d0. A more accurate consideration of the temperature factor is possible when calculating the heat of advective mass flows produced by the ejection effect of the flame and the heat released during fuel combustion.
We present experimental data for the operation of boilers E-210-140 and E-160. For the boiler, E-210-140 data are given at the nominal load of the boiler Dnom = 53 kg of steam per second, in addition, when burning coal with the heat of combustion of Kazakhstan 24,560 kJ/kg. Measured in the experiment parameters: speed w0 = 8 m/s, T0 = 1371 K, d0 = 0.9 m, T = 1721 K, l = 5.33. For the e-160 boiler, the data are given at the nominal load of the boiler Dnom = 44 kg of steam per second, in addition, when burning Kyrgyz coal with a combustion heat of 22,450 kJ/kg. The measured experimental parameters: velocity w0 = 9.2 m/s, T0 = 1389 K, d0 = 0.86 m, T = 1685 K, l = 4.75. After that, we calculate the desired speed of the flare continuum within the framework of the problem (Table 1).

4. Adapting the Model to Flame Characteristics

When adapting the relations of the multiple-size fuel flux to the aeromechanics of the flame continuum, it is necessary to define two asymptotic conditions that are associated with the initial and final flame sections. The formula ulim = (δlim − ξ)/(ξ·σ) is used to determine limiting minimum value of coordinate ulim, corresponding to the size of the fuel dust particle δlim, which does not yet participate in heating and ignition, and is not involved in calculation. If we place the flame cross-section l f = 0 at the point where the fuel flux starts to mix with the air and assume that δlim = 0 in this cross-section, then ulim = −ξ/(ξ·σ), that is, the location of the initial cross-section of the flame continuum depends on the characteristics of particles of the fuel dust ensemble.
The fuel particles of small thermal mass (the number Bi ≤ 0.1) are heated practically without temperature drop in the particle cross-section, therefore the dimensionless temperature Θ of a particle with the diameter δi can be described by the dependence
Θ = exp ( 3 B i F o ) ,
thermal conductivity, density, and mass heat capacity of the particle material in regulated units. Here Bi = αλ/δ, Fo = /δ2, where α is heat transfer coefficient, W/(m2K), λ is thermal conductivity, W/(m K), a is thermal diffusivity, m2/s, τ is time, s, δ defines the size of the surface, m.
The numerical analysis shows that the Bi number can exceed 0.1 when large fuel particles with a high heat transfer rate are heated. The effect produced by the thermal mass of particles on the temperature field (6) can be taken into account by introducing the thermal mass coefficient.
The thermal mass coefficient m of the fuel particle is determined by the relation of the heat-transfer coefficient on the fuel particle surface α to the heat-transfer coefficient k, including only a part of the internal thermal resistance and taking into account the particle shape factor ς
m = 1 + 0.5 B i / ( ς + 2 ) .
To a spherical fuel particle approximation ς = 3 and m = 1 + 0.1Bi with the maximum value of Bi = 1.0, the thermal mass coefficient in this case is m = 1.10. If we take the thermal conductivity coefficient of a coal particle equal to λ = 0.15 W/m·K, we can obtain Bimin = 0.03, m = 1.003 for a particle with a size of 50 μm = 0.05 × 10−3 m and a heat-transfer coefficient per particle surface of α = 100 W/m2·K, which refers to the field of thermally thin bodies. Under the same conditions, a particle of maximum size δ = 1000 μm = 1·10−3 m will correspond to the Bimax = 0.67, while the deviation from the calculation of the particle temperature field according to the laws of thermally thin bodies will reach 6.7% with mass ratio m = 1.067. The time of heating to the specified temperature of intensive volatilization should be calculated according to the formula obtained by taking the logarithm of dependence (7)
τ i n t = [ δ ρ c m / ( 2 α ς ) ] ln ( θ 0 / θ i n t ) ,
where θ0 = tgt0 is the difference in the temperatures between the flue gases tg and the initial temperature of the particle t0; θint = tgtint is the difference in the temperature between the flue gases and the temperature of intensive volatilization. As the first approximation, tint can be considered the ignition temperature of the combustible mixture. When tint increases, the temperature difference θint decreases and the required time τint increases; according to the results of [24], the drying time of the fuel particle can be ignored.
The gradual start of the fuel particles into the heat transfer processes in the flame and the additive dependence of the total heat content are defined by the integral of the distribution function of the deviation probability (3), and it is rational to set the integration limits from minus ∞ to the maximum fuel particle size δmax corresponding to the upper limit of integration U
Φ ( δ ) = 1 σ 2 π δ max exp 0.5 [ ( δ ξ ) / ( σ ξ ) ] 2 d δ = Φ u ( U ) .
Moreover, the normal distribution function Φ(U) is associated with the probability integral Φ(U) and the error function erf(U/20.5) by the dependence
Φ i ( U ) 1 / 2 = 0 , 5 e r f ( U / 2 ) = Φ ( U ) .
Thus, the normal distribution function Φi(U) can be defined by the probability integral Φ(U) or by the error function erf(U/√2), both dependences are tabulated [14,15]. But the practical adaptation of these dependences has a number of features that are related to the reliability of the dependences used. These features can be determined by the numerical analysis of the dependences (8) and (9) when changing the variability intervals of the arguments and their functions. Figure 5 shows the distribution function Фi(U).
The normal distribution law is characterized by a distribution function, which is determined in this case via the tabulated Laplace integral or via the tabulated error function. Figure 5 shows the normal distribution function through the tabulated error function, which is convenient for interpretation by scientists in the field of sieve analysis of particle sieving. The lower limit of integration in the definition of the error function is determined by the fractional composition of fuel particles in the region of small values of δ, and the upper limit is associated with the calculation of the incomplete combustion of fuel in a given length of torch. The specified range of values on the axes corresponds to the most common case of dust sieving with fractional composition R90.
e r f ( U / 2 ) = 2 / π 0 π exp ( [ U / 2 ] 2 ) d u .
The numerical analysis shows that when the probability integral (9) is used for the range of positive and negative values u = (δξ)/(ξ·σ), the difference between Φ(U) and Φ(δ) does not exceed 0.10% in the whole variation interval u. When calculating Φu(U) according to the dependence (10), the difference may exceed 10% in the range of negative values u, which is related to the lower limit of integration when obtaining the error function.
The resulting integral function is related to the known probability integral [14,15], which is given in tables in reference books; it makes it possible to calculate not only the total residue of the known particle size distribution, but also to define a number of properties of the multiple-size fuel flame. Thus, the reduced dependence (8) makes it possible to determine the start of a certain temperature regime depending on the particle size and the conditions of its motion with the multiple-size flame. In order to determine the burnout of fuel particles of different sizes, the adopted model should be complemented with the conditions under which the burning particle interacts with the flux of the furnace and flame medium.
If we assume that the interaction pattern of a burning fuel particle with the furnace atmosphere [23,24] has a surface model, which corresponds to a plane diffusion layer between a spherical particle and the furnace atmosphere, the burnout rate of the particle can be determined at Gc = const according to the equation dδ/dτ = −2McGc, where Mc = 12 kg/kmol is the molar mass of carbon, ρ ≈ 800 kg/m3 is the density of coke. The flow of burnable carbon in the flat diffusion layer, the most common model layer in simulation of moving burning particles of coal dust is assumed to be a constant to simplify the calculations, which does not contradict the actual conditions of combustion, because the boiler unit in a certain mode for a long time period supply a fuel-air mixture constant, and the fractional composition is determined by the operation of the mill devices and also with a certain mode of operation the boiler unit remains virtually unchanged.
The adopted pattern of surface one-dimensional burning in the primary reaction with CO2 generation has the first order in oxygen; combustion is considered quasi-stationary and self-similar in time. In this case, the total burnout time of a particle with the initial size δ0 is the following:
τ = 1 2 M c δ 0 δ ( ρ / G c ) d δ = 0.0417 δ δ 0 ( ρ / G c ) d δ
If we substitute the integral expression of time τ in (7) for the process associated with heating and burning of fuel in the flame, we can obtain the distribution function of heat release in the flame along its length. But this pattern requires refinement both in terms of assuming spherical surfaces to be flat, and in terms of the dependence of Gc on δ.
To take into account the polydispersity of the medium in the flame continuum, we need to apply the dependence for the single-particle burnout function associated with the local value of the kinetic and diffusive characteristics of the process. As a first approximation, we can neglect the effect of CO burning within the boundary layer, the reducing reactions, and the internal response [23,24]. Then, a one-dimensional model of the combustion process of the first-order oxygen reaction, with oxygen being insufficient, has the dependence for the burnable carbon flux Gc: Gc = ωCox, where ω is the kineticity parameter of the fuel burnout process, determined in accordance with Kirchhoff’s second electrical rule for a parallel circuit of two conduction bands. The oxidation function is determined by the amount of oxygen in the flame gas 0.21Pg/RT and its consumption for the combustion of volatile substances and carbon in stoichiometric ratios.

5. Dynamics of Particle Combustion

The dynamics of multiple-size fuel combustion is associated with the distribution of temperature and concentration of components along the length of the flame continuum; the application of electrothermal analogy and Kirchhoff’s electrical rules for nodes and circuits allow us to relate a number of factors. Marking with nodes the points on the burning particle surface without oxidizer Cox = 0 and those on the outer surface of the diffusive boundary layer, where the concentration of the oxidizer is equal to that of the flame continuum Coxf, we obtain the difference in oxidizer transport potentials ΔCox = Coxf 0 = Coxf. Arranging the distributed diffusive resistance 1D parallel to the lumped kinetic resistance 1/k1, we obtain the total diffusive-kinetic resistance to the combustion process RΣ = 1D + 1/k1 = 1, where the diffusive-kinetic conductivity ω = 1/RΣ = αDk1/(αD + k1).
The factor ω can also be referred to as a kinetic index of the fuel burnout process ωk for αDk1, when the diffusion processes do not limit the burning rate, ωk = 1.0; for αDk1, when the diffusion processes become limiting, ωk → 0. For αD = k1 ωk = k1/2, which corresponds to the intermediate regime; thus, applying the kinetic index of the combustion process ωk we can formally designate the boundaries of the kinetic and diffusion regimes including the intermediate region, which corresponds to the trans-regime representations in the dynamics of the solid-fuel burnout process.
The local concentration of oxygen in the flame continuum Coxf, kmol/m3, is determined by the initial volumetric concentration of Cox.i, associated with the combustion method and the burner with air blast Cox.i = 0.21, and with oxygen consumption for the combustion of volatiles and coke according to the oxygen balance
C o x f = C o x . i P T V k 0 R T V g [ ( α 1 ) V 0 V k 0 + K H p ] .
(13) takes into account that the stoichiometric air consumption V0, m3/kg, is equal to the air consumption for the combustion of coke and volatiles V0 = Vk0Kc + Vv0Vvc. The oxygen content according to (13) ensures the local burnout of the fuel particle with the initial size δi with respect to this dependence without taking into account the variability of the variability of the kinetic factor [23,24]
d δ i d τ = 2 M ρ k ω k C o k f .
The diffusion coefficient of mass delivery αD is generally determined by the dependence of NuD on numbers Re and Pr, where the former includes the relative velocity of a particle in the flux. When the flame flux develops, the relative velocity and the Re number reduce to zero, with NuD ≈ 2. Thus, the kinetic index of the process of fuel particle burning under these conditions is defined by the formula
ω k = k 1 N u D D / δ i / ( k 1 + N u D D / δ i ) = 2 k 1 D / δ i / ( k 1 + 2 D / δ i ) .
Here NuD is Nusselt diffusion criterion (dimensionless), δi is the size of the i-th particle, m, k1 is kinetic mass transfer coefficient, m/s, D is diffusion coefficient, m2/s, ω k is process kinetics index, m/s.
Taking into account (13)–(15), the combustion velocity of a fuel particle with size δi is determined by the formula
d δ i d τ = 2 M ρ k 2 k 1 D / δ i ( k 1 + 2 D / δ i ) C o k . H P T V k 0 R T V g [ ( α 1 ) V 0 V k 0 + K H p ] = ω k B T .
Depending on (16), parameter B includes all the constants of the specified task, the proportion of unburned coke Kur is assumed to be the same for the particle of any size, the temperature of the flame continuum T can be related to the abscissa X and the flame length, thus the effect of the particle size distribution consists in the factor ωk. This allows separation of the differentials and integration of the dependence (16) for obtaining the burnout function of fuel particles of size δi.
τ v y g = T B δ i 0 0 2 D + k 1 δ i 2 k 1 D d δ i = T B ( δ i 0 k 1 + δ i 0 2 4 D )
The change of sign in (17) is connected with the change in the limits of integration from 0 to δi0, which corresponds to the dynamics of fuel particle burnout. The temperature of the flame continuum, which affects the time of fuel particle burnout, is determined here by the approximation from Tint = tint + 273 to Tcon, corresponding to the final flame cross-section l f = 1.0
T = T i n t + l f ( T c o n T i n t ) .
The longitudinal dimensions of the flame continuum can be determined using different thermal physic techniques of diagnosis in the field of aeromechanics, temperature and radiation, concentration and analysis, as well as using techniques of mathematical statistics. The method for diagnosis should be selected in accordance with the main research task. In this paper, determination of flame characteristics is based on the dispersion analysis of fuel dust, so it is rational to determine the initial flame section l0 by the ignition time of the smallest fuel dust particles according to (8). In this case, l0 will differ from the distance to the burner port by the amount that is required to heat and ignite the smallest dust particles.
The main task under the stated conditions is to determine how much the initial velocity at the nozzle outlet F0 = πr02 is influenced by three main factors: jet expansion due to the involvement of the ICZ medium into the flux, change in velocity due to the change in particle size distribution of the flux caused by combustion and temperature change of the flame medium from T0 to Tf. Thus, the change in velocity of the flame medium from l f = 0 to l f = 1.0 during its motion and combustion process can be represented by the dependence
w f = w 0 k F k V k T ,
where kF, kV, kT are the influences of change in flame cross-section, change in volumetric flow in the flame, and its temperature change, correspondingly:
k F = 1 ( 1 + 0.21 l f / d 0 ) 2 ;
k V = 1 [ 1 + l f ( V g α V v α 1 ) ] ;
k T = T f / T 0 = 1 + Δ T g o r l f / T 0 .
The relation (20) is obtained geometrically on the basis of central angles of the heterogeneous flame divergence equal to 120 that are experimentally measured taking into account the advective mass fluxes. Relations (21) and (22) admit the linear approximation of the dependence Ф = f(X,) associated with fuel burnout and temperature rise in the flame, as Figure 3 shows. The factor lf/d0 is the distance from the nozzle in calibers. This factor is applied here because the flame jet advances up to the plane of fuel ignition initiation according to the aeromechanic laws and without internal heat sources. The factor ΔTcom takes into account the effect of flame temperature increases due to the combustion reaction of fuel.
After the fuel ignition at l f = 0, particle size related burnout of fuel and temperature increase of the flame continuum take place in accordance with the integral curve Ф = f(X). By taking into account this factor, we make the task iterative, therefore the linear approximation Ф = f(X) is a possible solution, while dependences (20)–(22) describe the second-order influences. In fact, these approximations result from the expansion of the V( l f ) function, as shown below.
The amount of gas phase in the flame continuum Vg, m3/kg, changes from αV0 at the flame starting point with l f = 0 before the ignition (α − 1) to V0 + Vgα=1 at the flame end when the combustion is completed, where l f = 1.0. If the function Vg( l f ) is expanded in the Taylor series, the first two addends of the series give
V g ( l f ) = α V 0 + l f V g α = 1 l f .
The authors determine the variables V0 and Vgα=1 that form part of Formula (23) by the reference data for a specific type of solid fuel [8] at 0 ≤ l f ≤ 1.0. Similarly, the authors determine the flame continuum enthalpy in its temperature change function, and the temperature changes can be determined in the first approximation by the enthalpy change with cr = const.
The determined flame boundaries make it possible to plot the starting and final points of the flame continuum on the X-axis, and, accordingly, to find the starting and final points of the flame on the Φ-axis, considering it as the heat source and taking into account the initial flame section. The corrections to the flame continuum velocity (19) are applied sequentially: only the correction kF is applied in the range from lf = 0 to lf = li, after that all three corrections up to the final flame cross-sections are applied.
The length of the initial section preceding the ignition initiation li depends on the initial velocity of the flame flux w0 and the heat transfer rate in terms of convection α, W/m2K, which gives the following if only kF is taken into account
l i 3 + 19.0476 l i 2 + 90.703 r 0 3 l i = 0.22576 10 3 w 0 r 0 2 / α ,
or if li is substituted with li,f = li/d0 and dimensionless similarity numbers Re = w0d0 and Nu = αd0 are applied
l i , f 3 + 9.5238 l i , f 2 + 5.6689 d 0 l i , f = 0.04 10 6 Re / d 0 N u .
Equation (24) can be solved either by Cardano’s method, or numerically; adaptation calculations show that the length of the initial section is 0.95 ... 1.30 calibers, which is consistent with the experimental data.
The transformations of the abscissa, make it possible to numerically determine the function Φ on the graph ordinate as an integral function of the increase in the flame continuum enthalpy, which can be related to its adiabatic temperature
t a d = Ф ( l f ) Q T T c T V g .
If the available heat of the fuel flux Qpp and the yield of combustion products per the fuel mass unit Vg with the heat capacity of the combustion products cp are substituted in the Formula (26), the local adiabatic temperature value of the flame continuum tad can be determined when the obtained result is multiplied by the normal deviation function Ф( l f ) for the given coordinate l f [25,26,27]. The transition from Φ(U) to Ф( l f ) is realized by calculating the time of the particle size related burnout of fuel in accordance with (17), which corresponds to the time determined by the actual velocity of particles in the flame τ = lf/wf according to (19). The obtained dependence is solved with respect to the relation between lf and u = (δi − ξ)/(ξ·σ).

6. Validation of the Results

The following results were obtained: the dependence of the parameters of the heterogeneous flame on the mathematical description of the fractional composition of coal dust was substantiated; the characteristics of the combustion process that were previously not taken into account in the theory of heat transfer were revealed; theoretically confirmed new methodological approaches to the theory of heat transfer in the furnace of the boiler unit.
The results of theoretical studies were checked by the authors at operating thermal power plants, in particular, when burning coal dust in E-160 boilers. In the laboratory of thermal power engineering of the scientific and technical center under the Ministry of Energy of the Kyrgyz Republic, comparative analyzes [28] and experimental studies were conducted on samples of coal dust of various compositions. The movement of coal dust through the sampling tubes of devices, as well as the degree of grinding in roller and hammer mills was taken into account. The results showed good convergence within the engineering error of theoretical and experimental data.
In addition, the possibility of treating the torch as a continual medium was theoretically confirmed; the prospects for identifying the distribution of the adiabatic temperature along the length and height of the torch continuum in the zone of intense combustion of the combustion chamber of the boiler unit are shown.
On boilers E-210-140, experiments were carried out to measure the length of the torch, in particular, its initial part, as well as the height of the zone of intense combustion. In the course of the experiments, techniques approved during the operation and measurement of parameters at high-temperature installations were used, for example, the errors of the primary and secondary measuring instruments were multiplied. The discrepancy between the results of theoretical and experimental studies on the measurement of the temperature and length of the torch tended to 3%, which could be explained by some error in conducting experiments at high temperatures in boiler installations, for example, re-radiation and high dust content of the torch in the combustion space.
Table 2 shows our results for the length l f and the temperature of the torch at the length of the torch l f during the operation of the boiler E-210-140. Data is given depending on the steam output of the boiler unit. The nominal load of the boiler is D nom = 53 kg of steam per second. In addition, the data are given when burning Kazakhstan coal with a heat of combustion of 24,560 kJ/kg.

7. Conclusions

The possibility of considering the torch as a single continuum, to which the authors refer flue gases in the zone of intense combustion and screen tubes, is theoretically confirmed. The prospects of determining the type of distribution of the adiabatic temperature along the length and height of the flare continuum in the zone of intense combustion of the combustion chamber of the boiler unit are shown.
The authors applied a new mathematical method to describe the behavior of the polydisperse flux ensemble of fuel particles in a combustion process, which allowed us to determine the function of particle size related burnout of fuel and the related enthalpy and local adiabatic temperature of the flame continuum.
Besides, the authors determined the length li of the initial flame section, where the fuel does not ignite yet. It depends on the intensity of heat supply to the initial section α, the accepted temperature tint of intensive volatilization, the velocity of flux w0 at the burner exit and its change due to the temperature rise and flame divergence.
The burnout time of δi-size particles of the multiple-size particle flux depends on the fuel characteristics and the kinetic parameter ωk of the combustion process, the temperature Tf of the flame continuum being taken into account.

Author Contributions

Conceptualization, E.T. and K.O.; Methodology, K.O. and S.A.; Validation, E.T., K.O. and S.A.; Formal Analysis, K.O.; Investigation, K.O. and S.A.; Resources, E.T. and K.O.; Data Curation, K.O.; Writing-Original Draft Preparation, K.O. and S.A.; Writing-Review & Editing, K.O. and S.A.; Visualization, K.O.; Supervision, K.O.; Project Administration, S.A.; Funding Acquisition, E.T., S.A. and K.O. All of the authors contributed significantly to the completion of this manuscript, conceiving and designing the research, writing and improving the paper. All authors have read and approved the manuscript.

Funding

The work was supported by the Government of the Russian Federation (Resolution No. 211 of 16 March 2013), contract No. 02.A03.21.0011.

Acknowledgments

The authors thank the editor and reviewers for their insightful comments. The authors thank South Ural State University (SUSU) for supporting.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Noble, A.; Luttrell, G.H. A Review Of State-of-The-Art Processing Operations in Coal Preparation. Int. J. Min. Sci. Technol. 2015, 25, 511–521. Available online: http://www.elsevier.com/wps/find/journaldescription.cws_home/727915/description#description (accessed on 2 November 2018). [CrossRef]
  2. Dong, L.; Zhang, Y.; Zhao, Y.; Peng, L.; Zhou, E.; Cai, L.; Zhang, B.; Duan, C. Effect Of Active Pulsing Air Flow On Gas-Vibro Fluidized Bed For Fine Coal Separation. Adv. Powder Technol. 2016, 27, 2257–2264. Available online: http://www.elsevier.com (accessed on 1 November 2018). [CrossRef]
  3. Peng, L.; Wang, Z.; Ma, W.; Chen, X.; Zhao, Y.; Liu, C. Dynamic influence of screening coals on a vibrating screen. Fuel 2018, 216, 484–493. [Google Scholar] [CrossRef]
  4. Peng, L.-P.; Liu, C.-S.; Song, B.-C.; Wu, J.-D.; Wang, S. Improvement for Design of Beam Structures in Large Vibrating Screen Considering Bending and Random Vibration. J. Cent. South Univ. 2015, 22, 3380–3388. Available online: http://www.springerlink.com/content/2095-2899/ (accessed on 4 November 2018). [CrossRef]
  5. Roy, P.; Dias, G. Prospects for pyrolysis technologies in the bioenergy sector: A review. Renew. Sustain. Energy Rev. 2017, 77, 59–69. [Google Scholar] [CrossRef]
  6. Han, J.; Elgowainy, A.; Dunn, J.B.; Wang, M.Q. Life Cycle Analysis of Fuel Production from Fast Pyrolysis of Biomass. Bioresour. Technol. 2013, 133, 421–428. Available online: www.elsevier.com/locate/biortech (accessed on 4 November 2018). [CrossRef] [PubMed]
  7. Moreno-Piraján, J.C.; Bastidas-Barranco, M.J.; Giraldo, L. Preparation of activated carbons for storage of methane and its study by adsorption calorimetry. J. Therm. Anal. Calorim. 2018, 131, 259–271. [Google Scholar] [CrossRef]
  8. Asotani, T.; Yamashita, T.; Tominaga, H.; Uesugi, Y.; Itaya, Y.; Mori, S. Prediction of ignition behavior in a tangentially fired pulverized coal boiler using CFD. Fuel 2008, 87, 482–490. [Google Scholar] [CrossRef]
  9. Mandø, M.; Rosendahl, L.; Yin, C.; Sørensen, H. Pulverized straw combustion in a low-NOx multifuel burner: Modeling the transition from coal to straw. Fuel 2010, 89, 3051–3062. [Google Scholar] [CrossRef]
  10. Inyang, H.I.; Bae, S.; Pando, M.A. Contaminant dust suppression materials: A cost-effectiveness estimation methodology. Meas. J. Int. Meas. Confed. 2016, 93, 563–571. [Google Scholar] [CrossRef]
  11. Yu, X.; Luo, Z.; Li, H.; Gan, D. Effect of vibration on the separation efficiency of oil shale in a compound dry separator. Fuel 2018, 214, 242–253. [Google Scholar] [CrossRef]
  12. Leimkulher, B.; Reich, S. Simulating Hamiltonian Dynamics; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  13. Hairer, E.; Lubich, C.; Wanner, G. Geometric Numerical Integration, 2nd ed.; Springer: Berlin, Germnay, 2006. [Google Scholar]
  14. Brugnano, L.; Iavernaro, F. Line Integral Methods for Conservative Problems; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
  15. Blanes, S.; Casas, F. A Concise Introduction to Geometric Numerical Integration; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
  16. Sanz-Serna, J.M. Symplectic Runge-Kutta Schemes For Adjoint Equations, Automatic Differentiation, Optimal Control, and More. SIAM Rev. 2016, 58, 3–33. Available online: http://epubs.siam.org/doi/pdf/10.1137/151002769 (accessed on 2 November 2018). [CrossRef]
  17. Hairer, E.; Wanner, G. Solving Ordinary Differential Equations II. Stiff and Differential-Algebraic Problems, 2nd ed.; Springer Series in Computational Mathematics; Springer: Berlin, Germany, 1996; Volume 14. [Google Scholar]
  18. Brugnano, L.; Frasca Caccia, G.; Iavernaro, F. Efficient implementation of Gauss collocation and Hamiltonian boundary value methods. Numer. Algorithms 2014, 65, 633–650. [Google Scholar] [CrossRef] [Green Version]
  19. Silva, I.R.; Maboudou-Tchao, E.M.; de Figueiredo, W.L. Frequentist–Bayesian Monte Carlo test for mean vectors in high dimension. J. Comput. Appl. Math. 2018, 333, 51–64. [Google Scholar] [CrossRef]
  20. Maboudou-Tchao, E.M.; Silva, I. Tests for mean vectors in high dimension. Stat. Anal. Data Min. 2013, 6, 578–598. [Google Scholar] [CrossRef]
  21. Robbins, H. Asymptotically subminimax solutions of compound decision problems. In Proceedings of the Second Berkeley Symposium on Mathematical Statistics and Probability; University of California Press: Berkeley, CA, USA, 1951; pp. 131–148. [Google Scholar]
  22. Robbins, H. An empirical Bayes approach to statistics. In Proceedings of the Third Berkeley Symposium on Mathematical Statistics and Probability; University of California Press: Berkeley, CA, USA, 1955; pp. 157–164. [Google Scholar]
  23. Spitzner, D.J. An Asymptotic Viewpoint on High-Dimensional Bayesian Testing. Bayesian Anal. 2008, 3, 121–160. Available online: http://ba.stat.cmu.edu/journal/2008/vol03/issue01/spitzner.pdf (accessed on 2 November 2018). [CrossRef]
  24. Kyprianidis, K.G.; Skvaril, J. Developments in Combustion Technology; InTech: Cranfield, UK, 2016. [Google Scholar] [CrossRef]
  25. Bartholomew, A.; Fenn, R. Quaternionic invariants of virtual knots and links. J. Knot Theory Ramif. 2008, 17, 231–251. [Google Scholar] [CrossRef]
  26. Carter, J.S. Classifying immersed curves. Proc. Am. Math. Soc. 1991, 111, 281–287. [Google Scholar] [CrossRef]
  27. Rodríguez-Nieto, J.G. Characterization of signed Gauss paragraphs and skew-symmetric graded matrices. J. Knot Theory Ramif. 2018, 27, 1850002. [Google Scholar] [CrossRef]
  28. Aliukov, S.V.; Alyukov, A.S. Analysis of Methods for Solution of Differential Equations of Motion of Inertial Continuously Variable Transmissions. Available online: https://doi.org/10.4271/2017-01-1105 (accessed on 2 November 2018).
Figure 1. Screen for screening coal dust.
Figure 1. Screen for screening coal dust.
Energies 12 00136 g001
Figure 2. Rotating coal (biofuel) vibrating screen.
Figure 2. Rotating coal (biofuel) vibrating screen.
Energies 12 00136 g002
Figure 3. The Gaussian distribution density of the normalized random variable y = φ(u).
Figure 3. The Gaussian distribution density of the normalized random variable y = φ(u).
Energies 12 00136 g003
Figure 4. Hodograph of the Gaussian distribution density of the normalized random variable y = φ(u).
Figure 4. Hodograph of the Gaussian distribution density of the normalized random variable y = φ(u).
Energies 12 00136 g004
Figure 5. Normal distribution function Ф of the random variable x = u = (δξ)/(ξ·σ) in the case of sieve analysis the sieve analysis of pulverized coal with fractional composition R90.
Figure 5. Normal distribution function Ф of the random variable x = u = (δξ)/(ξ·σ) in the case of sieve analysis the sieve analysis of pulverized coal with fractional composition R90.
Energies 12 00136 g005
Table 1. The desired speed of the flare continuum.
Table 1. The desired speed of the flare continuum.
BoilerE-210-140E-160
wl,T, m/s7.9649.165
Table 2. The length l f and the temperature of the torch at the length of the torch l f .
Table 2. The length l f and the temperature of the torch at the length of the torch l f .
Load of Boiler0.75·Dnom0.85·Dnom0.95·Dnom1.00·Dnom
l f , m2.323.163.914.9
T, K1648167316931721

Share and Cite

MDPI and ACS Style

Toropov, E.; Osintsev, K.; Aliukov, S. New Theoretical and Methodological Approaches to the Study of Heat Transfer in Coal Dust Combustion. Energies 2019, 12, 136. https://doi.org/10.3390/en12010136

AMA Style

Toropov E, Osintsev K, Aliukov S. New Theoretical and Methodological Approaches to the Study of Heat Transfer in Coal Dust Combustion. Energies. 2019; 12(1):136. https://doi.org/10.3390/en12010136

Chicago/Turabian Style

Toropov, Evgeniy, Konstantin Osintsev, and Sergei Aliukov. 2019. "New Theoretical and Methodological Approaches to the Study of Heat Transfer in Coal Dust Combustion" Energies 12, no. 1: 136. https://doi.org/10.3390/en12010136

APA Style

Toropov, E., Osintsev, K., & Aliukov, S. (2019). New Theoretical and Methodological Approaches to the Study of Heat Transfer in Coal Dust Combustion. Energies, 12(1), 136. https://doi.org/10.3390/en12010136

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop