3 9/17

1. For the Poisseuille flow of the previous question, derive the principal strain rates and the principal strain directions.

2. If you put a cup of coffee at the center of a rotating turn table and wait, eventually, the coffee will be executing a “solid body rotation” in which the velocity field is, in cylindrical coordinates:
   
   $$\vec v = {\widehat \imath}_\theta \Omega r$$
   
   
   where $\Omega$ is the angular velocity of the turn table. Draw samples of the streamlines of this flow. Find the vorticity and the strain rate tensor for this flow, using the expressions in appendices B and C. Show that indeed the coffee moves as a solid body; i.e. the fluid particles do not deform, and that for a solid body motion like this, indeed the vorticity is twice the angular velocity.

3. An “ideal vortex flow” is described in cylindrical coordinates by
   
   $$\vec v = \frac{C}{2\pi r} {\widehat \imath}_\theta$$
   
   
   where $C$ is some constant. Draw samples of the streamlines of this flow. In cylindrical coordinates, nabla is given by
   
   $$\nabla = {\widehat \imath}_r \frac{\partial}{\partial r} +... ...ial\theta} + {\widehat \imath}_z \frac{\partial}{\partial z}$$
   
   
   Evaluate
   
   $$\left\vert \begin{array}{ccc} {\widehat \imath}_r & {\wid... ...partial z} \\ v_r & v_\theta & v_z \end{array} \right\vert$$
   
   
   for this flow. Compare the answer with the vorticity $\omega=\nabla\times\vec v$, for which you can find the correct expressions in appendix B. Explain why the determinant does not give the correct result for the vorticity.

4. The “circulation $\Gamma$“ along a closed contour is defined as
   
   $$\Gamma \equiv \oint \vec v \cdot {\rm d}\vec r$$
   
   
   For both the solid body rotation of question 2, and the vortex flow of question 3, find the circulation along the unit circle in the $x,y$-plane. Next, only for the vortex flow, find the circulation along the closed curve consisting of the following segments:
   1. The part of the curve $y=\cosh(x)$, $z=0$ from $x=2$ to $x=-3$;
   2. vertically downward to the curve $y=-\cosh(\cosh(x))$, $z=0$ at $x=-3$;
   3. following the curve $y=-\cosh(\cosh(x))$, $z=0$ from $x=-3$ to the point $x=0$, hence $y=-\cosh(1)$, $z=0$;
   4. in a straight line along the $z$ direction to the point $x=0$, $y=-\cosh(1)$, $z=3.5$;
   5. in a straight line from $x=0$, $y=-\cosh(1)$, $z=3.5$ to the starting point $x=2$, $y=\cosh(2)$, $z=0$.
   Note that in cylindrical coordinates
   
   $${\rm d}\vec r = {\widehat \imath}_r {\rm d}r + {\widehat \imath}_\theta r{\rm d}\theta + {\widehat \imath}_z{\rm d}z$$
   
   

5. According to the Stokes theorem of Calculus III, you should have
   
   $$\oint \vec v \cdot{\rm d}\vec r = \int\nabla\times\vec v \cdot\vec n{ \rm d}S$$
   
   
   where the second integral is over the inside of the contour. So instead of integrating the circulation $\Gamma$ as you did in question 4, you could have integrated the component of vorticity normal to the circle over the inside of the circle. Show that if you do that integral using the vorticity that you found for solid body rotation in question 2, you do indeed get the same answer as you got in question 4. Fine. But now show that if you do the integral of the vorticity over the inside of the circle for the vortex flow of question 3, you do not get the same answer for the circulation as in question 4. Explain which value is correct. And why the other value is wrong.

6. Following Friday‘s lecture, you would of course love to do some integrals of the form $\int\vec{v}\cdot\vec{n}{ \rm d}{S}$ and $\int\vec{v}\cdot\vec{n}\vec{v}{ \rm d}{S}$. Here is your chance. Do question 5.1(b) and 5.1(d) and explain their physical meaning. Take the surfaces $S_{I}$, $S_{II}$, $S_{III}$, and $S_{IV}$ to be one unit length in the $z$-direction. (To figure out the correct direction of the normal vector $\vec n$ at a given surface point, note that the control volume in this case is the right half of the region in between two cylinders of radii $r_0$ and $R_0$ and of unit length in the $z$-direction. The vector $\vec n$ is a unit normal vector sticking out of this control volume.)