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

Concept


ID:(15487, 0)



Static and dynamic pressure

Description

When you have four columns with different cross-sections interconnected, the liquid will assume the same level in all of them. If you open the interconnecting channel, the liquid will start to flow towards the opening where the pressure is equal to the ambient pressure. In the first cylinder, the pressure is equal to the pressure of the water column plus the atmospheric pressure, so the difference with respect to the pressure at the exit is the pressure of the first column. The liquid begins to gain velocity while the dynamic pressure starts to decrease, which is evident in the increasingly smaller columns.

ID:(11092, 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 mmThese parameters are important to understand and analyze the process of column emptying and how the exit velocity varies with height.

ID:(9870, 0)



Model

Concept


ID:(15490, 0)



In liquid column

Model

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.

Variables

Symbol
Text
Variable
Value
Units
Calculate
MKS Value
MKS Units
$\tau_b$
tau_b
Characteristic time column with Bernoulli
s
$h$
h
Column height
m
$\rho$
rho
Density
kg/m^3
$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
$\rho_w$
rho_w
Liquid density
kg/m^3
$\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
$R$
R
Tube radius
m
$\Delta p$
Dp
Variación de la Presión
Pa
$\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 

Symbol
Equation
Solved
Translated

Calculations

Symbol
Equation
Solved
Translated

 Variable   Given   Calculate   Target :   Equation   To be used



Equations

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

$ dp = p - p_0 $



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)

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)

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$

(ID 14524)

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)

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)


Examples


(ID 15487)

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 mmThe initial liquid height is 25 cm.

(ID 11092)

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)

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)


(ID 15490)


ID:(1427, 0)