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

Applying Time-Parallelization to Turbulent Flows

Parallelization of the temporal domain is explored for the solution of turbulent flows. Multigrid reduction-in-time (MGRIT) is used to advance the large-scale fluid dynamics in time sequentially on the coarsest space-time grid but propagate the information in time parallel on all other levels. The goal of this process is to accurately and efficiently resolve the coarse-scale turbulence structure and use that to drive the fine-scales of the turbulent flow. The extra forcing from nonlinear multigrid facilitates the coupling and interaction between fine and coarse scales, through which the multiscale nonlinear physics is properly captured. Adaptive mesh refinement is employed to finely resolve only the regions with strong gradients, which provides further computational efficiency. The underlying computational fluid dynamics solver is a fourth-order finite-volume scheme with the standard 4-stage Runge-Kutta method. An advanced approach is devised and implemented to enable MGRIT to solve highly turbulent flows successfully. Furthermore, the method is applied to solve a Taylor-Green vortex problem and a doubleshear-layer turbulent mixing flow. Results are promising, validating that MGRIT with the filtering approach has the potential to efficiently solve general turbulent flows.

Computational Fluid Dynamics↗
DOE OSTI2022

Numerical investigation of the influence of shear and thermal stratification on the wind turbine tip‐vortex stability

Summary The interaction between wind turbine wakes and atmospheric turbulence is characterised by complex dynamics. In this study, two major components of the atmospheric boundary layer dynamics have been isolated, namely, the mean velocity profile shear and the thermal stratification, to examine their impact on the near‐wake development by undertaking a series of highly resolved large‐eddy simulations. Subsequently, instantaneous flow fields are extracted from the simulations and used to conduct Fourier analysis and proper orthogonal decomposition (POD) and compute the mean kinetic energy fluxes by different POD modes to better understand the tip‐vortex instability mechanisms. Our findings indicate that shear can significantly affect the breakup of the wind turbine tip‐vortices and the shape and stable length of the wake, whereas thermal stratification seems to only have limited contribution to the spatial development of the near‐wake field. Finally, our analysis shows that the applied perturbation frequency determines the tip‐vortex breakup location as it controls the onset of the mutual inductance instability.

17 WIND ENERGY↗
DOE OSTI2020

Evaluating turbulence models at high energy densities

The RESHOCK campaign at Lawrence Livermore National Laboratory has been working on increasing our understanding of the evolution of turbulent, unstable plasma interfaces which are applicable in various HED applications including ICF implosions. A common approach to model such interfaces involves the use of Renolds-Averaged-Navier-Stokes (RANS) models whose parameters have been constrained by theory and experiments. These parameter sets are not unique and the applicability of one set of tuned model parameters to systems with differing shock strengths or material densities has not been well tested. Our new NIF experiments, along with related studies of turbulent shear flows, provide data and model validation for plasma interfaces in the HED regime. Our experiments measure mixing-layer width at an unstable interface in the plasma regime and are specifically designed to challenge mix models of the Reynolds-Averaged-Navier-Stokes (RANS) type where turbulence is assumed to be fully developed, i.e. the flow exhibits a broad spectrum of length scales without memory of the initial condition.4 We utilize precise control of the initial interface conditions (densities on both sides of the interface, an initial rippled perturbation pattern consisting of equally weighted wavelengths between 10 and 20 microns, and initial material compositions) along with a repeatable drive history

36 MATERIALS SCIENCE↗
DOE OSTI2024

Exploring ship track spreading rates with a physics-informed Langevin particle parameterization

Abstract. The rate at which aerosols spread from a point source injection, such as from a ship or other stationary pollution source, is critical for accurately representing subgrid plume spreading in a climate model. Such climate model results will guide future decisions regarding the feasibility and application of large-scale intentional marine cloud brightening (MCB). Prior modeling studies have shown that the rate at which ship plumes spread may be strongly dependent on meteorological conditions, such as precipitating versus non-precipitating boundary layers and shear. In this study, we apply a Lagrangian particle model (PM-ABL v1.0), governed by a Langevin stochastic differential equation, to create a simplified framework for predicting the rate of spreading from a ship-injected aerosol plume in sheared, precipitating, and non-precipitating boundary layers. The velocity and position of each stochastic particle is predicted with the acceleration of each particle being driven by the turbulent kinetic energy, dissipation rate, momentum variance, and mean wind. These inputs to the stochastic particle velocity equation are derived from high-fidelity large-eddy simulations (LES) equipped with a prognostic aerosol–cloud microphysics scheme (UW-SAM) to simulate an aerosol injection from a ship into a cloud-topped marine boundary layer. The resulting spreading rate from the reduced-order stochastic model is then compared to the spreading rate in the LES. The stochastic particle velocity representation is shown to reasonably reproduce spreading rates in sheared, precipitating, and non-precipitating cases using domain-averaged turbulent statistics from the LES.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A two-fluid single-column model of turbulent shallow convection. Part II: Single-column model formulation and numerics

The two-fluid single-column model of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) is extended to include moisture and horizontal wind shear. Turbulent kinetic energy is introduced as a prognostic variable, dependence on a diagnosed boundary-layer height is removed, and subfilter fluxes are approximated using a two-fluid version of a Mellor–Yamada scheme. Three mechanisms for entrainment and detrainment processes are introduced, which represent entrainment of unstable air at the surface, forced detrainment of air at the top of the boundary/cloud layers, and turbulent mixing that relaxes the convective fluid to a reference profile. A semi-implicit Eulerian discretization replaces the semi-implicit semi-Lagrangian implementation of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) to improve numerical stability and conservation. The equations for the implicit time step are solved using a quasi-Newton method, which is shown to perform well in numerical tests for conservation and convergence. The two-fluid single-column model presented in this article will be applied to simulations of shallow cumulus convection in Part III.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Mixing Driven by Critical Reflection of Near-Inertial Waves over the Texas–Louisiana Shelf

Studies of internal wave-driven mixing in the coastal ocean have been focused on internal tides, while wind-driven near-inertial waves (NIWs) have received less attention in this regard. Here, this study demonstrates a scenario of NIW-driven mixing over the Texas–Louisiana shelf. Supported by a high-resolution simulation over the shelf, the NIWs driven by land–sea breeze radiate downward at a sharp front and enhance the mixing in the bottom boundary layer where the NIWs are focused because of slantwise critical reflection. The criterion for slantwise critical reflection of NIWs is ω=$\sqrt{f{^2}+s^{2}_{bot}N{^2}(1-s_{ρ}/s_{bot})}$ (where ω is the wave frequency, s bot is the bottom slope, and s ρ is the isopycnal slope), under the assumption that the mean flow is in a thermal wind balance and only varies in the slope-normal direction. The mechanism driving the enhanced mixing is explored in an idealized simulation. During slantwise critical reflection, NIWs are amplified with enhanced shear and periodically destratify a bottom boundary layer via differential buoyancy advection, leading to periodically enhanced mixing. Turbulent transport of tracers is also enhanced during slantwise critical reflection of NIWs, which has implications for bottom hypoxia over the Texas–Louisiana shelf.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

100 kHz OH PLIF of High-Pressure H2/NH3 and H2/CH4 Flames

The reaction zones in a multi-element, multi-stage, micromix injector are investigated using planar laser-induced fluorescence (PLIF) imaging of the hydroxyl (OH) radical at a 100 kHz repetition rate. The combustor is operated at a pressure of ∼10 bar. Multiple fuel blends of hydrogen, methane, and ammonia are considered to study the effect of changing fuel composition on the flame structure. Compact regions of high OH signal are observed in each flame for fuel mixtures with high hydrogen content, indicating rapid reactant consumption in the shear layers. Instantaneous images exhibit wrinkled flame fronts and the formation of unreacted pockets due to highly turbulent burning. The addition of methane is observed to reduce the reactivity and flame speed, resulting in more diffuse regions of OH between each flame. The shape of the reactant core is altered, and the length scale of the flame surface wrinkling is increased. Ammonia addition has a more severe effect on the distribution of OH and the local flame structure than methane addition. A time series of images shows the presence of OH vortical structures and significant turbulent fluctuations in the flames near the wall, suggesting the presence of a recirculation zone that could drive the observed disparity in flame lengths.

Hodge, Alexander J. [Purdue University]↗
DOE OSTI2024

Influence of Cross Perturbations on Turbulent Kelvin–Helmholtz Instability

Kelvin–Helmholtz instability has been studied extensively in 2D. This study attempts to address the influence of turbulent flow and cross perturbation on the growth rate of the instability and the development of mixing layers in 3D by means of direct numerical simulation. Two perfect gases are considered to be working fluids moving as opposite streams, inducing shear instability at the interface between the fluids and resulting in Kelvin–Helmholtz instability. The results show that cross perturbation affects the instability by increasing the amplitude growth while adding turbulence has almost no effect on the amplitude growth. Furthermore, by increasing the turbulence intensity, a more distinct presence of the inner flow can be seen in the mixing layer of the two phases, and the presence of turbulence expands the range of high-frequency motion significantly due to turbulence structures. The results give a basis for which 3D Kelvin–Helmholtz phenomena should be further investigated using numerical simulation for predictive modeling, beyond the use of simplified 2D theoretical models.

Sementilli, Mae (ORCID:0000000207125928)↗
DOE OSTI2025

Challenge Problem 1: Preliminary Results of the Direct Numerical Simulation of Transient Flows

This report presents the first direct numerical simulations (DNS) of transient mixed convection in an idealized downcomer-like channel (Challenge Problem 1, Phase II). Using the GPU-accelerated NekRS solver, we modeled a sudden decay in driving pressure, mimicking loss-of-flow events, and tracked the resulting evolution of Reynolds number, boundary-layer structure, turbulence statistics, and heat-transfer metrics. Key findings include the systematic thickening and eventual asymmetry of velocity and thermal boundary layers under buoyant deceleration; minimal “memory” lag in Reynolds shear stress and TKE profiles when sampled at matching Re, yet clear shifts of peak locations toward the cooled wall; overshoots in transient eddy-viscosity and eddy-diffusivity (and corresponding sub-unity turbulent Prandtl numbers) on the cooled side; and a pronounced transient Nusselt-number enhancement driven by wall-temperature inertia and residual eddy mixing. These effects combined to offer a temporary cooling margin above steady-state predictions during reactor LOF transients. Future work will extend this work to a more complex “Case II” geometry (90° turn + lower plenum) and generate multi-Re/Pr datasets for data-driven turbulence closures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

An enrichment wall modeling framework for spectral element methods

In the present work, a first-of-its-kind enrichment wall-model is developed within the spectral element method (SEM) framework for large-eddy simulations (LES) of wall-bounded turbulent flows. The method augments the polynomial solution in the wall-adjacent elements with an analytical law-of-the-wall enrichment function representing the mean velocity near the wall. In the solution representation, this enrichment function captures the large gradients in the boundary layer, which allows the polynomial modes to represent the turbulent fluctuations. The enriched solution is able to resolve the shear stress at the wall without any modification to the no-slip wall boundary conditions, which allows for greater accuracy in the near-wall region compared to traditional methods. The enrichment wall modeling approach is implemented in a high-order SEM computational fluid dynamics solver, Nek5000, and its performance is assessed in turbulent channel flow wall-modeled LES for a range of Reynolds numbers. It is demonstrated that the enrichment wall-model improves solution accuracy on under-resolved near-wall grids as compared to traditional shear stress wall-models.

42 ENGINEERING↗
DOE OSTI2020

Reduction of blob-filament radial propagation by parallel variation of flows: Analysis of a gyrokinetic simulation

Data from the XGC1 gyrokinetic simulation is analyzed to understand the three-dimensional spatial structure and the radial propagation of blob-filaments generated by quasisteady turbulence in the tokamak edge pedestal and scrape-off layer plasma. Spontaneous toroidal flows vary in the poloidal direction and shear the filaments within a flux surface resulting in a structure that varies in the parallel direction. Here, this parallel structure allows the curvature and grad-B induced polarization charge density to be shorted out via parallel electron motion. As a result, it is found that the blob-filament radial velocity is significantly reduced from estimates which neglect parallel electron kinetics, broadly consistent with experimental observations. Conditions for when this charge shorting effect tends to dominate blob dynamics are derived and compared with the simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

How eigenmode self-interaction affects zonal flows and convergence of tokamak core turbulence with toroidal system size

Self-interaction is the process by which a microinstability eigenmode that is extended along the direction parallel to the magnetic field interacts non-linearly with itself. This effect is particularly significant in gyrokinetic simulations accounting for kinetic passing electron dynamics and is known to generate stationary $E\times B$ zonal flow shear layers at radial locations near low-order mode rational surfaces (Weikl et al. Phys. Plasmas, vol. 25, 2018, 072305). Here, we find that self-interaction, in fact, plays a very significant role in also generating fluctuating zonal flows, which is critical to regulating turbulent transport throughout the radial extent. Unlike the usual picture of zonal flow drive in which microinstability eigenmodes coherently amplify the flow via modulational instabilities, the self-interaction drive of zonal flows from these eigenmodes are uncorrelated with each other. It is shown that the associated shearing rate of the fluctuating zonal flows therefore reduces as more toroidal modes are resolved in the simulation. In simulations accounting for the full toroidal domain, such an increase in the density of toroidal modes corresponds to an increase in the toroidal system size, leading to a finite system size effect that is distinct from the well-known profile shearing effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The role of RF-induced E×B flows in the mitigation of scrape-off-layer convective transport during ion cyclotron resonance heating

Abstract While multiple experiments have reported a decrease in intermittent fluctuations in the far Scrape-Off-Layer (SOL) during ion cyclotron resonance heating (Antar et al 2010 Phys. Rev. Lett. 105 165001, Li et al 2022 Nucl. Eng. Technol. 54 207–19, Antar et al 2012 Nucl. Fusion 52 103005), the physical mechanism behind this observation has not been fully established yet. In this work, we demonstrate, for the first time, a direct correlation between the amplitude of RF-induced E×B flows and turbulence suppression in the far SOL. Using the Gas Puff Imaging (GPI) diagnostic on Alcator C-Mod, we show again that Ion Cyclotron Range of Frequencies can significantly alter the flow in the SOL and introduce a shear layer in regions magnetically connected close to the antenna (Cziegler et al 2012 Plasma Phys. Control. Fusion 54 105019). With the 4-strap field-aligned antenna operated in dipole phasing, the ratio of the power coupled by the central two straps to the power coupled by the outer two straps was varied. The resulting RF-induced radial electric field magnitude thus varied, and we show that the impact on the far SOL turbulence correlates with the modified E×B velocity. We then apply a newly-developed blob tracking algorithm (Han et al 2022 Sci. Rep. 12 18142) to higher-resolution GPI videos in order to directly observe the process of blob shearing by RF-induced E×B flows. We show that the radially sheared poloidal flows act as a transport barrier by stretching, stopping, and destroying filaments, which is consistent with the observed difference in turbulence statistics.

Diab, R. (ORCID:0000000342065868)↗
DOE OSTI2020

Three-dimensional signatures of self-similarity in a high-energy-density plasma shear-driven mixing layer

A hydrodynamic shear mixing layer experiment at the National Ignition Facility had previously demonstrated Eulerian scaling of integrated, late-time quantities, including turbulent kinetic energy. In this manuscript, the experiment is repeated with new materials. In this work, using the new dataset, we demonstrate that Euler-number scalings hold not just for late time, but dynamically throughout the experiment, for measurements in all three spatial dimensions. Incorporating the dynamic scaling leads to an enhanced calculation that the heavier of the two scaled experiments has approached three generations of mergers of its primary instability's structures and a consistent observation of such a merger in action in the lighter of the two scaled experiments. Furthermore, the improved scrutiny of the time evolution of instability structures leads to sharper estimates of turbulent kinetic energy, including a demonstration of different behaviors correlating with surface roughness (quantitatively consistent with transitions between laminar and turbulent initial states), as predicted by a Reynolds-averaged turbulent model, which evidently correctly handles the differing shock-roughness interactions to drive its internal state of the model into different regimes. Altogether, a picture arises of the analytical improvements in treating these variations (of times, densities, and roughnesses) as a unified whole and of multiple ways by which deviations from the scaling could indicate an onset of non-hydrodynamic behavior. Such deviations were not expected for these experiments (which models correctly indicated would remain hydrodynamic) but could be introduced by, for example, imposing external fields or increasing drive energy to test conditions relevant to inertial confinement fusion or other high-energy-density experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Stochastic modal velocity field in rough-wall turbulence

Stochastically generated instantaneous velocity profiles are used to reproduce the outer region of rough-wall turbulent boundary layers in a range of Reynolds numbers extending from the wind tunnel to field conditions. Each profile consists in a sequence of steps, defined by the modal velocities and representing uniform momentum zones (UMZs), separated by velocity jumps representing the internal shear layers. Height-dependent UMZ is described by a minimal set of attributes: thickness, mid-height elevation, and streamwise (modal) and vertical velocities. These are informed by experimental observations and reproducing the statistical behaviour of rough-wall turbulence and attached eddy scaling, consistent with the corresponding experimental datasets. Sets of independently generated profiles are reorganized in the streamwise direction to form a spatially consistent modal velocity field, starting from any randomly selected profile. The operation allows one to stretch or compress the velocity field in space, increases the size of the domain and adjusts the size of the largest emerging structures to the Reynolds number of the simulated flow. By imposing the autocorrelation function of the modal velocity field to be anchored on the experimental measurements, we obtain a physically based spatial resolution, which is employed in the computation of the velocity spectrum, and second-order structure functions. The results reproduce the Kolmogorov inertial range extending from the UMZ and their attached-eddy vertical organization to the very-large-scale motions (VLSMs) introduced with the reordering process. The dynamic role of VLSM is confirmed in the –u'w' co-spectra and in their vertical derivative, representing a scale-dependent pressure gradient contribution.

42 ENGINEERING↗
DOE OSTI2021

A multispecies turbulence model for the mixing and de-mixing of miscible fluids

A de-mix capable model for turbulence in compressible, variable density flows is proposed. The model is based on the Besnard-Harlow-Rauenzahn (BHR) family of models (Besnard D, Harlow F, Rauenzahn R, et al. Turbulence transport equations for variable-density turbulence and their relationship to two-field models. NM (United States): Los Alamos National Laboratory; 1992 (Technical Report LA-12303-MS), but is extended to track the evolution of the turbulent fluxes and fluctuations in the material mass fractions for each species present. The new evolution equations are introduced without requiring additional closures or new empirically tuned coefficients relative to previous BHR models, and are shown to improve the model’s ability to reproduce the behaviour of simulations containing mixing layers that are constrained by a stabilising force. As a result, the model is tested in a range of canonical flows including Rayleigh-Taylor driven, shock driven, and shear driven turbulence, and is shown to produce reasonable agreement with simulations and experiments in these scenarios.

42 ENGINEERING↗
DOE OSTI2022

Evolving Interior Mixing Regimes in a Tidal River Plume

Microstructure profiling was utilized to estimate vertical mixing (via vertical turbulent buoyancy flux) during a tidal pulse in the interior Merrimack River plume in calm winds. Multiple stratified shear mixing regimes appear and evolve with time. Initially the plume acts as a nearfield jet, with mixing in the plume (plume layer mixing) and over the plume-ambient interface (nearfield interfacial mixing). As the plume grows, interfacial mixing is suppressed offshore of the nearfield as currents slow, diminishing turbulent exchange between plume and shelf. At the end of ebb, ambient tidal currents reverse direction below plume, initiating another mode of internal, interfacial mixing (coined here as tidal interfacial mixing), allowing exchange between plume and ambient waters offshore. This work highlights previously unreported tidally modulated mixing within the near and midfield of a river plume.

54 ENVIRONMENTAL SCIENCES↗