Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “state redistribution”

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 253 records · Page 14

Modeling High Current Pulsed Discharge in AA Battery Cathodes: The Effect of Localized Charging during Rest

During high current operation, substantial heterogeneity develops within battery cathodes, particularly when their thickness is large. Heterogeneity relaxation during subsequent rest is important for understanding battery performance under pulsed conditions. Localized charge balancing phenomena within batteries at zero net current are not well understood and merit investigation. In this work, the heterogeneity within cathodes of commercial alkaline Zn–MnO 2 batteries is measured during discharge and monitored during rest using energy dispersive X-ray diffraction (EDXRD). Significant gradients in protonation form during discharge and partially relax under rest. It is demonstrated that the proton gradient relaxation is through local redox activity at zero net current, where local (de)protonation works to redistribute charge across the cathode thickness. To support this redox-based relaxation, a fundamental kinetic study on prismatic MnO 2 cathodes is conducted to determine an appropriate model to describe both discharge and charge kinetics of MnO 2 . These kinetics are incorporated into a computational model to simulate the proton gradient formation and partial relaxation under identical discharge conditions as the operando EDXRD experiments. Model and experimental data are found to be in excellent agreement, correctly predicting localized charge balancing at rest.

Batteries↗

Understanding electronic structure tunability by metal dopants for promoting MgB 2 hydrogenation

Hydrogen is a promising energy carrier, but its onboard application is limited by the need for compact, low-pressure storage solutions. Solid-state complex metal hydride systems, such as MgB 2 /Mg(BH 4 ) 2 , offer high storage capacities but suffer from sluggish kinetics and poor reversibility. One avenue for improving reactivity is to introduce metal dopants to alter electronic and atomic properties, but the role of these chemical additives remains poorly understood, particularly for the hydrogenation reaction. In this work, we used density functional theory calculations on model MgB 2 systems to rationalize the potential role of metal dopants in destabilizing B–B bonding within the MgB 2 lattice. We carried out detailed electronic structure analyses for 28 different metal dopant adatoms to identify properties that contribute to a dopant’s efficacy. Based on the simulation results, we propose that an intermediate ionic and covalent character of the bonds between adatoms and B atoms is desirable for facilitating charge redistribution, disrupting the B–B bond network, and promoting H 2 dissociation and H atom chemisorption on MgB 2 .

08 HYDROGEN↗

Solute Redistribution During Aging in a Laser-Based Powder Bed Fusion-Produced Fe-18Ni-9Co-5Mo-Ti Maraging Steel

Solute distribution during the laser-based powder bed fusion of metals process is highly complex because of solidification. This study builds on a previous investigation of how solute heterogeneities across different scales influence phase distribution and aging physical phenomena in a Fe-18Ni-9Co-5Mo-Ti maraging steel (Santana et al. in Addit Manuf 104494, 2024). Previous work using energy-dispersive X-ray spectroscopy revealed heterogeneous Ni and Mo enrichment at intercellular regions. Using X-ray diffraction, 8 ± 3% retained austenite was quantified in as-built samples. Additionally, electron backscatter diffraction analysis exposed the location of retained austenite along cell boundaries. In this work, scanning transmission electron microscopy results revealed the presence of non-enriched retained austenite in the as-built condition. With increased aging temperature, Ni and Mo enrichment in austenite increased, especially at the austenite/martensite boundary. Solute diffusion was simulated in DICTRA Thermo-Calc® software to assess solute redistribution between the matrix and second phases and to comprehend whether elemental diffusion processes are interconnected. Three different scenarios were considered: (1) solute diffusion in the matrix; (2) diffusion with enriched austenite in the intercellular region; (3) diffusion with non-enriched austenite in the intercellular region. Results showed that Ni partitioning controls austenite growth and that boundary regions exhibit higher element enrichment than the core.

42 ENGINEERING↗

QMMM 2023: A program for combined quantum mechanical and molecular mechanical modeling and simulations

Combined quantum mechanical and molecular mechanical (QM/MM) methods play an important role in multiscale modeling and simulations. QMMM 2023 is a general-purpose program for single-point calculations, geometry optimizations, transition-state optimizations, and molecular dynamics (MD) at the QM/MM level. It calls a QM package and an MM package to perform the required single-level calculations and combines them into a QM/MM energy by a variety of schemes. QMMM 2023 supports GAMESS-US, Gaussian, and ORCA as QM packages and Tinker as the MM package. Four types of treatments are available for embedding the QM subsystem in the MM environment: mechanical embedding with gas-phase calculations of the QM region, electronic embedding that allows polarization of the QM region by the MM environment, polarizable embedding for mutual polarization of the QM and MM regions, and flexible embedding for both mutual polarization and partial charge transfer between the QM and MM regions. Boundaries between QM and MM regions that pass through covalent bonds can be treated by several methods, including the redistributed charge (RC) scheme, redistributed charge and dipole (RCD) scheme, balanced-RC scheme, balanced-RCD scheme, screened charge scheme that takes account of charge penetration effects, and smeared charge scheme that delocalizes the MM charges near the QM–MM boundary. Geometry optimization can be done using the optimizer implemented in QMMM 2023 or the Berny optimizer in Gaussian through external calls to Gaussian. Molecular dynamics simulations can be performed at the pure-MM level, pure-QM level, fixed-partitioning QM/MM level, and adaptive-partitioning QM/MM level. As a result, the adaptive-partitioning treatments permit on-the-fly relocation of the QM–MM boundary by dynamically reclassifying atoms or groups into the QM or MM subsystems.

97 MATHEMATICS AND COMPUTING↗

Phase Stability and Electrochemical Performance of La-Site-Doped Li6La3Zr0.5Nb0.5Ta0.5Hf0.5O12 High-Entropy Garnets

We investigate La-site substitution in the high-entropy garnet Li6La3Zr0.5Nb0.5Ta0.5Hf0.5O12 (LLZNTH) using Ba2+, Sr2+, and Sm3+ to elucidate how dopant governs phase stability, Li-site distribution, and electrochemical behavior. X-ray diffraction shows that Sr2+ is incorporated homogeneously into the garnet lattice, whereas the larger Ba2+ and smaller Sm3+ ions partially exceed the structural tolerance, generating secondary phases. Nevertheless, the Sm-doped composition (x = 0.05) exhibits the highest room-temperature ionic conductivity (2.7 × 10–4 S cm–1). Neutron powder diffraction reveals that Sm substitution drives a redistribution of Li+ from the tetrahedral 24 d sites into the higher-mobility 96 h positions, enhancing the connectivity of the three-dimensional Li-ion migration network. A Sm-doping series (x = 0.01–0.05) further shows that only sufficiently high Sm levels induce this redistribution, whereas lower concentrations retain Li arrangements similar to the undoped garnet. Critical current density measurements demonstrate that La-site dopants also influence interfacial stability against Li metal, underscoring a trade-off between bulk transport enhancement and mechanical robustness. Collectively, these findings reveal that in high-entropy garnets improved ionic conductivity can originate not only from phase-pure structures but also from targeted modification of the Li sublattice, even when accompanied by secondary phases, offering a compositional design principle for garnet electrolytes.

Li, Chang [Mechanical Engineering, School of Scien↗

Plasmonic Hot-Carrier Generation and Catalysis in Ti 3 C 2 O 2 from Real-Time TDDFT Simulations

Photoinduced hot electrons are central to plasmon-driven catalysis. Atomically thin Ti 3 C 2 O 2 , with high carrier density and broad optical absorption, offers a promising platform for plasmon-driven reactions. However, comprehensive investigations of its plasmon resonance, hot-carrier generation, and plasmonic catalytic performance remain limited. In this work, real-time time-dependent density functional theory (rt-TDDFT) was employed to study Ti 3 C 2 O 2 ’s plasmon excitation and hot-carrier generation from nonradiative plasmon damping. The temporal evolution of the dipole moment reveals plasmon resonance in Ti 3 C 2 O 2 , followed by strong plasmon damping that redistributes the stored energy to generate hot carriers. Ti 3 C 2 O 2 with low oxygen vacancy concentration (O v -Ti 3 C 2 O 2 ) exhibits plasmonic behavior resembling the pristine surface, and the plasmon-generated hot electrons can markedly reduce the dissociation barrier of CO 2 at the oxygen vacancy. These findings provide fundamental insights into the plasmonic properties of Ti 3 C 2 O 2 and how they drive its catalytic performance in surface reactions, which is valuable for advancing plasmon-driven catalysis.

CO2 reduction↗

Temporal and spatial tuning of optical constants in praseodymium doped ceria by electrochemical means

Abstract Temporal and spatial tuning of the refractive index of optical thin films is desired for flat optics applications. The redistribution of mobile ions in mixed ionic-electronic conductors (MIEC) has been demonstrated to serve as a viable means for achieving optical tuning down to the nanoscale. Here we studied the dynamic range of the optical tuning achievable in the refractive index, in the MIEC oxide – Pr x Ce 1− x O 2− δ (PCO), for x = 0.1, 0.2 and 0.4, at 500 °C, by in-situ spectrophotometry. Significant increases in the modulation of both the imaginary and real optical constants in the visible and the adjacent spectra were obtained for increased doping levels. Device employing an electrochemical titration method was implemented to modulate the oxygen concentration, and thereby the optical transmission of PCO. Incorporation of a patterned top electrode allowed for the demonstration of spatial control of PCO thin film properties by in-situ video imaging of the optical switching process. The electrochemically induced optical state is shown to remain non-volatile upon quenching the device to room temperature under applied bias.

36 MATERIALS SCIENCE↗

ARCANE Reweighting Demo

ARCANE reweighting is a Monte Carlo technique for tackling the negative weights problem in simulations in high energy physics. The technique works by redistributing (via an additive reweighting) the contributions of different pathways within an event generator that lead to the same final event. The technique is exact and does not introduce any biases in the distributions of physical observables. This repository contains the code and data associated with a demonstration of ARCANE reweighting, for the generation of $e^+ e^- \longrightarrow q\bar{q} + 1\,\mathrm{jet}$ events at next-to-leading-order accuracy in perturbation theory. ARCANE reweighting was shown to almost completely eliminate the negative weights problem for this process.

Shyamsundar, Prasanth [Fermi National Accelerator ↗

Evidence for Lignin–Carbohydrate Complexes from Studies of Transgenic Switchgrass and a Model Lignin–Pectin Composite

Lignin–carbohydrate complexes (LCCs) form through interactions of lignin with plant cell wall polysaccharides and are thought to be a significant source of biomass recalcitrance. In this work, we investigated LCCs formed between lignin and pectin homogalacturonan (HG). The structural changes in HG deficient transgenic switchgrass (GAUT4-knockdown, GAUT4-KD) after hot water pretreatment were compared to wild-type plants using small-angle neutron scattering (SANS), which showed that there were ~2.2-fold more lignin aggregates in GAUT4-KD biomass compared to the wild type. This demonstrated that decreased pectin resulted in more lignin redistribution and suggested that interactions between lignin and HG restrict lignin mobility in plant cell walls. To better understand the types of interactions between lignin and pectin, a model composite was prepared by polymerizing either protiated or partially deuterated coniferyl alcohol to form a dehydrogenation polymer (DHP) in the presence of HG. Small-angle X-ray scattering (SAXS) showed that the DHP and HG form a highly interconnected network structure that is not observed in a physical mixture of the individual polymers. Contrast matching SANS revealed the structure of DHP and HG in the composite and showed that the HG forms a swollen interconnected polymer network (power-law exponent, P = 1.5) interspersed with DHP particles (radius of gyration, Rg, 264 Å) that are composed of solvent-accessible DHP polymers (P = 2.3). Fourier transform infrared spectroscopy showed a unique ester absorption band in the DHP/HG composites. Solid-state nuclear magnetic resonance (NMR) analysis also supports interactions between DHP and HG. Overall, this study provides new insights into the relationship between primary and secondary cell wall polymers during cell wall synthesis and how LCCs formed between pectin and lignin could represent a previously unrecognized source of biomass recalcitrance. This knowledge may help develop new approaches to modulate cell wall properties to improve biofuel and bioproduct production.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Elucidation of Local Ordering and Atomic-Scale Structure in Polymer-Derived SiOC

Silicon oxycarbide (SiOC) is a versatile ceramic material with tunable microstructure and compositions that can be modulated through precursor chemistry and processing conditions. Though there are several noteworthy uses of SiOC across a range of application spaces, the difficulties in elucidating the short- to medium-range order within these materials have limited the maturation of strategies to precisely control SiC x O 4–x compositions for user-tailored applications. In this contribution, we implement a range of synchrotron scattering and spectroscopy methods coupled with stochastic modeling techniques to elucidate changes in local chemistry and structure associated with the pyrolysis of a commercially available SiOC polymer precursor. Stochastic modeling approaches provide valuable insights into decoupling local Si–O and Si–C environments while confirming predominate heterogeneous phases in materials. Using pyrolysis temperatures between 250 to 800 °C results in a heterogeneous material predominately composed of SiOC and amorphous SiO 2 domains. At 1100 °C, redistribution of Si–C pairs in the SiOC network and Si–O from the SiO 2 domains create a more ordered SiOC phase with local cubic SiC-like ordering. In addition, residual carbon leads to a detectable carbon phases at 1100 °C that persist at higher temperatures. These efforts address the difficulties of obtaining atomic-scale insights into the local structure and nanoscale heterogeneities in SiOC, providing pathways toward establishing structure–property relationships for future materials development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Multiscale Simulation Study of the Effects of Nantucket Island on Atmospheric Boundary Layer Flow

Multiscale atmospheric simulations were conducted to investigate how Nantucket Island modifies the marine atmospheric boundary layer (ABL) flow, how the mechanisms responsible are represented in mesoscale versus large-eddy simulation (LES) domains, and potential impacts on wind energy production downstream, in moderate and high wind speed cases. Here, to isolate the effects of the island, each case study was simulated (i) with and without island topography included and (ii) with observed versus spatially uniform sea surface temperature (SST) values, the latter removing the influence of spatial variability in SST, while also enabling the examination of how the island effects vary with near-surface thermodynamic stability over the surrounding ocean. Island effects were found to be strongest during daytime heating, which generates convective rolls that redistribute momentum vertically over the island and downstream. Colder SSTs and weaker winds strengthen these rolls and extend their influence. Flow statistics within a representative wind turbine rotor-swept area downstream show that while vertical momentum redistribution enhances near-surface wind speeds, those aloft decrease, resulting in minimal net change to rotor-equivalent wind speed and slightly reduced vertical shear. Comparisons of mesoscale and LES domains reveal that both capture the roughness-driven near-surface slowdown, but mesoscale models underestimate vertical momentum transport compared to LES, which explicitly resolves convective rolls. Increasing LES resolution beyond the grid spacing required to begin capturing the largest convective roll structures ∼ O(100) m yields minimal additional momentum transport but increases resolved turbulence kinetic energy, relevant for stress loading.

Atmosphere-land interaction↗

Distinct Surface and Bulk Thermal Behaviors of LiNi 0.6 Mn 0.2 Co 0.2 O 2 Cathode Materials as a Function of State of Charge

Understanding how structural and chemical transformations take place in particles under thermal conditions can inform designing thermally robust electrode materials. Such a study necessitates the use of diagnostic techniques that are capable of probing the transformations at multiple length scales and at different states of charge (SOC). In this study, the thermal behavior ofLiNi 0.6 Mn 0.2 Co 0.2 O 2 (NMC-622) was examined as a function of SOC, using an array of bulk and surface-sensitive techniques. In general, thermal stability decreases as lithium content is lowered and conversion in the bulk to progressively reduced metal oxides (spinels, rock salt) occurs as the temperature is raised. Hard X-ray absorption spectroscopy (XAS) and X-ray Raman spectroscopy (XRS) experiments, which probe the bulk, reveal that Ni and Co are eventually reduced when partially delithiated samples (regardless of the SOC) are heated, although Mn is not. Surface-sensitive synchrotron techniques, such as soft XAS and transmission X-ray microscopy (TXM), however, reveal that for 50% delithiated samples, apparent oxidation of nickel occurs at particle surfaces under some circumstances. This is partially compensated by reduction of cobalt but may also be a consequence of redistribution of lithium ions upon heating. TXM results indicate the movement of reduced nickel ions into particle interiors or oxidized nickel ions to the surface or both. These experiments illustrate the complexity of the thermal behavior of NMC cathode materials. The study also informs the importance of investigating the surface and bulk difference as a function of SOC when studying the thermal behaviors of battery materials.

25 ENERGY STORAGE↗

Intrinsic thermal conductivity of ZrC from low to ultrahigh temperatures: A critical revisit

Current phonon transport theory based on ground-state calculations has been successful in predicting thermal conductivity at room and medium temperatures but may misrepresent behavior at high temperatures. Here, in this work, we predict the thermal conductivity (κ) of ZrC including electronic and phonon contributions from 300 to 3500 K, by including high-order phonon scattering; lattice expansion; temperature-dependent (TD) second-, third-, and fourth-order force constants (2FC, 3FC, and 4FC); and interband phonon conduction by using first principles. For the phonon transport, we find that four-phonon scattering (4ph) significantly reduces the phonon thermal conductivity (κ ph ), by as much as ~75% at 3500 K. After including 4ph scattering and all other factors, κ ph shows a ~T –1.5 rather than ~T –1 dependence. TD 2FC decreases three-phonon scattering rates but increases 4ph rates by decreasing and increasing the scattering phase spaces, respectively. For 4ph phase space, the TD 2FC flattens phonon bands, and allows more redistribution-4ph processes (1 + 2 → 3 + 4) to happen. The combination effect of TD 2FC and TD 4FC reduces 4ph rates of acoustic modes but increases those of optical modes. The TD 3FC and 4FC decrease the phonon scattering cross section and increase the κ ph significantly (by 52% at 3500 K). The contribution from interband (Wigner) phonon conduction is small, even at ultrahigh temperatures. For electronic thermal transport, we find that it is sensitive to and can be changed by 20% by the TD lattice constants. The Lorenz number varies from 1.6 to 3.3 × 10 –8 WΩK –2 at different temperatures. The theoretical prediction in the literature overpredicts κph (e.g., ~28%) and underpredicts the κ el (e.g., ~38%), resulting in an overall underprediction of κ (~26% at 1500 K). The impacts of grain size and defects are found to be strong, leading to the lower observed thermal conductivity in experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Spatial inequities in access to medications for treatment of opioid use disorder highlight scarcity of methadone providers under counterfactual scenarios

Access to treatment and medication for opioid use disorder (MOUD) is essential in reducing opioid use and associated behavioral risks, such as syringe sharing among persons who inject drugs (PWID). Syringe sharing among PWID carries high risk of transmission of serious infections such as hepatitis C and HIV. MOUD resources, such as methadone provider clinics, however, are often unavailable to PWID due to barriers like long travel distance to the nearest methadone provider and the required frequency of clinic visits. The goal of this study is to examine the uncertainty in the effects of travel distance in initiating and continuing methadone treatment and how these interact with different spatial distributions of methadone providers to impact co-injection (syringe sharing) risks. A baseline scenario of spatial access was established using the existing locations of methadone providers in a geographical area of metropolitan Chicago, Illinois, USA. Next, different counterfactual scenarios redistributed the locations of methadone providers in this geographic area according to the densities of both the general adult population and according to the PWID population per zip code. We define different reasonable methadone access assumptions as the combinations of short, medium, and long travel distance preferences combined with three urban/suburban travel distance preference. Our modeling results show that when there is a low travel distance preference for accessing methadone providers, distributing providers near areas that have the greatest need (defined by density of PWID) is best at reducing syringe sharing behaviors. However, this strategy also decreases access across suburban locales, posing even greater difficulty in regions with fewer transit options and providers. As such, without an adequate number of providers to give equitable coverage across the region, spatial distribution cannot be optimized to provide equitable access to all PWID. Our study has important implications for increasing interest in methadone as a resurgent treatment for MOUD in the United States and for guiding policy toward improving access to MOUD among PWID.

59 BASIC BIOLOGICAL SCIENCES↗

Nano-Engineered Interfaces in Dual-Layer Electrodes for Protonic Ceramic Cells with Enhanced Stability and Kinetics

Enhancing interfacial stability and charge transfer in protonic ceramic cells (PCCs) remains a critical challenge, as structural degradation and interfacial resistance often compromise durability and efficiency. Here, we report a nanoengineered dual-layer oxygen electrode architecture designed to address these limitations by introducing a fine-grained nanoparticle interfacial contact layer beneath a porous catalytic backbone. The nanoscale powders, through enhanced sintering activity, densify into a robust interfacial layer that promotes strong chemical bonding, uniform adhesion, and continuous ionic/electronic pathways with the BCZYYb electrolyte. This hierarchical architecture mitigates delamination, redistributes mechanical stress, and establishes efficient charge and mass transport channels without relying on corrosive surface treatments. Electrochemical evaluation demonstrates that the dual-layer design markedly reduces interfacial polarization resistance and accelerates electrode kinetics. Compared to the single-layer counterpart, the architecture achieves a peel strength of 44.53 N/cm 2 , a 40% improvement in peak power density (0.96 W cm –2 at 600 °C), and a 130% enhancement in electrolysis current density (4.78 A cm –2 at 1.57 V). Faradaic efficiency remains as high as 88% under high steam concentrations, underscoring minimal charge loss during practical operation. Notably, the electrode retains stability across 450–600 °C and under transient voltage cycling, with impedance spectra confirming suppressed interfacial resistance growth over prolonged use. These results highlight nanoscale interface engineering as a powerful route to enhance both mechanical robustness and electrochemical kinetics in PCCs. The demonstrated scalability and durability of this architecture provide a versatile platform for advancing solid-state electrochemical systems, including reversible fuel cells and high-efficiency hydrogen production technologies.

Faradaic efficiency↗

Towards a NEAMS-based high-fidelity model of the MARVEL reactor

This report outlines the progress of Idaho National Laboratory in developing a high-fidelity and high-resolution model of the Microreactor Applications Research Validation and Evaluation reactor. The model was developed under the Nuclear Energy Advanced Modeling and Simulation microreactor application driver at Idaho National Laboratory. The overarching objective of this activity is the development of a high-fidelity multiphysics MARVEL model using NEAMS tools, and to verify and validate NEAMS tools against MARVEL reference simulation and experimental data, respectively. This is a unique opportunity to conduct multiphysics analysis on a soon-to-be-deployed microreactor. This multiphysics model developed under the NEAMS-funded INL microreactor application driver leverages three single-physics models coupled via the MOOSE’s MultiApp and Transfer systems. The latter systems enable in-memory data transfer between MOOSE-based and MOOSE-wrapped applications. The first single-physics model, that functions as main application, leverages Griffin to model the neutron transport in the core through the discontinuous finite element (DFEM) discrete ordinates solver (SN). Several optimization flags that were developed by the Griffin developer team were beta-tested to enhance the solver’s performance. These include the combined use of using_average_xs and update_averaged_xs_on that enable to avoid expensive on-the-fly cross sections evaluations at each linear iterations in favor of evaluations of the macroscopic cross sections at each Picard iteration. The second single-physics model uses BISON to handle solid heat transfer and asymptotic hydrogen redistribution analysis in the fuel. While the model returns consistent results for the temperature and hydrogen distribution in the fuel, a mismatch was noticed in the calculated temperature in the reflector due to the value of the gap conductance used in our model. Ongoing investigations are being performed to assess the origin of this discrepancy. Finally, the System Analysis Module (SAM) was used to model the flow of the sodium-potassium eutectic in the primary loop. A first verification was also performed showing good agreement in terms of mass flow rate and inlet temperature. All mesh files were generated using the MOOSE Reactor module, removing the need for external meshing tools. Notably, this workscope represents one of the initial applications of the MOOSE Reactor module for modeling highly irregular geometries. The use of the reactor module significantly streamlined the mesh generation process. The full multiphysics mode, that combines all the single physics models, was leveraged to conduct initial steady-state multiphysics simulations to compute power, and temperature distribution in the reactor. Initial testing was performed for transient simulations as well. In this case, the new checkpoint restart capability for eigenvalue calculations was tested showing the capability for streamlined restart of transient calculations. Future work will focus on improving the fidelity of the model by performing comprehensive code-to-code comparisons. For instance, the full-core Griffin neutronics model will be benchmarked against MCNP reference results, that were provided by the MARVEL design team. Additionally, the SAM T/H model will be verified against reference RELAP-5 results for selected accident scenarios. Besides code-to-code verification exercises, the model fidelity will be improved by replacing the single-channel SAM model with a more complex SAM-Pronghorn coupled model, in which the sub-channel capability is deployed to obtain radial temperature resolution in the coolant. This model will be developed in synergy with the NEAMS thermal hydraulics team.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

CNRS Molten Salt Reactor Benchmark Analysis Using Griffin-Pronghorn Coupled Multi-physics Code System

Molten salt reactors (MSRs) with flowing fuel have the unique feature of utilizing the fuel salt for heat generation and extraction at the same time, since the fuel salt is circulating through the whole primary loop. This movement of fuel salt results in a partial decay of the delayed neutron precursors (DNPs) outside the core and the corresponding redistribution of the DNPs in the active region of the core. To capture this effect, neutronics and thermal-hydraulics (T/H) codes need to be able to handle the movement of DNPs and their decay [1, 2]. Idaho National Laboratory (INL) is actively working on developing the neutronics and T/H codes to model the MSRs with flowing fuel. As one of the efforts, the neutronics code Griffin [3] and the T/H code Pronghorn [4] that are built upon the MOOSE framework [5] have been extended to handle the flowing fuel involving the drift of DNPs, and this capability is being tested and verified with MSR benchmark problems including the CNRS benchmark problem [6]. Additionally, a multi-physics analysis of the molten salt reactor experiment (MSRE) is being conducted to verify and validate the codes against available experimental data [7]. This work presents the analysis of the CNRS benchmark problem using the Griffin-Pronghorn coupled multi-physics code system of INL. Both steady-state and transient problems are analyzed, and the results are compared with those of other participants presented in [6]. The following section provides a description of the benchmark, along with input parameters and observables for each step. Then, the results for each step are presented in comparison with other participants’ results, followed by summary and conclusion.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Regioisomeric Engineering for Multicharge and Spin Stabilization in Two-Electron Organic Catholytes

Developing multicharge and spin stabilization strategies is fundamental to enhancing the lifetime of functional organic materials, particularly for long-term energy storage in multiredox organic redox flow batteries. Current approaches are limited to the incorporation of electronic substituents to increase or decrease the overall electron density or bulky substituents to sterically shield reactive sites. With the aim to further expand the molecular toolbox for charge and spin stabilization, we introduce regioisomerism as a scaffold-diversifying design element that considers the collective and cumulative electronic and steric contributions from all of the substituents based on their relative regioisomeric arrangements. Through a systematic study of regioisomers of near-planar aromatic cyclic triindoles and nonplanar nonaromatic cyclic tetraindoles, we demonstrate that this regioisomeric engineering strategy significantly enhances the H-cell cycling stability in the above two new classes of 2e – catholytes, even when current strategies failed to stabilize the multicharged species. Density functional theory calculations reveal that the strategy operates by redistributing the charge and spin densities while highlighting the role of aromaticity in charge stabilization. The most stable 2e – catholyte candidate was paired with a viologen derivative anolyte to achieve a proof-of-concept all-organic flow battery with 1.26–1.49 V, 98% capacity retention, and only 0.0117% fade/h and 0.00563% fade/cycle over 400 cycles (192 h), which is the highest capacity retention ever reported over 400 cycles in a multielectron all-organic flow battery setup. We anticipate regioisomeric engineering to be a promising strategy complementary to conventional electronic and steric approaches for multicharge and spin stabilization in other functional organic materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗