Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BOUNDARY LAYER FLOW”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Advanced Turbulence Models for Large-Scale Atmospheric Boundary Layer Flows

We present high-fidelity large-eddy-simulation (LES) modeling approaches for the turbulent atmospheric boundary layer (ABL) flows. Wind energy is a prime example of an application driven by ABL. Generation of electrical energy from farms of wind turbines at night in the stable ABL is a particularly interesting situation. In this report, we consider the well-known GEWEX (Global Energy and Water Cycle Experiment) Atmospheric Boundary Layer Study (GABLS) stably stratified benchmark LES case. We use a 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. In our earlier ANL report, we demonstrated 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. In this report, we extend the range of our SGS modeling approaches in the context of the mean-field eddy viscosity (MFEV), to include the solution of an SGS turbulent kinetic energy equation (TKE). We demonstrate the model fidelity of Nek5000/RS in comparison to that of AMR-Wind, a block-structured second-order finite-volume code with adaptive-mesh-refinement capabilities, with which we studied scaling performance for both codes in comparison on DOE's leadership computing platforms.

17 WIND ENERGY↗
DOE OSTI2026

A Multiscale Simulation Study of the Effects of Nantucket Island on Atmospheric Boundary Layer Flow

Multiscale atmospheric simulations were conducted to investigate how Nantucket Island modifies the marine atmospheric boundary layer (ABL) flow, how the mechanisms responsible are represented in mesoscale versus large-eddy simulation (LES) domains, and potential impacts on wind energy production downstream, in moderate and high wind speed cases. Here, to isolate the effects of the island, each case study was simulated (i) with and without island topography included and (ii) with observed versus spatially uniform sea surface temperature (SST) values, the latter removing the influence of spatial variability in SST, while also enabling the examination of how the island effects vary with near-surface thermodynamic stability over the surrounding ocean. Island effects were found to be strongest during daytime heating, which generates convective rolls that redistribute momentum vertically over the island and downstream. Colder SSTs and weaker winds strengthen these rolls and extend their influence. Flow statistics within a representative wind turbine rotor-swept area downstream show that while vertical momentum redistribution enhances near-surface wind speeds, those aloft decrease, resulting in minimal net change to rotor-equivalent wind speed and slightly reduced vertical shear. Comparisons of mesoscale and LES domains reveal that both capture the roughness-driven near-surface slowdown, but mesoscale models underestimate vertical momentum transport compared to LES, which explicitly resolves convective rolls. Increasing LES resolution beyond the grid spacing required to begin capturing the largest convective roll structures ∼ O(100) m yields minimal additional momentum transport but increases resolved turbulence kinetic energy, relevant for stress loading.

Atmosphere-land interaction↗
DOE OSTI2023

Large eddy simulation of atmospheric boundary layer flow over complex terrain in comparison with RANS simulation and on-site measurements under neutral stability condition

Large eddy simulation (LES) of the atmospheric boundary layer (ABL) flow over complex terrain is presented with a validation using meteorological tower (met-tower) data through an improved neutral stability sampling approach. The proposed stability sampling procedure includes a condition based on the most-likely occurrence time-periods of the neutral ABL and reduces the variabilities of the conditional wind statistics calculated at the met-towers in comparison to our previous work. The ABL flow simulations are carried out over a potential wind site with a prominent hill based using the OpenFOAM-based simulator for on/off-shore wind farm applications by applying the Lagrangian-averaged scale-invariant dynamic sub-grid scale turbulence model. A low-dissipative scale-selective discretization scheme for the non-linear convection term in the LES governing equation is adopted implicitly to ensure both the second-order accuracy and bounded solution. The LES inflow is generated through a precursor method with a “tiling” approach based on the flow driving parameters obtained from a corresponding Reynolds-averaged Navier–Stokes (RANS) simulation. Overall, the averaged wind velocity profiles predicted by the LES approach at all met-tower locations show a similar tendency as the RANS results, which are also in reasonable agreement with the met-tower data. An obvious difference in wind speed standard deviation profiles is seen between LES and RANS, especially at regions downstream of the hill edge, where the LES shows under-predicted results at the highest measurement levels in comparison to the tower data. The computational costs of the LES are found to be about 20 times higher than the RANS simulations.

Energy & Fuels↗
DOE OSTI2021

Reynolds stress scaling in the near-wall region of wall-bounded flows

A new scaling is derived that yields a Reynolds-number-independent profile for all components of the Reynolds stress in the near-wall region of wall-bounded flows, including channel, pipe and boundary layer flows. The scaling demonstrates the important role played by the wall shear stress fluctuations and how the large eddies determine the Reynolds number dependence of the near-wall turbulence behaviour.

42 ENGINEERING↗
DOE OSTI2022

Alcove formation in dissolving cliffs driven by density inversion instability

We demonstrate conditions that give rise to cave-like features commonly found in dissolving cliffsides with a minimal two-phase physical model. Alcoves that are wider at the top and tapered at the bottom, with sharp-edged ceilings and sloping floors, are shown to develop on vertical solid surfaces dissolving in aqueous solvents. As evident from descending plumes, sufficiently large indentations evolve into alcoves as a result of the faster dissolution of the ceiling due to a solutal Rayleigh–Bénard density inversion instability. In contrast, defects of size below the boundary layer thickness set by the critical Rayleigh number smooth out, leading to stable planar interfaces. Furthermore, the ceiling recession rate and the alcove opening area evolution are shown to be given to first-order by the critical Rayleigh number. By tracking passive tracers in the fluid phase, we show that the alcoves are shaped by the detachment of the boundary layer flow and the appearance of a pinned vortex at the leading edge of the indentations. The attached boundary layer past the developing alcove is then found to lead to rounding of the other sides and the gradual sloping of the floor.

58 GEOSCIENCES↗
DOE OSTI2023

Numerical Investigation of Wall-Cooling Effect on Aero-Optical Distortions for Hypersonic Boundary Layer

Compressible wall-modeled large-eddy simulations of Mach 8 turbulent boundary-layer flows over a flat plate were carried out for the conditions of the hypersonic wind tunnel at Sandia National Laboratories. The simulations provide new insight into the effect of wall cooling on the aero-optical path distortions for hypersonic turbulent boundary-layer flows. Four different wall-to-recovery temperature ratios, 0.3, 0.48, 0.71, and 0.89, are considered. Despite the much lower grid resolution, the mean velocity, temperature, and resolved Reynolds stress profiles from the simulation for a temperature ratio of 0.48 are in good agreement with those from a reference direct numerical simulation. The normalized root-mean-square optical path difference obtained from the present simulations is compared with that from reference direct numerical simulations, Sandia experiments, as well as predictions obtained with a semi-analytical model by Notre Dame University. Here the present analysis focuses on the effect of wall cooling on the wall-normal density correlations, on key underlying assumptions of the aforementioned model such as the strong Reynolds analogy, and on the elevation angle effect on the optical path difference. Wall cooling is found to increase the velocity fluctuations and decrease the density fluctuations, resulting in an overall reduction of the normalized optical path distortion. Compared to the simulations, the basic strong Reynolds analogy overpredicts the temperature fluctuations for cooled walls. Also different from the strong Reynolds analogy, the velocity and temperature fluctuations are not perfectly anticorrelated. Finally, as the wall temperature is raised, the density correlation length, away from the wall but inside the boundary layer, increases significantly for beam paths tilted in the downstream direction.

42 ENGINEERING↗
DOE OSTI2022

Multi-fidelity information fusion with concatenated neural networks

Recently, computational modeling has shifted towards the use of statistical inference, deep learning, and other data-driven modeling frameworks. Although this shift in modeling holds promise in many applications like design optimization and real-time control by lowering the computational burden, training deep learning models needs a huge amount of data. This big data is not always available for scientific problems and leads to poorly generalizable data-driven models. This gap can be furnished by leveraging information from physics-based models. Exploiting prior knowledge about the problem at hand, this study puts forth a physics-guided machine learning (PGML) approach to build more tailored, effective, and efficient surrogate models. For our analysis, without losing its generalizability and modularity, we focus on the development of predictive models for laminar and turbulent boundary layer flows. In particular, we combine the self-similarity solution and power-law velocity profile (low-fidelity models) with the noisy data obtained either from experiments or computational fluid dynamics simulations (high-fidelity models) through a concatenated neural network. We illustrate how the knowledge from these simplified models results in reducing uncertainties associated with deep learning models applied to boundary layer flow prediction problems. The proposed multi-fidelity information fusion framework produces physically consistent models that attempt to achieve better generalization than data-driven models obtained purely based on data. While we demonstrate our framework for a problem relevant to fluid mechanics, its workflow and principles can be adopted for many scientific problems where empirical, analytical, or simplified models are prevalent. In line with grand demands in novel PGML principles, this work builds a bridge between extensive physics-based theories and data-driven modeling paradigms and paves the way for using hybrid physics and machine learning modeling approaches for next-generation digital twin technologies.

42 ENGINEERING↗
DOE OSTI2021

Low-Level Mesoscale and Cloud-Scale Interactions Promoting Deep Convection Initiation

Data from scanning radars, radiosondes, and vertical profilers deployed during three field campaigns are analyzed to study interactions between cloud-scale updrafts associated with initiating deep moist convection and the surrounding environment. Three cases are analyzed in which the radar networks permitted dual-Doppler wind retrievals in clear air preceding and during the onset of surface precipitation. These observations capture the evolution of (i) the mesoscale and boundary layer flow, and (ii) low-level updrafts associated with deep moist convection initiation (CI) events yielding sustained or short-lived precipitating storms. The elimination of convective inhibition did not distinguish between sustained and unsustained CI events, though the vertical distribution of convective available potential energy may have played a role. The clearest signal differentiating the initiation of sustained versus unsustained precipitating deep convection was the depth of the low-level horizontal wind convergence associated with the mesoscale flow feature triggering CI, a sharp surface wind shift boundary, or orographic upslope flow. The depth of the boundary layer relative to the height of the LFC failed to be a consistent indicator of CI potential. Widths of the earliest detectable low-level updrafts associated with sustained precipitating deep convection were ~3–5 km, larger than updrafts associated with surrounding boundary layer turbulence (~1–3 km wide). It is hypothesized that updrafts of this larger size are important for initiating cells to survive the destructive effects of buoyancy dilution via entrainment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

An Investigation of LES Wall Modeling for Rayleigh–Bénard Convection via Interpretable and Physics-Aware Feedforward Neural Networks with DNS

Abstract The traditional approach of using the Monin–Obukhov similarity theory (MOST) to model near-surface processes in large-eddy simulations (LESs) can lead to significant errors in natural convection. In this study, we propose an alternative approach based on feedforward neural networks (FNNs) trained on output from direct numerical simulation (DNS). To evaluate the performance, we conduct both a priori and a posteriori tests. In the a priori (offline) tests, we compare the statistics of the surface shear stress and heat flux, computed from filtered DNS input variables, to the stress and flux obtained from the filtered DNS. Additionally, we investigate the importance of various input features using the Shapley additive explanations value and the conditional average of the filter grid cells. In the a posteriori (online) tests, we implement the trained models in the System for Atmospheric Modeling (SAM) LES and compare the LES-generated surface shear stress and heat flux with those in the DNS. Our findings reveal that vertical velocity, a traditionally overlooked flow quantity, is one of the most important input features for determining the wall fluxes. Increasing the number of input features improves the a priori test results but does not always improve the model performance in the a posteriori tests because of the differences in input variables between the LES and DNS. Last, we show that physics-aware FNN models trained with logarithmic and scaled parameters can well extrapolate to more intense convection scenarios than in the training dataset, whereas those trained with primitive flow quantities cannot. Significance Statement The traditional near-surface turbulence model, based on a shear-dominated boundary layer flow, does not represent near-surface turbulence in natural convection. Using a feedforward neural network (FNN), we can construct a more accurate model that better represents the near-surface turbulence in various flows and reveals previously overlooked controlling factors and process interactions. Our study shows that the FNN-generated models outperform the traditional model and highlight the importance of the near-surface vertical velocity. Furthermore, the physics-aware FNN models exhibit the potential to extrapolate to convective flows of various intensities beyond the range of the training dataset, suggesting their broader applicability for more accurate modeling of near-surface turbulence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

“Gray Zone” Simulations Using a Three-Dimensional Planetary Boundary Layer Parameterization in the Weather Research and Forecasting Model

Abstract Generating accurate weather forecasts of planetary boundary layer (PBL) properties is challenging in many geographical regions, oftentimes due to complex topography or horizontal variability in, for example, land characteristics. While recent advances in high-performance computing platforms have led to an increase in the spatial resolution of numerical weather prediction (NWP) models, the horizontal gridcell spacing (Δ x ) of many regional-scale NWP models currently fall within or are beginning to approach the gray zone (i.e., Δ x ≈ 100–1000 m). At these gridcell spacings, three-dimensional (3D) effects are important, as the most energetic turbulent eddies are neither fully parameterized (as in traditional mesoscale simulations) nor fully resolved [as in traditional large-eddy simulations (LES)]. In light of this modeling challenge, we have implemented a 3D PBL parameterization for high-resolution mesoscale simulations using the Weather Research and Forecasting Model. The PBL scheme, which is based on the algebraic model developed by Mellor and Yamada, accounts for the 3D effects of turbulence by calculating explicitly the momentum, heat, and moisture flux divergences in addition to the turbulent kinetic energy. In this study, we present results from idealized simulations in the gray zone that illustrate the benefit of using a fully consistent turbulence closure framework under convective conditions. While the 3D PBL scheme reproduces the evolution of convective features more appropriately than the traditional 1D PBL scheme, we highlight the need to improve the turbulent length scale formulation. Significance Statement The spatial resolution of weather models continues to increase at a rapid rate in accordance with the enhancement of computing power. As a result, smaller-scale atmospheric features become more explicitly resolved. However, most numerical models still ignore the impact of horizontal weather variations on boundary layer flows, which becomes more important at these smaller spatial scales. To address this issue, we have implemented a new modeling approach, using fundamental principles, which accounts for horizontal variability. Our results show that including three-dimensional effects of turbulence is necessary to achieve realistic boundary layer characteristics. This novel technique may be useful for many applications including complex terrain flows, pollutant dispersion, and surface–atmosphere interaction studies.

Juliano, Timothy W.↗
DOE OSTI2022

On the Response of a Herschel–Bulkley Fluid Due to a Moving Plate

In this paper, we study the boundary-layer flow of a Herschel–Bulkley fluid due to a moving plate; this problem has been experimentally investigated by others, where the fluid was assumed to be Carbopol, which has similar properties to cement. The computational fluid dynamics finite volume method from the open-source toolbox/library OpenFOAM is used on structured quad grids to solve the mass and the linear momentum conservation equations using the solver “overInterDyMFoam” customized with non-Newtonian viscosity libraries. The governing equations are solved numerically by using regularization methods in the context of the overset meshing technique. The results indicate that there is a good comparison between the experimental data and the simulations. The boundary layer thicknesses are predicted within the uncertainties of the measurements. The simulations indicate strong sensitivities to the rheological properties of the fluid.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Analysis of an Arctic cold air outbreak during autumn and related air mass transformations forced by surface changes and advection in higher altitudes

During marine cold air outbreaks (MCAOs), cold and dry Arctic air masses are transported from the central Arctic southward across the closed sea ice and much warmer open oceans. They experience significant transformations including a rapid heating and moistening, often leading to cloud formation. While intense wintertime MCAOs have been analyzed widely, the air mass transformations during other seasons have been studied sparsely. We address this gap by investigating an MCAO case observed in September 2020. To study the transformation processes, we combine the fifth generation of atmospheric reanalyses of the global climate (ERA5), trajectory calculations, as well as shipborne and airborne measurements. In the central Arctic, observations acquired from aboard the research vessel (RV) Polarstern during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition characterized the initial state of the air mass over closed sea ice. Trajectories indicated the pathway the air mass took from RV Polarstern southward to the Fram Strait. For the first 24 h of the southbound drift, the air masses remained quasi-stationary. Then, still 15 h ahead of the marginal sea ice zone, differential advection across the boundary layer flow introduced humidity and clouds at higher altitudes between 1.5 and 2.5 km. ERA5-derived temperature and humidity tendencies indicated complex vertical interactions. Radiative cloud-top cooling, entrainment, and turbulence were significantly reduced in the lower and enhanced in the upper advected cloud layer. Eventually, the lower cloud deck dissipated. After this confluence of 2 different air masses, observations gathered by Polar 5 in Fram Strait as part of the MOSAiC Airborne observations in the Central Arctic campaign revealed cloudy, moist layers throughout the lowest 3.5 km and an increasing boundary layer height. Comparing the initial with the final state 48 h later, the largest net heating of +8 K was found close to the surface, yet the largest net moistening of +2.5 g kg -1 at an altitude of 1 km, as the initial profile was exceptionally dry here. We conclude that the observed air mass transformations were driven by the surface changes from sea ice to open ocean but additionally strongly impacted by the differential advection of clouds and moisture across the near-surface MCAO flow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Aerodynamic and production comparison of wind farms with downwind versus conventional upwind turbines

Ever-increasing turbine scales and their associated logistical challenges have reignited questions about the performance of downwind rotor configurations. A particular potential benefit of downwind rotor configurations is the farm-scale power increase that may be conferred by tilt-driven downward wake entrainment and associated wake recovery. In this work, a comprehensive aerodynamic analysis is carried out to understand the mechanisms for wake entrainment and recovery across a spectrum of velocity and inflow alignment conditions on a small, structured farm in order to understand the impact of downwind rotors on farm production. The results show that the benefits demonstrated previously in the literature for downwind-rotor farms in aligned flows are fragile, and, outside of strong farm/flow alignment conditions, power production benefits for small farms with downwind rotor configurations are significantly if not completely mitigated by misalignment effects. The work indicates that farm-scale benefits for downwind rotors must be realized either from large-scale entrainment benefits, with more exotic farm arrangements that can take advantage of the aerodynamic effects, or from beneficial fatigue impacts from entrainment of less turbulent outer boundary layer flows.

17 WIND ENERGY↗
DOE OSTI2021

Realism of Lagrangian Large Eddy Simulations Driven by Reanalysis Meteorology: Tracking a Pocket of Open Cells Under a Biomass Burning Aerosol Layer

An approach to drive Lagrangian large eddy simulation (LES) of boundary layer clouds with reanalysis data is presented and evaluated using satellite (Spinning Enhanced Visible and Infrared Imager, SEVIRI) and aircraft (Cloud-Aerosol-Radiation Interactions and Forcing, CLARIFY) measurements. The simulations follow trajectories of the boundary layer flow. They track the formation and evolution of a pocket of open cells (POC) underneath a biomass burning aerosol layer in the free troposphere. The simulations reproduce the evolution of observed stratocumulus cloud morphology, cloud optical depth, and cloud drop effective radius, and capture the timing of the cloud state transition from closed to open cells seen in the satellite imagery on the three considered trajectories. They reproduce a biomass burning aerosol layer identified by the in-situ aircraft measurements above the inversion of the POC. Entrainment of aerosol from the biomass burning layer into the POC is limited to the extent of having no impact on cloud- or boundary layer properties, in agreement with the CLARIFY observations. The two-moment bin microphysics scheme used in the simulations reproduces the in-situ cloud microphysical properties reasonably well. A two-moment bulk microphysics scheme reproduces the satellite observations in the non-precipitating closed-cell state, but overestimates liquid water path and cloud optical depth in the precipitating open-cell state due to insufficient surface precipitation. A boundary layer cold and dry bias occurring in LES can be counteracted by reducing the grid aspect ratio and by tightening the large scale wind speed nudging towards the surface.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Is There a Scalar Atmospheric Surface Layer Within a Convective Boundary Layer? Implications for Flux Measurements

Top‐down entrainment shapes the vertical gradients of sensible heat, latent heat, and CO 2 fluxes, influencing the interpretation of eddy covariance (EC) measurements in the unstable atmospheric surface layer (ASL). Using large eddy simulations for convective boundary layer flows, we demonstrate that decreased temperature gradients across the entrainment zone increase entrainment fluxes by enhancing the entrainment velocity, amplifying the asymmetry between top‐down and bottom‐up flux contributions. These changes alter scalar flux profiles, causing flux divergence or convergence and leading to the breakdown of the constant flux layer assumption (CFLA) in the ASL. As a result, EC‐measured fluxes either underestimate or overestimate “true” surface fluxes during divergence or convergence phases, contributing to energy balance non‐closure. The varying degrees of the CFLA breakdown are a fundamental cause for the non‐closure issue. These findings highlight the underappreciated role of entrainment in interpreting EC fluxes, addressing non‐closure, and understanding site‐to‐site variability in flux measurements.

eddy covariance fluxes of scalars↗
DOE OSTI2024

Swirls and scoops: Ice base melt revealed by multibeam imagery of an Antarctic ice shelf

Knowledge gaps about how the ocean melts Antarctica’s ice shelves, borne from a lack of observations, lead to large uncertainties in sea level predictions. Using high-resolution maps of the underside of Dotson Ice Shelf, West Antarctica, we reveal the imprint that ice shelf basal melting leaves on the ice. Convection and intermittent warm water intrusions form widespread terraced features through slow melting in quiescent areas, while shear-driven turbulence rapidly melts smooth, eroded topographies in outflow areas, as well as enigmatic teardrop-shaped indentations that result from boundary-layer flow rotation. Full-thickness ice fractures, with bases modified by basal melting and convective processes, are observed throughout the area. This new wealth of processes, all active under a single ice shelf, must be considered to accurately predict future Antarctic ice shelf melt.

54 ENVIRONMENTAL SCIENCES↗