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.

SourceResearch connectionFrom yearThrough year

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

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Numerical algorithms for water waves with background flow over obstacles and topography

Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

Ambrose, David M.↗

DOE OSTI2023

A New Paradigm for Observing and Modeling of Air-Sea Interactions to Advance Earth System Prediction

The protection of people, property, and environmental resources from extreme weather, seasonal patterns, and climate change drives the need for predictions of weather, ocean, and climate patterns that have skill and value at timescales longer than traditional 1-10-day forecasts, including outlooks spanning weeks to decades. Advancing Earth System Prediction (ESP) skill at this range of timescales requires improved observations, understanding, and modeling of the processes in the ocean boundary layer, the atmospheric boundary layer, and their interface. A new way of referring to this coupled system is the Air-Sea Transition Zone (ASTZ). The report that follows is framed by the paradigm that the ASTZ is a single entity that regulates the flow of energy and matter between the ocean and the atmosphere. The ASTZ is thus the medium through which the ocean and atmosphere respond to and influence one another across their often disparate scales of variability. ASTZ modeling, observing, and understanding needs are particularly acute because very few measurements exist over oceans, and even fewer span the entire ASTZ, even though oceans cover 70% of Earth’s surface and are the source of most of the rain and snow that falls on both the land and the oceans.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Grid resolution requirement for resolving rare and high intensity wall-shear stress events in direct numerical simulations

Turbulent signals are intermittent with large instantaneous fluctuations. Such large fluctuations lead to small Kolmogorov scales that are hard to resolve in numerical simulations [P. K. Yeung, K. R. Sreenivasan, and S. B. Pope, Effects of finite spatial and temporal resolution in direct numerical simulations of incompressible isotropic turbulence, Phys. Rev. Fluids 3, 064603 (2018)]. Here, the present paper follows the above basic logic, but instead of dissipation events in isotropic turbulence, we study wall-shear stress events in plane channel flow. Wall-shear stress fluctuations are increasingly more intermittent as the Reynolds number increases. Hence, one has to employ higher grid resolutions as the Reynolds number increases in order to resolve a given percentage of wall-shear stress events. The objective of this paper is to quantify effects of the grid resolutions on the rare and high intensity wall-shear stress events. We find that the standard grid resolution resolves about 99% of the wall-shear stress events at Re τ = 180. A slightly higher grid resolution has to be employed in order to resolve 99% of the wall-shear stress events at higher Reynolds numbers, and if the standard grid resolution is used for, e.g., a Re τ = 10000 channel flow, one resolves about 90%–95% wall-shear stress events.

42 ENGINEERING↗

DOE OSTI2022

ExaWind: Opportunity for Discovery

The predictive simulation of the dynamics of modern wind turbines and wind farms is a high-performance-computing (HPC) grand challenge. Modern wind turbines are the largest rotating machines in the world, with rotor diameters exceeding 200 meters, and with heights well into the atmospheric boundary layer. To address this grand challenge, the U.S. Department of Energy (DOE) Wind Energy Technologies Office and the Exascale Computing Project have been supporting the creation of the ExaWind modeling and simulation environment since 2016. ExaWind is composed of the incompressible-flow computational-fluid-dynamics (CFD) solvers AMR-Wind and Nalu-Wind and the wind-turbine-dynamics solver OpenFAST. In our hybrid CFD approach, Nalu-Wind unstructured-grid models resolve the near-body flow around the turbine; those models are embedded in, and two-way coupled to, an AMR-Wind structured-grid mesh through overset meshes. In this talk I give an overview of ExaWind software stack and describe how it is being used to address important problems in wind energy.

computational fluid dynamics↗

DOE OSTI2022

Relevance of E × B drifts for particle and heat transport in divertors

Radial electric fields up to ~4 kV m –1 are observed in the boundary between the private flux region (PFR) and the scrape-off layer (SOL) driving E × B drifts between the inner and outer targets at speeds up to 2.8 km s –1 in the Tokamak à configuration variable divertor. The resulting E × B fluxes, located in a narrow region ($\Delta {\rho _\Psi } < 0.012$ in normalized radius or $\Delta $ R – R sep <4 mm mapped to the outer midplane) are equivalent to around 20% of the total heat and particle flux to the divertor targets (inner + outer). At the peak E r , the E × B poloidal transport is equivalent to parallel flows with M ∥ ~ 3. In the snowflake divertor with a second X-point in the outer SOL, the drifts in the PFR-SOL boundary were equivalent to around 30% of the total heat and particle flux to the divertor targets and cover a region ~50% wider than in the single null ($\Delta {\rho _\Psi }$ ~ 0.018, $\Delta $R – R sep ~ 6 mm). Furthermore, the location of the PFR-SOL boundary drift shifts radially in the E ∥ × B direction when reversing the toroidal field direction. Peaks in density and electron pressure have been identified near the primary X-point along with large gradients in density, temperature, and potential, the latter resulting in a local electric field ~2.7 kV m –1 which drives a drift (1.9 km s –1 ) upwards towards the closed flux surfaces. Floating potential (V f ) magnitudes up to 75 V (~2 kT e ) were measured, indicating that V f and parallel currents should not be neglected when estimating plasma potential.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

A note on thermal history kernel for unsteady heat transfer of a spherical particle

When a particle is subjected to an unsteady ambient flow, in terms of either time-dependent relative velocity or time-dependent temperature difference, the net heat transfer from the particle cannot be calculated based on the quasi-steady heat transfer correlation alone. Due to unsteady evolution of the thermal boundary layer, there is also a history contribution to heat transfer. The history contribution to heat transfer is expressed as a convolution integral of past evolution of temperature difference between the particle and the surrounding. While Basset history force and its finite Reynolds number extension have been well studied, similar understanding of unsteady heat transfer and thermal history kernel is lacking. Here, we use existing particle-resolved simulation results to develop a finite Peclet number thermal history kernel, which when used with the convolution integral is demonstrated to accurately predict unsteady heat transfer over a range of Peclet numbers and particle-to-fluid heat capacity ratio.

42 ENGINEERING↗

DOE OSTI2024

Inverse-Velocity Transformation Wall Model for Reacting Turbulent Hypersonic Boundary Layers

The present study builds on prior work by taking advantage of the novel framework proposed by Griffin et al. (hereafter referred to as the GFM) as a baseline. The model is progressively extended to multicomponent reacting mixtures, accounting for differential diffusion and finite-rate chemistry in a similar fashion to Di Renzo & Urzay and Di Renzo et al. The accuracy of the present approach, as well the prior model of Di Renzo & Urzay, is assessed in an a priori sense for the first time in a turbulent reacting boundary layer, using the boundary-layer data of Williams et al. Five species are included in the present analysis, i.e., N s = 5, namely N 2 , O 2 , NO, N and O, a neutral mixture most representative of dissociation/recombination phenomena for temperatures below 6000 K. A schematic of a flow over a wedge representative of the described configuration is presented in Figure 1. The brief is organized as follows: In Section 2, the wall-model equations and the computational framework are presented. In Section 3, the a priori results of the proposed model are described and compared to the extended EWM. Finally, in Section 4, some conclusions are offered.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Morphological and nanomechanical changes in tungsten in high heat flux conditions

Abstract Morphological and nanomechanical alteration of tungsten in extreme environments, like those in edge localized modes in nuclear fusion environments, up to 46.3 GWm −2 heat fluxes were experimentally simulated using electrothermal plasma. Surface and subsurface damage to the tungsten is seen mainly in the form of pore formation, cracks, and resolidified melt instabilities. Mirco voids, rosette-type microfeatures, core-shell structure, particle enrichment, and submicron channels all manifest in the damaged subsurface. The formation of voids in the subsurface was determined to originate from the ductile fracture of hot tungsten by plastic flow but not developed to cracking. The voids were preferentially settled in grain boundaries, interfaces. The directionality of elongated voids and grains is biased to the heat flow vector or plasma pathway, which is the likely consequence of the thermally driven grain growth and sliding in the high-temperature conditions. The presence of a border between the transient layer and heat-affected zone is observed and attributed to plasma shock and thermal spallation of fractural tungsten at high temperature. Plasma peening-like hardening effects in tungsten were observed in the range of 22.7–46.3 GWm −2 but least in the case of the lowest heat flux, 12.5 GWm −2 .

36 MATERIALS SCIENCE↗

DOE OSTI2024

Sulfuric Acid New Particle Formation Field Campaign Report

This campaign aimed to measure gaseous sulfuric acid concentrations at the surface level using a newly developed sulfuric acid dimethylamine-reactive condensation particle counter (SAD-RCPC). Sulfuric acid is a key nucleation precursor and is the likely driving compound for new particle formation at the Atmospheric Radiation Measurement (ARM) user facility’s Southern Great Plains (SGP) observatory. The SAD-RCPC will be deployed in the future at ARM SGP (AFC010010: "Vertically Resolved Nucleation Precursors at SGP") scheduled for 2024. The campaign was in conjunction with Eleanor Browne's campaign (AFC010075: Boundary-layer gradients in new particle formation) and took place from 05/7/23 to 5/21/23. The measured sulfuric acid concentrations by the SAD-RCPC will be compared to concentrations measured by Browne's chemical ionization mass spectrometer (CIMS). The SAD-RCPC consisted of a glass flow reactor coupled to a 1-nm diethylene glycol condensation particle counter (DEG CPC). Sampled air was mixed with a high concentration of dimethylamine vapor (~1 ppbv). Sulfuric acid within the sampled gas nucleated with dimethylamine to form detectable 1-nm particles. The measured particle concentration was converted to sulfuric acid concentration using nucleation reaction kinetics. The pulse height distribution of particles detected by the DEG CPC was observed using a multi-channel analyzer. Large pre-existing particles were separated from the freshly formed particles in the SAD-RCPC from the pulse height distribution (PHD). During the short campaign, seven new particle formation events were observed. The SAD-RCPC detected an increase in sulfuric acid concentration, but this also coincided with an increase in pre-existing <3 nm particle concentrations. In addition, Browne's CIMS detect sulfuric acid concentrations in the range of 10 5 -10 6 cm -3 .

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Implications of Shear and Thermal Stratification on Wind Turbine Tip-Vortex Stability

The interaction between wind turbines in a wind farm through their wakes is a phenomenon that has been studied for decades and is still relevant today. Turbines clustered together in arrays will often operate in the wake of other upstream turbines which may lead to significant power losses and fatigue loads. For modern large-scale wind turbines, the mean shear velocity profile and thermal stratification are major components of the atmospheric boundary layer so it is important to understand their impact on near-wake development. Additionally, veer is present due to the rotation of the Earth. The impact of shear, thermal stratification and veer on the stable wake length of turbines with a dynamic control strategy is studied numerically in this work using a suite of highly resolved large-eddy simulations. Instantaneous flow fields are extracted from the simulations and used to conduct proper orthogonal decomposition (POD) and compute the mean kinetic energy fluxes by different POD modes to better understand the tip-vortex instability mechanisms. Our findings show that the dynamic pitch control scheme is able to shorten the stable wake length to about 1.5R in uniform flow. Shear can significantly affect the break up of wind turbine tip-vortices as well as the shape and stable length of the wake, whereas thermal stratification seems to only have limited contribution to the spatial development of the near-wake field. Veer causes the wake boundary to skew but has a limited impact on the wake length.

length thermal stratification↗

DOE OSTI2023

Doppler Lidar Measurements of Wind Variability and LLJ Properties in Central Oklahoma during the August 2017 Land–Atmosphere Feedback Experiment

Abstract Low-level jets (LLJs) are an important nocturnal source of wind energy in the U.S. Great Plains. An August 2017 lidar-based field-measurement campaign [the Land–Atmosphere Feedback Experiment (LAFE)] studied LLJs over the central SGP site in Oklahoma and found nearly equal occurrences of the usual southerly jets and postfrontal northeasterly jets—typically rare during this season—for an opportunity to compare the two types of LLJs during this month. Southerly winds were stronger than the northeasterlies by more than 4 m s −1 on average, reflecting a significantly higher frequency of winds stronger than 12 m s −1 . The analysis of this dataset has been expanded to other SGP Doppler-lidar sites to quantify the variability of winds and LLJ properties between sites of different land use. Geographic variations of winds over the study area were noted: on southerly wind nights, the winds blew stronger at the highest, westernmost sites by 2 m s −1 , whereas on the northeasterly flow nights, the easternmost sites had the strongest wind speeds. Lidar measurements at 5 sites during August 2017, contrasted to the 2016–21 summertime data, revealed unusual wind and LLJ conditions. Temporal hodographs using hourly averaged winds at multiple heights revealed unorganized behavior in the turbulent stable boundary layer (SBL) below the jet nose. Above the nose, some nights showed veering qualitatively similar to inertial oscillation (IO) behavior, but at amplitudes much smaller than expected for an IO, whereas other nights showed little veering. Vertical hodographs had a linear shape in the SBL, indicating little directional shear there, and veering above, resulting in a hook-shaped hodograph with height. Significance Statement Doppler-lidar measurements at the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) were used to quantify the variability of winds and low-level jet (LLJ) properties between five sites of different land use and wind regimes across this area. Knowledge of wind and LLJ structure and dynamics is important for many applications, and strong southerly LLJ winds at night are an important resource for wind energy. The analysis of multiyear (2016–21) summertime LLJ parameters provided insight into the LLJ climatology in this part of the Great Plains.

Meteorology & Atmospheric Sciences↗

DOE OSTI2024

Simulated Chain Flow between the Augusta and Columbia Metropolitan Areas

Urban heat islands (UHIs) are enhanced temperatures in urban regions (Oke, 1982). • UHIs are studied fairly well for larger cities, but there is a knowledge gap in the extent of UHI effects in smaller or moderately-sized cities such as Augusta, GA. • Four different types of UHI: • Surface • Canopy Layer • Hydrological • Boundary Layer

Wermter, Joseph↗

DOE OSTI2023

Evaluating Scaling Effects of Steam Condensation Physics-based Model for Small Modular Reactor

This study focused on assessing condensation heat transfer (CHT) models for enhancing the passive safety of small modular reactors (SMRs). It also encompassed a review of the theoretical and semi-theoretical CHT models and correlations. Condensation models such as the degradation factor, heat and mass transfer analogy (HMTA), and boundary layer models are used to study advanced nuclear reactor PCCS. The previous CHT models were developed for small tube geometries, which differ from the containment of SMR. This study examined the scaled relations of a physics-based CHT model for analyzing laminar and turbulent film condensation downward flow conditions inside scaled tube geometries.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2023

Presentation: Evaluating scaling effects of steam condensation physics-based model for small modular reactors safety

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2025

Spectral distortions to momentum and scalar exchanges by non-turbulent motion and patchy landscape variability

Modifications to the spectra of turbulent velocity and scalars and co-spectra of vertical fluxes of momentum and scalars due to patchy landscape heterogeneity and non-stationarity are explored for a Mediterranean ecosystem. About 9 months of high frequency measurements of the three velocity components, water vapor concentration, carbon dioxide concentration, and air temperature were analyzed for different seasons (spring/summer) and prevalent wind directions (southeast/northwest). The two wind directions sampled a contrast of clumped and patchy landscape comprised of olive trees (southeast) and wall bounded flow disturbed by the presence of few upwind trees (northwest). The measured spectra and co-spectra were also compared to theoretical scaling forms from stationary, planar homogeneous flow, in the absence of subsidence as derived from the Kansas experiment. To assess the role of low frequency non-turbulent motion on the spectral and co-spectral content, a 5-min Fourier cutoff was introduced and the analysis was limited to near-neutral conditions where the boundary layer depth is shallow compared to its unstable counterpart. It was shown that the velocity statistics were not appreciably impacted by the low-frequency motion causing non-stationarity. Moreover, the turbulent scalar fluxes were also shown not to be significantly impacted by such low frequency motion. The scalar variances were impacted, especially the water vapor variance and its concomitant spectral shape. When the non-turbulent motion was filtered, the scalar spectra at low wavenumbers followed expectations from the so-called attached eddy hypothesis (i.e. exhibited a $k^{-1}_x$ scaling with $k_x$ defining the longitudinal wavenumber) applicable for near-neutral conditions. For momentum co-spectra, the canonical shapes from the Kansas experiment appear to describe well the measurements here and in both dominant directions and seasons with some adjustment to the integral time scales based on wind direction. For the scalar co-spectra, deviations from the Kansas experiment were prevalent. The most noticeable and surprising deviations were their slow decay with increased sampling frequency at inertial subrange scales. This slow decay was shown not to contribute appreciably to the overall scalar fluxes. At those fine scales, predictions from local isotropy were expected to hold. The scalar co-spectral deviations from local isotropy were then discussed using a simplified co-spectral budget model where scalar–scalar co-spectra naturally emerged and the interplay between landscape heterogeneity and a scale-dependent pressure-scalar de-correlation time was postulated. It is also envisaged that the findings here offer a preliminary template for analyzing eddy-covariance data in situations that deviate from ideal conditions, especially regarding low-frequency modulations of scalar spectra and vertical scalar flux co-spectra.

Canopy turbulence↗

DOE OSTI2023

Exploring the Meteorological Impacts of Surface and Rooftop Heat Mitigation Strategies Over a Tropical City

Different heat mitigation technologies have been developed to improve the thermal environment in cities. However, the regional impacts of such technologies, especially in the context of a tropical city, remain unclear. The deployment of heat mitigation technologies at city-scale can change the radiation balance, advective flow, and energy balance between urban areas and the overlying atmosphere. We used the mesoscale Weather Research and Forecasting model coupled with a physically based single-layer urban canopy model to assess the impacts of five different heat mitigation technologies on surface energy balance, standard surface meteorological fields, and planetary boundary layer (PBL) dynamics for premonsoon typical hot summer days over a tropical coastal city in the month of April in 2018, 2019, and 2020. Results indicate that the regional impacts of cool materials (CMs), super-cool broadband radiative coolers, green roofs (GRs), vegetation fraction change, and a combination of CMs and GRs (i.e., “Cool city (CC)”) on the lower atmosphere are different at diurnal scale. Results showed that super-cool materials have the maximum potential of ambient temperature reduction of 1.6°C during peak hour (14:00 LT) compared to other technologies in the study. During the daytime hours, the PBL height was considerably lower than the reference scenario with no implementation of strategies by 700 m for super-cool materials and 500 m for both CMs and CC cases; however, the green roofing system underwent nominal changes over the urban area. During the nighttime hours, the PBL height increased by CMs and the CC strategies compared to the reference scenario, but minimal changes were evident for super-cool materials. The changes of temperature on the vertical profile of the heat mitigation implemented city reveal a stable PBL over the urban domain and a reduction of the vertical mixing associated with a pollution dome. This would lead to crossover phenomena above the PBL due to the decrease in vertical wind speed. Therefore, assessing the coupled regional impact of urban heat mitigation over the lower atmosphere at city-scale is urgent for sustainable urban planning.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Lidar-Based Evaluation of HRRR Performance in California’s Diablo Range

The performance of the NOAA High-Resolution Rapid Refresh (HRRR) model for capturing low-level winds near a wind energy production site during summer 2019 is evaluated. This study catalogs the ability of HRRR to predict boundary layer dynamics relevant to wind energy interests over complex terrain, which has presented challenges for weather and energy forecasting. Performance is evaluated by comparing HRRR output to wind-profiling Doppler lidars at Lawrence Livermore National Laboratory Site 300. HRRR captured the diurnal profile of horizontal winds in the observed 150-m layer, despite strong underpredictions (∼4 m s −1 ) during evening and nighttime hours. These underpredictions may be a result of local speedup flows observed by the lidars, which were unresolved in HRRR due to their small spatial extent. HRRR bias magnitude relative to observations was found to be minimal during days with synoptic-scale troughs and strong 850-hPa geopotential gradients, while bias magnitude was maximal during days with synoptic ridging and weak 850-hPa geopotential gradients. To translate wind speed predictions to energy forecasting, generic turbine models were used to estimate power generation for turbines characteristic of the nearby Altamont Pass Wind Resource Area. Results show that HRRR-based energy estimates predicted daytime power generation adequately relative to lidar-based estimates with an 18-h lead time (bias magnitude < 0.4 MW from 0900 to 1400 LT) but overpredicted power during the rest of the diurnal cycle (bias > 1 MW). These results demonstrate conditions under which HRRR performs well for wind energy applications in complex terrain, while highlighting biases that require further investigation to support usage of a high-resolution model for wind energy forecasts.

Boundary layer↗

DOE OSTI2021

Quincke Oscillations of Colloids at Planar Electrodes

Dielectric particles in weakly conducting fluids rotate spontaneously when subject to strong electric fields. Such Quincke rotation near a plane electrode leads to particle translation that enables physical models of active matter. In this Letter, we show that Quincke rollers can also exhibit oscillatory dynamics, whereby particles move back and forth about a fixed location. We explain how oscillations arise for micron-scale particles commensurate with the thickness of a field-induced boundary layer in the nonpolar electrolyte. This work enables the design of colloidal oscillators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗