## PHY454H1S Continuum mechanics. Problem Set 2. Two layer flow driven by moving surface. (revisited)

Posted by peeterjoot on March 26, 2012

# Problem Q3 (revisited).

I’d produced the following sketches. 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}).

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, I find with Mathematica the following solution

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.

I suppose it is cheating to use Mathematica and then say that the solution is easy? To make amends for being lazy with my algebra, let’s show that it is easy to do manually too. I’ll do the same problem manually, but generalize it slightly. We can do this easily if we just be a bit marter with our integration constants. Let’s solve the problem for the upper and lower walls moving with velocity and respectively, and let the heights from the interface be and respectively.

We have the same set of differential equations to solve, but now let’s write our solution with the undetermined coefficients expressed as

Now it’s super easy to match the boundary conditions at and (the lower and upper walls respectively). Clearly the integration constants are just the velocities. Matching the tangential component of the traction vectors at we have

and matching velocities at gives us

This gives us

Plotting this with sliders or animation in Mathematica is a fun way to explore visualizing this. The results vary widely depending on the various parameters. Such a Mathematica notebook can be found in https://raw.github.com/peeterjoot/physicsplay/master/notes/phy454/problemSetIIQ3PlotWithManipulate.cdf. A pair of cached animations with variation of the pressure gradient for , , showing the superposition of the shear and channel flow solutions can be found here

- With fluids with densities of mercury and water in the lower (1) and upper (2) layers respectively.
- With , this is plotted fluids with densities of water and mercury in the lower (1) and upper (2) layers respectively.

## Leave a Reply