Strojniški vestnik - Journal of Mechanical Engineering 53(2007)3, 158-164 UDK - UDC 536.27 Izvirni znanstveni članek - Original scientific paper (1.01) Energijska in eksergijska analiza sotočnih in protitočnih prenosnikov toplote z uporabo merilnih podatkov Energy and exergy analysis of a parallel and counter-flow heat exchangers using measured data Antun Galovič1 - Marija Živič2- Ahmet Can3 ('University of Zagreb, Croatia; 2JJ. Strossmayer University of Osijek, Croatia; 3Trakya University, Turkey) Prispevek podaja energijsko in eksergijsko analizo sotočnih in protitočnih ločilnih prenosnikov toplote z uporabo marilnih podatkov. Sestavili smo merilno progo, na kateri smo merili vstopne in izstopne temperature in masni tok vode. Za določitev brezrazsežnega razmerja izgubljene eksergije in prenesenega toplotnega toka v odvisnosti od brezrazsežnih parametrov prenosnika toplote: razmerja absolutnih temperatur na vstopu nT razmerja toplotnih moči n3 in števila enot prenosnika toplote n2, smo razvili analitični model. Za vse primere smo izračunali tudi učinkovitost prenosa toplote. Rezultate smo prikazali v ustreznih brezrazsežnih diagramih. © 2007 Strojniški vestnik. Vse pravice pridržane. (Ključne besede: prenosniki toplote, energijske analize, eksergijske analize, izmerjene vrednosti) This paper presents an energy and exergy analysis of a parallel and counter-flow recuperative heat exchangers using experimental data. An experimental rig was constructed to measure the inlet and outlet temperatures and the mass flow rates of streams. The analytical model was developed to obtain a non-dimensional relationship between the destroyed exergy and exchanged heat-flow rate as a function of the non-dimensional parameters of a heat exchanger: the ratio of inlet absolute temperatures, nT the ratio oj the heat-capacity rates, nv and the number of heat-transfer units, nT The effectiveness of the heat exchange is also calculated for each case. The results are shown in appropriate non-dimensional diagrams. © 2007 Journal of Mechanical Engineering. All rights reserved. (Keywords: heat exchangers, energy analysis, exergy analysis, measured values) 0 INTRODUCTION Energy-exergy analyses of heat exchangers have been the subject of much research over the past few decades, [1] to [7]. These analyses are based on the first and the second laws of thermodynamics. From such analyses the parameters for the improved operation of a heat exchanger can be obtained. The destroyed exergy or the lost available work of a heat exchanger is due to two factors: the transfer of heat across the stream-to-stream temperature difference and the frictional pressure drop that accompanies the circulation of fluid through the apparatus. Taking into account these two parameters it is possible to optimize a heat exchanger on an entropy-generation minimization or on a minimum destruction of exergy ([8] and [9]). In this work the exergy destruction due to the pressure drop is neglected, because from experimental data it was clear that the pressure drop was very small. 1 DESCRIPTION OF THE MEASURING RIG Figure 1. shows the experimental setup of the heat exchanger. The studied heat exchanger was “double pipe” type with only one passage of every stream. The streams were a hot-water stream (as the stronger stream) and a cold-water stream (as the weaker one). The mass flow rate of the weaker stream was kept constant and equal to 0.002 kg/s. The values of the mass flow rate of the stronger stream were 0.004, 0.006, 0.008 or 0.01 kg/s. As can be seen from Figure 158 Strojniški vestnik - Journal of Mechanical Engineering 53(2007)3, 158-164 Heating Tank h He N ter High Flow Meter Control Va lv e s Fig. 1. The experimental rig 1, the inlet and outlet temperatures of the streams were measured. In each of four cases (for the each value of the mass flow rate), the inlet temperature of the stronger stream was kept constant and the inlet temperature of the weaker stream was varied four times. From the obtained four sets of the measured data, for each of the four cases, both for a parallel and a counter-flow heat exchanger, two non dimen- 7l1 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 ¦ ^3=0.50 n3=04> * ^3=0.33 • ^3=0.25 /+ ¦ ^3=0.20 /^* f V3=0.50 sional parameters (p1 and p3) were calculated. The parameters are presented in the appropriate diagrams, as shown in Figures 2 and 3. The total heat exchanger area was 0.08 m2. For each point the exchanged-heat flow rate and the exergy destruction (entropy generation) are calculated by using the equations of the following mathematical model. 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 7l1 0.2 0.4 0.6 0.8 712 1.2 0.2 0.4 0.6 0.8 7l2 1.2 Fig. 2. The positions of the measured operation points for a parallel-flow heat exchanger Fig. 3. The positions of the measured operation points for a counter-flow heat exchanger N 0 0 Energijska in eksergijska analiza - Energy and exergy analysis 159 Strojniški vestnik - Journal of Mechanical Engineering 53(2007)3, 158-164 2 MATHEMATICAL MODEL The exchanged-heat flow rate between two streams can be calculated using the following equation ([10] and [11]): ß = C1 (^-^) = C2 (5'-4) (1). The heat exchanger effectiveness e is usually defined as: ¦$-f1' '_ Q 1 2 Ä^max (2). The exchanged-heat flow rate Q& can be rewritten as: ß = ^fimax=^1 (51'-4) (3). Since the pressure drop of the two streams is neglected, the entropy generation can be calculated from the following equation [11]: gen 1 T T C ln + Cln 2' ' (4). Using Equations (1) and (2), it is easy to transform Equation (4) into the form: ' ^ln 1-e 1 C ln Qe\T ' c r, 1+1M1-1|(5). If the following terms for two non-dimensional parameters c (6) are introduced into Equation (5), the equation for the entropy generation assumes the following form: 5gen=Qln(1-e(1-^T)) + C2ln 1+*,*—-1| (7). After multiplying the above equation by the environmental temperature, T0, the equation for the irreversibility or the exergy destruction is obtained as: ^^ =r0šgen =r0kln ( 1-f( 1-^T )) + C2ln 1+^J1 -1j (8). The above equation can be written in a non-dimensional form. For that purpose, it will be divided by the product of the heat capacity rate of the weaker stream, C1, and the environmental temperature, T0, as follows: 1=^ = ln(1-f(1-^T)) + -ln T0& 1+M3L\-------1 7T (9). Because the scope of this work is to show the ratio of the irreversibility and the exchanged-heat flow rate, it is useful to write Equation (1) in a non-dimensional form, dividing it by the product C1T1': Q (10). Finally, Equation (9) is divided by Equation (10), and a relevant ratio is obtained: ln(1-f(1-a-T)) + ln 1+7lM------1 E\1-71 (11). For a parallel heat exchanger the effectiveness e is obtained using the following formula [10]: -(1+»3K 1-e{ 1 + ^3 (12). and for a counter heat-transfer heat exchanger e is estimated as [10]: 1-e (1-»3 )»2 1 - jie (1-»3)»2 where: Q (13). (14). As can be seen, this ratio is a function of the non-dimensional parameters p3 and e = p1, which represent the operation points of the heat exchanger and which are obtained by the energy analysis of heat exchangers. The additional parameters for irre-versibility or exergy destruction are the input temperature ratia, pT, and the ambient temperature, T0. 3 CALCULATION RESULTS AND DISCUSSION The diagrams in Figures 4 and 5 present the values of i1 and i1/q1 for a parallel-heat-flow heat exchanger, and the diagrams in Figures 6 and 7 for a counter-flow heat exchanger. These values are calculated using Equations (9) to (11) and Equations (12) and (13) for the parallel and the counter heat exchanger, respectively. The black points in the diagrams represent the values of the operating points, which are presented in Figures 2 and 3. 2 160 Ga/ovic A. - Živič M. - Can A. Strojniški vestnik - Journal of Mechanical Engineering 53(2007)3, 158-164 0.012 0.01 0.008 0.006 0.004 0.002 0 0.5 0.6 0.7 0.8 0.9 a) ;r=0.50 0.7 0.65 0.6 E 0.55 0.5 0.45 0.4 0.35 0.3 0.012 0.01 0.008 0.006 0.004 0.002 0 0.5 0.6 0.7 0.8 0.9 0.7 0.65 0.6 e 0.55 0.5 0.45 0.4 0.35 0.3 7h b) tt3=0.33 0.012 0.01 0.008 0.006 0.004 0.002 0 0.5 0.6 0.7 0.8 m" 0.7 0.65 0.6 0.55 0.5 0.45 0.4 0.35 0.3 0.012 0.01 0.008 0.006 0.004 0.002 0 cK=0.25 0.7 0.8 0.9 d) ^=0.20 0.7 0.65 0.6 0.55 0.5 0.45 0.4 0.35 0.3 33 Fig. 4. The non-dimensional irreversibility i1 and the effectiveness e as a function of p2 and the parametric curves pT for a) p3=0.50, b) p3=0.33, c)p3=0.25 and d) p3=0.20 for a parallel flow heat exchanger 0.12 0.1 0.08 0.06 0.04 0.02 0 0.5 0.6 0.7 0.8 0.9 a) ;r=0.50 0.12 0.1 i °-08 — 1\ 0.06 0.04 0.02 0 3 0.12 0.1 i, 0.08 — * 0.06 0.04 0.02 Sj=1.07 0'5 0^6 0^7 0^8 0^9 ^=1.11 m *=1M b)?r=0.33 ^=1.17 0.12 0.1 0.08 h 0.5 0.6 0.7 0.8 0.9 ih 1\ 0.06 0.04 0.02 0 cW^O.25 0.5 0.6 0.7 0.8 0.9 7I2 d) ;t=0.20 Fig. 5. Relative non-dimensional irreversibility i1/q1 as a function of p2 and the parametric curves pT for a) p = 0.50, b) p = 0.33, a) p = 0.25, a) p = 0.20 for a parallel-heat-flow heat exchanger h h 1.1 1.2 1.1 1.2 h rc=1.07 e TT = 1.11 TT =1.14 TT =1.17 h h 1.1 1.2 0.5 0.6 1.1 1.2 h 1.1 1.2 1.1 1.2 1.1 1.2 1.1 1.2 Energijska in eksergijska analiza - Energy and exergy analysis 161 Strojniški vestnik - Journal of Mechanical Engineering 53(2007)3, 158-164 0.012 0.01 0.008 0.006 0.004 0.002 0 0.5 0.6 0.7 a) ^=0.50 0.7 0.65 06 E 0.55 0.5 0.45 0.4 0.35 0.3 0.012 0.01 0.008 0.006 0.004 0.002 0 0.012 0.01 0.008 0.006 0.004 0.002 0 0.5 0.6 0.7 0.7 0.65 0.6 0.55 0.5 0.45 0.4 0.35 0.3 T: =1.01 K=\.\\ K=\.\A 7l=\M 0.5 E 0.8 0.9 b) tt3=0.33 0.7 0.65 0.6 0.55 0.5 0.45 0.4 0.35 0.3 0.8 0.9 7h 0.012 0.01 0.008 0.006 0.004 0.002 0 0.5 0.6 0.7 0.8 0.9 0.7 0.65 0.6 e 0.55 0.5 0.45 0.4 0.35 0.3 c) p3 =0.25 d) p3=0.20 Fig.6. Non-dimensional irreversibility i1 and effectiveness e as a function of p2 and the parametric curves 12 for a) p =0.50, b) p =0.33, c)p =0.25 d) p =0.20 for a counter-flow heat exchanger o.i h 0.08 1\ 0.06 0.04 0.02 0 0.5 0.6 0.7 0.8 0.9 p2 a) p =0.50 ^ 0.06 0.04 0.02 0 0.12 0.1 h_ 0.08 11 0.06 0.04 0.02 0 ^=1.11 0.5 0.6 0.7 0.8 0.9 b) 71=033 0.5 0.6 0.7 0.8 cK=0.25 0.9 0.12 0.1 0.08 0.06 0.04 0.02 0 0.5 0.6 0.7 0.8 0.9 d) ?r=0.20 Fig. 7. Relative non-dimensional irreversibility i1/q1 as a function of p2 and the parametric curves pT for a) p =0.50, b) p =0.33, c) p =0.25, d) p =0.20 for a counter-flow heat exchanger ii h e 1.1 1.2 0.6 0.7 1.1 1.2 h ii H E 1.1 1.2 1.1 1.2 p 0.12 0.12 0.1 L 0.08 h 1.1 1.2 1.1 1.2 k=1.07 7TT=1.14 Kj=1.17 1 1.1 1.2 1.1 1.2 162 Galovič A. - Živič M. - Can A. Strojniški vestnik - Journal of Mechanical Engineering 53(2007)3, 158-164 The above diagrams clearly show the relationship between the dimensionless values of i1 and the dimensionless parameters p2, p3 and pT. For each diagram it is possible to directly read the values of i1 and e for every case. It is obvious that the given values of p2 and p3 do not have a significant influence on i1, but they have an influence on e. The pT ratio has the greatest influence on i1, but it has no influence on e. This can be quantitatively seen from Figures 4a and 4d, where the values p3 = 0.5 and 0.2 and parametric values pT =1.11; 1.14 and 1.17 have the same value of p2 = 0.76 and 1.09 respectively. For p3 = 0.5, all operation points have the same value of e, i.e., e = 0.455, with the exergy destruction increasing from 0.0028 to 0.0078. On the other hand, for p3 = 0.2 and the same parametric values of pT and p2 of 1.09, the value of e is equal to 0.61. The associated exergy destructions are 0.0040, 0.0067 and 0.0076 respectively. It is obvious that parallel- and counter-heat-flow heat exchangers have very small values of the ratio i1 and i1/q1 for the measured operation points. It is not possible to make a comparison with respect to exergy destruction and heat-transfer effectiveness of the researched cases, because they did not have 4 CONCLUSION The presented analytical relationship between the dimensionless exergy destruction and the heat-transfer effectiveness of a heat exchanger seems to be rather convenient, because it relates the dimensionless parameters (p3; p 2 and pT as additional parameter) relevant to the operation of a heat exchanger and the exergy destruction of a parallel and a counter-flow heat exchanger. It is possible to include the values of the measured operation points into the presented mathematical model and simultaneously calculate the heat-transfer effectiveness and the exergy destruction of the considered heat exchangers. It can be concluded that the exergy destruction for each case of both investigated heat exchangers is small. The main reason for such a conclusion is the fact that the operation values of pT close to 1. Furthermore, by introducing additional exergy destruction due to the pressure drop in the model, it is possible to develop an analytical model for the minimization of heat exchangers’ exergy destruction. the same non-dimensional p variables. 5 NOMENCLATURE A0 C Ex i1 k Q& q1 S& T T overall heat-transfer area, m2 heat-capacity rate of the stream, W/K exergy, W exergy destruction, W non-dimensional exergy destruction overall heat-transfer coefficient, W/(m2K) heat-transfer rate, W non-dimensional heat-flow rate entropy generation rate, W/K thermodynamic (absolute) temperature, K ambient temperature, K Greek Letters J pT p2 p3 Celsius temperature, oC ratio of inlet absolute temperatures, K number of heat-transfer units ratio of heat-capacity rates e, p heat-exchanger effectiveness Subscripts 1 weaker stream 2 stronger stream destr destruction gen generated Superscripts ‘ inlet ‘’ outlet 6 REFERENCES [1] Bejan, A., (1977) Second law analysis in heat transfer, Energy, 5, (1977), pp. 721-732 . [2] Rant, Z. (1956) Exergie, ein neues Wort für technische Arbeitsfähigkeit. Forschung Ing. Wesens 22(1956), pp. 36-37. [3] Rant, Z. (1964) Thermodynamische Bewertung der Verluste bei technischen Energieumwandlungen, BWK 16(1964), pp. 453-457. Energijska in eksergijska analiza - Energy and exergy analysis 163 Strojniški vestnik - Journal of Mechanical Engineering 53(2007)3, 158-164 [4] Can, A., Buyruk, E., Eryener, D., (2002) Exergoeconomic analysis of condenser type heat exchangers, Exergy- an International Journal 2(2002) , pp.113-118. [5] Galovič, A., Virag Z., Mudrinič, S., (2003), Non-dimensional entropy analysis of condenser and/or evaporator type heat exchangers, Transactions of Famena, 27(2003), pp. 1-9. [6] Gregorig, R. (1973) Wärmeaustausch und Wärmeaustauscher, Verlag Sauerländer, Frankfurt am Main. [7] Galovič, A., M. Živič, M. Andrassy, (2003) Entropijska analiza sotočnih prenosnikov toplote, Strojniški verfwl49 (2), pp.l00-110. [8] Bejan, A., (1988) Advance engineering thermodynamics, John Willey and Sons, New York. [9] Bejan, A. (1996) Entropy generation minimizaton, CRC Press, New York. [10] Bošnjakovič, F., K. F. Knoche, (1997) Technische Thermodynamik, Teil II, Steinkopf Verlag, Darmstadt. [11] Çengel, Y,Boles, M., A., (1994) Thermodynamics: an engineering approach, McGraw-Hill Book Company, New York Authors’ Addresses: Prof. Dr. Antun Galovič University of Zagreb Faculty of Mechanical Engineering and Naval Architecture Ivana Lučiča 5 10000 Zagreb, Croatia antun.galovic@fsb.hr Dr. Marija Živič J.J. Strossmayer University of Osijek Mechanical Engineering Faculty 35000 Slavonski Brod, Croatia marija.zivic@sfsb.hr Prof. Dr. Ahmet Can Trakya University Faculty of Mechanical Engineering and Architecture 22030 Edirne, Turkey Prejeto: Sprejeto: Odprto za diskusijo: 1 leto 7.2.2006 25.10.2006 Received: Accepted: Open for discussion: 1 year 164 Galovič A. - Živič M. - Can A.