How will Sea Level Rise affect the Storm Surge in Norfolk, Virginia?
In this dissertation, subgrid corrections have been added to the ADvanced CIRCulation (ADCIRC) model, a widely used, continuous-Galerkin finite-element based, shallow water flow model. This includes the full derivation of averaged governing equations, closure approximations, and subgrid implementation into the source code. Testing of this new model was first performed on 3 domains: an idealized winding channel, a tidally influenced bay in Massachusetts, and a regional storm surge model covering Calcasieu Lake in Southwestern Louisiana with forcing from Rita (2005). By pre-computing the averaged variables from high-resolution bathy/topo data sets, the model can represent hydraulic connectivity at smaller scales. This allows for a coarsening of the model and thus faster predictions of flooding, while also improving accuracy. The implementation permits changing a logic-based wetting and drying algorithm to a more desirable logic-less algorithm, and requires averaging correction factors on both an elemental and vertex basis. This new framework further increases efficiency of the model, and is general enough to be used in other Galerkin-based, finite-element, hydrodynamic models. It is shown that the flooding model with subgrid corrections can match the accuracy of the conventional model, while offering a 10 to 50 times increase in speed.
Next, higher level corrections to bottom friction and advection were incorporated into the subgrid model, and the framework was expanded and tested at the ocean-scale. It was hypothesized that by adding higher-level corrections to the model and applying them to ocean-scale domains, accurate predictions of storm surge at the smallest coastal scales can be obtained. To accomplish this, higher-level corrections were derived and implemented into the governing equations and extensive elevation and landcover data sets were curated to cover the South Atlantic Bight region of the U.S. Atlantic Coast. From there, the subgrid model was tested on an ocean-scale domain with tidal and meteorological forcing from Matthew (2016). The improvements in water level prediction accuracy due to subgrid corrections are evaluated at 218 observation locations throughout 1500 km of coast along the South Atlantic Bight. The accuracy of the subgrid model with relatively coarse spatial resolution (RMSE = 0.41 m) is better than that of a conventional model with relatively fine spatial resolution (RMSE = 0.67 m).By running on the coarsened subgrid model, we improved the accuracy over efficiency curve for the model, and as a result the computational expense of the simulation was decreased by a factor of 13.
Finally, subgrid corrections were systematically tested on a series of five ocean-scale meshes with minimum nearshore resolutions ranging from around 60 m on the highest resolution mesh to 1000 m on the coarsest mesh. This study aimed to find the mesh resolution that offered the best trade-off between accuracy and efficiency. The limitations of the subgrid model were explored and guidelines for future users were established. In all, it was found that the primary limitation to the subgrid model came from the aliasing of important flow-blocking features such as barrier islands in the coarsest resolution meshes. However, in areas without these features subgrid corrections can offer tremendous advantages while running on very coarse meshes.
The work completed in this dissertation moves the science of subgrid corrections forward by integrating the corrections into a widely used ocean-circulation and storm surge model. This work offers improvements to both hurricane storm surge forecasting and long-term design by allowing for reduced run-times and increased accuracy on coarsened numerical meshes.
Jack Voight was featured on social media in a video about his summer research in our REU program. He is running simulations of storm surge and coastal flooding as part of a project about total water levels at coastal infrastructure. Glad he is part of our team!
Johnathan starts his dissertation defense.
The inland propagation of storm surge caused by tropical cyclones depends on large and small waterways to connect the open ocean to inland bays, estuaries, and floodplains. Numerical models for storm surge require these waterways and their surrounding topography to be resolved sufficiently, which can require millions of computational cells for flooding simulations on a large (ocean scale) computational domain, leading to higher demands for computational resources and longer wall-clock times for simulations. Alternatively, the governing shallow water equations can be modified to introduce subgrid corrections that allow coarser and cheaper simulations with comparable accuracy. In this study, subgrid corrections are extended for the first time to simulations at the ocean scale. Higher-level corrections are included for bottom friction and advection, and look-up tables are optimized for large model domains. Via simulations of tides, storm surge, and coastal flooding due to Hurricane Matthew in 2016, the improvements in water level prediction accuracy due to subgrid corrections are evaluated at 218 observation locations throughout 1500 km of coast along the South Atlantic Bight. The accuracy of the subgrid model with relatively coarse spatial resolution (ERMS = 0.41 m) is better than that of a conventional model with relatively fine spatial resolution (ERMS = 0.67 m). By running on the coarsened subgrid model, we improved the accuracy over efficiency curve for the model, and as a result, the computational expense of the simulation was decreased by a factor of 13.
