15. Compressible Fluid Mechanics#

15.1. Compressible Inviscid Fluid Mechanics#

15.1.1. Shocks#

15.1.2. Quasi-1d flows#

If no shock occurs in the flow, Euler equations in differential form governs the dynamics of the flow

{tρ+(ρu)=0t(ρu)+(ρuu+PI)=0t(ρet)+(ρuht)=0

Quasi-1d model for steady flows is a simple model that provides good-enough results for flows delimited by streamlines that varies gently in the streamwise direction (or by solid walls, in the limit of the high-Reynolds flow without separations, where viscous effects are confined to a thin region - boundary layer - close to the walls).

This model is derived integrating over the sections of the stream tube, so that the physical quantities are functions of the streamwise coordinate x only. Integration over an elementary volume between sections at x and x+dx of the mass and momentum equations, in steady conditions, gives

{d(ρuA)=0d(ρu2A)+d(PA)+Slat(x,dx)Pnx=0 ,

where the last equation comes from the contribution of the lateral surfaces, that has non zero contribution in the streamwise component of momentum equation if sections is not constant, and thus the unit normal vector n^ on the lateral surface has non-zero x-component nx. This contribution on the elementary lateral surface (where pressure is assumed to be uniform), can be evaluated summing and subracting the contributions on the A(x) and A(x+dx) surface,

Slat(x,dx)Pnx+P(A+dA)PA=PVnx=0P(A+dA)+PA=PdA

Thus the equations become

{0=d(ρuA)=00=d(ρu2A)+AdP=ρuAdu+AdP ,

having used the mass equation to simplify the first term in the momentum equation, since

d(ρu2A)=d(ρuA)=0u+ρuAdu .

Now, from momentum equation a relation from changes in velocity and pressure holds

ρudu=dP .

If the flow is isentropic (i.e. negligible viscous effects, no heat conduction, no shocks), s=s, the definition of the speed of sound

c2(ρ,s)=(Pρ)s ,

can be used to write a relation between changes in pressure and density dP=c2dρ. Using this formula and mass equation d(ρuA)=0, it’s possible to write relations between physical quantities and the variation of the section of the stream tube. Mass equation can be recast as

0=dAA+dρρ+duu==dAA+1c2dPρ+duu==dAA+1c2dPρdPρu2

and thus, introducing Mach number M:=uc,

dAA=(1M2)dPρu2=(1M2)duu .

From this equation, it’s immediate to realize that:

  • for subsonic flows M<1, if section of the stream tube increases, dA>0, thus velocity descreases du<0, and pressure increases dP>0, and viceversa

  • for supersonic flows M>1, if section of the stream tube increases, dA>0, thus velocity increases du>0, and pressure decreases dP<0