High-resolution subgrid models: background, grid generation, and implementation
Tóm tắt
The basic idea of subgrid models is the use of available high-resolution bathymetric data at subgrid level in computations that are performed on relatively coarse grids allowing large time steps. For that purpose, an algorithm that correctly represents the precise mass balance in regions where wetting and drying occur was derived by Casulli (Int J Numer Method Fluids 60:391–408, 2009) and Casulli and Stelling (Int J Numer Method Fluids 67:441–449, 2010). Computational grid cells are permitted to be wet, partially wet, or dry, and no drying threshold is needed. Based on the subgrid technique, practical applications involving various scenarios were implemented including an operational forecast model for water level, salinity, and temperature of the Elbe Estuary in Germany. The grid generation procedure allows a detailed boundary fitting at subgrid level. The computational grid is made of flow-aligned quadrilaterals including few triangles where necessary. User-defined grid subdivision at subgrid level allows a correct representation of the volume up to measurement accuracy. Bottom friction requires a particular treatment. Based on the conveyance approach, an appropriate empirical correction was worked out. The aforementioned features make the subgrid technique very efficient, robust, and accurate. Comparison of predicted water levels with the comparatively highly resolved classical unstructured grid model shows very good agreement. The speedup in computational performance due to the use of the subgrid technique is about a factor of 20. A typical daily forecast can be carried out in less than 10 min on a standard PC-like hardware. The subgrid technique is therefore a promising framework to perform accurate temporal and spatial large-scale simulations of coastal and estuarine flow and transport processes at low computational cost.
Tài liệu tham khảo
Balzano A (1998) Evaluation of methods for numerical simulation of wetting and drying in shallow water flow models. Coast Eng 34(83):107
Bates PD (2000) Development and testing of a sub-grid-scale model for moving boundary hydrodynamic problems in shallow water. Hydrol Process 14(2073):2088
Bates PD, De Roo APJ (2000) A simple raster-based model for flood inundation simulation. J Hydrol 236(54):77
Berg P, Weismann Poulsen J (2012) Implementation details for HBM. DMI, Copenhagen. DMI technical report no. 12-11. www.dmi.dk/dmi/tr12-11.pdf
Casulli V (1990) Semi-implicit finite difference methods for the two-dimensional shallow water equations. J Comput Phys 86(56):74
Casulli V (2009) A high-resolution wetting and drying algorithm for free-surface hydrodynamics. Int J Numer Method Fluids 60:391–408
Casulli V, Cattani E (1994) Stability, accuracy and efficiency of a semi-implicit method for three-dimensional shallow water flow. Comput Math Appl 27(99):112
Casulli V, Cheng RT (1992) Semi-implicit finite difference methods for three-dimensional shallow water flow. Int J Numer Methods Fluids 15(629):648
Casulli V, Stelling GS (2010) Semi-implicit subgrid modelling of three-dimensional free-surface flows. Int J Numer Method Fluids 67:441–449
Casulli V, Walters RA (2000) An unstructured grid, three-dimensional model based on the shallow water equations. Int J Numer Methods Fluids 32(331):348
Casulli V, Zanolli P (2005) High resolution methods for multidimentional advection-diffusion problems in free surface hydrodynamics. Ocean Model 10(137):151
Doms GA (2011) Description of the non-hydrostatic regional COSMO model. COSMO Consortium. (http://www.cosmo-model.org). Accessed Sept 2011
Horritt DM, Bates PD (2001) Effects of spatial resolution on raster based models of flood flow. J Hydrol 253(239):249
Mortensen NG, Landberg L, Troen I, Petersen EL (1993) Wind atlas analysis, and application program (WASP). User guide, vol 2. Risø National Laboratory Roskilde, Denmark, p 133
Stelling GS (2012) Quadtree flood simulations with sub-grid digital elevation models. In: Proceedings of the Institution Of Civil Engineers ICE, vol 165, issue 10. Water management, pp. 567–580
Stelling GS, Kernkamp HWJ, Laguzzi MM (1998) Delft flooding system: a powerful tool for inundation assessment based upon a positive flow simulation. Hydroinformatics 449:456
Sweers HE (1976) A nomogram to estimate the heat-exchange coefficient at the air-water interface as a function of wind speed and temperature: a critical survey of some literature. J Hydrol 30(375):401
Umlauf L, Burchard H (2003) A generic length-scale equation for geophysical turbulence models. J Mar Res 61(31):235–265
Van’t Hof B, Vollebregt EAH (2005) Modelling of wetting and drying of shallow water using artificial porosity. Int J Numer Methods Fluids 48(1199):1217
Yu D, Lane SN (2006a) Urban fluvial flood modelling using a two-dimensional diffusion-wave treatment, part 1: mesh resolution effects. Hydrol Process 20(7):1541–1565
Yu D, Lane SN (2006b) Urban fluvial flood modelling using a two-dimensional diffusion-wave treatment, part 2: development of a sub-grid-scale treatment. Hydrol Process 20(7):1567–1583