Sivells
Sivells
Sivells
78-63 COPUE
LEVELEI1,
FOR SUPERSONIC AND
WIND TUNNELS
J. C. Sivells ARO, Inc., a Sverdrup Corporation Company
VON KARMAN GAS DYNAMICS FACILITY ARNOLD ENGINEERING DEVELOPMENT CENTER AIR FORCE SYSTEMS COMMAND
.ARNOLD
Prepared for ARNOLD ENGINEERING DEVELOPMENT CENTERIDOTR ARNOLD AIR FORCE STATION, TENNESSEE 37389
D D D
JAN 8 1979
79 0
BEST AVAIL.ABLE -COPY
"4.
NOTICES When U. S. Governiment drawings, specifications, or other data are used for any purpose other than a definitely related Government procurement operation, the Government thereby incurs no responsibility nor any obligation whatsoever, and the fact that the Government may have formulated, furnished, or in any way supplied the said drawings, specifications, or other data, is not to be regarded by implication or otherwise, or in any manner licensing the holder or any other person or corporation, or conveying any rights or permission to manufacture, use, or sell any 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 sense as an indorsement 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 National Technical Information Service (NTIS). At NTIS, it will be available to the general public, including foreign nations.
LThis
ROBERT W. CROSSLEY, Lt Colonel, USAF Acting Director of Test Engineering Deputy for Operations
UNCLASSIFIED
REPORT DOCUMENTATION PAGE
BEFORE COMPLETING FORM
(6
nal
epaDc
J.
a Sverdrup
16.
fi
1,SCRT
LS,(fti
eot
IApproved
distribution unlimited.
17, L3STRIOUTION STATEMENT (of tho abstract onteted ir Block 20, It different from Report)
Available in
DDC.
19. KLY WORDS (Contitnuo on reverse side It nce...ry and Identify by block "limber)
boundary layers
* *..hypersonic
A computer program is presented for the aerodyna-mic design of axisymmetric and planar nozzles for supersonic and hypersonic wind The program is the culmination of the effort expended at tunnels. various times over a number of years to develop a method of designing a wind tunnel with an inviscid contour which has continuous curvature and which is corrected for the growth of the boundary layer in a manner such that uniform parallel flow can be
DD FOR
1473
~~~ ~UNCLASS 7
IT
1
%j 01
2
0ow
2
1
UNCLASSIFIED
20, ABSTRACT (Continued)
The continuous curvature is which has first and second derivaexpected at the nozzle ex-t. * of -'elocity (or Mach number)
achieved through specification of a centerline distribution tives that 1) are compatibl* -Lth a transonic solution near the throat and with radial flow near the inflection point and 2. approach zero at the design Mach number. The boundary-layer growtb ia calculated by solving a momentum integral equation by nuws'lcal iaitegration.
YEJ LEV
i, ,oilii W
........... ............. IT ............... Slt. AVAIL
AI
D D C
JAN 8 1979
11B
UNCLASSIFPTD
,,,.
--
AE DC-TR-78.63
PREFACE
The work reported herein was conducted by the Arnold Engineering Development Center (AEDC), Corporation Company), V32A-P1A. Air Force Systems Command (AFSC). The results of the research were obtained by ARO, Air Force Station, Tennessee, Inc., AEDC Division (a Sverdrup
operating contractor for the AEDC, AFSC, Arnold under ARO Project Numbers V33A-A8A and Elton R. Thompson. 1978.
The author wishes to acknowledge the assistance of Messrs. W. C. Moger and F. C. Loper, ARO, Inc., for providing thr. basic subroutines L. Shope, ARO, Inc., provided for smoothing and spline fitting, respectively, which were adapted for use with the subject program. Mr. F. technical assistance in the preparation of this report. publication of this report, the author retired from ARO, Prior to the Inc.
'I
AEOC-TR-73-63
*l
CONTENTS
1.0
INTRODUCTION ...................
......................
2.0 *
4
.................. .........
.................
...
10 15
20
CENTERLINE DISTRIBUTION .................. BOUNDARY-LAYER CORRECTION ..... DESCRIPTION OF PROGRAM ....... SAMPLE NOZZLE DESIGN ........... ............. ................ ..............
....
24 34 38 41
SUMMARY ..............................
REFERENCES ............. ......................
...
....
41
ILLUSTRATIONS
Figure 1. A Foelsch-Type Nozzle with Radial Flow at the Inflection Point .............. 2. .................. 6 7
, ,
Nozzle with Radial Flow and a Transition Region to Produce Continuous Curvature ............
Nozzle Illustrating Design Method of Ref. 13 ,
3. 4. 5.
8 9
...............
Relationships Obtained from Cubic Distribution of Velocity from Sonic Point to Point E for Axisymmetrin Nozzle ...... ................. ..............
=
....
18 19 23 28 29
6. 7. 8. 9.
Limitations of Fourth-Degree Distribution of Mach Number from Eq. (39) ..... ... 1 . ... Characteristics Near Throat of Nozzle with R Variation of Wake Parameter, A, w.,ith Reynolds Number (Incompressible) ... .......... Variation of Skin-Friction Coefficient with Reynolds Number (Incompressible) ... ........ .... Variation of Velocity Profile Exponent with Reynolds Number Based on Boundary-Layer .................. Thickness ...... .........
3 h 4;
10.
...
30
A EO C-TR-78-63
TABLE
Page
......
. . . .
. . .
39
...
45 49
. . .. . . . . . . . . 0. CUBIC INTEGRATION FACTORS. .. . . . . . . . . . . . . C. INPUT DATA CARDS . ....... D. COMPUTER PROGRAM................ . . . . . . . . . . :"NOMENCLATURE .......................
.52
61 139
AE DC-TR-78-63
1.0
INTRODUCTION
two
Early supersonic nozzles (circa 1940) were planar for many reasons:
state of the art was new with regard to both the design and the fabrication; 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 used for all the Mach numbers; and the relatively low stagnation temperature and 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. Ref. 2, 1, laid the foundation for determining Foelach,
flow in the region of the inflection point was radial, as if the flow
1',-
came from a theoretical source as illustrated in Fig. 1. The downstream boundary of the radial flow is the right-running characteristic AC from the inflection point, A, to the point,
C,
the design Mach number is first reached. The flow properties along this characteristic can be readily calculated; and inasmuch as all leftrunning characteristics downstream of the radial flow region are straight lines in planar flow, the entire downstream contour can be determined analytically. curve. Upstream of the inflection point, it was assumed that the source flow could be produced by a contour which was a simple analytic In the Foelsch design the Mach number gradient on the axis is This discontinuity produces a disdiscontinuous at the juncture of the radial flow region and the beginning of the parallel flow region. the theoretical exit of the nozzle. continuity in curvature of the contour at the inflection point and at
.-
AE DC-,TR-78-63
S~Flow
Uniformr
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.
AEDC-TR-78-63
Figure 2.
Nozzle with radial flow and a transition region 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, dimensional stability becomes more difficult to maintain in a planar nozzle. Therefore, axisymmetric nozzles are used when elevated stagnation temperatures 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 particular Mach 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 Inasmuch as the assumption of method of characteristics (Ref. 8). 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 10). Such discontinuities have been eliminated by 11, and 12; here, an axial distribution thus Mach number gradient on the axis produces a discontinuity in curvature on the contour (Ref. the design methods of Refs. 10,
of Mach number (or velocity) between points B and C (Fig. 2) introduces a transition region between the radial and parallel flow regions, radial flow values at the beginning of the parallel flow. Fig. 3, the upstream boundary of the radial flow region is flow angle is the same at points G and A. gradually reducing the gradient and/or second derivative to zero from the As shown in a left-running The
characteristic from the inflection point, G, to the axis at point E. a general nozzle design. As described in Ref. 12,
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 distance were calculated by an adaptation of the transonic theory of Hall, Ref. 13, or Kliegel and Levine, Ref. the axis. Inflection Reglon 14. The upstream limit of these cal-
M *
SI
Figure 3.
6
F
-
.'*. -:
..'. U,,.J-
AEDC-TR-78-83
This characteristic is
retical location of the throat and the intersection of the branch line
with the contour was a region which was not calculated but which increased in size as the throat curvature increased. This gap in the contour has been eliminated by the method described herein which utilizes 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 With this latest improvement upon the method of throat characteristic). 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 standpoints. 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.
After the design method was developed for axisymmetric nozzles, was adapted for planar nozzles having a prescribed centerline tribution of Mach number (or velocity). dis-
it
The current design is considerably different from that of Ref. 5. method is incorporated into the computer program included herein. As an option in the program, a complete centerline Mach number distribution
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 the zero (i.et, the acceleration starts from zero at the beginning .-. 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 twodimensional 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) R
+
. R2 . ..
10
R3
.Y.
(2)
10
AEDC-TR -78-63
where
2 +
12
Lv (y,z)
L,( L.
v+(y,z)
2
+ v,(y,)V(
.(
3)
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 Kliegel and Levine in Ref. 14 extended the applicaaxisymmetric flow.
bility of Hall's axisymmetric solution to lower values of R essentially by making the substitution S-1 + 9-2 + S(2) and (3). (4) . .(5)
where S
R + 1, into Eqs.
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
n=
1+
_1_
3S d u dx/y.
(14y-7 5 ) 2 270S
4464)
-(32y
+ 87y
540S2
561)
4S
2 1-4y -57)
288S
(8)
Sd__+__ dX/yu
-(64y2
as
+ I 17y 2 1152S
1026)
(9)
;!:
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 On the axis, y
u
=
_____(10)
0, v
_i 6S
+
0,
- 7767
429y_ 4320S 2
,Y,"
36S
(,&
-z'
S1 -
(2y 2
33y + 9)/72 + .
(1)
5+
+& -1 t Y
++
()2
(.2~
(--A)' Yr,
(4Y2
57y
27)/144 +
(12)
Because the sonic line is curved for finite values of R, the mass flow through the throat is reduced by the factor CD (discharge coefficient), which is the rat 4 .o of accual mass flow to that which could flow if R were infinite and the sonic line were straight. For planar flow,
CD
I-
[i
(] 13)
I12
12
AEDC-TR -78-63
and,
y 96 S'
[1
24S
.754y2
28EIOS'
-757y +36U.
(14)
The flow which passes through the throat also passes through the sonic area of the source flow which is at a distance r
In planar flow,
0(
or oY 0 /r 1
=
,/C 1
(16)
(17)
yo/r-
2 sin (1/2)/cD
(18)
In
the
value at x - 0, y - yo
Eq.
(6),
The half-height'~
Inasmuch
or radius, y0 , is divided into 240 equally spaced values of y. as the characteristic is right running,
dy/dx
tan(
where
sin
=
I 1/M
(20)
Also
W
=
+ K-l
y+1 sin
=
Y+-
M2)
(21)
(22)
v/W
and
o + ddrdo
'""a
d
y
(23)
13
--
AEIDC-TR-78-63
Sde
dx/cos(o - t)-
dy/sin(o - 1)
(24)
(2I -
2 - tan-1 (I -1)
(25)
Equations (19) and (23) are the characteristic equations and are solved
by finite differences. If at point 2 are found (y is
S
the values
x2 -
tan(
(26)
A(Y2
-Y)2
+ (x2!
x1 )2
(27)
V0 2 = 01 + 01 - 0 2
At the starting point W is
VI
+ - V 2
(28)
0. Values of v2I
evaluating the term in brackets in Eq. (28), the ratio v/y is defined by the 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. axis of symmetry.
It istic
(23) on the
may be noted that the value of W as calculated from the character,value from Eq. (21) differs from the value (u2 + V 2)I12 calculated
from the transoutic equations, but the difference decreases wT:Uh increasing R. on the axis, For the final point of the throat characteristic which 1,us the value of d 3 u/dx3 from the transonic solutioa for the "corrected" to make u - W for the axisymmetric The correction is abvut 16 perccnt made
axial distribution is
This correction is
14
ln
n.
"
..
AE :)C-TR -78-63
so that v,'alues of du/dx and d udx solution I1 for later application. 'much
for low values of R, because the series expression for u is 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
Y+1
h.-y+ I+ M 2)qy
++
(29)
'
or
)1+a W-(Y+' I W2--
(30)
can
First, second, and third derivatives of W -itM with respect to r/r be obtained as described in Ref. measured from the source. by the same factor, r,, can be transormed by Eqs. 12. Along the axis x
-
r when x is
Inasmuch as all, coordinates ma.st be normalized (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 be y/r and d& - dr/cos v, evaluated dimply. Inasmuch as sin
H
but
.
ik9!
S21
tan
= -tan (M
-I
"15
AEDoC-TR-78-63
'for a
1,
dr (24M
2
2(1 +2 E:
2
M)
Mi
Thus
,tnan dr s r
=
(M -1) 2 2(1+-Z2 MN )
2
dM M
do
(I +Y-1
2
2 ) M
M
M
21
Similarly, for the left-
(31)
ad
2
(32) (33)
(34)
-G
and
G - OF, = (a + 1)
and, from the design values n and MB (and/or 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
[(Y + 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 degree
polynomial
(35)
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
2 M = 1) + D2 X + D 3 X + 0) 4 X + D 5 X4 + D6X"
(37)
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 D 6 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 constant 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
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. of constant values of the ratio *E/n.
h + 4,
nozzles lie in the range covered in this figure, and inasmuch as *F/n
values of MF can also be obtained. 2,8
1,4
2.4 -2, 22
"-2
1,8 -
12
.11.6
1,2 4
20
24
Figure 5.
Relationships obtained from cubic distribution of velocity from tonic point to point E for axisymmetric nozzle.
- -
- -
-'
j18
,. -
AEDC-TR-78-63
In
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. relationship is shown in Fig. 6. for a Mach number greater than about 3.4, 20
+
This
~18
16
14 ,, 12 10 Axisymmetric
Planar
0
0 Figure 6. 2
12
14
point B is
Using such a
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 CD along 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 the region ABCD starting at point B and progressing downstream. The equations for a right-running characteristic were given previously.
dy/dx
+
=
(19) (23)
where
de - dx/cos(
-It)
(24) (40)
'!2s_0nL,uiJ y
d4
(41)
where Aldo
:, Also
p)
(42)
2 (1 + 1.i M )
dM M
acot I
dW
(43)
[I
20
__
A__
AEDC-TR-78-63
Values of x, y,
*,and M are
right-running characteristic, C, and at the general point 2 on the leftrunning 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) (45)
(0 2 + A2
where A4 ( -c and
Y3 Y2
3 222
-3
" 2) sec 1
tan~ 0
tan (0
+ jA 3 ) + 1.tan
(46)
011
+ (3-
01)
4,
(Min0mi&AMn
+. sinl 6minp
(47
(48)
where
(3-
x1) sec. a
and
=tan X -X
Itan ((3
13
+1L
2
tan(~
(49)
+ S6
+ P2 +
)(50)
03
(01
+ 0,+
02
-(51)
p2
0 because a
-0
can be solved directly, 14 is obtained from *j by t~he inverse application of Eq. (25), and vNsin (1/M ). In axisymnietric flow, the equations must be solved by iteration. A useful first approximation for P 1 ) /2. ~3 )/2 and P2 and P2 is the radial flow values, P 1 1 -(p-
21
--
.. ..
AEDC.TR-7S-63
At all points except on the axis in axisymmetric flow, Eqs. (44) On the axis, the and (47) are defined because Y2 and y, are nonzero. /Y2 and sin 0,/y, are indeterminate with the form zero/zero. terms sin 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 * W 3. Equation (41) Cot g L and Eq. 23 can be written
cot A dW= W
can be written do +
'in 6 sin I dx
(52)
y cos ( + -)
-d9
(53)
as and
964
9S-Psino,
tan
= i__i..-3 2
=.
1
3
3
tan Y3
sin 0.
tan P
(X 3 - X2)
Y3
--
x2)/y 3
(55)
Similarly
Wa
- W S(Wl 3)
=3
x3)/y
(56)
-*
(57)
Adding Eqs.
and rearranging,
lia
y4O Y
=i
2
~ __ w
W dx
(58)
and
(M1 -)Y
(d'
(59)
22
,.
- .
,,
. ,
./
'. ,
2 2
AEDC-TR-78-63
(44)
when point 2 is
HilIl
and
st
410
Yl
(60)
when point I is
In
the region TIEG, point E becomes point 1 and the first axis point upsteam of point E becomes point 2. The complete left-running characteristic 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
Y
0.6 0.4 "-
0.2
01
0.2
0.4
0.6
0.8
1,0
x
1.2
1,4
1,6
1.8
2.0
Figure 7.
-.
'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 correction will be made for one design condition of stagnation pressure and temperature although it is theoretically possible to reshape a flexibleplate 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 circumference at each station equal to the periphery of the actual rectangular cross section.
The method of calculating the boundary-layer growth is based on
axisymmetric flow.
0I
dx M 2 -M 2 M El + (y.-_)
+H 2
M /2
dM + I
dx
rw
dr
dx
C
2
(61)
The term
just described.
24
,.k
AEDC.TR -78-63
as a first approximation.
with new values of rw and drw/dx - tan 0w obtained each time by adding vectorially the displacement thickness to the inviscid contour. The value of mcmentum thickness used in Eq.
PI
(61)
q
is defined by
-__d__(62)
0___
where z is measured normal to the wall. Also
P'I
'(
1--
)d
(2
iiO.
5 (,_ 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 values w a* and 0, obtained from mass-defect and momentum-defect considerations by
12
(64)
and
-- . 62
(65)
(64) may
Because rw-6 coo Cw + y, where y is the inviscid radius, Eq. w a~ be rearranged to give w
4--
w I
Y 2
SC 2
y sec 9W
(66)
*a sec
radius in order that no correction be made to the longitudinal location. The integrations of Eqs. using Gauss' velocity distribution
q/q,
25
.-.
.......
'..
AEDC-TR-78-63
ai
(68)
Tw] (q/q.) 2
(Taw
Tw)
(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 1' data obtained in hypersonic wind tunnels with which agrees better wit 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 Fc, introduced by Spalding and Chi, Ref. 18, Fc Cf =Crf
and Cf is related to an incompressible Reynolds number, R,
(70)
which is
fI
by a factor F, RO
c
(71)
19, is given by
(72)
'
L I/(p/pe)
which uses Eqs, (68) and (69),
was implied, but Eq. (72)
d (q/qe)]
26
-. Ni
Zi
A E DC-TR-78.63
upon which R
is
based is
The equation used herein for incompressible skin-friction coefficient is that of Ref,
If
20,
(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,
, (fI
(0,242)2 (log 110+ 1.1696) (log lie + 0.3010) Also as shown in Ref. 20, Eq. (75) agrees
(76)
21,
based on Coles'
S~L
the wake,
(2/C 1 ) 2 fn 15 + 0.5 Fn (C /2) + YC + 211
-
if
400 to a maximum
of 0.5885 at Re 107.
and then decreases to about 0.49 at R. to agree with Eq. (77), I1must continually The data shown in Fig. as shown in Fig. 8.
(76)
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 function 2 by a sine distribution, and integrating to obtain
____
(78)
j27
A ED C-TH -78-63
and _8*
C=i(~,.
(Y9)
2K
,'
"0.6
' 0 0
0.50
0.3
0.2
I~ l Ill I l J lll ] I I
ll
I I I f
ll
I I J I ll
103
104
105
106
107
Figure 8.
number based on the actual boundary thickness, not corrected by FR evaluated through the use of the kinematic momentum thickneis
q Ifrom which
Ok/6
dz
3N + 2)
(80)
(81)
N/(N 2
or
-a+
- 6\
+ 1I
(82)
28
A EDC-TR -78-63
0. 006
I
I Tl l r
1 TrTiJ
lfl
! I
I 11 I 1 [1 1
0,003 5
Ct,
O.004
0.003
0.002
Eq .(75)
0.L0
0 1 1 11 tl 111 ilj ,J
Eq. (7610
I J1 1 iIlll f ,I I I I ILill'
104 Re,
106
Ia7
Two options contained in the program subroutine for the boundary layer utilize Coles' law of corresponding stations (Ref. 23), (83) if Cf /Cf
i
. F
c
is
(72)
for a - 0 or a
1, then one 1
~ ~(84) ~/I~
29
AEDC,.TR-78-63
and
=
!S/!L.(86)
1?.2(Cf/2)2 a(Ta
T)
305(Cf,/2) [a(Taw
Tw) + Tw
T.] (87)
10
.1
N 6
104 R6
106
Figure 10.
Variation of velocity profile exponent with Reynolds number based on boundary-layer thickness.
(88)
where
CTWAW
(89)
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 temperature, T WT, and a nozzle-exit wall temperature, T WD, both of which are input to the program. of wall temperature,
I(A i/A*)'I A/A*
(90)
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 be calculated by the program, T
'I; T=
' T + 00 l)-1 I('l' )
(91)
h. + 0
where h a is the airside heat-transfer coefficient at the throat as calculated by Reynolds analogy from the throat skin-friction coefficient Sp
C p 2/ c2
(92)
i) 777,64885,
(93)
material, the cooling water, and the geometry and would be a constant if the properties were constant. The assumption is made that the bulk temperature 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 temperature, T
S31
AEDC.TH ,78.63
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 with Sutherland's equation used for viscosity. With the inputs of Tw_ and
.. :
TWD, the wall temperatures can also be calculated as functions of Mach number, although TwTmay need to be obtained by interation if the option 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 as used herein, the viscosity variation with temperature is assumed to be linear. The integration of Eq. (61) is started at the throat where it is Iteration is
involved at each point because Cf is a function of Reynolds number based upon 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 iterations converge within specified tolerances, the value of P is calculated a from 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 integrated numerically to obtain the increment in 0 to be added to a previously determined value of 8. The trapezoidal rule is used to determine the second point, the parabolic rule for the third point, and cubic integration for the fourth and subsequent pointa.
For convenience, Eq. (61) may be written 8' + 8P - Q. The general
=On3 + Gu_3 1
+ Gn-2
n-2 +I Gn
1+
Gn 0.
(94)
32
+ ++
"+"
..... ,
.................
..........
AF DC.TR-78-63
where th' C's are functiw~s of the spacings s, t, and u between the points and are given in Appendix B. Except tor e and ef, the other fn n values in Eq. (94) are known from previous calculations. Inasmuch as
0
Eq. (92)
P 0
(95)
(95)
On-. 3 + "'
+ "'n-1 0,1-1 +
p, lon
(97)
which becomes the Gn 3 for the next point to be calculated. The values of 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 the a initial values obtained by the trapezoidal or parabolic integration. The success of the above type of integration depends upon the spacing 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 characteristics solution are selected with care. After the values of 6* sec w are calculated, the values of a w d(6* sec 0 )/dx are obtained by parabolic differentiation and added to a the inviscid values of dy/dx to obtain dr /dx. This procedure is believed w to be more accurate than differentiating the value (* sac w + y) .a w because dy/dx is obtained directly from the characteristics solution and not 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.
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, 4 ,d to and 16 subroutines arranged so that the program can be overla conserve computer storage. The four overlays consist of AXIAL, CONIC, SORCE, and TORIC; PERFC; BOUND and HEAT; SPLIND and XYZ. The input data 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 coefficients to obtain ordinates at evenly spaced points along the axis or, in
the case of the planar nozzle, at discrete points along the surface of the flexible plate at which the supporting jacks are located.
Subroutine AXIAL. In this subroutine, cards are read in with the
parameters used to define the axial distributions of velocity and/or Mach ,Lumber and with integers which define the number and spacing of the
points on the axis and on the key characteristics and the sequence of
34
AEDC-TR-78-63
subsequent calculations.
If
the upstream end of the upstream distribution starts at the intersection of the throat characteristic and the axis. at the point where M to not use the throat characteristic and thereby start the distribution 1. This option would normally be used for a e.g. a planar nozzle, or nozzle with a large throat radius of curvature,
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 distribution 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 The stagnation
conditions are input, as are the parameters to describe the wall temperature distribution, the temperature distribution in the boundary layer, and the factors relating the compressible skin-friction coefficients to incompressible values. Subroutine CONIC. (29). Function CUBIC. This subroutine is used to obtain the smallest positive This subroutine is used within AXIAL to give the
1
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
35
.-
AEDC-TR-78-63
Subroutine HEMT.
This subroutine is
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, This subroutine is used within PERFC to calculate the properties at the intersection of a left- and a right-running charSubroutine OFELD. 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 characteristic. If desired, certain designated intermediate characteristics may
be printed out.
associated with the smoothing are read and the smoothing applied. Inasmuch as the wall angle is interpolated from mass-flow considerations, 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 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 applicable to these shapes for the number of terms retained in the series,
U,,
36
AEDC-TR-78-63
Eqs.
(2)
and (3).
is applied to obtain the inviscid coordinates in inches, and the abscissas 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 Three points at a slopes and abscissas to obtain second derivatives. time are used to establish the parabola, and the slope is the center point. but with less accuracy. Subroutine SORCE. from Eq. (30). This subroutine computes the coefficients of cubic the unevenly spaced points obtained from the charThe initial and final slopes are used together This subroutine is used within AXIAL to give the obtained at The slopes at the first and last point are also obtained,
acteristics solution.
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 This would produce such a gradient from the transonic theory used. 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 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. one plus a submultiple of 240. Subroutine TWIXT. This subroutine is used in PERFC and BOUND to interA fourIn PERFC, the flow properties along this key characteristic are used at the number of points specified as
polate the ordinate and other properties at a specified point. the specified point. Subroutine XYZ.
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 The points may be the same points as those a comparison is desired between the derivatives so arbitrary uneven intervals. input to SPLIND if
7.0
The design of a Mach 4 axisymmetric nozzle is selected to illustrate 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. degree polynomial. Leaving SFOA blank ou Card 2 specifies that the upstream axial velocity distribution is not a fifthLeaving FMACH blank cn Card 3 specifies that the The computed value of FMACH is value of FMACH will be computed by the program, and leaving IX blank on Card 4 specifies a cubic distribution. 3.0821543, which is greater than the value of BMACH specified on Card 3;
t38
A ED C-TF-78.63
therefore,
value of PP means that the inflection point will be 60 in. an arbitrary point. Leaving XC blank specifies
downstream of
the downstream axial and the positive value The values of MT,
and LR determine the number of points on the key chardivided into an even number of the contour points to be
acteristics and are all odd numbers because each includes both end points of each distribution which is increments.
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 (201.3 191.3). increment,
running characteristic will be printed for the upstream contour together with the right-running characteristic through Point E. The smoothing
AN 1716s563
ZO 1.
RO 0,896
VISC 2.26968E-8
VISM 198,7a
SFOA
XAL 10000
60 RC Ix IN 10
FMACH IG MO 41
3. BMACH NO 49 NF -61
4. CNC MP MG
-12.25 SF dA 1 Jx
60. PP JC 10 IT
MC LR -21 NX 13
CARD 9 NOUP NPCT NODO S,) 5S 85 CARO 6 PPQ 200. CARD 7 XSt 1000. TO 16J8* XLOW 46. TwT 900. XENO 172. IWAT 540. XINC 2. GFUN .36 ALPH IHT IR I1 1 CN LV 5
Hi
XMID
XINC2
39
-.
A EDC-TR-78-63
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 (90). ILeavi.ng IR 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. 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 Xcoordinates 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 inflection angle, the curvature ratio, and the Mach number at the point B. and 3.0821543 (computed), respectively, The selected values of 8.67 deg, 6, length of about 14 in. length of about 31 in. and nearly 120 in. of the nozzle.
are not necessarily optimum but result in a nozzle with an upstream from the throat to the inflection point, a from the inflection point to point 3 (see Fig. 3),
could be reduced to some degree although experience with Mach 4 axisymmetric nozzles is very limited.
40
.:,.40
"AEDC-TR-78-63
The number of points used on the key characteristics should be consistent 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 aerodynamic design of planar and axisymmetric supersonic wind tunnel nozzles., 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 distributions 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 characteristics to compute the inviscid nozzle boundary. The final wall contour is then obtained by adding to the inviscid coordinates a boundarylayer 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 A listing of the computer program and an input descripoutput data. tion are included.
REFERENCES
1. Prandtl, L., and Busemann, A. "Nahrungsverfahren zur zeichnerischen
Stodola Festschrift.
41
r
AEDC-TR-78-63
2. Foel2ch, K.
"A New Method of Designing Two Dimensional Laval Report NA-46-235-1, Downey, California, March 1946.
Nozzles for a Parallel and Uniform Jet." North American Aviation, Inc., 3. Riise, Harold N,
Variable Mach Number Supersonic Wind Tunnel Nozzle Design." AGARDograph 3, October 1954. 5. Sivells, J. C. "Analytic Determinatiton of Two-Dimensional Supersonic Nozzle Contours Having Continuous Curvature." AEDC-TR-56-11 6. (AD-88606), S. July 1956. Fluid Flow and Heat Transfer at "Design and Testing of a
Mach 4 Axially Symmetric Nozzle for Rarefied Gas Flows." July 1952, University of California, California. "The Aerodynamic Institute of Engineering Research, Berkeley, 7. Beckwith, I. E., Ridyard, H. W., and Cromer, N.
Design of High Mach Number Nozzles Utilizing Axisymmetric Flow with Application to a Nozzle of Square Test Section." NACA TN 2711, 8. June 1952.
Cronvich, L. L. "A Numerical-Graphical Method of Characteristics for Axially Symmetric Isentropic Flow." Journal of the Aeronautical. Sciences, Vol. 15, No. 3, March 1948, pp. 155-162.
9.
Foelach, K.
"The Analytical Design of an Axially Symmetric Journal of Vol. 16, No. 3, March 1949, pp.
Laval Nozzle for a Parallel and Uniform Jet." the Aeronautical Sciences, 161-166, 188.
"i' '
42 4
AEDC-TR -78-83
,*
10.
Yu,
Y. N.
* 11.
AGARDograph 35,
November 1958.
nmetric Nozzles Part I - Analysis and Coordinates for Test Section Mach Numbers of 8, Wright Air Development Center, Dayton, Ohio, October 1958. 12. Sivells, J. C. "Aerodynamic Design of Axisymmetric Hypersonic Journal of Spacecraft and Rockets, 1970, pp, 1292-1299. 11, Nov.
"Transonic Flow in Two-Dimensional and AxiallyThe Quarterly Journal of Mechanics 15, Pt. 4, November 1962,
N.
May, R. J.,
H. D.,
and Hoffman, J. D.
"Comparison AFAPL-TR-
16.
Edenfield, E. E. "Contoured Nozzle Design and Evaluation for Hotshot Wind Tunnels." AIAA Paper 68-369, San Francisco, California, April 1968.
17.
AEOC.TR 78.63
18.
Spalding,
D. B. and Chi,
S. W.
Boundary Layer on a Smooth Flat Plate With and Without Heat Transfer." Journal of Fluid Mechanics, pp. 117-143. Vol. 18, Part 1,
January 1964,
19.
Van Driest,
E. R.
20.
Sivells, J.
C.
Ludwieg Tube."
21.
Coles,
D. E.
Proceedings
AFOSR-IFP-Stanford 1968 Conference on Turbulent Boundary Layer Prediction. Vol. II, Edited by D. E. Coles and E. A. Hirst.
22.
Wieghardt, bei
K. and Tillmann, W. Zur Turbulenten Reibungsschicht Z.W.B., K.W.I., U&M6617, 1944, translated
Druckanstieg.
as "On the Turbulent Friction Layer for Rising Pressure." NACA-TM-1314, 23. Coles, D. E. Fluid." 1951. a Compressible September 1962.
AEDC-TH-78-63
Eq.
(2)
GT S
+x(i -aF
8S
+2X2 +,-/-(1-
+ GT +Gv''
2
-2 - V +
/
x,
.)- +
33
K +
...
2
+ 2++ 2S (2+ UxP2 Y
S
4 4y4 ... U2 U
2 + y2 Y
U63
3 +U 23 2y2 y6 . U 43 y4 +
_2__3
U0 PO
S y42
222
+ 2 \
(10 - 30)y\+
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 (782
- y)/270
US
2 + 3507 y + 7767)/272160
GT
av - 5 y/18
GK
-
(A-5)
33y + 9)/24
(2y 2
(A-6)
(A-7) (A-8)
U4 2 -
(y + 6)118
U 2 2 - y/9
45
,j
.1,,
AEDC-TR-78.63
U 23 "
Up 2
LU 0 = (26 y 2 + 51 y
For axisymmetric flow,
27)/144
GR - (15 - 10 y)/288
(A-14)
(A-15)
(A-16)
9)/1152
GV - (y + 0)/8
* GK - (4 y2 - 57 y + 27)/48
(A-17)
(A-18)
U4 2 U2 2
(2 y + 9)/24
(4 y + 3)/24
(A-19)
(A-20)
U6 3 - (556 y U4 3 -(388 U 23 U
+ 1737 y + 3069)/10368
y2 + 777 y + 153)/2304
P2 - (52
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'!
--lo o ... ") :- ig i:,, _ T3:I r"l~:IIrI ]1'
AEDC-TR-78-63
(3)
can be written as 4 2
NS
2(3 2-) )s +
6 v
42
y4 +V +v 3'
s2V 2 2
2
+V 0 2
+ V6 3
43
23
-V
03
x 6U 6
I
3
2y
+,
1 .5q
(2 X12x2
4 U
Up
PO+
'
-2
+ x3
"4
..
)/
(A-26)
V2 2 - (10 y + 15)/108
V0 2 V 63
-
(A-28)
(A-29)
V4 3 "
(2254 y2 + 6153 y +
-
V 2 3 - (5026 y2 + 7551 y
2
-03
4923)/77760
(7570 y
+ 3087 Y + 23157)/544320
(A-33)
47
44'i
i:
A EDC-TR.78.03
III
L,
,
v 4 2 " (Y + 3)/9
v, 22
(20 y + 27)/96
V 23
"0O3
48
I ,,
48i
A EDC-TR -78-63
F
as
a+
Fbs
b-+ Fc
cs
+dd Fds c +
(B-i)
Areab-c
(B-2) (B-3) -3
li
where
Aresc-d c- "F
+F Fdu d
tota
Gb
+ G
c + Gd d
(B-4)
S(s
+ 4t
(+-5)
",2
iij
2 2
SFbs
bs
+ Zs I + 4t + 2u
+--)(B-6)
12 t(t + u)
a 3 (s + 2t+(2B) T2tu(a + t)
cs
I::
's-
7)
Fds
1 (s + t + u)(t + u)u
s_(L+2t,) 2
(B-B)
at F
bt "
-
t 3t+2u
12s(s + t)(s + t + u)
2+ t
+
(-9)
2u
2 t
12s(t + u)
F ct
2-u(s + t)
49
_....__
_%.-_
__._
.,.
AEDC-TR-78-63
Fdt
F
--
(B-12)
+ t + U) 3)
+ t + U)
u3(2s + 2t +_y
rbu
-
3i
(B14
(B-14)
'
12st(t + u)
u (2s + 4t + u)
cu 2 12t(s + t)
(-5
du
R .. 2
(B-16)
-F
+ Fat + F
(B-17)
Gb
bs + Fbt + Fbu
(1-18)
Gc
(B-19)
(B-20)
then
(B-21)
Fds
-Fat
-dt F
-Fau
-h/24
(B-22)
F Co
-F-5h/24
(B-23)
bu
50
AEDC-TR-78-63
= F
19h/24
(B-24)
Fas
Fdu -9h/24
(B-25)
Fbt
Fct Gt
d
13h/24
(B-26)
3h/ 3h/8
(B-27)
G c
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
j'"" 2 2$WA,,
AEDC-TR-78 -3
Columns
ITLE JD
Lard 2 GAM
2-12 14-15
1-10
AR ZO
11-20 21-30
Gas constant, ft2/sec2 R. Compressibility factor for an axisymmetric nozzle, constant for entire contour. Or, for a planar nozzle, ZO is half the distance (in.) between the parallel walls, and the compressibility
factor is one. Turbulent boundary-layer recovery factor. Constant in viscosity law.
RO VISC
31-40 41-50
V1SM
51-60
If VISM is
lb-sec/ft 2
2"
lb-sec/ft2
If
; T <VISM.
-
SFOA
61-70
Used for nozzle with radial flow region if 5th-deg axial velocity distribution is desired. If positive, the distance, in inches, from the throat to Point A
52
'-
.~-
AEDC.TR-78-63
on the characteristic diagram. If negative, absolute value is distance from the throat to Point G. If Blank, 3rd- or 4thdeg distribution iu used depending on value of IX on Card 4. XBL 71-80 Station (in.) where interpolation is desired (e.g., the end of a truncated nozzle). If XBL1000., the spline fit subroutines are used to obtain values at increments evenly spaced in length. Card 3 ETAD 1-10 Inflection angle in degrees if radial flow region is desired. Two characteristic solutions are obtained, one upstream and one downstream of Point A. If ETAD = 60., the entire centerline velocity distribution is specified and only one solution is obtained and the inflection point must be interpolated. If ETAD - 60., IQ 1, IX - 0, 1 on Card 4. Ratio of throat radius of curvature to throat radius. Must be given if ETAD - 60. If FMACH is given, RC is or FMACH - 0. calculated. If LR - u, IX = 0 givws thirddeg equation betweev Mach 1 and EMACH, matching first and second derivations at each end. If LR 0 0, the value of RC found for LR = 0 is used with given value of FMACH to define a fourth-deg equation. If IX = 1 and FMACH is given, RC is calculated to define a fourth-deg equation. If LR 0, a new value of FMACH is found, compatible with the value of RC calculated for LR = 0. Nega-. Mach number at Point F if ETAD # 60. tive value specifies Prandtl-Meyee angle at Point F as IFMACHI *ETAD (usually around -7). If PMACH and RC are given, IX = 0 If and 4th-deg distribution is used. FMACH - 0 and IX = 0, a 3rd-deg distribuIf FMACH = 0, and IX = i, tion is used.
a 4th-deg distribution is used. FMACH is
RC
11-20
FMACH
21-30
If
ETAD = 60.,
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 Point 2 2 C. If ETAD 0 60, positive CMC giyes d M/dx 2 -0, and negative CMC gives d M/dx' 0 0. If ETAD - 60., CMC ts positive. Scale factor by which nondimension coordinates are multiplied to give dimensions If SF - 0, nozzle will have an in inches. inviscid throat radius (or half-height) of If negative, nozzle will have an I in. inviscid exit radius (or half-height) of
SF
51-60
ISF1 in.
PP 61-70 Station (in.) at Point A. PP - 0 gives coordinates relative to geometric throat. Negative PP gives coordinates relative to
Nondimensional distance from source to 1. 1 requires centerline Mach Point C. XC No. distribution from Point B to Point C to be read in as input data on Unit 9. Otherwise, positive XC gives 5th.-dog distribution if CMC positive and 4th-deg if CMC XC - 0 gives 4th-deg dietribution negative. if CMC positive and 3rd-deg if CMC negative. Negative XC and IN gives 3rd-deg distribution with d 2W/dx 2 not matching source flow at Point B. If ETAD - 60. and XC > 1, XC is ratio
of length, from throat to Point C, to throat Negative XC gives 3rd-deg distribuheight. tion in M; XC - 0 gives 4th-deg distribution; XC - 1. XC > 1 gives 5th-deg distribution. requires centerline Mach No. distribution to be read in as input data on Unit 9.
Card4 MT 1-5 Number of points on charact',ristis EG if Maximum ETAD ' 60. or CD if ETAD - 60. A zero or Use odd nuwjer. value about 125. negative value stops calculation after tonterline distribution is calculated if NT positive.
54
NT
6-10
Number of points on axis IE. Maximum value is 149-LR. Use odd number. A zero or negative value stops calculation before centerline distribution is calculated but after parameters and coefficients of distribution are calculated.
IX
11-15
Determines if third derivative of velocity distribution is matched. IX - 1 matches third derivative with transonic solution. IX w -1 matches third derivative with source flow value. IX - 0 does not match third derivative but gives constant third derivative if RC - 0 or FMACH - 0. Determines type of distribution from Point B to Point C, positive for Mach No. distribution, negative for velocity distribution. IN - 0 for throat only. If XC is greater than 1., the downstream value of the second derivative at Point B is 0,1* JINJ times the radial flow value. Similarly, If ETAD - 60., the second derivative at Point I is 0.1i*IN times the "transonic value. Zero for a complete contour if ETAD 0 60., 1 for throat only or if ETAD - 60., -1 for downstream only. Number of points on characteris'ic AB. Maximum value about 125. Use odd number. A zero or negative value stops calculation similarly to MW. Number of points on axis BC. M~aximum value is 150. A zero or negative value acts like NT. Absolute value is number of points on characteristic CD for ETAD , 60. Maximum value is 149
or 200 - ND - MP
-
IN
16-20
IQ
21-25
MD
26-30
ND NY
31-35 36-40
55
*.
. .............
|
. . . . . .
--
AEDC-TA-78-63 A
MQ
46-50
Number of points downstream of Point D if parallel inviscid contour desired. A negative value can be used to eliminate the inviscid printout. Positive number if boundary-layer calculation ir desired before spline fit. Negative number transfers control of program to JX. Absolute values greater than one are used to approximately halve the number of points on the upstream contour even though LR + NT - I points are calculated from charnetwork if LR > 2, or (NT + 1)
points if LR a 0.
rB
51-55
Iacteristic *
JX
56-60
Positive number calls for calculation of streamlines, zero calls for repeat of inviscid calculations requiring new cards 3 and 4, or, if XBL 1000., for spline fit after inviscid calcu1 lation, negative number calls for repeat of calculations requiring new cards 1, 2, 3, and 4. If not zero, calls for printout of intermediate characteristics within upstream contour if JC is positive and downstream contour if JC is negative. Characteristics are (NT - 1)/JC or
(ND - 1)/(-JC). Opposite running characteristic
JC
61-65
through Point E (or B) is also printed. IT 66-70 Number of points at which spline fit is desired if points are not evenly spaced, such as Jack locations for a flexible plate. Used only for a planar nozzle, inasmuch as a nonzero value calculates distance along curved plate surface. Positive value of IT tequires additional cards to be read in (8 points per card) after boundary layer is calculated. Absolute value is number of points on throat characteristic used in characteristics solution. Negative values give printout of transonic solution. LR - 0 gives M - 1 at Point I. Number from 10 to 20 determines spacing of pcints on axis for upstream contour. NX w 10 gives linear spacing. NX > 10 gives closer spacing of points at upstream end than at downstream end. 20. Ratio of downstream NX - 0 same as NX
LR
71-75
NX
76-80
56
A E DC-TR-78-63
increment to upstream increment is (NT - 1 )NX/10 (NT - 2) . Optimum values, usually 13 to 15, determined by trial and error for specific contour desired. Negative NX used with negative LR limits printout to transonic solution.
-NX/1O
NOTE:
A zero value of MT, NT, MD, or ND will allow a repeat of calculations for parameters specified by new cards Nos. 3 and 4. A negative value will allow a repeat of calculations for new cards Nos. 1, 2, 3, and 4.
Card 5 NOUP NPCT NODO Card 5 orinviscd 1-5 If smoothing is desired, negative NF. times upstream contour is smoothed. Smoothing factor in percent.
- NPCT/100.
Number of
6-10 11-15
Smoothing factor
Number of times doimstream contour is smoothed. If boundary-layer calculation is desired using points smoothing). calculated from characteristics
or
solution.
(No
'
Card 6 or Card 7
If boundary-layer calculation is desired using evenly spaced points interpolated from spline fit of points from characteristics solution. If boundary-layer calculation is desired using evenly spaced points interpolated from spline fit of smoothed points. 1-10 11-20 21-30 Stagnation pressure (psia). Stagnation temperature, Rankine.
PPQ TO TWT
Throat wall temperature, Rankine, if QFUN - 0. If TWT - 0, the wall temperature is assumed to be the adiabatic value. Wall temperature, Rankine, at Point D. For water-cooled wall, the bulk water temperature i.s assumed to be 150 lower than specified The cooled wall temperature distribuTWAT. tiou is assumed to be
TWAT
31-40
57
.. t=,..
......
!.
j.
"u' ..........
I
AEDC-TR-78-63
(TWT TWTWAT +
WT)
"' C A*
1)
where A/A* is the area ratio corresponding to local value of Mach number and Ac refevs to Point C. For negative IHT
TW - TWAT + ,TWT-IA x A A*
QFUN
41-50
is obtained by Reynolds analogy from the skinIf QFUth is specified, friction coefficient. input value of TWT is ignored and TWT is calculated from QFUN. ALPH 51-60 Parameter specifying temperature distribution in boundary layer. ALPH - I. uses quadratic distribution both in the calculazion of the reference temperature TP and the calculation of ALPH w 0 usas boundary-layer shape parameters. parabolic distribution in boti. calculatiovs.
ALPH -u-1. uses quadratic distribution for TP
and parabolic in the calculation of boundaryWithin bcundary layer, layer shape parameters.
T - Tw + a(Taw+ TW) Tw) (U/11) - Twj (U/Ue) 2
Tes-
(Taw-
where a
IHT
61-65
Integer which determines temperature distribuIf nonzero, IRiT determines tion (see TWAT). how often subroutine HFAT is called. An absolute value of IHT greater than KO, the number of points "on the upstream contour, will prevent HEAT from being called but will allow the choice of temperature distribution to be made.
58
111.01 11 p1 Ml
AE DC.TR -78-63
NOTE:
HEAT is a special purpose subroutine for determ.ning heat-transfer values for the upstream contour. The subroutine HEAT incorporated in this program is a dummy. 66-70 Integer, parameter specifying transformation from incompressible to compressible values. If YR w 2, Coles' transformation is used for Cf and Re If IR - 1, TP is calculated by a modification of the Spalding-Chi (Van Driest) method. If IR - 0, the Van Driest value of Re 6 is used, but if IR- -1, Colas' lav: of
RD*Ree,
axis). If ID - 0, these effects are omitted. Negative ID suppresses the printout of the boundary-layev calculations. LV 76-80 Integer. Absolute value, usually 5, deterJ
mines number of rimes boundary-layer solution
is iLerated so that radxus terms in momentum equation refer to viscid radius instead of inviscid radius. Value of 0 or absolute val-e of 1 uses inviscid radius, Positive LV repeats boundary-layer calculations for new set of parameters on a new card if XBL j 1000. Card 5 METAD 1-10 If streariflines are desired, JX positive, (No
smnoo thing.)
Inflection angle in degrees for streamline 60. on Card 3, use ETAD - 60 on this card.
desired if ETAD 0 60. for Card 3. If ETkD -
11-20
Fract ion of r'.ottour desir.-ed if ETAD - 60. Otherwise, QON - lTAD (in Card 5 divided by ETAD oa Card 3. Value to updite JX for subr.eq4oent calculation, JX XJ.
X,1
21-30
!*.;.
,,
2
'
a a'11d598Y
I,
AZ-DC.TR-78.63
Card 5
(JX zero or negative and JB zero or nega(No smoothing.) tive). If SPLIND usad after viscid contouir (,7B positive and LV zero or negative). No v'.oothing of inviscid contour. Or, if inviscid contour is smoothed before SPLIND is used. If inviscid contour is smoothed, boundary layer is added and SPLINE Ls desired. Station (in,) for throat value of X. If XST a 1000., program uses value previously determined by specifying PP on Card 3. OtherWise, value of XST is used to shift contour
points by desired increments for arbitrmry Station 0. XLOW XEND 11-20 21-30 Starting value for interpolation. of interpolated X a XLOW + XINC. Second value
End value for interpolation. If zeto, SPLIND is used to calculate slcpe and d 2 y/dx2 at same points as previously defined. Increment in X for interpolation. If zero, and BJ > 10, contour is divided into BJ increments. Value to update JB for subsequent calculation.
JB - BJ. If negative and XEND - 0, interpola-
XINC BJ
31-40 41-50
tion is made at diecreate points read in on subsequent cards similar to case when IT > 0. XMID 51-60 Intermediate value for interpolaticn. Distance (XMID-XLOW) is divided into increments defined by XINC, and distance (XEND-XMID) is divided into increments defined by XINC2. Inctements in X between XMID and XEND if different than XINC.
XINC2 CN
61-70 /1-80
Number of copies desired of final tabulation of coordinates if more than one copy is desired.
,j
I,..60
...........................................
,
AEDC-TR-78-63
lz xxzsxxzzxzslzzxzzzz.0p-UeoxxzzRdgzzz*A zx
;1~M~ W2eV% n n n n t5n n n n n n n n n nnV
4;
a,~
Iris.
0 1
O-cc
@0.f
9
i
. 6
S
M o.
IA
0
as
0.4'W..
4MR3LO
*1
EM
~~~~~ @@ c
WhiP). UhZ 4*1; :af
~
9X
4~
ts
.0
a.
wi
iE
C 4 .0%.J Uu 0uJ%~cV
03 SeI "I~ m
--
0 ru4-.J
*0
4T%
V"OO
vv vhuIwwu
~I
0*
04
.. 0
).%.kU.61I
AEOC-TR-78-63
zzzzzzzzxgzzzz zzzxxzugzzzzzzzgxxzzZZZZXzzzzzzzzxxzzz;
941-
4;
a-9!! *
%..O. -O - z J
. 0J j A
a0
*.d
---
xw-
moo1
WIMa
IIAC
..
po
34
tv
5-
k"I
*4
*Jj
hIg.*
624*1
c.A Trrzt
77-
AEDC-TR-78-63
dN*Jfli
-AIzz NNNNN
0zzmaff
zPzEE4E4EEE
zzx
44444444
Z11
C>.0 >.g
I-
-0
,o
I.a*
1 ;.
&A 0. v m 0
3.pc 0
jz >. IL IL z 0
0 >. 0 4L
ccza.'
No *0 It
-0
1. N 3 .I0q
a Ax
45
w>
2e -2
.,C _j
0K
-
x
3r
IL
CL
0
4LJ, O
_;
z-
0.
u 400
x wV %
0i;2 3 aJ .O0204
40 W42 U 0 0 0
3
01 ).
co
0*O 12 Lo0*
~
4.
.
gag;;
N~ -
..
*.
JI
4 3
o
.
2jz -
0 0
3
3. *
0
'
*~ P-
hi
LN
N .2
~ ~
*
0 0
**04
.4.-.0
..
0.0
44
0~0
) 3Ix -
~ v0
*
0 1
U . 4-
00
1 0 3N .-
3
3
*02 03JI
N
0 *
C..
x -0--t
co-
4.P.
32
*-
3-
N *3-4
***
3-4**P
P.3-
3-0
AEDC-TR-78-63
4;
f4
Wo
ICc
iNa
w N
owh
Ni0Jcn;
44
"o
hi 04~~ of 0 00Im 1! WC .
OJ C;
VI x9C 2oW
mu
3-
--
VVo
---
wm4oo *
u K
o 4a*o= N N KyN
.01 M00
* 64
A ED C-T R-78-63
S4~S P. @
..
N~
~
4
0.
9:!
~
4
44
~
44
@0
NM
44
04
4O
M~~
iS
E..4
~ 4I,
44
*I
~~
44
,.j
1U
44
44 4
44
44
44
44
44
4 z
w W UO!x= *v ow20
P-O-6 tU
VS UO 24.4.. CLt
S. >
Z* . 4 4 S**o% 40.4 Z 0 h
0
:-W 0
4.
% %
a 0 A4
;T
3000. 0
,z4 I
4 C4 . O
~S
Ai
.0
* 0
~
U.
VW
2 -
N -
0-
.00 on
-sO
P4% %t
w
C h4 000.h. N. N0 :00a44*I,.
-
. S *IlSUzNl-'. X.
".
. -5
hNO
WO
O.P..K
.1
4nn
~ % ~
N-0* -
f~ 4 g..T ;,
9** 0 0 000*
4 US -
2 *00 e N* Z
5;
U;0.Z0 OO MO S2.
US4KS00t.JS-
-VS.4
4.P~S. K.
5.
SOK
*K65*
AE DC-TR-78-63
I-.A.90.Ncome*00.
0:- 0- 0-
- - -e-
E3 ---
~400y
Im
0.e
000
x 4 x -i z . S v
cc
F" -0fl4 e i.n SA 0 f0Ad 00 *
4
*;,"I * .4
.P
2 u .
.. 4... * 4 -0 *S @ .
-* 10 **N
--
* 0 0 4 % I twoM ID.5S0 .0%.*US0. S4 0-4INCI 0.0%401 OPA* 0-0* ZO *UW -00A *0 M0t1-. 4 ** 54 00 fu N*4N *1rP *~1 *0
.
t If * Uu 01 1 U A0 & As
Do-
0-
c4 C! w
* 554* 0
~
-. zo::4A U 06 PU N x S U 0 o w~~X a uv u 2m 4 o a 4 & 2u v*@ J0 K ZWJ N W.U V iUCI Z"--WO * J4 -2.& S *D vu...v N
u. mg4 -= u x. 0. M S
00 4 me... 9.0 30 0
vuIC
C 1
NO 4
5-00vv
*PS.
040 -
-C
UUS664
A ED C-T R-78-63
00
z
z
-0
wu =--
o1
4 2
ur
31 hi.-0
-a b
-4t
a.o~04
- Il0 2IE 0 i
;
a.
- hi z -s 1-w 4 P- 4 0 U w W ;;Va 0
We%
Z
zza
0411
. 0
21
0.-nhi 0 00
a. 'a0 a * uo w
rn-
wo
rn-
-O
U INOM1XE 3
o-z:-
;;. 0-za
*-;;ihv - z hi0
KLU-tU
I-z w h0-hi
0.
--0.LI
*IO
z * mz 4U-
~U
hi-ww
a.!M
0U-UX
Oh
hi
411
-0~
4-
4-
1-0
-0
VI
67i
a.3
..
AEDC-TR-78-63
II
aa I.In 4I W
o&0.0 S.-
a 3l ':0. 4
ly
U00
IL.
0o-
l,
-
a
Z
w -2; w.. m . 0
a -
0.I
3 XIL3
C-. od .1 0
m I-t4L
a5 I-n I-
I CL
IL
-!;
I0
,' I~
;10. . A0 '
~~~~~~
00
A1
V-03
i)-
468
AEDC.TR-76-63
........ 4t
4*
4b
4 b.4.
1414M44-H*~
-----H
""H
M"1 4
!Si
i
1(344
~ ~ ~ ~ 4 4444
t
444444q44q4.4
O'".J.=
' =. .i =.
,z:sc.
%C *
.=a =..l'
? " "
!
.
. '
.
' aI .
! fo
.
*,, 5..
.. ... .....
.. ..
=. ...= .= * ill,
-"''..
UU,- 3
0 -
68
~~68
..
A E DC-TR--78-63
WN
06
'Nm
Nr
x u c' Dt C ft 7
Z-4 CW ft
to
5;
0 ft 0s 9ft -4
q i'.x NooftS-1D;
z IX 0h
*.
2 1" w 0.uj K1 w4
xt
22 0 WT
0-oK ;;
*go. a
22l.z Wz W 2
w -4 4-: 05U Z0 WK
J
5-
-1
ftw -15
o 0z
u
5-
2 *.W
Li
c*; N1 c ly
--o N 0v
P.S;;
0 rY 0 fW 5-
;h 0
-
0c
-K
ftoa
US _ju
v) 0 -c ry 41I W4
wz.5 5h;h -9
611444 wLP W
--
NW U 55
4 g
--
y-
5U-
Si .0
i-
.--
--
ft
04W69.
AE DC-TR-78-63
-a~1 ZO
4444AWW
31 U. a.
.1
Wa.
0-1
3 I-
.j
lz
-. A.
jI
AA
* 06 a-
;I
3.. W 3W
a * a.
tW W -ma
--
a 0
1.1 t -WNO~ --
0 U
: Z
t 1-
tJ. I1
ti a A
-
W -W O--
3 xA *3. sa.*l0
a vQW OW 1.... 3
4-.
~~~~ 00
. W. W*.W MO O
ell 1*--0
0n0
-~-
Lb.0. W OW.
44
M1-4
.0
IA
IA4W
*0
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.
ZO
m*
I
I0 a.
I-
301
N 'a
4s
ILZtoIx a.
o
. -a
CL CI
-w LA f0o
o* 0 4 ;
w
0m 6
0 www
w 0 o 2
IWO
06
a.
".1
0S - - - o n
hI 20
a. w
a. !
a. iD
$0N4
kI.o. 6 I*"We3-
v Z"
4*
-"Ix.
h a a -a s-s-W .
-2X6v
,.-4 W 0s
4;
V/1
4 54
2 ~-oxwI.. &I,0
1 0z
.
. 41Iaa
4 14 OX
IN : I
a ;
WI UK 31 11
* x. a.
b uWs S wJW. 0 UI2*
'o .2I .0
66 -
. . wO hIS s. .- aNs
3-
o a.
l 0
-. hI ~~0h
12 hJ
le
N *
6N
4 -N
h..
S-J31
s-
4664
7106
-. a
AEDC-TFI-78-63
IL
a -it
-
T 2
.j
;a
-j
V; lI;
*
M 0* V~
-IWK
-::
&a
'a. W *U-hl . .W 3W&&UNM 3~~W * * V
Vto~
*8 0 00
4mw
.6
qh
09000-z
-
-aM W"z~ M
VmflOiWJ W00
V.I
A.
u XX in
-I
0@3 A A4V z*
*W 4
Vi4
0.**
-.
ha.. mO
a.U %~K
ww
b0.00.
0472
AEDC-TR-78-63
a 4443.
0 I*EI
.310
.0
1.1 L%uZ
aj31a
6a
I o 0PDO . Z
vu
.UI. .
u.
" '0 -9 0-a
0.M cN wU i-
2Z
02 0
I *Ifom 0 -
-0.2- -% 00%
a3 a a
00 4 U4 vv 00 --
R
*
0
a3 z --. x-
s.
m wa w wU *-------
o .~ .U 0o w .. m
uuav
Z
-----, I..0 -
OZU.VV
in-.x
In
inS 3a-uunInE2-
i-S
S.
AEDC-TR-7843
II
a,
fn
emL0M0
ZU . 06 IN
a. a
4 30 IL
NISI
a.ou~
00 U z 0z
0
w
ZR
a. 0 CL
4. CY UCC o aV * ODI 1W 0
CL
U" U PC CC
5---:VUUK:'-e c
vo VI- U U I.-S ; 4: -
.
U
U a 5
-~~ x
U 4D -a
--
1 -'IJ0 *OU2
--- v 2 A R
hi .
a-h~
.- K U
UU4
lay
Id
44II.
Kh.S-5A4
-00
*.
474
AEDC-TR-78-63
~w
*
di
~ ~~ ~
idid 4 44 444 id 44
~
id 4
~
ifud 4 id 4
~
idid 4 4
~
id 444
~ p---N---@**La5e0o*.cce ~ ~ ~ ~ ~ 0---------'
00 id d id id. 44 444 4 44
0000000.........
SSSSSd i id
44444
id
4 4 w K 44 W;Z
44 4 44
t 4 4 40444 4
4 44
4 4
4~( 4
0~~~.
llZ 'o
; e - %
4U
I:
N a
-V Z 0
0 di
.0 *ul CL
La S a,
-1
w u a:, 0-
U
:!0~
ga
0 a* --
0
000 U.
0 Z. ON Wi 00 o
0-
I0
-. 0
sK U 3*S *S
a3*
0
11
-9
0 Ui MIA* -
i.;~
;U !w
2
XZ
M%
-rI 4 1 CL0
K M0 M -ONO Ni
0--46.
Z 0'21
04
S --
*U~
z; 1 :&II
*n
0 K',
m~l'
.. -
LA-
l..W
*-*
0 0 K~l.
m.0 S .to
hjO- 9-0-0 N
'
40-l-
0.~
4.K'
i75
AECT-86
I,
Iu
xN fu 'a2W 02 4 M4--C fn 119 X44UXO.W 1.49-Z U W 119 4. Lr .
a.~
0. 0.
z
C 4u M4 Fn
13 0, U. 20 OZ -om
0. U M4 X .ux
~60.
7 211 .40 4. --MW
-z a
314 Z in 0. 11k
10
400
-4493t', Vol4.j
6:,Joao.N N r
f3
440WW .N
. 2 4
4rs
xIT
.03Or 0.0.
2.-Z XX~
0 32(
4 *.
~~
0. O
1-
W 434 - 2
ft a90 0 zX ,
~
. 06
'I90 0I.
.-
n .
M U.. * . 9 1
11
IZ
4430.3
W 0!2VI4k.4N-044VIr(0 0 6 - -41
394
ID3 *
0...
40 N
...
40-
qpU-.4p4.-.
-.
42
1.JZ J-k.
030 -44
2VW~uw4
76
A ED C-T R -78-63
N400Nm4A4N400-NM4IA4N
0,..l
4114*400 oU.'%M4t40.I1 -N
NN04
OOO 404
000 4tOO
0004
000----------
.0* -Z In 0 vC 3N x 2
v 0
~
i
E 4 0. ZZI 4 1
j
-
N
-o
N -
cc w. -W 4 .1 NI :o NI- X i 4 c x v 3 U
44
*4 M,
4.m
4!. i
_;wo 0
N9 Z CC
4 .6 Aj I%*
4r
2 4r
2~ P 4 IA
1-
A-
-C
W 1 .91 ju CL 2
6 am I.
Iz
o tn
-9
to0.1 4
2-C P:
mmt11m
*4
42
-9
06
0o
4.
4-4~~0 14
4211 r *Ly1 2 4
N 42Lm
. *
W~l U, 0Z
-j U10
X 410 .1
6611 It wa6
CL 4j,!*-
--w. 6
0
-4
N 4 *6 -w
3'.94 -W
40 ..3.
4x* .
W
*
.2 '.2
-1.. N
3 2 0 '1 4
N,
-U 34J -N
2-2
~
... *1L K0
2~
.X 0 ILN4 4. I* .
I 0
~.
z - 'i CU w 44w
%2 * i0 0U ;* tx 6i1
.T-. 0 2 ZIi 4.
00..3V3 OONNOVM-41311--
zN Z;i~~ a, co0
O*j
"*'
02;zz.0g2 20200n
Z 3 0
O*.I
-.
*0L
14NN.0% 4 4 ...N4
*41W m. .. 44 .
3 42On4 m4
.42 .
Xl
0 2D
000 Z -0.. 4200212U2242-211222222~~w.0210-22222 6-----U1* 0~.. 1 0 0 m 11Z 111 mm4 c. 1N1 11n IN1LiI 000 NM 1 0a 4 & 0o N 414
0Nw-.x4
-
Nmm'w"..Li
-OOOOOOIOCI . 11 11
10
.o02
1000.000002oo
77
AE DC-TR-78-63
It.
00 to ic4% % ;N -MOI
I b
~ W .
t
% It
0I
.2.
.4
. IL
. .N..O
..
fu
u4-000 %6OPJS0b
* .N 0. 1 'D in F 0 *n " m -
*%
j , 4
;A*..4 29.ia_
0
N %
T4
_ . ). a
~
:4N A
~
;=I-;-
of a %
00
-o
zt* *JAON
P.
: 3W 4 0; a I 0 0:
z
00. 00409400 ga
N a
0 @ P f" w . 3U 2 0
x.ouP!O
hi
aoDoi. 2.--
n AJ*
0--
---S 3-
'
4N
Oigg
5%
t~
hO
078
A ED C-T R-78-63
za6O o -----
OEA0-
- ------
N M*IA.0 0.NIII0.Q
wo
cccoccooooco000000
1 000
, It
P1
10 0
%-t
.0. a0 o L t
~ C.0w
t100 Na 't 0
00
IL
~I
W -c
I lI 0 . a
W a 0
2 o I 4a3 10,
-0% 1 0D*
A..
A NN 0~~l
. ig 0
I wN *
I I I N6.'.I.-P 4 I
I . 4
)
--
oS .II 0N U-O . CC -o
tb
A.
A*A . t~t A.
IAX
OSAI-
t-.
P79I
AE DC-TR-78-63
-----------------
---
--------
--
------------
$- a
0MN
C~m
0-U
0I 0
kM
CO Z -d I 0 00, US.. U .
S.%0
-1t1 .
h-h CO
fu
ow4
80
A EDC-T R-78-63
AJN
NNN
NNN
IJNN
NAJ
NNN
NNNIJNN
00000000000000000000000000000000000000000000000000000000
-,
0r u U. In 4* t 0ole 00 0.
0m3 0 C 0
4yq
v.
0 4 0
-
U 0
8 Nd
NU M0
8rljzX
0. P.I 0- K
09-L
N
*o
o U* 04.. 0,
0 % 0 1O
0 -ig 6
NSS 0L
) 2 0x
0onNI
--
49 - A~u A 1
0 8 -
a 0-
9z"0!zxx" knIh
h. U
I-
-2h.
LA
P-
981
AEDC-TR-78-63
CmNNA o
tuut4NNNNNNN N
NNNNJfy
0
00 61-
ar
mUa.mam
wee
wwww mw3
ILA
P1
Z3. N,
lyu
Iz 0 1101 314
ox -wI
IO
4 0
4 0a
1.. D*
02 1
I4 am
IX-
0 92 3 4 *0 0-C 00-0
UZ N6 h. h li4 at 4
A0 * * X, L. Is- 00 ,"
0
~W.e 20
~~10-.
U.
30-3
atU0Oa --
hi 4 41 Cm3 z S 'A hS -N214 K imk0*11 2--. e~ 07a4NIz , wPPP tov 1 2~3 '0a W3NW3.,. *S0 31 P-X3 4.P 20
*I.-
TOXOM 2
2 a.. a-
IA
82
AEDC-TA-78-63
FYMNNMO N cc
'"r o
O MM
MM
MM 0 M
M 0 MM
MM
4M MM 00000 M
1M3M
00 MM 0
0 MM 0
00 MM 0
MMMM 0 0 0 W m
~~~~~w wot
mff
clV0
4-
30
Oft W,:01
r.1.
itU
0W 9 4
0 , . z-
iz
0 0
91 x
w W3 P-0
be r4
..
Inh
-I-
1all
AEDC-TR-78-63
J UA
SO w *
X hi Z
xLii 01 a-w. n
-Z
or1
LW
-, co
J
I s
0I
K 5 .
ON or0 W
Q0
(0 UIm
-
0")w0Jca
--
-zz
0 zo Z 0
M4 00
Liar
C; : Z: a~ w-.
W 0zK
Li U
U
-L
.5 -
5 Li
0:. .
-3 1) -.-. 0
*-
- --
U)
C Li Ko U
xZN )
5'-. vu i P
C
.
UUL
S
4p
5f
UU
In
1N
Li.
~ t
.-
-I-
aD'
-0
SI
KU84
AEDC-TR-78-63
*4444*******444 gg g
4*
..
..
. m
4.4
44
4 **444 I
000
0000000000000000000000000000000000000000000000000000
~m
~n~
x0
0~
00 I I.-
ZQ
!22a I-1-
X TI~
21
2-w
4
2 A -4
dM 4lw
-
WW
-o u~
3-X
W I
-00'
~~
0
.J 0 OK 00 z
L
4Z
0 SIT" NJ -M?5 Ow
4 3-
J.1 33
0-12
-4
K U 0a
00 UQ
& ThO
SIu 4w 034f11
.. 4m -=4 0L.0
ONO*--j-Z
US~ 5-~
_ .0
x N .. 8*I1
jU-
I- IJN XmT
-4 -l
C!JWK:;
.
4'S
*KK -.
I -0)
--
MM e
*IU LC 5S
33~~~C
-- * 04
5 -03* 1
N0
0 ~I .1 I
*IA030 .. 1*
*.LL*-. J3. M OI 3
W~I K
V(3
'0
wwuu h
%NI0
4q
*4 ~
.1
6 4I
IAIA.
US f
MO
0-
0SI IV
'I0 1010-M U I
&- 0K-:
II-IA
S~A 0I
01
II
A0
11
**
-..
.E
.- U~8544
AEDC-TR-78-63
'asp
ix 'v
cc W 0 La 0-
a 20
@A
06.
,,
$-
z W
a -j 0 -'! a x
a, z
It U. T r U.
v rw ... 'a N.m x 41
'. 'vo
W Z--j 00 - -N W0zWu0_j
0 0 OZ IWO 42 02 -9.jW 0 ; z W W z k. -i a W'.. "HAZO mmmw #W zo-ONO IZJAWJTOZO;zx Cal 4 .. W 2 cc W toty ;
w0
lux T U. 4 .. -M u un 0 0. a W al .0 'd *a
Wo -j LL = . UWO -9 cr
c U. X,.Wz ..
0 -j
bj
tj D 24
in W &L 90 W19Z
I -i tA " 4A -C IL U Ic 2.910
WW-
;am . W
Olt
Z% ltt
Ito
-2;: Lt U
X-z
Oz
U,
In 0 a O K L.) CP a 't,
W W
n i! on
.9
uzkvt
-ka
IL z Uj -Wz
t'
XV
x HlizomwiS Q 6" 4 a a, m 0,
; LL Z .94 . 0 a IL 04 x x 0 ac
17 :)4 2 %A r 0 Z w::W,4
04 = 1XIt N -
It . 11 xf
U. 1 . cyx;4;1
W ru - X W
xWX .. , CC 0 1 2 3. z W:;,z
W-N Z=u 0 cc 00 fy u z x Or K z
flf":f*n:4
XMXA x
xmlx0 -4 CL x ; uj 11 p
.0 nr1%
1
01
=W=0O-X%X-,Y.C-x . 27 MO TZ mc lax r
VT C 70: .- JT WX tt 0 rNT
x a 0 raw -
ZX cu-joi zw% -W , a - zz
x x x x WN Z:.O.w Z:10:x"uuu 0
WO-ONZ --c
Z-OZZWO
;I ft
CP
CID
cc
a-
86
AEDC-TR-78-63
40
Of4-*&-ta~ 0::O-E4ml*
00
0@ OF :0.4P e In4P4I4-
NM
a xf
*4 p.m W
cc
co
U0 ol
J -; -W
20
0 Z ZI NO
-v
2.
-0
0 --
*
WIa .
U09
4 4M
Lai
-:c
0I;
4 40
US s: . 60
.0: U mN
CIAu 0*1 I -1 Z3
A cl
J..l
00 00 22-
*z me 04*010
al9-s
ZZ
a, z Wo
*ma
K WW
87
AEDC-TR-78-63
-0
0 *
Z 1-.3 U 0, Geo. .9 u* W % mt 0 -0 0 w N
9
0; z 0 % 11! l Z2Z bc W
. m
2c 4 00 22 27I q24;
mC
a
M-
'
0 zN
-0 0v * .J x y m
N -P4f 0
ml
Ui 0
wIl
LI~f
11
* *
WI88 4
AEDC-TR.-78.83
it.
00
u t. US
2z
z N
zo
I L 14 ,c N-n x.
0C 0 -o".i 5o w ux o -
00.w
i _j
w wgn
& a
o 0*.~
N 71,
c o o MZ!7 - -;--
*00
'I'V
--
6.
2 ,
-a
-jw*
c3 I%
-C
a,
0L
02
%0 VM rW o 4 z 01 wI(SU1 a
*05-5S
os
Iti5 A &440" .0
. ai w4 US~ 0 .3
-. n
v
. _W * 0 hihia
w
0 W
(SJM
'd-
-**~Z
289
AEDC-TR-78-63
w~ooeeMM *O
Ooee9oN
*o
P.:::=4eno
P. o
mooooo
Poaoo;Nm
4eC,
ooC
Owlw
-~.
W()
fnK 7,
-0
r 0
W -i IC
X, 0 2 -*
-
4 2
42W 215-W
-ea
*
2
aNx .4 N0 I
x.-
2C -Af4 W 0 0 j W 0--i(
-.-4-m-
ZO 0 0
Zu
) _ift x2 -
a "
~~~ Na ~
x-
. i-3
-i -
z C ii.
.d ! W. 'IJS
I-
Z -. 0.. !- . -0 h 2 u
AS owl0-
4Uo I
tfI
Ou IO-Z W u1 In
000-'*A~~ 2 2
10-.
Kzh
fn)W--K.-6
J~ u
4i-.
466.
U a
6 a'S*-N
*0
06
20
4JI
i-90LsSW.--
.JJ~
A ED C-T R-78.63
2Z2
2Z2 22Z
222
Ez~oo~o33
0000000
000
in
0-
'M
a-
zo
-U
-4w n,
MIN 00
:4
.
0-0 -: T
0j
~~~
I
%
On
_j
4 0'.
.4.
u1
r 11 x
*Z 4-
i91
AEDC-TR-78-63
W;U
OW 20-2
2 0
I-IN
z
Z4 02 PI-I 14~~~2 IT, N @2 02 wn y0-.
US x .2 M- U -
0 4 VI
*
% z
*60
*40
.4
;;.
. 0*
Zzo
2!
!0U2 fU j
. IUn ;
U *0
qit
4
014
*. % I-
oh 0
04 U.2;;
4-~~~~~~2 223t42
~I t
01 -mo01
sz-~~0 &;
01 a M.0P. 2
*2 **
0 ! '. 0.
# s
--
00we 22 041
40 w A
hMMNO
s43. CW4
Ro e
.2V
=U v
a04 itMe 19 z
a
m ENZ
x0 0 U-0@ -
_j S
-MN
ee=!tvIQ
12..NO vfI*z~ 11
UK
I .
*1h.I~
NO
WO
040
AW A- *1.M.UZA
20
hio 2
0.
Wu1024oO* W UOU S. .. N S.
o '. .N
**0 -
-02
O .
-SZU I
~h 3V -I VI
'.
.'
N'
%'
hi
.'
UO92hi
A ED C-T R-78-63
00
VZP
WIiWWWIJWhIJhhWW
au&
CLa.
-000
4
w
W
0 -1
0Z
0.. * 0 00
4, a w jl -
u
t Ow ww.WIL
0
'a w ONClpmii . 71
M-
Z:44
xa
m
!! ! m on . j
m--
11a0a0-t j
z vu
-&: m ---
S 4 :! 0. LO 7!7 hi m
Ow
v -.
-W
a IL. l
--
20
1 "
0lI 44
W -_J hi i
zoo
_we,
--
-1
LIoW-
*l
9.I~U. v
*U *
4%*.
e~.-.I.-~O
22*fJ4I0
*~))93I
AEOC-TR-78-63
K
mum~a~mm~~gg~a~m or
%0 Wz
Ow
..
--
I-
00 -)h W
z
-
xI-0Z n
@01
0Wa t 11
*
-2
',~l
---
hi% U IL
0~
~~
.. a'
0 N9 n-9,
.- 0
w 0-oZ nw
0 -Colo!--
IDO'
0 -. K h
a ZIA-
*J 0 U 4 *21 00
~~
-
h
-
hxiw2
0
-J
I 0-9
N 0
0 Zhi 0 2
~~
-4
-i
AEDC-TR-78-63
010
0 ow
fy 0
W
'A XZ -3
W
Tip V! 0 0 w m
-W10
ILA
A; ama Z 4K 2 VIJ
0 aX
00 to
%
Z-
0;
--. a- Ws - to 0- 0 UVI OnA 1.120001 -*-Wl %fmI
*
0 0
W*
0.
U.
WISo 004.j
959
AEDC-TR-78-63
IA
go fy 31 6
:ZN A V).00-." 0
Oa U,
Ia
-IA
-3 in
cc
TM
U N am 0 :;w S
-o C
tAt S U z . 1!A4 . ** W .
-W
;4 0 ;l 0 . j i-.
; 1"
WO z
0*0S v 1A0 S
h ,
. -
-J,
.. I..
. . . .0.-
fu 6
1650
0a UIUI
ey
MA in
In -.
l-
AbU
4*
ty
f"
5-
* *% 40
ft
aS
O90
A E DC -T R-.78-63
P. P. -P P. P.
P.
uwwawwl!6
www~w
www w
www
wwwwwwI.JwwwwJ
sS1 4 c
3'
4 0
aI
-)).I !,
"=0 0 ; -
!i.
.e 2~w
-
0-
a*4 zi 0. . 40mo -
Lo t
..J mIz
Z .J r
P. I
C! Ij>
-,"
I
to,
*
Z-17
42 01
Nio Iq 00 -0
:0 oU'UNm 0
w o0 00 N N
-j
j-
j
-
si*us
~
*S
1 U
~
-
Z u123 -
.6
4 N-
97
AEDC-TR-78-63
4EmEEE
4EU1m
Ix u I Ux
9
4n
l
w o IL 4. -3 -x 4 1
at a.
o-
11 0
o
..lw w t0-.
.j~
co
411
I. _jc2
LI
*
I
It,
*
-1 2
0 -0 ow1
-'m -L
to 4 cc -m m '.2 D. .0 42 w4t * 24 w 1 0x m
-~~ I.* w a !&
1 2j t*'2 0 oc
*2.,u.:!I .4-
4 Z
2Nf
o-
-411 0 0 4 11 4
0 0
-j~
20 2 o0 3. 24 4A ccUU
-e~z .. 2 w 0 IA
o 04
0l00
00,
214
in 00n.
Lm.I4---
4..U
4-
4 -
00I-4.J
I-.
0.J
0028
A ED C-.TR-78 43
z
#A W LI *~~ 3:0 %f" In .-i 0 !U
Ue
4 hi U
*s
osi
.1-
W3a
zo
Ua ; . 09-5-
,.
Co
0 Wz 4 -
~ Zj
U.sP ZA0
a3 a 0
U. gi
-n
z
S0
~
4
P a IL 31.
.*
u~m OW 5
~~
K -.. -;m U . M. -.
J...*
I Dh P.0* a
~U
.0--0990-
AEC-TR-78-63
op.o.-raftam
06
a60
0 L&&&
L0
LC
9 LC
LI
.4 z 4Z
VI
N
-
cc
V) z
a N 0 so. LI WI
-
z 0 u 0 m!. 1 x wU 'S ! x. 0 0 -1 o0
C oI
40 0%
x , 4
z
-
l -
0-
a x #4-
-W~U %L-mz wx
m*
0,
or
a1 02'WO
OZ VI .z
04 !
Z q 0 m .J 12
0!
0
woo
4'
tin
~~~~
-.
~-
0 a- X0P-3 z0&0U .
VI. . C-ID
5P-I
IP
u IW
1, a
~5
524-
'S
100
AEDC-TR-7843
I~P9Ehi949-6eO..NP94.J94N6tnNm 9--9Ip-51e6@@64e
Ill
N 0@ 21@
0.
z M-j
0 N42
25 oc000-------N2NNw
4~~~
$. - i
.h
I.j UN 2* 9 9 0.L *0 9* *
z 0 ,21. 4& * I% 2.
4-.: 2
4 40
M4
U h hi -
44
44 21. m1
6%~ 1%
04
h:,
222
xJJ
,O
11 2 Zhi
-.
00 .J20E
0 00O
ft I
OI hiO a
-'
Z0 h-CL -%c.
_Z 1 03
3 Z
2.0
vu9-U4Z
~
.43~I
-.
412
IA *~0 MU 0-0.0 -
~ 50
1?42 0 S.
*hi h4 ..
*% P1.. r
.~
; ;:
. .. 0
0. it-%a2
i01.%.I PC *0
.OW '%l0 2
4
-.
44 ZU 0 2.0.9
2U. 4 U 4h %I-O
.J~04
hiW43--
~* - ~
.1&P
%406Z
34*
-i4
-952j
m4-4 42
-4
0m
. a. 2
O4j2
U0.44Jl
-h.09W
P.
~
1
Qh4IO
%9U
Ow
Tit:!~*~.-. -4%.2 10 hi
x 09929N
4UIUJI. O2
-0 0
44. 4-ON;-
0 24h. *K WCK
h.2 a1 Z zW
J ii:
4Z- *.... .N
U.3 .
-'
-'
4 ...492 X, *
44
4444944444444-4422
x-~,;
-!4o0 IA
K~*1* 'L
'
a -
5..0990,5
2 2 . % 2 U:
1
004
00-zc
22 Cl O* 2 22 o2O -02 2 -2 2e
-9N
iw
*9**9 e
OQW
Z- M
-h* h0
* @
/9*.
-m--------------------------------
AEOC-TR-78-63
06
I.
A %. La4 0 IfV . c w
0
"
I1 3a.
wx 404
-
~ 3IN
iM: ol -mwm U
-
0 LM --
Fn go
-40
~~~~ ~J 0
c o -oM0
Il
0 'A -
z -
a-a 4 4
P.--0
'0
*.0
-L x It: o6 0
-
0 o
U A
Vac..
6.
a.
V 4DXo0 mo -a-.
.0
! -
1,W2 * z H
VO M S ) 3 N N-
a..0Z
-Z
U %
x -~0
0I
Zw
o a a. V01
jc
I.4
w-a
a
4
-j
U
.2 M.i0 314n0 * -25 05
*aji 0. a..ao
2a--z aoz
C
r. a3 wq -!, 5
-CL
-am --
g owa l.w
-. 4
10
o a. '
a-
m2
.
a-4 0W~n 0
a- N-w
W--)
00-
*u*
a-Cup
25 0
I.
..
UUSS.
5a-
40.
U102
AEDC-TR-78-63
.J
jIah
j Aa
-J -j-
ij
P-
ee -a..
a... aeaea.
I a.. &a
a. a c.a
m.~
IA Co
it x t
lo
OW -,C ON
31*1
m O e cb w14a
N-W
--
2y
in0
.0
VI
0@
it 103
AEDC-TR-78-63
x 4 2~ hI
S.
J ( C,
.4 (Y0
-0o 0 0 -02z
..
N Z
ID
9.
00
0 0 1.
4 * .9c
4 0 0 0 0 -
00
*0.
* m U 0.* -1 000*9 .. 4 9 ; ;n
0 .8
-p
-1t
w. U.
kv -0. '4
;z T. 0, 0 LbN
0x
L4
4 NQ 0
b548or%
z~ :
M
9ko0z*99N4 0 *Mx
,,
U
IZ 401
0cc
'4
2-
N@P. -.
N:NI
0.
c * 0
0 Nz 2 U b-00 aE~
-r 0
O
059
48 4 % tI
0b -1. u . 0 *W .4%
. 04 0 0. 1 0
0 N
bNN
4 00N0
. bb-9Ia 0bi4ol Nm04 n N
iu.10
w. fm46 CP .
01
*I
Eo*imS
b-
*0
0WW -0000410010
000;
-I
Cob-C.
-SZ~b~bNS.bs% 3 0-
.E
0E0
EZ
0104*
AEDC-TR-78-63
~ ~00*"0000@N~
4444 444 444 4444 444
NU
0N
*E*
0A~ P 444 N
* NC 4444
N OS g 444
4 0
0 MCI 4 444 0
~~~E'1~~~~~~Ec 4
444 mm~mm~ gg~~ g~
444
~
I U'
0%0
z*.2 In n!
~
0
4o
0 ;T
00 10-
I 0
0CU
C -
0 V@ . *y P
002 ;3-1
'A 3 L Zx4.0
CL 0- It-gO I34 . 4-
ty a01In 4%4 0a
* *
0 0 0 0 3-.533
~
00 14 T
e a-
w-01
(Yw -
.0~
~4'-Ie~4 m4
*-.. g O
.4.~S IL g4 -- o
eq-
0x
0@*~Clw
Z-00000u
~h%1
Z 0 w 2i
-
*3-0Z0Z-
.8
-j a0
i S ea 0--0
0-3 N.CL4 20
0wNZZZZ ,3
z4x&
4*
0-.k
0.
0V 4 4
4W4MS5
*00
5514-9-44 *N
%S
I@3-0
)-
3.40
N10511
AEDC-TR-78-63
-.
:::
ma
30.3 maw-
4 a
44
C~I I
.1
0L
It. o X0
4 .3
3*.0
In
*~
~
0v
IM
0~
~~
A 05
--
Z
I-
I.4
.1 4L-
a.
S -
%AJ 1
11
A1
.
2 W..
KD
Lf LA Q-0 S.1 o
(a. U z
*comNx.Z0f
0 3-0 N I m Z;4 1 th ZU
S Z h ,
I1
InZ.
x4-
I10-0.. I _
NM
S- If
M IC .-
imam
C T 0 0
'a c
.a
106
AEDC-TR-78483
0 f0u
of
*c
%~~I o~
dc
*%4 1, N U. 0h.*
.
-C a %I.
'.
m.
aK
IW
I.
a 0
I9
~
&
a. &M 1 .4 I z. x 11% IW 0 Oz-P
20 M0C
2w
Ex 1 1 z
cca M- CL (9 1 19 m z 19 *~t
a *a . -u
0 I
0 *
31 3K
a6 a.
_ % 'I
'*o K0 a i.
x OK " 1- : z ..
" I.
a
4 * K
~ 0.0'~
300 a, &3 -44K 10 y I-.00 -J Z 31 at NM41-1-0
~ 0.
0. 0.a 0& 0*3 &30
W 6N -
0.
1(-
* mISO *1 .30c; 69
N 190 N - 3
2T
y 41-
0 ;4O.%Xi -
U10 .
T6m *6
41I I 40 0'
00T
34 C2
I. 1-
~~m
*3U-Z~0K -00
. Q1..9.
1-
(1-O
0 0ey N
".11 PC M
-0.0010. 0'.1100
3-
-..
-0.-03-
-32-
W-~~6113*
1-l107-
AEDC-TR-7843
3'.
ey
.J
10I wS
I
in
V IS 0%
V"M 0
w tm 46 a w
u
N C
CAM 00 m
.0 0
Slaw
- x toS
.0
_00
W.
0. MSEIA
N
*U4 *i:f
SX
* Uz -u .Y m M
-S
1..
4108
AEDC-TR-78-63
o00000
008
00
R00000
0000
400000000000
00000 oo**
0000
- -
-06
01-*
MOO
4p00 004 4
O.P.0t
0 O
Q ~~~ ~~4
N
'a OM
P.
N-
44 N
P.
M vC ~ '04 O
~(I
~
:(03
004M -~ N
70
4 01
00 0, 0.
000
IL
mp.O PN 6 OP 0P NO 0EP0 .K
Om NNP.P A CP
hp P. P. ol In4a 10 V. 00 Imp6
w -9 a,
w
0 0
ON406P.PCOO E C'M N NO
6.
~~ ~ ~6 ~ 0 0 ~
4~~~~I U~
. O.P *
4 P o
N60
044NP1
~ ~
toP0
PPNNNN
..0*P 4
.O
P.
lee * 0 0
'r
N 0 tl0
D 0,0,P0 6001
4P1 4
CwwmwN0O04
0 0 '04NM
4 4D A
xl
I-
4. 4 N O . E 04N0P f N 4 N Np -M ' l M
80000 P.
000000000 3 M O00-0 v
0 I wP.c:!
CA 44NC~'4J
66P44PPC6
Ii.N
....
064109P1
AEDC-TR-78-63
o 0 o0 *P. 0.
NO 0 InN 0 0 000000000 0 0 0
Vl oa
V00 c C:
, 0 0
0M 0 0
0000000000 0
0 0
T -I
14 T
0 a
*
0000o0 0000000000 N
0oaooa a
UP -&' N 0
0 _j
4; 2N 0o
4N
4o00
0 la 0 : 0 -Zo ov~
N 00000
0 w~& w N
oo c: '
w
m i Io z 6. y . IL o
It. -In
56
0 N 3 B 0 3o
c . 4 N o. o hi h c
A'2 0 N o3 N
o 0 N Mn
IL 3 j ar
-
a 0 0 PP M 0 Noo 1 n U 0. PP 0 4 0 p N u 0 -
O.
N . . . . . . -~C1.C4 - - - -pPI~n
0 *4
p.
4 UP~~~
-Q
& N
A.
A.
00...
N0
4 . n0
't In .";0p 0.
z
p.
a.
U 0K U m
-W
nmN
Ni -
f3 N
a: 0.on P c4 wn o. 00 0- M 00 - 0 rw 0000
~o2
cN pBU 0000 0
-y
N W 0
0 &
4 000000&000
w m
0 p
0, 0 oi
P. Nl
to p.
o0:
.4 NNNNNA
o0
N NN
N0
vN
.0 0 u j
2 A
0N~ 0cl0
0u
4i
0 0 -00 -0 -0 - - - - -000
IS
' 0 00
m1 0
&o
Inp o 0
I N
o 0
0 .0
*
p. 4
. .0'
moNN P 0P
POn
.0pI00
0
1
NMNE.Pe p
l
P. T. 4
B.
B. 5 N U
0 1
.. . . . 51 5 1
110
AEDC-TR-78-63
ONO0 N
!Oo
top-' 4M.424N
10O N
-N4~4~0N
04?m0-M4U9
44 Ic
ON 0
N04
aNN41'E --
44
4CN 44..4' 61
0N
"fyy9
ry ury
ry'y y
~N.O00000
- - -
o0020000000
0 oo-o0
00
000
M-4%44
4 N~I~N0PN
'COM4
-N.
o MC
~
0
OMOONCNGI- -
0000 000
00000000
O~~~
*N1.. M
..
e-
0444
It
C. flN04Cr 0
0~00!
C:
N.r 4
00-04~~ Om
Nov
E-0
C 00
00
00 0
1. 000000004 0 N 01
IL
* 4
N N
NSM~fO4 4N 4o
N-.4
4-r04 ..
AEDC-TR-78-63
190.lfO 4000.U~4
.0019moo 0a 44awo"M0190
C'N.
m1
ON.10-NNN go,41l~4 -.
m m
A,2 m 0 4
Oon
0000000e
e0
00 2
00
000 OOll~f~f
g 0
~ ~ ~ ~ ~ ~~~~J
.19,0 M4 0N' .m
~
00=O-
0.-19 01IJ -
O 0 *ro.
1U0N N..-0 a
Ow *440fl*010
-0-----
0000000000O
-
0000000000
0 a
o
39 c.
000000000o
00 19019
acO00oooao
o
m 9:
OaO
N O9 0 047nW.J04 199. c
ONI-4~I0 9 4 o
O0000000000 0 00 ,G
00e0oooeoe 0 Z.P9.9' 0,
a 0000
LI
OON
INO
-0109 4N0
*0
044%0
0919 -
0000000000
000000000
0000000000 00 0 0 AM0000000 0 0
* 00 0
0-00
nN yNN9
01 N
0...o o 44N 0N m
0II000000000 .0
oo
0 -
ccU09290 0 m09.9a090 a
000 00
00 000 00 000 0
~~
0 00 Cm 0 00000000O N0O0co-O 04 N .. 9
0 0
0 0
0 00
0 0
000 0 004 0
199.4000-040 4
-0N
O: 409-0.041N19004
0.N1:M0' NNNOn
NNNNN
M-ONOI-00114 NNN991111 ru
. IV
N 0
0,
:0l
o 1999999999 0
NN
NN
NNNNNMA41
I 5a. x v
112
A ED C-T R-78-63
.0M!4JNNNU, M !M ! . ---4P0D 0P V .. F4 . N
4 P.N a N M
-. 00.t@NCE - 4 P0 . M~MM444.44.
N.U M 4' # m
0P4 O
Pi.CN00 *~ 0
00~N.0 0-M AA
P.044 AA .0 *ON
0 00 c
~
0
' 0
0 0 0 u
-0 00 0 0
0 0
0 0 0 0 0 0 0 0 0 *000 0P-Pb0 *
0 0 0-N00
;gv;;V NP... *
~
oUPo!!o
SN
0 T
00
0 0
0 N
0 4 P M
N0,
Q.
0-U 0
C0 4N
NP0I
--
~C
44U. oCP 00
P-4.1 0
j00 0
f00 000
N NIc, 00
P-SO, 0000
0000000000
OOOOOO:o
Zc
j'.j
.4 N~ O000U0'~ P* . O
z 1 a0
00 700
0 0
0 00000 0 000 0 0 0 0
0 0 0 0 40-A 0 c
a.40 '~U-P
00000
000000
WOO
00-0000000-
M0
000000
0:0-
4 xt4 4
CI0-...
4 ez 0 4 *0In0I
a~NO0 0 0..I K In I
0N M ;mOOm .N0~040n~..
4m
0.0, O...N OC
w I
0 .P. 040
Ib0
O"
00--.
N~M
*lCO-113
AEDC-TR-78-633
-.
~~~~
at-
SO
00. .0
.**..wf...S..0 0 0
00
..
OO
Of
..
N0
Of
NC
AV
NON
~OO4
00 0000000
.M
*WC~4t
00.* O
O00
00000
00
00
61
Nt4t0
g
4P
2S~t T NA. ;v
CP
0 0
MO I
03~
tOM
it .
MP4
S0004e
*P-lut
.MNN 060
NtOW 0 0 0
. . 00 e
000*0
it;o, nr'tA,. m
IMW.C
JS-.'I44cc~
49
U0
0000U4a-EN
fJ-0C.4C't., 0 9*t04004
*4
00000 0~~#9.0000000
a00400000P 000000000
0 0
0000020 000 0
000000 0 0
0 0
*~M
cc
00
co
00
000
F.~onO
.)*1--N..N 4 a t-
--
--
---
-N.... ....
--
..
.4
P-NCOOO,'i a 40 o'NC
A 1 l~a0N P n 1ME-I -m
~
I~--O~l
. Vo0a
NCI t
-Nm-m
f...4.4N01
L I-
NC E
041)0
OON. cc
.-
E 04.I-C.'4 ' 44 m
Nm444
,g;
iji
em
114
A ED C-TR -78-63
4N4lN -00~.~N0
C 4 ~~0N~44N~P.ID 4. N00, 0 4 04
40.P
4NeJ44 N
600000000 0-.0'4'4r-20
4 I
000000000
Il44
4~.I NN004M 4~
c *
0 c NN
4D4.fiNN
0000 000
00 0n 40 ; g0~0..N a
i440
o 00SN *MiN
O~l
440 '-0)0
*I--.0S04A
04
..
* 2 v
0000
00
0 U) 4 I 4i
. .0.U
? U
'.0'0499 90-
U J4
4oct oo
N00 oe v w
o o
ol0 oo
00 000 0 0 N)U..4)U
0 00 0
0 0
0 0 0
z))
000-.ON- N4994U)O1145
AEOC-TR-78-63
In
N
~
O o o
92
N c
~
N00 0
~ No.NN
!2~ 0 00
0 o
0 c
~!C~0
ec 00
n0! o
2N c
-~~~M 0
!! eooooooeo
0 oo N
00
00 0 NMoONeo! 00000
00M o
0 N0000:
o t
: 4-41: tt
luig12 a c . -
,sm, cuu 0 0. -
0 ue-.I
.,
In
0,a
EN %0 4 .. 0%.C
0 0 a 0 0 0 0 0 0 0a 0 3 M v v&
0 0
0 0 0 00
0 0
0 0
0 0
00-0
NVNNNVV
4j
O-4
0-.
01
40M.-.
O-
OI4
QO
La
EM ....
ON O
.A~0N
0M
0C
%.-N4-.
4~.1164
AEDC-TR-78-63
Ill
44 *~
0N,4N:6N
0000000----
NNNNN
NM~MIMIMM
~~0IIP-:Q
--
*440P.1S2N4O N 0:704o
p.MM&
0N MA&N
NN
U*l44444
444444444
co
eo 00!
ec 0 000
~~
04N@O-.
NPI-C~~-
M__7 0I 0NC
P V40
40 OP.0*
0 M
,0N I4 44444
ovN
0000000000, 00
N4000000 0 00 0
00000 000000000
-0000000000
0000000000
MOO P4
(.0m P4.4,
0~4N
N... *4V
oc 0 4 04 4
U
@ ley..M I"AN
JP-O
.*G s
-. N0
000000000
~ PP ~ P-0
mJIP.40
00.
0000000000 2 0
U
OOOOOC0
~ 4 ~ ~A 4 ~ 0 M 00
C0IV~@#010N _yN4
~P4'0. 07 M@040
M4~,fJ4-0101XI 11
44
-c
2~ 00
~~~~~~
00 Mp-0444N~
10 z =m~NNNWKNNN'.INN
00 00 0 0 0 0 0 0 0 0
*I~~mgg i - -0 -00---0000000 0 0 0 0
U
eg
.,.
'r 4
M44 aPp.
44'44..C~t
0P
404I
4.~0
I~N
117~
AEDC-TR-*78-63
AIIALIMll 0
0.
l g
~
-
g~~~
~ ~0I.A40.6 0.~~--0.
O...... 044MN-OON-Nf,
00
40.eC0U04N0
-t40--
E~
.00 0
()
-o
a,-----0
40.0
4 v .044 NNNNNNae.JaAoEJNNNE'gN 0 0 0 0
NO
N..O,
P N0 g 0, a m
0 0 0 000 0N40.0-.0
0~.P-..0
0 0 1
. 0
0 0 0 0 0 0 0 :0 O 4 0..-
0 I0, 4 N
N0 0 0
,I 4 4 4 U 4i .0 J * .0 0 0 MG(
co
*coca
NeN.000Ua
C. 24 U
00a00000-0
1i *0 ~ 4 40 ~ N~- 4
E NN' m
'cc 4
cc i ~ N NOo~ OO
N40c4 700
0
0 0 -0 0 000 r.0 _ 0000000o 000 0 4 00000 00 0 00 0 0000000000 04f 00-44doe444;Oa; ... M0 C ID 4E' g .. C'NJ1 043 0 $ m0 0 C4~ .0 0NNEQ0
0 0g ,0 coo m d N .0 0
.ON K
0M4
M*eo
..
0M.&44%ON0,0 .0-..0IwN
0.00.
N--N0O N
N aK
~ ~ O 4~ C
04i0
~ .0
00
fl 0 K10.
~
00
~M 0 0
0N
0.0 C
A .
Z 0 I 4. .
aJj 0: .
# N N. N0O .04.0
000000000
00
00
0000000000 000000
ON 0 mw N ,C 0 N U
CN4000
v000000000 ID0~.0 f~
O4
004400=
Ln
. CL
00
-..-
N4q
.0NNNNNC
mmaoo--N
cN -
--
tN
AEDC-TR-78-63
2.4~4N04
00 N4 t 0004
4 11 4
fD a,.fN
00000
000 0*0
Al 0 00 A.
0 44.0 hi A F 9
4A
4.PF..FMMMM o~ @ @ *
~~~~ o@ @
mm;MMFMMNNN * 0 0 00
NNNNNONNNNN 00 0
00 000
1
NNNNNNNNNN 000 0
N 0
00000
--
~
.
~
.
000
om
6.0
In
4. Q....NFF 10U.0-.1. .
........
...
..
..
...
F. m0
x
00060-6S P .04.00Al40b % o0F.64MAi%06 "0 00F.66Uo NM1114 0000 & 00 0 z*Z !2 fZM *4AN-:0, ' oe o 0 0: oo :'Q 0 e O e 0 O~o
NVA!.lt~1
11 4 -0'
06
t'
06PF10FO
4
M ..
N,-Ii00.
44
'IifI406 -. G. .
M"-00F.%f04P=
1~~~ --
0 -
F.J4 ~--
o -
. . .60 . .. -------
.-
coo~~ - It.',.
~
- -
~
--.
.
.l . I
- -
4=:*N4 '* )
.'t
-~.4
0 a.0
S. A 000.
-a (P:;
~
4
o OO oaDT ..
o oe Nm ..
..
MMF.F.4 h0gq .. 4M M
04
"o
4o 06F.06..~~~~~fl66G.NN
9PM NNFS6.~rNM W
N~ -. 064N00mm .
0M0FsI0N0M 0F0--.
M.N0
60.
..
P-119.4
AEDC-TR-78-63
.00 eC2 ON N 'C N 0M 01 ~~ 0 ,,, 1 N AINN NNNN..:N .=am09 c~am maa C-I~g OO C O CCo' OO NtzT 'cm&OO4NN'C4 ItN 0!CO . .. L'.., . N
*0N N "j N 0N NNN'C1
1 0vZ 0
f=44N-N'C
ISJ
m_:ccc . 0 0
0 ..
04N
v o
M N0
LM o.
-
oo o0 co N LACco 19 o 91 co w CPm N N IU
c d)moa: N
o1 N N o
I0
06.,.
a
cC
Booleo 0 0 co
*~
-A
~~~ ~~
.
a 'C v a a. e 4 0. U~~~V
C 0
0 o N 0. 2
. o
m N0 N
"': 1 4 . 0C a. 4 w41
4 .
90
N/ Is z a 0. 1 4 o .J
.0''49 c0 * Nm m i a s
409C5'0 -00 a
~ 0.z
'
4 04 2 C a F.N40P)0C-''N
-
c oc C In N 4Do 0,1 NN N
000
a4
o u .
l
11
Pa.S
C,
mC
2 U 0.
N a) 0. 0
0a
N
N N
4 cc 2u acoao 0 a 0 ua ,0 . N -
.NN ..
N . .
N mNN
N C N
o -m 0j 0 Z 31 -o 0% la
4!
a,
It It t9 a,4 N I -
% .14~0
.c ammD
2 20 m
LL _j
w0
w9
0 'C
'C 0n
x 4
00
0
m
m3 N am 00
mm4
m 0
mmm mi 00 m00
Ne
o'a 0
S
am 0 a
P0 3
IN U a. KM~
4~~~~
-i c in *o
N w
M0"
C -0.4'N0 IL
4'U
O NC
CSSCSO *A0PS w () .. l a
N N ~~ ~ a ~ ~ ~ 0 ~ ~ ~ ~- ~
lb 0 o a
0'C
4 0
45C5440
'CO..PS'..
'
4444C
4'
'C
'.
4CNN.590'C
3'12.0
N4
N54
AEDC-TR-78-63
0000 2 24
--
MNO.NN P PM
44 4N f.E004 4
M444
0 0 0
~
N1Pmom.0
40Cllo 4-0.
a 00000 -
00
o N-~m
a .0
40
00 1PP144
0Ji4
i.000000000 4444 44
00
-NNNNI.Af
3c
~*l.0i~0 0 04.N644iiJ
Em4wlcPO6N
.000:00iS. :01-...64000 :.4fl0.0 0i0dII4I. 0404f iflfiNoiOMOO WI0 MN& 9tCON 4-.444U= tN e.mtNm-0 0 * 6!=4
0 0
o0090000000 *
0000000000
0009000000
0000000000
e 0000 01 0 0
00 0000 2
,00 0
00
~4
. .4N M nI
me
I4-
.-
MN 0
46 .
fl InN4 -aMC0
~.P
. .. 0000000000
0N4
.... 00!00000000 maa WOiO C0404 0000000000 oo001000000 ~ ~~ ~ ~ ~ 4..c4..P.4 V UNOM~ IV-M44-P N2g *4 * ,0, InMM 4 4 -
.. 000 0000000
..
. 0
~
4 *
...
~ ~ c00000 P-i ~ ~
8M4 04 .. 4 44M
PPN-PN-PP 0 4 NPP
P.O
01
c.. IV.P~4
SM40.PNJf0NN4..4y *U atf44 0 4 * 10 *
.064M0 OO
40040mP..40
4 04404 N00NQ~sP M
w M-Z
P 4 9 10 0* OJi
00N4P4 Mow0N644m 0.~4 P4*0. In44-2-PO 00 NMMM .44 a IV4 44 0 wa& 0-i2N NNNNN
0 M A64
4 6
WNNM
4 MOON
N -----N N
-----
19
2, 4N uA1
ow 40
00 0
-0 M
0-
000 a&
0000000
000000000c 04P4-04 44 N 4
00
00
0 4
440-It444
I.
66 .~MN.
NP-6411P
04
40604
40N
40N41214
AEDC-TR-78-63
~ ~ oo ~ 000000000 Oo I o~~~~~ o oc
11 4'b.0 0 M J, 0.0 O
00
0 0 o o 0
-. 00 0 o
0000000000
0000060000
OO 00 . 00
00 0 0 00 00I
o
-n
.0=1 @-N f N.
0~C
N90i4O*I0*
04
NNNNNNNN
NNNNNNNNNN
NNNMM
A~4 A ".
44f
N w o
*40
4P.NO4ONO w N N rm o*CVN.C0 4 4 4. il
!~N040
0nf.I44 wc AA c0&
V4. c 00L0 4 00
gON4p
00ooooeocc
ceccoaccec
1 01 . ~* .csO-Ztim UI.0
~ N4'1
-
0~gg~i 0NC A N 0U NU Ni N
~ ~
top 4 0 NNN
0-N40
t-A 0 N N w 04N..J N
4 0N N --
NNENNN
JNJNVNNN
NNNNNNNNNN
NNJ
22
AEDC-TR-78-63
%n
a N;Z&
me
~*f~.lCP00 0 (PI~~0
N3000 0 00ggNO-0
N~4f
NM
00000
000
0 0
000 0000
~~~ :aooooooooc ~
000 000
0. 0 : 0 0
0 0 0 0
0 0 0 0
ie
4 I
me I ' 1
IS
II
ve~a
*Em4 00,00000000
,0 4
0 z* *Ecloo 0000000000
0. 44EN
wc0 00 0, 4
000000000000000
S.o
0..iI
g 0
(P.
01 N
f. 4E 0,
:N * 44* 4 ,0,0
**E * 0, 444 (p
40
ooceooce
0 ** U I
MOM 400,0
S d
'
0 N~~M0S 044000P o M.
d4 O - N 0
. .
M0 W
4 4 Z
0 4 2 a m0 f&Dy * k ul 0
i I0.4
4 0 0 4 U N A 004-0 N uo nL &0 -0 ,k j
W : 4 ! - N.
N0004.05nk -0 10f.a
sseo
M M0"
t- U o
~9Ia
A -0 # 0 00 40 0000 'a
;NM40,0J7m
M0000 00 00
123
A EDC-TR -78-83
,,
42
AEDC-TR-78-63
7TOOU
,77707
*7
ONA * 00 u
b 00
MM 00
NNAJNNNNNN N 7 00 7 00 0 000000000
z N 00
NN oe
N OO
NN O
N 0
't 'r
c N
NNno
0000000000O
00200000000 .. .. NN
00000000
N.04MS.sIONN gl-~P NI * 4e&I O : 4550 4 001,IS 09 v p~.q4 U 0 .. W~ L~ ---A. 0---SI --. 0N . 0055 ON~b~d~~fl15coo o
40 00 7~~I~ a4-PP-04sAp.m N.EpNIN i4OII Z0 A.0UC .N.0...4~ ~eP.NU~N 51 * 0N.P .0sNM 9E U2jW.0.j0Q V7 goI60 4c 0 *r4CO c* -N ... .f I
~~~~~~
0
000000000
00
000
00
00
00
00
000 0
MMW.J
00
~~
M e 3CM ODO.M5 E . . 00
. 1ItE WC ItCM . 0O . N 0 . 0
t . . 00 . . 00 . 00
f607
; 4
~ ~~
.. . .. A0 ~ 0 ~ - ~ .0~'s'E0 2
-q7N4 MNcJ0.
0Jjf. .
5N75I4_
J * mmmmmmm0f 45~-0
J A.I4
2 49
144..0.0.4.0
N *NN e
* 0E~0NO.. 2 m l .y4,ZO
*N..........
NN0NNmo
V 4Nmm
mmmmmm-44
*.415
4.MC'V 4444444444
v4 P..0-P.c5NNSw *
N 0~ N N N-NM -. NI O ~40~ 4 ~ ~ N?0. N ~ ~ ~ ~4dI~1 ~~ ~ ~~ ~ 400.. ~ ~~ ~ ~~~~. C50C0 emA.., ooANNM 5LNNpn.% TMMMIM c %4 NK
It
m 0 _ !-
~
_ Te N
A..ANANN VUI
J, NO , __
~ ~
fN.N
12
AEDOC-TR-78-B3
iy
1. NNNNN 1.,NaN N
NO N
f. NN
eeafP, IM~4
e 00000
0000000000 -NN..O..40 Co4 w~w NN4 N NN44M 0 110018NN NNNNMOM 000000,1 000000 40N @4NN -00 N4 -. - 1184 -0N N l 0
40 4 0A 494440Ue4 Ny $4 0Z
M440.P004 omNe on in
:4.Om: C4&0.40~ 0 N
N.7-..SIA
0fA004 P, N
4114
44 6J48~4
~N
1111111 41414441 NNNNNNNN 0 0 00
N:NTNW:~#NVN 00 0 000
NN 00
NN 000000000000 00 0 000
mmOP3Immmmm** O O
0000000000
4W,0m00 m00 Ot.4vo0N*60 000N444040 *0609-.N.N0N._44NV
0000000000
a *UN004A0. O*N-04n:.40 N0ON40
0000000000
0NN4S 04 ONN-40,11 0 440400.-..4 0940 NN044
00001000000
0NN O0 0N' 404 N4 NN40mwoowN WN*4M44m.04
0000000000
4Ne . elf"NO 04N000.....N:44000030M *4M4N0OONON &
00000
tono; 0.i004 402N0
0000000000 22ogeoooo
0440N-mNg44 0.i0..0-. *N0~.440N op, mi--44 3,VA~ m.0N4 V b0NOF
0000000
0000000
000000 00I
000000 00000000 -4?-NNNN044 =*04-.44oz .0"0 N 44O&N000Omo _0 0-'eON4..04@4N0 V0:14171..4 0E~00 I.0 :01 mm
0000
oooooo *044O...
0000000000 N~e
0000000000 ~N0N40N~ ~0M40Eg..3 044N &M0440 4.ge00~ 0Nlo ~M M,.~N 0mv 44 :% zN40:zN0'lO
07
07000
m0 0 -0
4040a:440NN
4 0000000000 0000000000
40-N4 4. M-0M000m0m
Wm:!N00N N
N N N C
0 0
N1C000 00000
0:
Nm a0
MCP 0M
1,p04 ' @ 0
M0 *00
00000z-a 00 0
0000000000 000 _000000000 00 00 4 : 4 0*4 4 t 404t 00U54* 0000.0 0009: .. aN 0004 44-NN O. M40Nm4 ONNNNW@000 . . . .
0 0 00000 * 0 4
0 00 4 4 00
~ ~
0 .
00 44 0 .
OU
04 VN"m r 00 .g4..P e-
1 01-1 C4JN0. a0
N.
lOm 4N %
NCN-...)4mmm
e 0
4Nt4
4-4
440MN004 4 4 0 O -0 4 . NI
n n0 440 0.00.. N
Lnn
L N4
n in-.0 M 0
%nLnU t00 O 4
.. ,4
4 0
!!-g4
44
~0~~404
04
C~'f~0~
N5..N~0
N0
4-~'0N1264N
A EDC-TR-78-63
SS I
ace I
P.00
'O'S
*S0
49,0
ii24P 45 P.9-0D 0
0000 0000
=ococ
*,. 00OI=
"0400
4090 oooo
e e e @9,0,,
99120
AEDC-TR-78-63
.6 Jo
...
tN
WN '"-O COO 40 0
P. 'r
0N0N O
o o e oooe~ocoo~o .....
!M*Ni0~ -
4OeoeeeOOOO .....
*~N
a. .0 0.:a*
0 1 N 0N04
'A00, 0-4
M O0 ONMMM*4,nPN 0 . 4
0 0
0 0
4
0 *440
*0*000
*0-----
a~aM4&~
---------
N
-----
c.040--N*4
----------
NNN Z N
W*0 It.oo V,
0.
IC
*M
NP.440 o~ 4 N
*0 o
04
40
N o o000000 a 0
*.
. *
.0 MM 0 0
440 N 0 N
M0
~
c I 01... 0.0
oe -
IL
4 4
4 0
~
b, fAJ lzj O m
U~~2
.,00 -
.4.a4.. ...
-C.Ooooooo. Q--
....
...
00000N000000000N100 ...
...
...
fNNNN 4 M 0 M
.. t NNE
aca M0 0
a.
m . . . . . .... . .
. .
. .
. .
f"
Co&&o41
O0
wc N
:a1w
0 0c
N.4 c
~
-j
M *4 h
*4441N..0 Zoo
*OO..
4a
P4a
N:;0P
4a-Oc N
Z04~
ON4
t440 0.~
a m h ,,
aa
a AM a 49.a0N
AnN
mm Ma4
NNN O4
NN
NN
NA - - - NOO.Z4
- - Z
7 N a.NNN 000,4
9
NNNw NNwmccc
~~~
a 0 0, &V a 000 aV 0, 00 m.-. 4 04 4 N4 0 aoa
~~ 04MN
0000'r04 ) ONO44,a4Nv*ON - ao sa~a~ 9NN;4O1.;NN a
cl2 N0 w m
2 . .'MOMZ4 DC 00 C 44
~
... 0-4
4 4D 4) 0
aa 44 n ",.*4 44 '
Mo
0 .Jp-
O4 N N
4.-O N N'm
. N A
o0 N Z4
zi
Of 4; NONONOEl
AMN zz NNNNNNNNONNN
g.
a.-aaa a
NNN
NNN
NNNN---------00
0000
2128
AEDC-TR-78-63
44
* ~T. S*
*0*..*
0
ss
... O0t
...
0 4
...
... t0
....... ----
0 : ---n
0:00 00'a00000-.p.
on000400
0no00; 0 A0 *
M 4
""'.00 *0..N0 NN
00
N.
C -
4.. 0
00-N o00
46-.44
*0-NNN 0
0000N0t0, ....
04 0. ..
M.N-.OOUN.02 0000000000000 . .. 00 .
O4MMNNNNN 0 . ....... M:
0 00000000000000000 ..............
- - - ft~ N"
- - -
- - - -
- --
~ - ---.
-------
-----
----MIA -. n I Tp
----------
------------
n
0
*
0~.
z 0Zm
NNNNNNNNNN
NN.NN4NN50...N00N0...44
40000-M04990NONN NNNNN--n--
----
09000
00000-0
00
0-
0049N-0
0490N
N NNN-NNNN
& 04 c: c 0 m 0c m 0440909CN004
008 0
~ ~ ~
000
0p ~ 000 000
Na 0,
~ 005NN~ 00 ~
N4 N
~MNCC 02 ~
5 C 250
00 0
500
000
n0
0A
000N00
714
'00 n
-0 mo
o-N -. 4n t, j,0.0N0
N00
4C
---------
N NNN
0000004
00000000
00
0004444
44
444444
4440
99
000000.4004 wo
0...OOMcN
0N
04
m N4N0. o
940.00 CPO,-O~
010N-
.. . ..
. .
. .. .
AEOC-TR-78-63
W-9.0460N 4 Ea0a
~ .-
~
NM ZNNMN*
0_
00
eo
44.:a
NNN
op-Ea~a-oana aA-
a,~~a~~tt aE
tit-
N4mON0flJNNNN Mm
rOZ00999 0N
PMM m. 0
ME~ % 0' , i
94
'a 0'.N
4 N 4 4 N 9 N
- c o o
4a .
p.
a.
2 4;201
4004 9 9 9~00 0 N000N.:coca I00 E000:0
2 xzn
1100 's j 0 ai
99 4M N..0 ON4P.4.*I09 0.
~~
000
00 00
00 00
N000
00
00 ,N 0 . n
0N
00
00
00
00000010aN .. .4p
0 N .. .p.P~ . 000
M 9 9 . .*... '9E .
M .4yrZNN .. 1 . .
NNNN 0,0
0009999994 o o o o o o o
4064404 0
4 e O
i 1
MP 2 -
Ma
IPN 0
.NNN nt N N
NzN
000 1.
4 z.
It
40
00000999
U U 4 N P.
U....I~l4
e4
WW~TIM-.1-,.4
NMU.942
NN9
-9...MEU14
44'a
ata
-..
ON 0 N. . N l ..
040 . ..
14. 1 9 :
N1
C C r0= M N**44P
V 40 M af.N 4
04,0
44 M EN 0 N.NE
09,6 -N
NS44 N co
4o
10 NggF.
m. C NN99.o9.
4
.4. .o.
0N4N ..
XN
.4
C O0ON~4 4~004
A,0
N-
N0o0
eN0~v Na0o
!2_22_22
.-
-04 -1 -44-N-.-04-9-99 9 O
N~J~4N N 0P5N'N~JNJNPl
~~IN
~130h
AEDC-TR-78-63
to
~~
~
00 % 0a 00 000
0Z
~
0
:4
N o
*v-
0.aP IO N -N 0 0 00 0000O
0 W% 00 0 0
00000oo00 a-
'
2 w
so aW ao-
ma 0 o
00000000004
000s00 i 440004,000444
~~~ 7 ~ 0
0 1 ~. ~ ~ - l~~ P 0
o i
M 7 o 00 .
N I0 t
. -
0-i
- .2
0loo P.Ol7f o-
00
02 t'4. --
~~~~~~~~~~~~ tN
= rP. i0
0Z~zI~ooo~roo~ooaoooooo
ep
.0 I..
*I
~~ 10.M MP 01
N.0N 0itNO
o"".0
coo 0
-0
0.
NN
IL
2 a
.ii.
------
e0---------10, CP.71 .0 . N
0 . 0 OOo4l a----------------------0OC~c*coo
iM 4!
0 0tI4.
00.(PM~*&
O&M0 C.
IDi~fPOO. C 0N1O -
nkPI
0 100 #i
P.p.
*7.OitPV?
11liPNNNNNNO C n N Ii 4010 m I I
NNNNNN
NNNN
NNN
INNNO NNN
NNNNNNON N* N m1 04
~ .0 ~
W N0 -
UO M
~ a
*0N
IWOU
*~~~~~~~
~~~~~
A0PN070oO.0NarS itP v
hio * a US~~~
. 4 -4 C0F C Nt m 4vPvM00-P.i.4.PPN *0,N 0f~ ' O 0M a 0" D,0 0 v kn 0: 00 UN 01 0 Inknin 01 01 U1 4401 .0i .0. ...t. N-N-P1 -PMM.NN-.mm mmt
hi.J
VS -
~NI ~
Z0 . N
. 4t0
..
P . 0
..
....
vz~~
.i..0.0.. it
a.....
..
0... P-C
~N
fn
( N.00A a, Cc do0,4N4
NN.007. 0 0 P_1N~00 N W
CL
D . ' .-.. O* N W PC 0
. . . QWDCD
C0
. . . WV W -----
. W
. . C N
CD
. m N
*----- ------
NNNNN
.NN..
NN
NN
NNN
NNNm
vi
0101
AEOC-TR-78-63
orumfu% MMMon"I
0..N4CP.0 MAAM4*T*:Y
-10 M 70 .. M4C0.-M
JM CC.
IN:: .
-NC'.
CU, 0
UUCC 0,0000000...
CN
-CO
oz44..M '
M 0, N 0;U'&Cmlciom~ N4COM440OJ
omm C0
NN-m-m44U
Nm M0O,
0 CU' 4M4U
COO
P. 0
a_e.C
... O @4r
'.
_
0
0F
0 coMo 00
4~
0
..
0
n.
0
N
0 0 0 0
MC
0
M
0
400
M000000000
00
000 60
000
00
00
0
-4
0t .2%
4
0 10 -
MOM
NNNNN
NNNONNNN
NNNNN-
41No0 0 44 00000
&44MON44-'6
4 4
MO-MCN-.000
4 4 4
~00
0000
1 0
0000 .0 M~~~~~
0000000 MM
MMM
4-,00 MMM-M M
0 m0000 MNNNNNN
.0
VlAP4O0O 0 000 N
cm -
N N N1 N N N N N! N N:!
N0N000 , NNNn
0 NNN-W O
--------------V A, ,04N'
--------
--------------1
0 N
omaP44N
zM-
2.
C' 0L
4mNN M0 44
MMMM
NNNNN 04 MM
NNN O NNNNN NC
NN -6 0M -
N 6
00
N NNN N U'4M.0v m
C 0.00000000
N -0046.00.464N.. 00 - 0 - 00 0 -0
I ON
NbN N
N V
N,
NNNON
iN NNN
N NyNv
ON-N NN
NN
N 0NN N M 0 1-
N N NN U..OZ
NNNNNN NN NN N NN
NCmw Ncco'U'U'4 0 0 0
.... 0 0
44 MMMNN 0
.N 40
N...OOC
0 00 0 .............. N .
O0 0.6
0 0 Co 6 C
- ,0 040MC0m
-0U 0
~~c 0 0
-
0 0~ It
~
.
0
000m 4
C N
,.N4T00...M4.0NQO
040 0
O~x
0,
06CC.-M
NN'..NN.~t NPC
0:06
U4-N0
C
'44
066.-NCUC4MP--0
4 4*40 44 440.
4 4
6N MMM
046 N6 N4040NN
4N
4-
'.
40N
60
.0
40MCI0
-00
000M .
.0U A
4' 4
N`
o "'NC C6 0N
CP m NN
't v~J9 m m NN 4 M
8 Q
00& N6N770U'0 L
,0
wm
0z *..004600N --
% -'-'-
'
4 N'M
c L..:.0 LO66NNN 0 m N 0
- -
- - 00..-6-N.
N,04
C44 .
OC:TN 4 C 60 CP444440
0 NC UCCC z M .00j,0
..
UUU4 4CC
U40 0
CC
O4C
N
wa
44
-0000
4
44
40 NN
w.NM4U'4N600
NN
NN
apt
132
AFDC-TR-78-63
af
in aN
VOMM.C0Nvoci
*..$
~ *tP-P.44O.Mv
00 --.A ---
IP 2
IP
*N4Cif0M
O~
-00
-co
------M. 404t
OM
---.
w
N~CC9I-04O
-o
ccP~
Am~0
*~ *JINQ41NO4MO
00000000 ---
z40-04; -
--. ------
84
0oaC&
M200ONNN
CONNNNNNA*N
O0
0
00
-14
---0'
`-0N401M 0.0000000ft'Oa
aMIiiNN ell 00 0 000 ,cce eo 0 0
4-4N.C0,41 *0 MCe
%.cv NN 0 N0 0 0 2.eO 00 0 0
---------0o m
NN 0 0 0
-00
~
1 N0 0 0 O
0
00 0 O
00
0 0 CO
0 M
0 0 .
U'
-0 .0 0 N O N 0 O 0
-0 &'
0 0 0 0 OC
-0
4 U
:1~-00 4~
0 000 0 C 0
444#4 0~t
4
0 0
41
-
0-
00 .OC
0 C
OCN
CCIE WOO
M00-4a O
mc c a .......
.0o
x CCM 10 0 .i0ttt--. .4
QN
1 -M~
~N
U '
.. . W
t~f N
4 -
N.CCC~ 0 '
0 0
0s0 0
TMOM4604CC
M42 Z 4-------
12! -
4.4f0 ---------
0s.l -
N v fv s-CC~~ C-
N 161CC
0C~C~l4
0000.
00Cft,0
44141,f4' 3,00 IN
000
C62'#
0 c
4 !
coo
-0
244
~133
AEDC-TR-78-63
--
--
ooooo
ooe
O gs o ee
O e
0 e-i 0
a-
00000 N
000
~
0f 0
O oooo0ooe j~ gi~ eg
~ e
M Miiil 0 0
jk
cmfu
00
mN
I 0 0
N 0
N 0 N 0 N 0
0 N1NI.J 0 0 0 0 0 0 0 0 0 0 0
*~~~~ 0
-0
U~
0 A0 0 ~~ ~~~~~ Zin
4 . 4
C) A 4N O N
. 0 N
6 A NW , ' A -. 0 0 4 NN -
On O 4N -
0N ' 0 .
ID
o44 InS NN
oj
AA 0 W
IZ *I*Ol0
zz 5 NJ M
j 0
4~ S 44j9000400*01N. MNNNN N
54 N
N c
*0 .. e
..
. . _-M0
. .
4-Itti
.0004 V :V 404 Z
.1 4
1i 0
Ow
0 4~.
~
0
N5.
~P4ONIO. ~ 5A ~O I ~
4
N OS N.--
4M %A 0
_j
-- N
l -O
5El N
4.0 D
-e4,ul -
IA, U
LA M
44C0 00 M
t05aNA
cc
a0-
l MO 0'4 4410NW
N2( 4
0l c,
M.l0M
C 0
0 0.
ry
rue
a---
N 515 IA 2 W z 0 2 0 4 .4 A4A5-45 0. M 0 0a-N 0 4 00 000-404410 0 0 0 4 4% . CC:4J4NItN0 .'t . It 0.Ct4 ItN 1cj, 0000000000 00 'N 0 0 00 0 6N 000000
0
In 05-4.
40n ---
0 0 9.
O 000 0-0 00
o 00 IA1 0
co 0 v 00 f;;
z
An
M5
to
0 1 0c00000000Ca 0~ 0000000000.
:00-Nmoom'A44 ---0000000080
o 0000000000
c In
2
*
W 0'SN-400I ~ ~ aN 0 4 ~ ~ a N 444s0
'. )-
n,.V0 1 N
000 0
;!.DNO:1: o0-.05-2
4-VI44..e~ 0 000 00. 000
40441 151~51 00 00
a,
IA~-, t 00 I 0 %UI 00
_V0
0 000
000l
0000 0000
000 2
I-.0 000
15u W0o
"'
j 0 2 N J'0.0N5-40
c, 0 W 47 N 00000 * *
5 . . 04 C5-N4I M1
,IN!!N4
4 1 04D700
.A551000-c..-'NmNNN
c M,.05 0
N4 LA If 2
~0
o. a, InI
5-A.55-045
4USC
5-25-0
4.I4-04
0045C S
VM1'5C4500...15
In-ONN M 15
W~M0 I a
u2 0. r U.
00-45-1541 440505-
It mN0- M 0.M454.2!M
.E .002 0 4 ZZ.0740N.205
-55.a, 504544N
0y0 0.
4
In Nt o 00JA Oo04I15 '0 2-51 =2.CA~f 41 55.. :,l In040.5s5 -1 M 4v 0090 *I05-0 N - 0000000000 -!557-0 = o 005.M g1 10-5 Coco5 MI w 41414=141 0000
_. W 2
-
.5
0,00 c
oo 0 0 0
ct5
0
-2
U'-0.155IA50155 N0 a0 .r004.4N-50 ( i 0 0 0 0
~~~~0 0 0 *~ 4 IA
-
02 VA,JCID M 1, W 000000 n7 0n 00 I 0 00
a.a 2
2i 4444451W 0
N
'
504I075 le a 0000#0N-..C
50
4 .
. .04N0~0 .O-.1-
7 i-.4. 0 054154M4
N0
04 Z oLfut0
n 0Nm; 0
4~~~ 2n
0 4l
~ 5-450E4 ~ .4
M.
~
G0 0--
~
5 2
0 Z.;
0 7 0 70-N415. 0
041oOIAOIA 0--
7N
20"M vi.00
134
AEDC-TR-78-63
777.E45J . u;s s .
7.770P. . . . I e0 N5NnNI' . . .
s s s m
A4~~I0U'S-E4~ 00
'990-
N'aCNONm
~-~~~UI%
*eM0-60 0000000000 .4'a~Mfl 0P 10~ M0
-tV55
*0-Qof5 0000000000m
000000
W, N O 0 N 5 4 N~~6t .40
5oe m.f P
N V
Cf
.A'sI5 . .* . .0 .
. . .. M
MN .
N~N
~
F-
NMM
MN~ONIMJ54-----
CJo
NN
00 -----
-O00
*m
00
5-
-----
----0
NNlyNNNNA
cecooc0 000000 00000
N*o-p Nmmm
z7NP 440
NNNN
NNN
0 ~ONq4 40 0
NNN
c
N N'1mCE I Op 0 N'555", I ,,m meNN6 06c60606 .. wo NN wNmNNN NN N NN N NN NN 0 0 0 0 000 0000000 Ooeeoeo 0 0a 0 0:00000 00 0 0 000000000 0 OCCOCCOOOD M. ;m P-5EJ..p 70 C . 0000N00 c~ .
0N.66 0 M:: N NN NN 0 0 0 0
0 0
0 0 0000
00 0 0
4W 0. m 0-N , NN -n!:00 0
~E
4- --0000000000 00 00000 00 0
fUN
50.~4.-*#~0
E0600 N,44 -
101
0 55.N0SoOO sO -
~501
0 0 00 0 ON-5 o00..... G04004-0 .0
ll..4E5-EMm' --
IM5M --
700 0 0 0 0
0 0
0 0
0 0
00 0 0
000 0
I. 4 m
M04:6o 6M~vw!0m
104VN00 NO..06
"0 50
owwN-61
NP-
N.
*O *600NC *55
O 00 omON0 t. 0.4
00
00 0
00 0 0 0
000000 00 00.J05'0 o0 0 0
000 !0
*0.153'00J5o 04
M5N0!60-.
0f~ SNoom4
f~~r.,ooo~~~
oo
oo
Mooncco. c m-,,
Oeo4o
055N~
PECN-
omS
ENam-,#"O
N6040NO
060
N5CN-..-NP...-
45500
-5544
om0
d 0N50.s
.f0,01-5*5 Do
-00-f0 004N-0.00wg5cm
0 055
l00550
135
AEDC-TR-78-63
000 00
00
000
000
00
00000
000
0000000000
00000000a
000000C0000 N-t
00
000 I 0 t
0000000 .
.N
0:
A
...
&&azwQ
-OZOO,~f44
*MNN
N7777
00000000 00000000
.q -O4 4 00N
00000000
NNISJnl * 0-'4-.CN
00 00~
.~0:.j NO*N m 0 000 V N .AIQ
000000
"00
100000 N-c0 4
0V-
100
~
4N
N4 4m-.P.
-4
0 w w
UN.004~ pNJ,
MI
I-
-~.
N ItlA ...
Vt 04 M0.4.4gNo
co'.
N 00 0 1000 00 0
~C~ 0 01
00
a .M
M004NOCmc-.N
, 00
:1.
00
00 00000 %f0.-. .
00M 00000
.0
4..0!-tNt-OcUI 00000=0000
0n I
444
.NO
44404440
0000000
0000000000
00*NO*040
..
SOAN4.
NJ
0
0 0
140..O
004
12O
0004.. 01-0,0N IC M N
M-00
0
000
~04
00
00
-- MO.N0 000000000
-S., 0
0000000000ca
000
0000000000
M
O NN0 In'0 0 m~0 2
NM-.'0
;;:t~S.
1~
M CY
NN
w
-.-
n MN
004
N0
40.
0M-
'a
NC;OM0S~.2~
NMMMIM
00 0 00 0 00 00 0 0
'WMMn
0 00 0 0
NNN
00 00 0
N
0 0
NNNNN
0 00 0 0000000000I C : t 00 !
.NNN .
0 00 00 ! iI 0 0 00
N..
00 0
-------N N
00 004'0 00 00
NN 10
NNN
NO.00
MM 001S.
MNC-C
N4NNC..NNN0N0 0 0
NM-0
0 0
0000000000
1i J'
P.00M4IMCM0~
~ ~ 4N0S.OMNO ~
000 000
O'. .. 1
OOOQ4Jl4
~ ~~~~c -00 4.
4C. 000 0 C ..
.MJ.0
4NMNO: .N4
00M~4
. 1
-.
N0. -O.
- O 4'44..0 .
0000000
0C40.~M..00 OtN0~0004 4 MN 0
0000 00 0000000000
0000 C.. 4 . N
0000 0 M
00
000 -41U'0
600
00.0000000 , C 00 P
0NMN40CU~~~~~...M 44IlNNS.c0MC
4M 4 S. N4 0N C~C~t 4N.-0~j,0~.~ MMNN .C.-. .M N 141CO -0N S~N 0N N . MN.-0N4~.Ii S.NM.-CM. I NI 074 MlMeS..0 N 4.0N 40..S.-, S404..
0JN1 0 M -..
I-.".o...........000...........
40N00
.OOM
NN M
44 INNN
w0
Ir. N M M .........
a4N 00
DwNo
0o
sIC0C'-T~~~~~0O -NNNNNNNNNNmm -
4&0CS-... 4000...
NN
.........
~mm
NNoNNom~4
z4
N N
NN
40
dc1 -..-.-..-.-..-.-..-.
. I -. mMII 0
t0 0
:~PS.4 C 400, a 4 S-
S. LI I N Cp0...am. M u . 0 M f
. N .
%%.O .
NI' .
S.%0 . .
"1 114G.
I 1. fl : v Ct
=-N
0.
..
tN
toa
136
AEDC-TR-78-63
0000000000000000000
0000000000
000000000O0
00000
oooooooo
N!0UN
oeoacaaa
:0N2:~d~4
MON
ON:44.g
oooaoo
oooo4o
0O.4 ;!4O:4NM LA
M.4 0
ooooo
4, 0td! NY
. 0-
0
a m4
lo
44 044 e P4N.NE
A.4Na
mm
Inz
.Nlo v--,m
* .0
-N..N.NNNNNN I *# 1 ,6 a i
NvNNNNNNNv ut,,i~
---su1g
000000000
00 0 00 0
0c0oo20
00 00 00 0
c2
00
0
00
00
0
0
00
000000a0o
00 0 00
02
00
00000
0 00 00 0
MI -MOOM
**4 0 4)
mcU* D 0 Dit f0
ft'.A
~~
ic0.4tP44...-0M0o40 M *40u NN N
4-oom
09 f9S 55
9 O
9 P042
19 O
V o4.
9 4 z0 M~.
;;
IS0C c
SSS am~
ry
af-m s I
, , g g It .OO~O .O C
1O
NIg I S S I S~ g "tItI 1 OC %. O
g s
e , , OOO
s , ... O3w:~O
i ,
, O
g I
i O0C4
0
4
51 44 % hi N5N59445 w m
C- 15.10 N0 nu
CN O -
.-..-
s15009...4
.M 9.
N
N0UomNU
oNE
00 *N -4-0
I1 t0 ;4 4 5955N09.50
""
O .C
0..P.0 59
0N
-.
~P.*
*0
40
4.
50
00
59
013754
AEDC-TR-78-63
too 42; me o
m *SN.ANu0. I
0 nN P ..
4IL~4t
NM*700o0o
- coo
W000000000
00o
NNM.I-.e..4138
AEDC.TR-78-3
NOMENCLATURE A AC A* Area Exit area, inviscid contour Sonic area Sonic speed Factor in logarithmic skin friction law, Eq. (77) Coefficients, Eq. (35)
a* C
C1 , 2 ,3,4,5,
*
C.
R were
1Cf
C p D 1,2,3,4,5,
6
Cf
incompressible
Fn
Multiplying factors, Eq, (97) Ratio, R /R i c Multiplying factors, Eqs. (94) and (96)
FR Gn
139
APEDC-TR-78.63
Ratio, h K In log
5*/1
Heat-transfer coefficient Streamline curvature Natural logarithm (base e) Common logarithm (Base 10)
M rm N n
Mach number
Exponent in Eq. (90)
P1,2
Pn
in
Coefficient of 6 at nth point on contour Prandtl number Factor related to heat transfer, Coefficient in momentum equation Eq. (91)
Pr
Q
Qn
Velocity along streamline or, in boundarylayer equations, velocity within boundary layer
qe
140
-.
~-.-
T --
..
. ....
. ..
AEOC-TR-78-63
Rg
,R
R1
R0
Re
i Sr r
Distance from source
Distance from source where M 1, used to non-dimensionalize distances for inviscid calculations
rV
I
,9
s,t,u T T
Cubic integration increments, Appendix B Temperature within boundary layer Adiabatic wall temperature Reference temperature, Eq. (87) Free-stream temperature at edge of inviscid contour 141
aw
Tc T
I..
........ . ,
AEDC-TR-78-63
T T
wD TWT u Wall temperature at nozzle throat Axial component of velocity, normalized by a* v Normal component of velocity, normalized by a* W X x Velocity along streamline, normalized by a* Ratio in Eq. (36) or (38)
Axial distance, normalized by yo in transonic equations, normalized by rl in .inviscid calculations, not normalized in boundarylayer calculations
y Y y
Normal distance, normalized same as x Throat half height, used to normalize x and y in transonic calculations
y*
infinite
142
: ;:: ,,F L i ,= . ! '. : & :|!J . W- s~l ';.4.,.2iL... ,..... ;,.. ,:L..,.,;,',al,.1;,.;,,,ili, ;',' '"" ::' ' i . ., ' . , .,, I,,; ,
AEDC-TR-7843
Mean angle of right-running characteristic, or factor in temperature distribution in boundary layer Mean angle of left-running characteristic
S~Boundary-layer
thickness
6* a*
6t
e
6
,Momentum
c c 01
0k
143
-I-
-..
AEDC-TR-78-63
0 n
K
Constant in
I2
(L y + I)S Mach angle, sin PC I PW e (1/M)
Viscosity at value of T c Viscosity at value of T Viscosity at value of Tw e Viscosity at value of T Distance along right-running characteristic Wake variable in logarithmic skin-friction law
Density within boundary layer Density at edge of boundary layer Zero for planar flow, 1 for axisymmetric flow
Pe a
4)
w
,Flow
angle
Flow angle of viscid contour Prandtl-Meyer angle
'44
AE2DC-T R-78.83
SUBSCRIPTS
1 Values at point 1 on right-running characteristic Values at point 2 ona left-running characteristic 3 A,B,C,DOE, F,G,I,J,T Values at intersection of characteristics Variables evaluated at points on Figs. 1 through 4F
C(Y)
Coefficient of third-degree term if throat contour is a cubic Coefficient of third-degree term if integrated throat contour is a cubic
C(YI)
145
....
AEDC-Tn-78-63
C(YP)
D2A/DX2
Second derivative of boundary-layer correction evaluated at the throat Second derivative of corrected contour evaluated at the throat
D2R/DX2
(66)
HYP/YO
Value of hyperbola with same throat curvature ratio 10' [C(YI) - C(Y)] for Point 2
IcY
INT.Y/YO
146
. ............
AEOC-TR-7843
PAR/TO
PE/PO
R(IU)
RASS
RTHI
SKPP
SMPPP?
source flow
THETA
0 from Eq.
(62)
used in Eq.
(61)
WE
3)
W1
3)
WO
SWoFP
147
ANDC-TR.7-43
WRPPP
w0
fV
148