Tracer Advection Using Dynamic Grid Adaptation and MM5

Abstract

A dynamic grid adaptation (DGA) technique is used to numerically simulate tracer transport at meso-and regional scales. A gridpoint redistribution scheme is designed to maximize heuristic characteristics of a "good" grid. The advective solver used in conjunction with the DGA is the multidimensional positive definite advection transport algorithm (MPDATA). The DGA results for regional tracer transport are compared against results generated using the leapfrog as well as MPDATA advection schemes with uniformly spaced, static grids. Wind fields for all tracer transport algorithms are provided by the fifth-generation Pennsylvania State University-NCAR Mesoscale Model (MM5). A mesoscale-sized test case with idealized initial condition and wind field clearly shows qualitatively and quantitatively the advantage of using the dynamic adaptive grid, which is a marked reduction in numerical error. These results are further corroborated by more realistic test cases that used NCEP-NCAR reanalysis data from 6-11 March 1992 to set initial and boundary conditions for (i) a mesoscale-sized, 24-h simulation with an idealized initial tracer field, and (ii) a regional, 5-day simulation with water vapor field initialized from the reanalysis data but then treated as a passive tracer. A result of interest is that MPDATA substantially outperforms the leapfrog method with fourth-order artificial dissipation (central to MM5) in all of our test cases. We conclude that with dynamic grid adaptation, results with approximately the same accuracy as a uniform grid may be obtained using only a quarter of the grid points of the uniform grid MPDATA simulations. Compared to results generated using the leapfrog method on a uniform grid, the DGA does even better. ABSTRACT A dynamic grid adaptation (DGA) technique is used to numerically simulate tracer transport at mesoand regional scales. A gridpoint redistribution scheme is designed to maximize heuristic characteristics of a "good" grid. The advective solver used in conjunction with the DGA is the multidimensional positive definite advection transport algorithm (MPDATA). The DGA results for regional tracer transport are compared against results generated using the leapfrog as well as MPDATA advection schemes with uniformly spaced, static grids. Wind fields for all tracer transport algorithms are provided by the fifthgeneration Pennsylvania State University-NCAR Mesoscale Model (MM5). A mesoscale-sized test case with idealized initial condition and wind field clearly shows qualitatively and quantitatively the advantage of using the dynamic adaptive grid, which is a marked reduction in numerical error. These results are further corroborated by more realistic test cases that used NCEP-NCAR reanalysis data from 6-11 March 1992 to set initial and boundary conditions for (i) a mesoscale-sized, 24-h simulation with an idealized initial tracer field, and (ii) a regional, 5-day simulation with water vapor field initialized from the reanalysis data but then treated as a passive tracer. A result of int...

EULAG, a computational model for multiscale flows

“…No single, simple GA mechanism appears to work well all the time. For example, grid adaptation is relatively straightforward in more elementary applications such as in flows past ''irregular" boundaries [63,15,62,16,35,2,54]; free surface flows [34,22,65,33]; or in tracer transport [7,14,2,21]. In these examples focusing the GA on boundary shape or targeting various combinations of tracer magnitude, gradient and/or second derivatives of the tracer field have all produced good results.…”

MPDATA and grid adaptivity in geophysical fluid flow models

“…13 To see this, set R ≡ 0 in (14), and apply (12) to rewrite the Christo el symbols as derivatives of the Jacobian. The ÿrst and middle terms can be reorganized to eliminate derivatives with respect to x i and x n , respectively.…”

“…Such methods include nested grids, grid point insertion=deletion, and redistribution methods, among others [12]. In this study, the focus is on the redistribution of continuous coordinates.…”

“…Grid adaptivity describes a wide variety of methods for locally increasing grid resolution in computational models in regions of interest. Such methods include nested grids, grid point insertion=deletion, and redistribution methods, among others [12]. In this study, the focus is on the redistribution of continuous coordinates.…”

A Cloud-Resolving Model with an Adaptive Vertical Grid for Boundary Layer Clouds

“…Much of the research in this area has focused on the application of adaptive horizontal grids for global or large regional-scale modeling (e.g., Bacon et al 2007; Barros and Garcia 2004;Iselin et al 2005;Jablonowski et al 2006;Lauter et al 2007). While atmospheric models more commonly use fixed grids than adaptive grids, the potential benefits of adaptive grids have attracted some attention.…”