Time-independent finite difference analysis of fully non-linear and viscous fluid sloshing in a rectangular tank

  • Authors:

  • Bang-Fuh Chen;Roger Nokes

  • Affiliations:

  • Department of Marine Environment and Engineering, National Sun Yat-Sen University, Kaohsiung 804, Taiwan and Department of Civil Engineering, University of Canterbury, Christchurch, New Zealand;Department of Civil Engineering, University of Canterbury, Christchurch, New Zealand

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2005

Quantified Score

Hi-index 31.46

Visualization

Abstract

A novel, time-independent finite-difference method for analyzing complete two-dimensional sloshing motion (surge, heave and pitch) in a tank has been developed based on the primitive 2D Navier-Stokes equations. Both the fully non-linear free surface condition and fluid viscosity are included. The boundary of the tank is mapped onto a fixed square domain through proper mapping functions and stretched meshes are employed near boundaries in order to more accurately evaluate the large disturbance of fluid along the boundary. The sloshing displacement agrees well with previously published results. The maximum transient amplitude is much larger than that of the steady-state. Clear beating phenomenon can be found when the tank is excited by near resonance frequency. The frequency dependence and Reynolds number effects are studied. For a fixed forcing-function amplitude, the sloshing response is greatest near resonance. An analysis under coupled surge and pitch motions is also made. The coupling effect is significant and simultaneous surge, heave and pitch motions should be included in the tank sloshing analysis. A simple formula is derived to approximate the horizontal force coefficient, C"F, on the tank walls. The formula implies that C"F is dominated by the free surface displacement when the tank is excited by small surge frequencies. Whereas C"F is attributed to added mass effects when the tank is under higher surge frequency forcing. A power spectra analysis is made to analyze the time series of sloshing displacement. For lower frequency of excitation, the system presents two peaks corresponding to the forcing frequency and fundamental frequency of the system. For higher frequency of excitation, the system shows only one major peak at the fundamental frequency. The limitations of the proposed method are also discussed.