The Richardson's Law in Large-Eddy Simulations of Boundary Layer flows

Relative dispersion in a neutrally stratified planetary boundary layer (PBL) is investigated by means of Large-Eddy Simulations (LES). Despite the small extension of the inertial range of scales in the simulated PBL, our Lagrangian statistics turns out to be compatible with the Richardson $t^3$ law for the average of square particle separation. This emerges from the application of nonstandard methods of analysis through which a precise measure of the Richardson constant was also possible. Its values is estimated as $C_2\sim 0.5$ in close agreement with recent experiments and three-dimensional direct numerical simulations.Comment: 15 LaTex pages, 4 PS figure

Interactions between downslope flows and a developing cold-air pool

A numerical model has been used to characterize the development of a region of enhanced cooling in an alpine valley with a width of order (Formula presented.) km, under decoupled stable conditions. The region of enhanced cooling develops largely as a region of relatively dry air which partitions the valley atmosphere dynamics into two volumes, with airflow partially trapped within the valley by a developing elevated inversion. Complex interactions between the region of enhanced cooling and the downslope flows are quantified. The cooling within the region of enhanced cooling and the elevated inversion is almost equally partitioned between radiative and dynamic effects. By the end of the simulation, the different valley atmospheric regions approach a state of thermal equilibrium with one another, though this cannot be said of the valley atmosphere and its external environment.Peer reviewe

Evaluation of large-Eddy simulations via observations of nocturnal marine stratocumulus

Data from the first research flight (RF01) of the second Dynamics and Chemistry of Marine Stratocumulus (DYCOMS-II) field study are used to evaluate the fidelity with which large-eddy simulations (LESs) can represent the turbulent structure of stratocumulus-topped boundary layers. The initial data and forcings for this case placed it in an interesting part of parameter space, near the boundary where cloud-top mixing is thought to render the cloud layer unstable on the one hand, or tending toward a decoupled structure on the other hand. The basis of this evaluation consists of sixteen 4-h simulations from 10 modeling centers over grids whose vertical spacing was 5 m at the cloud-top interface and whose horizontal spacing was 35 m. Extensive sensitivity studies of both the configuration of the case and the numerical setup also enhanced the analysis. Overall it was found that (i) if efforts are made to reduce spurious mixing at cloud top, either by refining the vertical grid or limiting the effects of the subgrid model in this region, then the observed turbulent and thermodynamic structure of the layer can be reproduced with some fidelity; (ii) the base, or native configuration of most simulations greatly overestimated mixing at cloud top, tending toward a decoupled layer in which cloud liquid water path and turbulent intensities were grossly underestimated; (iii) the sensitivity of the simulations to the representation of mixing at cloud top is, to a certain extent, amplified by particulars of this case. Overall the results suggest that the use of LESs to map out the behavior of the stratocumulus-topped boundary layer in this interesting region of parameter space requires a more compelling representation of processes at cloud top. In the absence of significant leaps in the understanding of subgrid-scale (SGS) physics, such a representation can only be achieved by a significant refinement in resolution—a refinement that, while conceivable given existing resources, is probably still beyond the reach of most centers

Comparison of two- and three-dimensional simulations of miscible Rayleigh-Taylor instability

A comparison of two-dimensional and three-dimensional high-resolution numerical large-eddy simulations of planar, miscible Rayleigh-Taylor instability flows are presented. The resolution of the three-dimensional simulation is sufficient to attain a fully turbulent state. A number of different statistics from the mixing region (e.g., growth rates, PDFs, mixedness measures, and spectra) are used to demonstrate that two-dimensional flow simulations differ substantially from the three-dimensional one. It is found that the two-dimensional flow grows more quickly than its three-dimensional counterpart at late times, develops larger structures, and is much less well mixed. These findings are consistent with the concept of inverse cascade in two-dimensional flow, as well as the influence of a reduced effective Atwood number on miscible flow

Pollutant dispersion in a developing valley cold-air pool

Pollutants are trapped and accumulate within cold-air pools, thereby affecting air quality. A numerical model is used to quantify the role of cold-air-pooling processes in the dispersion of air pollution in a developing cold-air pool within an alpine valley under decoupled stable conditions. Results indicate that the negatively buoyant downslope flows transport and mix pollutants into the valley to depths that depend on the temperature deficit of the flow and the ambient temperature structure inside the valley. Along the slopes, pollutants are generally entrained above the cold-air pool and detrained within the cold-air pool, largely above the ground-based inversion layer. The ability of the cold-air pool to dilute pollutants is quantified. The analysis shows that the downslope flows fill the valley with air from above, which is then largely trapped within the cold-air pool, and that dilution depends on where the pollutants are emitted with respect to the positions of the top of the ground-based inversion layer and cold-air pool, and on the slope wind speeds. Over the lower part of the slopes, the cold-air-pool-averaged concentrations are proportional to the slope wind speeds where the pollutants are emitted, and diminish as the cold-air pool deepens. Pollutants emitted within the ground-based inversion layer are largely trapped there. Pollutants emitted farther up the slopes detrain within the cold-air pool above the ground-based inversion layer, although some fraction, increasing with distance from the top of the slopes, penetrates into the ground-based inversion layer.Peer reviewe