## PHY454H1S Continuum Mechanics. Lecture 11: Worked examples of Navier-Stokes solutions. Taught by Prof. K. Das.

Posted by peeterjoot on February 16, 2012

# Disclaimer.

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

# Navier-Stokes equation.

The Navier-Stokes equation (our fluids equivalent to Newton’s second law) was found to be

In this course we’ll focus on the incompressible case where we have

We watched a video of the rocking tank as in figure (\ref{fig:continuumL11:continuumL11fig1}). The boundary condition that accounted for the matching of the die marker is that we have *no slipping* at the interface. Writing for the tangent to the interface then this condition at the interface is described mathematically by the conditions

\begin{figure}[htp]

\centering

\includegraphics[totalheight=0.2\textheight]{continuumL11fig1}

\caption{Rocking tank velocity matching.}

\end{figure}

Referring to figure (\ref{fig:continuumL11:continuumL11fig2}) where the tangents and normals are depicted an example representation of the normal and tangent vectors for the fluids are

\begin{figure}[htp]

\centering

\includegraphics[totalheight=0.2\textheight]{continuumL11fig2}

\caption{Normals and tangents at interface for 2D system}

\end{figure}

For the traction vector

we also have at the interface we must have matching of

More explicitly, in coordinates this is

# Steady incompressible rectilinear (unidirectional) flow.

In this case we can fix our axis so that

where the velocity components in the other directions

are both zero. Symbolically, the steady state condition is

We start with the incompressibility condition, which written explicitly, is

or

This implies

so our velocity can only be function of the and coordinates only

The non-linear term of the Navier-Stokes equation takes the form

With incompressibility and conditions killing this term, and the steady state condition 3.13 killing the term, the Navier-Stokes equation for this incompressible unidirectional steady state flow (in the absence of body forces) is reduced to

In coordinates this is

\begin{subequations}

\end{subequations}

Operating on the first with an x partial we find

Since we have

we also have

so our pressure must be linear with position

as illustrated in figure (\ref{fig:continuumL11:continuumL11fig3})

\begin{figure}[htp]

\centering

\includegraphics[totalheight=0.2\textheight]{continuumL11fig3}

\caption{Pressure gradient in 1D system.}

\end{figure}

latex x = 0$} \\ p_L & \quad \mbox{} \end{array}\right.\end{aligned} \hspace{\stretch{1}}(3.24)$

we have

and

## Example: Shearing flow.

The flows of this sort don’t have to be trivial. For example, even with constant pressure () as in figure (\ref{fig:continuumL11:continuumL11fig4}) we can have a “shearing flow” where the fluids at the top surface are not necessarily moving at the same rates as the fluid below that surface. We have fluid flow in the direction only, and our velocity is a function only of the coordinate.

\begin{figure}[htp]

\centering

\includegraphics[totalheight=0.2\textheight]{continuumL11fig4}

\caption{Velocity variation with height in shearing flow.}

\end{figure}

For such a flow 3.19a simplifies to

with solution

## Example: Channel flow

This time our simplified Navier-Stokes equation 3.19a is reduced to something slightly more complicated

with solution

The boundary value conditions with the coordinate system in use illustrated in figure (\ref{fig:continuumL11:continuumL11fig5}) require the velocity to be zero at the interface (the pipe walls preventing flow in the interior of the pipe)

\begin{figure}[htp]

\centering

\includegraphics[totalheight=0.2\textheight]{continuumL11fig5}

\caption{1D Channel flow coordinate system setup.}

\end{figure}

One solution, immediately evident is,

so our solution becomes

a parabolic velocity flow. This is illustrated graphically in figure (\ref{fig:continuumL11:continuumL11fig6}).

\begin{figure}[htp]

\centering

\includegraphics[totalheight=0.2\textheight]{continuumL11fig6}

\caption{Parabolic velocity distribution.}

\end{figure}

It’s clear that this is maximized by , but we can also see this by computing

This maximum is

The flux, or flow rate is

Let’s now compute the strain () and the stress ()

stress

This can be used to compute the forces on the inner surfaces of the tube. As illustrated in figure (\ref{fig:continuumL11:continuumL11fig7}), our normals at are respectively. The traction vector in the direction is at is

\begin{figure}[htp]

\centering

\includegraphics[totalheight=0.2\textheight]{continuumL11fig7}

\caption{Normals in 1D channel flow system}

\end{figure}

so that

(here the directionality comes from the index of the stress tensor).

The total force is then

FIXME: ask in class. This is the tangential force at the boundary of the wall. What is it a force on? If it is tangential, how can it act on the wall? It could act on an impediment placed right up next to the wall, but if that’s the case, why are we integrating from to ?

## Leave a Reply