Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective core potentials”

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 451 records · Page 25

In Situ Oxidation Studies of High-Entropy Alloy Nanoparticles

Although high-entropy alloys (HEAs) have shown tremendous potential for elevated temperature, anticorrosion, and catalysis applications, little is known on how HEA materials behave under complex service environments. Herein, we studied the high-temperature oxidation behavior of Fe 0.28 Co 0.21 Ni 0.20 Cu 0.08 Pt 0.23 HEA nanoparticles (NPs) in an atmospheric pressure dry air environment by in situ gas-cell transmission electron microscopy. It is found that the oxidation of HEA NPs is governed by Kirkendall effects with logarithmic oxidation rates rather than parabolic as predicted by Wagner's theory. Further, the HEA NPs are found to oxidize at a significantly slower rate compared to monometallic NPs. The outward diffusion of transition metals and formation of disordered oxide layer are observed in real time and confirmed through analytical energy dispersive spectroscopy, and electron energy loss spectroscopy characterizations. Localized ordered lattices are identified in the oxide, suggesting the formation of Fe 2 O 3 , CoO, NiO, and CuO crystallites in an overall disordered matrix. Hybrid Monte Carlo and molecular dynamics simulations based on first-principles energies and forces support these findings and show that the oxidation drives surface segregation of Fe, Co, Ni, and Cu, while Pt stays in the core region. The present work offers key insights into how HEA NPs behave under high-temperature oxidizing environment and sheds light on future design of highly stable alloys under complex service conditions.

36 MATERIALS SCIENCE↗

Strategy to Develop a Control Scheme for Core Thermal Performance Optimization

This report outlines a three-year research plan for optimizing reactor core thermal performance through use of a Digital Twin (DT) model and a set of neutron detectors to correspondingly adapt the reactor control strategy. It identifies reactor features and achievable neutron measurements that factor into this optimization task. It considers spatial effects that are important and how they can be managed by a real-time algorithm that controls reactivity actuators. The report presents a set of tasks, along with corresponding methods, for accomplishing this objective. The final task in this set is to validate the combined optimization procedure in the Purdue University’s PUR-1 reactor, for which we have an agreement of understanding. Additionally, we describe an alternative approach for overcoming some of the limitations inherent in the above approach, which includes the computational costs of high-fidelity simulations and achievable integration of neutron measurements into the DT. The proposed algorithms lower the technology readiness level (TRL) of the overall approach as some additional analysis and the development of a new neutron detector concept are necessary. In particular, a sensor measuring the gradient of the neutron flux is required, and the prototype is expected to be validated in the PUR-1 reactor. The potential benefits and the opportunity to significantly push forward the state-of-the-art make this approach worthy of further exploration in the upcoming years.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Study of alkaline carbonate cooling to mitigate Ex-Vessel molten corium accidents

To mitigate adverse effects from molten corium following a reactor pressure vessel failure (RPVF), some new reactor designs employ a core catcher and a sacrificial material (SM), such as ceramic or concrete, to stabilize the molten corium and avoid containment breach. Existing reactors cannot easily be modified to include these SMs but could be modified to allow injectable cooling materials. Current reactor designs are limited to using water to stabilize the corium, but this can create other issues such as reaction of water with the concrete forming hydrogen gas. Here the novel SM proposed is a granular carbonate mineral that can be used in existing light water reactor plants. The granular carbonate will decompose when exposed to heat, inducing an endothermic reaction to quickly solidify the corium in place and producing a mineral oxide and carbon dioxide. Corium spreading is a complex process strongly influenced by coupled chemical reactions, including decay heat from the corium, phase change, and reactions between the concrete containment and available water. A recently completed Sandia National Laboratories laboratory directed research and development (LDRD) project focused on two research areas: experiments to demonstrate the feasibility of the novel SM concept, and modeling activities to determine the potential applications of the concept to actual nuclear plants. Small-scale experiments using lead oxide (PbO) as a surrogate for molten corium demonstrate that the reaction of the SM with molten PbO results in a fast solidification of the melt due to the endothermic carbonate decomposition reaction and the formation of open pore structures in the solidified PbO from CO 2 released during the decomposition. A simplified carbonate decomposition model was developed to predict thermal decomposition of carbonate mineral in contact with corium. This model was incorporated into MELCOR, a severe accident nuclear reactor code. A full-plant MELCOR simulation suggests that by the introduction of SM to the reactor cavity prior to RPVF ex-vessel accident progression, e.g., core-concrete interaction and core spreading on the containment floor, could be delayed by at least 15 h; this may be enough for additional accident management to be implemented to alleviate the situation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

New proposed ITER divertor design using carbon insert on tungsten to mitigate ELMs and secondary radiation effects on nearby components

Abstract Building a successful device for the magnetic fusion energy production is a great challenge. ITER is an international project of the tokamak based magnetic fusion design being developed for the demonstration of the feasibility of thermonuclear technologies for future realization of successful commercial fusion energy. A key obstacle to a successful magnetic fusion energy production is however, the performance during abnormal events including plasma disruptions and edge-localized modes (ELMs). A credible reactor design must tolerate at least a few of these transient events without serious consequences such as melting of the structure. This paper investigates and compares the performance of the current ITER tokamak design during two types of transient events, i.e., ELMs occurring at normal operation and disruptions during abnormal operation. We simulated the divertor components response using our integrated 3D HEIGHTS package. The simulations include self-consistent modeling of the interaction of the released core plasma particles with the initial solid divertor material, energy deposition processes, vaporization of divertor material, secondary plasma formation and MHD evolution, incident core particles collisions and scattering from this dense secondary plasma, photon radiation of secondary plasma, and the resulting heat loads on nearby components. Our simulations showed that using a small carbon insert around the strike point can significantly reduce the overall expected damage on the tungsten dome structure, reflector plates, and prevent tungsten vaporization and its potential core plasma contamination.

47 OTHER INSTRUMENTATION↗

Functional Data Analysis for Extracting the Intrinsic Dimensionality of Spectra: Application to Chemical Homogeneity in the Open Cluster M67

High-resolution spectroscopic surveys of the Milky Way have entered the Big Data regime and have opened avenues for solving outstanding questions in Galactic archeology. However, exploiting their full potential is limited by complex systematics, whose characterization has not received much attention in modern spectroscopic analyses. In this work, we present a novel method to disentangle the component of spectral data space intrinsic to the stars from that due to systematics. Using functional principal component analysis on a sample of 18,933 giant spectra from APOGEE, we find that the intrinsic structure above the level of observational uncertainties requires ≈10 functional principal components (FPCs). Our FPCs can reduce the dimensionality of spectra, remove systematics, and impute masked wavelengths, thereby enabling accurate studies of stellar populations. To demonstrate the applicability of our FPCs, we use them to infer stellar parameters and abundances of 28 giants in the open cluster M67. We employ Sequential Neural Likelihood, a simulation-based Bayesian inference method that learns likelihood functions using neural density estimators, to incorporate non-Gaussian effects in spectral likelihoods. By hierarchically combining the inferred abundances, we limit the spread of the following elements in M67: Fe ≲ 0.02 dex; C ≲ 0.03 dex; O, Mg, Si, Ni ≲ 0.04 dex; Ca ≲ 0.05 dex; N, Al ≲ 0.07 dex (at 68% confidence). Our constraints suggest a lack of self-pollution by core-collapse supernovae in M67, which has promising implications for the future of chemical tagging to understand the star formation history and dynamical evolution of the Milky Way.

79 ASTRONOMY AND ASTROPHYSICS↗

Emulation and Adversarial Analysis of EV Charging Networks

In the effort of decarbonization and evolution of the modern electrical grid, electric vehicles (EVs) play a key part to transform the grid. However, due to the rapid adoption of EVs and the demand of the charging infrastructure required to power said EVs, risk of a cyber-attack may impose serious consequences. There is a need to analyze and protect the charging ecosystem infrastructure from cyber threats before it reaches wide-scale deployment. In an effort to secure vehicle to grid (V2G) communications, standardization is necessary for continued reliable system operation. The protocol ISO 15118 outlines controls and practices that should be implemented for secure vehicle to grid (V2G) communications. The standard is gaining momentum for American markets as the demand for EV infrastructure grows. The adoption of ISO 15118 in American markets poses several challenges: the deployment of a public key infrastructure (PKI) as outlined within the standard, interoperability of charging different EVs with chargers from different manufactures using the PKI, and scaling the ecosystem to meet the demand while managing risks. This project was created to understand potential cyber and scaling challenges of PKI for EV infrastructure through utilizing a series of emulated components mapping to what exists in the EV ecosystem today, and the components of the PKI that are under development. The key nodes within the emulation that are under development are: electric vehicle (EV), electric vehicle supply equipment (EVSE), charge network operator (CNO), certificate authority (CA), and online certificate status protocol (OCSP) that must all interact using secure and trusted communications. With these emulated components and utilizing orchestration methods to rapidly deploy and scale the components, the ability to analyze risks of the ecosystem and address gaps before the PKI ecosystem is fully deployed to production should yield a more robust and mature production charging infrastructure. Our approach will use a modular architecture of virtual machines within an orchestration platform and will target scales of 100s, 1000s, and 10,000s of entities interacting. The core research questions trying to be answered with this scope of work are: what are the impacts of a rogue CA, what are the risks of certificate revocation list (CRL) management, what is the value of OCSP stapling, what components are vulnerable to DOS attacks, and what test effective payloads may impact the components.

charging ecosystem↗

Feasibility of raised inner strike point equilibria scenario in ITER for detritiation from beryllium co-deposits

Abstract In ITER, tritium retention primarily occurs through co-deposition with beryllium. To avoid exceeding the strict tritium inventory limit, efficient tritium recovery techniques are essential. Baking is the ITER baseline for tritium recovery, but its effectiveness in removing tritium from thick beryllium layers is limited. A raised strike point scenario is considered an alternative method for removing tritium from the ITER inner vertical divertor target by heating components via plasma flux. This paper presents SOLPS-ITER code simulations conducted under various conditions, assessing the divertor performance and tritium outgassing of the raised strike point scenario. As the strike point is raised, recycled neutrals are not efficiently baffled by the dome and scrape-off layer, significantly changing the neutral trajectory and ionization source distribution. This improves detachment accessibility but worsens core-edge compatibility compared to the baseline scenario. However, in the partially detached condition, the impact of raising the strike point, perpendicular transport, and q 95 on target heat flux is not significant, as it primarily scales with the input power. Target heat flux is translated to target surface temperature using a simplified heat transfer model that considers the 3D target monoblock geometry and active cooling condition, excluding Be layer thermal properties. For partially detached divertor conditions, the bulk tungsten monoblock surface temperature remains below the baking temperature, which is insufficient for efficient tritium outgassing under the actively cooled ITER divertor condition. However, considering the potential thermal contact resistance between the beryllium and tungsten layers, which may significantly impact temperature distribution, the temperature of the beryllium layer can be raised to a level sufficient for efficient tritium outgassing. Therefore, the raised strike point scenario can be considered as an alternative in-vessel tritium removal technique.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of green-colour-emitting pyrotechnics as a core for 3D temperature imaging sensors inside coal boilers

Efficiency in the control operation of the boilers for coal and coal with biomass can be further improved if the flue gas temperature distribution can be better characterized. This is very difficult in these harsh environmental systems, where spatially resolved measurements are nearly impossible with solid-state sensors. Here, in this work, we evaluate the development of pyrotechnic compounds that would serve as the basis for a novel optical mapping of the temperature inside coal boilers. For this purpose, various green-colour-emitting pyrotechnics using BaCl 2 · 2H 2 O and Ba(NO 3 ) 2 as the green light source were prepared, as this colour offers a distinct signal from the combustion-based background in the boiler. These pyrotechnics were characterized using thermogravimetric analysis (TGA) and X-ray diffraction (XRD) and tested using a flat-flame burner. Furthermore, the composition was varied to evaluate the effect of different metal fuels such as Sn, Co, and Mg, as well as various binders such as ethylcellulose, shellac, parlon, and PVC on green light emission. The emission intensity and the apparent ignition temperature were strongly dependent on the metal type, with Mg showing higher intensities. On the other hand, the effect of the binder showed that the ignition behaviour, emission intensity, and spectral purity were influenced by the nature and exothermicity of the binder. The addition of other potential green light-producing materials, such as boric acid, increased the intensity of emission by 17% for a BaCl 2 · 2H 2 O-based composition. This study identified prospective compositions with intense and bright green-colour emissions that have high spectral purities.

42 ENGINEERING↗

Composition dependence of atomic order in strain-relaxed, metastable GeSn alloys

Extended x-ray absorption fine structure (EXAFS) measurements of single-crystal Ge/GeSn radial heterostructure nanowires are used to examine the effects of composition on both short-range order (SRO) and longer-range disorder in GeSn alloys. GeSn has prompted significant interest because it can achieve a direct band gap for sufficient Sn concentrations beyond the equilibrium solid solubility limit in an all-group IV system. Short-range order in this material is particularly interesting as it has been predicted to affect the band gap independent of average composition or strain effects. By independently controlling the Sn composition and GeSn thickness during chemical vapor deposition of misfitting GeSn shells around ultrathin, elastically compliant, Ge core nanowires, the elastic misfit strain in the GeSn is minimized for Sn compositions over the studied range ≈Ge 0.96 Sn 0.04 to Ge 0.88 Sn 0.12 . The degree of SRO was found to decrease with increasing Sn composition. Additionally, damping of the EXAFS signal was observed as the Sn content increased, particularly for increasingly distant atomic shells about the absorbing atom, even for scattering paths not involving Sn atoms. This result is quantified as an increase in the mean-squared relative displacement parameters of the shells. These measurements reveal the accommodation of local strain due to the presence of the highly size-mismatched Sn atoms in the Ge diamond cubic lattice (≈14%), which may have effects on the band structure of the material in addition to the influence of short-range atomic order. Comparison among the nanowire samples allows for calculation of the topological rigidity parameter, a ∗∗ , for the first-neighbor bond lengths. Furthermore, these exhibit chemically distinct values for Ge-Ge, Ge-Sn, and Sn-Sn, and they are consistent with the value a ∗∗ = 0.75 ± 0.07 confirming the general applicability of the model to alloys with both large amounts of natural misfit strain and the potential for short-range order.

Crystal structure Semiconductors Transmission elec↗

Plasma Wall Interaction with 3-D Plasma Boundaries

The interaction of the edge plasma and the material surfaces is one of the most critical challenges on the path to harness fusion power as new, fundamental energy source. This challenge typically combines the thrust to reach high density, low temperature (detached) plasmas in front of the divertor target plates as well as understanding the plasma material interaction (PMI) in particular in this regime. The combination of both research thrusts represents an extraordinarily challenging subject encompassing spatial scales spanning nanometers to meters in all states of matter and across a broad energy range. Modeling capabilities, which help to interpret data from nowadays experiments and enable extrapolation to future devices are urgently required. This is in particular true for toroidal magnetic confinement devices with three-dimensional (3D) plasma boundaries. Such plasma boundary geometries occur in tokamaks, when small amplitude magnetic perturbations are used to stabilize the unruly edge plasma or in stellarators, that are inherently 3D plasma confinement devices. In this project, the impact of 3D plasma boundaries on the plasma material interaction (PMI) was assessed. This work focused on plasma boundary conditions, in which high-density conditions at the material surfaces yield mitigation of the otherwise immense heat and particle loads that these materials would see. These so-called high recycling and eventually detached plasma regimes are of great interest for future reactor operation. In the project, key features that are unique to 3D boundaries were explored in comparison to canonically assumed axisymmetric plasma edge situations in tokamaks. In particular, the relevance of the 3D boundary situation in the extrapolation to the plasma boundary solution at ITER, the next step fusion energy experiment under construction as a multi-national, world-wide large-science experiment in southern France, has been explored. The EMC3-EIRENE plasma edge fluid and kinetic neutral transport code has been advanced to cope with the challenging and unprecedented conditions in the ITER boundary plasma including 3D fields that are planned to be used to suppress harmful edge instabilities, the so-called edge localized modes. This is a vital integration challenge for ITER and the results from this grant have provide a leading capability for this assessment. It was shown that the detachment process in a 3D edge solution for ITER follows the recycling regimes that are known from axisymmetric solutions, but that multiple plasma exhaust channels connected to the material surfaces are established which feature individual recycling characteristics. Because these channels touch the material surfaces in the divertor in a 3D geometry, the compatibility with the plasma material interaction (PMI), including erosion and impurity generation has been found to be an important part of the integration challenge. To address this, the fully 3D plasma material interaction code ERO2 has been adapted to these ITER specific geometries and a homogeneous mixing model was implemented, that allows to consider the mixing of Be and Was used at ITER in the PMI modeling. This model enhancement has been used to study non-local migration of Be in the JET ITER like wall configuration and it has been shown that with this model such complex migration processes in ITER relevant plasma shapes and with ITER relevant plasma boundary conditions can be addressed. The combined modeling approach using EMC3-EIRENE as a plasma boundary transport code and the ERO2 specialized PMI model will be an asset for the continued preparations of ITER operation as well as for Fusion Pilot Plant efforts that have emerged in the U.S. during the evolution of this grant. The predictive capability of this numerical tool has been validated at the DIII-D US national fusion facility. Here, dedicated plasma edge diagnostics were implemented to measure the impurity household around a 3D edge plasma during ELM suppression by 3D fields. Dedicated experiments with local material probes using these diagnostics and the state-of-the-art suite of boundary measurements at DIII-D have shown that the 3D perturbation of the plasma edge that is excreted by such 3D control fields yield a perturbation of the plasma boundary flux structure and hence also of the resulting PMI. The 3D boundary plasma is composed out of helical magnetic flux channels that intersect the divertor targets at an angle relative to the main guiding field, i.e., the toroidal magnetic field component of the tokamak. A similar effect has been measured as well on limiter surfaces during the startup campaign at the new stellarator experiment Wendelstein 7-X. These experiments ad initial analysis with the ERO plasma material interaction model, suggested that the place of erosion for a given particle from the surface and its re-deposition can be different in such 3D field geometries yielding potentially a significant level of net-erosion. This is not the case for axisymmetric solutions, where it was shown in the past that the eroded particles are effectively re-deposited into gaps produced by erosion at the same position and hence the net-erosion levels are small. For ITER, the quest to suppress the ELMs and at the same time maintain the integrity of the divertor is an issue, which these fundamental findings will help to resolve. The coupling of this work to the extrapolation in the ITER program has been addressed by both the PI and the lead numerical scientist being ITER Science Fellows in the duration of the contract and forward. A second focus in the exploration of 3D boundary effects on tokamaks and stellarators has been set on the measurement of helium exhaust features with such 3D fields. This is important because He represents the ash of the fusion process and needs to be exhausted. It was shown that 3D field application compatible with suppression of ELMs yields an increase of the helium exhaust performance. The ratio of the effective helium confinement time over the energy confinement time was reduced by almost 50% which demonstrated that the impact of helium accumulation in the plasma core with respect to the confinement of energy to sustain the fusion reaction is significantly improved with such 3D control fields. It was shown that this is the case for tokamaks as well as stellarators. At the Large helical Device in Japan, a similar enhancement of the helium exhaust features when small amplitude additional 3D fields were applied was measured. This is an important additional function of 3D field application and its impact on ITER is presently being studied in combination with investigations of helium exhaust in 3D field geometries of stellarator devices.

3D plasma edge transport↗

Connecting suborganismal data to bioenergetic processes: killifish embryos exposed to a dioxin-like compound

A core challenge for ecological risk assessment is to integrate molecular responses into a chain of causality to organismal or population level outcomes. Bioenergetic theory may be a useful approach for integrating suborganismal responses to predict organismal level responses that influence population dynamics. In this work, we describe a novel application of Dynamic Energy Budget (DEB), theory in the context of a toxicity framework (Adverse Outcome Pathways, AOP) to make quantitative predictions of chemical exposures to individuals, starting from suborganismal data. We use early life stage exposure of Fundulus heteroclitus to dioxin-like chemicals (DLCs) and connect AOP Key Events (KEs) to DEB processes through “damage” that is produced at a rate proportional to the internal toxicant concentration. We use transcriptomic data of fish embryos exposed to DLCs to translate molecular indicators of damage into changes in DEB parameters (damage increases somatic maintenance costs) and use DEB models to predict sublethal and lethal effects of young fish. By changing a small subset of model parameters, we predict the evolved tolerance to DLCs in some wild F. heteroclitus populations, a data set not used in model parameterization. The differences in model parameters points to reduced sensitivity and altered damage repair dynamics as contributing to this evolved resistance. Our methodology has potential extrapolation to untested chemicals of ecological concern.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Colloid-facilitated transport of 238 Pu, 233 U and 137 Cs through fractured chalk: Laboratory experiments, modelling, and implications for nuclear waste disposal

The influence of montmorillonite colloids on the mobility of 238 Pu, 233 U and 137 Cs through a chalk fracture was investigated to assess the transport potential for radioactive waste. Radioisotopes of each element, along with the conservative tracer tritium, were injected in the presence and absence of montmorillonite colloids into a naturally fractured chalk core. In parallel, batch experiments were conducted to obtain experimental sorption coefficients (K d , mL/g) for both montmorillonite colloids and the chalk fracture material. Breakthrough curves were modelled to determine diffusivity and sorption of each radionuclide to the chalk and the colloids under advective conditions. Uranium sorbed sparingly to chalk (log K d = 0.7 ± 0.2) in batch sorption experiments. 233 U(VI) breakthrough was controlled primarily by the matrix diffusion and sorption to chalk (15 and 25% recovery with and without colloids, respectively). Cesium, in contrast, sorbed strongly to both the montmorillonite colloids and chalk (batch log K d = 3.2 ± 0.01 and 3.9 ± 0.01, respectively). The high affinity to chalk and low colloid concentrations overwhelmed any colloidal Cs transport, resulting in very low 137 Cs breakthrough (1.1–5.5% mass recovery). Batch and fracture transport results, and the associated modelling revealed that Pu migrates both as Pu (IV) sorbed to montmorillonite colloids and as dissolved Pu(V) (7% recovery). Transport experiments revealed differences in Pu(IV) and Pu(V) transport behavior that could not be quantified in simple batch experiments but are critical to effectively predict transport behavior of redox-sensitive radionuclides. Finally, a brackish groundwater solution was injected after completion of the fracture flow experiments and resulted in remobilization and recovery of 2.2% of the total sorbed radionuclides which remained in the core from previous experiments. In general, our study demonstrates consistency in sorption behavior between batch and advective fracture transport. Finally, the results suggest that colloid-facilitated radionuclide transport will enhance radionuclide migration in fractured chalk for those radionuclides with exceedingly high affinity for colloids.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Bright Silicon Nanocrystals from a Liquid Precursor: Quasi-Direct Recombination with High Quantum Yield

Silicon nanocrystals (SiNCs) with bright bandgap photoluminescence (PL) are of current interest for a range of potential applications, from solar windows to biomedical contrast agents. Here, we use the liquid precursor cyclohexasilane (Si 6 H 12 ) for the plasma synthesis of colloidal SiNCs with exemplary core emission. Here, through size separation executed in an oxygen-shielded environment, we achieve PL quantum yields (QYs) approaching 70% while exposing intrinsic constraints on efficient core emission from smaller SiNCs. Time-resolved PL spectra of these fractions in response to femtosecond pulsed excitation reveal a zero-phonon radiative channel that anticorrelates with QY, which we model using advanced computational methods applied to a 2 nm SiNC. Our results offer additional insight into the photophysical interplay of the nanocrystal surface, quasi-direct recombination, and efficient SiNC core PL.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

H 2 BZmacropa-NCS: A Bifunctional Chelator for Actinium-225 Targeted Alpha Therapy

Actinium-225 ( 225 Ac) is one of the most promising radionuclides for targeted alpha therapy (TAT). With a half-life of 9.92 days and a decay chain that emits four high-energy α particles, 225 Ac is well-suited for TAT when conjugated to macromolecular targeting vectors that exhibit extended in vivo circulation times. The implementation of 225 Ac in these targeted constructs, however, requires a suitable chelator that can bind and retain this radionuclide in vivo. Previous work has demonstrated the suitability of a diaza-18-crown-6 macrocyclic chelator H 2 macropa for this application. Building upon these prior efforts, in this study, two rigid variants of H 2 macropa, which contain either one (H 2 BZmacropa) or two (H 2 BZ 2 macropa) benzene rings within the macrocyclic core, were synthesized and investigated for their potential use for 225 Ac TAT. The coordination chemistry of these ligands with La 3+ , used as a nonradioactive model for Ac 3+ , was carried out. Both NMR spectroscopic and X-ray crystallographic studies of the La 3+ complexes of these ligands revealed similar structural features to those found for the related complex of H 2 macropa. Thermodynamic stability constants of the La 3+ complexes, however, were found to be 1 and 2 orders of magnitude lower than those of H 2 macropa for H 2 BZmacropa and H 2 BZ 2 macropa, respectively. The decrease in thermodynamic stability was rationalized via the use of density functional theory calculations. 225 Ac radiolabeling and serum stability studies with H 2 BZmacropa showed that this chelator compares favorably with H 2 macropa. Based on these promising results, a bifunctional version of this chelator, H 2 BZmacropa-NCS, was synthesized and conjugated to the antibody codrituzumab (GC33), which targets the liver cancer biomarker glypican-3 (GPC3). The resulting GC33-BZmacropa conjugate and an analogous GC33-macropa conjugate were evaluated for their 225 Ac radiolabeling efficiencies, antigen-binding affinities, and in vivo biodistribution in HepG2 liver cancer tumor-bearing mice. Although both conjugates were comparably effective in their radiolabeling efficiencies, [ 225 Ac]Ac-GC33-BZmacropa showed slightly poorer serum stability and biodistribution than [ 225 Ac]Ac-GC33-macropa. Altogether, these results establish H2BZmacropa-NCS as a new bifunctional chelator for the preparation of 225 Ac radiopharmaceuticals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Time dependent model for non-inductive ECH X-I current ramp-up for SHPD tokamak facility

Non-inductive (NI) plasma current start-up and ramp-up is an important research topic for spherical tokamak (ST) based reactors and fusion pilot plant (FPP). For a compact FPP, the OH flux availability is highly restricted due to its compact geometry. Efficient fundamental extraordinary mode (X-I) electron cyclotron heating (ECH) current start-up and ramp-up regime was identified for a reactor-like high toroidal magnetic field range which has more than a hundred times higher current drive efficiency compared to more conventional ECH methods for the relevant start-up temperature range. High current drive efficiency is possible due to the strong X-I fundamental ECH interaction only with unidirectional passing electrons constrained by the wave accessibility conditions. Here, we extend the X-I electron cyclotron current drive (ECCD) investigation to a time dependent model to simulate the non-inductive current ramp-up to 10 MA for the Sustained High-Power Density (SHPD) facility. As the X-I ECCD driven current I EC rises, due to the back EMF driven negative current, the net plasma current Ip rises more slowly with the current resistive time scale. For tokamak confinement time (both L-mode and H-Mode) which tends to rise with I p , a positive feedback results and even with constant applied ECH power, T e0 , I EC , and I p can continue to rise to very high values. However, in a realistic situation, the T e0 rise should saturate due to a number of factors such as enhanced core radiation and increased power loss at high temperature. To simulate this effect, we adopt a maximum T e0 model which would limit the temperature rise to certain T e0 . With this model, we investigated the current ramp-up for various ECH power levels and the maximum T e0 of 15, 20, and 25 keV. We find that while the power required is reduced with increasing T e0 limit due to increased current drive efficiency, the time to reach 10 MA tends to go up due to the reduced plasma resistivity for the higher T e0 limit. We also find for a given Te0 limit, the time to reach 10 MA tends to be reduced by increasing the applied ECH power by over driving the current ramp-up where I EC is driven at significantly higher level than 10 MA. While the current ramp-up time may not be an issue for the steady-state reactor systems, if it is desirable to minimize the current ramp-up time, it is prudent to have a sufficient ECH power for current over-drive and have some T e0 limiting tools such as impurity seeding for enhanced radiation. In conclusion, a well-controlled NI ECH start-up also has a potential of improving the tokamak start-up reliability and avoid run-away electrons while the NI off-axis current drive could enhance MHD stability and plasma performance improvements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Structure and Behavior of the Ni End-Member Schreibersite Ni 3 P under Compression to 50 GPa

To better understand the potential presence of light element alloys of Fe and Ni in the Earth’s interior, the crystal structure and compressional behavior of the Ni-P binary compound, schreibersite (Ni 3 P), have been investigated using synchrotron X-ray diffraction experiments. Both powder and two single-crystal samples of synthetic Ni 3 P (in different orientations with respect to the loading axis of the diamond anvil cell) were compressed up to approximately 50 GPa at ambient temperature. The compressional data obtained for Ni3P were fitted with a 3rd order Birch–Murnaghan equation of state. All data indicated that the c/a ratio of unit cell parameters remained approximately constant up to about 30 GPa but then increased progressively with pressure, exhibiting a second slight discontinuity at approximately 40 GPa. The changes in unit cell parameters at ~30 GPa and ~40 GPa suggested discontinuous changes in magnetic ordering. Moreover, the threshold of these subtle discontinuities is sensitive to the stress state and orientation of the crystal in the diamond anvil cell. This study is the first report on the compressional behavior of both powder and single-crystal schreibersite at high-pressure (up to 50 GPa). It offers insights into the effects of Ni 3 P components on the compressional behavior of the Earth’s core.

36 MATERIALS SCIENCE↗

Assessment of the Griffin Reactor Multiphysics Application Using the Empire Micro Reactor Design Concept

In late 2019, INL and ANL agreed to jointly develop the reactor physics code named Griffin based on the integration of the two code suites, MAMMOTH/Rattlesnake (INL) and MC2 - 3/PROTEUS (ANL). Griffin is being developed based on the MOOSE framework and MOOSE quality assurance procedures. This decision was made to be able to allow DOE-NE to efficiently invest funding to this area and to provide effective and timely support for existing and potential users; the latter includes industry and government organizations who are developing various types of advanced reactors in the near and long term. Since MAMMOTH/Rattlesnake has been developed based on the MOOSE framework, the INL/ANL Griffin development team agreed to build Griffin beginning with a merger of MAMMOTH and Rattlesnake into a single code and moving forward by implementing capabilities from the PROTEUS suite into Griffin. Moving forward, both ANL and INL efforts are equally invested in the Griffin project, with management support, to provide an advanced reactor multiphysics tool to assist in reactor design, optimization, and safety analysis. Much work remains in moving Griffin forward to migrate PROTEUS capabilities and to optimize performance to meet user needs. The main objective of this work is to assess the current status of Griffin capabilities in terms of performance and accuracy, to determine priorities for PROTEUS migration, and to identify capabilities and features to improve for supporting the code integration effort. For this assessment, the Empire micro reactor problem that was developed in the ARPA-E MEITNER program was selected as an advance reactor concept of interest to the technical community. The Empire reactor problem was expanded from its original incomplete specification to be a small heat-pipe-cooled micro reactor core with ~113 cm radius and 70 cm in height, composed of 18 fuel assemblies, 12 control drums, and beryllium radial and axial reflectors. In the current model, using 5 cm axial reflectors specified in the original Empire assembly model, more than 10% of neutrons leak axially and through the empty center safety hole, as well as through heat pipe channels in fuel assembly elements that extend through the top reflector region. Several calculation models of the core were defined for systematic assessment, including 2-D and 3-D fuel assemblies and whole cores with cylindrical boundaries. Cross sections were generated using Serpent 2, and meshes were produced using the Argonne mesh tool or the INL neutronics meshing tools combined with CUBIT. Cross sections and meshes were converted to the ISOXML and Exodus formats, respectively, so that Griffin and PROTEUS could use consistent data for solving the reactor problems. With the prepared cross sections and meshes, PROTEUS was run first to ensure that all input data were correctly generated and input options in terms of angle, mesh, and energy group were accurately determined. Comparisons against Serpent 2 solutions were made in terms of eigenvalue and pin power. The same calculations and comparisons were then conducted using Griffin. For the fuel assembly and whole core problems, the PROTEUS eigenvalues agreed well with reference Serpent 2 solutions within 100 and 30 pcm, respectively, and pin power differences relative to Serpent 2 were overall less than 2.2% and RMS 0.8% for the whole core models. This indicated that all input data were properly prepared. Using the same data, Griffin was run selecting the SAAF-CFEM SN solver with Legendre-Gaussian quadrature and NDA and DSA for acceleration. It was found that the SAAF-CFEM solver of Griffin required finer meshes to achieve eigenvalue and pin power solutions in good agreement with Serpent 2, consequently requiring more memory requirement and longer computation time. On the other hand, the SPH-Diffusion 2-D core calculations performed using Griffin were able to recover the exact eigenvalue from the reference Serpent 2 solutions, resulting in a pin-power distribution with an RMS of 0.6% and maximum absolute difference of less than 1.4%. The runtimes for SPH-Diffusion for the 2-D core were less than 3 minutes on 40 cores. During this evolution of this evaluation, many updates were made in Griffin by the Griffin development team of INL (focusing on software updates) and ANL (reviewing and supporting software updates) to complete this assessment. Observations from the code assessment are presented in the conclusion section of this report, followed by a discussion of recommendations for future work.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Atomic-Scale STEM Analysis Shows Structural Changes of Au–Pd Nanoparticles in Various Gaseous Environments

We report dynamical restructuring effects in free-standing Au 0.75 Pd 0.25 nanoparticles occurring in gaseous environments at elevated temperatures. The freshly prepared sample was found to have a core–shell structure with a Pd-rich phase on the surface. The evolution of sample composition and morphology under exposure to 1 bar of pure gases, namely O 2 , H 2 , air, CO, and CO 2 , at different temperatures was studied using in situ scanning transmission electron microscopy (STEM) and energy-dispersive X-ray spectroscopy (EDS). We observed sharper facets on the surface of the particles under O 2 or air at 400 °C. Small islands of Pd were present on the surface, although some Pd was redistributed inside the bulk when the temperature was increased under O 2 . Subtle changes in surface roughness were noted when O 2 was substituted with H 2 at 400 °C, an observation correlated to density functional theory (DFT) calculations. The particles lost clean surface facets when CO was introduced at room temperature and at 200 °C. No substantial changes could be observed after exposure to CO 2 at 250 °C. The adsorption of CO molecules on the surface modifies the surface of the particles and decreases the facet prevalence. Furthermore, these in situ observations show how gases can induce subtle modification of the surface of nanocatalysts, potentially impacting their chemical properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗