# The Fluidyne engine

###### Abstract

The Fluidyne is a two-part hot-air engine, which has the peculiarity that both its power piston and displacer are liquids. Both parts operate in tandem with the common working gas (air) transferring energy from the displacer to the piston side, from which work is extracted. We describe analytically the thermodynamics of the Fluidyne engine using the approach previously developed for the Stirling engine. We obtain explicit expressions for the amplitude of the power piston movement and for the working gas temperatures and pressure as functions of the engine parameters. We also study numerically the power and efficiency of the engine in terms of the phase shift between the motions of piston and displacer.

###### pacs:

05.70-a, 88.05.De## I Introduction

The Fluidyne engine was invented by Colin West in .Colin () It is a Stirling machine with one or more liquid pistons. Usually it contains air as a working gas, and either two liquid pistons or one liquid piston and a displacer. The Fluidyne engine operates at a low frequency, typically Hz, and close to atmospheric pressure. The most common application of the liquid-piston system is in irrigation pumping, particularly for irrigation or drainage pumping in places where electric power may not be available.Grupta () Nowadays, many commercial setups are available for specific applications in agriculture, building services, drinking water and sanitation.Tom ()

The basic principle of a Fluidyne engine is the fact that air expands when heated and contracts when cooled. It is possible to find many interesting videos that present the construction and operation of this type of engine for demonstration and teaching purposes, on the internet. However it is not easy to find theoretical material that explains its thermodynamics without an overly complicated technical discussion of the Stirling cycle. In the present paper we study the thermodynamics of a particular type of Fluidyne engine that can be seen as a simple Stirling engine with one free liquid piston and displacer.Walker (); Darlington (); Grupta ()

Recently, we developed an alternative theoretical modelalternativo () for the usual Stirling cycle.Zeman () The main characteristic of that approach is the introduction of a polytropic process,romanelli () for which , as a way to represent the exchange of heat with the environment. The assumption of a polytropic process allows to model different dependencies of the working gas temperature with its volume. This means that the polytropic index can be used as an additional degree of freedom that could be adjusted to the experimental data of a real operating engine. In particular, we remark that all the discussion presented remains unchanged if is replaced by because this last is a special case of the model.

Our alternative model provides analytical expressions for the pressure, temperatures of the working gas and the work and heat exchanged with the heat reservoirs. The theoretical pressure-volume diagram achieved a closer agreement with the experimental one than the standard analysis. Due to the generality of the analytical expressions obtained, they can be adapted to any type of Stirling engine. In the present paper we use the mentioned model to study the thermodynamics of the liquid-piston Stirling engine, “the Fluidyne engine”.

The paper is organized as follows. In the next section we present the liquid-piston Stirling engine and we study in detail the dynamics of the liquid piston with the help of the results of the thermodynamic model. In the third section we obtain the pressure-volume and the work-efficiency diagrams for the Fluidyne engine. In the last section we present the main conclusions.

## Ii Liquid piston Stirling engine

The Fluidyne engine consists of two U-tubes partially filled with liquid and connected with a tube of negligible caliber as shown in Fig. 1. One end of the tube containing the piston liquid is open to the atmosphere. The three connected sections also contain the working gas. Figure 1 shows a snapshot picture of the engine working with its two liquids in motion modifying the gas volume. In this figure the liquids in the left and right U-tubes are respectively the power piston and the displacer of the engine. While the gas volume changes due to the motion of the liquids it maintains a uniform pressure determined by the action of the heat reservoir, represented by the torch in the figure. Two zones can be distinguished by their temperature in the volume occupied by the gas and labeled with the numbers and in Fig. 1. Zone is the hot zone where the gas is in contact with the reservoir at the external temperature , and zone is the cool zone where the gas is in contact with the reservoir at the external temperature . Then, the working gas is never completely in either the hot or cold zone of the engine. However, the upper connection between the zones imposes the same pressure in both. (Note that the diameter of the upper connection is relatively small but still very large compared with the mean free path of the gas molecules, this condition ensures that there is no effusion process involved.)reif () Consequently, the gas density must be different in each zone to keep the same pressure with different temperatures. The left side of the liquid piston is the free zone, here the pressure is the atmospheric pressure and the temperature is also .

In order to start the engine, we must set in motion externally both the piston and the displacer. One way to produce this initial motion starting from the static situation, where both liquids are at rest, is that an external agent performs the following three steps: first it rotates the engine a small angle with respect to an axis perpendicular to the main plane of the engine (the plane of Fig. 1), second it reverses the rotation to return the tubes to the initial position and third it temporally manipulates the pressure over the left side of the liquid piston in such a way as to obtain the desired relative movement between displacer and piston and simultaneously going back to the pressure . After these steps the Fluidyne engine is in the situation shown by Fig. 1, where both liquids are moving inside their U-tubes, and the only external agents that interact with the gas are the heat reservoirs.

With the liquids in motion, the mechanism of energy transfer between the hot and cold zones works without any type of one-way valves and can be qualitatively understood as follows. Initially we focus on the movement of the displacer and suppose the piston still. When the displacer is at the center of the U-tube, suppose that the air volumes are more or less the same in the hot and cold zones. When the displacer liquid moves towards the cold (hot) zone air is displaced, through the upper connecting tube, towards the hot (cold) zone and a greater quantity of air increases (decreases) its temperature and consequently air pressure increases (decreases). Then the air pressure of the right side of the liquid piston (see Fig. 1) depends on the displacer’s motion but the left side remains at constant atmospheric pressure. If we now incorporate the motion of the piston, the net motion of air between zones 1 and 2 depends on the phase of the relative motions of the displacer and piston fluids, which results in a cyclical transfer of energy from the high temperature zone 2 to the piston liquid of zone 1 (with a small amount going into sustaining the oscillation of the displacer liquid against dissipative losses). The phase of this motion is a crucial parameter for the efficiency and power delivered by the engine, and is treated in the analysis that follows.

When the displacer liquid is set into oscillation in its U-tube, then the gas above this liquid is transferred back and forth between the hot and cold zones. Both sides of the displacer tube have the same pressure, even when the engine is working, therefore if the displacer behaves as an ideal fluid any initial oscillation never decreases. In what follows we assume that the displacer behaves as an ideal fluid subject to the restoring force of gravity, and any oscillation has a natural frequency given by

(1) |

where is the gravity acceleration and is the displacer’s length. Then the amplitude of the oscillation is given by

(2) |

where is the maximum amplitude and we take a vanishing initial phase. If the displacer is not an ideal fluid it is always possible to use part of the energy produced by the engine to compensate losses and to maintain the original oscillation. As we shall see will be the operating frequency of the engine.

To study the dynamics of the power piston we observe that the pressures on the right side and left side give rise to a net force on the liquid piston proportional to . We also consider an external force that dissipates all the useful power delivered by the working gas. In a real engine the power piston should be coupled to a crankshaft to transfer the useful power. We opt for this simplified model in order to obtain a closed analytical solution, however the real system can be addressed numerically. Therefore, the power piston dynamics is determined by Newton’s second law applied to the displacement of the liquid piston (see Fig. 1)

(3) |

where , and are the piston liquid density, cross-sectional area and length of fluid respectively and is a damping coefficient.

To obtain the working gas pressure we assume the following simplifying hypotheses for the working gas: (i) The gas behaves as a classical ideal gas. (ii) The gas has uniform temperature in each zone namely and . (iii) The mass of gas inside the thin connecting tubes is negligible, i.e. all the gas is inside the U-tubes where the displacer and the piston move. (iv) The expansion and compression processes are treated as polytropic processes. For an ideal gas with a constant number of moles, the polytropic process may be defined through the relation , where is the polytropic index with typical values such that with for dry air at room temperature. In such a context we have obtained in Ref.alternativo, analytical expressions for the pressure and temperatures of the working gas as functions of the gas volume, the external temperatures and the polytropic index which characterizes the heat absorption process. The gas pressure is given by

(4) |

and the temperatures and are

(5) |

(6) |

where is the ratio between the temperatures of the reservoirs

(7) |

is the gas volume of the cool zone with temperature , is the gas volume of the hot zone with temperature , is the total gas volume, , and are initial conditions of , and respectively, and is the initial condition of . It is interesting to underline that Eqs. (4, 5, 6) for the isothermal case coincide with the so-called Schmidt solution, published in .Schmidt () For the case the same equations describe the classical adiabatic solution.Bercho (); Formosa () Both examples show the ductility of polytropic processes to describe the thermodynamics of this type of engines.

Assuming that both U-tubes have the same cross section area , the volumes of the zones and can be expressed as (see Fig. 1)

(8) |

(9) |

where the displacements and are taken from the initial (equilibrium) position, i.e. give the initial volumes in Eqs. (8, 9). The total volume is then

(10) |

The explicit dependence of on and is obtained from Eqs. (4), (8), (9) and (10), then

(11) |

It is clear now that Eq. (3) has a non-linear dependence with both and . However the engine operates in a closed regenerative cycle, which presupposes only one characteristic frequency. Then, only the periodic solutions of Eq. (3), with as a characteristic frequency, will be useful. Here and we further assume that the engine geometry is such that both and are satisfied, and then these magnitudes can be treated as perturbations in Eq. (11). Expanding equation Eq. (11) around the unperturbed volumes and substituting in Eq. (3) the following linear equation is obtained

(12) |

where

(13) |

(14) |

and is given by Eq. (2).

Equation (12) describes a linearly damped oscillator with external forcing. The first term on the right-hand side produces a constant shift from the equilibrium position, which is equivalent to a redefinition of the volume . Therefore, from now on we assume in order to ignore it and we will concentrate on the periodic solution.

When the engine achieves the steady state motion, the solution of Eq. (12) that determines its dynamics is

(15) |

where

(16) |

and

(17) |

From Eq.(16) it is clear that the engine does not work if , that is, in order to function the engine needs heat reservoirs with , see Eq.(7). Moreover, when then tends asymptotically to a finite value. On the other hand, if in Eq.(17) the power piston is in resonance with the displacer, , and the amplitude is maximum.

## Iii Fluidyne thermodynamics

To obtain the pressure-volume diagram for the Fluidyne engine we use Eqs. (2), (11) and (15). This is shown in Fig. 2 for two characteristic values of that are to be explained below. The area inside such smooth closed curves represents the total work of the cycle, whose expression is

(18) |

where the last equality follows from Eq. (10). Here is the work available for overcoming mechanical friction losses and for providing useful power.

The differential equation for the heat absorbed by the working gas in the polytropic process was obtained in the general case of a Stirling engine in Ref. alternativo, as follows:

(19) |

where , , are defined by Eqs. (8), (9), (10) and

(20) |

is the quotient of the specific heats of the gas at constant pressure and constant volume . In this paper we take the polytropic index such that .

It is important to emphasize that Eq. (19) together with Eqs. (4), (5) and (6) determine the Fluidyne thermodynamics. These four equations depend on the parameter (the temperatures ratio) and they have finite asymptotic values when . This means that after a certain value of , no matter how much we increase the temperatures ratio, the pressure and the absorbed heat are bounded, and this in turn explains why the useful work and the efficiency in any Stirling engine are asymptotically bounded.

Let us call the heat absorbed by the gas in the cycle with the convention ; similarly we call the heat rejected. They can be calculated numerically using Eq. (19) and the definitions

(21) |

(22) |

where and refer to the paths where and respectively. As the internal energy change in the entire cycle vanishes, the total work verifies . Therefore, the efficiency defined as can be written as

(23) |

We have integrated numerically Eqs.(18), (21) and (22) using the standard Simpson’s rule. Figure 3 shows the relation between work and efficiency as functions of the phase . From this figure is clear that the phasing between the piston and the displacer plays a central role in the functioning of the engine; in particular we observe that the maximum efficiency phase does not coincide with the maximum work phase. The diagrams in Fig. 2 correspond precisely to the phases of maximum work and maximum efficiency.

When the piston is in resonance with the displacer and this situation is intermediate between maximum work and maximum efficiency; there the engine is in a good performance zone, where the engine operates in a good compromise between power and efficiency.

Figure 4 shows the behavior of work and efficiency when the polytropic index varies. This parameter depends on the gas initial thermodynamic conditions and in particular on its relative humidity.romanelli () The figure shows that the dependence with is stronger for the power than for the efficiency.

## Iv Conclusions

Research on Stirling engines is one of the lines that contribute both to the rational use of energy and to sustainable development. In particular the solar thermal conversion systems based on these engines are amongst the most interesting and promising research lines. tesis1 (); tesis2 (); tesis3 (); tesis4 (); tesis5 (); tesis6 () This paper presents an unusual application of heat engine thermodynamics combined with the physics of oscillators at a level appropriate for advanced undergraduate students.

From a strictly technical point of view we have proposed a theoretical model that describes the thermodynamics of a Stirling engine in a simple, precise and natural way.alternativo () The engine studied in this paper is a special type of the Stirling engine, therefore we have applied the theoretical model to investigate its thermodynamics and dynamical evolution. The exchange of heat with the environment is modeled as a polytropic process. We obtain analytical expressions for the piston amplitude, the working gas temperatures and pressure. These quantities are expressed as functions of (i) the phase difference between the power piston and the displacer, ; (ii) the ratio of the temperatures of the heat reservoirs, ; (iii) the polytropic index, and (iv) the geometrical characteristics of the engine. We also study numerically the engine power and its efficiency as a function of and and show that, when in resonance, the engine works a good performance regime. These results show the versatility of the thermodynamic model developed here to describe different types of Stirling engines and they encourage further work on this line.

In summary, the theoretical approach proposed describes analytically the Fluidyne engine thermodynamics and it contributes to an understanding of the subtle way through which this engine transforms heat into work.

## Acknowledgments

I acknowledge the stimulating discussions with Víctor Micenmacher and Raul Donangelo and the support from ANII and PEDECIBA (Uruguay).

## References

- (1) C.D. West, Liquid Piston Stirling Engines, (Van Nostrand Reinhold, New York, 1983).
- (2) A. Gupta, S. Sharma, S. Narayan, Liquid Piston Engines, (John Wiley and Sons, Hoboken, 2017).
- (3) T. Smith, C. Markides, Thermofluidics Ltd., https://www.thermofluidics.co.uk/.
- (4) G.Walker, J.R.Senft, Free Piston Stirling Engines, (Springer-Verlag, Berlin, 1980).
- (5) R. Darlington, K. Strong, Stirling and Hot Air Engines, (The Crowood Press, Marlborough, 2005).
- (6) A.Romanelli, “Alternative thermodynamic cycle for the Stirling machine,” Am. J. Phys. 85, 926-931 (2017).
- (7) M.W.Zemansky, R.H.Dittman, Heat and Thermodynamics, (McGraw-Hill, London, 1997).
- (8) A.Romanelli, I.Bove, F.González, “Air expansion in a water rocket,” Am. J. Phys. 81, 762-766 (2013).
- (9) F.Reif, Fundamentals of statistical and thermal physics, (McGraw-Hill, New York, 1965).
- (10) G.Schmidt, “Theorie der Lehmann’schen kalorischen Maschine,” Z. Ver. Dtsch. Ing. 15, 1-12 (1871).
- (11) I.Urieli, D.Berchowitz, Stirling Cycle Engine Analysis, (International Public Service, 1983).
- (12) F.Formosa, G.Despesse, “Analytical model for Stirling cycle machine design,” Energy Convers. Manage., 51, Issue 10, 1855-1863 (2010).
- (13) L.G. Thieme, S. Qiu, M.A. White, “Technology development for a Stirling radioisotope power system,” NASA TM-2000-209791 (2000).
- (14) R.G. Lange, W. P. Carroll, “Review of recent advances of radioisotope power systems,” Energy Convers. Manage., 49, 393-401 (2008).
- (15) Caleb C. Lloyd, “A low temperature differential Stirling engine for power generation,” Master of Engineering thesis, University of Canterbury (2009).
- (16) D. García Menéndez, “Desarrollo de motores de Stirling para aplicaciones solares,” PhD thesis, Universidad de Oviedo (2013).
- (17) I. M. Santos Ráez, “Estudio de un motor Stirling con absorbedor interno alimentado con energía solar,” PhD thesis, Universidad de Málaga (2015).
- (18) G. Barreto, P. Canhoto, “Modelling of a Stirling engine with parabolic dish for thermal to electric conversion of solar energy,” Energy Convers. Manage., 132, 119-135 (2017).