Generalize heat conuction transfer functions in two imensional linear systems of assigne geometry F. Marcotullio & A. Ponticiello Dipartimento i Energetica- Universit i L 'Aquila Monteluco i Roio - 67100 L'Aquila, Italy ABSTRACT The present paper eals with the search for transfer functions vali for the whole class of object of an assigne two imensional geometry. The basic criterion was to turn the governing conuction equation an the transient bounary conitions as well into imensionless variables. The conclusion was that the same the TFs hol for all systems provie the Nusselt number (Nu) an the Fourier moulus (Fo) are equal. An appropriate TF's coefficients table has been set-up for the central an the corner points of a homogeneous soli square profile. The imensionless time interval AFo was selecte accoring to the sampling theorem, i.e. to the cut-off frequency of the system which epens on Nu. In aition, the Nu incremental steps were set in orer to provie "equal" percent eviation (10%) an thus avoi the necessity of temperature interpolations. The generalize TF's coefficients enables a wier use of the TF approach which is very useful for thermal transient analysis. INTRODUCTION The Z-transfer function approach represents a very powerful tool for preicting temperature histories in multiimensional thermal fiels, e.g. Mascheroni [1]. It is particularly useful, for instance, in simulating the ynamic behaviour of linear systems uner ifferent transient bounary conitions, e.g. Stephenson an Mitalas [2]. The search for the transfer functions (TF) in multiimensional heat conuction by a finite-element technique was envisage by the authors [3]. The basic criterion was that of rearranging an conensing the matrix finite-element equation to reuce the egrees of freeom to only one. The conensing proceure will yiel, in general, aitional terms (higher orer time-enva-
36 Heat Transfer tives) in the resulting "lumpe" matrix equation. At this stage, the numerical solution in the time-omain iscretization was obtaine explicitly using an appropriate recurrence formula. The coefficients of the Z-transfer functions were thereby compute from the cite recurrence scheme by matching the system free response in orer to set the correct recurrence parameters. In this stuy, the matrix conuction equation is turne into imensionless variables reaching the conclusion that the same TFs hol for the whole class of objects having the same geometric shape (plates, ros, bars, cyliners...) provie the Nusselt number (Nu) an the Fourier moulus (Fo) are equal. This expecte but relevant result enables a wier use of the TF approach, otherwise relegate to rare applications since the TF's ientification presupposes a certain amount of skilfulness which is infrequent in common practice. So far, it will suffice to set-up, at preselecte points an for iffering geometric shapes of practical interest, appropriate tables listing the series of the TF's coefficients for pairs of the inepenent imensionless variables Nu an Fo BASIC CONSIDERATIONS The governing conuction equation in a uniform thermal conuctivity soli ii (x,y,z) free of sources is given by: x* Let our task be the search for the general temperature history within the soli unergoing the following initial an transient bounary conitions: To turn Equation (1) into imensionless variables let: ' - - where L<. is a characteristic length an Fo is the Fourier moulus. By substituting in Equation (1) we obtain:
Heat Transfer 37 with the bounary conition expresse as: T (A\ ul ~. XT,, /nr T \ = Q (5) vv < where Nu=hLJk is the Nusselt number. It follows that the analytical solution T will be a function of the imensionless variables %,T,^, Fo, Nu an will not epen explicitly upon the soli thermophysical properties which, instea, are inclue in Fo (iffusivity a) an Nu (conuctivity *). The same conclusions can be rawn by performing a numerical analysis base, for example, on the finite-element metho. In this case, the iscretization of the whole omain Q (%,T,Q in a certain number M of finite elements leas to the well-known matrix equation Ca+Ka-f = 0 (6) here a is the vector of the unknown noal temperature of the system an a the first time-erivative with respect to the Fourier moulus Fo. All the matrices are assemble from the element submatnces in the stanar manner with submatnces K, CT an f given by: K = " j\$ $ T) ri t c^q ; (7a) (7b) (7c) where a is the ambient flui temperature varying with time (forcing function) an Ni, N^ the interpolation an weighting functions. Now, since the
38 Heat Transfer transfer functions (TF) at any selecte noe (,T,Q of the fiel are obtaine from the matrix Equation (6) by means of the proceure escribe in [3], it follows that exactly the same TF's hol for any assigne geometry provie the interpolation functions N(,r,Q, the Nusselt number Nu an the imensionless time interval AFo=Aia/L/ are equal. The transfer functions are given as sets of coefficients each relate to its own forcing function represente in general by m inepenent surrouning temperatures. The time-series of the preselecte noal temperature # +, are linke to each timeseries of the input # +, by means of a convolution equation such as: m ^ Nu AND AFo SELECTION To set-up a table listing the series of TF's coefficients for pairs of the inepenent imensionless variables Nu an AFo to be use without interpolation, some appropriate criteria shoul be followe. Firstly, the Fourier moulus iscretization AFo is selecte at least equal to X/15 accoring to the sampling theorem, X being the ominant eigenvalue of the matrix K C obtaine from Equation (6). The above conition implies that the system frequency response to the unitary step-function of the ambient flui temperature is ampe out at approximately 2% of its maximum value. In fact, the temperature history at any arbitrary point is given by: which, as time progresses, can be aequately represente by: T o) The Laplace transfer function is: an by setting the Laplace parameter s=27wi the moulus IW(s)l takes the form:
Heat Transfer 39 Table 1 - TF's coefficients for a soli homogeneous square profile 2^ i Nu AFo xlo* L;= t L ^~ I 4-/Vi-y+ tf-/v,-/ = <> 2=0 /=0 0.2 0.6 1.0 1.4 2.0 3.0 4.8 8.3 17 47 10* 86.12 30.60 19.52 14.78 11.25 8.53 6.51 5.13 4.20 3.66 3.41 b* b~ b> b^ b^ b- n + 1 0.006121 0.077561-0.000538 0.125824-0.000136 0.165626-0.000174 0.185363-0.000067 0.225369-0.000425 0.282798-0.000113 0.348444-0.000257 0.468378-0.000318 0.653359-0.000163 0.875796 0.000020 1.000003 n 0.053744 0.000113 0.964670 0.026946-0.038861 1.273946 0.010017-0.112147 1.524058 0.004604-0.144038 1.768537 0.001363-0.245666 1.983495 0.002131-0.402290 2.197197 0.000615-0.619427 2.584926 0.001524.005411 2.766646 0.001805.619145 2.862545 0.000838-2.401356 2.881824-0.000167-2.875266 2.875257 n 0.002747-0.015061-0.027282 0.019283-0.056898-0.335979 0.024448-0.042214-0.623488 0.019962-0.083457-0.967305 0.014111-0.032631.301144 0.003669 0.075296.669013 0.001004 0.210193-2.498636-0.003283 0.625223-2.914466-0.004215 1.375844-3.145066 -.001759 2.450002-3.192349 0.000505 3.176002-3.175991 n-2-0.001775 0.013850 0.018118-0.003093 0.019971 0.068193-0.001186 0.080627 0.186106 0.002830 0.091245 0.327282 0.010261 0.083987 0.514082 0.009362 0.178669 1.123435 0.010728 0.014074 1.458082 0.008473-0.399616 1.655889 0.003562.122439 1.697446 0.000376.682773 1.682767 n-3-0.000973-0.016261-0.009571-0.002402-0.022483-0.025468-0.004666-0.028821-0.053237-0.004761-0.137980-0.233970-0.003752-0.129534-0.345824-0.000382-0.031762-0.417386 0.004002 0.208447-0.432853 0.006433 0.427227-0.427226 n-4 0.000084 0.026291 0.018053-0.000486 0.031745 0.031088-0.001612 0.025071 0.040267-0.002860-0.006830 0.042312-0.003506-0.041532 0.041532
40 Heat Transfer 1 (12) If a cut-off frequency v, is chosen so as to obtain for Iwl a 2% attenuation, then v,«50/27txi an, corresponingly, AFo=l/2v«.=A,i/15. As it clearly appears from the set of Equations (7), since A,, epens upon Nu then the above time interval selection criterion allows us to set-up the TF's coefficients table with only one entry (Nu). Another problem that arose was the search for the appropriate Nusselt incremental steps ensuring "equal" percent temperature variations. For general engineering purposes, a 10% sprea will suffice to hol the transient response within 5% maximum eviation an thus avoi teious temperature interpolations. As an application example, the TF's coefficients for the central an the corner points of an infinite homogenous plate 0.6m by 0.6m of a soli square profile subject to transient convective flow are reporte in Table 1. Assuming the characteristic length L<=L=0.6m, the unit surface conuctance h=5.8wm^k"'an the thermal conuctivity k=0.65wm"*k% then Nu=hLJ( ^=5.35 against 4.8, which is the nearest Nu value in Table 1. Temperature histories for both noes 1 an 2 were plotte by using the series of coefficients corresponing to Nu=4.8 an AFo=6.510^ an assuming (Figure 1) an ambient flui linear rise temperature (3 C/h) an (Figure 2) a triangular pulse 10 C amplitue an 2Ai uration (AT=AFo-L//a=1.42h). The plots were checke in terms of actual time intervals against the results eriving from the use of the exact sets of TF's coefficients (Table 2). It is confirme that all eviations are within the chosen 5% maximum error range. Table 2 - Exact FT's Coefficients ( Ai=1.42h) Nu n+1 n n n-2 n-3 n-4 535 b* -0.000145 0.373902 0.000755-0.659232 2.551457 0.000605 0.217875-2.428483 0.010882 0.193255 1.072345-0.005264-0.146653-0.218698 0.000075 0.027753 0.016471 CONCLUSIONS The following comments can be emphasize from the above stuy: the availability of the transfer function coefficients for the whole class of objects of assigne geometry in terms of imensionless variables such as the Nusselt number Nu an the time interval AFo, enables a wier use of the
Heat Transfer 41 328 323 318 313 308 =0.65 Wm'iK'1 = 1700kgm" =0.60m = 1.647 10" =5.8 Wm'^K Noe! - Nu=4.8 Noel-Nu=5.35 Noe2 - Nu=4.8 Noe2-Nu=5.35 283 278 273 427 712 996 12,81 15,66 18,5 21,36 5,69 ' 8,54 11,39 14,23 17,08 19,93 Time (h) Figure 1 - Temperature histories at noes 1 an 2 with an ambient flui linear rise temperature (3 C/h) '"t" k=0.65 Wnr'K = 1700kgm'* c=836jkg-'k-' L =0.60 m a =1.647 10 2,85 5,69 8,54 11,39 Time (h) 14,23 19,93 21,35 Figure 2 - Temperature histories at noes 1 an 2 with an ambientfluitriangular pulse 2Ar=2.85 h uration an 10 C amplitue
42 Heat Transfer transfer function approach. appropriate tables listing the TF's coefficients have been set-up by choosing for the Nusselt number incremental steps of "equal" percent eviation (10%). This criterion permits the irect use of the sets of coefficients by referring to the closest Nu value, without the necessity of teious temperature interpolations. Moreover, by choosing the time interval AFo=l/2Vc=X/15, which epens on Nu, the TF's coefficients table turns into a more suitable one-entry table. REFERENCES 1. Mascheroni, R.H., Sanz, P.D. an Dominguez, M. 'A New Way to Preict Thermal Histories in Multiimensional Heat Conuction: the Z- Transfer Function Metho' Int. Comm. Heat & Mass Transfer, Vol.14, pp. 81-90, 1987 2. Stephenson, D.G. an Mitalas, G.P. 'Calculation of Heat Conuction Transfer Function for Multilayer Slabs' ASHRAE Transactions, 1971 3. Marcotullio, F an Ponticiello, A. "Determination of Transfer Functions in Multiimensional Heat Conuction by Means of a Finite Element Technique' (Eite by R.W. Lewis, University of Swansea U.K.), Vol.VIII, part 1, sect. 2, pp. 226-236, Proceeings of the Eighth International Conference on Numerical Methos in Thermal Problems, Swansea, Wales, 1993 4. Salvaori, U.O., Reynoso, R.O. an Mascheroni, R.H. 'The Use of Z- Transfer Functions to Evaluate Temperature an Quality Changes in Refrigerate Storage' Proceeings of meetings of Comm. B2, C2, Dl, D2/3, Int. Institute of Refrigeration, Paris, pp. 695-700, 1990 5. ASHRAE Hanbook of Founamentals, 1977, Chap. 25 6. Zienkiewicz, O C The Finite Element Metho McGraw-Hill, Lonon, 1977 7. Huebner, K.H. The Finite Element Metho for Engineers J Wiley&Sons, New York, 1975 8. Rao, S.S. The Finite Element Metho in Engineering Pergamon Press, Oxfor, 1989 9. Schneier, P.J. Conuction Heat Transfer Aison-Wesley Publishing Company, Inc, 1974