Search NASA⌕ Search

SEARCH · Search NASA

Results for “Transition to 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Overview of the COMPASS results *

COMPASS addressed several physical processes that may explain the behaviour of important phenomena. This paper presents results related to the main fields of COMPASS research obtained in the recent two years, including studies of turbulence, L–H transition, plasma material interaction, runaway electron, and disruption physics: Tomographic reconstruction of the edge/SOL turbulence observed by a fast visible camera allowed to visualize turbulent structures without perturbing the plasma. Dependence of the power threshold on the X-point height was studied and related role of radial electric field in the edge/SOL plasma was identified. The effect of high-field-side error fields on the L–H transition was investigated in order to assess the influence of the central solenoid misalignment and the possibility to compensate these error fields by low-field-side coils. Results of fast measurements of electron temperature during ELMs show the ELM peak values at the divertor are around 80% of the initial temperature at the pedestal. Liquid metals were used for the first time as plasma facing material in ELMy H-mode in the tokamak divertor. Good power handling capability was observed for heat fluxes up to 12 MW m –2 and no direct droplet ejection was observed. Partial detachment regime was achieved by impurity seeding in the divertor. The evolution of the heat flux footprint at the outer target was studied. Runaway electrons were studied using new unique systems—impact calorimetry, carbon pellet injection technique, wide variety of magnetic perturbations. Radial feedback control was imposed on the beam. Forces during plasma disruptions were monitored by a number of new diagnostics for vacuum vessel (VV) motion in order to contribute to the scaling laws of sideways disruption forces for ITER. Current flows towards the divertor tiles, incl. possible short-circuiting through PFCs, were investigated during the VDE experiments. The results support ATEC model and improve understanding of disruption loads.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Transition in ICF Capsule Implosions

Longstanding design and reproducibility challenges in inertial confinement fusion (ICF) capsule implosion experiments involve recognizing the need for appropriately characterized and modeled three-dimensional initial conditions and high-fidelity simulation capabilities to predict transitional flow approaching turbulence, material mixing characteristics, and late-time quantities of interest—e.g., fusion yield. We build on previous coarse graining simulations of the indirect-drive national ignition facility (NIF) cryogenic capsule N170601 experiment-a precursor of N221205 which resulted in net energy gain. We apply effectively combined initialization aspects and multiphysics coupling in conjunction with newly available hydrodynamics simulation methods, including directional unsplit algorithms and low Mach-number correction-key advances enabling high fidelity coarse grained simulations of radiation-hydrodynamics driven transition. Our presentation includes discussion of the capsule initialization and implosion dynamics, analysis of the vorticity production budget, transition signatures, quantities of interest—late-time ion temperature and fusion-neutron yield, numerical uncertainty quantification, and comparisons with NIF data.

42 ENGINEERING↗
DOE OSTI2022

Evaluation of entropy driven jet symmetry transitions

Here we determine whether entropy drives planar turbulent jets into round turbulent jets. Determining when a jet flow transitions from one symmetry to the next is an important but incompletely resolved problem. The constructal view argues that the transition between symmetries of jet flows is governed by the minimization or maximization of entropy. Here we explore whether entropy increases with the transition of a planar turbulent jet into a round turbulent jet and whether entropy maximization (or minimization) predicts the same location of symmetry transition as velocity matching. We find that entropy considerations presented do not predict this transition.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Turbulence Structures in the Very Stable Boundary Layer Under the Influence of Wind Profile Distortion

Abstract In very stable boundary layers (VSBL), a “cocktail” of submeso motions routinely result in elevated mean wind speed maxima above the ground, acting as a new source of turbulence generation. This new source of turbulent kinetic energy enhances turbulent mixing and causes mean wind profile distortion (WPD). As a results, this transient distortion in the wind profile adjusts the classical log‐law. Addressing how WPD‐induced turbulence regulates flow structures, turbulent fluxes, and transitions in stability regimes across layers remains a challenge. Eddy covariance data measured at four levels on a 62‐m tower are employed to address these questions. It is shown that the WPD initiates large turbulent eddies that penetrate downward, leading to enhanced vertical mixing and comparable turbulent transport efficiencies across layers. As a consequence, turbulence intensity and fluxes are increased. As the WPD is intensified, turbulent fluxes and turbulent flux transport caused by large eddies are also enhanced, leading to a transition from very stable to weakly stable regimes. Due to the influence of WPD‐induced large eddies, the large‐eddy turbulent Prandtl number does not deviate appreciably from unity and the partitioning between turbulent kinetic and potential energies is linearly related to the gradient Richardson number.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

Improved representation of laminar and turbulent sheet flow in subglacial drainage models

Subglacial hydrology models struggle to reproduce seasonal drainage patterns that are consistent with observed subglacial water pressures and surface velocities. We modify the standard sheet-flow parameterization within a coupled sheet–channel subglacial drainage model to smoothly transition between laminar and turbulent flow based on the locally computed Reynolds number in a physically consistent way (the ‘transition’ model). We compare the transition model to standard laminar and turbulent models to assess the role of the sheet-flow parameterization in reconciling observed and modelled water pressures under idealized and realistic forcing. Relative to the turbulent model, the laminar and transition models improve seasonal simulations by increasing winter water pressure and producing a more prominent late-summer water pressure minimum. In contrast to the laminar model, the transition model remains consistent with its own internal assumptions across all flow regimes. Based on the internal consistency of the transition model and its improved performance relative to the standard turbulent model, we recommend its use for transient simulations of subglacial drainage.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Flow-shear destabilization of multiscale electron turbulence

The impact of sheared ${\mathbf{E} \hskip -1pt \times \hskip -1pt \mathbf{B}}$ flow on multiscale turbulence is studied with nonlinear gyrokinetic simulations. Simulations are based on DIII-D-like, high-confinement mode (H-mode) pedestal parameters in the regime of low ion temperature gradient drive, where there is a broad spectrum of electron temperature gradient (ETG)-driven turbulence. It is found that ${\mathbf{E} \hskip -1pt \times \hskip -1pt \mathbf{B}}$ shear can have a significant effect on ETG-driven electron transport, with an unexpected transition from a turbulence stabilization regime at moderate to large shearing rates $γ_E$ to a novel turbulence destabilization regime at low levels of $γ_E$. In the turbulence stabilization regime, the electron energy flux decreases monotonically with $γ_E$, even when $γ_E$ is small compared to electron mode growth rates. The stabilizing effect comes dominantly from the electron, not ion, gyrokinetic equation. In the novel destabilization regime, reduction of zonal energy results from the interaction of $γ_E$-modulated nonlinear drive in the zonal ion gyrokinetic equation, increasing the electron transport over a broad range of shearing rates. Finally, neither of these effects has been observed in previous electron-scale simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Observation of a reduced-turbulence regime with boron powder injection in a stellarator

Abstract In state-of-the-art stellarators, turbulence is a major cause of the degradation of plasma confinement. To maximize confinement, which eventually determines the amount of nuclear fusion reactions, turbulent transport needs to be reduced. Here we report the observation of a confinement regime in a stellarator plasma that is characterized by increased confinement and reduced turbulent fluctuations. The transition to this regime is driven by the injection of submillimetric boron powder grains into the plasma. With the line-averaged electron density being kept constant, we observe a substantial increase of stored energy and electron and ion temperatures. At the same time, the amplitude of the plasma turbulent fluctuations is halved. While lower frequency fluctuations are damped, higher frequency modes in the range between 100 and 200 kHz are excited. We have observed this regime for different heating schemes, namely with both electron and ion cyclotron resonant radio frequencies and neutral beams, for both directions of the magnetic field and both hydrogen and deuterium plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Zooming in: SCREAM at 100 m using regional refinement over the San Francisco Bay Area

Pushing global climate models to large-eddy simulation (LES) scales over complex terrain has remained a major challenge. This study presents the first known implementation of a global model – SCREAM (Simple Cloud-Resolving E3SM Atmosphere Model) – at 100 m horizontal resolution using a regionally refined mesh (RRM) over the San Francisco Bay Area. Two hindcast simulations were conducted to test performance under both strong synoptic forcing and weak, boundary-layer-driven conditions. We demonstrate that SCREAM can stably run at LES scales while realistically capturing topography, surface heterogeneity, and coastal processes. The 100 m SCREAM-RRM substantially improves near-surface wind speed, temperature, humidity, and pressure biases compared to the baseline 3.25 km simulation, and better reproduces fine-scale wind oscillations and boundary-layer structures. These advances leverage SCREAM's scale-aware SHOC turbulence parameterization, which transitions smoothly across scales without tuning. Performance tests show that while CPU-only simulations remain costly, GPU acceleration with SCREAMv1 on NERSC's Perlmutter system enables two-day hindcasts to complete in under two wall-clock days. Our results open the door to LES-scale studies of orographic flows, boundary-layer turbulence, and coastal clouds within a fully comprehensive global modeling framework.

Geosciences↗
DOE OSTI2023

Comparison of nanocomposite dispersion and distribution for several melt mixers

Breakup (dispersion) and distribution of nanoparticles are the chief hurdles towards taking advantage of nanoparticles in polymer nanocomposites for reinforcement, flame retardancy, conductivity, chromaticity, and other properties. Microscopy is often used to quantify mixing, but it has a limited field of view, does not average over bulk samples, and fails to address nano-particle hierarchical structures. Ultra-small-angle X-ray scattering (USAXS) can provide a macroscopic statistical average of nanoscale dispersion (breakup) and emergent hierarchical structure, as well as the distribution on the nanoscale. Here, this work compares several common mixer geometries for carbon black-polystyrene nanocomposites. Two twin-screw extruder geometries, typical for industrial processing of melt blends, are compared with a laboratory-scale single screw extruder and a Banbury mixer. It is found that for a given mixer, nanoscale distribution increases following a van der Waals function using accumulated strain as an analogue for temperature while macroscopic distribution/dispersion, using microscopy, does not follow this dependency. Breakup and aggregation in dispersive mixing follow expected behavior on the nanoscale. Across these drastically different mixing geometries an unexpected dependency is observed for nanoscale distributive mixing (both nano and macroscopic) as a function of accumulated strain that may reflect a transition from distributive turbulent to dispersive laminar mixing as the mixing gap is reduced.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The influence of the shock-to-reshock time on the Richtmyer–Meshkov instability in reshock

Experiments on the Richtmyer–Meshkov instability (RMI) in a dual driver vertical shock tube (DDVST) are described. An initially planar, stably stratified membraneless interface is formed by flowing air from above and sulfur hexafluoride from below the interface location using the method of Jones & Jacobs ( Phys. Fluids , vol. 9, issue 1997, 1997, pp. 3078–3085). A random three-dimensional, multi-modal initial perturbation is imposed by vertically oscillating the gas column to produce Faraday waves. The DDVST design generates two shock waves, one originating above and one below the interface, with these shocks having independently controllable strengths and interface arrival times. The shock waves have nominal strengths of $M_L=1.17$ and $M_H=1.18$ for the shock wave originating in the light and heavy gas, respectively, with these strengths chosen to result in arrested bulk interface motion following reshock. The influence of the length of the shock-to-reshock time, as well as the order of shock arrival, on the post-reshock RMI is examined. The mixing layer width grows according to $h\propto t^\theta$ , where $\theta _H=0.36\pm 0.018$ (95 %) and $\theta _L=0.38\pm 0.02$ (95 %) for heavy and light shock first experiments, respectively, indicating no strong dependence on the order of shock wave arrival. Volume integrated specific turbulent kinetic energy (TKE) in the mixing layer versus time is found to decay according to $E_{tot}/\bar {\rho }\propto t^p$ with $p_H=-0.823\pm 0.06$ (95 %) and $p_L=-1.061\pm 0.032$ (95 %) for heavy and light shock first experiments, respectively. Notably, the 95 % confidence intervals do not overlap. Analysis on the influence of the shock-to-reshock time on turbulent length scales, transition criteria, spectra and mixing layer anisotropy are also presented.

Ferguson, Kevin (ORCID:0000000166537482)↗
DOE OSTI2024

Anticathode effect on electron kinetics in electron beam generated E × B plasma

Electron beam (e-beam) generated plasmas with applied crossed electric and magnetic (E x B) fields are promising for low damage processing of materials with applications to microelectronics and quantum information systems. In cylindrical e-beam E x B plasmas, radial confinement of electrons and ions is achieved by an axial magnetic field and radial electric field, respectively. To control the axial confinement of electrons, such e-beam generated plasma sources may incorporate a conducting boundary known as an anticathode, which is placed on the axially opposite side of the plasma from the cathode. In this work, it is shown that varying the anticathode voltage bias can control the degree to which the anticathode collects or repels incident electrons, allowing control of warm electron (electron energies in 10–30 eV range) and beam electron population confinement. It is suggested that the effect of the anticathode bias on the formation of these distinct electron populations is also associated with the transition between weak turbulence and strong Langmuir turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Turbulence and E × B flow correlations across the L-H transition in DIII-D deuterium and hydrogen plasmas

The isotope dependence of the low- to high confinement-mode (L- to H-mode) transition power threshold P LH presents significant challenges for the initial (non-nuclear) hydrogen operations phase of the international thermonuclear experimental reactor (ITER). Here, we examine the isotope dependence of turbulence and E × B flow correlation properties in the L-mode edge plasma, leading up to the L-H transition. At marginal auxiliary power (near P LH ), turbulence is initially suppressed periodically during limit cycle oscillations (LCO) that precede the transition to sustained H-mode confinement. Here, we present evidence that the long-range (toroidal) correlation of the E × B edge plasma flow across the LCO phase is much weaker in hydrogen than in deuterium in the DIII-D tokamak, congruent with the higher threshold power P LH required to access LCO and H-mode in hydrogen. Concomitantly, the time required to initially quench edge turbulence via localized edge E × B flow shear is significantly longer in hydrogen (1–1.5 ms) than in deuterium (approx. 100 µs). No toroidal long-range correlation of the turbulence amplitude is observed, in agreement with expectations based on the relatively short poloidal turbulence correlation length. Radial edge turbulence and flow correlation lengths are longer in deuterium than in hydrogen plasmas as one would expect from ‘naïve’ gyro-Bohm isotope transport scaling, despite the substantially higher thermal flux across the last closed flux surface in hydrogen before the L-H transition. This article is part of a discussion meeting issue 'H-mode transition and pedestal studies in fusion plasmas'.

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

Assessing turbulence and mixing parameterizations in the gray-zone of multiscale simulations over mountainous terrain during the METEX21 field experiment

Multiscale numerical weather prediction models transition from mesoscale, where turbulence is fully parameterized, to microscale, where the majority of highly energetic scales of turbulence are resolved. The turbulence gray-zone is situated between these two regimes and multiscale models must downscale through these resolutions. Here, we compare three multiscale simulations which vary by the parameterization used for turbulence and mixing within the gray-zone. The three parameterizations analyzed are the Mellor-Yamada Nakanishi and Niino (MYNN) Level 2.5 planetary boundary layer scheme, the TKE-1.5 large eddy simulation (LES) closure scheme, and a recently developed three-dimensional planetary boundary layer scheme based on the Mellor-Yamada model. The simulation domain includes complex (i.e., mountainous) terrain in Nevada that was instrumented with meteorological towers, profiling and scanning lidars, a tethered balloon, and a surface flux tower. Simulations are compared to each other and to observations, with assessment of model skill at predicting wind speed, wind direction and TKE, and qualitative evaluations of transport and dispersion of smoke from controlled releases. This analysis demonstrates that microscale predictions of transport and dispersion can be significantly influenced by the choice of turbulence and mixing parameterization in the terra incognita , particularly over regions of complex terrain and with strong local forcing. This influence may not be apparent in the analysis of model skill, and motivates future field campaigns involving controlled tracer releases and corresponding modeling studies of the turbulence gray-zone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A unified theory of zonal flow shears and density corrugations in drift wave turbulence

A unified theory of zonal flow shears and density corrugations in drift wave turbulence is presented. Polarization and density advection beat excitation are studied in combination with modulational response. Noise is driven by two-time flux correlation. While the effective zonal flow eddy viscosity can go negative, the zonal diffusivity is positive definite. There is no inverse cascade of density corrugation. The connection between avalanches and corrugations is discussed. Here, the zonal cross-correlation is identified and calculated. Conditions for alignment of zonal shears and corrugation gradients are determined, and the implications for staircase structure are discussed. We show that the synergy of beat noise and modulational effects is stronger than either alone. Strong zonal flows can be excited well below the modulational instability threshold. In the context of L–H transition, zonal noise quenches turbulence overshoot by eliminating the threshold for zonal flow excitation. The power threshold for L–H transition is lowered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Generation of momentum transport in weakly turbulent β -plane magnetohydrodynamics

Magnetohydrodynamic turbulence on a β plane with an in-plane mean field, a system which serves as a simple model for the solar tachocline, is investigated analytically and computationally. We first derive two useful analytic constraints: We express the mean turbulent cross-helicity in terms of the mean turbulent magnetic energy, and then show that (for weak turbulence) the time-averaged momentum transport in the system can be expressed in terms of the cross-helicity spectrum. Here, we then complete a closure of the system using weak turbulence theory, appropriately extended to a system with multiple interacting eigenmodes. We use this closure to perturbatively solve for the spectra at lowest order in the Rossby parameter β and thereby show that the momentum transport in the system is O ⁡(β 2 ), thus quantifying the transition away from Alfvénized turbulence. Finally, we verify our theoretical results by performing direct numerical simulations of the system over a broad range of β.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Some observations concerning "laminarization" in heated vertical tubes

For low "turbulent" Reynolds numbers in strongly-heated vertical gas flow, fundamental results of existing direct numerical simulation (DNS) databases have been examined to deduce differences between cases which "revert" to turbulent and those that yield integral parameters which correspond to laminar flow (henceforth laminarizing or laminarized). The objectives of this paper are (1) to examine the streamwise evolution of near-wall flow structures and (2) to determine which, if any, proposed laminarization parameters could be used to discriminate between turbulent and laminarizing flows resulting. Views of streamlines in r-θ cross sections showed evidence of a ring of irregular-shaped streamwise vortices near the wall at all locations. The transverse advection of these vortices is hypothesized to provide a path of least resistance to the transport of streamwise momentum in the wall-normal direction. It is demonstrated that a steady streamwise laminar vortex can provide apparent turbulent momentum transport (as a consequence of the Reynolds decomposition) such that a reasonable mean velocity profile is predicted near the wall. For the present cases, the values of the non-dimensional radius (y c,w + = Re τ,w ) and of the modified wall Reynolds number did discriminate whether turbulent or laminarizing flow resulted.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2020

Theory of the Tertiary Instability and the Dimits Shift from Reduced Drift-Wave Models

Tertiary modes in electrostatic drift-wave turbulence are localized near extrema of the zonal velocity U(x) with respect to the radial coordinate x. We argue that these modes can be described as quantum harmonic oscillators with complex frequencies, so their spectrum can be readily calculated. The corresponding growth rate $γ$ TI is derived within the modified Hasegawa-Wakatani model. We show that $γ$TI equals the primary-instability growth rate plus a term that depends on the local U''; hence, the instability threshold is shifted compared to that in homogeneous turbulence. This provides a generic explanation of the well-known yet elusive Dimits shift, which we find explicitly in the Terry-Horton limit. Linearly unstable tertiary modes either saturate due to the evolution of the zonal density or generate radially propagating structures when the shear |U'| is sufficiently weakened by viscosity. The Dimits regime ends when such structures are generated continuously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗