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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Exascale Multiphysics Nuclear Reactor Simulations for Advanced Designs

ENRICO is a coupled application developed under the U.S. Department of Energy's Exascale Computing Project (ECP) targeting the modeling of advanced nuclear reactors. It couples radiation transport with heat and fluid simulation, including the high-fidelity, highresolution Monte-Carlo code Shift and the Computational fluid dynamics code NekRS. NekRS is a highly-performant open-source code for simulation of incompressible and low-Mach fluid flow, heat transfer, and combustion with a particular focus on turbulent flows in complex domains. It is based on rapidly convergent high-order spectral element discretizations that feature minimal numerical dissipation and dispersion. State-of-the-art multilevel preconditioners, efficient high-order time-splitting methods, and runtime-adaptive communication strategies are built on a fast OCCA-based kernel library, libParanumal, to provide scalability and portability across the spectrum of current and future high-performance computing platforms. On Frontier, Nek5000/RS has recently achieved an unprecedented milestone in breaching over 1 billion spectral elements and 350 billion degrees of freedom. Shift has demonstrated the capability to transport upwards of 1 billion particles per second in full core nuclear reactor simulations featuring complete temperature-dependent, continuous-energy physics on Frontier. Shift achieved a weak-scaling efficiency of 97.8% on 8192 nodes of Frontier and calculated 6 reactions in 214,896 fuel pin regions below 1% statistical error yielding first-of-a-kind resolution for a Monte Carlo transport application.

Hamilton, Steven P.↗
DOE OSTI2022

Turbulent Rotating Rayleigh–Bénard Convection

Rotation with thermally induced buoyancy governs many astrophysical and geophysical processes in the atmosphere, ocean, sun, and Earth's liquid-metal outer core. Rotating Rayleigh–Bénard convection (RRBC) is an experimental system that has features of rotation and buoyancy, where a container of height H and temperature difference Δ between its bottom and top is rotated about its vertical axis with angular velocity Ω. The strength of buoyancy is reflected in the Rayleigh number (~ H 3 Δ) and that of the Coriolis force in the Ekman and Rossby numbers (~Ω -1 ). Rotation suppresses the convective onset, introduces instabilities, changes the velocity boundary layers, modifies the shape of thermal structures from plumes to vortical columns, affects the large-scale circulation, and can decrease or enhance global heat transport depending on buoyant and Coriolis forcing. RRBC is an extremely rich system, with features directly comparable to geophysical and astrophysical phenomena. Here we review RRBC studies, suggest a unifying heat transport scaling approach for the transition between rotation-dominated and buoyancy-dominated regimes in RRBC, and discuss non-Oberbeck–Boussinesq and centrifugal effects.

42 ENGINEERING↗
DOE OSTI2023

Isotope physics of heat and particle transport with tritium in JET-ILW type-I ELMy H-mode plasmas

As part the DTE2 campaign in the JET tokamak, we conducted a parameter scan in T and D-T complementing existing pulses in H and D. For the different main ion masses, type-I ELMy H-modes at fixed plasma current and magnetic field can have the pedestal pressure varying by a factor of 4 and the total pressure changing from <?CDATA $\beta_{\mathrm{N}} = 1.0$?> β N = 1.0 to 3.0. We investigated the pedestal and core isotope mass dependencies using this extensive data set. The pedestal shows a strong mass dependence on the density, which influences the core due to the strong coupling between both plasma regions. To better understand the causes for the observed isotope mass dependence in the pedestal, we analysed the interplay between heat and particle transport and the edge localised mode (ELM) stability. For this purpose, we developed a dynamic ELM cycle model with basic transport assumptions and a realistic neutral penetration. The temporal evolution and resulting ELM frequency introduce an additional experimental constraint that conventional quasi-stationary transport analysis cannot provide. Our model shows that a mass dependence in the ELM stability or in the transport alone cannot explain the observations. One requires a mass dependence in the ELM stability as well as one in the particle sources. The core confinement time increases with pedestal pressure for all isotope masses due to profile stiffness and electromagnetic turbulence stabilisation. Interestingly, T and D-T plasmas show an improved core confinement time compared to H and D plasmas even for matched pedestal pressures. For T, this improvement is largely due to the unique pedestal composition of higher densities and lower temperatures than H and D. With a reduced gyroBohm factor at lower temperatures, more turbulent drive in the form of steeper gradients is required to transport the same amount of heat. This picture is supported by quasilinear flux-driven modelling using TGLF -SAT2 within Astra . With the experimental boundary condition TGLF -SAT2 predicts the core profiles well for gyroBohm heat fluxes <?CDATA ${\gt}15$?> > 15 , however, overestimates the heat and particle transport closer to the turbulent threshold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Comparison of private flux region instability in conventional and super-X divertor configurations

Understanding turbulence in the divertor leg of tokamaks is essential to predicting the heat deposition profile on the divertor plate. This in turn is important for evaluating advanced divertor configurations, such as the super-X divertor. Within the divertor region, the private flux region is of interest because it is relatively unaffected by turbulence extending from the outboard midplane, so instabilities in this region could have a particularly pronounced effect on transport. These instabilities are modeled using the Arbitrary Topology Equation Reader (ArbiTER) eigenvalue code. Eigenmodes are examined further by comparing physics models to determine the fundamental mechanisms behind their formation, and quantifying the effect of individual terms. This analysis is conducted on both conventional and super-X divertors to compare these effects. Here, the resulting analysis reveals the presence of a geodesic curvature driven instability that is significantly more pronounced in the super-X configuration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Density fluctuation statistics and turbulence spreading at the edge of L–mode plasmas

Long-wavelength density fluctuations ($k{\rho _i}$ <1) are studied using beam emission spectroscopy (BES) at the edge of DIII-D L-mode plasmas (ρ = 0.88–1.1) in scenarios with electron cyclotron heating (ECH) power ramp (P ECH up to 1.5 MW), neutral beam injection (NBI) power ramp (P NBI up to 2.5 MW), and injected torque scan (-1 < T inj <0.6 Nm). We find that broadband turbulent density fluctuations (ƒ ~ 20–120 kHz) have a non-Gaussian distribution. The skewness of $\delta n/n$ changes sign from negative at ρ < 0.95–0.97 to positive at ρ > 0.97, indicating the prevalence of density 'voids' at inner radii and density 'blobs' at outer radii and outside of the separatrix. The turbulence intensity flux $\left\langle {{{\tilde v}_{\text{r}}}{{\tilde n}^2}} \right\rangle$ is calculated to characterize turbulence spreading at the plasma edge. During ECH/NBI power ramps and at counter-I p injected torque, $\left\langle {{{\tilde v}_{\text{r}}}{{\tilde n}^2}} \right\rangle$ is directed inward inside the separatrix, which is evidence of inward spreading of turbulence intensity from the edge gradient region caused by the inner propagation of density 'voids'. Significantly weaker $\left\langle {{{\tilde v}_{\text{r}}}{{\tilde n}^2}} \right\rangle$ is observed with co-I p torque. A correlation between co-I p torque, turbulence intensity $\delta n/n$ at ρ = 0.97, and increased srape-off layer (SOL) heat flux decay length ${\lambda _q}$ is found in the torque scan scenario, showing that edge turbulence plays a material role in determining the SOL conditions and heat flux width.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Convection in Scaled Turbine Internal Cooling Passages With Additive Manufacturing Roughness

Abstract Additive manufacturing processes, such as direct metal laser sintering (DMLS), enable the creation of novel turbine cooling internal passages and systems. However, the DMLS method produces a significant and unique surface roughness. Previous work in scaled passages analyzed pressure losses and friction factors associated with the rough surfaces, as well as investigated the velocity profiles and turbulent flow characteristics within the passage. In this study, the heat transfer characteristics of scaled additively manufactured surfaces were measured using infrared (IR) thermography. Roughness panels were CNC machined from plates of aluminum 6061 to create near isothermal roughness elements when heated. Fluid resistance differences between the aluminum roughness panels and roughness panels constructed from ABS plastic using the same roughness patterns from McClain et al. (2020, “Flow in a Simulated Turbine Blade Cooling Channel With Spatially Varying Roughness Caused by Additive Manufacturing Orientation,” ASME Turbo Expo 2020, Turbomachinery Technical Conference and Exposition, Virtual Conference, Sept. 21–25, GT2020-16069) were investigated. Finally, the overall thermal performance enhancements and friction losses were assessed through the calculation of surface averaged “global thermal performance” ratios. The global thermal performance characterizations indicate results in-line with those found for traditional commercial roughness and slightly below traditional internal passage convection enhancement methods such as swirl chambers, dimples, and ribs. The passages investigated in this study do not include compressibility effects or the long-wavelength artifacts and channel geometric deviations observed by Wildgoose et al. (2020, “Impact of Additive Manufacturing on Internal Cooling Channels with Varying Diameters and Build Directions,” ASME Turbo Expo 2020, Turbomachinery Technical Conference and Exposition, Virtual Conference, Sept. 21–25, GT2020-15049). However, the results of this study indicate that, based on the roughness augmentation alone, artificial convective cooling enhancers such as turbulators or dimples may still be required for additively manufactured turbine component cooling.

Engineering↗
DOE OSTI2023

New scaling law for turbulent boundary layers with high surface mass transfer

A fundamental understanding of the phenomena occurring in the turbulent boundary layer in the presence of surface blowing is limited, and considerable disagreements persist even in describing primary quantities, such as the boundary layer profile. The theories based on the linear boundary layer equations show that the thickness of the sublayer increases in the presence of surface blowing; therefore, the viscous sublayer and law of the wall modify. In this study, direct numerical simulations (DNS) of turbulent boundary layers with uniform surface mass transfers are carried out in order to scale the velocity profile. Emphasis is placed on moderate to high mass transfer rates, which are relevant to the most common hybrid rockets configuration. DNS data are used to establish a functional law of the wall and a law of wake by means of the relation between the wall shear stress and surface mass transfer. Analysis of the mean kinetic energy budget shows that the magnitude of turbulent kinetic energy increases by surface mass transfer, and the production rate extends significantly in the inner layer as the injection rate increases. DNS data of various surface blowing are used to complete the closure of turbulence kinetic energy equation and develop an eddy viscosity model. The predicted turbulent kinetic energy and eddy viscosity agree with DNS data for moderate to high blowing rates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Mechanisms and models of the turbulent boundary layers at transcritical conditions

Computational models for high-pressure transcritical turbulence in wall-modeled large-eddy simulation typically rely on wall-function models to overcome the need for resolving the boundary layer structures. However, the mechanisms and models of turbulent boundary layer at transcritical conditions remain poorly understood since the near-wall flow and heat transfer are significantly affected by the enhanced fluctuations and steep gradients in thermodynamic properties. Here, to address this issue, we study the mechanisms of transcritical turbulent boundary layers and wall-attached models at transcritical conditions. It is shown that the real-fluid variable-property effects are associated with the coupling between the near-wall cycle and the wall-normal coherent motions in the outer layer, resulting in the amplification of turbulent energy in the outer layer and noticeable energy transfer between the log-layer and the outer layer; hence, turbulence in the log-layer is modulated by the outer layer. Based on this underlying physical principle, we propose the characteristic velocity and length scales for the attached eddy at transcritical conditions, and extend the attached eddy model to transcritical turbulent boundary layers by introducing a mixed scaling that incorporates both inner and outer scalings. We show that the new characteristic scales and extended attached-eddy model perform well in characterizing the structures in transcritical turbulent boundary layers.

Li, Fangbo↗
DOE OSTI2021

Characterization of hollow cathode plasma turbulence using coherent Thomson scattering

Turbulence in hollow cathodes used for space propulsion is believed to play an important role in anomalous electron transport and ion heating. In this work, the implementation of coherent Thomson scattering to identify and characterize MHz-frequency ion acoustic turbulence and kHz-frequency oscillations in the plume of a hollow cathode is achieved. In the presence of a background magnetic field of a Hall thruster, a number of unstable modes are observed. A directive ion acoustic mode propagating predominantly within a restricted angle around the magnetic field is found, exhibiting an energy scaling with wavenumber k of the form k -5.2±0.58 , which differs from the classic Kadomtsev k -3 scaling for unmagnetized conditions. Bi-directional ion acoustic mode fluctuations propagating over a range of angles with respect to the magnetic field have been measured, possibly signifying the existence of a large-amplitude plasma wave, similar to the Buneman instability. Finally, electron density fluctuations in the kHz-frequency range, a possible consequence of drift-driven instabilities in the plane perpendicular to the magnetic field, have also been identified. Furthermore, these results not only are an indication of the diversity of wave types that exist in hollow cathode plumes but also point to the key role played by the presence of, and the configuration of, the magnetic field in their appearance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Convectively Induced Secondary Circulations and Wind‐Driven Heat Fluxes in the Surface Energy Balance Over Land

Increased resolution has enabled kilometer-scale weather and climate models to partially resolve secondary circulations, including horizontal convective rolls (HCRs) and cold pool gust fronts. Although these circulations are ubiquitous in convective boundary layers over land, their impacts on the surface energy balance are largely unknown. Doppler lidar and surface observations were combined with DOE E3SM land model experiments, revealing increased surface winds (5 m/s) and heat fluxes (50 W/m 2 ) in convergent branches of HCRs. Larger wind-driven flux responses (up to 150 W/m 2 ) were found along gust fronts. Surface energy balance shifts to accommodate wind-driven fluxes, reducing ground heat conduction and longwave cooling. Our findings from the US Southern Great Plains are broadly relevant to modeling convective boundary layers. In particular, widely used subgrid wind gust parameterizations were found to be physically inconsistent with resolved secondary circulations and could worsen climate prediction biases at kilometer-scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Simulation of Channel Flow with Square Ribs for Blanket First-Wall Cooling: Geometry-Specific Tuning of k-ω Model Using Adjoint Method

Cooling of the plasma-facing first wall is challenging in the design of blanket components because of the high heat flux (on the order of 𝑀𝑊/𝑚2) from the plasma, especially when a low thermal mass medium like helium is chosen as the coolant. Therefore, heat transfer enhancement in which the convective heat transfer rate is augmented by the addition of turbulence-promoting structures becomes a key initiative for providing sufficient cooling capability with helium. Previously, computational fluid dynamics simulations had been performed on pipe flows with different transverse and longitudinal ribbed geometries at Oak Ridge National Laboratory to compare the enhancement performance among different ribbed geometries. Rib shape morphing had been conducted to obtain an optimized rib profile. In the work presented here, the adjoint method is adopted in the ANSYS Fluent solver for turbulence model augmentation, and the Generalized k-ω (GEKO) turbulence model is employed because of its ability of tuning the turbulence model. The Nusselt number and pressure drop obtained from the channel flow with bottom ribbed wall experiments are used as the targets. Sensitivity analysis provides information as guidance to improve the turbulence model accuracy. The augmented GEKO model is tuned for the studied ribbed channel geometry and flow conditions, providing improved predictive accuracy within this context. Extension to other configurations offers potential but may require additional tuning and validation.

Xu, Tracy [ORNL] (ORCID:0009000193700887)↗
DOE OSTI2022

Investigation of core impurity transport in DIII-D diverted negative triangularity plasmas

Abstract Tokamak operation at negative triangularity has been shown to offer high energy confinement without the typical disadvantages of edge pedestals (Marinoni et al 2021 Nucl. Fusion 61 116010). In this paper, we examine impurity transport in DIII-D diverted negative triangularity experiments. Analysis of charge exchange recombination spectroscopy reveals flat or hollow carbon density profiles in the core, and impurity confinement times consistently shorter than energy confinement times. Bayesian inferences of impurity transport coefficients based on laser blow-off injections and forward modeling via the Aurora package (Sciortino et al 2021 Plasma Phys. Control. Fusion 63 112001) show core cross-field diffusion to be higher in L-mode than in H-mode. Impurity profile shapes remain flat or hollow in all cases. Inferred radial profiles of diffusion and convection are compared to neoclassical, quasilinear gyrofluid, and nonlinear gyrokinetic simulations. Heat transport is observed to be better captured by reduced turbulence models with respect to particle transport. State-of-the-art gyrokinetic modeling compares favorably with measurements across multiple transport channels. Overall, these results suggest that diverted negative triangularity discharges may offer a path to a highly-radiative L-mode scenario with high core performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Eddy Covariance Theory: A Review

Eddy covariance (EC), the gold standard for measuring ecosystem scale gas and heat exchanges, has transformed our understanding of the breathing of the biosphere, and thus global change biology. Despite numerous methodological improvements and insights gained from the technique, the community faces persistent challenges that have been present since the first EC measurements. Here, we review the theoretical developments underpinning EC. We present theoretical developments in four important areas that have relevance to EC measurements of the net ecosystem exchanges (NEE) of gases and heat from a single tower: (i) measuring the total vertical flux density, (ii) flux attenuation, (iii) coordinate rotations, and (iv) energy balance closure. Persistent problems with EC measurements, such as the inability to close the energy budget, led us to identify two priorities for revisiting the theory underlying: (i) sensible heat flux calculations, and (ii) constraining the mean vertical wind velocity. We present a framework for improved calculation of sensible heat flux derived from first principles of fluid mechanics and thermodynamics that considers coupled heat and mass transfer so that conservation of both is obeyed. These refinements are motivated by the need for unbiased measurements of energy and mass transfer between the land surface and atmosphere for ecosystem research and to validate satellite observations and land surface models.

ecosystem fluxes↗
DOE OSTI2023

Representation of atmosphere-induced heterogeneity in land–atmosphere interactions in E3SM–MMFv2

In the Energy Exascale Earth System Model (E3SM) Multi-scale Modeling Framework (MMF), where parameterizations of convection and turbulence are replaced by a 2-D cloud-resolving model (CRM), there are multiple options to represent land–atmosphere interactions. Here, we propose three different coupling strategies, namely the (1) coupling of a single land surface model to the global grid (MMF), (2) coupling a single land copy directly to the embedded CRM (SFLX2CRM), and (3) coupling a single copy of land model to each column of the CRM grid (MAML). In the MAML (Multi-Atmosphere Multi-Land) framework, a land model is coupled to CRM at the CRM-grid scale by coupling an individual copy of a land model to each CRM grid. Therefore, we can represent intra-CRM heterogeneity in the land–atmosphere interaction processes. There are 5-year global simulations run using these three coupling strategies, and we find some regional differences but overall small changes with respect to whether a land model is coupled to CRM or a global atmosphere. In MAML, the spatial heterogeneity within CRM induces stronger turbulence, which leads to the changes in soil moisture, surface heat fluxes, and precipitation. However, the differences in the MAML from the other two cases are rather weak, suggesting that the impact of using MAML does not justify the increase in cost.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Effects of spatiotemporal plasma power distribution on the modeling of ignition kernel evolution in quiescent and turbulent methane/air mixtures

Abstract The present work improves a phenomenological plasma-assisted combustion model by integrating the spatiotemporal distribution of plasma power density, thereby considering the evolution of plasma streamers in the modeling, and subsequently, better predicting the ignition kernel evolution. The improved phenomenological model is validated against experiments representing the plasma discharge and post-discharge ignition kernel evolution. Specifically, the new model demonstrates a more accurate prediction of ultrafast gas heating and O 2 dissociation during the plasma discharge, compared to the original model. In addition, the new model is found to closely match the experimental pressure wave and heated channel profiles post-discharge without the need for tuning the energy deposition (unlike the original model), highlighting its accuracy of post-discharge ignition kernel dynamics. The improved phenomenological model is then employed to investigate ignition kernel evolution for a stoichiometric methane-air discharge across various discharge gap configurations. Simulations reveal a non-uniform temperature and streamer distribution progressing from the electrode tips toward the center, contrasting uniform cylindrical discharges previously described in the original model. Streamer propagation is observed to be faster for larger gaps when maintained at the same average electric field for different discharge gaps. The tendency of smaller gaps to produce detached toroidal ignition kernels is observed, while larger gaps promote cylindrical and attached ignition kernels. Interactions between successive ignition kernels from consecutive discharges varied significantly, with the smallest gap (1 mm) promoting the quenching of the preceding ignition kernel due to the initial kernel–kernel separation. The intermediate gap (2 mm) promotes detached kernel growth. In contrast, in the largest gap (4 mm), kernels consistently combine and expand attached to electrodes. The impact of homogeneous isotropic turbulence is also explored, showing the persistence of ignition kernels early on but eventually quenching due to enhanced radical and heat losses with pronounced turbulence intensity.

Johnson, Praise Noah↗
DOE OSTI2025

Enhancing Turbulent Mixing and Microphysical Uniformity in a Tall Convection‐Cloud Chamber Through Idealized Heterogeneity of Boundaries

A large convection cloud chamber has been proposed for exploring aerosol–cloud–drizzle interactions under well‐controlled turbulent conditions. Recent theoretical and numerical studies suggest that a convection cloud chamber with two heated and two cooled sidewalls can significantly enhance the liquid water content and thus benefit drizzle initiation. However, a chamber with such a sidewall configuration develops stable stratification and extremely weak turbulence therein. In this study, we conduct large‐eddy simulations of a tall convection chamber with five different sidewall configurations consisting of alternating warm and cold patches. For each configuration, the total surface area of warm patches equals that of cold patches, resulting in the same expected cloud‐free supersaturation based on a flux budget model. Results show that changing the sidewall configuration, while keeping all other factors constant, can substantially enhance turbulent mixing and improve the uniformity of thermodynamic and cloud microphysical properties in the bulk region of the chamber. In addition, turbulence strength is positively correlated with liquid water content and negatively correlated with cloud droplet number concentration, consistent with theoretical predictions. Our results highlight the advantage of building a large cloud chamber using modular patches with individually controllable temperature and humidity to achieve well‐mixed conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Alfvén eigenmode-driven zonal modes saturate and heat thermal ions by cross-scale interactions

In scenarios where a sustained energetic particle source strongly drives toroidal Alfvén eigenmodes (TAE), and phase-space transport is insufficient to saturate TAE, this novel theory of TAE-zonal mode (ZM)-turbulence—self-regulated by cross-scale interactions (including collisionless ZF damping) – merits consideration. Zonal modes are driven by Reynolds and Maxwell stresses, without the onset of modulational instability. TAE evolution in the presence of ZMs conserves energy and closes the system feedback loop. The saturated zonal shears can be sufficient to suppress ambient drift-ion temperature gradient (ITG) turbulence, achieving an enhanced core confinement regime. The saturated state is regulated by linear and turbulent zonal flow drag. This regulation leads to bursty TAE spectral oscillations, which overshoot while approaching saturation. Heating by both collisional and collisionless ZM damping deposits alpha particle energy into the thermal plasma, achieving effective alpha channeling. This theory offers a mechanism for EP-induced transport barrier formation, and predicts a novel thermal ion heating mechanism.

ITB↗