Frequency Characteristics of Sloshing Resonance in a Three-Dimensional Shallow-Water Rectangular Tank

Abstract

The frequency characteristics of free surface elevation time histories of shallow-water sloshing in a three-dimensional rectangular tank is presented. The numerical model for sloshing motion uses an accurate velocity potential Boussinesq-type equation. A particular solution is adopted to express the external excitations and a linear damping term is introduced to replace a friction effect produced by tank walls. The excitation amplitude of the forced movement of the tank is tiny, while the excitation frequency varies around the first resonance frequency. The numerical results show that the wave energy in the rectangular tank under resonance is concentrated in the integer multiple of excitation frequencies, the sum excitation frequencies, and the difference excitation frequencies.

1. Introduction

Sloshing and its effects have become a crucial part of the structural and navigational safety of liquid-carrying ships such as LNG or FPSO. In addition to preventing the hazards due to sloshing, TLD and anti-rolling tanks that take advantage of the sloshing effect are also concerns of researchers [1,2,3,4]. Different relative filling levels in tanks can lead to significantly different regimes of liquid sloshing resonance. For intermediate water depth or more filling level cases, the nonlinear standing waves, local impact loads, and wave breaking phenomenon have been studied by the experimental methods [5,6,7], the numerical models based on Navier–Stokes equations [8,9,10], and inviscid incompressible flow theories [11,12,13].

For the shallow-water sloshing resonance, new shallow-water equations were derived by [14] under the assumption of inviscid vortical flow. The Boussinesq-type depth-averaged equations were used in [15,16] for a two-dimensional shallow-water sloshing. The Fourier decomposition was adopted for the reconstruction of the Boussinesq-type equations, and an artificial viscosity term was used for dissipation. The time histories of free surface elevations under the coupled sway-heave-roll excitation case were compared with the results calculated by the Smoothed Particle Hydrodynamics method.

For the strong nonlinear shallow water sloshing motion usually accompanied by wave breaking and hydraulic jump, the methods based on the Navier–Stokes equations and experimental methods were widely adopted. The experimental study of nonlinear sloshing characteristics in different shallow-water tanks was given in [17]. The results show that the excitation amplitude and dispersion parameters can greatly affect the resonance frequency. The hydraulic jump phenomenon in the shallow-water rectangular tank was analyzed in [18]. The frequency range of occurrence of hydraulic jump in tanks with different aspect ratios was discussed, and it has been found that the liquid depth plays a significant role.

Shallow-water sloshing in a rectangular tank was investigated in [19,20,21] by the Boussinesq-type equations with velocity potential as the variable to be solved. The Fourier series was applied for the particular solution of velocity potential, and a linear damping term was added in the dynamic boundary condition. The free surface profiles in tanks under six degrees of freedom excitations were shown and the single-bump, double-bump, and triple-bump traveling waves can be found when the excitation frequency changes around the first natural frequency. When the external excitations around the z-axis were added, the swirling waves appeared and were coupled with the nonlinear traveling waves.

The free surface profiles of shallow-water sloshing have an extremely complex wave compositions. In this paper, the power spectral densities of free surface elevations in tanks are used to analyze the wave compositions of the complex free surface profiles. Firstly, the numerical model based on the Boussinesq-type equations is built, and the linear damping term is given. Then, the numerical model is validated by comparing some calculations with experimental data in the literature. Finally, the power spectral densities of free surface elevations in the tank which suffered from different forced movements are shown. The results shown in this paper are useful for further understanding the regimes of nonlinear shallow-water sloshing motions.

2. Mathematical Formulation

2.1. Statement of the Problem

Sloshing of liquid in a rectangular tank is considered, symbol l is the length of the tank, b is the width, and h is the filling level. As shown in Figure 1, a Cartesian coordinate system $oxyz$ attached to the tank is adopted. $xoy$ plane coincides with the still water plane and the origin o with the center of it. The sides and bottom of the fluid are confined by the tank walls, and the upper boundary is the free surface $z=\eta (x,y,t)$. Assuming potential flow, there is a velocity potential $\mathsf{\Phi }(x,y,z,t)$ in the fluid region.

Nonlinear free surface boundary conditions are adopted. With the Laplace equations in fluid domain and boundary conditions at tank walls, the sloshing problem can be described as:

(1)$\begin{array}{cc}\hfill {\nabla }^2\mathsf{\Phi }& =0\hfill \end{array}$

(2)$\begin{array}{cc}\hfill \nabla \mathsf{\Phi }·\mathbf{n}={\mathbf{v}}_0·\mathbf{n}+\mathbf{w}·[\mathbf{r}\times \mathbf{n}]& =\mathbf{v}·\mathbf{n}\hfill \end{array}$

(3)$\begin{array}{cc}\hfill {\eta }_t-\mathbf{v}·\nabla \eta -{\mathsf{\Phi }}_z+{\eta }_x{\mathsf{\Phi }}_x+{\eta }_y{\mathsf{\Phi }}_y& =0\hfill \end{array}$

(4)$\begin{array}{cc}\hfill {\mathsf{\Phi }}_t-\mathbf{v}·\nabla \mathsf{\Phi }+{\displaystyle \frac{1}{2}}\nabla \mathsf{\Phi }·\nabla \mathsf{\Phi }+g\eta & =0\hfill \end{array}$

The symbol ∇ is the three-dimensional gradient operator defined as ∇ = ($\partial /\partial x,\partial /\partial y,$$\partial /\partial z$). $\mathbf{r}$ is the position vector, and $\mathbf{n}$ is the normal vector at the tank wall, and points tank outside. g is gravitational acceleration and $\mathbf{v}={\mathbf{v}}_0+\mathbf{w}\times \mathbf{r}=[v_a,v_b,v_c]$ represents the forced motion of the rectangular tank.

The translatory velocity along $(x,y,z)$ coordinate axes is expressed as ${\mathbf{v}}_0=[v_1,v_2,v_3]$. However, for the convenience of computation, the rotation axes of the angular velocity $\mathbf{w}=[w_1,w_2,w_3]$ are different. The dotted lines in Figure 1 show the rotation axis of the three angular velocity components, which are $(x,0,-h)$, $(0,y,-h)$, and z-axes, respectively.

2.2. Numerical Techniques

A linear damping term for replacing friction effect produced by tank walls is applied. According to [20], the linearized velocity potential ${\mathsf{\Phi }}_n$ is approximated by

(5)${\mathsf{\Phi }}_n=-a_n{\varphi }_n(x,y,z)sin({\omega }_{n}t+\theta )$

n represents the sloshing mode, and the parameter $ϵ$ is estimated by

(6)$ϵ={\displaystyle \frac{g}{{\omega }_n}}\sqrt{{\displaystyle \frac{\nu }{2{\omega }_n}}}{\displaystyle \frac{\int {\int }_{S_C}\nabla {\varphi }_n·\nabla {\varphi }_{n}dS}{\int {\int }_{S_F}{\varphi }_n^{2}dS}}$

(7)$ϵ=2{\omega }_n\sqrt{{\displaystyle \frac{\nu }{2{\omega }_n}}}{\displaystyle \frac{1}{bl}}(l+b+b{\displaystyle \frac{{\lambda }_{n}l-2{\lambda }_{n}h}{sinh2{\lambda }_{n}h}}){\lambda }_n=n\pi /l$

Here, $a_n$ is the amplitudes and ${\omega }_n$ the natural frequency. $S_C$ and $S_F$ represent the tank wall area and linear free surface area, respectively. $\nu$ is the kinematic viscosity. The velocity potential is divided into two parts $\mathsf{\Phi }=\varphi +\phi$, where $\varphi$ is a disturbance potential and $\phi$ is a particular solution satisfying the Laplace equation and boundary condition at the tank wall and bottom. According to [20], the particular solution considering the 6-DOF motion of the tank can be constructed as follows:

(8)$\phi =x{v}_1+y{v}_2+z{v}_3+w_1{\phi }_1+w_2{\phi }_2+w_3{\phi }_3$

(9)${\phi }_1=-z(y-b/2)+\sum _{n=1}^{\infty }{\beta }_{n}cos\left({\gamma }_{n}y\right){\displaystyle \frac{2sinh{\gamma }_n(z-h)}{{\gamma }_{n}cosh\left({\gamma }_{n}h\right)}}$

(10)${\phi }_2=-z(x-l/2)+\sum _{n=1}^{\infty }{\alpha }_{n}cos\left({\lambda }_{n}x\right){\displaystyle \frac{2sinh{\lambda }_n(z-h)}{{\lambda }_{n}cosh\left({\lambda }_{n}h\right)}}$

(11)$\begin{array}{c}\hfill {\phi }_3=\sum _{n=1}^{\infty }{\alpha }_{n}cos\left({\lambda }_{n}x\right){\displaystyle \frac{sinh{\lambda }_n(y-b/2)}{{\lambda }_{n}cosh({\lambda }_{n}b/2)}}\\ \hfill +\sum _{n=1}^{\infty }{\beta }_{n}cos\left({\gamma }_{n}y\right){\displaystyle \frac{sinh{\gamma }_n(x-l/2)}{{\gamma }_{n}cosh({\gamma }_{n}l/2)}}\end{array}$

where ${\lambda }_n=n\pi /l$, ${\alpha }_n=-4l/\left(n^2{\pi }^2\right)$, ${\gamma }_n=n\pi /b$, and ${\beta }_n=-4b/\left(n^2{\pi }^2\right)$. Based on the superposition principle, the decomposition of velocity potential leads to a simple boundary condition $\nabla \varphi ·\mathbf{n}=0$, which contributes to solving $\varphi$ by the Boussinesq-type approach. More details can be found in [19].

With initial conditions known, the time stepping process of $\mathsf{\Phi }$ and $\eta$ can be advanced by Equations (3) and (4). The central difference scheme and fourth-order Runge–Kutta method are used for spatial derivative calculation and temporal iteration, respectively.

3. Validation

The numerical model including the linear damping term has been valid and applied for the shallow-water sloshing in [19,20,21]. Comparison is made between the present results with the mode n of the particular solution equal to 9 and the experimental data given in [22]. According to [22], the time histories of wave elevations were measured by the resistance wave gauges fixed on the tank made of lucite, and a wetting agent was added to the water for reducing the influence of surface tension.

The geometric parameters of the tank are $l=$ 0.6095 m, $b=$ 0.23 m, and $h=$ 0.06 m. The harmonic movement of the tank occurs just in the x-direction, such that $X_{tank}=Asin\left(\omega t\right)$. The amplitude A is set to 0.196 cm, and excitation frequencies $\omega$ are chosen to be around the first natural frequency ${\omega }_r$, which was obtained from the linear solution for sloshing in the tank. As shown in Figure 2, a uniform grid with 31 nodes in the length direction and 11 nodes in the width direction is employed for calculation. The time step is set to 0.06, and the simulation has been over 100 periods to ensure that the sloshing motion reached a stable state.

The maximum transient dimensionless amplitudes at the tank wall versus normalized excitation frequencies are compared in Figure 3, where Exp represents the experimental data and Num the numerical results. The excitation frequencies divided by ${\omega }_r=3.8924{\mathrm{s}}^{-1}$ ranged from 0.9 to 1.12. The maximum value of the response curve appears around $\omega /{\omega }_r=1.05$ in experiment data while 1.06 in numerical data. The curves rise when approaching the maximum value from the left, but there is a sudden short pullback around $\omega /{\omega }_r=0.97$. Moving away from the maximum value from the right side, the curve drops sharply at first and then tends to decrease slowly. Overall, there are two discontinuities in the response curve. These discontinuities are produced by the effects of nonlinearities which result in a rapid change of sloshing motion form. More details about the sloshing motions can be found in [20].

In general, the results of the numerical model agree well with the experimental results. However, there were inconsistencies in the results from $\omega /{\omega }_r=$ 1.03 to 1.07. In addition to the occurrence of local breaking on free surface, the sensitivity of sloshing motion to the excitation frequency accounts for it. However, the correct prediction about location of the two discontinuities still shows that our numerical model has a good ability to capture the change of sloshing form caused by the variation of excitation frequency.

4. Results and Discussion

Based on the former numerical model, the free surface profiles in the rectangular tank are simulated. The sway, roll, coupled sway and roll, coupled surge and sway, coupled surge-sway-roll-pitch, and coupled surge ($v_1$)-sway ($v_2$)-roll ($w_1$)-pitch ($w_2$)-yaw ($w_3$) excitations are considered.

4.1. The Sway Excitations

A three-dimensional square tank with $l=b=0.5875$ m and $h=0.03$ m is considered. A uniform grid with 31 nodes in the length direction and 11 nodes in the width direction is employed for calculation. The time step is set to 0.06, and the simulation has been over 100 periods to ensure that the sloshing motion reached a stable state.

The translatory external excitations along the y-axis are adopted here. The excitation amplitude is 0.001 m, and the excitation frequencies vary between 0.92 ${\omega }_1$ to 1.08 ${\omega }_1$. Here, ${\omega }_1=2.8881{\mathrm{s}}^{-1}$ is the lowest natural frequency. According to the previous research, the violent sloshing motions appear in the vicinity of the lowest natural frequency and the transformation of free surface wave patterns happens. The time–history curves of the relative free surface elevation of the tank wall under different cases are displayed in Figure 4 and Figure 5. The free surface elevations $\eta$ divided by the still water depth h give the ordinate values and T are the excitation periods. From 0.92 ${\omega }_1$ to 1.08 ${\omega }_1$, several typical free surface wave patterns are shown in Figure 4 and Figure 5, including linear standing waves, nonlinear single-bump traveling waves, nonlinear double-bump traveling waves, and nonlinear triple-bump traveling waves.

As shown in Figure 6 and Figure 7, the power spectral densities of time histories of free surface elevations at an observation point within 600 s are calculated. The Matlab function $pwelch$ is used, and the frequencies divided by the excitation frequencies give the values of abscissa. As shown in the figures, the wave energy is distributed in the integer multiple of excitation frequencies. It should be noted that the rule can be applied to all the excitation frequencies from 0.92 ${\omega }_1$ to 1.08 ${\omega }_1$. For the excitation frequencies 0.92 ${\omega }_1$ and 1.08 ${\omega }_1$, the linear standing waves lead to the single frequency energy distribution. For the nonlinear double-bump traveling wave (Figure 5a) and triple-bump traveling wave (Figure 4c), the wave energy can be found on the double, triple, and higher multiple of excitation frequencies, but the main wave energy still focuses on the excitation frequency. In Figure 7b, the wave energy located on the double excitation frequency is higher than that on the excitation frequency, which leads to the nonlinear single-bump traveling waves (Figure 5b). Different free surface wave patterns in Figure 4 and Figure 5 are due to the superposition of different wave components (Figure 6 and Figure 7). The frequency distribution rule of the wave components has been found, but the amplitude distribution rule should be studied further.

4.2. The Roll Excitations

The three-dimensional rectangular tank used in Section 4.1 is adopted here and the following sections. The excitation amplitude is 0.001 rad and frequencies vary between 0.92 ${\omega }_1$ and 1.08 ${\omega }_1$. According to the numerical results, the free surface wave patterns of different excitation frequencies are similar to that in the tank excited by the sway motions. In Figure 8 and Figure 9, the time histories of relative free surface elevations and the power spectral densities with different excitation frequencies 0.96 ${\omega }_1$, 0.98 ${\omega }_1$ and 1.02 ${\omega }_1$ are given. The power spectral densities of free surface elevations are likewise distributed in the integer multiple of excitation frequencies. In addition, the amplitude distribution rule is similar to that in Section 4.1.

4.3. Coupled Sway and Roll Excitations

The amplitudes of sway excitations are 0.001 m and of roll excitations are 0.001 rad. The sway and roll motions are excited at the same frequencies, which ranged from 0.92 ${\omega }_1$ to 1.08 ${\omega }_1$. According to the numerical results, the wave amplitudes excited by the coupled sway and roll motions are much lower than that excited by the single sway motions or single roll motions. Meanwhile, the resonance frequency interval is narrower and closer to the natural frequency ${\omega }_1$. The conclusion and details can be found in [20]. For cases excited by the combined motion of sway and roll, the distribution positions of wave energy are also on the integer multiple of excitation frequencies (Figure 10). As shown in Figure 10c, the wave energy distribution indicates the linear standing waves present in the tank when the coupled excitation frequency equals 1.02 ${\omega }_1$. However, the strong nonlinear traveling waves can be observed under cases excited by the single sway motions or single roll motions at the same excitation frequency.

Next, the excitation frequency of sway motion is valued differently from the excitation frequency of roll motion. The sway excitation frequency is set to 1.02 ${\omega }_1$ and the roll excitation frequency is selected to 0.9 ${\omega }_1$. The excitation amplitudes are equal to 0.001 m and 0.001 rad. The power spectral densities of free surface elevations in the case of a sway and roll excitations are displayed in Figure 11a. In addition, the frequencies divided by the lowest natural frequency ${\omega }_1$ give the values of abscissa. As shown in Figure 11a, the wave energy distribution is not only on the integer multiple of excitation frequencies but also on the sum excitation frequencies and difference excitation frequencies. Figure 11a is enlarged, and the details are shown in Figure 11b and Figure 12. In Figure 11b, the locations of integer multiple of 0.9 and 1.02 are marked by the red ∘ and *, respectively, which means the integer multiple (N = 1–5) of sway excitation frequency 1.02 ${\omega }_1$ and roll excitation frequency 0.9 ${\omega }_1$ are denoted. Parts of power spectral density curves predictably coincide with the locations of integer multiple of excitation frequencies. In Figure 12a, the sum excitation frequency locations (0.9 + 1.02) ${\omega }_1$, (0.9 × 2 + 1.02) ${\omega }_1$, (0.9 + 1.02 × 2) ${\omega }_1$, (0.9 × 2 + 1.02 × 2) ${\omega }_1$, (0.9 × 3 + 1.02) ${\omega }_1$, (0.9 + 1.02 × 3) ${\omega }_1$, (0.9 × 3 + 1.02 × 2) ${\omega }_1$, (0.9 × 2 + 1.02 × 3) ${\omega }_1$, (0.9 × 4 + 1.02) ${\omega }_1$, (0.9 + 1.02 × 4) ${\omega }_1$ are marked by the red +. In addition, parts of power spectral density curves coincide with the sum excitation frequency locations. In Figure 12b, the difference excitation frequency locations (1.02 − 0.9 + 1.02) ${\omega }_1$, (1.02 − 0.9 + 1.02 × 2) ${\omega }_1$, (1.02 − 0.9 + 1.02 × 3) ${\omega }_1$, (1.02 − 0.9 + 1.02 × 4) ${\omega }_1$ and (1.02 − 0.9 + 1.02 × 5) ${\omega }_1$ are denoted by the red ×. The rest of power spectral density curves are located on these difference excitation frequencies. In summary, when the excitation frequencies of sway and roll motions are different, the power spectral density distributions of free surface elevation are compound, including the integer multiple of excitation frequencies, the sum excitation frequencies, and the difference excitation frequencies.

4.4. Coupled Surge and Sway Excitations with Different Amplitudes

The surge and sway motions have the same excitation frequency 1.02 ${\omega }_1$ but different excitation amplitudes. The direction of external combined excitation is not along the diagonal of $oxy$ plane in the tank since the surge motion has a 0.001 m excitation amplitude, while the sway motion has a 0.002 m excitation amplitude. The power spectral densities of free surface elevations at four tank corners $(l/2,b/2)$, $(-l/2,b/2)$, $(l/2,-b/2)$ and $(-l/2,-b/2)$ marked by A, B, C and D are shown in Figure 13 and Figure 14. As shown in the figures, the power spectral density distributions at the corners A and D are the same, and the wave energy locates on the integer multiple of excitation frequency 1.02 ${\omega }_1$. The power spectral densities of the other two tank corners B and C are also the same. There is an angle between the direction of combined excitations and the diagonal of $oxy$ plane. The diagonal component of combined excitations is dominant. As shown in Figure 13 and Figure 14, the spectral densities located on the excitation frequency at the corners A and D are much higher than that at the corners B and C.

4.5. Coupled Surge-Sway-Roll-Pitch Excitations

In this section, the coupled surge-swing-roll-pitch excitation is researched. The excitation amplitude of the translational motion excitation (surge and sway) is 0.001 m and 0.001 rad of the rotational excitation (roll and pitch). The excitation frequency of the coupled motions is selected to be 0.98 ${\omega }_1$. The power spectral density of free surface elevations at the tank corner ($l/2$, $b/2$) is shown in Figure 15. The wave energy distribution in the tank can be found on the integer multiple of excitation frequency 0.98 ${\omega }_1$. Using the same excitation frequency and amplitude, the superposition of results under coupled sway-roll excitations and coupled surge-pitch excitations equal results of the coupled surge-sway-roll-pitch case. However, the free surface height at the corner of the tank cannot be obtained by the same superposition due to the combined effect of multiple excitations. Compared to the results in Figure 10a, the spectral density due to coupled surge-sway-roll-pitch excitations cannot be calculated by the double spectral density due to coupled sway-roll excitations.

4.6. Coupled Surge-Sway-Roll-Pitch-Yaw Excitations

The excitation frequencies of surge, sway, roll, and pitch motions are equal to ${\omega }_1$. For the yaw excitations, the excitation frequency is calculated by ${\omega }_r^2=g\nu tanh\nu h$ and $\nu =\sqrt{{(\pi /l)}^2+{(\pi /b)}^2}$. The excitation amplitude of the translational motion(surge and sway) excitation is 0.001 m and 0.001 rad of the rotational excitation (roll, pitch and yaw). The power spectral densities of free surface elevations at the tank corners ($l/2,b/2$; $-l/2,b/2$) are calculated. Because the lowest natural frequencies of yaw excitation ${\omega }_r=4.0672{\mathrm{s}}^{-1}$ and others are different, the abscissa is valued by the frequency.

As shown in Figure 16a, the integer multiple of excitation frequency (${\omega }_1$) is marked by red *, and the integer multiple of excitation frequency (${\omega }_r$) is denoted by red o. The sum excitation frequencies ${\omega }_1+{\omega }_r$ and $2{\omega }_1+\omega {·}_r$ are expressed by red +. The wave energy distribution of the free surface elevations at the tank corner A ($l/2,b/2$) can be found clearly. At the other corner B ($-l/2,b/2$), the wave energy is focused on the integer multiple of excitation frequency ${\omega }_r$ and excitation frequency ${\omega }_1$ doubling. The power spectral density at the excitation frequency ${\omega }_1$ can not be found because the corner is not located on the direction of coupled surge-sway-roll-pitch excitations.

5. Conclusions

The shallow-water sloshing resonance in the three-dimensional tank was researched by Boussinesq-type equations with the linear damping term. The power spectral densities of free surface elevations in tank under resonance were analyzed for different external excitations which include single degree of freedom and multiple degrees of freedom excitations. Research on the sloshing phenomenon of external excitation frequency near the first natural frequency shows that: (1) According to the Figure 7, Figure 9, and Figure 10, whether it is a single degree of freedom excitation or two degrees of freedom excitations with the same excitation frequency, the power spectral density distributions are mainly at the integer multiple of excitation frequencies. (2) Combining the curves in Figure 11 and Figure 12, when the excitation frequencies of two directions are different, the power spectral density is distributed in the difference excitation frequencies and sum excitation frequencies in addition to the integral multiple excitation frequency. (3) The coupling effect of external excitation on sloshing motion is highly nonlinear. By comparing Figure 15 and Figure 10a, it can be found that the spectral density due to coupled surge-sway-roll-pitch excitations cannot be calculated by the double spectral density due to coupled sway-roll excitations. (4) At present, the power spectral densities of the free surface elevations has been given, and the power spectral densities of the forces should be consistent with it. The wave loads on the tank walls should be further studied.

Understanding the energy spectrum characteristics of sloshing motion can serve to establish a reduced-order model in the future, in which we can focus on these frequencies of the main distribution of sloshing energy. The amplitude distribution rules of power spectral densities should be studied further.

Author Contributions

Conceptualization, Y.S.; validation, X.Y.; formal analysis, X.Y. and P.X.; writing and editing, X.Y. and Y.S. All authors have read and agreed to the published version of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures

Figure 1. The coordinate system and rotation axis in a tank model.

Figure 2. The uniform grid on the free surface in the three−dimensional tank.

Figure 3. Variation of the maximum transient amplitudes divided by the water depth with excitation frequencies around the first natural frequency. (a) table; (b) figure.

Figure 4. Dimensionless free surface elevations time histories of pure sway cases under different excitation frequencies: (a) [Forumla omitted. See PDF.] = 0.92; (b) [Forumla omitted. See PDF.] = 0.94; (c) [Forumla omitted. See PDF.] = 0.96.

Figure 5. Dimensionless free surface elevations time histories of pure sway cases under different excitation frequencies: (a) [Forumla omitted. See PDF.] = 0.98; (b) [Forumla omitted. See PDF.] = 1.02; (c) [Forumla omitted. See PDF.] = 1.08.

Figure 6. The frequency spectrum of pure sway cases under different excitation frequencies: (a) [Forumla omitted. See PDF.] = 0.92; (b) [Forumla omitted. See PDF.] = 0.94; (c) [Forumla omitted. See PDF.] = 0.96.

Figure 7. The frequency spectrum of pure sway cases under different excitation frequencies: (a) [Forumla omitted. See PDF.] = 0.98; (b) [Forumla omitted. See PDF.] = 1.02; (c) [Forumla omitted. See PDF.] = 1.08.

Figure 8. Dimensionless free surface elevations time histories of pure roll cases under different excitation frequencies: (a) [Forumla omitted. See PDF.] = 0.96; (b) [Forumla omitted. See PDF.] = 0.98; (c) [Forumla omitted. See PDF.] = 1.02.

Figure 9. The frequency spectrum of pure roll cases under different excitation frequencies: (a) [Forumla omitted. See PDF.] = 0.96; (b) [Forumla omitted. See PDF.] = 0.98; (c) [Forumla omitted. See PDF.] = 1.02.

Figure 10. The frequency spectrum of coupled sway and roll cases under different excitation frequencies: (a) [Forumla omitted. See PDF.] = 0.98; (b) [Forumla omitted. See PDF.] = 1.0; (c) [Forumla omitted. See PDF.] = 1.02.

Figure 11. The frequency spectrum of coupled sway and roll cases under different excitation frequencies. excitation frequency: (a) [Forumla omitted. See PDF.] = 1.02, excitation amplitude: 0.001 m (sway); excitation frequency: [Forumla omitted. See PDF.] = 0.90, excitation amplitude: 0.001 rad (roll); (b) integer multiple.

Figure 12. The frequency spectrum of coupled sway and roll cases under different excitation frequencies. excitation frequency: [Forumla omitted. See PDF.] = 1.02, excitation amplitude: 0.001 m (sway); excitation frequency: [Forumla omitted. See PDF.] = 0.90, excitation amplitude: 0.001 rad (roll). (a) sum frequency; (b) difference frequency.

Figure 13. The frequency spectrum of coupled surge and sway cases under different excitation frequencies:. excitation frequency: [Forumla omitted. See PDF.] = 1.02, excitation amplitude: 0.001 m (surge); excitation frequency: [Forumla omitted. See PDF.] = 1.02, excitation amplitude: 0.002 m (sway). (a) corner A; (b) corner B.

Figure 14. The frequency spectrum of coupled surge and sway cases under different excitation frequencies. excitation frequency: [Forumla omitted. See PDF.] = 1.02, excitation amplitude: 0.001 m (surge); excitation frequency: [Forumla omitted. See PDF.] = 1.02, excitation amplitude: 0.002 m (sway). (a) corner C. (b) corner D.

Figure 15. The frequency spectrum of coupled surge−sway−roll−pitch excitations cases at the tank corner under different excitation frequencies. Excitation frequency is 0.98 [Forumla omitted. See PDF.].

Figure 16. The frequency spectrum of coupled surge−sway−roll−pitch−yaw excitations cases at the tank corner under under different excitation frequencies. (a) corner A: ([Forumla omitted. See PDF.]); (b) corner B: ([Forumla omitted. See PDF.]).

References

1. Zhang, Q.; Shui, B.; Zhu, H. Study on Sloshing Characteristics in a Liquid Cargo Tank under Combination Excitation. J. Mar. Sci. Eng.; 2022; 10, 1100. [DOI: https://dx.doi.org/10.3390/jmse10081100]

2. Wang, T.; Jin, H.; Lou, M.; Wang, X.; Liu, Y. Motion Responses of a Berthed Tank under Resonance Coupling Effect of Internal Sloshing and Gap Flow. Water; 2021; 13, 3625. [DOI: https://dx.doi.org/10.3390/w13243625]

3. Xue, M.A.; Dou, P.; Zheng, J.; Lin, P.; Yuan, X. Pitch motion reduction of semisubmersible floating offshore wind turbine substructure using a tuned liquid multicolumn damper. Mar. Struct.; 2022; 84, 103237. [DOI: https://dx.doi.org/10.1016/j.marstruc.2022.103237]

4. Liu, X.; Zhao, W.; Wan, D. Optimization of the roll motion of box-shaped hull section with anti-rolling sloshing tanks and fins in beam waves. J. Hydrodyn.; 2021; 33, pp. 688-697. [DOI: https://dx.doi.org/10.1007/s42241-021-0067-z]

5. Karimi, M.; Brosset, L.; Kaminski, M.; Ghidaglia, J. Effects of ullage gas and scale on sloshing loads. Eur. J. Mech. B/Fluids; 2017; 62, pp. 59-85. [DOI: https://dx.doi.org/10.1016/j.euromechflu.2016.11.017]

6. Xue, M.; Zheng, J.; Lin, P.; Yuan, X. Experimental study on verticle baffles of different configurations in suppressing sloshing pressure. Ocean. Eng.; 2017; 136, pp. 178-189. [DOI: https://dx.doi.org/10.1016/j.oceaneng.2017.03.031]

7. Lee, J.; Ahn, Y.; Kim, Y. Experimental study on effect of density ratio and phase transition during sloshing impact in rectangular tank. Ocean. Eng.; 2021; 242, 110105. [DOI: https://dx.doi.org/10.1016/j.oceaneng.2021.110105]

8. Liu, D.; Lin, P. A numerical study of three-dimensional liquid sloshing in tanks. J. Comput. Phys.; 2008; 227, pp. 3921-3939. [DOI: https://dx.doi.org/10.1016/j.jcp.2007.12.006]

9. Jin, X.; Lin, P. Viscous effects on liquid sloshing under external excitations. Ocean. Eng.; 2019; 171, pp. 695-707. [DOI: https://dx.doi.org/10.1016/j.oceaneng.2018.10.024]

10. Xue, M.; Jiang, Z.; Hu, Y.; Yuan, X. Numerical study of porous material layer effects on mitigating sloshing in a membrane LNG tank. Ocean. Eng.; 2020; 218, 108240. [DOI: https://dx.doi.org/10.1016/j.oceaneng.2020.108240]

11. Frandsen, J. Sloshing motions in excited tanks. J. Comput. Phys.; 2004; 196, pp. 53-87. [DOI: https://dx.doi.org/10.1016/j.jcp.2003.10.031]

12. Zhang, C. Nonlinear simulation of resonant sloshing in wedged tanks using boundary element method. Eng. Anal. Bound. Elem.; 2016; 69, pp. 1-20. [DOI: https://dx.doi.org/10.1016/j.enganabound.2016.04.006]

13. Zang, Q.; Liu, J.; Yu, L.; Lin, G. Boundary element analysis of liquid sloshing characteristics in axisymmetric tanks with various porous baffles. Appl. Ocean. Res.; 2019; 93, 101963. [DOI: https://dx.doi.org/10.1016/j.apor.2019.101963]

14. Ardakani, H.; Bridges, T. Shallow-water sloshing in vessels undergoing prescribed rigid-body motion in three dimensions. J. Fluid Mech.; 2011; 667, pp. 474-519. [DOI: https://dx.doi.org/10.1017/S0022112010004477]

15. Antuono, M.; Bouscasse, B.; Colagrossi, A.; Lugni, C. Two-dimensional modal method for shallow-water sloshing in rectangular basins. J. Fluid Mech.; 2012; 700, pp. 419-440. [DOI: https://dx.doi.org/10.1017/jfm.2012.140]

16. Antuono, M.; Bardazzi, A.; Lugni, C.; Brocchini, M. A shallow-water sloshing model for wave breaking in rectangular tanks. J. Fluid Mech.; 2014; 746, pp. 437-465. [DOI: https://dx.doi.org/10.1017/jfm.2014.127]

17. Gurusamy, S.; Kumar, D. Experimental study on nonlinear sloshing frequency in shallow water tanks under the effects of excitation amplitude and dispersion parameter. Ocean. Eng.; 2020; 213, 107761. [DOI: https://dx.doi.org/10.1016/j.oceaneng.2020.107761]

18. Gurusamy, S.; Sanapala, V.; Kumar, D.; Patnaik, B. Sloshing dynamics of shallow water tanks: Modal characteristics of hydraulic jumps. J. Fluids Struct.; 2021; 104, 103322. [DOI: https://dx.doi.org/10.1016/j.jfluidstructs.2021.103322]

19. Su, Y.; Yuan, X.; Liu, Z. Numerical research of three-dimensional nonlinear sloshing in Shallow-water rectangular tank. Appl. Ocean. Res.; 2020; 101, 102256. [DOI: https://dx.doi.org/10.1016/j.apor.2020.102256]

20. Su, Y. Numerical model of energy dissipation and shallow-water sloshing motions in tank under different coupled excitations. J. Mar. Sci. Technol.; 2021; 26, pp. 1038-1050. [DOI: https://dx.doi.org/10.1007/s00773-021-00797-y]

21. Su, Y. Transformation of wave patterns in three-dimensional shallow-water tank due to vertical excitations. Appl. Ocean. Res.; 2021; 115, 102847. [DOI: https://dx.doi.org/10.1016/j.apor.2021.102847]

22. Lepelletier, T.; Raichlen, F. Nonlinear oscillations in rectangular tanks. J. Eng. Math.; 1988; 114, pp. 1-23. [DOI: https://dx.doi.org/10.1061/(ASCE)0733-9399(1988)114:1(1)]