Edward G. Patton

The Influence of Idealized Heterogeneity on Wet and Dry Planetary Boundary Layers Coupled to the Land Surface

Abstract This manuscript describes numerical experiments investigating the influence of 2–30-km striplike heterogeneity on wet and dry convective boundary layers coupled to the land surface. The striplike heterogeneity is shown to dramatically alter the structure of the convective boundary layer by inducing significant organized circulations that modify turbulent statistics. The impact, strength, and extent of the organized motions depend critically on the scale of the heterogeneity λ relative to the boundary layer height zi. The coupling with the land surface modifies the surface fluxes and hence the circulations resulting in some differences compared to previous studies using fixed surface forcing. Because of the coupling, surface fluxes in the middle of the patches are small compared to the patch edges. At large heterogeneity scales (λ/zi ∼18) horizontal surface-flux gradients within each patch are strong enough to counter the surface-flux gradients between wet and dry patches allowing the formation of...

The Effect of Mesh Resolution on Convective Boundary Layer Statistics and Structures Generated by Large-Eddy Simulation

AbstractA massively parallel large-eddy simulation (LES) code for planetary boundary layers (PBLs) that utilizes pseudospectral differencing in horizontal planes and solves an elliptic pressure equation is described. As an application, this code is used to examine the numerical convergence of the three-dimensional time-dependent simulations of a weakly sheared daytime convective PBL on meshes varying from 323 to 10243 grid points. Based on the variation of the second-order statistics, energy spectra, and entrainment statistics, LES solutions converge provided there is adequate separation between the energy-containing eddies and those near the filter cutoff scale. For the convective PBL studied, the majority of the low-order moment statistics (means, variances, and fluxes) become grid independent when the ratio zi/(CsΔf) > 310, where zi is the boundary layer height, Δf is the filter cutoff scale, and Cs is the Smagorinsky constant. In this regime, the spectra show clear Kolmogorov inertial subrange scaling...

Turbulence structure above a vegetation canopy

We compare the turbulence statistics of the canopy/roughness sublayer (RSL) and the inertial sublayer (ISL) above. In the RSL the turbulence is more coherent and more efficient at transporting momentum and scalars and in most ways resembles a turbulent mixing layer rather than a boundary layer. To understand these differences we analyse a large-eddy simulation of the flow above and within a vegetation canopy. The three-dimensional velocity and scalar structure of a characteristic eddy is educed by compositing, using local maxima of static pressure at the canopy top as a trigger. The characteristic eddy consists of an upstream head-down sweep-generating hairpin vortex superimposed on a downstream head-up ejection-generating hairpin. The conjunction of the sweep and ejection produces the pressure maximum between the hairpins, and this is also the location of a coherent scalar microfront. This eddy structure matches that observed in simulations of homogeneous-shear flows and channel flows by several workers and also fits with earlier field and wind-tunnel measurements in canopy flows. It is significantly different from the eddy structure educed over smooth walls by conditional sampling based only on ejections as a trigger. The characteristic eddy was also reconstructed by empirical orthogonal function (EOF) analysis, when only the dominant, sweep-generating head-down hairpin was recovered, prompting a re-evaluation of earlier results based on EOF analysis of wind-tunnel data. A phenomenological model is proposed to explain both the structure of the characteristic eddy and the key differences between turbulence in the canopy/RSL and the ISL above. This model suggests a new scaling length that can be used to collapse turbulence moments over vegetation canopies.

Large-eddy simulation of windbreak flow

A large-eddy simulation has been performed of turbulent flow around multiple windbreaks set within a wheat canopy under neutral stability conditions. The simulation is validated against a wind tunnel data set taken under similar conditions. Velocity profiles and second-order statistics are presented and compared to those found in the wind tunnel. From the numerical simulation, we discuss spatial distributions of instantaneous velocity fields and pressure statistics, which are important and telling features of the flow that are difficult to measure experimentally. We present a discussion of the momentum balance at various locations with respect to the windbreak, and similarly, we introduce the budget of a passive scalar. These discussions show the importance of the terms in each budget equation as they vary upstream and downstream of the windbreak.

Turbulent Kinetic Energy Budgets from a Large-Eddy Simulation of Airflow Above and Within a Forest Canopy.

The output of a large-eddy simulation was used to study the terms ofthe turbulent kinetic energy (TKE) budget for the air layers above andwithin a forest. The computation created a three-dimensional,time-dependent simulation of the airflow, in which the lowest third ofthe domain was occupied by drag elements and heat sources to representthe forest. Shear production was a principal source of TKE in theupper canopy, diminishing gradually above tree-top height and moresharply with depth in the canopy. The transfer of energy to subgridscales (dissipation) was the main sink in the upper part of the domainbut diminished rapidly with depth in the canopy. Removal ofresolved-scale TKE due to canopy drag was extremely important,occurring primarily in the upper half of the forest where the foliagedensity was large. Turbulent transport showed a loss at the canopytop and a gain within the canopy. These general features have beenfound elsewhere but uncertainty remains concerning the effects ofpressure transport. In the present work, pressure was calculateddirectly, allowing us to compute the pressure diffusion term. Wellabove the canopy, pressure transport was smaller than, and opposite insign to, the turbulent transport term. Near the canopy top andbelow, pressure transport acted in concert with turbulent transport toexport TKE from the region immediately above and within the uppercrown, and to provide turbulent energy for the lower parts of theforest. In combination, the transport terms accounted for over half ofthe TKE loss near the canopy top, and in the lowest two-thirds of thecanopy the transport terms were the dominant source terms in thebudget. Moreover, the pressure transport was the largest source ofturbulent kinetic energy in the lowest levels of the canopy, beingparticularly strong under convective conditions. These resultsindicate that pressure transport is important in the plant canopyturbulent kinetic energy budget, especially in the lowest portion ofthe stand, where it acts as the major driving force for turbulentmotions.

Canopy element influences on resolved- and subgrid-scale energy within a large-eddy simulation

Large-eddy simulations described here have included the effect of canopy drag through the depth of a plant canopy. Specifically, we have considered the simulation of flow through a forest. Drag forces enter the simulation with the inclusion of form and viscous drag forces in the momentum equation. In addition, we have carried a variable to represent the kinetic energy (KE) associated with the turbulent wakes behind canopy elements. Assigning typical dimensions to canopy drag elements and, hence, to the scale of wake turbulence, we have evaluated wake effects on the dissipation process and on subgrid-scale (SGS) energy arising from the cascade of resolved-scale energy. Despite the fact that the rate of conversion of resolved-scale kinetic energy to wake energy is large, and the observation that wake energy is comparable with SGS energy, an effective diffusivity for wake turbulence can be ignored when calculated from wake-scale kinetic energy and a length scale based on element dimension. Thus, it is unnecessary to carry a wake energy variable in the simulation. On the other hand, it is very important that the process of conversion of SGS energy to wake-scale energy be included in the simulation because the action of wakes is to enhance the dissipation of SGS energy. Viscous drag appears to be of little consequence, at least in our calculations where elements are considered as flat plates and fully exposed to the grid-averaged flow.

WRF-Fire: coupled weather-wildland fire modeling with the weather research and forecasting model

AbstractA wildland fire-behavior module, named WRF-Fire, was integrated into the Weather Research and Forecasting (WRF) public domain numerical weather prediction model. The fire module is a surface fire-behavior model that is two-way coupled with the atmospheric model. Near-surface winds from the atmospheric model are interpolated to a finer fire grid and are used, with fuel properties and local terrain gradients, to determine the fire’s spread rate and direction. Fuel consumption releases sensible and latent heat fluxes into the atmospheric model’s lowest layers, driving boundary layer circulations. The atmospheric model, configured in turbulence-resolving large-eddy-simulation mode, was used to explore the sensitivity of simulated fire characteristics such as perimeter shape, fire intensity, and spread rate to external factors known to influence fires, such as fuel characteristics and wind speed, and to explain how these external parameters affect the overall fire properties. Through the use of theoret...

Decaying Scalars Emitted By A Forest Canopy: A Numerical Study

A large-eddy simulation is modified to include multiple scalars emitted by a plant canopy. Each of these scalars is subjected to varying rates of chemical loss. Presented is a detailed comparison between conservedspecies and species undergoing first- and second-order chemical loss.Profiles of mean mixing ratio, mixing-ratio variance and vertical mixing-ratio flux reveal the influence of chemical reactivity. Distribution of thescalar source through the depth of the canopy is shown to locally reducethe reaction rate for second-order species. Transport efficiencies, diffusioncoefficients, and mean source heights also exhibit chemical dependencies.Budgets of mixing-ratio variance and flux elucidate the mechanisms throughwhich chemistry modifies each. Instantaneous fields show the existence ofintermittently occurring coherent structures that are thought to enhancespecies segregation.

The canopy horizontal array turbulence study

The Canopy Horizontal Array Turbulence Study (CHATS) took place in spring 2007 and is the third in the series of Horizontal Array Turbulence Study (HATS) experiments. The HATS experiments have been instrumental in testing and developing subfilterscale (SFS) models for large-eddy simulation (LES) of planetary boundary layer (PBL) turbulence. The CHATS campaign took place in a deciduous walnut orchard near Dixon, California, and was designed to examine the impacts of vegetation on SFS turbulence. Measurements were collected both prior to and following leafout to capture the impact of leaves on the turbulence, stratification, and scalar source/sink distribution. CHATS utilized crosswind arrays of fast-response instrumentation to investigate the impact of the canopy-imposed distribution of momentum extraction and scalar sources on SFS transport of momentum, energy, and three scalars. To directly test and link with PBL parameterizations of canopy-modified turbulent exchange, CHATS also included a 30-m profile ...

Vertical cross-spectral phases in neutral atmospheric flow

The cross-spectral phases between velocity components at two heights are analyzed from observations at the Hovsore test site and from the field experiments under the Cooperative Atmosphere-Surface Exchange Study in 1999. These phases represent the degree to which turbulence sensed at one height leads (or lags) in time the turbulence sensed at the other height. The phase angle of the cross-wind component is observed to be significantly greater than the phase for the along-wind component, which in turn is greater than the phase for the vertical component. The cross-wind and along-wind phases increase with stream-wise wavenumber and vertical separation distance, but there is no significant change in the phase angle of vertical velocity, which remains close to zero. The phases are also calculated using a rapid distortion theory model and large-eddy simulation. The results from the models show similar order in phasing, but the slopes of the phase curves are slightly different from the observations, especially for low wavenumbers.