# A Boundary Element Investigation of Liquid Sloshing in Coupled Horizontal and Vertical Excitation

Authors: De-Zhi Ning; Wei-Hua Song; Yu-Long Liu; Bin Teng
Journal: Journal of Applied Mathematics
 (2012)

---

## Abstract

Sloshing flows in a two-dimensional rigid rectangular tank under specified excitations in the coupled horizontal and vertical modes are simulated by using a higher-order boundary element method (BEM). The liquid sloshing is formulated as an initial-boundary-value problem based on the fully nonlinear potential flow theory. And a semi-mixed Eulerian-Lagrangian technique combined with the 4th-order Runge-Kutta scheme is employed to advance the solutions in the time marching process. A smoothing technique is applied to the free surface at every several time steps to avoid the possible numerical instabilities. Numerical results obtained are compared with the available solutions to validate the developed model. The parametric studies are carried out to show the liquid sloshing effects in terms of the slosh frequencies and excitation amplitudes in horizontal and vertical modes, the second-order resonance frequency, a bottom-mounted vertical rigid baffle, free surface displacement, and hydrodynamic forces acting on the tank.


---

## Body

## 1. Introduction

Liquid sloshing, displaying a free-surface fluctuation in partially filled containers under external excitations, is a known physical phenomenon existing in a variety of engineering applications such as liquid oscillation in large storage tank by earthquake, the motion of liquid fuel in aircraft and spacecraft, the liquid motion in tank, and the water flow on the deck of ship. Such liquid motions may lead to the unexpected instability and failure of engineering structural system. Especially when the excitation stroke is large or the excitation frequency is close to the natural frequency of the container, the caused damage may be a tremendous loss of human, economic, and environmental resources. Thus, it is still urged the need for further research on the understanding of the complex sloshing behavior and the technique for sloshing suppression.Numerous studies related to liquid sloshing have been carried out numerically, theoretically, and experimentally in the past several decades. For example, Faltinsen [1] derived the linear analytical solution of a horizontally excitation tank in 2D. Frandsen [2] obtained the linear and second-order analytical solutions of a fixed tank with initial free surface being cosine curve. Faltinsen et al. [3] performed sloshing experiments in a square-base tank. Akyildiz and Erdem Ünal [4] made sloshing experiments in a moving tank with different baffle arrangements. Okamoto and Kawahara [5] performed the experiments of liquid sloshing in rectangular and hexagonal tanks. Faltinsen [6] and Nakayama and Washizu [7] adopted the boundary element method to model large amplitude sloshing in a 2D rectangular container subjected to a horizontal excitation. Wu et al. [8] gave an account of both 2D and 3D sloshing problems based on finite element method. Frandsen and Borthwick [9] simulated the sloshing motion in fixed and vertically excited containers using a finite difference solver. Wang and Khoo [10], and Sriram et al. [11] numerically simulated 2D sloshing waves due to random excitation. Eswaran et al. [12] numerically and experimentally studied the sloshing motions in baffled and unbaffled tank subjected to horizontal excitation. Cho et al. [13] investigated the resonant sloshing response in 2D baffled tank by using FEM. Biswal et al. [14] used FEM to compute the nonlinear sloshing response of liquid in a rigid rectangular and cylindrical tank with horizontal baffle. Panigrahy et al. [15] experimentally researched hydrodynamic pressure developed on the tank with horizontal and vertical baffles subjected to horizontal excitation. Liu and Lin [16] adopted VOF method to simulate liquid sloshing in 3D rectangular baffled tank. Belakroum et al. [17] studied the vibratory behaviors of three different configurations of tanks equipped with baffles to predict the damping effect of baffles on sloshing in tank. In addition, Ortiz et al. [18–21] founded a powerful closed-form model to study the coupling of nonlinear liquid sloshing and the container motion, in which the fluid is modeled using a 2D BEM based on potential flow theory with modified Rayleigh damping, the container boundary can be arbitrarily curved, and the container motion can be induced by a flexible multibody system.As an extension of the previous studies, the liquid sloshing in a 2D tank with and without a vertical baffle, subjected to the coupled horizontal and vertical excitations, are numerically investigated by a fully nonlinear higher-order boundary element method. The effects of two excitations and bottom-mounted vertical baffle on the nonlinear liquid sloshing are further researched. The relating studies are still limited. This paper is organized as follows. The mathematical formulation is firstly given in Section2. Then, in Section 3, the present numerical results including wave surface and hydrodynamic pressure are compared with analytical solutions, published experimental data, and numerical results to validate the proposed model. In Section 4, liquid sloshing in a 2D tank without baffle subjected to the coupled horizontal and vertical excitations is firstly considered. The effect of excitation frequencies and intensities on liquid sloshing is studied. Next, liquid sloshing in baffled tank is considered to analyze the effect of baffle position and dimension on wave surface and wave forces. Finally, conclusions are given in Section 5.


## 2. Mathematical Formulation

Referring to Figure1, liquid sloshing in a rectangular tank with length L and water depth H is considered in the present study. Two Cartesian coordinate systems, including a space-fixed coordinate O0X0Z0 and a tank-fixed coordinate OXZ, are defined with Z and Z0 positive upwards. The two systems coincide with each other when the tank is at rest.Figure 1
Definition sketch.The displacement of tank is defined asXb=[xb(t),zb(t)] in coordinate O0X0Z0. Under the assumptions that fluid motion is irrotational and fluid is incompressible and inviscid, the fluid motion is, therefore, governed by Laplace equation:
(2.1)∇2ϕ=0,
where ϕ is the velocity potential. The following condition is satisfied on the interface of tank and liquid:
(2.2)∂ϕ∂n=U⋅n,
where U=dXb/dt is the velocity of the tank and n is the outward vector normal to the tank walls. On the instantaneous free surface z0=η0(x0,t), both the fully nonlinear dynamic and kinematic boundary conditions are satisfied:
(2.3)∂ϕ∂t+12∇ϕ⋅∇ϕ+gη0=0,∂η0∂t+∂ϕ∂x0∂η0∂x0-∂ϕ∂z0=0,
where η0 is free surface and g is gravity acceleration.By using the following (2.4) and (2.5), the (2.3) can be transformed into (2.6) and (2.7) in the moving system OXZ:
(2.4)∇x0z0=∇xz,(2.5)(∂∂t)x0z0=(∂∂t-dXbdt⋅∇)xz,(2.6)∂ϕ∂t-∇ϕ⋅dXbdt+12∇ϕ⋅∇ϕ+g(η+zb)=0,(2.7)∂η∂t+(∂ϕ∂x-dxbdt)∂η∂x-(∂ϕ∂z-dzbdt)=0,
where η=η0-zb is the free surface in coordinate system OXZ for fixed x. A semi-mixed Eulerian-Lagrangian technique is adopted to change the potential with time on the free surface as follows:
(2.8)δϕ[x,η(x,t),t]δt=∂ϕ∂t+∂ϕ∂z∂η∂t.
Thus dynamic boundary condition is written as
(2.9)δϕδt+∂ϕ∂z∂η∂t-∇ϕ⋅dXbdt+12∇ϕ⋅∇ϕ+g(η+zb)=0.
The velocity potential ϕ can be split as follows in the space-fixed system:
(2.10)ϕ=φ+xu+zv,
where u and v are the components of U at x and z directions, respectively. Substituting (2.10) into (2.1), (2.2), (2.7), and (2.9), we can obtain
(2.11)∇2φ=0,∂φ∂n=0,∂η∂t=-∂φ∂x∂η∂x+∂φ∂z,δφδt=∂φ∂z∂η∂t-12∇φ·∇φ-gη-xdudt-ηdvdt.
The initial condition ϕ=0 is written as
(2.12)φ(x,z,0)=-xu(0)-zv(0).
The initial free surface η(x,0) is fixed according to practical simulation.In addition, hydrodynamic pressure acting on tank can be obtained by Bernoulli equation in space-fixed coordinateO0X0Z0(2.13)P=-ρ(∂ϕ∂t+12∇ϕ⋅∇ϕ+gη).Through (2.4), (2.5) and (2.10), (2.13) is transformed in tank-fixed coordinate OXZ(2.14)P=-ρ(∂φ∂t+12∇φ⋅∇φ+gη+x∂u∂t+z∂v∂t).
Then, the force acting on tank can be obtained by integration and the time derivative of velocity potential, that is, first term in (2.14), is obtained by the acceleration potential method.The direct boundary integral equation is derived to solve the prescribed boundary value problem by using the second Green’s theorem [22]. Applying the Green function satisfying the impermeable condition on the sea bed, integral equation can be written as
(2.15)C(p)φ(p)=∫Γ(φ(p)∂G(p,q)∂n-G(p,q)∂φ(q)∂ndΓ),
where p and q are source and field points, and C is the solid angle which can be conveniently and economically computed by the indirect method in the present study [23]. Γ is liquid domain boundary including free surface boundary and solid boundary and G is a simple Green function and is written as follows:
(2.16)G(p,q)=12π(ln⁡r1+ln⁡r2),
where r1=(x-x0)2+(z-z0)2, r2=(x-x0)2+(z+z0+2H)2, and H is water depth.Then the boundary surface is discretized with a number of three-node line elements. The geometry of each element is represented by the quadratic shape functions, thus the entire curved boundary can be approximated by a number of higher-order elements. Within the boundary elements, physical variables are also interpolated by the same shape functions, that is, the elements are isoparametric. Then the discretized integral equation can be written as linear equation systems:(2.17)[A11A12A21A22]{φ|on⁡solidsurfaces∂φ∂n|on⁡freesurface}={B1B2}.The discretized integral equation is reformulated at each time step, as the nodes are moved vertically and the grid is updated by stretching. For the sake of better stability and convenience, the 4th-order Runge-Kutta scheme is used at each time step. The methodology for obtaining the velocity components of the nodes, the updated position of the free surface, and the updated velocity potential is the same as described by Ning and Teng [24]. As the proposed container is a rectangular one, it is quite easy to use the five-point smoothing technique on the instantaneous free surface at every several time steps to prevent the so-called saw-tooth instabilities when some steeper waves are simulated. The detail can also be referred to Ning et al. [25]. For the container with curved boundaries, it is advised to use the smoothing and volume-correction techniques under polar coordinate system described by Ortiz and Barhorst [19].


## 3. Numerical Validations

In this section, the present model will be validated through comparing the computed wave surface and hydrodynamic pressure with analytical solution, published experimental data, and numerical results. A tank with and without baffle will be considered, respectively. The natural modes of a rectangular tank can be obtained by solving a linearized natural sloshing standing wave problem. And the natural frequency can be expressed as(3.1)ωn=gnπLtanh(nπLH),
where n=1,2,3 and so on.


### 3.1. Validation of Wave Surface

Simulation of liquid sloshing in a rectangular tank without baffle is the first benchmark validation test. A tank subjected to pure horizontal excitationxb(t)=Asin(ωt) is considered. The dimensions of the tank are taken as L=2.0 m and H=1.0 m. The initial wave surface is supposed to be stable. Through numerical convergence test, the temporal and spatial steps are selected as Δt=0.01 s and Δx=Δz=0.05 m. Figure 2 presents the distributions of wave profile measured from the bottom of the tank at two dimensionless time t/(h/G)0.5=13.0667 and t/(h/G)0.5=15.725 for A=0.0186 m and ω=0.999ω1. The comparisons among experimental data [26], results from finite element method [8], and present numerical results are also shown. It can be seen that the numerical simulation gives good agreement with experimental data and FEM results.Distribution of free surface elevation at two times.
(a)
t
/
(
h
/
G
)
0.5
=
13.0667
(b)
t
/
(
h
/
G
)
0.5
=
15.725The second set of validation tests is concerned with a tank subjected to pure vertical excitation. The tank scale isL/H=2.0. The vertical excitation is zb(t)=Avcos(ωvt), where Av   and ωv   are the amplitude and frequency of vertical excitation, respectively. The initial wave surface is taken as η(x,0)=acos(πx/L). Parameter kv=Avωv2/g is used as a measure of the intensity of vertical excitation and  ε=aω12/g is defined as a measure of wave profile nonlinearity. Figure 3 shows the comparison of wave surface time history at the left tank wall between the present and finite difference [2] numerical results in the condition of ε=0.0014 (i.e., a=0.00049 m), ωv=0.8ω1,  and  kv=0.5 (i.e., Av=0.27 m). It is found that present results have good agreement with Frandsen’s [2] result.Figure 3
Time history of wave surface at the left tank wall due to vertical excitation.Further, liquid sloshing induced by coupled horizontal and vertical excitation is studied. The dimension of the tank is the same as the above case. Horizontal excitation and vertical excitation are taken asxb(t)=Ahcos(ωht) and zb(t)=Avcos(ωvt). Parameter kh=Ahωh2/g is used as a measure of the intensity of horizontal excitation. In this condition, initial wave surface can be taken as η(x,0)=0. Figure 4 shows the comparison of time histories of wave surface at the left tank wall between present and Frandsen’s [2] result in the condition of ωh=0.98ω1, ωv=0.8ω1, kh=0.0014, and kv=0.5 (i.e., Ah=0.0005 m, Av=0.27 m). It can be seen that there still is a good agreement between two methods.Figure 4
Time history of wave surface at the left tank wall due to coupled horizontal and vertical excitation.Another case for liquid forced sloshing in a tank with a vertical baffle is investigated. The tank (L/H=2.0) is assumed to move as X(t)=asinωt. Excitation amplitude and frequency is a=0.002 m and ω=5.29 rad/s (ω is near first-order natural frequency of tank without baffle). A vertical baffle with Hb=0.75H and B=0.01 m is arranged at the center of tank bottom. Figure 5 shows the comparison of the time histories of wave surface at the right wall and the comparison between the present method and VOF method [16]. From the figure, it can be seen that the results of both methods are in good agreement and wave surface does not resonate due to the vertical baffle. Baffle changes the natural frequency of the tank and makes present natural frequency away from external excitation frequency to restrain liquid sloshing. Figure 6 shows the comparison of the first-order natural frequency of tank with different baffle length between present result and the linear frequency-domain result [27]. It is found that they coincide with each other well and the longer the baffle is, the smaller the natural frequency becomes, that is, far away from that of the unbaffled tank.Figure 5
Time history of wave surface at the right tank wall with a vertical baffle.Figure 6
Distribution of the first-order natural frequency for different baffle length.


### 3.2. Validation of Hydrodynamics

Firstly, a tank without baffle (L/H=2.0) is forced to move as X(t)=asinωt and excitation amplitude and frequency is a=0.00186 m and ω=5.3 rad/s. Because movement amplitude a is very small, the liquid sloshing can be considered as a linear problem. Figure 7 shows the comparisons of hydrodynamic pressures between the present results and linear analytical solutions [8] at free surface and tank bottom of the left tank wall. It can be found that the present results and linear analytical solution have a good agreement and the hydrodynamic pressure trends to be resonant.Time histories of hydrodynamic pressure at the free surface and tank bottom of the left tank wall.
(a)
at the free surface(-L/2.0,η)
(b)
at the tank bottom(-L/2.0,-H)A tank withL=30 m and depth H=15 m, moving as X(t)=-acosωt (a=0.5 m, ω=0.628 rad/s), is further considered. The tank motion is stopped after 5 s and the inner liquid begins to slosh freely. Figure 8 shows the time histories of horizontal wave force acting on the tank with and without vertical baffle. The vertical baffle with Hb=9.0 m is arranged at center of tank bottom. In the two figures, it is obvious that the present results and those obtained from FEM [28] coincide with each other very well and the wave force is lower significantly due to the effect of the vertical baffle. It also can be seen from the figure that both the horizontal wave loads have a sudden decrease at t=5 s due to the tank motion stopped at the moment.total horizontal hydrodynamic pressure time history with and without vertical baffle.
(a)
without baffle
(b)
with baffleThrough the above examples, the present numerical model is verified to be accurate to simulate liquid sloshing due to coupled horizontal and vertical excitations in a tank with and without vertical baffle.


## 3.1. Validation of Wave Surface

Simulation of liquid sloshing in a rectangular tank without baffle is the first benchmark validation test. A tank subjected to pure horizontal excitationxb(t)=Asin(ωt) is considered. The dimensions of the tank are taken as L=2.0 m and H=1.0 m. The initial wave surface is supposed to be stable. Through numerical convergence test, the temporal and spatial steps are selected as Δt=0.01 s and Δx=Δz=0.05 m. Figure 2 presents the distributions of wave profile measured from the bottom of the tank at two dimensionless time t/(h/G)0.5=13.0667 and t/(h/G)0.5=15.725 for A=0.0186 m and ω=0.999ω1. The comparisons among experimental data [26], results from finite element method [8], and present numerical results are also shown. It can be seen that the numerical simulation gives good agreement with experimental data and FEM results.Distribution of free surface elevation at two times.
(a)
t
/
(
h
/
G
)
0.5
=
13.0667
(b)
t
/
(
h
/
G
)
0.5
=
15.725The second set of validation tests is concerned with a tank subjected to pure vertical excitation. The tank scale isL/H=2.0. The vertical excitation is zb(t)=Avcos(ωvt), where Av   and ωv   are the amplitude and frequency of vertical excitation, respectively. The initial wave surface is taken as η(x,0)=acos(πx/L). Parameter kv=Avωv2/g is used as a measure of the intensity of vertical excitation and  ε=aω12/g is defined as a measure of wave profile nonlinearity. Figure 3 shows the comparison of wave surface time history at the left tank wall between the present and finite difference [2] numerical results in the condition of ε=0.0014 (i.e., a=0.00049 m), ωv=0.8ω1,  and  kv=0.5 (i.e., Av=0.27 m). It is found that present results have good agreement with Frandsen’s [2] result.Figure 3
Time history of wave surface at the left tank wall due to vertical excitation.Further, liquid sloshing induced by coupled horizontal and vertical excitation is studied. The dimension of the tank is the same as the above case. Horizontal excitation and vertical excitation are taken asxb(t)=Ahcos(ωht) and zb(t)=Avcos(ωvt). Parameter kh=Ahωh2/g is used as a measure of the intensity of horizontal excitation. In this condition, initial wave surface can be taken as η(x,0)=0. Figure 4 shows the comparison of time histories of wave surface at the left tank wall between present and Frandsen’s [2] result in the condition of ωh=0.98ω1, ωv=0.8ω1, kh=0.0014, and kv=0.5 (i.e., Ah=0.0005 m, Av=0.27 m). It can be seen that there still is a good agreement between two methods.Figure 4
Time history of wave surface at the left tank wall due to coupled horizontal and vertical excitation.Another case for liquid forced sloshing in a tank with a vertical baffle is investigated. The tank (L/H=2.0) is assumed to move as X(t)=asinωt. Excitation amplitude and frequency is a=0.002 m and ω=5.29 rad/s (ω is near first-order natural frequency of tank without baffle). A vertical baffle with Hb=0.75H and B=0.01 m is arranged at the center of tank bottom. Figure 5 shows the comparison of the time histories of wave surface at the right wall and the comparison between the present method and VOF method [16]. From the figure, it can be seen that the results of both methods are in good agreement and wave surface does not resonate due to the vertical baffle. Baffle changes the natural frequency of the tank and makes present natural frequency away from external excitation frequency to restrain liquid sloshing. Figure 6 shows the comparison of the first-order natural frequency of tank with different baffle length between present result and the linear frequency-domain result [27]. It is found that they coincide with each other well and the longer the baffle is, the smaller the natural frequency becomes, that is, far away from that of the unbaffled tank.Figure 5
Time history of wave surface at the right tank wall with a vertical baffle.Figure 6
Distribution of the first-order natural frequency for different baffle length.


## 3.2. Validation of Hydrodynamics

Firstly, a tank without baffle (L/H=2.0) is forced to move as X(t)=asinωt and excitation amplitude and frequency is a=0.00186 m and ω=5.3 rad/s. Because movement amplitude a is very small, the liquid sloshing can be considered as a linear problem. Figure 7 shows the comparisons of hydrodynamic pressures between the present results and linear analytical solutions [8] at free surface and tank bottom of the left tank wall. It can be found that the present results and linear analytical solution have a good agreement and the hydrodynamic pressure trends to be resonant.Time histories of hydrodynamic pressure at the free surface and tank bottom of the left tank wall.
(a)
at the free surface(-L/2.0,η)
(b)
at the tank bottom(-L/2.0,-H)A tank withL=30 m and depth H=15 m, moving as X(t)=-acosωt (a=0.5 m, ω=0.628 rad/s), is further considered. The tank motion is stopped after 5 s and the inner liquid begins to slosh freely. Figure 8 shows the time histories of horizontal wave force acting on the tank with and without vertical baffle. The vertical baffle with Hb=9.0 m is arranged at center of tank bottom. In the two figures, it is obvious that the present results and those obtained from FEM [28] coincide with each other very well and the wave force is lower significantly due to the effect of the vertical baffle. It also can be seen from the figure that both the horizontal wave loads have a sudden decrease at t=5 s due to the tank motion stopped at the moment.total horizontal hydrodynamic pressure time history with and without vertical baffle.
(a)
without baffle
(b)
with baffleThrough the above examples, the present numerical model is verified to be accurate to simulate liquid sloshing due to coupled horizontal and vertical excitations in a tank with and without vertical baffle.


## 4. Results and Discussion

In this section, liquid sloshing subjected to coupled horizontal and vertical excitation, with and without a vertical baffle, will be systematically investigated. The width and liquid depth in the tank are defined asL=1.0 m and H=0.5 m, respectively. The horizontal and vertical excitations are taken as xb(t)=Ahcos(ωht) and  zb(t)=Avcos(ωvt), respectively. Ah and Av are excitation amplitude. ωh and ωv are excitation frequency. Parameters kh=Ahωh2/g and kv=Avωv2/g are used as a measure of the intensity of the horizontal and vertical excitations, respectively. The initial wave surface is always supposed to be zero.


### 4.1. Effects of Horizontal and Vertical Excitation Frequencies on Liquid Sloshing

As it is well known, liquid sloshing is most violent when external frequency approaches to natural frequency. Therefore, three cases are considered for the external excitation frequencies separately or both equal to natural frequency. Case A is that parametersωh=ωv=ω1,  kh=0.003, and kv=0.2 (i.e., Ah=0.001 m and  Av=0.069 m) are taken. Case B is that ωh=ω1,  ωv=0.7ω1, kh=0.003, and kv=0.2 (i.e., Ah=0.001 m and Av=0.142 m) are defined. Case C is that ωh=0.7ω1, ωv=ω1,  kh=0.0014, and kv=0.2 (i.e., Ah=0.001 m and Av=0.069 m) are chosen. Through numerical convergence test, temporal step is selected as Δt=0.012 s and spatial step is defined as Δx=Δz=0.05 m. Figure 9 shows time histories of wave surface at the left tank wall and the corresponding spectral charts obtained by FFT method. From the figure, It can be seen that the resonance phenomena occur as shown in Figures 9(a) and 9(c) when the horizontal excitation frequency is equal to natural frequency. Horizontal excitation plays a dominated role on the shape and peak amplitude of wave surface and can be always found in the spectral charts. However, the vertical excitation frequency does not exist alone in spectral chart and occurs as the forms of sum- and difference-frequency with horizontal excitation frequency. Meanwhile, the first-order natural frequency ω1 and second harmonic frequency 2ω1 are also invoked and obviously seen in the spectral charts as shown in Figures 9(a), 9(c), and 9(f).Time histories of wave surface at the left tank wall and the corresponding spectral charts.
(a)
Time histories of wave surface for case A(ωh=ωv=ω1)
(b)
Spectral chart for Case A(ωh=ωv=ω1)
(c)
Time histories of wave surface for case B(ωh=ω1,ωv=0.7ω1)
(d)
Spectral chart for case B(ωh=ω1,ωv=0.7ω1)
(e)
Time histories of wave surface for case C(ωh=0.7ω1,ωv=ω1)
(f)
Spectral chart for case C(ωh=0.7ω1,ωv=ω1)


### 4.2. Effects of Horizontal and Vertical Excitation Intensities on Liquid Sloshing

As the excitation intensity is directly related to the sloshing nonlinearity, the effects of horizontal and vertical excitation intensities will be investigated in this section. The excitation intensity is measured by the corresponding parameter, namely,kh=Ahωh2/g and kv=Avωv2/g, in the present study.Firstly, the effect of horizontal excitation intensity is studied. Givenωh=0.7ω1, ωv=1.2ω1, kv=0.2 (i.e., Av=0.048 m), the horizontal excitation intensity parameter kh is taken as 0.0014 and 0.021 corresponding to small and large excitation amplitude of Ah=0.001 m and 0.015 m, respectively. Figure 10 illustrates the comparison of time series of wave surface at the left tank wall. From the figure, it can be seen that the wave phases are same for two cases but the wave nonlinearity is more apparent for the larger kh(kh=0.021), namely, higher and narrower wave crests, and flatter and lower wave troughs. Then, the FFT method is adopted to analyze time histories of wave surface shown in Figure 10, and the corresponding spectral charts are presented in Figure 11. It displays two dominating frequencies of which one is equivalent to the horizontal excitation frequency ωh, and the other is first-order natural frequency ω1. Their values in Figure 11(b) are larger than those in Figure 11(a) due to larger excitation intensity kh. More additional frequencies related to small spectral value, such as ω1-ωh,ω1+ωh,2ω1 and  ω3 are also induced by the stronger nonlinearity as shown in Figure 11(b).Figure 10
Time histories of wave surface at the left tank wall forkh=0.0014 and 0.021.Spectral charts corresponding to the Figure10.
(a)
k
h
=
0.0014
(b)
k
h
=
0.021Then, in order to investigate the effect of vertical excitation intensity, two cases withωh=0.7ω1, ωv=0.8ω1, kh=0.0014,kv=0.01 and 0.1 (i.e., Ah=0.001 m, Av=0.005 m and 0.05 m) are chosen. Figures 12 and 13 show the time histories of wave surface at the left tank wall for these two cases and the corresponding spectral charts. From the figures, it is found that although vertical excitation amplitudeAv increases ten times from 0.005 m to 0.05 m, the intensity of liquid sloshing does not increase accordingly and the wave surface value is only a little different, which is also represented in the corresponding spectral charts. The horizontal excitation frequency and first-order natural frequency are still the dominating frequencies. Thus it can be obtained that the effect of vertical excitation to liquid sloshing is weaker than that of horizontal excitation. But it has to be noticed that the sum-frequency and difference-frequency between horizontal frequency and vertical frequency occur due to the increase of vertical excitation intensity kv.Figure 12
Time histories of wave surface at the left tank wall forkv=0.01  and 0.1.Spectral charts corresponding to the Figure12.
(a)
k
v
=
0.01
(b)
k
v
=
0.1


### 4.3. Second-Order Resonance of Liquid Sloshing

Generally, the resonant phenomena can be excited when the excitation frequency is near to the odd natural frequencies of the container, but it is opposite for the even natural frequencies. Most previous researches were still limited to the first-order resonance problem, except that Wu (2007) [29] analytically studied the second-order resonance of sloshing due to half excitation frequency, sum-frequency, or difference-frequency equal to the even mode natural frequency. In this section, the second-order resonance of sloshing, excited by half second-order natural frequency, is numerically simulated and the relation with water depth is studied. As an example, a 1 m-length tank motion is governed by xb(t)=Ahcos(ωht), where Ah=0.001 m and ωh is taken as 0.5ω2. Water depth H is defined as 0.2 m, 0.3 m, 0.4 m, and 0.5 m, respectively. According to (3.1), the corresponding second-order natural frequency ω2 can be obtained as 7.238 rad/s, 7.671 rad/s, 7.798 rad, and 7.835 rad/s, respectively. Figure 14 shows the time histories of wave surface at left tank wall and the corresponding spectrum distribution for these four conditions. From the figure, it can be seen that liquid sloshing tends to resonance by time even if the excitation frequency is off natural frequency, especially for shallow water depth. Besides the first-order natural frequency, the second-order natural frequency is also provoked under such outer excitation. Both natural frequencies ω1 and ω2 are dominated in liquid sloshing compared with the excitation 0.5ω2, which coincides with the phenomenon of liquid resonance. The more shallow the water depth is, the more contribution of the second-order natural frequency results in, and the more violent the sloshing resonance is.Time history of wave surface at left tank wall and the spectral charts for excitation frequency 0.5ω2.
(a)
Time history of wave surface at left tank wall forH=0.2 m
(b)
Spectral chart of Figure14(a)
(c)
Time history of wave surface at left tank wall forH=0.3 m
(d)
Spectral chart of Figure14(c)
(e)
Time history of wave surface at left tank wall forH=0.4 m
(f)
Spectral chart of Figure14(e)
(g)
Time history of wave surface at left tank wall forH=0.5 m
(h)
Spectral chart of Figure14(g)


### 4.4. Effects of a Bottom-Mounted Vertical Baffle on Liquid Sloshing

Baffle is a kind of passive method to restrain liquid sloshing and has been applied in engineering widely. Many researchers have investigated liquid sloshing in baffled tank, but most of them supposed the tank is subjected to pure horizontal excitation. In fact, the tank movement is complicated and can be considered coupled movement of some basic movements. Thus, the effects of a vertical baffle on liquid sloshing are investigated comprehensively. Wave surface and wave force are chosen to reflect the corresponding effects. The same tank dimensions and motion are taken the same as the above case. An example ofωh=ωv=ω1, kh=kv=0.0029 (i.e., Ah=Av=0.001 m) is considered in this section.The effect of baffle position is firstly investigated. A vertical baffle of 0.4 m length (i.e.,Hb=0.8H) is arranged at the tank bottom (as shown in Figure 1). The thickness of the baffle is taken as 0.01 m and is supposed to be rigid. Figure 15 shows time histories of wave surface at the left and right walls versus varying baffle position D. It is observed that the closer to the center the vertical baffle is, the smaller is the wave surface. In other words, vertical baffle at center is most efficient to retrain liquid sloshing, which coincides with the common sense. Figure 16 shows the time histories of the total horizontal and vertical wave forces acting on the tank. It is obvious that the closer to the center the vertical baffle is, the smaller the total horizontal wave force is, which coincides with the conclusion of Figure 13. However, the total vertical wave force remains unchanged for different positions of the baffle.Time histories of wave surface at the left and right walls forHb=0.8H.
(a)
Left wall
(b)
Right wallTime histories of total horizontal and vertical wave forces forHb=0.8H.
(a)
horizontal wave force
(b)
vertical wave forceThe natural frequency of the above proposed tank will be different due to the presence of vertical baffle at the bottom. Figure17 shows the time histories of wave surface at the left tank wall in condition of different baffle positions under the excitations of various first-order natural frequencies and the motion governed by xb(t)=Ahcos(ωht), in which Ah=0.001 m and ωh is defined as the first-order natural frequency ω1, the baffle length Hb is taken as 0.75H, and baffle position D is taken as 0, 0.2L, and 0.4L, respectively. Based on the linear frequency-domain theory, the corresponding first-order natural frequency ω1 can be obtained as 4.342 rad/s, 4.581 rad/s and 5.134 rad/s, respectively. From the figure, it can be seen that the resonance phenomena still occur under the excitation of natural frequency as that without vertical baffle. However, for the smaller D, the sloshing waves propogate more slowly and the resonance amplitudes are also smaller, which presents a better effect on restraining sloshing.Figure 17
Time history of wave surface at left tank wall under the excitation of natural frequency.Furthermore, baffle length is also an important factor to influence liquid sloshing. In this part, a vertical baffle with varying length is arranged at the tank bottom-center. Figure18 presents time histories of wave surface at the left tank wall for three conditions of Hb=0.4H, Hb=0.6H, and Hb=0.8H. When baffle length is 0.4H, liquid sloshing is affected a little and wave surface is still near to resonant status. For a baffle of Hb=0.6H, wave surface becomes smaller obviously. When baffle length is increased to 0.8H, wave surface stabilizes at a small value. Meanwhile, with the increase of baffle length, the period of liquid sloshing becomes longer. Figure 19 shows the time series of total horizontal and vertical wave forces acting on the tank. It is obvious that the longer the vertical baffle is, the smaller is the total horizontal wave force. The total vertical wave force is almost invariant for these three conditions, which is the same as that shown in Figure 16(b).Figure 18
Time histories of wave surface at the left tank wall forD=0.0.Time histories of total horizontal and vertical wave forces forD=0.0.
(a)
horizontal wave force
(b)
vertical wave forceNext, liquid sloshing in the tank with different baffle length at the corresponding natural frequency is investigated. The baffle is fixed at the middle of tank bottom and baffle length is taken as 0, 0.6H, and 0.8H. Based on linear frequency-domain theory, the corresponding first-order natural frequency of tank can be obtained as 5.316 rad/s, 4.757 rad/s, and 4.12 rad/s, respectively. The tank motion is same as that described in Figure 17. Figure 20 gives the time history of wave surface at left tank wall in condition of different baffle lengths. From the figure, the similar conclusions can be obtained as those in Figure 17, that is, for higher Hb, the sloshing waves propogate more slowly, the resonance amplitudes are smaller, and a better effect on restraining sloshing is obtained.Figure 20
Time history of wave surface at left tank wall  under the excitation of first-order natural frequency.


## 4.1. Effects of Horizontal and Vertical Excitation Frequencies on Liquid Sloshing

As it is well known, liquid sloshing is most violent when external frequency approaches to natural frequency. Therefore, three cases are considered for the external excitation frequencies separately or both equal to natural frequency. Case A is that parametersωh=ωv=ω1,  kh=0.003, and kv=0.2 (i.e., Ah=0.001 m and  Av=0.069 m) are taken. Case B is that ωh=ω1,  ωv=0.7ω1, kh=0.003, and kv=0.2 (i.e., Ah=0.001 m and Av=0.142 m) are defined. Case C is that ωh=0.7ω1, ωv=ω1,  kh=0.0014, and kv=0.2 (i.e., Ah=0.001 m and Av=0.069 m) are chosen. Through numerical convergence test, temporal step is selected as Δt=0.012 s and spatial step is defined as Δx=Δz=0.05 m. Figure 9 shows time histories of wave surface at the left tank wall and the corresponding spectral charts obtained by FFT method. From the figure, It can be seen that the resonance phenomena occur as shown in Figures 9(a) and 9(c) when the horizontal excitation frequency is equal to natural frequency. Horizontal excitation plays a dominated role on the shape and peak amplitude of wave surface and can be always found in the spectral charts. However, the vertical excitation frequency does not exist alone in spectral chart and occurs as the forms of sum- and difference-frequency with horizontal excitation frequency. Meanwhile, the first-order natural frequency ω1 and second harmonic frequency 2ω1 are also invoked and obviously seen in the spectral charts as shown in Figures 9(a), 9(c), and 9(f).Time histories of wave surface at the left tank wall and the corresponding spectral charts.
(a)
Time histories of wave surface for case A(ωh=ωv=ω1)
(b)
Spectral chart for Case A(ωh=ωv=ω1)
(c)
Time histories of wave surface for case B(ωh=ω1,ωv=0.7ω1)
(d)
Spectral chart for case B(ωh=ω1,ωv=0.7ω1)
(e)
Time histories of wave surface for case C(ωh=0.7ω1,ωv=ω1)
(f)
Spectral chart for case C(ωh=0.7ω1,ωv=ω1)


## 4.2. Effects of Horizontal and Vertical Excitation Intensities on Liquid Sloshing

As the excitation intensity is directly related to the sloshing nonlinearity, the effects of horizontal and vertical excitation intensities will be investigated in this section. The excitation intensity is measured by the corresponding parameter, namely,kh=Ahωh2/g and kv=Avωv2/g, in the present study.Firstly, the effect of horizontal excitation intensity is studied. Givenωh=0.7ω1, ωv=1.2ω1, kv=0.2 (i.e., Av=0.048 m), the horizontal excitation intensity parameter kh is taken as 0.0014 and 0.021 corresponding to small and large excitation amplitude of Ah=0.001 m and 0.015 m, respectively. Figure 10 illustrates the comparison of time series of wave surface at the left tank wall. From the figure, it can be seen that the wave phases are same for two cases but the wave nonlinearity is more apparent for the larger kh(kh=0.021), namely, higher and narrower wave crests, and flatter and lower wave troughs. Then, the FFT method is adopted to analyze time histories of wave surface shown in Figure 10, and the corresponding spectral charts are presented in Figure 11. It displays two dominating frequencies of which one is equivalent to the horizontal excitation frequency ωh, and the other is first-order natural frequency ω1. Their values in Figure 11(b) are larger than those in Figure 11(a) due to larger excitation intensity kh. More additional frequencies related to small spectral value, such as ω1-ωh,ω1+ωh,2ω1 and  ω3 are also induced by the stronger nonlinearity as shown in Figure 11(b).Figure 10
Time histories of wave surface at the left tank wall forkh=0.0014 and 0.021.Spectral charts corresponding to the Figure10.
(a)
k
h
=
0.0014
(b)
k
h
=
0.021Then, in order to investigate the effect of vertical excitation intensity, two cases withωh=0.7ω1, ωv=0.8ω1, kh=0.0014,kv=0.01 and 0.1 (i.e., Ah=0.001 m, Av=0.005 m and 0.05 m) are chosen. Figures 12 and 13 show the time histories of wave surface at the left tank wall for these two cases and the corresponding spectral charts. From the figures, it is found that although vertical excitation amplitudeAv increases ten times from 0.005 m to 0.05 m, the intensity of liquid sloshing does not increase accordingly and the wave surface value is only a little different, which is also represented in the corresponding spectral charts. The horizontal excitation frequency and first-order natural frequency are still the dominating frequencies. Thus it can be obtained that the effect of vertical excitation to liquid sloshing is weaker than that of horizontal excitation. But it has to be noticed that the sum-frequency and difference-frequency between horizontal frequency and vertical frequency occur due to the increase of vertical excitation intensity kv.Figure 12
Time histories of wave surface at the left tank wall forkv=0.01  and 0.1.Spectral charts corresponding to the Figure12.
(a)
k
v
=
0.01
(b)
k
v
=
0.1


## 4.3. Second-Order Resonance of Liquid Sloshing

Generally, the resonant phenomena can be excited when the excitation frequency is near to the odd natural frequencies of the container, but it is opposite for the even natural frequencies. Most previous researches were still limited to the first-order resonance problem, except that Wu (2007) [29] analytically studied the second-order resonance of sloshing due to half excitation frequency, sum-frequency, or difference-frequency equal to the even mode natural frequency. In this section, the second-order resonance of sloshing, excited by half second-order natural frequency, is numerically simulated and the relation with water depth is studied. As an example, a 1 m-length tank motion is governed by xb(t)=Ahcos(ωht), where Ah=0.001 m and ωh is taken as 0.5ω2. Water depth H is defined as 0.2 m, 0.3 m, 0.4 m, and 0.5 m, respectively. According to (3.1), the corresponding second-order natural frequency ω2 can be obtained as 7.238 rad/s, 7.671 rad/s, 7.798 rad, and 7.835 rad/s, respectively. Figure 14 shows the time histories of wave surface at left tank wall and the corresponding spectrum distribution for these four conditions. From the figure, it can be seen that liquid sloshing tends to resonance by time even if the excitation frequency is off natural frequency, especially for shallow water depth. Besides the first-order natural frequency, the second-order natural frequency is also provoked under such outer excitation. Both natural frequencies ω1 and ω2 are dominated in liquid sloshing compared with the excitation 0.5ω2, which coincides with the phenomenon of liquid resonance. The more shallow the water depth is, the more contribution of the second-order natural frequency results in, and the more violent the sloshing resonance is.Time history of wave surface at left tank wall and the spectral charts for excitation frequency 0.5ω2.
(a)
Time history of wave surface at left tank wall forH=0.2 m
(b)
Spectral chart of Figure14(a)
(c)
Time history of wave surface at left tank wall forH=0.3 m
(d)
Spectral chart of Figure14(c)
(e)
Time history of wave surface at left tank wall forH=0.4 m
(f)
Spectral chart of Figure14(e)
(g)
Time history of wave surface at left tank wall forH=0.5 m
(h)
Spectral chart of Figure14(g)


## 4.4. Effects of a Bottom-Mounted Vertical Baffle on Liquid Sloshing

Baffle is a kind of passive method to restrain liquid sloshing and has been applied in engineering widely. Many researchers have investigated liquid sloshing in baffled tank, but most of them supposed the tank is subjected to pure horizontal excitation. In fact, the tank movement is complicated and can be considered coupled movement of some basic movements. Thus, the effects of a vertical baffle on liquid sloshing are investigated comprehensively. Wave surface and wave force are chosen to reflect the corresponding effects. The same tank dimensions and motion are taken the same as the above case. An example ofωh=ωv=ω1, kh=kv=0.0029 (i.e., Ah=Av=0.001 m) is considered in this section.The effect of baffle position is firstly investigated. A vertical baffle of 0.4 m length (i.e.,Hb=0.8H) is arranged at the tank bottom (as shown in Figure 1). The thickness of the baffle is taken as 0.01 m and is supposed to be rigid. Figure 15 shows time histories of wave surface at the left and right walls versus varying baffle position D. It is observed that the closer to the center the vertical baffle is, the smaller is the wave surface. In other words, vertical baffle at center is most efficient to retrain liquid sloshing, which coincides with the common sense. Figure 16 shows the time histories of the total horizontal and vertical wave forces acting on the tank. It is obvious that the closer to the center the vertical baffle is, the smaller the total horizontal wave force is, which coincides with the conclusion of Figure 13. However, the total vertical wave force remains unchanged for different positions of the baffle.Time histories of wave surface at the left and right walls forHb=0.8H.
(a)
Left wall
(b)
Right wallTime histories of total horizontal and vertical wave forces forHb=0.8H.
(a)
horizontal wave force
(b)
vertical wave forceThe natural frequency of the above proposed tank will be different due to the presence of vertical baffle at the bottom. Figure17 shows the time histories of wave surface at the left tank wall in condition of different baffle positions under the excitations of various first-order natural frequencies and the motion governed by xb(t)=Ahcos(ωht), in which Ah=0.001 m and ωh is defined as the first-order natural frequency ω1, the baffle length Hb is taken as 0.75H, and baffle position D is taken as 0, 0.2L, and 0.4L, respectively. Based on the linear frequency-domain theory, the corresponding first-order natural frequency ω1 can be obtained as 4.342 rad/s, 4.581 rad/s and 5.134 rad/s, respectively. From the figure, it can be seen that the resonance phenomena still occur under the excitation of natural frequency as that without vertical baffle. However, for the smaller D, the sloshing waves propogate more slowly and the resonance amplitudes are also smaller, which presents a better effect on restraining sloshing.Figure 17
Time history of wave surface at left tank wall under the excitation of natural frequency.Furthermore, baffle length is also an important factor to influence liquid sloshing. In this part, a vertical baffle with varying length is arranged at the tank bottom-center. Figure18 presents time histories of wave surface at the left tank wall for three conditions of Hb=0.4H, Hb=0.6H, and Hb=0.8H. When baffle length is 0.4H, liquid sloshing is affected a little and wave surface is still near to resonant status. For a baffle of Hb=0.6H, wave surface becomes smaller obviously. When baffle length is increased to 0.8H, wave surface stabilizes at a small value. Meanwhile, with the increase of baffle length, the period of liquid sloshing becomes longer. Figure 19 shows the time series of total horizontal and vertical wave forces acting on the tank. It is obvious that the longer the vertical baffle is, the smaller is the total horizontal wave force. The total vertical wave force is almost invariant for these three conditions, which is the same as that shown in Figure 16(b).Figure 18
Time histories of wave surface at the left tank wall forD=0.0.Time histories of total horizontal and vertical wave forces forD=0.0.
(a)
horizontal wave force
(b)
vertical wave forceNext, liquid sloshing in the tank with different baffle length at the corresponding natural frequency is investigated. The baffle is fixed at the middle of tank bottom and baffle length is taken as 0, 0.6H, and 0.8H. Based on linear frequency-domain theory, the corresponding first-order natural frequency of tank can be obtained as 5.316 rad/s, 4.757 rad/s, and 4.12 rad/s, respectively. The tank motion is same as that described in Figure 17. Figure 20 gives the time history of wave surface at left tank wall in condition of different baffle lengths. From the figure, the similar conclusions can be obtained as those in Figure 17, that is, for higher Hb, the sloshing waves propogate more slowly, the resonance amplitudes are smaller, and a better effect on restraining sloshing is obtained.Figure 20
Time history of wave surface at left tank wall  under the excitation of first-order natural frequency.


## 5. Conclusions

In the present paper, a higher-order boundary element method has been applied for the fully nonlinear numerical simulation of liquid sloshing in a 2D rectangular tank. The proposed numerical scheme has provided good results in comparison with the analytical solution, published numerical results, and experimental data.On the base of validations, a tank subjected to the coupled horizontal and vertical excitation is investigated for the effects of excitation frequencies and excitation intensities in the horizontal and vertical directions and vertical bottom-mounted baffle on liquid sloshing. The following conclusions can be obtained. Firstly, the horizontal excitation plays a dominated role on the wave shape and wave nonlinearity relative to the vertical excitation in the coupled sloshing problems. However, the latter can invoke more other higher and lower frequencies among natural frequency and base frequencies. Secondly, with the increasing of excitation intensities, liquid sloshing becomes more violent and the interaction of frequencies is enhanced. Thus more frequencies are invoked including sum frequency and difference frequency. Last, the vertical bottom-mounted baffle arranged at tank center is most efficient to restrain liquid sloshing due to the natural frequency of the tank changed. Under the excitation of the first-order natural frequency of unbaffled tank, both wave surface and horizontal wave load decrease with the increase of baffle length. However, the vertical baffle does little effect on the vertical wave load.



---

*Source: 340640-2012-06-27.xml*