Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reactor Physics Method”

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 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Electrochemical Synthesis of Zeolite Coatings with Controlled Crystal Polymorphism and Self-Regulating Growth

Zeolite coatings are studied as molecular sieves for membrane separation, membrane reactors, and chemical sensor applications. They are also studied as anticorrosive films for metals and alloys, antimicrobial and hydrophobic films for heating, ventilation, and air conditioning, and dielectrics for semiconductor applications. Zeolite coatings are synthesized by hydrothermal, ionothermal, and dry-gel conversion approaches, which require high process temperatures and lengthy times (ranging from hours to days). Here, we report the first zeolite coatings synthesized via electrochemical deposition on a cathodic electrode, with controlled crystal polymorphism achieved within subhourly duration. We demonstrate this approach by developing sodium zeolite (e.g., sodalite (SOD), NaA (LTA), and Linde Type N (LTN)) coatings on a titanium electrode and extending the synthesis method to porous stainless steel. The coating morphology and crystallinity depend on the temperature, time, and applied current. The coating thickness is independent of the applied current, showing the presence of a self-regulating mechanism to ensure a uniform coating thickness across the metal surface. The electrochemical zeolite growth mechanism was elucidated with high-resolution transmission electron microscopy, and applications of the resultant zeolite coatings for oil/water separation and ethanol/water pervaporation were exploited. Electrochemical synthesis represents a novel, simple, fast, and environmentally friendly approach to preparing zeolite coatings. It can potentially be generalized for developing zeolite materials with diverse framework structures, morphologies, and orientations for substrates with complicated geometries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

A Machine Learning Method for the Forensics Attribution of Separated Plutonium

Plutonium (Pu) source attribution would be a powerful tool to support nuclear nonproliferation efforts. This capability to find the source of a Pu sample would act as a deterrent to smuggling efforts, and also help regulatory agencies verify declared nuclear activities. Work at Texas A&M University yielded a nuclear forensics methodology, which is capable of determining separated Pu’s reactor of origin, fuel burnup, and the time since irradiation (TSI)—three parameters of interest. The methodology used a set of ten intra-element isotopic ratios found in separated Pu, which was compared to a library of isotopic ratio values produced using neutronics simulations for reactors of interest. By calculating the probability that unknown Pu sample’s isotopic ratio set matched a set in the library, the methodology could predict the three parameters of interest of the sample. One shortcoming of this methodology was an inability to correctly attribute spoofed Pu, where Pu sourced from two different reactors or two different fuel burnup levels are mixed. A new methodology to rectify this vulnerability using machine learning (ML) technique is developed, instead of the maximum likelihood calculation previously used and the results are satisfactory. The ML approach leverages the existing simulated data for training the algorithm, but use them efficiently by only using intra-element isotope ratios that contribute to the attribution one of the three parameters at a time. Previously, all isotope ratios were used to attribute all three parameters together. The new methodology attributes the Pu parameters in three steps, one for each parameter, rather than resolving all of the three parameters simultaneously like the previous maximum likelihood approach. First, a support vector machine classifier with a set of seven isotopic ratios finds the reactor of origin and a set of regression models trained using gaussian process predicts the burnup with a different set of seven isotopes. Finally, TSI is calculated analytically using decay equations. Thus far, the new methodology is capable of attributing pure Pu samples and has been validated using experimental data. The next step will to be augment the classifier training data set with spoofed Pu data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Electrons to Reactors Multiscale Modeling: Catalytic CO Oxidation over RuO 2

First-principles kinetic Monte Carlo (1p-kMC) simulations for CO oxidation on two RuO 2 facets, RuO 2 (110) and RuO 2 (111), were coupled to the computational fluid dynamics (CFD) simulations package MFIX, and reactor-scale simulations were then performed. 1p-kMC coupled with CFD has recently been shown as a feasible method for translating molecular scale mechanistic knowledge to the reactor scale, enabling comparisons to in situ and online experimental measurements. Only a few studies with such coupling have been published. This work incorporates multiple catalytic surface facets into the scale-coupled simulation, and three possibilities were investigated: the two possibilities of each facet individually being the dominant phase in the reactor, and also the possibility that both facets were present on the catalyst particles in the ratio predicted by an ab initio thermodynamics-based Wulff construction. When lateral interactions between adsorbates were included in the 1p-kMC simulations, the two surfaces, RuO 2 (110) and RuO 2 (111), were found to be of similar order-of-magnitude in activity for the pressure range of 1 × 10 –4 bar to 1 bar, with the RuO 2 (110) surface-termination showing more simulated activity than the RuO 2 (111) surface-termination. Coupling between the 1p-kMC and CFD was achieved with a lookup table generated by the error-based modified Shepard interpolation scheme. Isothermal reactor scale simulations were performed and compared to two separate experimental studies, conducted with reactant partial pressures of ≤0.1 bar. Simulations without an isothermality restriction were also conducted and showed that the simulated temperature gradient across the catalytic reactor bed is <0.5 K, which validated the use of the isothermality restriction for investigating the reactor-scale phenomenological temperature dependences. The approach with the Wulff construction based reactor simulations reproduced a trend similar to one experimental data set relatively well, with the (110) surface being more active at higher temperaures; in contrast, for the other experimental data set, our reactor simulations achieve surprisingly and perhaps fortuitously good agreement with the activity and phenomenological pressure dependence when it is assumed that the (111) facet is the only active facet present. Lastly, the active phase of catalytic CO oxidation over RuO 2 remains unsettled, but the present study presents proof of principle (and progress) toward more accurate multiscale modeling from electrons to reactors and new simulation results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

On the Interpretation and Correlation of High-Temperature Ignition Delays in Reactors with Varying Thermodynamic Conditions

Ignition delay time (IDT) is a useful global metric for fuel performance screening and a major target for kinetic modeling. Measurements of IDT are conceptually straightforward, however interpretation could be complicated, especially for systems with changing temperature and pressure. Some experimental conditions in the high repetition rate miniature shock tube (HRRST) exhibit complex temperature and pressure state histories. To better interpret and correlate IDTs, especially those obtained in reactors with varying thermodynamic conditions, an inverse Livengood-Wu (L-W) integral technique is applied to de-convolve the constant condition IDTs from measured IDTs in the HRRST using information on the varying state history. In this paper, the approximate problem is demonstrated using only the measured pressure history leading up to ignition, where the temperature history is estimated based on an isentropic assumption. The IDTs of several fuels were first measured in the HRRST including isooctane and acetone to represent those fuels without strong low temperature chemistry effects. Based on the measured pressure and approximated temperature history, measurements of IDTs in the HRRST compare very favorably with those measured using more conventional techniques, including conventionally sized shock tubes, via the inverse L-W correlation with relaxed Newton iteration and Genetic Algorithm (GA). Furthermore, this study demonstrates the feasibility of using a high throughput ignition testing facility, like the HRRST to extract the constant state IDT measured in regular shock tubes, assisted by the inverse L-W correlation. It is expected with an independent temperature measurement available in the future, IDTs in a broader range of thermodynamic conditions and effects of heat loss could be better resolved.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Machine learning molecular dynamics simulations toward exploration of high-temperature properties of nuclear fuel materials: case study of thorium dioxide

Predicting materials properties of nuclear fuel compounds is a challenging task in materials science. Their thermodynamical behaviors around and above the operational temperature are essential for the design of nuclear reactors. However, they are not easy to measure, because the target temperature range is too high to perform various standard experiments safely and accurately. Moreover, theoretical methods such as first-principles calculations also suffer from the computational limitations in calculating thermodynamical properties due to their high calculation-costs and complicated electronic structures stemming from f-orbital occupations of valence electrons in actinide elements. Here, we demonstrate, for the first time, machine-learning molecular-dynamics to theoretically explore high-temperature thermodynamical properties of a nuclear fuel material, thorium dioxide. The target compound satisfies first-principles calculation accuracy because f-electron occupation coincidentally diminishes and the scheme meets sampling sufficiency because it works at the computational cost of classical molecular-dynamics levels. We prepare a set of training data using first-principles molecular dynamics with small number of atoms, which cannot directly evaluate thermodynamical properties but captures essential atomistic dynamics at the high temperature range. Then, we construct a machine-learning molecular-dynamics potential and carry out large-scale molecular-dynamics calculations. Consequently, we successfully access two kinds of thermodynamic phase transitions, namely the melting and the anomalous λ transition induced by large diffusions of oxygen atoms. Furthermore, we quantitatively reproduce various experimental data in the best agreement manner by selecting a density functional scheme known as SCAN. Our results suggest that the present scale-up simulation-scheme using machine-learning techniques opens up a new pathway on theoretical studies of not only nuclear fuel compounds, but also a variety of similar materials that contain both heavy and light elements, like thorium dioxide.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

A Subcritical Assembly for Training and Education Use at the Oak Ridge National Laboratory

A near-final design has been completed for a new subcritical assembly at Oak Ridge National Laboratory (ORNL). This subcritical assembly will support the Nuclear Criticality Safety Program (NCSP) training and education program to augment current subcritical assembly capabilities at Sandia National Laboratories (SNL) and at the National Criticality Experiments Research Center (NCERC). Other collaborations have been proposed with the Sunrise Consortium and others. This proposed subcritical assembly uses legacy AGN-201M research reactor fuel plates that are available from the Y-12 National Security Complex. This subcritical assembly contains approximately 620 grams of 235 U as UO 2 particulates distributed homogeneously in radiation-stabilized polyethylene. The fuel plates will have a graphite neutron reflector to obtain a core multiplication, M, from 10 to 50, corresponding to a k eff of 0.9 to 0.95, respectively. The subcritical assembly will be able to support at least four experiments for the training courses: (1) the addition of fissile material to increase neutron multiplication (mass), (2) a core separation experiment (interaction), (3) the effect of adding moderators to the core (moderation), and (4) the effect of adding neutron absorbers to the core (poison/absorption). The proposed ORNL subcritical assembly will provide additional student bandwidth to support NCSP training courses and to provide for some backup capacity if nuclear facility operations are disrupted at SNL or NCERC due to some unforeseen issue. The subcritical assembly will be designed to be an inherently safe subcritical assembly—subcritical under all normal and abnormal conditions—and it will provide the capability to conduct hands-on training to support NCSP and general nuclear criticality safety (NCS) staff training and qualification goals.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

Facile One-Pot Synthesis of Pd@Pt 1L Octahedra with Enhanced Activity and Durability toward Oxygen Reduction

A successful strategy for reducing the content of Pt without compromising the activity of a Pt-based catalyst is to deposit Pt as an ultrathin overlayer on the surface of another metal. In this paper, we report a facile one-pot synthesis of Pd@Pt 1L (1L: one atomic layer) core–shell octahedra using a solution-phase method. The success of this method relies on the use of metal precursors with markedly different reduction kinetics. In a typical synthesis, the ratio between the initial reduction rates of the Pd(II) and Pt(II) precursors differed by almost 100 times, favoring the formation of Pd–Pt bimetallic octahedra with a core–shell structure. The reduction of the Pt(II) precursor at a very slow rate and the use of a high temperature allowed the deposited Pt atoms to spread and cover the entire surface of Pd octahedral seeds formed in the initial stage. More importantly, we were able to scale up this synthesis using continuous-flow reactors without compromising product quality. Compared to a commercial Pt/C catalyst, the Pd@Pt 1L core–shell octahedra showed major augmentation in terms of catalytic activity and durability for the oxygen reduction reaction (ORR). After 10000 cycles of accelerated durability test, the core–shell octahedra still exhibited a mass activity of 0.45 A mg –1 Pt . We rationalized the experimental results using DFT calculations, including the mechanism of synthesis, ORR activities, and possible Pd–Pt atom swapping to enrich the outermost layer with Pd. Specifically, the as-synthesized Pd@Pt 1L octahedra tended to take a slightly mixed surface composition because the deposited Pt atoms were able to substitute into Pd upon deposition on the edges; ORR energetics were more favorable on pure Pt shells as compared to significantly mixed Pd–Pt shells, and the activation energy barriers calculated for the Pd–Pt atom swapping were too prohibitive to significantly alter the surface composition of the as-synthesized Pd@Pt 1L octahedra, helping sustain their activity for prolonged operation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

SCALE 6.3 Validation: Radiation Shielding

Safe and reliable use of scientific and engineering computer codes requires validation for the types of applications in which they will be used. An example in the nuclear reactor engineering and licensing field is radiation transport employed in shielding analyses. The validity of computer codes for shielding applications is demonstrated in this report for SCALE version 6.3.0. Representative benchmarks corresponding to shielding analyses are selected for the validation study. Typical measurement results analyzed from these benchmarks include neutron fluxes, detector count rates, detector energy response functions, neutron and gamma dose rates, neutron activation rates and activities, neutron leakage fluxes, and skyshine dose rates. Thousands of points of comparison between measurement and calculation are presented in this work. Other than rare outliers typically explained by either a lack of information or large uncertainties in the experiment conditions, material, or dimensions, the Monaco with Automated Variance Reduction using Importance Calculations (MAVRIC) radiation transport computer code with built-in variance reduction methods distributed with the SCALE computer code system agrees well with the measurement results. In selected benchmarks, MAVRIC is also compared to Monte Carlo N- Particle® (MCNP® ) 1 calculations. Both computer codes generally agree well within the estimated uncertainties. With the release of SCALE 6.3.0, Shift was integrated as an alternative transport solver in MAVRIC, denoted MAVRIC-Shift. Although the traditional MAVRIC using Monaco was used primarily in this validation study, many results have also been generated using MAVRIC-Shift. Agreement between MAVRIC-Monaco and MAVRIC-Shift is generally very good. The benchmarks presented in this report were obtained from reliable sources such as the International Criticality Safety Benchmark Evaluation Project Handbook, the Shielding Integral Benchmark Archive & Database, and other shielding validation work found in the literature. Additional datapoints and benchmarks will be added to future versions of this report to expand the shielding validation suite.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2019

Quantum-Classical Science for the Plasma-Material Interface in NSTXU (Final Technical Report)

Plasma facing surfaces (PFS) have a profound influence on plasma performance and choice of material for PFSs is a burning issue for tokamak reactors. Using lithium wall conditioning in NSTX has resulted in a beneficial influence on plasma behavior including control of hydrogen recycling and improvements in energy confinement, longer confinement times, edge density control, and the elimination of Edge Localized Modes. A major challenge for the NSTX-U (increased power flux, increased campaign time) as well as for other magnetic fusion devices is to extend high-performance plasmas for very long durations, and to integrate this high performance with plasma facing components (PFCs) that can withstand very high heat and particle fluxes while maintaining structural integrity with minimal retention of fusion fuel. This project directly contributed to filling up the knowledge gaps and data for the dynamics of the plasma-facing interfaces in the NSTX-U at nano and nano-to-micro scales. This was achieved by both deciphering the poorly understood nanoscale phenomenology of plasma-surface interactions and providing the connection and comparison with the experiments which provide understanding at the mesoscopic scales. We modelled the dynamic response and evolution of a mixed material surfaces (D, Li, C, B, O) to bombardment by plasma particles, and investigated microstructure changes, erosion, surface chemistry, deuterium implantation and permeation. The main objectives of the project are i) a comparison of Li and B effects when deposited on carbon, and ii) the role of oxygen and other impurities (e.g. carbon, boron) in the lithium performance. In addition to predicting and understanding the phenomenology of the processes, the project achieved deliverables on the i) plasma induced erosion rates of PFCs (outer wall, divertor, etc.), including chemical and physical sputtering yields at various temperatures (300-700K), ii) deuterium uptake/recycling rates, and iii) accurate probability rates for the processes and reactions in i-ii which were used to understand the experimental results of the J.P. Allain and B. Koel NSTX-U projects. Lithium was a distinctive focus of the project. It easily charges in contact with other fusion material atoms, giving partially or fully its external electron cloud to other atoms in its vicinity and polarizing the mixed material. Our research strategy used experimentally validated multiscale theory applied to the large, multicomponent material surface systems. With the use of US supercomputing facilities we simulated evolution of atomic charges, chemistry, retention of deuterium, sputtering and morphology by introducing the corrections to the classical molecular dynamics, like are electronegativity equalization method (EEM) and quantum-classical molecular dynamics (QCMD). The goal was to understand and quantify the dynamics of the creation and evolution of the plasma-material interface under irradiation by atoms and molecules at boronized carbon, lithiated carbon, and boronized-lithiated carbon for both solid and self-healing liquid-metal divertor.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Development of an Ultrahigh-bandwidth Phase Contrast Imaging System for detection of electron scale turbulence and Gigahertz Radio-Frequency Waves

The study of waves and turbulence is vital to the development of future reactor-grade plasma devices developed in the quest for fusion energy. These fluctuations are responsible for moving heat and particles across the magnetic field, and a predictive understanding of them is needed to achieve the density and temperature required to sustain a plasma fusion reaction. While many techniques have been developed for measuring waves and fluctuations, every measurement method has limitations. There are relatively few techniques for measuring very high frequency fluctuations, such as radio frequency waves injected to heat the plasma, unstable waves driven by suprathermal particles, or short wavelength electrostatic waves driven by electron temperature or density gradients. The present project builds upon the proven phase contrast imaging (PCI) technique to extend the response of the diagnostic by orders of magnitude in frequency and almost a factor of ten in spatial resolution. PCI provides a measurement of electron density based on small angle scattering of a CO-2 laser beam by using optical techniques to render a phase shift as an intensity change on a detector. Due to the telecommunications revolution, technological development by manufacturers has focused on components in the near -infrared, so that high- quality lasers and detectors at 1.55 µm are readily available. Shifting PCI design to a new, shorter wavelength has numerous advantages and challenges. Similar detector performance is available with room-temperature arrays with GHz bandwidth, while the detector arrays for 10.6 µm required liquid nitrogen cooling and were therefore limited to a bandwidth of about 1 MHz. Shifting to a shorter wavelength reduces the angle at which the laser beam scatters off of plasma waves, allowing more such scattered components to pass through the aperture of the vacuum vessel port and be collected by the PCI, which increases the spatial resolution. Concomitant with these benefits, various questions of performance arise. At shorter laser wave- length, the sensitivity to mirror and lens quality is increased, the contribution of the laser to the overall system noise is increased, and the sensitivity to vibrations is increased. The custom optical components at the heart of the PCI technique were required to be properly scaled for the shorter wavelength, so fabrication technologies needed to be explored. This project was designed to show that a low noise, high response PCI system at 1.55 µm could be constructed and operated, and then to quantify potential issues to allow extrapolation to a full-size production system providing physics measurements on a large plasma device. The first stage, producing the custom optical component called a Phase Plate, was successfully achieved using two methods. First, an easily reproduced masking and coating technique was able to produce good phase plates with the required parameters. Second, a nanofabrication technique was found to produce extremely high quality phase plates at a competitive price. The PCI constructed with the new phase plates and 1.55 µm laser was found to provide excellent wavelength measurements with the theoretically expected response. The sensitivity to optical surface quality was found to be in line with previous measurements at 10.6 µm. The observed signal-to-noise ratio was similar to the theoretically expected value. The effect of vibrations on PCI was studied with the first measurement of the effect of beam motion on PCI response and comparison to theory, allowing for a quantitative prediction of the effect of vibrations on a production PCI system and the requirements for improved beam stabilization. PCI is an extremely cost-effective method to provide a low noise, absolutely calibrated measurement of plasma fluctuations across a wide spatial scale. This project has shown that a 1.55 µm PCI using modern techniques and components is less expensive than the 10.6 µm PCI of a few years ago, with the largest savings in phase plate fabrication and the infrared detector array.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Plan to Verify and Validate Multi-Hazard Risk-Informed Margin Management Methods and Tools

This report describes the development plan for verification and validation (V&V) of multihazard risk informed margin management (RIMM) methods and tools. The plan will include using new INL developed tool sets coupled with the new multi-partner External Hazards Experimental Group (EHEG) coordinated by the Idaho National Laboratory (INL) within the Risk-Informed Safety Margin Characterization (RISMC) technical pathway of the Light Water Reactor Sustainability Program. The RISMC pathway is developing a suite of new tools and methods for advanced evaluation of facility risk. The RISMC toolkit, includes tools and methods focused on evaluation of risk from external hazards (e.g, seismic and flooding events), which have been shown to be dominant risk contributors in probabilistic risk assessments performed for operating nuclear power plants (NPPs). The external events activity (EEA) within the RISMC pathway is tasked with developing tools and methods focused on evaluation of multi-hazard external risk (e.g, seismic and flooding events). The EHEG will provide technical expertise and experimental large-scale testing data needed for development and validation of tools and methods in the RISMC toolkit for external hazard safety evaluations. The EEA within RISMC has two key elements: (1) an organizational and research framework provided by INL and (2) a coordinated group of university and national laboratory partners with the complementary expertise and experimental capabilities needed to conduct large-scale external hazard-focused experiments. This cooperative group will allow INL to leverage a range of existing capabilities to meet the unique needs of RISMC tool and method development. In addition, the capabilities of the EHEG could be used to address needs of other national laboratories. Currently, there is limited data available for development and validation of the tools and methods being developed in the RISMC Toolkit. The EHEG is being developed to obtain high quality, large-scale experimental data validation of RISMC tools and methods in a timely and cost-effective way. The group of universities and national laboratories that will eventually form the EHEG (which is ultimately expected to include both the initial participants and other universities and national laboratories that have been identified) have the expertise and experimental capabilities needed to both obtain and compile existing data archives and perform additional seismic and flooding experiments. The data developed by EHEG will be stored in databases within RISMC. These databases will be used to validate the advanced external hazard tools and methods. The numerical tools under development in RISMC are used to evaluate the impacts that external hazards have on nuclear power plants. These analysis tools will be used to advanced current risk calculation processes and reduce uncertainty in these calculations. To have confidence in the predictive capability of these numerical tools it is important to verify and validate them. To verify and validate numerical tools the user must understand what physics is being used to represent the problem (i.e. constitutive models used to represent the soil and structure), the numerical solver approach used (i.e. finite element), what data already exists that could be used to validate the tools, and what data is needed for validation.

31274↗
DOE OSTI2023

Multi-Sourced Collaboration for the Production and Refining of Rare Elements and Critical Metals (Final Technical Report)

The project objective was to develop a feasible and cost-effective method for recovering rare earth elements (REEs) and critical materials (CMs) from coal and coal byproducts, resulting in high-purity individually separated REEs and CMs. The targeted REEs included Y, Pr, Nd, Gd, Dy, and Sm, with a purity of over 99.5%, while the CMs included Co, Mn, Ga, Sr, Li, Ni, Zn, and Ge, with a purity of over 90%. The project aimed to design a prototype facility capable of producing 1-3 tonnes/day of high-purity REO mixes. The work was divided into four designated circuits: 1) REE extraction and concentration, 2) REE separation and purification, 3) RE metal production, and 4) CM production. To achieve these goals, the project involved 11 tasks, including technology reviews, research, process flow diagram development, mass balance estimation, and preliminary technical-economic analysis. The project team included researchers from the University of Kentucky, University of Alabama and Virginia Tech as well as process specialists from Argonne National Laboratory. MP Materials provided technical support regarding rare earth markets and processing while Alliance Coal performed resource assessment. The project included a market analysis for Nd/Pr, Tb, Dy, Gd, Y, Co, Mn, Li, Sr, Ga, Ni, Zn, and Ge. These analyses provided insights into the supply and demand trends as well as historic and future projections of market price relative to purity requirements for these elements. Two coal resources were selected for the project: the West Kentucky No. 13 (Baker) Seam and an undisclosed lignite resource in the Illinois coal basin. The estimated quantities of REEs in these resources were calculated based on production samples and drilling data. It was estimated that there is adequate supply for an operation producing one metric ton daily of higher purity mixed rare earth oxides (MREO) for approximately 20 years at a site located in western Kentucky. In Circuit 1, project data was obtained from a pilot heap leach and REE concentration facility. It was concluded that the existing circuit, which generated a MREO concentrate, two types of CM mixed products, and Li- and Sr-containing waters, would be suitable feed for circuits 2-4. Data from the first-of-its-kind coal coarse refuse heap leach pilot pad played a crucial role in estimating reliable elemental concentrations of the pregnant leaching solution (PLS). The average total REE concentration in the PLS was found to be 28.6 ppm. In Circuit 2, several concepts were explored including a novel process referred to as solvent-assisted chromatography (SAC). This concept involved a novel columnar reactor that incorporated multiple mixer/settlers, thereby enabling the operation of counter-flowing aqueous and organic phases. Unfortunately, due to project time constraints, a complete fundamental modeling analysis could not be completed to fully evaluate the technology. Molten salt electrowinning was considered as an alternative for circuit 3 following circuit 2 purification circuit utilizing the novel SAC process. A mass and energy balance of Nd reduction to metal in a fluoride containing molten salt electrolyte was conducted. Comparisons were made with the current state of Asian molten salt electrorefining, and potential improvements in siphoning rare earth metals (REM) from the reactor were presented. A cost estimate was performed for the production of 1 tonne per day, which yielded a total of $2.29 million for the nine electrowinning (EW) cells required. The selected option for circuits 2 and 3 was a plasma distillation process, which initially separates rare earth elements (REEs) from other elements. This is followed by selective electrowinning in various ionic liquids. The selection was made on the basis of thermodynamic modeling and experimental data previously published by a project partner. The combination offers an innovative approach to integrated refining and RE metal production. For Circuit 4, an extensive literature review was conducted for the processing of the CMs. The ultimate decision was to utilize a combined plasma and ionic liquid process as well to produce individual high-purity concentrates of Zn, Ni, Co, Mn, and Mg. A separate flowsheet for Li and Sr was recommended, which would yield carbonates of these elements. Due to the lack of suitable experimental data at this time, a process recommendation could not be provided but several methods have been proposed for consideration. Lastly, a techno-economic analysis (TEA) was conducted to assess the effectiveness of the proposed process for further investigation. The TEA results revealed a capital expense (CapEx) of $737 million and an annual operational expense (OpEx) of $220 million. Due to the selected elements, the hypothetical heap leach pad can produce 1 metric tonne per day of REO equivalent, but a conscious decision was made to only treat targeted REEs, resulting in the production of 0.4 metric tonne of REM. An estimated annual revenue of $90.87 million was projected based on standard market pricing information provided by the funding agency. During the TEA, ten different modules were evaluated for costing purposes. The precipitation circuit was identified as the largest single operational expense, followed by the Mg/Mn process due to the amount of treated metal. In terms of capital expenditures, the heap leach process incurred the highest cost, followed by the Mg/Mn process. The scalability of the plasma process is a crucial consideration since the reactors cannot be scaled beyond the largest demonstrated size due to their reliance on surface area of the slag and vapor phase. The purity estimate for the REEs are generally 98%±2% to produce a metal. The purity level being lower than the project objective was due to the lack of specific experimental data needed to tighten the tolerance of the estimates. Based on literature and previous experience, the CMs are estimated as follows; Ga (95%+, metal), Sr (95%+, carbonate), Li (95%+, carbonate), Ni (98%±2%, metal), Zn (95%+, metal sponge), Ge (95%+, metal), Co (98%±2%, metal), and Mn (98%±2%, metal).

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Q1 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

This report describes the work and activities carried out towards the completion of each of the following milestones in FY25 Q1: 1. Demonstrate workflow for generating self-consistent CESOL plasma profiles + first wall and divertor loading prediction and generate the CAT plasma and neutron loading needed for further engineering analysis: $/circ$ Initially run CESOL with SOLPS and map heat flux using simple HEAT method: ▪ Generate SOLPS grid for CAT reference case, ▪ Implement and apply ‘lower’ fidelity method of using HEAT-like analytic method to map SOLPS heat flux from charged particles to first wall, and ▪ Provide initial heat flux + neutron loading for further engineering analysis. 2. Demonstrate multiphysics magnet analysis: $/circ$ Demonstrate validation of the Elmer workflow by comparing induced stresses due to EM forces generated by TF coils using multiple codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Statistical effects related to low numbers of reacting molecules analyzed for a reversible association reaction A + B = C in ideally dispersed systems: An apparent violation of the law of mass action

Theoretical analysis and computer simulations (Monte Carlo and numerical integration of differential equations) show that the statistical effect of a small number of reacting molecules depends on a way the molecules are distributed among the small volume nano-reactors (droplets in this study). A simple reversible association A + B = C was chosen as a model reaction, enabling to observe both thermodynamic (apparent equilibrium constant) and kinetic effects of a small number of reactant molecules. When substrates are distributed uniformly among droplets, all containing the same equal number of substrate molecules, the apparent equilibrium constant of the association is higher than the chemical one (observed in a macroscopic—large volume system). The average rate of the association, being initially independent of the numbers of molecules, becomes (at higher conversions) higher than that in a macroscopic system: the lower the number of substrate molecules in a droplet, the higher is the rate. This results in the correspondingly higher apparent equilibrium constant. A quite opposite behavior is observed when reactant molecules are distributed randomly among droplets: the apparent association rate and equilibrium constants are lower than those observed in large volume systems, being the lower, the lower is the average number of reacting molecules in a droplet. The random distribution of reactant molecules corresponds to ideal (equal sizes of droplets) dispersing of a reaction mixture. Our simulations have shown that when the equilibrated large volume system is dispersed, the resulting droplet system is already at equilibrium and no changes of proportions of droplets differing in reactant compositions can be observed upon prolongation of the reaction time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Molecular dynamics simulations of reflection and sputtering behavior of boron under deuterium ion irradiation

Boronization is a commonly used method of wall conditioning in fusion reactors. The application of boron films to the plasma-facing materials results in enhanced plasma performance due to the reduction of intrinsic impurities. This is primarily driven by a reduction in oxygen content that is chemically trapped in the boron film. The reactive nature of these boron films also raises questions concerning interactions with hydrogen isotopes. In this work, boron-deuterium interactions were studied using molecular dynamics (MD). Reactive force field potentials were used to model the chemical interactions between B and D. An amorphous boron substrate was irradiated by D atoms at varying incident energies, 10 eV < E i < 150 eV and angles, 0° < α < 85°. The reflection probability was calculated and compared to results from the commonly-used binary collision approximation (BCA) method. This comparison found that the BCA underestimated the reflection probability at E i < 35 eV and α > 45 ∘ . The source of this discrepancy was found to be the surface binding energy model. The BCA calculation with an isotropic surface binding energy model was more closely aligned to the MD result. This, in combination with a correction function based on the MD results allows for corrections to the reflection probability of deuterium impinging on boron surfaces. The sputtering of the substrate material was also studied. While this study did not contain sufficient events to quantitatively describe the sputtering behavior, some qualitative results emerged: namely, chemical sputtering of B and D-containing molecules (BD, BD 2 , BD 3 ) at low ( < 20 eV) incident deuterium energies. This result suggests that chemical sputtering could be a significant factor in limiting boron coating lifetime when exposed to lower ion energies, such as those in detached plasmas. The results show that chemical interactions should be taken into account when modeling the interactions between D ions and B surfaces.

Boron↗
DOE OSTI2024

Photoinitiated Reactions of Molecules and Radicals in Molecular Beams

The UV photochemistry of organic molecules is a fundamental process that governs reactions in the atmosphere, synthetic chemistry, processing of organic aerosols, and biological damage in living tissues. In many organic molecules, the ensuing evolution involves pathways that are in competition, giving rise to different products, isomerization, coupling to other electronic states, and secondary reactive collisions. Because photodissociation is usually fast (picoseconds to microseconds), these processes are far from equilibrium and are controlled by kinetic competition and dynamical forces. The work accomplished was focused primarily on the photochemistry of alpha-keto carboxylic acids, a group of acids produced from natural sources, which are implicated in aerosol formation and biological processes. In the atmosphere, they are destroyed mainly by solar radiation. Studying their photochemistry has been surprisingly difficult because of the complexity of their excited electronic states, and the effect of collisions and secondary reactions. The second project investigated the lifetime of excited electronic states of pyrazine and picoline and had both experimental and theoretical aspects. The work focused on a model molecule, pyruvic acid (PA), because it was hypothesized that it should be a good source of the unstable carbene, methylhydroxycarbene (MHC). PA has an internal hydrogen bond that controls the evolution of its decomposition. The decomposition was studied following excitation to two of its lowest excited states reached by laser irradiation at 351 and 193 nm. The photodissociation dynamics in the absence and presence of collisions was monitored and compared. Understanding the UV photochemistry requires the use of complementary experimental approaches. Two methods were enlisted, which together generated a comprehensive and detailed set of results: (i) The time-sliced velocity map imaging (SVMI) instrument at USC was exploited to determine kinetic energy release of fragments, fast dissociation timescales, and internal state distributions of fragments for which Resonance Enhanced Multiphoton Ionization (REMPI) schemes exist; and (ii) The multiplexed photoionization mass spectrometer (MPIMS) setup developed at the Sandia Combustion Research Facility was used for product discovery, achieved by exploiting tunable narrowband VUV radiation at the Advanced Light Source (ALS), and to follow in real time their subsequent unimolecular and bimolecular reactions. The experimental conditions ranged from collisionless molecular beams to study nascent products to flow reactors of variable pressures to study subsequent bimolecular reactions of the products. The goals stated above were successfully accomplished. By exciting PA to its lowest excited state, MHC was identified as the only primary product and its isomerization to vinylalcohol and acetaldehyde was directly observed in real time. Moreover, we reported the first bimolecular reaction of MHC, which was with acetaldehyde, and identified its reaction product. This paves the way for the study of other reactions of this important carbene intermediate. When excitation was carried out at 193 nm, which imparted much higher energy to PA, many more products were directly observed in real time, some deriving from three-body dissociation. Quantitative branching ratios and secondary reactions of radical products were also determined and analyzed. The studies of pyrazine and picoline focused on the second excited state of these molecules and showed (via ionization studies) that their excited states were very short lived (<100 fs) because of efficient couplings to lower electronic states via vibronic coupling. Theoretical work on modeling the absorption spectrum and decay mechanisms of the excited states are in progress in collaboration with theoreticians.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Report on Data Quality Required for Modeling an MSR

Calculations were performed to assess how uncertainties in molten salt property values impact the results of behavior models of molten salt reactor systems under steady-state conditions. The thermophysical and thermochemical properties of salt mixtures being considered for use in molten salt reactors are not yet well characterized and methods used to measure property values to molten salts are not yet well established. The limited property data for salts of interest that are available in the literature are inconsistent and unreliable. Many methods used to measure property values remain developmental in the sense that not all variables affecting measured values have been identified or controlled during the measurement and not all aspects of the measurements affecting data quality are calibrated. A prior report described the repeatability of measurements made at Argonne and the reproducibility of measurements made with the same salts at several institutions is being evaluated. There is significant uncertainty in measured property values due to instrumental limitations and the effects of contaminants in the salts used to make the measurements that impacts predictions of salt behavior in a reactor. Most salts of interest are multicomponent eutectic or near-eutectic mixtures and both uncertainty in the salt composition and uncertainty in measured values. Furthermore, the salt composition evolves over time due to the uptake of impurities from the atmosphere and corrosion of containment vessels. More significantly, fission products will accumulate during operation that change the salt compositions and property values. MSR developers need to understand how this uncertainty in the basic property data of the fuel will impact modeling predictions for normal reactor operation as well was in transient scenarios that will impact the safety basis. The calculations and analyses in this report address the impact of uncertainty in measured property values under steady-state conditions through a generalized parametric study. The pooled uncertainty from contributions by a number of factors were considered and quantified, including varying composition due to impurities and fission products, imprecision in the measured values, and uncontrolled aspects of the measurement. In this way, we can begin to gain an understanding of the magnitude of the impact uncertainty will have in modeling MSR systems.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Deuterium retention characteristics in Li film by coating and during flowing liquid Li limiter operation in experimental advanced superconducting tokamak

Lithium (Li) is a promising low-Z material for particle recycling and impurity control to improve plasma performance in fusion devices. In the experimental advanced superconducting tokamak (EAST), Li coating has become a routine method for wall conditioning, and a flowing liquid Li (FLiLi) limiter has been successfully tested several times. Deuterium retention characteristics in the Li film coated on the international thermonuclear experimental reactor-like tungsten divertor and FLiLi during plasma discharges, which is important for the utilization of Li in future fusion devices, were investigated in EAST. It is found that the absorption of the fuel particles by Li coatings decreases gradually, and recycling gradually increases over a series of discharges. The maximum net amount of deuterium retained reached ~0.8 g, corresponding to 12% deuterium in the Li. Furthermore, this corresponds to a whole day's worth of shots with a total of 87 plasma discharges and a total of ~640s plasma time after 11.75 g Li was deposited. Compared to the Li coating, it is shown that FLiLi continuously traps fuel particles and achieves a higher deuterium retention ratio over both the short- and long-term, leading to lower recycling. Meanwhile, it is also observed that the fuel particle retention ratio increases when FLiLi is closer to the plasma. This result is likely due to higher plasma heating power and limiter temperature, which cause an increased Li efflux from FLiLi due to a higher limiter temperature, which then redeposits on the other plasma-facing surfaces and increases fuel particle absorption. It is estimated that ≥ 80% of the retained D particles are captured by the continual renewal of the Li redeposition film during the FLiLi operation. This investigation would also be useful for D/T retention in future fusion devices if Li is used as a plasma-facing component.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗