A parallel 3D free surface Navier-Stokes solver for high performance computing at the German Waterways Administration

This paper deals with the numerical modeling of the highly complicated interactions between topography, hydraulic structures and vessels with the two phases water and air. Standard 3D models for river hydraulics are not necessarily in the position to describe these phenomena. For the approximation, a two-phase Navier-Stokes solver is indispensable. The goal is accordingly to extend the parallel 3D two-phase flow source code NaSt3DGPF for applications in the near vicinity of hydraulic structures and vessels. The model includes all significant processes relevant to hydraulic structures at waterways. Providing the transition between water and air the Level Set Method is implemented, turbulence is modeled by using a LES concept. Regarding an automated meshing of real world CAD applications, the substantial advantages of using a Finite Difference Method based on a non-uniform structured grid will be shown. In view of an effective usage of the HPC resources, the limitations concerning the sizes of the 3D modeling domains, the possibilities of optimizing the grid spaces and a necessary coupling with 1D far field models are in the focus of interest. The validation against laboratory models and recent projects at German waterways with up to 10 million grid points will be presented.

[1]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[2]  George Bergeles,et al.  DEVELOPMENT AND ASSESSMENT OF A VARIABLE-ORDER NON-OSCILLATORY SCHEME FOR CONVECTION TERM DISCRETIZATION , 1998 .

[3]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

[4]  B. Engquist,et al.  Discretization of Dirac delta functions in level set methods , 2005 .

[5]  S. Osher,et al.  Regular Article: A PDE-Based Fast Local Level Set Method , 1999 .

[6]  A. Chorin.  Numerical Solution of the Navier-Stokes Equations* , 1989 .

[7]  B. P. Leonard,et al.  A stable and accurate convective modelling procedure based on quadratic upstream interpolation , 1990 .

[8]  P. J. Mason,et al.  On the magnitude of the subgrid-scale eddy coefficient in large-eddy simulations of turbulent channel flow , 1986, Journal of Fluid Mechanics.

[9]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[10]  M. Griebel,et al.  A PARALLEL LEVEL-SET APPROACH FOR TWO-PHASE FLOW PROBLEMS WITH SURFACE TENSION IN THREE SPACE DIMENSIONS , 2004 .

[11]  Ronald Fedkiw,et al.  Level set methods and dynamic implicit surfaces , 2002, Applied mathematical sciences.

[12]  J. Smagorinsky,et al.  GENERAL CIRCULATION EXPERIMENTS WITH THE PRIMITIVE EQUATIONS , 1963 .

[13]  S. Osher,et al.  A PDE-Based Fast Local Level Set Method 1 , 1998 .

[14]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods , 1999 .

[15]  P. Sagaut.  Large Eddy Simulation for Incompressible Flows , 2001 .

[16]  Iehisa Nezu,et al.  Turbulence in open-channel flows , 1993 .

[17]  Michael Griebel,et al.  Numerical Simulation in Fluid Dynamics: A Practical Introduction , 1997 .

[18]  禰津 家久,et al.  Turbulence in open-channel flows , 1993 .