Large-eddy simulations (LES) above forests and cities typically constrain the simulation domain to the first 10--20\% of the Atmospheric Boundary Layer (ABL), aiming to represent the finer details of the roughness elements and sublayer. These simulations are also commonly driven by a constant pressure gradient term in the streamwise direction and zero stress at the top, resulting in an unrealistic fast decay of the total stress profile. In this study, we investigate five LES setups, including pressure and/or top-shear driven flows with and without the Coriolis force, with the aim of identifying which option best represents turbulence profiles in the atmospheric surface layer (ASL). We show that flows driven solely by pressure not only result in a fast-decaying stress profile, but also in lower velocity variances and higher velocity skewnesses. Top-shear driven flows, on the other hand, better replicate ASL statistics. Overall, we recommend, and provide setup guidance for, simulation designs that include both a large scale pressure forcing and a non-zero stress and scalar flux at the top of the domain, and that also represent the Coriolis force. Such setups retain all the forces used in typical full ABL cases and result in the best match of the profiles of various statistical moments.
more »
« less
This content will become publicly available on April 1, 2025
Analytical Model Coupling Ekman and Surface Layer Structure in Atmospheric Boundary Layer Flows
We introduce an analytical model that describes the vertical structure of Ekman boundary layer flows coupled to the Monin-Obukhov Similarity Theory (MOST) surface layer repre- sentation, which is valid for conventionally neutral (CNBL) and stable (SBL) atmospheric conditions. The model is based on a self-similar profile of horizontal stress for both CNBL and SBL flows that merges the classic 3/2 power law profile with a MOST-consistent stress profile in the surface layer. The velocity profiles are then obtained from the Ekman momentum balance equation. The same stress model is used to derive a new self-consistent Geostrophic Drag Law (GDL). We determine the ABL height (h) using an equilibrium boundary layer height model and parameterize the surface heat flux for quasi-steady SBL flows as a function of a prescribed surface temperature cooling rate. The ABL height and GDL equations can then be solved together to obtain the friction velocity (u∗) and the cross-isobaric angle (α0) as a function of known input parameters such as the Geostrophic wind speed and surface roughness (z0). We show that the model predictions agree well with simulation data from the literature and newly generated Large Eddy Simulations (LES). These results indicate that the proposed model provides an efficient and relatively accurate self-consistent approach for predicting the mean wind velocity distribution in CNBL and SBL flows.
more »
« less
- NSF-PAR ID:
- 10508221
- Publisher / Repository:
- Springer
- Date Published:
- Journal Name:
- Boundary-Layer Meteorology
- Volume:
- 190
- Issue:
- 4
- ISSN:
- 0006-8314
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Motivated by the need for accurate determination of wall shear stress from profile measurements in turbulent boundary layer flows, the total shear stress balance is analysed and reformulated using several well-established semi-empirical relations. The analysis highlights the significant effect that small pressure gradients can have on parameters deduced from data even in nominally zero pressure gradient boundary layers. Using the comprehensive shear stress balance together with the log-law equation, it is shown that friction velocity, roughness length and zero-plane displacement can be determined with only velocity and turbulent shear stress profile measurements at a single streamwise location for nominally zero pressure gradient turbulent boundary layers. Application of the proposed analysis to turbulent smooth- and rough-wall experimental data shows that the friction velocity is determined with accuracy comparable to force balances (approximately 1 %–4 %). Additionally, application to boundary layer data from previous studies provides clear evidence that the often cited discrepancy between directly measured friction velocities (e.g. using force balances) and those derived from traditional total shear stress methods is likely due to the small favourable pressure gradient imposed by a fixed cross-section facility. The proposed comprehensive shear stress analysis can account for these small pressure gradients and allows more accurate boundary layer wall shear stress or friction velocity determination using commonly available mean velocity and shear stress profile data from a single streamwise location.more » « less
-
Baroclinicity adds another layer of complexity to the much-studied barotropic atmospheric boundary layers (ABLs) by modulating the pressure gradient in height. Despite the prevalence of baroclinicity in real-world ABLs, our knowledge of the interacting effects of baroclinicity and atmospheric stability is limited. In this talk, we aim to address this knowledge gap by systematically varying baroclinicity and stability using the large-eddy simulation (LES) technique. We will present how baroclinicity alters the friction velocity, Obukhov length, shear production, ABL height, and low-level jet (LLJ) in diabatic baroclinic atmospheric flows. It will be shown that while baroclinicity significantly impacts unstable, neutral, and weakly stable ABLs, its effects reduce with increased stratification in the ABL. In strongly stratified ABLs, the LLJ height, friction velocity, and Obukhov length converge to a constant asymptote independent of the baroclinicity regime. We will demonstrate that this behavior is attributed to the strong turbulence destruction in very stable ABLs that decouples the surface from higher elevations where baroclinicity is more important. Finally, two rescaling methods in the inner and outer layers of stable baroclinic ABLs will be presented to non-dimensionalize and collapse the wind profiles in baroclinic environments. The developed reduced model for different baroclinic wind profiles will be shown against the LES results. The findings of this research elucidate the underlying physics of baroclinic diabatic ABLs and are useful for characterizing the wind profiles in weather/climate models, field measurements, and various industrial applications.more » « less
-
null (Ed.)Abstract The dynamical processes governing the seasonal cycle of the Atlantic meridional overturning circulation (AMOC) are studied using a variety of models, ranging from a simple forced Rossby wave model to an eddy-resolving ocean general circulation model. The AMOC variability is decomposed into Ekman and geostrophic transport components, which reveal that the seasonality of the AMOC is determined by both components in the extratropics and dominated by the Ekman transport in the tropics. The physics governing the seasonal fluctuations of the AMOC are explored in detail at three latitudes (26.5°N, 6°N, and 34.5°S). While the Ekman transport is directly related to zonal wind stress seasonality, the comparison between different numerical models shows that the geostrophic transport involves a complex oceanic adjustment to the wind forcing. The oceanic adjustment is further evaluated by separating the zonally integrated geostrophic transport into eastern and western boundary currents and interior flows. The results indicate that the seasonal AMOC cycle in the extratropics is controlled mainly by local boundary effects, where either the western or eastern boundary can be dominant at different latitudes, while in the northern tropics it is the interior flow and its lagged compensation by the western boundary current that determine the seasonal AMOC variability.more » « less
-
null (Ed.)Abstract Large-eddy simulations are used to investigate the influence of a horizontal frontal zone, represented by a stationary uniform background horizontal temperature gradient, on the wind- and wave-driven ocean surface boundary layers. In a frontal zone, the temperature structure, the ageostrophic mean horizontal current, and the turbulence in the ocean surface boundary layer all change with the relative angle among the wind and the front. The net heating and cooling of the boundary layer could be explained by the depth-integrated horizontal advective buoyancy flux, called the Ekman Buoyancy Flux (or the Ekman-Stokes Buoyancy Flux if wave effects are included). However, the detailed temperature profiles are also modulated by the depth-dependent advective buoyancy flux and submesoscale eddies. The surface current is deflected less (more) to the right of the wind and wave when the depth-integrated advective buoyancy flux cools (warms) the ocean surface boundary layer. Horizontal mixing is greatly enhanced by submesoscale eddies. The eddy-induced horizontal mixing is anisotropic and is stronger to the right of the wind direction. Vertical turbulent mixing depends on the superposition of the geostrophic and ageostrophic current, the depth-dependent advective buoyancy flux, and submesoscale eddies.more » « less
