Research Papers: Fundamental Issues and Canonical Flows

Investigation of Wave Characteristics in Oscillatory Motion of Partially Filled Rectangular Tanks

[+] Author and Article Information
M. Ozbulut

Engineering Faculty,
Piri Reis University,
Istanbul 34940, Turkey
e-mail: mozbulut@pirireis.edu.tr

N. Tofighi

Department of Mechanical Engineering,
University of Victoria,
Victoria, BC V8P 5C2, Canada
e-mail: ntofighi@uvic.ca

O. Goren

Faculty of Naval Architecture and
Ocean Engineering,
Istanbul Technical University,
Istanbul 34469, Turkey
e-mail: ogoren@itu.edu.tr

M. Yildiz

Integrated Manufacturing Technologies Research
and Application Center,
Sabanci University,
Tuzla 34956, Istanbul, Turkey;
Composite Technologies Center of Excellence,
Sabanci University-Kordsa,
Istanbul Technology Development Zone,
Sanayi Mah. Teknopark Blvd. No: 1/1B,
Pendik 34906, Istanbul, Turkey;
Faculty of Engineering and Natural Sciences,
Sabanci University,
Tuzla 34956, Istanbul, Turkey
e-mail: meyildiz@sabanciuniv.edu

1Corresponding author.

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received November 23, 2016; final manuscript received October 15, 2017; published online December 4, 2017. Assoc. Editor: Samuel Paolucci.

J. Fluids Eng 140(4), 041204 (Dec 04, 2017) (11 pages) Paper No: FE-16-1765; doi: 10.1115/1.4038242 History: Received November 23, 2016; Revised October 15, 2017

Simulations of oscillatory motion in partially filled rectangular tanks with different tank geometries, fullness ratios, and motion frequencies are presented. Smoothed particle hydrodynamics (SPH) method is used to discretize the governing equations together with new velocity variance-based free surface (VFS) and artificial particle displacement (APD) algorithms to enhance the robustness and the accuracy of the numerical scheme. Two-dimensional (2D) oscillatory motion is investigated for three different scenarios where the first one scrutinizes the kinematic characteristics in resonance conditions, the second one covers a wave response analysis in a wide range of enforced motion frequencies, and the last one examines the dynamic properties of the fluid motion in detail. The simulations are carried on for at least 50 periods in the wave response analysis. It is shown that numerical results of the proposed SPH scheme are in match with experimental and numerical findings of the literature.

Copyright © 2018 by ASME
Your Session has timed out. Please sign back in to continue.


Lee, S. J. , Kim, M. H. , Lee, D. H. , Kim, J. W. , and Kim, Y. H. , 2007, “ The Effect of LNG Tank Sloshing on the Global Motions of LNG Carriers,” Ocean Eng., 34(1), pp. 10–20. [CrossRef]
Delorme, L. , Iglesias, S. , and Perez, S. , 2005, “ Sloshing Loads Simulation in LNG Tankers With SPH,” International Conference on Computational Methods in Marine Engineering (MARINE), Oslo, Norway, June 27–29, pp. 1–10. https://www.researchgate.net/publication/228797725_Sloshing_loads_simulation_in_LNG_tankers_with_SPH
Delorme, L. , Colagrossi, A. , Iglesias, S. , Rodriguez, R. , and Botia-Vera, E. , 2009, “ A Set of Canonical Problems in Sloshing—Part I: Pressure Field in Forced Roll Comparison Between Experimental Results and SPH,” Ocean Eng., 36(2), pp. 168–178. [CrossRef]
Servan-Camas, B. , Cercos-Pita, J. L. , Garcia-Espinoza, J. , and Souto-Iglesias, A. , 2016, “ Time Domain Simulation of Coupled Sloshing Seakeeping Problems by SPH FEM Coupling,” Ocean Eng., 123(7), pp. 383–396. [CrossRef]
Faltinsen, O. M. , and Timokha, A. N. , 2009, Sloshing, Cambridge University Press, Cambridge, UK.
Okamoto, T. , and Kawahara, M. , 1990, “ Two Dimensional Sloshing Analysis by Lagrangian Finite Element Method,” Int. J. Numer. Methods Fluids, 11(5), pp. 453–477. [CrossRef]
Wang, C. , Teng, J. , and Huang, G. , 2011, “ Numerical Simulation of Sloshing Motion Inside a Two Dimensional Rectangular Tank by Level Set Method,” Int. J. Numer. Methods Heat Fluid Flow, 21(1), pp. 5–31. [CrossRef]
Guorong, Y. , Subhask, R. , and Kamran, S. , 2009, “ Experimental Study of Liquid Slosh Dynamics in a Partially-Filled Tank,” ASME J. Fluids Eng., 131(7), p. 071303. [CrossRef]
Shin, H. R. , 2005, “ Unstructured Grid Based Reynolds-Averaged Navier-Stokes Method for Liquid Tank Sloshing,” ASME J. Fluids Eng., 127(3), pp. 572–582. [CrossRef]
Gotoh, H. , Khayyer, A. , Ikari, H. , Arikawa, T. , and Shimosako, K. , 2014, “ On Enhancement of Incompressible SPH Method for Simulation of Violent Sloshing Flows,” Appl. Ocean Res., 46, pp. 104–115. [CrossRef]
Faltinsen, O. , and Timokha, A. , 2001, “ An Adaptive Multimodal Approach to Nonlinear Sloshing in a Rectangular Tank,” J. Fluid Mech., 432, pp. 167–200. https://www.cambridge.org/core/journals/journal-of-fluid-mechanics/article/an-adaptive-multimodal-approach-to-nonlinear-sloshing-in-a-rectangular-tank/D520FD5608316B33C3B778922EAF0FFC
Faltinsen, O. , and Timokha, A. , 2002, “ Asymptotic Modal Approximation of Nonlinear Resonant Sloshing in a Rectangular Tank With Small Fluid Depth,” J. Fluid Mech., 470, pp. 319–357. [CrossRef]
Chen, Y. , Djidjeli, K. , and Price, W. , 2009, “ Numerical Simulation of Liquid Sloshing Phenomena in Partially Filled Containers,” Comput. Fluids, 38(4), pp. 830–842. [CrossRef]
Monaghan, J. , and Gingold, R. , 1977, “ Smoothed Particle Hydrodynamics: Theory and Application to Non-Spherical Stars,” Mon. Not. R. Astron. Soc., 181(3), pp. 375–389. [CrossRef]
Lucy, L. , 1977, “ Numerical Approach to Testing the Fission Hypothesis,” Astron. J., 82, pp. 1013–1024. [CrossRef]
Koshizuka, S. , and Oka, Y. , 1996, “ Moving Particle Semi-Implicit Method for Fragmentation of Incompressible Fluid,” Nucl. Sci. Eng., 123(3), pp. 421–434. [CrossRef]
Khayyer, A. , and Gotoh, H. , 2013, “ Enhancement of Performance and Stability of MPS Mesh-Free Particle Method for Multiphase Flows Characterized by High Density Ratios,” J. Comput. Phys., 242, pp. 211–233. [CrossRef]
Antuono, M. , Colagrossi, A. , Marrone, S. , and Molteni, D. , 2010, “ Free Surface Flows Solved by Means of SPH Schemes With Numerical Diffusive Term,” Comput. Phys. Commun., 181(3), pp. 532–549. [CrossRef]
Rafiee, A. , Cummins, S. , Rudman, M. , and Thiagarajan, K. , 2012, “ Comparative Study on the Accuracy and Stability of SPH Schemes in Simulating Energetic Free Surface Flows,” Eur. J. Mech. B, 36, pp. 1–16. [CrossRef]
Chowdhurry, S. , and Sannasiraj, S. , 2014, “ Numerical Simulation of 2D Sloshing Waves Using SPH With Diffusive Terms,” Appl. Ocean Res., 47, pp. 219–240. [CrossRef]
Cao, X. Y. , Ming, F. R. , and Zhang, A. M. , 2014, “ Sloshing in a Rectangular Tank Based on SPH Simulation,” Appl. Ocean Res., 47, pp. 241–254. [CrossRef]
Meringolo, D. , Colagrossi, A. , Marrone, S. , and Aristodemo, F. , 2017, “ On the Filtering of Acoustic Components in Weakly-Compressible SPH Simulations,” J. Fluids Struct., 70, pp. 1–23. [CrossRef]
Pakozdi, C. , 2008, “ A Smoothed Particle Hydrodynamics Study of Two-Dimensional Nonlinear Sloshing in Rectangular Tanks,” Ph.D. thesis, Norwegian University of Science and Technology, Trondheim, Norway.
Monaghan, J. , and Kos, A. , 1999, “ Solitary Waves on a Cretan Beach,” ASCE J. Waterw. Port Coastal Ocean Eng., 125(3), pp. 145–154. [CrossRef]
Shadloo, M. , Zainali, A. , Sadek, H. , and Yildiz, M. , 2011, “ Improved Incompressible Smoothed Particle Hydrodynamics Method for Simulating Flow Around Bluff Bodies,” Comput. Methods Appl. Mech. Eng., 200(9–12), pp. 1008–1020. [CrossRef]
Liu, M. , and Liu, G. , 2010, “ Smoothed Particle Hydrodynamics: An Overview and Recent Developments,” Arch. Comput. Methods Eng., 17(1), pp. 25–76. [CrossRef]
von Neumann, J. , and Richtmyer, R. , 1950, “ A Method for the Numerical Calculation of Hydrodynamic Shocks,” J. Appl. Phys., 21(3), pp. 232–247. [CrossRef]
Antuono, M. , Colagrossi, A. , Marrone, S. , and Lugni, C. , 2011, “ Propagation of Gravity Waves Through an SPH Scheme With Numerical Diffusive Terms,” Comput. Phys. Commun., 182(4), pp. 866–877. [CrossRef]
Ozbulut, M. , Yildiz, M. , and Goren, O. , 2014, “ A Numerical Investigation Into the Correction Algorithms for SPH Method in Modeling Violent Free Surface Flows,” Int. J. Mech. Sci., 79, pp. 56–65. [CrossRef]
Antuono, M. , Colagrossi, A. , and Marrone, S. , 2012, “ Numerical Diffusive Terms in Weakly-Compressible SPH Schemes,” Comput. Phys. Commun., 183(12), pp. 2570–2580. [CrossRef]
Meringolo, D. , Colagrossi, A. , Aristodemo, F. , and Veltri, P. , 2015, “ SPH Numerical Modeling of Wave-Perforated Breakwater Interaction,” Coastal Eng., 101, pp. 48–68. [CrossRef]
Meringolo, D. , Colagrossi, A. , Aristodemo, F. , and Veltri, P. , 2015, “ Assessment of Dynamic Pressures at Vertical and Perforated Breakwaters Through Diffusive SPH Schemes,” Math. Probl. Eng., 2015, p. 305028.
Nomeritae, N. , Daly, E. , Grimaldi, S. , and Bui, H. , 2016, “ Explicit Incompressible SPH Algorithm for Free-Surface Flow Modelling: A Comparison With Weakly Compressible Schemes,” Adv. Water Resour., 97, pp. 156–167. [CrossRef]
Shepard, D. , 1968, “ A Two Dimensional Interpolation Function for Irregularly-Spaced Data,” 23rd ACM National Conference, Las Vegas, NV, Aug. 27–29, pp. 517–524. http://citeseerx.ist.psu.edu/viewdoc/download?doi=


Grahic Jump Location
Fig. 1

The view of free surface and fully populated regions during the evolution of sway-sloshing motion at t ∼ 20 s (a) and the conservation of final volume at t = 25 s for the dam-break problem with the local velocity coefficient in the APD treatment (b)

Grahic Jump Location
Fig. 7

Positions where pressure time series are calculated; dimensions are given in millimeters

Grahic Jump Location
Fig. 2

The representative geometry for the sway-sloshing simulations

Grahic Jump Location
Fig. 3

Comparison of the three different particle resolutions with experimental results (a) time histories of the free-surface elevations and (b) corresponding frequency analysis at x=−0.815 m for case 1-a

Grahic Jump Location
Fig. 4

Comparison of the (a) time histories of the free-surface elevations and (b) corresponding frequency analysis at x= −0.815 m for case 1-a

Grahic Jump Location
Fig. 5

Comparison of the (a) time histories of the free-surface elevations and (b) corresponding frequency analysis at the left wall at x= −0.815 m for case 1-b

Grahic Jump Location
Fig. 6

The comparison of wave elevations near wall with D/L = 0.35, A/L = 0.05

Grahic Jump Location
Fig. 8

Comparison of the pressure time histories at P1: (a) experimental results reported in Ref. [13], (b) results of this study, and (c) numerical results reported in Ref. [13]

Grahic Jump Location
Fig. 9

Quantitative comparison of the consecutive hits at pressure points P2 and P3. Horizontal lines indicate the mean impact pressure values where straight, dash-dot and dashed line types stands for the results of model test, present study and computational results of [13], respectively. (a) Consecutive pressure beat at point P2 and (b) consecutive pressure beat at point P3.

Grahic Jump Location
Fig. 10

Comparison of the pressure time histories at P4: (a) experimental results reported in Ref. [13], (b) results of this study, and (c) numerical results reported in Ref. [13]

Grahic Jump Location
Fig. 11

The snapshots of the free-surface profiles and pressure field between the time interval 10.54–11.17 s. The pressure values are given in kPa: (a) free-surface profile and pressure field at t1 = 10:54 s, (b) free-surface profile and pressure field at t2 = 10:64 s, (c) free-surface profile and pressure field at t3 = 10:86 s, and (d) free-surface profile and pressure field at t4 = 11:17 s.




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In