## PHY454H1S Continuum mechanics. Problem Set 2. Rectilinear flow with pressure gradient. Two layer flow driven by moving surface.

Posted by peeterjoot on March 16, 2012

# Disclaimer.

This problem set is as yet ungraded.

# Problem Q1.

## Statement

Imagine a steady rectilinear blood flow of the form through a two dimensional artery. It is driven by a constant pressure gradient maintained by an external `heart’. The top and bottom walls of the artery are distance apart and the fluid satisfies no-slip boundary conditions at the walls. Assuming that the fluid is Newtonian,

\begin{enumerate}

\item Show that the Navier-Stokes equation reduces to

where is the viscosity of the blood.

\item Show that the velocity profile of the fluid inside the artery is a parabolic profile.

\item What is the maximum speed of the fluid? Draw the velocity profile to show where the maximum speed occurs inside the artery.

\item If due to smoking etc., the viscosity of the blood gets doubled, then what should be the new pressure gradient to be maintained by the `heart’ to keep the liquid flux through the artery at the same level as the non-smoking one?

\end{enumerate}

## Solution. Part 1. Navier-Stokes equation for the system.

The Navier-Stokes equation, for an incompressible unidirectional fluid , assuming that there is no dependence, takes the form

\begin{subequations}

\end{subequations}

With a steady state assumption we kill the term, and 2.2d kills of the x-component of the Laplacian and our non-linear inertial term on the LHS, leaving just

\begin{subequations}

\end{subequations}

With , we have , so 2.3a is reduced to

Finally, since we have an assumption of no z-dependence () and from the incompressibility assumption 2.2d (), we have

as desired.

## Solution. Part 2. Velocity profile.

Integrating 2.5 twice, we have

Application of the no-slip boundary value condition , we have

Adding and subtracting these, we find

\begin{subequations}

\end{subequations}

so the velocity is given by the parabolic function

## Solution. Part 3. Maximum speed of the flow.

It is clear that the maximum speed of the fluid is found at

The velocity profile for this flow is drawn in figure (\ref{fig:continuumProblemSet2:continuumProblemSet2Fig1}).

\begin{figure}[htp]

\centering

\def\svgwidth{0.3\columnwidth}

\caption{Velocity profile for 1D constant pressure gradient steady state flow.}

\end{figure}

## Solution. Part 4. Effects of viscosity doubling.

With our velocity being dependent on the ratio, it is clear, even without calculating the flux, that we will need twice the pressure gradient if the viscosity is doubled to maintain the same flux through the artery and veins. To demonstrate this more thoroughly we can calculate this mass flux. For an element of mass leaving a portion of the conduit, bounded by the plane normal to we have

Integrating this over a width , our flux through the plane is

Doubling the blood viscosity for our smoker, our respective fluxes are

Demanding equality before and after smoking we find

where is the magnitude of the pressure gradient before the bad habits kicked in. The smoker’s poor little heart (soon to be a big overworked and weak heart) has to generate pressure gradients that are twice as big to get the same quantity of blood distributed through the body.

# Problem Q2.

## Statement

Consider steady simple shearing flow with no imposed pressure gradient of a two layer fluid with viscosity

latex -h < y < 0,$} \\ \mu^{(2)} & \quad \mbox{}\end{array}\right.\end{aligned} \hspace{\stretch{1}}(3.15)$

The boundary conditions are no-slip at the lower plate and at . The top plate is moving with a velocity at and fluid is sticking to it. using the continuity of tangential (shear) stress at the interface ()

\begin{itemize}

\item Derive the velocity profile of the two fluids.

\item Calculate the maximum speed.

\item Calculate the mean speed.

\item Calculate the flux (the volume flow rate.)

\item Calculate the tangential force (per unit width) on the strip of the wall .

\item Calculate the tangential force (per unit width) on the strip at the interface by the top fluid on the lower fluid.

\end{itemize}

## Part 1. Derive the velocity profile of the two fluids.

With coordinates referring to figure (\ref{fig:continuumProblemSet2:continuumProblemSet2Fig2}), our steady flow for layers and has the form

\begin{figure}[htp]

\centering

\def\svgwidth{0.4\columnwidth}

\caption{Two layer flow induced by moving wall.}

\end{figure}

\begin{subequations}

\end{subequations}

as we found in Q1. Only the boundary value conditions and the driving pressure are different here. In this problem and the next, we have constant pressure gradients to deal with, so we really have just the pair of equations

to solve. For this Q2 problem we have , so the algebra to match our boundary value constraints becomes a bit easier. Our boundary value constraints are

\begin{subequations}

\end{subequations}

plus one more to match the tangential components of the traction vector with respect to the normal . The components of that traction vector are

but we are only interested in the horizontal component which is

So the matching the tangential components of the traction vector at the interface gives us our last boundary value constraint

and we are ready to do our remaining bits of algebra. We wish to solve the pair of equations

for the four integration constants and using our boundary value constraints. The linear system to solve is

With , we have

Subtracting these to solve for we find

This gives us everything we need

Referring back to 3.20 our velocities are

Checking, we see at a glance we see that we have , , , and as desired.

As an example, let’s add some numbers. With mercury and water in layers and respectively, we have

so that our velocity is

latex y \in [-h, 0]$} \\ -1.37 U \left(1 + 1.7 \frac{y}{h}\right) & \quad \mbox{}\end{array}\right.\end{aligned} \hspace{\stretch{1}}(3.35)$

This is plotted with in figure (\ref{fig:continuumProblemSet2:continuumProblemSet2Fig3})

\begin{figure}[htp]

\centering

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

\caption{Two layer shearing flow with water over mercury.}

\end{figure}

## Part 2. Calculate the maximum speed.

With , and linearly decreasing, then linearly decreasing further from the value at , it is clear that the maximum speed, no matter the viscosities of the fluids, is on the upper moving interface. This maximum takes the value .

## Part 3. Calculate the mean speed.

As linear functions the average speeds of the respective fluids fall on the midpoints . These are

Averaging these two gives us the overall average, so we find

## Part 4. Calculate the flux

We can calculate the volume flux, much like the mass flux (although the mass flux seems a more sensible quantity to calculate). Looking at the rate of change of an element of fluid passing through the plane we have

Integrating over the total height, for a width we have

So our volume flux through a width is

## Part 5. Tangential force on the wall.

From 3.32 the tangential components of our traction vectors are

and

We see that the tangential component of the traction vector is a constant throughout both fluids. Allowing this force to act on a length of the lower wall, our force per unit width over that strip is just

The negative value here makes sense since it is acting to push the fluid backwards in the direction of the upper wall motion.

## Part 6. Tangential force on the lower fluid.

Due to constant nature of the tangential component of the traction vector shown above, the force per unit width of the upper fluid acting on the lower fluid, is also given by 3.41.

# Problem Q3.

## Statement

If on top of the problem described above a constant pressure gradient is applied between the boundaries , describe qualitatively what type of flow profile you would expect in the steady state. Draw the velocity profiles for two cases (i) (ii) . Explain your result.

## Solution

We showed earlier that the Navier-Stokes equations for this Q3 case, where is non-zero were given by 3.17, which restated is

Our solutions will now necessarily be parabolic, of the form

with the tangential traction vector components given by

The boundary value constants become a bit messier to solve for, and should we wish to do so we’d have to solve the system

Without actually solving this system we should expect that our solution will have the form of our pure shear flow, with parabolas superimposed on these linear flows. For a higher viscosity bottom layer , this should look something like figure (\ref{fig:continuumProblemSet2:continuumProblemSet2Fig4}) whereas for the higher viscosity on the top, these would be roughly flipped as in figure (\ref{fig:continuumProblemSet2:continuumProblemSet2Fig5}).

\begin{figure}[htp]

\centering

\includegraphics[totalheight=0.3\textheight]{continuumProblemSet2Fig4}

\caption{Superposition of constant pressure gradient and shear flow solutions ().}

\end{figure}

\begin{figure}[htp]

\centering

\includegraphics[totalheight=0.3\textheight]{continuumProblemSet2Fig5}

\caption{Superposition of constant pressure gradient and shear flow solutions ().}

\end{figure}

This superposition can be justified since we have no term in the Navier-Stokes equations for these systems.

### Exact solutions.

The figures above are kind of rough. It’s not actually hard to solve the system above. After some simplification, we get

Should we wish a more exact plot for any specific values of the viscosities, we could plot exactly with software the vector field described by these velocities.

## Leave a Reply