Abstract
In this paper, a multiphase method based on the Level Set Method is employed to study the sloshing phenomenon of two-layer liquid inside a two-dimensional rectangular container subjected to horizontal excitation. Validation of the multiphase approach is conducted through a comparative analysis with existing studies. Results show a fair agreement between the numerical model and available numerical and experimental data. Initially, a series of simulations were used to compare the sloshing behaviour of a two-layer fluid with that of a single-layer fluid. Even under identical external excitation, layered fluids demonstrate different sloshing patterns compared to single-layer liquids. Furthermore, the influence of the periodic excitation frequency on the sloshing dynamics was examined. Analysis was also conducted to explore the effect of internal baffles on the oscillatory behaviour of layered liquid sloshing. The findings reveal that the baffles significantly mitigate the sloshing of the layered fluid.
1 Introduction
Owing to the growing global demand for oil and gas, offshore technology such as floating production storage and offloading (FPSO) has recently begun to receive increased attention to accelerate the exploitation, extraction, and processing of offshore gas and oil. Crude oil extracted from deep-sea wells is initially a mixture of water and oil, piped to oily water separators in large tanks inside these floating units for processing. These separators operate on the principle of gravity settling to separate water from oil and require stable conditions for optimal operation. However, during their operations, these platforms are exposed to extremely harsh weather and sea conditions, such as waves, currents, and wind. Consequently, the motion of the floating platforms can cause intense sloshing of the layered liquid, which disturbs the equilibrium and makes obtaining the appropriate separation performance difficult.
Studies on the sloshing of multiphase fluid have been challenging over the past decade. Compared to single-layer liquid sloshing, multilayer sloshing presents more significant challenges due to variations in natural frequency through different layers [1]. Layered fluids (heterogeneous liquids) exhibit distinct sloshing behaviours compared to one-layer liquids (homogeneous liquids), even when subjected to the same external forces. Therefore, layered fluid sloshing has characteristics and behaviours that do not apply to assumptions generated using single-layer sloshing [2]. Therefore, it is essential to understand the impact of non-uniform density on the dynamic response of the layered fluid.
Most previous studies about liquid sloshing have mainly focused on single-layer fluid sloshing, while two-layer fluid sloshing cases have yet to be studied and analyzed. La Rocca et al. [3] conducted a theoretical investigation on the sloshing of a two-layer liquid, and they compared the mathematical model with experimental results to validate their study. Veletsos and Shivakumar [4] analyzed the sloshing behaviour of layered fluids within rectangular and cylindrical containers. Their study revealed that for each additional layer of fluid, there are also additional horizontal natural oscillation modes. This suggests that the complexity of the fluid layers contributes to a greater variety of sloshing patterns. Alongside theoretical studies, the sloshing of layered liquid was also studied numerically. Kargbo et al. [5] studied layered sloshing waves using the ALE method to track the fluid domain's free surface and interfacial wave. Korkmaz [6] studied the impact of viscosity on the sloshing of the stratified liquid, using high-speed imaging to track the deformations of the free surface. Liu et al. [7] conducted an experimental technique to investigate the sloshing characteristics of a three-layer fluid system.
Numerous studies have examined various elements and mechanisms in the control of sloshing [8]. One of the most common and effective methods used to suppress and mitigate fluid sloshing is the installation of baffles within tanks. These baffles work to dissipate sloshing energy and reduce pressure on tank walls through several mechanisms, such as increasing friction with the fluid and redirecting the flow, which reduces the momentum of the fluid and enables it to reduce its movement. Various shapes and types of baffles have been used to mitigate single-layer fluid sloshing, including perforated barriers, porous screens, slatted screens, and others. Few studies have been conducted on using baffles to study multilayer liquid sloshing. Therefore, this study aims to explore the effect of internal baffles on the behaviour of the sloshing of the two-layer liquid. Qin et al. [9] examined the hydroelasticity in the structural reaction of baffles with various shapes, demonstrating its importance and influence, while Zhang [10] studied the sloshing-damping effects of deformable baffles. Poguluri and Cho [11] developed analytical and numerical methods to study the sloshing behaviour of a two-layer liquid with vertical porous baffles. Bellezi [12] researched perforated bulkheads for sloshing mitigation. Guan et al. [13] utilized a Boundary Element Method to study the influence of different shapes of baffles on the sloshing. Results showed that baffle shape significantly impacts the sloshing of liquids. Cho et al. [14] investigated the utilization of a porous baffle to reduce sloshing in a partially filled tank. The porous baffle dissipated sloshing energy, reduced pressure on tank walls, and significantly reduced sloshing compared to a central baffle.
Despite previous research on multilayer liquid sloshing, most of it was primarily focused on theoretical analyses that are not valid for viscous and turbulent flows. Furthermore, numerical studies ignore difficulties caused by interface dynamics, nonlinear effects, viscosity effects, and surface tension. To address the limitations of existing methods, we have developed an innovative multiphase flow approach based on the Level Set method. This novel method seeks to capture the complex interactions between different liquid layers and their effects on sloshing behaviour, offering a fresh perspective on the problem.
This study presents a numerical model based on the Level Set Method to study the sloshing of two-layer liquid in a 2-D rectangular container under horizontal excitations. A series of simulations were used to compare the sloshing behaviour of a two-layer fluid with that of a single-layer fluid. In addition, the analysis illustrates the dynamic interactions associated with free surface waves and interface waves and how these interactions affect the fluid motion and stability in the tank. Moreover, the influence of the frequency amplitude of the periodic excitation on the dynamics of sloshing was examined. Analysis was also performed to investigate the impact of two horizontal baffles on the sloshing of layered fluids.
The paper comprises four sections. Section 1 covers the introduction above. The second section presents the numerical model of two-layer liquid sloshing, the validation schemes, and the mesh analysis of the developed numerical model. Section 3 discusses the numerical results. Finally, the fourth section outlines conclusions and further perspectives.
2 Numerical methods and verification schemes
2.1 Govern equations
More information regarding the Level Set Method is available in the COMSOL Multiphysics Users Guide [15].
2.2 Numerical method
The sloshing behaviour of two-layer liquid in a 2-D rectangular container has been analyzed using a computational model developed with COMSOL Multiphysics software. The model uses a multiphase flow technique based on the Level Set Method, which is well-known for its ability to simulate complex multiphase flows. In this model, the two immiscible fluids are assumed to be incompressible, Newtonian and the interaction between the fluids is modeled using the Navier-Stokes equations. The system is subjected to horizontal oscillations with a movement function represented by
Geometry of the model (Own source)
Citation: International Review of Applied Sciences and Engineering 2025; 10.1556/1848.2024.00868
FEM functions satisfy non-slip boundary conditions (on
2.3 Validation of the Level Set Method
The validity of the Level Set model is verified by comparing it with the numerical/experiment study conducted by Xue et al. [17]. Their experiments involved filling a tank, measuring
Figure 2 illustrates the oil–water phase interface wave at the right wall of the container. The results obtained using the multiphase Level Set Method demonstrate fair agreement with the IVOF-based numerical and the experiment conducted by Xue et al. [17].
Comparison of interface elevation between the present model and the simulation/experiment by Xue et al. [17] (Own source)
Citation: International Review of Applied Sciences and Engineering 2025; 10.1556/1848.2024.00868
2.4 Mesh analysis
To analyze the effect of the mesh size on the simulation generated by the model, a mesh convergence study was performed for the container configuration shown in Fig. 3. The mesh characteristics are divided into Coarse, Fine, and Normal mesh; the specific mesh properties are shown in Table 1. The simulations were performed using the same numerical parameters for the test model validation. This convergence is also related to the time integration step; in this example, the time step was set to t = 0.01 (Fig. 4).
Discretization of the domain (Normal mesh) (Own source)
Citation: International Review of Applied Sciences and Engineering 2025; 10.1556/1848.2024.00868
Mesh characteristics
Mesh type | Domain elements | Boundary elements | Maximal element size (m) | Minimal element size (m) | Curvature factor |
Fine | 12,274 | 480 | 0.00525 | 1.5E-4 | 0.3 |
Normal | 7,892 | 385 | 0.00675 | 3.0E-4 | 0.3 |
Coarse | 4,948 | 301 | 0.01 | 4.5E-4 | 0.4 |
Own source.
Effect of mesh type on free surface fluctuations (Own source)
Citation: International Review of Applied Sciences and Engineering 2025; 10.1556/1848.2024.00868
These results show that a mesh consisting of 4,948 domain elements and 301 boundary elements (considered as a Coarse mesh) is sufficient to discretize the considered geometry. Refining to a Normal mesh consisting of 7,892 domain elements and 385 boundary elements (Fig. 3), or even to a Fine mesh consisting of 12,274 domain elements and 480 boundary elements, did not significantly improve much in the results. It can be seen that the Normal and Fine meshes have significant similarities in terms of the temporal height of the free surface. Moreover, despite its lower discretization, the Coarse mesh also shows a strong similarity to the Normal mesh, indicating that the results remain acceptable even if the mesh is not too fine. Therefore, using a Normal mesh is valid and sufficient for this study, providing satisfactory accuracy while saving computational time.
3 Results and discussion
3.1 Sloshing analysis of layered liquids in a partially filled 2-D rectangular container
In this section, a series of numerical simulations using COMSOL Multiphysics Software were conducted to analyze the sloshing behaviour of a two-layer liquid in a 2-D rectangular container, as depicted in Fig. 5. The container is subjected to a horizontal excitation in the form of
Sketch of the geometry of the model in COMSOL (Own source)
Citation: International Review of Applied Sciences and Engineering 2025; 10.1556/1848.2024.00868
Physical properties of the two-layer liquid
Layered fluid | ||
Oil | 846.4 | 0.45 |
Water | 1,000 | 0.073 |
Own source.
The dynamics of the two-layer liquid sloshing are analyzed using the Level Set Method. The findings obtained from this analysis are compared with those of a single-layer liquid. The numerical results for the two-layered liquid sloshing are presented in detail in Figs 6–8. These figures provide a detailed analysis of the sloshing behaviour of the layered liquids and illustrate the differences between the two-layer and single-layer liquid sloshing.
Temporal comparison of free surface elevation generated by sloshing of two-layer and single-layer liquids, with
Citation: International Review of Applied Sciences and Engineering 2025; 10.1556/1848.2024.00868
Comparisons of free surface profiles of the sloshing of two-layer and single-layer liquids at several instants, with
Citation: International Review of Applied Sciences and Engineering 2025; 10.1556/1848.2024.00868
Comparison snapshots of the two-layer and single-layer liquids sloshing at several instants, with
Citation: International Review of Applied Sciences and Engineering 2025; 10.1556/1848.2024.00868
Figures 6 and 7 present a comparison between the free surface profiles of a two-layer liquid generated using the Multiphase Level Set Method and those of a single-layer liquid. The results indicate a significant difference in wave elevation and dynamic response. Even under identical external excitation, layered fluids reveal different sloshing patterns compared to single-layer liquids. This behaviour is potentially caused by density stratification of the layered liquid, and interfacial effects at the interface.
Furthermore, Fig. 8 shows snapshots captured at various time instants of the simulation of the two-layer liquid sloshing. The results show that the free surface elevation in a two-layer fluid is greater than a single-layer fluid, which is a key hydrodynamic characteristic of multi-layer fluid sloshing. The sloshing dynamics of multi-layered liquids are more complex and distinct from single-layered liquids, which can have significant implications in various industrial applications. These findings highlight the importance of considering the proprieties of layered liquids in sloshing systems.
3.2 Comparison of sloshing waves of the phase interface and the free surface of the two-layer-liquid
In this section, the sloshing characteristics of the interface waves of the layered liquid are analyzed and contrasted with those observed at the free surface.
Figures 9 and 10 compare free surface and interface wave profiles, which are two different types of waves that occur when a stratified fluid sloshes. Interface waves that propagate along the boundary between two immiscible liquids have a lower height than the free surface waves that propagate along the top surface of the liquid. This elevation difference is an important fluid dynamic property affecting stratified fluids' sloshing behaviour. It affects the stability of the fluid, the frequency of the waves, and the overall dynamics of the system. Therefore, understanding this difference is important to develop effective sloshing controls, and avoiding other safety risks.
Time history comparison of sloshing waves of the phase interface, and the free surface of the two-layer-liquid, with
Citation: International Review of Applied Sciences and Engineering 2025; 10.1556/1848.2024.00868
Time history comparison of sloshing waves of the phase interface, and the free surface of the two-layer-liquid, with
Citation: International Review of Applied Sciences and Engineering 2025; 10.1556/1848.2024.00868
3.3 Impact of frequency excitation on the dynamics of the layered liquid sloshing
In the following paragraph, we study the impact of varying frequency excitation on the dynamics of a two-layer liquid sloshing. This analysis is conducted using similar dimensions of the model considered in Fig. 5, to obtain comprehensive insights into the sloshing liquid's reaction to various excitation frequencies.
Figures 11 and 12 display the impact of frequency excitation on the behaviour of two-layer liquid sloshing. The sloshing amplitude follows a sinusoidal pattern due to the sinusoidal nature of the applied excitation. Furthermore, as the frequency amplitude of the periodic excitation rises, both the free surface and the elevation of the interfacial waves increase accordingly. The frequency of the periodic excitation plays a crucial role in shaping the free surface waves of layered liquid sloshing.
Time evolution of the interfacial layer waves at different values of the frequency amplitude of the horizontal excitation (Own source)
Citation: International Review of Applied Sciences and Engineering 2025; 10.1556/1848.2024.00868
Time evolution of the free surface waves at different values of the frequency amplitude of the horizontal excitation (Own source)
Citation: International Review of Applied Sciences and Engineering 2025; 10.1556/1848.2024.00868
3.4 Effect of two horizontal baffles on the layered liquid sloshing
Baffles are one of the well-known methods used to mitigate liquid sloshing. To investigate the effect of internal baffles on the sloshing characteristics of the two-layer liquid, a numerical model similar to the one depicted in Fig. 5 was used, with the integration of two horizontal baffles rigidly fixed at the right side and the left side of the tank, as shown in Fig. 13. The dimensions of the two horizontal baffles are given in Table 3. The movement of the tank was induced by a sinusoidal excitation function represented as
Sketch of the geometry of the model with two horizontal baffles (Own source)
Citation: International Review of Applied Sciences and Engineering 2025; 10.1556/1848.2024.00868
Description of horizontal baffles
Horizontal baffles | Length (m) | Width (m) | Position (m) |
Baffle 1 | 0.3 | 0.01 | x = 0, y = 0.3 |
Baffle 2 | 0.3 | 0.01 | x = 1.7, y = 0.3 |
Own source.
Figures 14 and 15 illustrate the interfacial layer and the free surface elevation at the right wall after installing the two horizontal baffles. The baffles effectively reduce sloshing amplitude. Furthermore, they demonstrate a notable effect in decreasing both the interfacial and the free surface wave height.
Comparison of phase interface elevation of the sloshing of the layered fluid with and without baffles (Own source)
Citation: International Review of Applied Sciences and Engineering 2025; 10.1556/1848.2024.00868
Comparison of the free surface elevation of the sloshing of the layered fluid with and without baffles (Own source)
Citation: International Review of Applied Sciences and Engineering 2025; 10.1556/1848.2024.00868
To visualize the interaction of the two-layer fluid with the horizontal baffles, screenshots are presented for different periods to compare the layered liquid sloshing with and without the two horizontal baffles, as shown in Fig. 16.
Comparisons of the snapshots of the sloshing of the layered fluid with and without baffles (Own source)
Citation: International Review of Applied Sciences and Engineering 2025; 10.1556/1848.2024.00868
Figure 16 shows two-layer liquid sloshing simulations for both baffled and unbaffled cases. The baffles act as a solid wall and effectively reduce the sloshing of the fluid. The region around the horizontal baffles experiences higher fluid velocity compared to the case without baffles. The velocity distribution differs significantly between the two cases, with noticeable vortex shedding near the tip of the baffle in the presence of baffles. This indicates the dissipation of fluid energy that reduces wave amplitude. Baffles interact with layered fluid, reducing turbulence and damping fluid oscillations. Overall, the internal baffles significantly reduce the free surface's elevation and the layered fluid's phase interface. This reduces impact forces that act on the tank walls and improves the structural integrity and safety of the system.
4 Conclusion
A numerical multiphase flow model based on the Level Set Method was built in COMSOL Multiphysics software to study the sloshing of two-layer fluid in a baffled tank under horizontal excitation. Validation of the multiphase approach is conducted through a comparative analysis with the numerical/experimental study carried out by Xue et al. [17], demonstrating favourable agreement in terms of interfacial wave elevations. Non-layered and layered liquid sloshing was analyzed. The results indicate a significant difference in wave elevation. The sloshing amplitude follows a sinusoidal pattern due to the sinusoidal nature of the applied excitation. Furthermore, as the frequency amplitude of the periodic excitation rises, both the free surface and the elevation of the interfacial waves increase accordingly. Even under identical external excitation, layered fluids demonstrate different sloshing patterns compared to single-layer liquids. This behaviour is potentially caused by density stratification of the layered liquid, and interfacial effects at the interface. Moreover, the influence of internal baffles on the oscillatory behaviour of layered liquid sloshing was examined, revealing that the presence of baffles greatly mitigates the sloshing of the layered fluid. The Level Set Method has been demonstrated to be highly suitable for the numerical analysis of the sloshing characteristics of multilayer liquids.
References
- [1]↑
M. Luo, X. Su, P. Lin, A. Khayyer, and X. Zhao, “Investigation of two-layer liquid sloshing by using the consistent particle method,”, vol. 32, no. 1, pp. 7–15, 2022. https://doi.org/10.17736/ijope.2022.jc848.
- [2]↑
H. Essaouini, L. El Bakkali, and P. Capodanno, “Analysis of small oscillations of a heavy almost-homogeneous liquid–gas system,” Mech. Res. Commun., vol. 37, no. 3, pp. 337–340, 2010. https://doi.org/10.1016/j.mechrescom.2010.01.003.
- [3]↑
M. La Rocca, G. Sciortino, C. Adduce, and M. A. Boniforti, “Experimental and theoretical investigation on the sloshing of a two-liquid system with free surface,” Phys. Fluids, vol. 17, no. 6, May 2005, Art no. 062101. https://doi.org/10.1063/1.1922887.
- [4]↑
A. S. Veletsos and P. Shivakumar, “Sloshing response of layered liquids in rigid tanks,” Earthq Eng. Struct. Dyn., vol. 22, no. 9, pp. 801–821, Sep. 1993. https://doi.org/10.1002/eqe.4290220906.
- [5]↑
O. Kargbo, M.-A. Xue, and J. Zheng, “Multiphase sloshing and interfacial wave interaction with a baffle and a submersed block,” J. Fluids Eng., vol. 141, no. 7, Jan. 2019. https://doi.org/10.1115/1.4041988.
- [6]↑
F. C. Korkmaz, “Damping of sloshing impact on bottom-layer fluid by adding a viscous top-layer fluid,” Ocean Eng., vol. 254, 2022, Art no. 111357. https://doi.org/10.1016/j.oceaneng.2022.111357.
- [7]↑
D. Liu, T. Lu, W. Cai, and L. Lei, “Experimental study on the sloshing of a three-layer liquid system with a free surface,” Ocean Eng., vol. 279, 2023. https://doi.org/10.1016/j.oceaneng.2023.114422.
- [8]↑
S. Koley, “Augmenting efficacy of global climate model forecasts: machine learning appraisal of remote sensing data,” Int. J. Eng. Trends Technol., vol. 72, no. 6, pp. 442–502, 2024. https://doi.org/10.14445/22315381/IJETT-V72I6P139.
- [9]↑
H. Qin, L. Mu, W. Tang, and Z. Hu, “Numerical study on structural response of anti-sloshing baffles of different configurations in a sloshing tank considering hydroelasticity,” Ocean Eng., vol. 188, 2019, Art no. 106290. https://doi.org/10.1016/j.oceaneng.2019.106290.
- [10]↑
Z. Zhang, “Numerical and experimental investigations of the sloshing modal properties of sloped-bottom tuned liquid dampers for structural vibration control,” Eng. Struct., vol. 204, 2020. https://doi.org/10.1016%2fj.engstruct.2019.110042.
- [11]↑
S. K. Poguluri and I. H. Cho, “Effect of vertical porous baffle on sloshing mitigation of two-layered liquid in a swaying tank,” Ocean Eng., vol. 289, 2023. https://doi.org/10.1016/j.oceaneng.2023.115952.
- [12]↑
C. A. Bellezi, L.-Y. Cheng, T. Okada, and M. Arai, “Optimized perforated bulkhead for sloshing mitigation and control,” Ocean Eng., vol. 187, 2019, Art no. 106171. https://doi.org/10.1016/j.oceaneng.2019.106171.
- [13]↑
Y. Guan, C. Yang, P. Chen, and L. Zhou, “Numerical investigation on the effect of baffles on liquid sloshing in 3D rectangular tanks based on nonlinear boundary element method,” Int. J. Naval Architect. Ocean Eng., vol. 12, pp. 399–413, 2020. https://doi.org/10.1016/j.ijnaoe.2020.04.002.
- [14]↑
I. H. Cho, J.-S. Choi, and M. H. Kim, “Sloshing reduction in a swaying rectangular tank by an horizontal porous baffle,” Ocean Eng., vol. 138, pp. 23–34, 2017. https://doi.org/10.1016/j.oceaneng.2017.04.005.
- [15]↑
Comsol, “The CFD Module User’s Guide,” 2023. [Online]. Available: https://doc.comsol.com/6.2/doc/com.comsol.help.cfd/CFDModuleUsersGuide.pdf. Accessed May 09, 2024.
- [16]↑
M.-A. Xue, O. Kargbo, and J. Zheng, “Seiche oscillations of layered fluids in a closed rectangular tank with wave damping mechanism,” Ocean Eng., vol. 196, 2020, Art no. 106842. https://doi.org/10.1016/j.oceaneng.2019.106842.
- [17]↑
M.-A. Xue, Z. Cao, X. Yuan, J. Zheng, and P. Lin, “A two-dimensional semi-analytic solution on two-layered liquid sloshing in a rectangular tank with a horizontal elastic baffle,” Phys. Fluids, vol. 35, no. 6, 2023. https://doi.org/10.1063/5.0153071.