SEARCH · Search NASA

Results for “Gradient flow”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2025

A subgrid model for electron-scale turbulent transport in global ion-scale gyrokinetic simulations of tokamak plasmas

A subgrid electron–temperature–gradient (ETG) model is demonstrated here which averages local electron-scale turbulence over intermediate scales in space and time to include in global ion–temperature–gradient (ITG) simulations. This approach results in ion-scale equations, which incorporate the electron heat transport from ETG turbulence and the effects of electron-scale turbulence on the ion scale. Flux-tube ETG Cyclone Base Case simulations are performed at different radial locations and a kinetic form of the flux is added to global ion-scale (IS) simulations as a source term. Analytic radial profiles of ETG heat flux are constructed and compared to flux-tube results at multiple radial locations. Different ratios of ITG to ETG heat flux levels are considered and the results of capturing ETG heat transport in global ion-scale simulations are discussed. The potential coupling of the ETG streamer potential and intermediate-scale zonal flows to the ion scale is further addressed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Divide and conquer: Learning chaotic dynamical systems with multistep penalty neural ordinary differential equations

Forecasting high-dimensional dynamical systems is a fundamental challenge in various fields, such as geosciences and engineering. Neural Ordinary Differential Equations (NODEs), which combine the power of neural networks and numerical solvers, have emerged as a promising algorithm for forecasting complex nonlinear dynamical systems. However, classical techniques used for NODE training are ineffective for learning chaotic dynamical systems. In this work, we propose a novel NODE-training approach that allows for robust learning of chaotic dynamical systems. Here, our method addresses the challenges of non-convexity and exploding gradients associated with underlying chaotic dynamics. Training data trajectories from such systems are split into multiple, non-overlapping time windows. In addition to the deviation from the training data, the optimization loss term further penalizes the discontinuities of the predicted trajectory between the time windows. The window size is selected based on the fastest Lyapunov time scale of the system. Multi-step penalty(MP) method is first demonstrated on Lorenz equation, to illustrate how it improves the loss landscape and thereby accelerates the optimization convergence. MP method can optimize chaotic systems in a manner similar to least-squares shadowing with significantly lower computational costs. Our proposed algorithm, denoted the Multistep Penalty NODE, is applied to chaotic systems such as the Kuramoto-Sivashinsky equation, the two-dimensional Kolmogorov flow, and ERA5 reanalysis data for the atmosphere. It is observed that MP-NODE provide viable performance for such chaotic systems, not only for short-term trajectory predictions but also for invariant statistics that are hallmarks of the chaotic nature of these dynamics.

Chaotic dynamical systems↗

DOE OSTI2024

Extreme elevational migration spurred cryptic speciation in giant hummingbirds

The ecoevolutionary drivers of species niche expansion or contraction are critical for biodiversity but challenging to infer. Niche expansion may be promoted by local adaptation or constrained by physiological performance trade-offs. For birds, evolutionary shifts in migratory behavior permit the broadening of the climatic niche by expansion into varied, seasonal environments. Broader niches can be short-lived if diversifying selection and geography promote speciation and niche subdivision across climatic gradients. To illuminate niche breadth dynamics, we can ask how “outlier” species defy constraints. Of the 363 hummingbird species, the giant hummingbird (Patagona gigas) has the broadest climatic niche by a large margin. To test the roles of migratory behavior, performance trade-offs, and genetic structure in maintaining its exceptional niche breadth, we studied its movements, respiratory traits, and population genomics. Satellite and light-level geolocator tracks revealed an >8,300-km loop migration over the Central Andean Plateau. This migration included a 3-wk, ~4,100-m ascent punctuated by upward bursts and pauses, resembling the acclimatization routines of human mountain climbers, and accompanied by surging blood-hemoglobin concentrations. Extreme migration was accompanied by deep genomic divergence from high-elevation resident populations, with decisive postzygotic barriers to gene flow. The two forms occur side-by-side but differ almost imperceptibly in size, plumage, and respiratory traits. The high-elevation resident taxon is the world’s largest hummingbird, a previously undiscovered species that we describe and name here. The giant hummingbirds demonstrate evolutionary limits on niche breadth: when the ancestral niche expanded due to evolution (or loss) of an extreme migratory behavior, speciation followed.

Science & Technology - Other Topics↗

DOE OSTI2025

Development of an In Situ Fission Gas Release Instrument for Fuel Sample Irradiations in the High Flux Isotope Reactor

Experimental measurement of gaseous fission product release with respect to temperature and burnup is a critical aspect of understanding nuclear fuel performance, validating predictive models, and qualifying new fuels. To measure this phenomenon in real-time, Oak Ridge National Laboratory has developed an instrument for measuring in situ fission gas release from small-scale fuel samples irradiated in the High Flux Isotope Reactor (HFIR). The instrument uses a continuous flow of Heover the fuel samples to sweep gaseous fission products from a sealed capsule in the HFIR Be reflector to an instrument station adjacent to the reactor. The instrument station houses two high-purity germanium (HPGe) detectors that measure decay gamma rays from fission products passing through a room temperature dwell chamber placed over the detector crystal. The sealed capsules in the reactor are designed to modulate fuel sample temperatures between 700 and 1,100°C by changing the Ar/He gas mixture surrounding the capsules during irradiation. N-type thermocouples are incorporated into the capsule housing to record real-time fuel temperatures. The capsules are heated primarily by prompt gamma rays emitted from the HFIR core with minimal heat contributions from fission in the fuel samples to minimize temperature gradients in the specimens for separate-effects characterization of the material. This paper describes modeling of time-dependent nuclear heating and fission product formation in fuel samples, thermal characteristics of the in-core capsules, and expected gaseous fission product gamma spectra at the HPGe instrument station.

Mulligan, Padhraic L [ORNL] (ORCID:000000025826540↗

DOE OSTI2026

Organic colloid composition in variable-redox porewaters within a mountainous floodplain

Redox gradients, often driven by changes in sediment moisture levels in porous, heterogeneous groundwater systems, create dynamic conditions that may promote the production and transport of colloids within natural waters. While much research has focused on the inorganic composition of colloids, the organic composition remains less well understood. Organic matter (OM) in colloids may associate with minerals, complex metal ions, and serve as an electron donor for microbial respiration; therefore, its composition is of high interest. We examined the composition of porewater OM along a redox gradient in a riparian soil located along the Slate River in Crested Butte, Colorado, USA as a function of depth (90, 130, 200, and 350 cm below ground surface). All depths were oxic to suboxic, except 200 cm, where the products of iron and sulfate reduction were observed concomitant with an increase in dissolved and/or colloidal OM, pH, alkalinity, and conductivity. We investigated the composition of porewater using correlated scanning transmission X-ray microscopy and transmission electron microscopy. We observed a change in the OM chemistry from carboxylate-rich at the 200 cm depth (where it was also enmeshed with non-crystalline iron) to phenol- and substituted-aromatic-rich at other depths. Radiocarbon dating revealed carbon in porewater at 200 cm was younger than depths above and below. Soil porewater can flow down into the underlaying gravel bed during baseflow conditions, thus we speculate whether riparian porewater could transport iron- and carboxylate-rich organic colloids into surrounding surface waters through the gravel bed conduit.

Alpine floodplain↗

DOE OSTI2026

Data from Stewart et al. 2026 "Organic Colloid Composition in Variable-Redox Porewaters within a Mountainous Floodplain"

alpine floodplain↗

DOE OSTI2024

Identifying Controlling Variables for Mercury Vapors in Alpha-4 at Y-12: Two Year Data Collection Update

Multiple sensor packages were deployed at Alpha-4 by SRNL, in collaboration with United Cleanup Oak Ridge LLC (UCOR), to monitor mercury vapor concentrations and meteorological parameters. These sensors collected data, inside and outside of the legacy-use facility, for approximately two years. Though data gaps still exist, particularly in colder months, several controlling variables were identified that govern mercury vapor concentrations within Alpha-4. Temperature, barometric pressure gradients, humidity, and wind speed have been identified as controlling variables. A strong positive correlation was seen between mercury vapor concentrations and temperature which generally followed diurnal fluctuations. Temperatures below approximately 10 degrees Celsius did not show any spikes above the PEL, indicating more work can be performed at any time during the winter months – more data should be collected to confirm consistency in this finding. Additionally, spikes in mercury vapor concentrations above that of the permissible exposure limit (PEL; 100 µg/m 3 ) occurred primarily in late afternoon or evening/overnight hours (between 3 PM and 6 AM), which suggests D&D operations might be best scheduled during morning or daytime hours prior to the late afternoon. However, a limited number of spikes did occur outside of the identified window, although this may be attributed to disturbances in air flow and mercury vapor release from work activities performed inside of the Alpha-4 building. The analysis conducted allows for a strong predictive capability for estimating mercury vapor concentrations based upon accurate meteorological parameters. Still, additional data collection, particularly in the winter months, could help to strengthen the predictive power and validate the identified data trends. Further, increased temporal resolution could also help to better characterize the incipient stages of the increases and decreases in the mercury vapor concentration. Continued monitoring support by SRNL at Y-12 is underway at Alpha-4 to further close remaining data gaps and support deactivation and decommissioning work. Within a collaborative effort with UCOR, the SRNL team is collecting mercury vapor data to study the efficacy of a novel mercury suppressant, FerroBlack® which was recently deployed at Alpha-4. In addition, modifications to the current monitoring setup to increase measurement resolution is also being investigated.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Computationally inexpensive part-scale thermal history of additive friction-stir deposition

This study presents an analytical model for steady-state power generation and tool heat loss in additive friction-stir deposition (AFSD), developed to enable part-scale thermal simulation while remaining computationally inexpensive. The model predicts total generated power, yielding 3.7–4.7 kW across deposition temperature setpoints of 400–460 °C for the deposition of AA6061 with a Be-Cu tool. This corresponds to 90–95% of the reported spindle power. Tool heat loss is experimentally determined by calibrating a steady-state energy balance between the generated power, the substrate-deposition thermal gradient, and a temperature dependent tool heat loss term: q tool (T) = a + b (T - 400°C) with a = 2.7 x 10 6 Wm -2 and b = 9.5 x 10 3 Wm -2 K -1 . The calibration indicates that about 69% of the generated heat is conducted into the tool for this configuration, which is much higher than previously reported. The calibrated heat-source is implemented in finite element software (Adamantine) to simulate the transient thermal history of a 100 cm 3 representative build in 8 min on a standard desktop (at 0.635 mm build-height resolution). For the first three layers, the substrate temperatures between simulation and experiment are within 10% mean absolute percentage error. Sensitivity analysis indicates that uncertainties in average deposition temperature and deformation localization (stir-zone geometry, depth, and spatial dependance of strain-rate and flow stress) dominate model variance, motivating additional experimental verification.

Additive Friction-Stir Deposition↗

DOE OSTI2025

TEM turbulence in simulation and experiment in quasisymmetric and degraded-symmetry configurations of the HSX stellarator

The Helically Symmetric Experiment (HSX) has demonstrated reduced neoclassical transport and flow damping with quasisymmetry, but the difference of turbulent transport between quasisymmetric and degraded-symmetry configurations has not been definitively assessed. Profile analysis presented here shows that anomalous transport at the mid-radius is larger with degraded quasisymmetry when temperature and density profiles are matched. An initial comparison (${T_\mathrm{i}/T_\mathrm{e}} = 1$) between the quasisymmetric and degraded symmetry configurations finds that trapped-electron mode (TEM) turbulence can support larger heat fluxes in the degraded symmetry configuration despite smaller linear growth rates. However, further simulations at the experimental temperature ratio ($T_\mathrm{i}/T_\mathrm{e} = 0.2$) do not reproduce this effect. Experimental measurements of heat transport and density fluctuation amplitude are compared to nonlinear gyrokinetic simulations at the experimental temperature ratio. The heat flux and density fluctuation amplitude in simulation reproduce a substantial dependence on the density gradient, and the simulated heat flux matches measurements within experimental uncertainties, indicating that TEM turbulence is responsible for the anomalous transport at the steep density gradient in HSX.

HSX↗

DOE OSTI2024

Ion-temperature- and density-gradient-driven instabilities and turbulence in Wendelstein 7-X close to the stability threshold

Electrostatic gyrokinetic instabilities and turbulence in the Wendelstein 7-X stellarator are studied. Particular attention is paid to the ion-temperature-gradient (ITG) instability and its character close to marginal stability [Floquet-type turbulence (Zocco et al., Phys. Rev. E, vol. 106, 2022, p. L013202) with no electron temperature gradient]. The flux tube version of the δƒ code stella (Barnes et al., J. Comput. Phys., vol. 391, 2019, pp. 365–380) is used to run linear and nonlinear gyrokinetic simulations with kinetic electrons. The nature of the dominant instability depends on the wavelength perpendicular to the magnetic field, and the results are conveniently displayed in stability diagrams that take this dependence into account. This approach highlights the presence of universal instabilities, which are less unstable but have longer wavelengths than other modes. A quasi-linear estimate of the heat flux suggests they are relevant for transport. Close to the stability threshold, the linear eigenmodes and turbulence form highly extended structures along the computational domain if the magnetic shear is small. Numerical experiments and diagnostics are undertaken to assess the resulting radial localisation of the turbulence, which affects the interaction of the latter with zonal flows. Increasing the amplitude of the magnetic shear (e.g.through current drive) has a stabilising effect on the turbulence and, thus, reduces the nonlinear energy transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Global gyro-kinetic ion temperature gradient and trapped electron mode turbulence modelling in X -point geometry in negative and positive triangularity

Comparative modelling of ion temperature gradient (ITG)/trapped electron mode (TEM) turbulence in negative (NT) and positive (PT) triangularity plasma shapes was done using the nonlinear global full-f gyrokinetic particle code JOREK-GK in the realistic X-point tokamak geometry including the Scrape Off Layer (SOL) for TCV and DIII-D parameters. A comparison of JOREK-GK code with the gyrokinetic codes GS2 and GENE was done using NT/PT triangularity TCV L-modes parameters showing good agreement between codes in linear growth rates and clear beneficial effect of NT as compared to PT. Global non-linear modelling of the ITG/TEM saturated turbulence for realistic DIII-D NT pulses was done and compared with numerically constructed PT equilibrium with the same plasma profiles. Existence of longer correlation length of density fluctuations in PT compared to NT was demonstrated. Stronger and more sheared zonal flows are generated via Reynolds stress in NT compared to PT. These factors are stabilizing for TEM/ITG turbulence in NT and lead to smaller heat fluxes and heat conductivities in NT compared to PT configuration. Weak dependence of plasma confinement on collisionality and plasma rotation was found in modelling of DIII-D NT shots similar to the experiment. The confinement scaling with normalized ion gyro-radius ρ* was estimated both for NT and PT. Bohm-like scaling was obtained in both configurations, however with better confinement for NT compared to PT which could be favourable factor for reactor size machines.

ITG↗

DOE OSTI2024

Highly Ordered Eutectic Mesostructures via Template‐Directed Solidification within Thermally Engineered Templates

Template-directed self-assembly of solidifying eutectics results in emergence of unique microstructures due to diffusion constraints and thermal gradients imposed by the template. Here, the importance of selecting the template material based on its conductivity to control heat transfer between the template and the solidifying eutectic, and thus the thermal gradients near the solidification front, is demonstrated. Simulations elucidate the relationship between the thermal properties of the eutectic and template and the resultant microstructure. The overarching finding is that templates with low thermal conductivities are generally advantageous for forming highly organized microstructures. When electrochemically porosified silicon pillars (thermal conductivity < 0.3 Wm −1 K −1 ) are used as the template into which an AgCl-KCl eutectic is solidified, 99% of the unit cells in the solidified structure exhibit the same pattern. In contrast, when higher thermal conductivity crystalline silicon pillars (≈100 Wm −1 K −1 ) are utilized, the expected pattern is only present in 50% of the unit cells. The thermally engineered template results in mesostructures with tunable optical properties and reflectances nearly identical to the simulated reflectances of perfect structures, indicating highly ordered patterns are formed over large areas. This work highlights the importance of controlling heat flows in template-directed self-assembly of eutectics.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Evaluation of thermal energy storage for central receiver concentrating solar power plants under charging cycles through structural assessment

Electricity production by concentrated solar power (CSP) systems has stood out among energy transition alternatives due to their large generation capacity (between 30 MW and 400 MW) and high capacity factor (80%) when incorporating thermal storage systems (TES). Although TES tanks in operation worldwide are made of austenitic steels with remarkable mechanical and anticorrosive properties, operational failures have been reported due to low-cycle fatigue, creep, and the abrupt release of residual stresses generated during tank manufacturing. A critical operation for salt tanks is the initial daily charging operation, given the low fluid level and the thermal gradients between the inventory and the salt entering the tank. This study presents a structural simulation of a 39.6 m-diameter molten salt tank for central receiver CSP plants, accounting for temperature and pressure variations during charging and different sparger ring inlet configurations. The initial pre-stressed condition of the floor, resulting from the manufacture of welded plates, is used as a boundary condition. The fatigue and creep life analyses were performed in accordance with the ASME BPVC and API standards. It was found that the tank's floor is in a less favourable mechanical condition than the wall, due to its initial manufacturing condition and the maximum temperature differences during charging. A critical zone was identified in an area affected by residual stress on the floor and under the sparger ring. The baseline sparger ring configuration, with 52 2-inch-diameter orifices and a flow direction of 90 degrees , results in an average floor stress of approximately 25 MPa and creep damage after 1.2 years of tank operation.

14 SOLAR ENERGY↗

DOE OSTI2024

Plasma potential shaping using end-electrodes in the Large Plasma Device

We perform experiments in the Large Plasma Device at the University of California, Los Angeles, studying how different end-electrode biasing schemes modify the radial potential profile in the machine. We impose biasing profiles of different polarities and gradient signs on a set of five concentric electrodes placed 12 m downstream from the plasma source. We find that imposing concave-down profiles (negative potential radial gradient) on the electrodes creates radial potential profiles halfway up the plasma column that are comparable to those imposed on the electrodes and a few electron temperature in height, regardless of the biasing polarity. On the other hand, imposing concave-up profiles (positive potential radial gradient) leads to non-monotonic radial potential profiles. This observation can be explained by the current drawn through the electrodes and the parallel plasma resistivity, highlighting their important role in controlling the rotation of plasma. Concave-down plasma potential profiles, obtained by drawing electrons on the axis, are predicted to drive azimuthal drift velocities that can approach significant fractions of the ion sound speed in the central region of the plasma column.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Continued Testing and Modeling of the AirSquared All-Metal Tritium Scroll Pump

The implementation of fusion energy requires the processing of the deuterium-tritium mixture used to fuel the reaction. Many gas processing operations require either a gas-motive force or a pressure gradient across a system; both are achieved using vacuum pumps. However, because of the material challenges of working with tritium, specialized pumps are required that use only metal on wetted surfaces. This prevents tritium retention and degradation of materials such as oils and polymers and reduces the generation of tritiated waste. Typical gas handling systems for tritium operations involve metal bellows pumps, all-metal scroll pumps, and metal turbomolecular pumps, depending on the pressure requirements. With the industry standard Normetex Type 15 m 3 /h (typically referred to as the Normetex 15) pump no longer available, the fusion market requires viable alternatives. In this study, we test an American-made pump, the AirSquared V16H034A-C01, and compare its performance to the Normetex. A performance model of the AirSquared pump is also produced, which can be applied to general process models to generate an accurate description of pressures and flows through various unit operations.

08 HYDROGEN↗

DOE OSTI2024

Analysis of thermal and mechanical properties with inventory level of the molten salt storage tank in central receiver concentrating solar power plants

Molten salt thermal energy storage (TES) tanks ensure steady power output of concentrating solar power (CSP) plants; however, recent tank failures have highlighted the need for further analysis. Current studies primarily focus on analyzing the molten salt flow, heat transfer, and thermal efficiency. Additionally, research on the latest tank structures is limited and lacks newest experimental validation. This study measures temperature and molten salt inventory levels in the high-temperature tank at a 50 MW central receiver CSP plant, connected to the power grid in 2019. A multi-physics model was developed to evaluate thermal and mechanical properties of TES tanks by combining computational fluid dynamics and finite element modeling using real plant data. Heat loss, temperature, displacement, and stress distribution of the tank at different inventory levels were investigated. Results show that ambient air velocity near the tank roof reaches 2.14 m/s, much higher than 0.2 m/s near the wall. The temperatures of inventory fluid and tank are close, varying slightly at different levels due to thermal conduction and radiation. Because the heat loss strongly depends on temperature, the total tank loss remains nearly constant across inventory levels. Larger temperature gradients and thermal stresses are primarily localized along the tank floor edge and the air-salt interface. Notably, the maximum thermal stress at the tank edge is three times higher than that at the interface. The magnitude of total stress changes by less than 5 MPa with and without thermal load, indicating that high temperatures exert only a minor impact on tank stress. In contrast, thermal load significantly affects tank deformation, particularly at the roof edge, where values exceed 150 mm. Despite the large variation in molten salt levels, tank wall temperatures and displacements present a minor change, suggesting a weak correlation with inventory levels. In conclusion, the findings obtained in this study provide important insights on the TES tank that could be used to optimize tank design and operation strategies.

14 SOLAR ENERGY↗

DOE OSTI2025

Corrosion of 316 L stainless steel under the natural circulation of molten NaCl-MgCl 2 salt

The corrosion behavior of 316 L stainless steel (SS) was studied via the natural circulation of molten eutectic NaCl-MgCl 2 salt through a microloop. The post-corrosion tested 316 L SS microloop sections were characterized with microscopy techniques to determine the microstructural and microchemical changes that occurred at the alloy/salt interface. It was found that 316 L SS showed heterogeneous dissolution at the hot-leg, whereas deposition of corrosion products occurred at the cold-leg. For the first time, experimentally obtained molten salt flow-induced corrosion of 316 L SS results were combined with computational thermodynamic-kinetic models to validate the dissolution and deposition in terms of elemental compositional changes at the alloy/salt interface. The thermodynamic-kinetic modeling predicted that the heterogeneous dissolution of Cr and Fe from the hot-leg section of 316 L SS persisted throughout the salt circulation. The model also estimated that, despite Cr deposition starting earlier than Fe, the total redeposition of Fe is expected to be significantly greater than that of Cr over the circulation of salt. Furthermore, the modeling accurately predicted the subsurface enrichment of Mo which is attributed to the reduced Cr activity and the relatively higher diffusion rate of Mo within the alloy matrix. Here, the agreement between modeling and experimental results confirms that Fe chlorides dissolve at the hot-leg and subsequently deposit at the cold-leg due to activity changes driven by the thermal gradient. By contrast, Cr was not detected in the cold-leg deposits, which is attributed to its weaker temperature dependence on activity, limiting its redeposition under these conditions.

36 - MATERIALS SCIENCE↗

DOE OSTI2025

Mitigated Dissimilar Mass Transfer with Aluminized RAFM Steels in Flowing PbLi at 650°C

To evaluate dissimilar material interactions in flowing eutectic Pb-16at. %Li between fusion-relevant materials and to establish a maximum operating temperature of future blanket designs, a series of thermal convection loop (TCL) experiments with flowing PbLi have been conducted. Following a 1000-h, 550°C to 650°C TCL experiment, a thin reaction product was observed on the surface of monolithic, high-purity SiC. To identify the source of the observed dissimilar material interaction and to understand the reaction kinetics, an identical 2000-h TCL exposure was conducted, but less reaction of the SiC specimen was observed. Characterization of the FeCrAlMo (alloy APMT) loop tubing and the PbLi-exposed austenitic stainless steel (SS) assembly parts suggested that the formation of the reaction layer on the surface of the PbLi-exposed SiC did not originate from the coated reduced activation ferritic martensitic (RAFM) steel, but from the dissolved loop materials (liquid metal exposed APMT tubing and type 316 SS fittings). Here, the results, after 2000 h, suggest that the dissimilar material reaction between the Al-coated RAFM steel and the SiC has very slow reaction kinetics within a temperature gradient of 550°C to 650°C.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗