# Liquid vibrations in a container with a membrane at the free surface

## Neelam Choudhary1, Elena Strelnikova2

1Department of Mathematics, Bennett University, Greater Noida, 201310, India

2A. N. Podgorny Institute of Mechanical Engineering Problems, National Academy of Sciences of Ukraine, Kharkiv, Ukraine

1Corresponding author

Vibroengineering PROCEDIA, Vol. 37, 2021, p. 13-18. https://doi.org/10.21595/vp.2021.21996
Received 13 April 2021; received in revised form 22 April 2021; accepted 29 April 2021; published 7 May 2021

51st International Conference on Vibroengineering, Dublin, Ireland, May 7-8, 2021

Copyright © 2021 Neelam Choudhary, et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Views 92
Abstract.

In the paper, the vibrations of liquid in a rigid and circular cylindrical container have been investigated. The effects of an elastic membrane placed on the free surface of the fluid have been analyzed on sloshing frequencies. The cylindrical container is assumed to be partially filled with an incompressible fluid. The potential theory is used to formulate the mathematical problem and set up the boundary value problem. The solutions are obtained for velocity potential and the deflection of elastic membrane in form of Fourier-Bessel series. The results are compared for the proposed analytical approach and BEM.

Highlights
• An analytical approach is used to solve the boundary value problem.
• Sloshing in presence of membrane on liquid free surface is discussed.
• Vibrations modes are reported.

Keywords: sloshing, vibration, membrane, cylindrical container.

#### 1. Introduction

“Sloshing is known as the free surface motion in a tank which is filled partially” . Slosh motion is a potentially dangerous situation to engineering structures and environment and can lead to the failure of structural units and their stability. The dynamic behaviour of structures carrying fuel tanks or storage reservoirs is significantly affected by fluid structure-interaction and is very dangerous to the safety and stability of the structure. It is necessary to control the vibrations of fluid-structure interaction to maintain the stability of the structures used in various engineering applications such as transporting liquid, petroleum reservoirs and space vehicles. Control of fluid sloshing inside a container has always been a challenge while designing any tank due to uninvited vibrations which are dangerous for the stability of the system. The Gas-free liquid is a key function of propellant tanks/containers in spacecraft and offshore oil industry where floating roofed oil tanks are used Propellant sloshing in tanks also influences the rigid body motion of the spacecraft or launch vehicles, which then needs to be controlled by the reaction control system.

Sloshing is a worrisome situation in oil storage tanks as well as in natural gas reservoirs. In recent years, floating liquified natural gas (FLNG) has gained popularity as innovative technology when it comes to natural gas exploitation. These floating platforms as well as huge tanks filled with LNG/ Oil are vulnerable to the vibrations due to ground motion during an earthquake. Suppression devices (rigid/elastic) are considered most effective to dampen the impact loads due to sloshing in a tank as shown in the studies reported in [1-6]. The effects of rigid baffles on sloshing frequencies in circular cylindrical containers are investigated analytically in [1, 2]. The non-linear sloshing problem is reported in  and stress analysis is investigated. The effects of vertically placed baffles have been reported in  and . Some studied on floating roof (considered as membrane) effects on sloshing frequencies since Tokachi-Oki (2003) earthquake are reported in [7-9], and . The effects of perforated baffles on sloshing are reported in  and . The vibrations of liquid in cylindrical tanks with and without baffles on the free surface under horizontal and vertical excitation are reported in . The problem of sloshing in horizontal elliptical tanks with T-shaped baffles is investigated in . Many researchers have been investigating the membrane’s dynamic response due to sloshing. In this paper, the fundamental sloshing modes of liquid with a membrane placed on the free surface are reported.

#### 2. Problem statement

Fig. 1 shows a schematic diagram of sloshing in a container with elastic membrane placed at the free surface. A vertical circular cylindrical tank with radius $R$ containing a fluid free surface at rest is considered at height $h$ from the rigid bottom of the container. To model the sloshing problem, the following assumptions are made: (i) the container is rigid and flow field is considered irrotational, (ii) the fluid is considered incompressible, (iii) small displacement of free surface of liquid are considered. The velocity potential theory is used to formulate the problem. The boundary value problem in terms of the velocity potential is given as follows.

Fig. 1. Schematic diagram of the problem For an irrotational flow motion, the velocity $V$ can be written as a gradient of potential $\varphi \left(r,\theta ,z,t\right)$ such that $V=\nabla \varphi$. Based on above assumptions of the potential theory, velocity potential satisfies the Laplace equation:

(1)

The boundary conditions at the rigid boundary of the cylinder are at the cylinder wall $r=R$ and at the rigid bottom:

(2)

The boundary condition for motion of membrane at the free surface is given by:

(3)
$\begin{array}{c}\frac{{\partial }^{2}\stackrel{-}{w}}{\partial {r}^{2}}+\frac{1}{r}\frac{\partial \stackrel{-}{w}}{\partial r}+\frac{1}{{r}^{2}}\frac{{\partial }^{2}\stackrel{-}{w}}{\partial {\theta }^{2}}-\frac{\mu }{T}\frac{{\partial }^{2}\stackrel{-}{w}}{\partial {t}^{2}}=-\frac{p}{T},\end{array}$

where $\stackrel{-}{w}$ denotes the deflection of the membrane, $T$ the tension per unit length, $\mu$ mass/unit area of the membrane and $p$ denotes the liquid pressure on the membrane. The membrane deflection satisfies the following boundary conditions at the rigid wall of the container $r=R$:

(4)
$\frac{\partial \stackrel{-}{w}}{\partial r}=\stackrel{-}{w}=0.$

The continuity of normal velocity components at the free surface $z=h$ should be satisfied:

(5)
$\frac{\partial \varphi }{\partial z}=\frac{\partial \stackrel{-}{w}}{\partial t},$

and the dynamic boundary condition at $z=h$ gives:

(6)
$p=-\rho \frac{\partial \varphi }{\partial z}-\rho g\stackrel{-}{w}.$

So, the coupled boundary value problem given by Eqs. (1)-(6), is required to solve to find the unknown functions $\varphi$ and $\stackrel{-}{w}$. First, we describe the analytical approach. The separation of variables method is used to find the velocity potential in the following form:

(7)
$\varphi \left(r,\theta ,z,t\right)=\sum _{m=0}^{\mathrm{\infty }}\sum _{n=1}^{\mathrm{\infty }}{A}_{mn}{J}_{m}\left(\frac{{k}_{mn}r}{R}\right)\mathrm{c}\mathrm{o}\mathrm{s}\mathrm{h}\left(\frac{{k}_{mn}z}{R}\right)\mathrm{c}\mathrm{o}\mathrm{s}m\theta {e}^{i\omega t},$

where ${k}_{mn}$ are roots of ${J\mathrm{\text{'}}}_{m}\left(kR\right)=0$. Using velocity potential function for Eq. (6), have the expression for pressure can be written as:

(8)
$p={e}^{i\omega t}\left[-pi\omega \sum _{m=0}^{\mathrm{\infty }}\sum _{n=1}^{\mathrm{\infty }}{A}_{mn}{\mathrm{J}}_{m}\left(\frac{{k}_{mn}r}{R}\right)\mathrm{c}\mathrm{o}\mathrm{s}\left(m\theta \right)\mathrm{c}\mathrm{o}\mathrm{s}\mathrm{h}\left(\frac{{k}_{mn}h}{R}\right)\right]-\rho g\stackrel{-}{w}.$

The thickness of the membrane is considered negligible. Eq. (4) gives rise to:

(9)
$\frac{{\partial }^{2}{w}_{1}}{\partial {r}^{2}}+\frac{1}{r}\frac{\partial {w}_{1}}{\partial r}+\frac{1}{{r}^{2}}\frac{{\partial }^{2}{w}_{1}}{\partial {\theta }^{2}}+\left(\frac{\mu {\stackrel{-}{\omega }}^{2}-g\rho }{T}\right){w}_{1}=\frac{i\rho \stackrel{-}{\omega }}{T}\sum _{m=0}^{\mathrm{\infty }}\sum _{n=1}^{\mathrm{\infty }}{A}_{mn}{\mathrm{J}}_{m}\left(\frac{{k}_{mn}r}{R}\right)\mathrm{c}\mathrm{o}\mathrm{s}\left(m\theta \right)\mathrm{c}\mathrm{o}\mathrm{s}\mathrm{h}\left(\frac{{k}_{mn}h}{R}\right),$

which is a nonhomogeneous differential equation in $\stackrel{-}{w}$. Solution for the deflection of membrane is following:

(10)
$\stackrel{-}{w}\left(r,\theta ,t\right)={e}^{i\omega t}\left[\sum _{m=0}^{\mathrm{\infty }}\sum _{n=1}^{\mathrm{\infty }}{A}_{mn}\left(\frac{i\rho \omega {R}^{2}}{T}\right)\mathrm{c}\mathrm{o}\mathrm{s}\mathrm{h}\left(\frac{{k}_{mn}h}{R}\right)\frac{{J}_{m}\left(\frac{{k}_{mn}r}{R}\right)\mathrm{c}\mathrm{o}\mathrm{s}m\theta }{\left[{c}^{2}-{k}_{mn}^{2}\right]}\right$

Using Eq. (5) at $z=h$, we obtain:

(11)
$\sum _{m=0}^{\mathrm{\infty }}\sum _{n=1}^{\mathrm{\infty }}{A}_{mn}\frac{{k}_{mn}}{R}{J}_{m}\left(\frac{{k}_{mn}r}{R}\right)\mathrm{s}\mathrm{i}\mathrm{n}\mathrm{h}\left(\frac{{k}_{mn}h}{R}\right)\mathrm{c}\mathrm{o}\mathrm{s}m\theta$

Boundary condition 0, at $r=R$ gives:

(12)
$\sum _{m=0}^{\mathrm{\infty }}\sum _{n=1}^{\mathrm{\infty }}{A}_{mn}\left(\frac{i\rho \omega {R}^{2}}{T}\right)\frac{{J}_{m}\left({k}_{mn}\right)\mathrm{c}\mathrm{o}\mathrm{s}m\theta }{\left[{c}^{2}-{k}_{mn}^{2}\right]}\mathrm{c}\mathrm{o}\mathrm{s}\mathrm{h}\left(\frac{{k}_{mn}h}{R}\right)+\sum _{m=0}^{\mathrm{\infty }}{B}_{m}{J}_{m}\left(c\right)\mathrm{c}\mathrm{o}\mathrm{s}m\theta =0.$

To determine the unknowns ${A}_{mn}$, ${B}_{m}$ and unknown frequency $\stackrel{-}{\omega }$, given by the Eqs. (11) and (12) is tedious task. This problem was solved using method developed in [1, 2]. For simulations using BEM, first we consider the problem of free vibrations of membrane without interaction with the liquid. We have:

(13)
$\frac{{\partial }^{2}\stackrel{-}{w}}{\partial {r}^{2}}+\frac{1}{r}\frac{\partial \stackrel{-}{w}}{\partial r}+\frac{1}{{r}^{2}}\frac{{\partial }^{2}\stackrel{-}{w}}{\partial {\theta }^{2}}-\frac{\mu }{T}\frac{{\partial }^{2}\stackrel{-}{w}}{\partial {t}^{2}}=0,$

with $\frac{\partial \stackrel{-}{w}}{\partial r}=\stackrel{-}{w}=0$.

Supposing that $\stackrel{-}{w}={w}_{1}\left(r,\theta \right){e}^{i\stackrel{-}{\omega }t}$, one can obtain:

(14)

The solutions of this problem are own modes and own frequencies, ${\stackrel{-}{\omega }}_{mk}\text{,}$${w}_{1mk}={w}_{mk}\left(r\right)\mathrm{c}\mathrm{o}\mathrm{s}m\theta$.

These modes and frequencies are reported in . To obtain the solution of coupled problem the next series is used for each $m$:

(15)
$\stackrel{-}{w}\left(r,\theta ,z,t\right)=\mathrm{c}\mathrm{o}\mathrm{s}m\theta {\sum }_{k=1}^{{N}_{w}}{c}_{k}\left(t\right){w}_{mk}\left(r\right).$

Then the next representation for function $\varphi$ is obtained:

(16)
$\varphi \left(r,\theta ,z,t\right)=\mathrm{c}\mathrm{o}\mathrm{s}m\theta {\sum }_{k=1}^{{N}_{w}}{\stackrel{˙}{c}}_{k}\left(t\right){\stackrel{-}{\mathrm{\Phi }}}_{mk}\left(r\right).$

We obtain the system of ODEs of second order for each $m$ in the following form:

(17)
${A}^{m}\stackrel{¨}{C}+{B}^{m}C=0,$

where:

For receiving functions ${\stackrel{-}{\mathrm{\Phi }}}_{mk}$ the boundary element method is in used . If we suppose that ${c}_{k}\left(t\right)=\mathrm{e}\mathrm{x}\mathrm{p}\left(i\omega t\right)$, we get the eigenvalue problem. We will seek the harmonic functions ${\stackrel{-}{\mathrm{\Phi }}}_{mk}$ in the form of the sum of simple and double layer potentials .

We shall indicate a moistened surface of a shell through ${S}_{1}$. The free surface of the liquid ${S}_{0}$ coincides with the plane $xy\left(orr,\theta \right)$ in unperturbed state. So, we have the singular integral equations for receiving ${\stackrel{-}{\mathrm{\Phi }}}_{mk}$ in the following form:

(18)
$2\pi \stackrel{-}{\mathrm{\Phi }}\left({P}_{0}\right)+\underset{{S}_{1}}{\iint }\stackrel{-}{\mathrm{\Phi }}\frac{\partial }{\partial n}\frac{1}{\left|P-{P}_{0}\right|}d{S}_{1}=\underset{{S}_{0}}{\iint }w\frac{1}{\left|P-{P}_{0}\right|}d{S}_{0}.$

This equation is reduced to one-dimensional one as in . For numerical simulation, the boundary elements with constant density are applied. Here, ${N}_{0}$ is the number of boundary elements along the free surface radius; ${N}_{w}$ along the shell wall, and ${N}_{bot}$ along the shell bottom. Consider the rigid circular cylindrical shell with the radius and height as $R=$ 1 m, and 1 m, respectively. Table 1 below provides the numerical values of the natural frequencies of liquid sloshing for 1.

Table 1. Slosh frequency parameters ${\omega }_{mk}^{2}/g$ of the fluid-filled rigid cylindrical shell

 $m$ BEM 1 2 3 4 5 1 ${N}_{0}$ ${N}_{w}$ ${N}_{bot}$ 1 2 3 4 5 25 20 20 1.6590 5.3301 8.5385 11.7071 14.8684 50 40 40 1.6579 5.3297 8.5372 11.7082 14.8655 100 80 80 1.6573 5.3293 8.5366 11.7066 14.8635 Analytical solution 1.6573 5.3293 8.5363 11.706 14.8635

The results of Table 1 testify convergence of proposed BEM. In should be noted that the accuracy $\epsilon =$ 10-4 has been achieved here for ${N}_{0}={N}_{w}={N}_{bot}=$ 80.

#### 3. Results

Several numerical experiments have been performed to validate the accuracy of the proposed analytical method and determine the frequency. Eqs. (11) and (12) are solved to determine frequencies $\omega$.

Table 2. Validation of numerical and analytical results

 Frequency (Hz) 1 2 3 4 Method BEM, 40 Elements 0.95582 1.62832 2.05986 2.41284 Analytical Solution 0.95597 1.62777 2.05970 2.41198

The validation of numerical results using BEM with the deployed analytical approach is shown in Table 2. The results show a good convergence. For computations, the parameters of clamped Silicon membrane are considered as the following: $R=$ 0.5 m, $\rho =$ 998 kg/m3, thickness  0.001 m, material density ${\rho }_{m}=$ 2800kg/m3, Young modulus $E=$ 50 MPa, Poison’s ratio $v=$ 0.49, and $h=$ 1 m. Fundamental modes of sloshing of liquid vibrations are shown in Fig. 2.

Fig. 2. Fundamental slosh modes of liquid in a circular cylindrical container Fig. 3. Sloshing modes with membrane vibrations as Eva plastic roof on free surface To see the modes of membrane, Eva plastic membrane is used at the free surface. The radius of the membrane is taken with 0.5 m, thickness ${h}_{m}=$ 0.001 m, material density as 950 kg/m3. The sloshing modes of the membrane are shown in Fig. 3.

#### 4. Conclusions

In this paper, the sloshing in a right vertical circular cylindrical container in the presence of a membrane placed at the free surface is investigated. The fundamental sloshing modes of liquid with and without membrane at the free surface are reported. An analytical approach used to determine sloshing frequencies is validated using Boundary Element Method (BEM). A good agreement is shown between two, even for a small number of boundary elements. The main aim of this research was to validate both analytical and numerical methods using comparison of the results.

#### Acknowledgements

The authors acknowledge the support of the Indo-Ukraine Joint Research Project Grant funded by DST, India, and Education and Science Ministry of Ukraine.

1. Choudhary N., Bora S. N. Liquid sloshing in a circular cylindrical container containing a two-layer fluid. International Journal of Advances in Engineering Sciences and Applied Mathematics, Vol. 8, Issue 4, 2016, p. 240-248. [Publisher]
2. Choudhary N., Bora S. N. Linear sloshing frequencies in the annular region of a circular cylindrical container in presence of a rigid baffle. Sadhana-Academy Proceedings in Engineering Sciences, Springer, Vol. 42, Issue 5, 2017, p. 805-815. [Search CrossRef]
3. Falahaty H., Khayyer A., Gotoh H. Enhanced particle method with stress point integration for simulation of incompressible fluid-nonlinear elastic structure interaction. Journal of Fluids and Structures, Vol. 81, 2018, p. 325-360. [Publisher]
4. George A., Cho I. H. Anti sloshing effects of a vertical porous baffle in a rolling rectangular tank. Ocean Engineering, Vol. 214, 2020, p. 107871. [Publisher]
5. Gotoh H., Khayyer A. On the state of the art of particle methods for coastal and ocean engineering. Coastal Engineering Journal, Vol. 60, Issue 1, 2018, p. 79-103. [Publisher]
6. Goudarzi M. A., Danesh P. N. Numerical investigation of a vertically baffled rectangular tank under seismic excitation. Journal of Fluids and Structures, Vol. 61, 2018, p. 450-460. [Publisher]
7. Hosseini M., Goudarzi M. A., Soroor Reduction of seismic sloshing in floating roof liquid storage tanks by using a suspended annular baffle (SAB). Journal of Fluids and Structures, Vol. 71, 2017, p. 40-55. [Publisher]
8. Ibrahim R. A. Liquid sloshing dynamics: theory and applications. Cambridge University Press, p. 970, 2005. [Publisher]
9. Koh C. G., Luo M., Gao M., Bai W. Modelling of liquid sloshing with constrained floating baffle. Computers and Structures, Vol. 122, 2013, p. 270-279. [Publisher]
10. Molin B., Remy F. Experimental and Numerical Study of the sloshing motion in a rectangular tank with a perforated screen. Journal of Fluids and Structures, Vol. 43, 2013, p. 463-480. [Publisher]
11. Molin B., Remy F. Inertia effects in TLD sloshing with perforated screens. Journal of Fluids and Structures, Vol. 59, 2015, p. 165-177. [Publisher]
12. Sakai F., Inouse R. Some considerations on seismic design and controls of sloshing floating roofed oil tanks. The 14th World conference on Earthquake Engineering, 2008. [Search CrossRef]
13. Strelnikoiva E., Choudhary N., Denis V., Krivtchenko, Vasyl Gnitko I., Tonkonozhenko Anatoly M. Liquid vibrations in circular cylindrical tanks with and without baffles under horizontal and vertical excitations. Engineering Analysis with Boundary Elements, Vol. 120, 2020, p. 13-27. [Publisher]
14. Kolaei A., Rakheja S. Free vibration analysis of coupled sloshing-flexible membrane system in a liquid container. Journal of Vibration and Control, Vol. 25, Issue 1, 2019, p. 84-97. [Publisher]
15. Gnitko V., Degtyariov K., Karaiev A., Strelnikova E. Multi-domain boundary element method for axisymmetric problems in potential theory and linear isotropic elasticity. WIT Transactions on Engineering Sciences, Vol. 122, 2019, p. 13-25. [Search CrossRef]
16. Brebbia C. A., Telles J. C. F., And Wrobel L. C. Boundary Element Techniques. Springer-Verlag, Berlin and New York, 1984. [Publisher]

#### Cited By

 Vibroengineering PROCEDIA Narveen Kumar, Neelam Choudhary 2022 Neelam Choudhary 2021 Narveen Kumar, Neelam Choudhary 2021