Grid convergence studies for the prediction of hurricane storm surge

SUMMARY The focus of this paper is a systematic determination of the relationship between grid resolution and errors associated with computations of hurricane storm surge. A grid structure is sought that provides the spatial resolution necessary to capture pertinent storm surge physics and does not overdiscretize. A set of numerical experiments simulating storm surge generation over 14 grid discretizations of idealized domains examines the influence of grid spacing, shoreline detail, coastline resolution and characteristics of the meteorological forcing on storm surge computations. Errors associated with a given grid are estimated using a Richardson-based error estimator. Analysis of the magnitude and location of estimated errors indicates that underresolution on the continental shelf leads to significant overprediction of the primary storm surge. In deeper waters, underresolution causes smearing or damping of the inverted barometer forcing function, which in turn results in underprediction of the surge elevation. In order to maintain a specified error level throughout the duration of the storm, the highest grid resolution is required on the continental shelf and particularly in nearshore areas. The disparity of discretization requirements between deep waters and coastal regions is best met using a graded grid. Application of the graded gridding strategy to the hindcast of Hurricane Camille reinforces the necessity of using a grid that has high levels of resolution in nearshore regions and areas of complex coastal geometry. # 1998 John Wiley & Sons, Ltd. Int. J. Numer. Meth. Fluids, 26: 369‐401 (1998).

[1]  U. C. Mohanty,et al.  Simulation of storm surges using a three‐dimensional numerical model: An application to the 1977 Andhra cyclone , 1983 .

[2]  P. Roache Perspective: A Method for Uniform Reporting of Grid Refinement Studies , 1994 .

[3]  R. Flather A numerical model investigation of the storm surge of 31 January and 1 February 1953 in the North Sea , 1984 .

[4]  A. M. Davies,et al.  Quantitative skill assessment for coastal ocean models , 1995 .

[5]  A. H. Al-Rabeh,et al.  A hydrodynamic model for wind-driven and tidal circulation in the Arabian Gulf , 1990 .

[6]  W. G. Gray,et al.  Shallow water modeling in spherical coordinates: equation formulation, numerical implementation, and application , 1994 .

[7]  Meshing Requirements for Large Scale Coastal Ocean Tidal Models , 1994 .

[8]  J. Bennett,et al.  Accuracy of a finite-difference method for computing lake currents☆ , 1987 .

[9]  R. Reid,et al.  An investigation of hurricane-induced forerunner surge in the Gulf of Mexico / , 1985 .

[10]  On modelling geophysical flows having low rossby numbers , 1993 .

[11]  Joannes J. Westerink,et al.  General Spectral Computations of the Nonlinear Shallow Water Tidal Interactions within the Bight of Abaco , 1989 .

[12]  Ingemar Kinnmark,et al.  The Shallow Water Wave Equations: Formulation, Analysis and Application , 1985 .

[13]  C. Blain,et al.  Influence of Domain Size and Grid Structure on the Response Characteristics of a Hurricane Storm Surge Model. , 1995 .

[14]  P. Roache,et al.  Convergence studies with the Sandia Ocean modeling system , 1990 .

[15]  S. K. Dube,et al.  On the effect of bathymetry in numerical storm surge simulation experiments , 1983 .

[16]  William G. Gray,et al.  Progress in Surface Water Modeling , 1991 .

[17]  J. Garratt Review of Drag Coefficients over Oceans and Continents , 1977 .

[18]  L. Leslie,et al.  A storm surge model for the Australian region , 1990 .

[19]  V. Cardone,et al.  Unified Program for the Specification of Hurricane Boundary Layer Winds Over Surfaces of Specified Roughness. , 1992 .

[20]  R. Lardner,et al.  A COMPARISON OF SPATIAL GRIDS FOR NUMERICAL MODELLING OF FLOWS IN NEAR-COASTAL SEAS , 1992 .

[21]  William G. Gray,et al.  A wave equation model for finite element tidal computations , 1979 .

[22]  R. Luettich,et al.  Modelling tides in the western North Atlantic using unstructured graded grids , 1994 .

[23]  Joannes J. Westerink,et al.  Aspects of nonlinear simulations using shallow-water models based on the wave continuity equation , 1994 .

[24]  S. Dube,et al.  Numerical simulation of storm surges in Bangladesh using a bay-river coupled model , 1986 .

[25]  C. Provost,et al.  Some tests of precision for a finite element model of ocean tides , 1986 .

[26]  Joannes J. Westerink,et al.  The influence of domain size on the response characteristics of a hurricane storm surge model , 1994 .

[27]  Daniel R. Lynch,et al.  Progress in Hydrodynamic Modeling Review of U.S. Contributions, 1979-1982 (Paper 3R0064) , 1983 .

[28]  Domain and Grid Sensitivity Studies for Hurricane Storm Surge Predictions , 1994 .

[29]  Norman W. Scheffner,et al.  ADCIRC: An Advanced Three-Dimensional Circulation Model for Shelves, Coasts, and Estuaries. Report 1. Theory and Methodology of ADCIRC-2DDI and ADCIRC-3DL. , 1992 .