Arxiv:1812.04415V1 [Physics.Ed-Ph] 4 Nov 2018 Stration Stirling Engine

Home , Stirling cycle

Quantitative Stirling Cycle Measurements: P-V Diagram and Refrigeration Y. J. Lu,1 Hiroko Nakahara,1 and J. S. Bobowski1, a) Department of Physics, University of British Columbia, Kelowna, British Columbia, Canada V1V 1V7 (Dated: December 12, 2018)

I. INTRODUCTION

This paper describes simple modifications to a demonstration Stirling engine that allows one to make quantitative measurements of the Stirling cycle. First, we de- scribe measurements of the P -V diagram using an inexpensive pressure sensor and a common photogate. Then, as a supplement, the engine was run as a refrigerator by using a rotary tool to turn the flywheel of the engine. A small, but measurable, temperature difference developed across the body of the engine. Stirling engines have long been used as in-class demonstrations and,1,2 in some cases, P -V diagrams have be constructed using a commercial apparatus.3 It is also worth noting that, despite the fact that Stirling’s original patent was filed over 200 year ago, the Stirling cycle is still being actively studied.4,5 The main objective of our work is to make an accurate measurement of the P -V diagram of a high-quality Stirling engine using inexpensive materials and a data acquisition system that Figure 1. Schematic diagram of the displacer-style Stirling engine. The power piston and displacer move 90◦ out of phase is commonly found in undergraduate laboratories. with one another. A temperature difference with T2 > T1 will This work was inspired by an earlier measurement of result in a counterclockwise rotation of the flywheel. a Stirling engine P -V diagram by Nakahara.6 In that project, the P -V diagram was measured using Vernier’s Logger Pro data acquisition system. Pressure data When covering the topic of heat engines in an introduc- were acquired at a rate of 10 samples per second us- tory thermodynamics course, one typically first describes ing Vernier’s pressure sensor. The gas volume was mea- the Carnot cycle.7 In the Carnot cycle, the isothermal sured by mounting a reflector onto the piston of the en- compression of an ideal gas results in heat Q1 being gine and using an infrared (IR) transmitter/receiver to released to a cold thermal reservoir at temperature T1 track the piston’s position. During operation, the cal- whereas an isothermal expansion causes the gas to ab- ibrated output voltage of the IR sensor was monitored sorb heat Q2 from a hot thermal reservoir at temperature using Vernier’s differential voltage probe. T2 > T1. The efficiency ε of the cycle is determined by The advantage of Nakahara’s scheme is that it pro- the fraction of the heat removed from the hot reservoir duced a real-time measurement of the P -V diagram. that is converted into useful work W : There are, however, a number of limitations and disad- W T vantageous to this scheme. One is that the pressure sen- ε = = 1 1 (1) sor does not have the sensitivity to to precisely measure Q2 − T2 the small pressure changes associated with the demon- In addition to the isothermal processes, the Carnot cycle arXiv:1812.04415v1 [physics.ed-ph] 4 Nov 2018 stration Stirling engine. Second, the time between pressure measurements is long compared to the period of the requires an adiabatic compression and expansion to take Stirling cycle. Third, the calibration of the IR sensor is the gas from T1 to T2 and T2 to T1, respectively. The nonlinear such that the measurement sensitivity of the Carnot cycle is conceptually simple because there is no piston position is low when the reflector is furthest from heat exchanged between the ideal gas and the thermal the IR transmitter/detector. Finally, to mount the IR reservoirs during the adiabatic processes. The cycle is, sensor, a new top plate for the Stirling engine had to be however, difficult to implement in practice precisely be- machined. In this project, we developed a measurement cause one cannot easily isolate the gas from the thermal scheme to overcome these limitations which allowed us to reservoirs. produce high-quality pressure and volume measurements The Stirling cycle, on the other hand, is easier to imple- when using the same Striling engine as Nakahara. ment but more challenging conceptually since it requires the idea of a regenerator. In the Stirling cycle, the adiabatic processes of the Carnot cycle are replaced with isochoric (constant volume) processes. The isochoric parts of the cycle do no work, however, the gas must absorb a) Electronic mail: [email protected] (release) heat as its pressure and temperature increases 2

Pressure Sensor 2.0 data ﬁt Chopper Photogate 1.5

(kPa) 1.0 P ∆

Voltage probe 0.5

Figure 2. Photograph of the setup to measure the Striling cycle P -V diagram. A photogate and chopper are used to 0.0 sense the position of the power piston while a pressure sensor 3 4 4 5 and voltage probe monitor the gas pressure. Pressure Sensor Voltage (V)

(decreases) at constant volume. This heat exchange oc- Figure 3. Pressure sensor calibration data for the MPXV7002 pressure sensor. The measured temperature change of a fixed curs via a regenerator through which the gas is forced to volume of air was converted to a pressure difference while pass. The regerenator is typically a porous mesh with simultaneously recording the voltage output of the pressure a large heat capacity. It can absorb or release a small sensor. amount of heat without a substantial change to its temperature. An analysis of the ideal Stirling cycle shows that it has the same efficiency, given by Eq. 1, as the displacer-type Stirling engine. A list of the all of the Carnot cycle. equipment and materials used in our experiments is given The Stirling cycle is more practical because the con- in Appendix A. Where appropriate, possible vendors are stant volume processes can be implemented in a number suggested and estimates of costs are given. of different ways. The Stirling engine used in our project is made by Kontax Stirling Engines.8 The main elements of the engine, shown schematically in Fig.1, are a power A. Pressure measurement piston that sets the volume of the working gas, a displacer that is used to move the gas between the hot and cold re- To monitor the pressure change of the gas in the Stir- gions of the engine, and a flywheel. The positions of the ling engine, an inexpensive differential pressure sensor, piston and displacer are set by connecting rods that are having the optimal sensitivity for our application, was attached to a disk centered on the flywheel. As shown purchased. From Nakahara’s previous work, the pressure in the figure, the connecting rods are attached such that change was expected to be just over 1 kPa. Therefore, we the piston and displacer positions are 90◦ out of phase purchase the MPXV7002 pressure sensor manufactured with one another. In Fig. 1, a small counterclockwise ro- by NXP Semiconductors which has a voltage output that tation of the flywheel would result in an (approximately) varies from 0.5 to 4.5 V for pressure changes from 2 to isothermal expansion of the gas while it is in contact with +2 kPa. The datasheet supplied by NXP Semiconduc-− the hot reservoir. This would be followed by an (approx- tors states that the typical calibration of the output is imately) isochoric process that would decrease the gas given by: temperature from T2 to T1.

This displacer-style Stirling engine does not have a ∆P = (Vout 2.5 V) 1 kPa/V (2) porous regenerator through which the gas passes dur- − × ing the constant volume processes. However, the surface where ∆P is the pressure difference across the sensor of the displacer and the cylindrical wall separating the in kilopascals and Vout is the sensor output voltage in hot and cold thermal reservoirs together act as a crude volts. (Equation 2 assumes that the sensor is powered by regenerator. With such a poor regenerator, the work +5 V.)9 output per cycle will be low and the engine’s efficiency is To precisely calibrate our specific pressure sensor, we expected to be much less than that of the ideal Stirling developed a simple calibration apparatus. Rubber tub- cycle. ing was used to connect the pressure sensor to one end of a length of standard 1/4-inch copper tubing. (The copper tubing was approximately 1 m long, but the precise II. P-V DIAGRAM length is not critical and does not need to be known.) A type-T thermocouple was then inserted into the op- Figure 2 shows the experimental setup used to si- posite end until its tip was approximately in the center multaneously measure the pressure and volume of the of the tube. Finally, the open end of the copper tube 3

(through which the thermocouple was fed) was sealed using J-B Weld epoxy. The temperature of the copper tube was then very slowly raised above room temperature by about 6 ◦C using a T-shirt heat press. The heat press was plugged into a variac which was used to drop the supply voltage from 110 to 15 Vrms so as to limit the temperature increase. The change in pressure of the air trapped in the copper tube is related to the temperature change via: ∆T ∆P = P (3) T where P and T are the ambient pressure and temperature of the room, respectively. The calibration data obtained are shown in Fig. 3. A linear ﬁt to the data yielded the following relationship between ∆P of the air and the pressure sensor output voltage:

∆P = (1.0485V 2.833 V) 1 kPa/V (4) out − × which is in reasonable agreement with manufacturer Figure 4. The geometry used to convert flywheel position in specifications. radians to piston height ∆y. To connect the pressure sensor to the Stirling engine, a hole was drilled through the engine’s aluminum top plate to accommodate the tip of a Luer lock syringe. The where the approximate expression is valid when r `. syringe was epoxied in place using J-B Weld. As shown The change in volume is then given by ∆y times the in Fig. 2, the pressure sensor was then connected to the cross-sectional area of the piston. syringe tip using tight-fitting rubber tubing. The figure The angular position of the flywheel was monitored us- shows a modified top plate, however, all that is needed ing a Vernier photogate and Logger Pro. The software is an appropriately-sized hole in the original plate. We was set to record the times at which the photogate went used a modified top plate only because it was left over from a blocked to an unblocked state. Under these con- from Nakahara’s earlier project. During operation, the ditions, after passing through the θ = 0 point, the angle pressure sensor’s output voltage was recorded at a rate at the nth recorded time is given by: of 200 samples/s using Vernier’s Logger Pro software and differential voltage probe. 2π 1 θn = n (7) N − 4 0

B. Volume measurement where N0 = 100 is the number of chopper intervals (before blocking one of the transparent sections), and n is an integer that runs from 1 to N 1. To calculate the change in gas volume, the height of 0 − the Stirling engine’s piston was measured. A circular chopper, shown in Fig. 9 of Appendix B, and a photo- C. Results gate were used to track the rotation of the engine’s flywheel. The chopper was designed to have 100 uniform dark and transparent intervals. One of the transparent Figure 5(a) shows the measured pressure and volume intervals was intentionally blacked out. When the chop- change of the Kontax Stirling engine as a function of per was attached to the flywheel, as shown in Fig. 2, the time. For these data, the Stirling engine was placed on blackened segment was positioned such that it was at the top of a cup of hot water. The water was boiled in a top of the flywheel when the piston was at its maximum kettle before it was poured into the cup. The data ac- height. This arrangement allows one to use the chopper quisition was started after the Stirling engine reached data to set the angular position of the flywheel to θ = 0 a speed that was approximately constant (the speed of when the gas volume is a maximum. the engine slowly decreases as the hot water cools). The Figure 4 shows the geometry that relates the flywheel data are high resolution and the acquisition rate is high compared to speed of the engine ( 5.3 cycles/s). The rotation to the piston height. A careful analysis of this ≈ geometry reveals that the change in piston height is given data also convincingly show that the pressure and vol- by: ume changes are 90◦ out of phase and that the pressure leads the volume. ` r 2 In Fig. 6, we show a plot of Stirling engine’s P -V dia- ∆y = r cos θ 1 sin θ 1 (5) gram. The P -V diagram that corresponds to the data in " − r − ` − !# r Fig. 5(a) is shown using orange circles. Because Logger r r cos θ + sin2 θ (6) Pro records the pressure and photogate readings at dif- ≈ 2` ferent times, the mathematics software Maple was used h i 4

0.6 1 0.2 0.2 0.4

0.5 0.1 0.1 0.2 ∆ ∆ V V

cm cm (kPa) (kPa) 0 0 0 0 P P 3 3 ∆ ∆

0.2 − 0.1 0.1 − 0.5 − − 0.4 − 0.2 0.2 − − 1 − 0.6 − 0 0.2 0.4 0.6 0.8 1 0 0.4 0.8 1.2 1.6 2 (a) time (s) (b) time (s)

Figure 5. (a) Measured change in pressure (red circles) and volume (blue squares) versus time with a free-running ﬂywheel. (b) Measured change in pressure (red circles) and volume (blue squares) versus time with light friction applied to the ﬂywheel.

1.0 surement i to measurement i + 1 is given by:

0.6 ∆Pi+1 + ∆Pi Wi = (∆Vi ∆Vi) (8) 2 +1 −

0.2 Summing all of the Wi contributions from one complete (kPa) cycle gives the net work. The result of this calculation P

∆ 0.2 was a net work of 0.43 mJ. Multiplying the work by the − previously-determined cycle rate of 5.3 Hz gives an output power of 2.3 mW, or 3.0 micro-horsepower. A typical 0.6 car might have a 120 horsepower engine – a factor of 40 − million greater than our demonstration Stirling engine! unloaded Finally, the measurements of ∆P and ∆V were re- 1.0 loaded − peated after light friction was applied to the flywheel 0.2 0.1 0.0 0.1 0.2 of the Stirling engine. A small piece of sponge was at- − − ∆V (cm 3) tached to a retort stand and then brought into light contact with the flywheel while the engine was running. The Figure 6. Measured P -V diagram of the Kontax Stirling en- sponge position was adjusted until the equilibrium period gine with (green squares) and without (orange circles) friction of the engine’s cycle noticeably increased. The pressure applied to the flywheel. The added friction causes the engine and volume changes are plotted as a function of time in to do more work per cycle while simultaneously reducing the Fig. 5(b). These data show the same 90◦ phase shift be- power output. tween the pressure and volume. Notice, however, that the magnitude of the pressure change has approximately doubled while the period of the cycle has increased from to interpolate the volume changes at the times of the 0.19 s for the free-running engine to 0.43 s after applying pressure readings. The plot in Fig. 6 shows a set of 25 friction. consecutive cycles that are indistinguishable from one an- The P -V diagram that corresponds to the data in other. Fig. 5(b) is shown using green squares in Fig. 6. Again, 25 Clearly, the measured P -V diagram differs from that consecutive and nearly-identical cycles are shown. The of the ideal Stirling cycle. For truly isothermal processes, increase in enclosed area is an indication that the engine 1 one would expect ∆P ∆V − . Also, as expected from is doing more work per cycle when friction is applied to the design of the connecting∝ rods, there are no prefectly the flywheel. The calculated area is 0.75 mJ. Dividing isochoric processes in the measured cycle. Nevertheless, the work per cycle by the period gives a power output of the high-resolution data do allow for an accurate mea- 1.7 mW or 2.3 micro-horsepower. Although loading the surement of the area enclosed by the loop which gives flywheel causes the engine to do more work per cycle, the the work output by the engine per cycle. Maple was also longer period results in power output that, in this case, used perform this calculation. The work done from mea- was reduced by about 20 %. 5

Differential thermocouple 0.3

0.2

0.1 C) ◦ ( T Dremel ∆ 0.0

Variac 0.1 −

0.2 − Figure 7. Photograph of the setup used to operate the Kontax Stirling engine as a refrigerator. A Dremel tool was used to 0 40 80 120 160 200 rotate the flywheel and a differential thermocouple measured time (minutes) the temperature difference between the top and bottom plates of the engine. Figure 8. Measured temperature difference across the top and bottom plates of the Stirling engine versus time. The raw ∆T data are shown by grey data points. The black data III. REFRIGERATION points show ∆T after removing the background drift of the room temperature, shown by the dotted line. A positive ∆T indicates that the top plate is at a higher temperature than Today, the most common commercial application of the bottom plate. the Stirling cycle is in closed-cylce cryocoolers.10 In a cryocooler, mechanical work is done on the working gas to generate a temperature difference. Several compa- plugged into a variac which was set to approximately nies now make single-stage Striling cryocoolers capable of 25 Vrms. reaching 40 K and two-stage systems capable of reaching Figure 8 shows the difference in temperature between lower temperatures. For a review of modern cryocooler the top and bottom plates of the Stirling engine as a technology, see Ref. 11. function to time. The raw data are shown by the grey To run the Kontax Stirling engine as a cryocooler, it data points. Initially, the Dremel was off and the fly- was first necessary to thermally isolate the top and bot- wheel was stationary. The Dremel was started at ap- tom plates of the engine from the surrounding environ- proximately 50 min and allowed to run continuously for ment. A 2-inch sheet of high-density polystyrene foam 30 min. During that time, the temperature difference in- was cut using a hot 32-AWG nichrome wire. The wire creased and was beginning to plateau by the end of the was under tension and we found that a current of 1 A was 30 min run time. An increasing ∆T corresponds to the sufficient to make make clean cuts. As shown in Fig. 7, top plate warming relative to the bottom plate. After the polystyrene blocks were cut such that the Stirling engine Dremel was stopped, the temperature difference steadily was completely enclosed except for the upper half of the decreased. The Dremel was run for another 30 min inter- flywheel. val starting at approximately 110 min. However, for this The lower-right corner of Fig. 7 shows the differential run the rotation direction of the flywheel was reversed. thermocouple that was used to monitor the temperature As expected, this caused the top plate to cool relative to difference between the top and bottom plates of the Stir- the bottom plate. The Dremel was stopped at 140 min. ling engine. A single #4-40 tapped hole was drilled into Figure 8 clearly shows a nonzero ∆T that evolves lin- each of the top and bottom plates and screws were used early with time even when the Dremel has been stopped to clamp down on the junctions of a type-E thermocou- for an extended period of time. This is due to a slow ple. Fine thermocouple wires were used so as to limit heat drift of the room temperature. The polystyrene used to conduction along the wires as much as possible. The ther- isolate the Stirling engine from its surroundings has a mocouple voltage was measured using a Keysight 34401A hole in the top to accommodate the flywheel. Therefore, multimeter. The data were logged as a function of time the top plate of the engine is more sensitive to changes in using a simple LabVIEW program. room temperature than the bottom plate which results in The flywheel of the Stirling engine was rotated using a the observed drift in ∆T . The grey data from 0 to 48 min Dremel tool and a polishing felt attachment. Even when and from 160 to 190 min were fit to a straight line. The 3 on its lowest setting, the Dremel spins the flywheel too best fit line has a slope of 0.89 10− ◦C/min and is fast to produce reliable results. We found that when shown in Fig. 8 as the dotted− orange× line. The black data the flywheel is rotated too fast, the dominate effect is in the figure show ∆T after removing the drift caused by heating of the top plate by friction. The best results were the changing room temperature. The result is an ap- obtained when the flywheel was rotated at a rate between proximately symmetric ∆T that reaches 0.2 ◦C after 1 and 2 Hz. To reduce the speed of the Dremel, it was running the Dremel for about 30 min. ± 6

IV. SUMMARY 2. Refrigeration

Stirling engines are commonly used as a demonstra- Rotary Tool – We used a Dremel tool that can be pur- tion of a working heat engine. These demonstrations can chased from most hardware stores for approximately $50. be significantly enhanced by adding instrumentation to We had to reduce the supply voltage to the Dremel us- measure the pressure and volume changes of a running ing a variable transformer (variac). Many physics de- engine. We have described a simple and inexpensive way partments will already have such an instrument. They to instrument a high-quality Stirling engine designed for can, however, be purchased from many suppliers. For demonstrations. Our pressure versus time measurements example, they are available from Digi-Key for $370 are high resolution and reveal the 90◦ phase shift be- (https://www.digikey.com/). tween the pressure and volume changes. The P -V di- Digital Multimeter – The thermocouple voltage is agram constructed from these measurements allows one very small and the Vernier voltage probe does not to determine the net work by the engine per cycle and, have the required sensitivity. We measured the ther- therefore, the power output. By adding friction to the mocouple voltage using Keysight’s 34401A multime- flywheel, it was shown that the work per cycle increased ter which is no longer in production. Their replace- while the power output decreased. Finally, we used a ment is the 34461A multimeter which costs $1121 rotary tool to drive the flywheel of the Stirling engine (https://www.keysight.com/). The multimeter we while monitoring the temperature difference across its used was interfaced to a computer running LabVIEW body. A small but measurable temperature difference via a GPIB-to-USB converter made by National Instru- was observed. Reversing the flywheel rotation direction ments. One advantage of the new 34461A multimeter is reversed the sign of the observed temperature difference. that it can connect to a computer directly via USB. Miscellaneous – This measurement also requires thermocouple wires which can be purchased from Omega (https://www.omega.com/) and polystyrene insulation which can be found at most hardware stores. Appendix A: Parts and Suppliers

This appendix provides a list of the equipment required Appendix B: Chopper Design to reproduce the measurements described in this paper. Where appropriate, vendors and cost estimates are also provided. The chopper that was used to track the rotation of the Stirling engine’s flywheel is shown in Fig. 9. As described Stirling Engine – The Stirling engine used earlier, the thick dark band at the top was used to indi- was manufactured by Kontax Stirling Engines cate that the piston was at its peak height. In our work, (https://www.stirlingengine.co.uk/). We pur- this image was printed onto an acetate sheet (overhead chased the KS90 Blue model which costs $118. transparency) and then cut out along the inside and outside circumferences. The chopper was then taped to the flywheel of the Stirling engine using a few small pieces of cellophane tape. If printed without any scaling, the im- 1. P-V diagram age in Fig. 9 will fit onto the flywheel of the KS90 models of the Stirling engines manufactured by Kontax.

Data Acquisition – Our data was acquired using Vernier’s differential voltage probe ($54), photogate ($62), LabQuest interface ($203), and Logger Pro soft- REFERENCES ware ($339) (https://www.vernier.com/). These items are standard in many university and college physics de- 1R. D. Spence and C. L. Foiles, “Stirling engines for demonstra- partments. We also note that the data acquisition could tion,” Phys. Teach. 20, 38 (Jan. 1982). 2H. Richard Crane, “The Stirling engine–173 years old and run- also be implemented using alternative systems, like those ning,” Phys. Teach. 28, 252 (Apr. 1990). offered by PASCO (https://www.pasco.com/) and Na- 3Christopher G. Deacon, Richard Goulding, C. Haridass, and tional Instruments (http://www.ni.com/). Brad de Young, “Demonstration experiments with a Stirling engine,” Phys. Educ. 29, 180 (1994). Pressure Sensor – The MPXV7002 pressure sensor 4Raul A. Simon, “Stirling’s cycle and the second law of thermo- made by NXP Semiconductors was purchased from Digi- dynamics,” Am. J. Phys. 52, 496 (Jun. 1984). Key for $19 (https://www.digikey.com/). The pres- 5Alejandro Romanelli, “Alternative thermodynamic cycle for the sure sensor requires 5 V DC to operate. We used the Stirling machine,” Am. J. Phys. 85, 926 (Dec. 2017). 6A website describing the project by Hiroko Nakahara is avail- E3630A power supply from Keysight, but any DC power able at the following url: https://people.ok.ubc.ca/jbobowsk/ supply can be used. Stirling/introduction.html. 7 Miscellaneous – Also required are a kettle to boil water, D. V. Schroeder, An Introduction to Thermal Physics (Pearson Education, Inc., San Francisco, 1999), pp. 122–136. a coffee mug, a syringe, epoxy (J-B Weld is a good and 8We purchased the KS90 Blue Stirling engine from Kontax Stir- inexpensive option), and rubber tubing. Wires also need ling Engines. These engines can be purchased online from the to be soldered to the pressure sensor. company website: https://www.stirlingengine.co.uk/. 7

Figure 9. The chopper design used to monitor the ﬂywheel rotation. The outside diameter of the chopper is 4 inches.

9MPXV7002 Integrated Silicon Pressure Sensor On-Chip Sig- 10J. W. L. K¨ohler,“The Stirling refrigeration cycle,” Sci. Am. 212, nal Conditioned, Temperature Compensated and Calibrated. 119 (1965). MPXV7002. Rev. 4. NXP Semiconductors. Mar. 2017. 11Ray Radebaugh, “Cryocoolers: the state of the art and recent developments,” J. Phys.: Condens. Matter 21, 164219 (2009).