AP08_6.vp 1 Introduction Stirling’s engine (Fig. 1) is a volumetric engine, in which work is done by changing the volume, the pressure and the temperature of the working gas. The working gas is moved by pistons on the hot and cold side through the regenerator. The motion of the pistons is regular, and the pistons are mechani- cally connected. The pistons are frequently moved by a crank shaft. The power output is through the crank shaft of the engine. One cycle of the Stirling engine is one turn of the crank shaft. The regeneration of the heat inside the regenerator is perfect according to Schmidt’s idealization. If the working gas is flowing from the hot side to the cold side, the working gas is cooled in the regenerator from temperature TH to tempera- ture TC. Heat is saved to the matrix of the regenerator during cooling. If the working gas is flowing from the cold side to the hot side, the working gas is heated from temperature TC to temperature TH inside the regenerator. The heat is taken from the matrix of the regenerator during heating. This pro- cess is called regeneration of the heat. The temperature of the working gas on the hot side and on the cold side is not con- stant during the cycle in a real Stirling engine (Fig. 1). The amount of regenerated heat inside the regenerator of the Stirling engine is much greater than the heat put into the engine (Urieli and Berchowitz – simply adiabatic analysis [1]). However no an analytical computing method for the amount of regenerated heat inside the regenerator has existed untill now. The amount of regenerated heat has an elemental action on the thermal efficiency, the performance and the dimen- sions of the engine. Schmidt’s idealization [2] is a comparative cycle of the en- gine with outer transfer of the heat (Stirling engine), when isothermal processes are taking place in all volumes of the en- gine Fig. 3. Schmidt’s idealization has been the only analytic computing method for thermodynamic design of Stirling engine untill now. The heat flows, the work autput and the thermal efficiency of the cycle without wastage of energy are computed on the basis of the assumptions of Schmidt’s ideal- ization. The assumptions of Schmidt’s idealization are stated below. 2 Notation symbol description unit cp specific heat capacity for working gas at constant pressure J�kg�1�K�1 I enthalpy of working gas J m total amount of working gas in engine kg p pressure Pa Q energy balance (heat) J r individual gas constant of working gas J�kg�1�K�1 T absolute temperature K V volume m3 � phase angle, of the hot side volume variation to the cold side volume variations rad � thermal efficiency - 10 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 48 No. 6/2008 The Amount of Regenerated Heat Inside the Regenerator of a Stirling Engine J. Škorpík The paper deals with analytical computing of the regenerated heat inside the regenerator of a Stirling engine. The total sum of the regenerated heat is constructed as a function of the crank angle in the case of Schmidt’s idealization. Keywords: Stirling engine, regenerator, regenerated heat, Schmidt’s idealization. V H =max. V C =max. V D V HD V R V CD V H (�) V C (�) T C T H T R T [K] V(�) V HC (�) V CC (�) hot side cold side regenerator V HC,max V CC,max Fig. 1: Scheme of a Stirling engine with dead volumes and prog- ress of the temperature of the working gas for crank angle � � crank angle rad � adiabatic index - � temperature rate - subscripts C cold side Car Carnot CC volume of cylinder of cold side D death (death volume) H hot side HC volume of cylinder of hot side max maximum min minimum R regenerator reg regenerated 3 The amount of regenerated heat inside the regenerator The amount of regenerated heat inside the regenerator of the Stirling engine can be found from an energy balance cycle. This energy balance can be applied to the Stirling engine (Fig. 1) under the following assumptions: 1) The working gas is an ideal gas. 2) There is no pressure loss, and the pressure is the same on all the volume. 3) The engine is totally sealed. 4) There is no heat transfer between the matrix of the regen- erator and the structure of the engine. 5) Steady state conditions are assumed for overall operation of the engine so that the pressures, temperatures, etc. are subject to cyclic variations only. The energy balance of the cycle on a part of a cycle 1-2 in the T-s chart (Fig. 2) is described by the following this equation d H C RQ Q Q Q Q S S 1 2 1 2 1 2 1 2 1 2� � � � �� � � �, , , . (1) The thermal flows and the entropy are functions of the crank angle, so equation (1) is d H C RQ Q Q Q Q � � 1 2 1 2 1 2 1 2 1 2� � � � �� � � �, , , . (2) The energy balance of the working gas inside the regener- ator on a part of cycle 1-2 can be found using the equation (2) Q Q Q QR H Cd, , ,1 2 1 2 1 2 1 2 � � �� � �� � � . (3) The energy balance of the working gas on the hot side of the engine on a part of cycle 1-2 � �Q I V p I V pH H H H Hd d d, ( ) ( ) ( )1 2 1 2 1 2 1 2 1 2 � � � � �� � � � � � � � � � � � �� . (4) The energy balance of the working gas on the cold side of the engine on a part of cycle 1-2 � � Q c m T V p I V pC C C C C C d d, ( ) ( ) ( ) ( 1 2 1 2 1 2 1 2 � � � � � � �� � � � � � � � � � ) .dp � � 1 2 � (5) For the energy balance of the cycle on a part of cycle 1-2, we use the first law of thermodynamics � � d d d d Q c m T V p c m T V p p � � � � � � � � � � � 1 2 1 2 1 2 1 2 � � �� � � � ( ) ( ) ( ) p � � 1 2 � . (6) The energy balance of the working gas inside the regener- ator on a part of cycle 1-2 can be found by substituting equa- tions (4), (5), (6) into equation (3) � � � � Q c m T V p I V p pR H H d d , ( ) ( ) ( ) ( ) 1 2 1 2 1 2 1 2 � � � � � �� � � � � � � � � � � � � 1 2 1 2 1 2 � � � � � � � � �� �I V pC C d( ) ( ) , (7) whereas for a medium temperature working gas in the engine (from the state equation) is T p V r m ( ) ( ) ( ) � � � � . (8) The total volume is the sum of all working volumes of the engine V V V V( ) ( ) ( ).� � �� � �H R C (9) The energy balance of the working gas inside the regener- ator on a part of cycle 1-2 can be found by substituting equa- tions (8), (9) into equation (7) © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 11 Acta Polytechnica Vol. 48 No. 6/2008 S [J·K-1] T [K] 1 2 S 1 S 2 S 1 S 2 dQ cycle Fig. 2: General thermal cycle in the T-s chart � � � � � � Q p V V p I I R R H C , ( ) ( ) ( ) ( ) ( 1 2 1 1 2 1 2 1 2 � � � � � � � � � � � � � � � � � � � � �) .� � 1 2 (10) This equation is the energy balance of the regenerator on a part of cycle 1-2. The total sum of the regenerated heat inside the regener- ator in the state �1 0� can be found from equation (10) � � � � � � � � Q p V V p I I R R H C , ( ) ( ) ( ) ( ) ( ) 0 0 0 0 0 1� � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � 1 0 0 0 0 p V p V V p p I I ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) . � � � � � � R H C (11) The pressure p(�) and the volume of the engine V(�) for a specific crank angle is calculated or measured. Fig. 4 describes the progress of the total sum of the regenerated heat inside the regenerator (11) as a function of the crank angle. From equations (4), (5) and (11) we can calculate the heat input to the engine, the heat output from the engine and the regener- ated heat inside the regenerator on a part of cycle 1-2. The difference between the maximum and minimum val- ue of the function Q R, 0�� is the amount of regenerated heat inside the regenerator in the cycle Q Q Qreg � �� �R R, , max , , min0 0� � . (12) The amount of regenerated heat can be computed exactly if the change enthalpy of the working gas on the hot side and the cold side � �I IH C( ) ( )� � �� 0 is known (more about this problem in [3], [4]). If the thermodynamic processes are iso- thermal on the hot side and the cold side, then the enthalpy of the working gas does not change (� �I IH C( ) ( )� � �� �0 0) and we can using equation (13). The thermodynamic processes are not isothermal in real Stirling engines. These changes of the temperature (enthalpy) of the working gas on the hot side and the cold side cannot be computed without measure- ment. However the amount of the regenerated heat inside the regenerator can be computed from equation (13) for the assumption � �I IH C( ) ( )� � �� 0 � � � � � � � � � � � � � � � � 1 1 2 1 2p V V p( ) ( ) ( )R � � � � Q p V p V V p p R R , ( ) ( ) ( ) ( ) ( ) ( ) . 0 1 0 0 0 � � � � � � � � � � � � (13) The value of the regenerated heat computed from equa- tion (13) is greater than the value of the regenerated heat computed from equation (11). The regenerator designed by equation (13) is therefore greater. The derivation of equation (13) was obtained using dif- ferential equations by a prof. Uriel, but only for conditions: TH const� , TC const� . 4 A calculation of the amount regenerated heat inside the regenerator for a Stirling engine fulfilling the assumptions of Schmidt’s idealization From equation (11) or (13) we can compute the amount of regenerated heat inside the regenerator of the Stirling en- gine if computed or measured the functions p(�) and V(�) are known. This section presents a method for calculation the amount of regenerated heat inside the regenerator in the cycle of the Stirling engine, fulfilling the assumptions of Schmidt’s idealization. Schmidt’s idealization can be applied to the engine in Fig. 3, under the following assumptions: 1) The temperature of the working gas which flows from the regenerator on the hot side is TH, and the temperature of working gas which flows from the regenerator on the cold side is TC. 2) There is no pressure loss, and the pressure is the same throughout the volume. 3) The working gas is an ideal gas. 4) The engine is totally close. 5) Sinusoidal motion of the pistons. 6) The temperature of working gas on the hot side is con- stant and equal to TH, and the temperature of working gas on the cold side is constant and equal to TC. The mean temperature of working gas in the regenerator is constant and equal to TR. 7) Heat enters the engine only through the walls on the hot and cold sides of the engine. 8) There is perfect regeneration. The summary equations for Schmidt’s idealization are evolved from the following assumptions [2]: Pressure p r T m V A B ( ) ( cos( )) , ,max � � � � � � H HC 1 2 (14) A V T V T V T V T k T T � � � � � � � � � � � � 1 2 1 1 HC H HD H R R CD C R C ,max ,� � � � � � � � � � � � � T T T k V V T T B x z H C H CC HC H C ln , , , ,max ,max 1 2 � 2 1 1 1x k z k z x k � � � � � � � � � � � � � � cos sin , arctan arctan 1 11 � � � � � � � � � � � sin cos . k 12 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 48 No. 6/2008 Volume of the hot side V V VH HC HD( ) ( )� �� � . (15) Volume of the cold side V V VC CC CD( ) ( )� �� � . (16) Heat input to the engine Q r T m k B A A B H H� � � � � � � � � � � � � � � 1 2 2 2 1sin( ) . (17) Heat output from the engine Q r T m k B A A B C H� � � � � � � � � �1 2 2 2 1sin( )� � � . (18) Thermal efficiency of the cycle � � � � � A Q H car1 1 . (19) As is evident from equation (19), an engine the fulfilling assumptions of Schmidt’s idealization has thermal efficiency equal to the thermal efficiency of Carnot’s for temperature rate T TH C. The temperature of the working gas on the hot and cold sides of the engine is constant. Therefore the enthal- py of the working gas on these sides does not change � �( ( ) ( ) )I IH C� � �� �0 0 ) and we can use the equation (13). © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 13 Acta Polytechnica Vol. 48 No. 6/2008 V H (�) T [K] T H =constant T C =constant T R 0 � � V C =max. V C (�) 0 � V H =max. V(�) Fig. 3: Simplified scheme of the Stirling engine, fulfilling the assumptions of Schmidt’s idealization, and the progress of the temperature of the working gas in volumes of the engine Q R,0-� [J] 0 60 120 180 240 300 360 -1000 0 1000 2000 3000 4000 5000 Q reg � [deg] a b Fig. 4: Progress of the sum total of regenerated heat (13) as a function of the crank angle for a Stirling engine fulfilling the assumptions of Schmidt’s idealization (a) – Q R, 0�� diagram for the Stirling engine United Stirling V-160 (working gas - helium, � � 167. , r � 2077 22. J/kg K, TH � 900 K, TC � 330 K, pmean � 15 MPa, VHC,max � 160 cm 3, VHD � 191 cm 3, VCD � 110 cm 3, VR � 28 cm 3, � � 105°) (b) – Q R, 0�� diagram for the same Stirling engine but in the case of V V VHD CD R� � � 0 cm 3 Qreg amount of regenerated heat inside the regenerator in case (a). Other results of these cycles are shown in Table 1 5 Conclusion The amount of regenerated heat inside the regenerator of a Stirling engine can be computed analytically if functions p(�) and V(�) from equations (13) are known. The progress of the total sum of the regenerated heat Q R, 0�� as a function of the crank angle for a Stirling engine fulfilling the assumptions of Schmidt’s idealization is shown in Fig. 4a. This figure shows the amount of regenerated heat inside the regenerator Qreg. From equation (17) and (12), we can compute the ratio of the regenerated heat to the heat input to the engine. This ra- tio is 4.94 (for an engine with the technical specifications in Fig. 4), and it corresponds perfectly with the ratio computed by Urieli and Berchowitz [1]. The larger this ratio is, the greater the impact on the efficiency of regeneration of the heat input to the engine. The progress of the total sum of regenerated heat Q R, 0�� as a function of the crank angle for the Stirling engine fulfill- ing the assumptions of Schmidt’s idealization without dead volumes is shown in Fig. 4b. The ratio of the regenerated heat to the heat input to the engine (1.59) is smaller than for the previous case. Thus, if the dead volumes in the engine de- crease, the influence of the regenerated heat inside the regen- erator on the thermal efficiency of cycle also decreases. References [1] Urieli, I.: Stirling cycle engine analysis, Bristol: A. Hilger, 1984, ISBN 0-85274-435-8. [2] Walker, G.: Stirling Engines, Oxford University Press, 1980, ISBN 80-214-2029-4. [3] Škorpík, J.: Příspěvek návrhu Stirlingova motoru, VUT v Brně, Edice PhD Thesis, 2008, ISBN978-80-214-3763-0. [4] Škorpík, J.: Energetická bilance oběhu Stirlingova motoru, [online] www.oei.fme.vutbr.cz/jskorpik/190.html. Ing. Jiří Škorpík, Ph.D. phone: +420 541 142 575 fax:+420 541 143 345 e-mail: skorpik@fme.vutbr.cz Department of Power Engineering Brno University of Technology Faculty of Mechanical Engineering Technická 2896/2 616 69 Brno, Czech Republic 14 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 48 No. 6/2008 parameter a b QH [J] 969.1 3129.1 QH [J] �355.4 �1147.4 A [J] 613.8 1981.8 QH [J] 4784.1 4963.8 Q Qreg D [-] 4.94 1.59 [-] 0.63 0.63 Table 1: Other results of the Stirling engine cycle fulfilling the assumptions Schmidt’s idealization for the parameters presented in the legend of Fig. 4.