Today, simplified 2D flood models are the state of the art for predicting flood wave propagation, or how floods spread across land. A team at IFC, led by UI Professor George Constantinescu, is creating 3D non-hydrostatic flood models that can more accurately simulate flood wave propagation and account for the interaction between the flood wave and large obstacles such as dams or floodplain walls. These 3D models also can be used to assess and improve the predictive capabilities of the 2D models that government agencies and consulting companies use for predicting how floods will spread and the associated risks and hazards.
Using one of the world's most powerful supercomputers - Titan, the 27-petaflop Cray XK7 at the Oak Ridge Leadership Computing Facility (OLCF) - George Constantinescu's team performed one of the first highly resolved, 3D, volume-of-fluid Reynolds-averaged Navier-Stokes (RANS) simulations of a dam break in a natural environment. The simulation allowed the team to map precise water levels for actual flood events over time. RANS is a widely used method for modelling turbulent flows.
"Flood events, like those generated by dam breaks, can be computationally very expensive to simulate", George Constantinescu stated. "Previously, there wasn't enough computer power to do these kinds of time-accurate simulations in large computational domains, but with the power of high-performance computing (HPC) and Titan, we are achieving more than was previously thought possible."
The project was supported in 2015 and 2016 within the OLCF's Director's Discretionary user programme. The OLCF, a US Department of Energy (DOE) Office of Science User Facility located at DOE's Oak Ridge National Laboratory, provides HPC resources for research and development projects to advance scientific discovery.
The team's 3D simulations showed that commonly used 2D models may inaccurately predict some aspects of flooding, such as the time over which dangerous flood levels last at certain locations and the amount of surface area flooded. Simulation results also demonstrated that 2D models may underestimate the speed at which floods spread and overestimate the time at which flood waves reach their highest point.
When the water sources that empty into a river rise simultaneously, they can trigger one or more successive flood waves. Accuracy of the 1D, 2D, or 3D flood models that track how these waves move is crucial for predicting maximum flood depth, hazardous conditions, and other variables.
"We need to know what's going to happen for situations in which a dam breaks", George Constantinescu stated. "We need to know who's going to be affected, how much time they will have to evacuate, and what else might happen to the environment as a result."
Because 2D models make simplified assumptions about some aspects of the flow, they can't account for changes in the flow, such as when the flood wave moves around large obstacles, changes rapidly in direction, or fully immerses bridge decks. The team needed a leadership-class supercomputer to run the 3D simulations and accurately capture these changes.
Using a fully non-hydrostatic 3D RANS solver, the team performed the first simulations of the hypothetical failure of two Iowa dams: the Coralville Dam in Iowa City and the Saylorville Dam in Des Moines. Each used a computational grid of about 30-50 million cells and covered a physical area of about 20 miles by 5 miles.
The team used the state-of-the-art computational fluid dynamics software STAR-CCM+. This software features a volume-of-fluid method to track the position of the water's free surface - the areas where water meets the air. In a scalability study, the team determined the peak performance of the code for the dam break simulations. The researchers used 2,500 of Titan's CPU processors for peak performance in each simulation.
The researchers also computed the same dam break test cases using a standard 2D model commonly used by IFC. When they compared the 2D results against those of the 3D simulations, they found that the 2D model underestimated how quickly the flood wave moved across land and overestimated the time at which the maximum flood occurred. This finding is important because government agencies and consulting companies use 2D shallow flow models to predict dam breaks and floods, as well as to estimate flood hazards.
"By performing these 3D simulations, we provided a huge data set that can be used to improve the accuracy of existing 2D and 1D flood models", George Constantinescu stated. "We can also examine the effectiveness of deploying flood protection structures for different flooding scenarios." The team ultimately showed that HPC can be used successfully to answer engineering questions related to the consequences of structural failure of dams and related hazards.
George Constantinescu said that as computers become faster and more powerful, simulations of full flooding events over larger physical regions will be possible. Summit, the OLCF's next-generation supercomputer that is scheduled to come online in 2018, will unearth new possibilities for George Constantinescu's research.
"Advances in numerical algorithms, automatic grid generation, and increased supercomputer power will eventually make the simulations of flood waves over large durations of time possible using Titan, and even more so with Summit", George Constantinescu stated. "Eventually, things we previously had to do by hand, such as generating a high-quality computational grid, will just be part of the typical software package."
The George Constantinescu team collaborated on this research with Professor Witold Krajewski, director of IFC and research engineer at the Iowa Institute of Hydraulic Research (IIHR)-Hydroscience and Engineering. IIHR PhD student Daniel Horna-Munoz, whose research was sponsored by IFC, was tasked with developing and testing the 3D model, running the simulations, and analyzing the data.