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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Generation of mean flows in rotating anisotropic turbulence: The case of solar near-surface shear layer

Results from helioseismology indicate that the radial gradient of the rotation rate in the near-surface shear layer (NSSL) of the Sun is independent of latitude and radius. Theoretical models using the mean-field approach have been successful in explaining this property of the NSSL, while global direct or large-eddy magnetoconvection models have so far been unable to reproduce this. We investigate the reason for this discrepancy by measuring the mean flows, Reynolds stress, and turbulent transport coefficients under conditions mimicking those in the solar NSSL. Simulations with as few ingredients as possible to generate mean flows were studied. These ingredients are inhomogeneity due to boundaries, anisotropic turbulence, and rotation. The parameters of the simulations were chosen such that they matched the weakly rotationally constrained NSSL. The simulations probe locally Cartesian patches of the star at a given depth and latitude. The depth of the patch was varied by changing the rotation rate such that the resulting Coriolis numbers covered the same range as in the NSSL. We measured the turbulent transport coefficient relevant for the nondiffusive (Λ-effect) and diffusive (turbulent viscosity) parts of the Reynolds stress and compared them with predictions of current mean-field theories. A negative radial gradient of the mean flow is generated only at the equator where meridional flows are absent. At other latitudes, the meridional flow is comparable to the mean flow corresponding to differential rotation. We also find that the meridional components of the Reynolds stress cannot be ignored. Additionally, we find that the turbulent viscosity is quenched by rotation by about 50% from the surface to the bottom of the NSSL. Our local simulations do not validate the explanation for the generation of the NSSL from mean-field theory where meridional flows and stresses are neglected. However, the rotational dependence of the turbulent viscosity in our simulations agrees well with theoretical predictions. Moreover, our results agree qualitatively with global convection simulations in that an NSSL can only be obtained near the equator.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

On the enhancement of boundary layer skin friction by turbulence: an angular momentum approach

Turbulence enhances the wall shear stress in boundary layers, significantly increasing the drag on streamlined bodies. Other flow features such as free stream pressure gradients and streamwise boundary layer growth also strongly influence the local skin friction. In this paper, an angular momentum integral (AMI) equation is introduced to quantify these effects by representing them as torques that alter the shape of the mean velocity profile. This approach uniquely isolates the skin friction of a Blasius boundary layer in a single term that depends only on the Reynolds number most relevant to the flow's engineering context, so that other torques are interpreted as augmentations relative to the laminar case having the same Reynolds number. The AMI equation for external flows shares this key property with the so-called FIK relation for internal flows (Fukagata et al. , Phys. Fluids , vol. 14, 2002, pp. L73–L76). Without a geometrically imposed boundary layer thickness, the length scale in the Reynolds number for the AMI equation may be chosen freely. After a brief demonstration using Falkner–Skan boundary layers, the AMI equation is applied as a diagnostic tool on four transitional and turbulent boundary layer direct numerical simulation datasets. Regions of negative wall-normal velocity are shown to play a key role in limiting the peak skin friction during the late stages of transition, and the relative strengths of terms in the AMI equation become independent of the transition mechanism a very short distance into the fully turbulent regime. The AMI equation establishes an intuitive, extensible framework for interpreting the impact of turbulence and flow control strategies on boundary layer skin friction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Reduction of transport due to magnetic shear in gyrokinetic simulations of the scrape-off layer

We report the effect of varying magnetic shear on scrape-off layer turbulence and profiles is studied via electromagnetic gyrokinetic simulations of a helical scrape-off layer model. We develop a model helical geometry with magnetic shear and a corresponding field-aligned coordinate system, which is used for simulations with the Gkeyll code. We find that perpendicular transport is reduced in cases with stronger shear, resulting in higher peak particle and heat fluxes to the endplates. Electromagnetic effects slightly increase transport in strong shear cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2021

The diurnal cycle of entrainment and detrainment in LES of the Southern Ocean driven by observed surface fluxes and waves

Empirical rules for both entrainment and detrainment are developed from LES of the Southern Ocean boundary layer when the turbulence, stratification, and shear cannot be assumed to be in equilibrium with diurnal variability in surface flux and wave (Stokes drift) forcing. A major consequence is the failure of downgradient eddy viscosity, which becomes more serious with Stokes drift and is overcome by relating the angle between the stress and shear vectors to the orientations of Lagrangian shear to the surface and of local Eulerian shear over 5 m. Thus, the momentum flux can be parameterized as a stress magnitude and this empirical direction. In addition, the response of a deep boundary layer to sufficiently strong diurnal heating includes boundary layer collapse and the subsequent growth of a morning boundary layer, whose depth is empirically related to the time history of the forcing, as are both morning detrainment and afternoon entrainment into weak diurnal stratification. Below the boundary layer, detrainment rules give the maximum buoyancy flux and its depth, as well a specific stress direction. Another rule relates both afternoon and nighttime entrainment depth and buoyancy flux to surface layer turbulent kinetic energy production integrals. These empirical relationships are combined with rules for boundary layer transport to formulate two parameterizations; one based on eddy diffusivity and viscosity profiles and another on flux profiles of buoyancy and of stress magnitude. Evaluations against LES fluxes show the flux profiles to be more representative of the diurnal cycle, especially with Stokes drift.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Impact of atmospheric turbulence on performance and loads of wind turbines: knowledge gaps and research challenges

Wind energy harvesting from the atmosphere takes place in the atmospheric boundary layer. The boundary layer shear and buoyancy create three-dimensional turbulent eddies spanning a range of scales that form a continuous forward cascade of kinetic energy to the smallest scales of motion where energy is dissipated. Large-scale atmospheric circulations modulate the boundary layer turbulence, characterized by coherence and intermittency. As wind turbines grow in size and the integrated control of both turbines and wind farms spans greater distances, the relationship between the scales of atmospheric turbulence and the design and operation of wind energy facilities has entered new territory. The boundary layer turbulence impacts both wind turbine power production and turbine loads. Optimizing wind turbine and wind farm performance requires an understanding of how turbulence affects both wind turbine efficiency and reliability. While the characteristics of atmospheric boundary layer turbulence have been observed and studied in detail over the last few decades, there are still significant gaps in our understanding of the impact of turbulence on wind power resources and wind farm operations. This paper outlines the current state of turbulence research relevant to wind energy applications and points to gaps in our knowledge that need to be addressed to effectively utilize wind resources.

Kosović, Branko [Johns Hopkins Univ., Baltimore, M↗
DOE OSTI2021

Spatially and Temporally Resolved Velocimetry for Hypersonic Flows

The development of new hypersonic flight vehicles is limited by the physical understanding that may be obtained from ground test facilities. This has motivated the present development of a temporally and spatially resolved velocimetry measurement for Sandia National Laboratories (SNL) Hypersonic Wind Tunnel (HWT) using Femtosecond Laser Electronic Excitation Tagging (FLEET). First, a multi-line FLEET technique has been created for the first time and tested in a supersonic jet, allowing simultaneous measurements of velocities along multiple profiles in a flow. Secondly, two different approaches have been demonstrated for generating dotted FLEET lines. One employs a slit mask pattern focused into points to yield a dotted line, allowing for two- or three-component velocity measurements free of contamination between components. The other dotted-line approach is based upon an optical wedge array and yields a grid of points rather than a dotted line. Two successful FLEET measurement campaigns have been conducted in SNL’s HWT. The first effort established optimal diagnostic configurations in the hypersonic environment based on earlier benchtop reproductions, including validation of the use of a 267 nm beam to boost the measurement signal-to-noise ratio (SNR) with minimal risk of perturbing the flow and greater simplicity than a comparable resonant technique at 202 nm. The same FLEET system subsequently was reconstituted to demonstrate the ability to make velocimetry measurements of hypersonic turbulence in a realistic flow field. Mean velocity profiles and turbulence intensity profiles of the shear layer in the wake of a hypersonic cone model were measured at several different downstream stations, proving the viability of FLEET as a hypersonic diagnostic.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Flow regime and Reynolds number variation effects on the mixing behavior of parallel flows

The hydraulic single-phase mixing of three parallel rectangular channels is experimentally investigated at various Reynolds numbers (Re) and flow regime combinations. Particle Image Velocimetry results for seven mixing cases are presented and discussed with varying Re combinations ranging from 1,824 to 20,844. While all cases result in the same Re ratio of ~0.69 between the inner and outer flows, two cases represent multi-regime mixing with the inner-outer regime pair of laminar-transitional and transitional-turbulent, while the other 5 cases are all characteristic of turbulent mixing with varying levels of turbulence. The outer channels initially share characteristics with a backward facing step. The center channel is found to initially behave like a slot jet, but then sees a significant increase in velocity decay. This inner flow velocity decay increased dramatically in the laminar-transitional mixing case, whose centerline velocity decay was ~6 times larger than the decay in the turbulent mixing cases. Second order statistics revealed a consistent mixing layer thickness of ~0.1 hydraulic diameters for all the cases but showed more intense shearing in the multi-regime mixing cases. The combined point and thereby the mixing layer length is determined using centerline velocity decay profiles, which show a much more aggressive mixing in multi-regime flows. Multi-regime mixing demonstrated superior characteristics relative to turbulent mixing due to a more dramatic velocity decay in the inner flow and a shorter mixing length. The contributions of this work include communicating the benefits of multi-regime mixing and providing detailed characterization efforts that can serve future efforts for validating computational models. Here this research also lays the groundwork for future studies aimed at achieving high levels of mixing without a severe penalty in pressure drop.

42 ENGINEERING↗
DOE OSTI2025

Predictive turbulence-driven flux model of scrape-off layer widths across confinement regimes in tokamaks

Reliable scrape-off layer (SOL) profile decay lengths predictions are needed to design and operate future tokamaks. The present manuscript describes a new model based on turbulent transport that is able to predict SOL widths for both L-mode and H-mode plasmas. The model is based upon the sheared-spectral filament paradigm (Peret et al (WEST Team) 2022 Phys. Plasmas 29 072306), however, incorporating the effects of thermal transport in order to calculate the parallel heat fluxes. The effects of magnetic shear and ExB shear on the cross-field transport are crucial to explain the shorter SOL decay lengths found in H-mode. The model is validated against a database of thousands of DIII-D L-mode and H-mode SOL profiles. We also calculate SOL decay length predictions in terms of plasma and engineer control parameters, which are in agreement with the multi-machine empirical H-mode scaling (Eich et al (ASDEX Upgrade Team and JET EFDA Contributors) 2013 Nucl. Fusion 53 093031), however, with an additional device geometry dependence. ITER SOL width predictions by the model are 3 times higher than the empirical scaling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Adaptive mesh based combustion simulations of direct fuel injection effects in a supersonic cavity flame-holder

Here, we present high-fidelity reacting simulations of a supersonic cavity flame-holder configuration. The focus of this work is on flame stabilization brought about by varying the location of fuel injection in a cavity stabilized supersonic flow of air. Central to our approach is a compressible multi-species reacting flow solver that uses adaptive-mesh-refinement (AMR), enabling the resolution of flame, shock-waves, boundary-layers, and small-scale structures in the computational domain. Our analysis indicates that fuel injection closer to the ramp at the aft end of the cavity allows for greater mixing and lower peak temperatures compared to fuel injection upstream that is closer to the backward facing step of the cavity. This difference is mainly due to greater turbulent fluctuations generated from the shear-layer towards the cavity ramp, thereby enhancing the mixing of fuel and air. A low frequency oscillatory behaviour in heat-release and pressure was also observed for the upstream injection case while a much higher-frequency phenomena was observed in the near-ramp injection case. By identifying the important physical determinants of the combustion processes, this study illustrates a promising pathway to design and optimize direct fuel injection strategies in supersonic cavity flame-holders that can improve flame stability, combustion efficiency, and reduce emissions.

42 ENGINEERING↗
DOE OSTI2021

Influence of gas expansion on the propagation of a premixed flame in a spatially periodic shear flow

It has been previously demonstrated that thermal gas expansion might have a role in boundary layer flashback of premixed turbulent flames, inducing local flow-reversal in the boundary layer's low-velocity streaks on the reactants’ side of the flame and facilitating its upstream propagation. Here, we perform a two-dimensional numerical investigation of the interaction between a periodic shear flow and a laminar premixed flame. The periodic shear is a simplified model for the oncoming prolonged streamwise velocity streaks with alternating regions of high and low velocities found in turbulent boundary layers in the vicinity of the walls. The parametric study focuses on the amplitude and wavelength of the periodic shear flow and on the gas expansion ratio (unburnt-to-burnt density ratio). With the increase of the amplitudes of the periodic shear flow and of the gas expansion, the curved flame velocity increases monotonically. The flame velocity dependence on the periodic shear wavelength is non-monotonic, which is consistent with previous theoretical studies of curved premixed flame velocity. The flame shape that is initially formed by the oncoming periodic shear appears to be metastable. At a later stage of the flame propagation, the flame shape transforms into the stationary one dominated by the Darrieus-Landau instability.

42 ENGINEERING↗
DOE OSTI2021

Bounds on edge shear layer persistence while approaching the density limit

This paper details the theory of edge shear layer collapse as the density approaches the Greenwald density limit. It significantly extends earlier work, which was restricted in applicability. The zonal shear flow screening length is calculated for banana, plateau and Pfirsch–Schluter regimes. Poloidal field scaling persists in the plateau regime. Neoclassical screening and drift wave–zonal flow dynamics are combined in a theory, which is then reduced to a predator–prey model. Zonal noise, due to incoherent mode coupling, is retained. The threshold condition for edge shear layer collapse is computed, and linked to a critical value of the dimensionless parameter ${\rho }_{\mathrm{s}}/\sqrt{{\rho }_{\mathrm{s}\mathrm{c}}{L}_{n}}$. Here ${\rho }_{s}$ is the ion sound radius, ${\rho }_{sc}$ is the zonal flow screening length and ${L}_{n}$ is the equilibrium density scale length. The limiting initial edge density for shear layer collapse is derived and shown to scale favorably with the plasma current. Here, the results are discussed in light of the density limit and Ohmic phenomenology.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Near-cancellation of up- and down-gradient momentum transport in forced magnetized shear-flow turbulence

Visco-resistive magnetohydrodynamic turbulence, driven by a two-dimensional unstable shear layer that is maintained by an imposed body force, is examined by decomposing it into dissipationless linear eigenmodes of the initial profiles. The down-gradient momentum flux, as expected, originates from the large-scale instability. However, continual up-gradient momentum transport by large-scale linearly stable but nonlinearly excited eigenmodes is identified and found to nearly cancel the down-gradient transport by unstable modes. The stable modes effectuate this by depleting the large-scale turbulent fluctuations via energy transfer to the mean flow. This establishes a physical mechanism underlying the long-known observation that coherent vortices formed from nonlinear saturation of the instability reduce turbulent transport and fluctuations, as such vortices are composed of both the stable and unstable modes, which are nearly equal in their amplitudes. The impact of magnetic fields on the nonlinearly excited stable modes is then quantified. Furthermore, even when imposing a strong magnetic field that almost completely suppresses the instability, the up-gradient transport by the stable modes is at least two-thirds of the down-gradient transport by the unstable modes, whereas for weaker fields, this fraction reaches up to 98%. These effects are persistent with variations in magnetic Prandtl number and forcing strength. Finally, continuum modes are shown to be energetically less important, but essential for capturing the magnetic fluctuations and Maxwell stress. A simple analytical scaling law is derived for their saturated turbulent amplitudes. It predicts the falloff rate as the inverse of the Fourier wavenumber, a property which is confirmed in numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Development and Verification of an Improved Wake-Added Turbulence Model in FAST.Farm

We introduce a generalized wake-added turbulence (WAT) model in the multiphysics, multiturbine simulation tool FAST.Farm. The WAT model introduces additional small-scale turbulence that represents the breakdown of vortical structures and shear layers in the wake. The article describes the development, implementation, calibration, and verification of the model. The novelties of the model include support for wake asymmetry, buildup of WAT across the wind farm, and secondary effects of wake-induced turbulence (e.g., wake meandering) driven by smaller-scale turbulence structures that arise from wake breakdown. Large-eddy simulations were run to support the calibration of the WAT parameters and verification of the model. Previous studies hypothesized that the lack of WAT modeling was the source of underprediction of fatigue loads, in particular for cases at low turbulence intensities and/or stable atmospheric boundary layers. This study confirms that the newly implemented WAT model enhances the loads predictions in these cases.

17 WIND ENERGY↗
DOE OSTI2026

Validation of OpenPronghorn for Periodic Hill Flow Separation

OpenPronghorn is an open-source, MOOSE-based thermal-hydraulics simulation tool used for advanced reactor analysis. As an open-source code, it offers a transparent framework for validating governing equations, assumptions, and numerical methods against established benchmarks. This study evaluates OpenPronghorn's RANS turbulence model against the ERCOFTAC Case 81 periodic hill benchmark, a standard test case for separated turbulent flow featuring curved-wall separation, recirculation, shear-layer development, and reattachment. Streamwise velocity profiles predicted by OpenPronghorn were compared to reference LES data at multiple x/h locations. Results show that OpenPronghorn captures the overall trend of the velocity profiles, but the largest discrepancies occur in the separated-flow region, where turbulence is highly anisotropic and strongly affected by adverse pressure gradients and wall curvature—conditions that are inherently difficult for standard RANS models to resolve. Future work will test alternative k-e model variants and correction terms to improve prediction accuracy in this region.

42 - ENGINEERING↗
DOE OSTI2025

On Master-Length Scale Formulations for Stable Conditions in Turbulence Closure Models

Three formulations of the turbulence-length scales used in numerical modelling of atmospheric flows are compared. The comparison is made using the Mellor–Yamada–Nakanishi–Niino turbulence closure model within the stable boundary layer local similarity framework. With an appropriate choice of model constants, the model predictions are barely discernible and compare well with the empirical data obtained from the SHEBA campaign.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Environmental sensitivities of shallow-cumulus dilution – Part 2: Vertical wind profile

This second part of a numerical study on shallow-cumulus dilution focuses on the sensitivity of cloud dilution to changes in the vertical wind profile. Insights are obtained through large-eddy simulations of maritime and continental cloud fields. In these simulations, the speed of the initially uniform geostrophic wind and the strength of geostrophic vertical wind shear in the cloud and subcloud layer are varied. Increases in the cloud-layer vertical wind shear (up to 9 ms -1 km -1 ) lead to 40%–50% larger cloud-core dilution rates compared to their respective unsheared counterparts. When the background wind speed, on the other hand, is enhanced by up to 10 ms -1 and subcloud-layer vertical wind shear develops or is initially prescribed, the dilution rate decreases by up to 25%. The sensitivities of the dilution rate are linked to the updraft strength and the properties of the entrained air. Increases in the wind speed or vertical wind shear result in lower vertical velocities across all sets of experiments with stronger reductions in the cloud-layer wind shear simulation (27%–47%). Weaker updrafts are exposed to mixing with the drier surrounding air for a longer time period, allowing more entrainment to occur (i.e., the “core-exposure effect”). However, reduced vertical velocities, in concert with increased cloud-layer turbulence, also assist in widening the humid shell surrounding the cloud cores, leading to entrainment of more humid air (i.e., the “core–shell dilution effect”). In the experiments with cloud-layer vertical wind shear, the core-exposure effect dominates and the cloud-core dilution increases with increasing shear. Conversely, when the wind speed is increased and subcloud-layer vertical wind shear develops or is imposed, the core–shell dilution effect dominates to induce a buffering effect. The sensitivities are generally stronger in the maritime simulations, where weaker sensible heat fluxes lead to narrower, more tilted, and, therefore, more suppressed cumuli when cloud-layer shear is imposed. Moreover, in the experiments with subcloud wind shear, the weaker baseline turbulence in the maritime case allows for a larger turbulence enhancement, resulting in a widening of the transition zones between the cores and their environment, leading to the entrainment of more humid air.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

DNS of Flame Stabilization Dynamics of a Swirl-Stabilized Spray Burner Using Sustainable Aviation Fuels

Concern for emission reduction has motivated the development of new cost-effective alternative sustainable aviation fuels (SAFs). Drop-in SAFs with blends of the certified and the alternative fuels are beneficial because they do not require engine modifications for use in current aviation engines. Combustion characteristics are of particular concern when comparing a new fuel with Jet-A. Lean blow-off and flame stabilization have been identified to be governed by the fuel propensity to autoignition, i.e. derived cetane number. However, fundamental investigation is required to identify the impact of derived cetane number on the combustion mode of flame stabilization in a realistic combustor. The combustion mode has also a direct impact on turbulent combustion modeling closures. DNS simulations are performed in the low-Mach solver of the Pele Suite called PeleLMeX. Lagrangian multi-phase modeling is used to capture the liquid spray injection of Jet-A (reference fuel) and C1 as a representative of a low cetane number SAF. Adaptive Mesh Refinement (AMR) is used to enable a more efficient simulation of a more realistic domain size and embedded boundary treatment is used to model a bluff-body geometry. Local extinction and edge flame propagation were observed for both fuels. The edge flame propagation mode was quantified in terms of a Damkohler number defined as the ratio between progress variable reaction rate and its diffusive flux. Initial analysis suggests that a mixed mode combustion occurs for the edge flame propagation, with flame propagation assisted by ignition and autoignition co-existing for both Jet-A and C1 flames. The analysis of the extinction region shows a larger progress variable for Jet-A which can further increase the local displacement speed and present a faster reignition of the stoichiometric mixture fraction.

adaptive mesh refinement↗