# Shape of a rotating fluid surface

I was going through fluid dynamics and I came upon this section on the shape of a liquid surface in a rotating cylindrical vessel. The working in the book is something like this:

It considers a mass $m$ of water on the surface of the rotating liquid at a distance $x$ from its axis of rotation .It takes into account a pseudo force ($m \cdot x \cdot \omega^2$) radially outwards and its weight force ($mg$).It then proceeds to calculate the net force and the angle it makes with the vertical.From there through basic integration it gets the equation of the surface of liquid i.e. $y= \dfrac{(x^2 \cdot \omega^2)}{2g}$ .

I get this working.

But if I were to view the same rotation through an inertial frame,how would I proceed? I mean, there wouldn't be any pseudo force, rather a centripetal acceleration of same magnitude radially inwards.This is where I am facing certain doubts.

The forces on the liquid would be the centripetal force inwards and the weight force downwards.The net force would therefore be in a different direction,and consequently the slope of the surface at that point would be in a different direction.How is it possible for the same surface to have two different shapes in two different frames?

I would be obliged if someone could help me out by pointing out the flaw in my understanding.

• Are you sure about the direction of your centripetal force? – Jon Custer Jun 12 '15 at 19:28
• yes, in the inertial frame, it is supposed to be directed inwards from the water surface to the axis of rotation,right? – Daipayan Mukherjee Jun 13 '15 at 8:24
• If you use an inertial frame, your fluid is no longer at rest. You can solve it, but not with fluidstatics formulae. Rather, you have tos olve the entire equations of motion (conservation of momentum, Navier-Stokes). This is obviously harder. That's why we choose a rotating frame! – FGSUZ Feb 10 '19 at 21:13

Imagine a little piece of water on the top surface at radius $x$. The height of the water is $y(x)$.

By moving inward $\mathrm{d}x$, the water would reduce its energy by $m g \mathrm{d}y = m g y' \mathrm{d}x$. It therefore feels a force of magnitude $m g y'$ towards the center.

This force causes the water to accelerate towards the center. We know this to be the centripetal acceleration, so

$$m \omega^2 x = m g y'$$

from which we recover your original expression.

"The forces on the liquid would be the centripetal force inwards and the weight force downwards."

The forces on a bit of water are gravity downwards and a force from a pressure gradient. Without knowing the exact pressure at every point, we can still calculate the net force from energy, as I did above.

"THe net force would therefore be in a different direction"

The direction of the net force on any bit of water is horizontally towards the center. If you wanted to calculate all the force exactly, you would need the pressure everywhere in the fluid. Note that the pressure is not the usual $\rho g h$ that you may have learned in an introductory class.

Although Mark Eichenlaub's answer is perfectly fine, I will post this answer, so that the time I spent writing it won't go to waste, and someone might find its detailedness useful.

To understand what happens when you change back and forth between inertial and rotating reference frames, consider a conical pendulum. In an inertial frame you see that the object at the end of the string is rotating around with uniform circular motion. The condition for that is that the sum of all forces acting upon it is pointing inwards, and this sum is equal to $m*\vec{\mathbf{a}}_{cp}=\left(m*v^2/R\right)\vec{\mathbf{e}}_{in}=m\omega^2*\vec{\mathbf{s}}$ where $v$ is the tangential velocity and $R$ is the radius of the trajectory, and $\vec{\mathbf{e}}_{in}$ is a unit vector pointing inwards, $\vec{s}=R*\vec{\mathbf{e}}_{in}$. As a conclusion you can write the forces acting on the object as $$\vec{\mathbf{K}}+m*\vec{\mathbf{g}}=\vec{\mathbf{F}}_{cp}.$$

Now switch to a rotating reference frame rotating with the same $\omega$ angular speed as the object. In this frame all forces acting on the object are the same, that is still $\vec{\mathbf{K}}+m*\vec{\mathbf{g}}$, but the object seems to be in a static equilibrium, thus a force equal to the oppostie of the sum of the former two must act on it. This is the centrifugal force, and its obvious that it must equal in magnitude and be opposite in sign with the centripetal force that is the sum of the forces coming from the string and gravitation. So now the forces acting on the object are $$\vec{\mathbf{K}}+m*\vec{\mathbf{g}} + \vec{\mathbf{F}}_{cf} = \vec{\mathbf{0}}.$$

From this you can see that simply taking the centrifugal force to the other side and reversing its sign gives you the centripetal force.

For a liquid rotating, first take the equation of continuums: $$\rho \frac{\text{d}\vec{\mathbf{v}}}{\text{d}t} = \vec{\mathbf{f}}+\text{Div}\hat{\sigma},$$ where $\hat{\sigma}$ is the stress tensor of liquids, and in a frictionless case $\text{Div}\hat{\sigma}=-\vec{\bigtriangledown} p$. $\vec{\mathbf{f}}$ is the density of outer forces. In case of a rotating frame the liquid is in equilibrium, thus this equation reads as $\vec{\mathbf{0}}=\vec{\mathbf{f}}+\text{Div}\hat{\sigma}$, where $\vec{\mathbf{f}}$ can be calculated as follows: $$\vec{\mathbf{f}}_{grav} = \frac{\text{d}\vec{\mathbf{F}}_{grav}}{\text{d}V}=\frac{1}{\text{d}V}\text{d}m*\vec{\mathbf{g}}=\frac{1}{\text{d}V}\text{d}m\left[-\vec{\bigtriangledown}(gz)\right]=-\vec{\bigtriangledown}\left(\rho gz\right)$$ $$\vec{\mathbf{f}}_{cf} = \frac{\text{d}\vec{\mathbf{F}}_{cf}}{\text{d}V}=-\frac{1}{\text{d}V}\text{d}m\omega^2\vec{\mathbf{s}}=-\frac{1}{\text{d}V}\text{d}m\left[-\vec{\bigtriangledown}(\frac{1}{2}\omega^2(x^2+y^2))\right]=+\vec{\bigtriangledown}\left(\frac{1}{2}\rho\omega^2s^2\right)$$ $$\vec{\mathbf{f}}=\vec{\mathbf{f}}_{grav}+\vec{\mathbf{f}}_{cf}$$ And now simply the equilibrium condition $\vec{\mathbf{f}}+\text{Div}\hat{\sigma}=\vec{\mathbf{0}}$ tells us that $$\vec{\bigtriangledown}\left(-p-\rho gz+\rho\omega^2s^2\right)=\vec{\mathbf{0}}$$ $$-p-\rho gz+\rho\omega^2s^2=-p-\rho gz+\rho\omega^2(x^2+y^2)=\text{const},$$ which gives the parabolic profile (*).

If you want to use an inertial reference frame, then just simply use instead of the centrifugal force that the sum of the two forces on a piece of material with mass $\text{d}m$ and volume $\text{d}V$ must be the centripetal force to uphold circular motion. In this case the l.h.s. of $\rho \frac{\text{d}\vec{\mathbf{v}}}{\text{d}t} = \vec{\mathbf{f}}+\text{Div}\hat{\sigma}$ won't be zero, but the centripetal force density. Using that $\text{d}\vec{\mathbf{v}}/\text{d}t$ must be the centripetal acceleration $\omega^2\vec{\mathbf{s}}$, you get: $$\rho\omega^2\vec{\mathbf{s}}=-\vec{\bigtriangledown}p - \vec{\bigtriangledown} gz$$ $$-\vec{\bigtriangledown}\rho\omega^2s^2=-\vec{\bigtriangledown}p - \vec{\bigtriangledown} gz,$$ which is the same equation as before giving the parabolic profile.

Sidenote: to see that the l.h.s is indeed the centripetal force density: $$\vec{\mathbf{f}}_{cp}=\frac{1}{\text{d}V}\text{d}\vec{\mathbf{F}}_{cp}=\frac{1}{\text{d}V}\text{d}m\omega^2\vec{\mathbf{s}}=\rho\omega^2\vec{\mathbf{s}}$$

* How to obtain the profile:

At $z=h_{\text{liquid}}=h_{\text{liquid}}(x,y)$ pressure is $p_{\text{air}}$, so the const is simply $const=p_{\text{air}}$, and taking the restriction of the l.h.s. to be on the surface of the liquid we get $$-p\left(x,y,z=h_{\text{liquid}}(x,y)\right)-\rho gh_{\text{liquid}}(x,y)+ \rho\omega^2s^2=p_{\text{air}}$$ $$-p_{\text{air}}-\rho gh_{\text{liquid}}+\rho\omega^2s^2=p_{\text{air}}$$ $$\rho gh_{\text{liquid}}=\rho\omega^2s^2$$ $$h(x,y)_{\text{liquid}}=\frac{\omega^2}{g}(x^2+y^2)$$

The centripetal force is not a force that is required to exist, but rather a description of a force that would keep a particle in uniform circular motion. The only forces on the particle here are its weight and forces from other particles (modeled by the fluid density and the incompressibility assumption).

Given incompressibility assumption and steady flow;

Consider Euler's equation in a uniform gravitational field:

$$\dfrac{\partial{\vec{v}}}{\partial{t}} + (\vec{v}*{\vec\nabla})\vec{v} = -\vec\nabla{\dfrac{p}{\rho}}+ \vec{g}$$

Particle velocity field is solenoidal, $$\vec \nabla \cdot \vec v$$ = 0, $$x$$ and $$y$$ velocities can be expressed in terms of the stream function $$\Psi$$.

$$v_{x} = \dfrac{\partial{\Psi}}{\partial{y}}$$ = $$-y\omega$$

$$v_{y} = -\dfrac{\partial{\Psi}}{\partial{x}}$$ = $$x\omega$$

$$v_{z} = 0$$ (I am assuming the z-axis to be the axis of rotation)

Upon substitution of each component individually into Euler's Equation, we arrive at the centripetal acceleration in terms of surface pressure (inertial frame, as was requested);

$$x\omega^{2} = \dfrac{1}{\rho}\dfrac{\partial{p}}{\partial{x}}$$

$$y\omega^{2} = \dfrac{1}{\rho}\dfrac{\partial{p}}{\partial{y}}$$

$$\dfrac{1}{\rho}\dfrac{\partial{p}}{\partial{z}}+g = 0$$

Now, evaluate the general integral of these expressions;

$$\dfrac{p}{\rho} = \dfrac{\omega^{2}(x^2+y^2)}{2} - gz + C$$

Given $$p=constant$$ at a free surface and solving for z,

$$z = \dfrac{\omega^{2}(x^2+y^2)}{2g}$$

The surface adopts the shape of a paraboloid. There shouldn't be any confusion; the centrifugal force and centripetal forces are merely different tools used to approach the same problem. The centripetal force points in the horizontal direction, NOT upwards as you have stated in the problem description.