In relating the pressures, temperatures and densities immediately across
the shock, it is conventional to consider the gas motion in relation to the
shock front. In such a frame of reference, the gas
enters the shock at a relative velocity , and leaves with a relative velocity
. In shock tube jargon, the shock is
then considered to be at rest;
is the gas velocity measured in shock-fixed coordinates and
is that measured in laboratory-fixed coordinates. These two frames of
reference are related by:
(7–5) |
(7–6) |
where is the shock velocity. Gas conditions associated with the incident shock
in the two coordinate systems are shown in Figure 7.3: Laboratory-fixed and Incident-shock-fixed coordinate systems
. The
Rankine-Hugoniot relations for properties across the
incident shock front are
(7–7) |
(7–8) |
(7–9) |
Utilizing the equation of state (Ideal Gas Law) and Equation 7–7
to eliminate the velocity
and
from Equation 7–8
and Equation 7–9
results in the following
expressions for the pressure and temperature ratios across the incident shock:
(7–10) |
(7–11) |
Since we assume no change in gas composition across the shock,
is a function of temperature alone and, hence, Equation 7–10
and Equation 7–11
represent a system of two
equations in two unknowns. The solution gives
and
when conditions before the incident shock are specified. Knowing these,
is determined from the equation of state and
from Equation 7–7
.
An iterative procedure is employed to solve Equation 7–10
and Equation 7–11
for and
. Letting
and
be the temperature and pressure ratios, respectively, across the shock,
Equation 7–10
can be solved for
in terms of
(see Equation 7–13
) to yield
(7–12) |
This expression is then substituted into Equation 7–11
to yield one equation with one
unknown, . Within Ansys Chemkin, a routine called
ZEROIN [52], which finds the zeros of functions, is employed to determine the value of
that satisfies this equation. An initial guess for
is provided by assuming that the test
gas is ideal (
and
are constant and independent of temperature). For ideal gases
(7–13) |
where is the specific heat ratio and
is the Mach
number of the incident shock.
(7–14) |
Many times the experimentalist reports the incident shock speed,
, and temperature and pressure behind the shock,
and
, respectively. Before the experiment can be modeled, however, the gas
velocity behind the shock must be determined. Employing the equation of state in Equation 7–10
to eliminate
results in
(7–15) |
This equation and Equation 7–11
again represent two equations
in two unknowns. The solution gives and
and from these the density in region 1 is determined from the equation of
state. The velocity behind the shock,
, is determined from Equation 7–7
. The solution to Equation 7–15
and Equation 7–11
is analogous to that already
described for Equation 7–10
and Equation 7–11
.