(ProQuest: ... denotes non-US-ASCII text omitted.)
De-Zhi Ning 1,2 and Wei-Hua Song 1 and Yu-Long Liu 1 and Bin Teng 1
Recommended by Armin Troesch
1, State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116024, China
2, Key Laboratory of Renewable Energy and Gas Hydrate, Chinese Academy of Sciences, Guangzhou 510640, China
Received 24 January 2012; Revised 22 April 2012; Accepted 2 May 2012
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.
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 Section 2. 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 Figure 1, 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.
[figure omitted; refer to PDF]
The displacement of tank is defined as Xb =[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: [figure omitted; refer to PDF] where [varphi] is the velocity potential. The following condition is satisfied on the interface of tank and liquid: [figure omitted; refer to PDF] 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: [figure omitted; refer to PDF] 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 : [figure omitted; refer to PDF] 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: [figure omitted; refer to PDF] Thus dynamic boundary condition is written as [figure omitted; refer to PDF] The velocity potential [varphi] can be split as follows in the space-fixed system: [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] The initial condition [varphi]=0 is written as [figure omitted; refer to PDF] 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 coordinate O0X0Z0 [figure omitted; refer to PDF]
Through (2.4), (2.5) and (2.10), (2.13) is transformed in tank-fixed coordinate OXZ [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] 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: [figure omitted; refer to PDF] 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: [figure omitted; refer to PDF]
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 [figure omitted; refer to PDF] 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 excitation xb (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
[figure omitted; refer to PDF]
(b) t/(h/G)0.5 =15.725
[figure omitted; refer to PDF]
The second set of validation tests is concerned with a tank subjected to pure vertical excitation. The tank scale is L/H=2.0 . The vertical excitation is zb (t)=Av cos (ωv t) , 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ωv 2 /g is used as a measure of the intensity of vertical excitation and [straight epsilon]=aω1 2 /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 [straight epsilon]=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.
[figure omitted; refer to PDF]
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 as xb (t)=Ah cos (ωh t) and zb (t)=Av cos (ωv t) . Parameter kh =Ahωh 2 /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.
[figure omitted; refer to PDF]
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 omitted; refer to PDF]
Figure 6: Distribution of the first-order natural frequency for different baffle length.
[figure omitted; refer to PDF]
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,η)
[figure omitted; refer to PDF]
(b) at the tank bottom (-L/2.0,-H)
[figure omitted; refer to PDF]
A tank with L=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
[figure omitted; refer to PDF]
(b) with baffle
[figure omitted; refer to PDF]
Through 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 as L=1.0 m and H=0.5 m, respectively. The horizontal and vertical excitations are taken as xb (t)=Ah cos (ωh t) and zb (t)=Av cos (ωv t) , respectively. Ah and Av are excitation amplitude. ωh and ωv are excitation frequency. Parameters kh =Ahωh 2 /g and kv =Avωv 2 /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 )
[figure omitted; refer to PDF]
(b) Spectral chart for Case A (ωh =ωv =ω1 )
[figure omitted; refer to PDF]
(c) Time histories of wave surface for case B (ωh =ω1 , ωv =0.7ω1 )
[figure omitted; refer to PDF]
(d) Spectral chart for case B (ωh =ω1 , ωv =0.7ω1 )
[figure omitted; refer to PDF]
(e) Time histories of wave surface for case C (ωh =0.7ω1 , ωv =ω1 )
[figure omitted; refer to PDF]
(f) Spectral chart for case C (ωh =0.7ω1 , ωv =ω1 )
[figure omitted; refer to PDF]
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ωh 2 /g and kv =Avωv 2 /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 for kh =0.0014 and 0.021.
[figure omitted; refer to PDF]
Spectral charts corresponding to the Figure 10.
(a) kh =0.0014
[figure omitted; refer to PDF]
(b) kh =0.021
[figure omitted; refer to PDF]
Then, 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 amplitude Av 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 for kv =0.01 and 0.1.
[figure omitted; refer to PDF]
Spectral charts corresponding to the Figure 12.
(a) kv =0.01
[figure omitted; refer to PDF]
(b) kv =0.1
[figure omitted; refer to PDF]
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)=Ah cos (ωh t) , 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 for H=0.2 m
[figure omitted; refer to PDF]
(b) Spectral chart of Figure 14(a)
[figure omitted; refer to PDF]
(c) Time history of wave surface at left tank wall for H=0.3 m
[figure omitted; refer to PDF]
(d) Spectral chart of Figure 14(c)
[figure omitted; refer to PDF]
(e) Time history of wave surface at left tank wall for H=0.4 m
[figure omitted; refer to PDF]
(f) Spectral chart of Figure 14(e)
[figure omitted; refer to PDF]
(g) Time history of wave surface at left tank wall for H=0.5 m
[figure omitted; refer to PDF]
(h) Spectral chart of Figure 14(g)
[figure omitted; refer to PDF]
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 for Hb=0.8H .
(a) Left wall
[figure omitted; refer to PDF]
(b) Right wall
[figure omitted; refer to PDF]
Time histories of total horizontal and vertical wave forces for Hb=0.8H .
(a) horizontal wave force
[figure omitted; refer to PDF]
(b) vertical wave force
[figure omitted; refer to PDF]
The natural frequency of the above proposed tank will be different due to the presence of vertical baffle at the bottom. Figure 17 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)=Ah cos (ωh t) , 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.
[figure omitted; refer to PDF]
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. Figure 18 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 for D=0.0 .
[figure omitted; refer to PDF]
Time histories of total horizontal and vertical wave forces for D=0.0 .
(a) horizontal wave force
[figure omitted; refer to PDF]
(b) vertical wave force
[figure omitted; refer to PDF]
Next, 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.
[figure omitted; refer to PDF]
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.
Acknowledgments
The authors gratefully acknowledge the financial support from the National Basic Research Program of China (973 Program, Grant no. 2011CB013703), the National Natural Science Foundation of China (Grant nos. 51179028 and 50921001), and the Open Fund of Key Laboratory of Renewable Energy and Gas Hydrate, Chinese Academy of Sciences (Grant no. Y007ka).
[1] O. M. Faltinsen, "A nonlinear theory of sloshing in rectangular tanks," Journal of Ship Research , vol. 18, no. 1074, pp. 224-241, ..
[2] J. B. Frandsen, "Sloshing motions in excited tanks," Journal of Computational Physics , vol. 196, no. 1, pp. 53-87, 2004.
[3] O. M. Faltinsen, O. F. Rognebakke, A. N. Timokha, "Classification of three-dimensional nonlinear sloshing in a square-base tank with finite depth," Journal of Fluids and Structures , vol. 20, no. 1, pp. 81-103, 2005.
[4] H. Akyildiz, N. Erdem Ünal, "Sloshing in a three-dimensional rectangular tank: numerical simulation and experimental validation," Ocean Engineering , vol. 33, no. 16, pp. 2135-2149, 2006.
[5] T. Okamoto, M. Kawahara, "Two-dimensional sloshing analysis by Lagrangian finite element method," International Journal of Numerical Methods in Fluids , vol. 11, no. 5, pp. 453-477, 1990.
[6] O. M. Faltinsen, "A numerical nonlinear method of sloshing in tanks with two dimensional flow," Journal of Ship Research , vol. 22, pp. 224-241, 1978.
[7] T. Nakayama, K. Washizu, "Nonlinear analysis of liquid motion in a container subjected to forced pitching oscillation," International Journal for Numerical Methods in Engineering , vol. 15, no. 8, pp. 1207-1220, 1980.
[8] G. X. Wu, Q. W. Ma, R. Eatock Taylor, "Numerical simulation of sloshing waves in a 3D tank based on a finite element method," Applied Ocean Research , vol. 20, no. 6, pp. 337-355, 1998.
[9] J. B. Frandsen, A. G. L. Borthwick, "Simulation of sloshing motions in fixed and vertically excited containers using a 2-D inviscid σ -transformed finite difference solver," Journal of Fluids and Structures , vol. 18, no. 2, pp. 197-214, 2003.
[10] C. Z. Wang, B. C. Khoo, "Finite element analysis of two-dimensional nonlinear sloshing problems in random excitations," Ocean Engineering , vol. 32, no. 2, pp. 107-133, 2005.
[11] V. Sriram, S. A. Sannasiraj, V. Sundar, "Numerical simulation of 2D sloshing waves due to horizontal and vertical random excitation," Applied Ocean Research , vol. 28, no. 1, pp. 19-32, 2006.
[12] M. Eswaran, U. K. Saha, D. Maity, "Effect of baffles on a partially filled cubic tank: numerical simulation and experimental validation," Computers and Structures , vol. 87, no. 3-4, pp. 198-205, 2009.
[13] J. R. Cho, H. W. Lee, S. Y. Ha, "Finite element analysis of resonant sloshing response in 2-D baffled tank," Journal of Sound and Vibration , vol. 288, no. 4-5, pp. 829-845, 2005.
[14] K. C. Biswal, S. K. Bhattacharyya, P. K. Sinha, "Non-linear sloshing in partially liquid filled containers with baffles," International Journal for Numerical Methods in Engineering , vol. 68, no. 3, pp. 317-337, 2006.
[15] P. K. Panigrahy, U. K. Saha, D. Maity, "Experimental studies on sloshing behavior due to horizontal movement of liquids in baffled tanks," Ocean Engineering , vol. 36, no. 3-4, pp. 213-222, 2009.
[16] D. Liu, P. Lin, "Three-dimensional liquid sloshing in a tank with baffles," Ocean Engineering , vol. 36, no. 2, pp. 202-212, 2009.
[17] R. Belakroum, M. Kadja, T. H. Mai, C. Maalouf, "An efficient passive technique for reducing sloshing in rectangular tanks partially filled with liquid," Mechanics Research Communications , vol. 37, no. 3, pp. 341-346, 2010.
[18] J. L. Ortiz, A. A. Barhorst, R. D. Robinett, "Flexible multibody systems-fluid interaction," International Journal for Numerical Methods in Engineering , vol. 41, no. 3, pp. 409-433, 1998.
[19] J. L. Ortiz, A. A. Barhorst, "Large-displacement non-linear sloshing in 2-D circular rigid containers--prescribed motion of the container," International Journal for Numerical Methods in Engineering , vol. 41, no. 2, pp. 195-210, 1998.
[20] J. L. Ortiz, A. A. Barhorst, "Closed-form modeling of fluid-structure interaction with nonlinear sloshing: potential flow," AIAA Journal , vol. 35, no. 9, pp. 1510-1517, 1997.
[21] J. L. Ortiz, A. A. Barhorst, "Modeling fluid-structure interaction," AIAA Journal of Guidance, Control, and Dynamics , vol. 20, no. 6, pp. 1221-1228, 1997.
[22] C. A. Brebbia, S. Walker Boundary Element Techniques in Engineering , Newnes-Butterworths, 1980.
[23] D.-Z. Ning, B. Teng, H.-T. Zhao, C.-I. Hao, "A comparison of two methods for calculating solid angle coefficients in a BIEM numerical wave tank," Engineering Analysis with Boundary Elements , vol. 34, no. 1, pp. 92-96, 2010.
[24] D. Z. Ning, B. Teng, "Numerical simulation of fully nonlinear irregular wave tank in three dimension," International Journal for Numerical Methods in Fluids , vol. 53, no. 12, pp. 1847-1862, 2007.
[25] D. Z. Ning, J. Zang, S. X. Liu, R. Eatock Taylor, B. Teng, P. H. Taylor, "Free-surface evolution and wave kinematics for nonlinear uni-directional focused wave groups," Ocean Engineering , vol. 36, no. 15-16, pp. 1226-1243, 2009.
[26] T. Okamoto, M. Kawahara, "Two-dimensional sloshing analysis by Lagrangian finite element method," International Journal for Numerical Methods in Fluids , vol. 11, no. 5, pp. 453-477, 1990.
[27] Y. S. Choun, C. B. Yun, "Sloshing characteristics in rectangular tanks with a submerged block," Computers and Structures , vol. 61, no. 3, pp. 401-413, 1996.
[28] M. Arafa, "Finite element analysis of sloshing in liquid-filled containers," in Production Engineering and Design for Development (PEDD '07), pp. 793-803, Cairo, Egypt, 2006.
[29] G. X. Wu, "Second-order resonance of sloshing in a tank," Ocean Engineering , vol. 34, no. 17-18, pp. 2345-2349, 2007.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Copyright © 2012 De-Zhi Ning et al. De-Zhi Ning 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.
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.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer