SCHMIDT THEORY FOR STIRLING ENGINES

Tentative version on January 20, 1997

KOICHI HIRATA

e-mail: khirata@gem.bekkoame.ne.jp


1. INTRODUCTION

The Schmidt theory is one of the isothermal calculation methods for Stirling engines. It is the most simple method and very useful during Stirling engine development.
This theory is based on the isothermal expansion and compression of an ideal gas.


2. ASSUMPTION OF SCHMIDT THEORY

The performance of the engine can be calculated a P-V diagram. The volume in the engine is easily calculated by using the internal geometry. When the volume, mass of the working gas and the temperature is decided, the pressure is calculated using an ideal gas method as shown in equation (1).

(1)

The engine pressure can be calculated under following assumptions:
(a) There is no pressure loss in the heat-exchangers and there are no internal pressure differences.
(b) The expansion process and the compression process changes isothermal.
(c) Conditions of the working gas is changed as an ideal gas.
(d) There is a perfect regeneration.
(e) The expansion dead space maintains the expansion gas temperature - TE, the compression dead space maintains the compression gas temperature - TC during the cycle.
(f) The regenerator gas temperature is an average of the expansion gas temperature - TE and the compression gas temperature - TC.
(g) The expansion space - VE and the compression space - VC changes according a sine curves.

Table 1 shows symbols used the Schmidt Theory.

Table 1 Symbols


3. ALPHA-TYPE STIRLING ENGINE

Figure 1 shows the calculation model of Alpha-type Stirling engine.

Fig. 1 Alpha-type Stirling Engine

The volumes of the expansion- and compression cylinder at a given crank angle are determined at first. The momental volumes is described with a crank angle - x. This crank angle is defined as x=0 when the expansion piston is located the most top position (top dead point).
The momental expansion volume - VE is described in equation (2) with a swept volume of the expansion piston - VSE, an expansion dead volume - VDE under the condition of assumption (g).

(2)

The momental compression volume - VC is found in equation (3) with a swept volume of the compression piston - VSC, a compression dead volume - VDC and a phase angle - dx.

(3)

The total momental volume is calculated in equation (4).

(4)

By the assumptions (a), (b) and (c), the total mass in the engine - m is calculated using the engine pressure - P, each temperature - T , each volume - V and the gas constant - R.

(5)

The temperature ratio - t, a swept volume ratio - v and other dead volume ratios are found using the following equations.

(6)

(7)

(8)

(9)

(10)

The regenerator temperature - TR is calculated in equation (11), by using the assumption (f).

(11)

When equation (5) is changed using equation (6)-(10), the total gas mass - m is described in the next equation.

(12)

Equation (12) is changed in equation (13), using equation (2) and (3).

(13)

Now;

(14)

(15)

(16)

The engine pressure - P is defined as a next equation using equation (13).

(17)

The mean pressure - Pmean can be calculated as follows:

(18)

c is defined in the next equation.

(19)

As a result, the engine pressure - P, based the mean engine pressure - Pmean is calculated in equation (20).

(20)

On the other hand, in the case of equation (17), when cos(x-a)=-1, the engine pressure - P becomes the minimum pressure - Pmin, the next equation is introduced.

(21)

Therefore, the engine pressure - P, based the minimum pressure - Pmin is described in equation (22).

(22)

Similarly, when cos(x-a)=1, the engine pressure - P becomes the maximum pressure - Pmax. The following equation is introduced.

(23)

The P-V diagram of Alpha-type Stirling engine can be made with above equations.


4. BETA-TYPE STIRLING ENGINE

Similarly, the equations for Beta-type Stirling engine are declared. Figure 2 shows a calculation model of a Beta-type Stirling engine.

Fig. 2 Beta-type Stirling Engine

The expansion momental volume - VE and the compression momental volume - VC are described in the following equations, with a swept volume of a displacer piston - VSE, a swept volume of a power piston - VSC and a phase angle -dx between the displacer piston and power piston.

(24)

(25)

In the case of the Beta-type Stirling engine, the displacer piston and the power piston are located in the same cylinder. When both pistons overlaps there stroke, an effective working space is created. The overlap volume - VB in equation (25) can be calculated in the next equation.

(26)

Then the total momental volume - V is found in equation (27).

(27)

The engine pressure - P based the mean pressure - Pmean, the minimum pressure - Pmin and the maximum pressure - Pmax are described in the following equations like the Alpha-type Stirling engine.

(28)

Several ratios and coefficients are defined as follows.

(29)

(30)

(31)

(32)

(33)

(34)

(35)

(36)

(37)

(38)

The P-V diagram of Beta-type Stirling engine can be made with above equations.


5. GAMMA-TYPE STIRLING ENGINE

Figure 3 shows a calculation model of a Gamma-type Stirling engine.

Fig. 3 Gamma-type Stirling Engine

Similar calculation equations are made as the Alpha- and Beta-type engine. The expansion momental volume - VE and the compression momental volume - VC are described in the following equations with a swept volume of a displacer piston - VSE, a swept volume of a power piston - VSC and a phase angle - dx between the displacer piston and the power piston.

(39)

(40)

The total momental volume - V is described a next equation.

(41)

The engine pressure - P based the mean pressure - Pmean, the minimum pressure - Pmin and the maximum pressure - Pmax are found in the following equations.

(42)?

Now,

(43)

(44)

(45)

(46)

(47)

(48)

(49)

(50)

(51)

The P-V diagram of Gamma-type Stirling engine can be made with above equations.


6. INDICATED ENERGY, POWER AND EFFICIENCY

The indicated energy (area of the P-V diagram) in the expansion and compression space can be calculated as an analytical solutions with use of the above coefficients. The indicated energy in the expansion space (indicated expansion energy) - WE(J), based on the mean pressure - Pmean, the minimum pressure - Pmin and the maximum pressure - Pmax are described in the following equations.

(52)

The indicated energy in the compression space (indicated compression energy) - WC(J) are described in the next equations.

(53)

The indicated energy per one cycle of this engine - Wi(J) is demanded in the next equations.

(54)

Relations between Pmean, Pmin and Pmax are determined in the following equations.

(55)

(56)

The indicated expansion power - LE(W), the indicated compression power - LC(W) and the indicated power of this engine - Li(W) are defined in the following equations, using the engine speed per one second , n(rps, Hz).

(57)

(58)

(59)

The indicated expansion energy - WE found equation (52) means an input heat from a heat source to the engine. The indicated compression energy - Wc calculated by equation (53) means a reject heat from the engine to cooling water or air. Then the thermal efficiency of the engine - e is calculated in the next equation.

(60)

This efficiency equals that of a Cornot cycle which is the most highest efficiency in every thermal engine.


7. EXAMPLE OF CALCULATION


EXERCISE:

Make a P-V diagram and calculate the indicated power of an Alpha-type Stirling engine under following conditions.
Swept volume of an expansion piston: 0.628 cm3, swept volume of a compression piston: 0.628 cm3, dead volume of the expansion space: 0.2cm3, dead volume of the compression space: 0.2cm3, regenerator volume: 0.2cm3, phase angle: 90deg, mean pressure: 101.3 kPa, expansion gas temperature: 400degC, compression gas temperature: 30degC, engine speed: 2000 rpm.


A temperature ratio - t, a swept volume ratio - v and other dead volume ratio are calculated with the equation (6) - (10).

Each coefficient is calculated with equation (14) - (16) and (19).

Engine pressure is calculated with equation (20).

When crank angle - x=0deg:

Similarly, when x=10deg:

When x=20deg:

c

Next each momental volume is calculated with equation (2) - (4).
When crank angle, x=0deg:

When x=10deg:

When x=20deg:

c

Repeat above calculation to one complete cycle and plot the volumes - V and pressures - P on a graph paper.
An example of the P-V diagram is shown in Fig. 4.

Fig. 4 P-V Diagram

The indicated energy is calculated with equation (52), (53) and (54).

The indicated power of this engines is calculated with equation (59).

The indicated power of this engine is 0.689 W.


ACKNOWLEDGEMENT

I wish to thank J. H. de Baat for many helpful suggestions during the translation of this script.


REFFERENCES

1) G. Walker., Stirling Engines, (1980),17, Oxford Univ. Press.


Return to Stirling Engiene Home Page Academic Edition
Koichi Hirata
e-mail: khirata@gem.bekkoame.ne.jp (Home)
e-mail: khirata@srimot.go.jp (Office)