21.11.2014 Views

Untitled - Aerobib - Universidad Politécnica de Madrid

Untitled - Aerobib - Universidad Politécnica de Madrid

Untitled - Aerobib - Universidad Politécnica de Madrid

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

4.5. APPLICATIONS 101<br />

By substituting these values and that of the temperature given by Eq. (4.37)<br />

into Eq. (4.36), the following expression, for the Hugoniot curve, is obtained<br />

xy = q + γ 2 − 1<br />

(y − x) . (4.40)<br />

γ 2 + 1<br />

Therefore, the Hugoniot curve is an equilateral hyperbola, with an asymptote parallel<br />

to the axis y for x 0 = (γ 2 − 1) / (γ 2 + 1). This curve cuts axis x at the point x 1 =<br />

q (γ 2 + 1) / (γ 2 − 1).<br />

Let us see how the Chapman-Jouguet points are <strong>de</strong>termined. The equation of<br />

an isentropic for the burnt gases is<br />

yx γ2 = const. (4.41)<br />

Therefore, the Chapman-Jouguet points are obtained by solving the system formed<br />

by Eq. (4.40) and by the equations that express the tangency condition of the curves<br />

Eq. (4.40) and Eq. (4.41), namely,<br />

(γ 2 + 1) y + (γ 2 − 1)<br />

(γ 2 + 1) x − (γ 2 − 1) = γ y<br />

2<br />

x . (4.42)<br />

The solution of Eq. (4.40) and Eq. (4.42) gives for x and y the following equations<br />

y 2 − ( (γ 2 + 1)q − (γ 2 − 1) ) y + q = 0, (4.43)<br />

x 2 − 1 γ 2<br />

(<br />

(γ2 + 1)q + (γ 2 − 1) ) x + q = 0. (4.44)<br />

When solving these equations the root x 1 < 1 must be assign to the root y 1 > 1 and<br />

root x 2 > 1 to the root y 2 < 1.<br />

In virtue of Eq. (4.8), the propagation velocity is<br />

v 1<br />

= 1<br />

√<br />

y − 1<br />

√<br />

a 1 γ1 1 − x , (4.45)<br />

where a 1 is the sound velocity of the unburnt gases, and γ 1 the ratio of their heat<br />

capacities.<br />

Similarly, the velocity of the burnt gases with respect to the wave front is<br />

v 2<br />

=<br />

x<br />

√<br />

y − 1<br />

√<br />

a 1 γ1 1 − x . (4.46)<br />

The corresponding Mach number M 2 is<br />

√<br />

M 2 = √ 1 ( )<br />

x y − 1<br />

. (4.47)<br />

γ2 y 1 − x

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!