AD-825198 Simulation of Turbofan Engine - Part II - User's Manual and Computer Program Listing
AD-825198 Simulation of Turbofan Engine - Part II - User's Manual and Computer Program Listing
AD NUMBER
AD825198
FROM
Distribution authorized to U.S. Gov't.
agencies and their contractors; Critical
Technology; NOV 1967. Other requests shall
be referred to Air Force Aero Propulsion
Laboratory, Wright-Patterson AFB, OH
45433.
AUTHORITY
"AFAPL-TR-67-125
PART 11
1-OVEMBER 1967
-qo
mm
• • mm
• m m mmint•
s m •m • mm
mm m • •mm ,
NOTICE
When Government. drawings, specifications, or other data are used for any purpose
connection with a definAt31y
other than inStates
the United
related Government procurement operation,
Government thereby incurs no responsibility nor any obligation -
whatsoever; and the fact that the Government may have formulated, furnisched, or in
anywa supled th addaigeecificatlons, or other data, is not to bee regarded
by Implication or otherwise as in any manner Licensing the holder or Any other person
or corporatlon, or conveying any rights or permission to manufature, use. or sell any
patented invention that may in any way be related thereto.
Warm
FOREWORD
This report was prepared in the Components Branch (APTC). Turbine Engine Division. Air
Force Aero Propulsion Laboratory. Wright-Patterson Air Force Base, Ohio, under Project
3066, "Gas Turbine Technology," Task 306603, "Advanced Engine Studies." with Charles E.
Bentz as Project Engineer.
This report covers work conducted withdn the Components Branch in the time period between
July 1965 and June 1967 and was submitted by the author 31 August 1967.
ERNESTT
Chief, C. S ON
Turbine Engine Division
Air Force.Aero Prou.lsion Laboratory -
1~
l -+
.4
.. ,+ Iz. •
AFAPL-TR-67-125
Part U jABSTRACT
This report d&scribes a digital computer program titled •TO-E (Simzlaft-sn of Turbofan
Ecgne). SMOTE s a&computer program for balancing-cy-le turbofan engines capable of r'unnlng
both design and off-.eslgn points. The program is ,ritten In rortran IV language, -d was
dcsa.gned for use or- an IBM 7090 Digital Computer, altho,.-h It has also beer, run on an IBM
System 380. Performance maps (Block Data format) of '.- 6 major engine compoaients are
.required.
Information for seWtig up the Block Data and L-mit data is given In the report. A-30
i~luded In the report is a complete program listing w iwecription of e4cb subroutine and
Sa inple data pasck.
(Distribution oi the abstract is unlimited.)
l$
-if
-____ ___ I
i4
A FAPL-TR-67-ý.25
Part 111
TA.BLE OF CONTENTS-
SECTION
PAGE
I IWTRODUCTION i
u- BR.OCK DATA
2
mI INPUT DATA
3
1. Controlled Output 3
2. Namelist Lnput 3
IV SUBROUTINE DESCRIPTIONS
4 V
V PROGRAM LISTING
1:-
V.
• -
ri•_1
IAMB
"APAL-TR-67-125
Part n
* SYMBOLS
STAVTIN~ NUMBERS
2 fan entrance
21 f=n e.,dt/conmpreasor and duct entrance
* 3 com ssr exit/burner entr-=ce
4 bttrm'r exit/hi preasure turbine entrance
5 bi Press. turbine exit/bo press. fturbine etrbance
55 lo: parsý. -rblze i.x1
6 :.fterhzrnoz- enL-ance
7 ~aftartxuner exit
8 mn 'rethroat
23 &ductburner eaft-rance
24 dwct bmrnr exit
25 duct exit (if mized-flaw engine)
23 duct nozzle throat
29 dact nozzle elxit
TRIfMODNA1Uc PROPERTIES
T tatz temponmture (degrees -5)
H h~ltr 0022Ipy
P tOWa Peenre (atmospheres)
S total ent~ropy
T's static temperatur'e (degrees R)
PS' stati pressure (aftnos-pheres)
FAR twl-air rstic
AM Usch number
V 'Mekclty (fbot/secood)
COMPOKENr SYMBOIB
F Ian
C eeiresor
S lxln~er
C0M burner
THP hi jczessurs turbin
TLPI lo resxure turbine
A snerbizrner
-~~~ AFT frbre
DUC duct
OB overboard
T towa
El~NG SYMBOLS
I
A FAPL-TR-67- 125
Part II
SYMBOLS (CONT)
DP pressure drop
DT temperature increase
MISCELLANEOUS
A area
ALTP altitude
AM Mach number of aircraft
BYPASS bypass ratio
CF correction factor
CS ambient speed of sound
CV nozzle velocity coefficient
DEL delta degradIon coefficient
DS design value
DUM dummy (not used)
FG gross thrust
FGM nonleDtUM 1-rust
PIGP pressure thrust
FN net tbrust
FRD ram drag
NOTE: Some symbols may be truncated wben combined with other symbols due to six
character limit Imposed by Fortran Computer Language.
I -
* 4 -
i f
I.-
'it i
%I
A FAPL-TR-67-125
Part II
SECTION i
INTRODUCTION
The purpose of this report is to describe
a digital computer p-ogram for balancing-cycle
turbofan engines. The progra!m, titled
SMOTE (Simulation of Turbotan Engine).
in the Components Branch, Turbine Engine wmas developed
Division, Air Force Aero Propulsion Laboratory.
The program, as written and described
in this report, is interded for use on
Digital Computa'x. but it haE also been an IBM 70o0
modified and tested on an IBM System
360. It is written
in Fortran ITV language.
The pe--fo--Pae of the major engine
components is based on component maps
be converted into Block Data subroutines which must
for use by SY.MOTE. Presently included
are maps for the fan. compressor, combut-or, and both turblnts. .n the p-rog
I1
f
j AFAPL~- I-T7-61-2•
.
PartlP
SECTION I1
To provide the basic component performance requires the maps for the fan, compressor,
(Oombustor tud both turbines to be converted to Block Data subroutines. The Block Data pre-
selted in thia report is very general and does not represent any specific engine. For the follow-
Ing dissosion on setting up Block Data, refer to the program listing, SectiL.. V. The format
for all the maps is very similar, with a maximum of 15 curves and a maximum of 15 points for
each curve. "N" refers to the numberofcurvee and "NP" is an array of the number of points
on euc curve.
For the fpr and compressor mape, the corrected speed (CN) Is Input as a ratio of the actual
corrected speed curve to the design corrected speed. "1he pressure ratio (PrR). corrected air-
flow• (WAC), and efficleac (ETA) are input In groups, with a &-ufor each corrected speed
curve. The corrected speeds and pressure ratios must be loaded low to high. Note that there
=ar Itwo points per card.
For the comb~tor mp-p, curves of pressure ente-ring the combustor (P3) in PSI aire loaded
low to high, as wsll DELT
as delta-T (T4 - T3) points f~r each curve. The efluiencies (ETA) cor-
responding ta each.h are loaded in one array. Note that in this array, each columrn is
obtained from one pressure cur-ve.
* For the burbine maps, turbine flow function (TFF) curves are input low to high, aad the
corrected speed (Ct), delta-H (MII), and efficieicy (ETA) are input !n groups, with a group for
each TFF !in--. The corrected speednmustbe loaded lowtohigh. TFF, CN. and DH are defined as
TFF= PIN
PLN
-
CN = P
DH HIN -HOUT
Formats other than those presented here may be used for the component performance mnaps.
ifowever, some hange in !ogic would then be required.
[.o
It-
A FAPL-TR-67-125
Part 1,
SECTION III
INPUT DATA
The Input data Is dlv4ded into two sections; data cards for the controlled output and data
cards in Namelist format for running each point. For the following'divoussif on setting up
the input data, refer to the listing of sample data immediately following the program listing.
1. CONTROLLED OUTPUT
The variables that are to be output are selected by the first section of data cards. Any vari-
able that is in one of the main commons (DESIGN, FRONT, SIDE, or BACK) may be selected
for output by punching the name of the variable as it appears in the common (with trailing
blanks, if necessary) in Columns 1 through 6. Up to 102 variables (seventeen lines of six viri-"
ables) may be chosen for a particular run. During the output phase, the name of the variable is
printed out, with its value printed immediately below the name.
Another feature of the controlled output is the ability to change the name of a variable to be
output; for example, it may be desired to change a station designation to one more common to
a particular programmer. In this case, the variable name would be punched in Columns I
through 6 as described above, b•t in addition, the desired name would be punched in Columns 13
through 18. Special symbols, such as /, miay be used in the new name. The last card of the
controlled output mi'st be a card with "THEEND" punched in Columns 1 through 6.
2. NAMELIST INPUT
The normal data for running the dczired points follows the conirollee -otiptut data and is in a
Namelist format, where the name of the Namelist is DATAIN. Ustirdly the first set of data it
the design point, as shown in the sample input data. When the design point is'run ,(IDES-1i), all
map scaling or correction factors are printed out,as well as being retained in common. There-
fore, it is possible to run off-design points immediately following the desiga point by making
"-seof the values in common, or to begin running an off-design point immediately by inputting
the scalingor correction factors. The first method is usually easier, but the second method may
be desired if many points are to be run using the same engine parameters with no changes
except for power setting, Mach number, and altitude.
The controls which fix the type of turbofan, the mode of operation, the method of calculating
ram recoveryd etc., are explained in the listingof subroutine INPUT. Subroatine ZERO'deter-
mines what values in common will be zeroed between points. None of the'design values or
correction factors are ever zeroed.
A title card must be input immediately after the first point of the data pack, and ITITLE must
be set equal to 1 in the data for the first point. This 1i' because a title to always printed for
each point and must, therefore, be previously defined. The input format for the title It 12A6,
and the resulting 72 spaces are centered on the page when printed out. The title may be changed
by setting ITITLE = 1 and inserting a new title card after the Namelist data for the point.
When an afterburning cr duct-burning point is run, the exhaust nozzle areas are allowed to
float to obtain optimum expansion. This means that there can be no balancing it the point, and it
is necessary to prebalance the engine cycle in a nonaugmented mode. That Is, an identical
point, except that it is nonaugmnented, must be run before either afterburning or ductburning.
When either IAFTBN or IDBURN is gr6ater than zero, the program will automatically set
INIT = 1 and use the balanced values from the preceding point. The nozzle areas are returned
to their standard design values after completing an augmented point. Some examples of this
type of afterburning are given in the sample data listing.
3
Am ,
iLj A L-TR-67-125
Part ji
I ~SECTION IVI
SULBROUTINE DESCRIPTIONS
The following will be a brief description of vdmt each subroutine calculates or controls.
See legure I for a s'broutine flow chart
SMOTE Dummy main program to initiate the calculations and cause the L-Mut of the
controlled output varir.bleb. Because of the loopingbetweensubroutine., control
Is never trasferred back to this routine.
ENGBAL Main subroutine. Controls all engine balancing loops. checks tolerances and
E Bnumber of looce. and loan- natrixn
MATRIX Solves error matrix.
GUESS Deterno nes initial values of independent variables (PCNF, PCNC, &-d T4) at
. each point It may be desired to change these equations to suit a partlcular
Sengine. The closer the initial values are to the final values, the faster the
program will balance.
COFAN Uses Blck Data to perform fan calculations.
COCOMP Uses Block Data to perform compressor calculations.
SCOCOS Uses Block Data to perform combustor calculations. May use either T4 or
WFB as the main pa-rameter.
COBFPT Ues Block Data to perform high prssure turbine calculations. Calculates
EMIR() and ERR(2).
COLPITB Uses Block Data to perform low pressure turbine calculations. Calculates ERR(3)
and ERR(4).
FRTOSD Dummy routine to tranfer val-ue from common FRONT to common SIDE.
CODUCT Performs duct and duct-burning cmlculations. May use eltber T24 of WFD as the
main paramet-r for duct-burnin& Controls the duct nozzle and cc ulhea EPI=J5}
If in separate-flow mode.
FASTBK Dummy routie to transfer values fram commons FRO,•T and S.DE to aommon
BACK.
- ,
It
AFAPL-TR-67-125
Part U
II Ii
,) rwm
! I
II
FAR T]IMU NOZZLE
HAIM
U
PERIG-ME W.JlOflLrh
lI
S1- 11 1
I
I I
I!
7-T
II
AFAPL-TR-67-125
-- ii--pu P
COMIX Performs gas-mixing calculations if in mixed-flow mode. At design points it
calculates areas from either an input static pressure (P"55) or an input Mach
number (AM55) if PS55 = 0. At off-design points it calculates static pressures
and Mach numbers from the design areas. Calculates ERR(-5).
S COAFEN Performs after-burning calculations. May use either T7 or WFA as ihe main
parameter.
SYG Controls printing from UNiTOB. Throughout the program and particular.-y in
ENGB.A.L. certain messages. variables, and marix values arý. written on UNIT08
as an aid In determining why an error occurred or why a point did not balance.
These values ace nrintedcutlfsubroutineERRORis called and 1DUMP Is great-
e. than zero, or after a good point IDUMP= 2.
TAPES Defines UNIT08, which Is Just a "scratch" diskand does not require a $SETUP
"card. Normal inpat and output are on UNITO5 and UNIT06, respectively.
THCOMP Performs Isentropic calculations for compressors.
SEARCH General tao:e look-up and interpolation routine to obtain data from the Block
Data sutroutines.
MAPBAC U3ed 'when calcula""ons result in values not on the turbin~e maps. Charges the
map value and an Independent variable %PCNF,PCNC, or T4) in an attempt to
x-ectify the situ.ffon.
AFAPL-TR-67-125
PartI :
"I-i]
I:-.l
The Bolo
cknData forpet lisigow
pres subreutrines reuie tornMIT.T rel
I-
ISECTION V
° 1-2
P•ROGRAM LIST'ING i:
The following is a rcmplete listing of all subroutines re~jufred to run SMOTE. Inunedaiely
following the program listing Is a listing of a sample input data pack~. -r
!-
JI
! I
I_ - I
7!
Pift H
Pad iH
,P9
,XPZ5
PXFAR24*XX~PI
P7 ,147
,H9
tDPAFT ,V55
,XS25r
iDUMB
IS7
IS9
YV25
i
t
I[
6PS6 ,V6 ,AM6 ,1T57 fps? ,V7 tAX7 2AM25 t
7TSB ,PS8
,VJD
Iv3 ,AMS
,FGKD ,VJMf
vTS9 SPS9
#FGfN ,FGPD ,FGPM
'V9 tAR11!
EVA ,FRO
9FGi' ,FGP SWFT I,;GT ,FART ,PG ,Fit ,SFC
DIMENSION VAR(6hODEL(6),EfRB(6),DELVAR(6),ENAT(6,6),YMAT(6),
ilAJ4T 6)
DATA AWOS(O/6HENHG9AL/ [
CALL INPUT
IF(INIT.EO.)) GO TO 50
TFFW2=TFHPDS
TFFLP=TFLPOS
NHN4AP=0
IGO=2
VMAT(R )=0. i
AMAT(I )=0. -
DELVAR( 11-0.
I DU 2 L=1#6
3 UMJP6R=LDOPE~t*1
T CAIL C-OFAM
I 1F(LGUPERGT.lTRYS) GO TO 18
1FINO)'AP.GT.O! GO TOI
£ VAR(1)-ZF*100.
I VA512)=PChF
VAR!3J=ZC'*100.
I lF il4ODE.EQU.O.OR.KO0E.EQ.2) VAR(41=PCr4C
j IHPMODE.EO.1) VAR(4)T4/1O.
VAR (5)=TFFHP
~VAR(6 jzTFFLP
00 4 1=1,6
14 lF(A8BS(tE~R(J).GT.TOLALL) 60 'U0 5
1
CALL PERF
WI
*
5
CALL ER'JP.
IFILLtfi.6T.0)
MAPEDG2O
~MAPSET wO
U 6 1=116
GO TO 7
EIMB(l )=ERR( I
7 lFf(.APE1JG.E0.LJ) GO TO 70
IMAP6OG=0
MAPSET2-1
VAR (LO0P)=VAK(LOOPJ4+2.SOEL(LOOPI
GO TO 10
70 UU'MAPSET.EU.01 VA8ILOOP)=VAR(L0QPJ.DEL(LLJOP)
lF(ItAPSET.EO.Il VAk(LWJP)=VAR(LUOP)-DEL(LUOPJ
I4APSET=O
Vu 8 1=1,8
8 E-4AT(ILLEP)=(EkP.B(1J-ERRftfJ)/OE-L(LOjPI
19
I
LOOP=LOOP41
~IFtLOOP.GT.6) GOJ TO 11
j j VAR(LO0P)=VAR(LUOP)-UEL (LWUPJ
It ZF=VAktl)/100.
PCNF=VAR (2)
ZC=VAi(3)/100.
I jF('4EODF.EU.'J.Ok.FPjE.ElU.Z)
FtM.OVE.E0.11 Tv4ýVAR(41910.
TFFHP=VA8 (51
PCNC=VAgt4)
- Tlf:LP=VAR(6)
IF(ZF.LT.D.) ZFO.O5
IFIZC.LT.Gi.) ZC=O.05
GO TO (1,3!,160
* 11 W0 12 1=1.6
00 14 J~lt6
IZER0=0
shG 13 LUOP=1,6
13 JF(EMtATf!,LUJPJ-.Eg.O.) lZEHO=IZk:X!+I
J IlFIIZERJJ.LT.b) GO TO 14
IdRITE16,100P1
LEWlPERJ TRYS+10flJ
GO TO IS5
14 CONTMtUE
00 16 Ltt0Pz1*6
10
U9Tw
AFAPL~-TR-67-125
Paz's U
I ZHW)P=O
Lfli 15 1=1,6
15 F(EUAT(1,L(J0P).kL,~0.) IZI:R0=IZENW+j
IF!IZERiJ.LT.6) GO TO 16
WRITE(6,101)LOUP I
LO(WER=i TRYS+l00
GO TO 18I
16 CONTINUE
CALL 1itTRIX(EMAT,V?-ATtAMAr!
DO 17 LcJOP=i,6
DELVAR (LLX"JP,=0.6-*VMAT (LOOP)
IF(ABS(O)E1VAR{LO0P)).GT.0.05*VAR(LOUaP))
1OELVAR(LOOPh=0.05*VAR(LOOP) DELVARtLOOP)/ABS(DELVARtLEUUP)
DO 19 1=1,6
DELAVE=(ABS(DELVAk(1l))ABS(DFLVAR(2)l:ABS(DELVAR(3))+
I ABS(iELVAR(4)l+ABS(DELVAR(5))+ABS(OELVbRI6Jfl/6.
WRITE(8,104)ERRAVE,DELAVE ~
IF(LOOPER.LT.ITRYS) GO TO 10
CALL ERROR
RETURN
IG10 FORMAT(4HOROW112,16H IS ZERO IN EMAT)
101 FL3RHAT(HOCOLUMN,12,16H JIS ZERO IN EMAT)
102 FORMAT(8148 ERRB*28X23HERROR MATRIX AFTER LOUP,!4?29X4iV.MAT,
I6X6HDELVARIXi4HVARIABLESS$SSS)
103 FiJRMAT(1H0,F8.4,8X6F1O.4,IOXF10..4,Fll.4,4XF11.4,6H$S$$$4i: 1-
104 FORMAT(1HOF8.4,32XL4HAVERAGE VALUES,42XF11.4,6H$Si$SS)3
END
* ~APAP!,TR-A~7-i-A
Part H
ii
5 T(ITiJ,jJ=T(IFROHJ)+RM*PIV(J)
I T0=I M-1
6 IFROIM=1FRDM-1
IF(IFROM.GE.1) GO TO 4
DO 7 J=!Pl,7
7 TtI,J)=PIV(J)
D0 a 1=1,5
J=7-1
K=6-1
00 8 (.=Jw6I
a T(K,71='I(K,. T(K*L)*T(Lt7)
D09 1=1,6
9 V(I)=T(1,7)
RETURN
END
'11
*SEIM
"AFAL-TR-67- 125
Part U7
14V
U AF"APL-T21-67-125
Part 11
-
TZC=T2
P2QOP2
DO 11I=1,63
!12=P20 4
CALL SfYt1 -
RETURN U
Emir
I.[4 . .......
'C-
ASF.. 1 T .j1
rat
A. SIFCCIL PCr91t
M AL
t.S4-t?!Hr-PZFT
21
T21 ,qP21 IN2TCW'SrET1DS.T3CSDCO ,DTCOCFEPA3CI
3TF0SCY*P4 EH4 S4 75FFC #C~PSF '*IC ,T2D
1 2HP5 I
5CUFL VPSFPD ELOSF IVCFC IL1'C VA3TLPCF ,riO
#OHLPC
6CN? *PS PETAC
DUD;WAOOWAC3O ETAJS tDPCOUC4TAOF ,
BA55W A2 *ETYLt*
HCLADTF #9L8 *MG5 viA285 #,A29
WOATA WOIRD1COHT
AJ=778.26
G=32. 174049
ALT-ALTP*2.0855531-4.O7/(12.O855531 Eo07.ALTP)
CALL AT~flS(ALTT12XXJXX2,XX3,P1,CSXX4,iIIR;
IF'lIAXTP.EQ.2! T1-Tl*TZ
1 FtlA,4TPd4E.1) CALL RANCA14,EAR)
CALL PRO 531(FA5R,T1,C3.XX2,XX3,R1,P~ilthl)
SI*PhI 1-R1SALOG(PI)
?#2=Hl+( A14CS j**Z! (2. eA*G)
no 2 1=1710
CALL THiERMO(P2T:H2,T2TS2TAV%,0,0.O,1)
2 PT=P1*EP((Metlk.98t375)t((S2T-Sl;t(1.986375/AW1.ALOG(P2T/Pfl).
* T20SST2
5 ZFXZ$L3S
6 REIVRft
~-77
A FA:'L -TR-67-125
--
1,,I
17
>1:
~Lri -
"AP.-TP.-67-125
I
1962
I C
C
C
UaS. STAROART ATMOSPHERE UP TO 9g KILOMETERTt
CALLING SEOUENCE...
C CALL ATHOS 'ZFT9 TM, SIGMA, RHUt THETAt DELTA, CAt A.t), X.)
C ZFT GE•METRIC ALTITUDE (FEET)
C TM = OLECULAR SCALE TEMPERATURE IDEGREEs%RANKINE)
C S[GA = PATIO OF JENSITY TO THAT AT SEA LEVEL
C (HU = DENSITY (LB-SEC**2-FT*.(-4) ,R SLUGS-FT**3)
C THETA = RATIU OF TEPPERATURE TU THAT AT SEA LEVEL
DELTA = RATIO OF PRESSURE TU THAT AT SEA LEVEL
"•C CA = SPEED OF SOUND (FT/SECI
C AKU = VISCOSITY CUEFFICIENT (LB-SEC/FT,.2
I C
C
C
X1=INORMIAL
= 2 ALTITUDE
z 3 FLUAT•NG
LESS THAN -5000
POIh-T OVERFLOW
METERS UR GREATER THAN 90 KM
C
C ALL DATA A40 Fi.UAMEI4TAL CON4STANTS ARE IN THE 1.ETRIC SYSTEM AS
- : C THESE 0iJANTITIES ARE DEFIED AS EXACT IN THIS SYSTEM.
C
C THE RAI)IUS OU THf EARTH IREFT59i IS THE VALUE ASSOCIATED WITH THE
C 1959 AROC ATMOSPHERE SO THAT PROGRAMS CURRErdLY USING THE LIBRARY
C kUUTINF kILL NOT REUWIRE ALTERATIUN TU USE THIS ROUTIIIE.
DIMENSION -B' 10) tT.B 1IO DELTAB( 1) ),ALlo(101
•+ ~~~~~DA,1A
IMBI ]) 7mbf]), UELTABI I) , L i)]=101
A -. 3.0, 320.65, 1.75363E 00, -6.5,
b 0.0, 288.151, 1.OOOOE 00, -6.5,
C 11.0, 216.65, 2.23361E-01, 0.cl,
U 20.0, 21.65#, 5.403?8tE-02, 1.0,
E 32.0, 228.65, 8._6663E-03, 2o8,
F 47.0, 270.65t 1.09455E-03, 0.01,
G 52.0, 270.65, 5.82289E-04, -2.U,
H 01.07 •
2S2.65, 679718E-04, -4.0,
1 79.U, 16U.65, 1.024) E-059 0.0,
1 -8.74T3, 18b.65, 1.6223 E-069 0.01
DATA REFT5912.0855531E 07/, GZ /9.80665/,
-"A AnZ /Zb.9b44 IT HSTAR /8.31432/,
v :" FTT(IKM/->.uAaEU-O4 I, S 1110.4 I,
. I C AMUZ 11.i.Z24E-05 I, CAL 11116.45/,
p P.H1*1 /0.076474 I, GZEN• /32.1?41/
C CnVERT GEIWETRIC ALTITUDE TO: GEOPUTENTIAL ALTITUDE
HFT = fkEFT59/(KEFT54+ZFT])17FT
C C(iNVERT liFT AND, ZFT To KILtOETFRS
I = FTTOKM*ZFT
m = FTIOKhFT-
K= I1
Tps7 = TMH(2)
IF (H.LT.-5.O.DR.Z.GT.9D.0) GO TO 16
-ou 10 ?=,10
IF (H-aBI4.IM) 11,12,10
10 CONhTIUE
•"•:++
-+11 •-1
M
IF (ALKIMJ.FU.O.0) G41 TO 13
p_ +•.T.MK = T'(B(I)+ALMrl'nO•FLH
Z GRADIENT IS NUNe ZERO, PAGE 10, EQUATI1ON 1.2.10-13)
DrLTA =:ELTAb(,M)(;(T•MIMI/TPK).*(GZ=AMZ/lRSTAR'ALP.i•))))
TI,- +*+++I,
'+I.+=
+ + + + •+:: • ++ • ..
AFAPL-T-h-67-125
Part II
GO TO 14
13 TM
C GRADIENT IS ZERO, PAGE 10, EQUATION 1.2.10-14) 1TB
DELTA = DELTAB(M).*EXP(CGZ-AMZ-DELH/(RSTAR*TMB(M)JJ
14 THETA =TMK/THZ
SIGMA = DELTA/THETA
ALPHA = SORT(THETA*.;3)*((ThZ+S)/(TMKS1))
C CONVERSION TO ENGLISH UNITS
TM =1.8*TMK '
RHO =RJIOZv*SIGMA/GZENG
CA =CAZ*SISRT(THETA)
AMU1= AMUZ*ALP~lA/GZENG
* G T2VERFL(J)II
(1,1),CALL
15 K = K+2
GO TO 1T
16 K =2
17 RETURN
19-
!mos
-.
1
, - Ft ------
I~j
~-~-
AFPL-TR-67-125
Part 1
Al0
AFAPL-TR-.67-125
Part H
I2
AYAPL-TIR-67-425
IISIBFTC
*
CUFAN DCe9/,R
SUJBROUTLINE COFAN
~~COMMNOE/ AJ
421
WORD VIOES ,JUES ,KIDhS 1?"nE 9JNIT ,IDUMP ,IAMTP
21GAS;4XqIOBUAHIAFTBNj1DCD 1INCU .IOShOCrIMSM(JCNIJZLT,
31TRYS tLOOPERA~OMAP kUMM~AP,PtAPE9G,TtlLALL,ERRf6)
COMMOUN/DESIGN/
IPCNFGUiPCNCGUiT4GU DUMPD1 ,LUUM02 ,DELFG ULELFN ,DLLSFC,
ZZFDS iPCNFDStPRFQS #ETAFDS*"AFDS ,PRFCF ,ETAFCF,WAFCF
3ZCDS *PCIECOSPPRCL)S iETACDSiWACOS ,PRCCF #ETACCF,WACCF
'.T4DS ,W~tDS .OTCtJDSETABDSWA3COSDPCOOS,DTCOCF ,ETABCF,
5TFHPDSCNHPDS,ETHPDS,TFHPCFVCNHPCF,ETIIPCF ,OHHI'CF,T21JS,
6TFLPDSClULPDS 1ETLPDSTFLPCF,C4iLPC.F,ETLPCFdrnLPCFtT21WSv
7T240S ,WFDOS ,DTDUDS,ETAUDS,$tAý30S,0OUDIJDSrTDUCFtETADCF9
8T70S *WFADS ,DTAFOS, ETAADSIWG6CDS,DPAFDS,DTAFCF ,ETAACFv
9A55 ,A25 ,A6 ,A7 ,A8 ,A9 ,A28 #429
APS55 ,AM55 tCVDMZ,CVri4NJZ,AaS. V ,A9SAV vA28SAVA29SAV
COf4N0-N / FR.ONT/
JITI.
'14755
1P1 Hli ,Sl vT2 ,P2 ti2 vS2 I
2T21 9P21 iH2i ts21 IT3 9P3 #H3 vS3 I
314 #P4 9114 ,S4 ,T5 ,P5 ,h5 TS5
v^55 ,Ii55s IS55 ,BLF ,BLC ,BLDU tBLOB ,
5CNF ,PRF ,ETAF IWAFC 9WAF ,PWA3 ,WG4 ,FAR4
6CHC ,PR~C
I9AI
,ETAC WIACC ,WAC ,ETAB 9UPCOq ,DUMiF
7CP4HP ,ETATHP,O~iTCmP,DHTC ,BLHP ,WG5 ,IFAR5 ,CS 1
8CtiLP IETATLP,DHITCLP,DHITF ,BLLP IW5 ,FAR55 HIPEXT
,ALTP ,ETAR ,ZF IPCNF ,zC vPCNC WhFb
ATFFHP ITF;LP ,PCSLF ,PCBLC ,PCBLLNJPCBLOBtPCBLHPtPCBLLP
COAXN/ FAN/CNX(15 ),PRX(1 15Jl5)WACXf15,151,ETAXfi~v5.)I
DIMENSION WLH(2)
DATA AWORDWLN/61i COFAN,6H RLU) ,6H (HI)I
WDRDwA WORD
THETA=SURT(T2/!518.665 3
CNF=PCNFI( 100.#THETA3
IF(ZF.LT.O.) ZF=0.
IF(ZF.GT.l.) ZF=l.
CNFS-ChF
CALL SFARCII(ZFlCNFPRFWAFCETAF,
EACUF=ETA. jHE-TA*CF
WAFPCBL-WF.GTO.
W BF=CLW
WIFIJESEU0l)PRGO TAFFWAC7TD
JDES1T(IOA ESG93U PFF*I588 TFC~k58
18?1 WAFCF.E 5.blbG
H TO
6 ,I58
I RsRUF*PF1tI
-
VI
<1
__
£10.TEA*h:
22AFF*TA
_____________
AFAPL-TR-67-125
Part 11
IF(II)FS.Ev.1) GOI TU 4
IF(HMODE.NE.2) Gti iTO 2
T4=GIJESS(3,Yl,Y2,PCNFI-PCNFDS,WF-iiWF~rWSY7,Y~,T40S)
PCNC=GUESS(4,Y1,Y2,PCNFPCNFflSWFBWFB[)SY7sY8,PC"CD)SI
j GO Tfl 5
2 IF(MOOE.EQ.1) GO TU 3
PCNC%1=GUESS(5,T4,T4D)S,y3,Y4tY5,y6,T2lT21l)SPCN'CDSI
I
1 3 G
CCPND
T4=G-UESS(6,Yl,Y2,PCNCPC.\'CDS,Y5,Y6,T21 ,TZIOST4DS)
GO TO 5
T4=T4DSf
WFB=WFODS
T21DS=T21
67 CALL CUCOMP 4
IFINOIMAP.EU,71 PCI-IF=)UMID1
RETURN
END 1
Al-4
' ' 2
18H1 vACCF=tEI5.b#8iH T2IC%=,Fl5.8)
PRC=PRCCFs(PRC-1.)+g*1.
ETAC-ETACCFQETAC
MAC=VACCFiiAC
IF(PCBLC.GT*.o. I LC=PCBLCsWAC
- -- * 8LfVJPC5LOU*RLC
* BLOB-PCBLDS*BLC
6LHP-FCBLt4Pv8LC
L2
PArlFAPL-TR-67-125j
WRITE (8,2000)CNCS7CNC
2000 FORMAT(IOHOCNC WAS= ,E15.8#111H A-%;- NUW= iE15.8i
1.24H CHECK PCNC JN4PUTSsst~s)
CALL ERROR
4 CALL COCOMB
RETURN
END
25I
"74
L
Psrt H
~~l~i-2 ALL/-R-
IE(!orS.;EO.1) A361r;-IA3c
-~ I !Ff74.&T.30O0.) 74.=3000.I
Z UTCO2TI.-73
~ Gt; TU 3
3 oJTf-='TCUC~v0TC0
P3PSltUv?3PSl
CALL SEARCýHf-l . P3PSI NtTCO, ETAB, OUiWA,
iID~.~'.1)GO TO 5
fTAB2wETA8Cf.2ITMI
t~vz((tlU(-.4S9&31*?E-193T')-.2034i1&E-i5J'T4+.2783643E-1lI*T4
I i*[email protected],3JST4-.fi4332'96E-L-i 1.T'.3 1846537E+ii5
CALL TIPAT,7,22bOO~
FAR4c(2*A-f-3fjmftVCETA6J
I.F '~F AR .I(4T. . AX3
26
AFAPL"-'R-67-125
Part It
IF (Mflj)L.NiL.2) (.11 To li
EHHWm (wp5-WFHiX 3/WF0
UI~RmSUkT(WFbi/Wf-SX) J 0 )K
CALL Ag~ilR((Jl).Tg,,&RkwOo,2(*,uooo1,1,~~~I
GO TO) (6,p9,73. I6it
Tr4a'4T
GO ro 1
7 CALL ER4UR
11 wFfIUWFBX 4 01
91 CALL THFI&M((P4,MI.,J4,S4,XX2,1,A-K9
W64swlH4.WA3
IF(IIESEU.1) wR1TIe(6,11O) WA3CI)SEFTARCFL)TC(jCF
ILTAtucFsev'elRt
10) FIJMMAT(17HOCODIRUSTUR 0)FSIGNtXRH IIA3DSUothIS.R8H j
18IH 0TC(lCFwE15&8)
CAL.L C(,NPTB
RETUhRN
ENfl
27
AFAPL'-TR-67-125
:17CfuHP SCUF
6C1HC
8CNLP
9A"I
,PRF
rpiRC
,ETAF ,WAFC
,ETAC *WACC.
,ETATHP,W)ICmP,OHTC
.ETATLPD8TCLPvO$TF
vALTP viETAK 9ZF
jWAF
,WAý
.5185'
tBLLP
,PCNF
IWA3
,ETAS
,5G5
%WG!55
InC
SWG4
IDPCOMt
vfAR5a
,PCtC ,WF8
,FAR4,
ID1.Fli
vCS
PFAR55 HIPEX7
I
I
lTFFISPSTFFII5
- In!EQO GO TO I
RE TUR-ýNMCVNCSRj4
IA~reS=TFHPS-FP.E007FPJG
- ~A8iCJ~iPS-iwWkP
jF ).GT692OO0CN CNSiMAPGO=
?C
if.As -*4hqStflt?fli7L1,GT.G.%*VCftl5IMAG-
-fiTUEXT=0.?06705*!JPEXT
- 5mTECC(dTUI-XT4,.AC*ltK38.1Zl 3(WG'q*T4)
ZF~l0~a~. Co) TO 6
21 _____2S
AFAPL-TR-67--125[
T1FI-PCF=TFHPD)S/TFHCAL
i)HHPCF-=iJHTCC/ DHTCHP
ET11PCF=ETHPOS/FTATHP
WRITE (6,1O?7)CN4HPCF ,TFHPCFTHPCF,O1HPCF
1-02 F(1RMAT(20MLIH.P, TURBINJE 0ESIGN,5X71ICNHPCF=,E15.8,811 TFHPCF=PE15.8,
18H1 ETHPCF=,El5.8,1Z1 DI4HPCF=,E15.8)
6 TFHCAL=TFHPCPF*TFHCAL
'JHTCHP=DHH4PCF -0HTCHiP
ETATHP-,-E T$.PCF*E TATHP
OKTC 4DHTCC*T4
ER ( ) = (TFHCAL-TFF!P) /TFHCALj
ERR(2)=iDHTCC-DHTCHP )/DlITCC
CALL T&ITURB(DHTC, ETATI'PFAk4,H-4,S4,P4,T5,1a5,S5,P5)
IF(BLHP.LE.O.) GO To 7
iAR5=WFB/ (WA3+BLHP)
;=WG4+BLHP
Wi5
115= (BLHPrwH3+WG4*H5)/WG5
CALL THE:tHO(P5.h5,T5,S5t,'X2:11FAR5,1)
I
GO TO b
7 FAR5=FAR4
WG5=WG4
8 CALL COLPTB
RETURN
END
Il
29I
'A
I sIBFTC COLPTS DECKYNh94/29XR7 :)0E ;NT :OP
ill
1PCI4FGJPCNCGUyT4CIJ ,DUMD1 wDM2 DLF 01G
DL51C,
2ZFOS #PCN4F0SvPRF0S ItETAFOStVAFDS 9PRFCF ,ETAFCFrWAFCFt
3ZCOS ,PC?4COSsPRC1JS tETACOSvWACDS wPRCCF VETACCF,WACCF
4T40S ,WFSOS ,DTCOOS,2T&BDSWA3CDSDPCOOSOTCOCFETA8CF,
5TFIPPOSCtIIPOSE~ihPLS TFI4PCP ,CMHPCF EThiPCF ,OI*PCF ,T2OS
6TFLFOS2 CNELPDSgETLPDSiTFLPCFYrCNLP FtcTLPCFvDHLPCFtT210S
TT24DS tWFIOOS IDTDqJDSETAODSPWAZ3DSOPDIJS,0DTDUFETAD)CP,
ST70S ,IaFADS ,OTAF0SETAADSVWGiCOSOPiAFOSOTAFCF ,ETAACF,
I
AAZTZI
A SS
~CWON~q~
R5
9AS
f FROMW
IPZI
, CVD?(3lCV?4NDZvA8SAV
6;.1
,21 S21 4T3
,A9SAV
tP3
PAZSSAVA29'AV
,h3 2S3 I
A6CNC
3T4 tP4 &H4 ?S4 ,Ts 'p ,IH5 VS5
47S5 *P55 3H55 iSs5 98LF vSLC iBLDU .8L06 i
5C4 ,PRF sETAF ,IiAFC ,viAF ,WA3 WG'C4 ,FAR4I
,PRc ,ETAC *UAC(. vWAC tETAS DLPC.4 ,DUMF
7CthHP ,ETATHPj0HTCH~P,DH-,C B8LHP *WG5 ,FAR5 ,CS I
SC14LP tETATLPD04TCLPvDi4TF ,BLLP tUG55 ,FARSS ,HPExT
9AJ( 9ALTP 9EIAR ,ZF ,PCIIF ,zc ,PCNC ,iwFB
1NTFFSINPTTF'i15)
DATA A)IORDIIL0,WMI/6HC(R.PTB,6N (L.0) P61H (hIi)
WQO=~AWORO
IF(DES.E0Ol O) 7010
I ~CV.PCF-CULPDS*SQRT IT5)IPChFf
I C4LPvC~4LPCF*PCNF/SQRT(15)
C)*LPSzCNLP
TFFLPS-TFFLP
CALL SEARCH(-1. ,TFFLPCULPOhTCLPETATLP.
1TFFX11)vftTFPFStC1X(v1 1vDHTC(1 91) 9ETATMIL t1 )ftIPTTFF (Ili 15,15,1IGO!
j 1FfIGO.E0. 1,OR.IGO.E0.II.OR.IGU*E0.2I) V1'RTEf8q10oo3TiFFLPS,WL9
IFfIGO.E0. 2*OR.IGO.E0.12.OR.1CC.EO.22) lWR1TE(8,l0G0)TFFLPS,WHI
IFI IGO.E~k.10.OR.1GiEO.E0il.0R.IGO.E0.12) WRITF(8,2000) CUfLPSWLO
1IFIGO.EQ.Z0.OR.IGO.EQ.e!1.cR.1GG.Eo.22) WKITE(8v2000) CmI~tPSEH
I
1000 FCR1XAT( 19140*****TFFLP OFF ItAPFIC.4q2XA6qIIH*****SSSS5s)
2000 FUAMAT(1I9H0s**b CkLP OFF MA3P, FIO*4 2XA6911Hu....SSSSIS 1
IP(1G0.h4E.7) GO TO 3
2 CALL ERROR
RETuAN
3 RAPBsECLSOL)G..0oNP)IAG=
lF(AaS(TFFLPS-TrFFLPI.LE.O.o-fleTFFLPSI GOS TO 4
j 4
5
a1;(ABSSCNLPS-CftLP) .GT,.O.O01*C'iLPS) PXnPG~s2
1 ,I O2,OAP,~PFYi4APGO.GT.O)
CALL i4APPACt2,ftAPGOTFFLPS*TFFLPCftLPSCft.PPCN~v
I ~IFIjiNmAP.GT*Q3 RETUtfl
r
TFLCALRMwG5*SOP.T(15)1(14.696$F5I
DHTCaUWAF*(iI21-HZ)/V5;5
1Ft1!qS.*EQ.O) GO TO 6
fFL3PCFzTFLPDSfTFLCAL
I
30
APAPL-TR-67-125
Part H
DHL PCF=DHTCF/VjHTCL P
E-TLPCF=ETLPDStE TATIP
WRITE(6,lQ23CNLPCFTFLPCIETLPCFDHLPCF
102 FURMATC2OHOL.P. TURBINE DESIGN,5X7HCfiLPCF=,E!5.8,Sf TFLPCF=,Eb5.8,
1HIH ETLPCF=tEj5.8t8H DHLPCF=tEl5.8)
6 TFLCAL=TFLPCF*TFLCAL
DHTCL P=DHLPCF*DHTCLP
ETATLP=ETLPCF*PTATLP
DHTF=DIITCF*T5
ERR (3 )= TFLCAL-TFFLP) /TFLCAL
ERR (4)=(DHTCF-DHTCLIP )/D;'iTCi=
CALL THTUR8(DHTFETATLPPAR5,H5,S5,P5,T55 ,H55,S55,P551
IF(BLLP.LE-.O.) GO TO 7I
FAR55=WFB/ (WA3+BLHP4.BLLP I
H55=( BLLP*H3+WG5*K55) /WG55
CALL THERMO(P55,H55,T55,S55,XX2,1,FAR55,1)
GO TO 8
7 FAR55=FAR5
WG55=WG5
8 CALL FRIOSO
RE TU'RN
END
* 11
~7J AFAP-l-TR-67-125
Part 11
ITI IP1 t,
,SI T2 ,P2 #H S2 I
2T21 ,P21% IP21 ,S21 ,T3 9P3 JH S3
I3T4t IP4 4 -.S4 ,T5 'P5 ?H 5 I
4T55 vP55 ,H55 ,S55 t BLF ,BLC v,8LVU vBLUB i
5cNF ,PRF ,ETAF tWAKC ,WAF ,WA3 ,boG4 ,iAR4
I6CNC ,ORC ,ETAC ,WACC ,iIAC *ETAB vD PC UM ,DUMP
7CNHP wETATHPIDWFCtIPPDHTC tBLHP ,WG5 7F-AR5 ,CS I
8CNLP iETATLPi0HlTCLPwDiiTF tBLLP tWG5i5 t1AR55 *HP~EXT
S-9AM vALTP tETAR .,ZF t PCifF ,ZC ,PCNC IIF8 I
XT21=T21
C-1 XP21-P21
.1 XH2l=1121
Ij
I
CALL CUDUCT
ENDi
H
'32,
.1 AFAPIT-TR-G7-125
Paft 11
WAX-WAF-WAC-BLF
P23mP1I
H23a (BLDU8*13.MAXSH21 )/W4
CALL THERMO(PZ34123,T23,S23,XX2,1,O.O,1?
BYPASSu(WAF-MACI IWAC
hA Z3C=WAD*SQRT (T233/P23
IF(IDES*EQ*j) WA23DSaIIA23C
DPMXJ'%DPDUDS*(WA--3C/WA230S)
IF(DPOUC.GTo1.) DPOVIC-i*
P24-P23*(1.-0PDU3C)
IF(IGASMX.GT.O) IDWWRNO
IF(IDBURN.NE.O) GO TO 21
T24-T23
W10u0.
FAA24aO.
2O41E1)T4 .2783RoE.2
643 E-113T2 4A
RFfID3Ukif.NE.2) GO TO 6
ERP.~ (WFO-WFDX)/VFD
DXRoSQRTIWFDIXFDX)
CALL AFQUIR(Q(1),T24,ERRW.QO.,OD*.0OOO1,DIRT24TiPlGO)
6O TO (497*5)91GQ
A 33
AFAPL-TlR-67-125
4 T24=T24T
4:
S 6
CALL ERROR
WFD=WFDX
I
~
7
101
CALL THERJ'O(P24,H24,T24,S24,XX2,1,FAR24,OI
M624 =WFD+HAD
IF(IDES.EQo1) WRITE(6,1O1) WA23DS
FORUMAT(12HODUCT DESIGNIMOX8 WA23DS=,E15.8)
C *** IF DESIRED, ENTER OTHER LOSSES HERE
S25=S24
IF(IGASHX.GToOl GO T0%11
WORO=AWORD2
A28SAV=A28
A29SAV-A29
NOO-
IDO=
CALL5)COP2R(-P25)/2P25R S5FR4W2pllN~A8 2R
~PIT28
l0Z.EQ21)2RTE(6,100) A28,AIM28,A29,AM2
j 18 A92,15.2BH AM2CE5g8
CALL F=STOK
GO
RTORI4lO )OO
--
END9T2
H2I8 I
I
?29=I2 ______ I;
TS9=_2
P-9=a2
V29mVr
* A:-APL-TR-67-125
Part 11
35
AFAPL-TR-67-125J
RHO=CAPSF*P355/ (AjeEX55*TS55;
54~ PSZS=psss
-aS25=T2S*a PSZ5IP2S);3*O.286
DO 4 l,;
CALL PRO 15ýPARZ4,TS25,CS25,AK25,pCP25,REX25tP~lSZ5,HSZ251
PHIS2P8!Z25-REX25C-AUGIGP25/PS25)
0EL PHI =PHI3I-P;4S25
IF! ABS (ELPHI ).LE. 0.0001'tPiI S) GO TO 5
4 TSZ5TSZS-*XP(4.0*DELPi4X)
GO TO 2
5 IFUi25.GT.14SZ5l GO TO 55
W'.!TE (8,10ZiPZ5,PS25rT25rTS25,H25,HS25
I
102 FCAMAT(22HOSORT OF I4Z5-HS25 ?4EG v6E15.6r6inSZSU-S)
CALL ERROR
55 V25aS~qT(2.*G*AJv(H25-'-2?5))
RHO-CAPSF?5SZ5/ CAJ*REXZS2*TS?5 )
A25=.G24/i RHO*V25 ý
h3125=Y25/CS25
)QRITE(6,,1OO)4A.5vAM55v-'25vAM25
10
GO ýOP.MA~fZ5e.OTMS~IWEfOULT AREAh DE-SlG~j4,7 A~55`ýi15.Zv
18
ISES.88 A25'I:El5.8:8II AtAZ5-eE15v8!
C-0 TO 20
* C S$* CALCMA~TE P555 AND PSZ5
6 OA~siGf5'A355
C1.=P55'eSQRT (6!lT55tAJ) )*,CPSF
QQ j21 0.87*S
00DO8 I1=115
CALI-. PRCIoM!FAR55,TS55,CS55tAK55,9CP55rREXS55PHl-555.Ks55)
V55=.%155*rC.S55
HSCAL=-H55-Vi55.r-2I (2 40*AJ)
OELMS-FlSCAL-HSS55
IF(ABS(0ELHS).LIý.O.0005*H4SCAL) GO TO 9
8 TS55=TS554OELHSlCP55
GO TO 2
9 WQAT=C1*SQRT(AK55/REXS5)*AM55/(1..41AK55-1.)*AK554.Zi2.)..
1((AK5541.M/2.*(AK5S-*1.M)
AM;X-AX55
IGOG9=0
10 DIP :WQA/biAT
E' =(WCA-tEOAT)/WQA
CALL AF0UJIR (QO(1) tAMXtEhlC.v30., 0-00050IR AMXT ICONI)
GO TO (11,3.5v2),ICO.N
11 IF(ANXT.LE.1.Ol GO TO 13
ANXTr-0.7
IF(ACON.LE.1) GO TO 13
PCNF=Di3P41
WRITE (8,1O3)PC14F,AKXPS55PS5,,P25,PS25
103 FORMAT112H0CORIX PfNF~vF7..4,4H APIZFS.6v5ti P355,F-9.Sv
168 PSS5,9F9.5,51t P25cvF9.S,6H- PS23=,F9.5,6HSSSUSSI
AFA.PL-TR-67- 125
Parl il
PCNF-1 .O1*PCNF
DUKl UIPCNF
12 t]Oe4APa?
RETURN
Ili tISIMOGO.EC0.11 GO TO 14
AM55uAMXT
GO TO 7
14 AM25mAMXT
GO TO 16
15 IF(IGOGO.EQ41) GO TO 19
PS55uP55IEXP( EPHIS5-PHIS55) /REX55)
IF(IGASMX*LE*OJ GO TO 28
WQAnWG26 /A25
ClwP25*SQRTfG/(y25*AjJ) *CAPSF
MCOi4UO
00 42) .0.
QQ(3 )u0.
AM 25-0.0*25
TS25-0.875*T25
16 00 17 1-1,15
CALL PROCOM(FAR24tTS25,CS25,AK25,CP25,K"CX25,PHIS25tHS25)
V25mAM25*CS25
HSCALnH25-V25**2/ (2.**G*AJ)
OELHSn"SCAL-HS25
XF(AfiStOELHS).LE*0.0005*I4SCAL) GO TO 18
17 TS25-TS25+OELHS/CP25
GO TO 2
18 WQAT=C1*SQRTIAK25/REX25)*AM25/I1.,(AK25-1.)*AM250*2/2d*ý0
IC .&I25+I.)/t2.**AK25-1.*)
AM X-A1425
lGOGO-1
GC TO 10
19 PS25nP25/EXPU(PII25-PHIS25)/REX25)
23 WG6aWG24+WG55
ERR(5)m(PS25-PSS) /PS25
WF6*WFD+WFB
FAII6-UF6/ (WG6-WF6)
H6u (WG24*!M254WG55*H55 /WG6
CALL THERMO(1.,H6,T6,PI16tAMJC,1,FARcb,1)
CIsPS55*A55*(1.+AK55*AM55**2)+PS25*A25*(l.4ýAK25*AM2'i..2)
T56=O.833*T6
r'o 23 1-1,15
CALL PROCOM(FAR6tTS6,CS6,AK6vCP6,REX6,PHIS6,NS6)
C2aWG6*SQRTIAJ*REX6*T6RiAK6*G~)
C3'wC2/ (CAPSF.*C1)
f4fIAiK6-1. )/2.-(C3*AK6)**2
C5,1.-2.*AK6*C3**2
C6u1:5**2+4. *C4*C3**2
XF fC6) 21, 22i-23
21 CALL ERROR
RETUJRN
2? AIN62=-C5/12.*C 4 )
GO TU 24
23 AM62-(SQRT(C6)-C5:/(2.*C4)
24 F(At62.LE.o.3 GO TO 21
AM6uSQkTfCAM621
V6aAM6*CS6
HSCA~.mM6-V6**2/ (2.*G*AJ)
DZLHS-IISCAL-HS6
IF(ABS(DE~ltS)*LE3.OOO05*HSCAL) GO TO 20
25 TS6=TS6+OELIPS/CP6
38
UFg
PAtFAPL-TR-67-12.5
GO TO 21
26 IF(IGASMX.EQ,2) GO TO 27
A6=A25+A55
27 C7=SQRT(l.+(AK6-l.)*AM62!2.)
PS6=C2/ (CAPSFzA6*AM6*C7!
P6=zPS64EXP( CP~-I6-PHIS6)/REX6)
CALL THERHO(P6,H6,T6,S6,XXl,1,3AR6.-O)
S6AVE= CWG24*S25+WG55*S55) /WG6
IF(S6 3 C.E.S6AVE) GO TO 29
S6=S6AVE
W ~P6=EXP(A4MX*(?i-16-S6)/l .986375)
GO TO 29
28 T6=T55
I P6=P55
I
H6=H55
I
S6=s55
WG6=WG-55
PS6=PS557
V6=V55
AM6A1455
IF(ICASMX.EQ.O) A6=A55
29 CALL COAFBM
RETURNJ
END
-------
- --
* ~APAPL-ifl-61-i25
* Part n
J
9A55
I
,A25 *A6 ,A7 tA8 ,A9 I.8
A2-, t
AFPS55 ,AXS5 vCV0NUZ*CMfl0Z,A5SAV 7A9SAY *A~eSAV,A29SAV
-. OR / BACK/
tT5
p5s ,H55 1!5 PZ5
-T .1425 9S25 I
Z2dF8 ,MG55 vFAR55 d'F.F ,1424 PFAIR24 P1I ,CUm3
3T6 P~6 IH6 256 ti ,P7 vH7 vS7
.8 f
'Ps vi4m 5 ,T9 ,pq ~I49 's9
5ifGb IWFA tWG7 .2F- A7 TAA tDPAFT 9V55 #V25
6PS6 ,V& #A."' tTST *PS7 tV 7 ,A117 rAM(25p
7TS8 #PS8 vys MA1t tT59 ,PS9 ,v~t IAM9
SVA v
F~RO 7VJD OW FG~R
IV1 ,FiGP0 ,FGPX .4
9F:GM IFGP' ,WPF .i;T tFART FIG ,FHa ,SFC
DLIKENSW0l Q!91
DATA AWflRD/6N4COAFPN/
UMN~DAWORD
0123=0.4
Q(31-09
AJlL778.26
I
I ~CASFU2116aI2o
~G-32. 174049j 4
f 11F62-WFB
J C
IF (IGASP.X.GT,01 WF i=bI?&+WlFD
WA6-.ev6-WF6
'LOSS
DAY
* 2 Wr-6C'-WG6*SQRT (T6)/P6
j DPAFT=DPAFDS*fWG6C/WG6COS)
IFIDPAFToGT.1.) OPAFT=14
A7mA6
I ~FAR6zidF6fkA6
* '
CALL PR03C0M1FAR6,t6,XXlWXXZtXX3,XX4,PXil6,XX ;
6
WOAwWG6/Atf
A IC1=tP?*SQRT(G/('!-6*AJ;)*CAPSiF
:S7=O.B75*T6
20 00 22 la1,15
CALL PPOC014(Fý:6,TS7,-CS7,AK7,CP7taEX7,PHIs7,ý-S)
YVsA3.t?OCS7
I ~HSCAL=K46-V7..2/Q;.vGcAJ)
0DcHSUKSC.AL-HS7
If (ABSSWELHSJ.LEo.0.0005*HSCAL) GO TO 24.
22 TS7=TS7.OELH~S/CP7
40
AFAPL-T-.--67- 125
Part 11
DIR=WQAIWC.AT
EW= (WQ~l-WQATi,fQA
CALL
GO TO (26t28r-s)v GO
26 AI17-AT
GO TO 20
28 PS7=P7/EXP( tPt16-PMlS7)lRvEr7!
!FUiA-F-T8.GT.0) GO TL, 4
MG**
-AFT-BtJARlNG
WFA=0.0
FAR? =FAk6
60 To 11
C 0** AFTERI~8Rfl1G
4 F11IAFT8N.&l-i21 1-7=T6+2000.-
IF(T7.LE.T4,) SO TO 3
P565=PS7
V65=V7
Q(3)=O.
5 IF(T7.GT*4000., fl'4000.
C l*F DESIRE&, ENTER UlA ýIATI-UINS FOR STAA HERES
EPa(PS7-PS7A)/PS7
CALL AFQUIR(Q(1),VTEPO.,50,.,O005.DIRV7TIGO)
~~1 sV7=VT
IF(V7.LT*100.) V7=100.
I 12
GO TO (llt12v8)i!GO
P7=PS7*EXP(CPH17-PHIS7)/REX7)
~CALL PRQCOM(FAR7,TS?,CS7,XX2,XX3,XX4,XX5,XXC6)
Ati7*V7/CS7
END
.21A
PArt ll1
DATA AWORDI6H~wNOZZL1
NORD=AWUORD
*
IFtNOZFLT.ECo.1.Qs.N3ZFL'T.EQ.3) tvOZN=1
IF~~co;~,1)GO 70 1
CALL COftVR~iT7,w11:97,?S?,FAR7,WG;PPP!pP1UWzArl
GO TO (393,3,2)tICCGN
p7 R*
;
CALL C~tuDJVT71,7P?,S7,FA.R7,MG7,P1,rXN4OZA8,ArP7Rt
1$SH{C<=CCJU
SO TO f4r414#219Xcoh
2 CALL ERROR
3 T9;r3?
TS 9=TSS
PS9P8
I)ISHflCm1C0+3
4 ERR (6 1P7R--!27 i/P7R_
JFt1RNGZ.Eo.10 FI4TE16s100) AatAM8,49,AMIS
100 FOA:4AT(14I-IOOZZLE QES1GNv,)08H Ar~mE15w81 3H BX*Ei5*.8,
isji A.153~#
dx,S15.81
RE TURN
43
AFAP2L-TR-67-125
paft u
I
j
7TT2'4tS tWFLVOS ,DTDUO3,ETADGSWA23DSsDPDUDS ,DTOUCFETADCFt
8ST70S ,IFADS tOTAFDSpETAADStWG6CD-f-,DPAFDSeOTAFCP ,ETAACF,
9A55 #AZS #A6 IA7 *A$ IA9 A29 A29 I
I
IAPSSS ,9'(55 vCV-aWZCVKNfOZtA&SAV ,A9SAV ,.AZOSAVPA29SAV
COMON I FRONIT1
2721 IP21 .1421 ,S2I rT3 ,P3 r.H3 tS3 ,
374 tP4 9H44 so# ,T5 IP5 %145 ,S5
4TT55 vP55 INS$ ti55 paBLF tBLC t8LOLI iSLOS I
nnr- #PRF tETAF tKAFC iWA*4 rWA3 iWG4 tFAR4T
6CSK *PRS PETAC *bVACC tWAC ttETP- *UPCO- vDWIF I
XUR tCI
E-TAth'P vD1TCHPvC14TC B8114P ,weGS IFAR5IS
9CI&P vE'!ATLPvEffTCLP 20MTF ,BLLP 9WZ55 9FAk~5 ,HPEXT
9Am sALTP tETAR ,ZF PPCHFtZC pPcNC $WF8
A-PFHP *T-rLP *PCSLF ,PCBLC vPCLS~tV-8L0B:PC8LHPt.PC0LLP
CtioG$ / SIDE/
XXPI tXWAF tXWAC *X,XSF :XT-LVD'j Z."13 ,XJMiSI sDURS2a
XXT21 vXP21 .11421 vXS21 vT! ,P23 .1423 VSz3
3T24 FPZ4 2H424 IS24 9T-55 tP25 W"425 tS25 #
0283 ,P28 .128 IS28 'T-9 IF29 ,1429 PS29
SWAD ,IFD, IGI-44 PFAA24 IETAV PDPOL" tBYPA'SSDUiISD
6TS28 ,PSZS IV24 ,Am2a i1S29 #PS29 #V29 PA1429
ce-ACi
XXIS5 ,XP55 ~z1435 iXS55 vXT23 PXP25 PX1425 tXS25 v
X^-VFB trvtG55 ,.XFARr-5praF0 PXWG24 97FAq249XXPI ,OUM!B 9
VA WA-~l
*--
* ~ 6?I'CS~S'~-P)qA9
-rc I
r ~ 44
=W -
A FAPL-TR-67-125
FGPDG=CAPSFI (PS29-Pl1)*A29
1 FGH=FGMM+FGMD
FGP=FGPM+FGPD
FG=FGJI+FGP
FN=FG-FRD
SFC=36 00. tWFT/FN
I FG=DELFG*FG
*1
Ffl=DEL FN*FN
SFC =DEL SF C*SFC
CALL OUTPUT
CALL ERROR
RETURN
END
A' %
4 j KFIIDES*EQ.ll SO TO '.
IF(IOSUP)E.TC*) GO TO 2
I RT1sO)Mr
* )-(iAFTBN.61.O) GO TO I
*J IRITE (6,1OO)WORDA11IALT?,T4,E7Aft
IG TO 3
_ O-O
AFAPL-TR-67-125 -
2 WRITEC6,1i12)WOROA~hALTPT4,T24,ETAR
3 CALL CONOUIT(2)
4 WRITE(69104)(W(IHSHOCi ),I~1,4),FGFrNSFC
IF(JGASMX*GT,0) GO TO 5
5 WRITE(6t106)LOOPERI
WORD =AWORD2
WRITE (6,1O7)WORDZFPCNFtZCPCNC 1 T4.l4ODE
WRITE(6ilOS)
WRlTE(6,109) (ANSl(I )v,=lt80)
WRITE(69108) I
WRlTE(6,109)%'ANS2(! )lI=1T8O)
WRITE(6,108)
WRlTE(6,109)(ANS3(I )#jI.4S8)
WRITE(61108)
WRlTE(6,109)(ANS4(I1),1=1,72)
IF(IDES.EQ.1) GO TO 6I
* A8=A8SAY
A9=A9SAV
A28=A28SAV
A28=A28SAV
IF(IDUMP.NE.2) GO TO 6
WRITE (6 ,110)
CALL SYG(2
6 CALL ENGBAL
RE-TURN
16X?1H T4=,FB.2,25X7H ETAR=IF7.4)
101 FORMAT(iHBA6,i4X7H AMlF7.316X7H ALTP=IF7.0,
16.X71 T4=*F8.2,5X7H '17=,F8.2,5X7H ETARcpF7S.)
102 FOR'4ATiIIIBA6t14X7H AH=rFT.3,&X7H ALTP~vF7.,O
l6x-fH T4=,FB1 2,54XH T24=,PB.2q5X71I ETAR=IF7*4)
"704 FORHAT(6HOMAIIJ ..4A6t9X3HFG=9F9.L'tlSX3HF.%=tF9.2,it3X4HSFC~iFS95)
~5 FORMATi6I1 DUCT 1 4A6)
:~FOPJ4AT(16HlCDNYERGE0 AFTER0i4t6H LOGPS,/#lHB)I
j:WZ FORMATMH 7A6vXv5E15.6v14)
"S14 FORMA7,(!H)
1. FORMAT(1H iSE15.6)
END
471P
ma-
"FAPL-TR-37-125
j~~1 1PCHFGUvPCNiCGUtT4GU
~2ZFDS
3ZDiPNORO ~A0~1CSRC
,PC.NF0SPftFS tETAFOStWAFDS
tODrnDI .04.102 POELFG #DELV4I tOELSFC,
vPRFCF #ETAF~C~eWAFCF
EACAC
4T4!3S tirrSS ,OTCDDS.ETAIBDSWA3CD$,DPCODSDTCD)CF.ETA~CCP
5TFMtPS*CNPOSETHPOSTF14PCFCNH14FETHCFD1*4PCF2 TZDS
6TFLPDSCNLPDSETLPOSTFLPCF ,CNLPCF.ETLPCFD14LPCFT21OSv
j 7?TZ405 ,ikFDOS ,DTOUOSETAOOSWA23DS,DPDUDS.DTDUCF.ETADCF,
4 1C0ORRO
8TMtS ,IFADS ,TFS.TAS ~ ,~ ~I~,TAP
9A53 IA23 ?A& ,A? SAO A9 *A28 PA29
APS55 .51155 ,CVOWI)ZvC%-POZvAeSAV 9A9SAV tAZ21SAV9A29SAV
/ FR~ONT/
ITI 9,PI .11 *SI MT VP2 .142 ,SZ 9
i:TZX P21 .1"21 qSZi '73 TP3 .143 IS3 r
3T4 ,P4 .144 vS4 ,T5 'P5 ,N5 llS5
4T55 r-S5 71455 IS55 ,B!.F *BLC TFLDU ,Bias I
5CRF tPRF ETAF rWAFC vWAF *v1A3 tbiG4 .FAR~4
6CF4 ,PRC ,ETAC ,IiACC vWAC ,E'-Ae ,DPCOII SDUMlF I
JRhr' 9ETATMPtD~iTCHP*OHTC veLKP vWG5 FPAR5 vCS P
tCtILP *ETATLPv3iTCLPO14TF .SLLP ,tG55 iFAR55 91tPEXT,
9A11 #A' ,ETAR ,ZF tpcIF ,tzr 7PCI4C tliF8
AThi~t' qTFFLP ,PCBLF ,PCSLC ,PC5i.DtJ.FC--BLO8,PC3LlHP.PC8LLP
COMO"1 I SV1DE-1
XXP1 *X1WA-@ #XWAC #XEIF PXSLIYJ wX83 ,OUNSI tDIRSZ#
XXT21 tXP21 *X142i 9AS21 YT23 *P23 91423 vS23 f
3T24 9P24 9H424 ,S2'. ,T25 PW-5 T"425 SS25 9
4T29 fP28 91428 IS28 ,T29 'PZ9O .129 IS29 0
5)itA0 ,VFo 9UGZý tFAR24 tETAfI wDPDUC POYPASS.04jt1339
6TS28 lp...Z8 ,VZG .AMZS ,TS29 #PS29 vVkv rAA29
CO#2It I SAWI
X~XTSS *XP55 tX1455 ,XS55 ,XT2S ,XP25 vXHZ5 tXS25v
* XF5 *XIG55 tXFAR55tXW~D 9XWG24 tXFAP.24*XXPI .DUme ,
3Tf- P6 H86 *Sfi '1P7 WHY7 *S? I
'.18 'Pe .1 He as$ f9 ,p9 .149 T59 t
5WG6 WIFA ,WG7 ,FAR7 vETA .0PAFT #V55 ,V?5 v
6PSG ,V6 ,AM6 T1S7 9PS7 ,Y7 9AN? YAK5 91
TTSB VSB ,Ys ,AP18 9TS9 ,PS9 MV oa49 i
OVA FRO ,vJD vFGO qVjii #FGKM1 eFGPD tFGPPM
WCA tFGP .VFT 'WGT vFART ,FG ,Rh tSF-C
DIM1ENSION4 PAaA91(280OhWMDY(2'60),IOUT(1O3) ,AOUT(6).IIOUTC6)
EQUIVALENCE (P'AN.HPCPAFGU)
DATA (WOR.DY~fI),I.19sl
16i4PC5FGUI 96HPCNCGUj, 6l$T46U v6HDUXIW , 614)M02 v6140E1FG t6mDEL FN v
?6*M1ELSFCv.HZF0S ,6KCX-fOSt64PtPFD3 .vKETAFDS96HWAF05 6H1PIEFCF
36tlETAFCF .6'4AFCF ,6H1CDS v6H14PC3CDS-6HPRCOS o6:4ETACDSv,6-HWACDSv
~.4WRCCF v614ETACCF961i.ACCF ,6aiT~tvS i6H%'FBDS #6hDTCOOStWCETABDS,
56MfA3CrS , 6XtIPCODS ,6UOTCOC F 6K-TL2CF.614TFKPDS .68ýCNHPDS .61ET14POS,
I
66)ITFhCF ,64C*I4HPCF , L4E r9CF,6I~b*h.ýCF ,6HTZDS * 6HTPLPOS ,6"CKL PD5
T6#4ETLPDSv6"TFLFCF v6HCWLPCF 6#-ITLPCF 96ND1LP0 961IT2 IDS t6tiTZ405
8614WFODS .6IiDTDUDS, 6HETADDS ,61UA230S,6H0P0U0S76WTOTUCF ,61ETADCF.:
96,1T70S ~61f#F-ADS ,614TAFO1S,6ttETAADS*6trJWDS ,6HPAi 0! t6KOTAFCFi
16trETAACF 61t555 1614525 ,&14A6 16mAT .61458 t6149 p
W4AZM 06)iA29 v6E19S55 t614A155 v6W-YD'QZ,614y1W0lv611A8SAV
36NA9SA v6JA28$CAV.#44529SAVs61471 96#PPI N61*41 v683 t
WI4142 ,61HP2 961mc2 06KS2 96HT121 ,614pzl t6HH*I
I4
56NS21 stfT3 ,6HF'3 i 6HH3 96HS3 ,6HT4 96HP4 I
DATA (WOROY(1,1=99,189)i
6,$KH4 ,6HS4 t6HIT5 ,6HP5 t6KH5i 961S5 ,6HT55
TSHP55 96KH55 t6MS55 %%tMLF v6HSLC t6H1310UI t6KBLOB
86HCXF tbHpfiF ,6HETAF v6FVAFC ,61iWAF t6HWA3 ,6HW4&A
9654FAR4 vbHCJNC 06H4PRC 964ETAC ,6VMACC 96HILAC ,61IET^B
lI&MrCOg ,6HD!JJ`F 168C."*P ~6HETATHP,6HOHH~P,6h0DHTC 96HBLtiPt
26HWGS 1 6i4FA!5 V6HCS 96KCW.P ip6HETATLPv6HD)1TCLP 1,68OtTF -#
36H3LLP ,6IWG55 ~6KFAR55 4'HHPEXT 63iAH ,64ALTP t6H1FTARt
46HZF ,6HPC.NF 6H*ZC t614PC NC t6t*VB jidHTFFHP ,6HTFFLP
561IPCBLF ,.,l1PCSLC ,6HPCBLOIJ ,61IPGLOB, 6HPC&LSIP 6,6-C5LLP ,61LXP1
66$ltXWAF~ t6HX-WAC t61IXSLF *6HXSLOU *6HrXH3 t6H0Uj451 s6HDWMSZ
76HXT21 t6tO -21 ,6hXI121 t61IXSZI 6ilTZ3 96KO23 96NH23
8SAS23 v6H124 ,6MP24 v61*424 ,6HS24 614T25 y6IKP25 I
I36HAA128
16HP29 96M"i29 96HiS29 ,6MYADo 96HWFD s6RWGZ4 v6IFAR24v
26HETAD -6HDPOV. 96;48YP:SS96HGUMS3 ,6HTS28 9,61PS28 16HV28
.6HTS29 ,64PS219 6HV29 ,6HKA29 ,6HXT55 ,6HXP55
46HXH55 t6HXS55 t6#4XT25 96HXP25 96KXH25 v6HXSZ5 v6HXWFS 7
56HXW05c5 6KXFM5S,6H4X%3IF ,614XWG24 ,61XFAR.24t6HXXPI ,64DJXS
6614T6 ,6HP6 96.1*6 t6MS4& v6T7 w6)1P7 I16HH14
76HS7 t6HTS 96ws~ v6H88 .6HS8 p6HT9 T61'$9 ,
86HH9 16HS9 r6HWG6 v61YFA ,64MG? ,6HFARl7 96METAA ?
96HDPAFT s6HV55 5 6!4Vs ,6HPS6 ,6HV6 ,SI1A2I6 ,611TS7 I
16HPS? ,6hV7 v6tAX7 96?iAK25 t6)ITSS t6HPS8 ,6H'/8 gf
26i1A11 ,6HTS9 1 6+CPS9 ~6HV9 ,6HAN9 ,6HVA p6KFrZD
36)WJD tbHFGMID ,6WEvjm 96I*GMM1 96HFGPO v6hfýGPM 96)iF~g
46?fGP :,6edFT v6HWGT 16N1FART ,6HFG ,&IFN ,61SFC
t /
DAIA THE ,81SAF.K,LlMJT/&HTHEEN0,6H ,2210/
GO TO 1191219ICON4
C **INPUT SECTIONE
READ(Stl00)AJNCHAP4GE
IFIAIN.E-Q.TMEEt4!) SO TO 5
- 00 2 JxlLIMIT
1F(A11I.EQ.MRORY(JJ) GO TO 31
4 WRITlo*11)IN-
3 **IOUTPUT S'Tlir
12 IFiHU"Xi.NEO.1) TO 18A4
1;Xj VMO
k:
12 IFZWN.GT.6) C0O M16
JaN
13 NaI4-6
DO14 KzItJ
Nr-it-UT(L)
WOJT Ms:=WORDY (9)l
14 AOUT(K)P1AUM(1)
IM
- --
16 RETso
7-7
- - AWE
A~~ 74-
7-, ~ -.- ~--- ~ - - ~ - -.-- -
-
I'
-
DRIMENSO
(.0TR-ASHJ(Sorf IIH180J#~S3'8)TAH(2
WRITE(6#105)GR
WRITE(6,0)tlh2ssi'(J 1,1=1C,.803 4,KM
WRITE(6t103)
WI'1TE16,104)ITRASn,(x 3.1-1,48)
WRITEE692051
hRlTE(6vI0f.)(TRASH2(I 1,1-1,72)
WRITE(6t1031
I5
WRITE (&,106)LOOPER
CAL -, -'W - -'
- - z ~~ -RETURN- -
13FRMETURH0
104 FOMiATMK0,8E15.6)
105 FORMAT(IH1)
106 FORMAT125HOFAILED TO CONVERGE AFTERfI416H LOOPS)
~EN
I5
- ~ ~ sz
- ~ -- ~~ .- ~ - ~ ---- ~ -A
SIF2 --- 0EC-tA94
'. 9 --
IT N -UR - 4- 7 _- - -
j REWIND 8
RETURN
~ jC TERMINATE THE FILE
2 WRITEi8,500)
500 r-ORMAT( 12H$~$eSS$Ss$Ss)
END FILE 8
REWIND 8
c READ RECORD
5 READ(P8,501)(4ORDfIIII=1,j32)
501 FORMAT1132A'j
C CHECK:i
FOR;1::L::: ING:DOLLAR SIGNS
10 CONTINUE
RETURN
IF'WORD(I )-OtNEDOL)15#l2#15-
12 K=;+5
DO 13 J=I#K
IF (WORD(J)-ONiEOOL 15tl3v15
13 CONTINU'EINYG
C PRTNT LINE
20 1-1-1
WRITE(6v501 I (YORDC'!, ,M=1,11)
GODTO 5
ENL
53I
A
- ---- - - - . -,,* - - -
. --. - ,flt r- S - - - - - - - - - -
- -. - - -- - '-- --- -. - - .- r. -
- -
I $H4F TAP$ cK
eITRY *UNC$b
.tJftO8. PZ UNITO8
II I
I
I
Lfl4TO8 FiLE tC(1),BCDBLKZ2,REAOs', ING-UT
---
£
II
I
II I
I
I
f
I
I
-.
-
I
I 1-
f--I
II
I
j
-
S
- j -
1
* I
I j
-1
'-<I
II___________ 54
-
A
-g__________
77--.~ '
I
;IEBFTC THCOIP O)ECI-,3,M4/42,XR1 A
SUBPP0BJTIFIE THCOMP(PR,eTA,THSPTQHOSOPO)
P04' *PR
TP=T*PRZ*06 28572
DO 1 1=1#25
I ~CALL T!1ERitO(POIHPTPiSPrX1,Ox2,O)
OELS=SP-S
JF(ABS(OELSJ.LE.O.00005*S) GO TO 2
I. TP=TP/EXP(4.*DELS)
CALL ERROR
I
2 HO=H+C(HP-*4)/ETA)
CALL TH~eRMO(POKOdTO,SOX1,OX2,1)
RETURN
END
IA
51
Z IOPTC 1rHT4JIB O)ECK,994/2tXR7
SUBROUTINE THTUP.3(D.u,FTA1 FAP. ,NS P&TO,140,SO PO)
HOPmi4-DH/ETA
PTcP/2,
DO 1 1=1*25
CALL THERhJ(PTKOt)PTTISTrAMWT,1tFAR,l)
DESST-S
A
7
-6-M -1 - -N-
FAP-r -r
IL
It
57I
A AiL- TP -ý12f
TEA=30((1.602E-6TXLO552Z.*4.
TEW5RITE-5(TE84102O2E'O1)*E
GO TO 3
GOTF(OEXLE.4000.E
3IFFAR.GE.*0)GO TO05
4 CPA =(;(7ft267B7150E~-25*TEýX-1l.435b677E-zo'-.*TEX
1+1.O6212913E-18&)*;TEX-'4.52o2a25E-14)naEX-6.96a78376-io-4T)
2-L37119O1E9-Cb*EX+1.22X-5.586o-368*TEx.0.?*3EX16638E-o2-
HEA=(49(t.264ý;-4288E-.26-*TEx-b-9O752O949E-21)#TEýX
1+1-27021525E-17*9Tg.X-3.02565E-14-*TsEX4Z*-.467945E-1Jo-).EX
22-.58982633E-07 IPTE-X46.129315QE-O4 I15-TE-X+7.38r102005E-O2!)*E
3TX3-.75!81530E+O1
SEA+?-. 502005'E-013;:LOG (TEX)( it (C 1.4450767EwTE
1-2.4221128--2I)*ýzTEX+1.2.43532E-17)*-,t-EX--1 8ý04E15*
'3+4.5432309E-02
5 CEZ(CPARX.L X*CPF0
GO+TO 5
CSEZSQTt:((AKEX*REX*TE-*265OEXI.37) 94E21)*E
3 + 6 . 4-3 --- 0 - -
5 CPX=(CA+F.RX-CPF!-011*+F--X
i~-o Ri 2EtTooas
Asc;.cI
-
NZOiJE=01 A~ LO
c OrJcric 27 E4RPOP-
C NJGrm . Lo
CZO.
C F*1~ A
!F~.T~A~1i1G-- TV 2
i CONTIRUE
!F(IA.GT.~AX!9) NCODE=2
2 O TOr(IltýGT.1!
WCODE -1 GO TO 31
A=AXi13
?M-4;A-AX(;GI)I(AXIH)-AX(1fL) I
!F(P.GE.O.) GO TO 6
l~t.GEaL)GO TO 4
GO TO%5
IFSLk,' GO TO 5
pM=O.5
lf(BH.LT.BX(1H*J)) GO TZ 8
7 CO14r lINU
8 JL=JH-1
DO 9 K=ZvLIML1
KH=K
IF(BL.LTa;8XtILK?i3 GO TO 'to
BT.BL4PRH*(BH-BL)
CTi=CL+PRM*CCH-CL)
j ODTZDL4PRM*( Oll-DL)
IIF (P.GE*O) GO TO 13
D~IR=SQRT( B/BT)
CALL AFQU'ý(Q(l) ,PPER~,O.,25.,O.OO1,D!RPTJC0I4s
GD TO fl1,13t,13)ilCONI
IF (PP*LT.O.) PP=04
GO TO 6
NCODE-7
C=BT
D=DT
j END
AFAPL-TR.-67-125
Part 11
61
AFAPL-TR-67-125
Part 1U
TSS=O.833*T I
1 J3J41
CALL PROCONM(FARTSSCSStAKStCPtREXSPHISStHSS)
HSCAL=HI-CSS**2/ (2.*G*AJ I
DELHSnHSCAL-HSS
IF(ABS(OELHS)-O.O005*HSCAL)4t4t2
2 TSS=TSS.DELHS/CP
3 1C0N-4
RE TURN
4 IF(IDES)12t12,5
C **ISENTROPIC EXPANSION CALCULATIONS
5 J-6
TSImTl*(PA/PI )**O.206
6 JSJ+1
CALL THERMO(PAHSI,'(61,SSIXX1I,~FARO)
IFfABSISSI-SI )-0O00i1*SI 18,8,7
7 TS!=TSI/EXP(4.*(SSI-SI))
IF(J-30)696t3
a VISmSQRT(2.*G*AJ*IHI-HSIJ)
IF(VIS-CSS)9t1i,11
C **SUBSONIC DESIGNt CALCULATE AD
9 VonVIS
TSO.TSI
PS0-PA
CALL PROCOM(FAR.TSOCSOPXX2,XX3,REXtPHISOtHSO)
RHOuCAPSF*PSOI IAJ*REX*TSO)
AOoWG/(RHO*VOI
!AHO=VoCSO
PRxPI
ICON-i
10 TOuTI
HOmmI
P0-P I
SOUSI
RETURN
C.*** SONIC OESIGN, CALCULATE AD
IL VO-CSS
TSOuTSS
PSOuPI*CTSO/TI)**(AKS/(AKS-1. II
RHG=CAPSFOPSO/(AJ*REXS*TSO)
62
W-4
Pa~rtf11
AO=WG/(RHO*VUI
AMO=1#0
PR=PI
* ICGN=2
GO TO 10
12 VO=CSS
TSO=TSS
PSO=PA
RHO=CAPSF*PSO/ (AJ*REXS*TSOJ
AOCRIT=WG/ (PHO*VOI
AMO=1 .0
I ~~PR=PSO*(TI/TSO)**(AKS/(AKS-1
I IF(AO-AOCRXT)13t13,14
14 PSO=PA
J=O
TSO=0*833*TSO --
15 J=J+1
I
CALL PROCOM(FARttSOCS0,AK0,CPREXPHISOHSO!
RHO=CAPSF*PSO/ (AJ*REX*TSO)
~ ~VO=WGI
fRHD*AO I
*ISCAL=HIl-VO**2/ (Z,*G*Ai)
I ~DELHS=HSCAL-HSOJ
16 IF(ABS(DELN.SJ-O.OOO5StSCAL) 17,17,16 -
16 TSO=TSOiDELHS/CP
17 AMO=%YUICSO
I - PR=PSO*(TlITSO)s*IAKO/(AKO-1. l)
I
ICON1
GO TO 10 -
END
T FAPL-TR-67-US!
-, - - 032.174049
C e*SONIC CALCUJLATIONS
T!hSsO.S33*TZ
CALL PROC0N(FAATSSCSS.AKCPR*XSPHISSHSS)
DELH$.44SCAL-RSS
IF(A9S(DEtHSJ-O0oOG5*H CLI494t2
1 2 - TSSaTSS*OELHS/CP
IFtJ,15)1*1*3
4 IF(IOES)11,11,5
I C **
5~ VTm-CSS
TSTwTSS
-
M~iIC 6ESIGN, CALCULATE AT
P.STsPISITST/fl J*"AJC/(AK-1. 31
:HSCAPSF*PSTI (AJSREXS*TST)
ROU/HOn
-~
IISP
I ICOln 3
1 9. TO-TI
HA!I~
P0-fl
SO-SI
part 11
10 Tl-TI
HTwHI
I
PTwPI
ST-SI
RETURN
C * ASSUMlE SDMIC MhOAT AND ISEXTROPIC EXPANSION TO AD
i1~ VTaCSS
ART-1e.O
TSTxTSS
RHO-WG/IAT*YT)
PSToRHD*A.IREXS*YST/CAPSF
PlRwPST*;TIITST)**(AKI(AK-l.)J
IF(PST-PA)12t2*924
j 12 TS~xG.95*TI
I 13
~
CALL PROCOK (FAP .TSOPcSOAKCPREXPHMISOHSO)
~AKNiSQftT(2.*( (TJ/TSOII. J/(AK-1.))
AOCALU-iAT/ANO)*(2*1.*gl.fAK-1.)*ANOS*2/2.)/(CAK+l.))..((hiK+1Io/
I
-
1(2.*(AK-loJi3
EAw(AO-AOCALIIAO
.
DIRmSOfT(ACIAOCAL 3
.3 CALL AFgUIR(Q(l-,Ts0.EA,0.,1OOoO.0oIOIRTSOTsJC(N)
6O TO (14tl5,3),JCOM
14 TSOxTSOT
IF(TSO-TJ)l4O,13tt41
IF(TSO*GTaTSCJ GO TO 142
21 TSOMO.9S*TI
12 IF(Q(2hoLT.30.O.Of.ANO.LT.0.95.0ft.NANEQ.1) CO TO 13
GO 1A3 TO
IFISO-O.MP33*TI
16 L P~AMOCSONFRTOCOACPEXPSS)
TF(SS (DE3*LS)OOOSAL2201
19 TSsTOCVLSCW(APStCtAttRtPItSO
L
2AOVOiCSOHGAD
OPLI -SCa-SO*TiO1*1K(~1
5 1SDLS)OGSI"L2901
19 TI0mS+DLS
IFJ2181#
II0 AMw3/
t Pift n
TSTwTSO
21 CALL PftOCOM(FARtT$T;CSTAKtCP,REXepx;-,HST)
PSTsPIR*(TST/TI)**iAK/(AK-l.))
- ~HO4ST#CAPSFl(A4*AEX*TSTP
A VTxNdIRMOSAT)
MSCALwhI-VT*.2/t,*G*AJI
Eiis (I4SCAL-IIST /HSvCAL
I
DIR-i .*(KSCAL-HST I/(CP*TST)
CALL AFcUIaqQu3,vTs1EHtO,,zo.,o*oDO,:c,:RTTsTrJCONJ
IL
-~c 6O 70 (22,23,33,JCO9N
22 TSTMISTT
GO TO 21.
24 PSSOtPA
J60
TSD=TI*CPSO/PIRJ**,286
25 J=J4.
CALL PR0COgq(FAftTSOCSOA~,CPtREXgPHISO*HSOJ
PHICALaPHI I-fEXOALOG(PlR/PSOJ
DEaPHI=PHICAL-PIIISO
* IF(ASSIDELPHIJ41-O.OO0*Pki4CA.1)27,27226
4 26 TSO=TSO*EXP(4.O*DELPHI)
* 2i7 VOsSOftT(2.*G*AJ*(HI-,tSO))
AIIO&YO/CSO
-;AO1Dau(AT/AM3)*(2**(l.+IAK-1.)*AWI**2/2.lI(AK~l.J)).e(f AK.1. I
29 TSO=Da*33*TI
30 43.1+1
CALL PROCON(FARTSO.CSOAKCPREXPHISO,14501
AM AL 0 1AT/ANO1* (2a 1.+ (AK-I)*AHO**2/2.JI AK+I.l.3 1~AK+ I. I
DELAsAG-AOCAL
31 lTsOjTSb*SQRTfAOC.AL/AO
I AS(EA-.O~A)33Z3
IFlJ-50)30t30t3
32 IF(N134#34933
c *s. WE-IEXPANDEDt SHOKII OUTSIDE NOZZLE
'33 PS~inPIRSETSO/flI**IAK/(AX-l.l)
YOiMIO*CS
1 itl 64.
__ ___ __
A FAPL-I"R-47-12,5
RI
PartI
PSY=PSX*(2.*AK*AMO**2/(AK+,l)-(AK-1. )/(AK+l.))
IF (PA-PSY)35t36 .36
35 PSO=PSX
VO=AHO*CSO
GO TO 9 1
C **OVEREXPAtIDED, SHOCK INSIDE NOZZLE-1
1 36 PSOzPA2
J3=0
~TSOO.*833*TI
CALL PROCOI4(FARTSOCSOAKCPtREXPHISOHSO) I
RHO=CAPSF*PSO/ (AJ*REX*TSO)
VD=1EG/(RHD*Ag)
HSCAL=HI-VD**2/ (2**G*AJ)
DELHS=HSCAL-+ISO
I
* ~IF (ABS (OELHS)-DO.OO5sHSCAL)39 939,38
38 TSO=TSO+DELHS/CP
IF(J-15137t37t3
1
39 AMV=VOICSO
TD=TI
HO=HI
PO=PSO*(TO/TSO)**(AK/ (AK-).))
SOmPIIII-REX*ALOG (PD)
ICON=2
GO TO 101
END
67
-iki
AFAPIA-TR-67-125
IF(ANSI1,2,1
-I DEPuOEPEgIO-AXS
TOLANSaTOLSAMS
2 DEFPc$EPEMD
* - TOLMS'aTOL.
-3 If-ABS(9EP`3-TOLAftSl)%5r4
-~ - 1(21s0.
N CNS
X(2)m0.
9 IF(X(31)9*9912
REIM
C,** FIRST ýGUESS USSING DIR
1(9)-AIN0D#21
10 Y&DIR8AIND
60 7106I
60 TO 6
12 IFIX131-4. 113,13*26
C * LINEAR GUESS
13 X(3J'2o
X16 1-CEPý
XI 7) AIND
I4: 41XS-16)491
- -77~
"AFPJ-TR-67-125
14 IFCX(91-X(71)15t9#15
15 Aw(X(9)-Xf7))/(Xt81-Xf6J)
Y-X (9)-A*X (8)
IF(ABS(10**X (93)-ADS (Y)39s9,6
C** QUADRATIC GUESS
16 X(4)xDEP
17 XM5-AIND
IF(X(71-X (53)18,17,18
17 IF(Xf6)-X(411139913,1
Is IF(X(*6)-X(4))19I3t,9,
23
241
242
19 F(X(9)ilO,11,10OtZ
IF(X(83-X(43J24eZ122
XY"-X(7)
GO TO 3A
iF(91114(7)Z3,
F
243 IF(XCI2X4!)44245~*2
244 Fm(l6-(JO(1J-()
245 GS-CIA
GO TO 270
243 IF(C)245v2482 245
244 Y--hA
.1 GO TO 370
249 G--C/A*I2
246 E-SORT(1-D
27O TO 27
200 29 16,3
IF(DEIX-ABS(X( 31)2,29,
DEPH!NwABSIXfI ))
30 IF(ASS((K)TJ-BSI(X(K)-YY99113371
31 V-TV
K-J-4
32 IF(J-6)33134,34 -
7Z -Z
-F
AFAPL-TR-67-125
part n
KK=K+2
GO 'TO 35
314 JJxJ-2
4 KK=K-2
35 SLOPE=(X{IKK)-XiK))/(X(JJJ-X(JI)
IFISLOPE*X(J)*(XCKJ-Y.))36,36i37
36 Y=yy
37 X(9)=X (7)
X(83.uX(6)
147i xX (5
X16frX(4)
GO TO 6
END
'-I16
-~~~~ -
-cw__~K~
AFAPL-TR-67-i2B
Part if
G1.250 t 7 .0.7150
IJ 1-290
K 1.370
*78.0
,64.0
10.8130
t.0.7500
,1.330
/
,
IiX V 1.400
W 1.570
1.680
Y 1.756
.148.0
-.143.0
1137.0
.13Ge0
10.80c0
10.8500
.0.8800
V0.8500
I 1.500o
-. 1.638
v 1.716
., 1-830
DATA CPR( 7vJ).bIAC( 7-7J),ETAf s)Jl1~
-1145.6
.140.0
1134.5
*-122.0
-10&8250
-90.8700
-sO.870!)
3_.7990
-f
-
-
1
1 ~
J 1.338
~~
1.556
li0.0
2187.8
0.,7000
w.8000
.7 1.52
I~ 1.85D
-:3.
2~.
.0X._7s00
075~
-.-
c
5ýE:ý:
LX
AFAPL-TR-67-125
part 11
72
AFAPL-TR- 67-42
I COMMN ICOM.BPSIt15),DELT(15,15),ETA(15t,15t1NvkiF(1
DATA NoNP/SvB*l11,7*OI
DATA PSI/jo.,ZO.,3O.,40 5O 60 'o*,74(2.I - 0
IA15*1325.,15*1450g,15*155Oh1i
DATA DELTI15*8OO.,159q0.'
DATA ETA/
BO.6OOO.82O,.9OO,~
2
1 5 *9 75 'vi
6 8 5.S
5 *1
065.t,15*1175,152
*8'0*6**
~ 0.960,0.q500 ,O.9600,O.9680vT*G.p
I
97lO~oh9787
CO.8783,O.8996,00;9183,D.9347tD.9899O9615wO
.9,8290943,O.5T3O.9 lOo9694,O*9782Sb.985ov7lO)
DO.8895 0910e
Eo,9OOI0 .O.9200 ,O.9372 59 6 3 70.790950095t
2 OvO
,O-95 9
5 O$Ob
85t1 -
Fo.9o8oo.9212,O.94,4O*9587,o,9718b~o*830v0I98
GO9020281 40060o9700.SCO95i~95i*-
H09330941&9coO180709
H0989.295-480
END
-- 7
AFAPL-TR-67-125
PaA 1l
74
AFAPL-.TR-67-125
Part HI
16
bIO.365 t0.0070 t0.7310 .0.766 .0.0142 .0.8000- 9
.0.05
~e315 10.8930 t2.007 .0.0346 P0.8975
A2*171 t0.0376 10,8999 .2.280 90.0400 10.4000 9
82s372 10.0424 90.8980 t2.482 tO.0454 90.8937 .
CZ.536 0.0474 .008900 ,2*591. 000 089
02.609 t0.0527 t0.8710 /
DATA (CN( 8,J)#DHI 8,JJETA( 8,J)tJ11)
EO.365 9000077 90.7100 f0.602 P0.0125 -.
0.7450 t-
FO.785 10.0163 #0&7680 t0.967 ,0.0198 .0.*8000- i
G1.223 90.0245 90o8380 .1.442 t0.0284 .0.8600 1
113.624 .0.0321 .0.8712 t1.825 .0.0362 .0.8780 i
11.934 ,0.0387 t0.8800 92.062 .0.0418 90#8775 v
J29135 t0.0438 .0.89760 92.190 90*0456 .98722-
K2.244 w0.0479 90.8660 12.281 .0.0502 t0.8600 -,
7-7,
- -- - - -- --- ---- -
II- -..
PCNF
ZF
4 - KAFCý
-zc-
kACC
VAC
T21
P21
T3-
PCRLPF
PCSLOO
VSLO
ElI!rAS
Al 'DITC
TS
J -OTZLP
P55F
T55
'PCOLOU
24 P24-
125
AWPD
ZFA24
WIT
73
-
4 -- i.
"AAPL-TP.-P-125
part n
ETAF
ETAC
ETATHP
ETATLP
AN155
AR125
T6
P6
PS6
£116
TG?
WFAR
rG7
DPAFT
PS8
P59
£169
Y9
V528
PS29
A1629
V29
BYPASS
"HPEXT
liFT
MG?
VA
cvmioz
CVDWOZ
* H
~FGM
FGP
THEEND
S0AOAIN ITITLEa1itlDES3-N0fvP~OIVUNPs1,A1TPm'09
IGASMXZ1#IDCD.O, CwwDZL-tTYSZ~TLALO05
ZFOSz0.633333,PtNFDS1O0O ,PRFDSm2.00,ETAFDSu0.850,VAFDSu1S8l.,
j 2COSO0.314596,?CNCDSS100.,PAtCDS.6.0O.ETACDSUOoS30,VACO$S7S.0 t
T4OSuZ400.,OTCODOSIZ6O., ETASDS=0.w935,DPODSaO.05s
TFHPDS-2I.25,CIHPDSw2.040,ETHPDSQ.Co90,
TFLPDS-53.41,CISLPDSSZ.2B0,ETLPDSUO.900
DPDUDSnO.04tOPAFDSO.*049
AX-0.0,ALTPm0.0'0NPEXT-O.0,-PC8LFZ0.0sPCSLC'in0.05,
PCILOU30.2O,PC8LOSU0.OgMLP-SLU.S0, PSLLPZO.O,
PS55s1 .63 CVDNOZzO.965,CVIWIZsO.985,DELFG=1 .0 ,ELFUNa1.0DELSFCu1.0S
SMOTE DESIGN POINT-
SOATAUI T4=1?5O.,A~aOo0,ALTPe0O000,, TITTLEs1S
SLS7 IDLE --
79
AFP-TR-7-125 1
SLS TAKE-OFF
ýsoxATINT4=24OO..,AP.a].29 ALTPzO0650flITITLEulS
SET-;UP LOW ALTITUDE DASH
SOATAIN-T4=24O~o. AN=1o2,ALT~i'z0500.PI71 TLE=I.
lAFTBNsIT?m3200**ETAAmO.91S
j SUATAIU( T422400.,AMt1.6,ALTPU50000.,XTITLE-1,
IAFTBNk1 ,Tl=1700o*,ETAA=O.94S
SUPERSONIC PARTIAL AIRl
J
IMAMA? T4'=24CD.*A.fwl 4 6:At.TP-5OOOO.pITITLE~i,
IIAFT5JXnhT7z,320G.,-TAAxO.91S
'~
McK~inney, John S.
Ceptain, USAF
S.REPORPT DAYE 7*. YO~rAL NO. OF VAgSC 17 . WO. OF MCPE
November 1967 94 0
da. CONTRACT 0OR GRANT NO. so. onseNiAvobea mapoaT Numear(s)
This report describes a digital computer program titled SIEOTF (Mlmublsin of Tuzbofan
ngine). SMOTE is a computer program for balancing-cyclz turbohan andiu aal of
rinigboth design andi off-design points. The program to written In Fortran TV language
ud was designed for use on an IBM 7090 Digital Compiter. althxou It bas slso been Tun on
an IBM System 360. Performanc~e maps (Block Data format) of the major e*lz components
are required. Information for setting up the Block Data and Input data to lgiven In the report.
Also included in the report is a complete program listing with a descrIption of each stib-
routine and a sormple data pack.
DD FJAONt'e0473 UNCLASSIFID
Sec'ulty Clasefficatioa
UNC LASS! FLD
SecurityClassification ______ ____________
KEY WORDS
____________________________________________
-IKAL4K6LN
POLK WT POLK uT_ POLK WT
I. ORtIGINATING ACTIVTY-: Enter the some and addes a amp=,dby security clmsovelo~in, irring standard statements
of the contractor. subcontractor. grantee, Department of Do-. uh m
fen.. mctivity orother organization rcororat.nawthat) lomn (1) "QuaeUi ed requesteri umy obtain laof this
Ihs report. r-iport.**ou DDC. "'
2&. REPowr swuntry CLAMI7rICATVDN. 'ZAWTm the over. (2) *'Forio announcement and dissemination of this
all security clasedifcmtias of the report. Indicate whether
"*RestrictedData" Is Indal d.I arking is to be iiaccord- r*port by IIDC to not author!zed.
£1112 with app prate secuity regujlations (3) '"U. &. 11owertmer agencies may obtain ctpieu of
2b. GROUP- Automatic downgrading ios ammcfied Is MOD D)I- 'eep rtkolly from DDC~. Other qualified DDC
hisLW
restiver 52D. 10 and Armed Fore Iibat '* Maratal. EnterusrShl ovttrug
the group haimther. Also,. whi.n applicable. , how that optional
nucukings hay. been usend for Croup 3 £rml Group -asauthor- (4) 11U. &.mllltary.agancies may obtain copies of this
izr.d. report directly:16omDDC. Other qualified users
3. REPORT TITLEM: Enter the Complete repot -title in all *hall vequeet 'through
capiltil letters Titles in alt case% should be anclansified.
If a meaningful title casant be selected without classifica--
tion. show tifle claseification In all capitals 1-1 parenhesis (:) "All distribution of this report io controlled. Qual-
inmawdiataiy following the titis. Mind DDC usrns uhall reqouet through
4L DE9CR1PTIVE NOTES: If mept topi. ontor the type of_______________
rep.yrt
prgres.
". inerim
ummry. nnul, rlial.If the sswr has 'bean furnished to the Office of Technical
Giv te atsncusvebo aspeifc eprtngperiod in If.les Daprtet~fCommre formsal. to the publice, indi-
Covered. CL1i etadOWte
S.Ersr
A~l~(S)
te 'amee) f auhors) ti IL SUPPLEM[NTARY IXOTEB: Use for additional explana-
or in the report. Enies loet ame. first name, middlelfnitial. tory notes
If =Auitary, show rank and branch of service. 'The some of
the principgl 4athor imiia absolute tainimum requirement. 1M IPJM!i .L4A ULITARY ACTIVITY: Enter the name of
6. RPOR th daa o th r at s dy,
nteDAT_ then departmental project office or laboratory xponsoring (pay-
6.~~~ ~ ne ~ h ae ~ ~ ~
h eota ~ ~a ~~-n
REOTDT E or) the research and dmevoopmeaLt Include address.
mant. year. or month. yeasm If toe tiuman*o dafte pp-tars 13A TRC:ntra btctgvgabiesdfata
lb.c
on the report, use date of publication,1.ABRC:2ta iig re n ata
sumýe-y of the document indicative of the report, even though
7.m. TOTAL NUMBER OF PArlM The tatal page count It may also appear elsewhere In -dhe body of the technicalI re-
sh~ould follow normal !7aginaflon procedures. Le.. enter the pot. Ufadditional spc is reurd a continuation sheet shall
number 01 Paesm COCAIt ining InfformaItion. be attached.
7b. NU"IrER OF RE*ERENCEft Entr the total turroer of It is highly desirabla that the oabtract of clsassified reports
references cited in the report. be unclecsiftad. Each peragraph of the abetmect shall rend with
S&. CONTRACT ONt O~tANT NtUMBER_ If 4Mprprite, entear an indication of the military security clsassification of the In-
the applicable number of the Contract or Veunt under whifh I ormation in the paragraph, represented as (TB). (S.). (c). or (u.).
the report Was writte%. There Is no limiltation on the length of the abstract. Ho0w-
Sb. ik. fi ad. PROJECT NUMBER: Enter the appropriate Ieer, The suggesetd length is from 150 to 225 words.
millitary deparmen identification. sach %aProject mahr
subpioject nulatez. 575tate mnambo, took tnumber,, ot-o
14. KE WORDS: Key words are tachnically meaningful terms
short phrases that chearcterive a report and may be used as
94. ORPIGANATORIS REPXOrT i1SUBER(S): Eaffl-he 0ff& lafdhX ntuas , for cataloging the report. Key word 8Must be
cial repoit number by which the document Will b14Iet~lifisd selected so diat no security classification Is required. Idrlnti-
and controlled by the originating activity. This muwuber mswt flare, such as equipment model designation, trade name, military
b* unique to this report. pr O J 5 -t cod s ame, geographic location. may be used as key
Si. =1 M"R RVPORT YUtMBER(S,: If the reporL has I words but will be followed byvan indication of tachnical con-
aealisned any oilhsr report numbers (either by thenmrigirater tszt. ",assignment of links. Tulesa, and weights is optional.
o. by the saaonor), aico enter this number(s).
10. AVAflAILIwy/LIunTATIow wa0th. -mc Vater any lim-
itations anmfarther dissemination cif the recrst. -`Mwr thant thoeei
UNC LASSI ID
Security Classification