Nozzle Design Theory

Download as doc, pdf, or txt
Download as doc, pdf, or txt
You are on page 1of 7

THE DESIGN OF A SUPERSONIC NOZZLE BY THE METHOD OF

CHARACTERISTICS
Read Anderson Section 13.3, especially the example discussed in detail on page
591. A computer code (Matlab) Is avalable at our web site for the design of the
straightening section of a 2-D nozzle, and execute the program for two test section Mach
numbers. The flow in the test section must be uniform, parallel to the x- axis, and free
from any compression or expansion waves. A reasonable contour for the subsonic section
is assumed.
The shape of the expansion section is arbitrary (see p. 591 and the enclosed
figure). For simplicity here, take all of the expansion at a sharp corner just downstream
of the nozzle throat as illustrated in the sketch below, i.e. set the radius of curvature of the
circular arc R , in the figure to zero.
The Mach number at the throat is unity. It may be shown that, for a given Mach
number in the test section, the shortest nozzle occurs when the expansion section opens
up to an angle equal to one half of the Prandtl-Meyer function corresponding to the test
section Mach number. The shortest possible nozzle occurs for this condition when all the
expansion occurs just downstream of the throat. Thus, the calculation required here
corresponds to the design of a shortest possible length nozzle. Because of boundary layer
problems, wind tunnel nozzles are normally not designed with the full expansion at the
throat.
The nozzles are to be designed for the following test section Mach numbers: 2.0
and 4.0.
In the figure, the expansion region immediately downstream of the throat is
represented by 3 steps. The program must allow for any number of expansion steps, N (N
< 10) . Take 6 to 10 steps through the expansion turning angle (the more you take, the
more accurate the solution).
How does the program work?
1. Numbering System: Numbering the nodes using a two-subscript system shown in the
enclosed figure simplifies the programming. It is convenient to number the "wall" points
on the straightening section using a single subscript system. If you find a more
convenient way of numbering the nodes, you may certainly do so.
2. Expansion at a sharp corner:
(1,0)
(2,0)
(N,0)
R = 0
Referring to the figure above, the radius R of the expansion section is taken to be
zero. Thus points (0,0) , (1,0), (2,0) , ...and (N,0) collapse to a single point. The flow
slope at these points steadily increases until it reaches
TS
/ 2. That is, at these points,

i,0
= (i/N) *
TS
/ 2
i = 1, 2, ... N
where N is the total number of steps.
Using the theory behind Prandtl-Meyer expansion fans covered in AE 3004, or
using MOC, it is easy to show that at these corner points,


i,0
=


i,0
i = 1 , 2 , ... N
The x and y values at this corner may be non-dimensionalized based on the throat
semi-height. Thus,
x
i,0
= 0
y
i,0
= 1
i = 1 , 2 , ... N
3. Interior Point Calculations: At the interior points (i,j) in the figure below, the flow
angle , the Prandtl-Meyer function , the Mach number and the (x,y) coordinates may
all be found, once information at the points (i-1,j) and (i,j-1) are known, using the interior
point process. Thus, these points may be computed in a "double DO loop" in a recursive
fashion.
(i,j)
(i,j-1)
(i-1,j)
i) Compute first the flow slope and at (i,j) from the compatibility relations:

i,j
= [
i,j-1
+
i,j-1
+
i-1,j
-
i-1,j
] / 2

i,j
= [
i,j-1
+
i,j-1
-
i-1,j
+
i-1,j
] / 2
ii) Compute from
i,j
the Mach number and Mach angle at (i,j).
iii) Compute the (x,y) locations of the node (i,j) as the intersection of
characteristics originating at nodes (i,j-1) and (i-1,j).
dy/dx = ( y
i,j
- y
i,j-1
) / (x
i,j
- x
i,j-1
) = tan [ (
i,j-1
-
i,j-1
+
i,j
-
i,j
) / 2 ]
dy/dx = ( y
i,j
- y
i-1,j
) / (x
i,i
- x
i-1,j
) = tan [ (
i-1,j
+
i-1,j
+
i,j
+
i,j
) / 2 ]
4. Points on the axis: These points have a subscript (i,i) where i varies between 1 and N.
require a slight variation of the interior point process, and are best handled in the same
"Double DO loop" as the rest of the interior points. At these points, the flow angle is
zero. consider a characteristic that connects a node (i,i-1) and (i,i). All the information at
the node (i,i-1) is known.

(i,i-1)
(i,i)
Nozzle Axis

i,i-1
+
i,i-1
=
i,i
+
i,i
=
i,i
Once the Prandtl Meyer function values
i,i
at these axis points are known, the
Mach number at these locations, and the corresponding Mach angles may be found. The
(x,y) locations of these axial points are finally found by solving:
y
i,i
= 0
dy/dx = ( y
i,i
- y
i,i-1
) / (x
i,i
- x
i,i-1
) = tan [ (
i,i-1
-
i,i-1
+
i,i
-
i,i
) / 2 ]
5. Straightening Section: The straightening section is made of wall points sequentially
numbered from 1 to N+1 as shown below.
1
2
3 N+1
(N,2)
(N,N)
(N,1)
(N,0)
The first wall point #1 is the easiest to handle, because it lies at the same physical
location as the expansion corner. The properties , , , M at this wall point, and the
physical location of this first wall point are identical to that of the node (N,0), previously
computed under the section on expansion corner.
The second and subsequent wall points are handled as follows. For the sake of
illustration, consider wall point #2 shown above. A characteristic line originating at the
interior node (N,1) intersects the nozzle wall at this point. This characteristic will reflect
back as a compression wave if the slope of the wall at wall point #2 is smaller than the
flow slope at the node (N,1). If the slope of the wall is larger than flow slope at (N,1)
then a new expansion fan will be generated at the wall point. Therefore, to ensure no
expansion waves or compression waves are generated at the wall point, we require
Wall slope at wall point #2 =
N,1
Then, using compatibility relations, we can conclude that the Prandtl Meyer angle
, the corresponding Mach number M and the Mach angle are identical to
corresponding values at the node (N,1).
The (x,y) locations of the wall point #2 may finally be found as the intersection of
(i) a characteristic line drawn from the interior node (N,1) of slope tan( + ) computed
at (N,1), and (ii) streamline of slope tan[(
Wall Point 1
+
Wall Point 2
) / 2 ] , drawn from
wall point #1.
The wall points 3 , 4.... N+1 are similarly handled. The entire nozzle wall point
calculation may be automated in a DO-loop.
6. Computation of Mach number M for a given value of : Our nozzle design problem
frequently requires finding given M , and finding M for a given . The former is easy,
and requires programming equation (9.42) in the text. Finding M that corresponds to a
given , however, requires an iterative process because of the highly non-linear,
complicated nature of equation (9.42). The following iterative process, known as
interval halving procedure is recommended. It takes advantage of the fact that is a
monotonically increasing function of M.
a) Assume that the value of M for the given lies in an interval between two
extremes M
LEFT
and M
RIGHT
.
Good starting guesses are: M
LEFT
= 1 and M
RIGHT
= 10.
b) Find the mid-point of this interval, M
MID
= (M
LEFT
+ M
RIGHT
) / 2
c) Check the width of the interval M
RIGHT
- M
LEFT..
. If this interval is
less than some tolerance, say 0.001, then the Mach number you are seeking is M
MID
.
Exit the iteration loop.
d) Compare given with (M
MID
) .
If is greater than (M
MID
) then you need to search in the right portion
of the current interval. Therefore, set M
LEFT
to M
MID
.
If is less than (M
MID
) then you need to search in the left half of your
current interval. Therefore, set M
RIGHT
to M
MID
.
e) Go to step (b). Repeat until the criterion in step (c) is satisfied.

You might also like