3D unsteady RANS modeling of complex hydraulic engineering flows. I: Numerical model

A general-purpose numerical method is developed for solving the full three-dimensional (3D), incompressible, unsteady Reynolds-averaged Navier-Stokes (URANS) equations in natural river reaches containing complex hydraulic structures at full-scale Reynolds numbers. The method adopts body-fitted, chimera overset grids in conjunction with a grid-embedding strategy to accurately and efficiently discretize arbitrarily complex, multiconnected flow domains. The URANS and turbulence closure equations are discretized using a second-order accurate finite-volume approach. The discrete equations are integrated in time via a dual-time-stepping, artificial compressibility method in conjunction with an efficient coupled, block-implicit, approximate factorization iterative solver. The computer code is parallelized to take full advantage of multiprocessor computer systems so that unsteady solutions on grids with 10⁶ nodes can be obtained within reasonable computational time. The power of the method is demonstrated by applying it to simulate turbulent flow at R ≃ 10⁷ in a stretch of the Chattahoochee River containing a portion of the actual bridge foundation located near Cornelia, Georgia. It is shown that the method can capture the onset of coherent vortex shedding in the vicinity of the foundation while accounting for the large-scale topographical features of the surrounding river reach. Journal of Hydraulic Engineering