Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

4.6 The Hugoniot Relation

Chalmers University of Technology
Department of Mechanical Engineering
Division of Fluid Dynamics

The Hugoniot equation is an alternative normal shock relation based on thermodynamic quantities only. It is derived from the governing equations and relates the change in energy to the change in pressure and specific volume. The starting point of the derivation of the Hugoniot equation is the governing equations (Eqns (4.1) - (4.4)).

The continuity equation is rewritten and inserted into the momentum equation

u1=(ρ2ρ1)u2u_1=\left(\frac{\rho_2}{\rho_1}\right) u_2

Replace u1u_1 in Eqn. (4.2) using Eqn. (4.87)

ρ1(ρ2ρ1)2u22+p1=ρ2u22+p2\rho_1 \left(\frac{\rho_2}{\rho_1}\right)^2 u^2_2 +p_1=\rho_2 u^2_2 + p_2
u22(ρ1(ρ2ρ1)2ρ2)=(p2p1)u^2_2\left(\rho_1\left(\frac{\rho_2}{\rho_1}\right)^2-\rho_2\right)=\left(p_2-p_1\right)
u22((ρ2ρ1)(ρ2ρ1))=(p2p1)u^2_2\left(\left(\frac{\rho_2}{\rho_1}\right)\left(\rho_2-\rho_1\right)\right)=\left(p_2-p_1\right)
u22=(ρ1ρ2)p2p1ρ2ρ1u^2_2=\left(\frac{\rho_1}{\rho_2}\right)\frac{p_2-p_1}{\rho_2-\rho_1}

Eqn. (4.87) and (4.91) gives

u12=(ρ2ρ1)p2p1ρ2ρ1u^2_1=\left(\frac{\rho_2}{\rho_1}\right)\frac{p_2-p_1}{\rho_2-\rho_1}

Eqn. (4.91) and Eqn. (4.92) inserted in the energy equation (Eqn. (4.4)) gives

h1+12(ρ2ρ1)(p2p1ρ2ρ1)=h2+12(ρ1ρ2)(p2p1ρ2ρ1)h_1 + \frac{1}{2}\left(\frac{\rho_2}{\rho_1}\right)\left(\frac{p_2-p_1}{\rho_2-\rho_1}\right)=h_2 + \frac{1}{2}\left(\frac{\rho_1}{\rho_2}\right)\left(\frac{p_2-p_1}{\rho_2-\rho_1}\right)
h2h1=p2p12[(ρ2ρ1)(1ρ2ρ1)(ρ1ρ2)(1ρ2ρ1)]h_2-h_1=\frac{p_2-p_1}{2}\left[\left(\frac{\rho_2}{\rho_1}\right)\left(\frac{1}{\rho_2-\rho_1}\right)-\left(\frac{\rho_1}{\rho_2}\right)\left(\frac{1}{\rho_2-\rho_1}\right)\right]
h2h1=p2p12[ρ22ρ12ρ1ρ2(ρ2ρ1)]=p2p12[ρ2+ρ1ρ1ρ2]h_2-h_1=\frac{p_2-p_1}{2}\left[\frac{\rho^2_2-\rho^2_1}{\rho_1\rho_2(\rho_2-\rho_1)}\right]=\frac{p_2-p_1}{2}\left[\frac{\rho_2+\rho_1}{\rho_1\rho_2}\right]
h2h1=p2p12(1ρ1+1ρ2)h_2-h_1=\frac{p_2-p_1}{2}\left(\frac{1}{\rho_1}+\frac{1}{\rho_2}\right)

Now, replacing the enthalpies with internal energies using h=e+p/ρh=e+p/\rho gives

e2e1=p1ρ1p2ρ2+p2p12(1ρ1+1ρ2)e_2-e_1=\frac{p_1}{\rho_1}-\frac{p_2}{\rho_2}+\frac{p_2-p_1}{2}\left(\frac{1}{\rho_1}+\frac{1}{\rho_2}\right)

which after some rewriting becomes the Hugoniot equation

e2e1=p2+p12(1ρ11ρ2)=p2+p12(ν1ν2)e_2-e_1=\frac{p_2+p_1}{2}\left(\frac{1}{\rho_1}-\frac{1}{\rho_2}\right)=\dfrac{p_2+p_1}{2}(\nu_1-\nu_2)

To give an idea about how the normal shock relates to an isentropic compression (a flow compression process without losses) the change in flow density as a function of pressure ratio is compared in Figure Figure 4.10. One can see that the normal-shock compression is more effective but less efficient than the corresponding isentropic process.

comparison of a normal shock compression process and an isentropic compression

Figure 4.10:Comparison of a normal shock compression process and an isentropic compression

Introducing C as the massflow per unit area (which is a constant)

ρ1u1=ρ2u2=C\rho_1 u_1 = \rho_2 u_2 = C

Inserted into the momentum equation this gives

p1+C2ρ1=p2+C2ρ2p_1+\dfrac{C^2}{\rho_1}=p_2+\dfrac{C^2}{\rho_2}

or

p2p1ν2ν1=C2\dfrac{p_2-p_1}{\nu_2-\nu_1}=-C^2

which implies that all possible solutions to the governing equations must be located on a line in pνp\nu-space (the so-called Rayleigh line). If we add the Hugoniot relation to this we will find that there are two possible solutions, the upstream condition and the condition downstream of the normal shock and the flow cannot be in any of the intermediate stages. The normal-process is a so-called wave solution to the governing equations where the flow state must jump directly from one flow state to another without passing the intermediate conditions. If we add heat or friction to the problem we will instead get continuous solutions as we will see in the following sections. Figures Figure 4.11 and Figure 4.12 shows a normal shock process in a pνp\nu- and TsTs-diagram, respectively. Note that the flow passes the characteristic conditions as it is going through the shock, which means that the flow goes from supersonic to subsonic.

the normal-shock process illustrated in a Pv-diagram

Figure 4.11:The normal-shock process illustrated in a pνp\nu-diagram

The normal-shock process illustrated in a Ts-diagram

Figure 4.12:The normal-shock process illustrated in a TsTs-diagram