next up previous index
Next: Solution of the Governing Up: gasDynamics Previous: The effects of Real   Index


Normal Shock

In this chapter the relationships between the two sides of normal shock are presented. In this discussion, the flow is assumed to be in a steady state, and the thickness of the shock is assumed to be very small. A discussion on the shock thickness will be presented in a forthcoming section5.1.

Figure: A shock wave inside a tube, but it can also be viewed as a one-dimensional shock wave.
\begin{figure}\centerline{\includegraphics {cont/shock/cv} }
\end{figure}
A shock can occur in at least two different mechanisms. The first is when a large difference (above a small minimum value) between the two sides of a membrane, and when the membrane bursts (see the discussion about the shock tube). Of course, the shock travels from the high pressure to the low pressure side. The second is when many sound waves ``run into'' each other and accumulate (some refer to it as ``coalescing'') into a large difference, which is the shock wave. In fact, the sound wave can be viewed as an extremely weak shock. In the speed of sound analysis, it was assumed the medium is continuous, without any abrupt changes. This assumption is no longer valid in the case of a shock. Here, the relationship for a perfect gas is constructed.

In Figure (5.1) a control volume for this analysis is shown, and the gas flows from left to right. The conditions, to the left and to the right of the shock, are assumed to be uniform5.2. The conditions to the right of the shock wave are uniform, but different from the left side. The transition in the shock is abrupt and in a very narrow width.

The chemical reactions (even condensation) are neglected, and the shock occurs at a very narrow section. Clearly, the isentropic transition assumption is not appropriate in this case because the shock wave is a discontinued area. Therefore, the increase of the entropy is fundamental to the phenomenon and the understanding of it.

It is further assumed that there is no friction or heat loss at the shock (because the heat transfer is negligible due to the fact that it occurs on a relatively small surface). It is customary in this field to denote x as the upstream condition and y as the downstream condition.

The mass flow rate is constant from the two sides of the shock and therefore the mass balance is reduced to

$\displaystyle \rho_{x} U_{x} = \rho _{y} U_{y}$ (5.1)

In a shock wave, the momentum is the quantity that remains constant because there are no external forces. Thus, it can be written that

$\displaystyle P_{x} - P_{y} = \left(\rho_{x} {U_{y}}^2- \rho_{y} {U_{x}}^2 \right)$ (5.2)

The process is adiabatic, or nearly adiabatic, and therefore the energy equation can be written as

$\displaystyle C_{p} T_{x} + {{U_{x}}^{2} \over 2} = C_{p} T_{y} + {{U_{y}}^{2} \over 2}$ (5.3)

The equation of state for perfect gas reads

$\displaystyle P = \rho R T$ (5.4)

If the conditions upstream are known, then there are four unknown conditions downstream. A system of four unknowns and four equations is solvable. Nevertheless, one can note that there are two solutions because of the quadratic of equation (5.3). These two possible solutions refer to the direction of the flow. Physics dictates that there is only one possible solution. One cannot deduce the direction of the flow from the pressure on both sides of the shock wave. The only tool that brings us to the direction of the flow is the second law of thermodynamics. This law dictates the direction of the flow, and as it will be shown, the gas flows from a supersonic flow to a subsonic flow. Mathematically, the second law is expressed by the entropy. For the adiabatic process, the entropy must increase. In mathematical terms, it can be written as follows:

$\displaystyle s_{y} - s_{x} > 0$ (5.5)

Note that the greater-equal signs were not used. The reason is that the process is irreversible, and therefore no equality can exist. Mathematically, the parameters are $ P, T, U,$ and ρ, which are needed to be solved. For ideal gas, equation (5.5) is

$\displaystyle \ln {T_y \over T_x} - \left(k - 1\right) {P_y \over P_x} > 0$ (5.6)

It can also be noticed that entropy, s, can be expressed as a function of the other parameters. Now one can view these equations as two different subsets of equations. The first set is the energy, continuity, and state equations, and the second set is the momentum, continuity, and state equations. The solution of every set of these equations produces one additional degree of freedom, which will produce a range of possible solutions. Thus, one can have a whole range of solutions. In the first case, the energy equation is used, producing various resistance to the flow. This case is called Fanno flow, and Chapter (9) deals extensively with this topic. The mathematical explanation is given Chapter (9) in greater detail. Instead of solving all the equations that were presented, one can solve only four (4) equations (including the second law), which will require additional parameters. If the energy, continuity, and state equations are solved for the arbitrary value of the Ty , a parabola in the T-s diagram will be obtained. On the other hand, when the momentum equation is solved instead of the energy equation, the degree of freedom is now energy, i.e., the energy amount ``added'' to the shock. This situation is similar to a frictionless flow with the addition of heat, and this flow is known as Rayleigh flow. This flow is dealt with in greater detail in Chapter (10).

Figure: The intersection of Fanno flow and Rayleigh flow produces two solutions for the shock wave.
Image intersection
Since the shock has no heat transfer (a special case of Rayleigh flow) and there isn't essentially any momentum transfer (a special case of Fanno flow), the intersection of these two curves is what really happened in the shock. In Figure (5.2), the intersection is shown and two solutions are obtained. Clearly, the increase of the entropy determines the direction of the flow. The entropy increases from point x to point y. It is also worth noting that the temperature at M=1 on Rayleigh flow is larger than that on the Fanno line.



Subsections
next up previous index
Next: Solution of the Governing Up: gasDynamics Previous: The effects of Real   Index
genick 2007-11-21