Peeter Joot's (OLD) Blog.

Math, physics, perl, and programming obscurity.

PHY454H1S Continuum Mechanics. Lecture 12: Flow in a pipe. Gravity driven flow of a film. Taught by Prof. K. Das.

Posted by peeterjoot on February 17, 2012

[Click here for a PDF of this post with nicer formatting and figures if the post had any (especially if my latex to wordpress script has left FORMULA DOES NOT PARSE errors.)]


Peeter’s lecture notes from class. May not be entirely coherent.

Review. Steady rectilinear flow.


\begin{aligned}\frac{\partial {}}{\partial {t}} = 0\end{aligned} \hspace{\stretch{1}}(2.1)

Rectilinear is a unidirectional flow such as

\begin{aligned}\mathbf{u} = \hat{\mathbf{x}} u( x, y, z ),\end{aligned} \hspace{\stretch{1}}(2.2)

Utilizing an incompressibility assumption \boldsymbol{\nabla} \cdot \mathbf{u} = 0, so for this case we have

\begin{aligned}\frac{\partial {u}}{\partial {x}} = 0\end{aligned}


\begin{aligned}u = u(y, z)\end{aligned}

Note that Prof. Das called this a continuity requirement, and justified this label with the relation

\begin{aligned}\frac{d\rho}{dt} = \rho (\boldsymbol{\nabla} \cdot \mathbf{u}),\end{aligned} \hspace{\stretch{1}}(2.3)

which was a consequence of mass conservation. It’s still not clear to me why he would call this a continuity requirement.

\item Nonlinear term is zero. (\mathbf{u} \cdot \boldsymbol{\nabla}) \mathbf{u} = 0
\item p = p(x). Since \frac{d^2 p}{dx^2} = 0 we also have \frac{dp}{dx} = -G, a constant.

\item \mu \left( \frac{\partial^2 {{u}}}{\partial {{y}}^2} + \frac{\partial^2 {{u}}}{\partial {{z}}^2} \right) = G


Solution by intuition.

Two examples that we have solved analytically are illustrated in figure (\ref{fig:continuumL12:continuumL12fig1}) and figure (\ref{fig:continuumL12:continuumL12fig2})

\caption{Simple shear flow}
\caption{Channel flow}

Sometimes we can utilize solutions already found to understand the behaviour of more complex systems. Combining the two we can look at flow over a plate as in figure (\ref{fig:continuumL12:continuumL12fig3})

\caption{Flow on a plate}

Example 2. Fluid in a container. If the surface tension is altered on one side, we induce a flow on the surface, leading to a circulation flow. This can be done for example, by introducing a heat source or addition of surfactant.

This is illustrated in figure (\ref{fig:continuumL12:continuumL12fig4})
\caption{Circulation flow induced by altering surface tension.}

This sort of flow is hard to analyze, only first done by Steve Davis in the 1980’s. The point here is that we can use some level of intuition to guide our attempts at solution.

Flow down a pipe.

Reading: section 2 from [1].

Recall that the Navier-Stokes equation is

\begin{aligned}\boxed{\rho \frac{\partial {\mathbf{u}}}{\partial {t}} + \rho (\mathbf{u} \cdot \boldsymbol{\nabla}) \mathbf{u} = - \boldsymbol{\nabla} p + \mu \boldsymbol{\nabla}^2 \mathbf{u} + \rho \mathbf{f}.}\end{aligned} \hspace{\stretch{1}}(4.4)

We need to express this in cylindrical coordinates (r, \theta, z) as in figure (\ref{fig:continuumL12:continuumL12fig5})
\caption{Flow through a pipe.}

Our gradient is

\begin{aligned}\boldsymbol{\nabla} = \hat{\mathbf{r}} \frac{\partial {}}{\partial {r}} + \frac{\hat{\boldsymbol{\theta}}}{r} \frac{\partial {}}{\partial {\theta}} + \hat{\mathbf{z}} \frac{\partial {}}{\partial {z}},\end{aligned} \hspace{\stretch{1}}(4.5)

For our Laplacian we find

\begin{aligned}\boldsymbol{\nabla}^2 &= \left(\hat{\mathbf{r}} \frac{\partial {}}{\partial {r}} + \frac{\hat{\boldsymbol{\theta}}}{r} \frac{\partial {}}{\partial {\theta}} + \hat{\mathbf{z}} \frac{\partial {}}{\partial {z}}\right) \cdot\left(\hat{\mathbf{r}} \frac{\partial {}}{\partial {r}} + \frac{\hat{\boldsymbol{\theta}}}{r} \frac{\partial {}}{\partial {\theta}} + \hat{\mathbf{z}} \frac{\partial {}}{\partial {z}}\right) \\ &=\partial_{rr} + \frac{\hat{\boldsymbol{\theta}}}{r} \cdot (\partial_\theta \hat{\mathbf{r}}) \partial_r+ \frac{1}{{r}} \partial_\theta \left( \frac{1}{{r}} \partial_\theta \right)+ \partial_{zz} \\ &=\partial_{rr} + \frac{1}{{r}} \partial_r + \frac{1}{{r^2}} \partial_{\theta\theta} + \partial_{zz},\end{aligned}

which we can write as

\begin{aligned}\boldsymbol{\nabla}^2 = \frac{1}{{r}} \frac{\partial {}}{\partial {r}} \left( r \frac{\partial {}}{\partial {r}} \right) + \frac{1}{{r^2}} \frac{\partial^2 {{}}}{\partial {{\theta}}^2} + \frac{\partial^2 {{}}}{\partial {{z}}^2}.\end{aligned} \hspace{\stretch{1}}(4.6)

NS takes the form

\begin{aligned}\boxed{\begin{aligned}\rho \frac{\partial {\mathbf{u}}}{\partial {t}} &+ \rho \left(u_r \frac{\partial {}}{\partial {r}} + \frac{u_\theta}{r} \frac{\partial {}}{\partial {\theta}} + u_z \frac{\partial {}}{\partial {z}} \right) \mathbf{u} =  \\ &- \left(\hat{\mathbf{r}} \frac{\partial {}}{\partial {r}} + \frac{\hat{\boldsymbol{\theta}}}{r} \frac{\partial {}}{\partial {\theta}} + \hat{\mathbf{z}} \frac{\partial {}}{\partial {z}}\right)p + \mu \left(\frac{1}{{r}} \frac{\partial {}}{\partial {r}} \left( r \frac{\partial {}}{\partial {r}} \right) + \frac{1}{{r^2}} \frac{\partial^2 {{}}}{\partial {{\theta}}^2} + \frac{\partial^2 {{}}}{\partial {{z}}^2} \right)\mathbf{u} + \rho \mathbf{f}.\end{aligned}}\end{aligned} \hspace{\stretch{1}}(4.7)

For steady state and incompressible fluids in the absence of body forces we have

\begin{aligned}\left(\hat{\mathbf{r}} \frac{\partial {}}{\partial {r}} + \frac{\hat{\boldsymbol{\theta}}}{r} \frac{\partial {}}{\partial {\theta}} + \hat{\mathbf{z}} \frac{\partial {}}{\partial {z}}\right)p = \mu \left(\frac{1}{{r}} \frac{\partial {}}{\partial {r}} \left( r \frac{\partial {}}{\partial {r}} \right) + \frac{1}{{r^2}} \frac{\partial^2 {{}}}{\partial {{\theta}}^2} + \frac{\partial^2 {{}}}{\partial {{z}}^2} \right)\mathbf{u},\end{aligned} \hspace{\stretch{1}}(4.8)

or, in coordinates

\begin{aligned}\frac{\partial {p}}{\partial {r}}  &= \mu \left(\frac{1}{{r}} \frac{\partial {}}{\partial {r}} \left( r \frac{\partial {}}{\partial {r}} \right) + \frac{1}{{r^2}} \frac{\partial^2 {{}}}{\partial {{\theta}}^2} + \frac{\partial^2 {{}}}{\partial {{z}}^2} \right)u_r \\ \frac{1}{r} \frac{\partial {p}}{\partial {\theta}}&= \mu \left(\frac{1}{{r}} \frac{\partial {}}{\partial {r}} \left( r \frac{\partial {}}{\partial {r}} \right) + \frac{1}{{r^2}} \frac{\partial^2 {{}}}{\partial {{\theta}}^2} + \frac{\partial^2 {{}}}{\partial {{z}}^2} \right)u_\theta \\ \frac{\partial {p}}{\partial {z}}&= \mu \left(\frac{1}{{r}} \frac{\partial {}}{\partial {r}} \left( r \frac{\partial {}}{\partial {r}} \right) + \frac{1}{{r^2}} \frac{\partial^2 {{}}}{\partial {{\theta}}^2} + \frac{\partial^2 {{}}}{\partial {{z}}^2} \right)u_z\end{aligned} \hspace{\stretch{1}}(4.9)

With an assumption that we have no radial or circulatory flows (u_r = u_\theta = 0), and with u_z = w assumed to only have a radial dependence, our velocity is

\begin{aligned}\mathbf{u} = \hat{\mathbf{z}} w(r),\end{aligned} \hspace{\stretch{1}}(4.12)

and an assumption of linear pressure dependence

\begin{aligned}\frac{dp}{dz} = -G,\end{aligned} \hspace{\stretch{1}}(4.13)

then NS takes the final simple form

\begin{aligned}\frac{1}{{r}} \frac{d}{dr} \left( r \frac{dw}{dr} \right) = - \frac{G}{\mu}.\end{aligned} \hspace{\stretch{1}}(4.14)

Solving this we have

\begin{aligned}r \frac{dw}{dr} = - \frac{G r^2}{2\mu} + A\end{aligned} \hspace{\stretch{1}}(4.15)

\begin{aligned}w = -\frac{G r^2}{4 \mu} + A \ln(r) + B\end{aligned} \hspace{\stretch{1}}(4.16)

Requiring finite solutions for r = 0 means that we must have A = 0. Also w(a) = 0, we have B = G a^2/4 \mu so we must have

\begin{aligned}w(r) = \frac{G}{4 \mu}( a^2 - r^2 )\end{aligned} \hspace{\stretch{1}}(4.17)

Example: Gravity driven flow of a liquid film

(This is one of our Professor’s favorite problems).

Coordinates as in figure (\ref{fig:continuumL12:continuumL12fig6})
\caption{Gravity driven flow down an inclined plane.}

\begin{aligned}\mathbf{u} = \hat{\mathbf{x}} u(y)\end{aligned} \hspace{\stretch{1}}(5.18)

Boundary conditions

\item u(y = 0) = 0
\item Tangential stress at the air-liquid interface y = h is equal.

\begin{aligned}\boldsymbol{\tau} \cdot (\boldsymbol{\sigma}_l \cdot \hat{\mathbf{n}}) = \boldsymbol{\tau} \cdot (\boldsymbol{\sigma}_a \cdot \hat{\mathbf{n}}),\end{aligned} \hspace{\stretch{1}}(5.19)


We write

\begin{aligned}\boldsymbol{\tau} &= \begin{bmatrix}1 \\ 0 \\ 0\end{bmatrix} \\ \hat{\mathbf{n}} &= \begin{bmatrix}0 \\ 1 \\ 0\end{bmatrix} \\ \end{aligned} \hspace{\stretch{1}}(5.20)

and seek simultaneous solutions to the pair of stress tensor equations

\begin{aligned}\sigma_{ij}^l &= - p \delta_{ij} + \mu^l \left( \frac{\partial {u_i}}{\partial {x_j}} +\frac{\partial {u_j}}{\partial {x_i}}\right) \\ \sigma_{ij}^a &= - p \delta_{ij} + \mu^a \left( \frac{\partial {u_i}}{\partial {x_j}} +\frac{\partial {u_j}}{\partial {x_i}}\right).\end{aligned} \hspace{\stretch{1}}(5.23)

In general this requires an iterated approach, solving for one with an initial approximation of the other, then switching and tuning the numerical method carefully for convergence.

We expect that the flow of liquid will induce a flow of air at the interface, but may be able to make a one-sided approximation. Let’s see how far we get before we have to introduce any approximations and compute the traction vector for the liquid

\begin{aligned}\boldsymbol{\sigma}^l \cdot \hat{\mathbf{n}} &= \begin{bmatrix}-p & \mu^l {\partial {u}}/{\partial {y}} & 0 \\ \mu^l {\partial {u}}/{\partial {y}} & -p & 0 \\ 0 & 0 & 0\end{bmatrix}\begin{bmatrix}0 \\ 1 \\ 0\end{bmatrix} \\ &=\begin{bmatrix}\mu^l {\partial {u}}/{\partial {y}} \\ -p \\ 0\end{bmatrix}\end{aligned}


\begin{aligned}\boldsymbol{\tau} \cdot (\boldsymbol{\sigma}^l \cdot \hat{\mathbf{n}})=\begin{bmatrix}1 & 0 & 0\end{bmatrix}\begin{bmatrix}\mu^l {\partial {u}}/{\partial {y}} \\ -p \\ 0\end{bmatrix}=\mu^l \frac{\partial {u}}{\partial {y}}\end{aligned} \hspace{\stretch{1}}(5.25)

Our boundary value condition is therefore

\begin{aligned}{\left.{{\mu^l \frac{\partial {u^l}}{\partial {y}}}}\right\vert}_{{y = h}} ={\left.{{\mu^a \frac{\partial {u^a}}{\partial {y}}}}\right\vert}_{{y = h}}\end{aligned} \hspace{\stretch{1}}(5.26)

When can we decouple this, treating only the liquid? Observe that we have

\begin{aligned}{\left.{{\frac{\partial {u^l}}{\partial {y}}}}\right\vert}_{{y = h}} ={\left.{{\frac{\mu^a}{\mu^l} \frac{\partial {u^a}}{\partial {y}}}}\right\vert}_{{y = h}}\end{aligned} \hspace{\stretch{1}}(5.27)

so if

\begin{aligned}\frac{\mu_a}{\mu_l} \ll 1\end{aligned} \hspace{\stretch{1}}(5.28)

we can treat only the liquid portion of the problem, with a boundary value condition

\begin{aligned}{\left.{{\frac{\partial {u^l}}{\partial {y}}}}\right\vert}_{{y = h}} = 0.\end{aligned} \hspace{\stretch{1}}(5.29)

Let’s look at the component of the traction vector in the direction of the normal (liquid pressure acting on the air)

\begin{aligned}\hat{\mathbf{n}} \cdot (\boldsymbol{\sigma}^l \cdot \hat{\mathbf{n}}) = \hat{\mathbf{n}} \cdot (\boldsymbol{\sigma}^a \cdot \hat{\mathbf{n}}) \end{aligned} \hspace{\stretch{1}}(5.30)


\begin{aligned}\begin{bmatrix}0 & 1 & 0\end{bmatrix}\begin{bmatrix}\mu^l \frac{\partial {u}}{\partial {y}} \\ -p^l \\ 0\end{bmatrix}= -{\left.{{p^l}}\right\vert}_{{y = h}} = -{\left.{{p^a}}\right\vert}_{{y = h}}\end{aligned} \hspace{\stretch{1}}(5.31)

i.e. We have pressure matching at the interface.

Our body force is

\begin{aligned}\mathbf{f} = \begin{bmatrix}g \sin\alpha \\ -g \cos\alpha \\ 0\end{bmatrix}\end{aligned} \hspace{\stretch{1}}(5.32)

Referring to the Navier-Stokes equation 4.4, we see that our only surviving parts are


\begin{aligned}0 = -\frac{\partial {p}}{\partial {x}} + \mu \frac{\partial^2 {{u}}}{\partial {{y}}^2} + \rho g \sin\alpha \end{aligned} \hspace{\stretch{1}}(5.33a)

\begin{aligned}0 = -\frac{\partial {p}}{\partial {y}} - \rho g \cos\alpha \end{aligned} \hspace{\stretch{1}}(5.33b)

\begin{aligned}0 = -\frac{\partial {p}}{\partial {z}} \end{aligned} \hspace{\stretch{1}}(5.33c)


The last gives us p \ne p(z). Integrating the second we have

\begin{aligned}p = \rho g y \cos\alpha + p_1\end{aligned} \hspace{\stretch{1}}(5.34)

Since p = p_{\text{atm}} at y = h, we have

\begin{aligned}p_{\text{atm}} = \rho g h \cos\alpha + p_1\end{aligned} \hspace{\stretch{1}}(5.35)

Our first NS equation 5.33a becomes

\begin{aligned}0 = \mu \frac{\partial^2 {{u}}}{\partial {{y}}^2} + g \sin\alpha,\end{aligned} \hspace{\stretch{1}}(5.36)


\begin{aligned}\frac{\partial^2 {{u}}}{\partial {{y}}^2} = -\frac{g}{\mu} \sin\alpha\end{aligned} \hspace{\stretch{1}}(5.37)

Solving we have

\begin{aligned}u = - \rho g \frac{\sin\alpha}{2 \mu} y^2 + A y + B\end{aligned} \hspace{\stretch{1}}(5.38)


\begin{aligned}u(0) &= 0 \\ {\left.{{\frac{\partial {u}}{\partial {y}}}}\right\vert}_{{y = h}} &= 0\end{aligned} \hspace{\stretch{1}}(5.39)

\begin{aligned}u = \rho g \frac{\sin\alpha}{2 \mu} \left( 2 h y - y^2 \right) .\end{aligned} \hspace{\stretch{1}}(5.41)

This velocity distribution is illustrated figure (\ref{fig:continuumL12:continuumL12fig7}).

\caption{Velocity streamlines for flow down a plane.}

It’s important to note that in these problems we have to derive our boundary value conditions! They are not given.

In this discussion, the height h was assumed to be constant, with the tangential direction constant and parallel to the surface that the liquid is flowing on. It’s claimed in class that this is actually a consequence of surface tension only! That’s not at all intuitive, but will be covered when we learn about “stability conditions”.

Study note.

Memorizing the NS equation is required for midterm, but more complex stuff (like cylindrical forms of the strain tensor if required) will be given.


[1] D.J. Acheson. Elementary fluid dynamics. Oxford University Press, USA, 1990.


Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s

%d bloggers like this: