Column casting with Hagen Poiseuille
Storyboard
It is considered a column with water with a hole in its lower part. The emptying is monitored obtaining an output speed depending on the height of the column.
If the data is modeled with Bernoulli but the exit through the hole is modeled with Hagen Poiseville, which corrects the problem of the case in which it was assumed without viscosity.
ID:(1428, 0)
Viscosity measurement
Description
If a small sphere with radius $a$ is dropped into a medium with viscosity $\eta$, it will accelerate until the gravitational force,
$mg=\displaystyle\frac{4\pi}{3}a^3\rho_sg$
where $\rho_s$ is the density of the sphere material, is equal to the viscous force,
$6\pi \eta a v$
where $v$ is the velocity.
Hence, it is possible to estimate the viscosity by measuring the velocity, since
$\displaystyle\frac{4\pi}{3}a^3\rho_sg =6\pi \eta a v$
where the radius can be directly measured.
Here you can observe the behavior of the sphere:
ID:(9871, 0)
Column emptying experiment
Description
This means that as the column empties and the height $h$ decreases, the velocity $v$ also decreases proportionally.
The key parameters are:
• Inner diameter of the vessel: 93 mm
• Inner diameter of the evacuation channel: 3 mm
• Length of the evacuation channel: 18 mm
These parameters are important to understand and analyze the process of column emptying and how the exit velocity varies with height.
ID:(9870, 0)
Column emptying experiment: viscosity effect
Description
If we analyze the equation
$ v = \displaystyle\frac{ \rho g R ^2}{8 \eta \Delta L } h $ |
which describes the application of Hagen-Poiseuille, we observe that the curve only matches the experimental data under the following conditions:
The velocity is low (when the column is nearly empty).
The radius of the evacuation channel must be reduced from 1.5 mm to 0.6 mm.
This shows that the flow is primarily turbulent and that only at low velocity levels is the velocity low enough for the Reynolds number to be low and the flow to be laminar.
ID:(11065, 0)
Viscosity calculation
Image
If we observe the path that the little ball takes over time, we can see that it mostly moves at a constant speed of about 0.31 meters in 25 seconds, which is equivalent to 0.0124 m/s.
By rearranging the equality between the gravitational force and the Stokes' resistance force:
$\displaystyle\frac{4\pi}{3}a^3\rho_sg =6\pi \eta a v$
we obtain:
$\eta = \displaystyle\frac{2 a^2\rho_sg}{9 v}$
Considering that the ball has a radius of 2 mm and weighs 8 mg, we can determine its density to be approximately $\rho_s\sim 2.38 g/cm^3$. Therefore, the viscosity is estimated to be around $\eta\sim 1.67, Pa s$.
ID:(9881, 0)
Model
Concept
Variables
Parameters
Selected parameter
Calculations
Equation
$h = h_0 e^{-t/\tau_{hp}}$
h = h_0 *exp(- t / tau_hp )
$ S \displaystyle\frac{dh}{dt} = - \displaystyle\frac{\pi R^4}{\Delta L}\displaystyle\frac{\rho g}{8 \eta} h $
S *DIFF(h,t,1) = - pi * R ^4 * rho * g * h /(8* eta * DL )
$ S\displaystyle\frac{dh}{dt} = \pi R^2 \sqrt{2gh}$
S*DIFF(h,t,1) = pi * R ^2*sqrt(2* g * h )
$ \tau_{hp} = \displaystyle\frac{S \Delta L}{\pi R^4}\displaystyle\frac{8 \eta}{\rho g}$
tau_hp = (8* eta * DL * S )/( pi * R ^4 * rho * g * h )
$ v = \displaystyle\frac{ \rho g R ^2}{8 \eta \Delta L } h $
v = rho * g * R ^2* h /(8* eta * l )
ID:(15494, 0)
Characteristic time column with viscous liquid
Equation
If we examine the equation for the draining of a viscous liquid column:
$ S \displaystyle\frac{dh}{dt} = - \displaystyle\frac{\pi R^4}{\Delta L}\displaystyle\frac{\rho g}{8 \eta} h $ |
we can condense the constants into a characteristic time unit:
$ \tau_{hp} = \displaystyle\frac{S \Delta L}{\pi R^4}\displaystyle\frac{8 \eta}{\rho g}$ |
This value becomes a characteristic time for the draining of a column with cross-sectional area $S$ of a viscous liquid with viscosity $\eta," as it flows through a tube of radius $R$.
ID:(14521, 0)
Non-viscous liquid column height over time
Equation
For the case of a non-viscous liquid flowing in a laminar fashion, the pressure difference generated by the column is:
$ \Delta p = \rho_w g \Delta h $ |
This results in a velocity flow $v$ through a tube according to Bernoulli's principle:
$ \Delta p = - \rho \bar{v} \Delta v $ |
Given the velocity and the tube's radius, we can calculate the flow, which is related to the flow within the column through the law of continuity. In turn, this is connected to the variation in height $h," as described in:
$ S\displaystyle\frac{dh}{dt} = \pi R^2 \sqrt{2gh}$ |
Using Bernoulli's equation, we can analyze the case of a column of water that generates a pressure difference:
$ \Delta p = \rho_w g \Delta h $ |
and induces a velocity flow $v$ through a tube, in accordance with:
$ \Delta p = - \rho \bar{v} \Delta v $ |
Thus, we can estimate the velocity as:
$v = \sqrt{2 g h}$
This velocity, through a tube section of radius $R$, results in a flow:
$J = \pi R^2 v$
If the column has a cross-sectional area $S$, and its height decreases with respect to the variation in height $h$ over time $t$, we can apply the law of continuity, which states:
$ S_1 j_{s1} = S_2 j_{s2} $ |
Therefore, the equation that describes this situation is:
$ S\displaystyle\frac{dh}{dt} = \pi R^2 \sqrt{2gh}$ |
ID:(9882, 0)
Column emptying experiment: model with Bernoulli
Description
Let's consider the system of a cylindrical bucket with a drainage hole. When the plug is removed, the water starts to flow due to the existing pressure. According to Bernoulli\'s principle, inside the bucket ($v\sim 0$), the velocity is zero, and we have:
$\displaystyle\frac{1}{2}\rho v^2 + \rho g h\sim \rho g h$
while outside the bucket ($h=0$), only the kinetic component exists:
$\displaystyle\frac{1}{2}\rho v^2 + \rho g h\sim \displaystyle\frac{1}{2}\rho v^2$
Since both expressions are equal, we have:
$\displaystyle\frac{1}{2}\rho v^2=\rho g h$
which gives the velocity as:
$v=\sqrt{2 g h}$
To compare with the experiment, we can use this expression to estimate, with:
$ S\displaystyle\frac{dh}{dt} = \pi R^2 \sqrt{2gh}$ |
the range that the stream should have. If we plot it graphically, we observe:
where:
• the red dots correspond to the experimental measurements,
• the blue dots correspond to the calculated range using a factor of 0.11,
• the transparent dots correspond to the calculated range using a factor of 0.09.
Therefore, we can conclude that Bernoulli\'s model overestimates the velocity at which the bucket empties. This is because in the vicinity of the drainage hole, the effects of viscosity are not negligible, and therefore, the velocity is lower.
ID:(11063, 0)
Column emptying experiment: Model with Hagen Poiseuille
Equation
Given the model for the flow of a viscous liquid through a tube and considering that the height of the column determines the pressure, we can estimate the velocity as a function of the column height:
$ v = \displaystyle\frac{ \rho g R ^2}{8 \eta \Delta L } h $ |
If we consider that the drainage channel presents hydraulic resistance, we can model it using the Hagen-Poiseuille equation:
$ J_V =-\displaystyle\frac{ \pi R ^4}{8 \eta }\displaystyle\frac{ \Delta p }{ \Delta L }$ |
where the pressure difference is determined by the water column:
$ p = \rho_w g h $ |
and the velocity is obtained through the flow:
$ j_s = \displaystyle\frac{ J_V }{ S }$ |
In this way, we obtain the relationship for calculating the velocity as a function of height:
$ v = \displaystyle\frac{ \rho g R ^2}{8 \eta \Delta L } h $ |
ID:(11064, 0)
Viscous liquid column height over time
Equation
The laminar flow of a fluid with viscosity $\eta$ through a tube of radius $R$ is described by Hagen-Poiseuille's law:
$ J_V =-\displaystyle\frac{ \pi R ^4}{8 \eta }\displaystyle\frac{ \Delta p }{ \Delta L }$ |
The pressure difference is determined by the height of the column $\Delta h$:
$ \Delta p = \rho_w g \Delta h $ |
which decreases as the liquid flows out. By applying the continuity equation, we can demonstrate that the height decreases over time as follows:
$ S \displaystyle\frac{dh}{dt} = - \displaystyle\frac{\pi R^4}{\Delta L}\displaystyle\frac{\rho g}{8 \eta} h $ |
If the flow through the tube is described by the equation:
$ J_V =-\displaystyle\frac{ \pi R ^4}{8 \eta }\displaystyle\frac{ \Delta p }{ \Delta L }$ |
and the pressure difference $\Delta p$ is proportional to the height of the column $\Delta h = h:
$ \Delta p = \rho_w g \Delta h $ |
we can apply the conservation of flow $J_{V1}=J_V$ between the tube and the column $J_{V2}$:
$ J_{V1} = J_{V2} $ |
,
where the flow in column $J_{V2}$ with cross-sectional area $S$ is given by:
$ j_s = \displaystyle\frac{ J_V }{ S }$ |
Here, the flux density $j_s$ corresponds to the average velocity, which is equal to the rate of change of height over time:
$j_s = \displaystyle\frac{dh}{dt}$
In this way, we obtain the equation for the height of the column as a function of time:
$ S \displaystyle\frac{dh}{dt} = - \displaystyle\frac{\pi R^4}{\Delta L}\displaystyle\frac{\rho g}{8 \eta} h $ |
ID:(14520, 0)
Temporal evolution of the viscous liquid column
Equation
The equation that describes the evolution of the viscous liquid column draining is as follows:
$ S \displaystyle\frac{dh}{dt} = - \displaystyle\frac{\pi R^4}{\Delta L}\displaystyle\frac{\rho g}{8 \eta} h $ |
We can rewrite this equation using the characteristic time:
$ \tau_{hp} = \displaystyle\frac{S \Delta L}{\pi R^4}\displaystyle\frac{8 \eta}{\rho g}$ |
Upon integration, this yields:
$h = h_0 e^{-t/\tau_{hp}}$ |
If in the equation
$ S \displaystyle\frac{dh}{dt} = - \displaystyle\frac{\pi R^4}{\Delta L}\displaystyle\frac{\rho g}{8 \eta} h $ |
the constants are replaced by
$ \tau_{hp} = \displaystyle\frac{S \Delta L}{\pi R^4}\displaystyle\frac{8 \eta}{\rho g}$ |
we obtain the first-order linear differential equation
$\displaystyle\frac{dh}{dt}=\displaystyle\frac{1}{\tau_{hp}} h$
whose solution is
$h = h_0 e^{-t/\tau_{hp}}$ |
where $h_0$ represents the initial height.
ID:(14522, 0)