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 541 records · Page 30

Validation of edge turbulence codes against the TCV-X21 diverted L-mode reference case

Self-consistent full-size turbulent-transport simulations of the divertor and scrape-off-layer (SOL) of existing tokamaks have recently become feasible. This enables the direct comparison of turbulence simulations against experimental measurements. In this work, we perform a series of diverted ohmic L-mode discharges on the tokamak à configuration variable (TCV) tokamak, building a first-of-a-kind dataset for the validation of edge turbulence models. This dataset, referred to as TCV-X21, contains measurements from five diagnostic systems from the outboard midplane (OMP) to the divertor targets—giving a total of 45 one- and two-dimensional comparison observables in two toroidal magnetic field directions. The experimental dataset is used to validate three flux-driven 3D fluid-turbulence models—GBS, GRILLIX and TOKAM3X. With each model, we perform simulations of the TCV-X21 scenario, individually tuning the particle and power source rates to achieve a reasonable match of the upstream separatrix value of density and electron temperature. We find that the simulations match the experimental profiles for most observables at the OMP—both in terms of profile shape and absolute magnitude—while a comparatively poorer agreement is found towards the divertor targets. The match between simulation and experiment is seen to be sensitive to the value of the resistivity, the heat conductivities, the power injection rate and the choice of sheath boundary conditions. Additionally, despite targeting a sheath-limited regime, the discrepancy between simulations and experiment also suggests that the neutral dynamics should be included. The results of this validation show that turbulence models are able to perform simulations of existing devices and achieve reasonable agreement with experimental measurements. Where disagreement is found, the validation helps to identify how the models can be improved. By publicly releasing the experimental dataset and validation analysis, this work should help to guide and accelerate the development of predictive turbulence simulations of the edge and SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

Hybrid combustion modeling approach for turbulent jet ignition in natural-gas pre-chamber spark-ignition engines at high EGR

Here, this study presented a hybrid modeling approach for simulating turbulent jet ignition and combustion processes in a natural-gas pre-chamber spark-ignition engine operating under exhaust gas recirculation (EGR) diluted conditions. In-depth analyses of experimental data and simulation results from previous work [Chinnathambi et al., ICEF2021-67836; Kim et al., Fuel 409: 137815, 2026] revealed two key findings: (i) the magnitude of pressure difference between the pre-chamber and main chamber ($∆P_{PC-MC}$) was positively correlated with the combustion duration from the moment of $∆P_{PC-MC}=0$ to the point of 5% mass fraction burned, with larger $∆P_{PC-MC}$ associated with longer duration; and (ii) the turbulent combustion regime in the main chamber transitioned from the broken reaction zone to the corrugated flamelet regime, with the Karlovitz number exceeding 100 immediately after turbulent hot jets were ejected from nozzles, coinciding with observed local extinction events. To accurately simulate the entire combustion process, a hybrid approach was developed under Reynolds-Averaged Navier Stokes framework, combining the G-equation model for pre-chamber combustion with the multi-zone well-stirred reactor approach and a turbulence-chemistry interaction (TCI) submodel for main chamber combustion. The TCI submodel accounted for the attenuation of reaction rates due to turbulent strain and modeled local extinction by suppressing reaction rates under certain flow and flame conditions. When applied to three EGR rate conditions toward the dilution limit, the hybrid modeling approach accurately reproduced experimental data in terms of cylinder pressure, apparent heat release rate, and the observed positive correlation, including the delayed onset of main chamber combustion—a feature not captured by existing combustion models.

computational fluid dynamics simulation↗

CFD analysis of a solid pin-fueled small modular fluoride salt-cooled reactor

A computational fluid dynamics (CFD) study is performed for the pre-conceptual solid pin-fueled small modular fluoride salt-cooled reactor developed by the Oak Ridge National Laboratory (ORNL). The solid and fluid regions in a 1/12th section of a fuel assembly are modeled. The solid region includes fuel and non-fuel pins, with heat generation in the active region of the fuel pins. Two different power profiles, uniform and center-peaked, are considered. The FLiBe coolant flows from the bottom to the top of the core, parallel to the bank of fuel and non-fuel pins. The k-ω shear-stress transport (SST) model is chosen as the baseline turbulence model. The effects of grid refinement, inlet turbulence specification, and turbulence models on the temperature and pressure drop predictions are studied. Turbulence model sensitivity is investigated by comparing the results from the k-ω SST model with other two-equation models (k-ω baseline or BSL and k-ε realizable) as well as with anisotropic Reynolds stress transport models (linear pressure-strain and stress-BSL). Here the results show that the choice of turbulence model has a significant impact on the pin temperature and bundle pressure drop predictions. The results from the baseline turbulence model show good agreement between the fuel pin temperatures and bundle pressure drop values predicted using a subchannel model previously developed by the authors of the present study.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

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

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental investigation of turbulent wake flows in a helically wrapped rod bundle in presence of localized blockages

In nuclear sodium fast reactors, bundles of rods are tightly packed into a triangular lattice, enclosed in a hexagonal duct, and each pin is spirally wrapped with a thin wire. Flow blockages can potentially impact the local flow characteristics and heat transfer mechanisms in the bundle due to its small subchannel size. The effects of the blockage on the flow structures and heat transfer mechanisms are important aspects that require an accurate investigation. In this study, the flow-field characteristics in the vicinity of a blockage located in the exterior subchannel of rod bundles with helically wrapped wires were experimentally investigated. The velocity fields in the exterior subchannel were acquired by applying matched-index-of-refraction and time-resolved particle image velocimetry (TR-PIV) techniques for Reynolds numbers of Re 1 = 4000 and Re 2 = 17 000, i.e., equivalent to Rew 1 = 19 600 and Rew 2 = 83 200, respectively, based on the blockage width. The results from the TR-PIV measurements revealed an arch-shaped vortex with a large flow recirculation and a pair of counter-rotating vortices in the wake region downstream of the blockage, which is commonly observed in the wake flow of bluff bodies. The relative lateral distance and angle between the two vortices decreased when the Reynolds numbers increased. Profiles of maximum turbulence intensity along the shear layers illustrated the transition process including the growth, peak, and decay along the flow direction. From the spectral analysis of the turbulent velocities extracted at points along the shear layer, the Strouhal numbers ( St ) representing the vortex shedding frequency were found to be St = 0.25 and St = 0.56 for the left and right shear layers, respectively. Characteristics of shear layers generated by the blockage in the exterior subchannel were investigated via the two-point cross correlation of fluctuating velocities. The patiotemporal cross correlations of turbulent velocities, computed at points in the region where the left shear layer exhibited rolling effects and vortex breakdowns, were considerably wider and longer. The convection velocity U c was estimated to be ~0.82 U m to 0.93 U m . Proper orthogonal decomposition (POD) analysis was applied to the instantaneous velocity fields to extract the statistically dominant flow structures. It was found that POD modes 2–3 and 4–5 formed the pair modes when the corresponding POD temporal coefficients depicted sinusoidal shapes and exhibited nearly circular orbits in the phase space. Spectral analysis of the POD temporal coefficients confirmed the vortex shedding frequencies detected in the analysis of turbulent velocities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Experimental observations of bifurcated power decay lengths in the near Scrape-Off Layer of ST40 High Field Spherical Tokamak

The scrape-off layer parallel heat flux decay lengths measured at ST40, a high field, low aspect ratio spherical tokamak, have been observed to bifurcate into two groups. The wide group follows established H-mode scalings (ranging between 2 to 8 mm) while the narrow group falls up to 10 times below these scalings (between 0.2 and 0.8 mm), being comparable to the ion total Larmor radius rather than the ion poloidal Larmor radius. The heat flux profiles of the latter group can only be described by a multi-exponential function, rather than the single exponential function convoluted with a Gaussian. The onset of the narrow scrape-off layer width is observed to be associated with suppressed magnetic fluctuations, suggesting reduced electromagnetic turbulence levels in the SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigating the Diurnal Radiative, Turbulent, and Biophysical Processes in the Amazonian Canopy‐Atmosphere Interface by Combining LES Simulations and Observations

Abstract We investigate the atmospheric diurnal variability inside and above the Amazonian rainforest for a representative day during the dry season. To this end, we combine high‐resolution large‐eddy simulations that are constrained and evaluated against a comprehensive observation set, including CO 2 concentrations, gathered during GoAmazon2014/15. We design systematic numerical experiments to quantify whether a multilayer approach in solving the explicit canopy improves our canopy‐atmosphere representation. We particularly focus on the relationship between photosynthesis and plant transpiration, and their distribution at leaf and canopy scales. We found the variability of photosynthesis drivers like vapor pressure deficit and leaf temperature to be about 3 times larger for sunlit leaves compared to shaded leaves. This leads to a large spread on leaf stomatal conductance values with minimum and maximum values varying more than 100%. Regarding the turbulent structure, we find wind‐driven stripe‐like shapes at the canopy top and structures resembling convective cells at the canopy. Wind‐related variables provide the best spatiotemporal agreement between model and observations. The potential temperature and heat flux profiles agree with an observed decoupling near the canopy top interface, although with less variability and cold biases of up to 3 K. The increasing complexity on the biophysical processes leads to the largest disagreements for evaporation, CO 2 plant assimilation and soil efflux. The model is able to capture the correct dependences and trends with the magnitudes still differing. We finally discuss the need to revise leaf and soil models and to complete the observations at leaf and canopy levels.

54 ENVIRONMENTAL SCIENCES↗

Background-Oriented Schlieren Velocimetry of Helium Coolant Flow in Additively Manufactured Channels

High-pressure helium gas cooling is an attractive solution for thermal management of the fusion blanket first wall, as this coolant is chemically and neutronically inert and separable from hydrogenic species. However, due to the low thermal mass of helium, geometric optimization of these channels is required to provide sufficient cooling at manageable flow rates and pumping burdens. Increasingly, analysis and optimization of these coolant channels rely on computational fluid dynamics (CFD) simulations, and these require relevant experimental data for turbulence model validation. Toward this end, a high-pressure helium gas flow visualization system has been employed to image the flow of helium in flow channels with one-sided heating, mimicking the blanket first wall environment. Flow of helium at 4 MPa pressure and flow rates up to 68 g/s (Reynolds number 57 000) is supplied to rectangular channel test sections, with uniform heating applied to the bottom wall of the channel at heat fluxes varied between roughly 50 and 130 kW/m2. A high-speed camera is used to image index of refraction (IOR) gradients in the fluid via background-oriented schlieren (BOS), and temperature and pressure instrumentation are used to characterize thermal-hydraulic performance of each channel. Cross correlation of time-resolved BOS images is then used to calculate time-averaged 2-D helium velocity fields. Flow in additively manufactured (AM) channels is examined in this manner, including both featureless channels and those containing baffling as a heat transfer enhancement. The flow distribution seen in the featureless case differs significantly from that seen in prior simulations, whereas the flow in the baffled case shows the predicted behavior of flow forced along the heated wall. This augmented flow distribution is seen to increase the heat transfer coefficient in the baffled test section. Here, strategies are discussed for ongoing and future validation of these simulations, with the aim of model deployment for blanket cooling design and optimization.

Additive manufacturing↗

Multi-device Studies of Pedestal Physics and Confinement in the I-mode Regime

This paper describes joint ITPA studies of the I-mode regime, which features an edge thermal barrier but with L-mode-like particle and impurity transport, and without ELMs. The regime has been demonstrated on Alcator C-Mod, ASDEX Upgrade and DIII-D tokamaks, over a wide range of device parameters and pedestal conditions, including low collisionality. Access to I-mode is independent of heating method (neutral beam injection, ICRH and/or ECRH). Normalized energy confinement H98y,2 ≥ 1 has been achieved for a range of 3≤q95≤4.9 and scales favourably with power. Changes in turbulence in the pedestal region accompany the transition from L-mode to I-mode. The L-I threshold increases with plasma density and current, and with device size, but weakly with BT. The upper limit of power for I-modes increases with toroidal field, and is largest on Alcator C-Mod at B > 5 T. Issues for extrapolation to ITER and other future fusion devices are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantitative comparisons of electron-scale turbulence measurements in NSTX via synthetic diagnostics for high-k scattering.

Two synthetic diagnostics are implemented for the high-k scattering system in NSTX [D. R. Smith, E. Mazzucato, W. Lee, H. K. Park, C. W. Domier, and N. C. Luhmann Jr., Rev. Sci. Instrum., 79, 123501 (2008)] allowing direct comparisons between synthetic and experimentally detected frequency and wavenumber spectra of electron-scale turbulence fluctuations. Synthetic diagnostics were formulated in real-space and in wavenumber space, and were deployed to realistic electron-scale simulations carried out with the GYRO code [J. Candy and R. E. Waltz, J. Comput. Phys. 186 (2003) 545]. A highly unstable ETG regime (electron temperature gradient mode) in a modest beta NSTX NBI-heated H-mode discharge was chosen for the analysis. Mapping the measured wavenumbers to field aligned coordinates showed that the high-k system is sensitive to fluctuations that are closer to the spectral peak in the density fluctuation wavenumber spectrum (streamers) than originally predicted. Analysis of synthetic spectra show that the frequency response of the detected fluctuations is dominated by Doppler shift and insensitive to the turbulence drive. The shape of the high-k density fluctuation wavenumber spectrum is sensitive to the ETG turbulence drive conditions, and can be reproduced in a sensitivity scan of the most pertinent turbulent drive terms in the simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

EAST steady-state long pulse H-mode with core-edge integration for CFETR

Recent EAST experiment has successfully demonstrated long pulse steady-state high plasma performance scenario and core-edge integration since the last IAEA in 2018. A discharge with a duration over 60s with β P ~2.0, β N ~1.6, H 98y2 ~1.3 and internal transport barrier on electron temperature channel is obtained with multi-RF power heating and current drive. A higher β N (β N ~1.8, β p ~2.0, H 98y2 ~1.3, n e /n GW ~0.75) with a duration of 20s is achieved by using the modulated neutral beam and multi-RF power, where several normalized parameters are close or even higher than the phase III 1GW scenario of CFETR steady-state. High-Z impurity accumulation in the plasma core is well controlled in a low level by using the on-axis ECH. Modelling shows that the strong diffusion of TEM turbulence in the central region prevents tungsten impurity to accumulate. More recently, EAST has demonstrated compatible core-edge integration discharges in the high βp scenario: high confinement H 98y2 >1.2 with high β P ~2.5/β N ~2.0 and f bs ~50% is sustained with reduced divertor heat flux at high density n e /n GW ~0.7 and moderate q 95 ~6.7. By combining active impurity seeding through radiative divertor feedback control and strike point splitting induced by resonant perturbation coil, the peak heat flux is reduced by 20-30% on the ITER-like tungsten divertor, here a mixture of 50% neon and 50% D2 is applied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Calculation of the Heat Transfer Coefficient in the Outer Body for a Rotational Detonation

Unsteady heat transfer characterization on the combustion surfaces of Rotational Detonation Engines (RDE) is not well understood. It is generally thought that the complex nature of the unsteady, reacting, compressible fluid flow inside the combustion anulus of the RDE causes the convective heat transfer coefficient to be significantly higher than it is in other applications. Empirical models that have been used to analyze this strictly apply to steady flows where dimensionless groups can be employed. The reacting flow fields that are characteristic of RDEs are inherently compressible, three dimensional, unsteady, and turbulent, having properties that change by orders of magnitude throughout the flow field. They will therefore contain multiple length scales and time scales operating everywhere in the flow during all times. Because of this it is not likely that the RDE flow fields will lend themselves to explanation using simple dimensionless parameters. The dimensionless groups have meaning only in situations where length scales and time scales are singular and well defined. In spite of this it may be possible to get a relatively good idea of what the convection heat transfer coefficient is. In this work a numerical study is performed where the inside wall surface temperature distribution in the RDE outer body is systematically changed over a given range that would be characteristic of the start-up flows inside an RDE. For each case, temperature distributions inside the outer containment wall of the RDE was calculated and compared with experimental data. The closet match can then be used to directly calculate the convection heat transfer coefficient on the inside surface of the RDE.

VanOsdol, John↗

Development of a Benchmark Eddy Flux Evapotranspiration Dataset for Evaluation of Satellite-Driven Evapotranspiration Models Over the CONUS

A large sample of ground-based evapotranspiration (ET) measurements made in the United States, primarily from eddy covariance systems, were post-processed to produce a benchmark ET dataset. The dataset was produced primarily to support the intercomparison and evaluation of the OpenET satellite-based remote sensing ET (RSET) models and could also be used to evaluate ET data from other models and approaches. OpenET is a web-based service that makes field-delineated and pixel-level ET estimates from well-established RSET models readily available to water managers, agricultural producers, and the public. The benchmark dataset is composed of flux and meteorological data from a variety of providers covering native vegetation and agricultural settings. Flux footprint predictions were developed for each station and included static flux footprints developed based on average wind direction and speed, as well as dynamic hourly footprints that were generated with a physically based model of upwind source area. The two footprint prediction methods were rigorously compared to evaluate their relative spatial coverage. Data from all sources were post-processed in a consistent and reproducible manner including data handling, gap-filling, temporal aggregation, and energy balance closure correction. The resulting dataset included 243,048 daily and 5,284 monthly ET values from 194 stations, with all data falling between 1995 and 2021. We assessed average daily energy imbalance using 172 flux sites with a total of 193,021 days of data, finding that overall turbulent fluxes were understated by about 12% on average relative to available energy. Multiple linear regression analyses indicated that daily average latent energy flux may be typically understated slightly more than sensible heat flux. This dataset was developed to provide a consistent reference to support evaluation of RSET data being developed for a wide range of applications related to water accounting and water resources management at field to watershed scales.

54 ENVIRONMENTAL SCIENCES↗

Cloud-atmosphere impacts on the central Arctic surface energy budget (Final Report)

This project investigated how clouds and atmosphere-sea ice interactions shape the central Arctic surface energy budget using observations from the MOSAiC expedition and coordinated model experiments. Analyses focused on three themes. First, a new method applied to thermistor string data revealed that snow thermal conductivity is systematically higher than assumed in most models, strongly influencing conductive heat flux and sea ice growth, especially under clear skies. Second, a full annual cycle of surface energy budget observations, combined with detailed cloud microphysical retrievals, demonstrated how cloud regimes regulate radiative, turbulent, and conductive fluxes, driving seasonal contrasts in sea ice energy balance: in winter, radiative forcing elicits strong turbulent and conductive responses, while in summer excess energy is primarily partitioned into surface melt. Third, these observations were used to evaluate and improve regional and global forecast models. While some models captured key Arctic cloud-radiation states, most exhibited persistent deficiencies in representing liquid-containing clouds, snow-on-sea-ice processes, and flux responses to radiative forcing. Together, these results provide unprecedented benchmarks for understanding coupled Arctic processes, improve parameterizations of snow, sea ice, and cloud interactions, and highlight critical pathways by which atmospheric variability drives sea ice change in a rapidly evolving Arctic system. Lastly, the project was a resounding success, resulting in 38 peer-reviewed publications and dozens of presentations, while also supporting an early career scientist.

54 ENVIRONMENTAL SCIENCES↗

Describing the Mechanism of Instability Suppression Using a Central Pilot Flame With Coupled Experiments and Simulations

Pilot flames are commonly used to extend combustor operability limits and suppress combustion oscillations in low-emissions gas turbines. Combustion oscillations, a coupling between heat release rate oscillations and combustor acoustics, can arise at the operability limits of low-emissions combustors where the flame is more susceptible to perturbations. While the use of pilot flames is common in land-based gas turbine combustors, the mechanism by which they suppress instability is still unclear. In this study, we consider the impact of a central jet pilot on the stability of a swirl-stabilized flame in a variable-length, single-nozzle combustor. Previously, the pilot flame was found to suppress the instability for a range of equivalence ratios and combustor lengths. We hypothesize that combustion oscillation suppression by the pilot occurs because the pilot provides hot gases to the vortex breakdown region of the flow that recirculate and improve the static, and hence dynamic, stability of the main flame. This hypothesis is based on a series of experimental results that show that pilot efficacy is a strong function of pilot equivalence ratio but not pilot flow rate, which would indicate that the temperature of the pilot products as well as the combustion intensity of the pilot flame play more of a role in oscillation stabilization than the length of the pilot flame relative to the main flame. Further, the pilot-flame efficacy increases with pilot-flame equivalence ratio until it matches the main-flame equivalence ratio; at pilot equivalence ratios greater than the main equivalence ratio, the pilot-flame efficacy does not change significantly with pilot equivalence ratio. To understand these results, we use large-eddy simulation (LES) to provide a detailed analysis of the flow in the region of the pilot flame and the transport of radical species in the region between the main flame and pilot flame. Furthermore, the simulation, using a flamelet/progress variable-based chemistry tabulation approach and standard eddy viscosity/diffusivity turbulence closure models, provides detailed information that is inaccessible through experimental measurements.

42 ENGINEERING↗