Research Articles

A Finite Element Solution of the Unidimensional Shallow-Water Equation

[+] Author and Article Information
Ali Triki

Research Laboratory: Applied Fluid Mechanics,
Process Engineering and Environment,
National Engineering School of Sfax,
University of SFAX B.P.,
1173, 3038 Sfax, Tunisia
e-mail: ali.triki@enis.rnu.tn

Manuscript received July 20, 2010; final manuscript received July 25, 2012; accepted manuscript posted August 17, 2012; published online January 15, 2013. Assoc. Editor: Arif Masud.

J. Appl. Mech 80(2), 021001 (Jan 15, 2013) (8 pages) Paper No: JAM-10-1249; doi: 10.1115/1.4007424 History: Received July 20, 2010; Revised July 25, 2012; Accepted August 17, 2012

Based on the finite element method, the numerical solution of the shallow-water equation for one-dimensional (1D) unsteady flows was established. To respect the stability criteria, the time step of the method was dependent on the space step and flow velocity. This method was used to avoid the restriction due to the wave celerity variation in the computational analysis when using the method of characteristics. Furthermore, boundary conditions are deduced directly from the scheme without using characteristics equations. For the numerical solution, a general-purpose computer program, based on the finite element method (FEM), is coded in fortran to analyze the dynamic response of the open channel flow. This program is able to handle rectangular, triangular, or trapezoidal sections. Some examples solved with the finite element method are reported herein. The first involves routing a discharge hydrograph down a rectangular channel. The second example consists of routing a sudden shutoff of all flow at the downstream end of a rectangular channel. The third one deals with routing a discharge hydrograph down a trapezoidal channel. These examples are taken from the quoted literature text book. Numerical results agree well with those obtained by these authors and show that the proposed method is consistent, accurate, and highly stable in capturing discontinuities propagation in free surface flows.

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



Grahic Jump Location
Fig. 1

Schematic cross section of the channel

Grahic Jump Location
Fig. 3

Depth velocity and discharge evolution

Grahic Jump Location
Fig. 5

Comparison between the finite element method and the method of characteristics (Streeter and Wylie [4])

Grahic Jump Location
Fig. 6

Comparison between the FEM method and Streeter and Wylie [4]

Grahic Jump Location
Fig. 7

Depths and velocities evolution

Grahic Jump Location
Fig. 8

Profiles of free surface evolution

Grahic Jump Location
Fig. 9

Comparison between the FEM method and Streeter and Wylie [4]

Grahic Jump Location
Fig. 10

Depth velocity and discharge evolution




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