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.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

The Effect of Pressure on NO x Entitlement and Reaction Timescales in a Premixed Axial Jet-In-Crossflow

This paper investigates the pressure dependency of a lean premixed jet injected into a lean vitiated crossflow with an experimentally verified detailed chemistry computational fluid dynamics (CFD) model and 53 species considered. Experimental data were taken in an axially staged combustor with an optically accessible test section, allowing the use of particle image velocimetry (PIV) and CH* chemiluminescence techniques as well as point measurement of species concentration, temperature, and pressure. The experimental data cases at one, three, and five atmospheres were selected to describe the flame stabilization dependency on pressure and gain the required knowledge for an extrapolation to engine condition. Simulated exit nitrogen oxide levels were validated with experimental emission data, and a global emission trend for the NO reduction at elevated pressure and constant turbine inlet temperature level was defined. The nitrogen oxide benefit at elevated operating pressure was justified with the significantly smaller flame surface area: the analysis of the simulated spanwise and top-view profiles showed a relatively short receded core flame with nitrogen oxide production in the center at high pressure relative to a longer and larger shear layer flame at atmospheric condition that produced NO toward the inner and outer side of the flame. Furthermore, decomposition of the Damköhler number revealed the strong influence of the reaction timescales with higher reaction rates at elevated pressure, along with a moderate influence of the turbulent timescales, showing higher turbulence intensity in the lee-side recirculation zone at lower pressure.

03 NATURAL GAS↗
DOE OSTI2026

A buoyancy–shear–drag–scalar-based turbulence model for power-law acceleration-driven Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A previously developed phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing based on a general buoyancy–shear–drag model [O. Schilling, “A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing,” Physica D 402, 132238 (2020)] is extended to include active or passive scalar mixing and power-law acceleration-driven Rayleigh–Taylor mixing. The buoyancy–shear–drag equations are coupled to a scalar variance equation that is used to define the molecular mixing parameter θ m , and when the scalar is active, modifies the Rayleigh–Taylor and Kelvin–Helmholtz mixing layer growth parameters to depend on the asymptotic value of this parameter, θ mol . Here, the scalar variance equation is closed by algebraically or differentially modeling the scalar variance dissipation rate. Nonlinear analytical solutions of the model are obtained in the total and separate bubble and spike mixing layer width formulations with the algebraic scalar variance dissipation rate for each instability, which are then used to calibrate the mechanical and scalar equation coefficients to predict specific values of physical observables and molecular mixing parameters. Surrogate mechanical and scalar turbulent fields can be constructed by multiplying a presumed self-similar spatial profile by appropriate functions of the width and its time derivative, and of the scalar obtained by solving the ordinary differential model equations. The explicit modeling and solution of turbulent transport equations are not required. The bubble and spike mixing layer width and scalar variance equations are then solved numerically for constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing, confirming that the prescribed level of molecular mixing is correctly predicted and illustrating the spatiotemporal evolution of the scalar fields.

Buoyancy–drag↗
DOE OSTI2022

Ice-shelf ocean boundary layer dynamics from large-eddy simulations

Abstract. Small-scale turbulent flow below ice shelves is regionally isolated and difficult to measure and simulate. Yet these small-scale processes, which regulate heat and salt transfer between the ocean and ice shelves, can affect sea-level rise by altering the ability of Antarctic ice shelves to “buttress” ice flux to the ocean. In this study, we improve our understanding of turbulence below ice shelves by means of large-eddy simulations at sub-meter resolution, capturing boundary layer mixing at scales intermediate between laboratory experiments or direct numerical simulations and regional or global ocean circulation models. Our simulations feature the development of an ice-shelf ocean boundary layer through dynamic ice melting in a regime with low thermal driving, low ice-shelf basal slope, and strong shear driven by the geostrophic flow. We present a preliminary assessment of existing ice-shelf basal melt parameterizations adopted in single component or coupled ice-sheet and ocean models on the basis of a small parameter study. While the parameterized linear relationship between ice-shelf melt rate and far-field ocean temperature appears to be robust, we point out a little-considered relationship between ice-shelf basal slope and melting worthy of further study.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Revisiting the validity of eddy viscosity models for predicting airflow over water waves

In this study, we revisit the validity of eddy viscosity models for predicting wave-induced airflow disturbances over ocean surface waves. We first derive a turbulence curvilinear model for the phase-averaged Navier–Stokes equations, extending the work of Cao, Deng & Shen (2020 J. Fluid Mech. 901, A27), by incorporating turbulence stress terms previously neglected in the linearised viscous curvilinear model. To verify our formulation, we perform a priori tests by numerically solving the model using mean wind and turbulence stress profiles from large-eddy simulations (LES) of airflow over waves across various wave ages. Results show that including turbulence stress terms improves wave-induced airflow predictions compared with the previous viscous curvilinear model. We further show that using a standard mixing-length eddy viscosity yields inaccurate predictions at certain wave ages, as it fails to capture wave-induced turbulence, which fundamentally differs from mean shear-driven turbulence. The LES data show that accurate representations of wave-induced stresses require a complex-valued eddy viscosity. The maximum magnitude of this eddy viscosity scales as ∼𝑢 𝜏 ⁢𝜁 𝑖𝑛𝑛𝑒𝑟 , where 𝑢 𝜏 is the friction velocity and 𝜁 𝑖𝑛𝑛𝑒𝑟 is the inner-layer thickness, the height at which the eddy-turnover time matches the wave advection time scale. This scaling aligns with the prediction by Belcher & Hunt (1993 J. Fluid Mech. 251, 109–148). Overall, the findings demonstrate that traditional eddy viscosity models are inadequate for capturing wave-induced turbulence. More sophisticated turbulence models are essential for the accurate prediction of airflow disturbances and form drag in wind–wave interaction models.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Turbulent dynamo action and its effects on the mixing at the convective boundary of an idealized oxygen-burning shell

Convection is one of the most important mixing processes in stellar interiors. Hydrodynamic mass entrainment can bring fresh fuel from neighboring stable layers into a convection zone, modifying the structure and evolution of the star. Because flows in stellar convection zones are highly turbulent, multidimensional hydrodynamic simulations are fundamental to accurately capture the physics of mixing processes. Under some conditions, strong magnetic fields can be sustained by the action of a turbulent dynamo, adding another layer of complexity and possibly altering the dynamics in the convection zone and at its boundaries. In this study, we used our fully compressible S EVEN -L EAGUE H YDRO code to run detailed and highly resolved three-dimensional magnetohydrodynamic simulations of turbulent convection, dynamo amplification, and convective boundary mixing in a simplified setup whose stratification is similar to that of an oxygen-burning shell in a star with an initial mass of 25 M ⊙ . We find that the random stretching of magnetic field lines by fluid motions in the inertial range of the turbulent spectrum (i.e., a small-scale dynamo) naturally amplifies the seed field by several orders of magnitude in a few convective turnover timescales. During the subsequent saturated regime, the magnetic-to-kinetic energy ratio inside the convective shell reaches values as high as 0.33, and the average magnetic field strength is ∼10 10 G. Such strong fields efficiently suppress shear instabilities, which feed the turbulent cascade of kinetic energy, on a wide range of spatial scales. The resulting convective flows are characterized by thread-like structures that extend over a large fraction of the convective shell. The reduced flow speeds and the presence of magnetic fields with strengths up to 60% of the equipartition value at the upper convective boundary diminish the rate of mass entrainment from the stable layer by ≈20% as compared to the purely hydrodynamic case.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

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↗
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

Measurement-driven large-eddy simulations of a diurnal cycle during a wake-steering field campaign

Abstract. High-fidelity flow modeling with data assimilation enables accurate representation of the wind farm operating environment under realistic, nonstationary atmospheric conditions. Capturing the temporal evolution of the turbulent atmospheric boundary layer is critical to understanding the behavior of wind turbines under operating conditions with simultaneously varying inflow and control inputs. This paper has three parts: the identification of a case study during a field evaluation of wake steering; the development of a tailored mesoscale-to-microscale coupling strategy that resolved local flow conditions within a large-eddy simulation (LES), using observations that did not completely capture the wind and temperature fields throughout the simulation domain; and the application of this coupling strategy to validate high-fidelity aeroelastic predictions of turbine performance and wake interactions with and without wake steering. The case study spans 4.5 h after midnight local time, during which wake steering was toggled on and off five times, achieving yaw offset angles ranging from 0 to 17°. To resolve nonstationary nighttime conditions that exhibited shear instabilities, the turbulence field was evolved starting from the diurnal cycle of the previous day. These background conditions were then used to drive wind farm simulations with two different models: an LES with actuator disk turbines and a steady-state engineering wake model. Subsequent analysis identified two representative periods during which the up- and downstream turbines were most nearly aligned with the mean wind direction and had observed yaw offsets of 0 and 15°. Both periods corresponded to partial waking on the downstream turbine, which had errors in the LES-predicted power of 4 % and 6 %, with and without wake steering. The LES was also able to capture conditions during which an upstream turbine wake induced a speedup at a downstream turbine and increased power production by up to 13 %.

17 WIND ENERGY↗
DOE OSTI2023

Turbulent Kinetic Energy Budgets over Gentle Topography Covered by Forests

Large-eddy simulations of flow over a “horizontally” uniform model forest are used to investigate the effects of gentle topography on the turbulent kinetic energy (TKE) budget within the canopy roughness sublayer. Despite significant differences between simulations using idealized sinusoidal topography and real topography of the Amazon forest, results indicate that the effects of topography are located predominantly in the upper canopy and above, and are mostly caused by mean advection of TKE. The “horizontally” averaged TKE budget from idealized and real gentle topographies are almost identical to that for flat terrain, including a clear inertial layer above the roughness sublayer in which shear production is balanced by local dissipation. At topography crests, where observational towers are usually located, mean vertical advection of TKE can be as important as horizontal advection. Here, we propose the use on an approximate TKE balance equation to estimate mean advection from single tower measurements, and introduce a new advection index that can be used as a proxy to quantify the importance of the topography on the TKE budget.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Long‐Term Observations of Turbulence Vertical Velocity Spectra in a Convective Mixed Layer: Dependence on Land‐Surface Forcing in the U.S. Southern Great Plains

Doppler lidar vertical velocity retrievals were analyzed for the scale and structure of mixed-layer turbulence over a 7-year period, on fair-weather warm season days (442 cases) in the U.S. Southern Great Plains. Based on the hypothesized influence of land-surface forcing and updraft size on convective boundary layer clouds, a spectral analysis was performed to quantify effects of surface forcing (surface buoyancy flux, friction velocity, and evaporative fraction) on turbulence scale. Significant (order-of-magnitude) variations in spectral density were found in the energy-production subrange and mesoscale regimes. Integral scale decreased with increasing buoyancy flux (and Monin-Obukhov stability parameter), while spectral density in the energy-production subrange increased, implying a transition to buoyancy-driven cellular structures with narrower updrafts. The influence of stability parameter was limited to the neutral to convective transition, and could not explain the wide variation in spectral density in the mesoscale regime. However, high friction velocity (u* > 0.5 m s –1 ) was associated with larger integral scale (updraft width), and a greater portion of spectral density in the mesoscale regime. Lidar profiles and radar reflectivity for individual cases revealed evidence of roll and wave-like structures contributing to mesoscale variability on high friction velocity days, suggesting shear instabilities on days with wetter land surface conditions and smaller buoyancy flux. The role of friction velocity emphasizes the multivariate nature of surface influences on turbulence, and implies that variance-based turbulence closures may not be adequate for capturing effects of surface forcing on updraft size.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Composited structure of non-precipitating shallow cumulus clouds

The normalized distributions of thermodynamic and dynamical variables both within and outside shallow clouds are investigated through a composite algorithm using large-eddy simulations of oceanic and continental cases. The normalized magnitude is maximum near the cloud centre and decreases outwards. While relative humidity (RH) and cloud liquid water ( q l ) decrease smoothly to match the environment, the vertical velocity, virtual potential temperature (θ v ), and potential temperature (θ) perturbations have more complicated behaviour towards the cloud boundary. Below the inversion layer, θ' v becomes negative before the vertical velocity has turned from an updraft to a subsiding shell outside the cloud, indicating the presence of a transition zone where the updraft is negatively buoyant. Due to the downdraft outside the cloud and enhanced horizontal turbulent mixing across the edge, the normalized turbulent kinetic energy (TKE) and horizontal turbulent kinetic energy (HTKE) decrease more slowly from the cloud centre outwards than the thermodynamic variables. The distributions all present asymmetric structures in response to the vertical wind shear, with more negatively buoyant air, stronger downdrafts, and larger TKE on the downshear side. We discuss several implications of the distributions for theoretical models and parameterizations. Positive buoyancy near the cloud base is mostly due to the virtual effect of water vapour, emphasizing the role of moisture in triggering. The mean vertical velocity is found to be approximately half the maximum vertical velocity within each cloud, providing a constraint to achieve possible power-law distributions for some models. Finally, the normalized distributions for different variables are used to estimate the vertical heat and moisture fluxes within clouds. The results suggest that distributions near the cloud edge and variability of maximum perturbations need careful treatment. The fluxes are underestimated in the inversion layer because cloud-top downdrafts cannot be captured well.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Simulations of a hypersonic turbulent boundary layer over wavy surfaces

Here, we conduct large-eddy simulations of a Mach 5.84 cold wall turbulent boundary layer over one-dimensional wavy walls with varying amplitudes and wavelengths. Across all wall topologies, a series of alternating shock and expansion waves is shown to influence the entire boundary layer, and generate repeating wave patterns in the turbulent stresses, dispersive stresses, and turbulent kinetic energy budget. The series of alternating shocks and expansions imposes repeating adverse and favourable pressure gradients across the wavy wall, and at sufficient wall amplitude, triggers flow separation in the trough of the wave. Flow separation is demonstrated to influence the behaviour of wall pressure fluctuations over the wavy wall. In attached flows, the prominent frequencies are consistent with integral-scale boundary layer turbulence, whereas in separated flows, a two-decade frequency range is present, akin to two-dimensional shock–boundary layer interactions. Counter-rotating streamwise-oriented structures are observed on the windward side of the wave, which diminish over the wave crest. A conditional analysis demonstrates that these structures are present in the upstream boundary layer, and are amplified with increasing wall amplitude. An examination of the Görtler number and probability density function (PDF) of the fluctuating lateral wall shear stress demonstrates the strong correlation between a large Görtler number and growth of the PDF tail density, suggesting that the amplification of the counter-rotating streamwise-oriented structures are linked to centrifugal instabilities in regions of streamline concavity.

boundary layers↗
DOE OSTI2026

Suppression of ion temperature gradient modes by Alfvén activity above a drive threshold in DIII-D

A recent study demonstrates that the suppression of ion temperature gradient (ITG) modes can occur during the nonlinear evolution of toroidicity-induced Alfvén eigenmodes (TAEs) in the DIII-D tokamak [Du et al., Phys. Rev. Lett. 135, 265101 (2025)]. In that work, ITG was suppressed by the formation of a narrow, TAE-induced shear flow layer, whose shearing rate exceeded the ITG decorrelation rate. Here, the shear flow arises from an imbalance between Reynolds and Maxwell stress forces, as the TAE departs from the conventional shear Alfvén wave polarization. This follow-up paper systematically identifies the plasma conditions required for robust ITG suppression through a series of comparative experiments. The results show that TAEs routinely suppress ITG turbulence in plasmas with higher local safety factor (q), elevated fast-ion beta, and larger populations of fast ions on passing orbits, i.e., conditions consistent with the substantial TAE drive. Database analysis further reveals the existence of a threshold in TAE drive for ITG suppression. That is, once the drive exceeds a critical value within a favorable q window, the system undergoes the nonlinear bifurcation process, characterized by a sharp increase in TAE saturated amplitudes at a nearly fixed fast ion drive, suppression of ITG turbulence, and the formation of an internal thermal transport barrier at the localized radii.

Du, X. D. [General Atomics, San Diego, CA (United ↗
DOE OSTI2022

Scientific challenges to characterizing the wind resource in the marine atmospheric boundary layer

Abstract. With the increasing level of offshore wind energy investment, it is correspondingly important to be able to accurately characterize the wind resource in terms of energy potential as well as operating conditions affecting wind plant performance, maintenance, and lifespan. Accurate resource assessment at a particular site supports investment decisions. Following construction, accurate wind forecasts are needed to support efficient power markets and integration of wind power with the electrical grid. To optimize the design of wind turbines, it is necessary to accurately describe the environmental characteristics, such as precipitation and waves, that erode turbine surfaces and generate structural loads as a complicated response to the combined impact of shear, atmospheric turbulence, and wave stresses. Despite recent considerable progress both in improvements to numerical weather prediction models and in coupling these models to turbulent flows within wind plants, major challenges remain, especially in the offshore environment. Accurately simulating the interactions among winds, waves, wakes, and their structural interactions with offshore wind turbines requires accounting for spatial (and associated temporal) scales from O(1 m) to O(100 km). Computing capabilities for the foreseeable future will not be able to resolve all of these scales simultaneously, necessitating continuing improvement in subgrid-scale parameterizations within highly nonlinear models. In addition, observations to constrain and validate these models, especially in the rotor-swept area of turbines over the ocean, remains largely absent. Thus, gaining sufficient understanding of the physics of atmospheric flow within and around wind plants remains one of the grand challenges of wind energy, particularly in the offshore environment. This paper provides a review of prominent scientific challenges to characterizing the offshore wind resource using as examples phenomena that occur in the rapidly developing wind energy areas off the United States. Such phenomena include horizontal temperature gradients that lead to strong vertical stratification; consequent features such as low-level jets and internal boundary layers; highly nonstationary conditions, which occur with both extratropical storms (e.g., nor'easters) and tropical storms; air–sea interaction, including deformation of conventional wind profiles by the wave boundary layer; and precipitation with its contributions to leading-edge erosion of wind turbine blades. The paper also describes the current state of modeling and observations in the marine atmospheric boundary layer and provides specific recommendations for filling key current knowledge gaps.

17 WIND ENERGY↗
DOE OSTI2024

Direct numerical simulation of open-channel flow over a heterogeneous particle bed at low relative submergence

In this article, this study investigates turbulent open-channel flows over beds of irregularly arranged particles, using direct numerical simulations at a friction Reynolds number of Re T = 300⁠. Two distinct cases are examined: a polydisperse bed (P800) composed of multiple layers of randomly distributed spheres of varying sizes, and a monodisperse bed (M1015) formed by a random distribution of uniform sized spheres, with a bottommost single layer of varied-sized particles to introduce realistic randomness. Our investigation unveils a rich network of low- and high-speed streaks within the flow field, exhibiting distinctive behaviors in different bed configurations. The P800 case presents a poorly organized flow pattern induced by the varied particle sizes and arrangements, while the M1015 case shows a more regular flow pattern, marked by larger streaks. We also observe that total wall shear stress is substantially influenced by surface roughness-induced drag, extending beyond the effects documented in existing studies of open-channel flows. The present study reveals intricate secondary flow patterns over irregular particle beds. Large-scale circulations are discerned around particle crests in the P800 case and localized circulations with increased turbulence in the M1015 case. Furthermore, analysis of Reynolds stress tensor components indicates that roughness disrupts coherent turbulent eddies, consequently mitigating peak stress. We quantify correlations between drag force and local fluid velocity fluctuations. Notably, a larger deviation in drag is observed in the P800 case compared to M1015, accentuating the influence of particle size and distribution on fluid–particle interactions.

42 ENGINEERING↗
DOE OSTI2025

Freestream turbulence effects on unsteady wind turbine loads and wakes: An IDDES study

We investigate numerically the effects of freestream turbulence on the unsteady aerodynamics and wakes of the National Renewable Energy Laboratory Phase VI wind turbine rotor for increasing wind speed. Turbulence is modeled using the Improved Delayed Detached-Eddy Simulation (IDDES) method. As a first step, a detailed mesh resolution study is conducted with the decaying freestream turbulence model at turbulence intensity of 0.5%. Our blade-resolved IDDES simulations show that grid-independent average torque and thrust results can be achieved with relatively coarse meshes, whereas dramatically higher mesh resolution is required for grid-independent results for power spectral densities of thrust force, especially in the deep-stall regime. Comparing the loads with the Shear-Stress Transport model demonstrates the superiority of IDDES in predicting massively separated flows. The aerodynamic performance and wake predictions with the decaying freestream turbulence model are compared with the synthetic freestream turbulence model. Both models predict nearly the same loads, spectral energy content, and wake characteristics. The properties of both the near- and far-wake regions are then examined. Furthermore, we show that separated boundary layers accelerate turbulent mixing and entrainment of the external flow, which results in faster wake recovery. The effect of increasing turbulence intensity to 6% is investigated using the synthetic freestream turbulence model. In contrast with the fully attached boundary layer, higher freestream turbulence in deep stall does not significantly affect the loads and vortex-shedding characteristics. However, the turbulent mixing in the wake is enhanced, which further hastens the recovery of the self-similar velocity profile. In general, increasing the wind speed at high turbulence intensity shifts the recovery farther upstream and increases the wake width.

17 WIND ENERGY↗
DOE OSTI2023

How the birth and death of shear layers determine confinement evolution: from the L → H transition to the density limit

Electric field profile structure—especially its shear—is a natural order parameter for the edge plasma, and characterizes confinement regimes ranging from the H-mode (Wagner et al. 1982 Phys. Rev. Lett.49, 1408–1412 (doi:10.1103/PhysRevLett.49.1408)) to the density limit (DL) (Greenwald et al. 1988 Nucl. Fusion28, 2199–2207 (doi:10.1088/0029-5515/28/12/009)). The theoretical developments and lessons learned during 40 years of H-mode studies (Connor & Wilson 1999 Plasma Phys. Control. Fusion42, R1–R74 (doi:10.1088/0741-3335/42/1/201); Wagner 2007 Plasma Phys. Control. Fusion49, B1–B33 (doi:10.1088/0741-3335/49/12b/s01)) are applied to the shear layer collapse paradigm (Hong et al. 2017 Nucl. Fusion58, 016041 (doi:10.1088/1741-4326/aa9626)) for the onset of DL phenomena. Results from recent experiments on edge shear layers and DL phenomenology are summarized and discussed in the light of L → H transition physics. The theory of shear layer collapse is then developed. In this work, we demonstrate that shear layer physics captures both the well known current (Greenwald) scaling of the DL (Greenwald 2002 Plasma Phys. Control. Fusion44, R27–R53 (doi:10.1088/0741-3335/44/8/201); Greenwald et al. 2014 Phys. Plasmas21, 110501 (doi:10.1063/1.4901920)), as well as the emerging power scaling (Zanca, Sattin, JET Contributors 2019 Nucl. Fusion59, 126011 (doi:10.1088/1741-4326/ab3b31)). The derivation of the power scaling theory exploits an existing model, originally developed for the L → H transition (Diamond, Liang, Carreras, Terry 1994 Phys. Rev. Lett.72, 2565–2568 (doi:10.1103/PhysRevLett.72.2565); Kim & Diamond 2003 Phys. Rev. Lett.90, 185006 (doi:10.1103/PhysRevLett.90.185006)). We describe the enhanced particle transport events that occur following shear layer collapse. Open problems and future directions are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Pore-resolved investigation of turbulent open channel flow over a randomly packed permeable sediment bed

Pore-resolved direct numerical simulations are performed to investigate the interactions between streamflow turbulence and groundwater flow through a randomly packed porous sediment bed for three permeability Reynolds numbers, Re K = 2.56 , 5.17 and 8.94, representative of natural stream or river systems. Time–space averaging is used to quantify the Reynolds stress, form-induced stress, mean flow and shear penetration depths, and mixing length at the sediment–water interface (SWI). Here, the mean flow and shear penetration depths increase with Re K and are found to be nonlinear functions of non-dimensional permeability. The peaks and significant values of the Reynolds stresses, form-induced stresses, and pressure variations are shown to occur in the top layer of the bed, which is also confirmed by conducting simulations of just the top layer as roughness elements over an impermeable wall. The probability distribution functions (p.d.f.s) of normalized local bed stress are found to collapse for all Reynolds numbers, and their root-mean-square fluctuations are assumed to follow logarithmic correlations. The fluctuations in local bed stress and resultant drag and lift forces on sediment grains are mainly a result of the top layer; their p.d.f.s are symmetric with heavy tails, and can be well represented by a non-Gaussian model fit. The bed stress statistics and the pressure data at the SWI potentially can be used in providing better boundary conditions in modelling of incipient motion and reach-scale transport in the hyporheic zone.

turbulence simulation↗