skip to main content


Title: Data‐Driven Prediction of Turbulent Flow Statistics Past Bridge Piers in Large‐Scale Rivers Using Convolutional Neural Networks
Abstract

Prediction of statistical properties of the turbulent flow in large‐scale rivers is essential for river flow analysis. The large‐eddy simulation (LES) provides a powerful tool for such predictions; however, it requires a very long sampling time and demands significant computing power to calculate the turbulence statistics of riverine flows. In this study, we developed encoder‐decoder convolutional neural networks (CNNs) to predict the first‐ and second‐order turbulence statistics of the turbulent flow of large‐scale meandering rivers using instantaneous LES results. We train the CNNs using a data set obtained from LES of the flood flow in a large‐scale river with three bridge piers—a training testbed. Subsequently, we employed the trained CNNs to predict the turbulence statistics of the flood flow in two different meandering rivers and bridge pier arrangements—validation testbed rivers. The CNN predictions for the validation testbed river flow were compared with the simulation results of a separately done LES to evaluate the performance of the developed CNNs. We show that the trained CNNs can successfully produce turbulence statistics of the flood flow in the large‐scale rivers, that is, the validation testbeds.

 
more » « less
Award ID(s):
1823530
NSF-PAR ID:
10445078
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Water Resources Research
Volume:
58
Issue:
1
ISSN:
0043-1397
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We develop a wall model for large-eddy simulation (LES) that takes into account various pressure-gradient effects using multi-agent reinforcement learning. The model is trained using low-Reynolds-number flow over periodic hills with agents distributed on the wall at various computational grid points. It utilizes a wall eddy-viscosity formulation as the boundary condition to apply the modeled wall shear stress. Each agent receives states based on local instantaneous flow quantities at an off-wall location, computes a reward based on the estimated wall-shear stress, and provides an action to update the wall eddy viscosity at each time step. The trained wall model is validated in wall-modeled LES of flow over periodic hills at higher Reynolds numbers, and the results show the effectiveness of the model on flow with pressure gradients. The analysis of the trained model indicates that the model is capable of distinguishing between the various pressure gradient regimes present in the flow. To further assess the robustness of the developed wall model, simulations of flow over the Boeing Gaussian bump are conducted at a Reynolds number of 2 million, based on the free-stream velocity and the bump width. The results of mean skin friction and pressure on the bump surface, as well as the velocity statistics of the flow field, are compared to those obtained from equilibrium wall model (EQWM) simulations and published experimental data sets. The developed wall model is found to successfully capture the acceleration and deceleration of the turbulent boundary layer on the bump surface, providing better predictions of skin friction near the bump peak and exhibiting comparable performance to the EQWM with respect to the wall pressure and velocity field. We also conclude that the subgrid-scale model is crucial to the accurate prediction of the flow field, in particular the prediction of separation. 
    more » « less
  2. Near-wall flow simulation remains a central challenge in aerodynamics modelling: Reynolds-averaged Navier–Stokes predictions of separated flows are often inaccurate, and large-eddy simulation (LES) can require prohibitively small near-wall mesh sizes. A deep learning (DL) closure model for LES is developed by introducing untrained neural networks into the governing equations and training in situ for incompressible flows around rectangular prisms at moderate Reynolds numbers. The DL-LES models are trained using adjoint partial differential equation (PDE) optimization methods to match, as closely as possible, direct numerical simulation (DNS) data. They are then evaluated out-of-sample – for aspect ratios, Reynolds numbers and bluff-body geometries not included in the training data – and compared with standard LES models. The DL-LES models outperform these models and are able to achieve accurate LES predictions on a relatively coarse mesh (downsampled from the DNS mesh by factors of four or eight in each Cartesian direction). We study the accuracy of the DL-LES model for predicting the drag coefficient, near-wall and far-field mean flow, and resolved Reynolds stress. A crucial challenge is that the LES quantities of interest are the steady-state flow statistics; for example, a time-averaged velocity component $\langle {u}_i\rangle (x) = \lim _{t \rightarrow \infty } ({1}/{t}) \int _0^t u_i(s,x)\, {\rm d}s$ . Calculating the steady-state flow statistics therefore requires simulating the DL-LES equations over a large number of flow times through the domain. It is a non-trivial question whether an unsteady PDE model with a functional form defined by a deep neural network can remain stable and accurate on $t \in [0, \infty )$ , especially when trained over comparatively short time intervals. Our results demonstrate that the DL-LES models are accurate and stable over long time horizons, which enables the estimation of the steady-state mean velocity, fluctuations and drag coefficient of turbulent flows around bluff bodies relevant to aerodynamics applications. 
    more » « less
  3. Abstract

    In meandering rivers, interactions between flow, sediment transport, and bed topography affect diverse processes, including bedform development and channel migration. Predicting how these interactions affect the spatial patterns and magnitudes of bed deformation in meandering rivers is essential for various river engineering and geoscience problems. Computational fluid dynamics simulations can predict river morphodynamics at fine temporal and spatial scales but have traditionally been challenged by the large scale of natural rivers. We conducted coupled large‐eddy simulation and bed morphodynamics simulations to create a unique database of hydro‐morphodynamic data sets for 42 meandering rivers with a variety of planform shapes and large‐scale geometrical features that mimic natural meanders. For each simulated river, the database includes (a) bed morphology, (b) three‐dimensional mean velocity field, and (c) bed shear stress distribution under bankfull flow conditions. The calculated morphodynamics results at dynamic equilibrium revealed the formation of scour and deposition patterns near the outer and inner banks, respectively, while the location of point bars and scour regions around the apexes of the meander bends is found to vary as a function of the radius of curvature of the bends to the width ratio. A new mechanism is proposed that explains this seemingly paradoxical finding. The high‐fidelity simulation results generated in this work provide researchers and scientists with a rich numerical database for morphodynamics and bed shear stress distributions in large‐scale meandering rivers to enable systematic investigation of the underlying phenomena and support a range of river engineering applications.

     
    more » « less
  4. Abstract

    Marine hydrokinetic (MHK) power generation systems enable harvesting energy from waterways without the need for water impoundment. A major research challenge for numerical simulations of field‐scale MHK farms stems from the large disparity in scales between the size of waterway and the energy harvesting device. We propose a large‐eddy simulation (LES) framework to perform high‐fidelity, multiresolution simulations of MHK arrays in a real‐life marine environment using a novel unstructured Cartesian flow solver coupled with a sharp‐interface immersed boundary method. The potential of the method as a powerful engineering design tool is demonstrated by applying it to simulate a 30 turbine MHK array under development in the East River in New York City. A virtual model of the MHK power plant is reconstructed from high‐resolution bathymetry measurements in the East River and the 30 turbines placed in 10 TriFrame arrangements as designed by Verdant Power. A locally refined, near the individual turbines, background unstructured Cartesian grid enables LES across a range of geometric scales of relevance spanning approximately 5 orders of magnitude. The simulated flow field is compared with a baseline LES of the flow in the East River without turbines. While velocity deficits and increased levels of turbulence kinetic energy are observed in the vicinity of the turbine wakes, away from the turbines as well as on the water surface only a small increase in mean momentum is found. Therefore, our results point to the conclusion that MHK energy harvesting from large rivers is possible without a significant disruption of the river flow.

     
    more » « less
  5. In this study, we conduct a parametric analysis to evaluate the sensitivities of wall-modeled large-eddy simulation (LES) with respect to subgrid-scale (SGS) models, mesh resolution, wall boundary conditions and mesh anisotropy. While such investigations have been conducted for attached/flat-plate flow configurations, systematic studies specifically targeting turbulent flows with separation are notably sparse. To bridge this gap, our study focuses on the flow over a two-dimensional Gaussian-shaped bump at a moderately high Reynolds number, which involves smooth-body separation of a turbulent boundary layer under pressure-gradient and surface- curvature effects. In the simulations, the no-slip condition at the wall is replaced by three different forms of boundary condition based on the thin boundary layer equations and the mean wall-shear stress from high-fidelity numerical simulation to avoid the additional complexity of modeling the wall-shear stress. Various statistics, including the mean separation bubble size, mean velocity profile, and dissipation from SGS model, are compared and analyzed. The results reveal that capturing the separation bubble strongly depends on the choice of SGS model. While simulations approach grid convergence with resolutions nearing those of wall-resolved LES meshes, above this limit, the LES predictions exhibit intricate sensitivities to mesh resolution. Furthermore, both wall boundary conditions and the anisotropy of mesh cells exert discernible impacts on the turbulent flow predictions, yet the magnitudes of these impacts vary based on the specific SGS model chosen for the simulation. 
    more » « less