0
TECHNICAL PAPERS

# Application of a Discontinuous Galerkin Finite Element Method to Liquid Sloshing

[+] Author and Article Information
Martin J. Guillot

Department of Mechanical Engineering, University of New Orleans, New Orleans, LA 70148mjguillo@uno.edu

J. Offshore Mech. Arct. Eng 128(1), 1-10 (May 25, 2005) (10 pages) doi:10.1115/1.2151204 History: Received April 20, 2004; Revised May 25, 2005

## Abstract

A Runge-Kutta discontinuous Galerkin (RKDG) finite element method is applied to the liquid sloshing problem using the depth-averaged shallow water equations in a rotating frame of reference. A weak statement formulation is developed by multiplying the equations by a test function and integrating over a typical element. The basis functions are Legendre polynomials of degree one, resulting in formally second-order spatial accuracy. Second-order time integration is achieved using a second-order Runge-Kutta method. A minmod slope limiter is incorporated into the solution near discontinuities to control nonphysical oscillations and to ensure nonlinear total variation bounded stability. The method is first applied to the dam-breaking problem with zero rotation to validate the basic numerical implementation. Grid independence of the solutions is established and solution error is quantified by computing the $L1$ norm and comparing the estimated convergence rates to theoretical convergence rates. Stability is demonstrated subject to a Courant-Fredricks-Lewey restriction. Sloshing in a nonrotating tank with a prescribed initial water surface elevation is first investigated to demonstrate the ability of the method to capture the wave speed of traveling waves, followed by a tank undergoing sinusoidal rotation. Time histories of water surface elevation at selected locations, as well as pressure distribution on the tank walls and the corresponding moment about the tank centerline are computed and compared to experimental data and to previous computations. Finally, a limited parameter study is performed to determine the effect of varying roll angle, depth to width ratio, and forcing frequency on the resulting maximum moment about the tank centerline.

<>

## Figures

Figure 1

Tank geometry and nomenclature

Figure 2

Domain discretization nomenclature

Figure 3

Dam-break problem, grid convergence study, depth solution for Ho∕H1=0.5 at t=0.1s

Figure 4

Dam-break problem, grid convergence study, velocity solution for Ho∕H1=0.5 at t=0.1s

Figure 5

Dam-break problem, effect of slope limiting and comparison to first-order Godunov method: 1600 elements, Ho∕H1=0.5, t=0.1s

Figure 6

Time history of free oscillations in 0.5m wide rectangular tank at x=0.15m; initial average depth, h=0.025m, ω=0, initial water surface elevation inclined at 5deg

Figure 7

Pressure time history on left tank wall 6cm above bottom of tank: Alkyidiz and Unal (18) tank problem, h∕L=0.168, ω=2.0rad∕s, δo=4deg

Figure 8

Computed and experimental periodic depth solution at x∕L=0.0, L=1.2m, h=0.09m, δo=2deg, ωo=ω=2.46rad∕s

Figure 9

Computed and experimental periodic depth solution at x∕L=0.492, L=1.2m, h=0.09m, δo=2deg, ωo=ω=2.46rad∕s

Figure 10

Figure 11

Maximum computed moment about tank centerline, 0.25ωo⩽ω⩽1.75ωo, h∕L=0.1, L=1.2m, δo=2.5deg

Figure 12

Maximum computed moment, 0.05⩽h∕L⩽0.20, δo=2.5deg, ω=ωo

Figure 13

Maximum compute moment about tank centerline, 0.5deg⩽δo⩽6.0deg, h∕L=0.1, L=1.2m, ω=ωo

## Discussions

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 Proceedings Articles
Related eBook Content
Topic Collections