Steady state velocity profile of stirred cup of non-bottomless coffee (continued).
Posted by peeterjoot on April 30, 2012
Solving for the flow below the stirring point.
Previously, we found the functional form for an azimuthal flow that has -axis dependence. Attempting apply boundary conditions that including a cylindrical mixing device got us into trouble, since our no-slip condition would simultaneously require zero and non-zero velocity at the point of contact of the mixing interface and the bottom of the vessel (i.e. where the stir stick contacts the bottom of the cup). We can avoid this issue by constraining the mixing to only occur above the bottom of the cup, and then look at the flow that this induces below the mixing point.
Let’s attempt to solve the boundary value problem. We’ll position the mixing cylinder at a height above the bottom of the cup, and set the radius of that inner cylinder to as before, with the mixing occurring at an angular velocity of . We now want apply these boundary value constrains to the velocity function we’ve found for the steady state problem. Provided , this will have the form
where are the zeros of the order one Bessel function . Observe that our boundary value conditions of and are automatically satisfied. Note that because of the equality here, this form of solution is no good if we are mixing at the edge of the cup, so we require .
Our remaining boundary value condition , means that we have to solve
Having had the same sort of problem in our steady state bottomless coffee problem, we know how to solve this
So we find that below the point of stirring (), our steady state solution for the velocity should be
For a cup size of , stir radius of , stir depth , and angular velocity , we find
We can plot this velocity field as in figure (1), but it’s hard to see the radial dependence.
The radial dependence of the magnitude of the velocity can be seen in figure (2), which plots . We see the zero velocity at the edges of the cup as expected, and once we hit the stir height, matching of stir velocity. An animation showing the variation of the radial velocity profile at various depths up to the stir height is available at http://youtu.be/BS8XQdXljSk, and also below.
Observing this animation we see that the velocity is dominated by the first term in the Bessel series, essentially just scaled by the hyperbolic sine that multiplies it.