SEARCH · Search NASA

Results for “Extended Kalman Filter”

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.

139 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Dual atom catalysts for rapid electrochemical reduction of CO to ethylene

Strong CO adsorption and facile CO dimerization are the key challenges in electrochemical CO2 reduction towards multi-carbon (C2+) products. We recently showed that CoPc immobilized on a single-walled carbon nanotube can selectively reduce CO2 to methanol. This is enabled through molecular strain, which dramatically improves the CO adsorption energy to CoPc, which in turn facilitates methanol formation. We now examine the extended Phthalocyanine (PcEx) dual atom catalyst (DAC), which is intrinsically strained and contains two catalyst centers, making it a candidate for reducing CO to C2+ products. Using Quantum Mechanics (QM), we screened 20 elements embedded in the PcEx, seeking catalysts with weak hydrogen binding, strong CO binding, and facile CO dimerization. We identi>ied Fe, Ru, Co, and Ir as the best performers and subsequently evaluated the entire CO to C2H4 mechanism (9 steps) using each of these elements as catalysts. In terms of limiting potential and overall exergonicity, we identi>ied CoPcEx as the best catalyst, followed by IrPcEx. We then examined the full CO to C2H4 mechanism on the bimetallic IrCoPcEx catalyst using grand canonical QM to obtain the reaction energetics as a function of applied potential. We conclude that the bimetallic IrCoPcEx is most promising for ef>iciently converting CO to ethylene.

Musgrave, Charles B.↗

DOE OSTI2024

Cosmological probes of Dark Radiation from Neutrino Mixing

Abstract Models of stepped dark radiation have recently been found to have an important impact on the anisotropies of the cosmic microwave background, aiding in easing the Hubble tension. In this work, we study models with a sector of dark radiation with a step in its abundance, which thermalizes after big bang nucleosynthesis by mixing with the standard model neutrinos. For this, we extend an earlier work which has focused on the background evolution only until the dark sector thermalizes by deriving the full background and perturbation equations of the model and implementing them in an Einstein-Boltzmann solving code. We expound on the behavior of this model, discussing the wide range of parameters that result in interesting and viable cosmologies that dynamically generate dark radiation during a range of epochs. We find that for the strongly self-coupled regime, there is no large cosmological impact for a tight prior on the mass, whereas larger mass ranges allow a smooth interpolation between a behavior close to the ΛCDM cosmological standard model and close to an additional component of strongly self-interacting dark radiation. In the weakly self-coupled regime we find that we can accommodate a parameter space relevant for the neutrino anomalies as well as one relevant to easing the Hubble tension.

Astronomy & Astrophysics↗

DOE OSTI2025

Strong anharmonicity dictates ultralow thermal conductivities of type-I clathrates

Type-I clathrate solids have attracted significant interest due to their ultralow thermal conductivities and sub- sequent promise for thermoelectric applications, yet the mechanisms underlying these properties are not well understood. Here, we extend the framework of vibrational dynamical mean-field theory (VDMFT) to calculate temperature-dependent thermal transport properties of solids using a many-body Green’s function approach. When applied to a coarse-grained description of 𝑋 8 Ga 16 Ge 30 , where 𝑋= Ba, Sr, we find that nonresonant scattering between cage acoustic modes and rattling modes leads to a reduction of acoustic phonon lifetimes and thus thermal conductivities. Moreover, we find that the moderate temperature dependence of conductivities above 300 K, which is consistent with experimental measurements, cannot be reproduced by textbook perturbation theory calculations, which predict a 𝑇 −1 dependence. Therefore, we suggest that nonperturbative anharmonic effects, including four- and higher-phonon scattering processes, are responsible for the ultralow thermal conductivities of type-I clathrates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Dynamic control of quantum phases in two-dimensional materials via Floquet engineering

The dynamical engineering of quantum states through periodic optical driving, known as Floquet engineering, has emerged as a powerful frontier in condensed matter physics, offering a pathway to realize material properties inaccessible in static equilibrium. This review provides a comprehensive overview of recent theoretical and experimental advances in the optical manipulation of two-dimensional (2D) quantum materials. We begin by systematically reviewing the evolution of the field from its pioneering applications in graphene and twisted moiré superlattices, highlighting the experimental realization of the light-induced anomalous Hall effect (AHE) to the complex spin-valley physics in transition metal dichalcogenides (TMDs). Furthermore, we briefly examine recent advances in 2D magnetic materials, demonstrating how optical driving can actively compete with intrinsic magnetism to dynamically switch magnetic orders and topological invariants. Moreover, we discuss the emerging frontiers of multi-frequency driving, quantum optimal control theory (QOCT), and ultrafast lightwave electronics. We highlight how tailored waveforms, such as bicircular light fields, and sub-cycle attosecond control can selectively break spatial symmetries to generate novel nonlinear photocurrents, mitigate dissipation, and extend the boundaries of quantum control well beyond the perturbative steady-state regime. Finally, we summarize the key experimental challenges for Floquet engineering, including effects such as heating and scattering, which limit coherent quantum control.

Wang, Wenpeng [Northeastern University, Shenyang, ↗

DOE OSTI2026

Contrasting structural reversibility and magnetic correlations in isostructural honeycomb magnets CrCl3 and 𝛼−RuCl3

We report a comparative neutron single crystal diffraction study of the structural and magnetic properties of layered halides CrCl3 and 𝛼−RuCl3. They host a honeycomb arrangement of transition metal ions with distinct electronic configurations and undergo a first-order structural transition between high-temperature 𝐶⁢2/𝑚 and low-temperature 𝑅⁢‾‾‾3. Both compounds show a step-like change in the 𝑐-lattice parameter across the structure transition. In contrast, the in-plane lattice response is quite different: 𝛼−RuCl3 exhibits an abrupt hysteretic change across the transition accompanied by progressive crystalline degradation upon thermal cycling, whereas CrCl3 shows a smooth in-plane lattice evolution and remains structurally robust. Magnetically, CrCl3 orders into an A-type antiferromagnetic structure at 𝑇𝑁=14 K and exhibits pronounced diffuse magnetic scattering extending up to about 40 K. 𝛼−RuCl3 shows no observable magnetic diffuse scattering above its zigzag antiferromagnetic ordering temperature 𝑇𝑁=7.6 K. These results suggest that the contrasting structural responses arise from an interplay between interlayer sliding energetics, stacking-strain coupling, and elastic accommodation of the stacking transition. The distinct chemical bonding and electronic configurations of the two compounds provide a microscopic basis for their different lattice responses to the structure transition and magnetic correlations.

Morgan, Zachary [ORNL] (ORCID:0000000243625911)↗

DOE OSTI2026

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗

DOE OSTI2026

Evolution of the Antarctic Ice Sheet from 2000–2300 and beyond: model sensitivity and uncertainty analysis using MPAS-Albany Land Ice

We present a description of the Antarctic Ice Sheet model configuration submitted to the ISMIP6-Antarctica-2300 experiment using the MPAS-Albany Land Ice model, along with three new sets of simulations: (1) a set of extended simulations to 2500 for three forced experiments and to 2775 for the control experiment; (2) a sensitivity analysis of our model configuration to parameters controlling basal sliding and sub-shelf melt, and to model structural choices including the choice of the energy and stress balances; and (3) a 72-member ensemble run on graphics processing units (GPUs) and analysis of variance to determine the primary sources of uncertainty in our ice-sheet model projections. Our extended simulations predict rapid retreat beginning after 2300 for SSP1-2.6 forcing and after 2500 for present-day (control) forcing, primarily in the Amundsen Sea Embayment. We find that varying the sub-shelf melt parameter between the 5th to 95th percentile values for a mean-Antarctic calibration target results in an up to ∼ ± 40 % change in sea-level contribution relative to our baseline simulations that used the median value. Using a linear basal sliding law reduces sea-level contribution by 51 %–73 % relative to our baseline nonlinear sliding law with an exponent of 1/5. When using basal sliding law exponents of 1/3 and 1/10, the overall difference from our baseline simulations at 2300 is on the order of 10 %. The Amundsen Sea Embayment region displays a strongly non-linear dependence of mass loss on the sliding law exponent, with no discernible relationship between the sliding law exponent and the mass loss by 2300, while the sectors feeding the Ross and Filchner-Ronne ice shelves exhibit more mass loss with a more-plastic sliding law. Our model fidelity sensitivity experiments reveal a 9 %–31 % increase in sea-level contribution when using a depth-integrated stress balance approximation relative to our three-dimensional solver, while using a fixed-in-time temperature field increases sea-level contribution by 14 %–88 % relative to two thermomechanically coupled configurations. Our 72-member ensemble and analysis of variance show that the uncertainty in long-term projections is dominated by the choice of Earth system model forcing and the presence or absence of hydrofracture forcing, rather than uncertainty in sliding and sub-shelf melt parameters.

58 GEOSCIENCES↗

DOE OSTI2026

Reduced-order CFD modeling of cryogenic hydrogen isotope extrusion for pellet fueling

This study presents a reduced-order model (ROM) for computational fluid dynamics (CFD) simulations of cryogenic hydrogen isotope extrusions, focusing on protium (H₂) and deuterium (D₂) piston extruders. Using a 2D axisymmetric ROM in ANSYS-Polyflow, significant computational savings were achieved (runtime reduced from 9∼24 h to 3∼5 min), with extrusion force discrepancies between the 2D ROM and 3D models being on the order of 1%. Parametric studies identified optimal cutoff shear rates in the viscosity model (0.01/s for H₂ and 0.001/s for D₂), providing recommendations for future simulations. Finally, a comprehensive comparison of ROM results with experimental data was performed across varying geometries, cryogenic materials, temperatures, extrusion lengths, and piston velocities. Predictions at low extrusion temperatures met the objective of providing quick and efficient solutions with an acceptable extrusion force error of approximately 10% or less, validating the effectiveness of the 2D ROM approach. However, at high temperatures closer to the triple point, extrusion force error grows, which necessitates developing an improved model that accounts for temperature effects, e.g. melting. Nevertheless, the findings still represent a significant improvement in efficiency of CFD modeling of cryogenic hydrogenic extrusion. The ROM framework can also be extended to tritium (T2) and screw extruders, which will ultimately provide a fast and effective tool for optimizing pellet injector design for ITER and future reactor systems.

Fan, Joy [ORNL] (ORCID:0000000229751735)↗

DOE OSTI2026

Nanoscale Compositional and Strain Gradients Enable High‐Speed and Amplitude‐Resolved Pyroelectric Sensing

The frequency response of pyroelectric sensors is fundamentally governed by thermal time constant (τth, determined by thermal mass and thermal conductance) and electrical impedance arising from film capacitance and readout circuit. Conventional bulk LiTaO3 detectors are optimized for high responsivity at low modulation frequencies (0.1-10 Hz), possessing a large τth that thermally averages rapid temperature oscillations at elevated modulation frequencies, limiting fidelity in resolving dynamic varying thermal signals. Here, compositional and strain gradients are introduced into 100-nm-thick relaxor-ferroelectric films reducing τth to ≈2 µs and producing built-in potentials (≈1.45 V or 145 kV cm-1) that enhance the pyroelectric coefficient and suppress the dielectric constant. This enables complementary dual-mode operation by enhancing current-mode electrical responsivity and improving the voltage-mode figure of merit - advantageous for superior temperature resolution (ΔTmin ≈ 30 µK). The responsivity peak shifts to near 1 kHz (>2500-times higher than conventional bulk sensors), with measurable responsivity extending to a carrier frequency of 100 kHz and amplitude-resolved detection at modulation frequencies up to 15 kHz. These results establish nanoscale internal-field engineering can reshape electro-thermal trade-off in pyroelectric thin films toward zero-bias, high-thermal-sensitivity, and amplitude-resolved thermal sensing across a wide frequency bandwidth.

Lin, Ching‐Che↗

DOE OSTI2026

Correlated Ion Transport Governed by Dynamic Local Structure in High Concentration and Localized High Concentration Electrolytes

Understanding the dynamics of cluster formation and network percolation provides the mechanistic link between microscopic solvation structure and transport in concentrated electrolytes, including localized high-concentration electrolytes (LHCEs). Although recent studies have shown LHCEs to form micelle-like aggregates at specific compositions, a quantitative understanding of how solvation structures and transport properties depend on salt–solvent–diluent ratios remains limited. Here, we integrate molecular dynamics with Onsager transport analyses to chart the evolution of solvation microstructure and associated ionic transport in LiFSI/DMC/TTE, achieving good agreement with experimental conductivities across composition. We show that micelle formation arises from a dynamic instability of the cation–anion network, and that network percolation is the primary determinant of conductivity. Two compositional thresholds emerge: A critical network concentration (CNC) and a critical micelle concentration (CMC) that delineate transitions from extended percolating networks to micelle-like clusters and then to fragments. These structural transitions rationalize the nonintuitive conductivity decrease at intermediate dilution despite monotonically decreasing viscosity and provide composition-level design rules for LHCEs.

Mohanakrishnan, Rohith Srinivaas [University of Ca↗

DOE OSTI2025

Cold Spray Cobalt Magnetostrictive Electromagnetic Acoustic Transducers for High Temperature Structure Monitoring

The Department of Energy’s Advanced Sensors and Instrumentation program seeks to develop and qualify advanced sensors for the nuclear industry. Reliable high temperature and high radiation sensors for detection and characterization of structural flaws in pipes, vessels, and structurally critical components is a weakness for both conventional light water reactors with coolant T-hot approaching 350oC, and for advanced reactors with T-hot temperatures in excess of 500 to 800oC. Magnetostrictive Electromagnetic Acoustic Transducers using a cold spray cobalt coating have been proposed as a sensor design that can withstand these kinds of temperatures and radiation levels to serve as online sensors to detect flaws before cracks, pits, or erosion/corrosion flaws progress to through-wall failures. This report tests cold spray cobalt as part of a magnetostrictive EMAT for high temperature service. Cobalt is known to have strong magnetostrictive properties however the effect of cold spray application is not well studied. This program was surprised to discover that cold sprayed cobalt exhibited little or no magnetostrictive behavior until it was thermally annealed. Following annealing to 650oC however, cold spray cobalt did exhibit a magnetostrictive response. Work to date prior to this milestone report publication showed that magnetostrictive EMAT was successfully tested to 400oC with an alnico permanent magnet. The program plans to extend testing with an electromagnet to higher temperatures. This follow-on work will be reported under subsequent publications or as a revision to this report.

Glass, Samuel W.↗

DOE OSTI2024

Resonant Metasurfaces with Van Der Waals Hyperbolic Nanoantennas and Extreme Light Confinement

This work reports on a metasurface based on optical nanoantennas made of van der Waals material hexagonal boron nitride. The optical nanoantenna made of hyperbolic material was shown to support strong localized resonant modes stemming from the propagating high-k waves in the hyperbolic material. An analytical approach was used to determine the mode profile and type of cuboid nanoantenna resonances. An electric quadrupolar mode was demonstrated to be associated with a resonant magnetic response of the nanoantenna, which resembles the induction of resonant magnetic modes in high-refractive-index nanoantennas. The analytical model accurately predicts the modes of cuboid nanoantennas due to the strong boundary reflections of the high-k waves, a capability that does not extend to plasmonic or high-refractive-index nanoantennas, where the imperfect reflection and leakage of the mode from the cavity complicate the analysis. In the reported metasurface, excitations of the multipolar resonant modes are accompanied by directional scattering and a decrease in the metasurface reflectance to zero, which is manifested as the resonant Kerker effect. Van der Waals nanoantennas are envisioned to support localized resonances and can become an important functional element of metasurfaces and transdimensional photonic components. By designing efficient subwavelength scatterers with high-quality-factor resonances, this work demonstrates that this type of nanoantenna made of naturally occurring hyperbolic material is a viable substitute for plasmonic and all-dielectric nanoantennas in developing ultra-compact photonic components.

Chemistry↗

DOE OSTI2025

Decoupling Failure Pathways in Li-O 2 Cells Operated Under Lean Electrolyte Conditions

The operation of lithium-oxygen (Li-O 2 ) batteries under lean electrolyte conditions offers higher energy density but leads to rapid degradation and short cycle life. Although cathode passivation and electrolyte decomposition occur in all regimes, we show that under lean electrolyte conditions, failure is primarily driven by progressive electrolyte consumption at the lithium/solid electrolyte interphase (SEI), rather than by irreversible cathode passivation. Poor wetting of the lithium surface results in heterogeneous SEI growth and high local current densities, which accelerate electrolyte loss and cell failure. Strategies aimed at improving interfacial stability, including optimized wetting and SEI forming additives, significantly extend cycle life without compromising energy density. Our results establish anode-electrolyte interactions as the dominant degradation mechanism under lean electrolyte conditions, and emphasize the need to engineer a stable Li/SEI interface for long-lasting Li-O 2 batteries.

Córdoba, Daniel [Argonne National Laboratory (ANL)↗

DOE OSTI2025

Exploiting ν-dependence of projected energy correlators in HICs

We extend the recently derived factorization formula for energy-energy correlators to study the analytic structure of general ν-point projected energy correlators in heavy ion collisions. The ν-point projected energy correlators (or, ν-correlators) are an analytically continued family of the integer N-point projected energy correlators, which probe correlations between N final-state particles. By tracking the largest separation (χ) between the N particles, in vacuum, their structure is closely related to the DGLAP splitting functions and exhibits a classical scaling behavior ∼ 1/χ which is modified by resummation through the anomalous dimensions. We show that, in a thermal medium, the ν-correlators display non-trivial angular scaling already at the leading order in perturbation theory. We find that for non-integer values, particularly ν < 1, medium-induced jet function is enhanced compared to ν > 1. This is particularly manifested in the ratios of ν-correlators with respect to the two-point energy correlator, which encodes intrinsic angular information for ν < 1 when compared to large ν values. Moreover, for small-ν values, the ν-correlators appear to saturate at ν = 0.01 . We further confirm our leading-order numerical computations against simulated events from JEWEL for the parton-level production cross-section. Finally, we qualitatively discuss the effect of BFKL resummation for various values of ν.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Strong gradient neoclassical transport in the plateau regime

Strong gradient regions in tokamaks such as the pedestal or internal transport barriers are regions of reduced turbulence where neoclassical transport can play a dominant role. In pedestals, gradient lengths comparable to the ion poloidal gyroradius have been measured. Standard neoclassical theory can miss important strong gradient effects in these regions because it assumes that the gradient length scales of density, temperature and potential are larger than the ion poloidal gyroradius. We extend plateau regime neoclassical theory into regions of gradients of the order of the ion poloidal gyroradius to capture strong gradient effects on transport processes in the pedestal and internal transport barriers. The fundamental idea behind our new framework is to keep a scale separation between the orbit widths and the gradient length scales by performing a large aspect ratio expansion. In the plateau regime, strong gradients cause poloidal variation that is in–out as well as up–down asymmetric. We study two different test cases assuming either radial force balance or the absence of turbulence and show that strong gradient effects can enhance or reduce standard neoclassical theory predictions in the plateau regime in strong gradient regions.

fusion plasma↗

DOE OSTI2026

Equation of state for Hf, Ta, W, Re, Os, Ir, Pt, and Au to multi-terapascal pressures from density-functional theory

We present the zero-temperature equation of state (pressure dependence of compression) and phase stability predictions for the 5d-transition metals obtained from all-electron density-functional theory (DFT) calculations. The results compare favorably with experiments but extend beyond current experimental capabilities to 10 TPa. Our study reveals phase changes that are explained from the calculated electronic structure. The cubic face-centered and body-centered structures (fcc and bcc), together with two-, three-, and four-layered hexagonal structures, play major roles under compression. The results’ dependence on the electron exchange and correlation in the DFT approach is investigated, and it is shown that the impact of the choice, while significant at lower pressures, diminishes in the terapascal regime. We further illustrate that the normal parabolic trends in atomic volume and bulk modulus with atomic number, due to the occupation of bonding and anti-bonding 5d states, break down at TPa pressures, suggesting drastically different chemical bonding at these extreme conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

HTGR Multiphysics Application Drivers FY26 Updates

This report summarizes FY26 progress under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program's high-temperature gas-cooled reactor (HTGR) application driver work, covering a wide range of activities such as code validation and multi-physics code assessment. 1) A detailed SAM model of the High-Temperature Engineering Test Reactor (HTTR) was developed using a unique-block grouping approach, with an extended parallel thermal network method to capture block-to-block conduction and radiation heat transfer, and applied to steady-state simulations of the HTTR 30~MW and 9~MW cases. 2) In another activity, SAM's newly implemented multi-component gas flow model was validated against the Natural convection Shutdown heat removal Test Facility (NSTF) argon ingress experiment, correctly capturing the density-driven suppression and thermal recovery of natural circulation observed when argon is introduced into the air-cooled Reactor Cavity Cooling System (RCCS) loop. 3) For the OECD/NEA High Temperature Test Facility (HTTF) benchmark, we co-led the international benchmark activities as well as the OECD/NEA final benchmark report to be released at the end of this year. 4) Finally, the coupled Griffin-SAM modeling capability for pebble-bed HTGRs was advanced by verifying the Griffin neutronics solution against Serpent Monte Carlo for a realistic non-uniform temperature distribution, resolving several deficiencies in the SAM-to-Griffin temperature transfer scheme, and enabling distinct fuel kernel, moderator, and coolant temperatures for cross section feedback. These new features were demonstrated in a PBR load-following transient.

Lee, Alvin↗

DOE OSTI2026

Transition from Vehicular to Structural Ionic Transport in Electrified Alkali Aqueous Solutions

A molecular understanding of the solvation and dynamics of ions under static electric fields is crucial for modelling a wide range of natural and technological processes. Yet, traditional simulation methods suffer from a trade-off that has to be made between accuracy and statistical convergence. To bridge this gap, herein we extend our recently introduced Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) approach to investigate the solvation structures and ionic transport mechanisms of electrified alkali cationic solutions. We obtain ionic conductivities for Li+, Na+ and Cs+ from the field dependence of the ionic current density in good agreement with experiment. Surprisingly, the migration mechanism is found to be strikingly different for the three ions despite their similar ionic conductivities. While Li+ conducts predominantly through vehicular migration of a stable 4-fold coordinated ion at all field strengths, Cs+ conducts strictly through a structural diffusion mechanism, where 9-12 transient first shell water coordination bonds are broken and renewed. Notably, aqueous Na+ emerges as a “Goldilocks” ion: its ion-water interactions are strong enough to maintain distinct 5-6 fold coordination shells at zero field (unlike Cs+) yet labile enough to be strongly perturbed by electric fields (unlike Li+). As a consequence, we observe an electric field-induced transition from vehicular to structural ionic transport for Na+ that is accompanied by a marked increase in the ionic current density. Our results imply that the conductance mechanism of ions with moderate ion-solvent interactions can be effectively tuned by external electric fields.

Joll, Kit↗