An Analysis Model Combining Gamma-Type Stirling Engine and Power Converter

Waste heat is a potential source for powering our living environment. It can be harvested and transformed into electricity. Ohmic heat is a common type of waste heat. However, waste heat has the following limitations: wide distribution, insufficient temperature difference (∆T < 70 K) for triggering turbines, and producing voltage below the open voltage of the battery. This paper proposes an energy harvester model that combines a gamma-type Stirling engine and variable capacitance. The energy harvester model is different from Tavakolpour-Saleh’s free-piston-type engine [7.1 W at ∆T = 407 K (273–680 K)]. The gamma-type Stirling engine is a low-temperature-difference engine. It can be triggered by a minimum ∆T value of 12 K (293–305 K). The triggering force in the variable capacitance is almost zero. Furthermore, the gamma-type Stirling engine is suitable for harvesting waste heat at room temperature. This study indicates that 21 mW of energy can be produced at ∆T = 30 K (293–323 K) for a bias voltage of 70 V and volume of 103.25 cc. Because of the given bias voltage, the energy harvester can break through the open voltage of the battery to achieve energy storage at a low temperature difference.


Introduction
Many forms of waste heat are produced in our living environment, such as the heat of combustion from factories, ohmic heating of computers, and heat from compressors.These heat sources increase the ambient temperature and cause global warming if we not do suitably handle them.Common methods of waste heat recovery include turbocharging, turbocompounding, thermoelectrics, and cabin heating and cooling [1].Therefore, waste heat should be harvested and stored as electric power.The most suitable method for achieving this objective involves thermoelectric material use.However, the costs of semiconductor thermoelectric material are too high for typical thermoelectric power generation applications at ∆T < 135 • C [2].This study aims to convert heat into electric power.To achieve this target, there exist three procedures.The first procedure involves converting heat into mechanical energy, the second involves converting mechanical energy into electric energy, and the third involves a combination of first two procedures.The three procedures are discussed in the following text.
The Stirling engine, which was created by Robert Stirling in 1816, has high efficiency for transforming heat into mechanical energy [3].The Stirling cycle involves the cyclic expansion and compression of the working fluid by maintaining a temperature difference between the hot and cold plates.The heat energy produced can be converted to mechanical work.Stirling engines can be classified into three types: alpha, beta, and gamma machines [4].The alpha-type Stirling engine has two pistons for compression and expansion.Both pistons are connected to the same crankshaft.The beta-type Stirling engine has one piston and one displacer, which share the same air room.The piston and displacer are linked by a wheel with two crankshafts.In the gamma-type engine, the piston and displacer are housed within two air rooms.The gamma-type engine (low temperature difference (LTD) engine) is the most suitable option for triggering the Stirling engine with a temperature difference from the environment (∆T = 10-30 • C) [4].The gamma-type engine has a high rotating speed, is triggered at LTDs, and has a small torque.The rotating speed is inversely proportional to the torque [5].
To calculate the performance of the Stirling engine, researchers developed computer fluid dynamic modes in commercial software [6].Robson et al. created a dynamic model which forms with each motion of component [7].Egas and Clucas compared the net work values of five types of Stirling engines by using the thermodynamics model [4].Abovementioned studies have assumed that the instantaneous pressure in a Stirling engine is constant.However, a variety of instantaneous pressures should be considered in reality.Kim et al. used thermal flow to analyze high and low pressures in an LTD Stirling engine [8].
For transforming mechanical energy to electrical energy, a conversion method that affects the mechanism of harvester should be selected.Vibration energy harvesters (VEHs) can be classified into three types according to the principle of electron transfer: (i) electrostatic, (ii) electromagnetic, and (iii) piezoelectric VEHs.The electromagnetic method is unsuitable for microelectromechanical systems because the magnetic coil requires complicated manufacturing processes [9].The piezoelectric material is highly brittle and hence is limited to a small deflection range [10].In the electrostatic method, a bias voltage is provided between two metal plates so that they become a capacitance component.The capacitance component then releases current from the elevated potential energy through vibration [11].The electrostatic method is more convenient than the electromagnetic and piezoelectric methods for harvesting energy from the environment.The capacitance component can be manufactured in any geometric shape and size to maximize the capacitive value.The component is widely used to increase the overlap area between two electrodes in the design concept of the interdigitated electrode [12,13].Moreover, the designer adjusts the resonance frequency of the VEH according to the vibration source [14].In addition to changing the gap of the variable capacitance, Moon et al. build a reverse electrowetting device with changing overlap area [15].Hsu et al. proposed the bubbler structure, which causes air bubbles to move into the gap between two electrodes.The entry of air bubbles changes the dielectric coefficient of the variable capacitance [16].
Energy harvesting and storage can be achieved at room temperature.Some researchers select the vibrations of pavements, which are caused by human walking and car movement [17,18].The frequency of human walking is 1.41-2.13Hz [19], and the frequency of bridge vibration is 0-30 Hz [20].The vibration source should have a high frequency, provide continuous waves, and have a stable amplitude for VEHs.The vibration source above has a low frequency (resonance frequency of the capacitive VEH is 6-13 GHz [21]), is discontinuous, and is unstable.Consequently, the vibration source is unsuitable for the electrostatic method.Therefore, the gamma-type Stirling engine can provide a stable and high-frequency source for capacitive VEHs.The bias voltage from the harvester can break through the open circuit voltage of the battery (open circuit voltage of lithium-ion batteries is 2.7-3.6 V during the charging process [22]).The Stirling generator has many applications.To accelerate the free-piston Stirling engine for reaching the resonance state (9.2 Hz), Tavakolpour-Saleh et al. [23] installed an additional DC motor on the Stirling engine.The obtained output power was 7.1 W at ∆T = 407 K (273-680 K).The combined heat and power system is an efficient energy recovery system.It can generate electricity and useful thermal energy and provide warm bath water, warm air, and electricity (10-36 kW) in the house [24,25].The dish Stirling generator can obtain solar energy through a concentrator and reflector.The range of output obtained with this generator is 2-50 kW.The dish Stirling generator is built in the United States, Japan, and Russia [26].
A mathematical model is proposed for evaluating the performance of the gamma-type Stirling engine with a capacitive harvester (GSCH).The governing equation of the energy harvester is constituted according to the conservation of energy.The governing equation includes the heat transfer, dynamic model, and electric model.If a temperature difference is applied between the hot plate and cold plate in the GSCH (input), the heat energy from the hot plate enters the work system of the Stirling engine.The energy is assigned into three parts, namely the mechanical, air, and friction parts.The mechanical part provides mechanical energy to generate power (output).According to the simulation results, the wheel of the energy harvester has a rotation frequency of 3.35 Hz at ∆T = 30 K (298-328 K) and can produce 21 mW electric power under unit area at a bias voltage of 70 V.Moreover, the optimized size of the engine can be determined by changing the geometric shape and size of the gamma-type Stirling engine, which has maximum conversion efficiency in a specific temperature range.

Heat Transfer
We assume that the temperature function of the GSCH only depends on the distance of the z-axis.The heat transfer model can be simplified by considering an ideal cylinder.Consider that the thermal conductivity (ignoring convection and radiation) and temperature have a linear relationship.Therefore, the heat transfer rate is q (Figure 1).The heat transmitted through the hot plate, air room, and cold plate can be written as follows: where ∆z Hot , ∆z Air , and ∆z Hot are the thicknesses of the hot plate, air room, and cold plate, respectively, and A Dc the cross-sectional area of the air room. of the Stirling engine.The energy is assigned into three parts, namely the mechanical, air, and friction parts.The mechanical part provides mechanical energy to generate power (output).
According to the simulation results, the wheel of the energy harvester has a rotation frequency of 3.35 Hz at ΔT = 30 K (298-328 K) and can produce 21 mW electric power under unit area at a bias voltage of 70 V.Moreover, the optimized size of the engine can be determined by changing the geometric shape and size of the gamma-type Stirling engine, which has maximum conversion efficiency in a specific temperature range.

Heat Transfer
We assume that the temperature function of the GSCH only depends on the distance of the z-axis.The heat transfer model can be simplified by considering an ideal cylinder.Consider that the thermal conductivity (ignoring convection and radiation) and temperature have a linear relationship.Therefore, the heat transfer rate is q (Figure 1).The heat transmitted through the hot plate, air room, and cold plate can be written as follows: where  Hot z ,  Air z , and  Hot z are the thicknesses of the hot plate, air room, and cold plate, respectively, and ADc the cross-sectional area of the air room.The thermal conductivity of the hot (kHot) and cold (kCold) plates can be obtained easily.However, the rigorous theory for the determination of the thermal conductivity of the air room does not exist.A forced convection enclosed flow field occurs in the air room of the Stirling engine.In this study, the heat transfer contribution of the displacer was not considered.The empirical relations for free convection are provided in the enclosures [27].The thermal conductivity of the hot (k Hot ) and cold (k Cold ) plates can be obtained easily.However, the rigorous theory for the determination of the thermal conductivity of the air room does not exist.A forced convection enclosed flow field occurs in the air room of the Stirling engine.In this study, the heat transfer contribution of the displacer was not considered.The empirical relations for free convection are provided in the enclosures [27].

Kinematics
Figure 2 displays the major components and strokes of the gamma-type Stirling engine.The displacement functions of the piston and displacer are represented using the angular displacement of the engine shaft as follows: where z p and z are the vertical positions of the piston and displacer, respectively, and θ the rotation angle of the shaft.The displacer is in its highest position (z = L) for at θ = 0.The phase difference between the movements of the displacer and piston is π/2.
where zp and z are the vertical positions of the piston and displacer, respectively, and θ the rotation angle of the shaft.The displacer is in its highest position (z = L) for at θ = 0.The phase difference between the movements of the displacer and piston is π/2.The displacer and piston perform periodical vertical movement when the Stirling engine functions with an angular velocity ω.The relationships between the vertical movement of the piston and displacer and the angular velocity are given as follows: Therefore, the Stirling engine performs three functions: (i) gas compression and expansion, (ii) mechanical motion, and (iii) friction loss when the engine is running.According to the first law of thermodynamics, the work of engine is given as follows: where Wgas and Wmech are the work of the gas and mechanical component, respectively, and Wloss the work of friction loss.

Gas Work
The source of power of the Stirling engine is the pressure gradient in air room, which is formed because of the difference in the temperature of the cold and hot plates.In this section, we refer to the method of Kim et al. to calculate the gas work of the Stirling engine [8].Because of the low speed and small size of the gamma-type Stirling engine, the airflow in the air room can be reasonably assumed to be laminar.Moreover, at 10 RPS, the maximum Reynolds number is 985 as per the mean velocity and hydraulic diameter for the flow through the annulus between the displacer and air room wall.The air room comprises a compression space (VhD + Vhd) and an expansion space (VcD + Vcd), as displayed in Figure 3.The displacer and piston perform periodical vertical movement when the Stirling engine functions with an angular velocity ω.The relationships between the vertical movement of the piston and displacer and the angular velocity are given as follows: Therefore, the Stirling engine performs three functions: (i) gas compression and expansion, (ii) mechanical motion, and (iii) friction loss when the engine is running.According to the first law of thermodynamics, the work of engine is given as follows: where W gas and W mech are the work of the gas and mechanical component, respectively, and W loss the work of friction loss.

Gas Work
The source of power of the Stirling engine is the pressure gradient in air room, which is formed because of the difference in the temperature of the cold and hot plates.In this section, we refer to the method of Kim et al. to calculate the gas work of the Stirling engine [8].Because of the low speed and small size of the gamma-type Stirling engine, the airflow in the air room can be reasonably assumed to be laminar.Moreover, at 10 RPS, the maximum Reynolds number is 985 as per the mean velocity and hydraulic diameter for the flow through the annulus between the displacer and air room wall.The air room comprises a compression space (V hD + V hd ) and an expansion space (V cD + V cd ), as displayed in Figure 3.In Figure 3, VhD and VcD change with the displacer movement and Vp changes with the piston movement.These relationships are expressed as follows: where ADc and Ap are the cross-sectional area of the air room and piston, respectively.The rate of airflow from the compression space to the expansion space depends on the pressure gradient across the displacer and the viscous force on the lateral surface of the displacer.According to approximation, the flow rate of the flow passing the annulus of the displacer is given as follows: where μ is the dynamic viscosity of air.Through Equation ( 6), the pressure gradient in the displacer annulus written as follows: We assume that the pressure gradient in the displacer annulus can be approximated using the following equation: where pc and ph are the pressure of the cold and hot regions in the air room, respectively.The relationship between the cold and hot pressure is given as follows: ( ) The gas in the air room is assumed to be an ideal gas, and no air is assumed to leak.The value of pc depends on the geometry of the air room, temperature, rotation angle of the shaft, angular In Figure 3, V hD and V cD change with the displacer movement and V p changes with the piston movement.These relationships are expressed as follows: where A Dc and A p are the cross-sectional area of the air room and piston, respectively.The rate of airflow from the compression space to the expansion space depends on the pressure gradient across the displacer and the viscous force on the lateral surface of the displacer.According to approximation, the flow rate of the flow passing the annulus of the displacer is given as follows: where µ is the dynamic viscosity of air.Through Equation ( 6), the pressure gradient in the displacer annulus written as follows: We assume that the pressure gradient in the displacer annulus can be approximated using the following equation: where p c and p h are the pressure of the cold and hot regions in the air room, respectively.The relationship between the cold and hot pressure is given as follows: The gas in the air room is assumed to be an ideal gas, and no air is assumed to leak.The value of p c depends on the geometry of the air room, temperature, rotation angle of the shaft, angular velocity, dynamic viscosity of air, and gas constant.
Energies 2019, 12, 1322 6 of 18 where M is the total mass of air in the air room and a 1 -a 9 the coefficients of the Stirling engine (Table A2).Therefore, the gas work of the Stirling engine per cycle is given as follows: where W pp and W pd are the moving-boundary work produced by the piston and displacer per cycle, respectively.

Mechanical Component Work
The engine work comprises three components, namely the work of the (i) flywheel, (ii) shaft assembly, and (iii) piston and displacer systems.The work of the engine is represented as follows: where W wheel is the work of the flywheel, W shaft the work of the shaft, W dis the work of the displacer system, and W pis the work of the piston system.As shown in Figure 4a, there exist two holes in the first quadrant of the flywheel.The two holes cause the center of mass to be located in the third quadrant.Therefore, the work of the flywheel is written as follows: where M wheel is the mass of the flywheel, g the gravity, R new the distance between the center of the flywheel and mass of the center point, and Z wheel the vertical distance (when θ = 1.5π,Z wheel is maximum).Moreover, a bearing links the rotation of the shaft and flywheel (Figure 4a).The work of the shaft is written as follows: where M c and M s are the mass of the bearing and shaft, respectively, and R c and R s the radii of the bearing and shaft, respectively.(10) where M is the total mass of air in the air room and a1-a9 the coefficients of the Stirling engine (Table A2).Therefore, the gas work of the Stirling engine per cycle is given as follows: ( )  (11) where Wpp and Wpd are the moving-boundary work produced by the piston and displacer per cycle, respectively.

Mechanical Component Work
The engine work comprises three components, namely the work of the (i) flywheel, (ii) shaft assembly, and (iii) piston and displacer systems.The work of the engine is represented as follows: where Wwheel is the work of the flywheel, Wshaft the work of the shaft, Wdis the work of the displacer system, and Wpis the work of the piston system.As shown in Figure 4a, there exist two holes in the first quadrant of the flywheel.The two holes cause the center of mass to be located in the third quadrant.Therefore, the work of the flywheel is written as follows: ( ) wheel wheel wheel wheel wheel wheel new where Mwheel is the mass of the flywheel, g the gravity, Rnew the distance between the center of the flywheel and mass of the center point, and Zwheel the vertical distance (when θ = 1.5π,Zwheel is maximum).Moreover, a bearing links the rotation of the shaft and flywheel (Figure 4a).The work of the shaft is written as follows: where  For convenient calculation, all the components (displacer, piston, screws, and links) are set up in two mass systems: displacer and piston systems (Figure 4b).The work of the displacer and piston systems is then given as follows: where M 1 and M 2 are the mass of the displacer and piston systems, respectively.As indicated by Equations ( 12)-( 15), the work of the engine depends on the angular velocity of the engine.

Friction Loss
The work of friction falls into two broad categories: solid and gas friction loss.Among them, solid friction loss occurs at contact surfaces between the bearing and shafts.Furthermore, the volume of solid friction depends on the torque acting on the bearing.The work of friction in one cycle can be represented as follows: where µ k is the coefficient of kinetic friction, τ shaft the load, V tan the tangential velocity of the shaft, T per the period, f the frequency, and α the angular acceleration.The method of [8] is used to analyze the work of gas friction loss.The analytical method is applicable for (i) steady flow, (ii) incompressible laminar flow, and (iii) completely developed flow in the air room.The shear stresses on the lateral surfaces of the displacer, piston, and displacer rod are given as follows:  17) is applied to obtain the work of friction loss from the viscous drags in the displacer, piston, and displacer rod per cycle.
where F s , F sp , and F sDR are the shear forces on the lateral surfaces of the displacer, piston, and displacer rod, respectively.
By integrating Equations (18), the relationship between the temperature difference (∆T = T 1 − T 4 ) and the frequency of the gamma-type Stirling engine movement can be obtained: where c 1 , c 2 , and c 3 are functions such as the mass, geometry, angle, force, kinetic friction coefficient, and dynamic viscosity (Table A1).

Electrostatic Capacitance
The capacitive energy harvester consists of two electrode plates: piston and piston cylinder (Figure 5).According to the capacitive formula, the capacitance value is proportional to the overlap area of the two electrode plates and the frequency of engine rotation.Therefore, changes in the capacitance value depend on the Stirling engine movement.

Electrostatic Capacitance
The capacitive energy harvester consists of two electrode plates: piston and piston cylinder (Figure 5).According to the capacitive formula, the capacitance value is proportional to the overlap area of the two electrode plates and the frequency of engine rotation.Therefore, changes in the capacitance value depend on the Stirling engine movement.Krupenkin and Taylor proposed a method to simulate the energy harvester with periodic movement [28].The equivalent circuit is displayed in Figure 5a.The circuit has a bias voltage of Vbias, resistive load of R, and variable capacitance of C. According to Kirchhoff's circuit laws, the RC (resistance and capacitance) circuit is represented as follows: where Qc is the charge on the variable capacitor, C0 the maximum capacitance of the variable capacitor, Vbias the bias voltage, and t the time.The variation in the capacitance with time is expressed as follows: Krupenkin and Taylor proposed a method to simulate the energy harvester with periodic movement [28].The equivalent circuit is displayed in Figure 5.The circuit has a bias voltage of V bias , resistive load of R, and variable capacitance of C. According to Kirchhoff's circuit laws, the RC (resistance and capacitance) circuit is represented as follows: where Q c is the charge on the variable capacitor, C 0 the maximum capacitance of the variable capacitor, V bias the bias voltage, and t the time.The variation in the capacitance with time is expressed as follows: where A over is the maximum overlap area between the outer part of the piston surface and the piston cylinder during the Stirling engine cycle, h the dielectric thickness of the piston cylinder, k die the dielectric constant of the thin film, and ε 0 the vacuum permittivity.Equation ( 20) is transformed into a dimensional equation as follows: where Qch = Q h V −1 bias C −1 0 , a 0 = ωRC 0 , and t = ωt.The solution of Equation ( 22) is given as follows: and Equation ( 24) can be written with the original variables as follows: Therefore, the average generated power is given as follows:

Result and Discussion
The temperature difference is used to change the volume of the air room, which in turn induces a change in the piston position.As indicated by Equation ( 5), the sizes of the compression space (V hD + V hd ) and expansion space (V cD + V cd ) depend on the phase angle of the wheel.In a gamma-type Stirling cycle, the isothermal compression is θ = 13~3(π/8), constant volume heating is θ = 4 (π/8), isothermal expansion is θ = 5-11 (π/8), and constant volume cooling is θ = 12 (π/8).When θ = 4 (π/8), the air room has maximum volume.The pressure difference between the hot and cold regions varies with the sinusoid motion of the displayer at fixed temperature differences (Figure 6).The variations in p h and p c for θ = 1-4 (π/8), θ = 9-12 (π/8), and ∆T = 0-30 K are calculated using Equations ( 9) and (10), where ω = 2πf for each ∆T is the result obtained from Equation (19)  The pressure difference in the air room of the gamma-type Stirling engine is very small because the thickness of the air room is only 10 mm.Consequently, hot and cold air mix quickly in the air room.Thus, the Stirling engine works in a steady state, which only requires a pressure difference of less than 0.1 Pa in the air room.
The relationship between the temperature differences of the Stirling engine and the frequency of the wheel can be obtained using Equation (19).Assume that the cold and hot plates have an initial temperature of 293 K (ΔT = 0 K) and the hot plate is then heated to 323 K (ΔT = 30 K).In this scenario, Equation ( 19) has one real root and two nonreal complex conjugate roots for each temperature difference.The solid line (real root) in Figure 7 represents the simulated results obtained using Equations( 1)-( 19), whereas the dotted line represents the fitting results obtained from an experiment.For the experimental results, the wheel rotation is determined in revolutions per minute, whereas for the theoretical results, the wheel rotation is determined in revolutions per second (Hz).There exist 63 experiment points, and the fitting line is a cubic equation, which is  The pressure difference in the air room of the gamma-type Stirling engine is very small because the thickness of the air room is only 10 mm.Consequently, hot and cold air mix quickly in the air room.Thus, the Stirling engine works in a steady state, which only requires a pressure difference of less than 0.1 Pa in the air room.
The relationship between the temperature differences of the Stirling engine and the frequency of the wheel can be obtained using Equation (19).Assume that the cold and hot plates have an initial temperature of 293 K (∆T = 0 K) and the hot plate is then heated to 323 K (∆T = 30 K).In this scenario, Equation ( 19) has one real root and two nonreal complex conjugate roots for each temperature difference.The solid line (real root) in Figure 7 represents the simulated results obtained using Equations ( 1)-( 19), whereas the dotted line represents the fitting results obtained from an experiment.For the experimental results, the wheel rotation is determined in revolutions per minute, whereas for the theoretical results, the wheel rotation is determined in revolutions per second (Hz).There exist 63 experiment points, and the fitting line is a cubic equation, which is represented as follows: ∆T = −0.0003f 3 + 0.0105 f 2 + 0.0975 f .
The relationship between the temperature differences of the Stirling engine and the frequency of the wheel can be obtained using Equation (19).Assume that the cold and hot plates have an initial temperature of 293 K (ΔT = 0 K) and the hot plate is then heated to 323 K (ΔT = 30 K).In this scenario, Equation ( 19) has one real root and two nonreal complex conjugate roots for each temperature difference.The solid line (real root) in Figure 7 represents the simulated results obtained using Equations( 1)-( 19), whereas the dotted line represents the fitting results obtained from an experiment.For the experimental results, the wheel rotation is determined in revolutions per minute, whereas for the theoretical results, the wheel rotation is determined in revolutions per second (Hz).There exist 63 experiment points, and the fitting line is a cubic equation, which is represented as follows:  = − + + 32 0.0003 0.0105 0.0975 T f f f .As displayed in Figure 7, the triggered temperature difference of the Stirling engine is at least 12 K because the engine must overcome the friction and mass of the mechanical system (M1, M2, Ms, Mc, and Mwheel).The activated resistance force is [M1+μk(M2+Ms+Mc+Mwheel)]g = 0.1731 N. Furthermore, through Equation ( 10), the difference in the pressures at 293 (P0) and 305 (P12) K is calculated to be 27.4Pa for θ = 0.The pressure on the displacer is equal to 0.1743 N. As displayed in Figure 7, the triggered temperature difference of the Stirling engine is at least 12 K because the engine must overcome the friction and mass of the mechanical system (M 1 , M 2 , M s , M c , and M wheel ).The activated resistance force is [M 1 +µ k (M 2 +M s +M c +M wheel )]g = 0.1731 N. Furthermore, through Equation ( 10), the difference in the pressures at 293 (P 0 ) and 305 (P 12 ) K is calculated to be 27.4Pa for θ = 0.The pressure on the displacer is equal to 0.1743 N.
The work system of the Stirling engine is divided into the mechanical, gas, and friction parts (Figure 8a).The work of the mechanical part comprises W dis , W pis , W shaft , and W wheel .The motion of the displacer and piston is vertical simple harmonic motion through wheel rotation.Therefore, the lines of W dis and W pis depend on the position functions of the displacer and piston.Because the Stirling engine works in a steady state, the angular velocity is independent of time.Therefore, W shaft is a constant volume.Moreover, there exist two holes in the first quadrant of the flywheel.Consequently, the moment of inertia depends on the phase angle (Figure 8b).The work of the gas contains W pp and W pd .Both the lines displayed in Figure 8c are sine functions.When W pp is a positive volume, the fluid in the engine outputs work to the piston; otherwise, the piston compresses the fluid in the engine.The displacer only moves the fluid in the air room and cannot compress or expand the fluid.Therefore, the volume of W pd is positive.Two peak volumes exist for W pd .These peak volumes represent the displacer touching the clearance in the hot and cold areas (Figure 8c).At the friction work part (Figure 8d), when F displacer and F piston drive the shaft rotation, the work of friction (Ω shaft ) reduces the mechanical work of gas expansion.Ω sp and Ω sDR are produced when the piston and displacer rod friction fluid (negative value).Ω s is produced when the fluid is reflected from the cold plate inside friction the lateral surface of displacer (positive value, reverse direction).
expand the fluid.Therefore, the volume of Wpd is positive.Two peak volumes exist for Wpd.These peak volumes represent the displacer touching the clearance in the hot and cold areas (Figure 8c).At the friction work part (Figure 8d), when Fdisplacer and Fpiston drive the shaft rotation, the work of friction (Ωshaft) reduces the mechanical work of gas expansion.Ωsp and ΩsDR are produced when the piston and displacer rod friction fluid (negative value).Ωs is produced when the fluid is reflected from the cold plate inside friction the lateral surface of displacer (positive value, reverse direction).To compare the effect of the area of heated plate with the Stirling engine, a dimensionless length ratio is defined.The length ratio is equal to the diameter of the displacer divided by the diameter of the wheel (D/2Rw).The length ratio is proportional to the heat flux of the air room if the gap between the displacer and displacer cylinder is fixed.Figure 9a displays the relationship between the length ratio, temperature difference, and frequency of the wheel.When the temperature of the hot plate is 30 K higher than that of the cold plate, the wheel rotation frequency equals 3.35 Hz.However, when the length ratio is increased by six times, the rotation speed increases by 2.69 times.Figure 9b illustrates power required (Wwheel + Wdis + Wshaft + Wpis + Wpd + Wpp + Ωs + Ωsp + ΩsDR + Ωshaft) by the Stirling engine in the steady state.When the length ratio increases, the mass of the wheel increases.Therefore, the engine requires additional energy to push the wheel.To compare the effect of the area of heated plate with the Stirling engine, a dimensionless length ratio is defined.The length ratio is equal to the diameter of the displacer divided by the diameter of the wheel (D/2R w ).The length ratio is proportional to the heat flux of the air room if the gap between the displacer and displacer cylinder is fixed.Figure 9a displays the relationship between the length ratio, temperature difference, and frequency of the wheel.When the temperature of the hot plate is 30 K higher than that of the cold plate, the wheel rotation frequency equals 3.35 Hz.However, when the length ratio is increased by six times, the rotation speed increases by 2.69 times.Figure 9b illustrates power required (W wheel + W dis + W shaft + W pis + W pd + W pp + Ω s + Ω sp + Ω sDR + Ω shaft ) by the Stirling engine in the steady state.When the length ratio increases, the mass of the wheel increases.Therefore, the engine requires additional energy to push the wheel.
The aforementioned frequency (Figures 7-9) of the Stirling engine is simulated through Equations ( 25) and (26).The dielectric coefficient of acrylic (ε a ) is 3.7, and the dielectric thickness (h) is 1 µm.The energy and power produced can be calculated.As displayed in Figure 10a,b, different bias voltages (26, 44, 53, 61, and 70 V) are produced between the two plates of the capacitor.The variable capacitance produces power along the piston of the Stirling engine.As per the computational results, the energy harvester produces 6.2 mJ of energy per cycle per unit area at a bias voltage of 70 V and temperature difference of 30 • C. Thus, the engine produced 21 mW of energy per unit area.The aforementioned frequency (Figures 7-9) of the Stirling engine is simulated through Equations ( 25) and (26).The dielectric coefficient of acrylic (εa) is 3.7, and the dielectric thickness (h) is 1 μm.The energy and power produced can be calculated.As displayed in Figures 10a and b, different bias voltages (26,44,53, 61, and 70 V) are produced between the two plates of the capacitor.The variable capacitance produces power along the piston of the Stirling engine.As per the computational results, the energy harvester produces 6.2 mJ of energy per cycle per unit area at a bias voltage of 70 V and temperature difference of 30°C .Thus, the engine produced 21 mW of energy per unit area.The aforementioned frequency (Figures 7-9) of the Stirling engine is simulated through Equations ( 25) and (26).The dielectric coefficient of acrylic (εa) is 3.7, and the dielectric thickness (h) is 1 μm.The energy and power produced can be calculated.As displayed in Figures 10a and b, different bias voltages (26, 44, 53, 61, and 70 V) are produced between the two plates of the capacitor.The variable capacitance produces power along the piston of the Stirling engine.As per the computational results, the energy harvester produces 6.2 mJ of energy per cycle per unit area at a bias voltage of 70 V and temperature difference of 30°C .Thus, the engine produced 21 mW of energy per unit area.The conversion efficiency of heat energy into electric power is calculated using the formula ηhe=100(P/q).The conversion efficiency of engine power into electric power is calculated using the formula ηee=100(P/η).In this study, η=1-(T4/T1), Aover = ADc, h = 1 μm, and Vbias = 70 V.A dimensional number (ss) is delineated to adjust the size of the Stirling engine.For instance, an ss value of 2 indicates that the length, mass, and gap between the displacer and displacer cylinder of the Stirling The conversion efficiency of heat energy into electric power is calculated using the formula η he = 100(P/q).The conversion efficiency of engine power into electric power is calculated using the formula η ee = 100(P/η).In this study, η = 1 − (T 4 /T 1 ), A over = A Dc , h = 1 µm, and V bias = 70 V.A dimensional number (ss) is delineated to adjust the size of the Stirling engine.For instance, an ss value of 2 indicates that the length, mass, and gap between the displacer and displacer cylinder of the Stirling engine become two times the dielectric thickness.The efficiency of the heat engine power (η) is proportional to the temperature difference.However, the relationship between the frequency of the wheel and temperature is a concave curve.When ∆T > 10 K, the rate of increase decreases.Therefore, the conversion efficiency of engine power into electric power follows the trend of a convex curve (Figure 11a).Considering the dimensional efficiency, a large energy harvester is suitable for high temperature differences.For instance, η he (ss = 150) is twice η he (ss = 30) at ∆T = 60-90 K.In Figure 10, each curve of η ee has a maximum value [η ee,max (ss = 30) at ∆T = 7 K, η ee,max (ss = 60) at ∆T = 18 K, and η ee,max (ss = 90) at ∆T = 30 K], the corresponding temperature difference is the best temperature difference (T bt ) for the conversion efficiency of heat energy into electric power.Thus, η ee exhibits positive growth for temperature differences below T bt .
(b) The conversion efficiency of heat energy into electric power is calculated using the formula ηhe=100(P/q).The conversion efficiency of engine power into electric power is calculated using the formula ηee=100(P/η).In this study, η=1-(T4/T1), Aover = ADc, h = 1 μm, and Vbias = 70 V.A dimensional number (ss) is delineated to adjust the size of the Stirling engine.For instance, an ss value of 2 indicates that the length, mass, and gap between the displacer and displacer cylinder of the Stirling engine become two times the dielectric thickness.The efficiency of the heat engine power (η) is proportional to the temperature difference.However, the relationship between the frequency of the wheel and temperature is a concave curve.When ΔT > 10 K, the rate of increase decreases.Therefore, the conversion efficiency of engine power into electric power follows the trend of a convex curve (Figure 11a).Considering the dimensional efficiency, a large energy harvester is suitable for high temperature differences.For instance, ηhe (ss = 150) is twice ηhe (ss = 30) at ΔT = 60-90 K.In Figure 10d, each curve of ηee has a maximum value [ηee,max(ss = 30) at ΔT = 7 K, ηee,max(ss = 60) at ΔT =18 K, and ηee,max(ss = 90) at ΔT = 30 K], the corresponding temperature difference is the best temperature difference (Tbt) for the conversion efficiency of heat energy into electric power.Thus, ηee exhibits positive growth for temperature differences below Tbt.

Conclusion
In this study, a mathematical model that considers the conversion of heat energy to mechanical energy and the mechanical energy to electric energy is proposed and the relationship between the frequency and temperature difference is verified.The simulated results indicate that the frequency of the wheel is 3.35 Hz at ΔT = 30 K (293-323 K) and the triggered temperature difference of the Stirling engine is 12 K (293-305 K).A small pressure difference exists between the hot and cold regions in the air room of the gamma Stirling engine ( 0.1 Pa at ΔT = 30 K).When heat enters the work system of the Stirling engine, the energy can be divided into three parts: mechanical, air, and

Conclusions
In this study, a mathematical model that considers the conversion of heat energy to mechanical energy and the mechanical energy to electric energy is proposed and the relationship between the +2π 2 0.5M wheel R Friction loss from the viscous drags in the displacer.
Friction loss from the viscous drags in the piston.
Friction loss from the viscous drags in the rod. 6µL Potential energy of the M 1 , M 2 systems, and flywheel component.

Figure 1 .
Figure 1.One-dimensional heat transfer (conduction only) through the GSCH and electrical analog.

Figure 2
Figure 2 displays the major components and strokes of the gamma-type Stirling engine.The displacement functions of the piston and displacer are represented using the angular displacement of the engine shaft as follows:

Figure 1 .
Figure 1.One-dimensional heat transfer (conduction only) through the GSCH and electrical analog.

Figure 2 .
Figure 2. Nomenclature and major components of the gamma-type Stirling engine.

Figure 2 .
Figure 2. Nomenclature and major components of the gamma-type Stirling engine.

Figure 3 .
Figure 3. Volumes in thermal flow of the gamma-type Stirling engine.

Figure 3 .
Figure 3. Volumes in thermal flow of the gamma-type Stirling engine.
Mc and Ms are the mass of the bearing and shaft, respectively, and Rc and Rs the radii of the bearing and shaft, respectively.

Figure 4 .
Figure 4. Nomenclature and major components of flywheel power systems.(a) The shafts, bearing, and flywheel; (b) The M 1 and M 2 systems.

Figure 5 .
Figure 5. Schematic and equivalent circuit diagram of energy harvesters.

Figure 5 .
Figure 5. Schematic and equivalent circuit diagram of energy harvesters.

Figure 6 .
Figure 6.Pressure difference lines (∆p(θ) = p h (θ) − p c (θ)) in the air room at different θ values for different temperature differences between the cold and hot regions of the Stirling engine (T 4 = 293 K, T 1 = 293-323 K).

Figure 7 .
Figure 7. Relationship between the temperature difference of the Stirling engine and the frequency of the wheel.Small circles represent the experimental data, solid line the simulation data, and dotted line the fitting line with the cubic equation.

Figure 7 .
Figure 7. Relationship between the temperature difference of the Stirling engine and the frequency of the wheel.Small circles represent the experimental data, solid line the simulation data, and dotted line the fitting line with the cubic equation.

Figure 8 .
Figure 8. Distribution of work from the heat source in one cycle of the gamma-type Stirling engine.(a) Major components in the engine, (b) work values of the gas part, (c) work values of the mechanical part, and (d) work values of the friction part.

Figure 8 .
Figure 8. Distribution of work from the heat source in one cycle of the gamma-type Stirling engine.(a) Major components in the engine, (b) work values of the gas part, (c) work values of the mechanical part, and (d) work values of the friction part.

Figure 9 .
Figure 9. Relationship between the length ratio (D/2Rw) and frequency of the wheel at various temperature differences.(a) Frequency of the wheel with different length ratios (0-6) at various temperature differences (ΔT = 0-30 K).(b) Required power of the engine with different length ratios (0-6) at different temperature differences (ΔT = 0-30 K).

Figure 9 .
Figure 9. Relationship between the length ratio (D/2R w ) and frequency of the wheel at various temperature differences.(a) Frequency of the wheel with different length ratios (0-6) at various temperature differences (∆T = 0-30 K).(b) Required power of the engine with different length ratios (0-6) at different temperature differences (∆T = 0-30 K).

Figure 9 .
Figure 9. Relationship between the length ratio (D/2Rw) and frequency of the wheel at various temperature differences.(a) Frequency of the wheel with different length ratios (0-6) at various temperature differences (ΔT = 0-30 K).(b) Required power of the engine with different length ratios (0-6) at different temperature differences (ΔT = 0-30 K).

Figure 10 .
Figure 10.Results of the energy generated per unit area per unit rotating cycle at the wheel for various temperature differences.(a) The energy produced with the Stirling engine energy harvester at different bias voltages.(b) The power produced with the Stirling engine energy harvester at different bias voltages.

Figure 10 .
Figure 10.Results of the energy generated per unit area per unit rotating cycle at the wheel for various temperature differences.(a) The energy produced with the Stirling engine energy harvester at different bias voltages.(b) The power produced with the Stirling engine energy harvester at different bias voltages.

Figure 10 .
Figure 10.Results of the energy generated per unit area per unit rotating cycle at the wheel for various temperature differences.(a) The energy produced with the Stirling engine energy harvester at different bias voltages.(b) The power produced with the Stirling engine energy harvester at different bias voltages.

Figure 11 .
Figure 11.Results for the dimensional efficiency (dimensional number = 30, 60, 90, 120, and 150) of the Stirling engine energy harvester at various temperature differences.Conversion efficiency of (a) engine power into electric power for the Stirling engine energy harvester and (b) heat energy into electric power for the Stirling engine energy harvester.

Figure 11 .
Figure 11.Results for the dimensional efficiency (dimensional number = 30, 60, 90, 120, and 150) of the Stirling engine energy harvester at various temperature differences.Conversion efficiency of (a) engine power into electric power for the Stirling engine energy harvester and (b) heat energy into electric power for the Stirling engine energy harvester.
Thermal resistance in series of the Stirling engine components.

Ef
Energy of Stirling engine harvester output per cycle J F s Shear force on displacer lateral surface N F sp Shear force on displacer rod lateral surface N F sDR Shear force on piston lateral surface N Frequency Length of clearance in the cold end of displacer cylinder 2.5 × 10 −3 M l h Length of clearance in the hot end of displacer cylinder 2.5 × 10 −3 M l D Displacer length 8 × 10 −3 M l DR Length of the bushing of displacer rod 3 × 10 −2 M M Total mass or air room kg 2 wheel + M wheel + R 2 new (R con )

Table A3 .
Nomenclature and dimensions of Gamma-type Stirling engine components.