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 (NO
x) 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(−
bδ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
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.
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
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
ycalc ≡
y2, was for the description of the practical particle size distribution of fuel dust
ypr ≡
y1. Thus, it is rational to subject the dependence (3) to the logarithmic linearization
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.56
d0/0.15
l 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.73
w02d0/l, whence it follows that
wl = 1.93
w0(
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 = ρw2/ρ0w02, 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
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 = 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
thermal conductivity, density, and mass heat capacity of the particle material in regulated units. Here
Bi =
αλ/
δ,
Fo =
aτ/
δ2, where
α is heat transfer coefficient, W/(m
2K),
λ is thermal conductivity, W/(m K),
a is thermal diffusivity, m
2/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
ςTo a spherical fuel particle approximation
ς = 3 and
m = 1 + 0.1
Bi 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/m
2·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)
where
θ0 = tg −
t0 is the difference in the temperatures between the flue gases
tg and the initial temperature of the particle
t0; θint = tg −
tint 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
UMoreover, the normal distribution function Φ(
U) is associated with the probability integral Φ(
U) and the error function
erf(
U/20.5) by the dependence
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 R
90.
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τ = −2
McGc/ρ, where
Mc = 12 kg/kmol is the molar mass of carbon,
ρ ≈ 800 kg/m
3 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 CO
2 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:
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.21
Pg/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 1/αD parallel to the lumped kinetic resistance 1/k1, we obtain the total diffusive-kinetic resistance to the combustion process RΣ = 1/αD + 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 αD ≫ k1, when the diffusion processes do not limit the burning rate, ωk = 1.0; for αD ≪ k1, 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/m
3, 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
(13) takes into account that the stoichiometric air consumption
V0, m
3/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]
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
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, m
2/s,
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
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.
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
= 1.0
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
= 0 to
= 1.0 during its motion and combustion process can be represented by the dependence
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:
The relation (20) is obtained geometrically on the basis of central angles of the heterogeneous flame divergence equal to 12
0 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 = 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() function, as shown below.
The amount of gas phase in the flame continuum
Vg, m
3/kg, changes from
αV0 at the flame starting point with
= 0 before the ignition (α − 1) to
V0 +
Vgα=1 at the flame end when the combustion is completed, where
= 1.0. If the function V
g(
) is expanded in the Taylor series, the first two addends of the series give
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 ≤
≤ 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/m
2K, which gives the following if only
kF is taken into account
or if
li is substituted with
li,f =
li/d0 and dimensionless similarity numbers
Re = w0d0/ν and
Nu = αd0/λ are applied
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
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 Ф(
) for the given coordinate
[
25,
26,
27]. The transition from Φ(U) to Ф(
) 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
and the temperature of the torch at the length of the torch
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
= 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.