In liquid column

Storyboard

In the case of a column of liquid, Bernoulli's law can be applied along with the hydrostatic pressure term. However, it's important to note that when viscosity of the liquid is not considered, the reduction in level occurs uniformly. In this regard, it can be modeled using the continuity equation to determine the downward velocity of the cylinder.

For a column of liquid with an outlet at the bottom, the behavior is similar to what is estimated with Bernoulli's principle. Differences arise due to the formation of small vortices at the outlet, effectively reducing the outlet area and obstructing the flow. However, the flow of a low-viscosity liquid can be modeled in the zone without vortices using Bernoulli's principle.

>Model

ID:(1427, 0)



Mechanisms

Iframe

>Top



Code
Concept
Liquid column exit velocity

Mechanisms

ID:(15487, 0)



Liquid column exit velocity

Description

>Top


If there is a column height ($h$) of liquid with the liquid density ($\rho_w$) under the effect of gravity, using the gravitational Acceleration ($g$), the variación de la Presión ($\Delta p$) is generated according to:

$ \Delta p = \rho_w g \Delta h $



This the variación de la Presión ($\Delta p$) produces a flow through the outlet tube with the tube length ($\Delta L$), the tube radius ($R$), and the viscosity ($\eta$) of a volume flow 1 ($J_{V1}$) according to the Hagen-Poiseuille law:

$ J_V =-\displaystyle\frac{ \pi R ^4}{8 \eta }\displaystyle\frac{ \Delta p }{ \Delta L }$



Since this equation includes the section in point 2 ($S_2$), the flux density 2 ($j_{s2}$) can be calculated using:

$ j_s = \displaystyle\frac{ J_V }{ S }$



With this, we obtain:

$ j_s = \displaystyle\frac{ \rho_w g R ^2}{8 \eta \Delta L } h $



which corresponds to an average velocity.

To model the system, the key parameters are:

• Interior diameter of the container: 93 mm

• Interior diameter of the evacuation channel: 3.2 mm

• Length of the evacuation channel: 18 mm

The initial liquid height is 25 cm.

ID:(11092, 0)



Column emptying experiment: model with Bernoulli

Description

>Top


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: Height and Reach

Description

>Top


If the Tracker program is used, the height of the meniscus of the column and the range of the jet can be measured. The relationship between the two is shown in the following graph:

Reach [m] vs Height [m]

The recorded data, which can be downloaded as an Excel table from the following link excel table, are as follows:

Time [s]Height [m]Reach [m]
02.23E-011.89E-01
42.14E-011.86E-01
82.04E-011.82E-01
121.94E-011.77E-01
161.86E-011.72E-01
201.79E-011.68E-01
241.71E-011.66E-01
281.63E-011.62E-01
321.54E-011.58E-01
361.46E-011.52E-01
401.39E-011.48E-01
441.32E-011.44E-01
481.24E-011.39E-01
521.18E-011.35E-01
561.11E-011.31E-01
601.06E-011.27E-01
649.88E-021.23E-01
689.29E-021.18E-01
728.70E-021.15E-01
768.11E-021.12E-01
807.52E-021.06E-01
847.12E-021.02E-01
886.51E-029.69E-02
926.00E-029.42E-02
965.58E-028.94E-02
1005.09E-028.52E-02
1044.70E-028.13E-02
1084.34E-027.63E-02
1123.97E-027.22E-02
1163.49E-026.79E-02
1203.15E-026.28E-02
1242.91E-025.96E-02
1282.58E-025.33E-02
1322.23E-024.92E-02
1361.98E-024.31E-02
1401.71E-023.85E-02
1441.54E-023.38E-02
1481.28E-022.85E-02
1521.11E-022.23E-02
1569.17E-031.54E-02
1607.15E-037.95E-03

Note: E indicates scientific notation (ej. 1.2E+3 = 1.2x10^3 = 1200, y 1.2E-3 = 1.2x10^-3 = 0.0012)

ID:(11062, 0)



Model

Top

>Top



Parameters

Symbol
Text
Variable
Value
Units
Calculate
MKS Value
MKS Units
$\tau_b$
tau_b
Characteristic time column with Bernoulli
s
$\rho$
rho
Density
kg/m^3
$g$
g
Gravitational Acceleration
m/s^2
$\rho_w$
rho_w
Liquid density
kg/m^3
$\pi$
pi
Pi
rad
$R$
R
Tube radius
m
$\Delta p$
Dp
Variación de la Presión
Pa

Variables

Symbol
Text
Variable
Value
Units
Calculate
MKS Value
MKS Units
$h$
h
Column height
m
$r$
r
Disc radius
m
$v_s$
v_s
Flow speed
m/s
$j_s$
j_s
Flux density
m/s
$j_{s1}$
j_s1
Flux density 1
m/s
$j_{s2}$
j_s2
Flux density 2
m/s
$\Delta h$
Dh
Height of liquid column
m
$h_0$
h_0
Initial height of liquid column
m
$\Delta p_s$
Dp_s
Pressure difference
Pa
$S_1$
S_1
Section in point 1
m^2
$S_2$
S_2
Section in point 2
m^2
$S$
S
Section Tube
m^2
$S$
S
Surface of a disk
m^2
$t$
t
Time
s
$\Delta t$
Dt
Time elapsed
s
$\Delta s$
Ds
Tube element
m
$\Delta V$
DV
Volume element
m^3
$J_V$
J_V
Volume flow
m^3/s

Calculations


First, select the equation: to , then, select the variable: to

Calculations

Symbol
Equation
Solved
Translated

Calculations

Symbol
Equation
Solved
Translated

Variable Given Calculate Target : Equation To be used




Equations

#
Equation

$ \Delta p = \rho_w g \Delta h $

Dp = rho_w * g * Dh


$ \Delta V = S \Delta s $

DV = S * Ds


$ h = h_0\left(1-\displaystyle\frac{t}{\tau_b}\right)^2$

h = h_0 *(1-t/tau_b)^2


$ j_s =\displaystyle\frac{ \Delta s }{ \Delta t }$

j_s = Ds / Dt


$ J_V =\displaystyle\frac{ \Delta V }{ \Delta t }$

J_V = DV / Dt


$ J_V = S j_s $

J_V = S * j_s


$ S = \pi r ^2$

S = pi * r ^2


$ S\displaystyle\frac{dh}{dt} = \pi R^2 \sqrt{2gh}$

S*DIFF(h,t,1) = pi * R ^2*sqrt(2* g * h )


$ S_1 j_{s1} = S_2 j_{s2} $

S_1 * j_s1 = S_2 * j_s2


$ \tau_b = \displaystyle\frac{S}{\pi R^2}\sqrt{\displaystyle\frac{h_0}{g}}$

tau_b = (S /( pi * R ^2))*sqrt(h_0/g)


$ v_s = \sqrt{\displaystyle\frac{2 \Delta p_s }{ \rho }}$

v_s = sqrt(2* Dp_s / rho )

ID:(15490, 0)



Volume element

Equation

>Top, >Model


If we have a tube with a the section Tube ($S$) moving a distance the tube element ($\Delta s$) along its axis, having displaced the volume element ($\Delta V$), then it is equal to:

$ \Delta V = S \Delta s $

$S$
Section Tube
$m^2$
6267
$\Delta s$
Tube element
$m$
10291
$\Delta V$
Volume element
$m^3$
10290

ID:(3469, 0)



Mean volume flow

Equation

>Top, >Model


The volume flow ($J_V$) corresponds to the volume flowing ($\Delta V$) flowing through the channel at the time elapsed ($\Delta t$). Therefore, we have:

$ J_V =\displaystyle\frac{ \Delta V }{ \Delta t }$

$\Delta t$
Time elapsed
$s$
5103
$\Delta V$
Volume element
$m^3$
10290
$J_V$
Volume flow
$m^3/s$
5448

ID:(4347, 0)



Average flow density

Equation

>Top, >Model


The flux density ($j_s$) is related to the distance traveled in a time ($\Delta s$), which is the distance that the fluid travels in the time elapsed ($\Delta t$), as follows:

$ j_s =\displaystyle\frac{ \Delta s }{ \Delta t }$

$j_s$
Flux density
$m/s$
7220
$\Delta t$
Time elapsed
$s$
5103
$\Delta s$
Tube element
$m$
10291

ID:(4348, 0)



Volume Flow and its Speed

Equation

>Top, >Model


A flux density ($j_s$) can be expressed in terms of the volume flow ($J_V$) using the section or Area ($S$) through the following formula:

$ J_V = S j_s $

$j_s$
Flux density
$m/s$
7220
$S$
Section Tube
$m^2$
6267
$J_V$
Volume flow
$m^3/s$
5448

The volume ($V$) for a tube with constant the section Tube ($S$) and a position ($s$) is

$ V = h S $



If the section Tube ($S$) is constant, the temporal derivative will be

$\displaystyle\frac{dV}{dt} = S\displaystyle\frac{ds}{dt}$



thus, with the volume flow ($J_V$) defined by

$ J_V =\displaystyle\frac{ dV }{ dt }$



and with the flux density ($j_s$) associated with the position ($s$) via

$ j_s =\displaystyle\frac{ ds }{ dt }$



it is concluded that

$ J_V = S j_s $

ID:(15716, 0)



Surface of a disk (1)

Equation

>Top, >Model


The surface of a disk ($S$) of ($$) is calculated as follows:

$ S = \pi r ^2$

$r$
Disc radius
$m$
5275
$\pi$
Pi
3.1415927
$rad$
5057
$S$
Surface of a disk
$m^2$
10361

ID:(3804, 1)



Speed with respect to rest

Equation

>Top, >Model


When applying the Bernoulli equation relative to a stationary point in the flow, it is established that the flow speed ($v_s$) is associated with the pressure difference ($\Delta p_s$) concerning the pressure at this point. With the density ($\rho$), the following is observed:

$ v_s = \sqrt{\displaystyle\frac{2 \Delta p_s }{ \rho }}$

$\rho$
Density
$kg/m^3$
5342
$v_s$
Flow speed
$m/s$
5821
$\Delta p_s$
Pressure difference
$Pa$
10117

In this case, it can be assumed that the mean Speed of Fluid in Point 2 ($v_2$) represents zero velocity and the mean Speed of Fluid in Point 1 ($v_1$) corresponds to the flow speed ($v_s$). Therefore, for the speed difference between surfaces ($\Delta v$) the following is established:

$\Delta v = v_2 - v_1 = 0 - v_s = - v_s$



and for the average speed ($\bar{v}$) it is calculated:

$\bar{v} = \displaystyle\frac{v_1 + v_2}{2} = \frac{v_s}{2}$



Consequently, with the variación de la Presión ($\Delta p$), which equals the pressure difference ($\Delta p_s$), we obtain:

$ \Delta p = - \rho \bar{v} \Delta v $



resulting in:

$\Delta p_s = \displaystyle\frac{1}{2} \rho v_s^2$



leading to:

$ v_s = \sqrt{\displaystyle\frac{2 \Delta p_s }{ \rho }}$

ID:(15710, 0)



Pressure difference between columns

Equation

>Top, >Model


The height difference, denoted by the height difference ($\Delta h$), implies that the pressure in both columns is distinct. In particular, the pressure difference ($\Delta p$) is a function of the liquid density ($\rho_w$), the gravitational Acceleration ($g$), and the height difference ($\Delta h$), as follows:

$ \Delta p = \rho_w g \Delta h $

$g$
Gravitational Acceleration
9.8
$m/s^2$
5310
$\Delta h$
Height of liquid column
$m$
5819
$\rho_w$
Liquid density
$kg/m^3$
5407
$\Delta p$
Variación de la Presión
$Pa$
6673

If there is the pressure difference ($\Delta p$) between two points, as determined by the equation:

$ \Delta p = p_2 - p_1 $



we can utilize the water column pressure ($p$), which is defined as:

$ p_t = p_0 + \rho_w g h $



This results in:

$\Delta p=p_2-p_1=p_0+\rho_wh_2g-p_0-\rho_wh_1g=\rho_w(h_2-h_1)g$



As the height difference ($\Delta h$) is:

$ \Delta h = h_2 - h_1 $



the pressure difference ($\Delta p$) can be expressed as:

$ \Delta p = \rho_w g \Delta h $

ID:(4345, 0)



Continuity by section

Equation

>Top, >Model


The principle of continuity dictates that the flow at the first point, which is equal to the flux density 1 ($j_{s1}$) times the section in point 1 ($S_1$), must be equal to the flow at the second point, given by the flux density 2 ($j_{s2}$) times the section in point 2 ($S_2$). From this, it follows that:

$ S_1 j_{s1} = S_2 j_{s2} $

$j_{s1}$
Flux density 1
$m/s$
10288
$j_{s2}$
Flux density 2
$m/s$
10289
$S_1$
Section in point 1
$m^2$
5257
$S_2$
Section in point 2
$m^2$
5413

ID:(4350, 0)



Non-viscous liquid column height over time

Equation

>Top, >Model


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}$

$h$
Column height
$m$
5406
$g$
Gravitational Acceleration
9.8
$m/s^2$
5310
$\pi$
Pi
3.1415927
$rad$
5057
$S$
Section Tube
$m^2$
6267
$t$
Time
$s$
5264
$R$
Tube radius
$m$
5417

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 characteristic time with non-viscous liquid

Equation

>Top, >Model


If we examine the equation for the draining of a non-viscous liquid column:

$ S\displaystyle\frac{dh}{dt} = \pi R^2 \sqrt{2gh}$



we can condense the constants into a characteristic time unit:

$ \tau_b = \displaystyle\frac{S}{\pi R^2}\sqrt{\displaystyle\frac{h_0}{g}}$

$\tau_b$
Characteristic time column with Bernoulli
$s$
10086
$g$
Gravitational Acceleration
9.8
$m/s^2$
5310
$h_0$
Initial height of liquid column
$m$
10085
$\pi$
Pi
3.1415927
$rad$
5057
$S$
Section Tube
$m^2$
6267
$R$
Tube radius
$m$
5417

This value becomes the time it takes for the column to completely empty, and it depends on the initial height.

ID:(14523, 0)



Temporal evolution of the column of non-viscous liquid

Equation

>Top, >Model


The equation describing the evolution of the draining viscous liquid column is as follows:

$ S\displaystyle\frac{dh}{dt} = \pi R^2 \sqrt{2gh}$



We can rewrite this equation in terms of the characteristic time:

$ \tau_b = \displaystyle\frac{S}{\pi R^2}\sqrt{\displaystyle\frac{h_0}{g}}$



Upon integrating, we obtain:

$ h = h_0\left(1-\displaystyle\frac{t}{\tau_b}\right)^2$

$\tau_b$
Characteristic time column with Bernoulli
$s$
10086
$h$
Column height
$m$
5406
$h_0$
Initial height of liquid column
$m$
10085
$t$
Time
$s$
5264

If in the equation

$ S\displaystyle\frac{dh}{dt} = \pi R^2 \sqrt{2gh}$



the constants are replaced by

$ \tau_b = \displaystyle\frac{S}{\pi R^2}\sqrt{\displaystyle\frac{h_0}{g}}$



we obtain the first-order linear differential equation

$\displaystyle\frac{dh}{dt}=\displaystyle\frac{1}{\tau_b} \sqrt{h_0 h}$



whose solution is

$ h = h_0\left(1-\displaystyle\frac{t}{\tau_b}\right)^2$

Where $h_0$ represents the initial height.

ID:(14524, 0)