Peeter Joot's (OLD) Blog.

Math, physics, perl, and programming obscurity.

PHY454H1S Continuum Mechanics. Lecture 9: Newtonian fluids. Mass conservation. Constitutive relation. Incompressible fluids. Taught by Prof. K. Das.

Posted by peeterjoot on February 10, 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.)]


section 1.4 from [1]. FIXME: Probably more elsewhere too.


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

Review: Relative motion near a point in a fluid

Referring to figure (\ref{fig:continuumL9:continuumL9fig1})
\caption{velocity displacements at a fluid point.}

we write

\begin{aligned}d\mathbf{x}' = d\mathbf{x} + d\mathbf{u} \delta t\end{aligned} \hspace{\stretch{1}}(3.1)

or in coordinate form

\begin{aligned}\begin{aligned}dx_i &= dx_i + du_i \delta t \\ &= dx_i + \frac{\partial {u_i}}{\partial {x_j}} dx_j \delta t \end{aligned}\end{aligned} \hspace{\stretch{1}}(3.2)

We can now split the components of the gradient of u_i into symmetric and antisymmetric parts in the normal way

\begin{aligned}\begin{aligned}\frac{\partial {u_i}}{\partial {x_j}}&= \frac{1}{{2}} \left( \frac{\partial {u_i}}{\partial {x_j}}+\frac{\partial {u_j}}{\partial {x_i}}\right)+\frac{1}{{2}} \left( \frac{\partial {u_i}}{\partial {x_j}}-\frac{\partial {u_j}}{\partial {x_i}}\right) \\ &\equiv e_{ij} + \omega_{ij}.\end{aligned}\end{aligned} \hspace{\stretch{1}}(3.3)

The antisymmetric term (name?)


\begin{aligned}\boldsymbol{\omega} = \boldsymbol{\nabla} \times \mathbf{u},\end{aligned} \hspace{\stretch{1}}(3.4)

we introduce the dual vector

\begin{aligned}\boldsymbol{\Omega} = \Omega_k \mathbf{e}_k = \frac{1}{{2}} \boldsymbol{\omega}\end{aligned} \hspace{\stretch{1}}(3.5)

defined according to

\begin{aligned}\Omega_1 &= \frac{1}{{2}} \omega_{32} = \frac{1}{{2}} \omega_1 \\ \Omega_2 &= \frac{1}{{2}} \omega_{13} = \frac{1}{{2}} \omega_2 \\ \Omega_3 &= \frac{1}{{2}} \omega_{21} = \frac{1}{{2}} \omega_3 \end{aligned} \hspace{\stretch{1}}(3.6)


\begin{aligned}\omega_{ij}= \epsilon_{ijk} \partial_i u_j\end{aligned} \hspace{\stretch{1}}(5.26)

we can write

\begin{aligned}\Omega_k = -\frac{1}{{2}} \epsilon_{ijk} \partial_i u_j.\end{aligned} \hspace{\stretch{1}}(3.10)

In matrix form this becomes

\begin{aligned}\omega_{ij} = \begin{bmatrix}0 & -\Omega_3 & \Omega_2 \\ \Omega_3 & 0 & -\Omega_1  \\ -\Omega_2 & \Omega_1 & 0\end{bmatrix}.\end{aligned} \hspace{\stretch{1}}(3.11)

For the special case e_{ij} = 0, our displacement equation in vector form becomes

\begin{aligned}d\mathbf{x}' = d\mathbf{x} + \boldsymbol{\Omega} \times d\mathbf{x} \delta t.\end{aligned} \hspace{\stretch{1}}(3.12)

Let’s do a quick verification that this is all kosher.

\begin{aligned}(\boldsymbol{\Omega} \times d\mathbf{x})_i&=\Omega_r dx_s \epsilon_{rsi} \\ &=\left(-\frac{1}{{2}} \epsilon_{abr} \partial_a u_b \right) dx_s \epsilon_{rsi} \\ &=-\frac{1}{{2}} \partial_a u_b dx_s \delta^{[ab]}_{si} \\ &=-\frac{1}{{2}} (\partial_s u_i-\partial_i u_s) dx_s  \\ &=\frac{1}{{2}} \left(\frac{\partial {u_s}}{\partial {x_i}}-\frac{\partial {u_i}}{\partial {x_s}}\right) dx_s  \\ &=\frac{1}{{2}} \left(\frac{\partial {u_j}}{\partial {x_i}}-\frac{\partial {u_i}}{\partial {x_j}}\right) dx_j  \\ &=\omega_{ij} dx_j.\end{aligned}

All’s good in the world of signs and indexes.

The symmetric term (strain tensor).

Now let’s look at the symmetric term. With the initial volume

\begin{aligned}dV = dx_1 dx_2 dx_3,\end{aligned} \hspace{\stretch{1}}(3.13)

and the final volume written assuming that we are working in our principle strain basis, we have (very much like the solids case)

\begin{aligned}dV' &= dx_1' dx_2' dx_3' \\ &= (1 + e_{11} \delta t) dx_1+(1 + e_{22} \delta t) dx_2+(1 + e_{33} \delta t) dx_3\\ &=(1 + (e_{11} + e_{22} + e_{33}) \delta t) dx_1 dx_2 dx_3 + O((\delta t)^2) \\ &=\left(1 + \left(\frac{\partial {u_1}}{\partial {x_1}}+\frac{\partial {u_2}}{\partial {x_2}}+\frac{\partial {u_3}}{\partial {x_3}}\right)\delta t \right) dV \\ &=\left(1 + (\boldsymbol{\nabla} \cdot \mathbf{u}) \delta t\right) dV \\ \end{aligned}

So much like we expressed the relative change of volume in solids, we now can express the relative change of volume per unit time as

\begin{aligned}\frac{dV' - dV}{dV \delta t} = \boldsymbol{\nabla} \cdot \mathbf{u},\end{aligned} \hspace{\stretch{1}}(3.14)


\begin{aligned}\frac{\delta(dV)}{dV \delta t} = \boldsymbol{\nabla} \cdot \mathbf{u},\end{aligned} \hspace{\stretch{1}}(3.15)

We identify the divergence of the displacement as the relative change in volume per unit time.

Newtonian Fluids.

\emph{(Newtonian Fluids)}

A fluid for which the rate of strain tensor is linearly related to stress tensor.

For such a fluid, the constitutive relation takes the form

\begin{aligned}\boxed{\sigma_{ij} = - p \delta_{ij} + 2 \mu e_{ij},}\end{aligned} \hspace{\stretch{1}}(4.16)

where p is called the isotropic pressure, and \mu is the viscosity of the fluid.

For comparison, in solids we had

\begin{aligned}\sigma_{ij} = \lambda e_{kk} \delta_{ij} + 2 \mu e_{ij}\end{aligned} \hspace{\stretch{1}}(4.17)

While we are allowing for rotation in the fluids (\omega_{ij}) that we did not consider for solids, we now impose a requirement that the strain tensor trace is not a function of the fluid displacements, with

\begin{aligned}\lambda e_{kk} = \lambda \boldsymbol{\nabla} \cdot \mathbf{u} = -p.\end{aligned} \hspace{\stretch{1}}(5.26)

What is the physical justification for this? In words this was explained after class as the effect of rotation invariance with an attempt to measure the pressure at a given point in the fluid. It doesn’t matter what direction we place our pressure measurement device at a given fixed location in the fluid. Note that this doesn’t mean the pressure itself is constant. For example with a gravitational body force applied, our pressure will increase with depth in the fluid. Noting this provides a nice physical interpretation of the trace of the strain tensor.

Can we mathematically justify this explanation? We see above that we have

\begin{aligned}\boldsymbol{\nabla} \cdot \mathbf{u} = \frac{\delta \ln(dV)}{\delta t},\end{aligned} \hspace{\stretch{1}}(5.26)

so we are in effect making the identification

\begin{aligned}\ln dV = -p t /\lambda + \ln dV_0\end{aligned} \hspace{\stretch{1}}(5.26)


\begin{aligned}dV = dV_0 e^{-p t/\lambda}.\end{aligned} \hspace{\stretch{1}}(5.26)

The relative change in a differential volume element changes exponentially.


\begin{aligned}[\mu] = \frac{\text{M}}{\text{L}\text{T}}.\end{aligned} \hspace{\stretch{1}}(4.22)

Some examples

\item \mu_{\text{air}} = 1.8 \times 10^{-5} \frac{\text{kg}}{\text{m s}}
\item \mu_{\text{water}} = 1.1 \times 10^{-3} \frac{\text{kg}}{\text{m s}}
\item \mu_{\text{glycerin}} = 2.3 \frac{\text{kg}}{\text{m s}}

Conservation of mass in fluid.

Referring to figure (\ref{fig:continuumL9:continuumL9fig2})
\caption{FIXME: continuumL9fig2}

we have a flow rate

\begin{aligned}\rho \mathbf{u} \delta t ds\end{aligned} \hspace{\stretch{1}}(5.23)


\begin{aligned}\rho \mathbf{u} ds,\end{aligned} \hspace{\stretch{1}}(5.24)

per unit time. Here the velocity of fluid particle is \mathbf{u}.

\begin{aligned}\oint \rho \mathbf{u} \cdot d\mathbf{s},\end{aligned} \hspace{\stretch{1}}(5.25)

we must have

\begin{aligned}\frac{\partial {}}{\partial {t}} \int \rho dV =-\oint \rho \mathbf{u} \cdot d\mathbf{s}.\end{aligned} \hspace{\stretch{1}}(5.26)

\begin{aligned}dm = \rho dV\end{aligned} \hspace{\stretch{1}}(5.27)

\begin{aligned}\frac{dm}{dt} = \frac{d}{dt} (\rho dV)\end{aligned} \hspace{\stretch{1}}(5.28)

positive if fluid is coming in.
negative if fluid is going out.

By Green’s theorem

\begin{aligned}\oint \mathbf{A} \cdot d\mathbf{s} = \int_V (\boldsymbol{\nabla} \cdot \mathbf{u}) dV,\end{aligned} \hspace{\stretch{1}}(5.29)

so we have

\begin{aligned}-\oint \rho \mathbf{u} \cdot d\mathbf{s} = -\int \boldsymbol{\nabla} \cdot (\rho \mathbf{u} ) dV,\end{aligned} \hspace{\stretch{1}}(5.30)

and must have

\begin{aligned}\int \left( \frac{\partial {\rho}}{\partial {t}} + \boldsymbol{\nabla} \cdot (\rho \mathbf{u}) \right) dV = 0.\end{aligned} \hspace{\stretch{1}}(5.31)

The total mass has to be conserved. The mass that is leaving the volume per unit time must move through the surface of the volume in that time. In differential form this is

\begin{aligned}\boxed{\frac{\partial {\rho}}{\partial {t}} + \boldsymbol{\nabla} \cdot (\rho \mathbf{u}) = 0.}\end{aligned} \hspace{\stretch{1}}(5.32)

Operating by chain rule we can write this as

\begin{aligned}\frac{\partial {\rho}}{\partial {t}} + \mathbf{u} \cdot \boldsymbol{\nabla} \rho = - \rho \boldsymbol{\nabla} \cdot \mathbf{u}.\end{aligned} \hspace{\stretch{1}}(5.33)

To make sense of this, observe that we have for f = f(x, y, z, t)

\begin{aligned}\delta f &= \lim_{\delta t \rightarrow 0} \frac{\delta f}{\delta t} dt \\ &=\frac{\partial {f}}{\partial {x}} \frac{\delta x}{\delta t}+\frac{\partial {f}}{\partial {y}} \frac{\delta y}{\delta t}+\frac{\partial {f}}{\partial {z}} \frac{\delta z}{\delta t}+ \frac{\partial {f}}{\partial {t}} \\ &=(\boldsymbol{\nabla} f) \cdot \mathbf{u} + \frac{\partial {f}}{\partial {t}}\end{aligned}

so we have

\begin{aligned}\frac{\partial {\rho}}{\partial {t}} + \mathbf{u} \cdot \boldsymbol{\nabla} \rho = \frac{d\rho}{dt}\end{aligned} \hspace{\stretch{1}}(5.34)


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

Incompressible fluid

When the density doesn’t change note that we have

\begin{aligned}\frac{d\rho}{dt} = 0\end{aligned} \hspace{\stretch{1}}(5.36)

which then implies

\begin{aligned}\boxed{\boldsymbol{\nabla} \cdot \mathbf{u} = 0,}\end{aligned} \hspace{\stretch{1}}(5.37)

at all points in the fluid.


[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: