Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “THERMOCHEMISTRY”

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 37 records · Page 2

Tree-based solvers for adaptive mesh refinement code $\scriptsize{FLASH}$ – IV. An X-ray radiation scheme to couple discrete and diffuse X-ray emission sources to the thermochemistry of the interstellar medium

X-ray radiation, in particular radiation between 0.1 and 10 keV, is evident from both point-like sources, such as compact objects and T-Tauri young stellar objects, and extended emission from hot, cooling gas, such as in supernova remnants. The X-ray radiation is absorbed by nearby gas, providing a source of both heating and ionization. While protoplanetary chemistry models now often include X-ray emission from the central young stellar object, simulations of star-forming regions have yet to include X-ray emission coupled to the chemo-dynamical evolution of the gas. We present an extension of the $\scriptsize{TREERAY}$ reverse ray trace algorithm implemented in the flash magnetohydrodynamic code which enables the inclusion of X-ray radiation from 0.1 keV < E γ < 100 keV, dubbed $\scriptsize{XRAYTHESPOT}$. $\scriptsize{XRAYTHESPOT}$ allows for the use of an arbitrary number of bins, minimum and maximum energies, and both temperature-independent and temperature-dependent user-defined cross-sections, along with the ability to include both point and extended diffuse emission and is coupled to the thermochemical evolution. We demonstrate the method with several multibin benchmarks testing the radiation transfer solution and coupling to the thermochemistry. Finally, we show two example star formation science cases for this module: X-ray emission from protostellar accretion irradiating an accretion disc and simulations of molecular clouds with active chemistry, radiation pressure, and protostellar radiation feedback from infrared to X-ray radiation.

79 ASTRONOMY AND ASTROPHYSICS↗

ELECTRONIC STRUCTURE METHODS AND PROTOCOLS WITH APPLICATION TO DYNAMICS, KINETICS AND THERMOCHEMISTRY

Hydrocarbon combustion involves the reaction dynamics of a tremendous number of species beginning with many-component fuel mixtures and proceeding via a complex system of intermediates to form primary and secondary products. Combustion conditions corresponding to new advanced engines and/or alternative fuels rely increasingly on autoignition and low-temperature-combustion chemistry. In these regimes various transient radical species such as HO2, ROO·, ·QOOH, HCO, NO2, HOCO, and Criegee intermediates play important roles in determining the detailed as well as more general dynamics. A clear understanding and accurate representation of these processes is needed for effective modeling. Given the difficulties associated with making reliable experimental measurements of these systems, computation can play an important role in developing these energy technologies. Accurate calculations have their own challenges since even within the simplest dynamical approximations such as transition state theory, the rates depend exponentially on critical barrier heights and these may be sensitive to the level of quantum chemistry. Moreover, it is well-known that in many cases it is necessary to go beyond statistical theories and consider the dynamics. Quantum tunneling, resonances, radiative transitions, and non-adiabatic effects governed by spin-orbit or derivative coupling can be determining factors in those dynamics. Building upon progress made during a period of prior support through the DOE Early Career Program, this project combines developments in the areas of potential energy surface (PES) fitting and multistate multireference quantum chemistry to allow spectroscopically and dynamically/kinetically accurate investigations of key molecular systems (such as those mentioned above), many of which are radicals with strong multireference character and have the possibility of multiple electronic states contributing to the observed dynamics. An ongoing area of investigation is to develop general strategies for robustly convergent electronic structure theory for global multichannel reactive surfaces including diabatization of energy and other relevant surfaces such as dipole transition. Combining advances in ab initio methods with automated interpolative PES fitting allows the construction of high-quality PESs (incorporating thousands of high-level data) to be done rapidly through parallel processing on high-performance computing (HPC) clusters. In addition, new methods and approaches to electronic structure theory will be developed and tested through applications. This project will explore limitations in traditional multireference calculations (e.g., MRCI) such as those imposed by internal contraction, lack of high-order correlation treatment and poor scaling. Methods such as DMRG-based extended active-space CASSCF and various Quantum Monte Carlo (QMC) methods will be applied (including VMC/DMC and FCIQMC). Insight into the relative significance of different orbital spaces and the robustness of application of these approaches on leadership class computing architectures will be gained. Synergy with other components of this research program such as automated PES fitting and multireference quantum chemistry will be used to address challenges encountered by the standard approaches to computational thermochemistry (those being single-reference quantum chemistry and perturbative treatments of the anharmonic vibrational energy, which break down for some cases of electronic structure or floppy strongly coupled vibrational modes).

74 ATOMIC AND MOLECULAR PHYSICS↗

Towards the Development of Molten Salt Reactor Multiphysics Modeling Capabilities: Depletion, Thermochemistry, Reactor Dynamics, and Species Tracking

This report documents the simulation and modeling work performed under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, utilizing the tools and codes developed under the program. The focus in this Fiscal Year (FY) is on the depletion, thermochemistry, reactor dynamics, and species transport in MSRs. The fission products travelling in the entire primary loop is unique in MSRs. There are many physical/chemical mechanisms that can significantly affect the fission products concentration in the fuel salt, e.g., diffusion into gas bubbles, diffusion in the core moderators, deposition on the surface of structures/pipes, etc. The fission product concentration in the fuel salt can potentially affect the neutronics, system dynamics, chemical potentials, etc., so accurately predicting the fission products distribution in the primary loop is important. Under the main objective, different activities are pursued to assess the NEAMS tools applicability for a variety of areas on species transport, including the depletion analyses using Griffin (fission product source), the redox potential assessment at different depletion points, the assessment of xenon transport effect on the neutronics and reactor dynamics, a survey of correlations for calculating the mass transfer coefficient of different fission products, etc.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

High-Fidelity Multiphysics Assessment of the Molten Chloride Reactor Experiment (MCRE) Integrating Neutronics, Thermal-hydraulics and Thermochemistry

This work presents a MOOSE-based Multiphysics model of the LOTUS MSR reactor system. The model includes the coarse-mesh turbulent thermal-hydraulics including passive advection of delayed neutron precursors of the primary system through Pronghorn, the power density and neutronics calculation employing Griffin, and the chemical interactions through Gibbs energy minimization and redox potential with Thermochimica. Incorporating a plate-out model, the full multi-physics model is used to evaluate the effects of deposition and removal of solid compounds formed in the molten salt in the reactor surfaces, evaluating the concentrations and regions in which nickel, chromium and iron would be deposited. Results indicate an almost isothermal state of the fuel during operation at 25 kW that do not increase the plate-out effects due to low temperature gradients.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Thermochemistry of sodium rare earth ternary fluorides, NaREF 4

Sodium rare earth fluorides, NaREF 4 (RE = rare earth), are used as luminescent materials for light emission and biomedical applications and are important compositions for extracting and separating RE elements. Solution calorimetric measurements of a series of β-structured NaREF 4 (Na 1.5 RE 1.5 F 6 ) phases with various RE elements determined their heats of formation. Though the lattice contracts from light to heavy RE elements, NaREF 4 compounds show more exothermic enthalpies of formation from binary components with the decrease of RE 3+ radii, contrary to behavior seen in most RE oxide ternary compounds. Furthermore, by constructing Born–Haber cycles, the different slopes of lines relating lattice energies to lattice parameters in binary and ternary fluorides appear to be the reason for this reverse trend, which can be associated with changes in the coordination number of RE cations. These trends and the metastability of sodium light RE fluorides not only reveal the key role of ionic radius in RE compound stability, but also are significant for the design and synthesis of new materials and motivate the more effective utilization of RE.

36 MATERIALS SCIENCE↗

Coupled neutronics, thermochemistry, corrosion modeling and sensitivity analyses for isotopic evolution in molten salt reactors

This study presents a computational methodology for analyzing isotopic evolution and associated uncertainties in molten salt reactors (MSRs), focusing on both fluoride- and chloride-based fuel salts. The primary goal is to enhance the understanding of isotopic behavior in MSRs and provide data to support future experimental efforts. The methodology integrates transport-coupled depletion calculations using OpenMC, equilibrium thermodynamics modeling with Thermochimica, and a corrosion model. Sensitivity analyses are performed to evaluate the impact of power density, air ingress, and humidity content on isotopic evolution in MSR concepts. This study examines representative F- and Cl-based MSR designs, highlighting the dominant influence of power density on isotopic composition, which significantly affects isotope production and depletion rates, accounting for approximately 76% of the observed variance in element concentration. Air ingress and humidity content also affect the redox potential, solubility of heavier elements, and corrosion rates, thereby altering the expected isotopic evolution in the reactor. On average, air ingress accounts for around 17% of the variance in element concentrations, while humidity explains the remaining 7%. These variances differ significantly from element to element, depending on the element’s role in depletion, redox potential evolution, and galvanic corrosion. The findings indicate that power density, air ingress, and humidity content are all critical factors for optimizing reactor design and operational strategies. Furthermore, the study provides expected ranges for key impurities in the fuel salt, which are crucial for guiding future experimental studies and refining MSR designs. Finally, this study demonstrates the importance of modeling depletion coupled with the evolution of redox potential and chemical interactions in MSR fuel salts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Redox Defect Thermochemistry of FeAl 2 O 4 Hercynite in Water Splitting from First-Principles Methods

Solar thermochemical hydrogen (STCH) production is a promising route to produce fuels from sunlight via high-temperature water splitting. However, efficient and technologically viable implementations of this process only allow a narrow window of thermodynamic boundary conditions that can be used to cycle the system, thus limiting the design space for suitable metal oxide redox mediators. An oxygen defect redox mechanism can contribute a favorable reduction entropy to expand this window, and computational evaluation of materials with high oxygen defect entropies could play a pivotal role in guiding the discovery and design of suitable oxides. This study employs first-principles calculations to investigate the redox mediating defect mechanism of the STCH candidate material, hercynite (FeAl 2 O 4 ). We compare the results of total energy calculations from density functional theory (DFT) with beyond-DFT approaches, including hybrid functionals and the random phase approximation, which are among the most advanced methods currently feasible for supercell defect calculations. Using the predicted formation energies, we perform thermodynamic modeling of FeAl2O4 reduction and oxidation via free energy minimization that incorporates ideal gas, configurational, and vibrational entropy contributions evaluated within the quasi-harmonic approximation. Special attention is devoted to understanding interactions among co-existing defects, such as the association of pairs and complexes of O vacancies and cation antisite defects, and the effect of mutually compensating defect charges. Our results corroborate the notion that the details of defect interactions can be decisive for the viability of hydrogen production within the desirable STCH process window.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermochemistry of Layered and Two-Dimensional Niobium Carbo-Chalcogenides

Two-dimensional transition metal carbo-chalcogenides (TMCCs) represent a novel class of layered materials with tunable electronic structures and high chemical versatility, making them promising candidates for energy and environmental applications. In this work, the structural and thermochemical properties of a series of TMCCs are investigated. High-temperature oxidative solution calorimetry, conducted at 800 °C in molten sodium molybdate (3Na2O·4MoO3), was used to determine the enthalpies of oxidation and formation from the elements at 25 °C, enabling a comprehensive thermodynamic assessment. Comparison of these enthalpies reveals distinct stability trends across the series. Cu0.67Nb2S2C exhibits the most exothermic oxidation enthalpy, indicating a strong thermodynamic driving force for oxidation and reduced resistance to oxidative degradation. Its positive enthalpy of formation from the elements suggests lower thermodynamic stabilization relative to the other compositions, consistent with the energetic influence of Cu and residual Fe within the Nb–S–C framework. In contrast, delaminated Nb2S2C exhibits poor oxidative stability but a highly exothermic enthalpy of formation, consistent with a defect-rich, metastable structural state rather than enhanced lattice stability. These findings establish correlations between crystal structure and energetic stability, providing insights into the design of robust TMCC-based materials for advanced energy and environmental technologies.

Cassell, Nakeshma [Clemson University]↗

Route to Chemical Accuracy for Computational Uranium Thermochemistry

Benchmark spinor-based relativistic coupled-cluster calculations for the ionization energies of the uranium atom, the uranium monoxide molecule (UO), and the uranium dioxide molecule (UO 2 ) and for the bond dissociation energies of UO and UO 2 are reported. The accuracy of these calculations in the treatments of relativistic, electron correlation, and basis-set effects is analyzed. Furthermore, the intrinsic convergence of the computed results and the favorable comparison with the experimental values demonstrate the unique applicability of the spinor representation of quantum-chemical methods to open-shell uranium-containing atomic and molecular species with uranium oxidation states ranging from U(0) to U(V).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reducing the Cost of CCSD Basis Set Extrapolation in Ab Initio Computational Thermochemistry

Here, a series of approximations to CCSD contributions in computational model chemistries is presented in the context of kcal mol –1 , kJ mol –1 , and 20 cm –1 theoretical predictions of total atomization energies, benchmarked within the HEAT+CH 4 test suite. A specific set of circumstances where MP2, without empirical scaling, may be used as an effective intermediate in the first two of these accuracy ranges was determined. However, SDQ-MP4, a method long used in pursuit of kcal mol –1 accuracy but relatively unstudied in the subchemical accuracy community, offers significant improvement over the quality of MP2 as a basis-set intermediate at significantly reduced cost compared to CCSD. Given this, we argue for SDQ-MP4 as the de facto CCSD basis-set intermediate in sub-chemical accuracy calculations when CCSD in a desired basis set becomes unaffordable. We additionally report on a “CBS-like” scheme, where MP2 and SDQ-MP4 are used in conjunction to create a “cheap” three-part approximation of large CCSD basis set limits. The data for the CCSD approximation schemes are organized in such a way that model chemistry developers can locate an analog of their current approach for the CCSD basis set limit and explore alternative intermediates that either decrease computational cost or increase computational accuracy. We also show, for a handful of molecules, that SDQ-MP4 shows promise as an effective basis-set intermediate for harmonic and fundamental frequency computations, allowing for zero-point corrections of nearly CCSD(T)/ANO1 quality using simple composite methods that only require CCSD(T)/ANO0.

Thorpe, James H. [Argonne National Laboratory (ANL↗

Molecular Properties of Thorium Hydrides: Electron Affinities and Thermochemistry

In this work, high-level electronic structure calculations of the ground and low-lying energy electronic states for ThH x and ThH x – for x = 2–5 are reported and compared to available anion photoelectron detachment experiments. The adiabatic electron affinities (EAs) are predicted to be 0.82, 0.88, 0.51, and 2.36 eV for x = 2 to 5, respectively, at the Feller– Peterson–Dixon (FPD) level. The vertical detachment energies (VDEs) are predicted to be 0.84, 0.88, 0.81, and 4.38 eV for x = 2–5, respectively. The corresponding experimental VDEs are 0.871 eV for x = 2, 0.88 eV for x = 3, and 4.09 eV for x = 5. As for ThH, there is a significant spin–orbit (SO) correction for the EA of ThH 2 , and this correction decreases substantially for x > 2. The observed ThH 2 – photoelectron spectrum has many transitions as predicted at the CASPT2-SO level. The FPD bond dissociation energies (BDEs) increase from 67 to 75 kcal/mol for x = 2 to x = 4 at the FPD level. The BDE for ThH 5 is much lower as it is a complex of H 2 with ThH 3 . The hydride affinities for x = 2 to 4 are all comparable and near 70 kcal/mol. A natural bond orbital analysis is consistent with a significant Th + –H– ionic contribution to the Th–H bonds. There is very little participation of the 5f orbitals in the bonding and the valence electrons on the Th are dominated by 7s and 6d for the neutrals and anions except for ThH 2 – where there is a significant contribution from the 7p.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theory-Based Mechanism for Fluoromethane Combustion I: Thermochemistry and Abstraction Reactions

A new detailed chemical kinetic mechanism is presented for small fluorinated hydrocarbons. Ab initio electronic structure theory is used to provide heats of formation with subchemical accuracy. The ANL0 method is extended to include fluorine. The resulting heats of formation at 0 K are in excellent agreement with 36 benchmark species in the Active Thermochemical Tables, with a mean error of μ = –0.02 kJ/mol and a standard deviation of σ = 0.91 kJ/mol. The thermophysical properties for 92 small-molecule H/C/O/F species are computed. The rate coefficients for 40+ H-abstraction reactions involving H, O, F, OH, OF, HO 2 , and various methyl radicals with CH 4 , CH 3 F, CH 2 F 2 , CHF 3 , CH 2 O, and CHFO are discussed. Here, the computed rate constants are in excellent agreement with the available literature. Additionally, 30+ rate constants are provided for F abstraction, which are several orders of magnitude smaller than H abstraction. The thermophysical properties and rate constants are provided in a mechanism. This mechanism is the first in a series of theory-based investigations into the thermal destruction of per- and polyfluorinated species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗