Search NASA⌕ Search

SEARCH · Search NASA

Results for “shear layer turbulence”

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.

158 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A reduced-model (nSOLT) simulation of neutral recycling effects on plasma turbulence in the divertor region of MAST-U

The 2D scrape-off-layer turbulence code (nSOLT), which includes 1D kinetic neutral–plasma interactions, is applied to study effects of neutral recycling on plasma turbulence for parameters illustrative of the MAST-U divertor region. Neutral recycling is modeled by injecting a fraction of the parallel plasma flux to the divertor back into the simulation domain as a source of Franck–Condon-distributed neutrals. Stationary sources, concentrated at the magnetic separatrix, model plasma streaming into the divertor region from the upstream scrape-off-layer and sustain plasma turbulence absent neutral recycling. Starting from one such no-neutrals equilibrium, we initiate recycling in a numerical experiment designed to diagnose and identify the effects of various neutral–plasma interactions on the divertor plasma, divertor turbulence, and plasma exhaust. The onset of recycling triggers an initial burst of enhanced cross field plasma transport that is quelled by ionization cooling and charge–exchange (CX) friction, with growing neutral pressure, leading to a quiescent, turbulence-free state. Diagnosis of this transient burst reveals that (1) the sudden increase in plasma density due to ionization dominates the onset of the burst, (2) electron cooling due to ionization increases collisionality and disconnects blob filaments from the sheath, and (3) CX friction drives tripole polarization of a blob that can dominate the curvature-driven dipole polarization, leading to the stagnation of blob propagation and reduced radial turbulent transport. It is shown that CX friction is negligible compared to sheath physics in determining equilibrium mean flow shearing rates, for parameters considered herein (specifically a short connection length to the divertor target), while it can significantly reduce interchange-instability growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Large-Eddy Simulation of Flow Over Boeing Gaussian Bump Using Multiagent Reinforcement Learning Wall Model: Preprint

We develop a wall model for large-eddy simulation (LES) that takes into account various pressure-gradient effects using multi-agent reinforcement learning. The model is trained using low-Reynolds-number flow over periodic hills with agents distributed on the wall at various computational grid points. It utilizes a wall eddy-viscosity formulation as the boundary condition to apply the modeled wall shear stress. Each agent receives states based on local instantaneous flow quantities at an off-wall location, computes a reward based on the estimated wall-shear stress, and provides an action to update the wall eddy viscosity at each time step. The trained wall model is validated in wall-modeled LES of flow over periodic hills at higher Reynolds numbers, and the results show the effectiveness of the model on flow with pressure gradients. The analysis of the trained model indicates that the model is capable of distinguishing between the various pressure gradient regimes present in the flow. To further assess the robustness of the developed wall model, simulations of flow over the Boeing Gaussian bump are conducted at a Reynolds number of 2 x 10^6, based on the free-stream velocity and the bump width. The results of mean skin friction and pressure on the bump surface, as well as the velocity statistics of the flow field, are compared to those obtained from equilibrium wall model (EQWM) simulations and published experimental data sets. The developed wall model is found to successfully capture the acceleration and deceleration of the turbulent boundary layer on the bump surface, providing better predictions of skin friction near the bump peak and exhibiting comparable performance to the EQWM with respect to the wall pressure and velocity field. We also conclude that the subgrid-scale model is crucial to the accurate prediction of the flow field, in particular the prediction of separation.

boundary layer↗
DOE OSTI2025

Characterization of HRRR-simulated rotor layer wind speeds and clouds along the coast of California

Stratocumulus clouds, with their low cloud base and top, affect the atmospheric boundary layer wind and turbulence profile, thereby modulating wind energy resources. GOES satellite data reveal an abundance of stratocumulus clouds in the late spring and summer months off the coast of northern and central California, where there are active plans to deploy floating offshore wind farms at two lease areas (near Morro Bay and Humboldt). Since the fall of 2020, two buoys equipped with multiple instruments, including Doppler lidar, have been deployed for about 1 year in these wind farm lease areas to assess the rotor layer wind conditions in these locations. The objective of this study is to evaluate how well the High-Resolution Rapid Refresh (HRRR) model represents stratocumulus cloud characteristics and turbine-relevant rotor layer winds (surface to 300 m) by comparing HRRR simulations with buoy and satellite observations. We first find that the HRRR model reproduces the seasonal cycle of cloud top height reasonably well in these regions. However, during the warm season – especially at Morro Bay – the HRRR-simulated stratocumulus clouds tend to have lower tops by about 150 m and exhibit weaker diurnal cycles than satellite observations. Our analysis also shows that rotor layer wind speeds and vertical shear are stronger at Humboldt than at Morro Bay, and both are generally stronger under clear-sky conditions. Finally, the HRRR model bias in rotor layer wind speed is small under cloudy conditions but larger and dependent on observed wind speed under clear skies. Specifically, HRRR underestimates wind speeds at Morro Bay and overestimates them at Humboldt under clear-sky conditions.

17 WIND ENERGY↗
DOE OSTI2022

Similarity of the turbulent kinetic energy dissipation rate distribution in the upper mixed layer of the tropical Indian Ocean

Turbulence within the upper ocean mixed layer plays a key role in various physical, biological, and chemical processes. Between September and November 2011, a dataset of 570 vertical profiles of the turbulent kinetic energy (TKE) dissipation rate, as well as conventional hydrological and meteorological data, were collected in the upper layer of the tropical Indian Ocean. These data were used to statistically analyze the vertical distribution of the TKE dissipation rate in the mixed layer. The arithmetic-mean method made the statistical TKE dissipation rate profile more scattered than the median and geometric-mean methods. The statistical TKE dissipation rate were respectively scaled by the surface buoyancy flux and the TKE dissipation rate at the mixed-layer base. It was found that the TKE dissipation rate scaled by that at the mixed-layer base exhibited better similarity characteristics than that scaled by the surface buoyancy flux, whether the stability parameter D/|L MO | was greater or less than 10, indicating that the TKE dissipation rate at the mixed-layer base is a better characteristic scaling parameter for reflecting the intrinsic structure of the TKE dissipation rate in the mixed layer, where D and L MO are respectively the mixed-layer thickness and the Monin-Obukhov length scale. The parameterization of the TKE dissipation rate at the mixed-layer base on the shear-driven dissipation rate and the surface buoyancy flux was further explored. It was found that the TKE dissipation rate at the mixed-layer base could be well fitted by a linear combination of three terms: the wind-shear-driven dissipation rate, the surface buoyancy flux, and a simple nonlinear coupling term of these two.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Experimental and numerical study of the shielding effect of two tandem rough cylinders in flow-induce oscillation

The shielding effect of the downstream cylinder in flow induced oscillation (FIO) of two cylinders arranged in tandem is studied experimentally and numerically at Reynolds number 30,000 to 120,000. Both cylinders are in one degree-of-freedom, transverse-oscillations, and have turbulence stimulation in the form of selective surface roughness to expand FIO beyond vortex-induced vibration (VIV) into galloping. Shielding of the downstream cylinder has a negative effect on harnessing hydrokinetic energy. Further, to study its effect and mechanics, selective cases are studied both numerically and experimentally and discussed to demonstrate the shielding effect on the downstream cylinder and understand its cause. The main conclusions are: (1) The shielding effect for the downstream cylinder shows a strong relation to the damping ratio. As the damping ratio increases, the shielding effect is mitigated. Additionally, the oscillation of the rear cylinder becomes stable and shows stable frequency. (2) In the VIV region, as the stiffness and natural frequency increase, the shielding effect decreases substantially. (3) In the VIV region, the vorticity of the vortices shedding from both the upper and the lower sides of the downstream cylinder does not accumulate enough due to the attraction by the vortices shed from the upstream cylinder, thus resulting in partial suppression of the oscillation on the downstream cylinder. (4) In the galloping region, the shielding effect for the downstream cylinder depends on whether the vorticity near the downstream cylinder is strengthened by the vortices generated by the shear layers of the upstream cylinder or weakened.

42 ENGINEERING↗
DOE OSTI2020

Instability of an electron-plasma shear layer in an externally imposed strain flow

The E x B shear instability of a two-dimensional (2D) filament (i.e., a thin, rectangular strip perpendicular to the magnetic field) of magnetized pure electron plasma is investigated experimentally in the presence of an externally imposed strain flow. Data are acquired using a specialized Penning–Malmberg trap in which strain flows can be applied in 2D by biasing segmented electrodes surrounding the plasma. The E x B drift dynamics are well-described by the Drift-Poisson equations, which are isomorphic to the 2D Euler equations describing ideal fluids. Thus, the experimental results correspond to the Rayleigh instability of a shear layer in a 2D ideal fluid, where the electron density is analogous to the fluid vorticity. Shear layers are prepared by stretching initially axisymmetric electron vortices using a strong, applied strain flow. The data at early times are in quantitative agreement with a linear model which extends Rayleigh’s work to account for the influence of an external strain flow. In the presence of weak strain, the system approximately maintains a phase relationship that corresponds to an instantaneous Rayleigh eigenmode. The instability develops into the nonlinear regime later in time and at smaller spatial scales as the strain rate is increased. A secondary vortex pairing instability is observed, but it is suppressed when the strain-to-vorticity ratio exceeds roughly 0.025. In this way, vorticity transport perpendicular to the filament is diminished due to the applied strain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Convective and Turbulent Motions in Nonprecipitating Cu. Part II: LES Simulated Cloud Represented by a Starting Plume

The dynamic structure of a small trade wind cumulus (Cu) is analyzed using a novel approach. Cu developing in a shear-free environment is simulated by 10-m-resolution LES model with spectral bin microphysics. The aim is to clarify the dynamical nature of cloud updraft zone (CUZ) including entrainment and mixing in growing Cu. The validity of concept stating that a cloud at developing state can be represented by a parcel or a jet is tested. To investigate dynamical entrainment in CUZ performed by motions with scales larger than the turbulence scales, the modeled fields of air velocity were filtered by wavelet filter that separated convective motions from turbulent ones. Two types of objects in developing cloud were investigated: small volume ascending at maximal velocity (point parcel) and CUZ. It was found that the point parcel representing the upper part of cloud core is adiabatic. The motion of the air in this parcel ascending from cloud base determines cloud-top height. The top-hat (i.e., averaged) values of updraft velocity and adiabatic fraction in CUZ are substantially lower than those in the point parcel. Evaluation of the terms in the dynamical equation typically used in 1D cloud parcel models show that this equation can be applied for calculation of vertical velocities at the developing stage of small Cu, at least up to the heights of the inversion layer. Dynamically, the CUZ of developing cloud resembles the starting plume with the tail of nonstationary jet. Both the top-hat vertical velocity and buoyancy acceleration linearly increase with the height, at least up to the inversion layer. An important finding is that lateral entrainment of convective (nonturbulent) nature has a little effect on the top-hat CUZ velocity and cannot explain the vertical changes of conservative variables qt and θl. In contrast, entrained air lifting inside CUZ substantially decreases top-hat liquid water content and its adiabatic fraction. Possible reasons of these effects are discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Effects of Surface and Top Wind Shear on the Spatial Organization of Marine Stratocumulus-Topped Boundary Layers

Typically, the convective nature of Stratocumulus topped boundary layers (STBL) involves the motion of updrafts and downdrafts, driven by surface fluxes and radiative cooling, respectively. The balance between shear and buoyant forcings at the surface can determine the organization of updrafts between cellular and roll structures. We investigate the effect of varying shear at the surface and top of the STBL using Large Eddy Simulations, taking DYCOMS II RF01 as a base case. We focus on spatial identification of the following features: coherent updrafts and downdrafts, and observe how they are affected by varying shear. Stronger surface shear organizes the updrafts in rolls, causes less well-mixed thermodynamic profiles, and decreases cloud fraction and liquid water path (LWP). Stronger top shear also decreases cloud fraction and LWP more than surface shear, by thinning the cloud from the top. Features with stronger top than surface shear are associated with a net downward momentum transport and show early signs of decoupling. Classifying updrafts and downdrafts based on their vertical span and horizontal size confirms the dominance of tall objects spanning the whole STBL. Tall objects occupy 30% of the volume in the STBL, while short ones occupy less than 1%. For updraft and downdraft fluxes, these tall objects explain 65% of the vertical velocity variance and 83% of the buoyancy flux, on average. Stronger top shear also weakens the contribution of downdrafts to the turbulent fluxes and tilts the otherwise vertical development of updrafts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

High speed PLIF study of the Richtmyer–Meshkov instability upon re-shock

The Richtmyer–Meshkov instability (RMI) of a twice-shocked gas interface is researched using high-speed planar-laser induced fluorescence (PLIF) in the Wisconsin Shock Tube Laboratory’s vertical shock tube. The initial condition (IC) is a shear layer with broadband diffuse perturbations at the interface between a helium-acetone mixture and argon. This IC is accelerated by a shock of nominal strength M = 1.8, and then accelerated again by the transmitted shock that reflects off the end wall of the tube. An estimate of the light gas mole fraction is extracted from high-speed imaging using an iterative process that accounts for the nonlinear temperature dependence of the acetone’s fluorescence quantum yield (FQY) and absorption cross-section. A vorticity deposition model for the initial growth rate after reshock is compared with the Mikaelian model for re-shock. Previously used in literature, the number of generations is shown to naturally arise from a normalisation of the scalar transport equation. A self-similar analysis is then further performed using the mole fraction data to explore the evolution of the RMI after reshock and the higher order moments of the light gas mole fraction are compared with a proposed model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Comparison of edge turbulence characteristics between DIII-D and C-Mod simulations with XGC1

The physical processes taking place at the separatrix and scrape-off layer regions are crucial for the operation of tokamaks as they govern the interaction of hot plasma with the vessel walls. Numerical modeling of the edge with state-of-the-art codes attempts to elucidate the complex interactions between neoclassical drifts, turbulence, poloidal and parallel flows that control the physical set-up of the SOL region. Here, we present post-processing analysis of simulation results from the gyrokinetic code XGC1, comparing and contrasting edge turbulence characteristics from a simulation of the DIII-D tokamak against a simulation of the Alcator C-Mod tokamak. We find that the equilibrium ExB flux across the separatrix has a similar poloidal pattern in both discharges which can be explained by rB-drifts and trapped ion excursions. Furthermore, collisionality is noted to play a major role in the way that it prevents local charge accumulations from having more global effects in the C-Mod case. In both cases, turbulent electron heat flux is observed to be higher than the ion one and is possibly related to the need of electrons to maintain quasineutrality through the only channel available to them for exiting the confinement. By Fourier analysis, we identify turbulent frequencies and growth rates of the dominant mode in both simulations. In the case of C-Mod, these numbers point to the presence of a drift wave. In the DIII-D case, further linear simulations with the Gene code reveal a trapped electron mode. Furthermore, using a blob detection and tracking tool, we present the amplitude and size distributions of the blobs from both simulations. The amplitude distributions are in qualitative agreement with experimental observations while the size distributions are consistent with the fact that most of the blobs are not connecting to the divertor plates and suggest that they are generated by the shearing of the turbulent modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Structural uncertainty assessment for fire-engulfed objects in crosswind: Establishing credibility for a multiphysics wall-modeled large-eddy simulation paradigm

A structural uncertainty validation study for a large-scale, fire-engulfed, elevated object subjected to crosswind is presented to establish the credibility of a high-fidelity, low-Mach, turbulent reacting flow wall-modeled large-eddy simulation (WMLES) approach that includes multiphysics coupling to participating media radiation and conjugate heat transfer. To establish that WMLES can accurately predict surface quantities including drag and pressure coefficient in the low-Mach crosswind regime, a foundational elevated isothermal cylinder validation case is presented at a similar gap-to-diameter ratio of 0.25, spanning the subcritical to supercritical drag regime (Re 𝐷 = 1.1 × 10 5 and 4.3 × 10 5 , respectively). Here, this study exercised both static and dynamic coefficient LES (Smagorinsky and 𝑘 sgs ) with both local and exchange-based velocity sampling. Results showcase that the drag crisis (or the sudden drop in drag coefficient at increased Re 𝐷 ) is well captured when using an exchange-based dynamic coefficient WMLES methodology, while noting lack of mesh convergence and overall drag and pressure coefficient predictively when using a static coefficient, local velocity sampling WMLES. For the 𝒪⁡(10) m JP-8 liquid pool fire crosswind validation study presented, two experimental crosswind configurations (2 m/s and 9.5 m/s) are showcased for a fire-engulfed mock fuselage roughly 4 m in diameter. Using the best model-form practices identified in the isothermal study, dynamic coefficient 𝑘 sgs exchange-based WMLES fire validation findings demonstrate accurate peak irradiation and skin temperature predictions as a function of crosswind magnitude. Excessive yaw in the low-crosswind fuselage configuration, consistent with experimental findings, captured a significant predicted asymmetry in flame attachment and heat flux toward the downwind cylindrical cap—indicative of axial vortex structures transporting the flame along the upper and lower fuselage leeward surface. All fire mesh resolution simulations captured the experimental finding that as crosswind increased, predicted flame shape and peak irradiation magnitude onto the fuselage transitioned from a windward to a leeward cylinder location due to the migration of the upper- to lower-shear fuel/air mixing layer thereby demonstrating the novelty, significance, and credibility of this high-fidelity WMLES reacting flow framework.

Domino, Stefan Paul [Sandia National Laboratories ↗
DOE OSTI2023

CROCUS Sodar Measurements of Lower Atmospheric Wind Profiles at Argonne Testbed for Multiscale Observational Science (ATMOS) Site

The Scintec MFAS Sodar (Multiple-Frequency Acoustic Sounder) is an autonomous, ground-based acoustic remote sensing system designed to measure vertical profiles of horizontal wind speed, wind direction, and vertical velocity in the lower atmosphere. The instrument transmits sequences of acoustic pulses and detects the Doppler-shifted sound waves backscattered by small-scale temperature and velocity fluctuations caused by atmospheric turbulence. From these Doppler shifts, the system derives three-dimensional wind vectors by combining radial velocities from multiple beam orientations.The MFAS Sodar operates with a first usable range gate beginning at approximately 30 m above ground level and a configurable vertical resolution of 10 m. Under favorable acoustic conditions, the system provides wind profiles extending up to 600 m above ground level. Measurements are processed into 15-minute averaged profiles containing wind speed, direction, vertical velocity, and diagnostic quantities such as signal-to-noise ratio and echo strength.This dataset was collected at the Argonne Testbed for Multiscale Observational Science (ATMOS) facility in Lemont, Illinois, as part of DOE's CROCUS Urban Integrated Field Laboratory (UIFL) initiative. The purpose of these observations is to characterize the vertical wind structure and boundary-layer evolution across the urban–suburban gradient of the greater Chicago region. In particular, these data are intended to improve understanding of how local meteorology, such as lake-breeze penetration, nocturnal jets, and daytime mixing, varies between the densely built urban core and the suburban periphery. The MFAS observations provide critical context for evaluating high-resolution model simulations and for integrating with complementary lidar, radar, and in-situ meteorological measurements within the broader CROCUS UIFL network.All data are archived in NetCDF (Network Common Data Form) format and include wind and diagnostic parameters. The files can be accessed and analyzed using standard software that supports NetCDF, such as Python (e.g., xarray, netCDF4), MATLAB, R (e.g., ncdf4, raster), or Panoply (NASA’s NetCDF visualization application).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Overview of the recent experimental research on the J-TEXT tokamak

The J-TEXT capability is enhanced compared to two years ago with several upgrades of its diagnostics and the increase of electron cyclotron resonance heating (ECRH) power to 1 MW. With the application of electron cyclotron wave (ECW), the ECW assisted plasma startup is achieved; the tearing mode is suppressed; the toroidal injection of 300 kW ECW drives around 24 kA current; fast electrons are generated with toroidal injected ECW and the runaway current conversion efficiency increases with ECRH power. The mode coupling between 2/1 and 3/1 modes are extensively studied. The coupled 2/1 and 3/1 modes usually lead to major disruption. Their coupling can be either suppressed or avoided by external resonant magnetic perturbation fields and hence avoids the major disruption. It is also found that the 2/1 threshold of external field is significantly reduced by a pre-excited 3/1 mode, which can be either a locked island or an external kink mode. The disruption control is studied by developing prediction methods capable of cross tokamak application and by new mitigation methods, such as the biased electrode or electromagnetic pellet injector. The high-density operation and related disruptions are studied from various aspects. Approaching the density limit, the collapse of the edge shear layer is observed and such collapse can be prevented by applying edge biasing, leading to an increased density limit. The density limit is also observed to increase, if the plasma is operated in the poloidal divertor configuration or the plasma purity is increased by increasing the pre-filled gas pressure or ECRH power during the start-up phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Influence of simple terrain on the spatial variability of a low-level jet and wind farm performance in the AWAKEN field campaign

In wind energy research, scientific challenges are often associated with complex terrain sites, where orography, vegetation, and buildings disrupt flow uniformity. However, even sites characterized as simple terrain can exhibit significant spatial variability in wind speed, particularly during stable boundary layers (SBLs) and low-level jets (LLJs). This study investigates these terrain interactions using both simulations and observations from the American WAKe ExperimeNt (AWAKEN). We employ a multiscale Weather Research and Forecasting (WRF) model simulation, integrating mesoscale forcing in the coarse domains and representing three rows of turbines from the King Plains wind farm as generalized actuator disks (GAD) in the large-eddy simulation (LES) domains. During a nocturnal LLJ event on 3 April 2023, the downstream, wake-affected turbine rows outperformed the upstream, unwaked row by 25 %–51 %. This counterintuitive result arises from terrain-induced streamwise variations in hub-height wind speed of approximately 4 m s −1 over 5 km – equivalent to ∼50 % of the upstream reference speed. This enhancement outweighs the wake-induced reduction in mean wind speed (∼12 %) and global blockage effects reported in the literature (∼1 %–3.4 %). The multiscale simulations capture the intra-farm spatial variability in power performance observed in SCADA data. Terrain-induced vertical displacement of the LLJ, coupled with large wind shear below the jet maximum, drives the substantial streamwise acceleration within the wind farm. These findings underscore the importance of accounting for spatial variability related to terrain, even in simple landscapes, particularly during LLJ conditions. Incorporating such effects into reduced-order modeling frameworks for wind farm design and control could significantly enhance their effectiveness.

17 WIND ENERGY↗