← 2017 Paper 2

UPSC 2017 Maths Optional Paper 2 Q8c — Step-by-Step Solution

15 marks · Section B

Potential flow · Mechanics & Fluid Dynamics · asked 10× in 13 yrs · Read the full method →

Question

If the velocity of an incompressible fluid at the point (x,y,z)(x,y,z) is given by (3xzr5,3yzr5,3z2r2r5), r2=x2+y2+z2\left(\dfrac{3xz}{r^5},\dfrac{3yz}{r^5},\dfrac{3z^2-r^2}{r^5}\right),\ r^2=x^2+y^2+z^2, then prove that the liquid motion is possible and that the velocity potential is zr3\dfrac{z}{r^3}. Further, determine the streamlines.

Technique

Continuity q=0\nabla\cdot\mathbf q=0 for “motion possible”; verify q=ϕ\mathbf q=-\nabla\phi for the stated potential; streamlines from dx/u=dy/v=dz/wdx/u=dy/v=dz/w — the x,yx,y pair gives meridian planes, the Stokes stream function ψ=ρ2/r3\psi=\rho^2/r^3 gives the in-plane curves.

Solution

Setup. Velocity components, with r=x2+y2+z2r=\sqrt{x^2+y^2+z^2}:

u=3xzr5,v=3yzr5,w=3z2r2r5.u=\frac{3xz}{r^5},\qquad v=\frac{3yz}{r^5},\qquad w=\frac{3z^2-r^2}{r^5}.

Step 1 — Motion is possible (incompressibility: q=0\nabla\cdot\mathbf q=0)

For an incompressible fluid the motion is possible iff the equation of continuity q=ux+vy+wz=0\nabla\cdot\mathbf q=u_x+v_y+w_z=0 holds. Using rx=xr\dfrac{\partial r}{\partial x}=\dfrac{x}{r} etc.:

ux=x(3xzr5)=3zr5+3xz(5)r6xr=3zr515x2zr7,u_x=\frac{\partial}{\partial x}\Big(3xz\,r^{-5}\Big)=3z\,r^{-5}+3xz\,(-5)r^{-6}\frac{x}{r}=\frac{3z}{r^5}-\frac{15x^2z}{r^7}, vy=3zr515y2zr7,v_y=\frac{3z}{r^5}-\frac{15y^2z}{r^7}, wz=z[(3z2r2)r5]=(6z2z)r5+(3z2r2)(5)r6zr=4zr55z(3z2r2)r7.w_z=\frac{\partial}{\partial z}\Big[(3z^2-r^2)r^{-5}\Big]=\frac{(6z-2z)}{r^5}+(3z^2-r^2)(-5)r^{-6}\frac{z}{r}=\frac{4z}{r^5}-\frac{5z(3z^2-r^2)}{r^7}.

Adding and using x2+y2+z2=r2x^2+y^2+z^2=r^2:

q=3z+3z+4zr515z(x2+y2)+5z(3z2r2)r7=10zr515z(r2z2)+15z35zr2r7.\nabla\cdot\mathbf q=\frac{3z+3z+4z}{r^5}-\frac{15z(x^2+y^2)+5z(3z^2-r^2)}{r^7}=\frac{10z}{r^5}-\frac{15z(r^2-z^2)+15z^3-5zr^2}{r^7}.

The numerator of the second term is 15zr215z3+15z35zr2=10zr215zr^2-15z^3+15z^3-5zr^2=10zr^2, so the second term is 10zr2r7=10zr5\dfrac{10zr^2}{r^7}=\dfrac{10z}{r^5}, giving

  q=10zr510zr5=0.  \boxed{\;\nabla\cdot\mathbf q=\frac{10z}{r^5}-\frac{10z}{r^5}=0.\;}

The continuity equation is satisfied, so the liquid motion is possible.

Step 2 — Velocity potential ϕ=zr3\phi=\dfrac{z}{r^3}

Flow has a velocity potential if it is irrotational; with the convention q=ϕ\mathbf q=-\nabla\phi we verify ϕ=zr3\phi=\dfrac{z}{r^3} reproduces the velocity. With x(zr3)=z(3)r4xr=3xzr5\dfrac{\partial}{\partial x}\Big(\dfrac{z}{r^3}\Big)=z(-3)r^{-4}\dfrac{x}{r}=-\dfrac{3xz}{r^5}:

ϕx=3xzr5=u,ϕy=3yzr5=v,-\frac{\partial\phi}{\partial x}=\frac{3xz}{r^5}=u,\qquad -\frac{\partial\phi}{\partial y}=\frac{3yz}{r^5}=v, ϕz=(1r3+z(3)r4zr)=1r3+3z2r5=3z2r2r5=w.-\frac{\partial\phi}{\partial z}=-\Big(\frac{1}{r^3}+z(-3)r^{-4}\frac{z}{r}\Big)=-\frac{1}{r^3}+\frac{3z^2}{r^5}=\frac{3z^2-r^2}{r^5}=w.

All three match, so

  q=ϕ,ϕ=zr3.  \boxed{\;\mathbf q=-\nabla\phi,\qquad \phi=\frac{z}{r^3}.\;}

(Since ϕ\phi exists, the flow is irrotational; one checks 2ϕ=0\nabla^2\phi=0, consistent with Step 1 — this is a doublet along the zz-axis.)

Step 3 — Streamlines

Streamlines satisfy dxu=dyv=dzw\dfrac{dx}{u}=\dfrac{dy}{v}=\dfrac{dz}{w}, i.e.

dx3xz/r5=dy3yz/r5=dz(3z2r2)/r5.\frac{dx}{3xz/r^5}=\frac{dy}{3yz/r^5}=\frac{dz}{(3z^2-r^2)/r^5}.

First relation (from the xx- and yy-ratios):

dxx=dyy  lnxlny=const  yx=const.\frac{dx}{x}=\frac{dy}{y}\ \Longrightarrow\ \ln x-\ln y=\text{const}\ \Longrightarrow\ \frac{y}{x}=\text{const}.

So each streamline lies in a meridian plane y=(tanβ)xy=(\tan\beta)\,x through the zz-axis (axisymmetric flow).

Second relation. Introduce the cylindrical radius ρ=x2+y2\rho=\sqrt{x^2+y^2}. In a meridian plane the in-plane components are uρ=3ρzr5u_\rho=\dfrac{3\rho z}{r^5}, w=3z2r2r5w=\dfrac{3z^2-r^2}{r^5} (r2=ρ2+z2r^2=\rho^2+z^2). The Stokes stream function ψ\psi with uρ=1ρψz, w=1ρψρu_\rho=-\dfrac1\rho\psi_z,\ w=\dfrac1\rho\psi_\rho is

ψ=ρ2r3(i.e. sin2θr in spherical (r,θ)),\psi=\frac{\rho^2}{r^3}\qquad\Big(\text{i.e. }\frac{\sin^2\theta}{r}\text{ in spherical }(r,\theta)\Big),

which one checks reproduces uρ,wu_\rho,w. Streamlines are ψ=\psi=const:

Answer

  yx=constandρ2r3=x2+y2(x2+y2+z2)3/2=const,  \boxed{\;\frac{y}{x}=\text{const}\quad\text{and}\quad \frac{\rho^2}{r^3}=\frac{x^2+y^2}{(x^2+y^2+z^2)^{3/2}}=\text{const},\;}
We post more of this — worked solutions, CSAT trap breakdowns, guide chapters — a few times a week on Telegram. Free, no sign-in. Join

This solution is part of the Maths Coverage Map — 13 years, mapped. Get the take-away PDF free.