The shock standoff distance ahead of a general rounded stagnation point has been estimated under the assumption of a constant-density-shock layer. It is found that, with the exception of almost-two-dimensional bodies with very strong shock waves, the present theoretical calculations and the experimental data of Zakkay and Visich for toroids are well represented by the relation Delta-3D/R(s) = ((Delta-ax sym)/(R(s))/(2/(K+1))) where Delta is the shock standoff distance, R(s),x is the smaller principal shock radius, and K is the ratio of the smaller to the larger of the principal shock radii

Reshotko, Eli

English

NASA Technical Reports Server

lo
o
n
I--
<
<
Z
NASA TN D-z050
, - / /
// • y /
TECHNICAL
D-1050
NOTE
ESTIMATE OF SHOCK STANDOFF DISTANCE AHEAD OF A
GENERAL STAGNATION POINT
By Eli Reshotko
Lewis Research Center
Cleveland, Ohio
NATIONAL AERONAUTICS AND SPACE ADMINISTRATION
WASHINGTON August 1961
https://ntrs.nasa.gov/search.jsp?R=19980227349 2020-06-15T23:05:24+00:00Z

8
NATIONAL AERONAUTICS AND SPACE ADMINISTRATION
TECHNICAL NOTE D-lOS0
ESTIMATE OF SHOCK STANDOFF DISTANCE AHEAD OF A
GENERAL STAGNATION POINT
By Eli Reshotko
SUMMARY
The shock standoff distance ahead of a general rounded stagnation
point has been estimated under the assumption of a constant-density
shock layer. It is found that_ with the exception of almost-two-
dimensional bodies with very strong shock waves, the present theoretical
calculations and the experimental data of Zakkay and Visich for toroids
are well represented by the relation
where _ is the shock standoff distance_ Rs_x is the smaller principal
shock radius, and K is the ratio of the smaller to the larger of the
principal shock radii.
INTRODUCTION
In recent years much attention has been given to the problem of the
inviscid flow about blunt shapes, particularly about bodies of revolu-
tion and cylinders. The theories proposed range from those that con-
sider the shock layer to be of constant density to exact numerical inte-
grations of the compressible flow equations. (A summary of these tech-
niques and an extensive list of references are given in ref. i.) The
inviscid flow about general blunt shapes (finite bodies with unequal
principal curvatures) has_ however_ received very little attention.
Hayes (ref. 2) has derived an expression for the shock standoff distance
ahead of a general stagnation point for a constant-density shock layer.
This expression includes centrifugal corrections to the pressure and
velocity distributions. It was, however_ not evaluated in reference 2,
perhaps because the centrifugal effects depended on shock shape away
from the stagnation point and a representative three-dimensional shock
or body shape is difficult to choose.
The present analysis retreats from that of Hayes (ref. 2) in that
centrifugal effects are neglected. In this lespect it is the three-
dimensional analog of an earlier analysis by Hayes (ref. 3). The shock
location is a function only of the density rstio across the shock and
the ratio of the principal radii of curvatur6 of the shock. It is inde-
pendent of whether the body is, for example, an ellipsoid or a toroid.
Because of the manyassumptions and similifications that are made
in the present analysis, it is not reasonable to expect the theory to
yield precise absolute values for the shock standoff distance ahead of a
general stagnation point. It will be shown, however_to give reasonably
good results for the ratio of this standoff distance to that for an
axially symmetric body. Then_ this ratio together with exact solutions
for axially symmetric bodies such as those o_ Van Dyke and Gordon (ref.
4) maygive reasonably good results for genelal bodies.
The main analysis is preceded by a very approximate calculation(OVERSIMPLIFIEDMETHOD)which_ in spite of ils crudity, yields a result
that is in general agreementwith the main aralysis. The present re-
sults are also comparedwith the experimental data of Zakkay and Visich
(ref. 5).
b_
!
GO
C
Cx, Cz
K
k
P
R
Rx, R z
U_ V, W
X
SYMBOLS
shock-layer velocity gradient, eq. (_)
shock-layer velocity gradients in principal directions
ratio of smaller to larger principal radius of shock wave,
Rs_x/Rs, z
ratio of free-stream density to that behind a normal shock_
P /Pz
free-stream Maeh number
static pressure
radius of curvature
principal radii of curvature
velocities in x_ y, and z directions, respectively
3x coordinate in shock layer identified with smaller principal
radius of curvature of shock
y direction normal to body
z coordinate in shock layer identified with larger principal ra-
dius of curvature of shock
r ratio of specific heats
shock standoff distance
dummy variable in shock surface in z direction
dummy variable in shock surface in x direction
p density
Subscripts:
ax sym
b
S
i
2D
3D
oo
axially symmetric
body
shock
conditions in shock layer
two-dimensional
three-dimensional
free-stream conditions
OVERSIMPLIFIED METHOD
In this approximation the constant-density shock layer is of uni-
fom thickness in the neighborhood of the stagnation point; which is the
region under consideration. The free-stream Mach number is very large;
the shock layer is very thin A << i; -- _ , and the pressure distri-
Rs
bution is assumed to be Newtonian and constant along a normal to the
body. The velocity in the shock layer is also assumed constant along a
normal to the body and equal to the velocity adjacent to the body.
Within these assumptions the shock-layer velocity maybe expressed
Ul _Cx. For Rb _ M_]-- _ I and _ - i) _ _, the velocity gradient isRs
% .y_--- i P_ ]" - i i + k
• =- _'_-- or y _- __k_ soC = _ss But in this limit k Pl Y + i i
the shock-layer velocity gradient can be wrJtten
C =_ss l+k
Three cases are now considered:
Case I: Axially Symmetric Flow
Shock
From continuity considerations (see sketch),
P_u_X2 = PlUl 2_x fiax sym
Using equation (i), the shock-layer thickness becomes
kRs _+k_x sym _- 2 2k
Case II: Two-Dimensional Flow
The continuity equation per unit span i3
p_u_x = PlUl A2D
(1)
, 5
o
I
so that the shock-layer thickness becomes
A2D _kRs _+
k
2k
which is just twice that for the axially symmetric case.
Case III: Three-Dimensional Flow
P_ lj
u
Rs_ z
Body surface
lhock surface
Let Rs_x be the smaller of the principal shock radii and let
K = Rs,x/Rs_ z. For this case, the continuity equation in each quadrant
can be written
p u xz = _lUl z ASD + 01WlX _3D
or
In terms of
p U
_3D = DI(C x + Cz )
__ + kk 2k
1 1
+--
Rs_x Rs_z
Rs_x and K_
_+kkRs, x 2k
A3D= i +K
6Note that in all cases the standoff di_tances are of the form
(A/Rs,x) _ _/_, which is not appropriate to rounded bodies. This incor-
rect form results from assuming that at a g_yen x and z the veloc-
ities in the interior of the shock layer ar_ equal to those at the body
surface_ so that in effect all the fluid in the shock layer entered
through the normal portion of the shock waw as did that fluid at the
body surface. Nevertheless, the ratio of the preceding results for a
given density ratio k_ namely_
(Rs_xlSD 2
I+K
(z)
will turn out to well represent the results of the following analysis.
ANALYSIS
The analysis for the three-dimensional case will now be somewhat
improved. In particular, the variation of velocities u I and vI
across the shock layer will be taken into account.
Consider the flow in a quadrant of the constant-density shock layer
ahead of a three-dimensional stagnation point. The coordinate system is
considered locally Cartesian.
Shock # l
surfacei / \_
J
.... "_ d_
/
- 7
a7
oG
r{
!
The pressure di.stributiDn is taken to be Ne_dson_an. It ::,s _nsldored q
function c,f x ,_md _ cn]f. i:,h_tsindependent of the n nr.',_ ,c_ord£n:_!,£
y. In the nei__hbsrhoc,£ ,:" the stasnRtion point it m%' %< ',.a'[:t,<n-':r
b_< >> ], k << t
2/ _'_<° R 2 R 8 +
With the aiQ of Bernc,ul]5 's reNa-tion_ the velocities a% _ ;gc:[nt in the
shock layer can be estimated. C,tnsider the streamline that r,r_sses
through the shock at the point (!;_I}). Bernoulli's re]ation states:
1 + Pl 2 :: _+ Ol 2 ( )
Frcm the oblique shock relations ]n the limit M_ >> ]-_ k << i
£
+ v + 'k %,
12
loo
•4- • )
Upon substituting relations S) and (5) into equation (4),
i/x,z : { { ,_2
CO \&x 2- +S, Z
(i_)
For the portion of the shock layer being considered, the principal
shock radii are assumed constant. From equation (_)_ the pressure _rs-
dient in a given principal direction is then independent of 8he other
principal coordinate so that uI : Ul(X,{ ) and wI = Wl(Z_). Since,
in addition: the normal velocity v/u_o is of order k (th_s can be ver-
ified a posteriori)_ the shock-layer velocities in the x and z direc-
tions are_ respective!y_
S--_- + _k-_-- + _(k2) (7)
-_ = + _k + _(k2)
Uc_O Z
(8)
Ul x
- = _ _ _ -_ and
For the stagnation streamZine (<1 : { 0)_ u_ Rs_ x
Wl z
%]lese are the velocities o-)tained upon assuming as in
L½ Rs, z
the aforepresented OVERSIHPLIFIED METHOD th_Lt all the flow enters the
shock layer through tile normal portion of t]_e bow shock wave. Currently_
however, the shock-layer velocity uI at :: depends also on the height
[ ab which the streanO_ine in question traversed the curved bow shock
wa,e. The \,,elocity w! at z depends sim:'.]_arly on {.
_e shock displacement distance is now to be determined. From con-
bi_n:dty, the mass f!c,w into the shock-layer quadrant through face ABCD
2i1S13 ,_ _"
_q_t_,__the strutof the flows leaving thr_Jugh the top face BCEF and the
s [de ['ace ABFG.
, hocl
./
B_<t since the normal velocities through the top and side faces depend on
and _, respectively, the division of outflow between the top and side
faces must be determined to properly carry cut the mass balance. In other
wcrds_ the curve BHD must be found such tha_ the inflow through area BCDH
ex]ts the top face of the control volume while the inflow through ABHD
exits the side face of the control vol<une. The line BHD is the inter-
section of the stream surface through the dividing line BF with the shock
wave and is found by tracing streamlines back from line BF_ which has the
coordinates (x_z). Once this is done_ the shock displacement distance
can be _ritten either as
fz /z _Bm_(_)d_
: ---
_3D -- 0lxw1(t,z
(_,z)
(9a)
or
A3D = f x
k _Bm)( _)d_
_---
(_,x)
(gb)
The equation of a streamline is
dx dz
Ul(X,_) : w1(z,_) (loa)
Upon substitution of the velocities from equations (7) and (8) and
integrating between the prescribed limits
f fzRs, x dxV_2(1- 2k)+ 2_ 2 RSj Z dz_/_2(1 - 2k) + 2kz 2 (lOb)
the equation of line BHD in the form suitable for equation (9b) is found
to be
sinh[l - K)sinh-l_+ K sinh-l__]
where
From equations (9b) and (ii) the expression for standoff distance
is
d_
sinhEl- K)sinh-l_ + K sinh-l_] V_2 + _2x 2
(13)
i0
By letting X = t/x, equation (13) can be written
_d.X
K)sinh-l_ _ K sinh-i _]_X 2 + _2
(i_)
Equation (14) is exactly the expression that would be obtained by neg-
lecting centrifugal effects in equation (4.5.7) of Hayes' analysis
(ref. 2).
In general, equation (14) must be integrated numerically. However,
in the special cases of two-dimensional and _xially symmetric flow,
closed-form expressions can be obtained. Th(,se are:
K = O, two-dimensional flow:
b_
!
Co
kRs'x l dX kRs,_:: sinh_ I _ - 2k
A2I) =-V_-- 21<_._X 2 +_2 -"v'T- 2k 2k
(is)
K = i, axially symmetric flow:
l
kRs,x X dX = kRs,x (16)
These two special results were obtained by H_.yes (ref. 3).
Equations (15) and (16) are now compare( with exact solutions. In
figure l(a), it is seen that equation (16) a_irees rather closely with
the calculations of Van Dyke and Gordon (ref. 4) for k _ 0.i, while at
higher density ratios (k _ 0.3) it overestimates the exact solutions by
no more than i0 percent. Equation (iS) for _wo-dimensional flow is less
satisfactory in that it overestimates the results of Van Dyke (ref. 6),
Belotserkovskii (ref. 7), and Uchida and Yaslhara (ref. 8) by 20 to 50
percent. This poor agreement is reflected also in the ratio of two-
dimensional to axially symmetric standoff di_:tances as shown in figure
l(b).
ii
CO
t--
!
For the three-dimensional case (0 < K < i) with density ratio k
greater than zero, equation (i¢) was numerically integrated on a desk
calculator using Simpson's rule. The results are given in table I and
are also plotted in figure 2. The results shown for k = 0 (table l(b),
fig. 2(b)) were obtained by evaluating the integrand of equation (14)
for _ _ 0 and then integrating. The result is
_ 1 in (17)
s,x/k_<) i - K K
The curves of figure 2 are seen to be quite regular, which is per-
haps not very surprising. However, the resulting ratios of three-
dimensional to axially symmetric standoff distances given in table II
and plotted in figure S show a more interesting result, namely, that
most calculated points show surprisingly good agreement with the crudely
derived expression (2) from the OVERSIMPLIFIED METHOD. The exceptions
are for bodies approaching the two-dimensional, for example, for
K _ 0.2 at density ratios k _ 0.I.
Before proceeding to a comparison with experiment, it must be real-
ized that the present analysis yields no information regarding the vari-
ation of shock-layer thickness about the body and leaves the body shape
unspecified. It is therefore not suited to determining the shock stand-
off distance ahead of a given body. Considering some of the results for
a sphere (ref. 2), it is doubtful whether much is gained by considering
centrifugal effects and higher order pressure terms in a constant-density
approach. It seems rather that, if the constant-density solution is to
be improved on, the compressible flow equations should be solved exactly
as done in references A, 6, 7, and 8 for axially symmetric and two-
dimensional bodies.
COMPARISON WITH EXPERIMENT
The only pertinent experiments known to the author are those at
M_ = 3 and M_ = 8 by Zaltkay and Visich (ref. 5). The three-dimensional
body tested was a toroid. The vital statistics of the experiments and the
theoretical comparison are given in the following table (shock radii un-
fortunately had to be measured from the schlieren photographs presented
in ref. S).
12
Experiment - (ref. 5)
Rb,x/%, z
K = Rs,x/Rs_ z
A/Rb, x
£v/Rs _x
Theory
k
_/Rs_x (from fig. 2)
(eq. 18)
(eq.is)
0 255
4.5
4.08
194
0 259
208
187
186
Moo :8
O. 255
• 34.
• 267
.152
0.180
.169
.157
.156
The experimental data for _/Rs_ x are compa:'ed with three theoretical
estimates. The first is that value taken di:'ectly from figure 2(a) for
the pertinent density ratio and shock radius ratio. The second calcula-
tion is according to the relation
F fa3D ]
]RS_X \ s /Van uyKe[_- _ss 7 (18)
where the ratio of three-dimensional to axial ly s_mmetric standoff dis-
tances for the proper density ratio k is t_Lken from figure 5. The
third calculation uses equation (2) for the _forementioned ratio:
Rs,x Rs /Van I (19)
For both equations (18) and (19) the axially ssqmmetric standoff distance
is taken from Van Dyke and Gordon (ref. A), _rhose data are partially
shown in figure i.
The theoretical estimates agree well wi;h the experimental standoff
distances_ in fact, the agreement is better -,han might be expected_ con-
sidering the author's ability to measure sho_k radii from photographs•
The results using equation (iS) or (19) are _dthin 5 percent of the
measured values. The theoretical estimate f:'om figure 2 is somewhat
higher_ which indicates primarily that the p:'esent constant-density ap-
proximation overestimates the absolute standc)ff distance for all bodies.
This has already been shown for two-dimensioz_al and axially symmetric
bodies in figure I.
bd
!
Oo
13
CONCLUDING REMARKS
The shock standoff distance ahead of a general rounded stagnation
point has been estimated under the assumption of a constant-density
shock layer. It is found that many of the present theoretical calcula-
tions as well as the two experimental points of Zakkay and Visich are
well represented by the relation
A3D
Rs,x \ RS /
where K is the ratio of the smaller to the larger principal shock ra-
dius. The exceptional cases are bodies approaching the two-dimensional
(K % 0.2) with shock layers whose density is much larger than that of
the free stream (k < 0.i0). In comparing with experiment, the axially
symmetric standoff distance was taken from the exact solutions of Van
Dyke and Gordon.
Unfortunately, a constant-density theory gives no information re-
garding the body shape corresponding to a given shock wave and is there-
fore not suited for obtaining the standoff distance ahead of a given
body. In fact, from the experience with the two-dimensional and axially
symmetric problems, about the only solutions that adequately relate the
flow field to the body are the exact compressible flow solutions. The
present results may nevertheless serve as a semiquantitative guide to
the phenomenon in the absence of exact solutions.
Lewis Research Center
National Aeronautics and Space Administration
Cleveland, Ohio, June 5, 1961
REFERENCES
i. Hayes, Wallace D._ and Probstein, Ronald F.: Hypersonic Flow Theory.
Academic Press, Inc., 1969.
2. Hayes, Wallace D.: Constant Density Solutions. Ch. IV of Hypersonic
Flow Theory, Academic Press, Inc._ 1959, pp. 162-165.
3. Hayes_ W. D.: Some Aspects of Hypersonic Flow. Ramo-Wooldridge Corp.,
1955.
4. Van Dyke_ Milton D., and Gordon, Helen D.: Supersonic Flow Past a
Family of Blunt Axisymmetric Bodies. NASA TR R-I, 1959.
'-standoff distances.
Figure I. - Comparison of exact solutions for two-
dimensional and axially symmetric shock standoff
distances with Hayes' simple theory (ref. 3).
16
18
X
C)
-rq
q_
_H
O
0J
4_
og
O
O
,--4
O
52
2_
24
/
/i
16 / //1
i I I /Shock radius
ratio, K /
0i
/
•IC/
/ /
/ /._/
/ / // /
/" 9 /
.// / ./
/
f
0 .04 .08 .12 .16 .20 .24
Density rat o_ k
.2S
(a) ASD/R_;,x.
i Figure 2. - Shock standoff distance for three-dimensional bodies.
19
cO
h-
I
5.2
2.0
v
A
1. G
1.2k
.8 _'-
.4
0
Figure 2. -
bodies.
I
Shock radius
rat io_ K\ I
_ .50
i. 00
•04 .08 .12 .IG .20 .24 .28
Density ratio, k
(b) A3D/(kRs, x).
Concluded. Shock standoff distance for three-dimensional
20
\
!.0
0 .2 .4 .6 .8
Shock radius ratio_ K
Figure 3. - Ratio of three-dLmensional to axially s_mmetric
shock standoff distances.
NA_ - L_gley Fle_, Va. E-1278
bJ
!
(x)


Estimate of Shock Standoff Distance Ahead of a General Stagnation Point

Abstract

Similar works

Full text

Available Versions

NASA Technical Reports Server