[Physics] Spin-up time for fluid in an open cylinder with no endwalls

fluid dynamicsrotational-dynamics

I'm interested in how quickly (if at all) the fluid in an open cylinder would reach solid-body rotation if the cylinder is suspended vertically in an unbounded fluid.

I've found plenty of work that considers spin-up in a cylinder with a single end-wall, but can't seem to track down any information on whether the end-wall is essential for driving the flow. Even a pointer to the correct keywords to search for would help.

Best Answer

This problem can get as complicated as you like, but let me show you how it can be approached in the most simple case. The Navier-Stokes equations describe the flow:

$\rho\left(\displaystyle\frac{\partial {\bf u}}{\partial t}+({\bf u}\cdot\nabla){\bf u}\right)=-\nabla p + \mu \nabla^2 {\bf u}$. $\qquad$ (1)

Assume that the motion of the fluid is slow enough so that the non-linear term $({\bf u}\cdot \nabla){\bf u}$ is negligible. This will be true as long as viscosity forces are large compared to inertial forces from acceleration. Ugh, now that it is a linear equation, something can be done analytically.

Consider the problem in cylindrical coordinates with $z$, $r$ and $\phi$. Let us assume $\partial/\partial z=0$ and $\partial/\partial \phi=0$ for all quantities due to the axial symmetry of the problem. Then the continuity equation will offer additional simplifications. In cylindrical coordinates, the continuity equation reads

$\displaystyle\frac{\partial \rho}{\partial t}+ \frac{1}{r}\displaystyle\frac{\partial(\rho r u_r)}{\partial r}+ \frac{1}{r}\frac{\partial\rho u_{\phi}}{\partial \phi} + \frac{\partial (\rho u_z)}{\partial z}=0$

With $\rho=\mathrm{const}$ (incompressibility assumption), the first term is zero, and last two terms are equal to zero due to the axial symmetry assumption. The remainder yields $u_r=0$. Hooray! Let us choose the problem so that $u_z=0$, too (true if the cylinder is not moving with respect to the fluid in $z$-direction) Then we are left with just $u_{\phi}$, and let us denote it $u\equiv u_{\phi}$.

If the cylinder is suspended vertically, i.e., the axis of the cylinder is along the gravity force, then the pressure gradient $\nabla p$ can also be neglected. This is because the only component of $\nabla p$ is in z-direction due to gravity, and we already assumed $\partial {\bf u}/\partial z$=0. And there may be no $\phi$-component of $\nabla p$, because it would violate our $\partial /{\partial \phi}=0$ assumption.

All of the above leaves us with the only equation for $u_\phi\equiv u$:

$\displaystyle\frac{\partial u}{\partial t} = \mu \frac{1}{r}\frac{\partial}{\partial r}\left(r\frac{\partial u}{\partial r}\right)$. $\qquad$ (2)

To understand where it came from, look up the explicit form of equation (1) in cylindrical coordinates at this Wikipedia page, take the $\phi$-component of it and use the assumptions outlined below to get rid of $\partial/\partial z$ and $\partial/\partial \phi$ terms, and of the non-linear and pressure terms.

The cylinder is spinning, which makes the boundary conditions for equation (2)

$u(r=0, t)=0$ and $u(r=R, t)=R\omega_0$,

where $R$ is the radius and $\omega_0$ is the angular frequency of the cylinder. To understand the second one, look at the well-known solution of a similar problem at this Wikipedia page or this paper referenced in the Wiki page. The initial conditions are

$u(r, t=0)=0$.

Equation (2) is a very well-known equation. It is the diffusion or the heat equation in disguise. Well, no surprise about it, because viscosity in is nothing more than diffusion of momentum. At least, this is how it is expressed in the viscosity term of equation (1). The solution is readily obtained using separation of variables and Fourier analysis. I will not solve it here, because it is a well understood mathematical problem, with which graduate students are routinely tortured with great success. This Wikipedia page describes it.

And the back-of-the envelope estimate I mentioned in my comment comes from dimensional analysis of equation (2). If the typical velocity scale of the problem is $U$, spatial scale is $R$ and time scale is $T$, then (2) can be approximated as

$\displaystyle\frac{U}{T} \approx \mu \frac{1}{R} \frac{1}{R} R \frac{U}{R}$,

or

$T \approx R^2/\mu$,

which is a crude estimate of your spin-up time.

Related Question