Paper 4, Section II, B

Waves
Part II, 2020

(a) Show that the equations for one-dimensional unsteady flow of an inviscid compressible fluid at constant entropy can be put in the form

(t+(u±c)x)R±=0\left(\frac{\partial}{\partial t}+(u \pm c) \frac{\partial}{\partial x}\right) R_{\pm}=0

where uu and cc are the fluid velocity and the local sound speed, respectively, and the Riemann invariants R±R_{\pm}are to be defined.

Such a fluid occupies a long narrow tube along the xx-axis. For times t<0t<0 it is at rest with uniform pressure p0p_{0}, density ρ0\rho_{0} and sound speed c0c_{0}. At t=0t=0 a finite segment, 0xL0 \leqslant x \leqslant L, is disturbed so that u=U(x)u=U(x) and c=c0+C(x)c=c_{0}+C(x), with U=C=0U=C=0 for x0x \leqslant 0 and xLx \geqslant L. Explain, with the aid of a carefully labelled sketch, how two independent simple waves emerge after some time. You may assume that no shock waves form.

(b) A fluid has the adiabatic equation of state

p(ρ)=AB2ρp(\rho)=A-\frac{B^{2}}{\rho}

where AA and BB are positive constants and ρ>B2/A\rho>B^{2} / A.

(i) Calculate the Riemann invariants for this fluid, and express u±cu \pm c in terms of R±R_{\pm} and c0c_{0}. Deduce that in a simple wave with R=0R_{-}=0 the velocity field translates, without any nonlinear distortion, at the equilibrium sound speed c0c_{0}.

(ii) At t=0t=0 this fluid occupies x>0x>0 and is at rest with uniform pressure, density and sound speed. For t>0t>0 a piston initially at x=0x=0 executes simple harmonic motion with position x(t)=asinωtx(t)=a \sin \omega t, where aω<c0a \omega<c_{0}. Show that u(x,t)=U(ϕ)u(x, t)=U(\phi), where ϕ=ω(tx/c0)\phi=\omega\left(t-x / c_{0}\right), for some function UU that is zero for ϕ<0\phi<0 and is 2π2 \pi-periodic, but not simple harmonic, for ϕ>0\phi>0. By approximately inverting the relationship between ϕ\phi and the time τ\tau that a characteristic leaves the piston for the case ϵ=aω/c01\epsilon=a \omega / c_{0} \ll 1, show that

U(ϕ)=aω(cosϕϵsin2ϕ32ϵ2sin2ϕcosϕ+O(ϵ3)) for ϕ>0U(\phi)=a \omega\left(\cos \phi-\epsilon \sin ^{2} \phi-\frac{3}{2} \epsilon^{2} \sin ^{2} \phi \cos \phi+O\left(\epsilon^{3}\right)\right) \quad \text { for } \quad \phi>0