## Continuum mechanics fluids review.

Posted by peeterjoot on April 25, 2012

# Motivation.

Review of key ideas and equations from the fluid dynamics portion of the class.

# Vector displacements.

Those portions of the theory of elasticity that we did cover have the appearance of providing some logical context for the derivation of the Navier-Stokes equation. Our starting point is almost identical, but we now look at displacements that vary with time, forming

We compute a first order Taylor expansion of this differential, defining a symmetric strain and antisymmetric vorticity tensor

\begin{subequations}

\end{subequations}

Allowing us to write

We introduced vector and dual vector forms of the vorticity tensor with

\begin{subequations}

\end{subequations}

or

\begin{subequations}

\end{subequations}

We were then able to put our displacement differential into a partial vector form

# Relative change in volume

We are able to identity the divergence of the displacement as the relative change in volume per unit time in terms of the strain tensor trace (in the basis for which the strain is diagonal at a given point)

# Conservation of mass

Utilizing Green’s theorem we argued that

We were able to relate this to rate of change of density, computing

An important consequence of this is that for incompressible fluids (the only types of fluids considered in this course) the divergence of the displacement .

# Constitutive relation

We consider only Newtonian fluids, for which the stress is linearly related to the strain. We will model fluids as disjoint sets of hydrostatic materials for which the constitutive relation was previously found to be

# Conservation of momentum (Navier-Stokes)

As in elasticity, our momentum conservation equation had the form

where are the components of the external (body) forces per unit volume acting on the fluid.

Utilizing the constitutive relation and explicitly evaluating the stress tensor divergence we find

Since we treat only incompressible fluids in this course we can decompose this into a pair of equations

\begin{subequations}

\end{subequations}

# No slip condition

We’ll find in general that we have to solve for our boundary value conditions. One of the important constraints that we have to do so will be a requirement (experimentally motivated) that our velocities match at an interface. This was illustrated with a rocker tank video in class.

This is the no-slip condition, and includes a requirement that the fluid velocity at the boundary of a non-moving surface is zero, and that the fluid velocity on the boundary of a moving surface matches the rate of the surface itself.

For fluids and separated at an interface with unit normal and unit tangent we wrote the no-slip condition as

\begin{subequations}

\end{subequations}

For the problems we attempt, it will often be enough to consider only the tangential component of the velocity.

# Traction vector matching at an interface.

As well as matching velocities, we have a force balance requirement at any interface. This will be expressed in terms of the traction vector

where is the normal pointing from the interface into the fluid (so the traction vector represents the force of the interface on the fluid). When that interface is another fluid, we are able to calculate the force of one fluid on the other.

In addition the the constraints provided by the no-slip condition, we’ll often have to constrain our solutions according to the equality of the tangential components of the traction vector

We’ll sometimes also have to consider, especially when solving for the pressure, the force balance for the normal component of the traction vector at the interface too

As well as having a messy non-linear PDE to start with, our boundary value constraints can be very complicated, making the subject rich and tricky.

## Flux

A number of problems we did asked for the flux rate. A slightly more sensible physical quantity is the mass flux, which adds the density into the mix

# Worked problems from class.

A number of problems were tackled in class

- channel flow with external pressure gradient
- shear flow
- pipe (Poiseuille) flow
- steady state gravity driven film flow down a slope.

## Channel and shear flow

An example that shows many of the features of the above problems is rectilinear flow problem with a pressure gradient and shearing surface. As a review let’s consider fluid flowing between surfaces at , the lower surface moving at velocity and pressure gradient we find that Navier-Stokes for an assumed flow of takes the form

We find that this reduces to

with solution

Application of the no-slip velocity matching constraint gives us in short order

With this is the channel flow solution, and with this is the shearing flow solution.

Having solved for the velocity at any height, we can also solve for the mass or volume flux through a slice of the channel. For the mass flux per unit time (given volume flux )

we find

We can also calculate the force of the boundaries on the fluid. For example, the force per unit volume of the boundary at on the fluid is found by calculating the tangential component of the traction vector taken with normal . That tangent vector is found to be

The tangential component is the component evaluated at , so for the lower and upper interfaces we have

so the force per unit area that the boundary applies to the fluid is

Does the sign of the velocity term make sense? Let’s consider the case where we have a zero pressure gradient and look at the lower interface. This is the force of the interface on the fluid, so the force of the fluid on the interface would have the opposite sign

This does seem reasonable. Our fluid flowing along with a positive velocity is imparting a force on what it is flowing over in the same direction.

# Hydrostatics.

We covered hydrostatics as a separate topic, where it was argued that the pressure in a fluid, given atmospheric pressure and height from the surface was

As noted below in the surface tension problem, this is also a consequence of Navier-Stokes for (following from ).

We noted that replacing the a mass of water with something of equal density would not change the non-dynamics of the situation. We then went on to define Buoyancy force, the difference in weight of the equivalent volume of fluid and the weight of the object.

# Mass conservation through apertures.

It was noted that mass conservation provides a relationship between the flow rates through apertures in a closed pipe, since we must have

and therefore for incompressible fluids

So if we must have .

## Curve for tap discharge.

We can use this to get a rough idea what the curve for water coming out a tap would be. Suppose we measure the volume flux, putting a measuring cup under the tap, and timing how long it takes to fill up. We then measure the radii at different points. This can be done from a photo as in figure (1).

After making the measurement, we can get an idea of the velocity between two points given a velocity estimate at a point higher in the discharge. For a plain old falling mass, our final velocity at a point measured from where the velocity was originally measured can be found from Newton’s law

Solving for , we find

Mass conservation gives us

or

For the image above I measured a flow rate of about 250 ml in 10 seconds. With that, plus the measured radii at 0 and , I calculated that the average fluid velocity was , vs a free fall rate increase of . Not the best match in the world, but that’s to be expected since the velocity has been considered uniform throughout the stream profile, which would not actually be the case. A proper treatment would also have to treat viscosity and surface tension.

In figure (2) is a plot of the measured radial distance compared to what was computed with 11.41. The blue line is the measured width of the stream as measured, the red is a polynomial curve fitted to the raw data, and the green is the computed curve above.

# Bernoulli equation.

With the body force specified in gradient for

and utilizing the vector identity

we are able to show that the steady state, irrotational, non-viscous Navier-Stokes equation takes the form

or

This is the Bernoulli equation, and the constants introduce the concept of streamline.

FIXME: I think this could probably be used to get a better idea what the tap stream radius is, than the method used above. Consider the streamline along the outermost surface. That way you don’t have to assume that the flow is at the average velocity rate uniformly throughout the stream. Try this later.

# Surface tension.

## Surfaces, normals and tangents.

We reviewed basic surface theory, noting that we can parameterize a surface as in the following example

Computing the gradient we find

Recalling that the gradient is normal to the surface we can compute the unit normal and unit tangent vectors

\begin{subequations}

\end{subequations}

## Laplace pressure.

We covered some aspects of this topic in class. [1] covers this topic in typical fairly hard to comprehend) detail, but there’s lots of valuable info there. section 2.4.9-2.4.10 of [2] has small section that’s a bit easier to understand, with less detail. Recommended in that text is the “Surface Tension in Fluid Mechanics” movie which can be found on youtube in three parts \youtubehref{DkEhPltiqmo}, \youtubehref{yiixltf\_HKw}, \youtubehref{5d6efCcwkWs}, which is very interesting and entertaining to watch.

It was argued in class that the traction vector differences at the surfaces between a pair of fluids have the form

where is the tangential (interfacial) gradient, is the surface tension, a force per unit length value, and is the radius of curvature.

In static equilibrium where (since if ), then dotting with we must then have

## Surface tension gradients.

Considering the tangential component of the traction vector difference we find

If the fluid is static (for example, has none of the creep that we see in the film) then we must have . It’s these gradients that are responsible for capillary flow and other related surface tension driven motion (lots of great examples of that in the film).

## Surface tension for a spherical bubble.

In the film above it is pointed out that the surface tension equation we were shown in class

is only for spherical objects that have a single radius of curvature. This formula can in fact be derived with a simple physical argument, stating that the force generated by the surface tension along the equator of a bubble (as in figure (3)), in a fluid would be balanced by the difference in pressure times the area of that equatorial cross section. That is

Observe that we obtain 13.52 after dividing through by the area.

## A sample problem. The meniscus curve.

To get a better feeling for this, let’s look to a worked problem. The most obvious one to try to attempt is the shape of a meniscus of water against a wall. This problem is worked in [1], but it is worth some extra notes. As in the text we’ll work with axis up, and the fluid up against a wall at as illustrated in figure (4).

The starting point is a variation of what we have in class

where is the atmospheric pressure, is the fluid pressure, and the (signed!) radius of curvatures positive if pointing into medium 1 (the fluid).

For fluid at rest, Navier-Stokes takes the form

With we have

or

We have , the atmospheric pressure, so our pressure difference is

We have then

One of our axis of curvature directions is directly along the axis so that curvature is zero . We can fix the constant by noting that at , , we have no curvature . This gives

That leaves just the second curvature to determine. For a curve our absolute curvature, according to [3] is

Now we have to fix the sign. I didn’t recall any sort of notion of a signed radius of curvature, but there’s a blurb about it on the curvature article above, including a nice illustration of signed radius of curvatures can be found in this wikipedia radius of curvature figure for a Lemniscate. Following that definition for a curve such as we’d have a positive curvature, but the text explicitly points out that the curvatures are will be set positive if pointing into the medium. For us to point the normal into the medium as in the figure, we have to invert the sign, so our equation to solve for is given by

The text introduces the capillary constant

Using that capillary constant to tidy up a bit and multiplying by a integrating factor we have

we can integrate to find

Again for we have , , so . Rearranging we have

Integrating this with Mathematica I get

It looks like the constant would have to be fixed numerically. We require at

but we don’t have an explicit function for .

# Non-dimensionality and scaling.

With the variable transformations

we can put Navier-Stokes in dimensionless form

Here is Reynold’s number

A relatively high or low Reynold’s number will effect whether viscous or inertial effects dominate

The importance of examining where one of these effects can dominate was clear in the Blassius problem, where doing so allowed for an analytic solution that would not have been possible otherwise.

# Eulerian and Lagrangian.

We defined

- Lagrangian: the observer is moving with the fluid.
- Eulerian: the observer is fixed in space, watching the fluid.

# Boundary layers.

## Impulsive flow.

We looked at the time dependent unidirectional flow where

latex t < 0$} \\ U & \quad \mbox{for }\end{array}\right.\end{aligned} \hspace{\stretch{1}}(16.76)$

and utilized a similarity variable with

and were able to show that

The aim of this appears to be as an illustration that the boundary layer thickness grows with .

FIXME: really need to plot 16.80.

## Oscillatory flow.

Another worked problem in the boundary layer topic was the Stokes boundary layer problem with a driving interface of the form

with an assumed solution of the form

we found

\begin{subequations}

\end{subequations}

This was a bit more obvious as a boundary layer illustration since we see the exponential drop off with every distance multiple of .

## Blassius problem (boundary layer thickness in flow over plate).

We examined the scaling off all the terms in the Navier-Stokes equations given a velocity scale , vertical length scale and horizontal length scale . This, and the application of Bernoulli’s theorem allowed us to make construct an approximation for Navier-Stokes in the boundary layer

\begin{subequations}

\end{subequations}

With boundary conditions

With a similarity variable

and stream functions

and

we were able to show that our velocity dependence was given by the solutions of

This was done much more clearly in [4] and I worked this problem myself with a hybrid approach (non-dimensionalising as done in class).

FIXME: the end result of this is a plot (a nice one can be found in [5]). That plot ends up being one that’s done in terms of the similarity variable . It’s not clear to me how this translates into an actual velocity profile. Should plot these out myself to get a feel for things.

# Singular perturbation theory.

The non-dimensional form of Navier-Stokes had the form

where the inverse of Reynold’s number

can potentially get very small. That introduces an ill-conditioning into the problems that can make life more interesting.

We looked at a couple of simple LDE systems that had this sort of ill conditioning. One of them was

for which the exact solution was found to be

The rough idea is that we can look in the domain where and far from that. In this example, with far from the origin we have roughly

so we have an asymptotic solution close to . Closer to the origin where we can introduce a rescaling to find

This gives us

for which we find

# Stability.

We characterized stability in terms of displacements writing

and defining

- Oscillatory unstability. .
- Marginal unstability. .
- Neutral stability. .

## Thermal stability: Rayleigh-Benard problem.

We considered the Rayleigh-Benard problem, looking at thermal effects in a cavity. Assuming perturbations of the form

and introducing an equation for the base state

we found

Operating on this with we find

from which we apply back to 18.103 and take just the z component to find

With an assumption that density change and temperature are linearly related

and operating with the Laplacian we end up with a relation that follows from the momentum balance equation

We also applied our perturbation to the energy balance equation

We determined that the base state temperature obeyed

with solution

This and application of the perturbation gave us

We used this to non-dimensionalize with

latex d$} \\ t & \quad \mbox{with } \\ \delta w & \quad \mbox{with } \\ \delta T & \quad \mbox{with }\end{array}\end{aligned} \hspace{\stretch{1}}(18.112)$

And found (primes dropped)

\begin{subequations}

\end{subequations}

where we’ve introduced the Rayleigh number and Prandtl number’s

\begin{subequations}

\end{subequations}

We were able to construct some approximate solutions for a problem similar to these equations using an assumed solution form

Using these we are able to show that our PDEs are similar to that of

where . Using the trig solutions that fall out of this we were able to find the constraint

which for , this gives us the critical value for the Rayleigh number

which is the boundary for thermal stability or instability.

The end result was a lot of manipulation for which we didn’t do any sort of applied problems. It looks like a theory that requires a lot of study to do anything useful with, so my expectation is that it won’t be covered in detail on the exam. Having some problems to know why we spent two days on it in class would have been nice.

# References

[1] L.D. Landau and E.M. Lifshitz. *A Course in Theoretical Physics-Fluid Mechanics*. Pergamon Press Ltd., 1987.

[2] S. Granger. *Fluid Mechanics*. Dover, New York, 1995.

[3] Wikipedia. Curvature — wikipedia, the free encyclopedia [online]. 2012. [Online; accessed 25-April-2012]. http://en.wikipedia.org/w/index.php?title=Curvature&oldid=488021394.

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

[5] Wikipedia. Blasius boundary layer — wikipedia, the free encyclopedia [online]. 2012. [Online; accessed 28-March-2012]. http://en.wikipedia.org/w/index.php?title=Blasius_boundary_layer&oldid=480776115.

## Leave a Reply