Past exam of the fluids and waves course of the University of Oxford / 2021 / 1

The velocity of a fluid element is the time derivative of its position with its Lagrangian coordinates $X,Y$ held fixed. Direct differentiation givesReplacing the Lagrangian coordinates by the current Eulerian coordinates therefore gives the velocity fieldThe distinction between following a fluid element and differentiating at fixed Eulerian coordinates is what connects the given deformation map to this velocity field.

For a planar incompressible flow, the divergence of the velocity field must vanish. HereEquivalently, the Jacobian determinant of the deformation map isso every material area is preserved. However, the out-of-plane vorticity isSince $\Omega >0$, this vorticity never vanishes. The motion is an incompressible flow but is not an irrotational flow; its local spin is clockwise.

Adopt the streamfunction convention $u={\psi }_y$, $v=-{\psi }_x$. Integrating the first component of the velocity field in $y$ gives $\psi =-xy+\frac{1}{2}\Omega e^{-2t}{y}^2+f(x,t)$. The second component then requires $f_x=\Omega e^{2t}x$. Thus, omitting an arbitrary function of time,At a fixed instant, the streamlines are the connected curves $\psi =K$. To classify this quadratic form, introduce the invertible scaled and rotated coordinatesThenFor $\Omega >1$, both coefficients are positive. Every $K>0$ therefore gives an ellipse, and $K=0$ is the stagnation point at the origin. For $0<\Omega <1$, the quadratic form is an indefinite quadratic form: each $K\ne 0$ gives a hyperbola. The level $K=0$ consists of the straight separatriceswhose rays approach or leave the origin; the origin itself is a stagnation point. At the transition $\Omega =1$, $\psi =d^2$, so the nonstationary streamlines are parallel straight lines and $d=0$ is a line of stagnation points. These are instantaneous streamlines; the time dependence means that they need not coincide with pathlines.

Use $w=\varphi +i\psi$ for the complex potential of a two-dimensional potential flow, so its complex velocity is $q=w^{\prime }=u-iv$. The incident uniform flow and line source haveThe coefficient of the complex logarithm is one because a line source of strength $Q$ has complex potential of a two-dimensional potential flow $(Q/(2\pi ))\log (z-b)$, and here $Q=2\pi$. Applying the Milne-Thomson circle theorem to the unit circle givesSince $1/z-b=-b(z-1/b)/z$, locally consistent branches of the complex logarithm giveChanges of branch only alter the constant in a local expression; $w^{\prime }$ is single valued away from its singularities. The added terms represent a two-dimensional potential doublet at the origin, an image line source at $1/b$, and a compensating line sink at the origin. The image line source and line sink have equal and opposite strengths, so no net volume flux enters the cylinder.

On $|z|=1$, the Milne-Thomson circle theorem makes $w=f(z)+\overline{f(z)}$ up to a constant. Consequently the streamfunction is constant there, establishing the no-penetration boundary condition. All line source strengths are real, and the two-dimensional potential doublet has zero circulation, so the circulation around the cylinder is zero. There is no additional line vortex term.

Let $F_x,F_y$ denote the force per unit length exerted by the fluid on the cylinder, and let $\rho$ be the constant mass density. For a counterclockwise contour $C$ around the unit circle, the Blasius force theorem givesAn efficient use of the residue theorem finds the sum of the enclosed residues by subtracting the external pole at $b$ from the sum of all finite residues. At large $z$,The sum of all finite residues of $q^2$ is therefore $2U$. At $b$, write $q=(z-b{)}^{-1}+H(z)$, with $H$ a holomorphic function near $b$. ThenThusUsing $a=1/b$ and $1/(b-a)-1/b=1/[b(b^2-1)]$ givesThe vanishing vertical force also follows from reflection symmetry about the real axis. This force on a unit cylinder due to a uniform flow and an external source uses the counterclockwise orientation in the Blasius force theorem to fix the sign of the horizontal force.

Since $b>1$, the denominator inis positive. The positive root of $U{b}^2-b-U=0$ isThe other root is negative, so the horizontal force points in the positive $x$ direction for $1<b<b_{\ast }$, vanishes at $b=b_{\ast }$, and points in the negative $x$ direction for $b>b_{\ast }$. In particular, the nearby line source pulls the cylinder towards itself as $b\downarrow 1$, while at large $b$ the horizontal force is $F_x\sim -2\pi \rho U/b^2$.