Past exam of the viscous flow course of the University of Oxford / 2023 / 1 / b

For the unidirectional flow $\mathbf{u}=(0,0,w)$, incompressibility gives ${\partial }_{z}w=0$. Hence $w=w(x,y,t)$, and the convective acceleration $w{\partial }_{z}w$ vanishes. The transverse components of the Navier-Stokes equation give $p_x=p_y=0$, while the axial component givesA steady-flow assumption is implicit in the question. Neither unidirectional flow nor a time-independent pressure gradient alone eliminates ${\partial }_{t}w$: a decaying transient can be superposed on the final pipe profile. The native PDF does not explicitly say steady flow in this part. For the intended steady flow, $w_t=0$, and thereforeTaking the pipe walls to be stationary, the no-slip boundary condition isThe no-penetration boundary condition is already satisfied by the assumed velocity field. This is the Poisson equation defining rectangular duct flow.

For a profile depending only on $y$, the Poisson equation becomes $w_1^{\prime \prime }=-G/\mu$. Integrating twice givesThe no-slip boundary conditions $w_1(b)=w_1(-b)=0$ imply $C_1=0$ and $C_2=G{b}^2/(2\mu )$. ThusThis is plane Poiseuille flow between the two horizontal walls. It is positive in the interior for $G>0$, so the velocity points down the pressure gradient. It does not satisfy the side-wall no-slip boundary conditions, which is precisely what the correction in part (iii) repairs.

Subtract the particular plane Poiseuille flow profile from the full rectangular duct flow. The remainder is a harmonic function:The data are even in $x$ and $y$. The boundary value problem is unique: if two solutions differed by $v$, Green's first identity and the zero Dirichlet boundary condition would giveThus $v$ is constant, and its zero boundary values force $v=0$. Reflection of a solution in either coordinate axis therefore gives the same solution. We may seek even separated eigenfunctions.

Use separation of variables withEach summand solves the Laplace equation, vanishes at $y=\pm b$, and has the required evenness. The hyperbolic cosine ratio equals one at either side wall. The functions $\cos (k_{n}y)$ have orthogonality relationThey are the even Dirichlet eigenfunctions on $[-b,b]$, so their Fourier series represents the even boundary data. Projecting that data onto the $n$th eigenfunction yields the Fourier coefficientFor completeness, two integrations by parts give the quadratic arch expansion in half-integer cosine modes: the displayed integral is $4(-1{)}^n/k_n^3$, agreeing with the identity supplied in the paper.

Adding the particular solution givesThe coefficients are absolutely summable and the hyperbolic cosine ratio is at most one throughout the closed rectangle, giving uniform convergence and the prescribed boundary values. On compact subsets of the interior, the hyperbolic cosine ratio supplies exponential decay in $n$, justifying differentiation of the series and the Laplace equation. The side-wall correction slows the plane Poiseuille flow enough to satisfy all four no-slip boundary conditions.

The volumetric flow rate is the integral of axial velocity over the cross-section:The plane Poiseuille flow contribution isFor each correction eigenfunction,Uniform convergence permits termwise integration. Multiplying these two integrals by $A_n$ givesConsequentlyEvery summand in the subtracted quantity is positive. Thus the side walls reduce the flux below the infinite-width plane Poiseuille flow prediction. Its dimensions are volume per unit time, as required for a volumetric flow rate.