0% found this document useful (0 votes)
37 views91 pages

AD-825198 Simulation of Turbofan Engine - Part II - User's Manual and Computer Program Listing

Uploaded by

Marek Cel
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
37 views91 pages

AD-825198 Simulation of Turbofan Engine - Part II - User's Manual and Computer Program Listing

Uploaded by

Marek Cel
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 91

UNCLASSIFIED

AD NUMBER
AD825198

NEW LIMITATION CHANGE


TO
Approved for public release, distribution
unlimited

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 ltr, 12 Apr 1972

THIS PAGE IS UNCLASSIFIED


S~ii

"AFAPL-TR-67-125
PART 11

SIMULATION OF TURBOFAN ENGINE


I0 PMT 11. USa's MANUAL A
COMPUTE PROGAM USTING

uJ JOHN S. .cKINNEY. CAPTAIN. USAF

"TIFCHXICAL REPRT .APL-TR-67-125, PART U

1-OVEMBER 1967

T-1= document is subj-ct to special eport controls and each transmittal


"toforeltg go.errments or fcreign nationals rna, be made onh. with prior
approhal of the Air Fjr-e Acro Propulsion iLa;ioratory, wnght.Patter.
mmAFB, Ohio.

AIR FORCE AERO PROPULSION LABORATORY -•

RESEARCH AND TECHNOLOGY DIVISION


AIR FORCE SYSTEM-iS COMMAND
WR!GHT-PATTERSON AIR FORCE BASE, OHIO

-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

Co~ls o otths e rportshold


rturnd xies reurn!a rqr, bysecrit
bliatios,
consdertion, cntzactmll
r nticeon spc~e ocuent
Z:

100 - J=nz-v Z-a6 - OD455 - 19-413


SIMULATION OF TURBOFAN ENGINE
PART It. USEWS MANUAL AND
COMPUTER PROGRAM LISTING

JOHN S. Vc]UVY.EY, CAPTAIN, USAF

"Tnsdomunmet is subjecft to spea m t controls and each trans~ml


to ioeign govenenits or foreign dtionals may be m=de only with prim -
approval of the Air Fore Aero Propulhion Labo=tory, Wright-Patter- "
son AFB, Ohio.j
AFAPL-TR-67-125
Part 11

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.

This technical report has been reviewed and is approved.

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

GU !nitial or guess welues


ITRYS number Gf loops thru engine before quitting
SFC specific fuel consumption
TOLALL tolerance
VA velocity of aircraft
Svj jet velocity

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

SMOTE uses a Namelbt input and!a controlled


be printed for each rur, output, Whereby only selected vaxiables will

Part I of thi's report describes the method


of engine calculations and 1l-e balancing tecnique
and gives some sample results. part _q is intentded
for setting up and running the program, as as a tserle
well as d programmamial and L'cludes instructions
listing. The part-s mry
icdependently of one another.. be_ ted

I1
f

j AFAPL~- I-T7-61-2•
.

PartlP
SECTION I1

.> BLOCK DATA

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.

INPLT RasZl fmelist data and title, rrints title.

ZERO Zeroes c-mmon and certain con-trts.

COINLT Determines ram recovery vnd performs inlet calcuiations.


ATMO62 1922 .RDC Atosphere Tables.

-RAM Ca--Cl zctes ram recovery defined by MIL-E-50(08B Specifications.

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.

4 1 COMNOZ Controls the main nzzme and calculates ERP(6).

PERF Calculates performance after the engine is bala_-ed.


UT Prints output except for controlled output. Prints the main commons in a close
format after each point.

I CONOUT Controls and prints the controlled output variables.


ERROR Controls all pr.lutoLts ifanerroroccurs. Prints name of subroutine where error
occurred and also prints the values oi all variables in the main commons.

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.

. TURB Perf-rms isentrople calculatlozs for tm-bines.

S THERMO Provides thermodynamic conditions using PROCOM.


i PROCOM Calculates thermodynamic gas properties for either air or a fuel-air mi•xtre,
based on JP-4.

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.

"CO.NYlR Perorms wozzle calculations for a convmrgent nozzle.

CONDIV1 1Perbm muzle calcula-io fr aconvergeni-divergent nozzle.

AI'R Gnrlqudai nepoainruie


U! -I

AFAPL-TR-67-125
PartI :

FANDAT Block Data for fan. -


CMPDAT Block Data for compressor. I
CMBDAT Block Data for combustor.
HPTDAT Block Data for high pressure turbine."
LPTDAT Block Data for low pressure turbine. I

"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!

-,.-=• -- _ : ----- • L-• • •"o4• =:


APA -- 2 - V7
-1

Pift H

%II8PTC SMOTE DECKM94/2,XR7


- ~COMMON /POINT/IDATPTj
IDATPT=O
CALL CONOUT(11
CALL ENGBAL
STOP
END
MSAPL-TR-67-125 4)

Pad iH

$SrFTC ENGBh4. DECKrM'94/2rxk(?

IwOrO 91DES ,JOES vKOE-S ,WUOE ,INIT *IDUMP qIAhiTP


21GAS-4X,IDBtA~N,IAFTBN.1VrD ,IMCD ,IOSIDC*IMSWOCtNOZFLTi
31TR~YS ,LOOPERvNH0AP ?4UKMIAPtl4APE0GvTOLALL ,EJRR(6)
CU1OXO?4 /DESIGN/
!PC14FGt),PCNCGUtT4GU *DUM01 ,DU14D2 ,DELFG vOELFH 9DELSFCt
2ZFDS *PCNiFOS*PRFDS jiETA'DSWAFDS ,PRFCF E~TAFCFtWAFCFv
3ZCDS iPCttCOStPRCOS gETACDS,kACDS ,PRCCF ,ETACCFWACCIP
4T40-~ WFBDS ,DTCIJDS, ETABDSvWA3CDS*DPCODSvOTCOCF*CITABCF,
5TFHPDSCUHPDSETHPOStTFI4PCF ,CNHPCFETIIPCF,0*iHPCFT2DS
6TFLPDSC?4LPDS,ETLPOSTFLPCFCHiLPCF.ETLPCFDHLPCFT21OSo
TT2!-OS ,IIFDOS ,DTDUDSETADIJSIAZ3OSDPUDUSDTOI3CFETAOCF,
ST70S ,WFAOS ,OTAFDS, ETAAOSWG6CDSOPAFOSDTAFCFETAACF,
9A55 ,A25 rA6 ,A? ,AS 9A9 PA28 1 A29 f
APS55 vA9155 *CVDNOZ,CVHtN0ZA8SAV ,A9SAV tA28SAVvA29SAV
IT! ,Pi IHI ,SI. *TZ ,P2 iftZ ,S2 I
2T21 9P21 *H21 vS21 ,T3 tP3 ,143 rS3 t
3T4 ,P4 H34 ,S4 IT5 TP5 ,H5 ,S5 I
4755 tP55 tHS5 tS55 E6LF ,8LC wBLDU tBLOS
5CNF ,PRF ,ETAF WIAFC I 1AF ,WA3 ,WG4. ,FAR4 TI
bCN- ,PRC vETAC WKACC ,WAC rlETAB 9DPCOJX ,DUMF v
7CUPP 9ETAT.HP,01TCt-.PtDtTC ?BtH? fWGS IFAR5 .CS t
8-CNLP VETATLP,OikTCLP,OHTF till-'P ,WG55 PFAR55 ,HP]EXT ,
9A14 ,ALTP ,ET.AR ,ZF ,PCNP ,ZC ,PCIC ,IdFB
ATFFHP ,TFFLP ,PCBLF ,PCBLC ,P BLO ,PC8LOSPCBLI .
£?-r'":
COMM~ON / SID3E/f
XXP1 ;XgAF 9XVAC ,XSLF fXBLDU ,XH3 sDUMSI vDUMS2 t
XXT21 vXP;1 rXH2i vXS2i *T23 9P23 91423 ,S23
3T24 VP24 H124 #S24 -PT25 tP25 ,1425 vS25 f
4T28 ,P23 ,112a PS28 ,T29 ~P29 11429 IS29 0
5WA0 ,WED tUC.24 PFAR24 tETAV *DPDUC ,BYPASStDUMFS3 i
6TS28 ,PSZ5 iV28 ,AM28 vTS29 9PS29 sV29 tAP129
XXT55
XXWFR
3T6
4`8
5;G6
,xp55 9'H55 ,XS55
*XWG55 vlrAR55*X1IPD
,P6
'P8
,WFA
tH6
,1a
WbG7
VS6
VSB
vFAR7
,XT25
tXWG24
,T7
,T9
rETAA
iXP25

,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

1Ii8FTC MATRIX D)ECK qM94/2 Xr(7


SUB~ROUTINE MATRIX(E,VtA)
DImENSION E(6,6),V(b) ,A(6.1,PIV(7),T(6,7)
DO 1 1=116
T(197)=A(1)
DO 1 j=1,6
D0 7 1=1,4
-* TEMP=O.
00 2 J=it6
IFfTEM1P.GT.A8%(T(JI))) GO TO 2
TEMP=ABStT(J*JH)
IPIV=J
2 CONTINUE
* - IP1=I+1
DO 3 J.:jPlt7
3 PlV(J)=TlIPIVqJI/Tt!PlVj,
j l!FROI4Zb
IT9=6
4IF(IFP.UM,.Eu.IPIV) 6O TO 6

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

SIBFTC INPUT D!ýCKqK94/29XR7


SUBKOUTINý INP4UT
COMMO4N /PUlt4T/Il;ATPT
CO"~ON / ALL'
1)WORD tiDES tJVES #KVES t4UDE
,JINJT ,1DU1P ,IAMTPj
2IGASRXIOBURNvIAf:Tkb,IV3uO 91MCDISC!S~Y~JFT
31TRYS ,LOOPERINUMAP ,NUMMAP,MAPEVG,T0lLALLER.~f6)
COM419OR IDES IGN/
1PCtNFGU,PCNCGUIT4GU ,DLRIMDI ,DUMO2 ,DELFG ,L)ELFN UE~LSFCt
ZZFflS PCt4FDS,PRFIJS 9ETAFDSWAFIS *PRFCF 9ETAFC.FWAFCF-I
3ZCflS IPC.ICDSPSRCDS ,ETACDSvWACDS rPRCCF ,ETACCFIWACCFI
4T4D WFLDS TC9DS,ETABOS,WA3CUS.!9PCIJDS,I)TCOr-FiFTABCFI
6TFLPOS,C.dLPDS,E-TLPDS.TFLPCF,CNLPCF,ETLPCF,Dh1LPCF,T21D)SI
7T24DS ,WFODD ,DTDUDS, ETADVS, WA23DSVOPDui)S CTDUCF ,ETADCF,
BT70S ,WFADS ,DTAFSiETAADSWG6COSvICfAFDS,DTAFCF,ETAACF,I
9A55 *A25 *A6 oA7 I AR ,Aq IA28 tA29 I
*APS55 ,AMS5 ,CVDNGD!?CVmNUZ,A8SAV jA9SAV ,AZBSAVYA29SAV
C014)13?' / FRONT/
M~ ,Pi ,11 -psi *T2 ,PZ $H2Z ,S2 1
2T21 fpZ1 ,iH21 tS21 vT3 ,P3 ,H3 ,S3 1
3T4 ,P4 H4 ,S4 vT5 ,P5 SH5 IS5
4T55 tP55 .H5, ,S55 vBLF SBLC ieLDJ v8LW I
5CJF ,pRF tETAF W~AFC ,WAF ,WA3 ,WG4 *FAR4
7CNiHP ;-ETATH-P,DHTCHP,DH1TC .SLHP W1G5 iFAR5 #CS I
SCkLP .ETATLP,DHITLLPtOHTF ,BJ.LP ,WG55 IFAR55 ,HPEXTI
9AM .ALTP ,E7-R ,ZF .PCNF InC *pcwc twFS f
ATFFHP ITFFLP ,PCBiF ,PCSLC ,PCBLOU*PCBL03,PC.-.4PvPCBLLI'
COP414GN f SIDE/
XXP1 rXkAF ,.XWAC iXBLF ,XBLDU J ,XF0 VDUMSl tOUS2 I
XXT21 ,XP21 ,XH21 ,XS21 ,T23 ,P23 0-23 2S23 9
3T24 9P24 HI24a ,S24 ,TZ5 tP25 "I25 OS25
4T28 ,P28 ,H28 IS28 ,T24; P29 ,1429 *S29 I
5WAD ,WFD kG24 ,FAR24 IFTAD IDPDIJC ,BYPASS9DtJlqS3I
6TS28 vPS28 vV23 ,hM28 ITS29 ,PS29 9V29 ,AM29c
COMMON / SACK/
XXT55 9XP55 ,XH55 ,XS55 iXT25 *XP25 tXH25. tXS25 5 *
XXWFS ,XI4;,55 ,XFAR55*xwFD ,XwG24 ,XFAR24,XXF1 ,DUX3 I
3T6 ,pb ,6 ,S6 -T7 9P7 ,147 S57
418 ,pa H8l 1SB 9T9 VP9 ,IN9 VS9
54G6 ,wFA .WG7 tFAR7 ,ETAA ,OPAFT ,V55 .V25
6P56 ,V6 ,AM6 ,TS7 ,PS7 ,V7 ,API7 IAlK25
7TSS ,pS8 Iva IAMB ,TSq ,PS9 1 vq *AM9 I
8VA trFRD %vj0 Fr-RD 9VJ I F&K vFOIF r4
9FGM. :=Gp ,b--r .WGT f-ART ,FG ,FN ,SFC
NAMELIST IOATAJNI
AITITLEJDES KItDOE ,INIT ,IWM.MP ,IAMTP PIGASHiX,
BI;BuR" IFTSNIDCD tIPZD ,1.CZLT-oITRYS ,TOLALL,
CZFIJS .?ýC--vFDS9PRFDS ,ETAFDSWAFDS ,PRFCF gETAFCFtlWAFCF
OzCDS ,PCXCDS,PRCODS ,ETACDS,IIACDS ,PRCCF ,ETACCF,eIACCF ,
ET4DS ,IWFSDS, Dt-CODSirETAB!IS,ltA3CDSiDPCOOSOTCEGCFETAZCIFI
FTFPDS.CNPDS,ETH)OS,TFHPCF,CNHPCF,E-ThPCF,DHl*CF,T2!;S
GTFLPDSCNLPDSETLPDSTFLPCFCtfL.PC1%-ýTLPCFO?iLPCFTZIOSI
HT240S ,WFDDS ,OTOUDSETADDStWA23D5,D"DLDSLbTDXCF,ETAD%-FI
lT7DS vWFADS ,DTAFDS,ETAADSWGC-6CE5DPAF9DSDTAKFCE-T4ACFt
J.A5 ,A25 ,A6 ,A7 ,AS A9 ,A28 4A29 2
tT2 ,P2 7T4 ,ZP *PCNF *ZC FPCNC I
LwFB ?TFI-ttP ,TFFLP ,AM ,ALTP wETAR hFPEX7 I
I'PCBLF ,BLF vPvCBLC B~LC ,PCSLDIJSLDU ,lPCBLOBi8LOe
* NPCbLHPBLS4P vPCBLLP,8LLP
P~art n

* 0PS55 tA1455 ,74 ETAV ,WFD I


jPT7 ,ETAA tWFA ?CVDNDZtCVtAN()ZiDELF5. ,)E'LFN IDEL;SFC
C.#*ITITLE=l WILL READ IN TITLE
C IDES =1 FOR~ CALCULATINf, DESIGN POIN,
C,~ MODE =0 FOR CONSTANT 74
r * MOD- =1 FOR CONSTANT PCNC
C '*MODE =2 FOR COI-STANT WI-5
~ 1 C * INIT =1 WILL NOT INITIALIZE PUINT
c* 10UR'P =2 v:I&LL OWJ,,i LCOP INGc wRITE-OUTS AFTER EVERY PUINS
IAi~rP =0 WiLL USE INPUT AM IND M-IL SPEC f-TARP
**v IA14'P =1 WhILL USE INPUT AIM ANID INPUT ETAH
C JAKTP =2 WILL 'JSE t12 AS 71=11+T2 AM4 STANO,6RD P1
C IAMIP =3 WILL USEL 02 AUD STANDAPRD T!
C IGASAX=-l SEPA..ATE FLOW, JINP.i A6
C. lGiASMX=O SEPARATE FLOW, A6=A55
C I~GASMX=I WILL ?Ii!; UCT 4ANt) MAIN STREAMS, A6=A254-A55
C IGASMX=2 WILL MIX DUCT AND -MAIN STREAMS, !kP(T A6
C *$10BURN=1 FOR DUCT BVP.PING, INPUT T2,4
C IDBURN=2 FUR DIRT B UR Nt 4% I NPU;
I WFD
C !AFTBN=; cOR AFTERBURkiNG, I NPUT T7
IRAFTBN=2 LOR AFTERBtiN~iNG INPUT WFA
IID Cu = 1 DUCT HOZZLF WIiLLB C-D
C 1MED =1 M4AIN NOZZLE WIL.L B~E -
C ;#*# NOZFLT1l FO'R FLOUATING MAINH N-OZZLE:
-C *;"UtZPLT=2 FUR FLfATINZ- tDUCT NOZZLE
C NtzHIQFLT=3 FOR FLOATING -MAIN AND DUCT NDZ2?LFS
C *q ITRYS =N NUMBEF OiF PASSES TKRU ENGINE BEFURE WUITTING
DItKENSIUN TITLE(12J
DATA Awi)RD/.f6H !NPUTI/
It0,ATP71I)A7PT-1
9- R!NT '100, .-
i*.AT~rI
CALL ZERO
REAO(5tVATA1NW

IF(I`ITsf.EsU.1) REAW50!'1O1 TITLE


TITi E=0
WR1TFf6t1OZi TITLk
lFihOrr-E.E0.O) v9RI TFi8tl03; 1IDEStA311ALTPiT4 JT241T7
t IFCMUOE.EfJ,1) W'RITEtS01O40 ItDES-AMjALTPrPCNC,T24tT7
IFIMGI)E.-EQ.2) hRISTEC5,1I05) lVESv4AP.ALT~,wFt$ ,T24,T7
CALL CUINLI
RETURN
100 FI)RHiAMIH fDMA P~jIiuTjl3)
101 FURMATtI244ý
102 FtlRMlhTII~il30)'12A6)
i93 FORMATi w*,7ti10t~fS=,3-10X714 Af4=vF7*3v6X7&I ALT04=,F7.O,
-36X7.4 s:-vF8.2,5X7H -l24=.F8.!,5X7H T=F.rls,-ýS
=104 PURi4ATf1H~x,7H I0ES=,13?1OXTM AM~,F7.3t6X7H ALTP=PF7.Ot
W6tI8 PCNC~qF8.3,5X7H T2':-,F8.2s5o^? TT1,F8.-2r6HSSSSS)
105 FOR-MATt1)I0,7H 19%..S.=11311OX74. AfM=,7.3t6X7H MlTP,'!-7.G,
16X?H wF8=,F8#..5,A7H T24-,F8.2,5X7m T7=,F8.2,6H:ý.,iS$SS)

14V
U AF"APL-T21-67-125
Part 11
-

SIbj:TC ZERO DEU~,'h94/2rXH7


SUBROIUTINtE ZEW -V
COAMON I ALL/I
1WOP.O %IDES IJ'JES ,KOES ?l~D ,IKIT PIDUP~P ,1A.%TP!5
21GASMXv1DBURNtIAFTShJ,0CG ,IkCG tbDSH0Ct1MSHOEC~ahtZFLTI
31TRYS PLOOPER,NUMAP xftUhAPqAPAPEDC-1 T~s!ALLfE~k~t6)
COMMON. / F~RONT/ S ,T .P ?2 ,2
ITI fps til s 7 P21q S
3T 24b4*4 7 P H 5
44R5 0P55 t,855 %,555 9BLF YBLC ,iBLDUI tSLIU4 v
5CNF P~RF tET4F ,WAFC tWAF ocA3 vNG4 FA1(4
1
6CNC ,P C tETAC ,WACC WIAC ttTAb 10OPCZ GS 154F
7CNHP ,EITA
THP,ERaTSHP ,04TC ,BLHP ,hG5 FEARS l,;S I
KhtL jETATL9,DHTCLPtD,4TF O~L1P IWGSS ,FAR3ib ,HExT
9AM wALTP ,ETAR ,ZF 1?PCh ZC ,PCNC ,k18 I
ATFFtIP fTFFLP ,PCL t:!FPCBLC P!tilDIJ,PCBLGBtP'C8LiFPYPf-LLP
CUMIG'N / ~ISIDEI/
*xp XWAP
7XZ~ ,XWAC *XBLF #XBtiJv
,XH3 ,JIMIMS ,OUIMS2 I
XXXTZ ~XP21 ,XH21 tXS21 *T23 1 P23 9"23 PS23 I
ýýT24 sP24 .*H24 #S24 IT25 orP25 1,
t25 *S25 ,
4T26. P28s H28 i:s2S '129 IP29
1 9,12-? .S29
5WAD 01FD swG24 ;:AR24 E1TAD ,PDUC tCypASS,Gt;1--S3
6TS28 ,PS28 ,V28 ,bM28 rTS29 .PSzq %V29 ,AM29I
C09MM4D / 5ACJ(/
XXT55 IXP55 Ixft55 *XS55 ,XT?5 rXP25 ,XH25 ,XS25r
XXiWFS ?XWG55 ,XFARS5,XWFO ,XWG24 ,XFAR24,XXPI ,iJUMi
4TS P,P 048 ,19
ISIS VP9 ,H9 PS9
SWC-6 ,WFA vwG? ,E7AP .OPAFT tV55
FAhR7 $V25
6ZS6 IV6 ?StI4 ITS?
,PS7 - ,V7 ,AM7 If-R25
-fse 7S ',Pss ,Ve lapis
sT&9 ,PS9 5 V9 --.
A#'.9
13VA ,FPD tYjD .FGHl 1 vj1I F5GK"a .FGP!) FIPP
9F G. ,FGP ,wFT VGT ,FART ,FGC dFN 04tC
D1'MENSWIN 2j(S3;,LZI4S)tz3i721
ECUIVALENCE tZ1*Ti;*QfXP
JOESO0
JDESO
x~),(Z-!,TX-55;
1
j
16SURN=0
1AFTSU=0

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

!--iPRI-rDS_#_D___IOE_ tMOD 9194-1 vIUP s'AA__

t.S4-t?!Hr-PZFT

':D tPFh :EAF-12FC

9AS5 SA25 #A6


8~IDC .aIC l*
IDSO 2s6SOCvtOZL
I 3TRS ,1 55rCVDN0PCMOt8V
,H$,sA?.&PEOGT0LeLLERRI6

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

ATFF5 tP55 ,115 #PBC


S5 ,LF .8DUPULCBtCLHLJ .8L0
DATAC AMORD6COTC / WC ,A ,LB ,OC$

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).

3 1Ff 1A1TP.Q.~3*OR*!AMTP.&O.4) ETAR*P2iP2T


- - - P2=ETAR'PZT _

IFflAJETPaNE*4? CALL TPEERM0(22tli cltXX5,,OO.0,11


IFCI!-1TP SO.A.) CALL 'THE-RJ4[P2,HZT2tS2,XX!),0,O.0,0)
IFjjf.IT,EO.1) 60 TO 6
IF6!OESA.1O 4 ) GCV TO
PCFG5(1.4~T!CP~SW~(BST,2SPfP
--II PCHFGljcfFt"
GO TO 5
- PCM~zFpCtFOS

* T20SST2
5 ZFXZ$L3S
6 REIVRft
~-77

A FA:'L -TR-67-125

--
1,,I

I' .. • ... .% . + •.. . . .. ' . . .. . . . .

17
>1:

~Lri -

"AP.-TP.-67-125

S.-$FTC ATt062 DECKh94/2,XR7


SUJBROUJTINaE ATHOS IZFTt TM, SIGMA, RHO, THETA, DELTA, CA, AMU, K)
"C THIS IS A SUBRO0UTINE TO CUMPUTECERTAIN ELEMENTS UF THE

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

SIBFTC RAM DECKtX94/2*XRI


SUBRlOUTINE RAM(AMETAR)
IF(At."GTolo) GO T0 2
ETAR1l.
I
1 RETURN
2 IFtAM.GT.5s) GO TO 31
ETAR=1.-0*075*((AM-i.)**1.35)
GO TO1
3O TETAR=i500*/((AH**4)+935..)
END

Al0
AFAPL-TR-.67-125
Part H

sIBFTC GUJESS DE~h419R


FUNCT JUN GU)ESSU(, TTLJP, PDWWD,D,D)D, VO)
IF(M.EO.O) GUESS=VD*((T/TD4)**1.60)*((DD/D,)**O.50)
IF(m.EQ.13 GuESS=VD*((P/-PD)*$1.80)*((DD)/D)**O.33)
IFIM.E(J.2) GUESS=VD*((W/WD)*sO.33)ý((DD/D)**1.00)
IF(PI.E',i.3) GUESS=VD*((W/HD)**O.OO)*((P/PD)**O.5O)
IFtM.EU0.4) GUESS=VD*U(W/1dD)**O.00))*((P/PD)#**.50)
IF(M.EU.5) GUESS=VD*((T/TD)I**1.1O)*((DD/D)**O.6O1
IFWM.EU.6) GUESS=VD*C(.P/PO)**1.OO)*c( (DIDD)3;*O.25)
RETURN
j ~~END ~1

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,

!F I(CNF-CNFS) .GT.0.0005*CNF) MAPEDG=1


IF(lGG0,EQ.1OR.IGO.EO.2) WRITE(8,1000i CNFS,WLH( IGUR
1000 FORMAT(19140 ** CkF OFF I-AP,FIO.4,2XA6?11Wý a *SsssssJ
RAF=WAFC*P2/ THETA
IF(IDES.NE,13 GO TO I
PRFCFw(PkFD'3-1.1/(PRF-1.) PFF,1.,NEAC,:58

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

ZI8FPC COtI.')MP f)ECX,P94/2,XR7


t ~~SWJROULTINE
CUOCMPp
CORW4IN / ALL,
- -WAORD ,IDES ,JUES ,KOES NW'UE ,INIT ,IDlJhP ,JAMTP ,
ŽIGAS!IXZDBIJINgIAFThNjDCD l1ItDC ,IOSHUC 2 II4SHUCtNUZFLT,
31TRYS ,LOOPER,V.-UKAP vNUMAPqkAPEUG,TOLALLvERR(6)

A 1?t~'NFECJtPCNCG1,T4GU tDUMtl'J 9DUkIJ2 s,[ILLFG PUELFU ,thLLSFCC


2ZFDS tPCNdFS,PftFUS sETAFUI-ISAFLJS ,PRI-CF vETAFC)-sWAFcI-
32C0S #PCICDS*PRCDS ,ETACOStWACOS ,PRCCF vkTACCFtWACCf-
4T4DS tiwFBUS ,0TCUOSETABDSIjA3CDSUPCWIDSIDTCL)~vETABCI.,
5TFHPOSECN'HIDS,ETHPDS,TFHPCFCNHPCF,ETHiPCFW..1.PCF,T i)S,
6 2
TRiP0SC.dLP0SETLPDS,TFLPCFCNLPCFETLPCF?0HLPCF-T2lJS,
J 7T2Z'0S ,WFDljS IOTO1JOSETA0OS,IIA23DS,0PDIUDS,0T0UCFETA0CF,
T70S -WFAOS ,i)TAFDSETAAI3S,t6CDS,DPAI-DS,tUTAFCF,ETAACF-
9A55 :A25
:IAPS55 ,A6 ,A? VAS Ihg ?AZ8
,AN'55 PCVONOZCVI4NOZAiSAV ,A9SAV tAZMSAVIA29SAV
COM~MON / FRONT/'
,A29

ITI ,P! thl ,SI 1T2 ,P2 ,H2 ,S2


I2TZI ,P21 v"21 IPS21 ,T3 ,P3 ,H3 ,S3
3T4 vP4 ,114 vS'. *T5 ,P5 :.Is IS5
41'T55 ,P55 ,H55 ,S55 vBIF *BiLC vBLDIJ *BLOJBv
5CWF IPRF vETAF ,4AFC tWAF ,WA3 ,IqG4 ,FAR4t
6CNIC 9PRC tETAC ,WACC ,hAC ,ETAS ,DPCOM- ,CUMF,
7CNHP vETATHPDllTCHPit)HTC ,BLPHP 9WG5 sFARS ICS f
SCNLP 9ETATL?,DHTCL?,DHTF ,BLLP ;WG5z ,FAR55 PHPEXTI
9Ak ,ALTP ,ETAK vZF ,PCWF ,ZC ,PCpiC IWFB
AFFhP ,TFFLP vPCHILF ,PCBiC ,PCBLUU,P%-EL(JbvPCSLHPvPCBLLP
COMJON / COPIP/CNX(151,PRX(15u15),VACXi15,I5),ETAX(lb,,lb,,
- 1NCNPdPT(15)
DATA AORDiLH/HU2) M~ (LOJ) ,6H (HI)/
W!JRDAWORO
THETAuSSWTCT21151as.6681
* CNC=PCNC/(100.*THETA)
IF (ZC.LT.0.) ZC&-O.
IF(ZC.61..1.) ZC=1.

f ~CALL SEARC1I(ZCC.NC*PRC, WACC,ETAC,


IF(HODE.E0.1) GO TO I
* ~IF(I(CNC-CNCS.GT.D.O0005*CNC) HAPEOG=l
IF(IGO.EU.1j.ft.IGO.E0.2) WRITE-ttit1000) CNCSICIH(!C4J3
1000 F!JRMA7(V*Ios* v CNC OFF IMAPvF10.4,2XA6,11Hl * *S~S4S$S
wAC-WACCtrP21/ThETA
1Ff WES.NE.1) GO TO 2
PRCCF=(PRClDS-1. )/PRC-1 -1
ET CC F=E TACz)S/E TA C
WACCF"~ACOS/LgAC
WdRITE(6,100)PRCCFETACCFWACCF,T4-1DS
100 FORM1AT(1IgHOCOP.PRESS~J ESGI6igPCf-1.ds tACFEis

' ' 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

OLL P=PC BLLP*FsLC UT4

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-

- Y:3IThYS ,LOOPtRvhUKWAP UA4AAEGTJLERt6


!PCNFGU,PCN'CGU,T4CU t,01M401 cOUIMfj #EL;ýG tOEE1:, ,OMSe-Ce
2ZF0C ,PCNFDS,PRFOS vETAFLJSiWAPUS tPRKCF ,ETAFCFtUAFCF v
13Zen) ,PCNui)SpkcDS ,ETACDStWACOS 1 PRCCF tETACC',WACCF,
4T40)S INFaDS %DCOfTO~W3ODCMI)CCt!3P
6TFPOSC#ILPFSEThPDSTFiiPCFCtiLPZr,ETHLPCFF4PCF,T2'-0S
- ~ ~ 7 TZ40S ,WFDOS Dit)UUDSIE7AI-DSIWA23i)S,DPD*,DStDTULCF,ETAUCI,~
8T7DS #qF4S TAiii- S VE TA4.-S - 6CjS . PA FOS, 0AFC F,!: IsAAZrF,
9rICAS5 IA25 46 ,A7 .AR ti9 VA28 vA2
IAP-S55 AN5,I) *V50 *CVMx! 2A8A ,A9SAV 2FSv29A
M~ vel1 NJ1 ,si vTZ P2 1'-:2 ,S2 f
z T'2 I P21 01ij s521 ,T3 *P3 2'~3 ,S3

4755 IPS5 ,H5


1"5 2555 PSLF TBLC 08 I-n'U aBL
%, Z
SLNF ,PF.F E~TAF ,WAFCr ,'A w.lA3 tw;4 *=AR4 I
6C1,c oar WACCý , ýAC ,ETAS 9DPC(Pi 9,OUAF I
7CftHP A f'g *PLDiTC -601P ,ber5 tFAR5 zCS
8CNLP iFTATLPiMTq-,LP,[,PTF %8LLP WtG55 ,FAR55 ,4P~xf ,
9AM VALTP ,ETAK rZ ,Lc?,F vzC tpcthC *5'IB
ATFFdIP ITFFLP ,PCBLF #PCSLC tPC8LDUiPC3LC5,PC8LhtPiPCbLLP
DII4ENSItJN fJ(91,DUI4SDI5t153
DATA AWG'RflI6t.CUdCGRBI

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

SEBPTC COPTH OECX,X9412,XR7

IWORO vIDES ~JOES ,iKDES 04W sINJT *IDUHP sIA14TP


4 ~ZIGASMXv1O8URIMJAFT6hslDCD vNCO ,I0SHMCIIISHOCPJZFLTt
31TRYS ,LUUPh~qh6GNAP 1 FhU.MAP,PAPEOGT0LALL9ERR(6)

IPC.NFGLOPCI4CGUT4GU ,OWMD1 ,O'JM0O2 YOELFG tOELFN tOELSFC9


r3ý!CfS vP':C9StPX.ZOS tETACDSgbEACDS tPACCF ,ETACCFiWACCF,
- 5HP0S.:K.PDSETKiPDSTFHPCFC?4HPCFETHPCF ,DIhsCF .T2OS

7T!4nS x~FOUS ,DIEXJDS,ETAOOSI'A230SDPDOODSOTDUCFsETaIJF,


STIMS tWnIAVS ,DTAFQS:ETAADSIWC6COSDPAFOSPt)TAFCFETAACF?
91.5 A25 ,A& A7 as ,A9 .A28 ,A29
I A?355~ ,?~55 ,CV~-OZCVI4ROZABSAV 9A9SAV ,ZSVA9A
J:ýI"W / FK04-
-2 T
2T1' tP21 9H21 9S21 T&3 2P3 "43 vS3 I
314, *P4 ,M4 ,S4 ,T5 ,ps h85 ,55 1
4T55 IP55 #H55 IS55 -BLF ,BLC v8100 tB5113 t

: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=

20ti FU)ABs(j4HtPS-=&Ch$P).T.OFF ACjI0.4*2,A)


?c.&P*S"SI
A CF(,ALL JGTG CALEA6C(1,.KHjTFORTFHNPCtPP1C
- JF~:~AP~T.0RETURNf

?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

CI4NROd ODES IGStPD


vJE INT91./l)~T

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

SJ8TC FRTOSD 0ECKjM94/2iXsR7

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

ATFF1IP ,TPFFP tPCBLF ,PC8LC vPC8LVUvPCBLtJ8rPCBl.HPtPCBUAý


commoUn / SIDEi
XXPI ,XWA!- r-W4C ,X9LF .XBLDU 9XH3 I DIUMS ,DUMS2
XXT21. tXP21 vX'!.21 9XS21 sT23 1,
Pz2 ,H23 ,S23 I
4T28 tP2H IH28 9S28 J129 9P29 1?42-9 IS29 I
4
5WA Wi0 vWG24 9 FAR-24 i ETAV tDPfD'C vBYPASS IDLMS3
6TS2ti vPS28 ,vaa TA14-28 *TS2, 9PS29 IV29 ~AM29
zp1zp
XWAF=WAF
XWAC=WAC
XBLF=BLF
XB L = BL D

XT21=T21
C-1 XP21-P21
.1 XH2l=1121

Ij
I
CALL CUDUCT

ENDi
H

'32,
.1 AFAPIT-TR-G7-125
Paft 11

SIBFTC COOD1CT BECKvM941 2tXR7


SUBROUTINE CODUCT
COMM4ONI ALL!
IWOR.D ,IDES tJDIES ADOES MOD0E 9INIT t!DUJ4P PIAMTP9
21GA:nEIO3URNYlAFTBNtXDCD 91WCD ,IOSHDCI*KSHOCtNOZFLTP
31TRYS ,LOOPER9NOHAP fNUXMAPftMAPE~t*TOLALLv2Rft(6)
COMMON /DES IGN/
1PCNFGUtPCNCGUtT4GU wM'UI~ tDUMD2 ,DELFG 9OELFN ,DELSFCv
ZZFDS ePC14FDStPRFDS tETAFDSIWAFDS PPR-FUi 9ETAFCFw(AFCF v
3ZCDS tPlIICOSiPRCDS ,ETACDSPRACDS ,PRCCF vETACCFIEACCF
4T4DS ,WFBDS ,DTCOOS ,ETABOS ,MA3CDSDPCOOS ,OTCOCIETABCFt
STFHPDSCWIPDSETHPDSTFHiPCFCNHPCFETHPCFDHH~PCFTZDS
6TFLPDSCNLPDSETLPEISTFLPCFCt4LPCFETLPC-FOXL.PC~wT21DSv
7T24DS *WFDDS ,DTUDUSETAODSVA23DSDPOUOSDTDUCFETADF,
8TTDS tWFADS ,DTAFDSETAADS,)IG6CDSDPAFOSDTAFCFETAACF,
9A55 sA25 VA6 ,A7 oA8 ,A9 9A28 ,A29
APS55 9AN55 PCVDNDZtCVM*4OZgA8SAV vA9SAV vA2BSAVtA29SAV
COMMON / SIDE/
iPI tWAF ,KAC ,5LF PBLDU vH3 vDUMS1 vDUKSS2v
2721 vP21 ,H21 ,S21 IM2 PPZ3 *H23 ,,S23 1
3T24 IP24 tH24 vS24 vTZ5 ,PZ5 SH25 9S25
4T25 tP28 ,H28 PS28 tT29 IP29 *H29 IS29
5VIAO RFD vWG24 ,FAR24 ,ETAO tDPDtX ,BYPASSDUfS3
6TS28 iPS2B tV28 9AM28 iTS29 ,PS29 vY29I rAM29
flIPENSION Q(9)
DATA AWORDI ,AWORDZI6HCOOUCT,61PJOZZLI -
IWORDAWORDI
Q(23 uO,

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

1.2051501E-07 I*T24-.2453116E-O3)I*T24-.9433296E-O1 1.724+. 1845537E405


CALL THERJ4(P2-tHA2T24,XXlXX2,~O,OtOO)

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

$IBFTC FASTBK DECKtM94/2,XR7


SUBROUTINE FASTBK
COMMO4N / FRONT/
1TI tp1 ,Hl :S1 vT2 TP2 VH2 lpS2
2T21 PP21 sH21 9S21 9T3 IP3 tH3 ,S3
3-14 ,P4 9H4 VS4 qT5 ,P5 vH S5
4T55 IP55 ,H55 #S55 v81P ,BLC ,BLDU ,BLOBv
5CNF 9PRF ,ETAF vWAFC ,WAF ,WA3 IWG4 ?FAR4I
6CNC tPRC #ETAC tWACC ,WAC ,ETAB ,DPCOI4 ,OUMF
7CNHP ,ECTATHP,DHTCH?,DHTC iBLHP tWG5
8-CNLP ,ETATLPrDfHTCLPvDHTF vBLLP PWG55- ,FAR5 vCS
iFAR55 ,f$PEXTr
I
9A1M ,ALTP ,ETAR vZF ?PCNF ,ZC rPCNC ,WFB
ATFR-IP tTFFLP ,PCBLF tPCI3L'% iPCBLDUtPCBLOBPCBLHPPCBLLP
* COMMON / SIDE/
XXP1 ,XWAF ,XWAC ,XBLP tXBLOU ,XH3 vDUMS1 tDUHS2
XXT21 ,XP21I XM21 ,AXS21IT23 tP23 t,123 iS23 I
3T24 IP24 ,RZ24 IS24 pT25 vP25 ,H25 ,S25 #
4T28~ IP28 ,1128 ,S28 ,F29 IJP29 .1129 vS29 I
5WAD tWFD ,WG24 ,FAR24 ,ETAD ,DPDUC tBYPASSvDtJMS3
6TS28 PPS28 ,V28 9AM26 iTS29 tPS29 tV29 ,A1429
COM4MON / BACK/
XXT55 ,XP55 ,XH55 X)S55 ,XT25 ,XP25 ,XH25 vXS25
XXWFB IXWG55 IXFAR55,PXWFD PXWG24 sXFARZ4tXX)Pl ,DWIB t
3T6 7P 14 S6 ,T7 tD7 v7 tS7
4T8 InB t ,S8 *T9 .P9 I89 TS9
SWG6 ,WFA ,WG7IFAR7 *ETA,- vDPAFT vY55 sV25 I
7"SS8 tPSe ,V8 ,AM8 ?TS9 ePS9 ,V9 YAH9 I
8VA 1,FRO tVJD ,FGMO *V..1Z tFGIM ,FGPD tFGPIM 7
9FGM ,FGP ~WFT -.WGT ,FART ,PG vFN vSFC
XT55=T5q
XP55=P55
XHi55=H55
XS55=S55
XT25=T25
XP25=P25
XH25=H25-
j XS25=S25
XWFB=14FB1
XWG55=WG55L
JXFAR55=FAR55
XWFD=WF1O
XWG24=WG24t
XF AR24=FAR24
-
xxPI*=p1
CALL COMiIX
RETURN
END
4
3.5-
"AAPL-TR-67-125
Pnrt 9

SISFiC COIQX DýC~vK94129XR7


SUIBROUTIME COMIX
COMON I ALL/
IWX 1DOES vJDES tKOES tMD0E ,IN!T wlZDUNP v1&MTP I
,7 2IGASMXvIDb8'R.1RJAFT8NJDCD PINCO tWSi*3CtSMS0CvNOZFLT .
2
31TRYS ,LOOPft~P4h0AP vNM4~APi MAPEDGTOLALL ,ERR (6)
COMON /DES IGN/
IFCNFGU, PCNCG*J, T4GU -,OUXD1 90UN02 tDELFG ,DELFA POELSF-Cv
2ZFOS iPCRFDSvPRFOS tETAFDSKAFOS ,PRFCF -,U.AFCFvWAFCF i
3ZCDS ,PCNCDSPRCDS ,ETACDStWIACVS #PRCCF ,E-TACC-Fti(ACCFv
4T4DS ,jWFBDS ,lTCODS,ETABIDS,3COS,DKCOM,OTCOC:FETAMCF,
5yFnftiS, NHPDS ,EThPDSTFHPCF,CNHPCF, EThiPCF ,DXPCF ,T20S v
6TRFPDSCHLPOSETLPDS;TFLPCFCULPCFPETLPCFDMiLPCFtT210Sv
712405 vWFO0)S ,DTDeJDS.ETAOOSWA23DSDOPUDSDTDUýCF7ETADCFt
BT7DS ,WFAOS ,DTAFDSETAADSWG6COStOPAFO5,OTAI#:CFrETAMCF,
9A55 ,A25 tAb ,A7 *A8 ,Ag 9X28 tA29 f
APS55 wAX55 tCVDNO!,ClMNOZA6SAV vA9SAV ,A28SA~tA29SAV
COIQCbG /' SACK/
1135 vP55 ;?-55 vS55 #T25 vP25 J1I2s 425 1
2VFB PWG55 wFAXSS *WFO ,W62'. ,FAR24 IPI ,D0Uie I
316 ,PtH6 .6 9T7 PP7 ,147 ISI ,

418 ,P8 ,Hs t.58 tT9 -P


IN9 2S9 f
SW(G6 ,bFA ,WG7 tFAR7 YETAA *DPAF'r ,V55 ,V25
6PS6 tV6 *A14& vTS7 ?PSr vV7 vAK? ,AKZS t
7158 ,PSS ,vs tAms I159 ,PS9 VV9 thX~9 ,
4YA ,FRD ,VJD ,FGKiD ,VJA *FGMP ,:GPD PFGPM
9FGH vFGP ,WFT vWGT ,FART VFG ,FN ,SFC
DATA AUrORD/61 COMJX/
OWKENSIG 0QQ9)
VOWDSAWORD
AJ=778.26
LAPSFu21*6.2170
G=32.174049
CALL PRZOCOqI~FAR55,T55,Ex1,XX2,XX'3,xX4,PHI55,XX5)
CALL PROCO'(FAR24tTZ&S,XX1,XX2,XX3,XX4,PH1259XX5)
lIfUDES*EO.O) GO TO 6
C * CALCEXATE A55 ANiD A25 WITH PS2S=PS55
j IF(PSSS4 E-Q.o,) GO To 50
TS55mT55*(PS55/P55) *sO.286
00 1 Iw1,15
CALL PRMOVK(F~AR5,tYS55,CS5S.AX55,CP55,REX55,PiiIS55 7 ,HS55
PHISsPHI15S-REX55*AWOG (P55/PS-r55 1
fF(5S(ELF~l.E*o.O0i*PHtIS) GO TO 3
TS55wTS55*EXPt4.0VDELPHi)
*2 CALI. R.RM
RETURN
so rS552O.875St55
DO 51 lu1,15
CALL PROCOHFAR55,TS55,CS55,AIC55,CP53.REX55,PHIS5,lHS553
V35mA1I55*CS55
ftSCAL-455-V55**Z/(2.*G*AJi)
OELHSPIKSCAL-4fSSS
1F(ABS(C6LHS$.LE.0o0005*14SCAL) GO TO 52
a 51 TS55mTL55+OELfiSCP55
GO TO 2
S2 PS55uP55iEXP( (P~ii55-PHISS5UIREX55)
3 1P0455.GT.HS5S' GO TO 53
val TEIS .101 Pp!5.PSSS,.Tss 7 T"s55,HsHS55
101 FORNATtZZHOSQRT OF H55-Hi$S5 NEG ,6S15.616HISUSSS)
CALL. ERIOR

35
AFAPL-TR-67-125J

RHO=CAPSF*P355/ (AjeEX55*TS55;

I A55c-WG55/ (RHO' V55)


~AN55 =V5~ICS55
'x(IGASt4'.GT.Cl GO TO 54
VRI TE (6,104)A551AIR55
104 FORMAT(20tITUfRB!NE AREA f'ESIGJ,6X614 A55,tE16*8,8H AF.55=*Ei5.8)

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

SUM&T COkr-SH DE'CKdWt4,f2,.XR-


SUBROUTINE CO3AFM

10AtM-S ~os KDES ,l400E 0Uth ,JQMUP PIA14TP


21GASX~tItbWU14?1AF;8NjIDO ,JfCD IJDSHOCtlftH0CtNOZFLT,
31TAYS fU OPERIN0IAP INMJM)APtHAPEOGT0LALLERa(j)

IPCNFCGJPPCNCGUtT46U tDUMD1 *DUN2 *OELFG ,Df'LFN ,OELSFC,


ZZFDS IPCNZDSPRFDS ,ETA90S,WAFOS ,PRF:(F ,ETAl' -FIACF
-,-CBS iPMDCStPRfCDS tETACnUsdiACOS ,PRCCF ,ETACZFvWAr F
AT40S tkFSDS tTC00SiTA8OýýrA3CDS,0PCODS,VýTC0I'FETABCFt
STPDt-2HDgE IDiFP-P2-?Co-HýFDH.4PCFT2DS
- ?724S 9WF*OS ,DTDS'',rA'ýODSV.A23DSPDPDJDSD7G)FETADC'2
$T7DS ,WFAOS ,D)TAFD$,E-TAAOSWC.-KS,DPAFDSIDTAFCFPCTAACF#

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

1 24 WQAT CI*SQfRT(AK7/REX7)::Aj4?/(1.,IX7i1.)wAW1..Z:Ž 3..

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

CALL THERRO(P7,HAT73 XA,4 XQXZ.-(*DO1 01D


FAR7=f(iA-t.6)/1(H8V*ETAA !
1IF(FAR7oGTo.0. GO -1O 6
T'7-T6
Go TO 5
6 WFAX-FAR7*U~G6
IM(AFTNrIEO.!z GO TO 9
ERRW= (WFA--wFAV
it
OIR=SQRT (WFA/WPA.0¶
CALL
GO TO (.7,10,8)tlGO
7 T-7-TiT
GOTO 5
a CALL ER.RO
9 WFA=WFAX
10 FAR.7=(WF6+UA~fli&6
WG7=WIG6iWFA
C MOMCUENfTUM LOSS
CALL PROCO.flF7R-T74(iX1,XX-7rXX3,REX7,PHI7,HT)
AHO7-CAPSF*P7IýAj*REX7wT7 1
V7-WGlI (RJ.'.07*A7)
Q(2) =0.
0!3)=D.
11 ;LG7TbiG7/I V7*A7 i
IIS7cU7-YT$*2/ (2.WGzAJ)
CALL -ttERMtOU.0,HS7,;TS7f~i.P57,XX2,!,FAR?,1)
lF(TS7.Gcuo301.1 GO T-2 110
CALL THC-R1D(I0t10HS7,A0O'.,PHIS7,iXS2,1,FArc7,0)
V7=$QRT2.*'GSAJ*fH7-HS7) I
GO TO 11
110 PS7s-RM~7sAJoitEXT*TS7/CAPSF
PS37. -P565+ (RK1065sV65 **2- RFDV*.z 11ctGCAPS
ps
0 1R SQ9 T fAS5 9 S7 P S7 A
Part 1

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

I 13 CALL THE MO(P7,H7tT?,S7,XX2..1,FAR7,O)


~IF(IOES.Egol) WRITE(6,00) WG6CDS
I 100 FORMAT(19HOAFTERBURNER DESIGNt5X8H WG6COS=,E15.8)
~CALL COM1NOZ
RETURN

END

.21A
PArt ll1

SIBFTC COW4UOZ DECK,1494/2,XR7


COMiMON I ALL/
14OaD #IDES 9JDES AGDES tMWOE ONIT PIDUJ4P PI1ATP
2IGAS#NX:JDBLJR14IAFTBNIDCDJIMCD tlDSHOCvJHSHDCHNOZFLTs

IPCNFGU*PCUCGUT4GU PDUHD1 *DUMOZ ,OELFG vDELFX PDELSFC,


2ZFOS ,PCNFOSPRFOS ,ETAFDSWAFDS ,PRFCF iE-TAFCFWAFCF
3ZCOS sPCNCDSvPRCDS tETACOSPiEACDS ,PRCCF tETACCFYACCF
4T40S aWFeDS ,DTCQODSETAB8DSWA3CDSDPCOOS ,OTCOCFETABCFv
5TFtIPOSCNHPDSETHPOS ,TFHPCF ,CflPCF, ET$4PCFDHJHPCFT.2CSv
6TFLDS0sCNLPDSETLPOSTFLPCF,CRLPCF.ETLPCFDHLPCF--21DSt
7TZ4DS w4FOOS ,DTOUOSETA0DStWA23DSDPDU0S.DTDLCFcEt-AOCF,
8T37S ,WFADS ,DTAFDSETAADSI566COSOPAFOc ;9TAVCF ,ETAACF,
9A55 IA25 ,A6 *A? f~ ,A9 tA28 VA29 91
APS55 ,AM55 ,CVOruOZiCV1INOZtA8SAV ,A9SAV ;A2SSAYA29SAV
/ BACK/
CM".4.O
I.T55 ,P55 IIH55 tS5 5 ?T25 ,P25 1,425 WSZ5
2WFS ~1C5 FAR55 tWc-; ,WG24 i FARZi- ,PI foDtsfl I
3T6 ,P6 ,H6 ,S6 .T7 ,P7 ;117 9S7?
!.TS PP_ INS i.58 tTq 09 IH9 25
5WG6 ,ufFA 9WG7 F~AR7 iETAA zDP-'T ,V55- ,V25
6PS6 SV6 ,AP6 ,1S7 tPS7 fY7 ijp- ,AM25t

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

SIBFTC PERFt DECKv194/2vXR7


SUBROUTIN~E PERF
CW4'OMMON ALIL/
lUOROo *IDES ,jDES, ADES wADDE ONIfT tIOUMP *IAMTP
RlASHXvIO3,JtHJAFTaN21OCD WIMCO vIDSHOCIMSHOCWNOZFLT,
MAY tLOOPER9NOXAP ;wiMHAPtHAPED:,TOLAU-.vF~RR16)
5
lPCtFU,~ tHCGU.T4G0 ,fAJMDl AUD2 9D
.ELFG rD7-LFN v0ELSFCv
2ZFD tPNFD9PRFS
PTAF)StVFDStPRCF ETAFCFtiAPCf
-*IFOS pPCIFCCr~PRC0S ,EAVvADSiRC
ETACCFeiIACCF v
4T05r~fES
*TC00oETBOS:A3C~sOPODSOTCOCFETASCFt

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

4`8 IPS &H38 ,9


WS ,P9 Ing9 159
5I alW*~A I~
VW PF7 tETAA PDPAFT *V55 vV25 t
Wls'. JV6 9A!X6 ITST" PS7 tV7 VAIA7 ,AK25 ,
7758 ZPS8 'vz TAN8 ,TS9 rPS9 *V9 tA149
V.ýA ,i-., IVJ3D fm ,VJm F1G.'M *FGPO vFGPN I
Q ,FCK diFT .1161 *FART PfQ ,FN ,S
t FC
DA.TA AVOaDfl& PERFI

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' %

SI5FTC OU)TPUT DECK,?(94/2vXR7


SuiBROUITINE OUTPUT
CWOI'W4 ALLI

IV=~ ,3DES ,NOOE ,INIT 91OUNP ,1AH'TPv


,JDES AKDES
1 21GASX~IX~iD-UMIAFTBNtIOCO ,1IecO ,IDSHOC,114S'0CNCZFLT,
~31TRYS tLOOPEa~ftOMAP ,NUMRANKAPEOGiTCLALLvERR (6)
lPC*-UI-1PCNCGUtT4GU ,cRJMD1 ,DJM2 *oELrFG *DELFU4 ,lfiLSFC,
2ZFDS *PCNiFOS*P$FOS vETAF9SIVAFDS tPRFCF tETAFCFtWAFCFI
3ZCOS *PCJ4COS#PRCOS *ETACOSPblACOL ,PRCCF ,ETACCF#WACCFI
I
4's405 ,bIFBDS ,OTCOf3SETABD5,WA3CDSDPCODSOTCOý-F,ETABCF,
5T.FHDSPCNHP&SE-THPDS4FHPCFCN4IPCF.ETHPCF,0t*IPCF ,T2OS fi
6'tFLPOSCNLPDSiETLPOSTF.PCFCtMPCFETLPCFDHLtPCFtT21DS,
7TZ49S ttlFOOS ,ODTOJDSETADOS,Rh23DS,DPIDUDS,OTDUCFETt4DCF,
ST70S ,bIFADS ,DTAFDS, ETAAOSAG6COSDPAFE;S DTAFCF ,ETAACI

9A55 ,A25 9A6 fAT ,AS vA9 ,AZ8 ,A29 9


APSS5 ,A%5 9CYDNOZ*CVi4NOZYA8SAV ,A9SAV ,A28SAVtAZ9SAV
COMM4O IFRONiT/ ,4 4 P5 , S
li
IT P1 .tq1 iS1 912 ,P2 tH2 IS2
2T21 ,PZ1 P"21 ,S21 J3 ,P3 tK3 tS3I
4TS5 tP55 qH55 ,S55 ,BLF *BLC ,8LDU B8LOB
5CNF ,PRF PETAF MWAFC tMAF ,WA3 twG4 *-AR4
6CtiC *PRC *ETAC rVACC PWAC ,ETAB ,DP-COMiIXMF
7C~iKp sETATHP,DttTCNP,OHTC t6Lt4P W~GS tFAR.5 PCS
8CNiLP vE'ATL?,D~iTCLPvD4ITF ,BLLP *WG55 tFAR55 ýHPEXT
9AM PALTF S~TAR ,ZF ,PCNF lZC aPChct,wF-6
ATFFtIP tTFFLP tPCBL.F ,PCSLC tPCBLv-WePC8LOB,PCBLHP,PCBLLP
lIP.). ?X~ SIDE/BF ,BLI XH OilS.,OiS
COMMON XWA 9^ iDBFtXLUsH3 tUSwlS
XXT21 ,XP21 9XH21 *XSZI ,TZ3 ,P23 ,H23 ,SZ3 I
3T2&-#;P24 ,H24 fS24 J25 IP25 H825 9S25 i
4T28 ,P2a IN28 tS2a ,T29 'P29 'H29 *S29
SWAB 99FD ,WG24 vFARt24 tETAD .DPDUC B6YP;%SS,DUftS3
6TS28 tPSZ8 ,VZG 9Afl28 iTS29 vPSZ9 ,V29 ,Ar-zq
COIRKON / BACK/
IXXT55 tXPSS 4,455 tXSS5 PXT25 ,XP25 1Xt125 ipXS25
I XI'F8 vXWG55 #XFAR55S.XWFD ,XWG24 SXFAit24,XXPI ,DU-NB
M ,6
9F I16 ,S6 T-t IP7 *H? 9S7?
4TS 'pa 9H8 ,so J9 ,p9 419 ,S9
SW"6 ,WFA ,WGT ,FAR? ,ETAA ,OPAFT ,V55 *V25I
6PS6 tV6 #A-46 ,TS7 tPST 'V7 ,Au7 2AQ5 1
1738 tPS8 ?V8 ,AM$ ,TS9 ,PS9 SVC SAP19 v
OVA ,FRD !VJD tFGHO ,VJ14 tFGMI ,FGPO 9FGPP?'
95GM 9FGP #WFT ,WGT ,FART ,FG IF14 ISM~
DIMERI4SOI WI5v,,rAN1I(80),Ak5Zt80),AJ-S3f48J),A~S4(7Z)
EDUTAL&C
Aa1AftD2/6 UFUTPLJT,624CO )4t*J4/ XIIA-I4,T
D~JIATA ACDvt(A.IS1,P
n TP!T4FGU),1NS2,T1)ttSX1fAS4T'5
DATA (Wf1,I)9!2X4i/fHSUBSQti,&HJC C-Gr6?i KiZLt6HF
DAITA (W(2#I)IzI124)/6i4SHOCK #6HltiSIOE*6H C-a 6!#-4Ok-ZZLEI
IDATA (Nt3,l),Ill1,4)l&HSM4DCK #6HOrJTS1GiHE C-C 9,6tftOZZLE/
DATA t(W4,I )9Islv4)l6ISUSSONv6HIC COtit6HYERF. v6WrOZZLEI
DATA ('4(5,1),1*1t4)/6tiSCNIC ~6W-OMVER,6%rGEAT Xio6fiCZ2LE/

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

SHER'- COe40UT 0ECK~ft4jZXR7I


ZiLSAOUTIME COMOUT(ICON)
CWA1Wk I ALL/
iWORDO ,1CES s.OES ,KOES 190CE lINIT .ili;P tIAMTPv
2IGAS1~vJD8UR)I.IAFT5NvIDCO ,1MCD vIDSHOCIKSHOC,?4OZFLT,

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

t %HHQ25 6r.325 9614T28


~DATA IWOROYL! 3.I190,23O)I
,6HP28 ,6K428 ,68S28 611T-n

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
- --

104 FO~4A,*l )Z EI.516)


ENDfi~tU

16 RETso
7-7

- - AWE
A~~ 74-
7-, ~ -.- ~--- ~ - - ~ - -.-- -
-

lP~F;JP;C~T.jt ICUf4DI ,40002~ #PkLFG tOELFlI PDELSFC.


Zr--S ,PCNFtDS,VIIFOS tETAFOS1 IeAIFS .PFCF #ETAFCFP1AFCF
3ZCDS tPCIKCIS,PRCDS tETACE,.k*ACOS ,PACCF $ETiACCFquAc..CF
4T40S tWFMS ,DTCODSETA9DSWA3COS,DPCICJSDTCOCFETABCF,t
57NDfHESEi Df~ict~:C,=-THPCF9oHHPCF ,TZOS f
712405 tWFDDS WOTDUDStETADDSWA2?,DSPE;P=-iSOIDT.Ft-TA.-Ff,
8T70S 9WFADS DOTAFDS,ErAADS:IiC-V-C5,DPAFDSPDTAFCP)ETAACF,
I9A55 *A25 vA6 YA7 ,A8 9A9 VA23 ,A29
APS55 i AH55 vCYONeZ CVHNHZ,A8SAV ,ASSAV iAZGSAV#A29SAV

I'
-

111 ,P1 :Hl ,Sl PTZ jP2 t312 O52


3T4 tP4.S4 IP ,T5 ,P5 PH5 ,S5
4155 qP55 .1155 tS55 vBLF tSLC v8LDU ,BLOB r
5CflF #PRF tETAF ,WAFC RWAF ,WA3 tWG4 rFAR4p
6CNC tYRC ,ETAC .IIACC tWAC iETAS ,DPCO.4 .OIJHF p
7CtIHP ,ETATiip,DfTCP8P,0H!C 9BLIIP vW;5 rFAR5 ,CS
P
SC ?-LI ETATLPlD.4TCLPPDHTF vBLLP ,UG55 #FAR55 .4?PEXT *~-
9A14 tALTP #ETAQ rF *PCHF 12C #PCHC rWFB ,
ATFFHP ,TFFLP *PCBLF ,PCBLC iPCSLDUPCBLOSPCOLHiPPPCBLLP
XXPI ,XWAF ,XWAC. tX8it )X~LDU oXH3 ,IXJ8S1 ,OWIS2
X,(T21 tXP22 .XH2I .XS21 J123 ,P23 tH23 r525
3124 rPZ4 .1124 #S24 J125 9P25 .1125 tSZ5
4T23 tP28 *828 9S20 PT2-9 fP29 IH29 I529
5VAZ tWFO ,V74? -zFAR24 iETAD tDPNXJ tEYPAS57DUMS3
6TIS2t tPS28 ,V28 ,AXZ6 JT529 #PS29 rV29 rAM29
COK4'0H / BACK/
XXT55 PXP55 ,XU-55 ,XS55 *XT2S ,XP25 ,XH25 mXSk5

0I~4T8lpa 9HS81(0)TRSe (PJT9 S83 11),TR SH4(7

WRIT FRO1 TWO IKHZFPCg tFZCFGPDT4,Kg

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 _- - -

ENDR FTIE SYUCN

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

C CHECK FOR 6 TRAILING DOLLAR SIGNS

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. -

#s.r- - - - - t- ----- -- ,-- 2 - - - - - -

- -

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

-~~ CALL EORORORP)

TSO Xtl 1.FAR 1A)


CALT5'IEPjMO (PD t4 ,O
RE-TURN
END

A
7

-6-M -1 - -N-
FAP-r -r

SIBFTC THERMO 0ECK ,H9412tXR7


SUBROUTINiE 7Ii"ERHtO(PX,HX,TX,SX,AMXLFARK
FX=O.
I
I IF(L.EQol) FX=FAR
Iz:(K-EQ.1) GO TO I
CALL PROCOHCFXiTXiCSAKCPtRtPHIHX)
GO TO 3
1 TX=4.*HX
00 2 1=1,15
CALL PROCOM(FXrTX9CStAKtCPtRtPHI9H)
DELffrHX-H
lF(ABS(OELH)oLE.0.00001*HIX) GO TO 3
2
100
TX=TX+4o*DELHI
WRITE8t!00
FORMAT(31HONf4 CONVERGENCE IN THERMOS$SSSW
3 SX=PHI-R*ALOG(PX')
AMX=1.9863751R
RETURN
END

IL

It

57I
A AiL- TP -ý12f

S18FTC PROM-1' VECKM94/2tXR?


SUBROUTINE PROCOMiPARX,TEXLCSEX,AKEX,CPEX-,REX,Pl4! 04-'E4
IF(FARX.LEO.067623) GO TO1
FARX=0.067623
I IFTEX.GE300) GO TO 2

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

102 =OMAT(1H 1 35KPRC02*L0


INPTE~ f IERtUR BELOWf 30ý2500.,6HSSS$S
103 o222T11HO,3&$IPROCM INPUT T26ERA7)TUR ABOVE 4000.6HSSSS$$)
104 2282OE10)4AT(6.8F9053HRCMIPTFEL-AIR RTEXIO BELOW ZERO,6HSSSSSSJ

3 + 6 . 4-3 --- 0 - -

5 CPX=(CA+F.RX-CPF!-011*+F--X
i~-o Ri 2EtTooas
Asc;.cI

IF WijT ie-PtA, P StWST ý-%JAL -14

-
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

6 BH=?PvP*UjX (IHLIM*I-BX (IH,1) +BX (it*I1)


BL=PPD(BXUI&L,L1Mt.#-SX(ILI )14BX(IL,1)3I
DO 7 J'-ZLJA'4

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

Pk-a I, )8W (BX(1


X4 JH BX( 1HiJLl)~
part 1I

lSi4 CXiIH,,Jt)-PR *CCXtlHtJH)-CX(IKJLI1


DH= DX(IfI,JLi-PR *(DX(I#4,JH)-DX(IHJL))I

IIDL- PR=(BXCILrKL)-BLJ/CBX (IL,.KH)-BX( 1LKL))


CL- CX(!L*KL)-PR *(CX(ILKH)-CXI*LKL))
OX(IL,KL)-PR 4(OXcIL9KH)-DX(lLt;Zt)I

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

SIBFTC MAPBAC DECK,1494/2tXR7


SUBROUTINE MAPBAC(MAPMAPGOTFFSTFFCNSCNPCNTNODE, !GONUM)
DATA WHWLtWTWS/6H H*P. ,6H L.P. ,6H TFF ,6HSPEED /
1 FORMAT(1HOPA6912HTURBINE MAP tA6,4HWAS-vel3o6tlOH AND tdOWtE13*6,
16H$$S$$$)
2 FORMAT(1HOrA6tA6922HWAS ALSO CHANGED FROM ,E13.695H TO tE13o6t
16HS$$S$S)
IF(NUM.GT.0) GO TO 3
NUMH=O
NUML=O
3 IGO=MAPGO+?*(MAP-1)
GO TO (100r200,300,400,500t600),IGO
C *** HIGH1 PRESSURE TURBINE
100 TFF=TFF+O.1*(TFF-TFFS)
WRITE(891 )WHWTtTFFSTFF
RETURN
200 CN=CNO.O5*(CN-CNS)
IF (MODE*NE.1) PCN-PCN*(CN/CNS)
IF(MODE.EQ.1J T =T *(CNS/CN)**2
WRITE (891)WH#WStCNStCN
IF(HUJ4H*GT*2) GO TO 210
14UM1l
NUMH=NUMH+1
RETURN
210 DELCN=CN-CNS
IF(DELCN.GE.O.) RETURN
TFF=TFF*( 1.+DELCNICN)
WRITE (8,2) WHWTtTFFStTFF
RETURN
300 TIFFTFF+Ool*(TFF-TFFS)
WRITE(8,1 )WHWTTFFSTFF
GO TO 200
C *** LOW PRESSURE TURBINE
4-30 TFF=TFF+O.1*(TFF-TFFS)
WRITE(8,1 )WLWTTFFSTFF
RETURN
500 CN=CN+0.05*(CN-CNS)
PCN=PCN*( CN/CNS)
WkITE(891 )WLWSCNSCN
IF(NUML.GT.2) GO TO 510
NUM1l
NUML =NUML+ 1
RETURN
510 DELCN=CN-CNS
IF(DELCN.GE.Oo) RETURN
TFF=TFF*( 1.+DELCN/CN)
WRITE (8,2)WLWTTFFS,TFF
RETURN
600 TFF=TFF+Ool*(TFF-TFFS)
WRITE(8,1)WLWTTFFStTFF
GO TO 400
END

61
AFAPL-TR-67-125
Part 1U

SIBFTC CONVRG OECK,N94/29XRT


SUBROUTINE CONVRG(TI ,HI PIt.51 FARWGPA .JOEStADPRv
IT09HO tPDSOtTSO#PSO9VOAMo, ICON)
C ICONSI SUBSONIC, COMPARE PI WITH PR
C ICON=2 SONIC, COMPARE P1 WITH PR
C ICON.4 ERROR
AJs-77S26
CAPSFs21 16.217
Go32.*114049
CALL PROCOM(FARTIXX1,XX2tXX3tXX4,PH1IXX6)
C **SONIC CALCULATIONS

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

C **NON-DESIGN, CALCULATE CRITICAL CGNDZTIONS

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

C **NON-DESIGN, CRITICAL AND SUPERCRITICAL CONDITIONS i


13 PSO=PSO*AOCRIT/AO
PR=PR*AOCRI1T/AO
ICON=2 -
GO TO 10

G **NON-DESIGN, SUBSONIC CALCULATIONS

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!

4- 'BFTC C90WIV OECKPM94/2tXR7


.? - SOUSOUTINE CoeONIVITII4I,PlSIFARMGPAIOESArAOPIft,

C ICW6-3 S0NICs SHOC OUTSIDE NOZZLE, COMPARE PIP WITH PI


C CON&4 EWRR
3~ DIMENSION 0191
Q(23m0.

-, - - 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

C ~*IDEAL EXPANSION i)ESISN, CALCULATE AO


I PSO&PA
j TSOzTI*(P$D/PI J**o2S6
6 JzJ~1
I ~~CALL PRDCOK(FAR*7O,5 ,A0C9EXP1SHO
PHICALmPHIsI-REX*ALOG (PI/PSOl)
6EA tI-PHICAL-PHISO
IFCASS(OELPHI ;-O.0001.SPHICALIS,8,7
7 7SO-TSO#EXP(i4.sDELPHI)I
-IF(J-ol536t6t3

-~
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)

fJSUBSONC pFLPOM iJETO

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.

C SUPER.CRITICAL FLO`Wv ISENTROPIC EXPANSION TO PA

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

21C SUPERCRITICAL FLOWt ISENTROPIC EXPANSION 70 AO

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

-Cr~* -OVEREXPAISOEV# FIND SHOIC% POSITI'ON


34 PSXcPR*R(TSOITIJ**(AKlgAK-l.;)

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

C *** OVEREXPANDEDI, SHOCK OUTSIDE NOZZLE 4

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

S1SFTC AFOUIR DECKM94/2tXft7


SUSK.OITINE AFQLIIR(XAINDDEPENO,AI4SAJTOLDIRANiEW#ICON)
DIMNEISION X(9)
X(1)lxMNE-OF ARRAY TO USE
t AIND'INOEPENDANT VARIABLE
:C DEPEND 6EPENDANT VARIABLE
C A*SwANSKWER UPION WHICH TO CONVERGE
C Al-MAX NUMIBER OF TRYS
Z ftO ý7ACESfT OLERANCE FOR C-OIVERGIENCE
:C OIRaOIREC71ONi A%9 PERCENTAGE FOR FIRST GUESS
C ANEWuCACULATED, VALUE OF hlEXT TRY AT INDEPEMDANT VARIABLE
C -ICOS~aoWTROL -1 GO TMRU LOOP AGAIk
f: -2 YOU HAVE REACHED THE ANSWER
ýc *3 COUNTER HAS HIT LIMITS
-CX12)-COWMTER STORAGE
C.X(3iwCHOOSES METHOD OF CO3NVERGENCE
DEEDVAR
AC X~iSECOND IRQ VAAR
C -X7I-)SECCWC IND VAR
C XID)=FIRST DEMON VAR
C XlffsIKST IND VAR FRTETYT OTN
C XM3 MUST BE ZERO LNFRTETYT OTN

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

sIBFTC FANDAT DECKt,494/29XR7


BLOCK VATA
COMMON / =AN/CNU15),PR(IS,15),WAC(15915)IETA(15el5)i;Ni-flP(1T
DATA N9NPI1Ot6i3*7i5*lOt8v5*0/ 5
DATA CN/0.3.0.4,0.5.Oo6.0.7,0.890.9,l.0,1.1,l.2, '04-
DATA fPR( l,J),WACI 1,J).E-A( 1.J),J1,v 61/
A 1.000 *71.0 ,0.7300 1 1.030 , 69.0 v0.7350
B 1.070 60.0 ,0.7400 , 1.096 2 50.0 t.0.300 1
C 1.112 *40.0 .0.70O0 t 1.120 , 26.0 C).6-200 /
DATA (PR( 2,JItWACC 2*J)vETA! 21J)vjalv 7)/
D 1.000 ,86.0 t0.7300 1 1.050 1 81.0 40.7500 1

G1.250 t 7 .0.7150

IJ 1-290
K 1.370
*78.0
,64.0
10.8130
t.0.7500
,1.330

/
,

DATA (PRf 4sJ)2WAC( I-vJ)iETA( 4.J),J1v -r)/


72.0 v 0
.
8 0 0 0
t

L 1.000 .1115.0 g,0.7200 1 1.092 9112.8 t0.75D0


M 1.220 il07.5 16.8000 r ".310 11132.0 .0.*8250 -

-N 1,400 ,94.0 -90.8500 -91.474 1 63.0 .0.8000


0 1.488 f.80.0 .0.7800 /
PATA EPR( 5,J3-,WAC-i 5,J)tETA( 5,J)#J=1,10G)/
p 1.000 -9132.0 10.7000 1.160 v131-0 10.7500 -,

01.296 v12G.5 10.8000 1 1.B70 .9126.2 -10*8250 1J


R~ 1.460 .122.0 .0.8500 1 1.524 7118.0 10.8.700 v
VT 1-610 .11.5 .0.-500 1.668 *102.8 -.0.8000 /
DATA A(PR( 6,j),WACt 6,J),ETA( 6tJ).*J=1%1O)I
(i1.000 .150.0 ,0.6650 ,1.25 19. .070

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
-
-

2 1.000 .170.0 10.620C 1.2190 .170.0 -*0.70mD0 .


A 1.380 v!70.4) rO.7500 .1.548 4-63'0 10.11000 if
B-1.650 1166.0 V.0.8300 -9 1.124 -,1mo -1.50
C 1-825 1158,6 .008700 -, 1.fi04 312..8 -100;S5D0 -3

D 1.978 .145.0 Oaasai - 2.000 .1I42.2 _0117-900.


DATA (PR( B.JIWACC 8q.)),ETAf 8,J),J=,l,;0)/
E1.000 1190-0 .7,05sco 1.110 .9150%0 -iD_620=

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

L 2.000 -,21).0 U~t.77.V50 -. 2.129 12D9.5 -. 0.7fsc*


5 0
H 2.25e *Zt;S.) -. 0.77750 20~34 -pZD - -0SQ .

2 .42LIF z209,0 ýA.TZZD .,~40 ,9~~~738 i -

v),3 -Cý q 1.350 &_22-5.O -.O.7-50 - - -r.

4j ;t Z.4-20 Z.24.5 ,(k..7043- -r .2.65D -0600 ~


3f1' -

5ýE:ý:

LX
AFAPL-TR-67-125
part 11

SIBFTC CMPOAT DECKtN94/2#XR7


BLOCK DATA
CORMON / CONP.'CN(15).PR.(15,15).bIAC(15915),ETA(15,15)tNNP(15)
DATA MNSP/10t2*6t2*8.4*10,2*8,5*0/
DAfA CN/0.5,0.6r0.7,O.611 0,89,pOe92,0.95.1.O ,1.l,1.15,540./
DATA (PR( 1.J),WAC( 1,J)tETA( 19JtJJu1 61/
A 1.000 , 21.4 90.5700 t 1.599 9 21.1 90.6000 t
5 2.018 9 20.8 .0.6200 9 2.335 1 20.5 90e6300 1
C 3.024 9 19.6 .0.6500 , 3.533 * 18.7 v0.6200 /
DATA (PR( 29J),WACI 2,J3,ETA( 29J)vJ~lt 61/
D 1.0o60 t 24.9 .0.5100 t 1.689 t 24.9 t,0.-201i ,
E 2.497 9 24.7 9O.670u t 3.096 t 24.3 90.7000 t
F 3.63S t 23.8 90o7700 9 4.234 v 23.1 .0.7000 /
DATA (PR( 3..J),WAC( 39J3,ETA( 3tJhvJo1. SI/
G 1*000 , 29o4 90.5650 9 1.599 , 29.4 .0*6200 9
H 2.198 9 29.4 g0.6600 , 2.719 t 29.2 90*7000 t
1 3.515 9 26.9 90.7500 , 4.174 9 28.5 t0.7650 t
J 4.796 9 27.9 .0.7500 9 4.946 9 27.6 90.7430 /
DATA (PR( 4%JitWAC( 49JJETA( 4,J),J1,t 83/
K 1.000 9 35,6 90.5600 , 1.719 9 35.6 ,0.6200 t
L 2o737 v 35.6 .0.7000 t 3.575 , 35.6 90.7500 ,
N 4.353 v B5.3 90.7800 9 5.012 9 35.0 .0.8000 ,
N 5.671 9 34*3 90.7800 9 5o970 t 33.8 90*7650 /
DATA (PR( 59J~aWACI 59J)tETA( 5@J)..J-1,10)/
0 1.000 9 42,7 *0.5520 9 2,084 t 42.7 ,0.6200 ,
P 3.132 9 42.7 .0.7000 9 4.054 t 42.7 90.,1500 9
0 4.922 9 42.5 .0.7800 t 5.431' 9 42.4 .0.8100 t
R 6.000 9 42.0 *0.8300 # 6e407 't41o8 .P0.8100 ,
S 6.988 9 41o2 .0.7800 t 7.138 9 40.9 9C~o7780 /
DATA (PR( 6.JIWAC( 6,J),ETA( 6tJ),J*1910)/
T 1.000 v 45.4 .0.5500 1 2.323 , 45.4 90.6200 9
Ui 3.353 9 45.4 .0.7000 ,i 4o293 9 45.4 90*7500 t
V 5o240 9 45.2 90.7800 r 51731 * 45.0 .0.8100 t
W 6o389 44.8
, .0.8300 't 6*808 v 44.5 v9*8100 t
X 7.377 9 43o9 #0&7800 # 7.*587 p 43.7 T0.7760 I
0AYA (PR( 79JJ.WACI 7tJ)#ETA(ý 79JI.Jm191OI/
Z 1.000 9 48.1 .0.5400 v 2.617 r 48o", t0.6200 ,
A 3o683 9 48.1 .0.7000 9 4.653 t 48.1 ,0.7500 .
B 50611 9 48.1 .0.7800 , 6.222 t 48.0 .0.6100 9
C 6o868 * 47.7 .0.8200 v 7.168 t 47.,6 .0.8100 ,
D 7.796 t 47.3 00.ldOO v 8.133 , 47,.2 .0.o7700 /
DATA (PR( 8,JlNAC( 8,J)PETAI 8,J)PJ=l,101/
E 1,000 , 51.6 ,0.5200 9 1.491, , 51,6 tOO550O ,
F 3n036 , 51.6 .0.6200 t 4.192 , 51.&6 ,0.7000 0
G 5.192 9 !~1*6 .0.7500 9 6.263ý 1 51.,6 .0.7800 ,
H 6.689 * 51.5 .0.7950 v 70347 , 51.3 ,0.8100o ,
1 3e078 9 51.1 .0.7800 9 8.725 t 50.7 .0.7500 I
DATA IPR( 9,JIWAC( 99J).ETA( 9,J).PJ-1, 8)/
J 1.000 t 58.7 ,0*4600 , 3.395 9 58s7 9095800 ,
K 5.731 p 58.7 ,0.7000 # 7.587 , 558*7 ,0.7500 I
L 8.186 s 58.7 .0.7550 9 0.784 q 5306 .0.7500 ,
M 9o683 9 58.5 .0.7000 ,10.042 * 58.4 ,0.6750 /
DATA (Pft(10,,'),WAC(1OJ),ETA(10,J),Ju1, 8)/
N 1*000 v 61.4 10.4500 , 2.976 v 6,ý*4 .0.5500 It
0 4.916 *61.4 .0.6200 , 6e838 9 61.4 .00.7000 ,
P 7*SS7 *61.4 .0.71250 f 8*485i r 61o4 .0.1250 ,
Q 9.'365 *61.4 .0.7000 .10.581 1,61.4 ,0.6200 I
END

72
AFAPL-TR- 67-42

$I8FT C14BOAT DECK,14941-r-XR7


BLOCK DATA 5 )

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

S18F1'C HPTDAT OECKPN94/29XR7


BLOCK DATA
COMMON / HTIJRB/TFF(15),CN(15.15),DH(15.15).ETA(15915) ,NNP(15)
DATA N9NP/10,9*15.12t5*0/
DATA TFF/ 16.86 1.18&27 920o17 t20e66 920*90 t
121.08 t21.25 t21.43 921o64 t21.92 .5*0.0
DATA (CN( 1,J)9DH( 1.J),ETA( 1,J)tJ1,1l51/
AO.191 90.0029 10.6150 90.344 9000051 .0.7000 9
80.526 .0.0076 90.7780 90.727 .0.0097 ,008000 t
CO.957 .0.0120 .0.8000 91.167 t0.0137 .0.7875 9
01.340 .0.0148 90a7692 .1.569 .0.0157 t0.7340 .
E1.,761 .0.0161 .0.7000 91.971 90.0158 .0.7550 ,
FZ.193 .0.0150 q0o6000 .2.454 vfoe0 3 0 .0.5250 ,
G2e641 .0.0108 t0.4720 .2.842 90.0074 .0.4000 ,
"H30005 .0.0031 90.3000 /
DATA (CNI 2,J),DH( 2#J)tETA( 2,J),J*1,15)/
104191 90*0034 .0.6000 #0.402 90*0072 .0.7000 t
JO*593 90.0102 .0.8000 .0.727 .0.0122 90.8200 t
KO*861 90.0140 *0.8270 .1.000 .0.0158 .0.8300 ,
11.129 .0.0173 .0.8275 ?1.301 .0.0191 9,04210 t
P01.474 90*0205 90.8162 .1.741 90 .0223 .0.8000 t
R2e009 .0.0234 .0.7610 92.316 .0.0235 .007000 9
02.751 #0.0217 90.6000 .3.158 90o0169 .0.5000 9
P3.445 .0.0115 90.4150 /
DATA (CN( 3,JJDH( 3,J),ETA( 3..J),Jn1,15)/
Q00.91 .0.0041 90*5700 90.445 .0.0090 ,097000 .
RO.670 .0.0130 .0.8000 .0.890 90 U166 90o8400 1
S10091 .0.0194 .0.8448 .1l.263 r0 .0216 90.8420 ,
T1.493 .0e0241 ,0.8400 .1.722 40.0263 ,0.8315 9
UJ2,009 90.0234 90.8170 92s2583 ,ooC298 .0.8000 F
VZ.603 ,0*0310 .0.7495 .2.861 .0.0311 90*7000 9
W3*100 90.0306 90*6578 93.330 90.0292 90o6000 9
X3*445 .0.0281' #0.5800 /
DATA (CNC 4#J)9DH( 4,JJ.ETA( 4,Jl9J.1,15)/
yo.191 ,0*0047 s0.5580J 90.260 .0.0061 90.6000 9
ZO*488 .060108 .0.7000 90.708 .0.0148 .0.8000 9
AO.933 90.0184 .0.8400 .1.167 .0.0220 90.8500 t
31.416 10*0252 ,0.8500 ,l@593 .0.0274 .0.8480 9
C1.837 ,0.0302 t0.8440 .2.019 90a0320 ,0.8400 t
D2.325 .0.0349 90.8270 .2.564 v0s0371 90.8170 9
EZ.890 .0.0397 .0.8000 .3.205 .0.0'S5 90.7710 ,
F3.445 .0.0445 90.7500 I
DATA (CN( 59J),DH( 5,JJETA( 59,,tJJ-1,1)/
GO.191 .0.0050 .0.5500 90*306 ,0.0079 90.6000 9
110.536 90.0130 .0.7000 ,0.765 ,0*0173 .0.8000 t
10099:; ,0.0212 .0.8400 .1.301 .0.0259 .0.8600 ,
J1.512 .0.0289 .0.8599 .1.799 90.0324 ,0*8566 .
K2,086 v0.0360 .0.8519 ,2.383 90.0400 90.8460 9
L2.698 .0.0446 r0684 2 6 .2.928 .0.0486 .0*8400 9
M3*138 90.0536 .0.8400 .3.215 .0.0576 .0.8437 ,
N3*225 .0.0595 P,08475 /
DATA (CN( 6tJ)tDHt( 6tJ)tETA( 6..~),J1,15)/
00.191 .0.0061 .0.5250 .0a364 9000108 .0.6000 I
P0.632 .0.0173 90.7000 .0.880 90.0227 g,0'000 9
Q1-115 ,0.0270 .0.8400 .1.311 .0.0306 04.8600 ,
R1.531 .0.0346 90.8721 l1.722 t0.0379 .0.8800 9
S1,952 .0.0425 .0.8840 .2.167 90*0'72 P068870 ,
T2.316 90*0508 .0.8875 .2.471 .0.0551 .0.8876 t
U2*545 #0.0576 ,0.8868 ,2o588 90a0601 ,0*8838 f
V2*607 .0.0628 .0.8798 /
DATA (CM( ?,JI.D4If 7,J),ETA( .7,J),Ju19151/

74
AFAPL-.TR-67-125
Part HI

I'IOM9 .0.0072 PO.5005 .0.440 *0*0148 90.6000


XOe698 .0.0212 .0.7000 t0.976 P0.0277 '00SO0O0
Y1.225 P0.0335 t0.8400 .1.411 .0.0374 i0.8600
Z1.541 t0.0403 9
.0.8700 .1.651 .0.0428' .0.6800
A1.780 .0.0459 ,0.8855 .1.899 90.0490 90*8900
81.995 .0.0518 O.0.691 P2.040 .0.0540 9008900
C2,086 90.0562 90.88C0 92.1.24 10.0594 t0.8826-
02*11'3 90-0630 .0.8695 /
DATA (CN( BtJ),D*4(
8,J3,ETA( 8tjtjJ=1,153/
EO.191 .0.0079 .0.4995 ,0.493 #0*0176 .0.600'0
FO*746 .0.0245 1067000 ,0.899 .0.0284 .0.7580
GI*043 #0.0320 .0.8000 ,1*167 .0.0353 10.82600
H1.306 ,0*0389 .0.8400 P1.397 .0.0414 .0.8500F
11*493 90*0439 .0.8600
J1.708 ,0.0504 ,0.8750 s1.627 .0.0475 .0.8710-
;1*780 .0.0536 .008750 ,-
K1.837 v0.0576 ,008690 ,1.852 .0.0598 .0.8600
11.,856 tOtO624 90.8495 /
DATA (CN( 9tJJDH( 9,J1,ETA( 9,J)#J=1.15)i
MO.191 .0.0084 90.4855 .0.344 .0.0143 .0.5320-
I00*823NO.545
P1.033
Q01.244
.0.0209
v0.0297
rO0.360
90*0432
90,6000
.0.7000
,0.7690
,0.8100
i0.689
,0e938
tle148
.1.340
0.0O256
90*0331
.0,4396
90*6500
.0.7380
.0.'8000 p
.0.0472 to0i8210
R1.397 t0.0500 .0s8255 ,1*435 #040,522- .0.6
g270-
S1.474 .0.0551 .0.8230 ,1.493 .0.0576 .0.68150 p
T1.497 t0.0601 90.8000 /
DATA (CNI10,J)',DH(2OJ),ETA(j~.J),J.1 ,1 2)/
130.191 .0.0119 ,0.4210 .0.287 .0.,0162- 40.4495
VO.388 .0.0205 ,0.5-000 90*478 .0.0241 iO,&5300--
Ii0.574 P0.0283 .0.5620 ,061 .. 31 057
XO*703 t0.0342 .0.6000 .0.751 ,0031. ,06110
Y00 75 .0.0396 ,0.6170 ,0*1823 .0.0428. .0.240
ZOe842 .0.0454 ,0s6195 ,0.847 .0.047? 90o6050'
ENDj /
AFAPL-TR-67-125
Part 11

SISFTC LPTDAT DECK,9194/2#XR7


BLOCK DATA
COMMON / LTLMRBITFF(151.CN(15,15i.DHI1S153 ,lETA(15,15),N.NP(15)
DATA NtNP/11,9915912,994*0/
DATA TFFI
131.50 .36&60 .o41.60 p4tv0.%0 #50.22 951089
253.41 .54o60 .55.69 .56.89 .50.10 04-0.0
DATA (CM( 1,Jlt0N( 1,t3.gETA( 1,J),Jmlk151/
AO.365 .0.0029 .0.7120 #0.529 .0.0042 .0.7300 0
80.730 f000056 .0.7472 90.967 1p0.0070 .0.7300 9
C1.204 .000081 P007140 v1.405 .000089 .0.7000 f
01.606 ,0.0045 .0.6850 .1.752 90.0098 t.0.6730 9
E2.007 00.0100 .0.6452 92*2ii3 .0.0098 90.6200 ,
F2.445 90.0092 .0.6000 p2o591 t000085 g0.5750 .
G2.792 .0.0070 .0.5310 .2o920 ,0.0056 .005000 f
H3*285 .0.00001 .0.3850 /
DATA ICU( 2.J).DHI 2.JltETA( 2.JlJml,159/
10.365 lo0.o-'42 .008000 90.547 t0o0063 .008100 .
JO.785 .0*0086 .0.8200 91.058 9000111 t0.8300 I
1(1.277 .0.0128 .0*8330 .1.432 O.900139 .0.8290 t
L1,679 .0.0153 .0.8100 t.1.79 t0.0162 .0.8000 9
912.04'. .0.0167 V007850 .2.208 ,0o0171 .0.7.600 .
P42.354 .0.0173 .0.7450 .2*600 ,0*0170 ,0.7000 f
02ý847 .0.0162 .0.680f) 03.029 .0.0151 .0.6450 t
P3,285 .0.0123 ,0.59c0
DATA (CNI 3tJ)tD~4( 3*JltETA( 3,.J),J-1,15)/
00.365 .0.0050 .008000 .0&586 .0.0061 .0.8300 9
R0.858 .0.0114 .0.8600 .1.067 .0.0199 ,008630 t
SI.241 90.0159 g0.8670 #1&423 .00.0178 .0.8700 ,
Tl.606 .0.0195 .0.8720 .1.825 #0.0215 .0.8720 ,
U2,007 90.0229 .008700 t2.226 t0.0243 v0.8670 9
V2*409 90o0252 90*8600 .2.513 vO.0259 .0.8500 t
W2.755 #0.0265 .0.8300 .3.029 .0.0268 90.8000 .
X3*285 .0.0263 .0.7600 /
DATA iCN( 49J)PDH( 4tJ3,ETA( 4#J9.Jultl5)/
YO*365 t.0.003 .0.7995 9C*420 90.0061 .0.8000 9
Z0.675 .0.0098 .0.8400 90.876 .0.0125 .0.8600 9
AI.m095 .0.0153 .0.86830 *1-277 .0.0176 .0.8*730 f
81.496 90.0201 .0.8800 .1.733 .0.0226 .0.8830 9
C2,,025 .0.0254 .0.8835 .2.190 ,0.0266 .0.8830 9
02.354 90.0279 p0.8800 92.582 .0.0293 .0.8740 t
E2.847 #0.0305 V308300 ,3o102 90.0312 #0-8350
F3.285 ,0o0315 .0.8200 /
DATA ICN( 59J)tDH( 5tJ),ETA( 5,J),tfrl,15)/
GO*365 90*0057 90.7750 v0.502 .000078 .0.08000 v
HOo730 90o0113 .0.8480 .0*967 V000148 90*8600 t
11.277 .060191 .0o8750 .1.553. 90.0226 .0.8900 ,
.1.71.5 .0.0248 .0*8912 #1@952 ,0.0276 90s8940 t
1(2.080 .0.0290 ,008955 .2.244 .000307 .0.8970 9
12e412 ,0.0323 .0.8961 ,2.646 90*0343 ,008900 ,
K2.920 90o0362 .0.8790 .3.157 .0&0376 ,0.8671 9
M43.;.85 .0*0382 .0.8600 /
DATA (CN( 6,J),D*4f 69J~tETA( 6,J)tJo1,15lI
00.36!7 .0.0061 t0.7600 90.611 .0#0103 .008000 9
P0.830 .0.0139 ,0*8450 .)1*095 .0.0181 ,0*8600 9
Q1*277 .0.0208 .0.8730 .1.515 .0.0243 90*8900 t
1
R1.r7 3 90.0273 .0.8950 .1.934 .0o0300 .0.9000 9
:ý2.194 #0.0334 .0.9005 P2&409 .0.0362 ,0.9010 I
T2.628 pe.0390 .0.9004 ,2o837 #0#0415 .0.9000 ,
U.3.111 ,0.045-7 .Oos80 t3.230. .0.0485 9008800 t
V3.285 40.0510 r0.8735 /

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 -,

L29285 t0.0514 v0.8480 1


DATA (CN( 9tJ),0E4( 9.J)vETA( 9tJ),J=1.15)/I
.0.*012-8 4007000,
H10,365 90,0086 #0.6780 900547
NO.657 t0.0153 90.7125 .0.821 ,0.0,190- .0.733#0 .
01,004 P0.0226 t0.7690 .1.159 .0.025'6 ,0.8000 p

P1.223 10.06270 t0.8060 91.369 .0.0301 '0`8225


Q1.5i5 90#0334 v0.8395 t1.606 .0.0357 ,00.856.,
R1.733 .0.0390 t.0.470 91.825 .0.e0421 ,0.84_45
S1.879 90.0446 90,8330 110898 .0.0463 v0.8235
T1,907 .0.0485 .0.8080 /
DATA CCN(10,J).DH(1oJ),ETA(10.i).J'1,12Z)/
110.365 10.0098 .0.6380 .0.474 t.0.01Z t.0.6550 4~
VO.639 .0.0167 t0.6700 .0.748 20.0195 .0.6850 ,-
WO.85b .0.0223 t0.7000 .0,967 .0.0251 .0.l7110 3

X1.095 90,0290 .0.7180 S1.191 .0.0321 -i017180- _9.


Y1.259 90.03248 .0.7170 .1.296 .0.0ý368 SOP'?7140 ,

Z1.325 90,0390 .0.7000 .1.329 #0'00401 $0.-6690.1


DATA (CN!11,J),DH1I(11J3tETA(11tJ)*JxlI 9)/
AO*365 .0.0111 .0.6000 90.438 .1000137 .0.6000 .
80547 .0.0169 .0s6120 .0.639 .0.01-97 SIPa-61i70
CO.730 .0.0226 .0.6210 .0.821 ,0i.0254 .0.6258- v:
00.876 .0.0275 g0.6250 90.931 .0.0029" .0.6230
EO*963 90.032.2 90*6009 I
END

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 --

SOATAIN Tl~s24O0.,AN.O.OAITPm0000.. ITITLE--1S


SLS 141LITARtY
SUATAIN T4s2400.eA)Sa0.O.ALTPi~O0o0, TITLtsZ.
IAFTBR.1*T1E3Z0o. ETA00.91S

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

Al LOW ALTITUDE DASH


SW~AIN T4x21OO. ,A~sOe.S ALTP=36100., XFITLE=IS
SUBSONIC CRUISE
SDATAIN 14in24OO.,Ahz1.2,ALTP=5OOOO.,ITITLEz1S
SUPERSONIC AT MILITARY POWER
ioATAistE 4z2400. .AMI&6,ALTPZ50000.,ITITLEZI
- SET-UP SUPERSONIC WITH AFTERBURNER

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
'~

4 SUPERSONIC FULL A/B


iEOF
UNCLASSIFIED
Security Classification
DOCUMAENT CONTROL DATA - R&D
(Security claevillcatemn atfihlft, body of abstract and Indesing an~notation viust he .nta...d . ovrepall#P
hdim~Me~U~ i. abfltdo
I O"IGINATING AC1~vI'Y (Corporate author) go. RZPwA1T £UCuNRITV C LASOIPICATION

Air Force Aero Propulson Laboratory Uzolaalfled


Wright-Patterson Air Force Base, Ohio 45433 11 GOU
3 REPORT TITLE

SIMULATION OF TURBOFAN ENGINE


Part EL. User's Manual and Computer Program Listing
T
4 OESCRIPTIVZ NOTES ( p" of epeart and Inclusive dadva)

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)

b6PROJECT No. 3066 APAPLTiR~t67-1.25. Part 11


Task No. 306603 11. . = O- .00
W(An 0%aouwf)mg eo.s

a-A VAIL AUILITY/LIMITATIONNOGTICES This document is Eubject to special;;xpv-nt ontrols a~ c


transmittal to foreign govenet or foreign nationals may be made only Withprlorý
of the Air Force Aero Propulsion Laboratory, WrIght-Patterson AME Ohilo.
II - SUPPLEMENTARY NOTES 12 SPONOPUNG 0ILITARY ACTSWIT'

13. ASTRAC Force Aero P3ropulsion X~bOrmtor

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.

(Distribution of the abstract ie unlimited.)

DD FJAONt'e0473 UNCLASSIFID
Sec'ulty Clasefficatioa
UNC LASS! FLD
SecurityClassification ______ ____________

KEY WORDS
____________________________________________
-IKAL4K6LN
POLK WT POLK uT_ POLK WT

Tudbln Engine Simulation


Turbtirw Engine Computer Programs
Turbine r-agine, CywIe Analyas

0~~ R-I-I- -. '

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

You might also like