AEDC-TR-.78-63
COPUE LEVELEI1,A COMPUTER PROGRAM FOR THE AERODYNAMIC
DESIGN OF AXISYMMETRIC AND PLANARONOZZLES FOR SUPERSONIC AND
SHYPERSONIC WIND TUNNELS
J. C. SivellsARO, Inc., a Sverdrup Corporation Company
VON KARMAN GAS DYNAMICS FACILITYARNOLD ENGINEERING DEVELOPMENT CENTER
AIR FORCE SYSTEMS COMMAND.ARNOLD AIR FORCE STATION, TENNESSEE 37389
LULDecember 1978
Final Report for Period December 1975 - October 1977
L Approved for public release; distribution unlimited.
Prepared for D D D
ARNOLD ENGINEERING DEVELOPMENT CENTERIDOTR JAN 8 1979ARNOLD AIR FORCE STATION, TENNESSEE 37389
D
79 0BEST AVAIL.ABLE -COPY
"4.
NOTICES
When U. S. Governiment drawings, specifications, or other data are used for any purpose otherthan a definitely related Government procurement operation, the Government thereby incurs noresponsibility nor any obligation whatsoever, and the fact that the Government may haveformulated, furnished, or in any way supplied the said drawings, specifications, or other data, isnot to be regarded by implication or otherwise, or in any manner licensing the holder or anyother person or corporation, or conveying any rights or permission to manufacture, use, or sellany patented invention that may in any way be related thereto.
Qualified users may obtain copies of this report from the Defense Documentation Center.
Reterences to named commerical products in this report are not to be considered in any senseas an indorsem*nt of the product by the United States Air Force or the Government.
This report has been reviewed by the Information Office (01) and is releasable to the NationalTechnical Information Service (NTIS). At NTIS, it will be available to the general public,including foreign nations.
APPROVAL STATEMENTI LThis report has been reviewed and approved.
Project Manager, Research DivisionDirectorate of Test Engineedng
Approved for publication:
FOR THE COMMANDER
ROBERT W. CROSSLEY, Lt Colonel, USAFActing Director of Test EngineeringDeputy for Operations
UNCLASSIFIEDREPORT DOCUMENTATION PAGE BEFORE COMPLETING FORM
(6 W. COMPUTER.PROGRAM FOR THE AERODYNAMIC nal epaDc 7
J. C./Sivellsl ARCO, Inc., a Sverdrup
Air Force Systems Command ogram Element 65807?
ArnldAiFrcSatonTenese_338__
16. AGERSLI N C S NAT MET Ao Ihi. OREpo(f dfeetfolotoln fi 1,SCRT LS,(fti eot
* IApproved for public release; distribution unlimited.
17, L3STRIOUTION STATEMENT (of tho abstract onteted ir Block 20, It different from Report)
IS. SUPPLEMENTARY NOTES
Available in DDC.
19. KLY WORDS (Contitnuo on reverse side It nce...ry and Identify by block "limber)
wind tunnel design boundary layerstransonic nozzles
* J supersonic nozzles* *..hypersonic nozzles
exhaust nozzle performance computer program2 0. ISTYAA CT (Con linue on .nre r.o .1 do I I nteo eeary and Identify by block numbe r)
A computer program is presented for the aerodyna-mic design ofaxisymmetric and planar nozzles for supersonic and hypersonic windtunnels. The program is the culmination of the effort expended atvarious times over a number of years to develop a method of de-signing a wind tunnel with an inviscid contour which has contin-uous curvature and which is corrected for the growth of theboundary layer in a manner such that uniform parallel flow can be
DD FOR 1473 EDITION OF I NOV 65 IS OBSOLET E
~~~ ~UNCLASS 7O 1 2 2IT %j 01 0ow 1
UNCLASSIFIED
20, ABSTRACT (Continued)
expected at the nozzle ex-t. The continuous curvature isachieved through specification of a centerline distribution
* of -'elocity (or Mach number) which has first and second deriva-tives that 1) are compatibl* -Lth a transonic solution nearthe throat and with radial flow near the inflection point and2. approach zero at the design Mach number. The boundary-layergrowtb ia calculated by solving a momentum integral equationby nuws'lcal iaitegration.
AILEV YEJi, W ii,oil- D D C
JAN 8 1979
IT ............... ........... .............
Slt. AVAIL
11B
UNCLASSIFPTD
-- ,,,.
AE DC-TR-78.63
PREFACE
The work reported herein was conducted by the Arnold Engineering
Development Center (AEDC), Air Force Systems Command (AFSC). The results
of the research were obtained by ARO, Inc., AEDC Division (a Sverdrup
Corporation Company), operating contractor for the AEDC, AFSC, Arnold
Air Force Station, Tennessee, under ARO Project Numbers V33A-A8A and
V32A-P1A. The Air Force project manager was Mr. Elton R. Thompson.
The manuscript was submitted for publication on September 12, 1978.
The author wishes to acknowledge the assistance of Messrs. W. C.
Moger and F. C. Loper, ARO, Inc., for providing thr. basic subroutines
for smoothing and spline fitting, respectively, which were adapted for
use with the subject program. Mr. F. L. Shope, ARO, Inc., provided
technical assistance in the preparation of this report. Prior to the
publication of this report, the author retired from ARO, Inc.
'I
AEOC-TR-73-63
*l CONTENTS
1.0 INTRODUCTION ................... ...................... 5
* 2.0 TRANSONIC SOLUTION ........... .................. ... 10
3.0 CENTERLINE DISTRIBUTION .................. ......... 154
4.0 INVISCID CONTOUR ........... ................. .... 20
5.0 BOUNDARY-LAYER CORRECTION ..... ............. ..... 24
6.0 DESCRIPTION OF PROGRAM ....... ................ ..... 34
7.0 SAMPLE NOZZLE DESIGN ........... .............. .... 38
8.0 SUMMARY .............................. ... 41
REFERENCES ............. ...................... .... 41
ILLUSTRATIONS
Figure
1. A Foelsch-Type Nozzle with Radial Flow at the
Inflection Point .............. .................. 6
2. Nozzle with Radial Flow and a Transition Region
to Produce Continuous Curvature ............ 7
3. Nozzle Illustrating Design Method of Ref. 13 , , , 8
4. Nozzle Throat Region ............ ............... 9
5. Relationships Obtained from Cubic Distribution
of Velocity from Sonic Point to Point E for
Axisymmetrin Nozzle ...... ................. .... 18
6. Limitations of Fourth-Degree Distribution of
Mach Number from Eq. (39) ..... .............. ... 19
7. Characteristics Near Throat of Nozzle with R = 1 . 23
8. Variation of Wake Parameter, A, w.,ith
Reynolds Number (Incompressible) ... .......... ... 28
9. Variation of Skin-Friction Coefficient with
Reynolds Number (Incompressible) ... ........ .... 29
10. Variation of Velocity Profile Exponent with
Reynolds Number Based on Boundary-Layer
Thickness ...... ......... .................. ... 30
3
h 4;
A EO C-TR-78-63
TABLE
Page
1. Input Cards for Sample Design ...... . . . . . . . 39
APPENDIXES
A. TRANSONIC EQUATIONS .................. . .. ... ... 45
0. CUBIC INTEGRATION FACTORS. .. . . . . . . . . . ... . 49
C. INPUT DATA CARDS . ....... . . . . . . . . . . . . .52D. COMPUTER PROGRAM................ . . . . . . . . . . 61
:"NOMENCLATURE ....................... 139
4
AE DC-TR-78-63
1.0 INTRODUCTION A
Supersonic and hypersonic wind tunnel nozzles can be placed in two
general categories, planar (also called two-dimensional) and axisymmetric.
Early supersonic nozzles (circa 1940) were planar for many reasons: the
state of the art was new with regard to both the design and the fabrica-
tion; the expansion of the air - the usual medium - was in one plane
only, thereby simplifying the calculations and requiring two contoured
walls for each test Mach number and two flat walls which could be usedfor all the Mach numbers; and the relatively low stagnation temperatureand pressure requirements did not create dimensional stability problems
in the throat region. Dimensional stability would in later years become
a primary factor in the development of axisymmetric nozzles.
Prandtl and Busemann, Ref. 1, laid the foundation for determining
the inviscid nozzle contours by the method of characteristics. Foelach,
Ref. 2, simplified the calculation of the contour by assuming that the
flow in the region of the inflection point was radial, as if the flow
came from a theoretical source as illustrated in Fig. 1. The downstream
boundary of the radial flow is the right-running characteristic AC from1',- the inflection point, A, to the point, C, on the axis of symmetry where
the design Mach number is first reached. The flow properties along thischaracteristic can be readily calculated; and inasmuch as all left-
running characteristics downstream of the radial flow region are straight
lines in planar flow, the entire downstream contour can be determined
analytically. Upstream of the inflection point, it was assumed that the
source flow could be produced by a contour which was a simple analytic
curve. In the Foelsch design the Mach number gradient on the axis is
discontinuous at the juncture of the radial flow region and the begin-
ning of the parallel flow region. This discontinuity produces a dis-
continuity in curvature of the contour at the inflection point and at
the theoretical exit of the nozzle.
5
.- .
AE DC-,TR-78-63
UniformrS~Flow
Figure 1. A Foelsch-type nozzle with radial flowat the inflection point.
As the state of the art progressed, it became desirable to cover a
range of Mach numbers without fabricating different nozzle blocks for
each Mach number. A limited range of Mach numbers could be covered by
using blocks with unsymmetrical contours which could be translated
relative to each other to vary the mean Mach number in the test section.
The widest range of Mach numbers with acceptably uniform flow in the
test section has been obtained in wind tunnels in which the contoured
walls consist of flexible plates supported by jacks which can be adjusted
to vary the contour to suit each Mach number. Inasmuch as the curvature
of a plate so supported must be continuous, methods of calculating
contours with continuous curvature were developed (Refs. 3, 4, and 5)
by introducing a transition region, A B C J, downstream of the radial
flow region (see Fig. 2). The shape of the wall between points A and J
was controlled to give continuous curvature. The contours used for the
von K~rmgn Gas Dynamics Facility 40- by 40-in. Supersonic Wind Tunnel
(A) at AEDC were obtained by the method of Ref. 5. Not only is a
continuous-curvature contour easier to match with a Jack-supported plate,
but it also satisfies the potential flow criterion for zero vorticity,,
dq/dn - Kq (1)
where q is the velocity measured along a streamline of curvature K and n
is the distance normal to the streamline. Inasmuch as the inviscid
contour is a streamline, this criterion implies that the flow will be
disturbed where a contour has a discontinuity in curvature.
6
AEDC-TR-78-63
D
B C
Figure 2. Nozzle with radial flow and a transitionregion to produce continuous curvature.
The usual wind tunnel criterion concerning temperature is that the
constituents of the gas should not liquefy during the expansion process
required to reach the test Mach number. For the usual pressure levels
involved, ambient stagnation temperatures can be used up to a Mach
number of about five. As the stagnation temperature is raised, dimen-
sional stability becomes more difficult to maintain in a planar nozzle.
Therefore, axisymmetric nozzles are used when elevated stagnation tem-
peratures are involved. Axisymmetric nozzles have also been used for
low-density tunnels (Ref. 6) because their boundary-layer growth is more
uniform than that of planar nozzles, which inherently have transverse
pressure gradients on the flat walls. The obvious disadvantage of
axisymmetric nozzles is that each one must be designed for a particularMach number. Moreover, disturbances created by imperfections in the
contour tend to be focused on the centerline.
Before the advent of high-speed digital computers, it was extremely
time consuming (Ref. 7) to calculate axisymmetric nozzle flow by the
method of characteristics (Ref. 8). Inasmuch as the assumption of
source flow saved time in designing a planar nozzle, it was logical to
use source flow as a starting point in the design of an axisymmetric
nozzle. In Ref. 9, Foelsch develops an approximate method of converting
the radial flow to uniform flow. Beckwith et al., Ref. 7, show that
Foelsch's approximations were quite inaccurate but utilized the idea of
7
./.
AE DC.TR-78-63
a region of radial flow followed immediately on the axis by uniform
flow, as in Fig. 1. As in the case of planar flow, the discontinuity in
Mach number gradient on the axis produces a discontinuity in curvature
on the contour (Ref. 10). Such discontinuities have been eliminated by
the design methods of Refs. 10, 11, and 12; here, an axial distribution
of Mach number (or velocity) between points B and C (Fig. 2) introduces
a transition region between the radial and parallel flow regions, thus
gradually reducing the gradient and/or second derivative to zero from the
radial flow values at the beginning of the parallel flow. As shown in
Fig. 3, the upstream boundary of the radial flow region is a left-running
characteristic from the inflection point, G, to the axis at point E. The
flow angle is the same at points G and A. Both are shown to illustrate
a general nozzle design. As described in Ref. 12, the contour upstream
of the inflection point can be calculated for an axial distribution of
velocity in the region between points I and E, which makes the transition
from sonic values to radial flow values. On the axis, the sonic values
of first and second derivatives of velocity with respect to axial distancewere calculated by an adaptation of the transonic theory of Hall, Ref.
13, or Kliegel and Levine, Ref. 14. The upstream limit of these cal-
culations was the left-running characteristic from the sonic point on
the axis.
Inflection Reglon
M * 6°
SI F -
Figure 3. Nozzle illustrating design method ofRef. 13,
8
-: • .'*. ..'. U,,.J-
AEDC-TR-78-83
This characteristic is also called a branch line. Between the theo-
retical location of the throat and the intersection of the branch line
with the contour was a region which was not calculated but which in-
creased in size as the throat curvature increased. This gap in the
contour has been eliminated by the method described herein which util-
izes a right-running characteristic originating at the throat as shown
in Fig. 4 (where point I has been moved from the sonic line to the
throat characteristic). With this latest improvement upon the method of
Ref. 12, contours can be designed which have throat radii of curvature
of the same order of magnitude as the throat radii although such an
extreme curvature would not normally be recommended from other stand-
points. A recent (1975) design of a Mach 6 nozzle utilized this method
with a throat radius of curvature of about 5.5 times the throat radius.
Sonic Line Branch Line
YO! Throat Characteristic
Figure 4. Nozzle throat region.
After the design method was developed for axisymmetric nozzles, it
was adapted for planar nozzles having a prescribed centerline dis-
tribution of Mach number (or velocity). This approach to such a design
is considerably different from that of Ref. 5. The current design
method is incorporated into the computer program included herein. As an
option in the program, a complete centerline Mach number distribution
9
AEDC-TR-78-63
can be used which does not include a radial flow region. Parts oZ the
computer program are subroutines for computing the boundary-layer correction
to the inviscid contour, for smoothing the contour, and for interpolat-
ing points at even axial positions by means of a cubic spline fit of the
contour.
2.0 TRANSONIC SOLUTION
In many early nozzle designs, it was assumed that the flow at the
throat was uniform (M - 1) and parallel. This assumption impliey that
the wall curvature is zero and that the acceleration of the fl, Is
zero (i.et, the acceleration starts from zero at the beginning .-. the
contraction, reaches a maximum in the contraction but is reduced to zero
again at the throat, and must be increased again in the beglaning of the
supersonic contour and reduced to zero at the nozzle exit). A nozzle so
designed therefore becomes considerably longer than one in which the
flow reaches its maximum acceleration in the vicinity of the throat,
where it is approximately proportional to the reciprocal of the square
root of the radius of curvature. The above argument indicates the
fallacy of some so-called "minimum length" nozzles, although some
designers have combined a contraction having a relatively high throat
curvature with the supersonic section having zero throat curvature.
For a throat with a finite radius of curvature there have been many
transonic solutions. Hall, Ref. 13, developed a small perturbation
transonic solution for irrotational, perfect gas flow, in both two-
dimensional and axisymmetric nozzles, by means of expansions in inverse
powers of R, the ratio of the throat radius of curvature to the throat
half-height, or radius. His solution gives the normalized (with the
velocity at the sonic point) axial and normal velocity components in the
form
(y,Z) + . . .. + + .Y. . . (2)R R2 R• 3
10
10
AEDC-TR -78-63
2
- + 12 Lv (y,z) v+(y,z) + v,(y,) V(L,( 1 2 + L. .( 3 )
where
and x and y are coordinates normalized with the throat half-height or
radius, yo. The value of a is zero for two-dimensional flow and one for
axisymmetric flow. Kliegel and Levine in Ref. 14 extended the applica-
bility of Hall's axisymmetric solution to lower values of R essentially by
making the substitution
- S-1 + 9-2 + S- . .(5)
where S - R + 1, into Eqs. (2) and (3). In the method used herein, the
same substitution is made in Eq. (4) for two-dimensional flow as well as
for axisymmetric flow and therefore becomes a special case of the general
transonic solution described in Ref. 15. The complete general equations
in terms of S are given in Appendix A.
At the throat, x - 0, y - yo' v 0 0, for planar flow,
n= 1+ _1_ - (14y-7 5) + (274Y2
- 861y,+ 4464)3S 270S
2 1701003
d u . -(32y2 + 87y 561) +
dx/y. S 540S2
and, for axisymmetric flow,
+ 1-4y -57) + (2364y 2 - 3915y + 14337) + .. . (8)4S 288S 2 82944S3
Sd__+__ = _ -(64y2 + I 17y - 1026) + (9)dX/yu as 1152S
2
;!: 11
AEDC-TR-78-63
where the derivatives are with respect to x nondimensionalized by the
throat half-height or radius, respectively, and
2
A _____(10)
On the axis, y - 0, v - 0, for planar flow,
u = 1 - _i + Y-5 782y 2 + 3 5 07y - 7767 +6S 270, 2
2721 60S'
+ + 34y2+ 429y_ 12+].4Y" 4320S 2 1
,Y," 6 36S
(,& -z' (2y 2 - 33y + 9)/72 + . (1)
and, for axisymmetric flow,
S1 - _. + IOy-. 1L5 _ 2708y2 + 2079y + 2 1 1 5+4S 288S 2 82944S 2
+& -1 -t I + 92y2 + 180X - 9 ++Y \ 8S 115282
()2 (.2~ + +
(--A)' (4Y2 - 57y + 27)/144 + (12)Yr,
Because the sonic line is curved for finite values of R, the mass
flow through the throat is reduced by the factor CD (discharge coef-ficient), which is the rat 4 .o of accual mass flow to that which couldflow if R were infinite and the sonic line were straight. For planar
flow,
CD I- [i - 4y24 + 334•2 - 457y + 4353 + (] 13)
I12
i 12
AEDC-TR -78-63
and, for axisymmetric flow,
C[ 1 - y [1 - - .754y2 -757y +36U. + (14)96 S' 24S 28EIOS'
The flow which passes through the throat also passes through the
sonic area of the source flow which is at a distance r from the source.
In planar flow, 0(
•! or
oY0/r1 = ,/C1 (16)
where the inflection angle, n, is in radians.
In axisymmetric flow,
7y. = ry~o CD =21 r• 1 -cs (17)
or I
yo/r- 2 sin (1/2)/cD (18)
In the calculation of the throat characteristic used herein, thevalue at x - 0, y - yo Eq. (6), is the starting point. The half-height'•~
or radius, y0, is divided into 240 equally spaced values of y. Inasmuch
as the characteristic is right running, its slope at each point is
dy/dx tan( (19)
wheresin I = 1/M (20)
Also
W = M + K-l M2) (21)y+1 Y+-
sin • = v/W (22)
and
drdo + d- o '""a d• (23)
y
13
- -
AEIDC-TR-78-63
Sde dx/cos(o - t)- dy/sin(o - 1) (24)
'The term i is the Prandtl-Meyer angle in two-dimensional flow,
y an. (2I - - tan-1 (I2 -1) (25)Y-1 y+l
Equations (19) and (23) are the characteristic equations and are solvedby finite differences. If all values are known at point 1, the valuesat point 2 are found (y is known a t both points) by
S + 2(y-y) (26)x2 - 1 tan( l-A, + tan
A (Y2 -Y)2 + (x2! - x1 )2 (27)
V02 = 01 + 01 - 0 2 + VI + - V 2 (28)
At the starting point W is the value of u because v - 0. Values of v2I
are calculated at each point (x 2 , y2 ) from the transonic solution,
and Eqs. (26) to (28) are iterated until convergence is reached. For
evaluating the term in brackets in Eq. (28), the ratio v/y is defined by 4the transonic solution even on the axis where both v and y are zero.This fact eliminates thbi general problem ia axisym-etric characteristics
solutions of evaluating the indeterminate sin */y in Eq. (23) on the
axis of symmetry.
It may be noted that the value of W as calculated from the character,-
istic value from Eq. (21) differs from the value (u2 + V 2)I12 calculated
from the transoutic equations, but the difference decreases wT:Uh in-
creasing R. For the final point of the throat characteristic which 1,us
on the axis, the value of d 3 u/dx3 from the transonic solutioa for the
axial distribution is "corrected" to make u - W for the axisymmetric
case fur values of R less than 12. The correction is abvut 16 perccnt
for R I and decreases rapidly as R increases. This correction is made
14
•
F • ' " i I i ' i I ii i i i ln n. . " . . ..
AE :)C-TR -78-63
so that v,'alues of du/dx and d udx 2 can be calculated from the transonic
solution for later application. The correction is believed to be JuatifiedI1•'•much as the accuracy of the transonic solution -s limited, particularly
for low values of R, because the series expression for u is truncated
after the x• term.
3.0 CEN'TERLMN DISTRIBUTION
In the radial flow region, the distance r, measured from the sourne,
is related to the locAl. Mach number by
Lr)+ _ (_L + h.-y+ I+ M 2)qy (29) i 'r Y+1 ++
or
)1+a W-(Y+' I W2-- (30)
First, second, and third derivatives of W -it M with respect to r/r can
be obtained as described in Ref. 12. Along the axis x - r when x is
measured from the source. Inasmuch as all, coordinates ma.st be normalized
by the same factor, r,, the transonic equation in terms of A/yo and y/yo
can be transormed by Eqs. (16) and (18), after which the distance froom
the source to the throat station must be taken I.-ntc account. This
latter distance is generally unknoun until after the distance from point
I to point E is determined.
In radial flow, the term on the right-hand side of Eq. (23) can beevaluated dimply. Inasmuch as sin y/r and d& - dr/cos v,
H ik9! S21 tan fy r
but
. -tan = (M2
-I ) 2
"15
AEDoC-TR-78-63
and, from Eq. (29) 'for a - 1,
. dr (2- 4M
r 2(1 +2 E: M2) Mi2
Thus,tnan s dr = (M 2 -1)
2 dM
r 2(1+-Z2 MN2 ) M2
I
From Eq. (25), do M(I +Y-1 M2 ) M
2
"therefore, Eq. (23), in radial flow, becomes
dot/ + do4 - 2- dot (31)21
which applies for characteristic AB or GF. Similarly, for the left-
running chaivacteristic EG,
do - do - ad (32)Thc.:efore, 2
13 - A + -G (33)
and G - OF, = (a + 1) (34)
and, from the design values n and MB (and/or ), MA, M, M' WE, and
the necessary derivatives can be calculated.
Within the accuracy of Eqs. (11) and (12), the second derivative of
velocity ratio at the sonic point is negative for values of R loss than
11.767 for planar flow and 10.525 for axcisyrmetric flow. The second
derivative of Mach number at the sonic point is positive for all values
of R. Inasmuch as the second derivative of either W or M is negative
for source flow, it seems better to use a velocity distribution rather
thau a Mach number distribution between points I and E. On the other
hand, a Mach number distribution between points B and C is preferable
16-IT
A E DC-T R -78-63
because the velocity ratio approaches the constant value of [(Y +
1)/Y _ 1)] 1/2 as the Mach number increases to infinity; therefore, the
change in velocity between points B and C becomes small relative to the
change in Mach number.
I.
The velocities and their first and second derivatives at points I
and E are used to determine the coefficients of the general fifth degreepolynomial
.,W C1 + C2 X + C3 X2 + C4 X3 + C'5X4 + C6X5 (35)where
X - (x - xl)/(XF - X[) (36)
Similarly, the Mach numbers and their first and second derivatives at
points B and C are used to determine the coefficients of the polynomial
M = 1) + D2 X + D3 X2 + 0)4 X + D5 X4 + D6X" (37)
where, in this case,X .: (x - XI)/(X. - xv) (38)
and the first and second derivatives at point C are usually set equal to
zero.
In these equations, the lengths (xE-xI) and (xC-xB) must be specified,
but can be determined by the conditions that C6 and D6 equal zero,
thereby reducing the polynomials to fourth-degree ones. If the velocity
at point E is determined by iteration, the third derivative at point I
or E can be included as a criterion for the fourth-degree polynomial;or, by setting C5 - 0, one can find a third-degree polynomial with a con-
stant third derivative. In either case, the Mach number at point B is
found from Eqs. (33) and (34) after the value at point E is found. All
"of these options are included in the program, but unless there are other
factors involved, the preferred options are the cubic between points I
and E and the quartic between points B and C.
17
- .'i- T-n i n
I
AEDC-TR-78-63
* For the cubic distribution for axisymmctric flow, the Mach number
at point E is related to the radius ratio as shown in Fig. 5 for y
1.4 for various values of inflection angle. Cross plotted are lines
of constant values of the ratio * E/n. Such values for most axisymmetric
nozzles lie in the range covered in this figure, and inasmuch as *F/n
h + 4, values of MF can also be obtained.
2,8
2.41,4 "-2-2, 22
1, 8 - 12
.11.6
1,2 4
0 4 8 12 16 20 24Radius Raik (R)
Figure 5. Relationships obtained from cubicdistribution of velocity from tonicpoint to point E for axisymmetricnozzle.
j18- - - - -' ,. -
AEDC-TR-78-63
In determining the length of the segment between points B and C,
using the fourth-degree polynomial distribution, there is a minimum
value of the Mach number at point B for the design Mach number at point
C. As given in Ref. 12,M= ,; + 0.715 M " /M (39)
where the primes indicate derivatives with respect to r/rI. This
relationship is shown in Fig. 6. For an axisymmetric nozzle designed
for a Mach number greater than about 3.4, the minimum Mach number at
20
+ ~18 /
16
14
,, 12
10 Axisymmetric
Planar
0 I . , I
0 2 4 6 8 10 12 14Minimum Mlach Number at Point B
Figure 6. Limitations of fourth-degree distribution
of Mach number from Eq. ?39).
point B is about two-thirds of the d2sign Mach number. Using such a
value visually causes the length to be •xcessive, and more realistic
19
zzzizI,7,.: ,. 71
* AEDC-TR-78-63
values of B are 75 to 80 percent of M It is important, however, as
illustrated in Ref. 16, that the distance between points B and C be
sufficient to allow for accurate machining of the contour between
points A and J, which lie on the characteristics through points B and C,
respectively.
4.0 INVISCID CONTOUR
The flow properties are determined at a desired number of points
along the key characteristics (i.e., the throat characteristic, TI, as
described earlier (a sub-multiple of 240 is used for subsequent calculations),
the characteristics EG and AB bounding the radial flow region by Eqs.
(33) and (34) for equal increments in n, and the final characteristic CDalong which the Mach number is constant and the flow angle is zero).
The flow properties are also determined at axial points from Eqs. (35)
and (37). The network of characLeristics is then calculated in the
region TIEG starting at point E and progressing upstream and in theregion ABCD starting at point B and progressing downstream.
The equations for a right-running characteristic were given previously.
dy/dx = tan(o - p,) (19)
+ ,. d (23)
wherede - dx/cos( -It) dy/sin(r/ - p) (24)
For a left-running characteristic, the equations are
dy/dx tan (0 + p) (40)
do- =b '!2s_0nL,uiJ d4 (41)y
where
Aldo dx/cos (0 + p) = dy/sin(0 + p) (42):, Also
do - .~J- dM acot I dW (43)(1 + 1.i M2) M W
2
[I 20
A__ __
AEDC-TR-78-63
Values of x, y, *,and M are known at the general point 1 on theright-running characteristic, C, and at the general point 2 on the left-
running characteristic, C. The characteristics intersect at the general
point 3 where the values are calculated by numerical integration of Eqs.
(23) and (41) along the respective characteristics.
S(si 6,mn A + , 02. n P~2 ) (44)
whereA4 c -( " 2) sec 1 (45)
andY3 Y2 tan~ 0 tan (0 3 + jA3) + 1. tan (02 + A2 (46)
3 222
-3 011 + (3- 01) = =
4, (Min0mi&AMn +. sinl 6minp (47
where (3- x1) sec. a (48)
and
=tan a Itan ((3 13 +1L tan(~ (49)X -X 2
Adding, substracting, and rearranging gives
03 (0 + - 0 + S6 + P2 + )(50)2
03 -(01 2 + 0,+ 02 + --(51)2
In planar flow, P1 p2 0 because a -0 and Eqs. (50) and (31)
can be solved directly, 14 is obtained from *j by t~he inverse applica-14-1
tion of Eq. (25), and vNsin (1/M ). In axisymnietric flow, the equa-
tions must be solved by iteration. A useful first approximation for P1
and P2 is the radial flow values, P1 -(p- 1)/2 and P2 ) ~ 3 /2.
21
.-- . .. ..
AEDC.TR-7S-63
At all points except on the axis in axisymmetric flow, Eqs. (44)
and (47) are defined because Y2 and y, are nonzero. On the axis, the
terms sin /Y2 and sin 0,/y, are indeterminate with the form zero/zero.
These indeterminates can be evaluated by assuming that the general
points 1 and 2 on the axis are very close together and that U f P s 13
and W OW W2 * W3 . Equation (41) can be written
Cot g L do + 'in 6 sin I dx (52)W y cos ( + -)
and Eq. 23 can be written
cot A dW= -d9 + sin 0 idnp dx (53)W y cos( S -)
as 964 9S-Psino, O±IA-.±ts
and tan 3 = i__i..-- =. 33 2 1 3
In finite-difference form,
co int an• IA.. (X3 -x2t a (W3 - W2) = 3 + . 3 3 2)
0 tan X sin 0. tan P (X3 - X2)
Y3 Y3
-- 2 sin 03 tan M3 (x3 - x2)/y 3 (55)
Similarly
S(Wl - W3) =3 + sin 0, tan /I3 (x1 - x3)/y 3 (56)Wa
-* 2 sin 03 tan p 3 (xl - x 3)/y 3 (57)
Adding Eqs. (55) and (57) and rearranging,
lia • =i °~ __ w (58)y4O Y 2 W dx
and s"' 0'2 Min (M1 -)Y (d' (59)
2 2
. , ,. . - . , ,, . , ./ '. , . . . , 2 2
AEDC-TR-78-63
for use in Eq. (44) when point 2 is on the axis, and
HilIl H st 410Yl 2W W \I (60)
for use in Eq. (47) when point I is on the axis.
In starting the calculation of the network of characteristics in
the region TIEG, point E becomes point 1 and the first axis point up-
steam of point E becomes point 2. The complete left-running character-
istic approximately parallel to EG is calculated, and the point on the
contour is determined from mass flow considerations as described in Ref.
17. The flow properties along this characteristic are then used to
calculate the next left-running characteristic, again starting on the
axis. This process is repeated until point I is reached, after which
the starting point for each left-running characteristic is a point on
the throat characteristic as illustrated in Fig. 7. The process in
region ABCD is similar except that right-running characteristics are
calculated for each point on the contour.1, 4
UI
1.0
0.8 -
Y0.6 "-
0.4 -
0.201
0 0.2 0.4 0.6 0.8 1,0 1.2 1,4 1,6 1.8 2.0x
Figure 7. Characteristics near throat of nozzlewith R = 1.
23
-. 'IL,
AEDC-TR-78-63
5.0 BOUNDARY-LAYER CORRECTION
To each ordinate of the inviscid contour must be added a correction
for the boundary-layer growth to obtain the viscid or physical contour
of the nozzle. Except for very low stagnation pressures, the boundary
layer is assumed to be turbulent. Generally, the boundary-layer cor-
rection will be made for one design condition of stagnation pressure and
temperature although it is theoretically possible to reshape a flexible-
plate type of planar nozzle to account for different boundary-layer
thicknesses corresponding to different stagnation conditions. The
correction for a planar nozzle is usually applied to the contoured walls
only, but the correction also allows for the growth of the boundary
layer on the parallel walls in order to maintain a constant Mach number
along the test section centerline. Therefore, the correction applied is
greater than the displacement thickness on the contoured walls, and the
flow in the test section is diverging in the longitudinal plane normal
to the contoured walls. In the longitudinal plane normal to the parallel
walls, the flow is converging because of the boundary-layer growth;
moreover, there is a tendency for the boundary layer to be thicker on
the wall centerline because of the transverse pressure geadients present
on the parallel walls. Although these physical effects make a true
correction impossible for a planar nozzle, the calculations described
herein are made as if the cross section were circular, with the cir-
cumference at each station equal to the periphery of the actual rec-
tangular cross section.
The method of calculating the boundary-layer growth is based on
obtaining a solution to the von Kfrmfn momentum equation written for
axisymmetric flow.
0I M 2 -M 2 +H dM + I dr C (61)dx M El + (y.-_) M2 /2 dx rw dx 2
The term [(I/rw)(drw/dx)1 becomes an effective one for planar flow as
just described. For either type of nozzle, the inviscid value is used
24
• , . ,.k
AEDC.TR -78-63
as a first approximation. The entire solution is iterated several times
with new values of rw and dr w/dx - tan 0w obtained each time by adding
vectorially the displacement thickness to the inviscid contour.
The value of mcmentum thickness used in Eq. (61) is defined by
PI -__d__(62)
0___ P'I '( 1-- q )d (2
where z is measured normal to the wall.
Also5
iiO. (,_ dz (63)
The quantities P* and 0 may be considered to be the displacement and
momentum thicknesses when the boundary-layer thickness is small with
respect to the radius, r . These values are related to tota2 valuesw
a* and 0, obtained from mass-defect and momentum-defect considerations
by12
- (64)
and -- .• 62 (65)
Because rw-6 coo Cw + y, where y is the inviscid radius, Eq. (64) mayw a~ wbe rearranged to give
4-- w I Y 2 SC2 2 - y sec 9W (66)
For the final correction, the value *a sec w is added to the in-viscid
radius in order that no correction be made to the longitudinal location.
The integrations of Eqs. (62) and (63) are performed numerically
using Gauss' 16-point formula, with the assumption of the power-law
velocity distribution
q/q, (67)
25 1
.-. ....... '..
AEDC-TR-78-63
ai and
pl/pý " TT (68)
whereT = T (Taw - Tw) q/q, + [Te - a - Tw) - Tw] (q/q.) 2 (69)
which is Crocco's quadratic temperature distribution if a - 1. 'However,
as shown in Ref. 12, a value of a - 0 gives a parabolic distribution
which agrees better wit1' data obtained in hypersonic wind tunnels with
water-cooled walls. The same distribution is obtained if T - Taw'w &
which is likely to be the case for planar, flexible-plate nozzles.
Before using the Gaussian integration, one must replace the values of z
and dz with 6(q/q dN and N6 (q/qe)N-1 d(q/qe), respectively, in order to
avoid the infinite slope, dq/dz, when q and z equal zero.
The value of the compressible skin friction coefficient, Cf, in
Eq. (61) is assumed to be related to an incompressible value, Cfi
by a factor F c, introduced by Spalding and Chi, Ref. 18,
Fc Cf =Cr f (70)
and Cf is related to an incompressible Reynolds number, R, which isfI
relateý to the compressible value, Re , by a factor F,c
FIa RO, RO (71)
The factor Fc, also used by van Driest, Ref. 19, is given by
L I/(p/pe)2 d (q/qe)] (72)
which uses Eqs, (68) and (69), In Refs. 18 and 19, a value of a ' 1
was implied, but Eq. (72) is used herein with a - 0 also, to give
a "modified" value of V . The factor Fc may be considered to be thei ~cratio of a reference temperature to the free-stream temperature. The
factor FR , as used by van Driest, isV F68 = ,/JAW (73)
26
-. Ni
Zi
A E DC-TR-78.63
The compressible momentum thickness, 6 upon which R is based is
the flat-plate valueC
li P% q P q" 'z 74
because the values of Fc and FR were developed to correlate flat-
plate data.
The equation used herein for incompressible skin-friction coef-
ficient is that of Ref, 20,
If (log H0 + 4.56]) (log t0o - 0.546)
This equation is believed to agree with experimental data slightly
better than the von Kgrm~n-Schoenherr equation,
C (0,242)2, (fI (log 110 + 1.1696) (log lie + 0.3010) (76)
at high Reynolds numbers. Also as shown in Ref. 20, Eq. (75) agrees
with the equation, Ref. 21, based on Coles' law of the wall and law of
the wake,S~L(2/C1 )2 fn 15 + 0.5 Fn (C /2) + YC + 211 (77)
if 1I varies as shown in Fig. 8 from about 0.41 at R6 - 400 to a maximum
of 0.5885 at Re - 50,000 and then decreases to about 0.49 at R.7i107. in order for Eq. (76) to agree with Eq. (77), I1 must continually
increase with increasing R0 as shown in Fig. 8. The data shown in Fig.
8 were computed by Coles iniRef. 21 from Wieghardt's flat plate data,
Ref. 22. A comparison of friction coefficients from Eqs. (75) and (76)
is shown in Fig. 9 together with Wieghardt's values as recomputed by
Coles. The constants K and C are 0.41 and 5.0, respectively. The
relationship between e and 6 is obtained from the logarithmic velocity
profile by neglecting the laminar sublayer, representing the wake function2
by a sine distribution, and integrating to obtain
____ (78)
j27
A ED C-TH -78-63
and
_8* C= i(~,. 2 179 11 + 1.5 n2) (Y9)
8 8 2K
•,'• I"0.6 ' 0 0
0.50
U From Eas. (76) and (77)'-Iata Tabulated In Ref. 21;
Identifled as Wloghardt Flat Plate Flow
0.3
0.2I~ l ll J I I l lll ] I I ~ ll I I I I I f ll I I I J I ll
103 104 105 106 107
Figure 8. Variation of wake parameter, n, withReynolds number (incompressible).
The value of N in Eq. (67) is assumed to be a function of Reynolds
number based on the actual boundary thickness, not corrected by FR
and is evaluated through the use of the kinematic momentum thickneis
q I- dz (80)from which
Ok/6 N/(N2 + 3N + 2) (81)
or
N -a+ - 6\ + 1I (82)
28
I
A EDC-TR -78-63
0 . 00 6 r I I I I [1 I T l l 1 TrTiJ I l f l I ! I I 1 1 1 1
0,003 5
O.004 Wleqhardt's Flat Plate DataTabulated In Ref, 21Ct,
0.003
0.002 Eq .(75)
0. L0 Eq. (7610
0 1 1 11 tl 1 1 f 111 ilj ,J I J iIlll ,I I I I ILill'
104 1 5 106 Ia7Re,
Figure 9. Variation of skin-friction coefficient withReynolds number (incompressible).
The value of 0k/ 6 is obtained from Eq. (79), where the value of nI isevaluated from Eqs. (75) and (77) with ek used instead of 0i. The re-sulting variation of N with Rd is shown in Fig. 10.
Two options contained in the program subroutine for the boundary
layer utilize Coles' law of corresponding stations (Ref. 23),
(83)
if Cf /Cf . F is calculated from Eq. (72) for a - 0 or a 1 1, then onei c
option gives
F 1 1 8 ~ ~/I~0''~ ~(84)
The second option divdes Eq. (83) into the two parts,
Cfl/( =f I' /T,, P11 (85)
29
AEDC,.TR-78-63
and
= !S/!L.(86)
where pis evaluated at the temperature
'T', . + 1?.2(Cf/2)2 a(Ta - T) 305(Cf,/2) [a(Taw - Tw) + Tw T.](87)
10
.1 8
N
6
104 106R6
Figure 10. Variation of velocity profile exponent with Reynoldsnumber based on boundary-layer thickness.
Still another option defines the incompress~ible skin-friction
coefficient as
Cf 0,0888(log 118 + 4.6221) (log 118 1.4402) (88)
where
CTWAW (89)
and F Cis calculated from Eq. (72).
30
V t.
AEDC-TFR-78.63
The wall temperature in the above equations can be the adiabatic
wall temperature or can be allowed to vary between a throat wall tem-
perature, T WT, and a nozzle-exit wall temperature, T WD, both of which
are input to the program. Two options are available for the variation
of wall temperature,
I- (A i/A*)'- I A/A* (90)
where m can be 1/2 or 1, A/A* is the area ratio corresponding to local
Mach number, and Ac/A* is the area ratio corresponding to the design Mach
number at the nozzle exit. Equation (90) is used in lieu of more
accurate values and approximates the way the heat transfer decreases as
the Mach number increases from 1 at the throat to the design value at
the exit. For a water-cooled throat, the value of T can also becalculated by the program, T
T ' + 00 I('l' l)-1 )'I; T= (91)
h. + 0
where h is the airside heat-transfer coefficient at the throat asacalculated by Reynolds analogy from the throat skin-friction coef-
ficientSp p C 2/ c2 (92)
with a constant specific heat based on the thermochemical BTU
T (y -. i) 777,64885, (93)
and Q is an input which is a function of the properties of the throat
material, the cooling water, and the geometry and would be a constant if
the properties were constant. The assumption is made that the bulktemperature of the water is 15*F less than T and that p2 /3 is the
square of the recovery factor used to obtain Vhe adiabatic wall tempera-
ture, T •
S31
AEDC.TH ,78.63
For "he integration of Eq. '61), the values of x, y, dy/dx, M,
and d&/dy .. : obtained from the iviscid contour at unevenly spaced
points as a result of the characterio:t:cs solution. With the inputs of
stagnation pressure and temperature, ga,• constant, and recovery factor,
the unit Reynolds number and static and adiabatic wall temperatures can
be calculated at the same points as functions of Mach number withSutherland's equation used for viscosity. With the inputs of T andw_
TWD, the wall temperatures can also be calculated as functions of Mach
number, although T may need to be obtained by interation if thewToption to input a value of Q is exercised. Sutherlandts equation is
also used with wall temperatures to obtain the viscodities at the wall.For any static temperature below the Sutherland temperature, 198.72*R asused herein, the viscosity variation with temperature is assumed to be
linear.
The integration of Eq. (61) is started at the throat where it is
assumed that dOidx w 0 in order to obtain a value of e. Iteration is
involved at each point because C f is a function of Reynolds number basedupon 8, and the relations e/6 and 8*/8 depend upon the value of N,
which is a function of Reynolds number based upon 6. After all itera-tions converge within specified tolerances, the value of P is calculated
afrom the value of 6*, and the values of B and dO/dx are used in the
calculation at subsequent points. The values of dO/dx are integratednumerically to obtain the increment in 0 to be added to a previouslydetermined value of 8. The trapezoidal rule is used to determine the
second point, the parabolic rule for the third point, and cubic integra-tion for the fourth and subsequent pointa.
For convenience, Eq. (61) may be written 8' + 8P - Q. The general
integration for the nth point is
On =On3 + Gu_3 1 n 3 + Gn-2 0 n-2 +I Gn 1+ Gn 0. (94)
32
S.... .. . : ': : .... . . .-- • " • -- • ' .... .. • -.. . . . ..• " + + + "+" , i ................. ..... ..........
AF DC.TR-78-63
where th' C's are functiw~s of the spacings s, t, and u between thepoints and are given in Appendix B. Except tor e and ef, the other
n fnvalues in Eq. (94) are known from previous calculations. Inasmuch as
0 Q P 0 (95) 1Eq. (92) can be rearranged to g(ve
on. (011_,1 111 0.,, : 3 1. "_• + G t + G; 0' 4.- G Q. .(6
(0 + Gri )
After convergence of the iterations, Eq. (95) is used to obtain de/dx.
Inasmuch as Eq. (94) depends upon the knowledge ot en-3, the value of
a2 is calculated by
on-2 = On-. 3 + "' n• - + + "'n-1 0,1-1 + p, lon (97)
which becomes the Gn3 for the next point to be calculated. The valuesof the Fls are also given in Appendix B. The values of 02 and e3 obtained
from Eq. (95) are used in the calculation of P* and 6P instead of thea
initial values obtained by the trapezoidal or parabolic integration.
The success of the above type of integration depends upon thespacing of the points. The values of the increments s, t, and u must
be of the same order of magnitude, although t is usually larger than s
and smaller than u if the parameters itnvolved in the characteristicssolution are selected with care.
After the values of 6* sec w are calculated, the values ofa w
d(6* sec 0 )/dx are obtained by parabolic differentiation and added toathe inviscid values of dy/dx to obtain dr /dx. This procedure is believed
wto be more accurate than differentiating the value (* sac w + y)
.a wbecause dy/dx is obtained directly from the characteristics solution andnot by differentiating y with respect to x.
In general, the boundary-layer correction at the throat will have
a gradient such that the viscid throat will be slightly upstream of the
33
AEOC-TR-78,63
inviscid throat. This displacement and the value of the viscid curve-
ture at the throat are calculated using the assumption that both the
inviscid throat and the boundary-layer correction are parabolic in
shape.
6.0 DESCRIPTION OF PROGRAM
The computer program is written in Fortran 1V for use with the IBM
370/165 Computer. The program consists of a main section, three functions,
and 16 subroutines arranged so that the program can be overla4,d toconserve computer storage. The four overlays consist of AXIAL, CONIC,
SORCE, and TORIC; PERFC; BOUND and HEAT; SPLIND and XYZ. The inputdata cards are described in Appendix C, and a listing of the program is
given in Appendix D.
Program MAIN. MAIN calls for the various overlays. The title card is
read in with the designation as to whether the nozzle is planar or
axisymmetric. A card defining the gas properties and a few pertinent
dimensions is then read in. The first subroutine called is AXIAL, in
which the upstream axial distribution is defined. PERVC Is called to
calculate the upstream contour. AXIAL is recalled to define the downstream
distribution, and PERFC is recalled to calculate the downstream contour.
BOUND is called to calculate the boundary-layer growth. SPLIND is
called to determine the coefficients of cubic equations to fit the
unevenly spaced points along the contour, and XYZ uses these coeffici-
ents to obtain ordinates at evenly spaced points along the axis or, in
the case of the planar nozzle, at discrete points along the surface ofthe flexible plate at which the supporting jacks are located.
Subroutine AXIAL. In this subroutine, cards are read in with theparameters used to define the axial distributions of velocity and/or
Mach ,Lumber and with integers which define the number and spacing of thepoints on the axis and on the key characteristics and the sequence of
34
AEDC-TR-78-63
subsequent calculations. If the throat characteristic is called for,
the upstream end of the upstream distribution starts at the intersection
of the throat characteristic and the axis. An option can be exercised
to not use the throat characteristic and thereby start the distribution
at the point where M - 1. This option would normally be used for a
nozzle with a large throat radius of curvature, e.g. a planar nozzle, or
if it were desired to repeat a calculation as in Ref. 13. Another option
is to avoid a radial flow section altogether by using a polynomial dis-
tribution from the throat to the beginning of the test cone or rhombus.
Other options will be described in Appendix C when the input cards are
discussed.
Subroutine BOUND. This subroutine is used to calculate the turbulent
boundary-layer correction to the inviscid contour. The stagnation
conditions are input, as are the parameters to describe the wall tem-
perature distribution, the temperature distribution in the boundary
layer, and the factors relating the compressible skin-friction coefficients
to incompressible values.
Subroutine CONIC. This subroutine is used within AXIAL to give the
derivatives of Mach number with respect to r/r 1 in radial flow from Eq.
(29).
Function CUBIC. This subroutine is used to obtain the smallest positive
root of a cubic equation.
Function FMV. This subroutine determines the Mach number for a given
Prandtl-Meyer angle.
Subroutine FVDGE, This subroutine is used within PERFC in conjunction
with NEO Lo smooth the inviscid coordinates as desired.
35
, .-
AEDC-TR-78-63
Subroutine HEMT. This subroutine is a dummy called by BOUND but is
included so that with a more elaborate subroutine a heat balance can be
made to determine the wall temperature if the material conductivity is
specified and the cooling water passage geometry and quantity of flow
are specified.
Subroutine NEO, This subroutine is used with PERFC in conjunction with
FVDGE to smooth the inviscid coordinates as desired by modifying the
ordinate such that the second derivative is more nearly linear after
smoothing than beforehand,
Subroutine OFELD. This subroutine is used within PERFC to calculate the
properties at the intersection of a left- and a right-running char-
acteristic.
Subroutine OREZ. This subroutine is used to make all values of an array
equal to zero prior to a new calculation.
Subroutine PERFC. In this subroutine, the properties along the key
characteristics are first calculated to go with those along the axis.
The intermediate characteristics are then calculated and the contour
points obtained by integrating the mass flow crossing each character-
istic. If desired, certain designated intermediate characteristics may
be printed out. If smoothing of the ordinates is desired, the inputs
associated with the smoothing are read and the smoothing applied.
Inasmuch as the wall angle is interpolated from mass-flow considera-
tions, independently of the coordinates, the wall slopes are integrated
from the inflection point toward the throat for comparison with the
"interpolated ordinates. Parabolic integration is used for this purpose
as well as for the mass flow. Also calculated for comparisorn are the
ordinates of a parabola and a hyperbola which have the same radius
ratio, R, inasmuch as the transonic solution should be equally applic-
able to these shapes for the number of terms retained in the series,
36U ,,
AEDC-TR-78-63
Eqs. (2) and (3). Finally, the scale factor, the value of rI in inches,
is applied to obtain the inviscid coordinates in inches, and the abscis-
sas are also shifted as desired.
Subroutine PLATE. This subroutine is also a dummy to allow additional
calculations to be made for a flexible plate contour after the coordinates
at each Jack location have been interpolated by SPLIND and XYZU
Subroutine SCOND. This subroutine is used in BOUND, NEO, and PERFC for
parabolic differentiation of coordinates to obtain the slopes, or of
slopes and abscissas to obtain second derivatives. Three points at a
time are used to establish the parabola, and the slope is obtained at
the center point. The slopes at the first and last point are also obtained,
but with less accuracy.
Subroutine SORCE. This subroutine is used within AXIAL to give the
derivatives of velocity ratio, W, with respect to r/rI in radial flow
from Eq. (30).
Subroutine SPLIND. This subroutine computes the coefficients of cubic
equations that fit the unevenly spaced points obtained from the char-
acteristics solution. The initial and final slopes are used together
with the coordinates to determine the cubic coefficients.
Function TORIC. If the velocity gradient is known at the axial point
where M - 1, this function gives the value of radius ratio, R, which
would produce such a gradient from the transonic theory used. This
function is used in AXIAL if the option is exercised of specifying the
Mach number at point F but not specifying the value of R. It is also
used to determine the value of R for calculating streamlines other than
the contour itself.
37
A E DC-TR-78-63
Subroutine TRANS. This subroutine calculates the throat characteristic
from the transonic theory. In AXIAL, at the point where the throat
characteristic intersects the axis, the derivatives of velocity and
Mach number are used to determine the coefficients of the polynomial
describing the axial distribution. In PERFC, the flow properties along
this key characteristic are used at the number of points specified as
one plus a submultiple of 240.
Subroutine TWIXT. This subroutine is used in PERFC and BOUND to inter-
polate the ordinate and other properties at a specified point. A four-
point Lagrangian interpolation is used with two points on either side of
the specified point.
Subroutine XYZ. This subroutine uses the cubic coefficients obtained in
SPLIND for calculating the ordinate, slope, and second derivative at
specified values of the abscissa read as inputs in the MAIN section of
the program. The points may be at even intervals in the abscissa or at
arbitrary uneven intervals. The points may be the same points as those
input to SPLIND if a comparison is desired between the derivatives so
determined and those obtained el8ewhere in the program.
7.0 SAMPLE NOZZLE DESIGN
The design of a Mach 4 axisymmetric nozzle is selected to illus-
trate use of the computer drogram. The input cards for the sample
design are given in Table 1. An axiiymmetric nozzle is specified by
leaving JD blank (JD - 0) on Card 1. Leaving SFOA blank ou Card 2
specifies that the upstream axial velocity distribution is not a fifth-
degree polynomial. Leaving FMACH blank cn Card 3 specifies that the
value of FMACH will be computed by the program, and leaving IX blank on
Card 4 specifies a cubic distribution. The computed value of FMACH is
3.0821543, which is greater than the value of BMACH specified on Card 3;
t38
4
A ED C-TF-78.63
therefore, BMACH also becomes 3.0821543. The negative value of SF
means that the inviscid exit radius of the nozzle is 12.25 in. The
value of PP means that the inflection point will be 60 in. downstream of
an arbitrary point. Leaving XC blank specifies the downstream axial
distribution will be a fourth-degree polynomial, and the positive value
of IN on Card 4 specifies a Mach number distribution. The values of MT,
NT, MD, ND, NF, and LR determine the number of points on the key char-
acteristics and are all odd numbers because each includes both end
points of each distribution which is divided into an even number of
increments. The negative value of NF specifies the contour points to be
smoothed according to Card 5, and the negative value of LR specifies
that the transonic distribution be printed as the first page of the
sample output. The NX value of 13 specifies the spacing of the axial
points between points I and E to be close together near Point I with the
last increment about 3.17 times as large as the first increment,
(201.3 191.3). The JC value of 10 specifies that every 10th left-
running characteristic will be printed for the upstream contour together
with the right-running characteristic through Point E. The smoothing
incegers on Card 5 are used to control the emoothing subroutine.
Table I. Input Cards for Sample Design
CARD IITLE JU
MACH4
CARD 2GAM AN ZO RO VISC VISM SFOA XAL144 1716s563 1. 0,896 2.26968E-8 198,7a 10000
CARO 38.67 60 3. 4. -12.25 60.
FTAD RC FMACH BMACH CNC SF PP MC
CARD 4MT NT Ix IN IG MO NO NF MP MG dA Jx JC IT LR NX41 21 10 41 49 -61 1 10 -21 13
CARD 9NOUP NPCT NODO
5S 85 S,)
CARO 6PPQ TO TwT IWAT GFUN ALPH IHT IR I1 LV200. 16J8* 900. 540. .36 1 5
CARD 7XSt XLOW XENO XINC Hi XMID XINC2 CN1000. 46. 172. 2.
39
-.
A EDC-TR-78-63
For the boundary-layer calculations for stagnation conditions of
200 psia and 1638R, the value of QFUN of 0.38 overrides the specified
throat temperature of 900R and produces the throat temperature of 866R
as indicated on the output. Leaving ALPH blank causes the temperature
distribution in the boundary to be parabolic for both the calculation of
the boundary-layer parameters and the calculation of the reference
temperature. Leaving IHT blank causes the longitudinal distribution of
wall temperature to vary as a square-root function of the area ratio
* corresponding to the local Mach number; m a 1/2 In Eq. (90). ILeavi.ng IR
:' blank causes the transformation from incompressible to compressible
values of skin friction coefficient to be calculated using a modified
Spalding-Chi reference temperature and a Van Driest reference Reyuolds
number. Specifying ID - I takes into account that the boundary-layer
thickness is not negligible relative to the radius of the inviscid core,
* and its positive value causes the boundary-layer calculations to be
printed for the first and last iteration; the number of iterations is
specified by the absolute value of LV (LV m 5 fcr the example).
For the final coordinates, interpolated at even intervals, speci-
fying XST - 1,000 (the same value as XBL on Card 2) keeps the X-coordinates consistent with the location of the inviscid inflection
point at 60 in. downstream of an arbitrary point.
The main parameters selected for the sample problem were the inflec-
tion angle, the curvature ratio, and the Mach number at the point B.
The selected values of 8.67 deg, 6, and 3.0821543 (computed), respectively,
are not necessarily optimum but result in a nozzle with an upstream
length of about 14 in. from the throat to the inflection point, a
length of about 31 in. from the inflection point to point 3 (see Fig. 3),
and nearly 120 in. from the inflection point to the theoretical end
of the nozzle. Such dcwnstream lengths are probably conservative and
could be reduced to some degree although experience with Mach 4 axisym-
metric nozzles is very limited.40
.:•,.40
"AEDC-TR-78-63
The number of points used on the key characteristics should be con-
sistent with the number of points used in the axial distributions in
order that the individual nets in the characteristics network should not
become too elongated (e.g., see Fig. 7). The spacing of the points on
the final contour should also progress in an orderly manner. Several
trials may be necessary to optimize the various inputs to the program.
8.0 SUMMARY
A method and computer program have been presented for the aero-
dynamic design of planar and axisymmetric supersonic wind tunnel noz-
zles., The method uses the well-known analytical solution for radial
source flow and connects this radial flow region to the throat and test
section regions via the method of characteristics. Continuous curvature
over the entire contour is attained by specifying polynomial distribut-
ions of the centerline velocity or Mach number and matching various
derivatives of these polynomials at the extremities of the radial flow
region, the test section, and a throit characteristic. The inviscid
contour is obtained by initiating characteristics outward from the
centerline and then integrating the mass flux along these character-
istics to compute the inviscid nozzle boundary. The final wall contour
is then obtained by adding to the inviscid coordinates a boundary-
layer correction based on displacement thickness computed by integrating
the von KArmAn momentum equation. To illustrate the method, a sample
design calculation was presented along with the associated input and
output data. A listing of the computer program and an input descrip-
tion are included.
REFERENCES
1. Prandtl, L., and Busemann, A. "Nahrungsverfahren zur zeichnerischen
Ermittlung von ebenen Stromungen mit uberschall Geschwindigkeit."
Stodola Festschrift. Zurich: Orell Susli, 1921.
41
rAEDC-TR-78-63
2. Foel2ch, K. "A New Method of Designing Two Dimensional Laval
Nozzles for a Parallel and Uniform Jet." Report NA-46-235-1,
North American Aviation, Inc., Downey, California, March 1946.
3. Riise, Harold N, "Flexible-Plate Nozzle Lesign for Two,.Dimensional
Supersonic Wind Tunnels." Jet Propulsion Laboratory Report"No. 20-74, California Institute of Technology, June 1954.
4. Kenney, J. T. and Webb, L. M. "A Summary of the Techniques of
Variable Mach Number Supersonic Wind Tunnel Nozzle Design."
AGARDograph 3, October 1954.
5. Sivells, J. C. "Analytic Determinatiton of Two-Dimensional Super-
sonic Nozzle Contours Having Continuous Curvature."
AEDC-TR-56-11 (AD-88606), July 1956.
6. Owen, J. M. and Sherman, F. S. Fluid Flow and Heat Transfer at
Low Pressures and Temperatures: "Design and Testing of a
Mach 4 Axially Symmetric Nozzle for Rarefied Gas Flows."
Rept. HE-150-104, July 1952, University of California,
Institute of Engineering Research, Berkeley, California.
7. Beckwith, I. E., Ridyard, H. W., and Cromer, N. "The Aerodynamic
Design of High Mach Number Nozzles Utilizing Axisymmetric Flow
with Application to a Nozzle of Square Test Section."
NACA TN 2711, June 1952.
8. Cronvich, L. L. "A Numerical-Graphical Method of Characteristics
for Axially Symmetric Isentropic Flow." Journal of the Aero-
nautical. Sciences, Vol. 15, No. 3, March 1948, pp. 155-162.
9. Foelach, K. "The Analytical Design of an Axially Symmetric
Laval Nozzle for a Parallel and Uniform Jet." Journal of
the Aeronautical Sciences, Vol. 16, No. 3, March 1949, pp.
161-166, 188.
42"i•' ' 4
AEDC-TR -78-83
,* 10. Yu, Y. N. "A Summary of Design Techniques for Axisymmetric
*• Hypersonic Wind Tunnels." AGARDograph 35, November 1958.
11. Cresci, R. J. "Tabulation of Coordinates for Hypersonic Axisym-
nmetric Nozzles Part I - Analysis and Coordinates for Test
Section Mach Numbers of 8, 12, and 20." WADD-TN-58-300,
Wright Air Development Center, Dayton, Ohio, October 1958.
12. Sivells, J. C. "Aerodynamic Design of Axisymmetric Hypersonic
Wind-Tunnel Nozzles." Journal of Spacecraft and Rockets,
Vol. 7, No. 11, Nov. 1970, pp, 1292-1299.
13. Hall, I. M. "Transonic Flow in Two-Dimensional and Axially-
* Symmetric Nozzles." The Quarterly Journal of Mechanics
and Applied Mathematics, Vol. 15, Pt. 4, November 1962,
pp. 487-508.
14. Kliegel, J. R. and Levine, J. N. "Transonic Flow in Small
Throat Radius of Curvature Nozzles." AIAA Journal, Vol. 7,
No. 7, July 1969, pp. 1375-1378.
15. May, R. J., Thompson, H. D., and Hoffman, J. D. "Comparison
of Transonic Flow Solutions in C-D Nozzles." AFAPL-TR-
74-110, October 1974.
16. Edenfield, E. E. "Contoured Nozzle Design and Evaluation for
Hotshot Wind Tunnels." AIAA Paper 68-369, San Francisco,
California, April 1968.
17. Moger, W. C. and Ramsay, D. B. "Supersonic Axisymmetric Nozzle
Design by Mass Flow Techniques Utilizing a Digital Computer."
AEDC-TDR-64-110 (AD-601589), June 1964.
AEOC.TR 78.63
18. Spalding, D. B. and Chi, S. W. "The Drag of a Compressible Turbulent
Boundary Layer on a Smooth Flat Plate With and Without Heat
Transfer." Journal of Fluid Mechanics, Vol. 18, Part 1,
January 1964, pp. 117-143.
19. Van Driest, E. R. "The Problem of Aerodynamic Heating."
Aeronautical Engineering Review, Vol. 15, No. 10, October
1956, pp. 26-41.
20. Sivells, J. C. "Calculation of the Boundary-Layer Growth in a
Ludwieg Tube." AEDC-TR-75-118 (AD-A018630), December 1975.
21. Coles, D. E. "The Young Person's Guide to the Data." Proceedings
AFOSR-IFP-Stanford 1968 Conference on Turbulent Boundary Layer
Prediction. Vol. II, Edited by D. E. Coles and E. A. Hirst.
22. Wieghardt, K. and Tillmann, W. Zur Turbulenten Reibungsschicht
bei Druckanstieg. Z.W.B., K.W.I., U&M6617, 1944, translated
as "On the Turbulent Friction Layer for Rising Pressure."
NACA-TM-1314, 1951.
23. Coles, D. E. "The Turbulent Boundary Layer in a Compressible
Fluid." RAND Corporation Report R-403-PR, September 1962.
AEDC-TH-78-63
APPENDIX ATRANSONIC EQUATIONS
When Eq. (5) is substituted into Eqs. (2),(3) and (4), Eq. (2)
can be written as:1 GR GS
u i 3-- cy)S -S2 S3 "S GT
+ x(i -aF + GT +Gv'' x,8S 2 /
+2X2 -2 V + )- + 33 K ++,-/-(1 - -• .- ...
2 S 3y2 y4 U Y2 U y6 U y4 +U 2 y
2+ 2+ 4 + 4 ... U2 + 63 . 43 + 23
2S _2__3
S y42(2+ UxP2 Y U0 PO )
222
+ 3 - (10 - 30)y\+2 \ 4S + A i
where the coefficients are written in the terminology of the program
and x and y are normalized with respect to yo. For planar flow,
GR - (15 - y)/270 (A-2)
US (782 2 + 3507 y + 7767)/272160 (A-3)
GT - (134 y2 + 429 y + 123)/4320 (A-4)
av - 5 y/18 (A-5)
GK - (2y 2 - 33y + 9)/24 (A-6)
U4 2 - (y + 6)118 (A-7)
U2 2 - y/9 (A-8)
45
- ,j .1,, i
AEDC-TR-78.63
63 = (362 y2 + 1449 ,y + 3177)/12960 (A-9)U43 y
43 ' (194 y + 549 y - 63)/2592 (A-10)
U2 3 (854 yL + 807 y + 279)/12960 (A-11)
" Up2 (26 y + 27 y + 237)/288 (A-12)
2LU0 = (26 y + 51 y - 27)/144 (A-13)
For axisymmetric flow,
GR - (15 - 10 y)/288 (A-14)
GS - (2708 y 2 + 2079 y + 2115)/82944 (A-15)
GT w (92 y2 + 180 y - 9)/1152 (A-16)
GV - (y + 0)/8 (A-17)
* GK - (4 y2 - 57 y + 27)/48 (A-18)
U42 (2 y + 9)/24 (A-19)
U2 2 (4 y + 3)/24 (A-20)
U6 3 - (556 y + 1737 y + 3069)/10368 (A-21)
U4 3 -(388 y2 + 777 y + 153)/2304 (A-22)
U Y2U2 3 (304 y + 255 y - 54)/1728 (A-23)
P2 - (52 + 51 y + 327)/384 (A-24)
Up0 - (52 y2 + 75 y - 9)/192 (A-25)
The first part of Eq. (A-i), which is independent of y, can be recognized
as Eq. (11) for planar flow or Eq. (12) for axisymmetric flow inasmuch
as x and y are normalized here with the value of yo.
46
own'!
--- ... :-• •ig lo o ") i:- ,, _ T3•:I• r"l~:IIr•I ]1'
AEDC-TR-78-63
In a similar manner, Eq. (3) can be written as
4 2
2-) 42 V22 +V 0 2NS• 2(3 •)s + s2
6 v y4 +V 2+ V6 3 - 4 3 +v 2 3 y -V 0 33'
2+ x I (2y +- 12 - 3+)y 2y 1 .5q
(9 - 3o)S
6U6 3 y4 -4 U4 3 + 2 U2 3 +,
S2
( 4 U2 '
X12x2 (2 Up PO+
2 S -2
+ x3 "4 y .. )/ (A-26)
For planar flow,
V4 2 - (22 y + 75)/360 (A-27)
V22 - (10 y + 15)/108 (A-28)
V0 2 - (34 y - 75)/1080 (A-29)
V6 3 - (6574 y + 26481 y + 40C59)/181440 (A..30)
" V4 3 - (2254 y2 + 6153 y + 2979)/25920 (A-,31)
2
V2 3 - (5026 y2 + 7551 y - 4923)/77760 (A-32) '
-03 (7570 y + 3087 Y + 23157)/544320 (A-33)
47 ,
44'i
IIIi: t A EDC-TR.78.•03
L, , t
, For axisymmetric flow,
v4 2 " (Y + 3)/9
v,22 (20 y + 27)/96(A-35)
V0 2 " (28 y - 15)/288 (A-36)
V6 3 - (6836 y2 + 23031 Y + 30627)/82944 (A-37)
V4 3 " (3380 y2 + 7551 Y + 3771)/13824 (A-38)
V2 3 - (3424 y2 + 4071 y - 972)/13824 (A-39)
"0O3 - (7100 Y2 + 2151 Y + 2169)/82944 (A-40)
48
* I , ,
48i
A EDC-TR -78-63
APPENDIX :CUBIC INTEGRATION FACTORS
If a curve through four points with ordinates a, b, c, and d,spaced at uneven increments in abscissa, s, t, and u, is definedby a cubic equation, the area under each section of the curve canbe found in the following manner:
Area b F a+ b-+ Fc +dd (B-i)a-b as Fbs cs Fds
Areab-c at a + Fbt b + F c + d (B-2)F Aresc-d "F a + Fbu b + F c +F d (B-3)
c- u u cu Fdu -3li Itota G a Gb + G c + Gd d (B-4)
where
2S(s + 4t (+-5)
•",2 I Zs + 4t + 2u
iij a + 2SFbs +--)(B-6)bs 2 12 t(t + u)
a 3 (s + 2t +(2B)cs T2tu(a + t) 's- 7)
I:: s_(L±+ 2t,)Fds 12 (s + t + u)(t + u)u (B-B)
t 3t+2u (-9)at 12s(s + t)(s + t + u)F - 2+ t + 2u
(bt 2 12s(t + u)
t 2 (2s + t -2B1)F " ct 2-u(s + t) (B-11)
49
_ _....__ _ _ _ _ _ __._ _ %.-- - .,.
AEDC-TR-78-63
t 3 (2s + t) (B-12)Fdt -- 12u(t + u)(s + t + U)
F u u3(2t + u)F = 32 )(B-1 3)
au 12s(s + M)s + t + U) '3iu3(2s + 2t +_y (B14
rbu - 12st(t + u) (B-14)
u (2s + 4t + u) (-5cu 2 12t(s + t)
..R u2 (2s + 4t + 3u) (B-16)du 2 12(t + u)(s + t + u)
G -F + Fat + F (B-17)
Gb bs + Fbt + Fbu (1-18)
Gc F + + F u (B-19)
G -F + F + F (B-20)d ds dt du
If all increments are equal, then
s t - u h (B-21)
Fds -F at F -dt -Fau -h/24 (B-22)
F -F-5h/24 (B-23)Co •bu
50
AEDC-TR-78-63
F = F - 19h/24 (B-24)
Fas Fdu -9h/24 (B-25)
F bt Fct - 13h/24 (B-26)
Gt 3h/
d 3h/8 (B-27)
G - c G 9h/8 (B-28)
The values of G's in Eq. (96) correspond to those in Eq. (B-4).The value of F's in Eq. (97) correspond to those in Eq. (B-i).
* '5 1
2 2$WA,,
j'""
AEDC-TR-78 -3
APPENDIX CINPUT DATA CARDS
Input Columns
Card 1
ITLE 2-12 Title
JD 14-15 Blank (0) for axisymmetric contour,-1 for planar.
Lard 2
GAM 1-10 Specific heat ratio.
AR 11-20 Gas constant, ft2/sec2 R.
ZO 21-30 Compressibility factor for an axisym-metric nozzle, constant for entirecontour. Or, for a planar nozzle, ZOis half the distance (in.) between theparallel walls, and the compressibilityfactor is one.
RO 31-40 Turbulent boundary-layer recovery factor.
VISC 41-50 Constant in viscosity law.
V1SM 51-60 Constant in viscosity law. If VISM isequal to or less than one,
VISC* T lb-sec/ft 2
If VISM is greater than one,
1.5VISC*T 2"
11 VI- T- lb-sec/ft2 . IfS"T + VISM
T is greater than VISM,
VISC* T ; T <VISM.2 VISM1/
2 -
SFOA 61-70 Used for nozzle with radial flow regionif 5th-deg axial velocity distributionis desired. If positive, the distance,in inches, from the throat to Point A
52
'- .~-
AEDC.TR-78-63
on the characteristic diagram. If nega-tive, absolute value is distance from thethroat to Point G. If Blank, 3rd- or 4th-deg distribution iu used depending on valueof IX on Card 4.
XBL 71-80 Station (in.) where interpolation isdesired (e.g., the end of a truncatednozzle). If XBL1000., the spline fit
subroutines are used to obtain values atincrements evenly spaced in length.
Card 3
ETAD 1-10 Inflection angle in degrees if radial flowregion is desired. Two characteristicsolutions are obtained, one upstream andone downstream of Point A. If ETAD = 60.,the entire centerline velocity distributionis specified and only one solution isobtained and the inflection point must beinterpolated. If ETAD - 60., IQ 1 1, IX - 0,on Card 4.
RC 11-20 Ratio of throat radius of curvature tothroat radius. Must be given if ETAD - 60.or FMACH - 0. If FMACH is given, RC iscalculated. If LR - u, IX = 0 givws third-deg equation betweev Mach 1 and EMACH,matching first and second derivations ateach end. If LR 0 0, the value of RC foundfor LR = 0 is used with given value of FMACHto define a fourth-deg equation. If IX = ±1and FMACH is given, RC is calculated todefine a fourth-deg equation. If LR ý 0,a new value of FMACH is found, compatiblewith the value of RC calculated for LR = 0.
FMACH 21-30 Mach number at Point F if ETAD # 60. Nega-.tive value specifies Prandtl-Meyee angleat Point F as IFMACHI *ETAD (usually around-7). If PMACH and RC are given, IX = 0and 4th-deg distribution is used. IfFMACH - 0 and IX = 0, a 3rd-deg distribu-tion is used. If FMACH = 0, and IX = ±i,a 4th-deg distribution is used. FMACH is
calculated if not given. If ETAD = 60.,Point F is not defined.
53
. . . .
•'•'• •'• • '••'• .... .. ' • •±• ,"• •P, --- • "-- --'; '• " • .. ... "• . .............. .. ........ ............. ..
AEDOCTR-78-63
BMACH 31-40 Mach No. at Point B if ETAD • 60.
(MC 41-50 Absolute value is design Mach No. at PointC. If ETAD 0 60, positive CMC giyes d 2M/dx 2
-0, and negative CMC gives d2 M/dx' 0 0. If jETAD - 60., CMC ts positive.
SF 51-60 Scale factor by which nondimension coordi-nates are multiplied to give dimensionsin inches. If SF - 0, nozzle will have aninviscid throat radius (or half-height) ofI in. If negative, nozzle will have aninviscid exit radius (or half-height) ofISF1 in.
PP 61-70 Station (in.) at Point A. PP - 0 givescoordinates relative to geometric throat.Negative PP gives coordinates relative tosource or radial flow (ETAD 0 60.).
XC 71-80 Nondimensional distance from source toPoint C. XC 1 1. requires centerline MachNo. distribution from Point B to Point Cto be read in as input data on Unit 9.Otherwise, positive XC gives 5th.-dog dis-tribution if CMC positive and 4th-deg if CMCnegative. XC - 0 gives 4th-deg dietributionif CMC positive and 3rd-deg if CMC negative.Negative XC and IN gives 3rd-deg distributionwith d2W/dx 2 not matching source flow atPoint B. If ETAD - 60. and XC > 1, XC is ratioof length, from throat to Point C, to throatheight. Negative XC gives 3rd-deg distribu-tion in M; XC - 0 gives 4th-deg distribution;XC > 1 gives 5th-deg distribution. XC - 1.requires centerline Mach No. distribution tobe read in as input data on Unit 9.
Card4
MT 1-5 Number of points on charact',ristis EG ifETAD ' 60. or CD if ETAD - 60. Maximumvalue about 125. Use odd nuwjer. A zero ornegative value stops calculation after tonter-line distribution is calculated if NT positive.
54
=l . . .. .. i I I III . •ll lN.• . . . . . f , ,
AED-TR - -63AEOC T-78*8
NT 6-10 Number of points on axis IE. Maximum valueis 149-LR. Use odd number. A zero or nega-tive value stops calculation before center-line distribution is calculated but afterparameters and coefficients of distributionare calculated.
IX 11-15 Determines if third derivative of velocitydistribution is matched. IX - 1 matchesthird derivative with transonic solution.IX w -1 matches third derivative with sourceflow value. IX - 0 does not match third deriva-tive but gives constant third derivative ifRC - 0 or FMACH - 0.
IN 16-20 Determines type of distribution from Point Bto Point C, positive for Mach No. distribution,negative for velocity distribution. IN - 0 forthroat only. If XC is greater than 1., thedownstream value of the second derivative atPoint B is 0,1* JINJ times the radial flowvalue. Similarly, If ETAD - 60., the secondderivative at Point I is 0.1i*IN times the"transonic value.
IQ 21-25 Zero for a complete contour if ETAD 0 60., 1 forthroat only or if ETAD - 60., -1 for downstreamonly.
MD 26-30 Number of points on characteris'ic AB. Maxi-mum value about 125. Use odd number. A zeroor negative value stops calculation similarlyto MW.
ND 31-35 Number of points on axis BC. M~aximum value is150. A zero or negative value acts like NT.
NY 36-40 Absolute value is number of points on character-istic CD for ETAD , 60. Maximum value is 149or 200 - ND - MP - jMQj - number of points onupstreem contour. Negative value calls forsmoothing subroutine.
MP 41-45 Number of points on conical section GA ifFMACH 0 BMACH. Use lialue to give desired in-crements in contour -, usually not known forinitial calculation.
55
* . |-- . ............. . . . . . .
A AEDC-TA-78-63
* I MQ 46-50 Number of points downstream of Point D ifparallel inviscid contour desired. A nega-tive value can be used to eliminate theinviscid printout.
rB 51-55 Positive number if boundary-layer calcula-tion ir desired before spline fit. Nega-tive number transfers control of program toJX. Absolute values greater than one areused to approximately halve the number ofpoints on the upstream contour even thoughLR + NT - I points are calculated from char-
* Iacteristic network if LR > 2, or (NT + 1)points if LR a 0.
JX 56-60 Positive number calls for calculation of stream-lines, zero calls for repeat of inviscid calcula-tions requiring new cards 3 and 4, or, ifXBL 1 1000., for spline fit after inviscid calcu-lation, negative number calls for repeat of cal-culations requiring new cards 1, 2, 3, and 4.
JC 61-65 If not zero, calls for printout of intermediatecharacteristics within upstream contour if JCis positive and downstream contour if JC isnegative. Characteristics are (NT - 1)/JC or(ND - 1)/(-JC). Opposite running characteristicthrough Point E (or B) is also printed.
IT 66-70 Number of points at which spline fit is desiredif points are not evenly spaced, such as Jacklocations for a flexible plate. Used only fora planar nozzle, inasmuch as a nonzero valuecalculates distance along curved plate surface.Positive value of IT tequires additional cardsto be read in (8 points per card) after boundarylayer is calculated.
LR 71-75 Absolute value is number of points on throatcharacteristic used in characteristics solution.Negative values give printout of transonic solu-tion. LR - 0 gives M - 1 at Point I.
NX 76-80 Number from 10 to 20 determines spacing of pcintson axis for upstream contour. NX w 10 giveslinear spacing. NX > 10 gives closer spacing ofpoints at upstream end than at downstream end.NX - 0 same as NX 20. Ratio of downstream
56
* , I I
A E DC-TR-78-63
increment to upstream increment is (NT - 1 )NX/10 -
-NX/1O(NT - 2) . Optimum values, usually 13 to 15,determined by trial and error for specific con-tour desired. Negative NX used with negative LRlimits printout to transonic solution.
NOTE: A zero value of MT, NT, MD, or ND will allow a repeat of cal-culations for parameters specified by new cards Nos. 3 and 4.A negative value will allow a repeat of calculations for newcards Nos. 1, 2, 3, and 4.
Card 5
NOUP 1-5 If smoothing is desired, negative NF. Number oftimes upstream contour is smoothed.
NPCT 6-10 Smoothing factor in percent. Smoothing factor- NPCT/100.
NODO 11-15 Number of times doimstream contour is smoothed.
Card 5 If boundary-layer calculation is desired usingorinviscd points calculated from characteristicsor solution. (No smoothing). '
Card 6 If boundary-layer calculation is desired usingevenly spaced points interpolated from spline
or fit of points from characteristics solution.
Card 7 If boundary-layer calculation is desired usingevenly spaced points interpolated from splinefit of smoothed points.
4 PPQ 1-10 Stagnation pressure (psia).
TO 11-20 Stagnation temperature, Rankine.
TWT 21-30 Throat wall temperature, Rankine, if QFUN - 0.If TWT - 0, the wall temperature is assumed tobe the adiabatic value.
TWAT 31-40 Wall temperature, Rankine, at Point D. Forwater-cooled wall, the bulk water temperaturei.s assumed to be 150 lower than specifiedTWAT. The cooled wall temperature distribu-tiou is assumed to be
57
t=,.. .. ...... . !. j. • "u•' ..........
I
AEDC-TR-78-63
(TWT C WT) "' A* 1)TWTWAT +
where A/A* is the area ratio corresponding tolocal value of Mach number and Ac refevs toPoint C.
For negative IHT
TW - TWAT + ,TWT-IA x A A*
QFUN 41-50 Heat-transfer function at the throat.
QFUN - ha(Taw - TWT)TWT - TWAT + 15
where ha has dimensions of BTU/sec/sq ft/R andis obtained by Reynolds analogy from the skin-friction coefficient. If QFUth is specified,input value of TWT is ignored and TWT is calcu-lated from QFUN.
ALPH 51-60 Parameter specifying temperature distribution inboundary layer. ALPH - I. uses quadratic dis-tribution both in the calculazion of the refer-ence temperature TP and the calculation ofboundary-layer shape parameters. ALPH w 0 usasparabolic distribution in boti. calculatiovs.ALPH -u -1. uses quadratic distribution for TPand parabolic in the calculation of boundary-layer shape parameters. Within bcundary layer,
T - Tw + a(Taw- TW) (U/11)
+ Tes- (Taw- Tw) - Twj (U/Ue) 2
where a 1 for quadratic dist.
a -0 for parabolic dist.
IHT 61-65 Integer which determines temperature distribu-tion (see TWAT). If nonzero, IRiT determineshow often subroutine HFAT is called. An absolutevalue of IHT greater than KO, the number of points"on the upstream contour, will prevent HEAT frombeing called but will allow the choice of tempera-ture distribution to be made.
58
11 p 1 Ml 1 11.01
AE DC.TR -78-63
NOTE: HEAT is a special purpose subroutine for determ.ningheat-transfer values for the upstream contour. Thesubroutine HEAT incorporated in this program is adummy.
66-70 Integer, parameter specifying transformationfrom incompressible to compressible values.If YR w 2, Coles' transformation is used forCf and Re If IR - 1, TP is calculated bya modification of the Spalding-Chi (Van Driest)method. If IR - 0, the Van Driest value ofRe6 is used, but if IR- -1, Colas' lav: of
corresponding stations is used.
Cf Cf * TE/TP, Re - RD*Ree,i f
ID 71-75 Integer. If ID - ±1, axisyminetric effectsare included in momentum equation and in cal-culation of boundary-layer parameters (6 nornegligible relative to coordinate normal toaxis). If ID - 0, these effects are omitted.Negative ID suppresses the printout of theboundary-layev calculations.
LV 76-80 Integer. Absolute value, usually 5, deter-mines number of rimes boundary-layer solutionis iLerated so that radxus terms in momentumequation refer to viscid radius instead of Jinviscid radius. Value of 0 or absolute val-eof 1 uses inviscid radius, Positive LV repeatsboundary-layer calculations for new set ofparameters on a new card if XBL j 1000.
Card 5 If streariflines are desired, JX positive, (Nosmnoo thing.)
METAD 1-10 Inflection angle in degrees for streamlinedesired if ETAD 0 60. for Card 3. If ETkD -60. on Card 3, use ETAD - 60 on this card.
* 11-20 Fract ion of r'.ottour desir.-ed if ETAD - 60.Otherwise, QON - lTAD (in Card 5 divided byETAD oa Card 3.
X,1 21-30 Value to upd•ite JX for subr.eq4oent calcula-tion, JX XJ.
,, 2 !*.;.a a'11d598Y '
I,
AZ-DC.TR-78.63
Card 5 If SPLIND used after inviscid calculation* (JX zero or negative and JB zero or nega-
or tive). (No smoothing.)
Card 6 If SPLIND usad after viscid contouir (,7B posi-tive and LV zero or negative). No v'.oothingof inviscid contour. Or, if inviscid contour
or is smoothed before SPLIND is used.
Card 7 If inviscid contour is smoothed, boundary layeris added and SPLINE Ls desired.
XST 1-10 Station (in,) for throat value of X. IfXST a 1000., program uses value previouslydetermined by specifying PP on Card 3. Other-Wise, value of XST is used to shift contourpoints by desired increments for arbitrmryStation 0.
XLOW 11-20 Starting value for interpolation. Second valueof interpolated X a XLOW + XINC.
XEND 21-30 End value for interpolation. If zeto, SPLINDis used to calculate slcpe and d2y/dx2 at samepoints as previously defined.
XINC 31-40 Increment in X for interpolation. If zero, andBJ > 10, contour is divided into BJ increments.
BJ 41-50 Value to update JB for subsequent calculation.JB - BJ. If negative and XEND - 0, interpola-tion is made at diecreate points read in on sub-sequent cards similar to case when IT > 0.
XMID 51-60 Intermediate value for interpolaticn. Distance(XMID-XLOW) is divided into increments definedby XINC, and distance (XEND-XMID) is dividedinto increments defined by XINC2.
XINC2 61-70 Inctements in X between XMID and XEND if differ-ent than XINC.
CN /1-80 Number of copies desired of final tabulation ofcoordinates if more than one copy is desired.
ý,j
I,..60
........................................... ,
AEDC-TR-78-63
lz zx xxzsxxzzxzslzzxzzzz.0p-UeoxxzzRdgzzz*A
W2eV% ;1~M~ - t5-
n n n n n n n n n n n n n nnV
W 4;
a,~Iris. O-cc j 9 .
0 1 ~ f @0. i 6
S 0 IA
as 0.4'W..M o. *1 EM 4MR3LO
~~~~~ N @@ c ~
4*1; :af WhiP). UhZ
9X ts .0
4~ a.
wi M iE C *O40 03
-- 0 T-R t
*0 . * 4T% 4 "I~ SeI
ru4-.J M-0 V"OO .0%.J ** m
vv vhuIwwu 0 Uu 0uJ%~cV .
~ 0 I * o 04 .. 0 ).%.kU.61I
AEOC-TR-78-63
zzzzzzzzxgzzzz zzzxxzugzzzzzzzgxxzzZZZZXzzzzzzzzxxzzz;
941-
j 4;
0 a-
* -9!!
W %..- . 0J j *.dA a0
O. J -O - - z
-- -- i xw- moo1 ..WIMa po 34 IIAC U
tv 5- k"I *4 *Jj - hIg.* -
624*1
c.A Trrzt = 77-
AEDC-TR-78-63
-AIzz dN*Jfli NNNNN 0zzmaff zPzEE4E4EEE 4 * zzx 44444444
Z11
C>.0>.g
I I-
N - -0 - ,o
I.a* 4 1 ;.
&A 0- 3.pc 0 -
v . 0 >. ILm 0 jz IL - z 0 ccza.'
No 1. 45 N >. 0 4L w>*0 3 a 2e Ax x 0. 0 _;
.I0q ý -- 0K 400 3r_j 2 .,C IL CL 4LJ, O z- u - x w- V %
It 0i;2 01 4 40 co 0*ýO-0 12 Lo0* - 3
aJ W42 0 0 ). ~.O0204 U JI 0 3 - 4. 0 *- - - .. *. * 4 3 3 0 3. ~ ) . *02gag;; jz v 0 2 - 0 I C .4.-.0 ) * 0 0 1 03JIo m f 0 0 0 3Ix U 00 3
N~ L N - N N 0 - - - ~ - 0 U . 3N 3 N- . hi *** OJ. .2 ~ **- 0.0 0~0 2)0 . 4 4- .- C.. co-*~ ~~~~~~~~ u- m 0-iJ-4 x0 ' 04 44 .. 1 *404, x 4 0 0 -0--t
a 0 P- 4.P. I W ~ * * 32 2 2 *- U*. * W * . . W *
3- *3-4 N *** 3-4**P P.3- . . 3-0 0 3 . - 3 * . 3 3 . . 2
AEDC-TR-78-63
J 4;
f4 Wo
ICc o iNa w N
owh Ni0Jcn; 0
44
"o OJ - VI -
hi of 0 Im 1! WC 2oW . I C; ý x9C
04~~ 00-
mu -- 3- VVo - --- wm4oo w u o
* K 4a*o=
.01 M00 KyN N N
~ ~ * 64
A ED C-T R-78-63
S4~S P. @ .. N~ ~ 4 0. N 44 @0 NM A 44 04 S 4O 0 ~ E *I ,.j
9:! ~ ~ N M~~ M M iS E..4 4I, ~~ 1U
44 4 44 4 44 4 44 4 4 44 4 44 4 44 4 4 44 4 44 4 44 4 4
4 W UO w !x= -
z *v
VS 4 Z* . 4 :-W
ow20 0 UO Z S*-tU P-O-6 40.4 *o%
CL t 24.4.. Z
S. a 4. % % 0 0 ;T 3000.
,z4 I N -> - 0 *C4 * 3 4 . O ~ S h 0 A4 Ai .0 0
- ~ W 00. V -4 2 4 N 0- .00 -sO
4 P4% %t b2 "'* -*.000Z - w on
U. ~ ~ ~ ~ C 000.h. ". S0.AS f~ -a~~ ~~~ * N. N0 :00a44*I,.
X. S . Z.. s I NO 4 - T ;, -5 9* **IlSUzNl-'. 4 h h . ~ 4 0
WO O.P..K .4.. U..Z .1 4nn n A % ~ g..- -
4 0.0.KznW...0 2 4 0 - 5; - ý U;0.Z0 N-0* 0
US U.0.. U .0.. OO MO *00 US4KS00t.JS- 000*
- 0* U i. .. e N* Z 4 - i S S2. - - 0
0 0 -VS 4.P~S. .4 K. 0 - 5. 5 5 4 SOK *K65*
AE DC-TR-78-63
I-.A.90.Ncome*00. E30:- 0- 0- - - - -e- ---
a 0 ~400y
Im 0.e P4-
to" *NM 0
0 000 w%% 10 v-: w x I
0.~~ 0d 0 UUH.0 * cc 4 C055 4
*Eý "I *;,- x z
V F" f0Ad 00 * f .P 2 -i 2ONO -0fl4 * .4 u z 06.19
N@* e i.n toON .. 4... . S . 0
%11 .. SA 0 *0@N. * @ Iwo v *
-0 *1rP 0 0 4 % * 4 . tI 4- twoM ID.5S- 0 -0 *S & U 19 As
.0%.*US0. ; S4 0-4INCI I- *0.- 1CO 0%401 U*5S z Z Z. WO ).
to. OPA* 0-0* 000 Uu wZO *UW -00A *0 01 f 1 0 0
-* 10 M0t1-. 4 o06 0 c4 C!**N !-j -1 ** 54 0- 1 2
-- 00 fu N*4N Do- A0 0 U 0 w .5
*~1 *0 OTPSPSR mg4 * ~ u.
EMS 00 4 . 554* -. 0 4 V Z"--WO vuIC -=es, me... zo::4A U 0 o o a iUCI
*4*9 9.0 w~~X 06 PU u x. 0. 4 * J4
04x0 30 0 C NO S *D N & x a u & J- ZW- 0 2m 2u M K N -2.- S
5-00vv vu...v uv v*@ 1 4 N S U J W.U
*PS. 040 - -C i U . 5 1 UUS664
A ED C-T R-78-63
00
z -0
z w u =-- w
o1 ur -a
2 4 hi.-0 31 b
-Il0 a. We%
-4t S 2IE - 1-w 4 4 z -s zza -o~04 0 i hi P- .a. ; - w W 0 U 0 ;;Va 0411 0 21 Z
0.-n- a. rn-
hi wo 0 'a0 00 a uo w *
rn- Q -O U - -z:- 0-za o ;;. *- u -;;ih v KLU-tU I-z 0.
IN- OM1XE * ---- i a.!M z mz z h- 0U-UX3 2 0.LI *IO 4U- hi0 0 ~U hi-ww w 0-hi Oh hi
o 411 -0~ 4- 4- 1-0 -0 VI i a.3 0 3 W Z ..67
AEDC-TR-78-63
IIaa
I.In
4I W
o- S.- &0.0a a ly U00
IL. 0o- 3l
l, 4 ':0.
a - a 2 0.I 3
C-. oz . - VIwhZ of *I Fo- * m w. I-t4La5 0 ow ':U; XIL3
od .1 - 0 w -2; . 0 I- 4.-
(D- a- ). -n x I i Sw.. m CL I- -1 IL 06I
I0 -!; ,' I~ ' . .;10. A0 o.C A '.
* ~ ~~~~~~ V-03 . 00 6 A1
* -. O - i)- W 0 468
AEDC.TR-76-63
........4t - 4* 4b 4 b.4. 1414M44-H*~ ----- H M ""H M"1 4
!Si
1(344 ~ ~ ~ ~ ~ t 4444 4 7 444444q44q4.4
i
68
.. ... ..... ....
O•'°".J.= %C • * !
•' =.• .i =. . .=a =..°l' .• . ! . fo=. .= * .. .= ,z:sc. ••? " " . ' ' a I . *,, ..5 •
ill,
-"''.. UU,- 3 -0
- ~~68 ..
A E DC-TR--78-63
WN 0 06 O 'Nm Nr
0 t
x Z-
c' Dt ft 7 ftu -4 CW C
to 0
5; z z 2 22 0 22l.zxt ft *. 0-oK Wz
0s 0h WT W IX9ft -4 4 1" w w ;; 2 JD; q i'.x Noo ftS-1- 0.uj *go. a w -4 05U 4-: 5-
-1 ft- o u . K1 Li 2 Z04 Z e w -15 0z 5- I *.W WK
U~mw -5 a -- v) fto 611444 - U --c*; 3 -o ma P.S;; 0 ;h w - 0 - 5h;h wz.5 NW N 4 wLPN1 -W1 N - 0 n - 0c US u _j- -9 -c
c zS' rY 0 -K W4 W U - 5U-ly fW 0v 5- u ry 41I -- 55 - g y
Si i- .0 0 .-- S 0 -- ft W 04W69.
AE DC-TR-78-63
. ~ -a~1 ZO 4444AWW 4 a
31
U. a.
.1 Wa. 0-1
3 I-
.j -. lz A. jI
o AA W L .
* 01.. 4W W - ;I 0 3..
06 a- I- - Iai.. P I P-o 3W W
W W t 1 .1 vQW a W 3a a t 1 0 : I1 t tJ. ti a I
* -ma xA -- U 1- A -W OWW- Z *3- 1....
a. -- 0 ell O-- - 1*--0
NO~ ~~~~ -Jl0- -- 0n0 . 34-. ~ -4 -- 00 -~- E -.0 sa.*l0
W. . W*.W Lb.0. W OW. *LM. 1- I.0. WW ILI
44 M1-4 MO O .0 IA IA4W *0
4 MW 4O- Al4Al~0 33~-3-701
A ED C-T R-78-63
*.NI4~t.0~Ae.d .~0.N*A -00~003 JM U.C.OA. M
ZO
m* -
301I I-
I- N
4s 0 'a .ILZtoIx -a
a. -w LA CL CIa. 0o o f
o- w a.* 0 4 0 0m IWO 06 iD
a. a. w www 6 $0N 4".1 a. ; ! kI.o. 6
w 0 o - 2 v Z" I*" We 3-
4 * 1 - I N : 4 54 1 . 4hI 0S - ,.-4 4 1 ; V/1 14 - 2
-"Ix. h 20 W 0s 0z ~ oxwI.. OX I- - - o n a . -2X6 v &I,0 . 41Iaa
a ; a x. -a * a. 'o 2I .0 . . . wO hISWI 3 s-s-W S UI- 2* 66 s. .- aNs
UK 31 11 b 0 uWs wJW. - 3-
o l a. 0 -. le 12 N N -N o- - ~~0h hI . 3 I hJ * 6 4 6 h..
3 - - S-J31 I s- 4664 7106 s C -. a
AEDC-TFI-78-63
ILa -it
- 2 T.j 4 ;a
-j
V; lI; 1%4 4 -IWK 2* M 0* .- F.
OhaW OUL .W4.W X c - - -
V~ -:: ct40 &a*8 0 - .Z 311 'a. W
.. h00 Se*a- 4- -0 o4mw -24 4 2 *U-hl -zz
00 .6 1d C: uu W . .W V U
-00.z N*03 3W- 14 Vto~09000 40 -aM W" z~ M qh -z W . Y A
am&&-Go VmflOiWJ &&UNM * *ih - .JW00ý -I0 I IL. 1a .oo *-" V.I -j
3~~W A. -I WU0 0@3 - A A4V z*
*W 4 Vi4 * 0.** * C S -. in Ip XX u in
ha.. %~K mO a.U ww b0.00. 0 ~ 0472
AEDC-TR-78-63
a 0
4443. I*E I
w *
.310
or. 0 .0 1.1
-iCLCL..CLCLmm w aj- owai L- 6a.a. - 31a
CL21 -W Lz %uZ Imm-- vu u.
MOM*: w " '0 -9 0-a o DO0,Uo x L .UI. 0P- .
0.43.- ow5 cN . 0.M2 -0, wU i- O Z 2Z Z
39IL -0.- -C -.. 02 -- I - z 0 -0.2- --.V.. f~l 0 *Ifom 0 - 00% %
-- a3 0 .a a .~ .U R s o 0-
0U303 vv uuav 00 00 a3 * z m o m
mm-a 32 --.4 -- x- wa Z w wU w .. OZU.VV1IDlIl ------ U4 ------, -I..0 *------- in-.x
In M - inS 0 3a-uunInE2- 0 i-S S.
AEDC-TR-7843
II
U a,
fn 0 I
emL0M0 a. aZU . 06 4
IN 30 IL
NISI a.ou~
0 ZR U a.w 00 z CL 0z 0
CL K5- U U" U CY 4.
1 VC 0U PC CC UCC
OD M 0w -- 5--- vo oKe Z T: M :VUUK:- VI- U aV
10 .5O '-e - 4: ; * 1W 0
:;z. Z c OD I . U
4 5-Uit! I.-S U - U a 5
-~~ U -'IJ K - 1 0 - hi -- v 2 .- K
x 4D -a 4 -- *OU2 A R . U
U UU4 lay 0 a-h~ U Id
4 44II. Kh.S-5A4 -00 . a *. . ~ 474
AEDC-TR-78-63
~w
* ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ' p---N---@**La5e0o*.cce 0----------
di SSSSSd i id id idid id id ifud id idid id id d id id. 00 0000000.........
4 4 44 4 444 44 44 44 4 4 4 4 4 444 4 4 44 444 44 4 t 4 4 4 4 4~(44444 44444 4444 u w K 0444 44 4 44 4 44 4 4 44 4 4
W;Z
0 0
0~~~. llZ 0 ; e
'o 4U X -
-V I: ý% -
N ~ Z 0 .0
a 0 di *ul CL
La w u a:, 0 U
S a, :! 0- a* 0~ ---1 ga 0 -
0 I0 000 sK U. -. 0
Z. 00 ON Wi o a3- U0- * 0 0 0 - 3*S *S
i.;~ ;U !w -rI 0--46. XZ11 -9 - Ui 1 4 CL0
2 0 9 Z 04 S MIA* M% K M *U~
0'21 - -ONO Ni M0 -- z; :&II 1
0 m~l *n m.0 ' - 0 U ~ 0 0 S 00 .. LA- K', - - l..W *-* l K~l. 0 .to
hjO- N 9-0-0 4 ' 40-l- 0.~ 4.K' i75
AECT-86I,
Iu
xN fu 'a2W 02
13 0, U 0. M4 4 .a.~ z x X .- ux
0. U. 20 -om M4--C
4u C X44UXO.W0. M4 - OZ fn 119 1.49- Z 119
~60. Fn U W 4. -Lr
& 7 1 -z 2~ i.0 314 Z 10 400
it a 11 2- . 2- in -4493t-
30 4 0 On .40 Vol4.j--- 4. .36 *4 11k 0 . ',
X021 w MW *O 42140me39 V. 0 a. .44 . 24 U X"NN &. 0 0- 6- -. 1 U. 2U f3 U 0
NO 0 v :,Joao.N N 0-.Z~l.N44 0. 0.-. WO N r 0 .-,3* 1
1.W . 4 0- TO11 440WW .N K 24 o32 . 4rs
xIT .03Or 2.-Z 0 4 w x* F2 X 04--..... 2 40.0. v m0 I a *. 44MX1*.30N1.. 1-
XX~ ~~ 32( 434 - U * . O - 4 I W - 2 - 4 1 9
ft a90 0. O ~ f .- n . v M U.. 11 4 IZ - .zX 0 , 'I90 - - 0I. 0 * . 9 1 . .
N@In* .040 W -4430.3 .J9IwI - 06 0 6 -- -0!2VIý4k.4N-044VIr(0 *0 .44
*~~1 11 11391 . 4 %* 41 A N -ý ON90.1.. . . .13- I3.I
ID3 394 * 0... N ... 40 .G 40- -. 9 42 qpU-.4p4.-.
1.JZ 030 J-k. -44 2VW~uw4 4 4 76
A E D C-T R -78-63
N400Nm4A4N400-NM4IA4N 0,..l 4114*400 -N oU.'%M4t40.I1
NN04 404 4tOO OOO 000 0004 000----------
.0* - E 4 j N cc N - N
-Z ~ v -o 44 4.m N9 w. 2~ -ý NI vIn 0 0 4 0. - i *4 4!. 4r P - -W :o U
vC 1- A- -C W - M, am 6 I. Z CC X NI- i 3 -3N _;wo ZZI Aj I%* 0 4 .6 2 4 N 4 4
x 2 i 4 1 .91 1 2 ju CL Iz 4r IA .1 c x tn o
4 4 2 -9 to0.1 *4 42 -C P: 0o mmt11m -9 06
4. r 4211 N W~l U, -j 6611 --- 0 It w- 0. UZ4-4~~- 0 *Ly1 a6 0Z U10 oN.2.-N-. n 4 Zol0Z-
2 42Lm X - CL 410 4j,!*- 6 -4 N w. 40 *4 Cc LA i14 4 . * -1.. .1 40 N 4- W w11 4 WO .2 x .e*
*z -U 34 -N 3 J .- .3. . * x* *6 -w '.2 N * 04 *
0. N4' Z -9 N1 *L304I 61 140 4 . .8 3'.94 -W .X 04.!x. Is 02 2 m-2 - 0 -C -C w ILN4 0 - 4. 4-.l214. 2 z-W 0 2-2 0.w42.1 ~ 2~ - . I* I f0J.i4
Ur.2 0 z - 'i '1 041L~-. 30ý *4. .11 v1. N23.4. 4
2.N ~. 0 CU -- S 4 0 -1 *N.3 ... %2 0 *.UIV N44w w 4 0U4. 011 *%Dgl0X4-. * .T-. 00..3V3
zN Z;i~~ N a, , 1 o0Ji0 0. 1 K0 *.2 *1L i0 ;* 0 OONNOVM- O*j
"*' 02;- Z Z 042 0 *WCWCb.0. 1.-o6wx- tx ZIi 2 -41311--co0 m33 3 4301 U 041 1-. 0U 6i1 O*.Izz.0g2 20200n 0 IUC N v -. 314 *N2 X4Z 4. 't -. *0L
14NN- Z0.0.% .2, *114 Xl *41W 3 N O .42 0 Jl4* 411442 41-a422 m. .. . 42On4 . 2D Li 41 04
.0% 4 41.44. a411 44 Z4-441Z1.24 * m4 '0.
4 N4 ...- 0 Z -0.. 000 0Nw-. x4 ZNM004w I-W . 4N
4200212U2242-211222222~~w .0210-22222 4004i- - - - - - - - - - - - - - - - 2 o 2oUi00 2
----- -- 6- ---- ---- N'.% Nmm'w"..Li
T..Zz m c. mm4 44M4-44 0 4 -OOOOOOIOCIU1* 11Z 1N1 IN1 111 11n LiI 414 x~x~r.Czz 4 . 11 11 10
*0 40 0- 0 0~.. 000 0a 0 & o .o02 ,00000000 H.2 0 1000.000002oo
CL 1 4P 11 V. - 0 1 0 - N M 1 4 N w o 0 0u N 11 411
77
AE DC-TR-78-63
It.00 to
ic4- I I
% % ;N b
-- ~ t - 0I .2.
MOI W . % It
.4 .IL . . .N..O.. fu
u4-000 *% j%6OPJS0b 4 , 0 O _ T4
* " m : of a ;A*..4 N 29.ia_ _.N 1 ~ ~ 0 % - % -o . ). a
0. 'D in F :4Nzt* *JAON P. 00 ;=I-;-
0 A :
*n 3W 4 0; a I 0 0:
z N a x.ouP!O aoDoi. hi N -2.- a n - - - 0--
00. 00409400 ga 0 3U - - ---- AJ*
@ P f" w . 2 0 0 S 3- * ' J
4N Oigg 5% % t~ hO£ J ~ 078
A ED C-T R-78-63
NPI4IA4I.*[emailprotected] 4~.*00-.NM.0II-Ia 0.NIII0.Q OEA0- 0 - N M*IA.0~~~P.~~~-~~Np-6SU2. za6O e .00. k aa e e o ----- - ------
wo cccoccooooco000000
4 Z
- 1 , It 1 1 P1 10 I I
* 000 a0 .0. 0 00 2 0o ~ ~ ~ ~ I IL I i L C.0w W o I
t %-t t t- W -c lI0 a . a 0 4-
a3 -0%
'0~~1 A. 10, 1 . A . I 100 I I I
i z N NW w ig w NN 't N Na oS N6.'.I.-P I )0~~l * 0D* 4 . --
P4- 0 P Op0 0 0 ~ ~ I 4 tb .II
P~.CI-hi i Os-X 0N CC -o
A.A 0. 0. A.N A. 1.A. U-O -% * .
in AI@ O P. eel t!C0 *0 -21!2IIA
A. . t~t A*A A. A 0 P IAX OSAI- C t-. C I P79I
AE DC-TR-78-63
----------------- - --- -------- -- ------------
$- a 3 1 .
0MN -5
0I kM 0z - C~m 0 0-U N
1~. CO S.%0
# O's Z 00, US.. U . -1t1 . o o.30 0-Vico -d
PO-:l;0;0- IXX 0 I
h-h CO fu ow4 V y f
80
A EDC-T R-78-63
- AJN NNN NNN IJNN NAJ NNN NNNIJNN
00000000000000000000000000000000000000000000000000000000
a -,
0r 0m3
u 0 0
U. C 0In 0 v. 4yq
4* 4 0 0 U zz 9
-0- t Iole 0. 8ý 00 - 8rljzX
x U.0 Nd NU M0
0. 09-L N o U* I.%
P.- 1- !%, *0W *o 04- -- U k. 0 % )I -T4- 4-3 4 .. 0, .-. 0-W O 0 0 NSS 2
0- I- K -O~ 1- O -ig 0L 0 x
w u1 - 0 -0 0- 49 - A~u A 1 0 6 a 9z"0!zxx- "
LAN onNI h- - -- 8 - 0- knIh
h. h U - I- ~ -2h. LA P- 0 981
AEDC-TR-78-63
CmNNA o N tuut4NNNNNNN NNNNJfy
0 mewww Moco me. ea wee wwww ar mUa.mam mw3
00
61-
ILA
P1Z3. N,
lyu IO
Iz 0 1-
4 4 1.. 02 1101 D* 1 00 0a 314
ox -wI I4 am
.*UU -0 L O =k. cc 0. 1 0
1-2 09 . 0CiI 4O -- 0 0 IX- or - - - *0 - Z -g jI 0 AlW N04 z z 0 0 . - UZ 4 301IA I .A
It 10 a 3 4 3IAIRA 0 92 0 at N6 i.1- 0-- NP11 h 0 h. - b -* .0 M444 0IJ 1 A0 3 3 li- 4 4 101 ~ 20 K *=2 0-W as * *0 * 3 .1- I 0I-P-I IAA 1 0-C 4 0 -!~ .a 0 0 4
IIIx 4% 0X * 00-0 - 0- 3 0
X, ~- Is- 00 ," .i 0 c0
0~ 4i 0-- hi 4 41 Cm3 zNPIM~c 2Pl- O2NM 03%.i -0 6 0 'A - S hS -N- 214 K*~~~~~~~~~' It *-AQ ** ,-4i *0 *-x . imk0*11 2--. - e~WWO V W.e 20 L. .. i.OW0 02,L y- *OIP1 O P4 07- , a4NIz wPP- 1 P tov*~~~~ 0 '0a W3 NW3.,.4 24 31 P-X3 2~3 4.P *S0 20 *I.-
atU0Oa O Is.4 TOXOM IAI.W 3- -- 2 -- a.. a- 2 2~~10-. x0..2.2 0-3
U. WIj WI 0
82
AEDC-TA-78-63
FYMNNMO O - - - - - - N N MM M M MM 4 E '"r 4M
cc N o 00 0 0 0 0 0 00 0 0 0 0 00000 MM 1M3M MM M MM MM M MM MM MMMM MM M
~~~~~w wot W m ~ mff
cl-
V0
4- 30
r.1.Oft itU 0W
W,:01 9 4
.In * 0 iz U -Z.
lN g oz , , j W14.0a-z Il WhO .
0 U0 -wn-u z-W -3 - ON 0'0 0 w K
h.Us NWO w W3 be --W0--n 0 0 0 * 0
w . X 0 0. 0 x 91 x P-0 r4
WU~~ U(n04pK~ Inh ..
-I- 1all
AEDC-TR-78-63
UA J
SO X
w hi
* Tx- Z . ý-Z
Lii n
01 a-w.or1 LW c 0 0 e
-, co J 0I K
I s ON 5 . (0 C
0")w0Jca -- Q0cam -zz or0 W - - Liar
KU0 C; : Z: 0 WU. '0~ 0 a~ w-. UIm 0zK
aL 05- n0 zo xZ- CNa- l OU WN Z 0 Li N )
%A i S e Li .5 - - -- Ko,e Li~i U -L -3 1) U) U CLi UK U "0 0:. M4 -.-. 0 vu . 5'-.U)Z 00 5 . 00 *- ~ UUL i P
00t 4p Li In U 5f 1N S UU U
Li. t U ~ .- U - -I- aD' -0 - S SI KU84
AEDC-TR-78-63
*4444*******444 4* .. .. 4 . 4.4 4... **444 4 4 44 4 4
gg g m ~m z mmD :)m ~n~ I000 0000000000000000000000000000000000000000000000000000
0~ 0 x0
- 00 ZQ !22a
I I- .- I-1-
X I~ T - WW ~ 3-X I W c21 4 dM 4lw -00-2-w 2 A - -4 - - -o u~ '
.-.JUZ - .J LM ~~ OK 0 4Z
W UM -0 0 W0 44 0 ý CCL * £00 :*x 00 '03 z K00 mw
m * .- I c W- - - P-Z O-xw
fy 40 0 4 J.1 w00-b - SIT" 33 V! .J. . -£
ow NJ -M ?5 Ow UN T-.0 hI-I-W
~~l o 3- .J 40 3- U.0 1 ".-USz 03. M.UK 0-12 z..J T 4U K 00
U'1 & MOI O £0 & -4 X 3 0 02r Uý -4 0a UQ
vA 0- z 0 ; ThO ( 0-.31 1 3 SI- .. 4m).. 1-M10 1U 1400 0 -1.4 N- M u 4w -=4
T.S- ONO*-- 4mm 03- 0££L.0LlI -j-Z _ x jU N - I- IJN Xm T 4f11 U, C!JWK:;
US~ ~ -4 .0 .. 8*I1 * . W 4'S *KK
5-~ -- U.JU. *IU *0_ -l -.N I aM1 0 3 MM LC .mk3.
-0) - 4 M . e 5S 4* C - N 0
33~~~C V(3 wwuu NI0 5 0 *IA030 *.LL*-. % *4 6-- - ~I -* -03* .1 - .. 1* J3. W~I '0 h 4q ~ 4I
- ~ f 0SI IAIA. US &- 0K-: M- 0- OI K 'I0N0 IV 1010-M 04 1 I - O M 3 U I .1
- II- S~A IA 0I 01 II O A0 11 ** -.. .E .- U~8544
AEDC-TR-78-63
'asp
ix a 20 6 0 06.
'v @A a ,, $- 1 z 4 . zW
cc u a
-jW La a0 0- 0 -'! '. W x
a, z It 'vo
W Z-ý-j 00 - -N U. T r U. w 0W0zWu0_jx 41 lux
0 0 OZ IWO 42 0- al .0ty v T U. 0
-9.jW 0 ; 2 ; 4 .. 0.z W W z k. -i a -M a
W'.. u un W 'd *a"HAZO Wo rw
mmmýw #W -j LL ... 'a czo-ONO = . N.ým U. 0 -j
IZJAWJTOZO;zx UWO ;am X,.Wz Z% -2;: X-z a OzW 2 Cal 4 .. -9 cr . W .. Olt ltt Lt Ucc W to- tj D in W (Y.J a Ito 0 In 0 a
bjý - &L 9- Czw to ýO K -L.) CP W n a24 0 W 19 Z OXO U, a 't,
Wýý W i! on
IL I -i Ic.9 uzkvt a z tA " 2.910 - - i P.-ka Uj 4A -Cýj tAzw t' C; OD M XV -09V u
X Ujz:.j0(LW -Wz IL U WW- x 0"A -tuzeu
W" x W-N; LL Z 17 0 , !q:)- 04 0 01 1 Z=uwoo 0 .94 . = z04 4 2 1XIt It It . 10 0 cc 000 ý a IL xWX .. ,x HlizomwiS ý%A ýr N - 0 N 11 - U. 1ý . W ru CC 0Q 6" 4 a a, x x 0 a- 0 - X W 1 2 3. 0 to fym 0, c Z w::W,4 u z
xf ýcyx;ý4;1 z W:;,z 2 ý- a x K zI" Orflf":f*n:4 XMXA xmlýx- X-r -wxi VT C Z X x x ý x x0 &L tn 0 - LL P) cu-joýi Z:.O.wW 70: .- JT WN-4 ; CL x uj 11 WX tt 0 zw% -W Z:10:x"uuu
x p 1 xlý02014x- rNT , a - 0=W=0O-X%- zz
X-,Y.C-x nr- WO-ONZ.ý .0 1% 01 QD 0 X Z-OZZWO -- -- ;--ýXo ý% ý
1.6" emno ZXNNCP c27 TZ r z x r x X 2 MO zwwm-mvw--ý I ftI z a raw N U.
MIMZ -a ýx Ll Cc x a L, a a M OD x P-- ý x x x z ujMO mc lax mo 0 0 . . 0 0 Xoý 0 -"MO -- 0 0 0 O.K xxxxa--N:X-z
J, U.WLw atli,04 U.Wo Wx'ý-Llmm W ý-uCCxxxm4wwuuumONM -N - - - -
CPCID 0 0 W cc a- u
86
AEDC-TR-78-63
40 - Of4-*&-ta~ 0::O-E4ml* 00 In4P4I4- OF :0.4P e 0@ NM
a 0 *4xf p.m W
ol- cc
1 -0 cowoI
eyUx x
U0 J ol 020 --M . *0 2. -; -W -0
Z ZI 4M 4 4 WINO U09 -:c 40 0I; a-v Lai .
US .0 : CIAu
s: . I J..l 0 0*160 U 00 mN
-1 00 Z3*z al9-s ZZ 0 a, z
me 22- 04*010 A cl Wo*ma
K WW 87
AEDC-TR-78-63
-0
0 * -Z 1-.3 0;
bc W
.U 0 z 2c
0, Geo. 0 Z2Z.9 - 4 q24; % 11!u* % m W 27I l
0 mt 00 mC a '
-0 M- 0 22 0
-0 0v* .J x w N zN
ml 1ý Ui 0 N 9 - P4ýf ý 0 x x wIlm y --
LI~f v N 11
* 4 4 * WI88
AEDC-TR.-78.83
it. 00
u t.
US 2 O2z - 0
z N zo
I 0C L 00.w14 0 -o".i ux
,c 5o w o - i _j
N-n 0 o 0- c o ox. x w w o w ý& - N - -;-- - 2 MZ!7
.5~ wca, w wgn *.~ 71, -jw* ,z Zh.- a 6. * -a 00
%0 4 hi - 02c3 a, 0L I% -C VM
4 z 'I'V *05-5S wI(SU1 a- 5- D Z; os -- rW o 01
It A .0 -.i5 .&440" .US~ 0 ai w4 * hihia n _W 0
.3 0 v 2 W w
(SJM ~ 'd- -**~Z 289
AEDC-TR-78-63
w~ooeeMM *O Ooee9oN *o P.:::=4eno P. o mooooo Poaoo;Nm 4eC, ooC
Owlw
-~. a t W()
fnK W 42W7, - r -i X, 4 215--
-0 0 IC 2 0 W
-ea 2 x.-a- 2 * - -.-- ZO
2C 0 a --* - 4-m- Zu "
-Af4 N- x .4 N- 0W 0 W a 0 z 0 I- Z -
a j 0--- - ! -.x- . -i !- . ..d 0.. AS 4Uo tfI _ift )
* ~ ~~~ Na - - - W. i( C I -0 IOu owl 0- x2 10-. - Kzh
i-3 IO-Z ii. 'I- h 000-'*A~~
fn)W--K.-6 W u1 In J~ u JS 2 u 2 2 4i-. *
466. U 6 *-N J a a'S *0 06 20 4JI i-90LsSW.-- .JJ~
A ED C-T R-78.63
2Z2 22Z 2Z2 222 Ez~oo~o33 3 0000000 000
in 3
0- 0
'M a-
zo j . -U M
:400 MIN 0-0 -4w
n n , . -: T
_j 0 On 0j I~~~ 4 0'. %
.4. r 11 x
-u- 1
*Z -4- i91
AEDC-TR-78-63
W;U OW U20-2 2 0
I-IN 0
z US 0 *
@2 Z4 IT, x VI 4 %
02 N 02 .2 wnPI-I U -M- z
14~~~2 * - y0- - 60 *40
.4 Zzo 2! !0U2 *0 fU j U qit 014;;. 0* . -. IUn ; 4 *.
o oh 0 0 -00 0- t ;.0 . 0 ZZW *.U5 UI- I- %
0 04 f;N N 0 1M. . 1
U.2;- ; 0-MMO1 ~ 0
4-~~~~~~2 223t42 V24 011
sz-~~0 4w 244;A~ *2 .000- ohijZ2 t ~I 01 01 *..3UlO IRE NO ** 0 ! # s
-- &; a M.0P. UUO-_0 '.mo0 1 2 ONuf IA4 C 1 0.
0 a- m x0 U .-w A 04 -- om S. wN 40 x..1 60 E I1M- '.'00we itMe 4 N.JMo mw mo *. U-0@ I . . 0 i122 ~~~~- - ENZ - -*, .2V . 1 0 . O N h h u 1 2 *
NN *0 hMMNO *-04* 12..NO N4W. V O- ox04 1 S -MN _j s4 3. CW4 11 *N =U v *l 19 z vfI*z~ z !*1h.I~ ee=!tvIQ Ro e NO Z a UK 040 K 0- evo~ o ~
AW A- *1.M.U ZA WO 20 hio 0. Wu1024oO* o
**0 02 - - - O . M 2 ~ 2 W UOU S. .. N S. '. .N
-- I ~h -I SZU 3V VI E - - - '. .' N' %' hi .' UO92hi
A ED C-T R-78-63
MI 00 VZP
aW~ D W W W W W W W W W W . W Wh WW h W WIiWWWIJWhIJhhWW&2209 &a;R~a~ac w w w.
a au& CLa. C
z
-000
4 T W
w -1 0 0Z
0 4, 0.. *
0 00 a w -jl
u 0 Z:44 mt Ow ww.WIL Z 4 m-- O a x
'a ~ ~~~ ~ ~ ~ 11a0a0-t S -&: !!-j-J t~~tm2uz z m 4 :! .
Ow v w -W vu 0. LO 7!7 ! -a IL. hi j m m on j
-. ONClpmii - . --- . .71 -- 20 w v l O o 1 " 0lI W -_ J i hi
"d M- - 44 _we,=OZ A; zoo jzzEzzj-- d0*V -0 0 2-VIwz=-O -1
LIoW- *l U toJUU VI.J4--' 9.I~U. v
* * 4%*. .*U e~.-.I.-~O 22*fJ4I0 2 *~))93I
AEOC-TR-78-63
K C
mum~a~mm~~gg~a~m or
%0
Wz - Ow .. -
I- Z -- -)h z 00 x-
@01 -2 ',~l -~ - I-0Z nND it -W
0Wa .. 0 W"4 IDO'ý 0 0 -. Kt * 11 -- a' z w N9 n-9, .- 0 w w
-- ~~ -Colo!-- hhi% U 0~ ~ ILl 0 -oZ a *J
IL 3~w nw ZIA- - h 0 0£~~ .U hxiw2 N Zhi U0 J -- I - 0 0 4 00
- ~ ~ ~ ~ I 0-9 2 *21
9 -4 -i - n -
AEDC-TR-78-63
010
0 fy
ow 0
W W -W10ILA
'A XZTip 0 -3 w m
A; ama N 0
4K Z VIJ 2 V! 0 aX
to 0 00 % 0; * W* U.Z- - --. -*-Wl 0 0. %fmI
- a- Ws - to 0- 0 0 WISo
UVI 1.120001 004.jOnA
959
AEDC-TR-78-63
a -
IA a
go U, ý:ZN Oa Ia * 0 Z -
fy 31 0 6 A V).00-."
-- 4 -3 in cc 0
IA U N am tAt ; ýTM :;w 0 S S -W ;4
0*0S - U z . W . 0 h .WO v -o 1A 0 1! C .A4 ~ - - i ;l i-. 1"
S z . j , -J, ** 0.. I.. . -. . .0.
1650 0a ey MA In ty f" - 4- Ifu 6 in UIUI AbU -. l- *
5- *% * 40 0 ~ ft * 0 aS O90
A E DC -T R-.78-63
Pý P. P. P. P. P. -P F
uwwawwl!6 wwwwwwwww www~w www w www wwwwwwI.JwwwwJ
IL~~~ 6.,04
sS1 3'
4 4
c 0
z z j
I a I f
" - -)).I
=- 0 0 ; !, !i.
4 - - P. Ia* 0- Lo t ..J .J - Z -
.e 2~w zi 40mo mIz r C! j Ij>
0. - -. -,"
I Z-17 Nio Iq :0 00 w o0 ! -j o j- j
oU'UNm N N Zto, A si*us 42 00 -0 -0 -u123
* ~ ~ 1 01 ~W .6 *S N- -U 4
- 97
AEDC-TR-78-63
m m00 - - - - .1.mNm1Jfn MMmmmmmmmfImmm 4EmEEE 4EU1m
I-
x I
uUx
9 4n
c 2 l
w at
o 4. -3 -- a.
IL 4 x 1
o- 11 2 z 411 co 0o LI I I. _jc2
0 * It, -1 0
to 4 ..lw w * 2 -0
-m cc m - m 1 - *2.,u.:!I t0-. 4-~~ ow1 w Z .- 4-
I.* w a - 0x !& 2j t*'- o- -'m .j~ - -L 2Nf
3 -411 2 * N4~4 N _j!-jcaaL
0 0 '.2 0 -j~ O -Oj j -o w p 20 0 -- o0 4 0 D. .0 oc onLj--wI .0 o0 2 24 2 3. e~z w ..
4 11 4- w4- w 1 0 001 IA 4A 04 2
2 t * w 0 * ý ccUU o0l 00 24 00, in 214 00n. Lm.I4--- 4 .
4..U 4 4- - S 00I-4.J I-. 0.J 0028
A ED C-.TR-78 43
z Ue
#A 4
W hi
LI !U U*~~ -i
3:0 0 .
%f" In .--
0 Wz *s -.1-
4 - U. -n z ~ giosi W3 a zo ,. Co Zj a3 a 0 S0
U ~ P a Ua a4C K J..-IL 31. -* ?W ; 09-5- u~m ~ -;m OW -.. U
9 M- 4 ZA0 . . .* 5 ~ . . U.sP .*
fn)4 P.- -. I DhID~- * 0 a
~ U : ~ - * .0--0990-
AEC-TR-78-63
op.o.-raftam 06 a60 6 0 L&&& 0 L0 6 LC L ( 6 9 LC L C LI
.4z
4Z
C VI -
N cc V) z
- -a4 0 N z
so. I 0 0 u
4 0 4 0m- LI - x
W- 0% !. 1 ,
z l a - 0- x 0 -1 4a wU 0 o0 OZ
a - 'S -- VI ! C oI .z-ý x W~U %L-mz - - x. or a1#4- wx m* 0, 04 - ! w
2 0! 02'WO Z 0 12
0 *P.pýL 4 0-4WVU zxwv -.q - m .Jwoo 0 00-~ I
4 -30m-04 10 W! LIO03 - 8I 1 a 0 VI.
0 ~~~~ ~- a- I. '0 X0P-34' tin i WE U. 10' ME Ic 414 ." ý0 z0&0U . . C-ID
* N I 5P-I . 4 u IP IW 1, a
~ 52 5 4- 'S - 0 0 100
AEDC-TR-7843
I~P9Ehi949-6eO..NP94.J94N6tnNm Ill z N E 0 0. N 4 4 0 M-j 0
9--9Ip-51e6@@64e 21@ 25 0@ oc000-------N2NNw N42
4~~~ $. 9 - i .h
I- 4-.: , z- .j 9 9 9- 0 ,2 1. 2
4 M4 UN 0.L * U 4-
40 2* *0 & * h *
44 21. I% - -6%~ 1% -hi 222
44 m1 04 h:, xJJ 2.
00 0 ,O -' O Z0 h-CL 3 - 2.0-11 00 hi- OI -%c. vu9--
2 .J20E O 1 ft a _Z Z U4Z
Zhi ~ ~ I - 03 . .- 0. 2U.-. 50 OW ; ;: '%- 44 U 4
412 IA *~0 1?42 *hi *% i01.%.I l0 4 - 4h -MU -0-0.0 r PC P1.. *0 it-%a2 ZU -. 0 %I-O
.4 ~ 3~I 0 h4 .~ %406 2 Z 2.0 .9 -h.09W P.-. ~* - S. . .. .. 2ý . a. O4j2 A A
.J~04 34- -i- -952j -4
hiW43-- .1&P m4-4 ~ * 0 4 42 0m Qh4IO U0.44Jl
%9U 1 Tit:!~*~.-. -4%.2 4UIUJI. -0 Owmom124.Ow 10 hi x 09929N O2 * 0 03 W4
0 44. * 24h. 09/I 44 - U..4-1420Azo*. 4Z- *.... 0.; -' -' 0
4 -ON;- *K 02 NNW a1 h.2 *' -h.* .N 4 l:4 -
0- ~ - J Z 4 hj. 1 U.3 . ...49- 2 .434N v1
WCK l.w ii: zW X, * ,*:44044 0 x
44 4444944444444-4422
a x-~,; -K~*1* 'L a 5..0990,5 -9N 1 0 00-zc
iw m09 4 -!4o0 0mov0 . U: Z- M * O- -- -- 2ee *9**9 IA OQW ' 6ID 2 4 2 % 2 004 22 2 22 02 2 2 Cl o2O
* -h* * h0 @ /9*. * -m--------------------------------
AEOC-TR-78-63
06 0
2 0 0
I. A 0 0
.%. IfV "
C S wx I- 404 -La4 c w 0 1 0 -
0 3IN LM --3a. ~
~~~~ iM: a-a z 4 P.--040 - I - I a l
Fn ol -L x It: *.0~J 0 0 - 4 '0 -mwm o6 0 0U A z 'A - 6. V 0 U 4DXo
go Vac.. c o a. mo -a-. - .0 - o
! - -oM0 a 1,W2 VO a..0Z -Z % U x
z M 0ýI Zw 4 * U 2- -~ 0a j c 0 -j w-a c o a a. a--z
a..ao *aji m - -O S. S r. a3 wq0. V01 Mn-' 2 OWON H ) 3 N N- .2 * M.i0 314n0 -CL-am g l.w o ' -- a 7, owa a. -- a*2 m2 a- 4 a- N-w aoz -!,I.4 4 4 -. 10 a- a-10- - -25 0W~n 5
W--) 00- 4 *u* . 0 05 C .
a-Cup 25 0 I. .. UUSS. 5a- 40. U102 *
AEDC-TR-78-63
jI ah .J j j Aa j -J -j- ij P- 1a... aeaea. ee -a.. I a.. &a c. a. a a m.~
Co IA
it x t
lo OW -,C ON
31*1
m e cb O w14a
N- -- 0 -x 0
-W 2y .0 VI in0
0@ - . S it -103
AEDC-TR-78-63
x J
4 -0 ( o
2~ S. .. 0hI 0 C, -02z
hi~P (Y- .4 0 0N Z SD 0 ID 00
-62 1: 9. .44
- 4 4 4 4 00 0U 0 0 * 0 *0. a4 *
0 -p m 0.* 0 84 4O U.8 -1t 0 .- 0 04 - 70* 0 P- -1 0-
1. 9c 0 - 0 ** I F,04 00*9w. kv -0. 0 N-.. 0 .: . 4 .. 9
U. '4 ;z T. 0, -- 0 0% 00 ; ;nN '*404. zoNZN@00L4ý z~ : 9ko0z*99N4 N:NI0 0 *Mx N@P. 4 NQ n
,, 401 IZ '4 N Lb- 0 2- 0bb 0N -.
U 0cc 0. 0 x 0 b- 548or% b 0 M S~-r 9 48 05 0b -1. u . 0 0 4
iu. c 10 ZN Nz 2 4 0 .4% 04 m x0 4. a .w .fm46 CP bb-9Ia0 0 05 U 00 . *W 0 * W 0 01 N 00N0 0bi4ol
O 4J4 45 b-00 aE~ 0. % tI 148 *I bNN . Nm04 n N
z 0o%4 Eo*imS b- *0 5 1 0 0WW -0000410010 000;
22OZO3 *:1 or5 *N jo US *.. 4% ZO) 1.jl- m * a 405 400a Ob~ v ZX 0 -I 05.aI IS5%5 Z* mIOE O * " *0 N N N NI# .D0 00*W0 Cob-C. z3 Do 000b-ObbSbb. 46 00 0EKw1. a -UI E5.. a -W S.4 -N 15-I Sb- S -100 NN S -*905-.8.4 *OW E %..- 248.2 * -S o *-.) =D9-4E4 M DD99OZ~b,
- 3 0- SZ~b~bNS.bs% .E 0E0 EZ * 0104*
AEDC-TR-78-63
~ ~00- *"0000@N~ NU 0N 4U0400 0A~ * A N 0 N MCI 4 0
~~~E'1~~~~~~Ec *E* *E 4 4 4 4 4 4 4 P N NC OS g 0 4 0
4444 444 444 4444 444 444 4444 444 444 4444 444 444 444
mm~mm~ g~ gg~~ ~
I U'
0%0
4o 0 00 oa I - -C
z*.2 ~ ;T 10- N000N 4 4- 0 0CU
IV 004 *In 0 N9 pi40y cn! *y V@ - 3 . P -oI D "I
002 0 ;3-1 % 0 # 'A CL
NNO% *;o 0- 4 It-gO 3
ty 1 4- * a0 0 In 3 - Z- I34 L
% - * 4 0 0 0 0 00 0% ~ .
3- a .533 ! a- x4.0 4-0a w-- 14 a 01 T 00 -(Yw :
.0~ *-.. .4.~S IL g4 eq- 0x~4'-Ie~4 * 3 .- g m4 e -- o
0@*~Clw Z-00000u O ~h%1 N 4
*3-0Z0Z- Z 0 w .8 4 2i - W -j a0
3-2,00 00 N ea i 3W N. 2. CL4 0. 0-3 ,3 S 0wNZZZZ z z4x& o 0 V 34000 4 *. ,S.. 0--0 20 4 4- 4 0-.k N 4
* 00 S00141414000S1g3-00% * 4W4MS5 4
*00 5514-9-44 *N %S I@3-0 )- 3.40 N10511
AEDC-TR-78-63
- --.
30.3 4 C 44::: ma maw- a *
C~I It. o 3*.0 4I .1 0L X0 .3
In A Z o 05 11 A1 Na-0.0~ ~~ ~ I- W - -- S Z.*. .*~ 4L- a. NIJ .1 - *...*- ~ ~ I .4 %AJ ...- 40000f
- 2 Lf LA KD o- 1 *comNx.Z0f - SW. Q-0 . 0v (a. 3fd :0. u- N*. 0 3-0 0
I 1 IM If M IC U N I S- .- imam Z
I K S.1 InZ. z x4- m Z;- ~ 5 SZ **4a C T hZU Z .,04 . -4 1 th ,
o I10-0.. I ~ NM 4 LP MNW-L NN00-1~005f. 0
'a a c . _ .:. 50.0..3.3 54..Z- 0
106
AEDC-TR-78483
0 f0u of *c -
%~~I 1, *%4 .o~ C N U. -
dc a %I.
'. 0h.* aK
m. I. IW
a ~ 20
0 & M0C -I9 a. 11% 2w Ex 1cc- I a x OK
a &M .4 I 0 IW - 1 z1 z. x Oz-P I " 1-ý
M- CL a6 .*a : a
(9 0 31 3K a. _ . -u z 4
* ~ ~ ~ I i % ' . .U-. 0 I. *1 0.0' 0. '*o O =-W' ..
m 19 300 0. K0 a . 6-% 3K% K
z a, &3 0.a 040 C 4 ft CC0 -L s " Wý 0.
19 10 -44K 0& 0. -2 :- t1 .*h.
y I-.00 -J 00A 0- 6N31 at Z &30 *3 0.14 X-3 * *1 -mISO
*~t NM .V 41-1-0 1(- 2 K 4 * 69 .30c; U
N 190 2T o m 4.-. &-0 IL m-4 .- 0 T6m U10 *6N y 41- IL3 K '0 Nap z K *K 4 -w ;4O.%Xi 41I
- 3 00 34 0. T 30.0 C2 T t0 &04. . 40 0' I
I. 1- I.-00..K*K 14* I. w 1- 98.101~~m ".11 PC *IJ.1 3.2 . .. 9 Z 0 1- .. NN1-.00K *3U-Z~0K K Q- (1-O ... 0ey
* 0'.1100 -00 -0.0010. 6311 1..9. 0.1.1. 3 N M
3- -.. 0.-03- -- -32- W-~~6113* 1-l107-
AEDC-TR-7843
3'. '. '. .. . 3 2'.3'.V.. 3. . V- I. I. V.3..3. 1. - b:xx MK Ku MW
ey 0
10I inwS
4 .J L a I
Vý IS
0% w uV"M 0 N C
CAM .0tm w
00 046 a m
Slaw
o - x
toS
.0 4
_00 N W. S
4 *U40.
*i:f MSE-
IA SX
* 9 -- , AUz M S
-u .Y m
C * 0
1.. 4108
AEDC-TR-78-63
o00000 0 008
00 R00000 0000
00000E00400-
K ON ton 40 7A iMO 064 41,N.6
400000000000
00000 0000 oo** - - - -06 01-* MOO 4p00
0 004
P. OM 'a N- 0 C 44 4
0 ~~~ 004M O 4 0641 O MQ ~ ~ ~ ~ ~ ~ ~ ~ -~ vC '04 N P. (I N O.P.0t -DMk*10t00000. N v 70 &1 4ý 0 01
0,0 00 0, 0. :(03 ' 000
40*00000000Oc0ee000-0
mp.O PN O- 6 Om NNP.P
IL P 0P NO P. ol P. hpImp6 In4a 10 V. 00
.K 0EP0 A CP
w L -9 ON406P.PCOO N60 M . 4a, 0 NO N E C'M N C PC O.Pw in446 6. 00. . 044NP1 P
0 ~ ~ ~ ~ ~ ~ ~ ~ * o toP0 6 0 0 4P .~0E0 PPNNNN
4~~~~I 000000 000 000 0*P .. 4 .O P. C
U~ 4 4.
lee 04N0P 4 N O . E'r 4P1 4 OOD N.0 , CwwmwN0O04 4
* 0 N D ' a g cc.0 40 C N0 C 0, 0 0 4D I- N 4 N Np -M ' f0 0 0,0,P0 &0 6, 0: .00 4 6 -. '04 NM A xl x M ~ l
tl0 6001 M00004 .00060
1 80000 000000000 00P. 3 6 o~oooo ooo
w 0 I o:OPC O00-0 16N 64 . 0. 0 -9 0 Nc o0004NN6P.c:! M v 000-ý zC-4-N. 4P4~64 0 0 00 3 PN64CIO1l
CA 44NC~'4J 66P44PPC6 I Ii.N .... 064109P1
* AEDC-TR-78-63
oý NO 0 Vl V00 , 0M
o0 oaý 0 c 0 0 0
* InN P. 0 000000000 0000000000 0
0 0. C: 0 0 0 0 0 0
4 4 x T 0 14 0000o0 0000000000 0oaooa
-I T a 0000000000 0000000000 a
&. N. N , *
UP -&' 0 4; la : N 4o00 0 0 0 0 0 0 0 0_j N Vpm.p -Zo 4N N w~& oo
2N 0o ý ov~ Uoo w c: 00000 00000
0 w It. 56 -In 0 .0In.PPp- OS .N N 'm i z 0 N c A'2 % 0.P.40P UN0PMP
6. y 3 hi 0 go:0ýg
o I . B . h N 4IIn~ tl.PP PPN - -
4 IL o 4 K oB.c~t o4 : mc o Mn
c c :4K!= t :=*ftcf 0 00000 00000o. N 3o o3 i4 v mz0 wcPo 00 0 0 00 0 00 0 0
o o m N N
IL .0 . .. . .n .0 .i .nP O. . . . . . . N3 j - -0 - ---- - - - - - - -pPI~n P-t -~C1.C4 -
ar a M 0 N u nN414P OPIPONP
- 0 *4 0 0 0 0In4p.0030- -N N co o aNc
p. 4 1 PP moo Noo & m 4UP~~~ ~ c n 0 - - - - N N N N N N
-Q -n N0 U p
A. A. 0. 00... . n0 't .";0p N0 4 In 0. . . P-0
z a. PP m nmN .U Ni f3 - N y N 0
p. 0K a: 0.on ~o2 P
c4 U cN pBU
m wn o. -W M - rw -00 0 00 0 0000 0000 0
N -y c & 0 4 w 000000&000 c-00000a0m 0 0
W 0 m Ii N1UEI41n LP0n .nn0 p
to 0 0, P. p. o0: .4 1
oi o0 r N N N0 Nl N vN NNNNNA vr
.0 2 4i 0N~ 0cl- o ' 0 0u m1 &o IS
u 0 j A 0 - - - - -0 -00 00 0 -0 0 00 0 -0
In p 0 p I . .. 0 p. ..- 0P3n .0' NMNE.Peo o o0 m .0 4 N 0 0ppIn0 POn . 0 p
0 1. moNN N 0 * .0n4 0P P .- 0pI00 1 l
P. .. . . .B. B. .0
T. 4 5 N U 1 51 5 1
U 0 * 110
AEDC-TR-78-63
ONO0 !Oo top-' 4M.424N 10O N
N -N4~4~0N 04?m0-M4U9
- 44 N04 Ic 44z ON 4CN 44..4'
0 0 - --aNN41'E 61
0N "fyy9 ry ury ry'y y ~N.O- - - - 4
o0020000000 00000
0 oo-o0 00 000 M-4%44
4 N~I~N0PN 'COM4 -N.
o MC ~ OMOONCNGI- -
0000 000 0 0 00000000
O~~~ *N1.. M P .. e- N 4 0444 It
C. 0 flN04Cr 0~00! 4 P C: N.r
00-04~~ Om Nov E-0
C 00 00 00 0 000000004 1. 0 01 N
* N NSM~fO4 4N N-.4 4-r04IL 4 4 N 4- o ~ .. ~
AEDC-TR-78-63
190.lfO 0 a .0019moo 4004NON.O m1 ON.10-NNN m4000.U~4 44awo"M0190 NOM 0019 C'N. -. go,41l~4 m
z A,2 m 0 Oon 0000000e 0 e0 000 00 00 000 0 0 0 g 0
4 ~ ~ ~ ~ ~ ~~~~J *Io00U Ily00,0I o0II00 2 OOll~f~f
~ 0.-19 01IJ M4 .19,0 0N' .m 9 N ~ O 1U0N
00=O- - - -40 NNN19191WOOO Ow *440fl*010 0 *ro. N..-0 a
--------- -0-----
0000000000O 0000000000 00000O000000 0000000000 0
- oaooo~o oo~0o000 000 00:000 000:00:000 a
o 000000000o acO00oooao o OaO o 009000
00 0O IOP-4-0 ON.1-N019O ONI-4~I0 fJfO,01.09 939 19019 m.19P.0 m : N O9 N G909904ý 0 N04 a9. *9 0 904 4 m oc . I 04 7nW.J04 9- 04 199. c 0U1019Nlove
O0000000000 00e0oooeoe 0000000000 0000000000 a
0 ý00 ,G Z.P9.9' 0 0, 129--0P0 *0,v000a2&a, 090 0 9.0,0, l92~4 0 0000
LI OON INO -0109 4N0 *0 1 - . 1 044%0 0919 - 0
0000000000 0000000000 0000000000 0000000000 0
000000000 0000000000 00 0 0 0 0 0 0 0 00m* 00 0-00 000 000000M-0 AM0000000 Won0:0000 m
0 44 002OU 02 N119 190 0 0 400 49-9NO.0M=f201 4'29. 0 1 1 0 040. 90 0 m Zo 0 or N1 20.9 01
I~~~~~~~~N 00 - .9 9 02 2 2 0 0 0 0 0 0 - - - nN N
m-ninmm----Mm- mmm mmN N9UN N rNN N NN AlN N yN N9 N
0000000000 00100000000 0II000000000 00:0000000000.ooo .o .oa o 0...o oo .0 - 00 -000 000 0
00 1000000c:.. 009040Z0000 ccU09290 000N010000-.N42 m 0190 N O 44N 0 N m 0 m09.9a090 a 10N101 -o 4Conm a m
000 000 00 00 0 0 0 0 0 0 0 0 0 0 0 000 0 Cmu 00 00 000 0 00000000O 0 0000W0000 0 00 0 0 004~~ N0O0co-O O ý0OO0M0O 0
199.4000-040 04 N -0N N04N01 N..904
4 ~ ~ ~ ~ O: O O;M 041N19004 019 .. 9 409-0.- N 4 . 0 4
404.9.22 0.N1:M0' M-ONOI-00114 0.0 9 0 m0aM-- ------ NNNOn NNNNN NNN991111 . 99999 IV 4In
ru .0 00a
040 IV0
N 0, :0l o
0 ~ ~ NN M NN N NNNNNMA41 1999999999 0
I 5a.
x v
112
A ED C-T R-78-63
NWoo-OMNP-o .0M!4JNNNU, 4!OwO*U;: 4 0 -. N.U
NAM aI N. P. M !- M D 4P0 0P mpNU. P.N 00.t@NCE: OM04 I Iý ! . V .. F4 . a 0 F. M : M N a - 4 P0 . M 4 ' # m
00000000-- ---- N PPM M~MM444.44.
EeNP.-u 0P4 Pi.CN00 N...P. 00~N.0 P.044
WO4 ,V a ,Nf NM4.0~ q N.0- O *~ 0 4 P. 0-M AA AA .0 *ON
0 ~ ~ ~ ' 00 0 0 0 0 0 0 0 0 0 0 0000 000 0 0 0
c -0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-N00 0 0 0 0 0 0 00 00-- - - -00 0 0
*N4..00P-0 *000 0P-Pb0 * 0N0. 0-..4 0 0 *fi.- u
00 00.. 000 ;g4Uý *M~o.UN'CON ;gv;;V ~ 0 - -0oOo a ogC PU',4aI-P- EP 0ofoP- NP... oUPo!!o 00 N 00
SN NJoumom.. * NOC ANP-PP-0 2
0000000000 0000000000 0000000000 0000000000mv-C 4 0000000000 000 00 00 00 0 0 0 0 000000000
0 O Mý ý* -0Mu0N0N*O ! 700;m" 4 'a a N 0,
T NOWOOM40M 0P'I 0 N P *0u~ mu
Q. 0-U 0 C0 P 4N ~ -- ~C -44.4N M NP0I
P-4.1 j00 0 fý00 N NIc, 44U. oCP P-SO, .0 . J N
0 0000000000 000 00 00 0000 0000 OOOOOO:o
Zc M U 14 ol to M M t-0 - LAOP44a .4 M N~4~M N~ 000U N~C0 O 4 N C4 O N O - 0P.2P-jCZ.E -flN m irlo .fm~ 0'~ jC90 .~1 0-4 0'0 2M ol j'.j 0 N M O~%w 0, NM -~4 N .~0 P* .
M M 44444, 4*0%.~M~E WOO ~ ~ I , 0 z 1 a0 O uM
00 0 0 00 0 0 0 0 0 0 0 0 000 00000 0 0* 000000 700 0000000000 0 0 0 0 0 0 0 0 0 00 M00ACOA04 Ltoo P0.40 i0U404.. ON 40-A
2 -t0-:& c tno.'v 4 nz: 2201N *N4PN. N0-CP.400 a.C0N 0-P0Coop- ~ ~ f~0 PO N-.N2r 40N00U 44rO Mc 40M aoNP.04POgoo N0~owM4~N: 0a.40 0C~ '~U-P 94- o:0NC- l0cO4MPNM*M.. ý QIMt -1 ý c
- WOO 00000 000000 00-0000000- M0 000000 0:0-
4 4 ez a~NO0 0 .P."N4 ZNNNMN....O C00 0N M 0.0, wxt4 0 4 0..I I CI0-... 00 NP NN NA0 ;mOOm 4m O...N4 *0In0I 04. 040 0CMfIiI OC
K In 4NP0NtN -'00N -. C..P- 4W~.NPNP.N.N .N0~040n~.. I
U 0 .P. 040 0 Ib0 O" ~ 00--. N~M *lCO-113
AEDC-TR-78-633 M l
-.~~~~ .0 ..N N N .NN.M.P.f.M .l.t . .
at- * N SO
0 00. 0 0 .0 0 00 0 .. 0 OO O .**..wf...S..0
Of .. N0 Of NC 0 0
NON AV ~OO4 *WC~4t - 00.* .M O
O00 00 0 0000000 00000 00 00
2 CP61 Nt4t0 4P 0 0 0 03~ it n MP4
NA. 2S~t MO I tOM g T ;v *P-lut S0004e .ýMNN
4,1 .4N4444 M,*4 *tdC . . .NtOW
o oso @@oo~ 0000 000*0 0 0 0 00 e 060
it;o, t mnr'tA,.
JS-.'I4- 0 fJ-0C.4C't.,
K IMW.C 4cc~ 49 U0 0000U4a-EN 9*t04004 *4
00000 a00 400000P 0 0000020 000000 00~~#9.0000000 000000000 0 000 0 0 0 0 0 0
*~M F.~onO cc 00 co 0 00 000
.)- -- - -- --- - - - -
-- -- - -
*1-- ....N..N ..N.... .4
4 a a 1 Pt- P-NCOOO,'i A l~a0N ~ .
40 o'NC n m Vo0a
-Nm-m 0 f...4.4N01 1ME-I -- I~--O~l - NCI t N
cc .-E L NC 041)0 OON. 04.I-C 5 E .'4 N
I- Nm444ý 0 ,g; m ' 44
iji
U em
114
A ED C-TR -78-63
~~0N~44N~P.ID 40.P 4 C4N4lN - 0 N00, 4.
--00~.~N0 4 04
- N 4NeJ44 600000000 4
£000000000 0-.0'4'4r-20 I
4~.I NN004M 0N..400..z Il44 4~ N~4UNOMM;!
M44NC.M*-O~c 4 '- 0N4004mc' 40J&N m 0f
c 444c --M 4a0 0 c .~i 4IiU
* 4 4 4 4 4 4 4N4 NN . N F. N .No*4** 4*4* 4 44 4 4 4 4D4.fiNN
000 0000 0000 00
0, 0 0 0 0 0 0 0 0 0 04 0 0 0 0 000 0n
m N4.4Ii N 0 NN0-4MO4N2 44 Ow 4Nr.SM 404 P....44 . .~ .......... 4
C f~ N4.N C4 244404#44 ;
v9;;;z ;999;gg;0Ng g0~0..N
0000000000 0000000000 a
- ~ ~ o * *MiN i440 00SN O~l
440 '-0)0 *I--.0S04A S
IlN4Nw9!:Mcolmo to4N4N43=02ONZ 04
:- w ,O N 00_,, .O NSN N
.. N ~ .-. .. ...... 0 ....
--------- N ----------
* 0 0000 00 0 0 0 0 0 0
2 0N9910400 4 4P.o NO c U)4 9N0.9.0z 4 N m4N4#40 " 4
v a-, 0 . . .0.U . ..- 94-.9N4.. 4i I J
? 44 N9N4 '.0'0499 U 4oct
U *40CO NOC 9- 0- J4ý - oo
N00 00000 ol0 00 0 0 0oe o oeX 000 00 0 0
w ~ o ~ oo 0 0 0 0 0v z)) N004. N)U..4)U
S 000-.ON- N4994U)O1145
AEOC-TR-78-63
InN ~ 92 c No.NN N ~ M E ~
0 0 0 0 !2~ ~ !ýC~ N00 n0! 2N00 o o c ec o o c O o o o o 0 0 0
-~~~M 00 0 0 0 00 0 00 00 0 00M 000 0!! N NMeoo oo 0000000cg NMoONeo! o 000000M N0000: N
eooooooeo o00000 00000 00000
m~~m~ o m : luig12 cuu ,sm, 0 ., In -&a. 0,a a
.NuN.0O . . t 4-41: tt c .- 0 0. - ue-. I EN 4 .. 0%.C %0 -
a 0 0 0 0 0 0 0a 0 0 ý 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 00-0 0
3 M v v- &NVNNNVV
4j O-4 0-. O 01 40M.-. O- O M N O 0 OI4 QO
La .... ON EM O .A~0N 0M 0 0C %.-N4-. N 4~.1164
AEDC-TR-78-63
Ill ~ *~ lom -4:00 44 0I:gPd 0N,4N:6N
0000000---- NNNNN NM~MIMIMM
* *440P.1S2N4O *40-P p- 04U4 0N~~0IIP-:Q N 0:704o ý.MM& &:QO MA&N
NN -- U*l44444 444444444 co
o 0000000 0000000000 eo ec e
* 0 00 00 00 0 0 000 00! 0 000
4 *aIU 04-2~ orN:M0Ghflt 0 04N@O-.
2 ~ ~ ~ ~ ~ ~ ~ ý NPI-C~~- **NNC.hcN~ ~
M__7 0I C0N4r-C P G 40 0 ,0- N ovN0NC --- :N4N V40 OP.0* M I4 44444
N4000000 000000000000000, 00 0 00 0 000000000
-0000000000 0000000000 000000000
0 ~ ~ M MOO - - P4 'J-I INGPM.
(.0m 0~4N P4.4, N... *4V @ JP-O .*G soc e ley..M I"A N lycAly404 -.0 N0
0 e 0000000000A 000000000 00000000
4 4P*tION 4 ocI0JZOM N 4 N 44004jmm04 0014N m Q L M.@4M: ,0-! vMN00 mJIP.40U ~ ~ t 0*44 C.. 4I4N404M Op~4 PP P-0 0o0.. . .N 4 . ,. . . 0 00.
0000000000 OOOOOC0 0000ý000
2 ~ ~ ~ ~ 07 M 0 4 4 0 P4'0. ON N t .P-NM N c.~A 00 M@040 N4N NO 0 44maJ'
U C0IV~@#010N M4~,fJ4-- _NVr.M
I 44 _yN4 0101ýýXI -44 4 M - - - - -- - - -
11 01 0 '3
10 z-c =m~NNN WKNNN'.INN2~ ~~~~~~ 00 00 00 00 0 0 0 0 0 0 0 0
*I~~mgg i *g,,e g eg m .,.- 0 0 0 0 -0 0 00 0 -00---0000000 --- ---
U Mp-0444N~ 'r 4 M44 aPp.
4 44'44..C~t 0P 404I 4.~0 I~N 117~
AEDC-TR-*78-63
0 0 0 0 0. f l AIIALIMll
,, :s-N4!-g 'ON, g ~ , g~~~
----- 4 N4 0 ~ ... N Ng 0 - O......-- -&Uý0 - 0.~~--0. ~0I.A40.6 f" mNNN -O- - 044MN-OON-Nf,
0 . 0 00 0 0 40.eC0U04N0 U -t40-- E~ 0 .00 0 0 0 0 () 0
0 -- a,-----0 40.0 4 v .044 NO N..O, P * N0 g 0, a m 0o NNNNNNae.JaAoEJNNNE'gN
4~0 1 .0 0 0 0 0 0 0 00 000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 N0 0
,I co c a* *coca C. 0 N40.0-.0 ~coo 0 'cc cc , I0 0
00a00000-0 NeN.000Ua 0~.P-..0 ~oU~M-4 4 4 N N 404 N 04 oP-.0 40 1i *0 24 - ~ . * i 0 4 4 ~ i 4
4 a -: 4tI 0 d ~ 4 ~ 4 U E f 4 ~ C . 0 0..- N NU No4 N4 ýN0 N~- NN' Nr~40, : .0....n~4. m *a1, O N40c4
4i O O M - ý-0 00 0 eo o oo o 0 0 0 000 r.0 _ NOo~ OO.0 m 00000000 0000000000 4 00000 000 0 0000000o 700 0
J * 0 000- n00 # _O 00 0 00000 00 0 00 0 0000000000.0 MG( N 0 MOO4..C 4D 04f 00-44- doe444;Oa; 0
0 K .ON .004.. 0M4 .0. .. ... 0 $ g .. 0g
M*eo c 4D @40-4 - C'NJ1 N4 M0 043 m 0 .0 0NNEQ0 ,0 m d .0.- 4MM4&'N M 0M.&44%ON 0,0 0 C ID coo N 0
I 44mý;0OC*4 .0-..0IwN 4E' ýC4~ 0.00. N-- N0O N
N ~ ~ ~ ~ ~ ~ ~ C A Z 0 aJj O 4 C .0 fl K10. - 0N 0.0 0 N N. N0O #
aK ItI .0.N 04i0 .N.t0 ~ M 00 . 4. . I 4. . 0: . .04.0
000000000 00 00 00 0 0 00 0 0 0 0000000000 0
ON CN4000 v000000000 000000 000000
mw ,C 0 N 0 N U ID 0~.0 f~ 01W2N U-.0 4N0 O4 004400= Ln
. 00 -..- N N4q .0NNNNNC mmaoo--N MCL
- - - - - - - - --cN tN
AEDC-TR-78-63
0000 00 00 00000000 0009000,000 0toNA:N07 _44:4m 4 m4 -6F # #Ný OAN9.062.4~4N04 11 S4N4%mo m 4ly. ~ U 0,
- a,.fN N4 t 4 mw fD 0 0004
* zzm.imF.F.F. FMF.FMMMNNN NNNNNNNNN.N NNNNNNNNNN Al
000000000 000T0000000 0000000000 0000000000 0,,,m,,,SI~~ t~,Sg NS,,,g *,0u N
S 00000 000 000 00 000 00000ý0000 4000 00 00
0 44.0 0*0 ZMAi4.-.0..OSflN 0fF-3~07-WOV A.
hi A F 9 4A M 4 - lo 0 - MM 00 , it I nI,,,t
o 4.PF..FMMMM mm;MMFMMNNN NNNNNONNNNN NNNNNNNNNN No~ ~~~~ @ @ @ o@ * 0 0 0 0 00 0 00 0 000 0 0 0
* 00 00000 000 000 ý000000000 0-- ~ ~ ~ 1 ;N.ONo0.N Md. om 6.0 In
4. Q....NFF 10U.0-.1. . . . ........ ... .. .. ...
00060-6S P .04.00Al40b :*N =N0F. 4 00002M NF 11% o0F.64MAi%06 "0 00F.66Uo Ao0Omfm6ONoO mA4-0NO!00' 4
F. & 0000 NM1114 !2 *4AN-:0- mIUoo.-. 6 d 6PO..-04z*Z fZM , ' a Z 000 00 0 ;:.llF4U 61N~ 0N (0. ON Ily1- cknAO 4a,~ 0&0 -0'
m0 x oe o oo e O e O~o 00...A N 0 0: :'Q 0 0 06NVA!.lt~1 t'
06PF10FO M N,-Ii00. 44 'IifI406 M"-00F.%f04P= 4 I
4 0 ac0. F.J4 .. 4.4N... 00004M -. NlIII.~.F01~~~ ~ . .oo ~ o .-... . . . . . .60 . .. .- G. . O..F .N0P. l:~9 1-- - -- - -- - -- - ------- oo oo oo---N M --------- 6.0
coo~~ ~ ~ N .It.',. .l . I
- - - - - - - - - - - - - - - - - - --.
4=:*N4
* 6Al4F.0Ad 04FO0.6... A-U--4011 P.-.619Ma9F4'- 0 6N440-.4'-G 0P-4NON:04' 4.F.04..4M *F9,.6..4
* ) 00 60. 4:... . :t. t . . .N. .4. 4P . .04 .6. -q% T-. ! . ..'t 0N - - - -.4 P. 4.'6.- O.- - - - -~64 - - - -~.4 -
olN4 ; ;Z* 00 0 00 0 : N M4'4 0.6 4 0 m .94 0 NIP-'.4.0-,4 a.0
S. A O ~ OO o oo o oe e 80 - - - A MMF.F.4 04-a (P:; oaDT ýNm C - M0 00 *d C h0gq
000. .. .. .. .. .. .. .. ..
N4N N 4M M
I N 4 4- o N~"o .06F.06..~~~~~fl66G.NN NNFS6.~r- - 9PM W NM -. 064N00mm
0M0FsI0N0M 0F0--. M.N0 4 0 60. .. P-119.4
AEDC-TR-78-63
.00
,,, eC2 1ON 1 ISJ N N AINN NNNN..ý:N
01 'C C-I~g maa c~am .=am09N N ~~ m_:ccc 0vZ CCo' OO C O OO
0M 0 f=44N-N'C 'cm&OO4NN'C4 NtzT. .. . ., 0 0 0 0 ItN 0!CO . ..L'.-
m 04N *0N N "j N 0N NNN'C1
v M LM oo o0 c o1LACco co o. -
o ý co CPm o d)moýa:- I- N0 19 0 N 91 w N N IU N N N o
a 06.,. a so m ll o Booleo 0 0 00 0 0 0 0cC c cC o 10E c'Q-.1N Q " oc co
*~ ~~ ~~~ m "': . a C 0 . N N0 5 1 4 4 . 90 0 N C C 4c--A o o 'C N 0C a. w41 . 4 'D.NN0'P mw'
U~~~V V: . .C . .1015-. .'- .''. .
N/ v N 0 .0''49 409C5'0 1'N'
Is a a e c0
4 ~ 0. * 0 N- - - - - -00 N'CCC4 9N.Pz 0. a. 0. - 2 m m i a s a m m m a a a
1 ' 4 c oooeeo00 000 000
4 o .z oc o o c c ococ ococo 000.J 0 4 F.N40P)0C-''N eC'09NN N N
2 C - C 40 , AwA900 0 P1N N0 .. N4-ND:-
a4 In N- 4m N 0 a N
4Do N N m m 0,1 a~v
o l Pa.S a C, '
u . 11 mC N ;oC4C'N~~eC5N9 'J
4 m N o3S N 0a a) 4 N mNN N N N .NN N N N N N
2 0. N Ncc .. N . . .V
U 0. 0 C N 2- - - 5J9
o acoao o u 0 0 0 0 0 00 0 00 0 000003
-m ua a Nm me 0 m90"a4d4N4 'D00549 wN'C
31 -o a, N N N. ,0 . N 1 0 00 4 N C 5 10 4 P0j Z 0 0% 4! No, oS Nm ,00 -. 4 1 waS4). . 0 -,A0 la It It t9 .4.C.' 1! . % .14~0 0 0 1,%. 1: . 1 9.
a,4 .c D za% rN- - -5 N I amm ma
2 w 'C 0 m3 N mmm am m a m mi a o' am20 m w9 0 0 0n 00 0 00 0 m00 00 0 0 a
LL _j x n m mm4 4 M0" e 'C In m m904CC'' m~'~ N S ' N 4
P- IN -i N .o ~ C 4'U O CSSCSO
3 0 c in w N -0.4'N0 NC N0N40'C *A0PSN M. N a. p. *o w ()
U ~ ~ ~ ~ ~ ~ N I- Om KM~ N 5555/5455 ..4~~~~ lb cc0 0 0 N a 4D a,ms a a a m a a a a l a
o ~ ~ ~ ~ ~ L a os- - - -- - - -I a, 3
4 0'C ~ 4 45C5440 ' 'NNN.PN mmNO 4444C- 0 4 0 'CO..PS'.. s-C '-'C-'C0MN
0 4 3 ' 'C * '. 4CNN.590'C 3'12.0 N4 N 0 N54
AEDC-TR-78-63
0000 0 -- - -MNO.NN 44 M M M444 441MP-41 e444MOMOO4* 0
2 ~ ~ ~ P f.E004 NN0 4N 4 N4 0l00 bO P.u4O 0 0
24 N1Pmom.0 PM 40Cllo 4-0. *0304ý0 = &.0 ON *i eP0 0
a 00 N-~ o m a 40 0 00 0Ji4 i.0ý00000000 00
00000 - - -NNNNI.Af .0 1PP144 4444 44
[emailprotected].. :01-...64000 .000:00iS. ::0-.::-:00a 04p:.0IMO oina 0i0dII4I. :.4fl0.0 cNFJ0QM:0aM 0
0 ~ 04.N644iiJ iflfiNoiOMOO - 0404f MN0AJ06000 03- *l.0ýi~0 4-.444U= tN 9tCON WI0 MN& 66O OCAN~ E 'a
c Em4wlcPO6N m-0 * 6!=4 0 C44 e.mtN- C !Eý04N
o0090000000 0000000000 0009000000 0000000000 e
* 00 0000 0 0 ,00 0 0 0 000000000& 000 0 0000 01
2 0000000 0000000 00000000ag m0!eOOoOeO 0. . . .. .. .. .. tN P4 00P.t 0P .4N t . . .4N 0
.0 # I.0- t 4#4o 4*#44g0 4 4 .. ~4 00 M IfmM mmmm nI
me o I . W .- 0 0 . -aMC0 4- InN4 0 MN 46 fl ~.P 4 0N4
.... . .. .. . .. .
00!00000000 0000000000 000 0000000 0000000000 0
maa WOiO oo001000000 c00000 v00000- ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ 0 4 NPP %p m.. NP.P 0. c.. ... P-i PPN-PN-PP
- P.O C0404 01 4..c4..P.4 ..8M4 04 in IV.P~4
SM40.P *U NJf0NN4..4y V UNOM~ IV-M44-P N2 g - - -
atf44 * * *4 4 * 4 4 * 4 44M 10 InMM M ,0,
0 OO .064M0 40040mP..40 j 0 N i 00N4P42 4 404 04 *46-P-4-4M P-0M:006N Mow0N644m 4
4 ~ ~ ~ : w a w O In44-2-PO NMMM 0.~4 U 4 P 00 P4*0. 0N00NQ~sP 7 9 10 * a M-Z 0 IV4 .44 0 A64 44 M 6
NO, 020 - M f",4~ -. N4~n OJýi 0 wa& 0-i2N
0 WNNM 4 MOON WN N N N 0NN ------ ----- NNNNN N
19 W-4-MNPM. 0-M6N N4 'm M..44 NO O 40 0 m2, 4N P 2 m4 M00tNt4 NO 44 N In4 ý!4w" m0uA1 on-3.00No MOM.4f 0N.6 4M 04 N64IP44 4 N-m -0
ow AO4M 444 m 444I44iý - o qIP-A 44PN.CoN m tomoN M 0-40 c4M. a"N 4 w 0P-4 Nt01 0i4 wO 4N.P4 MN m4,
00 0 000 a& 0 0- 0 0 0 0 0 0 0 0 0000000 0
000000000c 00 0 00 0 000000 0000 0000 00 04P4-04 NNNNP-P4 Nm mn60mmmm 4
44 N 4 440-It444 4-4 4 N -4
I. .~MN. 66 NP-6411P 04 40604 40N 40N41214
AEDC-TR-78-63
*::::0-*coo 0 0000: 0 0 :*Coco -. 000 f0'00 N0000a-----0~ 0 000 00 00 ý000 0 0 0000 00 00 00 000 00 OO:eO OO 00ýh. 00000000 000000000 000000000co0 a000a 00 0c0c0 00 0
- ~ ~ ~ I Oo oo oo %~o~ o 0000000000 0 00 0 0 0 0000000000 0o~~~~~ %co oc o C:ooo o Itooo o C:ooo e 000 00 00 00 . C .
0000000000 0000060000 0000000000 0000000000 0000000000 00I
.0=1 @-N 11 cf- P 4'b.0 3,N ~-NAP fS 0.00~e. N90i4O*I0* 04
o f 0~C O 0 M ** o ,C-Iý MO-O 0N'4 ---- in
-n J, 0.0 a0- a 0 N f -N11~.O- 05'11 4 0I
N. I NNNNNNNN NNNNNNNNNN NNNMM A~4 A ". m 44f
~S ~44~~-F1~J~...Ii *40 N02N..0 4P.NO4ONO !~N040 0nf.I44 V4.
N 4NON W' N w N w'NX&0 ~N.C M'v40MN- .4gN w o*ýCVN.C0 gON4pw NNNM"' E4im lroc4N--mec 0, a .-.. N N N rm 4 4 4. il 00L0 wc c0& c
o ftftrNNNft ('. NftSN lINNNN N NNS~ NN90mmmmm AA 4
ccocoeacoco 0oooooeoco ooooaccoac 00ooooeocc ceccoaccec 00
. 1 01 . .. 11M.-. .--. 4 ; ~ 0~gg~i ~ ~ 0-N40 o-o t-A 4 9 -4~* csO-Ztim N44wU Nt4 N4'1 0NC 4:90O~~ 0 N c*0 0
M4N40N 4-,I0.0 4 0N A N 0U N top w 04N..J J 0 0UI.0 -N~ oc. O - o NO a0..i~ N NU Ni 4 0 4 N N m O
--- --- -- --- --- - - -- NNENNN NNN JNJNVNNN NNNNNNNNNN NNJ
22
AEDC-TR-78-63
I %n a a0 0JME04'S0 me o c w ý0 (PI~~0 ~*f~.lCP00 N~4f
N;Z& 90ý: P.A- 0 t m = NM N3000 0 00ggNO-0 00000
~~~ 11d0 .0N~0 J MN 0~ 0.00 0. o : 0 0 0 00 . 0 000 0 000 :aooooooooc a000@0@ 0000 0 0 00000 00 0 00
000 0000 0000000000 000000 000009 : 0 0 0.000000 0000000000 00 00 0 00 0 00 0 00 0 0 0 0 0 0 0 .00 :0.0.0 0.000000
000 000 0000000 0000000 0000 0 0 0000 0000c 000000.00000 000 00 000 00 0000000000 00000C00000 000004;0000 00000000
ie ~ ~ 4 me ' I I 1 I IS II I I
o-PJI~~~af~~.IEn W.0I0O .---,0 wc0 0 0. 0 0 0, 00 o0000000* ve~a M** mo* *EEE E4E* *Em4 4 z* *Ecloo 44EN MM*NE*44 4 44404444
000000000000000 00,00000000 0000000000 000000000C 00000000
*415100n S.o E 0..iI N 0 0 g M-ffn0 0 (P. 01 f. N :Ný **E
WOO* M-~J0 MI-0 0ftE1- -. Mo00NSM0 ýW000.-..-NNM $E4 4E * 4 *-.- NNXAMM* *. 40MO&C.OO ý-P-P-000 0 WWW0000000 0100000M01 0, ,0,0 0, (p
V444040 4 40 444*44o : 40*4 444440 404*444404 44* 44400000ooocoooo ocacocoooo occoccococo cococcooco ooceooce
0 * I S d ' * U d4 0 O ~ - N . . M 0 4 4 4 2 a m0 0 i 4 0 0 4 U N A W : 4 ! - N .
MOM 400,0 N~~M0S 044000P 0 W Z f&Dy I A 004-0 0.4 N &0 -0 N0004.05-
sseo M M0" M. l 4 .0I~~o o 44.fJ4~ # SIt * k ul 0 t- o uo U nL j ,k nk - -0 10f.a
~9Ia MOM~Mf 4f4444 *4 40,0I 0 ;NM40,0J7m ; M0000
A 0 - # 0 00 00 40 00 0 1 0 'a 00 00 0
123
A EDC-TR -78-83
,,
o
a o
42
AEDC-TR-78-63
7TOOU 7 ,7OO.O . 7 ,77707T77
10001,cANiiAi~A!A
No oo
* ONA b MM NNAJNNNNNN NN N NN N z 't N 'r c NN no NNNNNI't
*7 7 7 N N 00 00 00 0 00 00 00 oe OO O 0 0000000000O 0000000000
u 00000000 000000000 00200000000 00000000 0000000
NN I.. i4OII Z0 40 00 7~~I~ a4-PP-04sAp.m N.EpNIN N.PM PK - N.0 gl-~P 4MS.sIONN .. ~eP.NU~N .N.0...4~ A.0UC
: 4 4 M M -0 550 *_Vc4 O 4e &I * NI .. 9E .0sNM * 0Ný.P 51 c0U M0
U 001,IS p~.q4 09 v 4mmP1Q~0 U2jW.0.j0Q N 0. 0U OW~ ~~~~~~ V7 .. . .- .40 .. 0 .4N1U .~M4P~ .. P.-10
-- --SI 0---- ---- A. L~ 3a.N'''~0
o ON~b~d~~fl15coo 0N . 0055 . goI60 4c 0 *r4CO -N f I c*
000000000 0 00 000 00 0 00 0 00 0 00 000 0 0 0
I ~~ .S0~ M e 10ZN ý. 1 ItE .012.60J It~N4I- . . N.5. It 0MMW.J ODO.M5 E 3CM CWWW WC ItC M 0 054V t WX M IDWVNC0 co co c 4 oo4%k
00 . . . . . . . . . . . . . . . . . . . . . . . . . .-------- -------- ------- 7
* -ý c * O N 0 0 0 00 00 0 0 0 0 a:& SO OC 00 00 00 0 00 00 00
; 4 6N C 0 .. 2 I0-1 -q7N4 2 4.0'V 4 M-.fU.flfJfU JII~f607 00 W00 ,N MNcJ0. S . -MS25-4O fU0 5N75I4_ N ,N0 Z 0 M J
.o0 * . .*-~ . .. . . .. .. -. 40 .0 4#ON..% .0Jjf. .0 E .0Q.f. A.I43 ~ ~ ~ ~ ~ ~ ~ ~ N N fm m mmmmmmm0f - A0 0 5 O I EO 4 J * 44~-7ml02IJ
2 144..0.0.4.0 .0~'s'E0 MN-00. 45~-0 M%0U.MUNN*
49 Cmftoot2o NS6 2 2 2 * 2O OPP - 4 4 4 4 IPt1
* 2 0E~0NO.. 00.rG.00NMY~E~~ *N.......... V *.415 4.MC'VN e *NN m l .y4,ZO N N N '-2-.N m -N4 NN0NNmo 4Nmm mmmmmm-44 4444444444
v4 NNSw NI P..0-P.c5 -. O N?0. N N 0~ N-NM 40~ N It m OaCmN ~ ~ ~ ~ ~ ~ ~ ~ ~ m 4 0NC.~M g..p.S47 f.O..00*MS0 S~ 0 N
* ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~~~~. J, NO 4dI~1 400.. C50C0 0 - fN.N ~ M~*emA.., , _ _ Te ~ ~ 7 ý ý l
A..ANANN 5LNNýpn.% ooANNM TMMMIM c 0MMMa __ !- N , MMIIVUI %ý4 P .
NK
12
AEDOC-TR-78-B3
1.1 14 1. 1.1 iy 1. 1.,NaN N NNNNNNN N NNNNN NO N f. NN eeafP, e IM~4
cococo00oaoooooocooaooocooaoocao 0000000000 00000
20 NC4;;. 0A 4 40 M440.P004 -NN..O..40 40N N l 0 :4.Om:
!v 4rv N 494440Ue4 omNe 0fA004 Co4 w~w @4NN C4&0.40~ 0
a a on in 44 P, Ny 0Z NN4 ~N N N NN44M 0 4 N 0- -00 N -. N - N4114 4 6J48~4 $4 1111111 41414441 110018NN - - --1184 N.7-..SIA
N:NTNW:~#NVN NN NNNNNNNN NN NNNNMOM mmOP3Immmmm** 4044000 0 000 00 0 00 0 00 0 000 000000000000 000000 000000,1 @OIOO O O 000 0
0000000000 0000000000 0000000000 00001000000 0000000000 000004W,0m00 m00 a *UN004A0. 0NN4S 04 0NN O0 0N' 4Ne . elf"NO & tono;Ot.4vo0N*60 O*N-04n:.40 ONN-40,11 0 404 N4 04N000.....N- 0.i004000N444040 N0ON40 440400.-..4 NN40mwoowN :44000030M 402N0
*0609-.N.N0N._44NV 0940 NN044 WN*4M44m.04 *4M4N0OONON
0000000000 0000000 0000000 000000 000000 0000
22ogeoooo oooooo 0000000000 0000000000 00000000 07 07000*044O... N~e ~0M40Eg..3 ~N0N40N~ 044N 00I I.0 0E~00 40.NN0
0440N-mNg44 '.m04 ;:Om;tJ ~M 0Nlo &M0440 4.ge00~ -4?-NNNN044 =*04-.44oz0.i0..0-. 44 M,.~N 0mv 44 :% zN40:zN0'lO .0"0 N :01 g*44
*N0~.440N @4444044N 4U 2N 0:404N M N.0N 44O&N000Omo NNM!2.0,op, C 4p4 004 N 40. e ~ N0~0 mw-.N= 0404 _0 0-'e- mm 4004N
mi--44 3,VA~ ON40210JON .04 -e 0000 Z-mm 400-.v ON4..04@4N0 0a~m.0N4 V b0N OF 44N-&wo mN a, f" M -' 4 0.1:1111!% V0:14171..4 ... M4
0000000000 0000000000 0000000000 000e000000 0000000000 000000000000 0000000 0000000 m0 000 000000000 000
0000 0000 00 0 00 0 0000000000 00 00 00 0 00000000 00004000-tNN -04 040~04 40O~NN.4 450 40M4 4 4 0.0V4e N4-NO:.@N44EO N 4 N0- - -~~CU~00 4 4 NON N04.. N 0N4moo.! 44400
4.4 NN 00 - S-0N -0 Moe04M4400- 0mOZ44004.. 4040a:440NN NN &
Wm:!N00N 40-N4 N 0 N1C000 0: Nm a0 MCP 0M 1,- ' p04 M M0 00000z-a4 4 4 4. N N N C 0 00000 M-0M000m0m @ 0 0 0 *00 00 0
0000000000 0000000000 0000000000 000000000,0000000 0 0 0 000000000000 400000 000 000 00 00 _000000000 : 0 0 00 0 0 0 000000 0 I *4 *4 * 0 . 4 0 0194It#40*-.39 t t 404 00U54* 4 4 0*4 4 1 . 4N 4 .N 0 .0 *4 * 4 4 404 0.00ao.4000N 0.00.0054000. 0000.0 ..0009: 9:::0~0 0M00 :: 0 00
0N44NNOCON ~ ~ ~ 4 400aO0 aN 0004 0-- 404 04 004N cc~~~20233-3,J3, N.4 N 4a 44-NN O. N-400 4-N 4 N0 00
NN4O0#0 _MN=N0 M40Nm4 000.-NN M 44 444 4400NNjI~O *4~4i44 ONNNNW@000 00 0 0 0 0 0 0 0 00 0 0
. . . . . . . . . . . . . . . . . . . . .
OU 04 VN"m 1 01-1 lOm NCN-...)4mmm e 0 4Nt4 -4 4 440MN004 4 440 n n0 Lnn L n in-.0 %nLnU t00
r 4N % .. 0 4 C4£JN0. a0 4 O -0 0.00.. O 400 .g4..P N. e- ~ ,4 g 0 !!-g4 44 0 4 . NI N ~ N4 M 0
~0~~4 04 ~ 04 C~'f~0~ N5..N~0 N0 4-~'0N1264N
A EDC-TR-78-63
SS I
'O'Sace IP.00
*S°0
49,0ii24P45 0
P.9-0D
00000000 •
=ococ*,. 00OI=
"04004090oooo
e e e@9,0,,
99120
AEDC-TR-78-63
P. tN '"-O WN COO 40 Ni~ 'r 0N0N V V.6 c.. &a !M*Niý0~ 0N4,z;aO: 4- - -
o~c a 0 0 0 e o o O e - - - - - NNN N NaJ0 0 0000000 oooe~ocoo~o 0000000000 OeoeeeOOOO C0000
... Jo *t to..... ..... ..... ..... . .-.
*~N 0 1 'A0 0, M O0 0 0 0 3 0
r a. .0 0.:a* N 0N04 0-4 ONMMM*4,nPN 0 . 4 0 4
*440
*0*000 *0- - - N a~aM4&~ c.040--N*4
----- ----- ----- ----- ---------- NNN
W*0 IC *M NP.440 *0 04 40 N *. . 4 .0 440 * ! M0 Z N ~
It.oo o~ o ~ oe o o 000000 * MM N
V , IL 0. c I 1..
4 N 0.0 - a 0 4 4 0 0 0 N 4 O 4 0 - - 0 .
- ~ .,0ý0 Q-- - - - -C.Ooooooo. 0 0 00000N000000000N100 0U~~2 .4.a4.. ... .... ... ... ... . ...
k0otlmoc*00a40Oo.t 4ICOM
b, lzj O ON:OC c: C!c: P 40 N 04 4
t O ... t
fAJ m fMMM MM fNNNN4 4yN NNNaNNNN r00O40O4444440
4jrq NNE
M0 0 0 aca M4aM M M O a a a aa a0~
a. m
. . . .. . . . . . . . . . . . . . . . .... . . .
wc :a1w m N.4 0 4 0 mf" N Co&&o41 O0 4 N N 0c c c
-~ ~ M h U *4 Zoo *4441N..0 *OO.. N:;0P N 4a P4a Z04~ 0 4a-Oc ON4
m m AM a mm aa a a AnN NNN NN NN NA a.NNN N 7 NNNw NNwmccc-j t- - - - - - - Z
4 h 40 ,, 0 0 49.a0N Ma4 O4 NOO.Z4 000,4
0.~ ~~~ ~~ 04MN cl2 .'MOMZ4 N N0 2 . ~ a Z 4 aa 9
0 &V w m C DC c a a aV ) 0, 0, 0 000 00 0000'r04 44 00 ... 4D 4) 0 n 44 ",.*4 44 '
Mo 0 m .-. 4 04 N4 ONO44,a4Nv*ON 0-4 a 0 O4 4.-O
-- - -N -.- tat-nJNN -N - - - 000w0000w0w4 0 aoa ao - sa~a~ a . o0 N N N N N NýJ#J N N N N'm
zi .Jp- -
Of 4; NONONOEl 9NN;4O1.;NN A N N Z4 z N AMN zz NNNNNNNNONNN
- g. a a a a.-aaa N N NNN NNN NNNN---------00 0000
2128
AEDC-TR-78-63
S* * * *0*.. ... ... ... ... ....... 044 T. ~
0 0 ss O0t
0 4t0 ---- : ---- n
0:00 00'a00000-.p. a.
N M 4
E E E E E N E o 4 4 0 00ooo oo0 N00 0 0 N N P PN44 c~~j gzcv
on000400 000000040in 000- 00 04 a 00 100000000 00000000000000 0 00000000000N
0no00; 0 A 0 00 00000 00,00 00 0 000 . 0000000000000000 "000 080000.00 08 0* p 0 N 1-0 M4 &- a
""'.00 0 . . .MNUov 40 0. 00 00 C 4.. 00-N*0..N0 NN NN NN N o00 00 N. 4 024mMN - - 0 0
0 U U 0
46-.44 *0-NNN 0000N0t0, 04 M.N-.OOUN.02 O4MMNNNNN NNNNNNNNMN
0 00000000000000000 0 0000000000000 00 0. 0000000 0.............. .... .. . .. . . ....... M: t.....0
- - - - - - - - - - - - - - - ---- - ~ - - - - - - - - --NN 0 0 0 - - - - - - - - - - 4 4 4 - - - - - --N Noco o 000 00 0 0 . 00 000000 . 90 0 9 0 0 000 0 0000 V 00 00 0
ft~ 0 N" NN mMM M99009 #990 N 4ý co t 002N! N NNNNNNN F NNNNNN - ý N WW NmeNN owNN
------- ----- ----- ---------- ------------ n
* 0Zm z , 1 MIA -. n I Tp . 0 o '
0~. NNNNNNNNNN NN.NN4NN50...N00N0...44 ONN NNN09 9NN - ----------
40000-M04990NONN NNNNN--n-- ---- 09000 00000-0 00 0-
0049N-0 0499 40400N0044 '9-OO S O 0490N .000 O Q M
N0 0 y 4.r 0~uj0N00. N04 .... ... ...N0N... N0 N 9 0 . N N . . 0 0 0 4 0 N . 0
FU F F tf NN( w t y0 NNN09 ow40N A f NN O 4.jf. N N fu m N 009 0 NNr NNeim0 N 0 yN N y NN t N00Nc044- r 0-oNO n y CeC0 I e0 .5 ý x 09NL t C j-C .4oo xNo # 0 & c0N & 00z
N000-0000000900 Ng ~ 4444 0 0 90NJNM m% N- N - N N N 0 NNlN N,ýN V - NN-N-N~~~~~~~- NNN-NNNN -N--NN-ONNNNNONNNNNONNNN-NN
09. N 0 N ac 't -N p 04 0 0c & c: c m m 0p Na 0, 0ý a 0ý (p0 0ý 4I -C40 N4 N 00 MNCC C w0"0000440909 ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ 1 cN--00C CN00ý4 ~MfO.~.0 005NN 02 4NN
008 0 0 000 000 00 000 000 000 000 00 500 000 5 C 250 CC0
) . n0 0 A A 0A M
000N00 714 '00 n -0 mo D o-N -. 4n t, j,0.0N0 N00 0 4C 8
0 --------- 0
N NNN 0000004 00000000 00 0004444 44 444444 4440 99
wo 000000.4004 0...OOMýcN m o m c CPO,-O~ 0N 04 N4N0. 940.00 010N-
. . . . . . . .. . . . . . . . . . . . . . . . . . . . . . .
AEOC-TR-78-63
W-9.0460N ~ ~ ~ 0_ 00 eo 4 Ea0a .- 44. :a NNN op-Ea~a-oana a,~~a~~tt aA- aE
tit- NM ZNNMN* N4mON0flJNNNN rOZ00999 PMM ME~
Mm 0N m. 0 % 0' , i
94 'a 0'.N
9 4 N N 4 4 N N - - c o o 4a . N 4 9 a p. 9 n 0 o . . a N 0 . a a a. -
4;201 xzn 2 2
4004 9 9 9~00 N000N.:coca 0 I00 E000:0 1100 's 0
j ai000 00 00 N0 00 00 00 00 00 0 0N 00 00 00 00000010aN
99 ~ ~ ~ ~ ~ ~ ~ ,N N ENN MU 4M N..0 M9 4N 09.4 0 N 0 0 .. M 9 9 M .4yrZNN .4pON4P.4.*I09 0. . .p- NE M-I' .P~ . . ~ '9E 1 .*... . .. . . ..
N N NNNN _aN NM 4V n 000 0009999994 4064404 4 i n
0,0 o o o o o o o 110 0 e O 1
MP 2 - Ma M 0 IPN nt 0 .NNN NzN a'01000 0000ca"N 0 000 -
N N 4 4ID 1.
N 00000999 z.It 9 9 4 4 4 C 4 4 4 4 4 4 4 40 4 4 4 4 4 U U 4 N P.
U....I~l4 e4 WW~TIM-.1-,.4 NMU.942 NN9 -9...MEU14 44'a ata
ON 040 .4 .0 14.00 OA 4 M-0N 1: 1
N a ! - -- - - - - t N. -- 1:4 - !_ .o!-.- 2t 0 l N . 9S -.- N~ 9 N1. . . . . . ..9N - . ... .. .. .. 4 N N. . .. .. .. . 1V4ý l :
C C r0= C, N V M t'.NN4 40 M af.N4
04,0 44 EN 09,6 NS44 M 0 -N**44P 4O4M040.E4N N.NE -N 1 N N 4 co o 10
0 C m. C NggF.
N44M NN4 .9NM . * NN99.- 9N~N09 0N4N 4NN 44N N4M~atfjU)94N4 9NN 4N4 0N~444 o9. .4. oO~-ON .. o N
4 O.4.. . "N'O04N4N~0N.4oNo4.,
4 .o. CNOMW 4 N 42 O -IIN
XN 0
. . . . . . . . . . . . . . . . . . . .4
C A,0 N-
O0ON~4 4~004 N0o0 o eN0~v Na0o 4 N M 0 O 0
!2_22_22 .- -0 4 9 9 O 0 - - N -1 4 -44-N-.-04-9-9-
N~J~4N ~~IN N 0P5N'N~JNJNPl ~130h
AEDC-TR-78-63
- ~ ~ ~ ~ :4 N 0Z 7 IO ON0l.0@ 00i0.0.0 0'*...0N. *v- 0.aP 011 0,
to W% % 0 a 00 00 0 0 c-ooo o 0-0- - - - - - - -a -N N 0 NNNNp 10 0 00 000 00000oo00 ' e o oo O O 00 0 0 00 0 0
a- 0N0 a . . .0 .000.e.oe.0.0000 d0OO0 0000O OONQN
0 a a,14V
2 w 0 aW so ma 000ace*0N 00000000004 i 000s00 440004,000444 000
ao- 0 a *-0 02P .4700N~~~~~~- -0 -0 -0 -t7ii~ P-.IN - .i -0 --.. l~- -O M 7 N .- -
o ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ o 0 0 I t 7 7 0 0 - l I P . - i i . P 0 - 1 . 0 i . o 0 - C * . i
-i .2 0 - p CO Ng 0loo Ai 00 02
S.-~ ~~~~~~~~~~~~ . .l .-va.- .i ..tN .-.P.Ol7f .0OUPON . t'4.
0Z~zI~ooo~roo~ooaoooooo o- --------- --
*I = I.. rP. N ep .0 10.M MP i0 ~~ 01 O N.0N 0itNO o"".- 0 4 coo 0 -0 0. NN
IL a 2 - v - ------ e0---------- 00----------000-------
a 10, ag OMNt; , .ii. CP.71 1S
~~~~~~~ 0 .10 0 . 0 ~ t I . C: .0 . N N ~ 0 C 0 a i ~ t~ . . .OOo4l a----------------------0OC~c*coo *coo 000 00 000
m O2
0 c 4! 0 iM 0tI4. 00.(PM~*& C. O&M0 0N1O IDi~fPOO. C - , , 0 D nkPI
a N. N0 0 11 0 100
a V vi 0 a1 0.0I 0 #i
ly .i . .a. .*. I.0M *7-Ct O0. .0t4l . . .%tui. .. . .iP- P.p. *7.OitPV?
m
-1 NON rye P.P 11liPNNNNNNO NNNN NNNNNN NNNN NNN INNNO NNN NNNNNNON N
a W UO N N* m1 *0 IWOUvi ~ ~ ~ n Ii .0 M a C N m I I Z z oL" 10.00 itP A0PN070oO.0NarS N 04
*~~~~~~~ ~~~~~ . . .. .0 .t . ..4010 .71PN-~ v .00~oP'.v ..1PM 7 N. .P0C:
N0 - 1I
4 -4 C0F C Nt . - m P. 00 mvhio *0,N 0f~ ' O 0M a * 0" D,0 4vPvM00-P.i.4.PPN 0 nP 010 o m0r-07.- 0NviPp .0.Ci
a v kn 0: UN Inknin U1 4401 00 01 0 01 01 .NO. .0i .0. .. P .t. P No D.0PNUS~~~ mm mmt N-N-P1 -PMM.NN-. - -00-0.OOQoaavia-1N,.-..P.P-
hi.J
VS ~ . . .. .NI Z0 4t0 P t.-.a...... .. . . .... . .i. vz~~ .0.0.. a..... .. 0...
-~ ~ ~ ( A a, Cc N N N.00 0 a7.P..0nCN iEm i NN.007. it 0 0 P_1N~00 N W P-C
fn do0- ,4N4
w CL D ' . .-.. . . . . . . . . . . . . . . . . . . . . . .O* N W PC 0 QWDCD C0 WV W W C N CD m
*----- ------ ----- NNNNN N .NN.. NN NN NNN NNNm
vi 2 0 0 0101
* AEOC-TR-78-63
-10 70 CC. M JM .IN:: a -NC'. 0, 0 ~ M N ;U'&Cmlciom~ ommorumfu% 0..N4CP.0 .. M4C0.-M oz44..M ' N4COM440OJ C0ý C
MMMon"I MAAM4*T*:Y CU, 0 UUCC CN -CO 0,0000000... NN-m-m44U
Nm 0 - a.
M0O, CC ~ M N00 CU' 4M4U C6O NN4C 0N~'~.0 0C6
COO
P. 0 a_ @4r e.C O ... '. _ C40 0F 4~ U'a O O .. . N _WCN MC M 400 n
eNNN emmO-MMM444.444444CC0CN60000 N 000 00 0 0 0 0 00 0 00 0 0 0 0 0 0 0 0 0 M000000000 coMo
0000060000000000 . ... 00 0 000 .2% 000 00 00 0 0 0 0 0 0t 0 0 0 0 0
0 10 60 -4 4
NNMMM MMM MMMMnM MOM NNNNN NNNONNNN NNNNN- - - - N
41No0 &44MON44-'6 MO-MCN-.000 --. m#.CON.46 N604U MCP00
0 44 4 4 .UU4 44 4 4 4 MMMM MM MOM' MOM MOmMO0000 *0 0000000 .0 000 000 ~00 000 00 .0 0,000 0 0*0 0 0 0
00000 1.00 0 0 .0 0 0 0 0 1 0000 0000 1 00 0 0 0 0N10 0 0 0 '
0000 .0 0000000 4-,00 0 m0000 .0 0 000 VlAP4O0O cm N0N000 0M~~~~~ z MMM-M MM MMM M MNNNNNN N - N N N N N! N N:! N N , N1 O NNNn NNN-W
--------------- -------- ---------------
V -A, C' ý ,04N' 1 zM- 2.0 i 0 N omaP44N 0L
4mNN NNNNN NNN NNNNNNNNN2ON NN N N N N NNN N N
M0 44 04 O NC CM.06 -6 0M 6 U'4M.0v m -0046.ý00.464N..
MMMM MM NNNNN N N-- -- - 00 C 0.00000000 00 0 00 0
I ON - -- - - -- 0
N i NbN V N N N N N, NNNON iN NNN N N yNv ON-N NN N N N N N N N N NVNY NNNNNN
Ny N N N N. 0 N N Nw N Ny N N N N N 0 N N N N N N N Nz Nz N 0; N N N N NZ N N NN N 0NN NN NmN NNQNN0N NN NN NN N N N
MO4NCZ CC4OC-OU'X6C NCmw .... 44 N .N 40 M U.ý.OZ N UN....4
.4MM N...OOC Ncco'U'U'4 MMMNN 1- 0 O0 0 V, , 4 .0 m 4 MM000000 00 0 00 0 a00 0 0 0 0 00000 0.6 0 6 C Co 666C 666C 0
.............. ........... :...........
- - -- - - 0 0 ,0 0 0 0 ~ 0 0O0,0 0 N00000000044N 000000 000m0~ It 00,N
0 0 ~~c a'0 U 0 0 N . O 4 N 4 N O N U . N 0 N N , N CN N4 U ' - C 0 0 C M N 0 0 - 4 C N N N
040MC0m o40---NNMMM44NMC 4-'''CCNO6600.... MM .J ,.N4T00...M4.0NQO
O~x - 0, 0 0 040 0
06CC.-M 0:06 U4-N0 0'* P004.-4MN.....00 066.-NCUC4MP--0 6NNN'..NN.~t NPC C C C CCA0j.UA C C A '44 4 4*4 0 44 440. 4 4 4 MMM
046 N4 NN N6 040 4N 4- '. 40N 60 .0 40MCI0 -00 000M C . .0U A O . N 6
4' 4 4 N- o ; CP 't m .v~J9 m m ' 8 00& ,0 wm 0z %
MM.C m666UUC n m M ` "'NC C6 NN NN 0N N'M 4 M 4 Q N6N770U'0 *..004600N
4PN.-0NUCM 00 L..:.0 LO66NNN c L LU' LL MNNLN.L. L - - - - 00..-6-N- -- -'-'-
0 O t O C ' N 0 m 0 00 U mw0 N,04 . C44 4 . .
OC:TN 0 NC UUU4 40000NNN- - - -0000 0 000V01066666 = NNNNNC Cv4 4 C UCCC C 4CC CC U w444CC C ' C O4 C wa 44 4 44 40 a,'U''''UUU VU %r..00 AU''U apt n % rPL "LPu
CP- ý 60 z .. M U40 0 600cý 7NM14U'4N60 w.NM4U'4N600 x ZNM40C'4N6m00 r .)M 4 0 a 40 0a444440 .00j,0 0.0.00 C0.0 ICt COC N NN NN NN 6666666660O m .00 0000000c
132
AFDC-TR-78-63
4 af in aN
~ *tP-P.44O.Mv VOMM.C0Nvoci *..$ 2 IP M ýIP *N4Cif0M O~ 4
00 -- -- -- -- -- -- ------- OM ---- w -o
Am~0 4 co .A M. N~CC 00 9I-04O ccP~ 404t .
*JINQ41NO4MO *~ ON4N N .. 0O:0MC z40-04; --. N 84
00000000 -- -- --- - - - - ------ 0oaC&
M200ONNN CONNNNNNA*N O0 0 0 4 0 1 N
- - - - -- - - ---- ---------- 00 0 0 0 -- - -- - 0000000----0' 0o &' U' U 41 0-
`-0N401M 0 14 4-4N.C0,41 m ~ 0 00 0 :1~-00 4 0~t 4 444#40.0000000ft'Oa *0 MCe I4U4l I~CS~ M I0fN 4~ - 0 -- I
aMIiiNN ell NN NN N- - - -0 0 0 0 00 0 0 0 -%.cv 1
00 0 0 00 0 00 0 0 0 0 0 00 0 0 .0 0 0 000 0A0 0 0
00 0 0 N0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 000000000nz
M00-4a Uý,cce eo OCN 2.eO CCIE WOO .OC O O CO . O O OC C C CC 0 C O N N
mc c a 10 x CCM N.CCC~ t QN 1 ~ t~f 4 0 0s- .. ........ .0o ~ 0 0ttt-- O C 1 1 I 4 4 .4 -M~ N N - 0 0 0 0 .i . W ' '
TMOM4604CC NxoN.. Z M42 4 12! NZMC I -Eg - 4.4f0 0s.l -
aO~O00 -- 000000000 4------- - ------- --------- -
N fv N v U
s-CC~~ 161CC C-
0C~- C~l4 0000. 0 00Cft,0 44141,f4'5
IN 3,00 000 C62'# -
c 0 c 0 coo4 !
- -0
244 ~133
AEDC-TR-78-63
-- - - -- - - - - - - - - - - -
a- N 00000 000 ooooo ooe O oooo0ýooe ~ O o O 0 0~ ~ ~ eg j~ gi~ e gs ee e e-i
jk 0 *'0 00* ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ 0f cmfu mN N 444 40 .. &00..eiMM~ MA 0M.~M M~ M Miiil I N N
*~ ~~~ 0 N1NI.J N N0 "2 0 0 0 0000 o 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
*~~~~ 0 -0 0 . 0 A0 0 C ) A 4N O W , . 0 4 N- 0 ' A -.6 A ¶ 4N O n O 0 N 'U~ ~~ ~~~~~ .0 o44 4 . 4 InS N N NN - - - - 0Zin ID NN
oj 4 zz j 44j9000400*01N. .. .-4~ S 54 *0 . . . . . 4-Itti .0004 V :V .1
AA W ~ 0 NJ M 5 M 0 MNNNN N N e Ow _-M0 (4p040 5-5 ON.a. : -404 4
IZ N 0 c *I*Ol0 Z 1i 0
0 ~ ~ ~ ~ ~ ~ N 4 M P4ONIO. OS AO0I 44 Ml5.flom...0NI wM-00404 00 5El4~. T1 . .. . .540 5A O I .0O00 54 N. .N..0. . . . It .% . N4.M60
-- -- - - -4A 4M N l N. 4 0 5- 0 - 0 4 A0 0 - -- -- -- -- NrYN u N N 4.0 f
_j N5. 4 LA M 00 M 44C0 t05aNA 9y.0 cc 4 N2( 0 M.l0M 0U C 0 . 0 ry rue D%A 0 0 IA, 0 a0- - 0'4 4 l 0l a--- 0 0 0.-e4, O 0 ul MO c, UN In
- 0 .4 ~ ~ ~ 0. 4 A4A5-45 M 0 0a-N 0 4410NW 00045.5 05-4. 10600 ftl... N515 0 2 0 4 00 000-404410 0 0 0 4 0 t"-00N I 40n --- -- n.0.'t . .4% CC:4J4NItN0 It 0.- j, 9. Ct4 0-' 54 ItN N, 0 4 1c IA0.....54 .....
IA W 0000000000 000000 O C CC O O o co c oo C c2 z- 0 N 00 0 00 0 0 0000 000 0 00 0 0 0 00 0 00 0 0000LA v f;; l
'N 0 1 N N0 In a45- 6 0 :00 0- -Nmoom'A44 0 IA1 50,0,0 0-..-.z M5 to W 0 0c00000000C- 0000 0--- -- -- --- --- N 4141N1N
a 0~ o0000000000. 0000000000 0000000080 :Oo c 0000000 0000
An c 2 W
N 0 0'SN-400I n,.V- 0 M ;!.DNO:1: 0 10 0515InO tfY M - 15-.0 04M tm* ~ ~ ~ ~ 0 '. N45-02.O. 40441 N000A 5. N0410 N a, 'n. 00In a a N )- 4 444s0 1 o0-.05-2 44-.05IA0154 mm-00 0 M0a 00
IA~-- _V0 N -VI44..e~ 151~51 4- -0000000 0000000000 000l, I t %UýI 0 00 00 000 00 0 000 00. 0 o0 0000000 0000000000 0000
* 2 I-.0 0 000 000 00 000 000 000 000 00 0000000000 0000
c, 0 W iii
Dm0 ,IN!!N4 i w . c cý nj 0 0- M,.05 47 V01 N5- 04 4 1 04D700 0 1604"' N N c ONNC440 g!'0, - 40 - 00575--0 10- .A55- 0.
15u 2 J'0.0N5-40 M (PO a0N0400 -4r -. _ -I0-155.0 00-c..-'NmNNN 4N N Nr W0o * 00000 0 0 0000
* 0 5 .
C5-N4I 4 4USC 5-25-0 4.I4-04 = 0045C S VM1'5C4500...15 CNCO4N4 M1 - o. a, 4 . ~0 InI 5-A.55-045 In-ONN 0 W~M0 &M
LA If u2 - 00-45-1541ý It mN0- M .E .002 0 4 0 15 M -55.a, I a OIOI2 r U. 0. N 440505- 0.M454.2!M O ZZ.0740N.205 ..504544N 4005
In Nt 0y0
'0 -0. _. a 5.4 .5 0,00 2-51 =2.CA~f In040.5s5 -1 M !557-0 o M g1W co 0 2 40-100 c oo ct5 o Oo04I15 00JA 41 55.. 4v :,l 0090 = 005.- 10-5
4 0 - IaIA05.00.-.s U'-0.155IA50155 0050.004000 4IAM1515515440 Coco5
-2 02 VA,JCID M 1, W N0ý .r 0 004.4N-50 ( i 07000--41451 N *I05-0 1010- :N MI w~~~~0 0 0 n7 a- 000000 vOo~ 00- -- - -- 1144 41414=141ý*~ I 0n 4 IA 0 00 0 00 0 0 0 0 0 0 0 0 0 000 0000000000 0000
- a ~ ! a. .2i .04N0~0 '5
504I075 50 7NU 4 . .O-.1- 7 i-.4.
0 04 2 4444451W le a 0000#0N-..C LP. 1560.Se 054154 0 M4Z oLfut0 0 n
N !a N0
0 0Nm;
0 ~ ~ ~ ~ 0 Z.; 0 N4-.00. N4050-. M M41IN0 4 4 U 044400CD4M. 041o-4~~~ 0 .4 2n 4l 5-450E4 n-452.0 u,0rsNS-. 045.q0.-5. OIAOIA t n nLnu
7N M. G0 0 20"Mvi.00 0-- 5 2 0 7 0 70-N415. 4A A .56 0 0--
134
AEDC-TR-78-63
07 TT 00.N..*?T?. 00OO4MMVT4 777.E45J 7.770P-
. . . . . . . . . . . . . . . . . . . . . . . .NNft =CmacN NNEENNý NNfg N I'4N N 0N5%~0 0' 0* * N 5 ~ I~~ , , s s s , * s u s e s i u;s s s s s m I e0.n'JN A ~ N%~N~ 0 aJ"~ AU N1AIf !NJUNNN . :10 N5NnNI'
A4~~I0U'S-E4~ W5.0000-N N'aC NONm ~-~~~UI%00 104* 4 -tV55 MMS5mAJNNNAJPJ NNNAJNt!
'990- 0~444000040 000 00N0 00o *0-Qof5 *eM0-60
000000 0000000000 eoOcooocao 0000000000 0000000000m 0000000000
W, N O m.f N Cf Mom 5oe @1.0ef..0 0E .O. .4'a~Mfl 10~ 0P
N 0 5 4 P V NMf560$U 29M -0~. *o!....0t _0.N010 6~M M0N~~6t .40 mm5 I0U .. U 6U M.N6-i0M *000
.A'sI5 . ~ M MN N .* i -. .0 N~N ~ . . . .. . % I
NMM MN~ON CJo IMJ54- 00 5 NN - -O00 F- *m 00 5-
----- ----- ----- ----- 000000000
N*o-p z7NP ~ONq4 0 0 N'555- N N'1mCE I Op 0 0N.66 00-N A 00 NNmmm 440 40 0 ", I ,,m meNN6 06c60606 .. wo 0 M:: a 0OO00 00
NNlyNNNNA NNNN NNN NNN NN wNmNNN NN N NN N NN NN N NN NN N N m N N N N m mcecooc0 ec e e c 0000000 0 000 0 0 0 0 0 0 0 0 0 0 0 0 0
000000 0000000000 0:00000 0 0a 0 Ooeeoeo 0000 0 0 0 0 0 0 0 0
00000 0 0000000 000000000 0 0 0 00 0 0 0 0 00 0000 0000
Oc-.00. 005S0NoONSN OCCOCCOOOD Oacc660600 00O00moacc 00000' N0600F0 N Pe0 mCy 4W -05~~ 0. m ~E M. O5-e.5 ;m P-5EJ..p 70 NW o OU%%9zm
N10-.U%~ *400'-ON 0-N , . c~ .NN f~a 0 C 0000N00
0.00 000000 -n!:- 4- --- fUN NNNNNNNNN- NNNNNNNNN N NN iNfUem0000 000000000 0000000000 00M00000 0 0000000000 .00000000
00 00 0 0 00 00 00 00 0 0 0000000 00 0 00 0 0000000000
50.~4.-* #~0 E0600 00000 101 0 -05 N00 0 55.N0SoOO 20-0 0002m!-
~501 N,44 ll..4E5-EMm' IM5M gm i 0NN N- sO O O v06 0 5--a- - - - - - -- - -- - -- - - - -0 000000000
0 0 0 0 0 .0 0 700 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 000 0 0 0 0 0 0 0
ON-5 G- c CN..c.0fU6 4-9010=1~44 0r9tZ-000.P. 6-00ommo 0-5o 5-.oommou04004-- . 0.isja *NM#4 -0-0N.6 000005-Co !N.o.C.5 ON N00
t . .O 0 5 0 5. 0o .- 0 ' 0 4 . . %. . NN--0. I I. . .1 .. 066:0 4 OO M I f C!o00..... IV* 150 va-Z-CaI Pt0.. t;*J 0. 0 N 50;0; W .t 0 0 -NM0 fl(P a. 0fU(Y. 505a, m
M04:6o 104VN00 cc N cc "0 LO N.. -N-C SO O *05P-10- omON0 N. *O O 006M~vw!0m owwN-61 -. 0-NION 50 .ON5OO~60 NP- 0 *600NC t.
NO..06 150.NOONN¶0 000..-...oC0 N550N0SSUS0ow0S0N *55 0.4
000000 00 00 00 0 00000000 000000 000 00
000000 000 0 0 0 00 0r0 0 00 0 0 !0 20 0 00 0 0000n000000&
., 00-.. *0.153'00J5o mVC 0~0 00.J05'0 ,J0*, 0f~ 000 00
..NN O P.N 04 0 7 cov. o 0 M5N0!60-. 6.0N SNoom4 oOOfl
f~~r.,ooo~~~ Mooncco. m-,, oo o oo c 0 0 z i 0 Oeo4o
055N~ PECN- omS ENam-,#"O N6040NO 060 N5CN-..-NP...- 0 45500
-55 44 4 om0 d 0N50.s Do .f0,01-5*5 -00-f0 004N-0.00wg5cm 0 055 a l00550
135
AEDC-TR-78-63
o 00 000 00 0 000 000 00 00000 000 0000000000 00000000a
000 000C0000 00 000 0000000 00 a "7 No494 A 0000000I 0I t .N-t . .~~.. 0 .0N --in.1! .C N - . 0: . 1 ... .N4 .0.-.N
&&azwQ -OZOO,~f44 *MNN N7777 W04N.-.ým0c
00000000 00000000 00000000 00 000000 100000 100
.q 4 -O4 N N4 NNISJnl * 0-'ý4-.CN .AIQ "00 -4 N-c0 0V- 0 0 w w00N ~ 4m-.P. UN.004~ pN- 00~ N V 4 ~ J, MI
% ItlA 4N -~. N ... 4 Vt .~0:.j f.4o1 .1:oINi1011 co'. I-
0 ~C~ 00 a M004NOCmc-.N 04 NO*N 4CSm*%OO0 0,4 4Ný 44.0 01 .M , M0.4.4gNo m 'n 4 On6 :1. v 0,-S 1 4f N 000
0000000000 w 00000000"' 00 00 0 00 0 00 0 00M 00 00 00 N
00000 0000 0 00000000 00000 000 0000 00000 0 0 0 0 0
:0N .C E G -C 0 .-- ' .0 ~ V .%f0.-. M . 4..0!-tNt-OcUI t :ItI!O 0n
444 .NO 44404440 0000000 0000000000 00000=0000 I
.. SOAN4. NJ 140..O ý12O 0004.. MO.N0 - -- -S.,00*NO*040 004 0 . 01-0,0N IC N M ~04 M-00 0
0000000000ca 0 0 0 000 0000000000 000 00 00 000000000 0
M ;;:t~S. NM-.'0 ý M CY N N w n MN ýM- N0 4 . 40. 0 'aNN0 O 0 2 In'0 1~ m~0 -.- 004 NC;OM0S~.2~
NMMMIM 'WMMn NNN N NNNNN .NNN . N.. ------- NN N N00 0 00 0 0 0 00 0 0 0 00 0 00 0 00 0 00 004'000 0 00 0 00 0 00 0 0000000000I C : t ! ý ! ýiI 00 0 00 0 00 00 00 10
NNN NO.00 001S. MM M MNC-C N4NNC.. NNN0N- NM-0
0000000000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
I!C!1!1!1! 1 1i 'C C C J' 0 j-0.. OOOQ4Jl4 C 00 -.M~4 N0. - O - ¶ -S.-00 O'. .. C .M 40N 0ý4 .N4 1 J.0 4NMNO: . 1 -O. 4'44..0 .
P.00M4IMCM0~ ~ ~ ~ ~ ~ ~~~~c OtN0~0004 4C. 0C40.~M..00N0N4'~-.~ 4N0S.OMNO W -00 4. 4 MN - 0 0000000 0
0000 000 0 00 000 000 0000 000 0 00 000 600 00.0000000
0000000000 ~ C.. 0 M .. C -41U'0 , C 00 P
0NMN40CU~~~~~...M 44IlNNS.c0MC O4Iv0-.N40VvC 0JN1 M4.-rI M .N O4M S. 4 N4 0N C~C~t 4 40..S.-, N.. IM00CC04N.-0~j,0~.~ MMNN .C.-. . 0 M 0 N . . 04 - -0 N0 *
0N N . .M N -0N S~N 141CO S- -.. c,-. 0 mMN.-0N4~.Ii I S.NM.-CM. K0% 404.. .N0d0CNI 074 MlMeS..0 N 4.0N N &~~444 MrJ0J~~J .f0 ý00flJ 0
I-.".o...........000...........
40N00 .OOM INNN Ir. a4N .S-g 44ý 00 DwNo 0o
sI¶C0C'-T~~~~~0O w0 NN M 4&0CS-... 4000... NNoNNýom~4- - --NNNNNNNNNNmm NmmNN N NN N M M ~mm......... ......... .........
00-'NNM4~ ~~ N30. r,-~~00 N-- - mNNNIJ N NNN N N N NN
0-----------N NoNN oNNNNN z4 NopggN "
.40 I c0.o " .It. , . %0 .44 1 4 t0 :~PS.4 0 C S. LI . %%.- I' N O N S.%0 .dc1 mMII 'D. -. 0 0 ,0 0 N c4000 Osc.-. & 400, N M u a I Cp0...am.
-..-.-..-.-..-.-..-. . . . . . .
"1 114G . I fl C -N 0. 0 .. N O0 M O 0 1. 4 S- 0 M f : v t = týoa t
136
AEDC-TR-78-63
0000000000000000000 0000000000 000000000O0 00000
oooooooo oeoacaaa MON oooaoo oooo4o ooooo
: Ný!0UN 0N2:~d~4 ON:44.g 0O.4 ;!4O:4NM 4, 0td! LA NY 0
N lo M 44 044 0 o A.4Na mm Inz lo 4 .N M.4 0 0- . a m4
* .0 e P4N.NE -N..N.NNNNNN NvNNNNNNNv v--- ----1 6 *# I , a i ut,,i~ ,m g , su1g s g
000000000 0c0oo20 c2 0 00 0 000000a0o 02 0000000 0 00 0 00 00 00 0 00 00 0 00 00 0 00 00 0 00 00 0
* hCO U:,ft4~~NO40 mcU* **4 ic0.4tP 44...-0M0o40 0400 4-oomCCIDM 13; MI -MOO M 0 4) D Dit 0 f0 ft'.A ~~ *40u M 4MM NN 1,0--0s--
4 09 9 9 19 9 ; ; 9 z ;; f9S Nf~N. NO-4 4 0 '.4 P N I. .Y No N Nf NN NN.S~s N N vN N NO-..0Nx cc 55 O P- O V Co.lsi M~. 770.f 4 4 s - N 4 . .. E 700..04P0
o4. IS0C 042 c am~ SSS NNNN--- -- INry
l af- m s I g g , , NI I S S I S~ g g g s e , , m s , i , g , , g I iIt .OO~O C .O ~ OC O 1O OOýO O3w:~O "tItI 1 O %. ... O0C4
51 % C- 44 15.10 CN N~ 4. 4 - 0 5N4 .-..- s15009...4 N0UomNU oNEhi ~ N5N59445 O N.N- C'm om 9loN04 04.:.O590EN- co0'c4-oNý
0 w m N0 -4N 15151 .1U5 4.59 NN. -. 04.44N 0 -0m15 044U49f% '!-- - - - - - - - - - - - -00 .0PP C f.U4 .M 9 . Nf...ti ...N
4 nu n4o1 .!! -------- N NNNNmvIz
00 o
"" t0 C 1 I1 ý C :0% l: V;4 4 m= 4*N -- - 4-0 5955N09.50 9 N-4-4 4 1452-9A0-0 4Nv5
O - 59 .C 0..P.0 C 0N -. ~P.* *0 40 4. 50 0 00 59 0 013754
AEDC-TR-78-63
faoooaocoa too
N 4 N4o4Veco 42;ON 0Q44E5@ me
10 -. 9m4EmmiE4 o
m I *SN.ANu0. N 0 N
4IL~4t nN P ..
O NM*700o0o - coo
W000000000 00o
NNM.I-.e..4138
AEDC.TR-78-3
NOMENCLATURE
A Area
AC Exit area, inviscid contour
A* Sonic area
a* Sonic speed
C Factor in logarithmic skin friction law,Eq. (77)
C1, 2 ,3,4,5, 6 Coefficients, Eq. (35)
* C. Ratio of actual mass flow to that if R were
* infinite
Cf Skin friction coefficient, compressible
1Cf Skin friction coefficient, incompressible
C Specific heat at constant pressure
p
D1,2,3,4,5, 6 Coefficients, Eq. (37) J
Ratio, Cfi/Cf
I Fn Multiplying factors, Eq, (97)
Ratio, R /RFR i c
Gn Multiplying factors, Eqs. (94) and (96)
139
APEDC-TR-78.63
Ratio, 5*/1
h Heat-transfer coefficient
K Streamline curvature
In Natural logarithm (base e)
log Common logarithm (Base 10)
M Mach number
rm Exponent in Eq. (90)
N Velocity profile exponent
n Distance normal to streamline
P1,2 Factors in axisymmetric characteristicsequations
Pn Coefficient of 6 at nth point on contourin
Pr Prandtl number
Q Factor related to heat transfer, Eq. (91)
Qn Coefficient in momentum equation
q Velocity along streamline or, in boundary-
layer equations, velocity within boundary
layer
qe Velocity at edge of boundary layer
140
-. • w ~-.- T -- ...- .~, - T-- T ' -.. .. . .... . . ..
AEOC-TR-78-63
R Ratio of throat radius of curvature to
throat radius (half height, a - 0)
2 2Rg Gas constant, ft /sec R
,R Reynolds number based on 6, compressible
R1 Incompressible Reynolds number
R0 Reynolds nuitiber beed on 0c, compressiblec
Re Incompreauible Reynolds numberi
Sr Distance from source
r Distance from source where M 1, used tonon-dimensionalize distances for inviscid
calculations
rV Radius of viscid contour
I
SR+ I
,9 s,t,u Cubic integration increments, Appendix B
T Temperature within boundary layer
T Adiabatic wall temperatureaw
Tc Reference temperature, Eq. (87)
T Free-stream temperature at edge ofe
inviscid contour
141
•I.. ... ..... . ,
AEDC-TR-78-63
T Wall temperaturew
T Wall temperature at nozzle exitwD
TWT Wall temperature at nozzle throat
u Axial component of velocity, normalized
by a*
v Normal component of velocity, normalized by
a*
W Velocity along streamline, normalized by a*
X Ratio in Eq. (36) or (38)
x Axial distance, normalized by yo in transonic
equations, normalized by rl in .inviscid
calculations, not normalized in boundary-
layer calculations
y Normal distance, normalized same as x
y Y Throat half height, used to normalize
x and y in transonic calculations
y* Theoretical throat height if R is infinite
z 1Function of x in transonic equations, or
distance normal to contour in boundary-layer
calculations
142L •,= •!• •: :•:• : & •:|!J • W •- s~l i .• '•.• ; ,,F • .• ';.•4.,.2iL.•.. , ..... ;,.. ,:L•..,•.,;,'•,al,.1;,.;,,,ili,• ;',' •'•"" ::' ' i . ., ' . , .,, I,,; ,
AEDC-TR-7843
a Mean angle of right-running characteristic,
or factor in temperature distribution in
boundary layer
Mean angle of left-running characteristic
A Prefix to indicate increment in value
Specific heat ratio
S~Boundary-layer thickness
6* Displacement thickness in boundary layer
a* Displacement thickness when boundary layer
is large relative to rw
6t Incompressible displacement thickness in boundary
layer
Distance along left-running nharacteris.ic
n Inflection angle, radians
e ,Momentum thIckness in boundary layer
6 Momen'tum thickness whon boundary layer
is large relative to r
c Compressible 0 for flat platec
0 1 Incompressible value of 0
0 k Kinematic rementum thicknesa
143
-I- ~ -. -..
AEDC-TR-78-63
0 Value of 0 at nth point on contourn
K Constant in logarithmic skin-triction law
I2(L y + I)S
Mach angle, sin (1/M)
PC Viscosity at value of T c
Viscosity at value of TI PW Viscosity at value of Tw
e e
Viscosity at value of T
Distance along right-running characteristic
Wake variable in logarithmic skin-friction law
p Density within boundary layer
Pe Density at edge of boundary layer
a Zero for planar flow, 1 for axisymmetric flow
4) ,Flow angle
•w Flow angle of viscid contour
Prandtl-Meyer angle
'44
AE2DC-T R-78.83
SUBSCRIPTS
1 Values at point 1 on right-running
characteristic
2 Values at point 2 ona left-running
characteristic
3 Values at intersection of characteristics
A,B,C,DOE, Variables evaluated at points on Figs.
F,G,I,J,T 1 through 4F
second-, and third-order approximations,
respectively
OTHER NOTATION
4/ dx:
OUTPUT NOMENCLATURE
BETA Pressure gradient parameter
2 8* d ci d
YM2 P~Ci
C(Y) Coefficient of third-degree term if throat
contour is a cubic
C(YI) Coefficient of third-degree term if integrated
throat contour is a cubic
145
....
AEDC-Tn-78-63
C(YP) Coefficient of third-degree term determined
from slope of contour
D2A/DX2 Second derivative of boundary-layer correction
evaluated at the throat
D2R/DX2 Second derivative of corrected contour
evaluated at the throat
DA/DX Slope of boundary-layer correction
DELR(IN) Boundary-layer correction to inviscid contour
DELTA* 6* from Eq. (66)a
DELTA* - 1 6* from Eq. (63)
FMY Bracketed term in Eq. (61)
HYP/YO Value of hyperbola with same throat curvature
ratio
IcY 10' [C(YI) - C(Y)] for Point 2
INT.Y/YO Value of Y/YO obtained by integrating contour
slopes starting at inflection point
KCF 1000 Cf
KCFI 1000 Cf
KCFS KCF sec cw
146
. ............ x
AEOC-TR-7843
lO 1000 dO/dx
MSS Result of mass integration along characteristic
EG or AB (measure of accuracy of numerical
integration)
PAR/TO Value of parabola with same throat curvature
ratio
PE/PO Ratio of static to stagnation pressure
R(IU) Ordinate of viscid contour
R•ASS C1/0+0)
RTHI Incompressible Reynolds number based on
momentum thickness
SKPP Second derivative of Mach number in source
flow evaluated for BMACH
SMPPP? Third derivative of Mach number in source flowevaluated for BMACH
THETA - 1 0 from Eq. (62) used in Eq. (61)
WE Velocity ratio at Point E (Fig. 3)
W1 Velocity ratio at Point I (Fig. 3)
WO Velocity ratio on axis at throat
SWoFP Third derivative of throat velocity distribution
147
BEST AVAILABL[ COPY
ANDC-TR.7-43
WRPPP Third derivative of velocity ratio in source
flow evaluated at WE
w0 Velocity ratio on vail at throat
fV
148BEST AVAILABLE COPY'