Mit
Mit
Mit
by
RAHMATALLAH SHOURESHI
DOCTOR OF PHILOSOPHY
at the
MASSACHUSETTS INSTITUTE OF TECHNOLOGY
June, 1981
A -. 0% 'r
Signature of Author
Department of Mechanical Engineering
IN
Certified by
Thesis Co-supervisor
Thesis Co-supervisor
111
Accepted by
Chairman, Department Committee on Graduate Students
Arcfilven
MASSACHUSETTS INSTITUTE
OF TECHNOLOGY
OCT 29 1981
uBwPlES
ANALYSIS AND DESIGN OF
STIRLING ENGINES FOR WASTE-HEAT RECOVERY
by
Rahmatallah Shoureshi
TABLE OF CONTENTS
Page
ABSTRACT ii
LIST OF TABLES vii.
LIST OF FIGURES ix
ACKNOWLEGEMENTS xii
MOMENCLATURE xiv
CHAPTER I:
CHAPTER II:
CHAPTER IV:
LIST OF FIGURES
ACKNOWLEDGEMENTS
F Correction factor
Lp displacer length
N engine speed
P pressure
AP pressure drop
S stroke
t thickness
T temperature
TC cold temperature
M
TH hot temperature
Tr temperature ratio (TC/TH)
V volume
VD dead volume
Greek letters:
a aspect-ratio
y ratio of Re/M 2
6 temperature ratio
p density
1- viscousity
t time constant
Subscripts:
am ambient
c compressor, cooler
C cooler, cold
for dead volume
e expansion space
EQ equivalent
h heater
H Heater, hot
in input
m mean value
min minimum
max maximum
out output
regenerator
T for temperature
1.1- Introduction
ter a patent on the invention was issued in his name. The eng-
gine originally used hot air as a working fluid and was, there-
from competing with first the steam engine and later with the
engines were made and used in a wide variety of sizes and shapes
in Britain, Europe, U.S.A., and other parts of the world. They were
blood pumps and many other services. Since most of these appli-
source, and low output power level design. As shown in next sec-
tion, there has not been adequate research along these li.nes;-
tive. Much process heat above 200.F was (and still is) rejected
which might have some hope is the Rankine cycle; the following
sections (6-1-1) outlines efforts for using this cycle for waste-
this ratio might become before the engine could no longer gene-
traded off for lower cost, i.e. where the heat source itself is
at lower temperature, and thus use of conventional or less ex-
water and such systems were first promulgated more than 150
years ago. Almost a century ago, the principle was related spe-
drives a gas turbine which provides the plant power. Then the
candidate fluids.
of the cycle showing the waste heat source and the pertinent
portion of this hot acid for the heat recovery system. The
condenser.
steam and Freon, and is projected to recover 500 KW. The steam
are only a few research efforts in that area based on the Stir-
lirg cycle.
ting temperatures [15]. The engine which has been used is the
and continues with ideal cycle analyses (i.e, the Schmidt equa-
between the hot and cold spaces. Chapter V shows the special
kine engines.
Fic. (10)
Analysis of
alternative
4,
Thermofluid analysis
required for using
Determination
of different
forms of S.E. at low AT and alternatives
regenerator for. low power-level for minimum
units friction
V
Usi ng scal ing 1aws .to model
this class of engine
and at low-temperature-ratios.
gine should be more practical and not taken as the ideal case.
entrance to the pump and at exit from the expander, the flow
TI
H2
] Fig. (2-l)
17
is a mixture of liquid and vapor and would create some design
5
Fig. (2-2)
MOMMMahmm
For the absolute temperature-ratio ranging from 1.0 to 2.0
with the maximn engine pressure less than 500 psia a proper
TC = T = 55 F
TH/TC=l. 4
Td = T3' = 260 F
State 1:
T = 55 F
P = 33.23 Psia
H = 21.86 Btu/lb
S = .0477 Btu/lb R
V = 1/71.24 Ft /lb
State 2:
S = .0477 Btu/lb R
P = 413.8 Psia
H = 22.85 Btu/lb
S
State 3:
P = 413.8 Ps ia
T = 260 F
H = 133 Btu/lb
S = .22 Btu/lb R
State 4:
P = 120 Ps ia
H = 122.5 Bt u/lb
S = .22 Btu/lb R
State 5:
P = 120 Psia
H = 151 Btu/lb
S = .264 Btu/lb R
State 6:
P = 33.23 Psia
H = 135 Btu/lb
S = .264 Btu/lb R
Wout
Efficiency = 18.4%
Qin
adiabatic work transfer processes (in ideal case they are isen-
of 85%, the maximum power is only 49% of the power which could be
T>
WonJ 10 4
4
"WOMOMMMOMMAMb
Fig. (2-3)
Qin = CP (Tj - T)
T (TS /T - 1)
T11
C001.0 ItsG.
I1
WT _ _S
Fig. (2-4)
T1
TH TC
TC 4
TH T
S
Fig. (2-5)
losses which seem most important and push the engine farther
TR - TC
TH - TC
L - TC/TH
Efficiency =
I + (1-E)(1-TC/TH)/((k-1)*(ln PgL/P ))
that Stirling and Rankine engines are more efficient than Bray-
that not all elements of the working fluid pursue the same ther-
modynainic cycle. Since the ideal cycle analyses have been pre-
sented in closed form solutions, then they are useful for pre-
process (4-1). Then the gas flows through the regenerator and
gains enough heat to reach the hot temperature, TH. This heat-
ing is such that the same volume of gas which enters to the re-
hot space, i.e. the volume variations of hot and cold spaces
moved from the working fluid and stored for its subsequent
Where
A =\f" /[TC/TIH+VC/VH+4X*TC/(TC+TH)]
B = (TC/TH) +(VC/VH) + 2TC*VC*COSC7/(TH{*VH)
resulted.
Instantaneous mass:
Hot Space:
MH = 1/2VH*Pm*(l-A ) *(l+COSwt)/[R*TH*(l+ACOS(wt- 6 )]
Cold Space:
MD = X*VH*Pm(l-A 1 )/[R*TD(l+ACOS( w t- )]
Total Mass:
VH*Pm(l-A )
MT CO
TC/TH+2X*TC/(TC+TH)+VC(1+COSq)/(2VH)]
R*TC(l+A*COSG)
ters. To get the optimum dead volume ratio (X VD/VH) the deri-
quadratic equation:
ature ratio and swept-volume ratio. Walker [ 25] has found the
Fig. (1-4) shows his results for two different temperature ratios,
The above equation would give the same results directly and
quickly.
engines.
30
less than the heat which was extracted from the gas by the re-
way to the hot space, the gas leaves the regenerator with a
temperature less than the hot temperature of the cycle (TH) and
and hot space in real engines. Fig. (15) shows a schematic diag-
ram of a real engine, for this engine based on hot and cold tem-
tor as follows.
TR = (TH-TC)/(ln TH/TC)
1 - TC/TH
Efficiency =
1+Cv/R*(1-E)(1-TC/TH)/(lnP, /P )
Efficiency
(E-16A)
Ideal Eff. l+(l-E)/(k-1)(l-TC/TH)/(lnP3
ten. These dead volumes would push the engine farther fron ideal
heater and gas cooler will not change the work realized per
cycle but will definitely increase the heat required per cycle.
'
Where:
t
TR/TH = (1-TC/TH)/(ln TH/TC) ,
(E-16B)
increasing dead volume of the system for case A=2 and different
generator and the presence of dead volume are for steady state
cases and are only first order, nevertheless they show how the
III-- -- -- M
32
real engine trends farther and farther from ideal case. Also
that not all elements of the working fluid pass through the
take into account the movement of working fluid into and out of
each one. Since the volume in the heat exchangers (dead volume)
analytical techniques.
In the ideal cycle there are only three components, i.e. there
position. Since the free flow area inside the hot and cold
components, then the pressure drop inside these spaces are neg-
exchangers.
35
nent and let L denote the length of this component. The pres-
tion factor fx as
this point that the hydraulic diameter d and the free flow area
mx = fx*AFR*Vx
=d(P) =1/2*fx*(L/d)*(mx)/(fx*AFR*AFR)*d(x/L)
x = Px/(R*Tx)
d(AP)I= 1/2*fx*(L/d)*(mx*R*Tx)/(AFR*AFR*Px)*d(x/L)
I I= 1/2 ffx(L/d)*(mx*R*Tx)/(AFR*AFR*Px)*d(x/L)
For Cooler:
2.
I&Pj= 1/2(L/d) *(R*TC)/(AFRC*AFRC) fx*mlx/Px*d(x/LC) (2-1)
For Heater:
For Regenerator:
Based on the given geometry, the free flow area (AFR) and other
fx = a/(Re)**b
cess, a time averaged mass flow rate can be used for each heat-
ten as:
For Cooler:
outlet conditions.
This approach can be used here, i.e. the effect of the pres-
VC = VC+ A VC
with AP
PC = PH+AP
WC PC*d(VC)= (PH+4dP)d(VC+AVC)
WC =fPH*d(VC)+ P*d*VC
pressure loss, and the second term is the power loss due to the
Wloss = AP*d(VC)
This integration can be done either by computer or by some ave-
cold space would push the engine farther from the ideal case
and generate entropy inside the system. Those losses which are
heat transfer loss, heat leakage from the system, transient heat
and heat pumping loss inside the gap between cylinder and pis-
sented.
This means that the temperature of the gas entering the hot and
cold spaces (cylinders) and the regenerator will not be the same
temperature.
39
equal to the work done by the gas in the adjacent space (cylin-
40
regenerator the temperature is Td. For a perfect regenerator
Ta = Td, furthermore AT& = Tc - Tb. The temperature of the
(Ts - Ta)*EXP(-NTU) = Ts - Tb
(Ts - Tc)*EXP(-NTU) = Ts - Ta
&T Ts - Tb = ( AT Ec )/EEXP(2*NTU)-1 =
TC = TC * [l-Wc*(k-l)/(mc*R*TC*k)/[EXP(2*NTUc)-l]] (2-8)
Effective Temperature in Hot Space:
loss which is taken out in the compressor. This means that the
Then by using equations (2-8) and (2-9) the following are resul-
ted.
VittA~kr"C-t if I e-1
Which means%
Wh*Wc*(k-1)/(mh*R*TH*k)/[EXP(2*NTUh)-l] (2-10)
temperature TH, it releases heat such that by the time that the
and the regenerator gains all of the heat. Then similarly, when
tor such that when the flow leaves the regenerator it should
E = QR/[m*Cp(TH-TC)] (2-12)
is used which results in the same amount that the Rios [22]
ture gradient between the two ends of the Stirling engine, then
the warm side of the system (heater) to the cold side (cooler).
solid material and also the gas inside the void volume of the
tance is now the contact between adjacent wires. Since the wires
Piston or Displacer
at T<TH
Cylinder Wall at
T4-yn
rTH
Piston or Displacer
at TH
When the displacer (piston) moves toward the top dead cen-
heat conduction.
Q = TT/8*Kg*S*(TH-TC)*(B/t)*(S/L)*(BET) (2-16)
C Kp/ Kg * * W/(2*oc p)
Where B and S are the bore and stroke of the hot cylinder, Kp
assumed.
Appendix ( F ) shows the stress analysis for the calcula-
t = (Pm*R)/5000 (2-17)
2 J7(TH - Tam)
OL (2-18)
1/(L*R*hi)+ln ((R+t)/R)/(L*K )+l/EL*(R+t)*ho
get:
7T(TH - Tam)
QL (2-19)
1/(3S*Kg)+Ln ((R+t)/R)/(S*K)+l/(3.5S*(R+t))
tant value, the mean (or constant) vaue can be used in equa-
tion (2-19) for heat leakage calculation.
the gas in the center and near the wall of the cylinder produces
a neat flow from the wall to the center and vice-versa, which
Their solution is used for the present work in the next section.
It can be summarized as
Pm = mean pressure
= VH/VC
been obtained.
An alternative for calculating this loss would use experi-
coefficient between the gas inside the cylinder and the walls.
bers which depend on the geometry of the flow system and on the
following form.
Where Tm and Pm are mean values for the flow temperature and
for heat transfer between the gas and the cylinder wall, then
dix ( G ) .
heat transfer between the gas and cylinder wall. The result be-
comes:
8 4
Q= .052 (Re) (Pr) *(k-1l)/k*Pl/Pm*COS(>1*Kg*Tm*Xm (2-24)
from the system into it' Thu Ite gas inside behaves like infinite
number of thermal capacitances with thermal resistances in bet-
shows circuit graph and bond graph of such a gas cylinder. The
Distributed To Lumped
Sub-layer Model Model
Tw/Tm=1/ (l+ S)
Substitution of io for S shows that TO/Tm=(l+1+T*iw)/(l+Tg*iw)
TN/TO=((l+ TS)/(l+T
0 S))N . Although the N-lump solution is more
have: P k/1-k=Const.
(2-25)
By taking the logarithmic derivative of equation (2-25) respec-
tively we get:
M- - -M ____- M-
temperature.
be rewritten as
calculate the basic heat input and the output work, would intro-
This means that in this way the transient heat transfer loss can
fer in the radial direction, but it can be used for axial direc-
tion as well. The same r-egion i.e.; sub-layer which exists on the
cylinder wall is present on the piston face too. The longer stroke
produces a thicker boundary layer which results to a higher tem-
2.3.2.7- Heat Pumping Loss Inside the Piston and Cylinder Gap:
Gap eV
the engine cylinder and the outer diameter (OD) of the piston
out of this gap. Since the closed end of the gap is cold,
extra heat must be added to the gas as it comes back from this
(2-38)
Where B is the bore, I is the gap size, Lp is the length of the
loss. Martini [13] has arbitrarily used 20% of the basic output
ship between fmep and engine speed, then the following corre-
fmep=.002N+l (2-29)
Then the power loss due to this friction effect, based on the
P=fmep*V*N*n/531
volume.
ted from the Schmidt power. Various heat losses are calculated
and added to the Schmidt heat input. All these engine processes
der methods are much more laborous, but since fewer assump-
to be more accurate.
Rios [22] second order model. Similarity of this model and the
idealized.
55
form distribution.
areas, but up to the present time they have not been of prac-
tical importance.
subsequent corrections.
matrix material which has a heat capacity higher than the gas
gas enters the cylinder from the heat exchanger. Generally, the
free flow area of the heat exchanger will be less than that of
when the jet of gas enters from the heat exchanger into the
helium and hydrogen are the best gases to use as working fluids
that the mass contained in the dead spaces may vary only with
meters are:
vcEVC/VC (2-33)
vh5Vh/VH (2-34)
hot cylinder.
P=p/Pmax
-=k *t*'c +A
-PO (]/r*Vc+W+ksV) (2-35)
dP,~ ~(rA44L/v
/4, > 0 >ka*4L
yrt fMcA
t4+k * p d (2-36)
(2-37)
VGi~ct
XPxr
d Mc~a, d M14 )e
M /(2-38)
The differential equations for the mass are:
a VC v.+aA (2-39)
AAWC >0
4c 40
tance of gas and solids (C), and fluid and thermal resistances
is taken for the case assuming a given form of the volume varia-
that one can pass from a first order analysis to higher orders
very easily simply by adding more bonds and elements for more
M____
inertance is used for the gas motion inside the cylinder. There
gas state equation such as the ideal gas law (P*V=MRT). Besides
derived.
0 M- M___
For the present work, where the forcing function is the volume
62
2.4.2-Real Engine Model:
from the engine are the effective hot and cold temperatures;
factor 1(TH/TC)/(TH TC )]
In order to calculate the real engine performance, the
.T"
J
~ {j (2-43)
I
PMRTG - ' f I --am 1, 1Lwt --.J (wt)
(2-44)
UTEGRA
11W"-,10,(2-44)
Re<2000 F=16/Re
Re>2000 F=.0457/Re.2
In regenerator-
Reg60 F=53.7/Re. 9 3
60<Re<l000 F=5.176/Re.3 6 5
Re>1000 F=1.035/Re. 1 2 5
Figures (23) and (24) show the data for thermal conduc-
tivity and Prand:il number with their assigned correlations.
Table (6) and Fig. (25) show the Ward pertormance aata
and the complete model predictions. There is a small dif-
ference between the experimental data and the model results.
This error would appear to come from two sources: (1) the
first one results from having no access to the real geometry
of heat-exchangers for the experimental engine; (2) the
second reason is that the cold temperature of the engine is
not specified so that 70 F (530R) has been used for TC. Since
the experimental data does not give the amount of input heat,
then in Fig. (26) the calculated input heat is plotted versus
the measured fuel mass flow rate. As shown, this plot is a
straight line which means its slope is the heat value of the fuel.
a phase shift of 112 between the cold and hot space displace-
ment. In the second one, this phase difference has been chang-
ed to 118 . In Figs,(28) and (30) the Carnot efficiency (corres-
put per cycle with increased speed, and the change in perfor-
mance with phase angle difference of the hot and the cold
cision.
Hot heat-exchanger:
Regenerator:
c- GPU-3 Engine:
ferent models. Two of the last model GPU-3 were preserved and
are now being used by NASA Lewis Research Center to obtain re-
TH=1860 R
TC=560 R
Displacement volume=2.175 in 3
Length =.89 in
Diameter =.89 in
Wire diameter =.0016 in
Mesh =213 in
Porosity =.714
Lewis Research Center (LeRC) has used a third order analysis and
model results in Tables (9) and (10), and Figs.(31) and (32).
the present complete model from the output power view point.
On the other hand, their heat input is more than what the com-
the engine, such as dead volumes inside the system; and that
mation.
M-
75
CHAPTER III
engine performance.
engines and its results are very close to the second order
late basic heat input and work output first for an engine with
sion processes inside the cylinders, but the real engine works
tors to the Schmidt basic work output equation (see section 2.2),
form:
it/
A Pm A *119 S109 F(I-Tc/To)/(,+ x'to) (3-1)
(3-1) Pm, VH, <p, and TC/TH are the terms which determine the
(TC/TH) the complete model has been run where helium as the
solution, the error between that solution and the Schmidt equa-
(3-3) we get:
ted and compared with the results of equation (3-2) and the
Y=.36167x-87.4476 (3-5)
Where W-W
Y--- ----*100
Wreal
and x=360+$
been varied, this has been done for different sets of constant
($)and (TC/TH) ard rigs.(35) and (36) show the results of F,(D)
versus the total dead volume in the system. These two figures
onto curves indicates that no other terms than F,(T), F(*) ,and
79
F(D) would appear to have any significant effect.
Following are linear correlations which has been used for F,(D).
.1143D+.72 ; if D43
For TH/TC>l.5 F (D)= .02D+1.02 ; if 3<D<6 (3-7B)
.004D+l.12; if D>6
Where D=(VDR+VDC+VDH)/VC
tion (3-2), we get the basic output for a real engine very accu-
rately.
then the new efficiency can be used and the calculations should
be repeated again. Usually, this has not been the case and
then it means what ever mass flows into the three heat-
81
vatives of the two volume variations give the net volume flow
VA (- SIA/(Wt))/2 (3-10)
(3-11)
I-c W) NCIVM
(3-12)
Where
TW~ot) =(V'N/vC+ CastV)/6IN4P)
Equation (3-10) shows the instantaneous volume flow rate
flow rate, Qmax , and the half period, tH.P., have the follow-
Where
V-C*0NW/(AI*9Cso(
Equation (3-12) is true for an ideal case, i.e. the maxi-
mum flow rate and half period can show the occupied volume
., - f cwoJe
xSIMSI9 NT+*)$2 #4( 3-14 )
P= T= TR (T ,Tc)/
(r./Tc)j
Heater:
(3-15)
Cooler:
Regenerator:
(3-17)
Where Fch, Fcc, and Fcr are the correction factors, which
factors can once again be applied here. This means that each
geometry and the phase angle are kept constant and with the
ratio of the mass flow rate determined from the complete
$ has been changed. Fig. (39) shows the results for f c(')
Heater- t
-.24 (TN/Tc.)+.127(TO/c)+. 8455
(T) j (3-21)
t4 -20 6 (T*/Tc.)+ .79445 ; T/tc)41.i
Cooler-
{ 1972 (Tr/Tc)-.20&$(TN/rc,)+{.310
(3-22)
-7 (To /Tc)4.3o To/-Tc <.S
Regenerator-
TH/-rc >~2
(3-23)
TY/TC(j 2
oI4 (TH/rc) +.423MT)/c+. f/;
-..
For the phase angle difference ():
Heater- fh $)=. 4 6 36 6 7 *1042-1.5516*10-2++.798 (3-24)
tions.
.1373*DH+.8864; DH<2
Heater- fh(DH)= (3-27)
ho .0408*DH+l.08 ; DH>2
-.333*DC+l.293 ; DC.6
Cooler- f (DC)= 3-28)
ca -.lll*DC+1.163 ;DC>.6
1.06 TH<l000 R
Regenerator-f (DR) .06 DR<l.17 ,(3-29)
rb 1/30*DR+l.072 ; DR>l.17&TH>l000R
Where:
DH=VDH/VC
DC=VDC/VC
DR=VDR/VC
(3-15) through (3-17), the mass flow rates inside the three
range .8 to 1.2).
86
(3-30)
Where f is the friction factor, A is the flow path area, L is
RH = AFR*L/(AHT)
For Heater:
A~$.DI 4* -~T CpS20.
(3-31)
For Cooler:
(3-32)
For Regenerator:
('3-334)
-p -o-~c (for sphere filling)
3 --- R
C3- 33 P,)
AFR=VDR/LR
es5
~~ *W CI- WS S(ot )+NS2Yos
I14TT R L(DVnRA Csaok * R*TR F+
(3-33)
t (2-11)
Power loss due to cooler AT= s &T * c
We=1/2 (Wout+Qin)
Wc=1/2 (Qin+Wout)
Wout is given by equation (3-2) and Qin is given by (3-9), mH
Where Pr is the Prandlt number and ReC and ReH are Reynold's
For the regenerator Martini [13] has shown that its in-
Ineffectiveness=2/(NTUV+2) (3-34)
S*cr# I TT-1c)
+
IV ToV+ -( 3-35)
ter for any part of the process. In the next section the sim-
the accuracy of the simplified model has been checked over the
CHAPTER IV
tire engine.
Inputs:
TH=hot source temperature
TC=cold source temperature
Pm=mean pressure of the system (given based on the engine
material)
Qin=available rate of input heat
VD=percentage of the dead volum (based on the available
space)
Performance of the optimum designed engine will be cal-
culated by either complete or simplified models. Both of them
require the complete data about the engine geometry. Deter-
mination of this geometry data is the optimization task which
will be shown in this chapter and is summarized as follows.
a)- By using equation (3-9) the required cold volume
(cold cylinder volume) can be calculated.
T-Th*)4 *1 (3-9)
Nu = 5.75 (4-1)
Qin=h*A*(TH-TR) (4-2)
A=NH*7r*DH*LH (4-3)
NH=(I.l*Qin)/[5.75rKg*LH*(TH-TR)J (4-5
Equation (4- 5) has LH, the heater tube length, but this
regenerator.
heat is less than Qin, then the cold space volume (VC) will
Fig. (59) shows a block diagram of how the optimum design model
works.
put V*&Qi.
4 Cfd i(2V()
(4-8)
Where C, 41itI'T/ceSk)s Fc29 [%%v'f 4. SIA(zA-3)* rf(cdt)]
OR2 . 3 *CP
z~.C 4 ,- Ito Tc. cOS.
(4-9)
Where ,
31.7
53.1 ~
~*)a (~L~r4.1 s~Lzd-ff) 4jzLtCf)
* .
TV
F0(h) 4-Ty) 5T(T)
) [ sr F, 1)]
Cp
S2
oc. Pr4s waslev
cP
ft T*/TC
4I
A' a to F)Z. (I .t.r)2/16186%it 41
~ ~~cw*(Trc,) (335)
Where
In * ) (4-11)
Where 2
'm(
C2. ~ ces 342.U62(k- **
Where
Stirling Engine
with Perfect
Effective heat Components
Leakage Due to Net Output Power
Imperfect Components
Power Loss Due
to Engine Compo-
nent Imperfection
In order to see how much power loss the engine has due
to the heat transfer losses and compare them with the other
power losses we can put an irreversible engine in the by
pass flow and convert those heat leakages to mechanical
power (work transfer) losses. This idea is shown in the
following diagram.
mout
Component Imperfection
103
(4-13)
For the first round calculation efficiency ( can be approxi-
mated as:
Ot.017A.0 .4 Y"gt [2 oS
(4-14)
(4-15)
For typical values of porosity (r) and Prandelt number Pr
i.e. 701* and .72, respectively, with the assumption of hydrogen
as the working fluid the values for C C2, and C can be cal-
culated as:
Ct ~47 I'~oq
CS (4-16)
(4-17)
tion (4-14).
L$)=8go -- 20
(4-18)
( .R+fL
L/)=(% R (4-19)
except that for the cooler there are two types of losses:
Power loss due to the pressure drop; and power loss due to
the temperature drop.
- Cfc) a MA ? +x-%)*Jywz4*cf
WT T3 c(Akts (0C044 TC C*s20
WJher'e
C:: ILyg .-. s6gbg/t) -+J SV(2d+4f)]
egg2Z( (4-21)
LetS define:
Pon DC
SV-t 'A
AR~T 0C2
A.jV'C
Then
COS(
I.
32 P., #DC
R itTC
)E,/ - 2
/.) F 6,Tv, V
00 T (4-22)
Where 3.8S
E I= 2 .04)7 ( su+ +(2 ) S3A11de 4P]
T 1P a'':
7T 4"/C.)IS(CVOjI'
I,Sr-+ s (--)+s u.,)}F c
NI = #CTR-Tc)]A4.1
(.L (4-23)
By using a standard tube diameter for cooler and apply equa-
tion (4-23) for the number of tubes, the dead volume of the
cooler would be constant and frozen for a given Qin and a
given temperature-ratio.
Appendix (L) shows how the power loss due to the tempera-
ture drop in cooler can be written in the following non-dimen-
sional form.
VOUT (4-24)
Where 22 s@ *S --
-041 (-Tev)
STo/T7C
Therefore, the total power loss of the cold heat-exchanger
would be:
OTWus.. (4-2522).
)cc
YD (4-26)
108
IT 1 #(4-27)
(4-30)
(4-31)
Where
VCD/r 1/2(VH+VC)=(Cooler Dead Volume + Cooler Ends Dead Volume/[
-3 = =-f2
= 2-7e > =-
TOT41= ).4
t .U 3 g~
4)e~ V.C 4_____A_(IAlone)_ti _
Woro/r~, 1T3T)&F(4)
(4-32)
WJhere 4 E 1 t(/)
.1 /14 19
t TH -T(4-36)
ill
Where
SVHD/ 1/2(VH+VC)J
=Heater Dead Volume + Heater Ends Volumej/[1/2(VH+VC)J
- 7
5= T12 (4-39)
.N=)-t (4-40)
-* 3 (4-41)
B=.321
v=550
6=1.83
=. 27
This difference in L/D calculation is seven units which has
Chapter III for the different power losses and the basic
given in Chapter II. Since the first three losses are the
dominant ones (as indicated in Fig. (64)), then only they will
113
the maximum net output power. This means that the cylinder
two cylinders (Ae, Ac). But these two areas are directly
basic output power would also increase and both of them have
Non-Uniform Temperature:
as:
A=VH/VC
6=TH/TC
(= (VDC+VDH*TC/TH+VDR*TC/TR)/VC
3 (4-45)
the complete model has been used. It has been run for two
follows.
6=1.8113
6
=1.859 =4 x=1.
6=8678
115
6=1.6226
=2.03 Wo A=1.3
a=9140
6=1.434
(=2.42 =: A=1.1
a=9900
2 =- 1-t E1340('-1)//3
4+ -. 02L(-') Id-5 . (4-46)
to it.
(4-47)
lar to the previous cases, the correction terms have been de-
termined as following.
(4-48)
= -e . 54 (4-50)
the three curves in that figure are all for 165 KW total
this figure also shows that helium can be used with small
tolerable efficiencies.
119
in Figs.(45) and (46). From these figures also we can see the
used and also that helium gives lower performance than hydro-
overall performance.
121
Miodel
optimum phase angle by using the complete model. The only teason
why that complete model can be used for this purpose is through
before differentiation.
Wout W
Then ' . ) - '(Teral Lors) (4-51)
-o Tz4/9
CI)( 4 2*TAjf -
(4-52)
' ~ ~ (Y t7K-l
[(( t )esA'--I- 87 [Cosqc? c S
C(Zo C) &A1 4I[18 (K? ti4),
L4(.44)]
(4-53)
Where / ()f/c42 ic*CScf
(4-54)
The equations for each of the K's are given in Appendix
state speed. This comes from the fact that for an engine,
ihe mean effective pressure decreases with increasing the
engine speed.
then:
Torque i:: 6
-- 84A
2
Power = Torque * "J :& r c -w
have effects only on the heat input not the output power.
fied model.
as follows:
'UT = (4- .
AT-Cooler = k1 1 W
4)eg (4-60)
written as:
A/,f
t) (rf (4-63
12 (4-64)
2.Kt .*)A+-. Y 4+04 Ks)' + 4C ~17+-=k' 7( -6)
some trial and error. In order to get the first guess we can
has a long stroke and small bore, then the temperature gra-
*t (4-67)
(4-68)
(4-69)
and
xm= 1 / 2 S
R (4-70)
8s277 et TPC
*Mwt
(4-71)
TT (2-16)
rewritten as:
3 / - 4' +
(4-72)
Where
K2=.28 Kl (4-74)
IK S
('=*e tT
- T-)*(LST 'f (4-75)
There is an example in Appendix(N) to show how equation
equation.
P-4 (4-76)
M -0
.L .' -
Data
ture range.
one KW.
In 1976, Philips Laboratories [151 started a study on
lysis.
131
were taken.
of 500 Psia and input heat rate 6850 Btu/hr or 2000 Watt have
closer to unity.
CHAPTER V
rator -eometry.
136
perature is 90 F (TC=90"F).
temperature.
TH/TC<l.4 .
Table (15) and Fig. (51) show the results of their ana-
engine.
tures.
139
CHAPTER VI
such a design model appears for the first time in the open
literature.
Finally, the results of the optimum model of Chapter IV
tive conclusions:
a)- Fig. (55) shows the variation of a regenerator
less.
c)- Table(12) and Figs. (45) and (46) show a comparison
optimum phase angle and engine speed, and for optimum bore,
Table (18) and Figs. (45) and (46) show the behavior of
engine can work efficiently with one working fluid. But for
terest are:
1- An experiment should be set up to see the behavior
ponent has been the major part of this engine and its elimi-
practically.
two phase flows then that engine can behave more efficiently
recovery application.
151
Dissip.
TH/TC Wout Energy Wcomp. Qin Ef f.
Btu/lb Btu/lb Btu/lb Btu/lb %
p p 'I
TH/TC=].2
TH/TC=1.6
TH/TC=2. 0
TH/TC-oo
TH/TC = 1.2
Wout/
Wout 100% 45.4% 20.3% 8.0% 2.2% 0.%
ideal
TH/TC = 1.6
Wout/
wout 100% 42.5% 17.4% 6.0% 1.2% 0.%
ideal
TH/TC = 2.0
Wout
wout 100% 42.2% 17.3% 6.0% 1.2% 0.%
ideal
158
TH N Pm Q m w Eff.
R RPM psi atgmi Watts %
Hydrogen Helium
Helium Hydrogen
1.1 - - - - --
1.2 40 5% 2% 118 5% 3%
(a)
(b)
WARMWATER INTAKE
ELECTRIC POWER OUTPUT
VENTI LATIO
DUCTS
WARM
WATER INLE
AFT COLD
WATER
TANKS
EXHAUSTS
FORWARD GATE
COLD VALVES
WATER WARM WATER
EXHAUSTS EXHAUSTS
COLD WATER
PLENUM CHAMBER
RANKINE CYCLE
HEAT RECOVERY
Ref. {12]
175
EXHAUST
STACK FLOW
MIZERI
ONO
EC
EVAPORATOR
SUPERHEATER STEAMLOOP
430*
40 PSIG
$2CF*r GENERATOR
PSIG
.2 PSIG
75 PSiG
CONDENSER COOLINO
WATER
FREONLOOP
>1 50
U
THERMIZER
40
CARNOT
30 /
oq /FOIL FINS
20 -
10!
10
O
0 200 400 600 800
TH 13
Fig. (7): Variation of Efficiency With Heater Temperature
G.M. Benson, Thermal Oscillators, Ref.[5]
177
0.6-
0.6
0.5 TT0*C
0.4 0.4
TH=850*C -H- -- 400*C
0.3 0.3
TH=400*C 0.2 TH 250*C
0.2 HELIUM
HELIUM
Tc = 100* C Tc 0.0* C
0.1
* T= P50*C
0.0 - 0.01
20.0
0.1 0.2 0.5 1.0 2.0 5.0 10.0 20.0 0.I 0.2 0.5 1.0 2.0 5.0 10.0
P(kWATT) -+ P (k WATT) ---
0.6 0.6
0.5 1 0.5
-. -. .TH 850*C
0.4 0.4
TH= 850*C
0.3 0.3 ok
TM= 400*C
0-0 - 0.0,-
0.1 0.2 0.5 1.0 2.0 5.0 10.0 20.0 0.1 0.2 0.5 1.0 2.0 5.0 10.0 20.0
P(kWATT) -+w P(k WATT) -*
Indicated
- - --Shaft
EPANDERIN.ET !FPFRArtURE,
OF
100 200 - in
20
TEMP.
CONIIENSING 1YL.E FID1
3000K (0 F)
"* - - R-111
3110K (000Fl i
wITH
REGENERATrO
0
322'K(120 F) ) -PM13
,VAPIKAIC TFMPERAmTIIE, F -
0 to !0 to 40 it) 60
I I I I I
10 IN1.ET
TjIM . 311., K (20' IJNSRTEP
.ESMNS0EK MITp.
S * 1IND.1 R COlL)
10 0 0 (S0F)
3110 10000F)
0
32? K tI20J'Fl
I
2)0 210 20
(VAIOkAIlKm 'II (OIUTO'Jk
TEMPURLA~IUKf. Cull.)
>1
S50
40
STIRLING
RANKINE
30
Ideal
BRAYTON
- .00e
20
9/f00' STIRLING .NGINE WITHOUT REG.
9O
10
0
1.2 1.4 1.6 1.8 2.0
TH/TC
.25
E-15 =4
a
TH /TC=2
10
0.01
0.0 0.5 1.0 1.5 2.0
DEAD VOLUME RATIO VD/VE
TH/TC=4
.1 6
.12
TH/TC=2
*08
-04 ( I
/
I
I
/
0.0
-50 1-0 1-5 2-0
Swept-Volume Ratio VH/VC
4.20 TH/TC=4
41
0
91 .16
.12
TH/TC= 2
00 wm00
P
woo~~
~ 0 d
0 wo-0
s
0000 -s-
.08
.04
60 70 60 90 100 110 120
Phase Angle #
TH 4 1 --- ,MTF- - P
SE . R - w i 1
S T ./ Q
TCT
R R
0
BOND-GRAPH R
T
SE
80
60
40
20
20 40 60 e0 100
Regenerator Effectiveness 1%]
*8
*o
/ 0
/ (D
CN/
/- 0
o 0 0 0 0
0 QD(0
TepjI'4noM/qnoM
186
Hot Cold
Cylinder Cylinder
QIN
WOUT
IDEAL ENGINE
QOUT
TH P F
SE -- 41. TF TF - .F
TC T
RL
SE
TA-AMBIENT TEMP.
TH
F
SFF- pTF k~
v0
P.
N MT FF- T
S.
1
a)- Simple
~1
SFL- F 1 F 1 MTF
v Q
R T
100
b) - Detaile d x 1P
SHUTTLE LOSS
1 0 R0
C R
QIN
T- S
DOMAI N
Entropy
F-V Generation --
Domain *in Working
OUT Mechanica P-v DOMAIN Fluid
< I ~Dissp.
Heat leakage
Dry Friction
QOUT
co
Fig. (20): Energy Flow of a Real Stirling Engine with Different Loss-Domain
0 5 S S S S
0
0 00 0
0 0 0 0 a
Calculations
Completion
for Geometry
of
1-4
S 0 0 0* 0 0 0 0
U)
10
16
14
12 K= (1.7704*T+521) /17300
10
6
500 700 900 1100 1300 1500
Temperature [R)
Fig. (23) : Thermal Conductivity of Helium and Assigned Correlation
Ref. [13]
21 80
'78
19
17 .76
'74
15
13 .72
11 .70
9
500 700 900 1100 1300 1500
Temperator [*R
Fig. (24) : Thermal Conductivity and Prandtl Number for Hydrogen
and Assigned Correlations,
Ref. {13]
193
6-T=2l2 R
EXPERIMENTj
-TH=1752 R
4Ji
520
COMPLETE MODEL
0
480
0
440
400
360
320-
COMPLETE MODEL
280 0 0
I I I
240
1000 1200 1400 1600 1800
Speed RPM
14001-
12001- I/
7
0-
1000k
is0
7
800 '-Si-z
6001 - I --- -1 O
0-EXPERIMENT
150
130
COMPLETE MODEL
110
90
70
1000 1500 2000 2500 3000
Speed (RPM]
O-EXPERIMENT
Z70
CARNOT
LU
-60
50
COMPLETE MODEL
40
0
30
20 F-
10
1000 1500 2000 2500 3000
SPEED CR PM)
O -EXPERIMENT
150 F-
0
130
0
110 COMPLETE MODEL
90
70
1000 1500 2000 2500 3000
Speed [RPM]
O -EXPERIMENT
80
70
CARNOT
60 -
0 0 0 0
401-
301-
20 -
10
1000 1500 2000 2500 3000
Speed CRPMI
LeRC MODEL
FOR HELIUM I-Pm=600
O-Pm=400 Psia
rP
110 >1
50 .
4-4
90
40
70 COMPLETE
MODEL
30
0
0
50
20
30
1000 1500 2000 2500 3000
Speed (RPM1
.0 70 F
60 0
60 .9-4
44I
1
COMPLETE MODEL 50 r"i
0
50 0
40 40
30
30
20
10 20
1000 1500 2000 2500 3000
Speed (RPM]
0 60
60
50
S4o
30
20
10
-lt
os2 os O.4 0-5 0-6 0-7
1-TC/TH
Fig. (34): Variation of Error between Schmidt and Complete Model
Outputs with Phase Angle
Simplified Model, Equation (3-6)
r-i
15
s4
10
-5
-10 '-
220 240 260 280
qr=360+cr
203
1.5
1.3
1.1
.90
.70
.50
1.0 1.5 2.0 2.5 3.0
Dead Volume Ratio
Fig. (35): Dead Volume Correction Factor for Schmidt
Equationp e M
Simplified modell5quation
, qato (3-7A)
(-A
204
F(D)
1.18
1.14
1.10
1.06
1.02
.9e
.94
.90
2 3 4 5 6 7 8 9 10
Dead Volume Ratio
Cold C3
.I
rr, Time
Half Period
f(T)
1.20
1.121
(T)
.8L
to 2.0 3.0 4.0
1.0
.90
.80
.70
-
-o0
.60
7
.50
.401L
70 80 90 100 110 120 130
yf [DEGI
Fig. (39): Phase Angle Correction Factors for Mass Flow Rates
in Heater, Cooler, and Regenerator
Simplified Model, Equations (3-24) to (3-26)
208
fh(DH)
1.25
1.20
1.15
1.10
1.05
1.0
.95-
o..9 -
0.6 1.2 1.8 2.4 3. 3.6 4.2 4.8
DH
Heater Dead Volume Ratio
Fig. (40): Dead Volume Correction Factor for Mass Flow Rate in
Heater
Simplified Model, Equation (3-27)
209
fc (DC)
1.16
1.14-
1.12 -
1.10 -
1.08
1.06 -
1s04 0
1.02 -
0 .2 .4 6 .8 1. 1.2 14
DC
COOLER DEAD VOLUME RATIO
Fig. (41): Dead Volume Correction Factor for Mass Flow Rate in
Cooler Simplified Model, Equation (3-28)
210
fr (DR
1.18-
1.15-
1.12
1.09
1.06
1.03
10
1 2 3 4
OR
REGENERATOR DEAD VOLUME RATIO
Fig. (42): Dead Volume Correction Factor for Mass Flow Rate in
Regenerator Simplified Model , Equation (3-29)
211
1600 24
20
1200
800
-18
400
-14
I I I
II'
I'
I I
TEMPERATURE RATIO(TH/TC)
60
650
15 00 7-
i 1000
40| -
2.660
1000 1600
400 RjPM
20 4- -4
225hp/ctrwDERWI
0 i I I I I
D 20 40 60 so
HP/LITER
800
ISOTHERMA
CARNOT Ho.
600
OX H%
0
zHe
400-
ADIABATIC
0
1.0 1.2 1.4 1.6 1.8 2.0 2.2
CARNOT
dP
>4
U
z
H 40
ISOTHERA Hi
30 -
20-
ADIABATIC
Temperature-Ratio (TH/TC)
TN
PMI
I.AI
Regenerator
Switch
Contrller - - - -1
r
riP
40-
z
CARNOT
I DIC ATED
P34
30-
200 O SHAFT
COMPUTED RESULTS 0 F
OPTIMUM MODEL
10-
0
1.0 1.2 1.4 1.6 1.8 2.0 2-2
TEMPERATURE-RATIO [%]
&.140 -
>I
z INDICATED
CARNOT
r4
30 -
20 -
/A
10 COMPUTED RESULTS
OF
OPTIMUM MODEL
0
1.0 1.2 1.4 1.6 1.8 2.0 2.2
TEMPERATURE-RATIO(TH/TC)
O. 'hermo-Electron Results
50
40
STIRLING ENGINE
301
RANKINE ENGINE
STIRLING ENGINE
WITHOUT REGENERATOR
0- BARBER-NICHOL'S RESULTS
50-
CARNOT
d 40-
z
W STIRLING ENGINE
H
r4
30 -
RANKINE ENGINE
20
10- --
STIRLING-ENGINEl
WITHOUT REGENERATOR:
0
1.0 1.2 1.4 1.6 1.8 2.0 2.2
TEMPERATURE-RATIO (TH/TC)
P-4
S1.4-
oU21A
U2 W2
W. 1.2
0.6
0.
0.8
0
04
1.6
1.4
0
C?
0
41.0
C- 1
06
0.4
0
02-
0
14
z
S1.2-
0
0
1.0-
rz
018
0.6
0.4
Cl)
0
.-
TEMPERATURE RATIO(TH/TC)
REFERENCES
19- Einar B. Qvale, and J.L. Smith Jr., "A Mathematical Model
for Steady Operation of Stirling Type Engines", Transac-
tions of ASME, Journal of Engineering for Power 1968.
20- Einar B. Qvale, and J.L. Smith Jr., "An Approximate Solu-
tion for the Thermal Performance of the Stirling Engine
Regenerator", AS4E Paper 68-WA-Ener-l, 1968.
APPENDIX(A)
APPENDIX (B)
233
it can have high fuel economy and low pollution and low noise.
Ford has not been able to meet emissions standards with the
vehicle engine. The results of other tests show that the engine
is about 50 lbs over weight objective. The slower start-up and
acceleration is attributed to a higher pressure drop through the
combustion side of the engine than was anticipated. However, for
the program to be successful the big development has to be to
lower cost.
Major technical problems encountered and resolved on Philips-
Ford program are :3]
- Swashplate surface galling.
- Crankcase failure.
234
- Preheater leakage.
- Preheater binding.
APPENDIX (C)
a)- TH/TC=1.l
TC=555 R
TH=610 R
State 1:
T=555 R
P=66.11 Psia
H=33.04 Btu/lb
S=.069 Btu/lb R
State 2:
P=ll2 Psia
H-33.166 Btu/lb
S=.069 Btu/lb R
State 3:
T=610 R
P=112 Psia
H=125 Btu/lb
S=.226 Btu/lb
237
State 4:
P=86 Psia
H=122.5 Btu/lb
S=.226 Btu/lb
T=595 R
State 5:
T=610 R
S=.235 Btu/lb
H=127.3 Btu/lb
State 6:
P=66.11 Psia
H=124.5 Btu/lb
S=.235 Btu/lb
Eff .=Wout/Qin=5.35%
b)- TH/TC=1.2
TC=550 R
TH=660 R
State 1:
T=550 R
P=60.99 Psia
H=31.58 Btu/lb
S=.0663 Btu/lb R
V=l/68.ll ft 3/lb
238
State 2:
T=640 R
S=.0663 Btu/lb
P=2.14.3 Psia
H=32 Btu/lb
State 3:
T=660 R
S=.225 Btu/lb R
H=131 Btu/lb
P=214.3 Psia
State 4:
P=115 Psia
H=124 Btu/lb
T=595 R
State 5:
H=137 Btu/lb
T=660 R
S=.24 Btu/lb R
State 6:
P=60.99 Psia
S=.246 Btu/lb
H=131 Btu/lb
Wout=12.58 Btu/lb
Qin=112 Btu/lb
Eff.=11.23%
c)- TH/TC=1.4
TC=515 R
TH=720 R
239
State 1:
T=515 R
P=33.23 Psia
H=21.86 Btu/lb
S=.0477 Btu/lb
V=1/71.24 ft 3 /lb
State 2:
P=413.8 Psia
H=22.85 Btu/lb
S=.0477 Btu/lb
State 3:
P=413.8 Psia
T=720 R
H=133 Btu/lb
S=.22 Btu/lb
State 4:
P=120 Psia
H=122.5 Btu/lb
T=600 R
State 5:
H=151 Btu/lb
T=720 R
S=.264 Btu/lb R
State 6:
P=33.23 Psia
H=135 Btu/lb
S=.264 Btu/lb R
M-
240
Wout=25.51 Btu/lb
Qin=l38.65 Btu/lb
Eff.=18.4%
d)- TH/TC=1.6
TC=450 R
TH=720 R
State 1:
TC=450 R
P=8.37 Psia
H=6.2 Btu/lb
S=.0146 Btu/lb
State 2:
P=413.8 Psia
H=7.18 Btu/lb
S=.-146 Btu/lb
State 3:
P=413.'8 Psia
T=,720 R
H=133 Btu/lb
S=.22 Btu/lb
State 4:
P=60 Psia
T=552 R
H=116 Btu/lb
State 5:
P=60 Psia
H=154 Btu/lb
T=720 R
241
S=.281 Btu/lb R
e)- TH/TC=1.8
TC=400 R
TH=720 R
State 1:
T=400 R
P=2.136 Psia
H=-3.74 Btu/lb
S=-.0093 Btu/lb R
V=1/80.27 ft 3/lb
State 2:
P=413.8 Psia
H=-2.79 Btu/lb
S=-.0093 Btu/lb R
State 3:
P=413.8 Psia
T=720 R
H=133 Btu/lb
S=.22 Btu/lb
State 4:
P=30 Psia
H=109.77 Btu/lb
T=510 R
State 5:
P=30 Psia
T=720 R
H=155 Btu/lb
S=.295 Btu/lb
242
State 6:
P=2.136 Psia
H=122 Btu/lb
S=.295 Btu/lb
Wout=55.28 Btu/lb
Qin=181.02 Btu/lb
Eff.=30.54%
f)- TH/TC=2.0
TC=380 R
TH=760 R
State 1:
P=1.118 Psia
T=380 R
H=-7.44 Btu/lb
S=-.019 Btu/lb
V=1/81.74 ft 3 /lb
State 2:
P=413.8 Psia
H=-6.5 Btu/lb
S=-.019 Btu/lb
State 3:
T=760 R
P=413.8 Psia
H=147.5 Btu/lb
S=.24 Btu/lb R
State 4:
P=21.5 Psia
T=562 R
H=122 Btu/lb
S=.24 Btu/lb R
State 5:
T=760 R
P=21.5 Psia
H-=165 Btu/lb
S=.318 Btu/lb
244
APPENDIX (D)
written as:
Mh=P*Vh/(R*TH) (D-4)
T2 2.R*TC0 A(To-tT)
(D-7)
nitions.
Tcm=() o+2C/ * S4)/[TC/ + (VC/V#-*C (D-8)
P=Const./(1+A*COS(wt-e)) (D-10)
245
Pmin=Const./(l+A) (D-ll)
Pmax=Const./(l-A) (D-12)
P=Pmaxm(1-A)/(l+A*COS(t-)) (D-14)
r
P 1~/020) .P *d (wt- 6) (D-15)
9*Erg (D-17)
(D-18)
result is:
Wve =lTt* ivi
4,%fAIc SIAC#)/(')+,'A (D-20)
(D-21)
\WO We -W/C PWIf
* NIt UPIT7+
(D-22)
246
Instantaneous Masses:
Hot space:
R1-TH R* ~T cos(ot-4)1
R+
(D-23)
Cold space:
1/2 * _ E+ CosA-9
MC
Ra *C R *-TC jt A C0sc t- El (D-24)
Dead space:
AD
R*LTw+TC) R* (-rct T) +A4cos t )
(D- 25)
Then the total mass in the system is the summation of the above
masses.
Wa4 r I
2(D-27)
VC/a & y TC/T# & 7
- - =
W (D- 29)
'bX *
247
(D-30)
4. a -P /(%+ T)
-r. VO
4-A Xar-+t%) (D-31)
V
X T'=
S I _-4A1_KAI-tTr) -
-tN
(D-32)
Since there is no X term in the numerator it means there is no
would decrease.
_-
A[ * +41(,t ] Te +1r+41Tr/(Te)
(D-35)
248
Assuming COS$%0
(2+.)[
tv-+ IT) [r434-2Tr) +IoXTe +3 +rT(8X+27)I = 6-w )*
IV +Tte4 41 oe6,+T")0
(D-36)
or if A assumed to be constant and first use Vr=1 to calculate A,
ratic equation:
11 I g
____,F I r (AZ+ A? -1+ 4x 2-r 2A IF,' 42T _ +1 A2V', I -1)_
_,eX 14r,,
I-..Tv
2
~ 2*#-" X r
S+'are
A2 A6 %-+ Tr )f SA2X il,1 *
T
-
I+ T,
(D-37)
Finally, to find the optimum phase angle difference $, the
derivative of output is taken respect to $:
V TY -V 10 -bt.* W
(D-38)
*
%TWOc
W+ Y + 2- .T . CS I
(D- 3 9)
""01/ Te AT.
roee 11 W
'+Y;' coy+ V.-rT1r. Su. NI
OIL jii[14 (1--f T+-4Xir.
249
-bq
- =-I
TA +
g[Tjit,-4x-rAtTy]
(D-40)
Or
CoCO-VV f[ {1;tv;2e ,41- F +__ _
I]T,
(D- 41)
This quadratic equation would give the optimum design $
APPENDIX (E)
held, then:
WSXe Jp ch'i (E-l)
=' Y"
(E- 2)
Where
TR=TH-TC/ln TH/TC
r~i (E-3)
WC-0V4
Aye(E-4) S
"T C 4)
(+ ) T (E-5)
(E-6)
IL Oz Yr I
U ~ ~ T. ~ ~ T i 44 T
(E- 7)
251
(E-9)
VMl2 /y-* =
(E-10)
Or,
-Vb 1
'rT n.,4
(E-12)
Z =X
Vi
*I
xr -+4 (E-13)
ra~c
IYT
I... X
44
14 A- (E-14)
____ I V
L. *I L (E-15)
M'iA 1TArc
If X=O then,
WO-7 ='e' A+T
X=P l/P2
Then,
WOU t+~
(TU/TC1)e (E-16)
TK T/Tc -f
e
&-TrqTc & T(/Tc
252
APPENDIX (F)
14
a$/Rl+ae/R2=Pm/t
For cylinder:
Rl=o , R2=R
Therefore,
ae=Pm*R/t , a#=Pm*R/2t
t= (Pm*R) /5000
253
APPENDIX (G)
tribution as:
k Pow e(G-1)
Where $ and $ are phase differences between temperature
displacement.
can assume:
S(G-2)
T 0 !').sl)3t (G-3)
This temperature difference can be used to calculated
a= R Ai 4T (G-4)
254
A =-T * D A,
X = 3-D -Xin ( 14 SLVWea)
(G-5)
0 3( e
=--=Nw P (G-6)
Therefore,
(G-8)
255
I I
capacitances I I
resistances
-a
would be:
. .... c A
4 Cz'
1 L 1 T T
The simplest model would be a two-lump one:
V_// 4r- a 2
Sf C
I*=
=.4
Q-(IT.
LiGLI)
256
Ct 4; rs) (G-9)
(G-10)
Ct I+%5)
Where t =RC ' t =RC, Ct=C +C ' SEd/dt
If we assume a one-lump system then its temperature would be
Q=C t*S*Tm
or
- (G-ll)
(G-12)
Calculation of t & t
................. T&
CwtOt Pr?"it
257
t 1 =RC 2 (G-13)
t =RC 1 C 2 / (C1 +C 2)
(G-14)
amplitude is used.
aT (G-15)
&( (G-16)
2& -T (G-17)
Annular
Internal Cylinder
258
=L +_-D-__ T2 I.
Ig (G-19)
4 ~ 2
or D (G-20)
=> 2
C. = +2
(G-21)
Therefore:
259
'D G8 4 (G-22)
'(G-23)
i (G-26)
; = -98 -<6
(G-27)
Example:
TH=860 R
k =.118 Btu/hr ft R
APPENDIX(H)
CaT = -- pT Z
JA a (H-2)
CV aT = --
P t CpT am m<] (H-3)
RT* k 9 T*(H4
262
and (H-5)
Where K is the specific heat ratio C /C and R is the gas
constant.
Equations (H-4) and (H-5) are valid for both the hot
cylinder, which will be denoted by the subscript h, or the
cold cylinder, which will be denoted by the subscript c.
Since it has been assumed that the heat-exchange components
introduce no pressure drop, then the pressure P will be uni-
form throughout the system at all times.
The mass in the heat-exchange components or dead space
is proportional to the pressure in the system. Therefore, we
get the following.
44 AR v T-) (H-6)
The quantity Vd represents the ratio of the mass contained
in the dead space to the mass contained in one half the volume
displaced by the warm piston at the same pressure and at tem-
perature TH. This mass ratio Vd may be called the reduced dead
volume, since it represents the effect of introducing dead space
in the system. The amount of working gas which must be moved in
an out of the dead space without actually moving through all
the heat-exchange components is proportional to Vd for a given
pressure ratio. This means that additional work must be trans-
ferred in and out of the system in order to pressurize the
dead space without actually increasing the net work.
If (H-6) is differentiated, then the change of mass in
the dead space is given by:
263
4g
YA= VA_ * V~4 I PP/CR *Ti H)7
MYA= Vc +C
t +M Mf Coert (H-8)
or
dMC
mMR40d(H-9)
Because working gas may be accumulated in the dead space,
combinations are:
C":> (H-10a)
dmr <0
(H-10b)
I nA.
>0(H-10c)
(H-10d)
'TC RTc R T4 RT Tl
(H-11)
If the variables
(H-12)
264
P,cT k 0 TC (H-15)
If the dimensionless mass variable
MA' (H-16)
is defined, then (H-15) may be written as
(H-19)
is defined, then (H-18) may be written as
dM . V+ -L A-
.A [jMA>] (H-20)
265
f c (H-21)
and
jd qIg ?
(H-22)
(H-23)
..t dAer + d + .p o0
t A:-- 7 A4 A<v]
(H-24)
rWhl tA4s
angs a
P
While the mass changes are given by
c LIC) [d.4o]
(H- 25)
and
[ j A, Zo , dA4# XD ]
L4 r; (H-27)
and
4 +~k--1
(H- 28)
266
Where
TC (H-30)
267
T -t~: (H-31
and
NTDfA~~
268
APPENDIX (I)
0.1% and 0.5% change. This point has usually been reached
and the pressures, works and masses are normalized with the
cylinder model.
X/L=X
DMRE=
2.T JL Mt (6
M__
273
XII= IlX
X12= I2X
X13= I x/I2X
APPENDIX (J)
C- Lt t (j-3)
-
-c QhT(N&
C .dL ("..- - ChST ) C549t' S14/wt
or
'VNT :\IC'C ' t if CDt + IVW t (j-5)
If we define a to be:
N t - Ces P( )IN
j--6C.
Q .e= , ( IWA - COS) +
2 Co O(L( j-6 )
275
(j -7)
2 CS a
The area under Qnet curve in half of the cycle shows an effec-
(j-8)
Where Tp is the period and TH. is half of the period.
Therefore, if we 2C'o
use a correction factor for effective
exchangers.
For Heater: /i I e)
For Cooler:
P=P ~-/'7i
T= TC Pm Ac.
For Regenerator: P
r PM1
Where TR=(TH-TC)/ln(TH/TC)
277
APPENDIX (K)
#g# A 0 n (k-1)
flow acceleration core friction
and f is the friction factor. For the full cycle the acce-
out of hot and cold spaces are the same. Therefore, the
2M P (k-2)
ev L~-
fCp So(16 N (k-3)
Based on the discussion in Chapter II, the power loss
VL b cLv
AP
V-c -=rc:+Sx )
278
~'*
is 4=
2
- ~L-
.t
ts) i1 ci(Sa _______
j Cw= st-)
0-4.. Cs& CLt'.e)
(k-5)
which
i.(
s a)Cssft-l)ldw C# - as (.wemI
Ther ef or e,
(x- C)
The only unknown in this equation is f, i.e. friction factor,
For Heater:
FR44)(k-7)
RT4 (k-8)
T 4
279
If ReH<2000
(k-9)
If ReH>2000
.067
A
tD
2
= (k-10)
(Re*)
InH is given by equation (3-15).
For Cooler:
4w~
Ke[)c /CAc (k-12)
If ReC<2000
fc=16/Rec (k-13)
If ReC>2000
fc=.0457/(Rec). 2 (k-14)
For Regenerator:
AP LA - (k-15)
fR=53.7/(ReR) 93 (k-16)
If ReR<60
regenerator.
(L/D)R=l/4*AHT/AFR (k-19)
T =_Iq V-TDR
(k-20)
Pb Y=bRA/LR
Therefore,
or
ND -- x-
AF (k-21)
This means:
(k-22)
as: 3
281
AT =Aim
A# rz 12 V Q
~. I
or LD) I2 L?(k-23)
This means:
3 -a- (k-24)
A Sample of SimpifdMoeCacltn
is:
TH=860 R
TC=530 R
VC=2.476*2in 3
VH=2.845*2in 3
Working fluid=helium
Phase angle=108
Heater:
LH=9 in
NH=50
DH=.07 in
VHD=l.3 in3
282
Cooler:
LC=6.8 in
NC=77
DC=.05 in
VCD=l. 22 in 3
Regenerator:
LR=.946 in
DRR=1.5 in
aT=.69
Wire diameter=.0016 in
VRD=l.038 in3
A /,P /Ty+
2
/+T--)(-vL (E+r ,)
Vp4M)ece+ +%
T(R Te -rc-
=:P90= 172'
283
F (T)=. 822
F ($)=. 9644
D= (VDH+VDR+VDC)/VC=1. 5085
F (D)=.8924
Therefore,
Losses Calculations:
Equation (3-12):
Therefore,
f =1.005
cc
fcr=.9324
284
JT3
-6v.6 %) 1-iaVU
( K *n 1.
Wloss=6 Watts
For Cooler:
pc=.0475
Rec 4t 1 135
IT(0c) Wc-/c4
-- S = *x0cY1 = ./
Wloss=13.75 Watts
c.,(#e. ap
For Regenerator:
ReR=38
fR=l.464
285
Wloss=148 Watts
Heater-Equation (2-10):
Mw.-R-To -I
We=l/2(Wout+Qin)=1534 Watts
Wc=l/2(Qin-Wout)=856 Watts
2
NTUH=4(LH/DH) (.023*ReH '*Pr-. 6 )=2.35
Wloss=.6 Watts
Cooler-Equation(2-11)
vie- R-Tc -
NTUC=2.15
Wloss=.6 Watts
Equation (3-35):
#TO rt 2
AIU.=
-74(3 Pr;'3 LAUc.) -7
NTUV=l27.3
QR=66. Watts
Qs=86 Watts
Heat leakage by equation (2-19)
QL=29 Watts
286
8
QFR= 6 Watts
Pumping loss in negligible because temperatures are not very
high.
[equation (2-31)].
QA=3 Watts
Final Results:
Wout=270.8 Watts
Complete model:
Qin=2257 Watts
Effic.ienc =270.8/2459=11%
Efficipncy =262/2257=11.6%
287
APPENDIX (L)
OPTIMIZATION METHOD:
(L-1)
The net output power of a Stirling engine is given by equation
(3-2):
Vv =TF*l*
Pn A x A *-STAU - TC/-*; ()h q)*. D
Nav \+{1
(L-2)
Therefore the normalized power loss would be:
(L-3)
Where:
(E) /
288
wJ6a (I $e
s (L-5)
=C TZW - (L-8)
R T ~*- 2CPk
289
Also
2 2 2 2
Therefore,
Lets define:
CR
ZLCp
b'=A Tg/Tc
cs = C; 1017 2(-0
(L-10)
Where
Foto) YZ
(16-Tc/r) E(T)
(L-11)
290
Or
(L-12)
Where:
0~7 '14
acvZ
[11 4 S/2 9-() +/I q t-~
2 T*t
PvrC)*4J 4V (*C
(L-14)
20 A* (rs
V-C1'NTOV
(L-15)
Where
C2.-CL
=. lj F *
coset "
3'if2.V52 (k-)
291
Where (
=C2 )
Heat transfer loss due to the regenerator axial conduction
is given by equation (2-15):
AR =}r (Ra)
Then the normalized loss would be:
V'VT T- t* PMWO
ICC
Or
devy (L-18)
Where
S*(L-19)
Ct D( - T b)a (L-20)
Carnot Efficiency
.- /0=.
T-.4 1 14.1-
-/D7 Cons .07
cIT (L-22)
293
Where r.)7(y
Lets define: P 3
(L-24)
( 'c*,4RC P-*TC
IAJ (L-25)
Or cr C
0 f c (L-27)
Where:
73 IS
8
(L-28)
(-TC/ITc) FTLT) V)
(L-29)
=V* IT*0 iL
=k'9gS.7 6 /D
AT =.-TR-TC (cooler)
A'T'=*T14 -'R (heater)
For Cooler:
1 c = 5'.76' *A/C *Tri LC x-(TR-TC)
or = aVCr
It/C
LC/(5.75'1T
f Z *LC * (R -Tc)
(L~-'-./-30)
295
wc- , (L-32)
fc * R*Tc_
T k
Where,
NTOC 4 A-023 ac
Since we can be approximated by,
wc s; Pm 6t - %
/--t) wco c t-c * (T H/Tr-.)
Then the normalized power loss would be:
s. A/j'-)(A Y~i)
-A Ftt) * (T .I L P0 (
Viv T 9Vll PFa0) *(T)*
f TC/TVlf*k 4j C .10) -o23 W2 .,v (
(L-33)
Or L-4 1 #
(L-34)
296
Where:
22 g
__ P ~'____
404's (CT/c) IAM?
=TR /TC
Therefore, the total power loss due to the cooler would be:
-r,-'i 2 (LI It--
(L-35)
In equation (L-35) Fc is independent of (L/D) because of
cooler frozen volume for calculation of NC. Therefore, Fc has
no effect in the following differentiation.
I ~V/Oe
D (L-37)
Where:
S2*
Yq
')/ ,V. R, D
VC(L-38)
if-Z4, T * 0Vr)(L-39)
297
= -26 , -16 , -6 , -3 , +4
f( ) = -19.2 * +6
or w
ss149g ($ i-(/g1L-1
Where:
'= TI)/TC
298
Equation (2-16) expresses the shuttle loss as:
--
%AI . (
4 LTI-TC)AX (BET)
Or
Where
P* (
W*t) 2 / +/(f(G)
as: 2K?
&-'~~
R + J4 PI,) (L-43)
(i +1
44 P:S)
Tqvr
(L-44)
Wei- 04 A. Toll
WIns=
V/ou ,. _ vV
T
-.5
9 1.xs 00 X113r AD50_
-+ __ _/
(L-45)
Or i+ g
-5
+2 ~- M.R~ = 0
(L-46)
299
(L-47)
Figures (L-1) through (L-6) show how the optimum geomet-
ries of different components vary with the non-dimensional
numbers, such as Mach number and. Reynolds number.
300
150
130
ClQ-
I110(
90
140
120.C
100
100
140
120
100
C .
Jo
0'
1.4 i
PIN
020
0oo
.80
APPENDIX (M)
-O-- (M-1)
to $.
AP-loss in Heater:
AP-loss in Cooler:
K4 V.6
@( 1e YLI (M-5)
Y Cu(M-6)
1155
The only other losses which have effects on the net out-
put power and include the phase angle are temperature drops
K1 =TA-
4L' Tor)/Loo*))-t ) (M-8)
I T&V.--
t,1 2._V; T 05) (r - 2 2M-9)
28 z bpa Is
.0_ Vca- (I*v'p)
+L.
T2 4-(8 3
4 8 8 8 (M-10)
-93 7 9 A7?
K7 = (M-14)
SIinc -+ CbS
Since
-4/vAod = V7 + IA
Cps
and
T44{.9) = ~cAq
~T+ C
Then :
CaS X =
lV * cs (
(M-15 )
and 0
V' CPCP
VLT" (M-16)
1T p=t WC (M-17)
is given:
WpIV'PLA~tA
(M-19)
y & .'
-o4 42.a+ FI -iA9/(t (M-20)
309
If we define M as:
"9I (M-21)
(M-22)
Therefore,
STQ cv(
(K, t c S~ IAVcf)- (c.-scp Vj [KX3+ 4)1A', 7
1" J
- +'+ ( s10- kt 7Cs')=o
(M-24)
TVV -+2Wr4)Yo ;
lo -= 2- P%-,V2rT, COS I
Therefore,
+ COS- - +3 kM-27)
Example:
We want to calculate the optimum $ for the previous
~ A.
K2=.3770
311
Y
= VASS
*3 1V %5W(2U-)+, 59039
K 4= 11.4
k 5=146
temperature.
Since. $q = 17 Uept
p
C 5-- 15.z-37
1,637-121I
A4 =i
KI = . 11-B 02
e3+4M + 1.07 ksit "3 174
re
opt =110 then we can see how well the above equation predict
the optimum $.
Heater:
(M-28)
Where - 0 7 .457'T . 2
-*6 057 D'A#o
qt 4
I. A A t I~
t- q VI
(M-29)
313
A~:'v
A, c ~DC~
- DeC 11 T A, 4-
T' S's4Ar
e
(M-30)
Regenerator:
(M-31)
CL ol
(M-32)
Therefore, pressure drop power losses can be written as:
Heater:
ess I (M-33)
Cooler:
l,-os3 -34 LA
(M-34)
Regenerator:
LS3 1 9 Y7
(M-35)
For Power loss due to the temperature drops in heat and
dependent are:
SS wU)
k& 'k/91001 (M-39)
Therefore, the net output can be written as
KW 9
2.avt uJ.... 2.+? 4 -4(jK7k.+
g).67 c4
(M-41)
315
AI100 RPA
WZ , 606P =2.'1rli1
J4Y
xt1o
VI
2. -0.27
W pt=ll20 RPM
Wopt=1080 RPM
I
316
APPENDIX (N)
follows.
C = , At 6T (N-1)
AT is given by (2-22).
2
A=Tr/4*B
( p(N-2) P)
This equation can be integrated over a cycle to find the
integration follows.
(N-3)
Y3 (N-4)
-2
(N-5)
fs4
-. ~ -Ot i
(N-6)
0=.0727 /0 /1
(N-7)
33
/7
-4 -53
Since the shuttle loss is very significant, then we have
- .. (2-16)
5 -T Tc
318
-I *
3s 0 (~q -TC 4~I %
a 7) it TTs t - (N-9)
L
Where
I4::9
sOP. ' h-I it -P-1- i CtS ( 5 (N-12)
Example:
N=l000 RPM
P=500 Psia
p=.0 6 8 lbm/hr ft
38925 d, 8t 4 . ott
8345
Pl/Pm. 1
Kg=.l18 Btu/hr ft
Q =17
.i3
V=2. 845 in
1=.005 in
319
L=1.5 in (assumed)
From equation (4-73) through (4-75) the K's have been cal-
culated:
K =107.3
K 2=30
K 3=162.3
APPENDIX (P)
. . 0
TI, I *
TC
'4T
P U 30 (P-3)
Fcr=1*
____SY f C~ LT,-TC~ 2, k r-
3 C&SdL (P-4)
or
6 2K* 2 PMrrc Cp %_____ T
=- Rol 0 "-0-7-) CS -R-Tc (P-5)
tor material.
322
APPENDIX (Q)
20000 MM=1
SHW=SHW+1.
NAH=1
PUSH=l.
OMG=3.1415*N/30
TR=(TW-TC)/ALOG(TW/TC)
IF(M.EQ.2.)GOTO 1635
MIUC=(196.14+.464*(TC/1.8-293)-.093*PM*.00689)/1000000*241.9076
MIUW=(196.14+.464*(TW/1.8-293)-.093*PM*.00689)/1000000*241.9076
GOTO 1636
1635 MIUC=(88.73+.2*(TC/1.8-293)+.118*PM*.00689)/1000000*241.9076
MIUW=(88.73+.2*(TW/1.8-293)+.118*PM*.00689)/1O00000*241.9076
1636 NMM=1
AKC=(1. 7704*TC+521) /17300
AKH= (1.7704*TW+521) /17300
IF(M.EQ.2. )AKC=(.0194444*TC+8. )/173
IF(M.EQ.2.)AKH=(.0194444*TW+8.)/173
1740 NW=WH*3.563/(AKH*LW*(TW-TR) )*1.1*ZAP
NC=AWC*3.563/(AKC*LC*(TR-TC) )*1.1*ZAP
IF(NW.LT.10) NW=10.
IF(NC.LT.10) NC=10.
30000 NC=NC*FACTC
NW=NW*FACT H
C OPTIMUM COOLER
1728 BC=((R*TC)**.5)*RC*RC/(OMG*VAC)*34.047
GC=PM* (RC**3) *83462.4/ (OMG*VAC*MIUC)
A=1.535*(NC**.8)
DLC=TW/TC
IF(SIG.EQ.1. )VCD=VCDO+.5*(VDEER+3.1415*DRR*DRR*LR/4)
IF(MM.GE.2)GOTO 1700
ZC=VCD/ (2*VAC)
GOTO 1702
1700 ZC=VCD/ (VAC+VAW)
1702 LC=RC*(DLC*A*(GC**.2)*(BC**.6)+192-59.2*DLC+.0275*GC
1-688.3*BC-19.2*ZC)
AFRC=NC*3.1415*RC*RC/4
C
C OPTIMUM HEATER
C
BH=( (R*TW)**.5) *RW*RW/(OMG*VAC)*34.047
GH=PM* (RW**3)*83462.4/ (OMG*VAC*MIUW)
DLH=TW/TC
B=1.535*(NW**.8)
IF(SIG. EQ.1. )VWD=VWDO+.5* (VDEER+3.1415*DRR*DRR*LR/4)
IF(MM.GE.2)GOTO 1703
ZH=VWD/(2*VAC)
GOTO 1704
1703 ZH=VWD/ (VAC+VAW)
1704 LW=RW*(DLH*B*(GH**.2)*(BH**.6)+244- 7 5.7*DLH+GH/100-356*BH-
116*ZH)
AFRW=NW*3 .1415*RW*RW/4
NMM=NMM+1
IF(NMM.LE.3)GOTO 1740
3241
C
OPTIMUM REGENERATOR
IF(SIG.EQ.1.)GOTO 1749
VEF=AFRW*LW+AFRC*LC+VDEER+VWD+VCD+3.1415/4*DRR*DRR*LR
FOR METAL C=.11 BTU/LB/R , RHO=488 LB/FT3
CAL=(VAW/VAC+COS(AFI) )/SQRT (1+(VAW/VAC)**2+2*VAW/VAC*COS(AFI))
DRR=(288/3*PM*VAC/(488*TR)*SHR/(SHR-1)*SIN(AFI)/(1-SIG)/.11
1*(TW-TC)/(TR-TC)/(CAL*778*LR))**.5
AL=ATAN( (VAW/VAC+COS(AFI) )/SIN(AFI))
A=47.6*SIN(AFI)/(CAL)/(SIN(AFI)+1/6*
1SIN(2*AL-AFI)+. 5*SIN(2*AL+AFI))
DEL=TH/TC
DR=SIG*D0/( 1-SIG)
Z=.11*778*(SHR -1)/(R*SHR)*OMG*DR*DR*488/
1( (MIUC+MIuW)/2)*3600/144
GA=(DR**3)*PM*83462.L4/(MG*VAC*(MIUC+MIUW))
BET=((R*TR)**.5)*DR*DR*((1-SIG)**2)*68.o94
1/(SIG*0MG*VAC)
IF(EFF.EQ.0 0 ) EFF=0.9-1/DEL
AREG=(.285*A*((ALOG(DEL))**1.66) *EFF*GA*
1GA/(BET*BET*(Z**. 66)) )**2. 66
ZET=(VDR+VDEER )/(VAW+VAC)
FZ=80*ZET-20.
ARN=AREG+FZ
LR=ARN*DO*SIG/(1-SIG)
1749 VDDR=SIG*(3.1415*DRR*DRR*LR/4)
AFRR=VDDR/LR
C
C OPTIMUM CYLINDERS
C
BCY=PM*OMG*VAC*VAC/(R*TC*MIUC*AC*AC)*3600*.047548
DLCY-TW/TC
ZCY=(AFRC*LC+VCD+TC/TW*(AFRW*LW+VWD)+VDDR*TC/TR)/VAC
FZCY=.357*ZCY-.78
FDL=.548*DLCY-.54
FBY=.525-8.8*BCY/100000
AAA=.11*(ZCY+1)*DLCY/(1+3400*(DLCY-1)/BCY-.02*(DLCY-1)*( (BCY
1/DLCY)**.5))
VAW=VAC* (AAA**. 5+FZCY+FBY+FDL)
WRITE(6,1743 )VAW,BCY, DLCY, ZCY
IF(SIG.EQ.1)GOTO 1730
IF(NCC.NE.1)GOTO 1730
DELA=(TC/TW)**2+2*TC*VAC/(TW*VAW)*COS(AFI)+(VAC/VAW)**2
DELB=TC/TW+VAC/VAW+2*TC/VAW*((AFRW*LW+VWD)/TW+(VDDR+VDEER)/TR
1+(AFRC*LC+VCD) /TC)
325
DEL=(DELA**. 5)/DELB
TET=ATAN(VAC*SIN(AFI) / (VAW*(TC/TW+VAC/VAW*COS(AFI))))
VAC=WOUT*(1+(1-DEL*DEL)**.5)*1.409*3/(OMG*PM*DEL*SIN(TET)*
1(1-TC/TW) ) /1. 5*AZ
IF(VAC.LT.O.9) VAC=0.9
NCC=NCC+1
GOTO 30000
1730 VAW=VAW
IF(ABS(VAW/AVM-1) .LT.0.05)GOTO 1725
AVM=VAW
MM=MM+1
GOTO 1728
1743 FORMAT(///////////,1OX, 'REAL VAW=', F1O.5, //, 5X, 'BCY=',F1O.5
1,//,5X,'DLCY=',F10.5,//,5X,'ZCY=',F10.5)
1725 WRITE(6,u1726)VAC,VAWLW, LR,LC,NWNC,DRR
1726 FORMAT(1H1,///,15X, 'OPTIMUM DESIGN RESULTS' ,///,5X, 'VAC=',
1F10. 5,/ /, X, 'VAW=', F10. 5, //, 5X, 'LH=' , F10. 5,//,5X, 'LR=', F10. 5
1,//,5X,'LC=',F10.5,//,5X,'NW=',F10.5,//,5X,'NC=',F10.5, //,
15X, 'DRR=',FlO.5)
IF(SIG.EQ.1)GOTO 5000
IF(SIG. EQ. 1)GOTO 5000
IF(CYL.EQ.1. )GOTO 4000
DR=SIG*DO/ (1-SIG)
GOTO 5000
4000 DR=2*SIG*D0/(3*(1-SIG))
5000 VDC= (LC*AFRC*TW)/ (VAW*TC)
VDEC=VCD*TW/ (VAW*TC)
VDEW=VWD/VAW
VDW=AFRW*LW/VAW
VDER=VDEER*2*TW/ (VAW*(TW+TC))
VDR=VDDR*2*TW/ (VAW*(TW+TC))
RVT=VAC*TW/(VAW*TC)
VD=VDC+VDEC+VDER+VDR+VDW+VDEW
NNN=NDIV
XND=NDIV
NDIV1 = NDIV + 1
DALF = 6.2831853/XND
NT = NDIV/4+2
NE = NT - 1
ALF = 4.7123889
NF = NDIV/4
CALL VOLC(DALFNFC,CIDCDCI,ZZCNDIV,SIFI, COFI, SALF, CALF)
100 FI = DALF*NFI
SFI = SIN(FI)
CFI = COS(FI)
IND(1,1) = 1
OCALL VOLW(W,WIDWDWI, CFISFI, ZZWNDIV,SIFICOFISALF,CALF,
DALF)
IND(1,2) =3
IND(2,1) =4
IND(2,2) =2
NN = 1
XMC = 0.
p = 1.
326 (a)
XMW = 1.-CFI
XMWS = 0.
PS = 1.
NO = 4
WW = 0.
WC = 0.
NITE = 1
NST = 1
NFIN = NFI
NLOP = 0
404 DO 102 I=NST,NFIN
VW = W(I)
vc = C(I)
VWI = WI(I)
VCI = CI(I)
DVW = DW(I)
DVC = DC(I)
DVWI = DWI(I)
DVCI = DCI(I)
GO TO (201,202,203,204),NO
601 NLOP = NLOP+1
IF(NLOP-6) 201,201,605
201 DP = -SHR*P*(RVT*DVCI+DVWI)/(RVT*VCI+VWI+SHR*VD)*DALF
S = P+DP/2.
DP = -SHR*S* (RVT*DVC+DVW)/ (RVT*VC+VW+SHR*VD)*DALF
DMW = S*DVW*DALF+VW*DP/SHR
DMC = -(DMW+VD*DP)/RVT
IF(DMW)302, 301, 301
301 K = 1
GO TO 303
302 K = 2
303 IF(DMC) 304,305,305
305 L = 1
GO TO 306
304 L = 2
306 NO = IND(KL)
GO TO (400,602,603,604),NO
602 NLOP = NLOP+1
IF(NLOP-6) 202,202,605
202 IF(XMC) 803,801,801
803 XMC = 0.0
801 IF(XMW) 805, 802, 802
805 XMW = 0.0
8020DP -SHR* (XMC*RVT*DVCI/VCI+XMW*DVWI/VWI) /
1 (XMC*RVT/P+XMW/P+SHR*VD)*DALF
DMC = XMC*(DVCI*DALF/VCI+DP/SHR/P)
DMW = -RVT*DMC-VD*DP
S = P+DP/2.
SMC = XMC+DMC/2.
SMW = XMW+DMW/2.
ODP = -SHR*(SMC*RVT*DVC/VC+SMW*DVW/VW)/
1 (SMC*RVT/S+SMW/S+SHR*VD)*DALF
326 (b)
DMC = SMC*(DVC*DALF/VC+DP/SHR/S)
DMW = -RVT*DMC-VD*DP
IF (DMW)3 12, 312, 307
312 K = 2
GO TO 308
307 K = 1
308 IF(DMC) 309,309,310
309 L = 2
GO TO 311
310 L = 1
311 NO = IND(KL)
GO TO (601,400,603,604),NO
603 NLOP = NLOP+1
IF(NLOP-6) 203,203,605
203 IF(XMC) 704,703,703
704 XMC = 0.
7030DP = -SHR*(P*DVWI+XMC*RVT*DVCI/VCI)/(VWI+XMC*RVT
1 /P+SHR*VD)*DALF
DMC = XMC*(DVCI*DALF/VCI+DP/SHR/P)
DMW = -RVT*DMC-VD*DP
S = P+DP/2.
SMC = XMC+DMC/2.
SMW = XMW+DMW/2.
ODP = -SHR*(S*DVW+SMC*RVT*DVC/VC)/ (VW+SMC*RVT
1 /S+SHR*VD) *DALF
DMC = SMC*(DVC*DALF/VC+DP/SHR/S)
DMW = -RVT*DMC-VD*DP
IF(DMW) 313,314,314
314 K = 1
GO TO 315
313 K = 2
315 IF (DMC) 316,316,317
316 L = 2
GO TO 318
317 L = 1
318 NO = IND(K,L)
GO TO (601,602,400,604),NO
604 NLOP = NLOP+1
IF(NLOP-6) 204,204,605
204 IF(XMW) 705,702,702
705 XMW = 0.
7020DP = -SHR*(P*RVT*DVCI+XMW*DVWI/VWI)/(RVT*VCI
1 +XMW/P+SHR*VD) *DALF
DMW = XMW*(DVWI*DALF/VWI+DP/SHR/P)
DMC = -(DMW+VD*DP)/RVT
S = P+DP/2.
SMC = XMC+DMC/2.
SMW = XMW+DMW/2.
ODP = -SHR*(S*RVT*DVC+SMW*DVW/VW)/(RVT*VC
1 +SMW/S+SHR*VD) *DALF
DMW = SMW* (DVW*DALF/VW+DP/SHR/S)
DMC= -(DMW+VD*DP)/RVT
IF(DMW) 319, 319, 320
319 K = 2
GO TO 321
320 K = 1
327
CMMAX = XLARGE(XMCX,NDIV)
WMMAX = XLARGE(XMWX,NDIV)
CMMAX = CMMAX/PMAX
WMMAX = WMMAX/PMAX
WRITE(6,19) WMMAXCMMAX
ARG = 2.*RP/(RP-1.)*WC/3.1416
IF(1.-ARG**2) 1607,1608,1608
1608 FIPV=ASIN(ARG)
WRITE(6,18) FIPV
XNDS = NDS
X = 0.
DX = 1./XNDS
WRITE(6,21)
NIN = NDS + 1
COR = PMAX**(XNHT-2.)*DALF**(XNHT-1.)
DO 854 I=lNIN
CALL PDINT(X, XDMWXDMC,RVT, DCNDIV,DMREPR,XINTDPRXI1, XI2, XNHT)
XINT = XINT/DALF/PMAX
AXINT(I)=XINT
DMRE = DMRE/PMAX/6.2832
ADMRE(I)=DMRE
XIl = XI1*COR/(1.5708*DMRE)**(1.-XNHT)
X12 = XI2*COR/(1.5708*DMRE)**(2.-XNHT)
AXI2(I)=XI2
X13 = XI1/XI2
AXI3(I)=XI3
AX(I)=X
WRITE(6,22) XXINT,DMRE,XI1,XI2,XI3
X = X+DX
854 CONTINUE
1607 IF(NWR) 1509,606,1509
1509 DO 509 I=1,NDIV
PR(I) = PR(I)/PMAX
XMCX(I) = XMCX(I)/PMAX
509 XMWX(I) = XMWX(I)/PMAX
WI(NDIV1) = WI(1)
CI(NDIV1) = CI(I)
DO 1001 I=1,NDIV
IF(XMCX(I)) 1003,1003,1002
1002 TEC(I) = PR(I)*CI(I+1)/XMCX(I)
GO TO 1006
1003 TEC(I) = 0.
1006 IF(XMWX(I)) 1004,1004,1005
1005 TEW(I) = PR(I)*WI(I+1)/XMWX(I)
GO TO 1001
1004 TEW(I) = 0.
1001 CONTINUE
TEW(NDIV1) = TEW(1)
TEC(NDIV1) = TEC(1)
PR(NDIV1) = PR(1)
XMCX(NDIV1) = XMCX(1)
XMWX(NDIV1) = XMWX(1)
329
TWDM = 0.
TCDM = 0.
DO 573 I=1,NDIV
DMW = XMWX(I+1)-XMWX(I)
IF(DMW) 574,575,575
574 TMPW = (TEW(I)+TEW(I+1))/2.
TWDM = TWDM+( TMPW-1. ) *DMW
575 DMC = XMCX(I+1)-XMCX(I)
IF(DMC) 576,573,573
576 TMPC = (TEC(I)+TEC(I+1))/2.
TCDM = TCDM+(TMPC-1.)*DMC
573 CONTINUE
TWDM = TWDM*SHR/(SHR-1.)
TCDM = TCDM*SHR/(SHR-1.)
DO 1021 I=1,NDIV
1021 XMT(I) = XMCX(I)*RVT+XMWX(I)+PR(I)*VD
WRITE (6, 51) TWDM , TCDM
WRITE(6,14) PR(NDIV),(PR(I),I=1,NDIV)
WRITE(6, 15) XMCX(NDIV), (XMCX(I) ,I=1,NDIV)
WRITE(6,16) XMWX(NDIV), (XMWX(I) ,I=1,NDIV)
WRITE(6, 1010) TEC(NDIV),(TEC(I),I=1,NDIV)
WRITE(6,1011) TEW(NDIV), (TEW(I),I=1,NDIV)
WRITE(6,1022) (XMT(I),I=1,NDIV)
606 NDIV=O
IF(NDIV) 511, 511, 100
511 GO TO 10000
110FORMAT (22HlSPECIFIC HEAT RATIO =F8.3,10X,18H (VC/VW)*(TW/TC) =F8.
13//17X,5H FI =F8.3,22X,5H VD =F8.3//11X,11H DIVISION =I8//)
120FORMAT (10X,12H COLD WORK =F8.3//10X,12H WARM WORK =F8.3//5X,
117H PRESSURE RATIO =F8.3///)
13 FORMAT (10X,10,11H ITERATIONS )
14 FORMAT (24H1ARRAYS START AT MC = 0. ///12H PRESSURE /(10F10.4))
15 FORMAT (12HlCOLD MASS /(1OF10.4))
16 FORMAT (12HlWARM MASS /(10F10.4))
17 FORMAT (SHOLOOP 3110)
18 FORMAT (20H P-V ANGLE IN RAD = F10.4//)
19 FORMAT(16H MAX WARM MASS = F10.4, 16H MAX COLD MASS = F10.4//)
21 FORMAT (/// 23H PRESSURE DROP INTEGRAL)
220FORMAT (6H X/L = F6.2,11H INTEGRAL = E12.4,7H DMRE = E12.4,
1 6H XII = E12.4, 6H X12 = E12.4,6H X13 = E12.4/)
45 FORMAT (2110,10X,6E10.3)
510FORMAT (16H INTEGRAL (H*DM) /9H WARM END F10.4,10X,9H COLD END
1 F10.4,1H1)
1010 FORMAT(12HlCOLD TEMP /(10F10.4))
1011 FORMAT(12H1WARM TEMP /(10F10.4))
1022 FORMAT (liHlTOTAL MASS/(10F10.4))
10000 HR=SHR
IF(SHR.LT.1.4.AND.M.EQ.2. )HR=1.4
IF(SHR.LT.1.4.AND.M.EQ.4.)HR=1.66
NDIV=NNN
PMIN=2*PM/ (1+RP)
PMAX=RP*PMIN
IF(SIG.EQ.1)DR=DRR
330
FCR=10.**(ARl)
GO TO 716
715 AR2=. 714-.365*ALOG10(RECR)
FCR=10. **(AR2)
716 QCRP=(OMG**3) *PMAX* (VAC**3 ) *(LV/DR)*XINTCR*FCR
1/(TC*AFRR*AFRR*R*699212.9672*.73756)
C REGENERATOR HOT END
REWR=RECR* (MIUC/MIUW) * (DMREWR/ DMRECR)
IF(REWR.LE.60. )GOTO 717
IF(REWR.GT.60.AND.REWR.LE.1000)GOTO 718
AR4=. 015-.125*ALOG10(REWR)
FWR=10.**(AR4)
GO TO 719
717 AR5=1.73-.93*ALOG10(REWR)
FWR=10.**(AR5)
GO TO 719
718 AR6=.714-.365*ALOG10(REWR)
FWR=1 0.* * (AR6)
719 QWRP=QCRP*(XINTWR/XINTCR)*(FWR/FCR)*(TW/TC)
C MIDDLE OF REGENERATOR
REAM=RECR*(MIUC/(MIUW+MIUC)*2)*(DMREMR/DMRECR)
IF(REAM.LE.60.)GOTO 720
IF(REAM.GT.60.AND.REAM.LE.1000)GOTO 721
AR7=.015-.125*ALOG1O(REAM)
FMR=l0. ** (AR7)
GO TO 722
720 AR8=1. 73 -. 93*ALOG10(REAM)
FMR=10.**(AR8)
GO TO 722
721 AR9=.714-.365*ALOG10(REAM)
FMR=10. **(AR9)
722 QMRP=QCRP*(XINTMR/XINTCR)*(FMR/FCR)*(TW+TC)
1/(2*TC)
WRITE(6, 10006)REAVC, REAVW, RECR,REWR, REAM,FC, FWFCRFWR,FMR
10006 FORMAT(//,5X,5F13.5,5F7.3)
C TOTAL ENERGY DROP IN REGENERATOR
ORP=(QCRP+QWRP+4*Q4RP) /6
C HEAT TRANSFER LOSSES
C COOLER
VDR=VDDR+VDEER
ANUC=.023*(REAVC**.8)*(PRC**.4)
ANTUC=4*LC*ANUC/(RC*REAVC*PRC)
AMAC=1.5708*DMRECL
AC=EXP(2*ANTUC) -1
COOLOSS=WDC*WC*(HR-1) /(HR*AC*2*AMAC* .73756)
C HEATER
ANUW=.023*(REAVW**.8)*(PRW**.4)
ANTUW=4*LW*ANUW/ (RW*REAVW*PRW)
AMAW=1. 5708*DMREHT
AW=EXP(2*ANTUW)-l.
HEATLOSS=WDW*WW* (HR-1)/(HR*AW*2*AMAW*.73756)
C REGENERATOR
C COLD END
IF(SIG.EQ. 1. )GOTO 24000
33~3
REC=l. 5708*RECR
STC=1. 1*SIG/ (REC** .405)/ (PRC** .667)
ANTUCR=4*LR/DR*STC
C HOT END
REH-l. 5708*REWR
STW=1.l*SIG/((REH**.405)*(PRW**.667))
ANTUWR=4*LR/DR*STW
AMR=(DMRECR+DMREMR/2)/ (DREWR+DMREMR/2)
ANPH=.63662*VDR*(HR-1) /(VAC*HR* (DMRECR+DMREMR
1/2) * (TW/TC-1) )
ALANDR=(1+ANPH* (XI3CR+XI3WR*AMR)/2)/ (ANTUCR/
1XI2CR+ANTUWR*AMR/XI2WR)
IF (SIG-l)2 50000 24000 25000
24000 ALANDR=l.
25000 RHEATLOSS=ALANDR*HR/(HR-1) * (TW-TC)/TC*PMAX*VAC
1*OMG*DMRECR/ (24*M*.73756)
DLPC=8*FC*LC* ((PMAX*VAC*N*DMRECL/(RC*RC*NC ))**2) / (28980*RC*PM
1*R*TC)
DLPH=8*FW*LW*((PMAX*VAC*N*DMREHT/(RW*RW*NW))**2)/(28980*RW*PM*
1R*TC) *TW/TC
DLPR=8*FMR*LR*((PMAX*VAC*N*DMREMR/(SIG*DRR*DRR) )**2 )/(28980*DR*
lPM*R*TC) *TR/TC
WRITE(6, 26000)ALANDRDLPC iDLPHDLPR
26000 FORMAT(///, lOX, 'LANDA=', FlO.6,///
1,5X, 'PRESSURE DROP IN COOLER EPSI]=' ,F8.3,//,5X, 'PRESSURE DROP
1 IN HEATER EPSI]=',F8.3,//,5X, 'PRESSURE DROP IN REG. [PSI]=,
1F8.2)
C SHUTTLE HEAT TRANSFER
C ASSUMING CONSTANT:DENSITY=488 LB/FT3 , C=.1l BTU/LB.R
C FOR AKGAKM THE UNITS ARE WATT/CM.K
AKG=(1.7704*TR+521) /1000000
IF(M.EQ.2. )AKG=( .0194444*TR+8)/l0000
AKM=.01111*( (TR/1.8)** .45425)
ALP=AKM/ .928664
1624 BS=BS
B=(8*VAW*BS/3. 1415)**. 3333
S=(8*VAW/(3 .1415*BS*BS))** .3333
AL1=AKM/AKG* ( (OMG*GL*GL/ (2*ALP) )**. 5) *5
BET1=(2*AL1*ALl-AL1)/(2*AL1*AL1-1.)
SHUTLOSS=AKG*S*(TW-TC)*B/(2*GL)*BET1*.554142
C REGENERATOR AXIAL CONDUCTION
AMM=(1+AKM/AKG) / (1-AKM/AKG)
VDR=VDDR
AKMG=AKG* (AMM+S IG-1. )/(AMM+1. -S IG)
AXILOSS=AKMG*DRR*DRR*(TW-TC)/LR*1.108284075
C DRY FRICTION LOSS
FRLOSS=N*( .002*N+1)*(VAW+VAC) /265.5*SW1
C SPRING GAS LOSSES
AE=(8.*VAW*VAW*3.1415/BS)**.3333
PA=(PM**.5)*(RP-1)/(2*(RP**.5))
AKL=(1.7704*TC+521) /1000000
334
IF(M.EQ.2. )AKL=(.0194444*TC+8)/10000
AC=(8. *VAC*VAC*3.1415/BS) ** .3333
ALANDA=VAW/VAC
VETA=(AFRC*LC+VCD+TC/TW* (AFRW*LW+VWD)+VDDR*TC/TR)/VAC
AKW=(1.7704*TW+521) /1000000
IF(M.EQ. 2. )AKW=(.0194444*TW+8)/10000
FE=3.84*PM*N*( (2*VAW)**.6667)/(100000*TW*MIUW)+2.22
FC=2.1*PM*N*( (2*VAC)**.6667)/(100000*TC*MIUC)+2.95
SGL=.35381*((OMG*(((SHR-1)/SHR)**3))**.5)*PA/(1+ALANDA
1*TC/TW+VETA)*(TC/TW*((AKW*TW)**.5)*FE*AE+( (AKL*TC)**.5)
1*FC*AC)*((COS(FIPV)-SIN(FIPV))*(ALANDA+COS(AFI))-(SIN(AFI))
1* (SIN(FIPV)+COS(FIPV)))
TIC=PM*((VAW/3.1415)**.3333)/5000
ARI=B/2
ARO=TIC+ARI
AB=12/ (1089.54*S*AKG)+ALOG(ARO/ARI) *12/(3.1415*S*10)+
1144/ (21.99*S*ARO)*2
QLEAK=(TW-TC)/AB*778/ (3600*.73756)
GK=AKW/ .0173
QP=2*OMG/(3*3.1415)*(300/778*OMG*B/(2*GK))**.6*((12*(PMAX-PMIN)
1*GL/TR)**1.6)*GL*(TW-TC)*S*HR/(144*.73756*(HR-1))
C FINAL RESULTS
C TCP=TOTAL COOL CYLINDER SHAFT POWER [WATTS]
C THP=TOTAL HOT CYLINDER SHAFT POWER [WATTS]
C THC=TOTAL HEAT FLOW IN COOLER [WATTS]
C THH=TOTAL HEAT FLOW IN HEATER [WATTS]
C SGL=SPRING GAS LOSS [WATTS]
TCP=WDC/.73756-QCP-QWP-QRP-COOLOSS-HEATLOSS
THC=TCP-RHEATLOSS-SHUTLOSS-AXILOSS-SGL
THH=WDW/.73756+RHEATLOSS+AXILOSS
C BORE-STROKE OPTIMIZATION
AK=(1.7704*TW+521) /17300
AKK1=.026*((PM*MG*3600/(R*TW*MIUW))**.8)*(PRW**.4)*(SHR-1)
1SHR*PA/PM*AK*TW*COS(FIPV)*( (4*VAW/3.1415)**.9)
IF(M.EQ. 2. )AK=( .0194444*TW+8)/173
AKK2= *28*AKK1
AKK3=3.1415/8*AK/(S*GL)*(TW-TC)*(BET1)*4*VAW/3.1415*778/12/3600
RLG=.4*AKK1*(BS**2. 6)+3*AKK3
RIG=2. 6*AKK2* (BS**3.6)
IF(ABS(RIG-RGL)-10 )1610,1611,1612
1610 BS=BS+.05
IF(BS.GT.2.5)GOTO 1611
GOTO 1624
1611 WRITE (6,1626)BS
GOTO 1625
1612 BS=BS-.05
IF(BS.LT.O.8)GOTO 1611
GOTO 1624
* 1626 FORMAT(////,1OX, 'OPTIMUM BORE-STROKE RATIO=',F5.2)
1625 IF(NAH.GT.3)GOTO 28000
NAH=NAH+1
WH=ABS(THH)
AWC=ABS (THC)
335
WRITE (6,1741)WH,AWC
1741 FORMAT(/////, 5X, 'WH= , F10. 5,//, 5X, 'WC=',F10. 5)
SHR=SH
GOTO 1740
28000 IF(SIG.NE.1.)GOTO 27000
NUM=NUM+1
IF(NUM.NE.2)GOTO 27000
FACTC=FACTC*THC/ (THC+RHEATLOSS)
FACTH=FACTH*THH/(THH-RHEATLOSS)
SHR=SH
GO TO 30000
27000 THP=WDW/.73756-QLEAK-QP
OUTPUT=THP+THC+RHEATLOSS+AXILOSS-FRLOSS
DLV=2*VAW+VAC* (1+COS(AFI))
BMEP=0UTPUT*.73756*12/ (DLV*OMG)
PANOR=BMEP*OMG*DLV/(12*.73756)
THPN=THP/PANOR
TCPN=TCP/PANOR
WRITE(6, 1726)VACVAWLW,LRLC,NWNCDRR
THCN=THC/PANOR
WCC=WDC/ .73756
WRITE(6, 10011)PMAX, PMIN, TCTW,N,SIZE, WCC
10011 FORMAT(1H1,//,10X, 'MAX. PRESSURE EPSIA]=',F7.2,//,1OX
1, 'MIN. PRESSURE EPSIA]=', F7.2, //,10X, 'COLD TEMP. [R]='
lF6.1,//,10X, 'HOT TEMP.[R]=',F6.1,//,10X,'SPEEDERPMJ=', F6.1
1, //,10X, 'SIZE=', F6.3, //, 10X, 'COMPRESSION WORK=', F14.5)
IF(SIG.NE.1. )GOTO 7013
WRITE(6,10010)FACTC, FACTH
10010 FORMAT(///,20X, 'ALL OF THE FOLLOWING QUANTITIES ARE IN WATTS',
1/////,5X, 'FACTOR FOR CHANGING COOLER=' ,F7.3,//,5X, 'FACTOR FOR
1 CHANGING HEATER=',F7.3)
7013 THHN=THH/PANOR
QCPN=QCP/PANOR
QWPN-QWP/PANOR
QRPN=QRP/PANOR
QM=QP/PANOR
COOLOSN=COOLOSS/PANOR
HEATLOSN=HEATLOSS/PANOR
RHEATLOSN=RHEATLOSS/PANOR
SHTOSN=SHUTLOSS/PANOR
AXILOSN=AXILOSS/PANOR
FRLOSN=FRLOSS/PANOR
SGLN=SGL/PANOR
OUTPUTN=OUTPUT/PANOR
TOR=OUTPUT/OMG*. 73756
TORN=TOR/( (PMAX-PMIN) *2*VAW) *12
BMEPN=BMEP/ (PMAX-PM IN)
EFFN=100*OUTPUT/(THH* (1.-TC/TW))
TOTLOSS=QCP+QWP+QRP+COOLOSS+HEATLOSS+RHEATLOSS+SHUTLOSS+AXILOSS
1+FRLOSS+SGL+QLEAK+QP
IF(SIG.EQ.1.) TOTLOSS=TOTLOSS-RHEATLOSS
336
TOTLOSN=TOTLOSS/PANOR
EFF=OUTPUT/THH*1 00
QL*OLEAK
QN=QLEAK/PANOR
WRITE(6, 90)QCP, QWP, QRP, COOLOSS, HEATLOSS, RHEATLOSS, SHUTLOSSIQL
1,QOP
1, SGL, AXILOSS, FRLOSS, TCP, THP,THC,THH,OUTPUT, EFF,TOR, BMEP,
1TOTLOSS
90 FORMAT(//, 25X, 'PRESSURE LOSSES:',//, 5X, 'COOLER=', F10.5, /
1,5X,'HEATER=',F10.5,/,5X,'REGENERATOR=',F10.5,//,25X,'
1HEAT TRANSFER LOSSES:' ,/,5X, 'COOLER=' ,F1O.5,/,5X, 'HEATER
1=',F1O.5,/,5X,'REGENERATOR=',F13.5,/,5X, 'SHUTTLE LOSS=',
1F10.5, /,5X, 'HEAT LEAKAGE=', F10.5, /5X, 'PUMPING LOSS=',
IF10.5,/,5X,'SPRING GAS LOSS=',F1O.5,/, 5X,'AXIAL LOSS=',F1O.5,
1/, 5X, 'COULUMB FRICTION LOSS=',F1O. 5,//, 5X, 'TOTAL COOL CYLINDER
1 POWER=',F12.5,/,5X,'TOTAL HOT CYLINDER POWER='',F12.5,/,5X,
1'TOTAL HEAT FLOW IN COOLER=',F14.5,/,5X, 'TOTAL HEAT FLOW IN
1 HEATER=',F12.5,//,10X,'NET OUTPUT POWER=',F12.5,/,10X,'EFFIC
1IENCY=',F8.4, '%',///,10X, 'TORQUE [LB-FT]=', F12.5, //,10X, 'BMEP
1 PSIAt=',F12.5,//,10X,'TOTAL LOSS=', F12.5)
IF(ABS(OIN-THH) .LE.50. )GOTO 1738
IF (QIN-THH )1627 ,1738, 1628
1627 VAC=VAC-SAR
IF(VAC.LE.O.7.OR.SHW.GT.(VAX/SAR) )GOTO 1738
SHR=SH
GOTO 20000
1628 VAC=VAC+SAR
IF(VAC.GT.VAX.OR.SHW.GT. (VAX/SAR))GOTO 1738
SHR=SH
GOTO 20000
1738 WRITE(6, 10013)PANORN, WDWN, WDCN
10013 FORMAT(////,20X, 'THE NONDIMENSIONAL OUTPUTS IN PERCENTAGE' ,////
1, 10X, 'NORMALIZE QUANTITY: ',//, 10X, 'PM*VC*N [WATTS)=' ,F11.3,
1////,1OX,'WARM WORK=',F10.5,///,10X,'COLD WORK' ,FlO.5)
WRITE (6, 90)QCPN,QWPN,QRPN, COOLOSN,HEATLOSN, RHEATLOSN, SHTOSN,QN,
1M,
1SGLN,AXILOSN, FRLOSSN,TCPN, THPN,THCN,THHN,OUTPUTN,EFFN, TORN,
1BMEPN, TOTLOSN
STOP
END
337
FUNCTION SMALL(X,NDIV)
DIMENSION X(720)
SMALL = X(1)
DO 507 I=2,NDIV
IF(SMALL-X(I)) 507, 507, 508
508 SMALL = X(I)
507 CONTINUE
RETURN
END
NC= 90.00000
DRR= 1.98000 in
FI = 252.000 VD = 4.628
DIVISION = 360
2 ITERATIONS
MAX WARM MASS = 1.6278 MAX COLD MASS 1.3421
SPEED[RPM] =1200.0
PHASE ANGLE=108