Development of a Beta-Type Moderate-Temperature-Di ﬀ erential Stirling Engine Based on Computational and Experimental Methods

: Stirling engine is a favorable technique in the application of waste heat recovery or cogeneration system. This paper aims at developing a beta-type Stirling engine which is operated at moderate heating temperature (773–973 K). Rhombic drive mechanism is utilized to make coaxial motion of displacer and piston. Based on the proposed dimensions, a theoretical model combining thermodynamic and dynamic analysis is built to predict the performance of the Stirling engine. Thermodynamic analysis deals with variations of properties in each chamber while dynamic analysis handles the resultant shaft torque produced by the Stirling engine. Furthermore, a prototype engine is manufactured, and experimental test is carried out to validate the simulated results in this research. Under heating temperature of 973 K, charged pressure of 8 bar, rotation speed of 1944 rpm, shaft power of 68 W is obtained from the prototype Stirling engine. Power density is calculated to be 1.889 W-c.c. − 1 by theoretical prediction and 1.725 W-c.c. − 1 by tested result. The impact of the geometrical dimensions is investigated to survey the optimal piston diameter which is related to compression ratio and swept volume.


Introduction
As an external heat engine, the Stirling engine features availability to most heat sources, such as solar, biomass, geothermal, nuclear, and waste heat energy. In recent years, capacity and efficiency of Stirling engines have been remarkably lifted and lots of effort were paid to implement the practical utilization of this technology [1]. Atmosphere gas like air, nitrogen, helium, and hydrogen can be adopted as working medium in Stirling engines [2,3]. Three main kinds of structure are exploited in Stirling engines, which are alpha-, beta-, and gamma-type. According to first-order analysis on design parameters, beta-type Stirling engine has the highest dimensionless work output [4]. Coaxial movement of piston and displacer is the characteristic of beta-type Stirling engine, which results in a more compact configuration and better power density. Cinar et al. [5] designed an air-charged beta-type Stirling engine with crank mechanism. At hot-source temperature of 1000 • C generated from a radioactive heater, the test engine can generate 5.98 W at 208 rpm. Aksoy et al. [6] adopted a slotted lever to implement the displacement of displacer and piston. The engine was tested with halogen lamp and power of 127.17 W was achieved at heating temperature of 873 K. In contrast to crank mechanism, Meijer [7] demonstrated a rhombic drive mechanism firstly which successfully operated to achieve the scale of 3 W-8.95 kW. The advantage of rhombic drive mechanism is that piston and displacer bear less lateral force due to the symmetrical geometry, and hence vibration and noise are remarkably reduced. Erol and Caliskan [8] compared performance of four mechanisms of beta-type Stirling engines with isothermal and kinematic analysis. At the same swept volume of 365 cm 3 and charged mass of 7.16 × 10 −4 kg, net

Theoretical Model
In beta-type Stirling engine, a displacer and a power piston are used to generate volume variations in expansion and compression chamber with a certain phase angle. Expansion and compression chamber are set at high and low temperature regions, respectively, and working gas will shuttle between these two chambers. Figure 1 shows the schematic of studied beta-type Stirling engine where three heat exchangers are also incorporated. A heater and cooler are introduced to enhance the heat transfer in hot and cold ends. They are composed of numerous tubes and immersed in the heat sinks. Between the heater and cooler, regenerator is an internal heat exchanger that serves as thermal inertia and provides a great amount of heat transfer area. The regenerator is usually composed of porous material and benefits maintaining the axial temperature gradient. While high temperature working gas passes through the regenerator, its internal energy will be preserved inside the material. The energy is afterward recycled as working gas returns from the compression chamber. Stainless steel (SS304) mesh that is composed of crossed metal wire is employed as regenerator for the present Stirling engine. Below power piston, the back chamber functions as a buffer for the working zone and it is taken isolated in the present model. Rhombic drive mechanism is adopted to drive a pair of gears where the Energies 2020, 13, 6029 3 of 14 right gear serves as the main rotor to endure shaft load. In the figure, y p and y d are displacements of the piston and displacer while θ is the crank angle about the main shaft. D d , D p , and D b are the diameters of displacer, piston, and back chamber rod, respectively. L ce and L cc are clearance lengths of expansion and compression chamber.
Energies 2020, 13, x FOR PEER REVIEW 3 of 14 the present Stirling engine. Below power piston, the back chamber functions as a buffer for the working zone and it is taken isolated in the present model. Rhombic drive mechanism is adopted to drive a pair of gears where the right gear serves as the main rotor to endure shaft load. In the figure, In this paper, a theoretical model is used to simulate the transient process via thermodynamic analysis for working gas and dynamic analysis for mechanical components. As indicated on the schematic, the working zone consists of expansion, heater, regenerator, cooler, and compression chambers. In thermodynamic analysis, properties including volume, mass, pressure, and temperature in each chamber are calculated by ideal-gas, continuity, momentum, and energy equations. Volume variations of expansion and compression chamber are calculated once the displacements of piston and displacer are updated with crank angle. After that, average pressure inside engine ave p is obtained and mass in each chamber can be evaluated as / In the above equation, mi, Vi, and Ti are mass, volume, and temperature in each chamber. Rconst is gas constant whose value is 2077 J-kg −1 -K −1 in the present case that helium is used. Subscripts e, h, r, k, c stand for expansion, heater, regenerator, cooler, and compression chamber.
Mass flow rate across the boundaries can be evaluated by continuity equation. As working gas flows through heater, regenerator, and cooler, it suffers from viscous friction on the channel surface. Therefore, pressure drops in these three elements should be considered. Friction factors f regarding Reynolds number can be described as [22] 64 / Re where Reynolds number is determined in terms of the hydraulic diameter. As a result, pressure drop across three heat exchangers can be estimated with the help of Darcy-Weisbach equation. In present model, average pressure is set at the regenerator, and the pressures in other chambers are determined as In this paper, a theoretical model is used to simulate the transient process via thermodynamic analysis for working gas and dynamic analysis for mechanical components. As indicated on the schematic, the working zone consists of expansion, heater, regenerator, cooler, and compression chambers. In thermodynamic analysis, properties including volume, mass, pressure, and temperature in each chamber are calculated by ideal-gas, continuity, momentum, and energy equations. Volume variations of expansion and compression chamber are calculated once the displacements of piston and displacer are updated with crank angle. After that, average pressure inside engine p ave is obtained and mass in each chamber can be evaluated as In the above equation, mi, Vi, and Ti are mass, volume, and temperature in each chamber. R const is gas constant whose value is 2077 J-kg −1 -K −1 in the present case that helium is used. Subscripts e, h, r, k, c stand for expansion, heater, regenerator, cooler, and compression chamber.
Mass flow rate across the boundaries can be evaluated by continuity equation. As working gas flows through heater, regenerator, and cooler, it suffers from viscous friction on the channel surface. Therefore, pressure drops in these three elements should be considered. Friction factors f regarding Reynolds number can be described as [22] where Reynolds number is determined in terms of the hydraulic diameter. As a result, pressure drop across three heat exchangers can be estimated with the help of Darcy-Weisbach equation. In present model, average pressure is set at the regenerator, and the pressures in other chambers are determined as Energies 2020, 13, 6029 4 of 14 p r = p ave (6) p k = p r − (∆ p r + ∆p k )/2 (7) p c = p k − ∆p k /2 (8) where ∆p h , ∆p r , and ∆p k are pressure of heater, regenerator, and cooler chamber, respectively. For the back chamber, it is assumed that there is no mass leakage in the chamber and its temperature remains constant. Thus, the pressure in the back chamber (p b ) can be evaluated with ideal-gas equation as where back chamber volume V b is related to the position of piston. In order to determine temperature variation of the working gas, the energy equation of each chamber is used. That is In the above equation, .
Q is the heat transfer input and .
W is the work output. Work output is evaluated by the boundary work in the expansion and compression chambers. Since volumes of heater, regenerator, and cooler chamber remain unchanged, no boundary work is exerted. Heat transfer input depends on thermal resistance between the working gas and the heat sink. For the working gas flowing in the channel of heater, regenerator, and cooler, convection heat transfer dominates. The convective heat transfer coefficient can be determined from Nusselt number [23] as In thermodynamic analysis, volume, mass, pressure, and temperature are calculated. Subsequently, pressure of working zone is brought into dynamic analysis for predicting the resultant shaft torque of the engine. Yang, et al. [24] presented the dynamic relations of rhombic drive mechanism for a 1-kW Stirling engine, in which friction of each joint is considered. The friction caused by the ball bearing is relatively small compared to the dry friction at the piston-cylinder interface and the resistance of the shaft seal. Besides, gravitational and inertial energies of components are demonstrated to be trivial compared to the driving power as well. Cheng et al. [25] conducted experimental measurements for the mechanical loss of rhombic drive mechanism. It is found that the overall friction power can be represented in polynomial form as a function of the rotation speed. In this study, the analysis in [24] is simplified by neglecting the effects of joint friction, gravitational, and inertial energies. All the frictions involved in the engine are summarized in terms of angular speed. The angular acceleration α of the main shaft can be evaluated in the following equation once the pressures in the expansion, compression, and back chambers are determined. . .
In the above equations, . W dri is the driving power induced by the gas pressures of expansion, compression, and back chambers; . W f ri is the frictional power represented in a second-order polynomial . W sh is the shaft power defined by the product of exerted shaft torque τ sh and angular speed ω. v p and v d are linear velocities of the piston and displacer. In addition, c 1 and c 2 are two coefficients determined according to the experimental data. For the present case, c 1 = 2.025 × 10 −1 W-s, and c 2 = 1.125 × 10 −4 W-s 2 . Figure 2 shows the flow chart of the computation with the proposed theoretical model. An initial rotation speed is applied and after a starting period the engine is able to reach the steady operation state. At each time step, thermodynamic and dynamic analyses are performed simultaneously and iteratively until the properties, m, p, T, and α converge to the solutions. Table 1 lists the design parameters of the studied Stirling engine.
In the above equations, dri W  is the driving power induced by the gas pressures of expansion, compression, and back chambers; fri W  is the frictional power represented in a second-order polynomial form; sh W  is the shaft power defined by the product of exerted shaft torque sh τ and angular speed ω. p v and d v are linear velocities of the piston and displacer. In addition, c1 and c2 are two coefficients determined according to the experimental data. For the present case, c1=2.025 × 10 −1 W-s, and c2=1.125 × 10 −4 W-s 2 . Figure 2 shows the flow chart of the computation with the proposed theoretical model. An initial rotation speed is applied and after a starting period the engine is able to reach the steady operation state. At each time step, thermodynamic and dynamic analyses are performed simultaneously and iteratively until the properties, m, p, T, and α converge to the solutions. Table 1 lists the design parameters of the studied Stirling engine.

Prototype Engine and Experimental Setup
In this paper, a prototype beta-type Stirling engine based on the above-mentioned parameters was manufactured to validate the raised model. Figure 3a shows the picture of the designed Stirling engine prototype. Components of the Stirling engine require high precision machining and accurate assembly. The piston and cylinder were made with fine surface roughness and matched in tight tolerance. Hermetic seals were applied at the junctions of parts to avoid gas leakage. Copper pipe was used as the heat exchanger in the heating zone for its good thermal conductivity.
For the purpose of acquiring performance of the designed Stirling engine, a setup of experimental apparatus shown in Figure 3b was constructed. Kuehl [26] put forward empirical equations for estimating the transport properties of gases. Dynamic viscosity and thermal conductivity of the available working gases in Stirling engine are listed in Table 2, where the values were calculated based on the equations in [26] and the conditions were set at 873 K and 6 bar based on Table 1. It can be found that dynamic viscosity of hydrogen is roughly half of air, nitrogen, and helium. In terms of thermal conductivity, helium and hydrogen perform 6-13 times better than air and nitrogen. Despite the fact that conductivity of helium is 44 % of hydrogen, helium was chosen as the working gas due to its availability and safety. Interior air should be drawn out firstly by vacuum pump and high purity (99.999%) helium is then filled into the engine. Charged pressure was controlled by a pressure regulator whose operational range is 0-10 bar and resolution is 0.01 bar. In experiments, an electric heater was used to supply thermal energy for the engine, and the heating temperature was adjusted with a PID controller. The engine shaft was connected to a torque and speed sensor and a hysteresis brake. In this manner, the shaft power can be determined in terms of the torque and rotation speed measured. Variation of heating temperature, engine torque, and rotation speed were collected and recorded by a data acquisition system.

Prototype Engine and Experimental Setup
In this paper, a prototype beta-type Stirling engine based on the above-mentioned parameters was manufactured to validate the raised model. Figure 3a shows the picture of the designed Stirling engine prototype. Components of the Stirling engine require high precision machining and accurate assembly. The piston and cylinder were made with fine surface roughness and matched in tight tolerance. Hermetic seals were applied at the junctions of parts to avoid gas leakage. Copper pipe was used as the heat exchanger in the heating zone for its good thermal conductivity.
For the purpose of acquiring performance of the designed Stirling engine, a setup of experimental apparatus shown in Figure 3b was constructed. Kuehl [26] put forward empirical equations for estimating the transport properties of gases. Dynamic viscosity and thermal conductivity of the available working gases in Stirling engine are listed in Table 2, where the values were calculated based on the equations in [26] and the conditions were set at 873 K and 6 bar based on Table 1. It can be found that dynamic viscosity of hydrogen is roughly half of air, nitrogen, and helium. In terms of thermal conductivity, helium and hydrogen perform 6-13 times better than air and nitrogen. Despite the fact that conductivity of helium is 44 % of hydrogen, helium was chosen as the working gas due to its availability and safety. Interior air should be drawn out firstly by vacuum pump and high purity (99.999%) helium is then filled into the engine. Charged pressure was controlled by a pressure regulator whose operational range is 0-10 bar and resolution is 0.01 bar. In experiments, an electric heater was used to supply thermal energy for the engine, and the heating temperature was adjusted with a PID controller. The engine shaft was connected to a torque and speed sensor and a hysteresis brake. In this manner, the shaft power can be determined in terms of the torque and rotation speed measured. Variation of heating temperature, engine torque, and rotation speed were collected and recorded by a data acquisition system.
The major parameters measured in the experiment were charged pressure, heating temperature, shaft torque, and rotation speed. Uncertainty analysis proposed by Moffat [27] is performed in Table  3 where the parameters are estimated with 95% confidence. Relative uncertainty of shaft power is calculated to be 3.7%.    The major parameters measured in the experiment were charged pressure, heating temperature, shaft torque, and rotation speed. Uncertainty analysis proposed by Moffat [27] is performed in Table 3 where the parameters are estimated with 95% confidence. Relative uncertainty of shaft power is calculated to be 3.7%.

Results and Discussion
In the computation, once the steady operation was reached, an increment of shaft torque was added. The shaft torque was increased to alter the rotation speed as well as the shaft power. Figure 4 displays the dynamic simulation by plotting the variations in the average rotation speed and the shaft power in response to shaft torque change. The charged pressure is 4 bar and the heating temperature is 723 K. As shown in Figure 4a, the Stirling engine is accelerated from the initial condition without shaft torque and it takes about 64 s to attain the steady operation regime. While the resisting shaft torque is applied, engine speed will be gradually declined to another level. The variation in the shaft power is plotted in Figure 4b. In this figure, data of the friction power, driving power, and net heat input are also provided. The net heat input of the Stirling engine, . Q net , is the sum of heat transfer rates in all chambers. It is clearly observed that larger heat input is required into the engine in the first few seconds. At each steady level, net heat input rate is balanced by the driving power. Note that the driving power is consumed by the friction power and the shaft power in the steady operation regime with α = 0, as indicated in Equation (13).

Results and Discussion
In the computation, once the steady operation was reached, an increment of shaft torque was added. The shaft torque was increased to alter the rotation speed as well as the shaft power. Figure 4 displays the dynamic simulation by plotting the variations in the average rotation speed and the shaft power in response to shaft torque change. The charged pressure is 4 bar and the heating temperature is 723 K. As shown in Figure 4a, the Stirling engine is accelerated from the initial condition without shaft torque and it takes about 64 s to attain the steady operation regime. While the resisting shaft torque is applied, engine speed will be gradually declined to another level. The variation in the shaft power is plotted in Figure 4b. In this figure, data of the friction power, driving power, and net heat input are also provided. The net heat input of the Stirling engine, net Q  , is the sum of heat transfer rates in all chambers. It is clearly observed that larger heat input is required into the engine in the first few seconds. At each steady level, net heat input rate is balanced by the driving power. Note that the driving power is consumed by the friction power and the shaft power in the steady operation regime with α = 0, as indicated in Equation (13).  Performance curves of the designed Stirling engine under different charged pressures and heating temperatures are illustrated in Figure 5, where both numerical and experimental results are provided. For the present engine the heating temperature varied from 500 to 700 °C (i.e., 773-973 K). Thus, the difference between heating and cooling temperatures ranges within 473-673 K. For 4-bar charged pressure shown in Figure 5a, the engine can be actuated at heating temperature of 773 K and peak shaft power of the prototype is 16 W at 1108 rpm. As the heating temperature is increased to 973 K, shaft power is lifted to 47 W in numerical simulation and 44 W in measurement. For 6-bar charged pressure, the operated heating temperature is between 823 and 973 K as displayed in Figure  5b. Shaft power ranges from 36 to 64 W in theoretical model and from 29 to 60 W in experimental data. In case of 8-bar charged pressure, the lowest heating temperature for starting engine is observed to be 873 K in testing. It can be found that the lowest heating temperature depends on the charged pressure of the designed engine. Since more working gas is filled in the engine at larger charged pressure, higher heating temperature is required to supply sufficient heat energy to actuate the engine. At 8-bar charged pressure and 973-K heating temperature, shaft power is predicted to be 74 W in simulation and verified to be 68 W in experiment. Performance curves of the designed Stirling engine under different charged pressures and heating temperatures are illustrated in Figure 5, where both numerical and experimental results are provided. For the present engine the heating temperature varied from 500 to 700 • C (i.e., 773-973 K). Thus, the difference between heating and cooling temperatures ranges within 473-673 K. For 4-bar charged pressure shown in Figure 5a, the engine can be actuated at heating temperature of 773 K and peak shaft power of the prototype is 16 W at 1108 rpm. As the heating temperature is increased to 973 K, shaft power is lifted to 47 W in numerical simulation and 44 W in measurement. For 6-bar charged pressure, the operated heating temperature is between 823 and 973 K as displayed in Figure 5b. Shaft power ranges from 36 to 64 W in theoretical model and from 29 to 60 W in experimental data. In case of 8-bar charged pressure, the lowest heating temperature for starting engine is observed to be 873 K in testing. It can be found that the lowest heating temperature depends on the charged pressure of the designed engine. Since more working gas is filled in the engine at larger charged pressure, higher heating temperature is required to supply sufficient heat energy to actuate the engine. At 8-bar charged pressure and 973-K heating temperature, shaft power is predicted to be 74 W in simulation and verified to be 68 W in experiment.  Performance curves of the designed Stirling engine under different charged pressures and heating temperatures are illustrated in Figure 5, where both numerical and experimental results are provided. For the present engine the heating temperature varied from 500 to 700 °C (i.e., 773-973 K). Thus, the difference between heating and cooling temperatures ranges within 473-673 K. For 4-bar charged pressure shown in Figure 5a, the engine can be actuated at heating temperature of 773 K and peak shaft power of the prototype is 16 W at 1108 rpm. As the heating temperature is increased to 973 K, shaft power is lifted to 47 W in numerical simulation and 44 W in measurement. For 6-bar charged pressure, the operated heating temperature is between 823 and 973 K as displayed in Figure  5b. Shaft power ranges from 36 to 64 W in theoretical model and from 29 to 60 W in experimental data. In case of 8-bar charged pressure, the lowest heating temperature for starting engine is observed to be 873 K in testing. It can be found that the lowest heating temperature depends on the charged pressure of the designed engine. Since more working gas is filled in the engine at larger charged pressure, higher heating temperature is required to supply sufficient heat energy to actuate the engine. At 8-bar charged pressure and 973-K heating temperature, shaft power is predicted to be 74 W in simulation and verified to be 68 W in experiment. Clearance length is defined to describe the least distance that causes dead volume for a chamber. In the present design of beta-type Stirling engine, clearance lengths of expansion and compression chambers are 0.003 and 0.0165 m, respectively. If the dead volume of the compression chamber can be reduced, the compression ratio is increased and engine performance can be improved due to the more efficient compression of working gas. Figure 6 shows the effects of the clearance length and piston diameter on the engine performance under 8-bar charged pressure and 973-K heating temperature. In Figure 6a, clearance length of the compression chamber is decreased from 0.0165 to 0.0005 m, which makes an increase of compression ratio from 1.242 to 1.300. As a result, shaft power is raised from 74 to 82 W, indicating 11 % enhancement.
Influence of the piston diameter is displayed in Figure 6b. Compression ratio will be elevated to 1.463 while the diameter is increased to 0.07 m. It is noted that performance curve is narrowed, and the maximum rotation speed is lowered as the diameter gets larger. Stirling engine is accelerated by the driving torque which comes from the normal pressure force on the piston. However, pressure of back chamber which also exerts on the cross area of the piston will balance the resultant force on the piston. As the cross section is larger, the force induced by the back chamber is increased under the same charged pressure, which drops the maximum speed of the engine. It can be seen that shaft Clearance length is defined to describe the least distance that causes dead volume for a chamber. In the present design of beta-type Stirling engine, clearance lengths of expansion and compression chambers are 0.003 and 0.0165 m, respectively. If the dead volume of the compression chamber can be reduced, the compression ratio is increased and engine performance can be improved due to the more efficient compression of working gas. Figure 6 shows the effects of the clearance length and piston diameter on the engine performance under 8-bar charged pressure and 973-K heating temperature. In Figure 6a, clearance length of the compression chamber is decreased from 0.0165 to 0.0005 m, which makes an increase of compression ratio from 1.242 to 1.300. As a result, shaft power is raised from 74 to 82 W, indicating 11 % enhancement.
Influence of the piston diameter is displayed in Figure 6b. Compression ratio will be elevated to 1.463 while the diameter is increased to 0.07 m. It is noted that performance curve is narrowed, and the maximum rotation speed is lowered as the diameter gets larger. Stirling engine is accelerated by the driving torque which comes from the normal pressure force on the piston. However, pressure of back chamber which also exerts on the cross area of the piston will balance the resultant force on the piston. As the cross section is larger, the force induced by the back chamber is increased under the same charged pressure, which drops the maximum speed of the engine. It can be seen that shaft power is lifted to 85 W in the case of 0.06 m piston diameter, where 15 % power improvement is obtained. For 0.065 m piston diameter, shaft power will slightly decrease with a lower speed. When the piston diameter is 0.07 m, shaft power declines to 80 W.  Figure 7 shows the effect of piston diameter on the compression ratio and power density under 8-bar charged pressure and 973-K heating temperature. Compression ratio is increased with larger piston diameter, as plotted in the figure. Power density is defined as the ratio of peak shaft power and swept volume of the piston. In the proposed Stirling engine, power density is predicted to be 1.889 W-c.c. −1 and tested to be 1.725 W-c.c. −1 . On the other hand, it is noticed from the figure that an optimal value of power density exists. Highest power density is 2.020 W-c.c. −1 in the case of 0.042 m piston diameter. It can be deduced that piston diameter of 0.06 m leads to the largest shaft power and piston diameter of 0.042 m yields the best power density. For the purpose of more comprehensive analysis of the Stirling engine, design parameters including piston stroke, displacer diameter, displacer stroke, and phase angle should be studied simultaneously and optimized through optimization method.  Figure 7 shows the effect of piston diameter on the compression ratio and power density under 8-bar charged pressure and 973-K heating temperature. Compression ratio is increased with larger piston diameter, as plotted in the figure. Power density is defined as the ratio of peak shaft power and swept volume of the piston. In the proposed Stirling engine, power density is predicted to be 1.889 W-c.c. −1 and tested to be 1.725 W-c.c. −1 . On the other hand, it is noticed from the figure that an optimal value of power density exists. Highest power density is 2.020 W-c.c. −1 in the case of 0.042 m piston diameter. It can be deduced that piston diameter of 0.06 m leads to the largest shaft power and piston diameter of 0.042 m yields the best power density. For the purpose of more comprehensive analysis of the Stirling engine, design parameters including piston stroke, displacer diameter, displacer stroke, and phase angle should be studied simultaneously and optimized through optimization method. A comparison of developed beta-type Stirling is shown in Table 4. It is noted that the present Stirling engine yields higher power density than those presented in [6] and [9] at nearly equal pressure and temperature. In comparison to a high-pressure system (29.6 bar) demonstrated in [10], the present engine is lower only by 0.075 W-c.c. −1 .
Note that the present study aimed to develop a moderate-temperature-differential Stirling engine. From design point of view, the β-type Stirling engines are well suitable for development of the compact engines. Therefore, both modeling and prototyping phases were combined in this study so as to fulfill the aim. The compact property of the engine is a criterion to be competitive in personal power generation applications.

Conclusions
In this study, a beta-type Stirling engine for heating temperature from 773 to 973 K was developed. Dynamic and thermodynamic analyses were built to simulate the characteristic of the engine behavior. A prototype Stirling engine was manufactured, and experiments of different charged pressure and heating temperature were performed. Conclusions obtained from this research are listed as below.
(1) Dynamic behavior of the Stirling engine in response to varied shaft torque was demonstrated in the paper. Time-dependent variations of rotation speed and energies were illustrated. Conservation of energy was also satisfied to certify the presented method.
(2) At 4-bar charged pressure, the prototype Stirling engine generates 16-44 W shaft power in the temperature span of 773-973 K. Meanwhile, numerical results according to the proposed model show 17-47 W within the same temperature. As the working conditions are elevated to 8-bar charged pressure and 973-K heating temperature, the shaft power of the designed engine is able to achieve 68 W while predicted power is 74 W. Close agreements between the experimental results and the numerical predictions were found. (3) By means of parametric study on the compression chamber, shaft power is increased with smaller clearance. It was observed that optimal piston diameters corresponding to maximum