An automatic mesh generator for coupled 1D/2D hydrodynamic models

Abstract. Two-dimensional (2D), depth-averaged shallow water equation (SWE) models are routinely used to simulate flooding in coastal areas — areas that often include vast networks of channels and flood-control topographic features/structures, such as barrier islands and levees. Adequately resolving these features within the confines of a 2D model can be computationally expensive, which has led to coupling 2D simulation tools to less-expensive, one-dimensional (1D) models. Under certain 1D/2D coupling approaches, this introduces internal constraints that must be considered in the generation of the 2D computational mesh used. In this paper, we further develop an existing automatic unstructured mesh generation tool for SWE models, ADMESH+, to sequentially (i) identify 1D constraints from the raw input data used in the mesh generation process, namely, the digital elevation model (DEM) and land/water delineation data, (ii) distribute grid points along these internal constraints according to feature curvature and user-prescribed minimum grid spacing, and (iii) integrate these internal constraints into the 2D mesh-size-function and mesh-generation processes. The developed techniques, which include a novel approach for determining the so-called medial axis of a polygon, are described in detail and demonstrated on three test cases, including two inland watersheds with vast networks of channels and a complex estuarian system on the Texas, US coast.