Paper 2, Section II, B

Waves
Part II, 2017

Show that, for a one-dimensional flow of a perfect gas (with γ>1\gamma>1 ) at constant entropy, the Riemann invariants R±=u±2(cc0)/(γ1)R_{\pm}=u \pm 2\left(c-c_{0}\right) /(\gamma-1) are constant along characteristics dx/dt=u±c.d x / d t=u \pm c .

Define a simple wave. Show that in a right-propagating simple wave

ut+(c0+12(γ+1)u)ux=0\frac{\partial u}{\partial t}+\left(c_{0}+\frac{1}{2}(\gamma+1) u\right) \frac{\partial u}{\partial x}=0

In some circumstances, dissipative effects may be modelled by

ut+(c0+12(γ+1)u)ux=αu\frac{\partial u}{\partial t}+\left(c_{0}+\frac{1}{2}(\gamma+1) u\right) \frac{\partial u}{\partial x}=-\alpha u

where α\alpha is a positive constant. Suppose also that uu is prescribed at t=0t=0 for all xx, say u(x,0)=u0(x)u(x, 0)=u_{0}(x). Demonstrate that, unless a shock develops, a solution of the form

u(x,t)=u0(ξ)eαtu(x, t)=u_{0}(\xi) e^{-\alpha t}

can be found, where, for each xx and t,ξt, \xi is determined implicitly as the solution of the equation

xc0t=ξ+γ+12α(1eαt)u0(ξ)x-c_{0} t=\xi+\frac{\gamma+1}{2 \alpha}\left(1-e^{-\alpha t}\right) u_{0}(\xi)

Deduce that, despite the presence of dissipative effects, a shock will still form at some (x,t)(x, t) unless α>αc\alpha>\alpha_{c}, where

αc=12(γ+1)maxu0<0u0(ξ)\alpha_{c}=\frac{1}{2}(\gamma+1) \max _{u_{0}^{\prime}<0}\left|u_{0}^{\prime}(\xi)\right|