Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “initio molecular dynamics”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Lower-length-scale modeling of chemical additions, corrosion and fission products in select MSR base salts

Molten salts play a crucial role in advancing Generation IV nuclear energy technology, with chloride-based salts like NaCl-UCl 3 garnering significant attention due to their distinctive properties. However, the corrosive nature of molten salts causes the dissolution of chromium (Cr), leading to the formation of CrCl 2 and CrCl 3 species in molten chlorides. Moreover, the radioactive decay of nuclear fuel gives rise to fission products, including Cs, Sr, and I chlorides. The first part of this report presents a comprehensive study utilizing ab initio molecular dynamics (AIMD) simulations to investigate the properties of eutectic NaCl-UCl 3 molten salt in the presence of corrosion products (CrCl 2 and CrCl 3 ) and fission product (CsCl). The study focuses on essential structural and thermophysical properties such as density, mixing energy, coordination numbers (CN), and Radial Distribution Functions (RDF) of the salts with varying compositions of corrosion products and fission product, covering a range from 0% to 13.5%. The results offer valuable insights into the behavior of corrosion and fission products in uranium-based molten salts, providing essential data that can be used as input to the MTDB-TC and MSTDB-TP property databases being developed by the NEAMS program. Due to their favorable characteristics such as low melting points, high boiling points, and low costs, MgCl 2 +NaCl+KCl (MNK) eutectic salts have recently attracted attention as high-temperature heat transfer fluids. The incorporation of LiCl into MNK salts can further reduce their melting points and increase their specific heat capacities, which is de- sirable for high-temperature heat transfer applications. The second part of this report presents the development and validation of a new shell-model potential for the MgCl 2 +NaCl+KCl+LiCl system, which captures the polarization of Cl anions. The extensive comparison with experimental data and AIMD simulations demonstrates the accuracy and reliability of the potential. Furthermore, using this potential, we elucidate the intricate network structure in MNK eutectic salts that is formed through polyhedron sharing. This research contributes to a fundamental understanding of the atomic structures and thermophysical properties of multi-component molten salts, which is crucial for future development of heat transfer fluids for applications in molten salt reactors.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Design of CO 2 Selective Type 3 Porous Liquids Through Porous Host Morphology

Direct air capture (DAC) of CO 2 is a negative emission technology under development to limit the impacts of climate change. The dilute concentration of CO 2 in the atmosphere (~400 ppm) requires new materials for carbon capture with increased CO 2 selectivity that is not met with current carbon capture materials. Porous liquids (PLs) are an emerging candidate for carbon capture and consists of a combination of solvents and porous hosts that creates a liquid with permanent porosity. The fundamental mechanisms of carbon capture in a PL are relatively unknown. To uncover these mechanisms, PLs were synthesized consisting of three different zeolitic-imidazolate framework (ZIF-8, ZIF-67, or ZIF-69) porous host in a water/glycol/2-methylimidazole solvent. The most stable composition was based on ZIF-8 and exhibited carbon capture following exposure to CO 2 . Density functional theory identified a three-step carbon capture mechanism based on (i) reaction of OH- with ethylene glycol in the solution followed by (ii) formation of 2-hydroxyethyl carbonate, which (iii) further react with OH- to form a carbonate species. This mechanism was validated with experimental nuclear magnetic resonance spectroscopy (NMR) to identify the dissolved carbonate phases and the decrease in the pH during CO 2 exposure. Deuterated samples of the ZIF-8 PLs were synthesized and analyzed via neutron diffraction at the Spallation Neutron Sources at Oak Ridge National Laboratory. Results identified differences in diffraction for PLs pre- and post-CO 2 exposure that will be combined with ab initio molecular dynamics data of the same PL composition to identify how the presence of a solvent-porous host interfaces results in carbon capture.

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

FY24 Progress Report on Viscosity and Thermal Conductivity Measurements of Nuclear Industry Relevant Chloride Salts: An Experimental and Computational Study

As presented in this report, experimental and computational techniques were performed to assess the viscosity and thermal conductivity of key alkali and actinide chloride mixtures for molten salt reactor developers. These mixtures were pure LiCl, NaCl-KCl, LiCl-NaCl, LiCl-KCl, LiCl-NaCl-KCl, and NaCl-UCl 3 . Experimental measurements of viscosity were performed with a rolling ball viscometer, whereas experimental measurements of thermal conductivity were performed with a variable gap apparatus. Additional benchmarking work was performed using both property measurement systems to prepare for x-ray radiography in stainless-steel crucibles for viscosity and to ensure that calibration methods were accurate for thermal conductivity before assessing the NaCl-UCl 3 system. Validation data for the NaCl-UCl 3 in literature are minimal. Details on the calibration methods, salt measurement processes, and sources of error and uncertainty are discussed in detail for both property measurements. The computational methods described herein involved ab-initio molecular dynamics (AIMD) calculations using CP2K. The calculations were performed for the LiCl-KCl-NaCl and NaCl-UCl 3 systems. These calculations not only provided thermophysical property estimations for comparison to experimental data, but they also allowed for the determination of diffusion coefficients, coordination numbers, and radial distribution functions to provide insight into ion mobility and local coordination environments, which is linked to macroscopic property trends.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

An Overview of the Molten Salt Thermal Properties Database–Thermophysical, Version 4.0 (MSTDB-TP V.4.0)

A central repository of thermophysical and thermochemical properties of molten salt compositions of relevance to molten salt reactors (MSRs) is vital in supporting the broad community of MSR developers, who are at various stages of developing and deploying their reactor designs. In general, these MSR designs differ significantly from developer to developer (e.g., with respect to the hardness of the neutron spectra, level of fissile loading, target multicomponent temperatures and power levels, and moderating capabilities). Therefore, the fuel and coolant salts being considered vary greatly: they may be chlorides or fluorides, they utilize different actinides at different ratios, and the cations in the melt are selected based on perceived advantages and disadvantages. Considering the general need for thermal properties, and the vastness of the array of potential candidate salt mixtures, the Molten Salt Thermal Properties Database (MSTDB) was initiated in 2018 with the goal of providing thermophysical and thermochemical characterization of key molten salt compounds and mixtures across their temperature and compositional domains. The MSTDB is thus divided into the thermophysical arm (MSTDB-TP) and the thermochemical arm (MSTDB-TC). The MSTDB is an effort funded by the Department of Energy, Office of Nuclear Energy (DOE-NE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, and the MSR Campaign. This report provides an overview of the MSTDB-TP v4.0 in terms of the data contained within, the state of the tools used to access the data, the availability of predictive models that leverage the raw data in the database, the preliminary status of developmental efforts that are currently underway, and an account of future goals for MSTDB-TP. The primary goal for the update from MSTDB-TP v.3.1 to v4.0 was the incorporation of surface tension data into the database; this property is important for thermal hydraulics modeling and species transport in other tools that have been developed under the NEAMS program. A breakdown of the surface tension data that have been added into MSTDB-TP v4.0 is provided herein, and the manner in which the quality of the data has been assessed is also documented. For MSTDB-TP v4.0, newly published thermophysical property data—primarily from collaborative experimental efforts under the MSR Campaign—have been incorporated into the database, and the resulting expansion is documented here. Because of the size to which MSTDB-TP has grown, the raw data format has now been recast into JavaScript Object Notation (JSON) format for easier connection with the MSTDB-TP application programming interface (API). Saline; the pre-existing comma-separated value (CSV) format has been deprecated but is still maintained, accessible, and up to date. As a final effort in packaging the MSTDB-TP v4.0 update, the graphical user interface (GUI) for MSTDB has been updated to allow full accessibility to the density and viscosity predictive models, which are based on Redlich-Kister expansions of MSTDB-TP raw data. Some other major aspects of this report, in terms of preliminary and future work, include: (1) documentation of the formalism and preliminary testing of a kinetic theory model that may act as a predictive model for thermal conductivity; (2) documentation of the candidate predictive models that may be considered in the future for surface tension, making use of the surface tension data now in MSTDB-TP v4.0; (3) a preliminary account of a data collection process that will enable the filling of additional gaps within MSTDB-TP, namely with data which have been collected computationally (e.g., through ab initio molecular dynamics).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Atomistic Simulations for Thermophysical Properties of Uranium-Containing Halide Molten Salts

Characterizing the thermophysical properties in both fuel and coolant salts are critical in modeling, developing, process optimizing and utilizing molten salt reactors (MSRs), as these properties directly relate to operation metrics and can inform on the selection of candidate salts. The demand for consistent, accurate and publicly available thermophysical property data has become more apparent in recent years as interests have increased from molten salt reactor developers. There are a number of challenges in experimentally measuring properties such as thermal conductivity, viscosity, density and heat capacity , which have led to sparse and often times conflicting data points or molten salts in general. Additionally, there are a number of hazards to consider when synthesizing, storing, using, treating and disposing of molten salts. With the advances in computational capabilities over the last 10 years, the use of atomistic simulations can be implemented to support these efforts. The primary objective of this work is characterize the thermophysical transport properties in a number of molten chloride salts, and in particular NaCl-UCl 3 using ab-initio molecular dynamic (AIMD) simulations. In this binary salt the UCl 3 acts as the primary fissile material and NaCl acts as a carrier salt due with its’ high solubility for actinides A number of studies on the thermophysical properties of NaCl-UCl 3 have been published but there is not a vast amount of viscosity data for this system. In 1975, Desyatnik, et al published a study reporting dynamic viscosities that were calculated from kinematic viscosity measurements, and using the coefficients provided the viscosity in a 70:30 NaCl:UCl 3 mixture is 2.29 cP and 2.88 for a 60:40 mixture. Termini et al. recently reported viscosities in the range of 2.75 – 3 cP for the 63:37 NaCl-UCl 3 mixture in the same temperature range using rolling ball viscosity measurements. Computational viscosity of a similar mixture (64:36) can be obtained from the work Andersson et al. using the reported diffusion coefficients, and the hydrodynamic radius from the pair-radial distribution functions (RDFs). Using Eq (1) (vida infra), the viscosity would be 2.50 cP at 1100K. This is not to say that these values are incorrect due to the varying reported values, but aims to highlight the necessity of this work. The data reported in this ongoing work are computations on a 64:36 mixture of NaCl-UCl 3 at 987K. This work is likely to be expanded into varying concentrations of this mixture along with the inclusion of other salt candidate mixtures.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Development of a High-Rate Lithium-Air Battery Using a Gaseous CO 2 Reactant

Li-air batteries are considered a potential alternative to Li-ion batteries for transportation applications due to their high theoretical specific energy. Most works in this area focus on use of O 2 as the reactant. However, newer concepts for using gaseous reactants (such as CO 2 , which has a theoretical specific energy density of 1,876 Wh/kg) provide opportunities for further exploration. The main objective of this project was the development of a novel strategy that enables operation of Li-CO 2 batteries at high-capacity and high-rate, with a long-cycle-life. The team was able to: (1) Synthesize two novel transition metal chalcogenide (TMC) catalysts that work in synergy with ionic liquid-based electrolytes to enhance the efficiency of reactions during discharge and charge processes; (2) Fabricate high-porosity cathode electrodes with 3D printing to increase electrode surface area and gas permeability; (3) Develop a multiscale modeling framework that integrates Density Functional Theory (DFT), Ab-Initio Molecular Dynamics (AIMD), classical Molecular Dynamics (MD), and Finite Element Analysis (FEA) to investigate atomic and cell-level properties of Li-CO 2 batteries; (4) Assemble a stackable Li-CO 2 pouch-cell able to deliver a capacity of >200 mAh. These achievements were realized through an integrated approach based on materials synthesis, testing, characterization, analysis, and computation. This project produced a thorough understanding of key chemical, electronic, and kinetic parameters that govern the operation of Li- CO 2 batteries in realistic conditions. The methodologies employed, and the insight generated, will be valuable beyond advancing the field of Li-CO 2 batteries

25 ENERGY STORAGE↗
DOE OSTI2026

Hydrated Metal and Metal-Nitrate Complexes in Water: Full Lanthanide(III) Series plus Miscellaneous Metal Ions

This is a dataset of hydrated metal complexes and metal–nitrate hydrated complexes intended for public use, reproducibility, and downstream structural analysis. A key feature is coverage across the full lanthanide(III) series (La–Lu), enabling systematic comparisons of coordination motifs and bonding trends across the entire lanthanide sequence. In addition to the lanthanides, the dataset also includes other metal ions such as UO2(VI), Fe(II), and Fe(III). The dataset provides optimized geometries for hydrated and nitrate-containing hydrated complexes, together with representative ab initio molecular dynamics (AIMD) trajectories saved in standard XYZ formats. The accompanying NWChem input decks enable reproduction of the reported calculations and provide a starting point for extending the simulations to related coordination environments. Computationally, DFT calculations employ the B3LYP functional with DFT-D3BJ dispersion corrections and a COSMO continuum solvent model (dielectric constant 78.4) to represent solvation beyond the explicitly treated first hydration shell. AIMD simulations are performed with the NWChem qmd module at 298 K, integrating nuclear motion with the velocity-Verlet algorithm and controlling temperature using a Nosé–Hoover thermostat. Trajectories are approximately 4.8 ps in length and are used primarily to assess short-time stability of candidate coordination motifs, including (for lanthanides) differences between 8- versus 9-water coordination and comparisons between nitrate-bound and nitrate-free hydrated complexes.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2024

Functionalization of Cathode–Electrolyte Interface with Ionic Liquids for High-Performance Quasi-Solid-State Lithium–Sulfur Batteries: A Low-Sulfur Loading Study

We introduce a quasi-solid-state electrolyte lithium-sulfur (Li–S) battery (QSSEB) based on a novel Li-argyrodite solid-state electrolyte (SSE), Super P–Sulfur cathode, and Li-anode. The cathode was prepared using a water-based carboxymethyl cellulose (CMC) solution and styrene butadiene rubber (SBR) as the binder while Li6PS5F0.5Cl0.5 SSE was synthesized using a solvent-based process, via the introduction of LiF into the argyrodite crystal structure, which enhances both the ionic conductivity and interface-stabilizing properties of the SSE. Ionic liquids (IL) were prepared using lithium bis(trifluoromethyl sulfonyl)imide (LiTFSI) as the salt, with pre-mixed pyrrolidinium bis(trifluoromethyl sulfonyl)imide (PYR) as solvent and 1,3-dioxolane (DOL) as diluent, and they were used to wet the SSE–electrode interfaces. The effect of IL dilution, the co-solvent amount, the LiTFSI concentration, the C rate at which the batteries are tested and the effect of the introduction of SSE in the cathode, were systematically studied and optimized to develop a QSSEB with higher capacity retention and cyclability. Interfacial reactions occurring at the cathode–SSE interface during cycling were also investigated using electrochemical impedance spectroscopy, cyclic voltammetry, and X-ray photoelectron spectroscopy supported by ab initio molecular dynamics simulations. This work offers a new insight into the intimate interfacial contacts between the SSE and carbon–sulfur cathodes, which are critical for improving the electrochemical performance of quasi-solid-state lithium–sulfur batteries.

25 ENERGY STORAGE↗
DOE OSTI2021

Free Energies of Hydrated Halide Anions: High Through-Put Computations on Clusters to Treat Rough Energy-Landscapes

With a longer-term goal of addressing the comparative behavior of the aqueous halides F−, Cl−, Br−, and I− on the basis of quasi-chemical theory (QCT), here we study structures and free energies of hydration clusters for those anions. We confirm that energetically optimal (H2O)nX clusters, with X = Cl−, Br−, and I−, exhibit surface hydration structures. Computed free energies, based on optimized surface hydration structures utilizing a harmonic approximation, typically (but not always) disagree with experimental free energies. To remedy the harmonic approximation, we utilize single-point electronic structure calculations on cluster geometries sampled from an AIMD (ab initio molecular dynamics) simulation stream. This rough-landscape procedure is broadly satisfactory and suggests unfavorable ligand crowding as the physical effect addressed. Nevertheless, this procedure can break down when n≳4, with the characteristic discrepancy resulting from a relaxed definition of clustering in the identification of (H2O)nX clusters, including ramified structures natural in physical cluster theories. With ramified structures, the central equation for the present rough-landscape approach can acquire some inconsistency. Extension of these physical cluster theories in the direction of QCT should remedy that issue, and should be the next step in this research direction.

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

Mechanochemical Synthesis of Sustainable Ternary and Quaternary Nanostructured Cu2SnS3, Cu2ZnSnS4, and Cu2ZnSnSe4 Chalcogenides for Thermoelectric Applications

Copper-based chalcogenides have emerged as promising thermoelectric materials due to their high thermoelectric performance, tunable transport properties, earth abundance and low toxicity. We have presented an overview of experimental results and first-principal calculations investigating the thermoelectric properties of various polymorphs of Cu2SnS3 (CTS), Cu2ZnSnS4 (CZTS), and Cu2ZnSnSe4 (CZTSe) synthesized by high-energy reactive mechanical alloying (ball milling). Of particular interest are the disordered polymorphs of these materials, which exhibit phonon-glass–electron-crystal behavior—a decoupling of electron and phonon transport properties. The interplay of cationic disorder and nanostructuring leads to ultra-low thermal conductivities while enhancing electronic transport. These beneficial transport properties are the consequence of a plethora of features, including trap states, anharmonicity, rattling, and conductive surface states, both topologically trivial and non-trivial. Based on experimental results and computational methods, this report aims to elucidate the details of the electronic and lattice transport properties, thereby confirming that the higher thermoelectric (TE) performance of disordered polymorphs is essentially due to their complex crystallographic structures. In addition, we have presented synchrotron X-ray diffraction (SR-XRD) measurements and ab initio molecular dynamics (AIMD) simulations of the root-mean-square displacement (RMSD) in these materials, confirming anharmonicity and bond inhomogeneity for disordered polymorphs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Violation of the T –1 Relationship in the Lattice Thermal Conductivity of Mg 3 Sb 2 with Locally Asymmetric Vibrations

Most crystalline materials follow the guidelines of T –1 temperature-dependent lattice thermal conductivity (κ L ) at elevated temperatures. Here, we observe a weak temperature dependence of κ L in Mg 3 Sb 2 , T –0.48 from theory and T –0.57 from measurements, based on a comprehensive study combining ab initio molecular dynamics calculations and experimental measurements on single crystal Mg 3 Sb 2 . These results can be understood in terms of the so-called “phonon renormalization” effects due to the strong temperature dependence of the interatomic force constants (IFCs). The increasing temperature leads to the frequency upshifting for those low-frequency phonons dominating heat transport, and more importantly, the phonon-phonon interactions are weakened. In-depth analysis reveals that the phenomenon is closely related to the temperature-induced asymmetric movements of Mg atoms within MgSb 4 tetrahedron. With increasing temperature, these Mg atoms tend to locate at the areas with relatively low force in the force profile, leading to reduced effective 3 rd -order IFCs. The locally asymmetrical atomic movements at elevated temperatures can be further treated as an indicator of temperature-induced variations of IFCs and thus relatively strong phonon renormalization. The present work sheds light on the fundamental origins of anomalous temperature dependence of κ L in thermoelectrics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Suppressed paramagnetism in amorphous Ta 2 O 5 − x oxides and its link to superconducting-qubit performance

Amorphous-oxide layers in thin-film capacitors are linked to reduced transmon-qubit T 1 coherence times. Ta -based capacitors outperform Nb -based ones, suggesting that amorphous Ta 2 O 5 − x is less lossy than Nb 2 O 5 − x . We investigate the microscopic features of these amorphous oxides using ab initio molecular dynamics and density functional theory, revealing the origins of the superior performance of Ta 2 O 5 − x . We establish that oxygen deficiency is less likely to occur in amorphous Ta 2 O 5 − x than in Nb 2 O 5 − x for 0 ≤ x ≤ 0.25 and that for a given oxygen deficiency x , metal Ta — Ta bond formation is enhanced. Such bonds, which are accommodated by structural flaws in the amorphous network, capture electrons better than in amorphous Nb 2 O 5 − x . These thermochemical differences quench or highly suppress magnetic moments in amorphous Ta 2 O 5 − x and eliminate a potential source of quasiparticles and magnetic flux noise. We also show that hyperfine couplings between Nb nuclei and local magnetic moments in Nb 2 O 5 − x can form “two-level systems” (TLSs) or “two-level fluctuators” with energy splittings of 100–1000 MHz or higher. This reveals a TLS mechanism in amorphous Nb 2 O 5 − x oxide layers that is likely inactive in Ta 2 O 5 − x . Our work provides a fundamental understanding of the materials chemistry and limitations imposed by native oxides of superconducting qubits that can be used to guide materials selection and processing.

Pritchard, P. Graham [Northwestern U.] (ORCID:0000↗
NASA NTRS2016

Li-Doped Ionic Liquid Electrolytes: From Bulk Phase to Interfacial Behavior

Ionic liquids have been proposed as candidate electrolytes for high-energy density, rechargeable batteries. We present an extensive computational analysis supported by experimental comparisons of the bulk and interfacial properties of a representative set of these electrolytes as a function of Li-salt doping. We begin by investigating the bulk electrolyte using quantum chemistry and ab initio molecular dynamics to elucidate the solvation structure of Li(+). MD simulations using the polarizable force field of Borodin and coworkers were then performed, from which we obtain an array of thermodynamic and transport properties. Excellent agreement is found with experiments for diffusion, ionic conductivity, and viscosity. Combining MD simulations with electronic structure computations, we computed the electrochemical window of the electrolytes across a range of Li(+)-doping levels and comment on the role of the liquid environment. Finally, we performed a suite of simulations of these Li-doped electrolytes at ideal electrified interfaces to evaluate the differential capacitance and the equilibrium Li(+) distribution in the double layer. The magnitude of differential capacitance is in good agreement with our experiments and exhibits the characteristic camel-shaped profile. In addition, the simulations reveal Li(+) to be highly localized to the second molecular layer of the double layer, which is supported by additional computations that find this layer to be a free energy minimum with respect to Li(+) translation.

molecular dynamics↗
NASA NTRS2020

Molecular Simulations of Ultrafast Radiation Induced Melting at Metal-Semiconductor Interfaces

Understanding radiation induced ultrafast melting at material interfaces is essential in designing robust electronic devices for aviation/space applications and in laser machining. While it is difficult to achieve the spatial and temporal resolution required to quantify the phenomenon experimentally, simulations can provide the detailed mechanisms of the structural changes that happen during phase transition. In this work, we use molecular simulations to study the effect of radiation damage on silicon carbide (SiC) - tungsten (W) interfaces which is of interest in high power electronics. A multi-scale approach is involved wherein the reactions at the interfaces are quantified using ab-initio molecular dynamics (MD) simulations and classical MD simulations are employed to understand the structural and diffusional changes across the material interface. Finally, coarse-grained Lennard-Jones type models are used to study the larger scale mechanisms and structures obtained due to the induced damages. We show that the response of the material to radiation damage depends on factors such as energy of the incident radiation, thermal properties, and molecular structure of the material.

Ravichandran, Ashwin↗
DOE OSTI2022

Synthesis and Thermophysical Property Determination of NaCl-PuCl3 Salts

Currently, a knowledge gap exists in the available data and understanding of thermophysical properties relating to fresh fuel salts, especially those containing plutonium. Determination of this data is necessary for the design and construction of test reactors, as well as the licensing of future commercial molten salt reactors. Thermophysical properties such as melting temperature, salt stability, density, and heat capacity were determined on synthesized eutectic NaCl-PuCl3 and a more sodium rich composition containing 25 mol% PuCl3. These measurements document the baseline properties of the salt as a function of temperature for future experiments on irradiated fuel salt which will provide a holistic perspective on the change of thermophysical properties during reactor operations. It was determined that the NaCl-PuCl3 ingot synthesized for this study contained 63.4 mol% NaCl, 36.3 mol% PuCl3, and was 99.7% pure. Upon heating it was shown that the NaCl- PuCl3 eutectic was stable at temperatures of 800°C. The onset of melting occurred at 541°C and the enthalpy of fusion was determined to be 140.7 ± 8.4 J/g. Specific heat capacity measurements showed a slightly decreasing trend with respect to temperature in the liquid phase ranging from 0.67 to 0.57 J/g·K, with an average value of 0.637 ± 0.03 J/g·K (104 ± 5 J/mol·K) between 500 to 720°C. Three independent trials of the molten NaCl-PuCl3 eutectic salt found the density to be ?(T) = 3.8589 – 9.5342·10-4 T(°C), validated between 500 to 800°C. In addition to salt synthesis and experimentally determining thermodynamic properties, Ab Initio molecular dynamic (AIMD) simulations were used to calculate density and heat capacity values.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Demonstrate Improved Ag Diffusion and Describe the Basis for Pd Penetration Modeling in SiC

In past work, an effective diffusivity coefficient was determined for Ag transport through the silicon carbide layer of a tristructural isotropic fuel particle. The effective diffusivity coefficient accounts for the microstructure of the silicon carbide and includes both bulk diffusion and grain boundary diffusion of Ag. In this report, the model has been improved by accounting for the enhanced concentration of vacancies in the bulk due to irradiation, which substantially influence bulk diffusivity at low temperatures. To improve the BISON model and make it fission rate dependent, effective diffusivity calculations have been performed that incorporate the radiation modified bulk diffusivity. The microstructure and irradiation-dependent effective diffusivity has also been implemented into BISON, and its predictions for Ag release from tristructural isotropic fuel have been successfully compared to AGR-1 post irradiation measurements. Moreover, a new feature has been developed in the Multiphysics Object-Oriented Simulation Environment (MOOSE) to account for different grain boundary types. The Ag diffusivity in 5 (210)/[001] grain boundaries has been computed and was found to be greater than in random high-angle grain boundaries. The presence of the fission product Pd can also have an important effect on the properties of the silicon carbide layer in tristructural isotropic particles. The penetration of Pd into the silicon carbide layer causes a corrosion reaction that can lead to the failure of the silicon carbide layer; however, this corrosion reaction is not well understood. To enable an improved understanding of the mechanism, ab-initio molecular dynamics simulations of Pd interaction with bulk silicon carbide have been performed. The improved understanding of the reaction will form a basis for future improvements to the BISON’s Pd penetration failure model

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Predicting Hugoniot equation of state in erythritol with ab initio and reactive molecular dynamics

Erythritol has been proposed as an inert surrogate for developing theoretical and computational models to study aging in energetic materials. In this work, we present a comparison of mechanical and shock properties of erythritol computed using the ReaxFF reactive force field and from ab initio calculations employing density functional theory (DFT). We screened eight different ReaxFF parameterizations, of which the CHO parameters developed for hydrocarbon oxidation provide the most accurate predictions of mechanical properties and the crystal structure of erythritol. Further validation of the applicability of this ReaxFF parameterization for modeling erythritol is demonstrated by comparing predictions of the elastic constants, crystal structure, vibrational density of states, and Hugoniot curves against DFT calculations. The ReaxFF predictions are in close agreement with the DFT simulations for the elastic constants and shock Hugoniot when the crystal is loaded along its c axis but show as much as 30% disagreement in the elastic constants in the ab plane and 12% difference in shock pressures when shocked along the a or b crystal axes. Last, we compare thermomechanical properties predicted from classical molecular dynamics with those calculated using the quasi-harmonic approximation and show that quantum mechanical effects produce large discrepancies in the computed values of heat capacity and thermal expansion coefficients compared with classical assumptions. Combining classical molecular dynamics predictions of mechanical behavior with phonon-based calculations of thermal behaviors, we show that predicted shock-induced temperatures for pressures up to 6.5 GPa do not exceed the pressure-dependent melting point of erythritol.

Hu, Jing↗