# An Approximate Analytical Solution of Sloshing Frequencies for a Liquid in Various Shape Aqueducts

Authors: Yuchun Li; Zhuang Wang
Journal: Shock and Vibration
 (2014)

---

## Abstract

An approximate analytical solution of sloshing frequencies for a liquid in the various shape aqueducts is formulated by using the Ritz method. The present approximate method is, respectively, applied to find the sloshing frequencies of the liquid in rectangular, trapezoid, oval, circular, U-shaped tanks (aqueducts), and various shape tuned liquid dampers (TLD). The first three antisymmetric and symmetric frequencies by the present approach are within 5% accuracy compared to the other analytical, numerical, and experimental values. The approximate solutions of this paper for the various shape aqueducts are acceptable to the engineering applications.


---

## Body

## 1. Introduction

Sloshing is a common physical phenomenon which occurs in moving tanks with contained liquid masses, such as liquid cargo carriers, rockets, aircrafts, and the seismically excited storage tanks, dams, reactors, and nuclear vessels. The sloshing frequencies of contained liquid are essential in the analysis and design of the liquid tanks and the associated structures.The modal analysis of the linearized free liquid surface in a rigid container has been widely studied by many scholars [1–4]. The various methods of modal analysis were well documented in the books by Ibrahim [3] and Faltinsen and Timokha [4]. These approaches are, respectively, variational formulation (Lawrence et al. [5]), integral equation (Budiansky [6]), conformal mapping (Fox and Kuttler [7]), bipolar coordinate transformation (McIver [8]), and so forth. The analytical and approximate methods to estimating the sloshing frequencies were developed by Evans [9] and Evans and Linton [10]. The numerical methods [11] in common use are the finite difference, finite element, and boundary element. The finite element may be the most commonly used approach in the sloshing analyses. On the basis of the small displacement assumption, the irrotational fluid can be modeled by a degenerate solid finite element [12, 13] with no shear resistance. This displacement-based finite element formulation, which had been introduced into the commercial ANSYS Code [14], was usually used for simulating the sloshing problems [15, 16] in certain engineering fields.The sloshing frequencies of a liquid in an arbitrary-section aqueduct can be accurately extracted by the above mentioned methods. However, it is not easy work. For many engineering problems, the engineers only need the approximate previous order frequencies. The purpose of this paper is to provide a simple and quick method to evaluate the sloshing frequencies of a liquid in the various shape aqueducts.


## 2. Basic Equations of Two-Dimensional Free Sloshing

An arbitrary-section aqueduct is shown in Figure1, where the 2-D Cartesian coordinates oxz are fixed to the aqueduct and the symbols Ω, ∂Sf, and ∂Sw, respectively, denote the static liquid region, the free liquid surface, and the wetted boundary. The x axis is in the plane of free surface and y axis points vertically upwards through the midpoint of the free surface. The width of free surface is 2a. The aqueduct is presumed to be rigid.Figure 1
An arbitrary-section aqueduct.On the basis of the linear potential theory [1–4], the fluid in absence of surface tension on the free surface can be assumed to be inviscid, incompressible, irrotational, and of small displacement; the two-dimensional free sloshing of a liquid inside a motionless aqueduct can be described by the following boundary value problem:
(1)∇2Φ(x,z,t)=0,(2)∂Φ(x,z,t)∂n|∂Sw=0,(3)[∂Φ(x,z,t)∂z-∂h(x,z,t)∂t]|∂Sf=0,(4)[∂Φ(x,z,t)∂t+gh(x,z,t)]|∂Sf=0,
where Φ(x,z,t) is the function of velocity potential, h(x,z,t) is the wave-height function on the free surface, t is the time, g is the acceleration of gravity, n is the normal drawn outwardly on the wetted boundary (see Figure 1), and ∂/∂n denotes differentiation in the direction of the normal.


## 3. Approximate Analytical Solution

In order to find the solutions of (1)–(4), we may proceed directly by assuming trial solutions:
(5)Φ(x,z,t)=iωφ(x,z)exp⁡(iωt),h(x,z,t)=f(x,z)exp⁡(iωt),
where i=-1, ω is the sloshing (circular) frequency and φ(x,z) and f(x,z) are, respectively, the amplitude functions of Φ(x,z,t) and h(x,z,t). Substituting (5) into (1)–(4) and merging (3) and (4) into a single equation, we arrive at the following eigenvalue problem:
(6)∇2φ(x,z)=0,(x,z)∈Ω,(7)∂φ(x,z)∂n|∂Sw=0,(8)[∂φ(x,z)∂z-ω2gφ(x,z)]|∂Sf=0.
The sloshing frequency can also be expressed as [1, 3]
(9)ω2=g∬Ω[(∂φ(x,z)/∂x)2+(∂φ(x,z)/∂z)2]dxdz∫∂Sfφ2(x,z)ds.
From (8) and (9), we further obtain
(10)ω2=g(∂φ(x,z)/∂z)|z=0φ(x,z)|z=0,(11)∬Ω[(∂φ(x,z)/∂x)2+(∂φ(x,z)/∂z)2]dxdz∫-aaφ2(x,z)dx=(∂φ(x,z)/∂z)|z=0φ(x,z)|z=0.For thejth frequency ωj (j=1,2,3,…) and mode function φj(x,z) (j=1,2,3,…), they both satisfy (10) and (11). It can be seen from (10) that the sloshing frequency ωj exclusively depends on the mode function φj(x,z). In order to obtain an approximate frequency ωj, the function φj(x,z) can be approximately assumed to be a product of two individual functions ξj(x) and ηj(z) of the coordinates; that is, φj(x,z)=ξj(x)·ηj(z). In the physical sense, the function ξj(x) presents the jth mode shape of free liquid surface. On the basis of the Ritz method, ξj(x) can be approximately expressed as a known trigonometric function sinαjx for the antisymmetric mode or cos⁡βjx for the symmetric mode, in which αj, βj are the known constants. Inserting φj(x)≈sinαjx·ηj(z) or φj(x,z)≈cos⁡βjx·ηj(z) into (6), we can find the solution ηj(z)=Aj(cosh⁡αjz+Bjsinh⁡αjz) or ηj(z)=Cj(cosh⁡βjz+Djsinh⁡βjz), in which Aj, Bj, Cj, and Dj (j=1,2,3,…) are the undetermined coefficients. Because the coefficients Aj and Cj have no effect on the final outcomes of the frequencies, they are ignored here. Finally the mode function φj(x,z), corresponding to the jth antisymmetric and symmetric modes, can be proximately presented as
(12)φj(x,z)≈{sinαjx·(cosh⁡αjz+Bjsinh⁡αjz)iiiiiiiiiiiiiiiantisymmetricmodecos⁡βjx·(cosh⁡βjz+Djsinh⁡βjz)iiiiiiiiiiiiiiiiiiiiiisymmetricmode(j=1,2,3,…),
where
(13)αj=(2j-1)π2a,βj=jπa.
Thus, substitution of (12) into (10) yields the antisymmetric and symmetric sloshing frequencies
(14)ωj2≈{BjαjgantisymmetricDjβjgsymmetric(j=1,2,3,…).
It is noted that the function φj(x,z) has only one unknown coefficient Bj (or Dj), which can be uniquely determined by (11). Substituting (12) into (11), we, respectively, acquire the two quadratic equations with the unknowns Bj and Dj:
(15)aAjBj2+bAjBj+cAj=0,(j=1,2,3,…),aSjDj2+bSjDj+cSj=0,(j=1,2,3,…),
where
(16)aAj=∬Ω(sin2αjx+sinh2αjz)dxdz,bAj=2∬Ωsinh⁡αjz·cosh⁡αjz·dxdz-aαj,cAj=∬Ω(cos⁡2αjx+sinh2αjz)dxdz,iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii(j=1,2,3,…),aSj=∬Ω(cos⁡2βjx+sinh2βjz)dxdz,bSj=2∬Ωsinh⁡βjz·cosh⁡βjz·dxdz-aβj,cSj=∬Ω(sin2βjx+sinh2βjz)dxdz,iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii(j=1,2,3,…).
In a real physical system, there should be one root for (15). On the basis of an analysis in Section 4.1, where the sloshing frequencies of a rectangular tank are evaluated, the authentic roots of (15) can be determined as follows:
(17)Bj=-bAj-bAj2-4·aAj·cAj2aAj,(j=1,2,3,…),Dj=-bSj-bSj2-4·aSj·cSj2aSj,(j=1,2,3,…).
The rest of the roots of (15) are the unwanted solutions, which are neglected here. For an arbitrary-section aqueduct, the coefficients (aAj, bAj, cAj) and (aSj, bSj, cSj) can be obtained by using (16) through integrals in the static liquid domain Ω. Then, the coefficients Bj and Dj can be figured out by means of (17). Finally the antisymmetric and symmetric sloshing frequencies can be calculated by (14).


## 4. Numerical Examples

The analytical solution of this study is suitable for the fluid in an arbitrary-section aqueduct (tank). To validate this solution, the present method is applied to solve the sloshing frequencies of liquid in various shape tanks, and the frequency results of this study are compared with the other theoretical and experimental values.


### 4.1. Rectangular Tank

Figure2 shows a rectangular tank, in which 2a is the width of liquid and H is the depth of liquid. On the basis of (13)–(17), the sloshing frequencies of liquid inside a rectangular tank can be derived as
(18)ωj2={(2j-1)πg2atanh[(2j-1)π2aH]iiiiiiiiiiiiiiiiiiiiiiiiiiantisymmetricjπgatanh⁡(jπaH)iiiiiiiiiiisymmetric(j=1,2,3,…).Figure 2
A rectangular tank.It can be found that (18) is just the exact solution of sloshing frequencies. That is because the trigonometric functions in (12) are exactly the mode shape functions of the free fluid surface.In this example, solving (15) can yield an extraneous root, which is apparently not the authentic solution of the liquid system and is discarded here.


### 4.2. Trapezoid Tank

Gardarsson [17] investigated the sloshing characteristics for a sloping bottom tank (isosceles trapezoid tank). In his PhD thesis, Gardarsson proposed a finite difference method for calculating the sloshing frequencies of fluid in a trapezoid tank. The tank configuration used for the experimental and theoretical analyses is sketched in Figure 3, in which the water depth H, half-width of the free surface a, and slope angle θ of side walls are fully presented in Table 1.Table 1
Comparison of the fundamental frequencies by Gardarsson [17] and present method.
H (m)
a (m)
θ (degree)
Gardarsson [17] theory (Hz)
Gardarsson [17] experiment (Hz)
Present method (Hz)
0.04
0.1145
30°
1.2577
1.2812
1.2454
0.07
0.166
30°
1.0875
1.0625
1.0889
0.10
0.218
30°
0.9614
0.9375
0.9663Figure 3
A trapezoid tank.The present approximate method is applied to calculate the fundamental frequencies of the three trapezoid tanks. The results by Gardarsson and present approach are listed and compared in Table1. It can be seen from Table 1 that the results of this study are within 3% accuracy compared to the analytical and experimental values of Gardarsson. The solutions of this study have a good agreement with those by Gardarsson.


### 4.3. Semiellipse Tank

Hasheminejad and Aghabeigi [18] used a conformal transformation approach for studying the sloshing frequencies of fluid in half-full horizontal elliptical tanks. In the light of [18], the baffled and unbaffled ellipse container geometries are shown in Figure 4. A liquid fills the tanks to their half capacity, with the major and minor semiaxes and baffle length “a′,” “b′,” and “L,” respectively. The symbol “a” is the half width of still-liquid surface. For the un-baffled tank, the baffle length L is zero (i.e., a′=a).Figure 4
Baffled and unbaffled elliptical tanks.One baffled and one unbaffled elliptic tanks from [18] are selected as the numerical examples. The dimensions of the tanks are indicated in Table 2. The first three normalized antisymmetric and symmetric sloshing frequencies (Ωj=a′b′ωj2/g (j=1,2,3)) by present study and [18] are displayed and compared in Table 2.Table 2
The first three normalized antisymmetric and symmetric sloshing frequencies by [18] and present study Ωj=a′b′ωj2/g,  (j=1,2,3).
Unbaffled tank (a′=a=1.0m,b′=0.2m)
Baffled tank (a′=1.0m,b′=0.2m,L=0.05m)
Antisymmetric
Symmetric
Antisymmetric
Symmetric
Present study
[18]
Present study
[18]
Present study
[18]
Present study
[18]
Ω
1 0.19866
0.19542
0.67487
0.66059
1.46170
1.43313
3.22321
3.21824
Ω
2 1.34934
1.29511
2.12996
2.05458
4.83482
4.91081
6.44640
6.56741
Ω
3 2.94933
2.84342
3.76968
3.58803
8.05800
8.21077
9.66961
9.84664It can be observed from Table2 that the maximum relative error of the normalized frequencies between [18] and present study is approximately 5%; that is, the maximum relative error of the frequencies is around 2.5%. The approximate results of this paper agree well with the analytical solutions of [18].The solution of present paper not only applies to the semiellipse section but also to the ellipse tanks with arbitrary liquid depth.


### 4.4. Various Shape Tuned Liquid Dampers (TLD)

Idir et al. [19] developed an equivalent bottom-wet-length method for calculating the fundamental frequencies of the various shape tuned liquid dampers (TLD). The three TLD geometries from [19], as shown in Figure 5, are presented here as the illustrative examples. The first antisymmetric (fundamental) sloshing frequencies by present approach and [19] are exhibited and compared in Table 3. Via comparison, the maximum relative error of the fundamental frequencies between [19] and present study is less than 3%. The two results match well.Table 3
Comparison of the fundamental frequencies by Idir et al. [19] and present study.
Bottom shape
Idir et al. [19] theory (Hz)
Idir et al. [19] experiment (Hz)
Present study (Hz)
Flat
1.0270
1.0050
1.0257
V-shaped
1.0270
1.0310
1.0279
Arc-shaped
1.0270
1.0550
1.0558Various shape tuned liquid dampers (length unit: mm).
(a)
Flat bottom
(b)
V-shaped bottom
(c)
Arc-shaped bottom


### 4.5. Circular Aqueduct

Figure6 shows a circular-section aqueduct, in which D is the inner diameter (R is the inner radius), H is the liquid depth, and a is the half width of the free liquid surface. Both the finite element and present approaches are used for computing the sloshing frequencies. The ANSYS Code [14] provides a fluid element (Fluid80) for simulating the sloshing problem of contained fluid. In this example, the parameters of fluid are chosen as mass density ρ=1000 (kg/m3), bulk modulus Ev=2.067×109 (Pa), and radius of aqueduct R=3.0 (m). The viscosity coefficient of the fluid is assumed as zero (without damping). The interface between liquid and the rigid aqueduct is treated as a slip boundary condition, such that the normal relative displacement of the fluid and rigid container on the interface is forced to be zero and the tangential relative displacement on the interface is not restrained. Many obtained results using ANSYS Code have been verified by the previous investigations [15, 16]. For example, in [16] the fundamental frequencies in the rectangular tanks by ANSYS Code and the exact solution [3] were compared with each other, and the results showed the two solutions have a good agreement.Figure 6
A circular-section aqueduct.When the liquid depthH varies from 0.35D to 0.70D, the corresponding sloshing frequencies are both extracted by the finite element and present methods. Figures 7(a) and 7(b) show the variation of ω(R/g)1/2 with liquid-depth ratio H/D. It can be drawn from Figure 7 that the maximum relative error between present study and finite element is not greater than 5.0%. The two results match well.Variationω(R/g)1/2 with liquid-depth ratio H/D (circular aqueduct).
(a)
Antisymmetric frequency
(b)
Symmetric frequency


### 4.6. U-Shaped Aqueduct

Figure8 shows a U-shaped aqueduct, in which R is the inner radius of the bottom semicircle and H (H=h+R) is the liquid depth. The computing methods and parameters of the fluid are the same as in Section 4.5. The liquid depth above the bottom semicircle is fixed as h=1.5 m. When the liquid-depth ratio h/R varies from 0.1 to 1.0, the corresponding sloshing frequencies can be acquired by ANSYS Code and present method. Figures 9(a) and 9(b) show the variation of ω(R/g)1/2 with liquid-depth ratio h/R.Figure 8
A U-shaped section aqueduct.Variation ofω(R/g)1/2 with liquid-depth ratio h/R (U-shaped aqueduct).
(a)
Antisymmetric frequency
(b)
Symmetric frequencyIt can be observed from Figure9 that the results of this study have a good consistency with those by ANSYS Code (finite element). The maximum error between present study and finite element is less than 5.0%.It can be noted from Figures7 and 9 that the results of present study are slightly greater than those of finite element (ANSYS Code). This may be attributed to the compressibility of the fluid. A compressible fluid was applied in ANSYS Code, while an incompressible liquid was used in present study. The compressible case can produce lower natural frequencies. According to [16], the sloshing frequencies in a rectangular tank by ANSYS Code are slightly lower than the exact solutions using (18).


## 4.1. Rectangular Tank

Figure2 shows a rectangular tank, in which 2a is the width of liquid and H is the depth of liquid. On the basis of (13)–(17), the sloshing frequencies of liquid inside a rectangular tank can be derived as
(18)ωj2={(2j-1)πg2atanh[(2j-1)π2aH]iiiiiiiiiiiiiiiiiiiiiiiiiiantisymmetricjπgatanh⁡(jπaH)iiiiiiiiiiisymmetric(j=1,2,3,…).Figure 2
A rectangular tank.It can be found that (18) is just the exact solution of sloshing frequencies. That is because the trigonometric functions in (12) are exactly the mode shape functions of the free fluid surface.In this example, solving (15) can yield an extraneous root, which is apparently not the authentic solution of the liquid system and is discarded here.


## 4.2. Trapezoid Tank

Gardarsson [17] investigated the sloshing characteristics for a sloping bottom tank (isosceles trapezoid tank). In his PhD thesis, Gardarsson proposed a finite difference method for calculating the sloshing frequencies of fluid in a trapezoid tank. The tank configuration used for the experimental and theoretical analyses is sketched in Figure 3, in which the water depth H, half-width of the free surface a, and slope angle θ of side walls are fully presented in Table 1.Table 1
Comparison of the fundamental frequencies by Gardarsson [17] and present method.
H (m)
a (m)
θ (degree)
Gardarsson [17] theory (Hz)
Gardarsson [17] experiment (Hz)
Present method (Hz)
0.04
0.1145
30°
1.2577
1.2812
1.2454
0.07
0.166
30°
1.0875
1.0625
1.0889
0.10
0.218
30°
0.9614
0.9375
0.9663Figure 3
A trapezoid tank.The present approximate method is applied to calculate the fundamental frequencies of the three trapezoid tanks. The results by Gardarsson and present approach are listed and compared in Table1. It can be seen from Table 1 that the results of this study are within 3% accuracy compared to the analytical and experimental values of Gardarsson. The solutions of this study have a good agreement with those by Gardarsson.


## 4.3. Semiellipse Tank

Hasheminejad and Aghabeigi [18] used a conformal transformation approach for studying the sloshing frequencies of fluid in half-full horizontal elliptical tanks. In the light of [18], the baffled and unbaffled ellipse container geometries are shown in Figure 4. A liquid fills the tanks to their half capacity, with the major and minor semiaxes and baffle length “a′,” “b′,” and “L,” respectively. The symbol “a” is the half width of still-liquid surface. For the un-baffled tank, the baffle length L is zero (i.e., a′=a).Figure 4
Baffled and unbaffled elliptical tanks.One baffled and one unbaffled elliptic tanks from [18] are selected as the numerical examples. The dimensions of the tanks are indicated in Table 2. The first three normalized antisymmetric and symmetric sloshing frequencies (Ωj=a′b′ωj2/g (j=1,2,3)) by present study and [18] are displayed and compared in Table 2.Table 2
The first three normalized antisymmetric and symmetric sloshing frequencies by [18] and present study Ωj=a′b′ωj2/g,  (j=1,2,3).
Unbaffled tank (a′=a=1.0m,b′=0.2m)
Baffled tank (a′=1.0m,b′=0.2m,L=0.05m)
Antisymmetric
Symmetric
Antisymmetric
Symmetric
Present study
[18]
Present study
[18]
Present study
[18]
Present study
[18]
Ω
1 0.19866
0.19542
0.67487
0.66059
1.46170
1.43313
3.22321
3.21824
Ω
2 1.34934
1.29511
2.12996
2.05458
4.83482
4.91081
6.44640
6.56741
Ω
3 2.94933
2.84342
3.76968
3.58803
8.05800
8.21077
9.66961
9.84664It can be observed from Table2 that the maximum relative error of the normalized frequencies between [18] and present study is approximately 5%; that is, the maximum relative error of the frequencies is around 2.5%. The approximate results of this paper agree well with the analytical solutions of [18].The solution of present paper not only applies to the semiellipse section but also to the ellipse tanks with arbitrary liquid depth.


## 4.4. Various Shape Tuned Liquid Dampers (TLD)

Idir et al. [19] developed an equivalent bottom-wet-length method for calculating the fundamental frequencies of the various shape tuned liquid dampers (TLD). The three TLD geometries from [19], as shown in Figure 5, are presented here as the illustrative examples. The first antisymmetric (fundamental) sloshing frequencies by present approach and [19] are exhibited and compared in Table 3. Via comparison, the maximum relative error of the fundamental frequencies between [19] and present study is less than 3%. The two results match well.Table 3
Comparison of the fundamental frequencies by Idir et al. [19] and present study.
Bottom shape
Idir et al. [19] theory (Hz)
Idir et al. [19] experiment (Hz)
Present study (Hz)
Flat
1.0270
1.0050
1.0257
V-shaped
1.0270
1.0310
1.0279
Arc-shaped
1.0270
1.0550
1.0558Various shape tuned liquid dampers (length unit: mm).
(a)
Flat bottom
(b)
V-shaped bottom
(c)
Arc-shaped bottom


## 4.5. Circular Aqueduct

Figure6 shows a circular-section aqueduct, in which D is the inner diameter (R is the inner radius), H is the liquid depth, and a is the half width of the free liquid surface. Both the finite element and present approaches are used for computing the sloshing frequencies. The ANSYS Code [14] provides a fluid element (Fluid80) for simulating the sloshing problem of contained fluid. In this example, the parameters of fluid are chosen as mass density ρ=1000 (kg/m3), bulk modulus Ev=2.067×109 (Pa), and radius of aqueduct R=3.0 (m). The viscosity coefficient of the fluid is assumed as zero (without damping). The interface between liquid and the rigid aqueduct is treated as a slip boundary condition, such that the normal relative displacement of the fluid and rigid container on the interface is forced to be zero and the tangential relative displacement on the interface is not restrained. Many obtained results using ANSYS Code have been verified by the previous investigations [15, 16]. For example, in [16] the fundamental frequencies in the rectangular tanks by ANSYS Code and the exact solution [3] were compared with each other, and the results showed the two solutions have a good agreement.Figure 6
A circular-section aqueduct.When the liquid depthH varies from 0.35D to 0.70D, the corresponding sloshing frequencies are both extracted by the finite element and present methods. Figures 7(a) and 7(b) show the variation of ω(R/g)1/2 with liquid-depth ratio H/D. It can be drawn from Figure 7 that the maximum relative error between present study and finite element is not greater than 5.0%. The two results match well.Variationω(R/g)1/2 with liquid-depth ratio H/D (circular aqueduct).
(a)
Antisymmetric frequency
(b)
Symmetric frequency


## 4.6. U-Shaped Aqueduct

Figure8 shows a U-shaped aqueduct, in which R is the inner radius of the bottom semicircle and H (H=h+R) is the liquid depth. The computing methods and parameters of the fluid are the same as in Section 4.5. The liquid depth above the bottom semicircle is fixed as h=1.5 m. When the liquid-depth ratio h/R varies from 0.1 to 1.0, the corresponding sloshing frequencies can be acquired by ANSYS Code and present method. Figures 9(a) and 9(b) show the variation of ω(R/g)1/2 with liquid-depth ratio h/R.Figure 8
A U-shaped section aqueduct.Variation ofω(R/g)1/2 with liquid-depth ratio h/R (U-shaped aqueduct).
(a)
Antisymmetric frequency
(b)
Symmetric frequencyIt can be observed from Figure9 that the results of this study have a good consistency with those by ANSYS Code (finite element). The maximum error between present study and finite element is less than 5.0%.It can be noted from Figures7 and 9 that the results of present study are slightly greater than those of finite element (ANSYS Code). This may be attributed to the compressibility of the fluid. A compressible fluid was applied in ANSYS Code, while an incompressible liquid was used in present study. The compressible case can produce lower natural frequencies. According to [16], the sloshing frequencies in a rectangular tank by ANSYS Code are slightly lower than the exact solutions using (18).


## 5. Conclusion

An approximate analytical solution of sloshing frequencies was formulated for a liquid in the arbitrary-section aqueducts. This approximate method has been, respectively, applied for finding the sloshing frequencies of the liquid in the various shape tanks (aqueducts). The first three antisymmetric and symmetric frequencies by this study in the various shape tanks are within 5% accuracy compared to the other analytical, numerical, and experimental values. The present approximate solutions for the various shape aqueducts are acceptable to the engineering applications. This method can be a simple and convenient way to estimate the sloshing frequencies of the liquid in various shape aqueducts for practical usage.



---

*Source: 672648-2014-02-13.xml*