Abstracts
The dynamics of an ideal fluid contained in two communicating vessels is studied. Despite the fact that the static properties of this system have been known since antiquity, the knowledge of the dynamical properties of an ideal fluid flowing in two communicating vessels is not similarly widespread. By means of Bernoulli's equation for non-stationary fluid flow, we study the oscillatory motion of the fluid when dissipation can be neglected.
Keywords:
Bernoulli's equation; non-stationary fluid; fluid dynamics
Estuda-se a dinâmica de um fluido ideal contido dentro de dois vasos comunicantes. Embora as propriedades estáticas desse sistema sejam conhecidas desde a Antiguidade, são bem menos conhecidas as propriedades dinâmicas do fluxo de dois fluidos ideais contidos vasos comunicantes. Utilizando a equação de Bernoulli para um fluxo não estacionário, nós estudamos o movimento oscilatório do fluido na ausência de dissipação.
Palavras-chave:
equação de Bernoulli; fluxo não estacionário; mecânica dos fluidos
1. Introduction
The topic of communicating vessels is often adopted as a common example in physics teaching [1[1] D. Halliday, R. Resnick and J. Walker, Fundamentals of Physics (Wiley and Sons, New York, 2005), 7th ed.]. The static property of this system is also used in other subject areas as, for example, sociology and economics [2[2] F. van Waarden, in: Advancing Socio-Economics: An Institutionalist Perspective, edited by J.R. Hollingsworth, K.H. Müller and E.J. Hollingsworth (Rowman and Littlefield Publishers, Lanham, 2005), chapter 9.] and in metaphoric expressions. In fact, it is a widespread common knowledge that a fluid in adjacent containers reaches the same height, measured with respect to a common reference point, independently of the shape of the vessels. However, reference to the dynamical properties of these types of physical systems in physics teaching books is rarely found. One of this rare examples can be found in ref. [3[3] P.A. Tipler, Physics (Worth Publishers Inc., New York, 1980), chapter 14.], where it is reported that the height of an ideal fluid contained in a U-shaped tube is seen to oscillate harmonically with a frequency given by:
being the total length of the fluid in the tube. Therefore, although the dynamics of these systems have been studied in quantum mechanical and classical contexts [4[4] R.J. Donnelly and O. Penrose, Phys. Rev. 103, 1137 (1956).–6[6] S.L.T. de Souza, I.L. Caldas, R.L. Viana, J.M. Balthazar and R.M.L.R.F. Brasil, Journal of Sound and Vibration 289, 987 (2006).], it is still not much studied at elementary levels. In particular, we would like to mention a specific technological application of this system in structural engineering [6[6] S.L.T. de Souza, I.L. Caldas, R.L. Viana, J.M. Balthazar and R.M.L.R.F. Brasil, Journal of Sound and Vibration 289, 987 (2006).]. In this specific field, U-tube-like containers, called tuned-liquid column dampers, are used to mitigate earth-quake induced vibrations in tall buildings.
In the present work we thus consider the dynamics of two communicating vessels of unequal cross-sectional areas, say and , as those shown in Fig.1. Therefore, the mechanical properties of the system shown in Fig. 1 are analyzed by means of Bernoulli's equation for non-stationary states. We write the nonlinear differential equations governing the motion of the system and point out that the result in equation (1) can be obtained by taking two equal sections; i.e., by writing . In this case a harmonic solution is obtained. We also notice, by a phase plane numerical analysis, that the solutions for the heights and is still periodic in the most general case . In this general dynamical system the fixed point is identified with the static equilibrium solution; i.e., the common height of the two fluids. Further generalizations are proposed as, for example, the possibility of allowing leakage from the system.
Communicating vessels containing an ideal fluid out of equilibrium. The first vessel has cross-sectional area and the fluid column reaches an height . The second vessel has cross-sectional area and the fluid in it is at an height . The connection channel is positioned at the base of the two vessels and is supposed to be sufficiently small.
2. Preliminary dynamical notions
We start by writing the continuity equation for the fluid [5[5] M.J. Smith, J.J. Kobine and F.A. Davidson, Proc. R. Soc. A 464, 905 (2008).]:
where is the flux flow rate in the connecting channel, as shown in Fig. 1. Summing the above equations, we have:
In equation (3) the static solution is implicitly defined. In fact, by integrating with respect to time, we get:
where is the equilibrium height.
In order to tackle the dynamical problem, we need to specify the type of fluid we are dealing with. In what follows we shall thus assume that the fluid is ideal and that Bernoulli's equation, written for the non-stationary regime of the fluid, holds. A derivation of this equation can be found in ref. [7[7] R. De Luca and P. Desideri, Eur. J. Phys. 34, 189 (2013).]. Naturally, different behaviour is expected when the fluid is not ideal. We therefore write, for two points in the fluids, one at an height , corresponding to the free liquid surface in the left container, the second at , corresponding to point in Fig. 1:
Setting and in equation (5), we have:
where we have assumed that the velocity profile does not appreciably vary with the height and where the dot stands for “derivative with respect to time”. Similarly, for the second vessel, by considering one point corresponding to the free liquid surface and the other at , we may set:
In this second case, we have considered and . By now combining equations (6) and (7), we obtain:
Because of equation (3), we have , with the obvious definition of the parameter . On the same token, by equation (4), we may also set
Therefore, equation (8) can be written only in terms of the variable as follows:
If we now go in horizontal from point to point , assuming the connection channel to be sufficiently small, equation (5) gives:
so that, because and , we may rewrite equation (11) as follows:
Finally, by combining equation (10) and equation (12), we have:
By now making the change variable , we may rewrite equation (13) as follows:
For we find the already known result [3[3] P.A. Tipler, Physics (Worth Publishers Inc., New York, 1980), chapter 14.]:
In this way, for , the liquid is seen to oscillate with a period . A trivial solution can be found for . For , on the other hand, we can proceed in normalizing the variable and as follows:
In this way, we obtain the following non-linear second-order ordinary differential equation:
where , and the prime stands for “derivative with respect to ”. The above equation can be cast in the form of two first-order differential equations as follows:
Equations (18a) and (18b), under the condition , are sufficient to describe the dynamics of the system. Notice that, for , the harmonic oscillator equation is obtained from (18a) and (18b).
3. Analysis of the dynamical system
By considering equation (18a) and (18b) we can start by describing the dynamical propertied of the system by means of a phase-plane analysis. In fact, by noticing that equation (15) gives a periodic solution of the problem, we expect that, at least for small values of the quantity , the solution of equation (17) to be periodic with an attraction point corresponding to . In order to detect this property, let us consider the phase-plane representation of the coupled first-order ordinary differential equations (18a) and (18b). By considering the initial conditions
we preliminary notice that the value of must fall within the interval . In fact, the dimensionless quantity represents the normalized deviation from the equilibrium position of the height of the liquid in the first column. In this way, we take the fluid to be initially at rest and completely contained in the right vessel or in the left vessel for and , respectively. As shown in Fig. 2a and 2b, obtained by means of the software Mathematica, the system shows periodical solutions, since the curves are closed. In Fig. 2a we report solutions for and for different values of , corresponding to different colours. In Fig. 2b, on the other hand, solutions for are shown for the same values of as in Fig. 2a. The black dashed line represents, in both Figs 2a and 2b, the harmonic solution to equation (15). We notice that, for decreasing values of , the vs. curves tend to collapse on the harmonic solution.
Phase-plane representation of the solution of the problem of communicating vessels. The initial conditions are the following: a) =1.0, v=0; b) =0.75, v=0. The colours magenta, red, orange, green cyan correspond to the following values of : 0.15, 0.25, 0.40, 0.55, 0.70. The dashed black circle corresponds to the solution for =1. Periodicity of the solution can be argued by the closing of all curves.
Another useful hint on the periodicity of the solutions of equation (17) can be given by perturbation analysis. In fact, by writing
where and are functions to be found, we may rewrite equation (17) as follows:
where second and higher order terms in have been neglected. By carrying out further algebra and by still retaining only first-order terms in , we write:
Therefore, the above equation splits into two different parts, one for the unperturbed solution only, one for the perturbed term , so that, we have:
Notice that the differential equation (23a) describes free harmonic oscillations, whose solution is
where the constants and are to be determined by means of the initial conditions (19) applied to the solution . However, notion of the solution to equation (23a) is necessary to solve equation (23b). In this way, by substituting equation (24) into equation (23b) we obtain a differential equation describing the motion of a forced harmonic oscillator:
where a higher harmonic term appears. A general solution to equation (25) can be found by writing it as the sum of the homogeneous and particular solutions as follows:
In this way, the approximate solution is, according to equation (20), the following:
where , , , and need to be evaluated according to (19), so that:
We thus find the following solution by separating the zero-order and the first-order terms:
Finally, we write equation (27) as follows:
By taking the derivative with respect to , we can also write:
We can thus rewrite these solutions in terms of the original variable :
By taking the derivative with respect to , we can also write:
The constant is defined in equation (19) and the quantities and can be known by means of equation (9).
The above solutions are periodic. In fact, we can now prove that the maximum deviation from a circular path in the phase plane of the functions and is a first-order term in . By considering equation (30) and (31) we can see that
Let us now consider the radial distance from the circumference in the phase plane and the curve in (34):
To first order in , we can write:
The maximum value of can be found by maximizing the time dependent term in (36), obtaining:
Therefore, the effects due to non-linearity of the ordinary differential equation (17) are rather weak when very small values of are considered. A comparison between the numerical solutions of equation (18a) and (18b) and the analytical solutions (30) and (31) obtained by means of the above perturbation approach for 0.1 and 0.2 is shown in Fig 3a and 3b. The latter figures, as for Fig. 2a and 2b, have been obtained by means of the software Mathematica.
Phase-plane representation of the numerical (dashed) and approximated (full line) solution. The latter solution has been obtained by means of a perturbation analysis to first order in the parameter . The following values of have been chosen: =0.10 (a) and =0.20 (b).
4. Conclusions
The celebrated system of communicating vessels, which also finds applications in other subjects than physics, has been seen to possess interesting dynamical properties. In the present work we have considered two cylindrical vessels containing a non-viscous fluid. The two cylinders are connected by a small horizontal channel placed at their base. In the symmetric case, in which the two containers are perfectly identical, the system behaves like a harmonic oscillator, whose angular frequency is given by the following simple relation:
where is the height, measured with respect to the bottom of the vessels, up to which the free surface of the liquid rises at equilibrium. Surprisingly, when we consider the asymmetric case, harmonic oscillations recede because of non-linear effects. The resulting non-linear differential equation preserves periodicity, as both a phase-plane analysis and a perturbation approach suggest.
As for the didactical value of the present work, we may affirm that the present study could be useful to students approaching the topic of nonlinear dynamical systems for the first time. In fact, the non-linear effects in the dynamical properties of the system can be seen either by means of a phase-plane analysis or by a perturbation approach. Both methods confirm periodic behavior of the system in the presence of non-linearity. Moreover, as a starting point students can rely on the harmonic oscillator mechanical equivalent system when the parameter , defining the degree of asymmetry between the sections and of the two cylindrical vessels, is zero. Nonlinear effects can be seen to gradually arise as increases. Finally, students may notice that, while analytic solution of the problem is possible by means of a perturbation approach for small values of the parameter , numerical analysis is necessary for value of approaching one.
Referências
-
[1]D. Halliday, R. Resnick and J. Walker, Fundamentals of Physics (Wiley and Sons, New York, 2005), 7th ed.
-
[2]F. van Waarden, in: Advancing Socio-Economics: An Institutionalist Perspective, edited by J.R. Hollingsworth, K.H. Müller and E.J. Hollingsworth (Rowman and Littlefield Publishers, Lanham, 2005), chapter 9.
-
[3]P.A. Tipler, Physics (Worth Publishers Inc., New York, 1980), chapter 14.
-
[4]R.J. Donnelly and O. Penrose, Phys. Rev. 103, 1137 (1956).
-
[5]M.J. Smith, J.J. Kobine and F.A. Davidson, Proc. R. Soc. A 464, 905 (2008).
-
[6]S.L.T. de Souza, I.L. Caldas, R.L. Viana, J.M. Balthazar and R.M.L.R.F. Brasil, Journal of Sound and Vibration 289, 987 (2006).
-
[7]R. De Luca and P. Desideri, Eur. J. Phys. 34, 189 (2013).
Publication Dates
-
Publication in this collection
2017
History
-
Received
22 Dec 2016 -
Reviewed
26 Jan 2017 -
Accepted
05 Feb 2017