1.II.17B

Fluid Dynamics
Part IB, 2008

Two incompressible fluids flow in infinite horizontal streams, the plane of contact being z=0z=0, with zz positive upwards. The flow is given by

U(r)={U2e^x,z>0U1e^x,z<0\mathbf{U}(\mathbf{r})= \begin{cases}U_{2} \hat{\mathbf{e}}_{x}, & z>0 \\ U_{1} \hat{\mathbf{e}}_{x}, & z<0\end{cases}

where e^x\hat{\mathbf{e}}_{x} is the unit vector in the positive xx direction. The upper fluid has density ρ2\rho_{2} and pressure p0gρ2zp_{0}-g \rho_{2} z, the lower has density ρ1\rho_{1} and pressure p0gρ1zp_{0}-g \rho_{1} z, where p0p_{0} is a constant and gg is the acceleration due to gravity.

(i) Consider a perturbation to the flat surface z=0z=0 of the form

zζ(x,y,t)=ζ0ei(kx+y)+st.z \equiv \zeta(x, y, t)=\zeta_{0} e^{i(k x+\ell y)+s t} .

State the kinematic boundary conditions on the velocity potentials ϕi\phi_{i} that hold on the interface in the two domains, and show by linearising in ζ\zeta that they reduce to

ϕiz=ζt+Uiζx(z=0,i=1,2).\frac{\partial \phi_{i}}{\partial z}=\frac{\partial \zeta}{\partial t}+U_{i} \frac{\partial \zeta}{\partial x} \quad(z=0, i=1,2) .

(ii) State the dynamic boundary condition on the perturbed interface, and show by linearising in ζ\zeta that it reduces to

ρ1(U1ϕ1x+ϕ1t+gζ)=ρ2(U2ϕ2x+ϕ2t+gζ)(z=0)\rho_{1}\left(U_{1} \frac{\partial \phi_{1}}{\partial x}+\frac{\partial \phi_{1}}{\partial t}+g \zeta\right)=\rho_{2}\left(U_{2} \frac{\partial \phi_{2}}{\partial x}+\frac{\partial \phi_{2}}{\partial t}+g \zeta\right) \quad(z=0)

(iii) Use the velocity potentials

ϕ1=U1x+A1eqzei(kx+y)+st,ϕ2=U2x+A2eqzei(kx+y)+st,\phi_{1}=U_{1} x+A_{1} e^{q z} e^{i(k x+\ell y)+s t}, \quad \phi_{2}=U_{2} x+A_{2} e^{-q z} e^{i(k x+\ell y)+s t},

where q=k2+2q=\sqrt{k^{2}+\ell^{2}}, and the conditions in (i) and (ii) to perform a stability analysis. Show that the relation between s,ks, k and \ell is

s=ikρ1U1+ρ2U2ρ1+ρ2±[k2ρ1ρ2(U1U2)2(ρ1+ρ2)2qg(ρ1ρ2)ρ1+ρ2]1/2.s=-i k \frac{\rho_{1} U_{1}+\rho_{2} U_{2}}{\rho_{1}+\rho_{2}} \pm\left[\frac{k^{2} \rho_{1} \rho_{2}\left(U_{1}-U_{2}\right)^{2}}{\left(\rho_{1}+\rho_{2}\right)^{2}}-\frac{q g\left(\rho_{1}-\rho_{2}\right)}{\rho_{1}+\rho_{2}}\right]^{1 / 2} .

Find the criterion for instability.