SEARCH · Search NASA
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.
Using neutron scattering dynamics measurements to validate molecular dynamic simulations, it is found that the rates of solvent exchange are not likely to limit the rates of calcite crystal growth and dissolution.
This study systematically investigates the impact of increased upper limit voltage in the reliability and degradation of a scaled vanadium redox flow battery over long-term testing (500 cycles).
Limit-cycle oscillations are studied for ion temperature gradient turbulence, which, in the absence of large diamagnetic (mean) shear flows, saturates through energy transfer from unstable modes to large-scale stable modes via zonal-flow intermediary modes. Oscillations of zonal flow and turbulence levels are strongly constrained by the reactive, largely non-dissipative character of the zonal flows. Since existing predator–prey models for observed oscillations in experiments do not include energy transfer through zonal flows to stable modes, low-order fluid models with this physics are constructed and investigated. A simple three-wave truncation produces low-amplitude zonal flows that slowly oscillate around a zero mean, with turbulence oscillations between coupled wavenumbers that exceed linear frequencies by orders of magnitude. This inconsistency with experimental observations is caused by the weak non-linear drive of zonal flows in three-wave systems and the lack of multiple-wavenumber turbulent interactions. A more comprehensive model that preserves multiple wavenumber interactions within the context of conservative zonal-flow-mediated energy transfer to stable modes accurately reflects observed dynamics when the phase between stable and unstable modes is occasionally randomized.
Advanced tokamak scenarios often feature equilibriums with zero and reversed magnetic shear. To isolate and investigate their impacts on the resistive wall mode (RWM) instability analytically, we construct a series of cylindrical limiter equilibriums with different magnetic shear in the core and toward plasma edge, as a prototype of the configurations in advanced tokamak scenarios. Equilibriums with zero β are considered so that the analysis focuses on the current-driven RWMs. Analytical solutions for the n=1 resistive wall mode, derived from the reduced ideal magnetohydrodynamics (MHD) equations, consistently demonstrate key insights for both the core and the edge: the growth rate of RWMs is inversely proportional to the averaged magnetic shear when the safety factor at the edge (qa) is held constant. Specifically, an increase in positive averaged magnetic shear results in a decreased growth rate, whereas an increase in negative averaged magnetic shear leads to an increased growth rate. Furthermore, the influence of nqa on the mode's stability significantly surpasses the effects of magnetic shear changes previously discussed. These analytical calculation results provide a verification of NIMROD simulations.
In this work, we present the development of robust surrogate models for Ion Cyclotron Range of Frequencies (ICRF) and High-Harmonic Fast Wave (HHFW) heating predictions in fusion plasmas. Building upon our previous efforts to achieve real-time capable models, we identify the cause of the outliers found using TORIC in certain HHFW heating scenarios. The outliers are observed to be spurious ion Bernstein wave (IBW)-like modes caused by a wavelength control algorithm designed to address challenging scenarios with high perpendicular wavenumbers. The effect arises from the modulation in the perpendicular susceptibility, which can induce sign reversal and IBW-like propagation for scenarios featuring normalized ion Larmor radius λ i ≫ 1. We use TORIC with this algorithm disabled to generate a novel HHFW-NSTX database that is free of outliers. Surrogate models trained on this database, including Random Forest Regressor (RFR), Multi-Layer Perceptrons, and Gaussian Process Regressors (GPR), demonstrate the ability to accurately predict HHFW heating profiles, with regression scores of R 2 ∈[0.93−0.99]. Additionally we demonstrate that it is possible to generalize predictions beyond training data by the use of both RFR and GPR models, enabling the prediction of scenarios previously limited to the original model. GPR models also provide uncertainty quantification, offering insights into model confidence. This work introduces a comprehensive Verification, Validation, and Uncertainty Quantification methodology for surrogate modeling, applicable not only to ICRF heating but also to other RF heating challenges and fusion physics problems. Beyond accelerated inference, these models show effective extrapolation capabilities, providing an alternative for addressing numerical challenges.
We show that trapped ions in virtual cathode potential wells can raise the transmitted current of emitted electrons into a plasma much closer to the full emission than is predicted by cathode sheath theories without trapped ions. The transmitted current is controlled by the well barrier voltage, which must adjust to balance the creation of low-energy ions within the well, and their loss. Our model considers the case of a plasma-facing cathode where trapped ions are created passively via charge-exchange collisions and lost passively via thermal leakage over the well. We quantify these rates and estimate the current in terms of system parameters for thermionic emission into a plasma with several cathode geometries. A general prediction is that the current as a function of emitted flux does not saturate at the traditional space charge limit (the onset of a well) but can reach far higher values until the trapped ion balance breaks down, causing instability. The maximum stable current depends on parameters but in principle can be arbitrarily high if active techniques are used to manipulate the trapped ion balance. We conclude that major improvements in plasma technologies with hot cathodes might be achieved by optimizing the current enhancement enabled by trapped ions.
Palladium (Pd) catalysts have been extensively studied for the direct synthesis of H 2 O through the hydrogen oxidation reaction at ambient conditions. This heterogeneous catalytic reaction not only holds considerable practical significance but also serves as a classical model for investigating fundamental mechanisms, including adsorption and reactions between adsorbates. Nonetheless, the governing mechanisms and kinetics of its intermediate reaction stages under varying gas conditions remain elusive. This is attributed to the intricate interplay between adsorption, atomic diffusion, and concurrent phase transformation of catalyst. Herein, the Pd-catalyzed, water-forming hydrogen oxidation is studied in situ, to investigate intermediate reaction stages via gas cell transmission electron microscopy. The dynamic behaviors of water generation, associated with reversible palladium hydride formation, are captured in real time with a nanoscale spatial resolution. Our findings suggest that the hydrogen oxidation rate catalyzed by Pd is significantly affected by the sequence in which gases are introduced. Through direct evidence of electron diffraction and density functional theory calculation, we demonstrate that the hydrogen oxidation rate is limited by precursors’ adsorption. These nanoscale insights help identify the optimal reaction conditions for Pd-catalyzed hydrogen oxidation, which has substantial implications for water production technologies. The developed understanding also advocates a broader exploration of analogous mechanisms in other metal-catalyzed reactions.
Fiber-optic bolometers (FOBs) based on a fiber-tipped silicon Fabry–Perot (FP) interferometric temperature sensor and a gold disk absorber have been shown to be an attractive alternative to conventional resistive bolometers for plasma radiation measurement in fusion devices. Either a high-finesse FP or a low-finesse FP can be used, each with trade-offs between noise performance and fabrication complexity. In this paper, we present an FOB design that overcomes these limitations by combining a low-finesse long silicon FP cavity with a large gold disk absorber to achieve enhanced sensitivity and noise performance without increasing the fabrication complexity and the time constant. We also demonstrated a fabrication method for the sensor head facilitated by thermal release tape and precise laser heating. Our FOB demonstrates a temperature resolution of 0.08 mK, a cooling time constant of 230 ms, and a noise equivalent power density of 0.015 W m −2 . This represents an eightfold improvement over previous high-finesse FOBs and 26-fold improvement over previous low-finesse FOBs with similar demodulation bandwidths and similar cooling time constants.
Abstract A hybrid Monte-Carlo molecular-dynamics method for determining solidus and liquidus compositions in multicomponent systems is presented that overcomes both the time limitations in conventional molecular dynamics that prevent the evolution of distinct solid and liquid compositions via diffusion and the complexity challenge that prevents use of thermodynamic assessment in systems of many components. This hybrid method is validated in the Cu–Ni system against an independent assessment of solidus and liquidus compositions based on the regular solution model. Strategies for efficient mapping of different phase diagrams, based on the thermodynamic parameter T 0 , the temperature at which two phases of the composition X 0 have equal free energies, are presented and then demonstrated for the copper-nickel fully miscible system and the gold–silicon eutectic system. A calculation of the solidus and liquidus sampled during the equilibrium melting of equiatomic CrMnFeCoNi is performed, indicating that this method has potential to be extended to the study of many component alloys.
Abstract We present a model for the particle balance in the post-disruption runaway electron plateau phase of a tokamak discharge. The model is constructed with the help of, and applied to, experimental data from TCV discharges investigating the so-called ‘low- Z benign termination’ runaway electron mitigation scheme. In the benign termination scheme, the free electron density is first reduced in order for a subsequently induced MHD instability to grow rapidly and spread the runaway electrons widely across the wall. We show that the observed non-monotonic dependence of the free electron density with the measured neutral pressure is due to plasma re-ionization induced by runaway electron impact ionization. At higher neutral pressures, more target particles are present in the plasma for runaway electrons to collide with and ionize. Parameter scans are conducted to clarify the role of the runaway electron density and energy on the free electron density, and it is found that only the runaway electron density has a noticeable impact. While the free electron density is shown to be related to the spread of heat fluxes at termination, the exact cause for the upper neutral pressure limit remains undetermined and an object for further study.
Big-bang nucleosynthesis (BBN) probes the cosmic mass-energy density at temperatures ~10 MeV to ~100 keV. Here, we consider the effect of a cosmic matter-like species that is non-relativistic and pressureless during BBN. Such a component must decay; doing so during BBN can alter the baryon-to-photon ratio, η, and the effective number of neutrino species. We use light element abundances and the cosmic microwave background (CMB) constraints on η and N ν to place constraints on such a matter component. We find that electromagnetic decays heat the photons relative to neutrinos, and thus dilute the effective number of relativistic species to N eff < 3 for the case of three Standard Model neutrino species. Intriguingly, likelihood results based on Planck CMB data alone find N ν = 2.800 ± 0.294, and when combined with standard BBN and the observations of D and 4 He give N ν = 2.898 ± 0.141. While both results are consistent with the Standard Model, we find that a nonzero abundance of electromagnetically decaying matter gives a better fit to these results. Our best-fit results are for a matter species that decays entirely electromagnetically with a lifetime τ X = 0.89 sec and pre-decay density that is a fraction ξ = (ρ X /ρ rad |10 MeV = 0.0026 of the radiation energy density at 10 MeV; similarly good fits are found over a range where ξτ X 1/2 is constant. On the other hand, decaying matter often spoils the BBN+CMB concordance, and we present limits in the (τ X ,ξ) plane for both electromagnetic and invisible decays. For dark (invisible) decays, standard BBN (i.e. ξ = 0) supplies the best fit. We end with a brief discussion of the impact of future measurements including CMB-S4.
We investigate two classes of dark matter (DM) candidates, sub-GeV particles and primordial black holes (PBHs), that can inject low-energy electrons and positrons into the Milky Way and leave observable signatures in the X-ray sky. In the case of sub-GeV DM, annihilation or decay into e + e - contributes to the diffuse sea of cosmic-ray (CR) leptons, which can generate bremsstrahlung and inverse Compton (IC) emission on Galactic photon fields, producing a broad spectrum from X-rays to γ-rays detectable by instruments such as eROSITA and XMM-Newton. For PBHs with masses below ∼10 17 g, Hawking evaporation similarly yields low-energy e ± , leading to comparable diffuse emission. Using the first data release from eROSITA and incorporating up-to-date CR propagation and diffusion parameters, we derive new constraints on both scenarios. For sub-GeV DM, we exclude thermally averaged annihilation cross sections in the range ∼ 10 -27 –10 -25 cm 3 /s and decay lifetimes of ∼ 10 24 –10 25 s for masses between 1 MeV and 1 GeV, with eROSITA outperforming previous X-ray constraints below ∼ 30 MeV. For asteroid-mass PBHs, we set new bounds on the DM fraction based on their Hawking-induced emission. Finally, we revisit previous constraints in 10.1088/1475-7516/2023/07/026JCAP 07 (2023) 026 and consequently 10.3847/1538-4357/ad41e0Astrophys. J. 968 (2024) 46 derived from XMM-Newton, finding that they were approximately four orders of magnitude too stringent due to the use of the instrument's geometric solid angle rather than its exposure-weighted solid angle. Upon using the exposure-weighted solid angle, we show that the revised XMM-Newton limits are slightly weaker than those from eROSITA.
At temperatures below the QCD phase transition, any substantial lepton number in the Universe can only be present within the neutrino sector. In this work, we systematically explore the impact of a non-vanishing lepton number on Big Bang Nucleosynthesis (BBN) and the Cosmic Microwave Background (CMB). Relying on our recently developed framework based on momentum averaged quantum kinetic equations for the neutrino density matrix, we solve the full BBN reaction network to obtain the abundances of primordial elements. We find that the maximal primordial total lepton number L allowed by BBN and the CMB is -0.12 (-0.10) ≤ L ≤ 0.13 (0.12) for NH (IH), while specific flavor directions can be even more constrained. This bound is complementary to the limits obtained from avoiding baryon overproduction through sphaleron processes at the electroweak phase transition since, although numerically weaker, it applies at lower temperatures and is obtained completely independently. We publicly release the C++ code COFLASY-C on GitHub (https://github.com/mariofnavarro/COFLASY/tree/COFLASY-C) which solves for the evolution of the neutrino quantum kinetic equations numerically.
We present a semi-automated algorithm for designing three-dimensional divertor or limiter plates targeting low heat loads. The algorithm designs the plates in two stages: firstly, the parallel heat flux distribution is caught on vertically-inclined plates at one or several toroidal locations. Secondly, the power per unit area is reduced by stretching, tilting and bending the plates toroidally. Heat transport is modelled using the EMC3-Lite code, which uses an anisotropic diffusion model. We apply this scheme to HSX, a medium-sized stellarator located at the University of Wisconsin–Madison. Starting from the current machine with an extended vessel wall, we construct plates which are able to effectively catch and spread the heat for three different magnetic configurations. The scheme has a computational cost in the order of tens of CPU-minutes, making it a powerful tool for semi-automated plasma-facing component design in three-dimensional environments.
We present a simplified energy-conserving approach to incorporate wave–wave nonlinear effects within the framework commonly used to describe wave–particle nonlinearities. In particular, the effects of zonal mode (ZM) generation on the determination of the saturation amplitude of energetic particle (EP)-driven Alfvénic instabilities is studied. The model assumes that the zonal perturbations grow at a rate twice that of the original (pump) wave, consistent with a beat-driven (or force-driven) generation mechanism. The evolution and saturation of the mode amplitude are investigated both analytically and numerically within our reduced model assumptions, in both the collisionless and scattering-dominated regimes. These studies underscore the crucial role of sources and sinks in capturing the impact and the role of beat-driven zonal perturbations on mode evolution. In the realistic case of saturation set by sources and sinks, we discuss the role of a finite amplitude ZM in reducing microturbulent particle scattering, thus limiting the energy source for the resonant mode. We then discuss comparisons between the model’s predictions and simulation results. The model reproduces key features observed in gyrokinetic simulations as the reduction in saturated mode amplitude and the onset of wave–wave nonlinear effects as functions of mode growth rate and amplitude. Thanks to its simplicity, it can be readily implemented into codes based on reduced models, thereby improving their predictive capability for strongly driven instabilities.
The ITkPixV2 chip is the final production readout chip for the ATLAS Phase 2 Inner Tracker (ITk) upgrade at the upcoming High-Luminosity LHC (HL-LHC). Due to the extraordinarily high peak luminosity at the HL-LHC of 5 × 1034 cm-1 s-1, ITkPixV2 must meet significant increases in nearly all design requirements compared to the current ATLAS Inner Detector (ID), including a 10× increase in trigger rate, a 7.5× increase in hit rate, a 3× increase in radiation tolerance, and a 12.5× decrease in pixel current draw per unit area, all while maintaining a similar power per unit area as present pixel detectors. Here we present the first measurements of the ITkPixV2 chip operated at the limits of the full chip design requirements, including in particular a measurement of the activity-induced current of the chip as a function of increasing hit rate.
China has large, estimated potential for direct air carbon capture and storage (DACCS) but its deployment locations and impacts at the subnational scale remain unclear. This is largely because higher spatial resolution studies on carbon dioxide removal (CDR) in China have focused mainly on bioenergy with carbon capture and storage. This study uses a spatially detailed integrated energy-economy-climate model to evaluate DACCS for 31 provinces in China as the country pursues its goal of climate neutrality by 2060. We find that DACCS could expand China’s negative emissions capacity, particularly under sustainability-minded limits on bioenergy supply that are informed by bottom-up studies. But providing low-carbon electricity for multiple GtCO 2 yr -1 DACCS may require over 600 GW of additional wind and solar capacity nationwide and comprise up to 30% of electricity demand in China’s northern provinces. Investment requirements for DACCS range from 330 to 530 billion dollars by 2060 but could be repaid manyfold in the form of avoided mitigation costs, which DACCS deployment could reduce by up to $6 trillion over the same period. Enhanced efforts to lower residual CO 2 emissions that must be offset with CDR under a net-zero paradigm reduce but do not eliminate the use of DACCS for mitigation. For decision-makers and the energy-economy models guiding them, our results highlight the value of expanding beyond the current reliance on biomass for negative emissions in China.
Permafrost peatlands store vast amounts of frozen carbon across northern landscapes. When ground ice melts, surface subsidence produces thermokarst landforms that expand wetlands at the edges of permafrost plateaus. Thermokarst represents an accelerating climate feedback, but uncertainties remain about how ground ice, hydrology, and vegetation interact to shape landscape change and carbon fluxes. We extended the process-based model ecosys to simulate thermokarst dynamics in laterally coupled 2D transects at a well-characterized boreal peatland site in Canada’s Northwest Territories. After benchmarking against site observations, we varied ground ice content and hydrologic boundary conditions across ranges typical near the southern permafrost limit. Simulations revealed distinct degradation regimes governed by the elevation difference between the frost table and the external water table. Rates of lateral retreat, the thaw-driven encroachment of wetlands into adjacent plateaus, ranged from 0 to >2 m yr −1 under identical weather forcing, consistent with observations and highlighting the strong role of WT and ground ice. Simulated vegetation dynamics indicate that black spruce mortality cannot be explained by anoxia alone, pointing to additional stressors such as root damage, pathogens, or physical destabilization. Despite large hydrologic shifts, net ecosystem CO 2 exchange remained a slight sink after collapse, while methane (CH 4 ) emissions rose by one to two orders of magnitude. As a result, lateral retreat substantially increases the greenhouse warming potential of permafrost peatlands (1.7 million km 2 in area), with simulated emissions of 0.1–10 Mt CO 2 -eq decade −1 depending on hydrology and retreat rates. These results underscore the need to account for both ground ice and hydrologic dynamics when assessing thermokarst-driven climate feedbacks.