Search NASA⌕ Search

SEARCH · Search NASA

Results for “Turbulent Heating”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Numerical investigation of two-dimensional buoyancy-driven eddies in liquid metal magnetohydrodynamic flows in breeding blankets *

Abstract Lead–lithium flows are key features in the design of tokamak breeding blanket concepts such as the dual-coolant lead–lithium (DCLL). Since they flow under magnetic fields, they are affected by magnetohydrodynamic (MHD) effects. The neutron flux originating in the tokamak plasma heats the breeding blanket channels in a non-uniform manner, inducing buoyancy forces in the liquid metal. Buoyancy may become a source of quasi-two-dimensional (Q2D) turbulence, and the appearance of eddies may affect the transport of heat and tritium across the blanket. Blankets characterized by high-speed liquid metal flows (such as DCLL) will need ceramic insulating walls to reduce the MHD-related pressure drop in the channels. In our simulations, we have used the Q2D model proposed by Sommeria and Moreau (SM82) which is especially suitable for modeling electrically insulating channel flows. On top of that, we have modelled buoyancy forces in the momentum equation using the Oberbeck–Boussinesq approximation. In this work, we include a validation of the implemented Q2D model in buoyancy-driven cases and identify a thermohydraulic configuration that promotes the generation of eddies and the accumulation of tritium. We also show the results provided by our post-processing tool based on the bi-dimensional fast Fourier transform for eddy detection and characterization. We complete our investigation by performing an initial assessment of how eddies can accumulate tritium in the breeding blanket. We conclude by discussing the relevance of preventing accumulation of tritium within the blanket and suggest a possible solution.

Suarez, Daniel (ORCID:0000000251740100)↗
DOE OSTI2020

Evidence and modeling of turbulence bifurcation in L-mode confinement transitions on Alcator C-Mod

Analysis and modeling of rotation reversal hysteresis experiments show that a single turbulent bifurcation is responsible for the Linear to Saturated Ohmic Confinement (LOC/SOC) transition and concomitant intrinsic rotation reversal on Alcator C-Mod. Plasmas on either side of the reversal exhibit different toroidal rotation profiles and therefore different turbulence characteristics despite the profiles of density and temperature, which are indistinguishable within measurement uncertainty. Elements of this bifurcation are also shown to persist for auxiliary heated L-modes. The deactivation of subdominant (in the linear growth rate and contribution to heat transport) ion temperature gradient and trapped electron mode instabilities is identified as the only possible change in turbulence within a reduced quasilinear transport model across the reversal, which is consistent with the measured profiles and inferred heat and particle fluxes. Experimental constraints on a possible change from strong to weak turbulence, outside the description of the quasilinear model, are also discussed. These results indicate an explanation for the LOC/SOC transition that provides a mechanism for the hysteresis through the dynamics of subdominant modes and changes in their relative populations and does not involve a change in the most linearly unstable ion-scale drift-wave instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Theory for transport in magnetized plasmas

We present a two fluid model for anomalous thermal, particle and momentum transport which has been derived from kinetic theory using a small parameter ε~10-2. This small parameter can come from either the ratio of wave frequency to cyclotron frequency or the ratio of potential to thermal energies. This fluid model is strongly nonlinear and has an exact closure within the accuracy ε. It can be used to simulate the full radius of a tokamak including internal and edge transport barriers at the same time without introducing any local assumptions for physics model or grid size. It involves a full transport matrix, thus giving possibility for both density and heat pinches. The model can be derived from the Vlasov equation, adding also close collisions and includes particle pinches of the type found in the magnetosphere and in the Massachusetts Institute of Technology elevated dipole machine. Pinches are also in agreement for the heat pinch experiments on Doublet III D including simultaneous particle and heat pinches. For tokamaks it gives the right scaling of confinement time versus heating power. Finally, we have also recovered the Turbulent Equipartition result that the density scales as inverse safety factor q in typical simulations. Thus this model is both fundamental and of high practical use.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Modeling electron temperature profiles in the pedestal with simple formulas for ETG transport

Abstract This paper reports on the refinement (building on (Hatch D.R. et al 2022 Phys. Plasmas 29 062501)) and application of simple formulas for electron heat transport from electron temperature gradient (ETG) driven turbulence in the pedestal. The formulas are improved by (1) improving the parameterization for certain key parameters and (2) carefully accounting for the impact of geometry and shaping in the underlying gyrokinetic simulation database. Comparisons with nonlinear gyrokinetic simulations of ETG transport in the MAST pedestal demonstrate the model’s applicability to spherical tokamaks in addition to standard aspect ratio tokamaks. We identify bounds for model applicability: the model is accurate in the steep gradient region, where the ETG turbulence is largely slab-like, but accuracy decreases as the temperature gradient becomes weaker in the pedestal top. We use the formula to model the electron temperature profile in the pedestal for four experimental scenarios while extensively varying input parameters to represent uncertainties. In all cases, the predicted electron temperature profile exhibits extreme sensitivity to separatrix temperature and density, which has implications for core-edge integration. The model reproduces the electron temperature profile for high η e = L n e / L T e scenarios but not for low η e scenarios in which microtearing modes (MTMs) have been identified. We develop a proof-of-concept model for MTM transport and explore the relative roles of ETG and MTM in setting the electron temperature profile.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Core electron temperature turbulence and transport during sawtooth oscillations in the DIII-D tokamak

Abstract Sawteeth are one of the concerning instabilities in ITER and future burning plasma experiments. Sawtooth dynamics and its interaction with broadband plasma turbulence has been a challenge for predictive simulations of core transport in future fusion devices. This study provides new observations of core turbulence behavior during sawtooth oscillations in DIII-D hydrogen L-mode neutral beam injection heated plasmas in an inner wall limited configuration. A strong correlation of electron temperature and density turbulence levels with the sawtooth oscillation phase has been observed at locations inside the T e inversion radius and/or safety factor q = 1 magnetic surface. The T e turbulence amplitude in the core during the sawtooth ramp exhibits a critical T e gradient behavior inside but not near the T e inversion radius/ q = 1 magnetic surface. The most unstable mode calculated from the trapped gyro-landau fluid turbulence simulations reveal a change from low- k ion-type to low-k electron-type modes from pre- to post- sawtooth crash time periods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Measurement of the Alfvén Wave Parametric Decay Instability Growth Rate

Alfvén waves, a fundamental mode of magnetized plasmas, are ubiquitous in space and laboratory plasmas. The nonlinear behavior of these modes is thought to play a key role in important problems in space plasma, such as the heating of the solar corona and solar wind turbulence. In particular, theoretical predictions show that these Alfvén waves may be unstable to various parametric instabilities, but space observations of these processes are limited. We demonstrate the first measurement of the Alfvén wave parametric decay instability (PDI) growth rate. Experiments are conducted on the Large Plasma Device at UCLA in which a high amplitude 𝛿⁢𝐵/𝐵 0 ∼ 0.7% pump Alfvén wave is launched from one end of the device and a smaller seed Alfvén wave is launched from the other side. When the frequency of the seed wave is chosen to match the backward wave expected from PDI, damping of the seed wave is reduced. We compare this reduction in damping to the theoretically expected PDI growth rate while accounting for acoustic mode damping. Results show agreement between measurements and theoretical predictions. As a result, this not only provides critical validation for PDI theories and simulations that could help interpret future space observations but also suggests a new way of studying similar nonlinear wave phenomena.

Alfvén waves↗
DOE OSTI2021

AmeriFlux US-Wwt Willamette Wheat

This is the AmeriFlux version of the carbon flux data for the site US-Wwt Willamette Wheat. Site Description - The site was established in summer 2014 and is part of the Oregon eddy covariance flux tower network . Meteorological variables such as temperature, humidity, solar irradiance and wind are measured at this site. Turbulent fluxes of water vapor, carbon dioxide and heat are observed and stored as 30 minute averages. The crop changes periodically and rye grass and fescue is grown alternately at this site located in Oregon's Willamette Valley with its relatively mild climate. The Wheat site is also part of the PhenoCam network operated by Harvard University and the University of New Hampshire.

Law, Bev↗
DOE OSTI2021

Morning boundary layer conditions for shallow to deep convective cloud evolution during the dry season in the central Amazon

Observations of the boundary layer (BL) processes are analyzed statistically for dry seasons of 2 years and in detail, as case studies, for 4 shallow convective days (ShCu) and 4 shallow-to-deep convective days (ShDeep) using a suite of ground-based measurements from the Observation and Modeling of the Green Ocean Amazon (GoAmazon 2014/5) Experiment. The BL stages in ShDeep days, from the nighttime to the cloudy mixing layer stage, are then described in comparison with ShCu days. Atmospheric thermodynamics and dynamics, environmental profiles, and surface turbulent fluxes were employed to compare these two distinct situations for each stage of the BL evolution. Particular attention is given to the morning transition stage, in which the BL changes from stable to unstable conditions in the early morning hours. Results show that the decrease in time duration of the morning transition on ShDeep days is associated with high humidity and well-established vertical wind shear patterns. Higher humidity since nighttime not only contributes to lowering the cloud base during the rapid growth of the BL but also contributes to the balance between radiative cooling and turbulent mixing during nighttime, resulting in higher sensible heat flux in the early morning. The sensible heat flux promotes rapid growth of the well-mixed layer, thus favoring the deeper BL starting from around 08:00 LST (UTC-4 h). Under these conditions, the time duration of morning transition is used to promote convection, having an important effect on the convective BL strength and leading to the formation of shallow cumulus clouds and their subsequent evolution into deep convective clouds. Statistical analysis was used to validate the conceptual model obtained from the case studies. Despite the case-to-case variability, the statistical analyses of the processes in the BL show that the described processes are very representative of cloud evolution during the dry season.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Comparative Thermal Performance of Downdraft and Updraft Forced Convection in a Representative 3x3 PWR Fuel Assembly

A plethora of theoretical, numerical, and experimental investigations have relied on updraft forced convection of operating fuel assemblies for light water-cooled fuel assemblies. The inertial scales of turbulence in updraft flow through PWR fuel assemblies are primarily influenced by upstream mixing patterns within the lower plenum and consequent acceleration of flow through the lower core and support plates into the bottom nozzles. In secondary heat transfer applications such as steam generation, recuperation and economization, downdraft flow is utilized to retrieve maximal sensible heat from the primary coolant. Sub-channel volumes aligned vertically could also benefit from force of gravity in developing turbulent flows. A large reduction in pumping power, associated cost of operation and maintenance, and improved balance of plant is suggested to be possible through downdraft forced convection in PWR-type reactors. This article explores the impact of downdraft forced convection within the well-studied 3x3 sub-channel within a range of 20 hydraulic diameters across a representative intermediate flow mixing grid. A range of inflow Reynolds number from 10,000 to 60,000 is simulated over a range of equivalent heat loads of 20% to 100% spanning 20 hydraulic diameters. The fuel rods are represented in true geometric detail, and a constant-value heat profile is assumed for the investigation. At each inflow condition, the turbulence intensity, peak velocity, local heat transfer coefficient, peak wall heat flux, and peak wall temperature are collected across the sub-channel and compared with an updraft configuration evaluated at the same conditions.

Rao, Vivek↗
DOE OSTI2020

Turbulence/flame/wall interactions in non-premixed inclined slot-jet flames impinging at a wall using direct numerical simulation

In this paper, three-dimensional turbulent non-premixed oblique slot-jet flames impinging at a wall were investigated using direct numerical simulation (DNS). Two cases are considered with the Damköhler number (Da) of case A being twice that of case B. A 17 species and 73-step mechanism for methane combustion was employed in the simulations. It was found that flame extinction in case B is more prominent compared to case A. Reignition in the lower branch of combustion for case A occurs when the scalar dissipation rate relaxes, while no reignition occurs in the lower branch for case B due to excessive scalar dissipation rate. A method was proposed to identify the flame quenching edges of turbulent non-premixed flames in wall-bounded flows based on the intersections of mixture fraction and OH mass fraction iso-surfaces. The flame/wall interactions were examined in terms of the quenching distance and the wall heat flux along the quenching edges. There is essentially no flame/wall interaction in case B due to the extinction caused by excessive turbulent mixing. In contrast, significant interactions between flames and the wall are observed in case A. The quenching distance is found to be negatively correlated with wall heat flux as previously reported in turbulent premixed flames. The influence of chemical reactions and wall on flow topologies was identified. The FS/U and FC/U topologies are found near flame edges, and the NNN/U topology appears when reignition occurs. The vortex-dominant topologies, FC/U and FS/S, play an increasingly important role as the jet turbulence develops.

42 ENGINEERING↗
DOE OSTI2024

Increasing Main Cooler Thermal Performance for sCO2 Power Cycles

Increasing sCO2 power cycle cooler performance will increase the efficiency and lower the cost of electricity for sCO2 power cycles. This work considered the heat transfer performance of two counter-flow, shell-and-tube heat exchangers using local measurements. The heat exchangers were constructed using commercial tubing. One heat exchanger was constructed with a conventional, smooth inner tube (7 mm nominal inner diameter). The other was constructed with an additively manufactured tube with a square cross section and angled rib turbulators. For mixed mean temperatures away from the pseudocritical temperature, the heat transfer coefficient matched the Dittus-Boelter correlation. As the mixed mean temperature approached the pseudocritical temperature, the heat transfer coefficient deviated from the Dittus-Boelter correlation due to transverse property variation and buoyancy. The local and average heat transfer coefficients for angled rib tubes were nominally 150% to 300% greater than the heat transfer coefficients for the smooth tube. Trends in heat exchanger effectiveness indicated that the ribs were most beneficial when the water flow rates were at the low end of the range considered: 0.016 – 0.126 kg/s. Using results from prior system modeling, it is anticipated that adding heat transfer enhancement features in a cycle cooler will increase the sCO2 cycle efficiency by 0.34 percentage points.

Searle, Matthew↗
DOE OSTI2021

Gyrokinetic benchmark of the electron temperature-gradient instability in the pedestal region

Transport from turbulence driven by the electron temperature-gradient (ETG) instability is likely a major source of electron heat losses through the pedestal. Due to extreme gradients and strong shaping, ETG instabilities in the pedestal are distinct from those in the core, having, for example, multiple branches (toroidal and slab) in different wavenumber ranges. Due to its importance for pedestal transport, and its rather exotic character, a rigorous multi-code benchmarking exercise is imperative. In this work, we describe such an exercise, wherein we have carried out a detailed comparison of local linear pedestal ETG simulations using three gyrokinetic codes, CGYRO, GEM, and GENE and testing different geometric parameters (such as circular, Miller, and equilibrium EFIT geometry). The resulting linear frequencies, growth rates, and eigenfunctions show very good agreement between the codes in the three types of employed geometries. A nonlinear benchmark between CGYRO and GENE is also described, exhibiting good agreement (a maximum of 20% difference in the heat fluxes computed) at two locations in the pedestal. This lays the foundation for confidently modeling ETG turbulence in the pedestal

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Global gyrokinetic simulation of neoclassical ambipolar electric field and its effects on microturbulence in W7-X stellarator

Global neoclassical simulations of a model equilibrium of the W7-X stellarator find an ambipolar electric field with either an ion root or an electron root solution depending on the temperature ratio between electrons and ions. The ambipolar electric field is then used as an equilibrium field in the turbulence simulations of ion temperature gradient (ITG) instability. The shear of the ambipolar electric field has modest effects on the ITG linear instability, nonlinear saturation amplitude, and turbulent transport in the ion root case. However, in the electron root case, the ambipolar electric field significantly reduces the linear ITG growth rate, turbulence intensity, and radial correlation length, resulting in an ion heat conductivity comparable to the neoclassical transport level in the strong shear region.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Dimits shift, avalanche-like bursts, and solitary propagating structures in the two-field flux-balanced Hasegawa–Wakatani model for plasma edge turbulence

We show that the recently introduced two-field flux-balanced Hasegawa–Wakatani (BHW) model captures the key features of drift-wave turbulent transport mediated by zonal flows observed in more complete and accurate gyrokinetic simulations, such as the existence of a nonlinear upshift of the threshold for drift wave turbulence driven transport, often called the Dimits shift, as well as non-local transport with avalanche bursts and solitary propagating structures. Because of the approximations made in the BHW model, these observations are made for the particle flux instead of the heat flux more commonly studied in ion temperature gradient (ITG) driven turbulence in fluid or gyrokinetic codes. Many of these features are not seen in other Hasegawa–Wakatani models, which confirm the critical role of the electron dynamics parallel to the magnetic field lines. To address questions regarding the role of boundary conditions on the drift-wave zonal flow dynamics, we apply our model to both a channel domain geometry and the more typical doubly periodic geometry. We only observe strong soliton-like solutions in the particle flux for the channel geometry, in the vicinity of the boundaries, where strong velocity shear and density gradients are generated, which are absent in the doubly periodic simulations. Changing the aspect ratio of the simulation domain also has a significant effect. In domains which are elongated in the radial direction, more complex multi-scale dynamics takes place, with multiple zonal jets interacting with each other, and large scale avalanches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Lagrangian and Eulerian Supersaturation Statistics in Turbulent Cloudy Rayleigh–Bénard Convection: Applications for LES Subgrid Modeling

Abstract Turbulent fluctuations of scalar and velocity fields are critical for cloud microphysical processes, e.g., droplet activation and size distribution evolution, and can therefore influence cloud radiative forcing and precipitation formation. Lagrangian and Eulerian water vapor, temperature, and supersaturation statistics are investigated in direct numerical simulations (DNS) of turbulent Rayleigh–Bénard convection in the Pi Convection Cloud Chamber to provide a foundation for parameterizing subgrid-scale fluctuations in atmospheric models. A subgrid model for water vapor and temperature variances and covariance and supersaturation variance is proposed, valid for both clear and cloudy conditions. Evaluation of phase change contributions through an a priori test using DNS data shows good performance of the model. Supersaturation is a nonlinear function of temperature and water vapor, and relative external fluxes of water vapor and heat (e.g., during entrainment-mixing and phase change) influence turbulent supersaturation fluctuations. Although supersaturation has autocorrelation and structure functions similar to the independent scalars (temperature and water vapor), the autocorrelation time scale of supersaturation differs. Relative scalar fluxes in DNS without cloud make supersaturation PDFs less skewed than the adiabatic case, where they are highly negatively skewed. However, droplet condensation changes the PDF shape response: it becomes positively skewed for the adiabatic case and negatively skewed when the sidewall relative fluxes are large. Condensation also increases correlations between water vapor and temperature in the presence of relative scalar fluxes but decreases correlations for the adiabatic case. These changes in correlation suppress supersaturation variability for the nonadiabatic cases and increase it for the adiabatic case. Implications of this work for subgrid microphysics modeling using a Lagrangian stochastic scheme are also discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Statistical analysis of magnetic divertor configuration influence on H-mode transitions

DIII-D plasmas are compared for two upper divertor configurations: with the outer strike point on the small angle slot (SAS) divertor target and with the outer strike point on the horizontal divertor target (HT). Scanning the vertical distance between the magnetic null point and the divertor target over a range 0.10–0.16 m is shown to increase the threshold power, $P_\mathrm{th}$, and edge plasma power, $P_\mathrm{Loss},$ for the low-to-high confinement (L–H) and H–L transitions respectively, by up to a factor of 1.4. The X-point height scans were performed at three L-mode core plasma line average electron densities, $\bar{n}_\mathrm{e} = $ 1.2, 2.2 and 3.6 $\times 10^{19}\,\mathrm{m}^{-3}$, to investigate the density dependence of divertor magnetic configuration influence on $P_\mathrm{th}$. The X-point height, $Z_\textrm{x-pt}$, was further extended across the range 0.16–0.22 m with the more open HT divertor configuration, for which a clear decrease in $P_\mathrm{th}$ with increasing $Z_\textrm{x-pt}$ is observed. The dependence of $P_\textrm{th}$ on divertor magnetic geometry is further investigated using a time-dependent probability density function (PDF) model and information geometry to elucidate the roles played by pedestal plasma turbulence and perpendicular velocity flows. The degree of stochasticity of the plasma turbulence is observed to be sensitive to the plasma heating rate. The calculated square of the information rate shows changes in the relative density fluctuations and perpendicular velocity PDFs begin 2–5 ms prior to the L–H transition for three plasmas; providing a crucial measurement of the dynamic timescale of external transport barrier formation. Additionally, both information length and rate provide potential predictors of the L–H transition for these plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Chiral effects and Joule heating in hot and dense matter

Initial states of dense matter with nonzero electron chiral imbalance could potentially give rise to strong magnetic fields through chiral plasma instability. Previous work indicated that unless chiral chemical potential is as large as the electron vector chemical potential, the growth of magnetic fields due to the instability is washed out by chirality flipping rate enabled by electron mass. We reexamine this claim in a broader range of parameters and find that at higher temperatures the hierarchy is reversed supporting a growing magnetic field for an initial electron chiral chemical potential much smaller than the electron vector chemical potential. Further, we identify a qualitatively new effect relevant for magnetized hot and dense medium where chiral magnetic effect (CME) sourced by density fluctuation acts as a powerful source of Joule heating. Remarkably, even modest chiral chemical potentials (keV) in such environment can deposit energy densities set by the QCD scale in a relatively short time of the order of a few milliseconds or seconds. We speculate how this mechanism makes CME-driven Joule heating a potentially critical ingredient in the dynamics of turbulent density fluctuation of supernovae and neutron star mergers.

Particle astrophysics↗
DOE OSTI2021

Diffraction of radio-frequency waves by plasma turbulence in the edge of a tokamak

The Mie-Lorenz-Debye theory of scattering of radio frequency waves (RF) by filamentary structures in the edge of a tokamak plasma has been established in prior publications. These results have been verified. Theory and computations show that RF can be reflected, refracted, and diffracted by turbulence. Furthermore, the filamentary structures lead to side scattering of waves and excitation of wave modes different from the launched RF. It is important to quantify the overall effect of turbulence on the propagation of RF since efficiency of heating and current drive by RF waves depends upon the wave characteristics. The commonly used RF ranges vary from the low frequency ion-cyclotron (IC) waves, to medium frequency lower-hybrid (LH) waves, to high-frequency electron cyclotron (EC) waves with wavelengths respectively, large, comparable, or smaller than the radial dimensions of the filaments. Currently, we discuss the scattering as well as the forces exerted by different RF in the edge plasma. Also, we present a code for 3D full wave scattering analysis of arbitrary tensor permittivity and plasma density profiles. The study is applicable to ITER-like plasmas, as well as to plasmas in medium sized tokamaks such as TCV, ASDEX-U and DIII-D and high field concepts like SPARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗