Simplified model of the heat exchange process in rotary regenerative air pre-heater Текст научной статьи по специальности «Физика»

URAL MATHEMATICAL JOURNAL, Vol. 2, No. 2, 2016

SIMPLIFIED MODEL OF THE HEAT EXCHANGE PROCESS IN ROTARY REGENERATIVE AIR PRE-HEATER1

Abdulla A. Azamov

Institute of Mathematics, National University of Uzbekistan, Tashkent, Uzbekistan, abdulla.azamov@gmail.com

Mansur A. Bekimov

Institute of Mathematics, National University of Uzbekistan, Tashkent, Uzbekistan, mansu@mail.ru

Abstract: A simplified mathematical model of a rotary regenerative air pre-heater (RRAP) is suggested and studied based on the averaged dynamics of the heat exchange process between nozzles and a heat carrier (i.e. air or gas-smoke mixture). Averaging in both spatial coordinates and time gives a linear discrete system that allows deriving explicit formulas for determining the characteristics of the air heater and establishing some properties such as periodicity, stability, ergodicity and others.

Key words: Heat exchange, Cyclic process, Averaging, Linear discrete system, Stability, Ergodicity, Inverse problem.

Introduction

A rotary regenerative air pre-heater (RRAP) is a special unit, usually applied in thermal power plants (TPP) in order to increase its efficiency by heating the air which is blowing into a boiler of the plant by means of an exhaust of hot mixture of smoke and gas (from now on simply gas) generated in fuel combustion. The RRAP essentially reduces the thermal pollution of atmosphere [6], [8-12].

Currently, several types of RRAPs are used in TPP. In the present paper, we consider the case of a unit with the main part consisting of a rotating cylindrical drum with metal nozzles of high thermal conductivity. The region of space occupied by the wheel of the RRAP is divided into two parts Ba and BG by a fixed conditional plane passing through the axis of the cylinder. During the work of the RRAP, the atmospheric air passes through the part Ba in one direction parallel to the axis of the drum, being heated by the nozzles; as a result, the temperature of the nozzles is reduced. Through the part BG, the gas flows in the opposite direction, being cooled by the heat output to the nozzles. The final transfer of heat from the hot gas to the cool air occurs due to the rotation of the drum around its axis.

Monitoring and control of the temperatures of air and gas leaving the RRAP and especially the temperature of the nozzles is the important problem for the effective exploitation of the RRAP [3], [15], [17-19]. The direct measurement of the temperatures of air and gas of both inlet and outlet is easily carried out; in contrast, the observation and control of the nozzle temperature require to use sophisticated measuring tools and represent an important problem. Therefore, the treatment for mathematical modeling of the heat exchange process in the RRAP is reasonable.

1This work was supported by Committee for Coordination Science and Technology Development Under Cabinet of Ministers of Uzbekistan (project no F4-FA-F014).

In recent decades, a variety of mathematical models of the heat exchange process in RRAPs were suggested [1-3], [5] and see also [8]. It should be noted that the thermodynamic process in RRAP nozzles can be in principle described by the classical heat conduction equation, but, without additional simplifying assumptions, it is difficult to write out the boundary conditions since the nozzle system has a rather complex geometry. On the other hand, modeling of the thermodynamic process of air heating and gas cooling passing through the drum of the RRAP is even more challenging problem pertaining to thermos-aerodynamics [8]. Furthermore, the necessity of considering the turbulence of flow that arises when air and gas pass through the RRAP drum also adds complexity. Because of these and other features, all mathematical models of the RRAP are built under essential simplifying assumptions.

In the present paper, we propose a mathematical model of the thermodynamic process in the RRAP based on averaging the quantities associated with the heat exchange process between the nozzles, air and gas in both spatial coordinates and time interval. As a result, we obtain rather simple linear discrete equation. This allows us to write out explicit calculating formulas for the current values of parameters, and find steady and periodical states, establish the ergodicity and other properties. Then, we provide appropriate calculations for the case when averaging over the time variable is carried out on a period of time equal to a half-cycle of rotation of the RRAP drum. In subsequent parts of the work, models with averaging performed over a small time interval are considered, a comparative analysis of numerical results and experimental data is fulfilled.

1. Derivation of equations

Let the RRAP drum be of the form of a cylinder x2 + y2 < R2, 0 < z < H. Assume that the parts Ba and BG are described by conditions y > 0 and y < 0, respectively. Let 6(t, x, y, z) denote the temperature at the point (x, y, z) of the drum occupied by the nozzles at time t, t > 0, and let T(t, x, y, z) denote the temperature of heat transfer (air or gas) at the point (x, y, z) of the drum outside the nozzles. The pair of quantities 6(t, x, y, z) and T(t, x, y, z) completely characterizes the heat change process in the RRAP drum. However, as mentioned above, due to the complexity of the configuration of the nozzles, the initial-boundary problem for corresponding system of equations of thermo-aerodynamics is too hard to be investigated by analytical methods. One of the ways to overcome such complexity is the method of averaging [4], [7]. To this end, we consider the corresponding average values

®A(t) = 1 I QA(t,x,y,z)dv, QG(t) = „A.x [ QG(t,x,y,z)dv,

V (BA) JB'A V (Bg)JB'G

Ta^ = T7Tw\I TA{t,x,y,z)dv, TG(t) = —^—[ TG{t,x,y,z)dv,

V (BA) J BA V (BG) J B'G

where BA = Ba n B, BG = Ba n B, B is the part of the drum occupied by the nozzles, BA = Ba\BA, BG = BG\BG, dv is the volume element, V denotes the volume of the corresponding part of the drum.

Next, we perform averaging over the time intervals I(n) = [nh, (n + 1)h) as well, where h is the half-turn time of the RRAP drum, n = 0,1,2,.... Denote the average temperatures obtained in this way by

Xn = ]- [ &A(t)dt, yn = i [ eG(t)dt,

h J I(n) h JI (n)

Un = \ [ TA(t)dt, vn = i [ TG(t)dt.

h J I(n) h JI(n)

Equations connecting these quantities are derived under the following simplifying assumptions:

1. for nh < t < (n + 1)h, the RRAP drum is motionless, the portion of air (respectively, gas) filling the part BA (respectively, EG) is also motionless, the heat exchange between the heat carrier and nozzles occurs in accordance with the linear Newton law [8];

2. at the time t = (n + 1)h, the heated portion of air with temperature un leaves the area BA, the gas portion cooled down to temperature vn leaves the region EG, the drum turns to 180o with jump (i.e., the part BA goes to the part EG and vice versa); then the part BA is filled in by a new portion of air from outward (or from calorifer in case of its connection to the RRAP, [12]) of temperature pn, and the part EG is filled in by a new portion of gas of temperature qn.

In accordance with the Newton law, we have the following relations:

Xn+1 = (1 - i3h)y,n + fihqn, yn+1 = (1 - ah)x,n + ahpn, (1.1)

Un = Pn + Yh(xn - Pn), Vn = qn + Sh(y,n - q,n), (1.2)

where a, ft, y, S are parameters depending on characteristics of the heat exchange process in the RRAP (the geometry, the heat capacity of the carcass of drum, the system of nozzles and their structure, the humidity of air and gas and their thermodynamic characteristics, the coefficients of heat conductivity and diffusion, parameters that characterize the heat exchange process on the contact surface of the nozzles with air and gas and etc.).

Relations (1.1) and (1.2) have been obtained as a result of extremely simplifying assumptions on the heat exchange process in the RRAP. Nevertheless, due to such simplification, system (1.1) allows a fairly complete analysis; therefore, it can serve as a basic model to describe the work of the RRAP.

2. Solution of the system

Introducing the vectors zn = (xn, yn)T and rn = h(ftqn, apn)T, where T is the transpose sign

to transform the row-vector to the column-vector, and the matrix

A = ( 0 1

- ah 0 Y

one can rewrite system (1.1) as follows:

Zn+1 = Azn + rn, n = 0, 1, 2, ... . (2.1)

All further reasoning is conducted under the assumption 0 < ah, fth < 1, called the physical realizability of the model. This condition implies the fact that the eigenvalues of the matrix A, which are where ^ = belong to the interval (-1, 1) and therefore all the

solutions of system (1.1) are asymptotically stable [13,14], [16].

The solution of equation (2.1) can be written by the Cauchy formula [13]

n—1

m„ i \ ^ An— 1—k

AnZ0 + An—1—krk. (2.2)

zn

k=0

Since An = /j,nE for even numbers n, and An = ¡in—1A for odd numbers n, then equation (2.2) can be transformed to the form

n/2—1 n/2—1

E jn—2—2j r2j + z j=0 j=0

Zn = jnZ0 + A Z jn—2—2j r2j + Z jn—2—2j r2j+1,

for even n, and to the form

(n-1)/2 (n-1)/2

zn = ^n-1Azo + A Y, ^n-1-2jr2j-i + E ^n-1-2jr2j j=1 j=0

for odd n.

Based on these relations, one can easily obtain explicit formulas for xn and yn. In practice, however, it is more convenient to calculate them, in case of need, directly from equations (1.1) that are accommodated to computer calculations. In the rest of the paper, some properties of solution of systems (1.1), (1.2) are established.

3. Steady state and periodic regimes

Let us first consider the case of a steady state, which occurs under the assumption of constant incoming flows. Such a state can be established in the periods of time measured in hours, if the parameters of the air and the energy load at the TPP remain virtually unchanged. Since h is measured in minutes, the work of the RRAP actually consists of long periods of steady states and relatively short transition intervals from one steady state to another. Therefore, it is important to determine the parameters of the RRAP under steady state conditions.

Thus, let

pn = p, qn = q n = 0,1,2,..., r = (Pf , ap)T. Then formula (2.2) can be simplified even more: for even n

Zn = ^nZo + (1 - ^n)(E - A)-1r , (3.1)

or in the component form

1 — un

xn = txnx0 + ——-—[p + q( 1 - f3h)],

a + p — aph

1 — an

Vn = l~inyo + ——5-—[q+p(l - ah)],

a + p — aph

and for odd n

zn = nn-lAz0 +---—(E - nnA)(E -

a + p — aph

or

= lT-\ 1 " Ph)yo + b(l - /Xra+2) + 9(1 - ph)( 1 - z/1"1)],

a + p — aph

Vn = ~ ah)x0 + 1 [q{ 1 - /xra+2) +p{ 1 - ah){ 1 -

a + p — aph

For the steady state, we obtain z = (E — A)-1f from the equation z = Az + r (the invertibility of the matrix E — A follows from the condition of physical realizability) or in the component form

1 (ap + f3q(l- f3h)),

a + p — aph

\ (3.2)

y=a + P-aPh Ofr + °P(1 ~ "*)) •

Due to the fact noticed above, steady state (3.2) is asymptotically stable. On time intervals measured in weeks, a periodic state for the RRAP can be formed.

Theorem 1. Let the sequence rn be periodic with the period m, m > 2, i.e. rn+m = rn, where n = n0, n0 + 1, n0 + 2, ... ; n0 is some fixed value of n. Then system (2.1) has a unique periodic solution of period m, which is asymptotically stable.

Proof. The uniqueness and asymptotic stability follows from the condition of physical real-izability. To establish the existence, we set

m— 1

z"n = (E — Am) — Am—1—krk (3.3)

k=0

(it follows from the condition of physical realizability that the matrix E — Am is invertible). Let zn be a trajectory with the initial point z0. Rewrite relation (3.3) as follows:

m— 1

* _ Am * I V"^ Am—1—k _ z* zn _ A z0 + / A rk _ zm-

k=0

Then,

m—1 n+m-1

zn+m _ An+mz0 ^ An+m—1—krk + Z An+m—1—krk. (3.4)

k=0 k=m

The sum of the first two summands in the right hand side of (3.4) is Anzm, the third one is equal to

n— 1 n—1

Y^An—1—krk+m An—1—k rk k=0 k=0

since rk is periodic. Thus,

Y^An—1—k rk = zn — Anz0-

k=0

Now, in accordance with (3.4), we obtain

* _ \n * I An * _

zn+m _ A zm + zn — A z0 _ zn,

and the proof is complete. □

4. Boundedness and ergodicity of the solutions

In time intervals of longer duration, in a changeable external environment and energy load at the TPP, the sequence rn is not periodic, a fortiori, it is not stationary. In this regard, we give two properties of the solution for more general classes of systems (2.1), which models more or less irregular regime of the heat exchange process.

Theorem 2. If rn is a bounded sequence, then each solution of equation (2.1) is also bounded.

This statement is also a special case of more general theorem [11].

Theorem 3. Let lim rn = l. Then each solution zn approaches the limit (E — A)-1l as n ^ to

n—^^o

independently of z0.

Proof. First assume that l = 0. Let M = max |rn|. Given any e > 0, choose n0(e) such that

n

|rn| < -(1 — ||A||)£ at n > no(e). Then choose N(no(e)) such that

ßn <

1 -

2M (1 - ||Ano I

for all n > N(no(e)). Clearly, N(n0(e)) can be chosen to satisfy N(n0(e)) > no(e). Then, for n > N (n0(e)), we have

|zn| < ||An|| jzöj +

no —1

£A

k=0

n— 1—k

rk

+

n— 1

^ ^ An—1—k

k=no

rk

Since 0 < ah, ßh < 1, the first summand tends to zero as n ^ to. Denote two other terms by S, and Sn1. Then, in view of ||A|| = ß € (0, 1), we get

Si < \\An—no

no —1

no 1

k=0 n— 1

Ano

— 1—k

no — 1

|rk | <M \\An—no \ Ê 11 Ano 1 k k=0

< Mß

1 -|An

e

1-IIAII <2'

SU ^ |An—1—k |rk| <

k=no

1-

\ Ak

k=0

1-

1-

e

<2-

Hence, |zn| < S^ + SI1 < e at n > N(n0(e)), i.e. zn 0.

Let now rn — l as n — to, where l = 0. We make the change of variables

Then

zn — zn + (E A) l rn — rn l

Zn+1 + (E - A)—1l — A[zn + (E - A)—1l] + rn + l — Azn +

□

where rn — 0. As proved, zn — 0. Therefore, zn — (E — A)

In general, the change of the parameters characterizing the state of air, as well as the variation of the load at the TPP is random with a hard-determinable distribution function.

Considering the work of the RRAP as a stochastic process, we leave it for the next part of the paper, and we now present another property, taking into account irregular characters of values of phase variables of system (1.1), (1.2).

A sequence an is called almost-periodic, if it can be represented in the form bn + cn, where bn is periodic and cn — 0 as n — to.

Corollary 1. If rn is almost-periodic, then each solution of (2.2) is also almost-periodic.

Definition. A sequence xn is called ergodic, if the sequence of Cesaro means

a„

X1 + X2 +----+ x,n

n

converges as n — to.

Theorem 4. If the sequence rn is ergodic, namely

r0 + r1 +-----+ rn-1 ,

Pn =--> I,

n

then each solution

Zn is also ergodic with

Sn = *!+** + •••+*» {E - A)-H.

n

e

1

2

2

Proof. Express SN as the sum of two terms: SN = £/ + £//, where

N N n— 1

Then

as n —^ to.

For £//, we have

N^ 0 JJ N

n=1 n=1k=0

1 N 1 Ej = - E Anz0 = -A{E - A)~\E - An)Z0 0

n=1

N—1 N N—1 /N—1—k \

= E E

k=0 n=k+1 k=0 \ n=0 /

N

—1

'N—1 N—1

£ rk — E AN —1—k rk k=0 k=0

= (E — A)'1 Pn - ^^ £ A^"1" V

N k=0

By hypothesis of the theorem, (E — A)—1 pn — (E — A)—1l as N — to. Assuming pk = 6 + l, where 6 — 0, we arrive at the equation

rk = (k + 1)Pk+1 — kpk = (k + 1)6+1 — k6 + l.

Thus,

N ^ k N

k=0

'N — 1 N — 1

E AN—1—k [(k + 1)6+1 — k&] — E AN—1—kI k=0 k=0

Obviously,

N — 1

iN—1—k;

AN—1—kl — 0

— k=0

as N — to. Let

N—1

5 = jj E

N k=0

Show that S — 0 as N — to. Indeed,

N—1 , N—1

\s\ < e I6i < E i&i • k=0 N k=0

By Theorem 3, the right-hand side of this inequality tends to 0 as N — to. Thus, finally we obtain

1N

n=1

which is our claim. □

5. Finding the coefficients of the system

As mentioned above, besides the temperature and velocity of incoming air and gas, the drum rotation speed, the heat exchange process in the RRAP depends on many parameters expressing the thermodynamic characteristics of air and gas, the material and geometry of the nozzles of the RRAP, thermal properties of the drum casing etc.

There are a number of works devoted to the mathematical modeling of the RRAP where formulas are given to calculate the values of the parameters above on the basis of molecular physics laws [1-3], [8], [9]. In addition, finding the numerical parameters characterizing the RRAP is only possible with a certain accuracy. From this point of view, it is much easier to define them with a satisfactory accuracy by solving the inverse problem for system (1.1), (1.2) on the basis of empirical data. With regard to model (1.1), (1.2), such a problem consists in calculating the values of a, ft, y, and 5 based on the results of measurements of temperatures of outgoing air un and gas vn. Here, there is a wide field of application of the least squares method and tools of mathematical statistics. Here, we confine ourselves to the simplest case, when a, ft, y, and 5 are found by measuring ui, u2, vi, and v2 assuming that incoming streams pn = p, qn = q are stationary. We can assume that uo = p, v0 = q, and the values of u1, u2, v1, and v2 are found by the direct measurements. As a result, we arrive at the problem of finding the unknowns a, ft, y, 5 based on the given values h, p, q, u1, u2, v1, v2, without measuring the values xn, yn (the nozzle temperature).

We have

x2 = (1 — fth)y1 + fthq, y2 = (1 — ah)x1 + ahp, (5.1)

u1 = p + Yh(x1 — p), v1 = q + 5h(y1 — q), (5.2)

u2 = p + Yh(x2 — p), v2 = q + 5h(y2 — q). (5.3)

Set uk = uk — p, vk = vk — q,k = 1,2, (these values have the clear physical meaning). Substituting the values of x2 and y2 from (5.1) into (5.3), we obtain the following system

Yh(x1 — p) = u1, 5h(y1 — q) = v 1, (5.4)

Yh[(1 — fth)y1 + fthq — p] = u2, 5h[(1 — ah)x1 + ahp — q] = v2, (5.5)

with 6 unknowns a, ft, y, 5, x1, y1. It is nonlinear and, in general, cannot be solved explicitly.

Therefore, we use the fact that, for the values of a, ft, y, 5, there is a priori estimate 0.1 ^ 0.6 and h < 1. This allows us to neglect the terms containing ftYh2 and a5h2.

Then equations (5.5) take the form Yh(y1 — p) = u2, 5h(x1 — q) = v2. As a result, for the intermediate unknowns x1 , y1 , we obtain the linear system

xi -p _ ui yi-g _v2 Vi-p Û2 xi-q vi

(5.6)

with the determinant equal to UiV i — u2V2. We call the quantity x = UV i — u2V2| the divergence coefficient of the RRAP. The deviation of x from zero is a characteristic of the RRAP that expresses how the rates of air heating and gas cooling differ. Further, we assume x = 0. It follows from (5.6) that

(Ui — U2)V2 p + (vi — V2)uiq (Ui — U2)Vip + (Vi — V2)u2q

X\ =-—-—-, y 1 =-—-—-.

UiVi — U2V2 UiVi — U2V2

Substituting the values xi, yi into (5.4), (5.5), we obtain the final formulas

xi — q_ yi — ^ yi — p_ xi — p_

ah =-vi--v2, fth =-u 1--u2

xi — p xi — p yi — q yi — q

jh = ———, =

xi - p yi - q

Table 1 shows posteriori values of parameters a, P, 7, 5 calculated by these formulas under the assumption that p = 32°, q = 282°, h = 0.25 [4, 5].

U\ Vi U2 Vi a ß 7 5

3.61 -2.58 6.75 -4.82 0.048 0.068 0.041 0.029

Table 1. Posteriori values of parameters a, ß, 7,

6. Conclusion

In the present paper, we have proposed the mathematical model of thermodynamic process of the RRAP, which is described by linear discrete equations. To obtain this, we have used averaging the quantities associated with the heat exchange process between the nozzles, air and gas in both the spatial coordinates and time interval. We have found steady and periodical states, established the ergodicity and other properties. Next, we have studied the cases when the time averaging is performed over the period of time equal to the half-cycle of rotation of the RRAP drum as well as when the time averaging is performed over a small time interval. Finally, we have provided the comparative analysis of the numerical results obtained and the experimental data.

Acknowledgements

Authors express their gratitude to G.I. Ibragimov for useful discussion and help.

REFERENCES

1. Alagic S., Kovacevic A. and Buljubasic I. A numerical analysis of heat transfer and fluid flow in rotary regenerative air pre-heaters // Strojniski Vestnik, 2005. Vol. 51, no. 7-8, P. 411-417.

2. Armin Heidari-Kaydan, Ebrahim Hajidavalloo. Three-dimensional simulation of rotary air pre-heater in steam power plant // Applied Thermal Engineering, 2014. Vol. 73. P. 397-405.

3. BoStjan Drobnic, Janez Oman, Matija Tuma. A numerical model for the analysis of heat transfer and leakages in a rotary air preheater //Int. J. of Heat and Mass Transfer, 2006. Vol. 49. P. 5001-5009.

4. Burd V.Sh. Method of averaging for differential equations on an infinite interval: Theory and Applications. Chapman and Hall/CRC, 2007. 360 p.

5. Chi-Liang Lee. Regenerative air preheaters with four channels in a power plant system // J. of Chinese Institute of Engineers. 2009. Vol. 32, no. 5, P. 703-710.

6. Grebennikov E.A. Method of averaging in applied problems. Moscow: Nauka, 1986. 256 p. [in Russian]

7. Kirsanov Yu.A. Cyclic thermal processes and the theory of thermal conductivity in regenerative air heaters. Moscow: Fizmatlit, 2007. 240 p. [in Russian]

8. Kovalevskii V.P. Simulation of heat and aerodynamic processes in regenerators of continuous and periodic operation. I. Nonlinear mathematical model and numerical algorithm // J. of Engineering Physics and Thermophysics, 2004. Vol. 77, no. 6, P. 1096-1109.

9. Kudinov A.A. The study of heat transfer processes in rotary regenerative air power boilers // J. of Energetics, 2012. No 6. P. 32-34. [in Russian]

10. Kudinov A.A., Ziganshina S.K. Energy savings in power and heat technologies. Moscow: Mashinos-troenie, 2011. 374 p. [in Russian]

11. Liu Fuguo, Zhou Xingang. Heat transfer model of tri-section rotary air preheater and experimental verification // J. of Mechanical Engineering, 2010. Vol. 46, no 22, P. 144-150.

12. Ismatkhodjayev S.K., Nuraliyev E.N., Abduraimov F.E., Mokrushev V.A., Azamov A.A., Bekimov M.A. A method for automatic adjusting the air temperature in air heater. Patent No IAP 05072. [in Russian]

13. Romanko V.K. The course of difference equations. oscow: Fizmatlit, 2012. 200 p. [in Russian]

14. Saber Elaydi. An introduction to difference equations. Third edition, Springer Science+Business Media, New York, 2005.

15. Shruti G., Ravinarayan Bhat, Gangadhar Sheri. Performance evaluation and optimization of air preheater in thermal power plant // International J. of Mechanical Engineering and Technology (IJMET), 2014. Vol. 5, iss. 9. P. 22-30.

16. Srochko V.A. Numerical methods: Lectures. Irkutsk: Irkutsk University Press, 2004. 205 p. [in Russian]

17. Wang H.Y. Analysis on thermal stress deformation of rotary air-preheater in a thermal power plant // Korean J. Chem. Eng., 2009. Vol. 26, P. 833-839.

18. Wang H.Y., Zaho L.L., Xu Z.G., Chun W.G., Kim H.T. The study on heat transfer model of tri-sectional rotary air preheater based on the semi-analytical method// Appl. Therm. Eng. 2008. Vol. 28, no. 14-15. P. 1882-1888.

19. Zhuang Wu, Roderick V.N. Melnik, Finn Borup. Model-based analysis and simulation of regenerative heat wheel // Energy and Buildings, 2006. No. 28. P. 502-514.