Search NASA⌕ Search

SEARCH · Search NASA

Results for “NEUTRON ACTIVATION”

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 235 records · Page 13

GRB 180128A: A second magnetar giant flare candidate from the Sculptor Galaxy

Magnetars are slowly rotating neutron stars that possess the strongest magnetic fields known in the cosmos (10 14 − 10 15 G). They display a range of transient high-energy electromagnetic activity. The brightest and most energetic of these events are the gamma-ray bursts (GRBs) known as magnetar giant flares (MGFs), with isotropic energies E iso ≈ 10 44 − 10 46 erg. Only seven MGF detections have been made to date: three unambiguous events occurred in our Galaxy and the Magellanic Clouds, and the other four MGF candidates are associated with nearby star-forming galaxies. As all seven identified MGFs are bright at Earth, additional weaker events likely remain unidentified in archival data. We conducted a search of the Fermi Gamma-ray Burst Monitor database for candidate extragalactic MGFs and, when possible, collected localization data from the Interplanetary Network (IPN) satellites. Our search yielded one convincing event, GRB 180128A. IPN localizes this burst within NGC 253, commonly known as the Sculptor Galaxy. The event is the second MGF in modern astronomy to be associated with this galaxy and the first time two bursts have been associated with a single galaxy outside our own. Here we detail the archival search criteria that uncovered this event and its spectral and temporal properties, which are consistent with expectations for a MGF. We also discuss the theoretical implications and finer burst structures resolved from various binning methods. Our analysis provides observational evidence of an eighth identified MGF.

Trigg, Aaron C. (ORCID:000900068598728X)↗

Reaching For New Physics With MeV-scale Reconstruction In The MicroBooNE LArTPC Neutrino Detector

Large neutrino liquid argon time projection chamber (LArTPC) experiments can broaden their physics reach by reconstructing MeV-Scale energy depositions, or blips, in their data. We demonstrate new calorimetric and particle discrimination capabilities at the MeV scale using reconstructed blips in MicroBooNE LArTPC data at Fermilab. A concentration of low-energy ($<$3 MeV) blips is observed around fiberglass mechanical support struts along the TPC edges, with spectral features consistent with the Compton edge of the 2.614 MeV $^{208}$Tl decay $\gamma$ ray. With these features we perform the electron energy scale calibration to few-percent precision and yield the specific activity of $^{208}$Tl in the struts, $(11.7 \pm 0.2 \text{(stat)} \pm 2.8 \text{(syst)})$ Bq/kg. Using cosmogenic blips above 3 MeV, we demonstrate the ability of large LArTPCs to discriminate low-energy proton and electron depositions. An enriched low-energy proton sample selected with this technique is smaller in data than in dedicated CORSIKA simulations, pointing to possible mismodeling in CORSIKA incident cosmic fluxes or Geant4 particle transport. These methods are applied to MicroBooNE's inclusive single-photon search, which reported a 2.2$\sigma$ excess below 600 MeV in shower energy for events with no reconstructed protons. By identifying and classifying blips near single-photon events selected by the WireCell reconstruction framework, a more comprehensive labeling of nearby hadronic activity is established: blips upstream of the shower axis indicate previously unidentified final-state protons, while elevated blip counts at wide angles signal final-state neutrons. Taken together with MiniBooNE's long-standing low-energy excess (LEE) and MicroBooNE electron-like and sterile neutrino searches disfavored as possible explanations of the MiniBooNE anomaly, this analysis motivates an expanded exploration of the single-photon channel in Fermilab's short-baseline LArTPC program. This thesis documents the current status of this enhanced analysis, which will form a key part of MicroBooNE's final low-energy-excess results.

Andrade Aldana, Diego Armando [IIT, Chicago (main)↗

Superionic-like diffusion in yttrium dihydride

For the next-generation high temperature microreactors, yttrium dihydride (YH 2 ) is an attractive solid state neutron moderator. Despite a number of recent investigations, the mechanism of hydrogen transport remains poorly understood. Experimental evaluations of diffusivity are inconclusive with large variations in diffusivities and activation energies. In this work, we perform ab initio molecular dynamics (AIMD) simulations on YH 2 for temperatures spanning 300 K to 1200 K. Our main finding is that YH 2 shows a superionic-like behavior with hydrogen atoms hopping from one native site to another above a characteristic temperature of 800 K. This correlated motion results in quasi-one-dimensional string-like displacements that enable the hydrogen atoms to diffuse rapidly. We confirm that the octahedral sites are mostly unoccupied, although channeling through them is the most favored pathway between lattice hops above 800 K. At the highest temperature of 1200 K, the string relaxation time is merely of the order of a few picoseconds, which indicates a liquid-like diffusive behavior. Based on the formation of spontaneous thermal vacancies, an order-disorder crossover temperature T α ~ 800 K is established for YH 2 with an activation energy of 0.83 eV for hydrogen diffusion in the superionic-like state.

Superionic-like Diffusion↗

Role of Nuclear Science User Facilities (NSUF) in Nuclear Energy Materials Research

The Nuclear Science User Facilities (NSUF) is one of a diverse number of U.S. Department of Energy (DOE) user facilities established to provide researchers with the most advanced tools of modern science. The NSUF is unique and represents a consortium of capabilities distributed across the U.S. at twenty-one institutions. The NSUF is centered at and managed from the Idaho National Laboratory (INL), where it was originally founded, but it coordinates activities at twenty “partner” institutions that include universities, the Center for Advanced Energy Studies (CAES), national laboratories, and a nuclear industry vendor. These institutions have capabilities that include neutron, ion, and gamma irradiation, hot cells, advanced materials characterization equipment, and high-performance computing resources. Many of these capabilities were beyond reach for most researchers before NSUF. The NSUF provides researchers to access these capabilities at no cost to nuclear energy researchers to produce the highest quality research results to increase understanding of advanced nuclear energy technologies important to DOE-NE and support national priorities by adapting to the needs of DOE-NE programs, industry, and new innovative concepts for sustainable nuclear future.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Measured Reaction Rate Ratios of 235 U, 238 U, 237 Np, 239 Pu Samples in the PFUNS 235 U Prompt Fission Neutron Spectrum Criticality Experiment

The Prompt Fission Uranium Neutron Spectrum experiment, an experiment to reduce uncertainties in the high energy tail of the 235 U prompt fission spectrum, achieved success by performing two separate irradiations measuring approximately 40 different IRDFF-II reactions total using over 20 different foil materials at the National Criticality Experiments Research Center in February 2024. The criticality experiment utilized a set of highly enriched uranium hemispherical shells of increasing diameters with a large void in the center where the samples were located. The focus of this work is the first of two PFUNS irradiations focused on irradiating two of each fission foils, one bare and one cadmium covered, along with metal activation foils containing reaction products with short half-lives, such as indium, iron, and gold along with nickel fluence monitors. This work focuses on presenting the initial reaction rate ratio results of the fission foils from the aforementioned first irradiation to assist in nuclear data validation of those species in a nearly pure 235 U prompt fission neutron spectrum and compares to the Lady Godiva and Flattop-25 historic experiments at the Los Alamos Critical Experiments Facility. Future work will combine fission foil and metallic activation foil reaction rate ratio results from both the first and second higher power irradiation and reaction rates for a final spectral adjustment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

New constraint on the Np 237 ( n , γ ) Np 238 integral cross section using the Godiva-IV critical assembly

Accurate knowledge of the 237 Np(n, γ) 238 Np cross section at fast neutron energies is important for applied nuclear science. The presently available experimental data has large disagreements in the fast neutron region. Perform a model-independent measurement of the 237 Np(n, γ) 238 Np integral cross section using a well characterized fast neutron source and compare the result with previous measurements and current nuclear data evaluations. Provide an integral measurement that can be used as a benchmark for current evaluations. Multiple samples of 237 Np were irradiated in the Godiva-IV critical assembly. Following the irradiation, the samples placed in a γ-ray counting setup and the γ-rays emitted from the decay of 238 Np were measured over a time period of approximately 7 days. Multiple γ-ray decay branches of 238 Np were observed. The observed activity of 238 Np was used to calculate the amount of 238 Np produced during the irradiation via the 237 Np(n, γ) 238 Np reaction and an integral cross section of 342(11) mb was measured for the Godiva-IV neutron spectrum. Further, the 238 Np half-life has been measured with a result of 50.31(5) hours. The 237 Np(n, γ) 238 Np integral cross section measured in this work is in agreement with overlapping 1σ error bands to ENDF/B-VIII.0. However, the measured value is 3σ away from the calculated integral cross section using JENDL-5. This measurement offers a reliable benchmark for future 237 Np(n, γ) 238 Np cross section evaluations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Li 21 Ge 8 P 3 S 34 : New Lithium Superionic Conductor with Unprecedented Structural Type

Abstract Lithium superionic conductors are pivotal for enabling all‐solid‐state batteries, which aim to replace liquid electrolytes and enhance safety. Herein, we report the discovery of an unprecedented lithium superionic conductor, Li 21 Ge 8 P 3 S 34 , featuring a novel structural type and a new composition in the Li–Ge–P–S system. This material exhibits high lithium ionic conductivity of approximately 1.0 mS cm −1 at 303 K with a low activation energy of 0.20(1) eV. It's unique crystal structure was elucidated using three‐dimensional electron diffraction (3D ED) and further refined through combined powder X‐ray and neutron diffraction analyses. The structure consists of alternating two‐dimensional slabs: one of corner‐sharing GeS 4 tetrahedra and the other of isolated PS 4 tetrahedra, enabling efficient lithium‐ion transport through a tetrahedrally interconnected network of 1D, 2D, and 3D diffusion pathways. This distinctive structural motif provides a novel design strategy for next‐generation solid electrolytes, broadening the structural landscape of lithium superionic conductors. With further advancements in compositional tuning and interfacial engineering, Li 21 Ge 8 P 3 S 34 could contribute to the development of high‐performance all‐solid‐state batteries.

Chemistry↗

Operando neutron imaging-guided gradient design of Li-ion solid conductor for high-mass-loading cathodes

High-mass-loading cathodes are crucial for achieving high energy density in all-solid-state batteries from the lab scale to industry. However, as mass-loading increases, electrochemical performance is significantly compromised due to sluggish kinetics. In this work, operando neutron imaging is deployed on a high-mass-loading NMC 811 cathode of 33 mg/cm 2 (5.0 mAh/cm 2 ) and directly visualizes the lithiation prioritization of the cathode active material (CAM) from the solid electrolyte membrane side to the current collector side. In addition to the tortuosity, another key limitation on ion transfer in the cathode arises from the mismatch between the uniform distribution of the solid electrolyte (catholyte) in the conventional composite cathode and the non-uniform Li + flux generated by the faradaic reaction of CAMs. Therefore, we engineer a gradient in the catholyte concentration to match the Li + flux distribution as a means of eliminating the ion transfer obstacle. This approach demonstrates enhanced rate performance, even with high-mass-loading cathodes. A LiCoO 2 composite cathode with 100 mg/cm 2 high-mass-loading exhibits an areal capacity of 10.4 mAh/cm 2 at a current density of 2.25 mA/cm 2 . This work provides insight into the ion-transport limitation in thick cathodes and demonstrates an effective gradient design to overcome the kinetic barrier and achieve high battery performance.

Batteries↗

Computationally designed coiled coil ‘bundlemers’ as model colloidal nanoparticles for solution assembly and materials design (Final Report)

As a collaborative team at the University of Delaware and the University of Pennsylvania, Kloxin, Pochan and Saven designed new biomimetic nanomaterials de novo, leveraging a variety of complementary areas of expertise: computational design of biopolymers (Saven at the University of Pennsylvania), and synthesis and characterization (Kloxin and Pochan at the University of Delaware). Overall activities included: sequence-specific peptide synthesis; covalent crosslinking; noncovalent assembly; site-specific functionalization; and nanostructural characterization using electron microscopy and solution-phase (x-ray and neutron) scattering. Using natural and non-natural amino acids, the team created modular, functional peptide building blocks for elaboration of new nanostructured materials. Ultimately, the development of robust peptide-based, building blocks provides tools for researchers to readily produce complex nanomaterial structures in a wide range of applications. The project had three, interconnecting goals in an effort to provide the broader scientific community with a new peptide-based paradigm for materials design and characterization. First, we further developed the coiled-coil bundle-based toolbox (otherwise known as the ‘bundlemer’ toolbox) via computational design with experimental bundle assembly verification. Second, we developed new uses of covalent interactions, in addition to desired physical (noncovalent) interactions, to assemble bundlemers into 1-D polymer chains with targeted chain rigidity, length, and dispersity. Thirds, we used the above designs to experimentally realize (physical or covalent) polymers to target the creation of liquid crystals or to realize interparticle assembly into nanoporous lattices. The close integration of the three groups was instrumental in success of the biomolecular materials design, formation, and understanding for future designs.

36 MATERIALS SCIENCE↗

Direct cross-section measurement of the weak r -process Sr 88 ( α , n ) Zr 91 reaction in ν -driven winds of core-collapse supernovae

About half of the heavy elements beyond iron are known to be produced by the rapid neutron capture process, known as the r process. However, the astrophysical site producing the r process is still uncertain. Chemical abundances observed in several cosmic sites indicate that different mechanisms should be at play. For instance, the abundances around silver measured in a subset of metal-poor stars indicate the presence of a weak r process. This process may be active in neutrino-driven winds of core collapse supernovae where (α,n) reactions dominate the synthesis of Ζ ≈ 40 elements in the expelled materials. Scarcely measured, the rates of (α,n) reactions are determined from statistical Hauser-Feshbach calculations with α-optical-model potentials, which are still poorly constrained. Further, the uncertainties of the (α,n) reaction rates therefore make a significant contribution to the uncertainties of the abundances determined from stellar modeling. In this work, the 88 Sr (α,n)⁢ 91 Zr reaction which impacts the weak r-process abundances has been probed at astrophysics energy for the first time; directly measuring the total cross sections at astrophysical energies of 8.37–13.09 MeV in the center of mass (3.8–7.5 GK). Two measurements were performed at ATLAS with the electrically segmented ionization chamber MUSIC, in inverse kinematics, while following the active target technique. The cross sections of this α-induced reaction on 88 Sr, located at the shell closure N = 50, have been found to be lower than expected, by a factor of 3, despite recent statistical calculations validated by measurements on neighboring nuclei. This result encourages more experimental investigations of (α,n) reactions, at Ν = 50 and towards the neutron-rich side, to further test the predictive power and reliability of such calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Dissecting Disorder: Defect-Driven Structural Complexity in Layered Li3InCl6 Solid Electrolyte

Halide solid electrolytes have emerged as promising candidates for solid-state batteries owing to their high oxidative stability and ionic conductivity. Among them, Li3InCl6 (LIC) has attracted significant attention. However, diffraction patterns of LIC synthesized via different methods exhibit distinct differences particularly at low-angle reflectionsindicative of underlying structural disorder. These variations are attributed to deviations from ideal crystallographic order, especially stacking faults, whose impact on structure and ion transport remains poorly understood. Here, we identify and quantify stacking faults in LIC samples prepared under different synthetic conditions. Using X-ray diffraction and time-of-flight neutron diffraction, we construct and refine stacking fault models that accurately reproduce the experimental diffraction features. LIC samples with higher degrees of stacking faults exhibit only negligible differences in ionic conductivities and activation energies. This indicates that stacking faults have a limited impact on altering the Li+ diffusion pathway along the c-axis, likely due to the high concentration of vacancies in the In layers, while Li+ diffusion remains nearly unchanged in the ab-plane. Our results account for the observed differences in diffraction patterns across samples and provide a quantitative assessment of faulting probabilities and stacking sequences. The insights gained from this study are expected to be broadly applicable to other layered halide solid electrolytes and contribute to a deeper understanding of the role of structural disorder in ion transport.

Liu, Jue [ORNL] (ORCID:000000024453910X)↗

Summary Report of the FY24 DOE Contributions to the GIF VHTR CMVB

The Generation-IV Forum (GIF) Very-High-Temperature Reactor-Computational Methods Validation and Benchmark (VHTR-CMVB) initiative, involving organizations from Korea Atomic Energy Research Institute (KAERI) (South Korea), Institute of Nuclear and New Energy Technology of Tsinghua University (INET) (China), U.S. Department of Energy (DOE) (U.S.), Joint Research Centre (JRC) (Europe), and Japan Atomic Energy Agency (JAEA) (Japan), is dedicated to the verification and validation of tools for High-Temperature Gas-Cooled Reactors (HTGRs) analysis, using data shared by Computational Methods Validation and Benchmark (CMVB) signatories. For FY24, the US DOE CMVB has committed to several critical activities. Under WP1, led by the US, the integration of the High Temperature Gas Cooled Reactor - Pebble-Bed Module (HTR-PM) Phenomena Identification and Ranking Table (PIRT) into the comparative PIRT is progressing, with a new draft of the comparison tables issued earlier this year and currently being utilized by INET for their contribution. Neutronic validation efforts under WP3 include the preparation of the burnup analysis benchmark, preliminary calculations, and the development of reference models and results. In WP2, a validation exercise for hot gas mixing in the lower plenum of HTR-PM is in progress, using experimental data from INET (China) to validate modeling approaches. A model of the experimental facility has been developed using StarCCM+, with initial calculations slated for presentation at the GIF CMVB meeting this fall. Another WP2 activity focuses on validating numerical models for air-cooled Reactor Cavity Cooling System (RCCS) with experimental data from the Wisconsin Madison RCCS facility. A high-fidelity model, developed using NEK-RS, is currently being validated with available data from a low power forced convection test. These efforts are aimed at enhancing and confirming the accuracy of HTGR analysis tools, ensuring their alignment with experimental data and regulatory requirements.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Melatonin-Induced Modulation of Cholesterol-Enriched Model Neuronal Membranes

Melatonin, a hormone primarily produced by the brain’s pineal gland, not only regulates circadian rhythms, but also influences the structural and biophysical properties of neuronal membranes. Its amphiphilic nature enables direct incorporation into lipid bilayers and preferential interactions with cholesterol-rich lipid rafts, critical hubs for cellular signaling and membrane organization. Despite increasing recognition of its membrane activity, the molecular basis of melatonin’s interactions with coexisting liquid-ordered (L o ) and liquid-disordered (L d ) phases remains unclear. Here, in this study, we combine small-angle neutron scattering (SANS) and all-atom molecular dynamics simulations to examine model neuronal membranes composed of DSPC, DOPC, POPC, and cholesterol. Our results show that melatonin preserves domain morphology while adopting distinct orientations within the bilayer and at the membrane interface, allowing both lateral and transmembrane bridging across lipid phases. These findings establish the molecular underpinnings of melatonin’s modulation of membrane heterogeneity and provide strong support for its receptor-independent actions.

atomistic simulations↗

Neutron capture measurements at UMass lowell research reactor

A new facility was designed around the thermal column beam port of the 1 MW Research Reactor at the University of Massachusetts Lowell. Thermal neutrons are collimated to a 1-inch diameter beam and incident on samples to induce the radiative neutron capture. New measurements of capture gamma rays are planned for Mn, Cu, Ni, Cr, and Gd samples in the next few years. The experiments will be carried out in close collaboration with the ENSDF (Evaluated Nuclear Structure Data File) evaluation group at Brookhaven National Laboratory. The gamma rays are measured using an array of high-resolution HPGe detectors. Some of the HPGe detectors use active Compton shields that were designed from scintillation detectors and improve significantly the signal-to-background in the measured gamma-ray spectra with HPGe. The experimental spectra will be validated by the GEANT4 simulations of the array and theoretical models of the emission of gamma rays from the compound nucleus. In conclusion, the new experimental results on the capture gamma ray intensities will be incorporated in future ENSDF evaluations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Determination of Site Occupancy in the M–Pd–Zn (M = Cu, Ag, and Au) γ-Brass Phase by CALculation of PHAse Diagrams Modeling and Rietveld Refinement

The Pd–Zn γ-brass phase provides exciting opportunities for synthesizing site-isolated catalysts with precisely controlled Pd active site ensembles. Introducing a third metallic element into the γ-brass lattice further perturbs the catalytic active site ensembles. Here, in this work, we introduce coinage metallic elements M (M = Cu, Ag, and Au) into the Pd–Zn γ-brass phase and investigate the site occupation factors of each element in the γ-brass lattice. The CALculation of PHAse Diagrams (CALPHAD) modeling approach supported by energetics predicted by the density functional theory and X-ray and neutron diffraction with Rietveld refinement were used to identify the SOF on each Wyckoff site for various M amounts alloyed into the Pd–Zn γ-brass phase. The present analysis unveils the strong preference for Pd occupying the outer tetrahedral (OT) site in the γ-brass lattice, while the coinage metallic elements tend to substitute for Zn on the octahedral (OH) site. The determination of site occupancy in the bulk M–Pd–Zn γ-brass phase provides opportunities to investigate and tailor potential catalytically active site ensembles in the γ-brass phase materials.

36 MATERIALS SCIENCE↗

Pervasive symmetry-lowering nanoscale structural fluctuations in the cuprate La 2−𝑥 ⁢Sr 𝑥 ⁢Cu⁢O 4

The cuprate superconductors are among the most widely studied quantum materials, yet there remain fundamental open questions regarding their electronic properties and the role of structural degrees of freedom. Recent neutron and x-ray scattering measurements uncovered exponential scaling with temperature of the strength of orthorhombic fluctuations in the tetragonal phase of La 2−𝑥⁢ Sr 𝑥 ⁢Cu⁢O 4 and Tl 2 ⁢Ba 2 ⁢Cu⁢O 6+𝛿 , unusual behavior that closely resembles prior results for the emergence of superconducting fluctuations, and that points to a common origin rooted in inherent correlated structural inhomogeneity. Here, in this study, we use neutron and x-ray diffuse scattering to provide further insight into structural fluctuations in La 2−𝑥 ⁢Sr 𝑥 ⁢Cu⁢O 4 . We perform measurements up to temperatures that approach the melting point of the undoped parent compound, and we investigate the effects of in situ in-plane uniaxial elastic stress and plastic deformation. Our neutron scattering results for the parent compound La 2 ⁢Cu⁢O 4 reveal that short-range orthorhombic fluctuations persist to the maximum experimental temperature of nearly 1000 K, i.e., to a significant fraction of the crystallization temperature. At this temperature, the spatial characteristic length extracted from the momentum-space data is still about three lattice constants. The neutron scattering experiment enables quasistatic discrimination and reveals that the response is increasingly dynamic at higher temperatures. We also find that neither compressive elastic stress along the tetragonal [110] direction nor irreversible plastic deformation activating the ⟨100⟩⁢[001] slip system significantly alter the robust orthorhombic fluctuations. Overall, these results support the notion that these fluctuations couple to subtle, underlying inhomogeneity that underpins the cuprate phase diagram. Finally, in La 1.8⁢ Sr 0.2⁢ Cu⁢O 4 , we uncover low-energy structural fluctuations at a nominally forbidden reflection that are distinct from the orthorhombic distortions. While the origin of these fluctuations is not clear, they might be related to the presence of extended defects such as dislocations or stacking faults.

Spieker, R. J. [Univ. of Minnesota, Minneapolis, M↗

Alkali Metal Cation Effects for Rapid C–H Activation by Iron(0) Complexes

C–H oxidative addition is a key reaction in organometallic catalysis, motivating efforts to accelerate it. Here, we examine an anionic beta-diketiminate-supported iron(0) species that was previously observed to activate C-H bonds with Na(15-crown-5), but not with K(18-crown-6) or Rb(18-crown-6). Though crown ethers are usually seen as beneficial due to their ability to solubilize alkali metal cations, we observe that removing the crown ether leads to rapid and complete oxidative addition of the C-H bond even by K, Rb, and Cs. The products are iron(II) phenyl hydride complexes that exist as dimers bridged by the alkali metals. Neutron crystallography of the cesium complex verifies the presence and location of the bridging hydrides. It is likely that the crown-free alkali metal cations have greater Lewis acidity that enables them to facilitate oxidative addition of the C-H bond. This system gives insight on how to control the rate and favorability of C-H activation through manipulation of the countercation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Initial Development of Fusion Magnet Simulation Capabilities for Performance and Safety Evaluation Using the MOOSE Framework

Fusion energy holds the promise of being a transformative technology as a carbon-neutral, sustainable source of energy. Whole device modeling and the development of fusion digital twins will be increasingly important for emerging fusion device concepts at both national laboratories and within the commercial fusion industry. However, meeting the challenge of whole device modeling of fusion energy devices requires robust, multiphysics, multiscale modeling and simulation technologies capable of running on large-scale supercomputers. Detailed analysis of individual systems at-scale is also required to ensure safe and efficient operation as well as provide the safety basis for future device designs and licensing activities. In a tokamak, toroidal and poloidal magnets confine and shape the fusion plasma to promote the fusion reaction. High plasma temperatures and high magnetic field requirements in modern design concepts (leading to high amounts of energy stored within each magnet) impose electrical, thermal, and mechanical loads on the magnet components, which in turn impacts the safety considerations of the magnet and their supporting systems. Idaho National Laboratory (INL) has a history of working in this space, including development and benchmarking of the Magnetic System Circuitry Analysis Program (MSCAP) and Magnet Arcing (MAGARC) codes to study magnet quench events; notably, MAGARC was used to study quenching during the ITER Engineering Design Activity. However, these legacy codes and capabilities are not parallel and scalable, and new tools are required for future advances in this area, which leads to the INL-developed Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Developed originally for fission reactor systems under United States Department of Energy, Office of Nuclear Energy modeling and simulation programs, the MOOSE framework has been well-suited to multiscale, multiphysics modeling and simulation needs for nuclear systems. The framework is open-source, well-tested, under continuous development and deployment, and developed to a Nuclear Quality Assurance, Level 1 software quality standard. MOOSE has also been used in the fusion space previously in several projects: INL’s Tritium Migration Analysis Program, Version 8 (TMAP8) for tritium migration, UK Atomic Energy Authority’s A Unified Resource for OpenMC (fusion) Reactor Applications (AURORA) code for fusion thermo-mechanical and neutronics analysis, and Argonne National Laboratory’s Cardinal for high-fidelity computational fluid dynamics and neutronics. However, to model superconducting magnets, several MOOSE enhancements are required: additions to the current MOOSE electromagnetic capabilities, new material libraries for superconductors of interest (such as YBCO), as well as fusion-specific models for thermo-mechanics. This talk will discuss initial development activities to build these capabilities in MOOSE, focusing on initial validation and benchmarking activities. Proposed coupling workflows and future work to support the simulation of fusion magnets and magnet structural assemblies for performance and safety evaluation in MOOSE will also be discussed.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗