19690025987
19690025987
19690025987
A. S. Predvoditelev, Editor
-0
-0-g
and Physics of Combustion
-
m
"W
a 2
Academy of Sciences of the USSR. Ministry of Power and Electrification of the USSR
G . M. Krzhizhanovskii Power Institute
A. S. Predvoditelev, Editor
HIGH-TEMPERATURE
PROPERTIES OF GASES
Physical Gas Dynamics, Thermodynamics, and Physics
of Combustion
(Svoistva gazov pri vysokikh teinperaturakh: fizicheskaya gazodinamika,
terinodiiiainika i fizika goreniya)
Copyright Q 1969
Israel Program for Scientific Translations Ltd.
IPST Cat. No. 5437
Translated by A. Barouch
Edited by P. Greeriberg
Available f r o m the
U. S. DEPARTMENT O F COMMERCE
Table of Contents
supersonic nozzles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
NABOKO. I.M. and R.G.NEMKOV. Study of the flow state behind a shock
o f a blunt body . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
71
KON'KOV, A. A., A. P. RYAZIN, and V.S.RUDNEV. Experimental study
hot nitrogen. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
90
BASHILOV, V.A. and V.P.SHEVELEV. Flow structure of ionized gas i n an
...
111
. . 127
PREDVODITELEV, A.S. Ignition theory of two-phase mixtures . . .
. . . 153
FEDOSEEVA, N.K. Apparatw for projecting stills . .. .. . .. .. .. . . .
. . 157
ZAITSEV, S.G., E. V. LAZAREVA, and E.I.CHEBOTAREVA. The boundary
. . . . . 160
process. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
180
in a supersonic stream , .. .
. .. ..... . . . . .. . .. . .
. . . . .. 196
VORONTSOV, Yu. N., V. M. EROSHENKO, A . V. hlORAR, and Yu. N. TERENT'EV
iv
I
ANNOTATION
vi
E . V . Samuilov
+
and Z = Z*/(1 Rx). Function (3) w i l l therefore replace (1)in what follows.
For R ravthe electron interacting with the particle may be considered
as having come f r o m infinity.
Following /4/, the motion of an electron in a field cp ( T ) reduces to the
motion of the electron in the effective field
where rg) i s the value of P for which 'pl (r) is a maximum in the interval
ro<r<w.
At given electron velocity t h e r e corresponds to each Z and 2 a maximum
value of the impact p a r a m e t e r bg) f o r which the electron becomes incapable
of reaching the particle surface. The values of b$ and r z ) a r e determined
for given electron velocity f r o m the s y s t e m of equations
Since 9 (0) = 1 and 'II, (00) = 00 equation (8) has no positive r e a l solutions
> <
f o r 11 ( p a ) 0 , while for g ( p z ) 0 two positive roots exist, of which p m ,
2
<
situated in the range p 1 p,,, Q p 2 , corresponds t o the maximum value of
>
the function cp;. and t h e other one pma p z corresponds t o its minimum
value for pml < p < w. When Z <
- 1 equation (8) has only one positive root
for pml > p a . This root corresponds t o the maximum of the function cp;.
The minimum value of the energy Eminfor which the electron is incident
on the particle is a s follows:
for Z >
0
E,," = 0,
for 2 Q - 1
where
By solving equation (8) relationships were derived for pnl and pmz a s
functions of Bmaxfor different Z . The values of p,, (( 1 . The dependence of
B,,,,, on E for different Z and finite ranges of values of E is approximated
by a linear function of the f o r m
TABLE 1.
-~ ~
Z =i
-__ - .
I
If, in expression (3), only the Coulomb interaction energy is considered,
(12) is replaced by
dNinc = f ( V ) V b d b d e d V , (15)
n, = 5 f(V)dV. (16)
Ma
B = - mRea , e2R
and E = - -mVa
2
dNinc= 5
f sin6d6dsde'dBdE. (21 1
2ma
4
With the aid of (21), (22), (17), one may derive by integration the total
number of electrons incident on the particle p e r unit time:
where rminis the minimum distance between the electron and the
particle during collision (the maximum positive root of the equation
b 2
1 - mV'V
- (y)z ( r ) = 0).
The angle of deviation for the potential ( 1 ) can only be calculated
numerically. For the potential ( 3 ) the expression for X, can be reduced
t o an elliptic integral. In t e r m s of the variables E, B , y = R / r , expression
( 3 ) can be used to transform integral ( 2 4 ) t o the f o r m
CD(y) = 0. (26 1
5
If only the first t e r m is taken into account in ( 3 ) , then
1
xz = 2 arc sin
f
The relationship
( 3 6 ) becomes
(38)
6
In t e r m s of the variables V 2 ,Ma, e, where mF= +cp(ro), (38) has the
form
.~
1
d V 2d M a de. (39)
4(rd2
e4
-dBdE de.
2 (romIa
dN,, = 4xae4
- 2-
e hTp
(y+al*-!+(r,))d B d E de;
E > TJ:(PI. (42 )
V b db de d V . (45)
where
where
for Z g - I
2
S (2,x ) = zexz.
The Boltzmann equation for the electron distribution function has the
f o r m 151
8
thermionic emission of electrons f r o m the s u r f a c e of p a r t i c l e s with effective
charge ( 2 - l ) e , inelastic collisions of electrons with particles of effective
charge Ze, and elastic collisions of electrons with particles of effective
charge z r .
Inelastic collisions of t h e electrons with atoms o r ions a r e not taken into
account, since t h e i r relative weight i s s m a l l f o r low values of E .
The following calculations a r e realized under t h e assumption T,= T i =
= T, =T.
Using ( 4 5 ) and ( 1 5 ) (according t o I S / ) , and taking into account that the
m a s s of particle, atom, o r ion p > m . we can write
In expressions (52), (54), ( 5 5 ) f ' = f ( V ' ) , where V ' i s the velocity of the
electron relative t o t h e p a r t i c l e a f t e r an elastic collision with it.
The solution of the Boltzmann equation ( 5 0 ) is taken in the f o r m
f = f ( 0 ) (1 -/- a).
Under the assumptions made relative t o the distribution function f , the
left m e m b e r of equation (50) can be written in the folowing form, using the
standard transformation procedure 15 1 of the kinetic theory of g a s e s :
where
9
The s u m of t h e integrals of the inelastic collisions in the left m e m b e r
of equation (50) is trazsformed t o the f o r m
@=-A a. -In- BT
ar
.aYo +
e. ar 12 Cd, (62)
c = w c (M.2), (63)
10
where
where
m ) X'X
(731
where z = eaJRkT.
11
5. ELECTRICAL CONDUCTIVITY
The c u r r e n t density is
-
j = n,eV, (74)
where v = -'
ne S fVdV. Using (57), (62), (67), we obtain f o r the ohmic
component of t h e c u r r e n t density
Emin
+ 1
oj
~e-xEII-cosxz] EdbdE
n o
Q:*')* is the integral of the collisions of the electrons with atoms o r ions of
type a having the usual f o r m 151.
We s h a l l calculate Q;,')', retaining only the f i r s t t e r m in (3). Using (14),
(32), (77), and (78), w e obtain
Bibliography
13
I. B . Rozhdestvenskii, V .P . Sheueleu, K . K. Olevinskii
14
end of the paper examples a r e considered of realizing this model in a t h r e e -
a d d r e s s computer .
where ek is a symbol designating the type of atom k , and !Ii,( is the number
of given atoms in the i-th component (in a l l t h e r e a r e k - 1 types of a t o m s ;
the l a s t index k always corresponds t o e l e c t r o n s ) .
It is likewise possible t o expand [zil not only with respect t o a l l the atoms
entering the mixture composition, but a l s o with respect t o k - 1 molecular
components selected in s o m e specific way. The molecular symbols a r e
designated by e; . Then
15
.. ... .. . .
l i n e a r algebra, where [xi]is regarded a s a vector in the space of base e;
or e*.
To that end w e select the matrix* 1 A , 11 containing initial information on
the compos ition.
In o r d e r to t r a n s f o r m to a general [I @G 11 one must find the t r a n s f o r m a
tion matrix 1 Bik 1, which is the inverse of the transposed matrix consisting
of the coefficients Ai,for any base e; in,the space of base e,. It is only
n e c e s s a r y that the molecules forming e k contain at l e a s t once any atom ek
and not be a l i n e a r combination of each other (such a s 0, and 0).
We shall give an example of such a transformation. Take t h r e e types of
a t o m s 0, N, C, and e-, and consider the combinations CO,, CO, N2,N,O+, C2N2,e;
The matrix )I A,, 1 w i l l have respectively the following f o r m for them:
ONCe
2010
11010
j 0220
0001
, I , ,
Take now a s el, e2, e,, e4 respectively COz, CO, N, and e-.
The transformation matrix 1 B,, 1 is then found f r o m the matrix
2010
I010
0200
0001 .
1 00 01 I.
The elements of @)ik for NzO a r e hence found a s
16
multiplying t h e s e vectors by the matrix IIBixII. The l a t t e r can easily be
found with the aid of the adjoint of the matrix \ l A i k l \ e k . i . e . , the matrix
formed of the coefficients A , in the space e k :
where IIAikIlel is the adjoint matrix, and d is the determinant of the matrix
.
I1 A ik Ile,
With the aid of the matrices IIAik[Io r !!(Dikll the fraction of an a r b i t r a r y
component xi is written a s follows on the strength of the law of m a s s action:
or
2 N~ - N = 0,
i
2 AijNi - cjN0 = 0,
i
17
where Ni is the number of moles of component i; N is the total number of
moles: A i j a r e the elements of the matrix of stoichiometric coefficients;
cj a r e the material constants. For the electrons cj = 0. The index j takes
the values 1 t o k . The components i in these relations a r e expressed, with
the aid of the law of m a s s action, in t e r m s of k a r b i t r a r y components, but
include necessarily each type of the k atoms and electrons.
For determinacy we shall express the components i in t e r m s of the
number of moles of atoms and electrons (cf. above):
For the quantities with index k f r o m the totality i we obtain, with the
aid of (13), the identity iV,<= N , and K,, = 1. Substituting (13) in (11)and
(12), we obtain a s y s t e m of k +
1 nonlinear algebraic equations for
determining the unknowns N , , N2,... N,.N. The s y s t e m ( l l ) , ( 1 2 ) w i l l be
completely determined i f the s e t of numbers c j . the equilibrium constants
KPi, and the matrix of stoichiometric coefficients Q, a r e specified.
Applying Newton's method t o ( l l ) , ( 1 2 ) and solving this s y s t e m with
respect t o the logarithmic corrections A In N,, A In N, we obtain
18
3 . CALCULATIONS IN DIFFERENTIAL
THERMODYNAMICS
Thus, it is easily seen that the determinants of systems (24), (25), and
(26), ( 2 7 ) coincide with that of s y s t e m (14), (15) in the l a s t iteration cycle.
When using the standard subroutines for solving-systems of l i n e a r
equations, the matrix of coefficients is not conserved in the c o u r s e of the
19
calculations and cannot be used f o r solving other s y s t e m s . It is therefore
advisable t o include the solution of the s y s t e m with respect to the derivatives
in every cycle of iterations of the s y s t e m f o r composition (14), (15). Since
the systems differ in t h e i r right sides, it is convenient t o u s e a solution
p r o g r a m based on the Gauss method with selection of the main element and
with many right s i d e s (cf. below). The solution of t h e s y s t e m with respect
t o the derivatives is obtained in this c a s e in the last cycle of iterations of
the composition system. The matrix of coefficients of the determinants and
the matrix of the right s i d e s of the s y s t e m s a r e
1
N- A;i IIT 2( O i k H k-If i) N i 2 (1 - 2 Oik)N i
1 r
I
Ci-xAijNi
.
4T~ A i j ( @ i x H , - H i ) N i ~ A i( j1 - 2 ( D i , < )
t
N,
20
t h e thermodynamic quantities and composition of t h e mixture corresponding
t o given T and p . It can a l s o be used as a subroutine in the p r o g r a m t o solve
the gas-dynamic problem. Both when used independently and a s a subroutine,
t h e p r o g r a m n e c e s s i t a t e s a s input data, in addition t o T and P, a l s o the
initial approximation .VI,N , , ..., N , , N . The p r o g r a m itself u s e s the l i b r a r y
of standard subroutines available in the IS-2 system, and a l s o a p r o g r a m
f o r solving a s y s t e m of l i n e a r equations by the Guass method with selection
of the main element. The advantages of this program, besides the accuracy,
a r e twofold: f i r s t , the selection of the main element, which makes it
possible t o avoid division by z e r o when dividing by an a r b i t r a r y element of
t h e s y s t e m determinant; second, by t h e Gauss method the expanded matrix
comprising the determinant coefficients and of s e v e r a l right s i d e s is
inverted, which enables the computing t i m e t o be considerably reduced.
Several comments a r e now made as r e g a r d s the possibility of using the
p r o g r a m with a computer having an internal memory of about 4000 t h r e e -
a d d r e s s c e l l s . We r e c a l l that the formation of the coefficients of l i n e a r -
s y s t e m m a t r i c e s is realized with t h e aid of t h e m a t r i x of stoichiometric
coefficients Ai,:. This matrix occupies ik cells. Since it is constantly in
use, it must be s t o r e d in the internal memory. For s m a l l values of iand
k (i- 100, k = I O ) the matrix & i s s t o r e d in such a w a y t h a t a cell corresponds
t o every m a t r i x element. For numbers i and k such that ik 2 1 0 0 0 , s e v e r a l
c e l l digits a r e assigned t o each matrix element, which is coded in binary
f o r m . Every vector i in space k is s t o r e d in two c e l l s of t h e memory, and
a s a r e s u l t the matrix occupies 2i c e l l s . Provision was made for modifying
the number of digits assigned t o each element Aikr depending on the maximum
content of the atom of type k in particle i. Independently of the method of
coding the matrix Ai,, t h i s unit r e a l i z e s the formation of the matrix of
determinant coefficients and t h e right s i d e s of l i n e a r s y s t e m s and t h e i r
distribution in a definite o r d e r for a r b i t r a r y i and k . Note that the maximum
o r d e r of the l i n e a r s y s t e m for the existing standard subroutines does not
exceed 60. We shall consider two optimum possibilities of the p r o g r a m .
In the first c a s e the m a t r i x Ai,, is s t o r e d in ik c e l l s , and t h e coefficients of
the polynomials cp?" a r e s t o r e d in the internal memory (i. e., the whole
p r o g r a m is s t o r e d in the internal memory). Possible i and k can then be
determined f r o m the approximation relation 15i ik 4-k + + 10k N 3200.
F o r k = 10, we obtain i 2 120; for k = 15, i N 90.
In the second c a s e the matrix A.is s t o r e d in 2i cells, and the coefficients
of the polynomials rp!"'are s t o r e d in the internal memory (i. e., the p r o g r a m
occupies the internal memory and one group of d r u m s ) . In this c a s e
possible i and k a r e determined f r o m the relation 7i k2 10k N 3200. + +
F o r k = 20, we obtain i N 370; for k = 40 i = 150.
Bibliography
1. S t u p o c h e n k o , E . V . , I . P . S t a k h a n o v , E . V . S a m u i l o v , A.S.
P l e s h a n o v , and I.B. R o z h d e s t v e n s k i i . Termodinami
cheskie svoistva vozdukha v intervale t e m p e r a t u r 1000- 1200K
(Thermodynamic P r o p e r t i e s of A i r in t h e T e m p e r a t u r e Range
1000- 1200K). - Sbornik "Fizicheskaya gazodinamika, I t Moskva,
Izdatel'stvo AN SSSR. 195 9.
21
2. S t u p o c h e n k o , E . V . , E . V . S a m u i l o v , A.S. P l e s h a n o v , and
1 . B . R o z h d e s t v e n s k i i . Termodinamicheskie funktsii
vozdukha p r i vysokikh temperaturakh (Thermodynamic Functions
f o r Air at High T e m p e r a t u r e s ) . - Zhurnal Fizicheskoi Khimii,
Vol. 34, No. 6. 1960.
3. H i r s c h f e l d e r , J . , C . C u r t i s s , and R. B i r d . Molecular Theory
of Gases and Liquids.-John Wiley and Sons, New York. 1949.
4. G u r v i c h, L. V. and others. Termodinamicheskie svoistva individual'
nykh veshchestv (Thermodynamic P r o p e r t i e s of Individual
Substances).- Moskva, Izdatel'stvo AN SSSR. 1962.
5 G u r v i c h , L. V. and N.P. R t i s h c h e v a . Analiticheskoe predstavlenie
tabulirovannykh znachenii termodinamicheskikh svoistv gazov
(Analytic Representation of the Tabulated Values of the Thermo
dynamic P r o p e r t i e s of Gases). - Teplofizika vysokikh t e m p e r a t u r ,
No. 6. 1964.
6. S a m u i l o v , E . V . and K . K . O l e v i n s k i i . Tablitsy s u m m d l y a
r a s c h e t a termodinamicheskikh svoistv gazov (Tables of Sums for
Calculating t h e Thermodynamic P r o p e r t i e s of Gases). - Moskva,
Izdanie VTs AN SSSR. 1963.
7. A r a v i n , G.S. and V . P . S h e v e l e v . Termicheskaya ionizatsiya i
elektroprovodnost' nekotorykh smesei i produktov sgoraniya
( T h e r m a l Ionization and Electrical Conductivity of Some Mixtures
and Combustion Products).-PMTF, No. 2. 1962.
a. P l e s h a n o v , A. S. Obshchii r a s c h e t sostava i termodinamicheskii
analiz proizvol'nykh reagiruyushchikh gazovykh s i s t e m (General
Calculation of the Composition and Thermodynamic Analysis of
General Reacting Gaseous Systems). - Sbornik "Fizicheskaya
gazodinamika, " Moskva, Izdatel'stvo AN SSSR. 1962.
9. P 1e s h a n o v , A. S. and S. G. Z a i t s e v. Sostav, termodinamicheskie
i gazodinamicheskie svoistva uglekislogo g a z a p r i vysokikh
temperaturakh (Composition, Thermodynamic, and Gas -Dynamic
P r o p e r t i e s of Carbonic Acid at High T e m p e r a t u r e s ) . - Sbornik
"Fizicheskaya Gazodinamika. " Moskva, Izdatel'stvo AN SSSR.
1962.
10. P 1e s h a n o v , A. S. Sostav, termodinamicheskie i gazodinamicheskie
svoistva azota p r i vysokikh temperaturakh (Composition, Thermo
dynamic, and Gas-Dynamic P r o p e r t i e s of Nitrogen at High
T e m p e r a t u r e s ) . - Sbornik "Fizicheskaya gazodinamika, " Moskva,
Izdatel'stvo AN SSSR. 1962.
11. E p s h t e i n , P. S. K u r s termodinamiki (A Course of Thermodynamics). -
Moskva-Leningrad, GITTL. 1948.
22
I . B . Rozhdestvenskii, K. K. Olevinskii, V . P . Shevelev
h, kca
2wl
zoo0
1200
800
2.2 2.4 2.6 2.8 3.0
s , kcal/kg. deg
FIGURE 1. h - I -diagram of the combustion products.
23
FIGURE 2. Composition of the combustion products as a function of the oxidizer
excess coefficient a t a temperature of combustion T = T, and a pressure in the
combustion chamber p = p c = 1 a m .
24
molecular weight p (g/mole), specific heats c p . c, (cal / g - deg), = c,/c,,
sound velocity n ( m / s e c ) , density p (g/cm3). The calculations were
conducted as a function of the oxidizer e x c e s s coefficient a for the following
p r e s s u r e s in the combustion chamber: p c = 1, 1.5, 2, 3, 4, 6, 8, 10, 15, 20,
30, 40 atm. The thermodynamic functions of the pure components (polyno
mial coefficients) for h and s were taken f r o m a handbook / 2 / . The
equilibrium constants f o r the reactions 02+e - - 0 2 - 2 0 and C + e - - - - Z 0
a r e taken from / 3 / . The dissociation or ionization energies and the enthalpy
r e f e r e n c e levels at 0K a r e given in the table for the main components.
0 IO 20 30 40 Pc 0 1 2 3
FIGURE 3. Combustion temperature a t a = 1 a s FIGURE 4. Combustion temperature p c = 1 atm
a function of the pressure in t h e combustion as a function of the oxidizer excess coefficient.
chamber.
Dissociation,
ionization Enthalpy at
energy, )OK, cal/mole
cal/mole--
CaHBO1+ 3 O ~ - Z C O f + 3 H 2 O 1 +W,
where t h e s u p e r s c r i p t s 1 and g r e f e r t o liquid and gas respectively.
25
T. "K
FIGURE 6. Specific-heat ratio y = c,/c, of the combustion
products a t p = 1 as a function of the temperature for dif
ferent a.
26
The heats of combustion of CzH60in liquid and gaseous s t a t e s for T =
= 298"K, following /4/, are equal respectively t o W(Cz%01)-326.66 f
f 0.10 kcal/mole, IV(CzH,0h')-336.78 f 0.10 kcal/mole. It follows that the
heat of evaporation of alcohol at 298K (q,::~.o) is equal t o 10.12 kcal/mole.
The products of combustion a r e in both c a s e s liquid water and gaseous
carbon dioxide. The heat of evaporation of HzO at 2 9 8 K ( q ~ ~is, ~equal
) to
10.503 f 0.006 kcal/mole.
We can thus write t h e following reaction equations:
+
These equations yield w - W ( C , I ~ ~ O ' ) qC&O - 3 qH,O = 305.271 kcal/mole
f o r the heat of combustion of gaseous alcohol into gaseous products. At the
enthalpy r e f e r e n c e levels adopted the enthalpy of Cz%Og at 298K is found
t o have the value 310.628kcal/mole (fIY&o = 310,628 c a l / m o l e ) .
The calculation w a s performed f o r 96 70 (by volume) alcohol. The quantity
It.'includes the heat of evaporation of the water contained in 96 70 alcohol.
The theoretical combustion t e m p e r a t u r e w a s determined by l i n e a r interpola
tion of the enthalpy values with simultaneous calculation of the h - s - d i a g r a m s
of the combustion products.
The range of the k - s diagram (cf. Figure 1 ) which is of practical
interest for calculating the efflux through a nozzle is the lower region,
lying beneath t h e horizontal line of the corresponding combustion t e m p e r
a t u r e ( h ," = h p r ) .
B i b 1i o g r a p h y
1. R o z h d e s t v e n s k i i , I . B . , V . P . S h e v e l e v , and K.K. O l e v i n s k i i .
Raschet sostava i termodinamicheskikh funktsii proizvol'nykh
reagiruyushchikh gazovykh sistem (Calculation of the Composition
and Thermodynamic Functions of General Reacting Gaseous
Systems).-In this collection, p. 14.
2. Gu r v i c h, L. V. and others. Termodinamicheskie svoistva individual'
nykh veshchestv (Thermodynamic P r o p e r t i e s of Individual
Substances).- Moskva, Izdatel'stvo AN SSSR. 1962.
3. P l e s h a n o v , A.S. and S . G . Z a i t s e v . Sostav, termodinamicheskie i
gazodinamicheskie svoistva uglekislogo g a z a dlya t e m p e r a t u r ot
1000 do 12,000"K i davlenii ot lo-' t o l o 3 a t m (Composition,
Thermodynamic, and Gas-Dynamic P r o p e r t i e s of Carbonic Acid
for T e m p e r a t u r e s between 1000 and 12,000"K and P r e s s u r e s
between lo-' and l o 3 atm). -Sbornik "Fizicheskaya gazodinamika,"
Moskva, Izdatel'stvo A N SSSR. 1962.
4. L e w i s , B. and G. von E l b e. Combustion, Flames and Explosions of
Gases.-Academic Press, New York. 1951.
21
N . N . Tsitelauri
28
r
L, p e c The given calculation of the p a r a m e t e r s
mi of the intermolecular interaction potential
of CO is based on the experimental r e s u l t s
summarized in 151. The authors d e t e r
mined the relaxation time z over a wide
101 temperature range f r o m the variation of
the infrared emission in the g a s behind a
shock wave. The inaccuracy of t h e i r
measurements (of t h e o r d e r of 15 70)was
JI;
due t o e r r o r s in the front velocity, the
temperature and p r e s s u r e behind the front,
the initial t e m p e r a t u r e and p r e s s u r e .
E r r o r s due t o impurities, and a l s o the
I effect of the boundary l a y e r and the self-
008 01 absorption, could be neglected. Their
r-b r e s u l t s 15 1 a g r e e satisfactorily with those
FIGURE 1. Vibrational relaxation time of of others, who measured the time of
CO niolccules as a function of temperature. vibrational relaxation by means of infrared
emission / 6 1. bv the line-reversal method
- I
IC! z = 75.8T"'-10.lf. (1 1
29
t h e internuclear distances a r e sufficiently l a r g e compared with t h e
r a d i u s of action u-l of the intermolecular forces, and t h e r e f o r e only the
interaction of t h e n e a r e s t a t o m s is important, while that of the second
atoms can be neglected.
The full intermolecular potential w i l l then equal the interatomic potential
given by
V = V,, exp (- aR),
V = V o exp (- a R ) - E . (3 1
- In A =( ~ T-'".
3 /a'/,) (5 1
30
V , = E exp [ u ( s - - L ) ] . (7 1
U o e x p ( - - a , , r , ) = E ' =L ~ ~ , (8)
I Nz-Nz
where u* = (4n2kTAE/~hp)':a is the velocity of
4
the molecules which a r e most effective
I in exciting vibrations 141.
I The whole procedure is then repeated
I
I with the new values of r , . This p r o c e s s
\ converges fairly rapidly.
\ As a r e s u l t the following values w e r e
\
',Lennard -Jones
potential
obtained f o r the interactions of C and 0
atoms contained in difierent colliding
molecules: u = 4.899A1', V, =
= 2.624. e r g ; the correspondoing values
<,co-co f o r the molecules a r e aN=4.629A-', lJo=
18 22 26 30 ~ 3'i = 3.158. e r g . The value of the c o r r e c
r ,A tion f a c t o r A in equation ( 4 ) is equal t o 0.2.
The effective interaction potential obtained
FIGURE 2. C u r v e s of rhe interaction f o r the CO molecules,
potential of the molecules.
li = 19 700 r x p ( -4.629 r ) cli (9)
31
..
Bibliography
33
excited in a stable manner behind the wave, and dissociation has taken place
(variant I); a l l the vibrational modes a r e stable behind the wave, and t h e r e
is no dissociation (variant 11); neither dissociation has taken place n o r have
a s y m m e t r i c valent vibrations been excited in the g a s behind the wave
(variant 11'). The graphs compare curves for a constant specific-heat r a t i o
y = 1.4 during transition a c r o s s the shock (variant 111). Variant I1 was used
f o r nitrogen, since the degree of dissociation of nitrogen molecules is low
f o r the corresponding shock velocities. The calculation of the values of
the enthalpy and molecular weight for variant I is described below. The
values of the enthalpy f o r variant I1 were taken f r o m the tables in / 6 / , and
the values of the enthalpy for variant 11' f r o m /7/.
The calculations were f i r s t conducted for the gas s t a t e in front of the
shock w i t h p i n = 1.68, lo-' a t m and Ti, = 293"K, and wave velocities between
1000 and 3000 m / s e c . Calculations were also performed for shocks under
higher initial t e m p e r a t u r e s and p r e s s u r e s . These initial s t a t e s a r e obtained
a s a r e s u l t of the passage through a gas (whose initial temperature is 293K)
of shock waves moving at one of the following velocities: u = 1.0, 1.3, 1.6,
1.9, 2.2, 2.5, 2 . 8 k m / s e c .
TABLE 1. C 0 2 . Variant I
~.
293 I 17i0 2185 1 2690
- 1.75 1
k
5
.
Q
-
ha
- b"
0.7 1.84 7.2 4.23 1.03 1.33 1.29 1.00 .2? . .22 .04 - -
0.8 2.08 10.2 4.92 1.07 1.73 1.61 1.03 .41 . .43 .23 .00f 1.14
0.9 2.36 13.1 5.61 1.11 2.20 1.97 1.06 .85 .74 .51 .02 1.35
1.0 2.62 16.1 6.16 1.15 2.72 2 36 1.08 t28 1.09 .84 .04 1.66
1.1 2.93 19.6 6.69 1.19 3.31 2.77 1.11 .76 :.45 .19 .06 1.96
1.2 3.26 23.4 7.17 1.22 3.96 3.2t 1.14 .29 I. 83 .55 .07 2.30
1.3 3.60 27.5 7.65 1.26 4.67 3.65 1.17 .88 1.24 .93 .09 2.63
1.4 3.97 32.0 8.06 1.31) 5.42 4.11 1.20 .50 :.65 .33 .tt 2.98
1.5 4.35 36.8 8.46 1.33 6.16 4.59 1.22 .20 .I3 .74 .I3 3.25
1.6 4.75 42.0 8.85 1.37 7.14 5 06 1.23 -94 .63 .16 .I4 3.71
1.7 5.19 47.5 9.16 1.41 8.08 5.53 1.27 .70 1.02 .59 .I7 4.12
1.8 5.62 52.3 9.45 1.45 9.08 6.03 1.29 .54 1.54 .02 .18 4.56
1.9 6.04 59.2 9.81 1.49 10.1 6.48 1.31 '.43 i.11 .41 .20 4.93
2.0 6.42 65.820.2 1.52 11.2 6.97 1.34 .34 i.52 - .22 5.34
2.1 6.83 72.7 10.6 1.56 12.4 7.47 - - - - -
2.2 7.23 80.010.9 1.59 13.6 7.95 - - - - -
2.3 7.58 87.711.4 1.63 14.9 8.42 - - - - -
2.4 7.88 95.711.8 1.66 16.3 8.89 - - - - -
2.5 8.14 04 12.3 1.69 17.7 9.38 - - - - -
2.6 8.45 13 12.7 1.72 19.2 9.88 - - - - -
2.7 P.70 22 13.2 1.75 20.7 10.3 - - - - -
2.8 8.94 31.513.6 1.78 22.2 10.7 - - - - -
2.9 9.16 41 14.1 1.81 23.9 11.1 - - - - -
3.0 9.39 51 14.6 1.84 25.6 11.6 - - - - - -
- I
~ . .-
The tables give the ratios of the p r e s s u r e pa, the temperature Tb, and
the density pa behind the wave to the initial p r e s s u r e pin, t e m p e r a t u r e T,n,
34
and density P , ~ ,respectively, a s a function of the wave velocity u . The
initial p r e s s u r e and temperature a r e indicated for variant I. For variants
I1 and 11' only the initial t e m p e r a t u r e I",, is indicated, since in this c a s e
the specific enthalpy of the gas behind the wave i s a function of t e m p e r a t u r e
alone, and the given ratios depend only on the velocity of the shock and the
initial temperature. The calculations have an e r r o r of 1-2 70. The data
can be used for determining the gas p a r a m e t e r s behind reflected waves in
a shock tube, and behind waves at obstacles, when a hot gas flows past
them. These r e s u l t s were used t o determine the gas p a r a m e t e r s for shock
diffraction a t a wedge (Tables 1-4).
T,,,O K
1390
__ -
51 TI;.
3 - 3
P
h a
a h
-P
35
5 h" 0
-
0
2810
u - 3_ _ - 1400 1800 !- -
2 h .5
,o
,Pc
< P
8
8 -
5 .a
1
5
- h
- a
36
. .. . .
X0Xc
-- KPCO
BCO - P '
where H , is the molar enthalpy of the components; p,, is the molecular weight
of the components; n is the number of components.
Such calculations w e r e conducted for T = 1100-4500K with increments
of 200" and p r e s s u r e s of 0.1, 0.5, 1, 5, 7, 9, 11, 1 3 atm.
Tables 5 - 11 give data for h, I.L,Xc0,, X C O ,XO,,X o , X c .
The numerical material of the tables is given in normalized form.
The enthalpy of the gas C 0 2 at To = 293.16"K was taken a s 43.485 cal/g.
The thermodynamic functions of the mixture components, given in / 6 /,
w e r e used t o calculate the composition and thermodynamic properties of
carbon dioxide.
The reference levels of the enthalpies of the components, and the values
of the dissociation energies of the main components, a r e given in Table 1 2 .
A comparison of the results obtained with those of 161, which allows for
a l a r g e number of possible reactions but gives values of the enthalpy for
temperature increments of lOOO", showed that the failure to allow for the
gas components figuring in very low concentrations is unimportant over
these temperature ranges. The discrepancy does not exceed 1 % . The
values of the enthalpies calculated in o u r case differ by 10-20 70 f r o m the
values obtained in / 6 / by linear interpolation, in particular over the
temperature range 2000- 3000K.
Consider the interaction between a shock wave and a wedge (Figure 1).
Let the shock wave SA be incident on a wedge of apex angle a,, and let the
wave velocity be u,,. An attached shock wave ON is then formed.
37
__
P. am
T. "K
-
0.1
I
0.5
1 A i A
2
T A
3
r
1
(
5
r ~
I 1
A , T I A P T
0.452
i.'a0.451 3 0.451
i3 l 0.520
0.451 3
1
1700 0.457 3 0.454 0.453 3 0.452 3 0.452 3 0.452 3
0.496 3 0.491 3 0.489 3 0.487 3 0.487 3 0.486 3 0 . 486 0.485 3 0.485 0.485 3
1800 3
0.800 3
2500 0.113 4 0.963 3 0.914 3 0.873 3 0.854 3 0.832 3 0.820 3 ' 0.811 3 0.805
0.151 4 0.124 4 0.115 4 0.108 4 0.104 4 0.100 4 0.980 3 0.965 3 0.953 3 0.944 3
0.161 4
0.222 4
3900
4100 0.405 4 1 0.388 4 0.373 4 0.352
n.379
4
4
0.338
0.368
4
4
0.320
0.351
4
4
0.308
0.340
4
4
0.299
n mi
4
4
0.292
0.324
4
4
0.286 4
0.318 4
4300 0.414 4 0.405 4 n.395 4
0.422 4 0.417 4 0.411 4 0.400 4 0.392 4 0.378 4 0.368 4 0:360 0.353 4 0.347 4
W
1700 ; 0.993 0 0.996 0 , 0.998 0 0.998 U n. w x 0 0.998 I) 0.999 [I,O.XI!I n
1800 0.986 0 0.992 0 0.995 0 ' n.995 0 [l.!197 0 0.997 0 0.997 0 0.9!17 0
lono 0.973 0 0.984 0 ' 0.990 n 0.991 0 0.903 0 0.994 0 ' 0.094 0 0.995
2100 0:jzi 0 0.953 0 0.9711 0 (1.974 n (I.9x0 0 ' 0.982 0 0.983 0 0.984 n
2300 0.815 0 0.886 0 0.927 0 0.Y3li 0 0.951 0 0.955 0 0.958 0 u.9tiu 1)
2500 0.649 0 0.773 0 U.S4!1 0 0.867 0 0. H98 0 0.905 0 0.911 0 0.916 0
2700 0.449 0 O.Gl7 0 0.734 0 0.7G3 II 0.814 0 0.827 0 ' 0.838 0 0.s45 n
2900
3100
0.266
0.138
0
0
0.442
0.283
n
11
0.58!1
n.435 0
I 0.(i'LS
11.470
11
II
0.7111
ll.%!I
0
0
0.721
0.594
0
0
0.736
0.613
0 11.748
iI U.UY!I
0
I1
3300 0.622 -1 0,164 o n.m 0 u 339 0 11.433 0 0.461 0 0.482 0 U.501 n
3500 0.265 -I 0.~7n 1) 0.187 0 0.224 11 n.xo 0 0.337 0 0.359 n 0.377 (1
n.xo 0 0.271 u
1
3700 0.112 -1 11.437 -1 0.111 0 0.139 U 0 0.234 0 0.254
3900 0.400 -2 0.214 -I (1.631 -1 11.827 0 (1.136 0.156 0 0.172 0 0.18U n
4100
4300
4500
0.227
1
0.112
0.587
-2
-2
-3
1 ll.10F
nso
(1.288
-1
-1
-2
,
I
0.,%0
0.193
0.108
-1
-1
0.477
u.271
-1 -. 1 0.155
--I
-1
-'i____
-1
0.855
11.524
0.318
1; 0.999 --I ' 0.112
0.626 -1
0.358
0.718 -I
u 0.124
~
0.801 -1
0.452 -1 10.511 -1
I1
- - ~
W
W
A r A r . A r1-A P I A r ( A r
1700
1800
0.451 -2
0.934 -2
0.123 -2
0.255 -2
0.110
0.228
-3
-2
I 0.102
0.210 -2 1 0.945 -3 10.894 -3
0.191; -2 0.18G -2
19ll(l 0.178 -1 0.489 -2 0.437 -2 0.402 -2 0.376 -2 ' 0.358 -2
2100 0.526 -1 0:314 -1 0.251 --2 0.200 -1 0.175 -1 0.148 -1 0.133 -1 0.122 -1 ' 0.114 -1 0.108 -1
2300 0.122 0 0.753 -1 0.607 --1 0.487 -1 0.428 -1 0.354 -1 0.3'26 -1 0.30.1 -I 0.282 -1 0.267 --I
2500 0.229 0 0.149 0 0.122 0 0.997 -1 0.88' -1 0.754 -1 0.679 -1 O.R'28 -1 0.590 -1 0,559 --1
2700 0.353 0 0.250 0 0.210 0 0.175 0 0.15ti 0 0.135 0 0.123 0 0.114 0 0.108 0 0.102 0
2900 0.455 0 0.357 0 0.312 0 0.268 0 0.243 0 0.214 0 0.196 0 0.lH3 0 0.174 0 0.1GG 0
3100 0.510 0 0.447 0 0.406 0 0.362 0 0.335 0 0.302 0 0.280 0 0.~5 0 0.253 0 0:~43 0
3300 0.5'24 0 0.502 0 ' 0.476 0 0.441 0 0.417 0 0.38li 0 0.3% 0 0.X7 0 0.334 0 0.323 0
3500
3700
0.519
0.510
3900 1 0.505
:'
0
0.524
0.524
0.517
0
0
0
0.514 n
0.5%
0.524 0
0
0.404
0.520
0.527
0
0
0
0.478
0.512
0.52G
0
0
0
0.453
0.407
0.520
0
n
0
0.43
0.484
0.514
0
u
0
0.410 o
0.473 0
0.507 0
0.407
0.464
0.501
u 0.396 I)
o 0.455 n
0 0.485 0
4100 0.503 0 0.510 0 0.517 0 0.523 0 0.526 0 0.527 0 0.53 0 0.524 0 0.521 0 0.518 0
4300 0.501 0 0.506 0 0.510 0 0.517 0 0.521 0 0.525 0 0.527 0 0.528 0 0.528 0 i 0.577 0
4500 0.501 0 0.503 0 0.506 0 0.511 0 0.515 0 0.520 0 0.523 0 0.525 0 0.527 0 0.527 0
P
0
- I A r I - r 1
A - A 1.1 A A r ) A r ! A r ) A 1.1 A I . ) A Y
1700
1800
1 0.GGG -5 I
0.261 -4
(1.28 -5
0.894 -5
I 0.144 -5
0.563 -5
I 0.906:
0.355 -5
-6 1
I
0.691 -6
0.2ii -5
0.492 -13
'
0.193 -5
I 0.393 -(i
0.154 -5
0.331. -8
0.131) -5
I 0.391 -6 I 0.2130-6
0.114 -5 0.102 -5
1900 0.882 -4 0.303 -4 u.ini -4 0.120 -4 0.920 -5 n.~55 -5 0.523 -5 0.443 -5 0.387 -5 0.348 -5
2100 0.703 -3 0.244 -3 0.15443 0.974-4 0.744-4 0.530-4 0.424-4 U.35!3--'t 0.614 -4 0.281 -4
2300 0.379 -2 0.134 -1 0.853 -3 0.541 -3 0.415 -3 ' 0.20(i -3 0.237 -3 0.2nl -3 0.lifi -3 0.158 -3
2500 0.149 -1 0.546 - 'J 0.351 -2 0.25 -2 0.173 --1! 0.124 -2 0.995 -3 0.844 - - 3 0.i40 -3 0.663 -3
2700 0.445 -1 0.173 -1 0.113 -1 0 734 -2 0.568 -2 0.412 -2 o.331 -Z 0.28: -2 0.248 _ ' I 0.2':! -2
2900 0.104 0 0.440 -1 0.300 -I 0.200 -1 0.154 -1 0.112 -1 O.!)lt -2 0.779 -2 0.G87 -2 O.ti18 -2
3100 0.196 0 0.932 -1 0.M8 -1 0.44'!-1 0.351 -1 0.2Iit - 1 0.213 -1 (1.183 -1 O.lli2 --I 0.t44T - 1
3300 0.302 0 0.16ti 0 0.121 n 0.855 -1 0.691 -1 0.523 -1 0.434 -1 0.9ili --Z 0.335 -1 0.304 -I
3500 0.391 0 0.254 0 0.195 n 0.145 0 0.120 0 0 . 9 2 ~ -1 0.781 -1 0.684 -1 n . w -1 0.560 -1
3700 0.446 0 0.339 n 0.258 0 0.217 n 0.184 0 0.148 n n.i2i; (I n.ii2 0 0.101 I 0.931 -1
3900 0.475 0 0.405 0 0.354 0 0.293 o 0.257 0 0.213 o 0.186 n o.iiii n 0.153 1 0.141 0
4100 0.488 0 0.448 0 0.411 0 0.360 n 0.326 0 0.281 o 0.250 n n.r2o 11 u . 2 ~ 1 0.198 0
4300 0.494 0 0.471 0 0.448 0 0.411 0 0.384 0 0.343 0 0.313 11 0.291 0 0.~73 1 0.258 n
45(M 0.497 0 0.484 01 0.470 n 0.446 0 0.425 o 0.393 0 I 0.368 u , o 347 0 O.:XIO 1 0.315 0
TABLE 11. XC=A.10'
T. 0.1 05 I 2 1 7 9 II 13
I 1
1700 0.488 -22 0.167 -22 0.105 -22 0.663 -23 0.506 -29 0.360 -23
1800 fl.186 -20 0.635 -21 n.c,iio -21 0.252 -21 0.193 -21 0.137 -21
1900 U.481 -19 0.165 -19 u. Ill4 -19 0.655 -20 0.500 -20 0.:35fi -20
2100 0.127 -19 11.436 -17 0.275 -17 0.174 -17 0.133 -I7 0.945 -18
2300 0.124 -14 0.434 -15 0.275 -15 0.174 -I5 0.133 -15 0.948 -16
2500 0.571 -1:3 0.203 -13 0.129 -I3 0.823 -14 0.631 -14 0.452 -14
2700 0.144 -11 0.526 -12 0.338 -12 0.217 -12 0.167 -12 0.120 -12
29nn 0.228 -IO 0.848 -11 u.551 -11 o.:m -11 0.276 -11 o.zno -11
3100 , 0.254 - 9 0.939 -10 0.614 -10 0.401 -10 0.312 -IO 0.227 -10
,3300 0.223 - 8 0.778 - 9 0.507 - 9 0.332 - 9 0.259 -9 fl.190 -9
3500 0.168 - 7 U.521 - 8 n.:m - 8 0.216 - 8 0.168 - 8 0.123 -8
3700 0.111 - 6 ii.:tun - 7 0.184-7 0.116 7 - 0.898 - 8 0.654 - 8
0.645 - 6 -7 -7
3900 0.155 - 6 11.898
0.399
-7
-6
-5
0.545 - 7
0.230 - 6
'
0.8Y5 - 6
' 0.414
0.171
0 645 -
-6
0.297
fl.119
0,437
-6
-6
0.162
0.604 - 5 0.322 - 5 1 0:226 -
6
SI n.148 -5
TABLE 12.
Dissociation Dissociation
energy at OX,
Component Reaction
1 energy at OX,
caymole 1
HT =O I' Component
I/ I
Reaction
1 cal/mole
HT = 11
+ 2c 2
+ co f ;02 3 0
(:2 --cz 0 143170
--CUI 66767 0
u - - 58986
--co + c + 0 2 0 255790 66767 c - - 263571
-02 + 20 2 0 117973 0
Designate by P , u. p . T,p, h the density, velocity, p r e s s u r e , temperature,
molecular weight, and specific enthalpy of the gas.
The subscripts 0, 1, 2, 3, 4, 5 r e f e r respectively to the g a s s t a t e ahead
of the incident wave, behind the incident shock wave, behind the reflected
wave, behind the Mach wave, ahead of the Mach wave (this s t a t e coincides
with s t a t e 0), and behind the attached wave.
42
.. . . ~.
..... .. .. . .. ..
TABLE 13. C02. Variant I
1770 20"06' 1.22 2nno 17"54' 1.76 2385 16"54' 2.22 2620
1830 23 06 1.63 2120 21 24 2.50 2450 20 06 3.12 2680
1890 26 48 2.18 2150 25 00 3.34 2520 23 48 4.22 2750
1970 31 00 2.86 2201) 28 48 4.32 2590 27 36 5.72 2830
2050 35 18 3.60 2270 33 24 5.60 2660 32 12 7.40 2910
2120 40 18 4.53 2340 38 00 7.04 2720 37 00 9.40 2990
2210 45 36 5.60 2420 43 00 8.60 2800 41 48 11.65 3070
35 00 54 12 5.26 2285 51 36 , 6.72 2490 48 48 10.5 2880 47 24 14.2 3150
'
40
41
44
00
30
00
58 24
41 48
-
j 6.40
7.10
-
2370
1
55 36
70 i
7.94
-
9.7 i
2570
1
52 24
-
- I
12.45
-
-
I -
1
2970
-
-
Id-
-
e
TABLE 14. C q . Variant I1
uo, kmlsec
1.9 2.2 2.5 2.8
O"O0' 21%' 0.71 1390 20"30' 0.98 1780 20"30' 1.33 19"48' 1.72 2730
5 00 25 00 O.YR 1420 25 00 1.30 1850 24 06 1.82 23 06 2.37 28511
10 00 28 36 1.24 1450 29 00 1.72 imn 27 48 2.44 27 00 3.m 29x0
15 on 33 00 1.54 1480 33 on 2.24 2025 31 48 3.14 31 12 4.20 3125
20 00 37 48 1.90 1540 37 42 2.86 21'0 36 24 4.00 35 48 5.40 3175
25 on 43 12 2.36 1610 42 48 3.55 zin 41 36 5.05 41 00 6.80 3430
30 00 49 06 2.86 1690 48 30 4.36 2320 47 00 6.20 46 18 8.33 3620
35 00 1 56 30 3.50 1790 55 48 5.35 2470 , 53 36 7.44 52 48 10.1 3830
9.20 1 61 00 1 12.3 4080
1
~ ~
0"OO 22"48'
1
11.99
I
1900 1 22"Ofl 1.34 2390 'Y"12' 1.74
I 2970
5 00 26 24 1.34 1990 1 25 12 1.80 2480 24 36 2 . :14 3110
10 00 3n 12 1.72 2080 29 00 2.:x 2610 28 34 3.10 3260
25 00 34 I? 3.20 2170 33 15 3.1111 2740 33 30 4 . IJ2 3420
20 00 38 48 2.76 2260 38 00 3.82 2890 37 00 5.ill xno
25 01) 44 12 :3.47 2.380 43 12 4.84 3040 42 30 6.4 3800
30 00 50 36 4.30 2500 411 34 6.00 3215 48 24 7.8R 4040
35 00 58 4 R 5.26 2680 57 on 7.28 3445 55 42 9.60 4300
37 18 68 00 6.10 2840 - - - - - -
38 06 - - - 67 I2 8.80 3680 - - -
0,
I 13 I 16
Ta." I i
__
10"OO 3140
15 00 3350
20 00 JSRO
25 00 3930
19 no - - - 65 00 1.311 -
22 12 -
24 36 -
25 36 4010
process. Thus, for avelocity ofthe incident shockwave equalto 1.9- l o 5 m / s e c
(Figure Z), the change in the specific-heat r a t i o y caused by exciting the
vibrational levels of the C 0 2 molecules causes the shock at the wedge apex
t o remain attached up t o wedge angles a. = Ib= 41"31'.
110
45"
50
, -13"36'
-14 96
92"24'
93 OG
2170
2130
2.50
2. I15
2730
2570
2.97
2.74
28'18'
33 06
55 -14 36 94 12 2100 1.83 2440 2.56 38 30
GO -12 42 95 90 2070 1.64 2320 2.41 43 30
65 -6 42 97 18 2050 1.48 2220 2.30 48 06
45
4. THE TRIPLE CONFIGURATION USING
SHOCK POLARS
where ti, is the gas velocity in the coordinate f r a m e linked with the incident
wave.
Construct the shock polar for the reflected wave in ( p . e) coordinates.
Let the reflected wave l i e a t an angle q1 t o the flow incident with velocity
6,; the normal component of the flow velocity will then be u,,, = El sin cp.
The value of p z behind the reflected wave i s determined from Tables 1-4
and 9 2 f r o m the continuity equation.
F o r different (p, a number of values of El" can be derived for correspond
ing e2. and p z , i. e., the shock polar in the p-8 plane. Construct now the
shock polar for the Mach wave. If the Mach wave l i e s a t an angle (p4 to the
incident flow, the normal velocity component t o the Mach wave w i l l be
U M = uo (sin cp4/sinq).
By assuming different values for uMr one can construct the shock polar
for the Mach wave in ( p , 8 ) coordinates, using Tables 1-4 and the continuity
equation.
In o r d e r t o find the required solution the compatibility condition must
be fulfilled, namely, equality of the p r e s s u r e s and parallelism of the flows
in regions 2 and 3; p z = p 3 and O3 = el- O2.
These conditions a r e fulfilled at the intersection point of the shock
polars f o r the reflected wave and the Mach wave.
The intersection point of the polars yields the required values of and
U1n.
All the angles of the Mach configuration a r e determined f r o m the
e x p r e s s ions
-
sin (p3 = "In, Oz=~z-el,
u1
sin (p4 = 5
uo
sin ml, 0 3 = qa.
Tables 17-20 give the r e s u l t s obtained using this method, and a l s o the
values of the gas-dynamic p a r a m e t e r s and the angles between the shock
waves forming the t r i p l e configuration a s a function of the velocity of the
incident wave.
46
P I = 0.46 atm
PI 0.706 atm
47
TABLE 19. Cq. Variant I1
uo = 1 . 9 !un/sec TI 1900" K
40" -8"36 92"24' 2220 2.40 3970 2.95 21'42'
45 -11 00 93 18 2125 1.97 3340 2.68 27 30
1
50 -12 I 8 94 18 2055 1.67 2900 2.47 32 54
55 -12 18 95 36 2010 1.46 2590 2.31 38 24
60 -10 00 97 06 1975 1.30 2360 2.18 43 30
65 4 30 99 12 1950 1.17 2170 2.07 47 54
2 km/sec TI= 2390' K
94012' 33042'
95 24 38 48
96 48 2.52 43 42
98 54 2770 2.40 48 36
u1 = 2.5 !un/sec Ti = 2965" K
60"
65 1 -11"12'
-6 06 1 96"4?'
98 42 [ 3090
3030 1 2.26
2.05 I 3760
3440 1 2.87
2.73
44024'
49 18
50" 3.25
55 3.03
60 2.85
40 50 60 70
UP
FIGURE 3. The reflection angle mt vs.
the incidence angle 0,:
1-ideal gas Y = 1.4; 2-unexcited asym
metric vibrations of the molecules and no
dissociation; 3,4-the curves overlap, the
vibrations are excited, there is no dissocia
tion, and complete equilibrium.
Bibliography
1. G v o z d e v a , L . G . and O . A . P r e d v o d i t e l e v a . Eksperimental'noe
issledovanie makhovskogo otrazheniya udarnykh voln p r i s k o r o s t i
1000-3000 m / s e k v uglekislom gaze, azote i vozdukhe (Experi
mental Study of the Mach Reflection of Shock Waves at Velocities
of 1000-3000 m / s e c in Carbon Dioxide, Nitrogen, and Air).-
Doklady AN SSSR, Vol. 163, NO. 5. 1965.
2. G v o z d e v a , L.G. and O . A . P r e d v o d i t e l e v a . Osobennosti
makhovskogo otrazheniya udarnykh voln, dvizhushchikhsya v
uglekislom gaze s o skorostyami poryadka 2000 m / s e k (Special
F e a t u r e s of the Mach Reflection of Shock Waves Propagating in
Carbon Dioxide and Nitrogen at Velocities of the O r d e r of
2000 m/sec).-Sbornik "Issledovaniya PO fizicheskoi gazo
dinamike, Moskva, Izdatel'stvo "Nauka. 1966.
3. Z e 1 d o v i c h, Ya. B. and Yu. P. R a i z e r. Fizika udarnykh voln i
vysokotemperaturnykh gidrodinamicheskikh yavlenii (Physics
of Shock Waves and High-Temperature Hydrodynamic Pheno
mena). - Moskva, Fizmatgiz. 1963.
49
50
G . N . Nikolaev, V .P . Popov
1. EXPERIMENTAL SET-UP
51
Figure 2 shows a nozzle with a wedge in its working position at the end
of the shock tube (the p r e s s u r e chamber is removed). Figure 3 depicts the
working section of the apparatus in assembled form.
The nozzles were axisymmetrical, the cross-sectional radius along the
axis obeying the hyperbolic law
ra = r2cr + (GIn) 9 9
52
-
. ...... ..-.
..... __ ... -. .. ..
where z is t h e distance f r o m the c r i t i c a l section; r is t h e c r o s s - s e c t i o n a l
radius corresponding t o the coordinate x ; K N = f i t g 6 ; 6 is the half-angle
of t h e asymptotic cone of t h e unparted hyperboloid of revolution r e p r e s e n t e d
by the nozzle in its supersonic p a r t .
. --
1i
- -
Cross-sectional Cross-sectional
Distance from Distance from
Nozzle expan- diameter, m m Nozzle expan- diameter, m m
throat to exit,
I -~ throat to exit,
sion ratio
8.1
26.6
t
10.15
h
I ",::1
r o ~ ~ m~m l
80
135
1 sion ratio
.55:: 1 throat
10.31
5.67
exit
76.0
44.7
mm
200
125
2. EXPERIMENTAL DATA
53
54
.. I
of NO cannot be completed by the time the g a s has flowed through the nozzle.
It was assumed in the calculations that the vibrational degrees of freedom
a r e permanently half-excited (the Lighthill model) / 2 3 / .
The experimental data presented in Figure 4 implies that the s t a t e of the
g a s flowing out of the nozzle not only differs noticeably f r o m the equilibrium
state, but does not correspond t o the conditions of frozen flow. Such a
p r o c e s s of nonequilibrium gas expansion, in which physicochemical
transformations proceed a t a finite r a t e , is described by a s y s t e m of
equations including the kinetic equations of the chemical reactions 11, 8 1 .
The calculations involved a r e very cumbersome and the data on the reaction
r a t e s is unreliable; an attempt has therefore been made t o simplify the
description of the given process.
The simplification reduces basically t o the following:
1 ) if the temperature is not too high, only one chemical reaction i s
important f r o m the point of view of the energy balance, namely, the
recombination of the atomic oxygen through t r i p l e collisions (the amount of
atomic nitrogen and NO is negligible);
2 ) the nonequilibrium effects a r e only connected with the recombination
of oxygen. A s t o the vibrational degrees of freedom, they a r e either in
equilibrium with t h e translational and rotational degrees of freedom or a r e
considered permanently half excited ;
3 ) by analogy with the r e s u l t s /8, 24-27 / for simple diatomic dissociating
g a s e s , the instantaneous freezing model of oxygen recombination is
considered sufficiently close t o the r e a l p r o c e s s . The possibilities and
conditions of application of the instantaneous freezing model of one reaction
t o the processes taking place in gaseous mixtures, and in p a r t i c u l a r in a i r ,
a r e examined in detail in 1 2 5 , 2 8 1 .
F o r temperatures up to 6000" a i r approximately s a t i s f i e s the above-
described model of physicochemical transformations at rapid expansions.
It is possible in t h i s c a s e t o calculate the variation in the p a r a m e t e r s of
the a i r expanding in the supersonic nozzle, assuming a n overall equilibrium
flow up to s o m e section of the nozzle, a f t e r which the flow proceeds
independently of any energy variation due t o the chemical activity of the
components.
The r e s u l t s of calculating the Mach numbers M by this approximate
method a r e shown in Figure 4 (curves 3, 4, 5 ) . They correspond t o the
interruption of the oxygen recombination, originating a t a section of
dimensionless a r e a 2 . 7 , 5.4, and 8.3.
A comparison of the obtained experimental data with these curves leads
t o the conclusion that a s h a r p slowing down of the oxygen recombination
occurs a t the nozzle section of a r e a equal t o 3-7 t i m e s the throat a r e a .
In our opinion the best representation of the experimental r e s u l t s is curve 4.
55
middle curve of Figure 4, corresponding t o the instantaneous freezing
model of the phys ic oc he mica1 t r a n s formations .
Table 2 presents the measured values of the static piessure given in
the published papers, together withthe calculated values for a nonequilibrium
flow based on measurements of the Mach number M. The expansion ratio,
relative t o which all the data of Table 2 r e f e r , has been s o selected a s t o
enable a comparison t o be made with the r e s u l t s of 131. Under t h e a i r
efflux conditions created by u s and the expansion r a t i o selected for
comparison, the p r e s s u r e t u r n s out t o be lower than the equilibrium
p r e s s u r e by a factor of 2 - 2 . 5 (in the c a s e of complete freezing it is lower
than the equilibrium p r e s s u r e by almost a factor of 4). With due allowance
f o r the inaccuracy in determining the freezing section from tlie data on RI,
our result can h e considered a s agreeing satisfactorily with that in 1 3 1 .
The computed data of 1 3 , 4 / , given in Table 3, have been measured f r o m
the curves given in these papers, and therefore the third figure in the
values given i s only tentative. The experimental values have an e r r o r of
the o r d e r of k 1 0 70. The values of the p r e s s u r e derived from measurements
of the Mach number 14 a r e less reliable, since the instantaneous freezing
section is determined with a considerable e r r o r . Thus, foi. freezing in the
section a t which A = 2.72, the s t a t i c p r e s s u r e in the section corresponding
t o A = 144 will be 30 % lower than the value indicated in Table 2 : nevetheless
the result a g r e e s well with the values obtained by direct measurements of
the s t a t i c p r e s s u r e .
TABLF 2 .
ior a fruzcii
To. R P,. atiii che iiiical
for c qu i I i br i ui i
sourct of thechpcriiiiciiral colilposlrloil
expa inion ixpc riimtiit
results from tlit start
(calculation)
of thc cxpaiision
(calcul~lt~o~iJ
* The results of the calculations for T,, =7000 and 4000K arc takvn from /4/, where it was assumed that
the chemical composition was inbaiiablc a n d thr excitation of the vibrations \\as in cquilibriuni. The
result for T , = 5700% presupposes that thc vibration5 are likewise frozen.
56
B ib l i o g r a phy
li
12. S i c h e 1, M. Leading Edge of a Shock-Induced Boundary L a y e r . -
P h y s . of Fluids, Vol. 5. 1962.
13. B u r k e, A. and K. B i r d. P r i m e n i konicheskikh i profilirovannykh
s o p e l v giperzvukovykh ustanovkakh (Use of Conic and Profiled
Nozzles in Hypersonic Installations). -Sbornik "Sovremennaya
tekhnika aerodinamicheskikh issledovanii p r i giperzvukovykh
skorostyakh, ' I Moskva, Izdatel'stvo "Mashinostroenie. '' 1965. t
14. H a y e s , W.D. and R.F. P r o b s t e i n . Hypersonic-Flow Theory.-
Acad. Press, New York. 1959.
15. V i d a v s k i i , A . , R. O s v a l d t , and J. H a r n . Eksperimental'noe
opredelenie konstanty rekombinatsii vodoroda (Experimental
Determination of t h e Hydrogen Recombination Constants). -
Raketnaya Tekhnika i Kosmonavtika, No. 12. 1962. *
16. S t u p o c h e n k o , E.V. and I . P . S t a k h a n o v . Oburavneniyakh
relaksationnoi gidrodinamiki (Equations of Relaxation Hydro
dynamics).-Doklady AN SSSR, No. 4, p. 134. 1960.
17. S t a k h a n o v , I . P . and E.V. S t u p o c h e n k o . O s t r u k t u r e linii
Makha v relaksiruyushchikh sredakh (Structure of t h e Mach Lines
in Relaxing Media).-Doklady AN SSSR, No. 5, p. 134. 1960.
18. W o o d, W. W. and F. R. P a r k e r. Structure of a Centered Rarefaction
Wave in a Relaxing Gas.-Phys. of Fluids, Vol. 1. 1958.
19. M o o r e , F.K. and W.E. G i b s o n . Propagation of Weak Disturbances
in a Gas Subject to Relaxation Effects.-J. Aerospace Sci.,
Vol. 27. 1960.
20. B r o e r, L. J. F. C h a r a c t e r i s t i c s of the Equations of Motion of a
Reacting Gas.-J. Fluid Mech., Vol. 4. 1958.
21. W i l s o n , J. An Experiment t o Measure t h e Recombination Rate of
Oxygen.-J. Fluid Mech., Vol. 15. 1963.
22. W e g e n e r , P. P. Supersonic Nozzle Flow with a Reacting Gas
Mudure.-Phys. of Fluids, Vol. 2. 1959.
23. L i g h t h i 11, M. Dynamics of a Dissociating Gas: Part I. Equilibrium
Flow.-J. Fluid Mech., Vol. 2. 1957; Part 11. Quasiequilibrium
T r a n s f e r Theory.-J. Fluid Mech., Vol. 8. 1960.
24. B r a y , K. N. C. Atomic Recombination in a Hypersonic Wind-Tunnel
Nozzle.-J. Fluid Mech., Vol. 6. 1959.
25. B r a y , K. N. C. Sudden-Freezing Analysis for Nonequal Nozzle Flows.-
A m e r . Rocket SOC.J., Vol. 31. 1961.
26. H a l l , J. G. and A. L. R u s s 0 . Proceedings of t h e First Conference
on Kinetics, Equilibria, and Performance of High T e m p e r a t u r e
Systems. -Butterworth Scientific Publications, London. 1960.
27 L o r d i, J. Comparison of -act and Approximate Solutions f o r Non-
Equilibrium Flows in Nozzles.- Raketnaya Tekhnika, No. 8.
1962. *
28. B r a y , K. Ninth Symposium (International) on Combustion.
N. Y. - London, Academic Press. 1963.
5437 58
I . M . Naboko, R . G . Nemkov
T =S / U -TI,
(1 1
59
the n e a r e s t generatrix of the half-wedge was measured before the experi
ment. Together with a knowledge of the distance between the a p e r t u r e s
a, p, 7 (which s e r v e d t o obtain reference l i n e s ) this made it possible t o
r e c o r d the r e a l s c a l e of the phenomenon on the film. The Mach number of
the flow was determined f r o m the relations
ST
I
"I
ti
60
Figure 3 shows a typical f r a m e - s c a n photograph of the flow past a
half-wedge (shock wave in COP,Mo=7.45 is Mach number ofthe shockwave).
The Mach line is indicated on Figure 3 by the l e t t e r M .
The photograph indicates the vibrations of the Mach line in the COz flow,
which point t o an inhomogeneity of the other flow p a r a m e t e r s . Such
vibrations a r e absent in an N2 flow.
~ ~~
61
The r e s u l t s a r e given in the table, where M corresponds t o s o m e
intermediate p a r t of the plug, in which the flow p a r a m e t e r s are almost
constant with t i m e . When calculating c , all the internal d e g r e e s of freedom
except the rotational ones were considered frozen, i. e., 7 = (c,/cI,) = 1.4.
CONCLTJSIONS
I
Bibliography
62
I
6. B a z h e n o v a , T. V. Izmenenie s k o r o s t i potoka gaza z a s k a c h k o m v
udarnoi t r u b e (Variation in the Gas Flow Velocity behind a Shock
in a Shock Tube).- Sbornik "Fizicheskaya gazodinamika i
teploobmen, ' I Moskva, Izdatel'stvo AN SSSR. 1961.
7. Z a i t s e v , S . G . , E. V . L a z a r e v a , and A.P. S h a t i l o v . Issledovanie
normal'nogo otrazheniya udarnykh voln v udarnoi t r u b e (Study of
t h e Normal Reflection of Shock Waves in a Shock Tube).-PMTF,
No. 4. 1964.
8. B u s y g i n , E.P. and G . K . T u m a k a e v . Izmerenie plotnosti gaza za
udarnoi volnoi v udarnoi t r u b e metodom elektronnogo lucha
(Measurement of the Gas Density behind a Shock Wave in a Shock
Tube Using an Electron Beam). - Zhurnal Tekhnicheskoi Fiziki,
Vol. 34, No. 1. 1964.
63
R . G . Nemkov
It was shown in 111 that in a relaxing gas the Mach line occupies a whole
region (the cone A O B in the figure). It was a l s o established that in the
c a s e s (x,y)+ 0 and (5,y) -> w the s t r u c t u r e of the disturbance region
(emanating f r o m a s o u r c e at x = y = 0 ) is such that a clearly expressed
disturbance maximum exists and is localized along the frozen Mach line
(cf. the line OA in Figure 1 ) and correspondingly along the equilibrium
Mach line (the line OB).
Thus, in these two limiting c a s e s the introduction of the concept of an
"effective" Mach line a s s o m e one-dimensional region, given by
An attempt is made in this paper t o find the profile of the effective Mach
line (as the curve describing the propagation of the wave packet center), i. e.,
64
the function (l), over the e n t i r e range of values of z. The s a m e assumptions
a r e used a s in 111, namely, the relaxing substance (gas) is described by
the Kneser model, and the Mach line corresponds t o such s m a l l disturbances
that it can be described by some linear acoustic equation. The effective
Mach line in dispersive media is considered.
65
Applying the convolution theorem t o the functions C-~(O) and O ( O ) , and
using the relation
T o the best of our knowledge equation ( 6 )has not been used in gas-dynamic
studies. Similar acoustic equations have however been used t o d e s c r i b e
viscoelastic media (according t o the Maxwell model) in 17, 8 1 .
+
Multiplying ( 6 ) by the operator [I i3 ( a / & ) ] , we obtain the t h i r d - o r d e r
partial differential equation
66
2. EQUATIONS O F THE EFFECTIVE CHARACTERISTICS
IN DISPERSIVE MEDIA
OA y = x/T/M&-i,
OB y=x/l/M~--l.
where O A (cf. the figure) is the frozen Mach line and OB the equilibrium line.
This fact can a l s o be derived intuitively. When x =: 0 the disturbance
a) (5,y), i. e. , (1/2n) (a)exp ( i [-~k (a)y
])do,is localized in a s m a l l spatial
region, and in this wave packet a l l the frequencies participate with equal
intensity (i. e., 0 (a)= const). Each frequency then propagates with its
region.
a = - agist,
k = VQ,
61
By analogy, the effective c h a r a c t e r i s t i c s of equation (8) a r e the solutions
of
-
is satisfied along them, where V, is determined f r o m
C Y = v(@M),
Rec(aM) = mas R e c ( o )
0 I
68
3 . PROFILE O F THE EFFECTIVE MACH LINE
IN A RELAXING GAS
We shall discuss the relation obtained. The fact that y (z)= 0 when z < 0
implies that the only liquid particles which interact directly with the s o u r c e
a r e those situated (for x < 0) on the line y = 0.
It follows further f r o m (25) that in the limiting c a s e s (z+ 0, Z > 0 , Z - 0 0 )
the direction of the Mach line coincides (asymptotically) with the direction
of the frozen ( 9 ) and equilibrium (10) Mach lines. The Mach line itself
tends asymptotically t o the line
69
As T 40 we obtain t h e equilibrium Mach line, and as T -+ 00 t h e frozen
line is derived. The limit of y i (T) a s r -+ 0,
i/W-j, x > o ,
lim yS,(r)=
2*0 x=o
Bibliography
71
I .
a plane shock front propagating along the tube channel and the corresponding
equations f o r the detached compression shock around the blunt body under
flight conditions.
The bending of the shock wave discontinuity around the blunt body is
neglected when establishing the correspondence between t h e stagnation
enthalpies of the gas in the shock tube (in the plug) and under flight conditions.
Thus, if u2 is the gas velocity in the plug for a laboratory coordinate
f r a m e , the stagnation enthalpy in the region before the stagnation point is
given by
V2
h' = 2
2
+ hz,
where h, is the enthalpy of the plug gas.
On the other hand, for a coordinate f r a m e linked with the shock front,
v? v? u?
h'=1+"--1
2 2 2 '
72
I , I I I II I1
In the absence of data on the gas velocity in the plug, this velocity can
be eliminated with the aid of the mass-conservation equation on the front,
namely,
Vf ,l ~ /sec
n
FIGURE 1. Correspondence between the stream
parameters in the shock tube and the supersonic-
flight conditions a t different altitudes with respect
to the stagnation parameters:
Shaded regions-simulation conditions; A -typical
trajectory of an earth satellite; B -gliding vehicle;
c-ballistic vehicle; v-maximum altitude of the
equilibrium boundary layer.
Hence
where psrV, a r e the density and velocity of the g a s in the vicinity of the
stagnation point behind the compression shock under flight conditions.
The model of supersonic plug-gas flow can be simulated by
PIUS =P J S .
73
Then
or
which, with the aid of 1 2 1 o r 151, shows that even for M, > 1 0 we can take
or
On the strength of (6) and (9), Figure 1 shows curves of altitude, initial
shock-tube p r e s s u r e , flight velocity and M, for the shock wave in the tube.
Several typical flying-vehicle t r a j e c t o r i e s a r e represented 17 1, in additional
t o the equilibrium boundary l a y e r n e a r the s u r f a c e of a body of radius
R =0.3m.
74
The operating range of the shock tube easily includes the t r a j e c t o r i e s of
different hypersonic vehicles at entry into the dense l a y e r s of the atmosphere.
2. EXPERIMENTAL TECHNIQUES
75
The method of measuring heat exchange by film t r a n s d u c e r s provides a
means f o r detecting minute variations in the s u r f a c e t e m p e r a t u r e . The
oscillogram in Figure 4 clearly shows the vertical jump of the film t e m p e r
a t u r e a t the moment the shock wave touches the model. The analysis of
/ 4 / and the many oscillograms f o r different r e g i m e s show that the t i m e
constant of the film t r a n s d u c e r does not exceed l p s e c .
FIGURE 4. Oscillogram of the pulse from a heat transducer near the stagnation p i n t of a
blunt body:
M s = 6.6; P, = atm; marker frequency 10 psec.
3. EXPERIMENTAL RESULTS
- to1 h,
[t (0,
=, 2.0.5642 w.
Figure 5 shows the dependence of the heat flux on the Mach number M,
of the shock wave and on the flight velocity, calculated f r o m the relations
governing the simulation conditions. Rose and Starks r e s u l t s a r e a l s o
presented; they w e r e obtained by means of c a l o r i m e t e r s on models of
12.7 d i a m e t e r 111. It is seen that the film t r a n s d u c e r s provide a means
for studying the heat exchange at lower M, and Re.
Up t o flight velocities V ,= 6 k m / s e c the experimental r e s u l t s coincide
satisfactorily with the theory of Fay and Riddell / l o / for the equilibrium
boundary layer. The calculation w a s realized in 141.
16
To obtain quantitative data on the heat flux f o r the leading p a r t of a
flying vehicle a correction must be introduced t o allow a correction for the
body geometry. Then
where r,,, is the model radius; R , i s the radius of the leading p a r t of the body.
qsI km/cm2 Y 6 8
8
e
2
I
0.8
6?t
0.2
D.1
008
0.09
2.0 20 g.0 $0
Vf. h/sec
FIGURE 5. Heat flux near the stagnation point of the model as
a function of rhe shock-wave Mach number M:
1- p,=1.3.10-2 atm. H = 22--28km; 2- ~ s = l O - ~ a u n ,
H e 35-4Okm; 3 - r e s U l t s of /l/; 4-the author's results.
Bibliography
77
3. O f f e n h a r t z , E . , H. W e i s b l a t t , and R. F l a g g . Stagnation-Point
H e a t - T r a n s f e r Measurements at Super-satellite Speed.-J. Roy.
Aeronaut. SOC.,Vo1.66, No. 1 . 1962.
4. P o l y a k o v , Yu. A. P r i m e n e n i e plenochnykh datchikov dlya izucheniya
teploobmena v dissotsiirovannom potoke gaza (Use of F i l m
T r a n s d u c e r s f o r Studying Heat Exchange in a Dissociated Gas
Flow). - Sbornik "Fizicheskaya gazodinamika, teploobmen i
termodinamika gazov vysokikh t e m p e r a t u r . I ' Moskva, Izdatel'stvo
AN SSSR. 1962.
5. R o z h d e s t v e n s k i i, I. B. Termodinamicheskie i gazodinamicheskie
svoistva potoka vozdukha z a pryamym skachkom uplotneniya s
uchetom d i s s o t s i a t s i i i ionizatsii vozdukha (Thermodynamic and
Gas-Dynamic P r o p e r t i e s of the Air Flow behind a Normal
Compression Shock allowing f o r Dissociation and Ionization of
t h e Air).-Sbornik "Fizicheskaya gazodinamika, ' I Moskva,
Izdatel'stvo AN SSSR. 1959.
6. 0 s t o s 1a v s k i i, I. V. and I. V. S t r a z h e v a. Dinamika poleta.
T r a e k t o r i i letatel'nykh apparatov (Flight Dynamics. T r a j e c t o r i e s
of Flying Vehicles).- Moskva, Oborongiz. 1963..
7. C h e n g , H. K. Recent Advances in Hypersonic Flow Research.-AIAA
Journal, Vol. 1, No. 2. 1963.
a. P o l y a k o v , Yu.A., I. M. N a b o k o , and Yu. V. M a k a r o v . Eksperi
mental'noe opredelenie rabochego v r e m e n i udarnoi truby metodom
teplovogo zonda (Experimental Determination of the Operating
T i m e of a Shock Tube Using a Heat Probe).-Teplofizika vysokikh
t e m p e r a t u r , No. 3. 1965.
9. L y k o v , A. V. Teoriya teploprovodnosti (Theory of Heat Conduction).-
Moskva, GITTL. 1952.
10. F a y , J. A. and F. R. R i d d e 11. Theory of Stagnation Point Heat
T r a n s f e r in Dissociated Air.-J. Aero. Sci., Vol. 25. 1958.
78
A . A . Kon'kov, A . P . Ryazin, V . S . Rudneu
1. EXPERIMENTAL SET-UP
19
2. SPECTRAL COMPOSITION O F THE RADIATION
AND DEPENDENCE O F THE ABSORPTION
COEFFICIENT O F AIR ON WAVELENGTH
80
using the reference lines of iron, the wavelengths of all the lines in the a i r
s p e c t r a w e r e interpreted and determined. The impurity lines in the a i r
radiation s p e c t r a were highly broadened, and it was consequently difficult
t o accurately determine t h e i r wavelength. W e therefore obtained /7/ the
radiation s p e c t r a of krypton and xenon at 20000-~5000"K, where the lines
of the impurities only appear in the absorption spectrum. The r e v e r s e d
lines a r e quite narrow, and the accuracy in finding t h e i r wavelength is
much higher than determining it f r o m the radiation s p e c t r a . If the mutual
spacing of the lines was the s a m e in the s p e c t r a of a i r ,
krypton, and xenon, and the wavelength coincided with
that of the impurity line, it was considered a s
established that the given line belongs t o the impurity.
The spectrograms were processed using hetero
chromatic photometry 16 /. The graphic variant of this
method was used: using an MF-4 microphotometer,
microphotograms of the s p e c t r a of the SI-8-200 tungsten
band lamp w e r e recorded on the s a m e plate with the
aid of a nine-step attenuator on the s a m e film a s the
spectrum being studied. The dependence of the
emission power of tungsten on the wavelength is well
known, and s o these microphotograms enable t r a n s
formations t o be made f r o m the blackening in the
spectrum of the gas studied t o the dependence of t h e
intensity in relative units on the wavelength.
This variant of the method f r e e s u s from the cumber
s o m e and time-consuming construction of the c h a r a c
t e r i s t i c curves of the film for every wavelength. As
indicated above, the relative dependence of the intensity
on the wavelength was normalized t o absolute values
with the aid of data f r o m a p h o t o e l e c t r i c channel. The
s p e c t r u m zone of width 1 0 A about the wavelength4850A
was selected f o r the normalization, because it is f r e e
of impurity lines, a s seen f r o m the s p e c t r o g r a m
(Figure 2). The emissive power was thus obtained a s
a function of wavelength. The dependence of the
absorption coefficient on the wavelength was derived
using the laws of Kirchhoff and Beer.
This method of finding the absorption coefficients
presupposes a known gas temperature. Under the
given experimental conditions, t h e r e a r e insufficient
grounds for taking the a i r temperature equal t o the
value calculated on the basis of the conservation laws;
accordingly, the accuracy of such a calculation was
corroborated by direct measurements of the gas
temperature, whose r e s u l t s a r e presented below.
Figure 3 (curve 1)gives the experimental relation
ship between the a i r absorption coefficient and the
wavelength for a t e m p e r a t u r e of 11500K and a p r e s s u r e
of 60 atm. The a i r absorption coefficients, obtained
FIGURE 2. Time-re-
solved air radiation
by directly measuring the s p e c t r o g r a m blackenings,
spectrum; the temper- a r e given here. The vertical lines give t h e values of
81
height of the lines corresponds t o the value of the absorption coefficient
measured directly a t the given point of the spectrum. It follows f r o m an
examination of both this curve and the s p e c t r o g r a m (cf. Figure 2 ) that v e r y
few zones a r e f r e e of impurity lines.
Kh.102. cm-1
U
3000 3500 'to00 5000 i
FIGURE 3. The air absorption coefficient as a function of the wavelength:
1-experimental data of the present paper, temperature 11500'K, pressure 60 arm; P-theore
tical values according to f2, 3 / , temperature 12000oK, pressure 50 a n n ; 3-theoretical values
according IO /4/9 temperature 12000"K, pressure 50 atm.
82
Figure 4 (curve 1 ) shows the experimental relationship between the a i r
absorption coefficient and the wavelength at t e m p e r a t u r e 10500K and
p r e s s u r e 50 atm. The impurity lines have been excluded f r o m consideration
h e r e and the given values of the absorption coefficient only correspond t o a i r .
The accuracy of the absorption coefficients is determined by the accuracies
of the heterochromatic photometry and the normalizing factor. The accuracy
of the heterochromatic photometry is determined by the graduation accuracy
of the standard source, and the measurement accuracy of blackening by a
microphotometer, and f o r m s *20 70. The measurement accuracy of the
normalizing factor is determined by the accuracy with which the standard
s o u r c e is s e t up, its graduation, the measurement of the potential by an
oscillograph, and the determination of the gas p a r a m e t e r s , and f o r m s rt 25 70.
83
L ________ -____ Z -I L..---l
W
P
< "K
I
a
1-
I4000
I3000
12000
ll000
IO000
I;
1
21 25 27 29N
FIGURE 8. Radiation intensity FIGURE 9. Temperature of the air heated
(a) and absorptive power (b) as a by a reflected shock wave as a function of
function of time. the Mach number of the incident shock wave.
85
can be calculated f r o m the data of the second one. After measuring the
r a t i o of sensitivities of the channel, one can, if necessary, suitably c o r r e c t
the values obtained f r o m the f i r s t channel.
The a i r temperature in the section situated 4 m m fromo the end of the
shock tubeb was measured. A s p e c t r u m zone of width 5 A was selected
n e a r 4238 A f o r measuring the temperature. The s p e c t r o g r a m (Figure 2 )
and the dependence of the absorption coefficient on the wavelength (Figure 3 )
show that the given s p e c t r u m zone is f r e e f r o m impurity lines and corresponds
t o the absorption maximum (0.1) of the f i r s t negative band of nitrogen.
The radiation intensity ( a ) and the absorptive power (b) have been plotted
in Figure 8 a s afunction of time. The curves show that after the a r r i v a l of
the reflected shock wave the radiation remains constant for a p e r i o d of
lG-30psec (at times considerably longer), then drops rapidly t o a relatively
s m a l l value, and remains at that level for a considerable time. The absorption
is constant for 10-30psec after the a r r i v a l ofthe reflected shock wave, and
t h e n s t a r t s t o i n c r e a s e t o unity. The s a m e shielding effect of the radiation
ofthe auxiliary s o u r c e bythe gas being studied after the reflected shockwave
and the contact surface meet, was observed e a r l i e r in 161. No reliable
interpretation of this effect has s o far been given.
The experimental values of the a i r temperature a r e represented by dots
in Figure 9. The initial a i r p r e s s u r e in the low-pressure chamber was
5 m m mercury. The accuracy of the temperature measurement was 10000". *
The continuous curve r e p r e s e n t s the dependence of the temperature of the
a i r heated by the reflected shock wave on the Mach number of the incident
wave. These data a r e taken f r o m /8 /, where the temperature was calculated
f r o m the conservation laws. The theoretical values of the temperature a r e
givenin / 8 / for Mach numbers up t o 24; the values corresponding to higher
Mach numbers w e r e obtained by extrapolation.
The theoretical a i r temperature is seen t o coincide with the experimental
value up to Mach number 23. F o r higher Mach numbers a systematic
discrepancy is observed; however, taking into account the accuracy of the
temperature measurement and the fact that the theoretical values in this
region a r e extrapolated, it is impossible to s t a t e with certainty that a
discrepancy does indeed exist.
86
effect of impurities, whose emission can make an additional contribution,
which makes it n e c e s s a r y t o consider t h i s value of the oscillator strength
of the second positive s y s t e m of nitrogen a s the upper limit. The above-
mentioned papers imply that t h e value n e a r e s t the t r u e value is 0.039
obtained in Ill/, where the lifetime of the excited s t a t e of the second
positive s y s t e m of nitrogen was measured. The oscillator strength of the
f i r s t negative s y s t e m of nitrogen - was taken a s 0.035 in 141, according- t o the
data of 7 1/. The adiabatic approximation is assumed c o r r e c t . The authors
of 1 2 1 assumed the s a m e value f o r
Re! at. units t h e oscillator strength of t h e f i r s t
negative s y s t e m of nitrogen, but
assumed in addition that the square
of the matrix element of the dipole
transition moment Re'depends on
the internuclear distance: the
assumed relationship is represented
by t h e dotted line in Figure 10.
Hence, the values 0.032-0.035
and 0.039 a r e the best approximations
of the oscillator strength of the f i r s t
4 000 a,a" negarive and second positive nitrogen
FIGLIREi n . Dcpcndcnce of thi. s w a r t ' of the matrix s y s t e m s respectively. Since no
~ l c t i i ~ of
n t the dipole 1 1 1 0 r i i c n 1 of Vlectron transition radiation increase was observed in
of the first negative systcm of nirrogcn on the wave the region of 3500A in the present
length .
study, the relationship represented
by t h e solid line in Figure 1 0 w a s
adopted a s expressing the s q u a r e of the matrix element of the f i r s t negative
s y s t e m of nitrogen a s a function of the wavelength.
It was therefore found advisable t o introduce the following changes in the
computed values of / 2 / and 141:
1 ) to exciude the contribution of the negative nitrogen atom f r o m the data
of 141;
2 ) t o take 0.039 a s the oscillator strength f o r the second positive s y s t e m
of nitrogen;
3 ) t o a s s u m e that the s q u a r e of the matrix element is given by the solid
curve in Figure 1 0 .
87
FIGURE 12. Dependence of the absorption coefficient of air on the wavelength:
CONCLUSION
B i b l i o g r a p hy
1. A r m s t r o n g , B . H . , I . S o k o l o f f , R . W . N i c h o l l s , D.H. H o l l a n d ,
and R. E. M e y e r o t t . Radioactive P r o p e r t i e s of High Temper
a t u r e Air.- J. Quant. Spectr. Radiative T r a n s f e r , Vol. 1,
No. 1/2. 1961.
88
89
Yu.A. Plustinin
SPECTROSCOPIC INVESTIGATIONS OF THE
EMISSION OF HOT NITROGEN
1. INSTRUMENTATION
5437 90
bore. The low p r e s s u r e compartment terminated in a s t a i n l e s s s t e e l s q u a r e
section 35 X35 mm, connected by an intermediate section with a compartment
of c i r c u l a r c r o s s section. The radiation was observed both along and
a c r o s s the tube axis f r o m corundum a p e r t u r e s at the end of t h e shock tube.
L a t e r a l a p e r t u r e s were used in most of the experiments; the observation
axis was 4 m m f r o m the end surface.
The d r i v e r gas was helium under a p r e s s u r e of 200atm, and s o the
intermediate compartment was filled with helium up t o 1 atm. The shock
tube was equipped with a s y s t e m t o purify the gas being studied (removal
of solid foreign particles, drying the gas, freezing out the C O z ) . The low
p r e s s u r e was p 1 = 25 m m mercury. The velocity of the normal shock wave,
*
measured by means of ionization s e n s o r s , was v = 3450 30 m / s e c . The
given p i and u corresponded t o a design t e m p e r a t u r e and p r e s s u r e behind
the reflected shock wave of T3 = 7100 f 30" K and pR= 34 f 0.4 atm. The gas
p a r a m e t e r s and composition of hot nitrogen w e r e determined by the usual
method 111, using the thermodynamic functions of nitrogen in 171. The
molar fraction of nitrogen was z=O.821. It follows f r o m / 8 / that the
design t e m p e r a t u r e of nitrogen -7000K coincides with the measured value;
the time-consuming pyrometric investigations of the working conditions of
our apparatus could therefore be dispensed with, On the strength of the
data in / 9 / we established that under our experimental conditions the
relaxation p r o c e s s e s behind the reflected shock wave ceased a f t e r 6 0 . l p s e c .
The equipment for the s p e c t r a l measurements consisted of a g l a s s
spectrograph ISP-51 with c a m e r a of f = 2 7 0 mm, an image converter with
multislit photocathode / 6 /, a f o u r - c h p n e l diffraction s p e c t r o m e i e r with
relative a p e r t u r e 1:7, dispersion 2 6 A / m m t o first o r d e r and 1 3 A / m m t o
second, and a condenser of f = 90 mm. The e n t r a c e slit of the s p e c t r o
m e t e r was 0.1 m m wide, and the exit slit 1m m wide. The radiation was
recorded in a l l four channels by means of FEU-38 photomultipliers and
OK17-M pulse oscillographs. The receiving-amplifier circuit was graduated
by means of a standard tungsten SI 8-300 lamp.
The beam of light f r o m the l a t e r a l a p e r t u r e was split by a s e m i t r a n s
parent m i r r o r into two p a r t s , one of which was directed t o the entrance slit
of the s p e c t r o m e t e r , and the other t o the entrance slit of the spectrograph.
With this apparatus, the radiation f r o m a givenvolume of gas couldbe studied
simultaneously by photoelectric and photographic methods. A n image-convert
e r photocathode was placed where the spectrum image formed i n the focal plane
of the spectrograph chamber. The spectrum was photographed f r o m the image-
converter s c r e e n o n highly sensitive f i l m . The spectrograph slit was 0.02 m m
wide. The image-converter had a delayed f r a m e - s c a n c i r c u i t , triggered by
a signal f r o m an ionization s e n s o r at the end of the shock tube.
2. MEASUREMENT RESULTS
91
The following conclusions w e r e reached. Atomic lines of nitrogen, the
band s p e c t r u m of the nitrogen molecules, and atomic lines of the impurities
a r e a l l excited in the spectrum. The following atomic lines w e r e identified:
Element Wavelength, A
N . . . 8216, 8185, 7468, 7442, 6723, 6645
0 . .. 8446, 7773
0
One notes the presence in the s p e c t r u m of the intense oxygen line 7773A,
due to the u s e in the experiments of industrial nitrogen, which contains
0.3-0.8 70 0.
Wavelength, a
7400 7600 78Q0 8#00 8ZQU 84Q0 86N
v i - I ,---T~-rl
I
a b
FIGURE 1. Frame scan of the nitrogen emission spectrum (1) for wavelengths 7 3 0 0 - 8 5 0 0 i ,
-
3
D
A f = 25 psec,
and microphorogram of the emission spectrum (b) ( the second frame):
1- theoretical curve.
Wavelength, a
68K ~ 7000 7200 7400 76?0 7800
Sequence d o - 2
a b
FIGURE 2. Frame scan of the nitrogen emission spectrum (a) for wavelengths 6800--7800R, A t = 80 wec,
and microphotogram of the emission spectrum (b) (the second frame):
1-theoretical curve.
92
Figures 1- 3 indicate that the number of impurity lines and t h e i r
intensity increase sharply with time. Thus, t r a c e s of K lines can be s e e n
on the s p e c t r o g r a m of the f i r s t f r a m e in Figure 1, while in subsequent
f r a m e s t h e i r intensity increases sharply. The lines of the other easily
excitable element Ca appear much l a t e r in the s p e c t r u m than the K lines.
Unlike the K and Ca lines, the intensity of the 0 lines is almost constant
from the t i m e the spectrum i s first photographed. The differences in
behavior of the radiation intensity of the K, Ca, and 0 lines a r e probably
due to the fact that the K and Ca atoms originate at the walls of t h e shock
tube, while oxygen is present in the gaseous phase throughout the radiating
gas. Additional observations of the K, Ca, and 0 lines from the end
a p e r t u r e corroborate this opinion. The emission s p e c t r o g r a m s established
that during the f i r s t 25psec a few weak impurity lines a r e present.
Sequence d u = 3
FIGURE 3. Frame scan of the nitrogen eniission spectrum (a) for wavelengths G 4 0 0 - 7 2 0 3 i , A f 80 vsec,
7
93
F i g u r e s I- 3 present microphotograms of emission s p e c t r a and t h e
r e s u l t s of calculating the blackening density, allowing f o r the s p e c t r a l
sensitivity of the converter radiation r e c e i v e r . The s p e c t r a l absorption index
of the f i r s t positive s y s t e m of nitrogen w a s calculated by a computer using
the formula /11/
where
- . . ..
a. K Intensity Intensity
__ .~ . - -_ _ . . _ _ _
1
~
5450
GOOll
4.16; 3.G
20.8; 23.5; 17.2;
16.G
30.6: 33.0 3R.0;
15.7; 1 7200
73i3
18.1; 19.3; 25.2; 26.9
l R . O : 27.2: 27.8: 20.6
94
The values of the absolute intensities obtained by the photoelectric method
a r e given in the table ( [Ihl=olw/cm3. p * s t e r a d ) . The s p e c t r a l range isolated
by the spectrometer was 2 6 A f o r wavelengths >
0 . 6 4 ~and 1 3 A for wave
lengths < 0.64~.
95
The formula holds f o r an optically thin gas l a y e r . To t e s t this assump
tion, the intensity was measured during recordings of the radiation a t the
end of the shock tube. The corresponding radiation oscillogram is shown
in F i g u r e 5. A l i n e a r f o r m is observed on the oscillogram for 25psec, a
fact indicating that the radiation intensity and t h e thickness of the radiating
gas a r e related linearly 113, 141.
The r e s u l t s of determining the s q u a r e of the m a t r i x element over the
wavelength range 5450-85OOA a r e given in Figure 6, together with values
ifRe'obtained f r o m the relative
distribution of the radiation intensity
when processing the microphotograms
of Figures 1 and 3. A relative relation
ship was normed with the aid 2f data
on Re'for the wavelength 7490 A, obtained
by the photoelectric method. It is seen
62
B i b 1i o g r aphy
96
5. K e c k , J . C . , R.A. A l l e n , a n d R . L . T a y l o r . E l e c t r o n i c T r a n s i t i o n
M o m e n t s f o r A i r Molecules.- J. Quant. S p e c t r . Radiative
T r a n s f e r , Vol. 3, No. 4. 1963.
6. B u t s 1o v , M. M. E l e k t r o n n o - o p t i c h e s k i e p r e o b r a z o v a t e l i d l y a
izucheniya s v e r k h b y s t r y k h p r o t s e s s o v (Image C o n v e r t e r s f o r
Studying U l t r a r a p i d P r o c e s s e s ) . - U s p e k h i n a u c h n o i fotografii.
Vol. 6, T r u d y S o v e s h c h a n i y a PO v y s o k o s k o r o s t n o i f o t o g r a f i i i
k i n e m a t o g r a f i i , M o s k v a - L e n i n g r a d , I z d a t e l ' s t v o AN SSR. 1959.
7 . G l u s h k o , V. P . , e d i t o r . T e r m o d i n a m i c h e s k i e s v o i s t v a individual'nykh
v e s h c h e s t v ( T h e r m o d y n a m i c P r o p e r t i e s of Individual S u b s t a n c e s ) . -
M o s k v a , I z d a t e l ' s t v o AN SSSR. 1962.
8. K u d r y a v t s e v , E . M . , N.N. S o b o l e v , L . N . T u n i t s k i i , and
F. S. F a u z i l o v . P i r o m e t r i c h e s k o e i s s l e d o v a n i e s o s t o y a n i y a
g a z a z a o t r a z h e n n o i u d a r n o i v o l n o i ( P y r o m e t r i c Study of t h e Gas
State h e h i n d a Reflected Shock Wave).-Trudy FIAN, Vol. 18.
1962.
9. L o g a n , D. G. R e l a x a t i o n P h e n o m e n a i n E y p e r s o n i c A e r o d y n a m i c s ,
C o l l e c t i o n of t r a n s l a t i o n s . - Mekhanika, No. 1. 1 9 59.
1 0 . A t l a s s p e k t r a r t u t i ( A t l a s of t h e M e r c u r y S p e c t r u m ) . - A l m a - A t a ,
I z d a t e l ' s t v o AN KazSSR. 1959.
11. P 1a s t i n i n , Yu. A. a n d G. G. B a u 1 a. S e c h e n i y a p o g l o s h c h e n i y a
e l e k t r o n n y k h s i s t e m d v u k h a t o m n y k h m o l e k u l N,, O,, N,, NO,
C,, CN. F i z i c h e s k i e j z m e r e n i y a v g a z o v o i d i n a m i k e ( A b s o r p t i o n
C r o s s S e c t i o n s of E l e c t r o n i c S y s t e m s of t h e D i a t o m i c M o l e c u l e s
N,, 0,. N,, NO, C z , CN. P h y s i c a l M e a s u r e m e n t s in G a s
D y n a m i c s ) . - C o l l e c t i o n of p a p e r s e d i t e d by A. S. P r e d v o d i t e l e v
" I s s l e d o v a n i e PO f i z i c h e s k o i g a z o v o i d i n a m i k e , " M o s k v a ,
I z d a t e l ' s t v o "Nauka. ' I 1966.
12. Z a r e , R.N., E . O . L a r s o w , a n d R . A . B e r g . F r a n c k - C o n d o n
F a c t o r s for E l e c t r o n i c Rand S y s t e m s of M o l e c u l a r N i t r o g e n . -
J . M o l e c u l . S p e c t r o s c . , Vol. 15, No. 2. 1965.
13. D i e c k e , G. H. and D. F. H e a t h. T h e First a n d S e c o n d P o s i t i v e
R a n d s of N,.-Johns Hopkins S p e c t r o s c o p . Rept. No. 1 7 . 1959.
1 4 . K e c k , J . C . , J . C . C a m m , B. K i v e l , a n d J. W e n t i n g . R a d i a t i o n
f r o m Hot A i r . - Ann. P h y s . , Vol. 7, No. 1. 1 9 5 9 .
15. T u r n e r , R.G. a n d R . W . N i c h o l l s . An E x p e r i m e n t a l Study of
Rand I n t e n s i t i e s in t h e F i r s t P o s i t i v e S y s t e m IT,. T h e T r a n s i t i o n
in t h e F i r s t Moment,-Canad. J . P h y s i c s , V o l . 32, No. 6. 1954.
16. N i c h o 11 s, R. W. F r a n c k - C o n d o n Factors t o High V i b r a t i o n a l
Q u a n t u m N u m b e r s N, a n d N,+.- J. R e s . Nat. B u r e a u S t a n d a r d s ,
Vol. 6 5 A . No. 5 . 1961.
97
V . A . Bashilov, V . P . Shevelev
the whole cone; 4--5psec l a t e r a glow appears in the region behind the
lOpsec it fills
98
behind the annular electrode approximately lOpsec after the discharge
commences. Near the tube walls a narrow, pale ejection appears, which
gradually shrinks t o the center of the tube, forming a cone with a glowing
vertex. The oblique shocks converging t o the tube center gradually become
distinct, and t h e general pattern of the plasma flow f r o m the conic chamber
during this period resembles the pattern of supersonic flow f r o m a nozzle.
The angle at the vertex of the luminous cone increases, the vertex itself
becomes blurred, and after -12 p s e c the regular geometrical configuration
disappears. After another 5 p s e c , the pattern is r e s t o r e d . This restoration
was found t o be linked with the second half-period of the discharge c u r r e n t .
F I G W 1. Rotating iiiirror photograph of the process in the discharge chamber. The axis of the camera
mirror lies a t right angles to the tube axis.
99
F o r higher initial tube p r e s s u r e s the plasma efflux f r o m the cone is not
accompanied by the formation of a distinct and stable pattern, and t h e r e
a r e no oblique shocks a t the exit f r o m the cone. AS the plasma moves
f u r t h e r along the tube, the leading front of the plasma glow c e a s e s t o have
a distinct boundary.
The investigation of the plasma flow pattern at l a r g e distances f r o m the
discharge chamber showed that for low initial p r e s s u r e s (0.1-0.5 m m
m e r c u r y ) of the d r i v e r gas the leading front of the plasma remains plane
up t o a distance of 1.5 m f r o m the discharge chamber (no investigations
were conducted at l a r g e r distances). F o r higher initial p r e s s u r e s the
leading front of the plasma flow did not have a distinct shape a t l a r g e
distances f r o m the annular electrode, and the formation of "tongues" and
diffuse ejections f r o m the brightly glowing region was observed.
V, kin /sec
3l
0,!un/sec
Zff
IO 0
F: mm mercury u v , k"
FIGURE 2. Glow front velocity as a function FIGURE 3. Glow front velocity as a
of the initial pressure in the shock tube: function of the voltage on the capac
itor battery, for air.
1-hydrogen: 2-air; 3 -argon.
IO0
FIGURE 4. Rotating mirror photograph of the propagation process of the glow front in the top part of the
shock tube over a portion enclosed by light filters, for p = 7 m m mercury.
The fluid was dyed t o aid studies of the front s t r u c t u r e . Sodium chloride
was introduced at the central electrode of the discharge chamber, and the I
0 2 4 E 8
of the shock tube f o r different initial
P,m m mercury discharge conditions (Figure 4). At initial
FIGURE 5. Distance between the glow p r e s s u r e s of 0.1-0.5 m m m e r c u r y and
front and the boundary of the sodium- voltages of 3.5-5 kv applied t o the
colored region, as a function of the discharge chamber, t h e r e was no differ
initial pressure of a i r ence in the recorded glow whether o r not
101
I I I I I Il 1 I I I
Some papers, such a s 171, established that in discharge shock tubes the
t i m e of flow past a model considerably exceeds the computed value f o r a
shock wave with corresponding velocity. This discrepancy is normal, since
the existence of two zones has been established in a plasma flow, namely,
that of the shock-heated gas and that of the gas-discharge plasma, the l a t t e r
exceeding the f o r m e r in its dimensions by a factor of 10-15.
The flow past differently shaped bodies was studied using the IAB-451
Schlieren unit with an S F R c a m e r a taking s e p a r a t e films. Satisfactory
Schlieren photographs of the flow of ionized gas in discharge tubes a r e
difficult t o obtain, in view of the strong plasma glow, To obtain a s a t i s
factory Schlieren photograph, the optical s y s t e m between the exit s l i t of
the Schlieren unit and t h e entrance slit of the SFR c a m e r a was adjusted,
The two-line framing inset together with the conversion s y s t e m produced
an image of the model in a s c a l e such that a working field of vision of 13 c m
diameter was represented on a f r a m e of 10" diameter. To that end RO-2
and h d u s t a r - 2 3 objectives with focal distances of 75 and 1 1 0 m m respectively
were placed between the shadow apparatus and t h e SFR c a m e r a . Two knife
edges, each 0.4" wide, were inserted at right angles t o the flow in the
IAB-451 instrument t o produce a slit of height 8 mm. These adjustments
enabled a satisfactory Schlieren photograph of the flow past the models t o
be obtained.
102
The window in the plexiglass section was made of polished optical g l a s s
of type K8. The g l a s s was washed a f t e r every experiment, and the film
of copper and soot removed.
The light s o u r c e was an IFK-300 pulse lamp with a bank of capacitors
of 200pf capacity and 3 kv voltage. To initiate breakdown, a specially made
s p a r k generator with delay circuit was used in the lamp. The g e n e r a t o r
was triggered by a pulse f r o m the SFR panel.
I03
Bibliography
I04
V .A . Bashilov
with the obviously overestimated value of 20000K for the temperature (at
a front velocity of 20 k m / s e c ) , we obtain for the hydrogen atom
A h D = 0.25 A:
3 ) shock broadening, appearing only at very high p r e s s u r e s .
These effects can be neglected for a halfwidth of HP- 40 A.
The s p e c t r a obtained on the spectrograph for different initial discharge
conditions were processed a s follows. The continuous s p e c t r u m of a
tungsten l a m p was photographed on the s a m e film through a nine-step
attenuator. The characteristic curves of the film were then plotted for
t h r e e points with wavelengths 6560, 4860, and 4 3 4 0 A (all the s p e c t r a
processing was conducted on an M F - 4 microphotometer). The value of
for each s p e c t r u m range was determined on the basis of these curves. The
H, and Hp lines of the s p e c t r u m under observation w e r e then measured
photometrically.
I05
L!%=
JO
0.5: Ig k
Jo
= Ig 0.5 = 0.301; Ig Jo,5= Ig J o - 0.301.
$2
8
6-
4- 3or
2-
S, cm
FIGURE 1. Electron concentration behind the shock FIGURE 2. The temperature behind the shock
front and in the gas-discharge plasma: front in hydrogen for p 0 = 4 m m mercury:
1-region of shock-heated gas; 2-region of gas- o -experimental points obtained on the basis
discharge plasma. of Sachs equation.
106
allowing f o r t h e electrons
Bibliography
107
T . V .Bazhenova, Yu.S . Lobastou
I08
This is probably linked with the fact that a distortion to the reflected wave
front only exists in a narrow region near the front 161. The increase in
radio-wave absorption during 6Opsec is due t o the establishment of
equilibrium ionization of nitrogen.
_._
__IC1
f,
109
reaction N 4- 0 t NO' 4-e (cf. Figure 3, curve 4), whose r a t e constant equals
3.10-ST-S/2cm-3. It was shown in / 7 / that the ionization t i m e of a i r for
T = 3000--5000K is mainly determined by the formation t i m e of the nitrogen
and oxygen atoms. The high values obtained in t h i s paper f o r the ionization
t i m e of nitrogen indicate that not only the formation of nitrogen atoms, but
the ionization p r o c e s s itself also represent a limiting process.
Two possible ionization reactions f o r nitrogen a r e analyzed in / 9 / :
2,,pet. a t m
!
I
68 MO
110
111
G . P. P e t r o v a . E . K . Chekalin
I12
The electrodes w e r e thoroughly insulated f r o m the whole apparatus and
cooled by distilled water t o prevent shunting of the electrodes in the
direction of the water cooling. The layout of the plasmatron with channel
and electrodes is shown in Figure 2. The screening casing with a 4 m m
thick w a l l w a s made of Armco iron, and prevented the magnetic field f r o m
influencing the plasmatron a r c c u r r e n t . The s a m e figure shows the layout
of the electrocies and poles of the electromagnet (view f r o m the channel end).
The electromagnet poles had a c r o s s section of 14x16 mm', and the
maximum s t r e n g t h of t h e magnetic field w a s 4500gauss. The magnetic
field s t r e n g t h in the gap was measured by a coil of 1.5 m m radius with
11 = 50 t u r n s . The field in the gap w a s found t o be homogeneous t o an
accuracy of 1 0 % .
FIGURE 1. Layout of the plasmatron with quartz chamber and exit cone:
FIGURE 2. Layout of the plasmatron vith quartz channel and electrode unit:
1-water-cooled arc chamber; 2-tungsten cathode; 3-nozzle anode; 4
screening casing; 5-quartz channel; 6- water-cooled quartz jacket; 7-brass
electrodes; 8 -insulating talc -chloride washers; 9 -poles of the electromagnet;
10-water; 11-argon; 12-distilled water.
113
20 30 4o %IC cham, '
2. E X P F R I M E N T A L RESULTS
The emf induced by the flowing plasma and the current in the electrodes
circuit w e r e both investigated for various magnetic fields and different
working conditions of the plasmatron.
The values of the emf were measured by an M-193 voltmeter and an
ENO-1 oscillograph with an input impedance of 0.5 megaohm. In the absence
of a magnetic field the plasma flow induces a potential difference on the
electrodes varying f r o m experiment t o experiment, and probably due t o
the nonuniform combustion of the plasmatron a r c or a disturbance in the
coaxial alignment of the nozzle and electrodes.
Tests of two types of electrode showed that the s c a t t e r of the initial
potential difference is maximum for thin tungsten electrodes. This is
probably due t o the fact that when electrodes a r e introduced in a f r e e jet
of relatively s m a l l diameter, even s m a l l deviations f r o m the s y m m e t r y
can cause changes in the temperatures and velocities of the probes (in a
jet of such s m a l l diameter the temperature and velocity gradients a r e
considerable). In addition, during a prolonged operation the wire probes
became very hot. F o r l a r g e water-cooled electrodes in the channel, the
value of the initial potential difference was much s m a l l e r and m o r e stable
than in the f i r s t variant. These electrodes w e r e capable of operating for
long periods (all the experiments w e r e conducted with the s a m e electrode
unit ).
Oscillograph m e a s u r e m e n t s established that the voltage on the electrodes
has a variable component, the amplitude of whose oscillations is of the s a m e
o r d e r of magnitude as the emf; the frequency of these oscillations is
I14
I
v;v
-----
I I
D I 2 3 t,sec
FIGURE 6. Typical oscillogram of the current in the circuit of cooled electrodes
115
III I l l l l l
emf. v a b
The emf and the c u r r e n t on the electrodes have been plotted in Figure 7
as functions of the magnetic field strength f o r two lengths of the quartz
channel, namely 58 and 38 mm. The emf and the c u r r e n t in the electrodes
c i r c u i t a r e s e e n t o d e c r e a s e f o r increasing channel length.
1 I6
In fact, the longer the channel, the s t r o n g e r the cooling and the lower
the average temperature of the plasma in the region between the electrodes.
F o r two channels of lengths lI and Lz respectively w e have
Bibliography
1 I7
E . K. Chekalin, V . S. Shumanov
1 I9
FIGURE 2. Radiation of the plasma of an exploding copper wire of 0.3mm dian~ etcrand l O O m m length for fro=3 0 k v and p o = 3. lO-'mm mercury:
a-motion-picture record along the wire axis: b-motion-picture record on one sidc of the wire; c-rotatins mirror photograp11 011 one sidc of the wire; d -rotating mirror
photograph of the emission spectrum.
the narrow horizontal slit. These moving glow zones a r e probably produced
by f l o w s of electrode material, whose formation is not completely c l e a r .
It is mentioned .in / 4 / that the j e t s of electrode m a t e r i a l a r e possibly due
to the sputtering of the electrodes m a t e r i a l in the liquid phase. The
actual sputtering is analogous to the p r o c e s s of the e l e c t r i c explosion
of metallic wires. The continuous s p e c t r u m a f t e r 18 and 32psec is
t h e r e f o r e apparently of the s a m e nature as the continuous s p e c t r u m at the
initial stage of t h e explosion.
The existence of an absorption s p e c t r u m indicates that t h e t e m p e r a t u r e
corresponding t o t h e glow intensity of the continuous s p e c t r u m is higher
than t h e electron t e m p e r a t u r e in the outer l a y e r s of the expanding plasma
cloud. No absorption s p e c t r u m was observed in t h e r e d region.
Thus, both the integrated radiation of the plasma of an exploding wire and i t s
s p e c t r a l distribution are inhomogeneous in time and space.
The appearance of electrode jets is probably t h e r e a s o n for the intense
continuous spectrum, against whose background is observed the absorption
s p e c t r u m of the lines due t o the metal in the wire.
These facts complicate somewhat t h e utilization of exploding w i r e s as
a radiation s o u r c e f o r practical applications.
B i b l i og r a phy
121
E . K . Chekalin, V . S . Shumanov
e
0.24 P ( t )R ( t )dt = M d T , (2 1
5437 122
I
By solving equations (1) and ( 2 ) simultaneously one can find the variation
with t i m e of the discharge current and the wire r e s i s t a n c e ; this a s s u m e s
a knowledge of the dependence of the wire resistance on temperature, which
up t o the instant of melting approximately satisfies the relationship
R ( T )= & ( I + PT), (3 1
1
R = ROrxp(At3), T = - [exp ( A t 3 ) - l ) ] ,
P
where
I23
. ... .. . . .. .
TABLE 1.
-
Material
aluminum
,. . * $ :- c
l i-
*I- r L-
.4-a-
* :*; :% ,1 4%
FIGURE 1. Oscillogram of the current I and voltage U during the explosion of a copper wire of 0.1 m m
diameter and 100mm length:
0. = 30 kv: time markings 0.1 w e c .
The variation with time of the discharge current in the intial stage of
the process is linear. The well-known current pause is absent in these
experiments. Only at the instant of wire explosion does one s e e a s m a l l
variation in d l l d t on the c u r r e n t oscillogram. Such a variation in the r a t e
of i n c r e a s e of the current is to be expected, since at the explosion time
considerable energy is transmitted to the wire f r o m the bank of
capacitors, leading t o a rapid r i s e in the temperature of the wire. This
I24
is accompanied by a s h a r p increase in the w i r e r e s i s t a n c e , and therefore
a d e c r e a s e in d l l d t , t h e r a t e of increase of the c u r r e n t .
L a t e r , a f t e r the breakdown of the s p a r k gap and the ignition of the a r c
discharge in the metal vapors, the c u r r e n t variation is mainly determined
by the e l e c t r i c p a r a m e t e r s of the circuit, since the resistance and inductance
of the discharge interval a r e s m a l l compared with t h e corresponding
p a r a m e t e r s of the discharge c i r c u i t .
The initial voltage drop on the w i r e is determined by the ratios of the
resistances and inductances of the w i r e and the whole circuit. The w i r e
resistance l a t e r increases with increasing temperature. The increase of
the w i r e r e s i s t a n c e , and the increase of the discharge current with time,
lead t o an i n c r e a s e in i t s voltage potential. The f i r s t point tl a t which the
curve slope v a r i e s rapidly corresponds t o the transition f r o m a solid t o a
liquid s t a t e , since the material cocductivity v a r i e s discontinuously during
the phase transition. At s o m e instant the w i r e is p a r t solid and p a r t liquid.
The second voltage discontinuity is observed at t i m e f2, the beginning of
the evaporation, a f t e r which the w i r e resistance i n c r e a s e s until the voltage
gradient along it becomes sufficiently high t o ignite an a r c in the metal
vapors a t time t 3 . Thus, f o r a copper wire of 0.31 m m diameter the
oscillogram yields /,,, ( t l - t u ) = 1.70psec, while the theoretical value is
t,,, _v 1.49psec (cf. the table). The s m a l l increase in the experimentally
obtained melting t i m e of the wire is due t o the lower discharge current at
the melting t i m e a s compared with i t s theoretical value.
FIGURE 2. Oscillogram of the current I and the voltage U during the explosion of a copper wire of
0.31mm diameter and 100" length:
U O = 30 kv; t i m e markings 0.2 psec.
I25
One can derive data on the electrical conductivity of metals for different
phase s t a t e s of the m a t e r i a l by analyzing the electric c h a r a c t e r i s t i c s of the
explosion. A study can a l s o be made of the thermophysical p r o p e r t i e s of
metals in liquid and vapor phases a t high t e m p e r a t u r e s and p r e s s u r e s .
126
..
.. .. - . .. . . ,
V . F . Vygovskii, E. K. Chekalin, V . S . Shumamv
DEFORMATION O F SOLIDS BY A N
EXPLODING WIRE
3 4
7 b
a
FIGURE 1. Experimental apparatus:
a -electric blockdiagram of the apparatus: 1-rectifier; 2-charge resistance; 3 -bank ofcapacitors;
b-layout before the explosion: 1-copper pipe: 2-plasticine: 3- wire; 4-steel ring: 5-forward
conductor
127
where a:., din a r e the internal d i a m e t e r s of the pipe after and before the
wire explosion.
Very strong pipe-ring joints w e r e obtained in the experiments.
y kj
Bibliography
128
A . S . Predvoditelev
dT
c,p x d v .
Suppose the burning r a t e is $(T,z ) and the specific heat liberated during
combustion is q ; the quantity of heat liberated in t h e burning zone will then
be given by
s Y 9 do.
1N w h do.
129
where (jlis the latent heat of vaporization and +,
is the r a t e of evaporation.
The heat liberated by the external heat s o u r c e is
or
where the b a r indicates the mean value taken over the burning zone.
Designate P by P,. We can then always alter s o m e mean values t o satisfy
t h e equalities
130
Taking these relations into account, the lasc equality can be written a s
follows :
I1
or
-
2. Rapid heating of the material leads t o the formation and motion of a
burning zone. The value of P , w i l l then satisfy inequalities of the f o r m
F P 8 > =;
3 LGPs>q.
Hence
p - AT
_ = _A-T ='c AT
-=--: A T An 1
'- Az r Az Trp An r r p '
P . = r gQ = . I (8 1
Q
V
-=a=ggt.
S
Therefore
AT
Q = TAT = I,-
An
Equality (10) now takes the form
$= (T)/(z) J
- _-
aq A
g: = a l j
p .
MzQ
g1 =
('Pi& + 891
41 F )EN
algF
QaT
Hence, formula (11)can be rewritten in the simple form
I33
I11
g=y+
y 3 - l2
y3 - 2 (2+-$ ) =o.
ys + 3py + 24 = 0.
Therefore, the following identities a r e valid in our c a s e :
alternatively
25 +$ >o.
i. (13a)
I 34
I
With the aid of this equality g1 takes t h e form
Suppose the heat source for the ignition yields exactly the amount of heat
n e c e s s a r y for the evaporation of the drop. Then
v l S ( T - T d )N V QS,
V
- = d = gt = g,
ill
.
S rl
Therefore
The factor before the brackets in this formula is identical with the first
t e r m of ( l o a ) , which is equal t o g; if AT is equal t o the difference ( 7 ' - T d ) .
These differences a r e not in fact equal in the general case. Jntroduce y,
defined by
or
I35
or
136
A . S . Predvoditelev
INTRODUCTION
where y i s the specific heat of the chemical reaction, and second by the
heat t r a n s f e r r e d f r o m the walls t o the g a s :
cps ( T , - T,)
9
Equation (1)and (2) do not f o r m a closed system, and do not give the
condition f o r the ignition of a volume of gas mixture. During ignition, as
mentioned above, a direct link must exist between the t e m p e r a t u r e and
the mixture composition at any point in the volume, i. e., a relationship of
the f o r m
T = f(z).
The t h r e e equations (1)- ( 3 ) must hold simultaneously f r o m the ignition
time .
138
When the f i r s t equation i s differentiated with respect to the t e m p e r a t u r e
T , and t h e second one with r e s p e c t t o the variahle
z, where
x e obtain
For slow heating, relation (4a) is valid at any instant. T o find the
ignition condition, the g a s p r e s s u r e P at the ignition t i m e must be known.
Consider s e v e r a l p a r t i c u l a r c a s e s . Suppose P s a t i s f i e s the conditions
139
chemical reaction. Bearing this in mind, we can w r i t e relation (4a) a s
The mean values of P,c,,p and a$faT a r e selected t o satisfy the equalities
dT 1
g = - --, dT dT
d t z gx=-Z*
dn
F r o m F o u r i e r ' s equation we obtain
and therefor e
where h r e p r e s e n t s h,/pc,.
Let AT be the t e m p e r a t u r e difference between the s u r f a c e s of the burning
zone. Hence
=VAT.
V Sld.
Relation (11)then becomes
141
A fundamental r e s u l t is now drawn f r o m this expression. Consider the
following relations f o r t h e burning zone:
d=gz, M = l i d t = % ,
0
- .-
2az
g= aTF- (12a)
A T 7 M
I42
lead under suitable temperature conditions t o the development of the
chemical process. The quantity PZ is the partial p r e s s u r e of the gas minus
the concentration of the active center. The following equality w i l l then hold:
P = P1 + Pz.
When the gas is not ignited its p r e s s u r e p fluctuates with t i m e about
s o m e mean value.
The Volterra equations can be written in the f o r m
_
dPl -
d l - PI (a1 - P I P ~ ) , = -Pa (az- ~ L P ~ ) .
Pz
dP1
a't + PI
dP2
= a1Pzp1 - azP1pz.
Now multiply the f i r s t equation by zz/pl and add the r e s u l t s . This yields
o r alternatively
z = E!..
P
The Volterra equations, expressed in t e r m s of z, will now have t h e f o r m
143
(which may be the mean molecular weight of any g a s mixture), and t h e
notation
d M = [ a l m l z -aazm4(1 - z ) ]p
dt +[ m & ( l - z)z - m & (1-z)]p2=
=F1(T, 4P+F2(T1,4PZ.
The following additional notation is introduced:
dT=u=o.
dM
5 u)
(I++ [(%-+ PFl) p + (S+ 2PFz) P q = 1.
When condition ( 3 ) is observed, it follows that
[(S
144
whole c o u r s e of our reasoning, leading t o the existence of two ignition limits
with r e s p e c t t o the p r e s s u r e , shows that this phenomenon is wholly
determined by second-order reactions. Therefore, whatever model of
chemical transformations we attempt t o construct, given the two ignition
l i m i t s in t e r m s of p r e s s u r e , the leading reaction must always be of second
order.
Following Semenov it can be shown that t h e r e exist ignition regions
defined by c u r v e s having t h r e e roots with r e s p e c t t o t h e p r e s s u r e .
According t o our reasoning, a t h i r d - o r d e r leading reaction must exist in
this case.
w i l l then hold for this zone. The mean values a r e chosen t o satisfy the
identities
This relation holds for the g a s ignition, i. e., for the formation about
the heat s o u r c e of a combustion front propagating throughout the g a s volume.
It is a l s o valid when a burning zone, which is incapable of propagating
throughout t h e gas volume, is formed about the heat s o u r c e . The l a t t e r
can be realized for very s m a l l values of P,. Since this magnitude (accord
ing t o equality (17))depends on the p r e s s u r e , it i s obviously possible t o
attain minimum values of P,. The displacement r a t e of the burning zone
I45
w i l l then be very low, and the whole g a s volume w i l l only ignite after a
l a r g e time interval. In fact, P , can be expressed by
-
VSAT = Q = -q -
d x
d t Sd = - 4 ( F l p + F2p2)Sd,
I p TSAT 9 (&P + F2PZ) .
d 8 - SdAM - AM
F1 ( T , Z ) =- + ( ~ z m+s ) = -b
~ 1 m lz
+ biz,
F,(T, z ) = (P2mz- (1 - 2 ) z = a ' ( l -
P1ml) 2) 2.
[- b + b l A M ) + b,;
"P = - a' [- A M - (1- 2 A M ) z-+ 91 '
146
The roots of the quadratic expression in t h i s equation a r e positive for
all values of A M , however small. We can therefore w r i t e
b
a*
- . b l L
a*
AM - f >0.
Formula ( 2 2 ) was corroborated experimentally.
147
G . D . Salamandra
148
The gas motion ahead of the flame front was calculated by the method
of c h a r a c t e r i s t i c s . The construction of the c h a r a c t e r i s t i c s for a plane
isentropic gas flow in a bounded space was t r e a t e d in detail in 14, 51.
The method of c h a r a c t e r i s t i c s was used in / 6 / t o calculate the gas s t a t e
ahead of the flame front, propagating in long tubes, when the combustion
is accompanied by the formation ahead
. of it of a shock wave o r a transition
d b into detonation. This method gave
r e s u l t s agreeing satisfactorily with
experimental r e s u l t s .
F o r combustion in s h o r t tubes, the
calculation is complicated somewhat
by the necessity of allowing f o r the
influence of the waves reflected
f r o m both the tube end and the
flame front. Figure 2 r e p r e s e n t s
a part of the wave diagram of the
p r o c e s s shown in Figure 1 , a . The
a b s c i s s a r e p r e s e n t s the distance x,
and the ordinate the t i m e t . The
flame is represented by the bold line.
The position of the walls is
indicated by vertical lines with
coordinates x = 0 and r = 1 6 7 mm. The
f i g u r e s I, 11, IIIdesignate the artificially
created inhomogeneities. The ( 5 ,t )
plane is covered by a grid of c h a r a c
t e r i s t i c s , constructed a s follows.
The slope of the characteristic
0-0, starting f r o m point 0 is
determined by the sound velocity in
the undisturbedgas, equalto 3 9 4 m / s e c .
This characteristic i n t e r s e c t s the wall
at Oo. The region 0, O,,x corresponds
t o the state of r e s t . To construct the
c h a r a c t e r i s t i c starting f r o m point 1
it is n e c e s s a r y t o know the gas velocity
u and the sound velocity a at this point.
The gas velocity ahead of the flame
front is determined f r o m the t i m e
r e c o r d of the process 161. The sound
velocity was determined f r o m the
c h a r a c t e r i s t i c in the ( u , a)-plane
(Figure 3). Due t o a suitable selection
of the s c a l e , this c h a r a c t e r i s t i c
FIGURE 1. Combustion o f a hydrogen-oxygen r e p r e s e n t s a straight line inclined at
mixture containing 67 qo 0, and 33 q n HI:
I
an angle of 45" t o the coordinate axes.
a-continuous time scan of the process: b
Knowing the value of u and a at 1, we
series of instantaneous photographs. Rate of draw a straight c h a r a c t e r i s t i c s t a r t i n g
photography 200000 f r a meshec; 1- flame; f r o m point 1 in the ( 5 , t )-plane (cf.
2-artificially created inhomogeneity. F i g u r e 2). Its equation is dxldt = u f a .
149
Through 0, we draw a straight c h a r a c t e r i s t i c ofthe second family dxidt = -a,,
where a, is the sound velocity in f r e s h a i r . The intersection of the two
straight c h a r a c t e r i s t i c s determines the point O1. The characteristic
s t a r t i n g f r o m 2 is constructed similarly. The point 0,lies at the i n t e r s e c
tion of the c h a r a c t e r i s t i c s starting f r o m points 2 and 0,. The slope of the
c h a r a c t e r i s t i c O,,O, is determined by the sound velocity at 0,, equal to the
sound velocity at 1. To find the values of u and a, necessary f o r construct
ing the c h a r a c t e r i s t i c 0,, 1, we draw the straight line u - u1 = - (2Iy - 1) x
X(a--a,) through point 1 in the ( u, a)-plane. This straight line, making an
angle of 90" with the characteristic 1, ... , 7, intersects the ordinate axis a t
l , , whose coordinates a r e the required values. The slope of the c h a r a c t e r
istic 0, , l , is determined by the coordinates of the point of intersection
of the c h a r a c t e r i s t i c s in the ( u, a )-plane, drawn f r o m points 1, and 0,.
The remaining c h a r a c t e r i s t i c s a r e plotted in the s a m e manner.
150
is given in Figure 4, where the a b s c i s s a axis r e p r e s e n t s t h e distance
traveled by the flame f r o m the ignition point, and the ordinate axis
r e p r e s e n t s the gas velocity. The vertical line indicates the flame
position at the instant under consideration. The horizontal line is
drawn at a level corresponding t o the flame velocity at the given time.
The points indicate the experimentally obtained gas velocities, with
which the calculation a g r e e s satisfactorily, Note t h e coincidence of
the oscillation phases at a l l sections along the tube; this is c h a r a c t e r
istic for a standing wave in a medium at rest.
The compression waves propagating in the burning medium affect
the shape of the flame front by altering the quantity of g a s burning
in it. Although the question of
2 the interaction between the
- a , m/sec
Y-1 c o m p r e s s i o n waves and t h e flame
front is of great i n t e r e s t ,
sufficient data does not yet exist
on the variation of the flame
front shape during vibrational
combustion.
It is known that Schlieren
photography visualizes regions
in which the first derivative of
the refraction index differs f r o m
zero. Since the gradient of the
refraction index is maximum in
the zone of heating, this zone
is always recorded in t h e
S c N i e r e n photographs. F r o m its
variation one obtains an idea
about the variation .in the flame
front shape (assuming that during
interaction with the acoustic waves
the burning zone is deformed in
the s a m e manner as the heating
zone). The photographs in Figure 5
(on p. 1 5 5 ) show that this is not so.
According to the f i l t e r a r r a n g e m e n t
in the lattice used in the experi
ments, the undisturbed field of the
Schlieren photograph is r e d [ g r a y
background], the region of s m a l l
gradients of the refraction index is
green [fringes], and the region of
l a r g e gradients is blue [band
between fringes]. When inter
preting color Schlieren photographs
one must b e a r in mind that due
t o the finite width of t h e l a t t i c e
0 20 40
s t r i p s every color includes
U,m/sec
certain angles of deviation. T h e
FIGURE 3. Characteristics in the ( Y, a 1-plane. wider t h e lattice s t r i p s , the
151
l a r g e r the region of deviations c o r
responding t o the s a m e color. W e
u s e d a l a t t i c e ofstripwidth0.5 mm.
The Schlieren photograph
shows that t h e development of
the p r o c e s s is accompanied not
only by a variation in t h e total
width of the zone recorded on
the photograph, but a l s o by a
redistribution of t h e regions with
different refraction-index gradients
inside this zone. The leading
boundary of t h e heating zone, which is
habitually used f o r outlining the shape
of the flame front, v a r i e s much m o r e
Distance, m m
slowlythanthe combustion zone. F o r
FIGURE 4. Distribution of the gas velocities between a detailed study of the flame front
the flame fmnt and the end of the tube at different s t r u c t u r e t h e n u m b e r of filters in
instances of time: the lattice must be increased.
1-3*10-4sec; 275*10-4;3--7*10-4; 4- 8.10-4; The following conclusions can be
5--10-'. drawn f r o m here.
1. For the vibrational propaga
tion of a flame in a tube, the gas ahead of the flame front is s e t in motion
by the disturbances generated by the front. The problem of the gas
motion ahead of the flame front can be considered as one-dimensional.
The calculations of the gas s t a t e ahead of the flame front by the method
of c h a r a c t e r i s t i c s a g r e e satisfactorily with experimental r e s u l t s .
2. F o r a m o r e detailed study of the flame front s t r u c t u r e in the c a s e of
vibrational combustion it is n e c e s s a r y t o u s e color Schlieren photographs.
Bibliography
I53
5437 154
FIGUREl. Schlieren
photographs of the
vibrational propaga
tion of the flame:
Rate of photography
about 3000 frames
per sec .
155
Bibliography
I56
N . K . Fedoseeva
157
FIGURE 1. Generdl view of the apparatus forprojecting films
158
Bibliography
159
S . G . Zaitsev, E . V . Lazareva, E . 1. Chebotareva
160
interfering beams). This condition, which is usually fulfilled for the main
s t r e a m , makes the relations given below inapplicable in the immediate
vicinity of a wall. A s a result, the region within 0.2-0.3" of a wall is
inaccessible to measurements.
Consider the c a s e when the direction of the bands coincides with the
direction of the density gradient in the s t r e a m . The interference pattern
then makes it possible t o determine in a z e r o approximation the variation
of the density field, given by
YI
Yt
161
We shall estimate the e r r o r in the density caused by the bending of the
light r a y due t o an inhomogeneity. Its magnitude is determined by the
variation in the refractive index along the r a y trajectory in the volume
studied, and the variation in geometrical length of the ray path AC-AB
(Figure 3).
Let the r a y propagate along the straight line A B in the absence of a
density gradient, in a gaseous medium of refractive index nl. The equation
of this line is
y - yo = ax.
(3 )
where
N
y-yo=-
2 (x+ 7)
fas (7 1
t o an accuracy of t e r m s of o r d e r (y' - d 3 .
The light r a y propagating along AC continues beyond point C along a
straight line coinciding with the tangent t o A C at C ; this line i n t e r s e c t s A B
at the point E . Therefore, in the approximation adopted, the inhomogeneity
considered does not disturb the plane in which the interference pattern is
localized. The angle between .4B and EC = ( W h y ) 1. F o r the regions of
the boundary l a y e r studied, y > 0.2 m m and y < When using interfer
ence fringes of width S equal to -0.2 m m (magnification of the interference
5431 162
.. . . - ....
._ .~ ~
. ___ .. . . .. __ . . .. . .. ..
\
A= n(y)C+mdx-
AC
\
AB
nll/l-+3dx.
AP = P1 (Yo)
where
2-3
Relation (10) enables one t o estimate the
e r r o r when t h e density variation is determined
f r o m the interferograms by relation (1).
In this paper, for y* 31 0 . 2 mm, max el 6 3 70.
Obviously, a s y v a r i e s f r o m y* t o 6 (the boundary
l a y e r thickness), the value of e, drops t o zero.
Consider the boundary-layer effect due t o
A-1/z the cover glass, on the density measurement
_L
163
In this relation, the density on segment CD is taken a s invariable in the
direction of r a y propagation and determined by the distance f r o m the wall.
On segments D B and A C the density v a r i e s due t o the effect of the s i d e walls
y = 2 1 1 2 . T o estimate t h e s e t e r m s we used the density distribution in the
boundary l a y e r obtained by analyzing the interferograms on the b a s i s of (l),
i. e . , disregarding the effect under consideration.
The density variation is represented in the f o r m
where
164
sharp, and the e r r o r in determining its position did not exceed 0.1 mm.
The boundary l a y e r behind the shock front was investigated by
interferograms with fringes oriented parallel (horizontal fringes) and
perpendicular (vertical fringes) t o the upper and lower walls of the chamber.
When vertical fringes w e r e used, the interferograms w e r e recorded on a
fixed film. Thus, each photograph gives an interference pattern of the
p r o c e s s in the measuring chamber at a given instant of time. To study the
variation in the boundary l a y e r along the plug, a number of interferograms
w e r e taken under the s a m e working conditions, but at different t i m e s during
which the light s o u r c e functioned. The reproducibility of the working
conditions was not lower than 0.1 M , and the light s o u r c e was a s p a r k of
duration l 0 p s e c .
When horizontal fringes were used, the interference fringes w e r e oriented
parallel to the upper and lower walls of the chamber. The film moved at
a speed of 0.117 m m / p s e c in the propagation direction of the shock wave,
and the a r e a of measuring chamber projected on it was r e s t r i c t e d by means
of an optical slit oriented parallel to the shock front. The photography thus
gives the variation with time of the interference pattern at a given section.
The light s o u r c e was an I F P - 8 0 0 flash l a m p of flash duration -800psec.
Since the density increases gradually in the boundary layer, the fringes
passing through the boundary layer region a r e inclined at an angle of
inclination which d e c r e a s e s toward the main s t r e a m . When the interfero
g r a m s were processed, we measured the displacement of the fringe center
at a distance y f r o m the wall relative t o the center of this fringe on the
shock front. A s I IL the method with vertical fringes, - the relationship pllp[y)
was established by means of formula ( l ) , and the boundary layer thickness
determined.
165
where 5 is the displacement of the fringe center at a distance y f r o m the
wall relative t o the fringe center in the main s t r e a m ; S is the fringe width,
measured in the main s t r e a m . The measured values of E w e r e c o r r e c t e d
t o allow f o r fringe distortions caused by g l a s s imperfections. The r e s u l t s
were represented in the f o r m of the relationship pl/p(y) = f(y) ( F i g u r e 5).
The thickness of the boundary l a y e r was taken as the distance f r o m the
wall corresponding t o pl/p(y) = 0.99.
166
d mm
167
viscosity in the s t r e a m ; Re6 is the boundary-layer Reynolds number,
determined by
S;mm
d mm
FIGURE 9. Variation of the boundary layer FIGURE 10. Variation of the boundary layer thickness
thickness with t i me, in nitrogen: with time, in nitrogen:
P.= 0 . 1 5 0 ~1 0 - 4 g/cm3,T o = 3009(, 0,=353 m/sec; p 0 = 0.150.10-4 g/cm3, T.=3009(; 1- p 0 = 50 mm
l-.p.= 10 mercury; M,= 5.1; u;= 1470 m/sec; mercury; M,= 2.1; u r = 690m/sec; T , = 690K; p l =
T, = 1670%; p , = 0 . 8 3 1 . 1 0 - ~ g / c m ~ P; ,= 310 mm = 2.70. 10-4g/cm3; p l = 4 1 2 m m mercury; P
mercury; ~.r=0.552.10-g/cm.sec; Reg= L046, lo6; = 0.318.10-~ g/cm.sec; ~ e 6 1= .557.106; 2- P o = 50;
2-~0=10; M,= 5.4; =;=1570; T,=1820; P, = M, = 3.6; u;= 980; Tr = 1010; P.= 3.25. p,=
= 0.863, ~ 1 ~ 3 5 !J=
0 ; 0.598. io-; Re6 = 730; P = 0.406.10-3; Re6=2.67.1OS; 3 - P,= 50;
1.100. lo6. M , = 3.1; u;= 1010; T,=1050; P I = 3.30. PI=
= 770; P = 0.416.10-3; Re6=2.79, lo7; 4- p o = 50;
M s = 4.0; u;= 1120; T,=1160; p t = 3.63.
PI = 930; P = 0 . 4 0 io-; Reg= 3.63. io7; 5 - P o =
= 50; M,= 5.6; u;= 1640; TI=1930; & = 4.38.
R = 1880; P= 0.620. lo-; Reg= 5.75.10.
168
variation pattern of the density in the plug remains the s a m e a s f o r hydrogen
or helium, but the fluctuation i n c r e a s e s . The fluctuation is affected by the
method of igniting the oxygen-hydrogen-helium mixture. If the ignition is
distributed (e. g., when a glow-lamp filament is placed along the whole high-
p r e s s u r e chamber), the fluctuations a r e a minimum.
CONCLUSION
It was shown that for boundary-layer Reynolds numbers below 1.1. lo6
and 0.3. l o 6 for the nitrogen and argon s t r e a m s respectively, the variation
of 6 ( f ) a g r e e s satisfactorily with the growth law of a fully l a m i n a r boundary
layer.
Bib1i o g r a p hy
169
7. M a r t i n , W. A. An Experimental Study of the Turbulent Boundary
L a y e r behind the Initial Shock Wave in a Shock Tube.-J.
Aerospace Sci., Vol. 25, No. 10. 1 9 5 0 .
8. S t e v e r , H., E. W i t m e r , and W. H e r m a n n . Development of the
Boundary Layer behind a Shock Wave. - In: "Boundary Layer and
Heat T r a n s f e r P r o b l e m s , " Gortler, H. and W. Tollmien, editors.
[Russian Translation 1960.1
9. M i r e 1s , H. Attenuation in a Shock Tube due t o Unsteady-boundary
layer Action.- NACA Tech. Note No. 3278. 1 9 5 6 .
170
R . G . Nemkov
T = T (u). (1)
Suppose we wish t o light a flash bulb at the exact instant that a shock
wave p a s s e s through a section c - c' (Figure 1 ) of the shock tube. The
simplest procedure would be t o place a probe a t section c - c'to t r i g g e r
the bulb. However, this is sometimes impossible, e. g . , when t h e t i m e t o
light the bulb (to its maximum luminosity), although stable, is too s m a l l
compared with the t i m e taken f o r the wave t o pass the section. A probe
located u p s t r e a m f r o m section c - c'(cf. Figure 1) is then used, and the
problem is t o delay the pulse generated by it for a t i m e
circuit:
1-control
Trigger pulse
171
Figure 2 r e p r e s e n t s the block diagram of this delay circuit. The delay
T~~ of the phantastron is a l i n e a r function of its controlling voltage (Figure 3),
s o that
U, = c,At,
where cy is constant.
It is s e e n f r o m the block diagram that
UDh(t) = ";hf uY ( t ) ,
I I
+At -i
8
I
8 n n
Q b
-
t
FIGURE 3. The phantastron delay T~ as a FIGURE 4. Diagram of the functioning
function of iu controlling voltage Uph. of the control circuit:
A-output pulse as a function of time:
B-probe pulses.
172
6N6P 6Kh2P TG3-0.1/1.3 6N1P 6N1P 6N6P
e
4
W
Bibliography
174
M . V . Gusev, V .A . Gordyushin, V . D . Lobanov
-
FIGURE 1. Principal schematic diagram of pulse-counting chronometer.
175
Channel I Channel I1 Piezoelectric t r a n s d u c e r s , which
.- w e r e designed in our laboratory 141,
a r e now being used as a pulse-counting
chronometer, and the pulse counter
used h e r e is t h e standard s c a l e r of
the PS-10000type.
The chronometer counts the pulses
with a frequency of 1MHz which a r r i v e
between two pulses f r o m the piezo
e l e c t r i c t r a n s d u c e r s ; it is a two-
channel electronic s h a p e r of the pulses
f r o m the piezoelectric t r a n s d u c e r s
into rectangular pulses, which a r e
equal in length t o the t i m e interval
between the leading edges of the pulses
f r o m these s e n s o r s .
The principal schematic design is
shown in F i g u r e 1. The voltage
diagrams (by blocks) of this design
(Figure 2 ) clarify the operational
procedures of the blocks and t h e i r
functions.
The signals f r o m t h e s e n s o r s
(Figure 2, a ) at the separating inputs
of Channels I and I1 open up the
thyratrons of the shaping blocks
(Figure 2, c ) by t h e i r leading edges
a f t e r amplification in the f i r s t cascade
(Figure 2, b). As a result of thyratron
FIGURE 2. Diagram of the output voltage shaping, the instrument effects a single
from the instrument blocks: operation during the discharge t i m e
a-shape and characteristics of the signals and pileup of the capacitance charge
of the piezoelectric pressure sensors; b C6, at the anodes T, and T,. This time
amplifying stages (T1 and T5 - 6E5P): is regulated by the proper selection of
c-pulse-shaping stages (T1and Ts - RB,SIIand C6, , s o that it is much
TGI-O.1/1.3); d-trigger (T3 and T4 longer than the t i m e for a signal f r o m
6N6P); e-regulated generator (T, and
t h e piezoelectric s e n s o r , i. e . , the
T8-6F1P); f-coincidence circuit (T9
6Zh2P).
t i m e for attenuation of the wave
p r o c e s s e s in the shock tube. On
a r r i v i n g a t the t r i g g e r (Figure 2, d), the signals f r o m the thyratrons
alternately open it up (Channel I ) and cut it off (Channel 11) with t h e i r leading
edges. A rectangular voltage pulse a r i s e s at the t r i g g e r output and opens
up the coincidence circuit (Figure 2, f). At this moment, the coincidence
circuit becomes passable f o r sinusoidal oscillations of frequences 1 MHz
coming f r o m t h e generator. This generator is regulated by a quartz
resonator of the K-17 type (Figure 2, e). Negative pulses of duration
0.4psec and repetition frequency 1MHz a r e obtained at the output of the
coincidence circuit and recorded by the PS-10000 s c a l e r .
176
MEASUREMENT ERRORS
177
the input pulses f r o m the leading edges, the delays in the actuation of the
first channel produce a negative e r r o r , while those of the second channel
produce a positive e r r o r .
The variations in the delay for both channels a r e determined, not only
by the duration t of the leading edges, but a l s o by the r a t i o k between the
output voltage of each stage and the triggering threshold of the subsequent
stage.
Thus, the absolute e r r o r , e. g . , for the shaping stage, is expressed
a s follows:
Sign of
Stage ChI ChII
v,,=1000m/sec I-o = 4000m/sec
+
0.25 0.25 m 1.0
_ _ ~
ATmmax =
+ (0.78S--O.G12) wsec.
- (0.287-0.112)
5437 178
Considering that the electronic tube (Stage f ) is open the e n t i r e t i m e
f r o m t h e beginning of the leading edge t o the end of t h e t r a i l i n g edge of a
t r i g g e r signal, as well as the fact that only t h e negative p a r t of the s i n e -
wave of t h e generator is transmitted in the PS-10000 device, we can expect
a difference between the measurable and measured t i m e of one pulse in the
+ 0.788
AT,,,,,<
1
s c a l e r , with a maximum random e r r o r of AT,,,^^ = -0.287 psec, since
T is the period of t h e oscillations of t h e generator.
Thus, f o r any m e a s u r a b l e time interval, the maximum possible random
e r r o r is f 1 psec. The relative e r r o r 6t in this c a s e is a value connected
l i n e a r l y with t h e measurable t i m e , and is equal t o 1 % for an interval of
1OOpsec.
Bibliography
179
M . V . Gusev, V . D. Lobanov, V . A . Gordyushin
180
kv
0
Lvasf constant of the circuit R,& (cf. Figure 1).
At the multivibrator output we obtain a
rectangular resolving pulse of controllable
width z (Figure 2,c), which is selected
t equal t o the duration of the p r o c e s s
FIGURE 2. Piagram of the output voltage studied. (By delaying the pulse f r o m the
from the instrument blocks: piezoelectric t r a n s d u c e r during a t i m e t
t h e leading edge of the rectangular pulse
a-amplifying stage (TL-6E5P); b-thy
rafron stage for shaping and delaying the
has already been made t o coincide with
pulse (T, and T3-TGZ-0.3/1.3): c-multi t h e beginning of the p r o c e s s phase being
vibrator (T1, 5 -GN1P); d-phase-shifting studied in the experimental section. )
stage (T6-6NGP); e-coincidencecircuit The resolving pulse is fed t o the
(D,, D,. R 3 ; f-thyratron stage of the coincidence circuit D,, D,, Ra6 (Figure l),
control pulse (T7. T,, T,-TGI-O.3/1.3). and when a pulse reaches it simultaneous
ly f r o m the pickup located on the SFR
m i r r o r (Figure 2, d), a pulse of positive voltage appears at the output of
the coincidence circuit (Figure 2, e), which is then used t o t r i g g e r t h e
thyratron block. As a result, a control pulse (Figure 2, f ) was only obtained
when the phase under observation coincided with the correSponding position
of the SFR m i r r o r . The pulse f r o m the pickup of the SFR m i r r o r was taken
f r o m the outlet synchronization on t h e SFR control panel during its functioning
under controlled conditions. To bring t h i s pulse in correspondence with
the resolving pulse of the t r i g g e r e d multivibrator, it passed through the
phase -shifting stage of the circuit (T6lh 6N6P).
After i t s transformation in a high-voltage t r a n s f o r m e r ( T p , ) , the control
pulse triggered the flash bulb of the IAB-451 instrument (bulb IFP-500). The
film in the SFR was exposed while the s o u r c e was luminous.
181
The flash of the IFP-500 bulb was selected t o satisfy t < ( U 2 ) T . since
the SFR m i r r o r was two-way; T is the period of rotation of the m i r r o r .
In the s c h e m e described the beginning of recording does not have t o be
synchronized with the f i r s t f r a m e , since the p r o c e s s studied required
10-15 f r a m e s when using one series of a two-series lens i n s e r t and an
SFR m i r r o r rotation speed of 30,000 r p m . This made it possible to
dispense with additional pulse shaping f r o m the SFR m i r r o r pickups and,
f o r a pulse width of 5Opsec, t o s t a r t recording the p r o c e s s f r o m any one
of the first six f r a m e s (for the above rotation speed).
The flash shows whether the image of the p r o c e s s was incident on the
film, and when n e c e s s a r y t o repeat the experiment without having t o
develop the film. This was impossible previously.
Bibliography
182
V .M . Eroshenko, A . V .Moray, Yu.N . T e r e n t ' e v
183
discontinuity s u r f a c e s in the medium s t a t e on which t h e initial discontinuities
a r e conserved during the e n t i r e development t i m e of the p r o c e s s . The front
s u r f a c e is c h a r a c t e r i s t i c f o r s o m e wave propagation, and its motion is
described by a Hamilton-Jacobi partial differential equation. The solution
is based on Huyghen's idea that a t any instant of t i m e the wave front
propagating f r o m a s o u r c e may be regarded a s t h e envelope of t h e fronts of
s p h e r i c a l waves originating f r o m different points of the initial surface.
The interpretation of the front a s the envelope of a given family of s u r f a c e s
l e a d s t o Monge's representation of the front s u r f a c e a s consisting of the
contact lines of the enveloping surface with the enveloped s u r f a c e s . The
equation of the enveloping surface, which moves in s p a c e and is the front
of s o m e propagation process, can be written a s
c , p ( a T / a t ) = kAT + q ( a M / a t ) ,
aM/at = / ( M ) . I
The first equation of this s y s t e m is the equation of heat conduction f r o m
a s o u r c e , where q designates the specific heat liberated o r absorbed
during the transformation of 1 g solid material; the second equation
controls the reaction r a t e of the m a t e r i a l , M representing the quantity
of t r a n s f o r m e d phase of the m a t e r i a l p e r unit volume. A combined p r o c e s s
184
is regarded a s one inwhich the basic quantities determining it a r e interrelated.
In o u r c a s e these a r e the temperature T and m a s s M . Therefore, in
combined processes they must be connected by the functional relationship
d M l d T , which is always assumed t o exist in these processes. Assumptions
relative t o higher o r d e r derivatives w i l l be made l a t e r . The first of
equations ( 3 ) now takes the f o r m
8T / a t = ( a l a ) AT, (4 1
P-(-/-...)
daT d T . (71
dt2 dS
d a M / d T 2= 0 . (10)
185
Hence, equation (6) becomes
g; + a H,, + --
a AV
--gc a dM = 0,
df
The minus sign before the radical has been discarded, s i n c e it is physically
meaningless under certain conditions. In fact, the phenomenon can exist
even if the derivative dfldM vanishes. When the sign before the radical is
negative the equality dfldM=O l e a d s t o the vanishing of the velocity gE a s
well, which is inadmissible.
We now a s s u m e that
G(g)2-&<0.25,
i. e., the second t e r m under the r a d i c a l can be neglected.
This inequality is satisfied if any one of the cofactors of (15) is s m a l l .
The value of each can be determined for the p a r t i c u l a r phenomenon
considered. Thus, if (15) is valid, then
a AV
4=----'
a H,
186
I
... ... - .. .
F o r t h e complete removal of the transformed m a t e r i a l by t h e g a s -
dynamic s t r e a m , no considerable e r r o r is introduced if we take
T -To
q-.
Too - To
2
H,=1; A V = -
Allowing f o r (22), equation (16) can be written as
gr = - ( Z a / a r ) . (23)
187
can however be derived theoretically. Consider a spherical body, f o r which
equation (14)(allowing f o r (22)) becomes
-
d t = --
dr ar (If v-;%).
Integrating and performing a number of elementary transformations, we
obtain
where
and
p = L L
a dM'
5'= So - G t , (351
188
where
4n = _1 =
_ i (37 1
So r2 (C/4n)-4(a/a) .
This equation implies that i f the phenomenon under observation develops
according t o Sreznevskii's law, t h e reaction r a t e of the m a t e r i a l
v a r i e s such that the Gaussian curvature of the s u r f a c e r e m a i n s constant.
This can be realized only i f t h e r e exists a plane dividing the s p h e r e into
two p a r t s , in which t h e m a t e r i a l i s transformed at different r a t e s . This is
observed in t h e presence of n a t u r a l convection, or when the m a t e r i a l
react+-ontakes place in a moving medium. The integration of ( 2 7 ) i s
possible only when n i s independent of the radius r of the spherical particle.
The Sreznevskii l a w is therefore m o r e general than t h e conservation law
of the Gaussian curvature of the reaction surface, and can be satisfied
in the f o r m of (26), while the Gaussian curvature of the s u r f a c e can vary.
T h e r e is one corollary f r o m t h e l a w s established, concerning the
geometry of t h e moving reaction s u r f a c e of the material. Gauss' famous
t h e o r e m a s s e r t s that the total (Gaussian) curvature of the s u r f a c e
depends only on the l i n e a r elements (and t h e i r f i r s t - and second-order
d e r i v a t i v e s ) on the surface, and t h e r e f o r e does not change with t h e defor
mation of the s u r f a c e , It follows that e v e r y s u r f a c e of constant Gaussian
curvature is applicable t o (can be deformed into) any other s u r f a c e of the
s a m e curvature. We r e c a l l that two s u r f a c e s S and S * a r e isometric, i. e.,
the s u r f a c e S i s applicable to the s u r f a c e S* (or S* r e p r e s e n t s a deformation
of S ) if t h e i r line elements a r e equal. It follows f r o m Gauss' t h e o r e m that
developable s u r f a c e s (surfaces enveloping the oo-family of tangent planes
t o a s p a c e c u i v e ) , cones, and cylinders a r e applicable t o a plane. Thus,
s i n c e the material transformation takes place under specific conditions,
fulfilling the Sreznevskii l a w and conserving the Gaussian curvature of the
surface, the motion of the reaction s u r f a c e can be represented a s
an uniform and continuous deformation of the whole s u r f a c e o r individual
p a r t s of it. If the solid p o s s e s s e s s p h e r i c a l symmetry, the surface
is transformed with t i m e in such a way that s y m m e t r y w i l l be p r e s e r v e d
for the front and back of the surface, subject t o different exchange conditions
with t h e surrounding medium. For bodies of conical, cylindrical, o r any
other developable s u r f a c e the dynamic transformation of t h e m a t e r i a l will
proceed such that even when Sreznevskii's law is obeyed and the Gaussian
c u r v a t u r e is conserved f o r t h e whole s u r f a c e o r its individual p a r t s , t h e
body shape can be altered substantially by deformation (e. g., transformed
into a plane).
189
The r e s u l t s listed r e p r e s e n t only s o m e of the advantages of the proposed
model, whose simplicity considerably facilitates t h e study of many pheno
mena of various combined p r o c e s s e s .
The theoretical model proposed does not contradict published data o r
our own experimental r e s u l t s , which a r e considered below.
One of the first papers emphasizing the importance of Sreznevskii's law
was 151, where the law was applied t o the evaporation kinetics of liquid
drops and volatile solid bodies. Zubkov /6/ observed the evaporation
kinetics of naphtalene and camphor s p h e r e s , over an air-velocity range of
0.5-lOm/sec, and air t e m p e r a t u r e of 15-120". These a i r t e m p e r a t u r e s
enabled both pure sublimation and evaporation through a liquid phase t o be
studied. The investigations showed that t h e r e is no difference between the
evaporation kinetics of naphtalene in liquid and solid phases, i. e., that
although the d e c r e a s e in m a t e r i a l p e r unit t i m e is different, the p r o c e s s is
the s a m e qualitatively. This is understandable in the f r a m e of the model
discussed. Figure 1 shows c m s e c u t i v e stills of an evaporating naphthalene
s p h e r e . The evaporation proceeds at different r a t e s on different portions of
the surface. A zone of minimum evaporation r a t e is c l e a r l y visible, and
it divides the s u r f a c e of the s p h e r e into front and back p a r t s . The zone of
minimum evaporation is situated where the impinging s t r e a m divides. By
observing t h e shape of the front and back parts, the naphtalene was s e e n
t o evaporate at different but comparable r a t e s . According t o Zubkov's data,
in a first approximation, both the front and back p a r t s p r e s e r v e a s p h e r i c a l
shape during evaporation. In addition, Zubkov's calculations of t h e r a d i i
of the front and back s u r f a c e s showeu that they vary linearly, i. e., following
Sreznevskii's law. Similar r e s u l t s w e r e obtained by Aksel'rud /7/ for
compressed s a l t p e t e r s p h e r e s and sodium chloride dissolved in a liquid
s t r e a m . Figure 2 r e p r e s e n t s successive stages in the dissolution of the
s a l t p e t e r s p h e r e s . A comparison of Figures 1 and 2 indicates that both
p r o c e s s e s obey the s a m e law.
FIGURE 1. Successive stills of the evaporation of naphthalene spheres for low subsonic speeds of the air stream.
190
0 1 2 3 4 5 6
- 2, sec
3-30? 4-40',
191
Sreznevskii's law, and conserves the Gaussian curvature of the melting
surface.
We conducted experiments on the melting of cylindrical ice models with
hemispherical heads, placed lengthwise in a water s t r e a m . Sreznevskii's
lzw w a s obeyed, and the Gaussian curvature of the hemispherical nose
surface and the cylindrical part of the model s u r f a c e was conserved.
The numerous experiments conducted at the Power Institute im.
Krzhizhansvskii on the melting of various metals and alloys in gas-dynamic
s t r e a m s corroborate the front c h a r a c t e r of the displacement of the
reaction surface. The shortening of the conic s a m p l e along its axis
obeyed a l i n e a r function of time, undisturbed by the variation in cone
angle. If the cone length versus time is plotted, a straight line is obtained
whose slope d e c r e a s e s with increasing cone angle (Figure 4).
Experimental characteristics of the dynamic melting and sublimation of
axisymmetrical bodies, data on the contraction of the models, and values
of the radius of the models' noses a r e given in / l o / . The varying geometry
of the reaction s u r f a c e of the material draws attention to the complex
interconnected p r o c e s s e s taking place. The approximate theoretical
calculation of the reaction surface does not, in our opinion, elucidate the
t r u e development mechanism of the investigated phenomenon. Stetson's
r e s u l t s s e e m t o indicate that the phenomenon of dynamic melting and
sublimation developed according t o the model of front dynamic processes.
a-sphere @I" in an oxygen stream: M = 1.3; T*= 820C. T i m e interval 0.5sec between the stills;
b-sphere @I" in an air stream; M = 3; T*= 870C. Time interval 5 min E t w e e n the stills.
192
We conducted an experimental study of the combustion of s p h e r i c a l and
cylindrical coal particles in a supersonic a i r and oxygen s t r e a m f o r M =
= 1 . 3 - 3 . 5 and s t r e a m stagnation t e m p e r a t u r e s between 780 and 950". The
r e s u l t s (Figure 5 ) show that t h e model of front p r o c e s s e s is a l s o applicable
t o combustion processes. The a r e a of the reaction s u r f a c e of the s p h e r i c a l
coal particle varied linearly with t i m e : it was a l s o established that the
Gaussian curvature is conserved on the front and back p a r t s . These
observations a r e illustrated in F i g u r e s 6 and 7. The r e s u l t s obtained a r e
qualitatively s i m i l a r t o those of Blinov and Khaikin on the combustion of
carbon in low-speed gas s t r e a m s . T h e i r r e s u l t s a r e discussed in detail by
Predvoditelev 1 3 1 f r o m the viewpoint of front p r o c e s s e s .
M 1 4; T * = 8 2 0 T .
S, c m 2 s.105 cm2
160
I20
0 2 4 6 T,sec D 20 40
a b
FIGURE I. Variation in the surface area of spherical models 6 I mm from coal particles
during combustion in a stream of air and oxygen:
a-oxygen; M = 2; T*= 82OOC; b-air: M = 2.5; 1- T * = 82OoC; 2- T*= 860; 3- T*= 950.
193
Of considerable i n t e r e s t a r e t h e experimental r e s u l t s on the interaction
between graphite s p h e r i c a l bodies and a high-enthalpy supersonic s t r e a m
(M = 2 . 4 ) of a i r , argon, and nitrogen / l l / . The paper presents data on the
stable s u r f a c e shape acquired by the cylindrical graphite bodies in a plasma
j e t of a i r and two other gases at stagnation t e m p e r a t u r e s above 8000" and
s t r e a m speeds about 3000 m / s e c . The authors established that graphite
rods rapidly acquired a stable conic s u r f a c e geometry with a s h a r p nose,
which remained invariable throughout the experiment. These r e s u l t s a r e
illustrated in Figure 8. They a l s o presented data on the contraction of the
s a m p l e and the removal of transformed m a s s a s a function of time, and
showed that both phenomeca a r e described by l i n e a r relationships. The
experiments w e r e therefore found t o obey both laws of motion of the
reaction s u r f a c e of the m a t e r i a l in combined p r o c e s s e s developed
by a front.
Bibliography
194
6. Z u b k o v , V. I. Ob isparenii sharikov tverdykh t e l v potoke gaza
(Evaporation of Solid Spheres in a Gas Stream).-Doklady AN
SSSR, Vol. 123, No. 5. 1958.
7. A k s e 1 r u d, G. A. K voprosu o kinetike rastvoreniya tverdoi
chastitsy, dvizhushcheisya v potoke zhidkosti (Dissolution
Kinetics of a Solid Particle, Moving in a Gas Stream).-Zhurnal
Fizicheskoi Khimii, Vol. 28, No. 10.
8. E r o s h e n k o , V.M. a n d Y u . N . T e r e n t ' e v . Analiz sublimatsii
tverdykh t e l v sverkhzvukovom Potoke (The Sublimation of Solids
i n a Supersonic Stream).- In this collection, p.196.
9. D o n a l d , J. The Melting of Ice in a Hot Humid S t r e a m of Air.-J.
Fluid Mech., Vol. 8, p a r t I. 1960.
10. S t e t s o n , N, Issledovanie ablyatsii konicheskikh modelei (Investiga
tion of the Ablation of Conic Models).- VRT, No. 2. 1966.
11. C h r i s t e n s e n , D. and R . D . B u h l e r . O n t h e Stable Shape of an
Ablating Graphite Body.-J. Aero. Sci., Vol. 26, No. 1. 1959.
195
V . M . Eroshenko, Yu.N . Terent' ev
THE SUBLIMATION OF SOLIDS IN A
SUPERSONIC STREAM
196
accuracy the coordinates of t h e model contour at different s t a g e s of its
evaporation. This is illustrated in F i g u r e 4. Besides t h e experimental
points, the figure shows c i r c l e s whose radius equals that of the initial
sphere, but relative t o a point situated at a distance f r o m t h e initial center
equal t o the displacement of the leading stagnation point of the surface. It is
s e e n f r o m Figure 4 that the experimental points characterizing t h e front
s u r f a c e of the model at different sublimation t i m e s are located n e a r t h e
a r c s of the c i r c l e of initial radius, plotted by t h e method described above.
13.5
._
FIGURE 1. Naphthalene samples FIGURE 2. Scheme of the mounting and introduction
of models in the stream :
1-model; %lever for introducing the model in the
stream; 3-refrigerator; 4-De Lava1 nozzle; 5-low
pressure chamber.
Stream directi
197
Thus, when a s p h e r e evaporates in a supersonic s t r e a m , t h e front
s u r f a c e of the model conserves a s p h e r i c a l shape with its initial radius.
In other words, t h e front s u r f a c e evaporates as if the c e n t e r of a s p h e r e
with constant radius (equal t o its initial value) moves downstream. A
s i m i l a r pattern was observed f o r models of 7 and 10 d i a m e t e r and other
stagnation t e m p e r a t u r e s . The radius of the front p a r t of the s p h e r e
remained unchanged up t o v e r y extensive evaporation. In many experiments
the displacement of t h e leading stagnation point reached 70- 85 % of the value
of the initial radius, without any change in the total curvature of the s p h e r i c a l
model surface. Experimental plots of Ah = f(z) a r e given in F i g u r e 5, with
Ah representing the distance between the leading stagnation point of the model
at the initial moment of the p r o c e s s and the stagnation point corresponding to
s o m e instant z during evaporation. The relationship At = f ( z ) is linear, i. e.,
for given shape and dimensions of the model the displacement velocity of
the leading stagnation point r e m a i n s almost constant for invariable s t r e a m
parameters.
Ah, mm
20
2.0
198
--
Ah. nun
I-
J
min
FIGURE 6. Displacement of the leading FIGURE 7. Schlieren photographs
stagnation point of naphthalene spheres of of the successive stages in the
different initial diameters under identical evaporation of conical naphthalene
flows conditions: models.
1-sphere $7"; 2-$ 10".
199
s p h e r i c a l front and back regions varied with t i m e . In addition, both the
total s u r f a c e of the s p h e r i c a l particle studied, and its front and back
s u r f a c e s , v a r i e d linearly with time.
Thus, the r e s u l t s of / 2 / a r e qualitatively s i m i l a r t o our r e s u l t s &d can
be explained by the front model of the motion of the evaporation s u r f a c e .
Ah. mm
process 1 'a
FIGURE 9. Contours of the successive sublimation stages
of naphthalene cones for stream stagnation temperatures
of 30" and 75C.
200
- -- '- -- --- -- - -
where gc is the n o r m a l displacement velocity of the reaction surface
of the material; a is the thermal diffusivity; H,, and AV a r e the first and
second s c a l e factors of the reaction surface V ; dfldM is the derivative
of the function controlling the reaction kinetics with respect to the
quantity of transformed material. The quantity a is defined a s
P9 = C"P ( T , - TW). (3 1
V = 2 + - br2
R* '
where r = v v ; R is the radius a t the cone base; b is the distance
f r o m the cone apex t o the base.
This means that the evaporation of the cone is determined by the velocity
of the reaction surface of the m a t e r i a l in the direction of the solid
residue. The first and second s c a l e factors of the surface V have the f o r m
dz
H" = 1; 4 = Zi
If the m a s s variation p e r unit volume obeys the l i n e a r l a w
aM
- = (7
at kM,
where k is a proportionality constant, then
20 1
For a cone of unit volume we have
_ddfM_ - - 3
36
&?EH,* (11
Allowing for (5), (6), and ( l l ) ,we can write equation (1)in the f o r m
ba
gefDgE=O. (12
d r = - _ ba
_
dt aR2 '
z = -
-t
ba
zRZ
+ const.
It is easily shown that the constant has the value b, if the coordinate
f r a m e is selected a s follows:
ctg 0
Ah = --ta
aR 2
or
202
. . .
This is the final f o r m of the equation which must be obeyed by t h e l a w
of variation of the cone length. To u s e it, w e must determine T , , and T , .
Suppose a c r i t i c a l t e m p e r a t u r e T,, exists below which no sublimation of
the m a t e r i a l is observed. This t e m p e r a t u r e corresponds t o t h e melting
t e m p e r a t u r e in metals. F o r sublimating bodies it a l s o depends on the
manufacture of the material. The m o r e friable the material, the lower t h e
sublimation temperature. Its value w i l l be highest f o r a compressed
sublimating m a t e r i a l with a polished s u r f a c e .
In o u r experiments t h e samples w e r e made of compressed naphthalene,
which hardly sublimated during many hours in a gas-dynamic s t r e a m at
t e m p e r a t u r e s around 25". This t e m p e r a t u r e was adopted as t h e c r i t i c a l
sublimation t e m p e r a t u r e . To calculate T, by formula (3), one must
determine the heat of m a t e r i a l transformation (i. e., the heat of sublimating
1 g m a t e r i a l ) a s follows. Before a molecule leaves t h e solid material, it
must do work in o r d e r t o approach the surface. This work is equivalent
t o the escape of a molecule f r o m the s p h e r e of influence of the other
molecules (namely, a distance f r o m the s u r f a c e l a r g e r than t h e radius of
influence of the cohesion forces), and i s equal t o t h e heat of fusion. After
the molecule has reached the surface, it must do work against the s u r f a c e
forces. Only a f t e r t h e s e forces have been overcome can it evaporate. This
work i s identified with the heat of evaporation o r sublimation.
Designate by Q the total latent heat of evaporation, by q r the sublimation
heat, and by q , the heat of fusion. The following equality holds:
TOO. E
30 75 0.00450 0.00475
50 85 0.00570 0.00572
203
Satisfactory agreement is observed between the theoretical and experi
mental r e s u l t s . Thus, in addition t o examining the qualitative laws of
sublimation in a gas-dynamic s t r e a m , we have established the suitability
of the model of front p r o c e s s e s f o r analyzing the complex phenomenon of
heat and m a s s exchange.
B i b l i og r a p hy
1. E r o s h e n k o , V . M . , A . V . M o r a r , a n d Y u . N . T e r e n t ' e v . Frontovoi
kharakter poverkhnosti prevrashcheniya veshchestva p r i vzai
modeistvii s vneshnei s r e d o i (Front Nature of the Reaction
Surface with the Medium).--In this collection, p. 183.
2. Z h u b k o v , V. I. Ob isparenii sharikov tverdykh tel v potoke gaza
(Evaporation of Solid Spheres in a Gas Stream). - DAN SSSR,
Vol. 123, No. 5. 1958.
204
Yu.N . Vorontsov, V .M . Eroshenko, A . V .Moray,
Yu.N . Terent ' e v
205
I
a b c
FIGURE 1. Successive stages in the motion of the isothermal
translucent surface of polytetrafluoroethylene in spherical
models of @ = 7 mm for M= 2.5 and T k = 770K:
a-brass holder; b, c-porcelain holder. The time interval
between the frames i s 10 sec.
206
It is seen f r o m F i g u r e 3 that for s t r e a m stagnation t e m p e r a t u r e s over
the range 600- 900"K, t h e t i m e of complete recrystallization of a spherical
polytetrafluoroethylene model v a r i e s according t o a near-exponential l a w .
At t e m p e r a t u r e s above 900K t h e recrystallization is independent of the
s t r e a m stagnation t e m p e r a t u r e , and is determined only by internal
p r o c e s s e s , i. e., the maximum recrystallization r a t e .
t. m m
t, sec
FIGURE 2. Velocity of motion of the isothermal FIGURE 3 . Total time taken by a polytetrafluoro
translucent surface i n a spherical model for Y=2.5: ethylene sphere to become translucent as a function
of the stagnation temperature of the stream.
1- Ti0-7630K; 2-793"; 3-813'; 4-833".
Bibliography
207
2. S r e z n e v s k i i , V. Ob i s p a r e n i i zhidkostei (Evaporation of Liquids).-
ZhRFKhO, vol. 14. 1882.
3. P r e d v o d i t e l e v , A.S. 0 frontovom p r o t s e s s e ispareniya i subli
m a t s i i sfericheskikh chastits (Front Process of Evaporation and
the Sublimation of Spherical P a r t i c l e s ).-Trudy Odesskogo
Gosudarstvennogo Universiteta, Seriya Fizicheskikh Nauk, No. 7.
1960.
208
\
Natioaal AWOIUL~~~CS
andSpace Administration POSTAGE AND FEES PAID
NATIONAL AERONAUTICS AND
-
WMHlNGTONo 0. C.
OFFICIAL BUSINESS
20516 SPACE ADMINISTRATION