A Potential Method to Predict Performance of Positive Stirling Cycles Based on Reverse Ones

: There are two kinds of working mechanisms for the Stirling cycle, i.e., the positive and the reverse cycles, and a Stirling engine (SE) can be operated as a Stirling refrigerator (SR). This indicates that a probable practical method for evaluating the performance of a Stirling engine is to run it as a refrigerator, which is much easier to operate. For this purpose, an improved Simple model for both the positive and the reverse Stirling cycles, considering the various loss mechanisms and actual operating conditions, is proposed and veriﬁed by a self-designed Stirling engine. As to the positive cycle with helium and nitrogen at 2.8 MPa, the model errors range from 5.4–11.3% for the indicated power, and 1–10.2% for the cycle efﬁciency. As to the reverse cycle with helium and nitrogen, the errors of the predicted input power range from 7.9–15.3% and from 2.5–10.9%, respectively. The experimental cooling temperatures can reach − 92.2 and − 53.6 ◦ C, respectively, for the reverse cycle with the helium and nitrogen at 2.8 MPa. This Stirling-cycle analysis model shows a good adaptability for both the positive and the reverse cycles. In addition, the p-V maps of the positive and reverse cycles are compared in terms of “pressure ratio” and “curve shape”. The pressure ratio of the reverse cycle is signiﬁcantly higher than that of the positive one at the same mean pressure. A method is proposed to predict the indicated work of the positive Stirling cycles using the reverse ones. A mathematical model to predict the indicated power of the positive Stirling cycles based on the reverse ones is proposed: W heat 2 − W cool 1 = A · ( T ge 2 − T gc 2 T gc 1 − T ge 1 (cid:17) B . The most critical issue with this method is to establish an associated model of the temperatures of the expansion and the compression space. This model shows a good adaptability for both the positive and the reverse cycles and can provide detailed information for deep discussion between the positive and the reverse cycles.


Introduction
It is imperative to search for other alternative sources of energy on account of the finite amount of fossil fuels in the world. Solar energy is an attractive option due to its availability in a large portion of the world. The Stirling engine is a reciprocating piston power machine driven by flexible heating sources, which has a broad potential for application in the fields of solar thermal power, waste-heat utilization, and combined cooling, heating, and power systems [1][2][3][4].
The Stirling-cycle analysis methods can be divided into five categories (zeroth-order, first-order, second-order, third-order, and fourth-order analysis methods). Recently, some models have been developed to accurately predict the performance of Stirling engines.
However, most of these models are not suitable for wide speed range, especially at a high rotational speed. Fawad et al. [5] analyzed the heat-transfer and flow-friction effects of the heater, cooler, and regenerator on the performance of the engine using a quasi-steady flow approach. Tlili [6] proposed an optimization method for an endoreversible Stirling heat engine with finite heat-capacitance rates of the external fluids in the heat source/sink reservoirs, which provided a finite time thermodynamic for the design, performance evaluation, and improvement of the Stirling heat engine. Hosseinzade, H. and H. Sayyaadi [7] established an analysis model called CAFS, which combined adiabatic analysis with finite speed thermodynamics. Stirling engines can be designed in different types and capacities, and most performance tests need flexible, stable heat sources, whose heat-exchange efficiency is relatively high, especially for some engines with a high power level, such as the V-160 [8], Cummins Power Generation [9], I-98 [10], I-365 [11], Mod I [8], and some double-acting Stirling engines [4,8]. It is noted that tests at high temperatures are not easy to fulfil. Çınar et al. [12] tested an α-type Stirling engine heated with an electrical heater within the range of a heater temperature of 800-1000 • C with 50 • C increments; however, the maximum output power was only 30.7 W at 1000 • C. Cheng et al. [13,14] developed and tested a beta-type 300 W Stirling engine whose results showed that the shaft power can reach 390 W at 1400 rpm with 32.2% thermal efficiency. They also studied a 1-kW class prototype Stirling engine heated by a 3-kW electric infrared heater whose results showed that the output power may reach a maximum value at a certain rotational speed when the control pressure of the working gas is not higher than 1 MPa. In addition, the Stirling engine with a low temperature difference can be realized easily, and many researchers have designed and tested different Stirling engines with a low temperature difference [15][16][17][18].
There are two kinds of working mechanisms for the Stirling cycle, i.e., the positive and the reverse cycles, and a Stirling engine can be also operated as a refrigerator based on a reverse Stirling cycle. Batooei et al. [19] optimized a γ-type Stirling refrigerator using a multi-objective optimization method and found out that the cooling capacity increased with the rotational speeds while the COP had a maximum value. Guo et al. [20] proposed a general analytical model for various types of Stirling refrigerators with a good agreement with the experiment. Hachem et al. [21] studied a β-type Stirling engine under different operating and geometrical parameters and investigated the effect of the geometric parameters. The results showed that regenerator porosity should be about 85% for the maximum refrigeration power, and the optimal values of the regenerator length and diameter for the prototype were about 60 mm and 22 mm, respectively. Tekin et al. [22] performed a thermodynamic analysis of a V-type refrigerator with air, H 2 , and He, and the instantaneous pressure and the COP were calculated. The results showed that the COP for the H 2 was higher than that for the He or air due to the low pressure drop and high specific heat. Cheng et al. [23] developed a β-type Stirling refrigerator that could reach 93 K with He at 1000 r/min and 3 atm. In addition, the total heat losses and the input power increased with the cooling temperature decreasing, while the COP decreased. Li et al. [24] built a free-piston Stirling refrigerator with an experimental cooling power of 350 W at 80 K and a relative Carnot efficiency of 26.8%. Smirnov et al. [25] developed an α-type Stirling refrigerator with a high-tolerance, lubricated piston-cylinder pair and a Ross-yoke transmission structure, finding that the loss of the ideal cooling capacity due to gas leakages amounted to 2.2%.
A Stirling engine usually works with a stable heating source of 500-800 • C to produce power, while a motor is required to drive a Stirling refrigerator in order to produce a cold source of usually less than −80 • C. Generally, the operation of a Stirling refrigerator is more convenient than that of a Stirling engine, whereby an operation can be easily regulated by a motor. It is probably a practical method to test the design and test characteristics of a Stirling engine run as a refrigerator, i.e., from a reverse cycle to a positive cycle.
However, most of the models were mainly used in a positive cycle or a reverse cycle. Most of the models of a Stirling refrigerator were more concerned with the acquisition and In the present work, a Simple model was improved to analyze the performance of both the Stirling engines and the refrigerators by modifying various power and heat losses. An experimental system based on a β-type Stirling prototype was built to study the positive and reverse Stirling cycles. The performance was further analyzed by experiment and model to understand the effects of the rotational speed, pressure, and working gases. The indicated power difference of the Stirling engine and the refrigerator is mainly caused by the upper and lower temperature limit, and it is possible to use the temperature ratio to establish the relationship between the positive and the reverse Stirling cycles, which means the performance of a Stirling engine can be predicted by a Stirling refrigerator. A novel method was developed to predict the performance of the positive Stirling cycles based on the reverse ones. Figure 1 illustrates a schematic diagram for a framework of an improved Simple model combined with various heat and power losses, which was developed based on an ISAM model [26,27] and can be applied to Stirling engines and refrigerators simultaneously. The Stirling-cycle analysis model mainly consists of four parts: the adiabatic analysis module, the loss analysis module, and the temperature and pressure check modules. Based on an ideal adiabatic model, the power/heat losses of engine components and the actual heat-transfer processes are fully considered. The control equations of the adiabatic analysis model [4] are shown in Appendix A. The heat/power losses considered in the model mainly include the heat conduction loss, regenerative heat loss, displacer shuttle heat loss, flow resistance loss of the heat exchanger, cross-section mutation power loss, gas spring hysteresis power loss, and seal leakage power loss.

Improved Simple Model Combined with Various Heat and Power Losses
(1) Heat conduction loss: heat conduction through the engine body.
(2) Regenerative heat loss: the temperature difference between the fillers and gas cannot be neglected for an actual regenerator, indicating that more heat is needed for a thermal cycle.
(3) Shuttle heat loss of the displacer: the displacer brings some heat from the hot chamber into the cold one, which means the heater should absorb more heat.
(4) Flow resistance power loss of the heat exchanger: pressure drops in the heat exchangers, especially in the regenerator, reduce the final power output.
(5) Cross-section mutation power loss: the velocity direction and pressure distribution of the gas flow, due to the sudden expansion or contraction of the pipe section, which makes a drop in the output power.
(6) Gas spring hysteresis power loss: for an ideal gas, the pressure/volume relationship is either isothermal or adiabatic. For real gas, this imperfect thermodynamic process dissipates some work.
The corresponding calculation equations of these losses are shown in Table 1. In the model, the actual mass of the working gas calculated by the Schmidt model was corrected, and the change of the physical parameters of the gas and the change of pressure and temperature were coupled. The convergence criterion is achieved when the temperature error is less than 1 • C and the relative error of the pressure is less than 0.1%. The smaller values (such as 0.1 • C and 0.01%) were adopted, and the convergence results were almost the same when calculating the convergence results. The above values (1 • C and 0.1%) were adopted in order to save calculation time.

Category Equation Description
Power losses flow resistance loss of heat exchanger (W) [31]  cross-section mutation power loss (W) [28,32] W gp : gas spring hysteresis power loss (W) γ: insulation factor T w : average wall temperature (K) P mean : average pressure (MPa) ∆V: volume amplitude (m 3 ) V B : mean volume of the gap spring cavity (m 3 ) A w : mean wetted area (m 2 ) seal leakage power loss (W) [33,34] W leak = q leak · C p · T leak q mleak = ρπD c As for the positive cycle, power losses directly affect the engine power whilst heat losses affect its ability to absorb heat directly, impacting the working gas temperature of the compression/expansion space and the cycle efficiency of the engine. However, for the reverse cycle, the power losses directly affect the input power, with heat losses affecting the actual cooling power of a cold head. This will also affect the temperature of the working gas and COP. The heat absorption/cooling capacity and cycle power equations of positive and reverse can be expressed as follows.
For the positive cycle: While for the reverse cycle: The heat transfer process is implemented to obtain the actual gas temperature in the hot/cold ends.
T gc = T wk − Q acco /h Aco (8) Nu = 0.023·Re·0.8·Pr n (heat : n = 0.4, cool : where Q ach is the actual heat input (W) of SE; Q adh is the adiabatic analysis heat input (W); Q acco is the actual cooling power (W); Q adco is the adiabatic analysis cooling power (W); W acipo is the actual cycle power (W) of SE; W adip is the adiabatic analysis cycle power (W); T ge /T gc is the working gas temperature in the expansion/compression space ( • C); T wh /T wk is the heater (heat absorber)/cooler wall temperature ( • C); A h /A co is the heater (heat absorber)/cooler convection heat transfer area (m 2 ); Q acc is the actual cooling power of SR (W); W acipi is the actual cycle input power (W) of SR; Q w is the heat conduction loss (W); Q rloss is the regenerative heat loss (W); Q sh is the displacer shuttle heat loss (W); W fr is the flow resistance loss of the heat exchanger (W); W fj is the cross-section mutation power loss (W); W gp is the gas spring hysteresis power loss (W); and W leak is the seal leakage power loss (W).

Experimental System
The performance of the Stirling engine and the refrigerator was studied to verify the accuracy and applicability of the improved Simple model. The three-dimensional structure of the engine and the measuring point arrangement are shown in Appendix A, with detailed parameters of the Stirling engine and the measuring instruments. Figure 2 shows the Stirling engine and the refrigerator-performance test system, consisting of three parts: a salt bath heating system, a β-type Stirling engine, and a data acquisition system [26]. The Stirling engine was heated by a salt bath whose temperature was maintained at 535 • C ± 2 • C. The internal working gas of a Stirling engine experiences a Stirling cycle, and the outputs work externally, driving a generator to create electricity. The Stirling refrigerator-performance test system also included three main parts: a servo motor, a refrigerator, and a data acquisition system. The measuring points of the Stirling engine and the refrigerator mainly included the gas temperature in the expansion/compression space, the gas pressure in the compression space, the temperature of the cooling water, the engine speed, and the crank angle of the flywheel.

Verification of the Model with GPU-3 Stirling Engine
The model was validated and analyzed by taking the GPU-3 Stirling engine, which was developed by General Motors [2,4], as a case. Some second-order analysis models [4,5,7,26,30] produced good results and exhibited good agreement with the experimental data when compared with the original Simple model proposed by Urieli and Berchowitz [4]. Timoumi's team [30] considered five kinds of energy losses, i.e., regenerator heat transfer loss, shuttle heat loss, flow resistance loss, gas spring hysteresis loss, and heat conduction loss. The finite speed thermodynamics principles were incorporated in the CAFS model [7]. The quasi-steady flow approach was adopted to analyze the heat-transfer and flow-friction effects of the heater, cooler, and regenerator on the performance of the engine in Fawad's optimization model [5]. Table 2 shows the model results compared with those of some previous methods under rated conditions (a heat source temperature of 977 K, a heat sink temperature of 288 K, an average pressure of 4.14 MPa, a rotational speed of 2500 r/min, and a total mass of working gas of 1.1362 g). Compared with the ISAM model [23], the cycle's average pressure was corrected, and the accuracy was greatly improved with the average pressure reduced from 4.26 MPa to 4.14 MPa and the mass from 1.255 g to 1.137 g. The model has a good applicability with a wide operating range. Figures 3 and 4 compare the indicated power and the cycle efficiency between the model values and the experimental ones of the GPU-3 Stirling engine when working at different conditions. The indicated power for H 2 increases, while the cycle efficiency increases at first and then decreases with the rise of rotational speed. The optimal rotational speed of H 2 is around 1500-2500 r/min. In general, the greater the pressure, the greater the corresponding indicated power and cycle efficiency. As to He, the indicated power increases at first and then gradually decreases with the increasing of the rotational speed, while the cycle efficiency shows an overall downward trend. This is due to the sharp increase in the flow resistance loss of the heat exchanger and the cross-section mutation power loss as the rotational speed increases. Indicatively, the predicted errors of the power and cycle efficiency are within 12% and 14%, respectively, when the rotational speed is in a range of 1000-3500 r/min with the H 2 at a constant pressure of 1.38 MPa and 2.76 MPa. The predicted errors of the indicated power and cycle efficiency are within 9% and 13%, respectively.

A Method to Predict the Performance of Positive Stirling Cycles Based on Reverse Stirling Cycles
The Stirling cycle, which consists of two isothermal processes and two isochoric processes, can be used as a prime mover, a refrigerator, and a heat pump depending on the direction of the cycle and the difference in the temperature of the hot and cold end. When used as the prime mover, the Stirling cycle goes along the clockwise direction and absorbs heat at a higher temperature. When used as a refrigerator, the direction of cycle is exactly the opposite. The refrigerator is driven by a motor and requires an external input of electric energy. Moreover, the cycle absorbs heat from the outside at a lower temperature. Figure 5 shows the p-V map and the T-S map of the ideal Stirling cycle, showing the cycle direction and temperature domain of the ideal Stirling cycle. The 1-2-3-4-1 positive cycle is the prime-mover cycle, and the 1-2-3'-4'-1 reverse cycle is the refrigerator cycle. Based on the ideal Stirling cycle, the positive and reverse cycles have the same compression work, but the expansion work of the positive cycle is much larger than the compression work, and the area of 1-2-3-4-1 is the indicated work of the Stirling engine. As to the refrigerator, the expansion work of the cycle is much smaller than the compression work, and the area of 1-2-3'-4'-1 is the indicated work input; that is the reason why the refrigerator needs a motor to drag. According to the ideal Stirling-cycle map, the cycle thermal efficiency of the prime-mover cycle and the coefficient of the performance of the refrigerator are shown in Equations (11) and (12), respectively. The efficiency of the Stirling cycle is the same as that of the Carnot cycle at the same temperature.
where η C is the cycle thermal efficiency of the Stirling engine, COP re f is the coefficient of the performance of the Stirling refrigerator, T C is the temperature of the gas in the compression space, T E is the temperature of the gas in the expansion space, and T ref is the refrigerator temperature. The internal working gas of the Stirling engine carries on the complex oscillating flow-heat and mass-transfer process, and in the applicability and accuracy of the Stirling cycle model similar problems also exist. Some problems of the performance evaluation of the Stirling engine caused by testing difficulties can be solved if the indicated work of the Stirling refrigerator can be used to predict the Stirling positive-cycle results. Figure 6 is the logic graph of the method to predict the performance of the positive Stirling cycles based on the reverse ones. The method mainly includes three parts: Carry out the positive and reverse Stirling cycles at the same low conditions according to the same Stirling engine working as both a prime mover and a refrigerator. II. Establish a mathematical model based on the positive and reverse Stirling cycles in part I. III. Carry out reverse Stirling cycles at high conditions and predict the performance of the positive Stirling cycles based on the mathematical model in part II. Specifically, the Stirling positive-and reverse-cycle experiments were carried out under the same low conditions (mean pressure and rotational speed) to obtain the upper and lower temperature limits of the Stirling cycle, including the gas temperature of the expansion and compression space as well as the cyclic p-V diagram, which can be used to calculate the indicated work. The most important part is to use mathematical models to establish the relationship between the work of the positive and reverse Stirling cycles. The mathematical model to predict the indicated power is shown in Equation (13). The difference between the cycle works of the positive and the reverse Stirling cycles is mainly affected by the cycle temperature difference in the case of controlling the same mean pressure and the rotational speed. Equation (13) is used to establish a dimensionless parameter based on temperature difference.
where W heat2 /W cool1 is the cycle power (W) of the positive/reverse Stirling cycle, T ge2 /T gc2 is the gas temperature in the expansion/compression space for the positive Stirling cycle, and T ge1 /T gc1 is the gas temperature (K) in the expansion/compression space for the reverse Stirling cycle, where A is the pressure term coefficient, related to the pressure. When A < 1, the control pressures of the positive/reverse Stirling cycles are equal, and the average temperature of the reverse Stirling cycle is lower, which means that the mass of corresponding working fluid is larger. B is the temperature term coefficient, which is assumed to be unchanged for the same Stirling engine. The Stirling engine cycle power results under the needed conditions are predicted by performing a refrigeration cycle under the required operating conditions and the indicated work relationship model.

Positive and Reverse Stirling Cycles' p-V Maps
As shown in Figure 7a, there is a "cusp phenomenon" in the p-V map obtained by the ideal adiabatic model, which is more remarkable in the experimental results (Figure 7b). The "cusp phenomenon" refers to the fact that the minimum volume part of the p-V map is sharper than the maximum volume part. Figure 8 shows the variation between the total volume of the expansion and compression spaces with the crank angle. The change rate of the volume near the minimum volume point is greater than that near the maximum volume point, which contributes to the "cusp phenomenon".   Figure 9a shows the p-V maps of the positive and reverse Stirling cycles at the same mean pressure. The distinction between the p-V maps is mainly reflected in the "pressure ratio" and "curve shape". The pressure ratio of the reverse Stirling cycle is significantly higher than that of the positive one, and the p-V map of the Stirling engine is asymmetric, while relatively harmonious, for the refrigerator. For the positive cycle, the expansion-space volume continues to increase while the compression-space volume decreases after the minimum volume point, and the temperature of the working gas increases overall; thus, the maximum pressure of the positive cycle is ahead of the minimum volume point. Similarly, the maximum pressure of the reverse cycle lags behind the minimum volume point, as shown in Figure 9b. The maximum pressure of the positive cycle is smaller than that of the reverse one, as shown in Figure 9b. In addition, the minimum pressure of the positive cycle is ahead of the minimum volume point; however, that of the reverse cycle lags behind the maximum volume point. The minimum pressure of the positive cycle is larger than that of the reverse one. The peak pressure occurs at the expansion progress for the positive cycle and the compression progress for the reverse cycle, while the trough pressure occurs at the compression progress for the positive cycle and the expansion progress for the reverse cycle. The effect of the minor loss caused by the sudden change of the cross-section area and flow resistance loss on the positive cycle is more remarkable than on the reverse cycle, which is a possible reason for the obvious "higher pressure ratio" phenomenon of the reverse cycle. Figure 10a,b shows the changes in the indicated work and power with the rotational speeds at a constant pressure of 2.84 MPa. Both the model and the experimental results demonstrate that the indicated work gradually decreases while the indicated power increases, with a relative error range of 6.2-11.3%, as the rotational speed increases with He. Figure 10c shows that the cycle efficiency for He increases with the rotational speed, while the slope of the increase gradually slows down. However, the cycle efficiency for N 2 and CO 2 increases at first and then decreases, indicating that there is an optimal rotational speed for a certain pressure. The value of the optimal rotational speed is ordered as He > N 2 > CO 2 and decreases with the increasing pressure, as shown in Figure 10d. The flow resistance loss of the heat exchangers and the regenerative heat loss of the regenerator increases sharply as the rotational speed increases, resulting in a decrease for the indicated power and a dramatic increase in the heat-absorption rate. The model's predicted errors are within 10.2%, 7.5%, and 13.7% for the He, N 2 , and CO 2 , respectively. The indicated work and cycle efficiency are mainly influenced by temperature ratio, pressure ratio, regenerative heat loss, flow resistance loss, and minor loss. The leakage loss is independent whilst the flow resistance loss increases exponentially along with the increase in the rotational speed, as shown in Figure 11; thus, the cycle power increases slowly and even decreases. The total flow resistance loss increases from 1.0 W to 18.9 W, from 4.42 W to 51.2 W, and from 8.26 to 63.1 W for the He, N 2 , and CO 2 at 2.84 MPa, respectively. The flow resistance for CO 2 is larger than that for He and N 2 , and the difference is more noticeable with the increase in rotational speed. The regenerator is the largest flow-resistant component in the Stirling engine and has a great influence on the cycle efficiency. About 50% of the power and heat losses are related to the regenerator [35]. Gheith et al. [36] proposed that the regenerator is a key component of the Stirling engines. Furthermore, the flow asymmetry should be considered for a deep analysis of the internal flow and temperature fields. The researcher put much focus on the inner structure of a regenerator, filling porosity, or material [37,38]. Here, we focus on the internal heat and mass transfer with different working gases, as shown in Table 3.

Performance of Stirling Refrigerator
The gas temperature of the expansion space reaches a relatively stable state after a certain period when the Stirling engine is connected with a servo motor. Both the model and the experimental results show that the cooling temperature for He decreases as the rotational speed increases, while the cooling temperature for N 2 decreases at first and then maintains a flat level with even a slight increase, as shown in Figure 12. The minimum cooling temperatures for He are smaller than those for N 2 at same mean pressures. The minimum cooling temperatures are −92.2 • C and −53.6 • C for He and N 2 at 2.84 MPa, respectively. The regenerative heat loss of He is much smaller than that of N 2 , which leads to a significant decrease in the cooling power of N 2 . Furthermore, there is an optimal rotational speed for the minimum cooling temperature, which is related to the type of working gases and pressures. It can be seen from Figure 12a that the optimal rotational speed for He is higher than that for N 2 and decreases as the mean pressure increases, as shown in Figure 12b. In addition, the cooling temperatures for He and N 2 are similar at lower rotational speeds (500 r/min), and the cooling temperature for He is lower compared with that of N 2 at higher rotational speeds, which indicates that He and other small molecule refrigerants have a greater refrigeration capacity in the Stirling refrigerator. Figure 13 shows the temperature variation of the gas in the compression space and the wall of cold head at different rotational speeds. The mass-flow rate increases as the rotational speed increases, resulting in an overall increase in the heat-exchange rate, and the temperature of the working gas in the compression space also increases. The minimum wall temperatures of the cold head are −87.2 • C and −46.1 • C for He and N 2 , respectively.  It can be found from Figure 14 that the pressure ratio of the p-V map for He increases with the rotational speed, while for the N 2 it is almost invariant, which is closely related to the temperature ratio. The pressure difference for the He increases with an increase in the rotational speed, while that for the N 2 increases first and then decreases, as shown in Figure 15a. The indicated work reduces from −4.64 J to −7.20 J, and the cycle input power increases from 38.6 W to 144.0 W for He when the rotational speed increases from 500 to 1200 r/min at 1.96 MPa, as shown in Figure 15b,c, while for the N 2 , it reduces from −3.59 J to −4.68 J (660 r/min) and then increases to −4.35 J when the rotational speed increases from 360 to 1200 r/min at 1.96 MPa. The cooling power for the He increases when the rotational speed increases, while the cooling power for the N 2 increases at first but then gradually decreases, as shown in Figure 15d. The predicted errors of the cycle input power and cooling power are within 10.9% and 14.5%, respectively, for the He and N 2 at different mean pressures.  It can be found from Figure 14 that the pressure ratio of the p-V map for He increases with the rotational speed, while for the N2 it is almost invariant, which is closely related to the temperature ratio. The pressure difference for the He increases with an increase in the rotational speed, while that for the N2 increases first and then decreases, as shown in Figure 15a. The indicated work reduces from −4.64 J to −7.20 J, and the cycle input power increases from 38.6 W to 144.0 W for He when the rotational speed increases from 500 to 1200 r/min at 1.96 MPa, as shown in Figure 15b,c, while for the N2, it reduces from −3.59 J to −4.68 J (660 r/min) and then increases to −4.35 J when the rotational speed increases from 360 to 1200 r/min at 1.96 MPa. The cooling power for the He increases when the rotational speed increases, while the cooling power for the N2 increases at first but then gradually decreases, as shown in Figure 15d. The predicted errors of the cycle input power and cooling power are within 10.9% and 14.5%, respectively, for the He and N2 at different mean pressures.  The decrease in cooling power is mainly caused by an increase in heat conduction loss and regenerative heat loss. Figure 16 shows the variation of the cycle input power and the heat/power losses in the model at different rotational speeds. With the increasing of the rotational speed, the displacer shuttle heat loss is practically unchanged, while the heat conduction loss and the regenerative heat loss increase significantly. The regenerative heat loss for N2 is greater for a lower thermal conductivity and a higher heat capacity compared with that of He. This indicates that the performance of the Stirling refrigerator is worse for N2 at a high rotational speed. The COP coefficient (defined as: COP = ) gradually decreases with the rotational speeds with a relative error The decrease in cooling power is mainly caused by an increase in heat conduction loss and regenerative heat loss. Figure 16 shows the variation of the cycle input power and the heat/power losses in the model at different rotational speeds. With the increasing of the rotational speed, the displacer shuttle heat loss is practically unchanged, while the heat conduction loss and the regenerative heat loss increase significantly. The regenerative heat loss for N 2 is greater for a lower thermal conductivity and a higher heat capacity compared with that of He. This indicates that the performance of the Stirling refrigerator is worse for ) gradually decreases with the rotational speeds with a relative error of <4.5%, as shown in Figure 17, which is a combined result of an increase in the cycle input power and an overall decrease in the cooling power.

Relationship between the Indicated Work of Stirling Engine and Refrigerator of the 100 W Stirling Engine
The partial experimental results of the Stirling engines and refrigerators at 1.96 MPa and 800 r/min are listed in Tables 4 and 5. According to the numerical model given by Equation (13), detailed relations between the cycle work of the 100 W beta-type Stirling engine and the refrigerator are obtained according to Figure 18, as shown in Equation (14).  Table 5. Stirling engines and refrigerators results with mean pressure (800 r/min, He).   Furthermore, combined with the characteristics under different pressures, the pressure coefficient A is obtained, referring to Table 6. The reason for the pressure coefficient A < 1 is due to the mass of the working gas difference on account of the difference of the temperature at the same mean pressure.  Figure 19 shows the experiment and model results for He at the constant pressure of 2.84 MPa. The error between the model and the experimental value is at the range of −0.5-7.4%. Based on the results of the reverse cycle, the cycle power of the Stirling engine is predicted relatively accurately (the cycle-power result error is the same as with the cycle-indicated work on account of the same rotational speed).

Conclusions
In order to develop a practical method to evaluate the performance of a Stirling engine by it running as a refrigerator, an improved Simple model was proposed at first for positive and reverse Stirling cycles simultaneously. A self-designed β-type Stirling engine was tested and the experimental results verified the good applicability and accuracy of the model. As to the positive Stirling cycles, the model errors of the indicated power and the cycle efficiency were within 11.3% and 13.7%, respectively. The indicated power maximum was 189.3 W, with a cycle efficiency of 14.7% at 1406 r/min for He. As to the reverse cycles, the predicted errors of the cycle input power and cooling power were within 10.9% and 14.5%, respectively. The experimental cooling temperatures could reach −92.2 • C and −53.6 • C, respectively, for He and N 2 at 2.8 MPa.
The p-V maps of the positive and reverse Stirling cycles were obtained at the same control pressures. The distinction between the positive and reverse Stirling cycles' p-V maps was mainly reflected in the "pressure ratio". The pressure ratio of the reverse cycle was higher than that of positive one at the same mean pressure. The peak pressure occurred at the expansion progress for the positive cycle and the compression progress for the reverse cycle, while the trough pressure occurred at the compression progress for the positive cycle and the expansion progress for the reverse cycle.
In addition, a method was proposed to predict the indicated work of the positive Stirling cycles based on the reverse ones. The most critical issue with this method was to establish an associated model of the gas temperatures in the expansion/compression space. A mathematical model to predict the indicated power of the positive and reverse Stirling cycles was proposed: The error between the values of the model and the experiment with He at 2.8 MPa was at the range of −0.5-7.4%. These results indicate that the model can provide much useful information for studying a Stirling engine and a refrigerator simultaneously and will probably bring a practical method for forecasting a positive cycle through a reverse one.   Appendix A Table A1. Equations of the ideal adiabatic model.
dm c = (pdV c + V c dp/γ)/(RT ck ) mass accumulations dm e = (pdV e + V e dp/γ)/(RT he ) dm k = m k dp/p dm r = m r dp/p dm h = m h dp/p m ck = −dm c mass flow m kr = m ck − dm k m he = dm e m rh = m he + dm h if m ck > 0 then T ck = T c , else T ck = T k conditional temperature if m he > 0 then T he = T h , else T he = T e dT c = T c (dp/p + dV c /V c − dm c /m c ) temperatures dT e = T e (dp/p + dV e /V e − dm e /m e ) Figure A1. β-type Stirling engine measuring point arrangement and 3D structure. 2-heating tube, 3-regenerator, 5-gas piston, 6-cooler, 7-buffer chamber, 8-power piston, 10-flywheel, 11-diamond transmission. A-working gas temperature of expansion space, B-working gas temperature of regenerator near expansion space, C-working gas temperature of regenerator near compression space, D-inlet cooling water temperature, E-cooling outlet water temperature, F-working gas temperature of compression space/working space pressure measuring point.       Figure A2. Equipment photograph.