High-Resolution Large-Eddy Simulations of Scalar Transport in Atmospheric Boundary Layer Flow over Complex Terrain

This paper presents high-resolution numerical simulations of the atmospheric flow and concentration fields accompanying scalar transport and diffusion from a point source in complex terrain. Scalar dispersion is affected not only by mean flow, but also by turbulent fluxes that affect scalar mixing, meaning that predictions of scalar transport require greater attention to the choice of numerical simulation parameters than is typically needed for mean wind field predictions. Large-eddy simulation is used in a mesoscale setting, providing modeling advantages through the use of robust turbulence models combined with the influence of synoptic flow forcing and heterogeneous land surface forcing. An Eulerian model for scalar transport and diffusion is implemented in the Advanced Regional Prediction System mesoscale code to compare scalar concentrations with data collected during field experiments conducted at Mount Tsukuba, Japan, in 1989. The simulations use horizontal grid resolution as fine as 25 m with up to eight grid nesting levels to incorporate time-dependent meteorological forcing. The results show that simulated ground concentration values contain significant errors relative to measured values because the mesoscale wind typically contains a wind direction bias of a few dozen degrees. Comparisons of simulation results with observations of arc maximum concentrations, however, lie within acceptable error bounds. In addition, this paper investigates the effects on scalar dispersion of computational mixing and lateral boundary conditions, which have received little attention in the literature-in particular, for high-resolution applications. The choice of lateral boundary condition update interval is found not to affect time-averaged quantities but to affect the scalar transport strongly. More frequent updates improve the simulated ground concentration values. In addition, results show that the computational mixing coefficient must be set to as small a value as possible to improve scalar dispersion predictions. The predicted concentration fields are compared as the horizontal grid resolution is increased from 190 m to as fine as 25 m. The difference observed in the results at these levels of grid refinement is found to be small relative to the effects of computational mixing and lateral boundary updates.