Third M.I.T. Conference on Computational Fluid and Solid Mechanics June 14–17, 2005  

Three-dimensional numerical modelling of curved open channel using nonhydrostatic turbulent finite element solver for free-surface flows

Célestin Leupia,*, Edie Migliob, Mustafa S. Altinakarc
aISE-STI-LIN, Ecole Polytechnique Fédérale de Lausanne, Lausanne 1015, Switzerland bPolitecnico di Milano, Dipartimento di Matematica 'F. Brioschi', Piazza Leonardo da Vinci 32, 20133 Milano, Italy cNCCHE, The University of Mississipi, Carrier Hall Room 102 University, MS 38677, USA

  Full Text
ABSTRACT
A novel three-dimensional model for a free-surface without the conventional hydrostatic pressure assumption has been applied to simulate the curved channel flows. Using the unstructured grid the finite element solver adopts the socalled Raviart-Thomas finite element in the horizontal plane and the conventional linear galerkin finite element in the vertical direction. The fixed strata system is used in the vertical plane and this procedure allows accurate grid resolution at the bed and the free-surface. A time-marching scheme is achieved by using the fractional three-step semi-implicit method while the characteristic method is used for the computation of convective terms. A new algorithm is used for the k − ∈ turbulence solver, and the related algebraic eddy coefficients are modified to account for some anisotropic and the related secondary effects appearing in the curved open channels flows. The model is applied successfully to the three-dimensional unsteady curved open channel flow for which experimental data are available for comparison.

Keywords:  Anisotropic turbulence; Curved open channel; Secondary currents; Nonhydrostatic pressure; Lagrange-Galerkin method Raviart-Thomas finite element

* Corresponding author. Tel.: + 41 (21) 652 2507; Fax: + 41 (21) 652 3646; E-mail: Celestin.leupi@epfl.ch