Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boundary layer”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

Evaluation of Atmospheric Boundary Layer Height From Wind Profiling Radar and Slab Models and Its Responses to Seasonality of Land Cover, Subsidence, and Advection

In this study, we evaluated the effect of land cover, atmospheric subsidence, and advection on the annual dynamics of atmospheric boundary layer (ABL) height from two contrasting sites. The first site is the Walker Branch forest, a deciduous forest of temperate climate, complex topography, and cloudy summers. The second site is the Sacramento-San Joaquin River Delta, a site of Mediterranean climate, flat terrain on a local scale, and clear summers. After testing a new algorithm to calculate ABL heights from 915 MHz radar wind profilers, we evaluated a hierarchy of three slab models to recreate the diurnal and annual patterns of ABL growth. We found that the lower ABL heights in the Delta, particularly during late summer, are driven by the combined effects of increased atmospheric subsidence and marine air advection. In both sites, the annual pattern of ABL height was strongly correlated to total daily incoming radiation, and in the Delta, the annual pattern of ABL height closely followed the seasonal patterns of sensible heat flux from a mosaic of different land covers. A land composite of latent and sensible heat fluxes obtained through a meso-network of eddy covariance measurements and the ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) mission resulted in higher model skill, thus showing that land cover heterogeneity is an important driver of ABL growth. To conclude, model simulations show that in the Delta, restoring agricultural land to wetlands with large open water areas could result in a reduction of ABL height during those months with low subsidence and advection.

54 ENVIRONMENTAL SCIENCES↗

Low ozone dry deposition rates to sea ice during the MOSAiC field campaign: Implications for the Arctic boundary layer ozone budget

Dry deposition to the surface is one of the main removal pathways of tropospheric ozone (O 3 ). We quantified for the first time the impact of O 3 deposition to the Arctic sea ice on the planetary boundary layer (PBL) O 3 concentration and budget using year-round flux and concentration observations from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) campaign and simulations with a single-column atmospheric chemistry and meteorological model (SCM). Based on eddy-covariance O 3 surface flux observations, we find a median surface resistance on the order of 20,000 s m –1 , resulting in a dry deposition velocity of approximately 0.005 cm s –1 . This surface resistance is up to an order of magnitude larger than traditionally used values in many atmospheric chemistry and transport models. The SCM is able to accurately represent the yearly cycle, with maxima above 40 ppb in the winter and minima around 15 ppb at the end of summer. However, the observed springtime ozone depletion events are not captured by the SCM. In winter, the modelled PBL O 3 budget is governed by dry deposition at the surface mostly compensated by downward turbulent transport of O 3 towards the surface. Advection, which is accounted for implicitly by nudging to reanalysis data, poses a substantial, mostly negative, contribution to the simulated PBL O 3 budget in summer. During episodes with low wind speed (<5 m s –1 ) and shallow PBL (<50 m), the 7-day mean dry deposition removal rate can reach up to 1.0 ppb h –1 . Our study highlights the importance of an accurate description of dry deposition to Arctic sea ice in models to quantify the current and future O 3 sink in the Arctic, impacting the tropospheric O 3 budget, which has been modified in the last century largely due to anthropogenic activities.

54 ENVIRONMENTAL SCIENCES↗

Coastal Urban Boundary-layer Interactions with Convection (CUBIC)

The CUBIC project added three extra weather‑profiling stations across Houston during the 2021 TRACER campaign to better understand how sea breezes, urban heat, and local circulations affect air quality and help trigger storms. These observations, combined with drone data and model simulations, give a clearer picture of how the boundary layer behaves in a coastal city.

54 ENVIRONMENTAL SCIENCES↗

Learning the structure of wind: A data-driven nonlocal turbulence model for the atmospheric boundary layer

In this work, we develop a novel data-driven approach to modeling the atmospheric boundary layer. This approach leads to a nonlocal, anisotropic synthetic turbulence model which we refer to as the deep rapid distortion (DRD) model. Our approach relies on an operator regression problem that characterizes the best fitting candidate in a general family of nonlocal covariance kernels parameterized in part by a neural network. This family of covariance kernels is expressed in Fourier space and is obtained from approximate solutions to the Navier–Stokes equations at very high Reynolds numbers. Each member of the family incorporates important physical properties such as mass conservation and a realistic energy cascade. The DRD model can be calibrated with noisy data from field experiments. After calibration, the model can be used to generate synthetic turbulent velocity fields. To this end, we provide a new numerical method based on domain decomposition which delivers scalable, memory-efficient turbulence generation with the DRD model as well as others. We demonstrate the robustness of our approach with both filtered and noisy data coming from the 1968 Air Force Cambridge Research Laboratory Kansas experiments. Using these data, we witness exceptional accuracy with the DRD model, especially when compared to the International Electrotechnical Commission standard.

17 WIND ENERGY↗

Summertime Marine Boundary Layer Cloud, Thermodynamic, and Drizzle Morphology over the Eastern North Atlantic: A Four-Year Study

Abstract Summertime remote sensor and in situ data from 2016 to 2019 collected at the ARM Eastern North Atlantic (ENA) Observatory are combined with aircraft measurements from the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) campaign to quantify marine boundary layer (MBL) cloud, thermodynamic, and drizzle morphology in the region. A radar reflectivity–rainfall rate relationship ( Z – R ) is developed from aircraft data and 6-h cloud morphological regimes are identified from ENA data using a k -means algorithm driven by three independent inputs quantifying cloud thickness, drizzle intensity, and cloud field geometric complexity. Four separate MBL structural regimes representing non- or weakly drizzling single-layer stratocumulus, drizzling stratocumulus and cumulus-coupled stratocumulus, deep convection, and broken clouds embedded in northerly flow are identified. Single-layer stratocumulus is indicated when weak subtropical anticyclones are significantly west of the ENA site, and the MBL is cooler and drier than when drizzling and cumulus-coupled stratocumulus and broken clouds are observed. Drizzling and cumulus-coupled stratocumulus clouds are observed on the eastern flank of strong subtropical anticyclones in deep warm moist air masses with wind speeds exceeding 7 m s −1 and strong near-surface wind shear. Broken clouds exhibit strong wind shear near the inversion, while single-layer stratocumulus clouds have lower wind speeds and minimal shear. Net latent heat fluxes in the subcloud layer resulting from a combination of the ocean surface heat flux and evaporating drizzle average near zero over long periods in drizzling and cumulus-coupled stratocumulus. The ECMWF reanalysis version 5 (ERA5) is found to accurately represent single-layer stratocumulus properties, while producing significant discrepancies when drizzling stratocumulus and cumulus-coupled stratocumulus are observed.

54 ENVIRONMENTAL SCIENCES↗

Deep Convection‐Driven Downward Transport of Trace Gases and Aerosols From the Free Troposphere to the Boundary Layer

Deep convective clouds regulate Earth's energy and moisture budgets, yet their impact on the atmospheric boundary layer (BL) composition remains underexplored. Using long-term observations from three mid-latitude sites, we show that deep convection (DC) consistently enhances nighttime surface ozone and is often accompanied by modest increases in ultrafine particle concentrations. Within the BL, the condensational growth of these transported ultrafine particles may contribute up to 60% of total cloud condensation nuclei (CCN). Mass flux calculations suggest that short-lived convective cores (∼30 min) account for ∼2% of total vertical air mass transport relative to steady entrainment, increasing to ∼13% when the trailing stratiform regions are included. These results show that DC provides an episodic but efficient pathway linking the free troposphere and BL, influencing oxidant budgets, CCN variability, and climate forcing. Accurately representing this process in climate models may help reduce uncertainties in climate projections, under both preindustrial and present-day conditions.

54 ENVIRONMENTAL SCIENCES↗

A Moving-Wave Implementation in WRF to Study the Impact of Surface Water Waves on the Atmospheric Boundary Layer

Abstract While numerous modeling studies have focused on the interaction of ocean surface waves with the atmospheric boundary layer, most employ idealized waves that are either monochromatic or synthetically generated from a theoretical wave spectrum, and the atmospheric solvers are typically incompressible. To study wind–wave coupling in real-world scenarios, a model that can simulate both realistic meteorological and wave conditions is necessary. In this paper we describe the implementation of a moving bottom boundary condition into the Weather Research and Forecasting Model for large-eddy simulation applications. We first describe the moving bottom boundary conditions within WRF’s pressure-based vertical coordinate system. We then validate our code with idealized test cases that have analytical solutions, including flow over a monochromatic wave with and without viscosity. Finally, we present results from turbulent flows over a moving monochromatic wave with different wave ages, and demonstrate satisfactory agreement of the wave growth rate with results from the literature. We also compare atmospheric stress and wind parameters from two physically equivalent cases. The first specifies a wind moving in the same direction as a propagating wave, while the second involves a stationary wave with the wind adjusted such that the wind relative to the wave is the same as in the first case. Results indicate that the velocity and Reynolds stress profiles for the two cases match, further validating the moving bottom implementation.

17 WIND ENERGY↗

Simulating Atmospheric Boundary Layer Turbulence with Nek5000/RS

We present large-eddy-simulation (LES) modeling approaches for the simulation of atmospheric boundary layer turbulence that are of direct relevance to wind energy production. In this report, we study a GABLS benchmark problem using high-order spectral element code Nek5000/RS, which is supported under the DOE’s Exascale Computing Project (ECP) Center for Efficient Exascale Discretizations (CEED) project, targeting application simulations on various acceleration-device based exascale computing platforms [1, 2]. We demonstrate our newly developed subgrid-scale (SGS) models based on high-pass filter (HPF), mean-field eddy viscosity (MFEV), and Smagorinsky (SMG) with no-slip and traction boundary conditions, provided with low-order statistics, convergence and turbulent structure analysis. The model fidelity and scaling performance of Nek5000/RS on DOE’s leadership computing platforms in comparison to those of AMR-Wind, a block-structured second-order finite-volume code with adaptive-mesh-refinement capabilities, are discussed in [3].

17 WIND ENERGY↗

Formulation, Implementation and Validation of a 1D Boundary Layer Inflow Scheme for the QUIC Modeling System

Recent studies have highlighted the importance of accurate meteorological conditions for urban transport and dispersion calculations. In this work, we present a novel scheme to compute the meteorological input in the Quick Urban & Industrial Complex () diagnostic urban wind solver to improve the characterization of upstream wind veer and shear in the Atmospheric Boundary Layer (ABL). The new formulation is based on a coupled set of Ordinary Differential Equations (ODEs) derived from the Reynolds Averaged Navier–Stokes (RANS) equations, and is fast to compute. Building upon recent progress in modeling the idealized ABL, we include effects from surface roughness, turbulent stress, Coriolis force, buoyancy and baroclinicity. We verify the performance of the new scheme with canonical Large Eddy Simulation (LES) tests with the GPU-accelerated FastEddy"Equation missing" solver in neutral, stable, unstable and baroclinic conditions with different surface roughness. Furthermore, we evaluate QUIC calculations with and without the new inflow scheme with real data from the Urban Threat Dispersion (UTD) field experiment, which includes Lidar-based wind measurements as well as concentration observations from multiple outdoor releases of a non-reactive tracer in downtown New York City. Compared to previous inflow capabilities that were limited to a constant wind direction with height, we show that the new scheme can model wind veer in the ABL and enhance the prediction of the surface cross-isobaric angle, improving evaluation statistics of simulated concentrations paired in time and space with UTD measurements.

54 ENVIRONMENTAL SCIENCES↗

Large-eddy simulation of helical- and straight-bladed vertical-axis wind turbines in boundary layer turbulence

Here, turbulent wake flows behind helical- and straight-bladed vertical axis wind turbines (VAWTs) in boundary layer turbulence are numerically studied using the large-eddy simulation (LES) method combined with the actuator line model. Based on the LES data, systematic statistical analyses are performed to explore the effects of blade geometry on the characteristics of the turbine wake. The time-averaged velocity fields show that the helical-bladed VAWT generates a mean vertical velocity along the center of the turbine wake, which causes a vertical inclination of the turbine wake and alters the vertical gradient of the mean streamwise velocity. Consequently, the intensities of the turbulent fluctuations and Reynolds shear stresses are also affected by the helical-shaped blades when compared with those in the straight-bladed VAWT case. The LES results also show that reversing the twist direction of the helical-bladed VAWT causes the spatial patterns of the turbulent wake flow statistics to be reversed in the vertical direction. Moreover, the mass and kinetic energy transports in the turbine wakes are directly visualized using the transport tube method, and the comparison between the helical- and straight-bladed VAWT cases show significant differences in the downstream evolution of the transport tubes.

17 WIND ENERGY↗

Validation of Actuator Line and Actuator Disk Models with Filtered Lifting Line Corrections Implemented in Nalu-Wind Large Eddy Simulations of the Atmospheric Boundary Layer

Turbine generator power from simulations using Actuator Line Models and Actuator Disk Models with a Filtered Lifting Line Correction are compared to field data of a V27 turbine. Preliminary results of the wake characteristics are also presented. Turbine quantities of interest from traditional ALM and ADM with the Gaussian kernel (..epsilon..) set at the optimum value for matching power production and that resolve the kernel at all mesh sizes are also presented. The atmospheric boundary layer is simulated using Nalu-Wind, a Large Eddy Simulation code which is part of the ExaWind code suite. The effect of mesh resolution on quantities of interest is also examined.

Actuator Disk Models↗

Direct numerical simulations for hybrid rocket boundary layers: Performance modeling and scaling

This paper presents a comprehensive performance and scaling analysis of direct numerical simulations for reacting boundary layers, focusing on slab burner configurations. Using a PETSc-based finite volume CFD framework, the study evaluates the scalability and computational cost of flow, chemistry, and radiation evaluations across 2D and 3D simulations. Polymethyl methacrylate (PMMA) is the fuel with pure O 2 as the oxidizer, modeled using a detailed chemical kinetics mechanism with 113 species and 660 reactions. A ray-tracing-based radiation solver, designed for distributed memory applications, is implemented to model radiation heat transfer. Parallel scalability is analyzed for the coupled flow, chemistry, and radiation heat transfer processes. Weak and strong scaling studies are conducted on up to 15,000 computational ranks, revealing robust performance when flow cells exceed 200 per rank. Chemistry evaluations dominate the computational cost in large 3D simulations, accounting for approximately 40% of the total runtime, while flow processes contribute around 35%, and radiation solver contributions remain below 10% due to reduced evaluation frequencies. GPU accelerated chemistry evaluation, implemented with Zero-RK, demonstrates significant promise, achieving up to a 4x speedup for workloads exceeding 30,000 cells per GPU. However, diminishing returns are observed for smaller workloads due to CPU-GPU communication overhead. This study identifies key challenges, including memory bottlenecks and the effects of domain partitioning on flow scalability, while highlighting the potential of GPU-accelerated chemistry to reduce computational costs. In conclusion, these findings provide realizable run configurations for 2D, 3D, and GPU-accelerated cases, offering insights for optimizing reactive flow solvers.

CFD Scalability↗

Morning boundary layer conditions for shallow to deep convective cloud evolution during the dry season in the central Amazon

Observations of the boundary layer (BL) processes are analyzed statistically for dry seasons of 2 years and in detail, as case studies, for 4 shallow convective days (ShCu) and 4 shallow-to-deep convective days (ShDeep) using a suite of ground-based measurements from the Observation and Modeling of the Green Ocean Amazon (GoAmazon 2014/5) Experiment. The BL stages in ShDeep days, from the nighttime to the cloudy mixing layer stage, are then described in comparison with ShCu days. Atmospheric thermodynamics and dynamics, environmental profiles, and surface turbulent fluxes were employed to compare these two distinct situations for each stage of the BL evolution. Particular attention is given to the morning transition stage, in which the BL changes from stable to unstable conditions in the early morning hours. Results show that the decrease in time duration of the morning transition on ShDeep days is associated with high humidity and well-established vertical wind shear patterns. Higher humidity since nighttime not only contributes to lowering the cloud base during the rapid growth of the BL but also contributes to the balance between radiative cooling and turbulent mixing during nighttime, resulting in higher sensible heat flux in the early morning. The sensible heat flux promotes rapid growth of the well-mixed layer, thus favoring the deeper BL starting from around 08:00 LST (UTC-4 h). Under these conditions, the time duration of morning transition is used to promote convection, having an important effect on the convective BL strength and leading to the formation of shallow cumulus clouds and their subsequent evolution into deep convective clouds. Statistical analysis was used to validate the conceptual model obtained from the case studies. Despite the case-to-case variability, the statistical analyses of the processes in the BL show that the described processes are very representative of cloud evolution during the dry season.

54 ENVIRONMENTAL SCIENCES↗

Effects of Vegetation and Topography on the Boundary Layer Structure above the Amazon Forest

Observational data from two field campaigns in the Amazon forest were used to study the vertical structure of turbulence above the forest. The analysis was performed using the reduced turbulent kinetic energy (TKE) budget and its associated two-dimensional phase space. Results revealed the existence of two regions within the roughness sublayer in which the TKE budget cannot be explained by the canonical flat-terrain TKE budgets in the canopy roughness sublayer or in the lower portion of the convective ABL. Data analysis also suggested that deviations from horizontal homogeneity have a large contribution to the TKE budget. Results from LES of a model canopy over idealized topography presented similar features, leading to the conclusion that flow distortions caused by topography are responsible for the observed features in the TKE budget. These results support the conclusion that the boundary layer above the Amazon forest is strongly impacted by the gentle topography underneath.

54 ENVIRONMENTAL SCIENCES↗

Using ground-based lidar data to investigate the water–vapor budget in the daytime atmospheric boundary layer

The moisture advection term in the water–vapor budget equation is investigated with a combination of a vertically-staring water–vapor lidar and Doppler lidar systems. These instruments make it possible to get the mean profile of moisture tendency and the latent heat flux (LHF) divergence. We use data of the Land–Atmosphere Feedback Experiment (LAFE) at the Atmospheric Radiation Measurement (ARM) Program’s Southern Great Plains (SGP) site, Oklahoma, USA, collected on 30 August 2017 between 15 and 24 UTC, which corresponds to 09 to 18 LT. The lidars provide turbulence resolving profiles of moisture and vertical wind fluctuations. The LHF profile is derived from the covariance of these moisture and vertical wind fluctuations. The mean boundary layer height z i is determined from the peak of the moisture variance. The results demonstrate that the combination of two remote sensing instruments can be applied for determining the dominant water–vapor budget terms, namely moisture tendency, latent heat flux divergence and moisture advection.

Advection↗

A Vertically Resolved MSE Framework Highlights the Role of the Boundary Layer in Convective Self-Aggregation

Abstract Convective self-aggregation refers to a phenomenon in which random convection can self-organize into large-scale clusters over an ocean surface with uniform temperature in cloud-resolving models. Previous literature studies convective aggregation primarily by analyzing vertically integrated (VI) moist static energy (MSE) variance. That is the global MSE variance, including both the local MSE variance at a given altitude and the covariance of MSE anomalies between different altitudes. Here we present a vertically resolved (VR) MSE framework that focuses on the local MSE variance to study convective self-aggregation. Using a cloud-resolving simulation, we show that the development of self-aggregation is associated with an increase of local MSE variance, and that the diabatic and adiabatic generation of the MSE variance is mainly dominated by the boundary layer (BL; the lowest 2 km). The results agree with recent numerical simulation results and the available potential energy analyses showing that the BL plays a key role in the development of self-aggregation. Additionally, we find that the lower free troposphere (2–4 km) also generates significant MSE variance in the first 15 days. We further present a detailed comparison between the global and local MSE variance frameworks in their mathematical formulation and diagnostic results, highlighting their differences.

54 ENVIRONMENTAL SCIENCES↗

Boundary Layer Structures Over the Northwest Atlantic Derived From Airborne High Spectral Resolution Lidar and Dropsonde Measurements During the ACTIVATE Campaign

The Planetary Boundary Layer height (PBLH) is essential for studying PBL and ocean-atmosphere interactions. Marine PBL is usually defined to include a mixed layer (ML) and a capping inversion layer. The ML height (MLH) estimated from the measurements of aerosol backscatter by a lidar was usually compared with PBLH determined from radiosondes/dropsondes in the past, as the PBLH is usually similar to MLH in nature. However, PBLH can be much greater than MLH for decoupled PBL. Here, in this study, we evaluate the retrieved MLH from an airborne lidar (HSRL-2) by utilizing 506 co-located dropsondes during the ACTIVATE field campaign over the Northwest Atlantic from 2020 to 2022. First, we define and determine the MLH and PBLH from the temperature and humidity profiles of each dropsonde, and find that the MLH values from HSRL-2 and dropsondes agree well with each other, with a coefficient of determination of 0.66 and median difference of 18 m. In contrast, the HSRL-2 MLH data do not correspond to dropsonde-derived PBLH, with a median difference of -47 m. Therefore, we modify the current operational and automated HSRL-2 wavelet-based algorithm for PBLH retrieval, decreasing the median difference significantly to -8 m. Further data analysis indicates that these conclusions remain the same for cases with higher or lower cloud fractions, and for decoupled PBLs. These results demonstrate the potential of using HSRL-2 aerosol backscatter data to estimate both marine MLH and PBLH and suggest that lidar-derived MLH should be compared with radiosonde/dropsonde-determined MLH (not PBLH) in general.

54 ENVIRONMENTAL SCIENCES↗

A North Sea in Situ Evaluation of the Fitch Wind Farm Parameterization Within the Mellor-Yamada-Nakanishi-Niino and 3D Planetary Boundary Layer Schemes

Wind resource assessments and wind power forecasts that account for wind farm wakes are sensitive to the choice of planetary boundary layer (PBL) scheme. This work compares the one-dimensional Mellor-Yamada-Nakanishi-Niino (MYNN) PBL scheme with a three-dimensional PBL (3DPBL) scheme, evaluating predictions made with both schemes against two sets of North Sea in situ observations of wind farm wakes. The optimal PBL scheme varies based on the observations (FINO1 tower vs. aircraft), the quantity of interest (wind speed vs. turbulence kinetic energy [TKE]), and the error metric (bias, centered root mean square error [cRMSE], R2, and earth mover's distance [EMD]). Whereas 3DPBL wind speeds outperform MYNN wind speeds with respect to the cRMSE at the FINO1 site located at a single point within the turbine rotor layer, 3DPBL TKE bias is larger than MYNN TKE bias when compared to aircraft observations taken 100 m above a wind farm. Wind speeds in the aircraft region are ambiguous with regard to which PBL scheme is optimal. Aircraft MYNN wind speeds outperform 3DPBL wind speeds with respect to R2 and cRMSE but underperform with respect to bias and EMD. Future evaluations across broader temporal and spatial scales may offer further insight into model differences.

17 WIND ENERGY↗