Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermodynamics of computing”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

First-principles modeling of chemistry in mixed solvents: Where to go from here?

Mixed solvents (i.e., binary or higher order mixtures of ionic or nonionic liquids) play crucial roles in chemical syntheses, separations, and electrochemical devices because they can be tuned for specific reactions and applications. Apart from fully explicit solvation treatments that can be difficult to parameterize or computationally expensive, there is currently no well-established first-principles regimen for reliably modeling atomic-scale chemistry in mixed solvent environments. We offer our perspective on how this process could be achieved in the near future as mixed solvent systems become more explored using theoretical and computational chemistry. We first outline what makes mixed solvent systems far more complex compared to single-component solvents. An overview of current and promising techniques for modeling mixed solvent environments is provided. We focus on so-called hybrid solvation treatments such as the conductor-like screening model for real solvents and the reference interaction site model, which are far less computationally demanding than explicit simulations. We also propose that cluster-continuum approaches rooted in physically rigorous quasi-chemical theory provide a robust, yet practical, route for studying chemical processes in mixed solvents.

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

Nonintrusive manufactured solutions for ablation

The study of heat transfer and ablation plays an important role in many problems of scientific and engineering interest. As with the computational simulation of any physical phenomenon, the first step toward establishing credibility in ablation simulations involves code verification. Code verification is typically performed using exact and manufactured solutions. However, manufactured solutions generally require the invasive introduction of an artificial forcing term within the source code such that the code solves a modified problem for which the solution is known. In this paper, we present a nonintrusive method for manufacturing solutions for a non-decomposing ablation code, which does not require the addition of a source term.

42 ENGINEERING↗
DOE OSTI2021

Predicting large-scale pool fire dynamics using an unsteady flamelet- and large-eddy simulation-based model suite

A low-Mach, unstructured, large-eddy-simulation-based, unsteady flamelet approach with a generalized heat loss combustion methodology (including soot generation and consumption mechanisms) is deployed to support a large-scale, quiescent, 5-m JP-8 pool fire validation study. The quiescent pool fire validation study deploys solution sensitivity procedures, i.e., the effect of mesh and time step refinement on capturing key fire dynamics such as fingering and puffing, as mesh resolutions approach O(1) cm. A novel design-order, discrete-ordinate-method discretization methodology is established by use of an analytical thermal/participating media radiation solution on both low-order hexahedral and tetrahedral mesh topologies in addition to quadratic hexahedral elements. Here, the coupling between heat losses and the flamelet thermochemical state is achieved by augmenting the unsteady flamelet equation set with a heat loss source term. Soot and radiation source terms are determined using flamelet approaches for the full range of heat losses experienced in fire applications including radiative extinction. The proposed modeling and simulation paradigm are validated using pool surface radiative heat flux, maximum centerline temperature location, and puffing frequency data, all of which are predicted within 10% accuracy. Simulations demonstrate that under-resolved meshes predict an overly conservative radiative heat flux magnitude with improved comparisons as compared to a previously deployed hybrid Reynolds-averaged Navier–Stokes/eddy dissipation concept-based methodology.

42 ENGINEERING↗
DOE OSTI2022

Quantum Computing by Coherent Cooling

Interesting problems in quantum computation take the form of finding low-energy states of (pseudo)spin systems with engineered Hamiltonians that encode the problem data. Motivated by the practical possibility of producing very low-temperature spin systems, we propose and exemplify the possibility to compute by coupling the computational spins to a quantum coherent bath that serves as a heat sink. The quantum tunneling effect provides additional cooling channels to accelerate the cooling process. Here, we demonstrate both analytically and numerically that this strategy can achieve a quantum advantage in the unstructured search problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Structure and Dynamics of Domains in Ferroelectric Nanostructures – Phase-Field Modeling (Final Technical Report)

Domain pattern formation is one of the most common phenomena in nature and is a topic of immense interest in many fields ranging from materials science and physics to chemistry and biology. This DOE sponsored research project explored the basic science concerning the thermodynamic stability of mesoscale polarization domain patterns and their temporal evolution mechanisms during formation and subsequent switching in ferroelectric nanostructures and heterostructures. The project employed the computational phase-field method in combination of microelasticity and electrostatic theories. The research was carried out in close collaboration with a number of experimental groups who used High Resolution Transmission Electron Microscopy (HRTEM), in situ TEM with Scanning Probe Microscopy (SPM), or Piezoresponse Force Microscopy (PFM) to characterize the domain structures and dynamics of ferroelectric thin films and heterostructures and who grow high-quality ferroelectric and multiferroic thin films using advanced growth techniques such as Molecular Beam Epitaxy (MBE), Pulsed Laser Deposition (PLLD), and sputtering. The research efforts of the project help establish the phase-field method as the most powerful method for understanding and predicting domain structures in ferroelectric thin films and nanostructures. The findings of the project led to the basic understanding of stability and switching mechanisms of ferroelectric domains under different mechanical boundary conditions and under either homogeneous capacitor configurations or local fields using metallic probes as electrodes and the emergence of charged domain walls during domain switching. The project predicted the spatial length scales, temperature ranges, and electromechanical conditions for different polar states in heterostructures and guided the discovery of both transient and stable novel polarization states containing vortex lattices in oxide superlattices. The project resulted in 134 journal publications and 6 PhD theses with all the PhD graduates currently working in either academia or industry within the United States. The basic understanding on the stability of mesoscale polar states and pattern evolution achieved by the project improved our ability to control and engineer properties of ferroelectric thin films and heterostructures for potential applications in nanoscale electronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Development and Validation of Physics-Based Sub-Models of High-Pressure Supercritical Fuel Injection at Diesel Conditions Insert

The project sought to develop and validate accurate, physics-based, numerical submodels and implement them in computational fluid dynamics (CFD) software codes to enable accurate predictions of diesel sprays at pressures higher than the current systems. The combined experimental-computational effort has four main objectives: (1) acquire spatially and temporally resolved scalar measurements of fuel-oxidizer mixing in the near field of the fuel jet, (2) develop software modules to calculate thermo-physical properties of diesel surrogates and their mixtures with oxidizer of interest in real-gas operating regime, (3) integrate the real-fluid model into a commercial and an open-source CFD code to simulate fuel-air mixing in high-pressure diesel sprays, and (4) assess the robustness, accuracy, and uncertainty of the integrated CFD solvers in computing diesel sprays. Project accomplishments are consistent with the original goals of the project. Although we acquired spatially and temporally resolved scalar measurements, they were limited to refractive index difference and not to thermodynamic properties such as temperature and equivalence ratio. However, we successfully developed the theoretical background needed to acquire such measurements at diesel conditions and demonstrated the capability of the RSD for scalar measurements in a canonical turbulent jet flow replicating features of the diesel jet. Additional work will be needed to implement these procedures at high-pressure, high-temperature engine-like conditions of the CPFR. Mainly, we did not anticipate the need to find the liquid boundary and its statistical features or the state relationships between refractive index and thermodynamic properties, assuming the widely known Dale-Gladstone equation would suffice. Still, efforts devoted to these tasks were essential to fulfil the goals of the project and serves as the foundation on which many of the conclusions and key contributions of this project are based. Further assessment and application of the developed CFD codes is recommended to fully realize the benefits of this research.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Thermodynamic Modeling of Complex Solid Solutions in the Lu-H-N System via Graph Neural Network Accelerated Monte Carlo Simulations

Metal hydrides are important across diverse applications, such as hydrogen storage, batteries, gas sensors, nuclear reactions, and high-temperature superconductivity. Previous computational studies of metal hydrides under extreme pressures, e.g., 𝑂⁡(10 2 ) ⁢GPa, usually treat them as stoichiometric compounds without considering interstitial lattice disorder. As pressures become more moderate in the 𝑂⁡(10 0 ) ⁢GPa and below range, hydrogen disorder at interstitial lattice sites becomes prominent, e.g., in the N-doped Lu hydride that was recently claimed superconducting near 1 GPa. Further adding compositional complexity from alloying and/or multielement interstitial occupation makes elucidating pressure- and temperature-dependent observables intractable by first-principles calculations alone. We therefore propose a lattice graph neural-network surrogate modeling approach to predict configuration- and pressure-dependent equation-of-state properties. Their efficiency permits Monte Carlo simulations to calculate Gibbs energies and pressure-dependent phase diagrams, thereby revealing insights into the synthesis conditions required for achieving desired phase equilibria. We demonstrate this concept for the compositionally complex cubic Lu(H,N,Va) 3 system where three constituents (hydrogen, nitrogen and vacancy) have disordered multielement interstitial occupancies and insights into pressure-dependent phase equilibria are critically needed, e.g., N-doping levels can significantly lower dehydrogenation temperatures and provide a new strategy to optimize hydrogen-storage alloys. This work can improve the thermodynamic understanding of the Lu-H-N system and help rational synthesis of N-doped Lu hydrides, but more generally demonstrates an efficient approach to model pressure-dependent thermodynamics of multicomponent solid solutions.

Monte Carlo methods↗
DOE OSTI2024

2D-EFICACY: Control of Metastable 2D Carbide[1]Chalcogenide Heterolayers: Strain and Moire Engineering

The experimental isolation of graphene led to the discovery of an entirely new world of two-dimensional (2D) materials in which the 2D nature often leads to emergent behaviors not seen in bulk systems. 2D transition metal dichalcogenides (TMDs) exhibit physico-chemical properties that depend on the transition metal, polymorph, thickness, and presence and type of defects. Recently, a group of thin (10-100nm) transition metal carbides (TMCs), such as Mo2C, has been synthesized that exhibit a thickness-dependent superconducting critical temperature (Tc). These thin TMCs are different from MXenes, another class of 2D materials consisting of few layers of nitrides or carbides (<5nm) produced by chemical etching and delamination. The goal of this renewal proposal is to combine experiment and computation to synthesize and elucidate the guiding principles that control the growth, orientation and strain of heterostacks of thin TMCs and TMDs composed with Nb, Ti and W. We expect to stabilize metastable hybrid phases of TMCs sandwiched between TMDs (H-TMD/Cs) with unprecedented physico-chemical properties. As part of previous DOE-funded work by the Terrones/Sinnott groups, thin (10-100 nm thick) Mo2C flakes were successfully synthesized by chemical vapor deposition (CVD). By subsequently exposing Mo2C to H2S, partial chalcogenization was achieved, resulting in heterostacks of MoCx phases and MoS2. The formation of MoS2 led to a deficiency of Mo atoms in the underlying Mo2C, resulting in an inhomogeneous phase change from α-Mo2C to γ’-MoCx and then to γ-MoC. The γ’-MoCx is a strained metastable phase and the heterostack of all three phases demonstrated an increased Tc relative to that of α-Mo2C, from 4 to 6K; its interleaved layered structure consisting of superconducting and semiconducting phases is ideal for future studies of Josephson junction series arrays. Moiré patterns in these heterostacked systems could result in new phenomena, as moiré patterns in bilayer graphene showed unconventional superconductivity and moiré excitons have been observed in twisted TMD heterobilayers. The scientific hypothesis of the proposed synergistic computational and experimental research is that orientation and strain control within confined thin metastable TMCs, sandwiched by stable phases of TMCs and layered TMDs, will depend on kinetic and thermodynamic “knobs” that include fast temperature changes, chalcogen diffusion through preferred crystallographic planes, reaction times, pressure, reactive atmosphere, precursors, and surfactants, which will also tailor properties such as superconductivity, magnetism, ferroelectricity, piezoelectricity, and catalytic performance. We will develop the guiding principles for the synthesis and stabilization of metastable H-TMD/Cs based on Nb, Ti and W. In order to validate the hypothesis, four tasks are proposed: The first task will synthesize ultra-thin TMCs based on Nb, W and Ti, by: 1) adapting the CVD method used for Mo2C, 2) plasma assisted CVD, 3) defect-mediated CVD processes, and 4) cryo-milling of carbide powders. The second task will accomplish the synthesis and basic physico-chemical characterizations of H-TMD/Cs by chalcogenization of the materials synthesized in task one, and by carbonization of TMDs. H-TMD/Cs will also be investigated for their suitability in energy conversion applications such as supercapacitors, Li and multivalent ion batteries, and electrocatalysts, topics of interest to DoE. These tasks will be carried out in close conjunction with density functional theory (DFT) calculations with insights into energetics, lattice parameters, stability, phase diagrams, band structures, and density of states of H-TMD/Cs. The third task will characterize and evaluate strain and moiré patterns at the interfaces of different H-TMD/Cs by high-resolution scanning transmission electron microscopy (HR-STEM), scanning tunneling microscopy (STM), and conductive tip atomic force microscopy. Nudged elastic band calculations with DFT will be performed to understand the chalcogen diffusion process, which will provide insights into the interfaces between different phases of TMCs and TMDs. The fourth task aims at quantifying the stability and dynamics of H-TMD/Cs by in-situ TEM and Raman studies under heating, strain, and electrical biasing. Phonon calculations using DFT will provide a basis for interpreting Raman spectra. This coherent framework involving synthesis, characterization, and computation will result in a broad scientific impact for energy related applications. The ability to develop new H-TMD/Cs will enhance a range of applications that include batteries, catalysts, switches, sensors, quantum computing components and smart coatings.

2-Dimensional materials↗
DOE OSTI2024

2D-EFICACY: Control of Metastable 2D Carbide Chalcogenide Heterolayers: Strain and Moire Engineering

The experimental isolation of graphene led to the discovery of an entirely new world of two-dimensional (2D) materials in which the 2D nature often leads to emergent behaviors not seen in bulk systems. 2D transition metal dichalcogenides (TMDs) exhibit physico-chemical properties that depend on the transition metal, polymorph, thickness, and presence and type of defects. Recently, a group of thin (10-100nm) transition metal carbides (TMCs), such as Mo 2 C, has been synthesized that exhibit a thickness-dependent superconducting critical temperature (Tc). These thin TMCs are different from MXenes, another class of 2D materials consisting of few layers of nitrides or carbides (<5nm) produced by chemical etching and delamination. The goal of this renewal proposal is to combine experiment and computation to synthesize and elucidate the guiding principles that control the growth, orientation and strain of heterostacks of thin TMCs and TMDs composed with Nb, Ti and W. We expect to stabilize metastable hybrid phases of TMCs sandwiched between TMDs (H-TMD/Cs) with unprecedented physico-chemical properties. As part of previous DOE-funded work by the Terrones/Sinnott groups, thin (10-100 nm thick) Mo 2 C flakes were successfully synthesized by chemical vapor deposition (CVD). By subsequently exposing Mo 2 C to H 2 S, partial chalcogenization was achieved, resulting in heterostacks of MoCx phases and MoS 2 . The formation of MoS 2 led to a deficiency of Mo atoms in the underlying Mo 2 C, resulting in an inhomogeneous phase change from α-Mo 2 C to γ’-MoCx and then to γ-MoC. The γ’-MoCx is a strained metastable phase and the heterostack of all three phases demonstrated an increased Tc relative to that of α-Mo 2 C, from 4 to 6K; its interleaved layered structure consisting of superconducting and semiconducting phases is ideal for future studies of Josephson junction series arrays. Moiré patterns in these heterostacked systems could result in new phenomena, as moiré patterns in bilayer graphene showed unconventional superconductivity and moiré excitons have been observed in twisted TMD heterobilayers. The scientific hypothesis of the proposed synergistic computational and experimental research is that orientation and strain control within confined thin metastable TMCs, sandwiched by stable phases of TMCs and layered TMDs, will depend on kinetic and thermodynamic “knobs” that include fast temperature changes, chalcogen diffusion through preferred crystallographic planes, reaction times, pressure, reactive atmosphere, precursors, and surfactants, which will also tailor properties such as superconductivity, magnetism, ferroelectricity, piezoelectricity, and catalytic performance. We will develop the guiding principles for the synthesis and stabilization of metastable H-TMD/Cs based on Nb, Ti and W. In order to validate the hypothesis, four tasks are proposed: The first task will synthesize ultra-thin TMCs based on Nb, W and Ti, by: 1) adapting the CVD method used for Mo 2 C, 2) plasma assisted CVD, 3) defect-mediated CVD processes, and 4) cryo-milling of carbide powders. The second task will accomplish the synthesis and basic physico-chemical characterizations of H-TMD/Cs by chalcogenization of the materials synthesized in task one, and by carbonization of TMDs. H-TMD/Cs will also be investigated for their suitability in energy conversion applications such as supercapacitors, Li and multivalent ion batteries, and electrocatalysts, topics of interest to DoE. These tasks will be carried out in close conjunction with density functional theory (DFT) calculations with insights into energetics, lattice parameters, stability, phase diagrams, band structures, and density of states of H-TMD/Cs. The third task will characterize and evaluate strain and moiré patterns at the interfaces of different H-TMD/Cs by high-resolution scanning transmission electron microscopy (HR-STEM), scanning tunneling microscopy (STM), and conductive tip atomic force microscopy. Nudged elastic band calculations with DFT will be performed to understand the chalcogen diffusion process, which will provide insights into the interfaces between different phases of TMCs and TMDs. The fourth task aims at quantifying the stability and dynamics of H-TMD/Cs by in-situ TEM and Raman studies under heating, strain, and electrical biasing. Phonon calculations using DFT will provide a basis for interpreting Raman spectra. This coherent framework involving synthesis, characterization, and computation will result in a broad scientific impact for energy related applications. The ability to develop new H-TMD/Cs will enhance a range of applications that include batteries, catalysts, switches, sensors, quantum computing components and smart coatings.

2-Dimensional Materials↗
DOE OSTI2021

Solubility Product of Vanadinite Pb 5 (VO 4 ) 3 Cl at 25 °C—A Comprehensive Approach to Incongruent Dissolution Modeling

Although vanadinite (Pb 5 (VO 4 ) 3 Cl occurs in abundance in various terrestrial geochemical systems of natural and anthropogenic origin and is seriously considered as a potential nuclear waste sequestering agent, its actual application is severely limited by a lack of understanding of its basic thermodynamic parameters. In this regard, the greatest challenge is posed by its incongruent dissolution, which is a pivotal hurdle for effective geochemical modeling. Our paper presents an universal approach for geochemical computing of systems undergoing incongruent dissolution which, along with unique, long-term experiments on vanadinites’ stability, allowed us to determine the mineral solubility constant. The dissolution experiments were carried out at pH = 3.5 for 12 years. Vanadinite has dissolved incongruently, continuously re-precipitating into chervetite (Pb 2 V 2 O 7 ) with the two minerals remaining in mutual equilibrium until termination of the experiments. The empirically derived solubility constant K sp,V,298 = 10 –91.89 ± 0.05 of vanadinite was determined for the first time. The proposed modeling method is versatile and can be adopted to other mineral systems undergoing incongruent dissolution.

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

Force Field X: A computational microscope to study genetic variation and organic crystals using theory and experiment

Force Field X (FFX) is an open-source software package for atomic resolution modeling of genetic variants and organic crystals that leverages advanced potential energy functions and experimental data. FFX currently consists of nine modular packages with novel algorithms that include global optimization via a many-body expansion, acid–base chemistry using polarizable constant-pH molecular dynamics, estimation of free energy differences, generalized Kirkwood implicit solvent models, and many more. Applications of FFX focus on the use and development of a crystal structure prediction pipeline, biomolecular structure refinement against experimental datasets, and estimation of the thermodynamic effects of genetic variants on both proteins and nucleic acids. The use of Parallel Java and OpenMM combines to offer shared memory, message passing, and graphics processing unit parallelization for high performance simulations. Overall, the FFX platform serves as a computational microscope to study systems ranging from organic crystals to solvated biomolecular systems.

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

Combining Theory and Experiment to Map the Atomic-Level Structure–Energy Pathways of Adsorbate-Mediated Phase Changes in a Cooperatively Flexible Metal–Organic Framework

An important subclass of metal–organic frameworks (MOFs) exhibits cooperative flexibility, wherein individual crystallites undergo global structural phase changes in response to external stimuli. Where cooperative flexibility results in reversible changes between crystalline states of distinct accessible porosity, these frameworks can exhibit rare yet desirable behaviors that cannot be explained by local dynamics alone. Yet, the chemical and structural origins of cooperative flexibility and how frameworks undergo these reversible phase changes at the atomic level remain poorly understood. Deliberate design for specific applications is therefore exceedingly difficult, and there is great impetus to develop a fundamental understanding of this phenomenon. Here, an effective and widely accessible computational approach is developed, which is designed to provide microscopic resolution via direct comparison to experimental data along the desorption-guided pathway. The strategy is applied to explain the desorption-induced phase change in an experimentally well-characterized framework, CdIF-13 (sod-Cd(benzimidazolate)2), where experiment alone was unable to resolve the atomistically detailed phase change landscape. Our findings reveal that the cooperative phase change pathways are adsorbate dependent with thermodynamics of intermediate structural states dictated by a nuanced interplay of ligand orientation, skeletal symmetry, and modes of surface adsorption. The results reveal that this isotropically flexible framework is “chaperoned” through a complex energy landscape by specific adsorbates, revealed by the reported computational approach with atomic-level insight and validated by experimentally determined structures. Thus, this work facilitates both understanding and future design of flexible materials for applications in gas storage, transport, delivery, and separation technologies.

03 NATURAL GAS↗
DOE OSTI2022

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↗
DOE OSTI2020

Thermomechanical conversion in metals: dislocation plasticity model evaluation of the Taylor-Quinney coefficient

Using a partitioned-energy thermodynamic framework which assigns energy to that of atomic configurational stored energy of cold work and kinetic-vibrational, in this study we derive an important constraint on the Taylor-Quinney coefficient, which quantifies the fraction of plastic work that is converted into heat during plastic deformation. Associated with the two energy contributions are two separate temperatures – the ordinary temperature for the thermal energy and the effective temperature for the configurational energy. We show that the Taylor-Quinney coefficient is a function of the thermodynamically defined effective temperature that measures the atomic configurational disorder in the material. Finite-element analysis of recently published experiments on the aluminum alloy 6016-T4 [1], using the thermodynamic dislocation theory (TDT), shows good agreement between theory and experiment for both stress-strain behavior and temporal evolution of the temperature. The simulations include both conductive and convective thermal energy loss during the experiments, and significant thermal gradients exist within the simulation results. Computed values of the differential Taylor-Quinney coefficient are also presented and suggest a value which differs between materials and increases with increasing strain.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Non-intrusive data-driven model reduction for differential–algebraic equations derived from lifting transformations

In this paper we present a non-intrusive data-driven approach for model reduction of nonlinear systems. The approach considers the particular case of nonlinear partial differential equations (PDEs) that form systems of partial differential–algebraic equations (PDAEs) when lifted to polynomial form. Such systems arise, for example, when the governing equations include Arrhenius reaction terms (e.g., in reacting flow models) and thermodynamic terms (e.g., the Helmholtz free energy terms in a phase-field solidification model). Using the known structured form of the lifted algebraic equations, the approach computes the reduced operators for the algebraic equations explicitly, using straightforward linear algebra operations on the basis matrices. The reduced operators for the differential equations are inferred from lifted snapshot data using operator inference, which solves a linear least squares regression problem. The approach is illustrated for the nonlinear model of solidification of a pure material. The lifting transformations reformulate the solidification PDEs as a system of PDAEs that have cubic structure. The operators of the lifted system for this solidification example have affine dependence on key process parameters, permitting us to learn a parametric reduced model with operator inference. Numerical experiments show the effectiveness of the resulting reduced models in capturing key aspects of the solidification dynamics.

42 ENGINEERING↗
DOE OSTI2021

A Machine Learning Approach for the Prediction of Formability and Thermodynamic Stability of Single and Double Perovskite Oxides

Perovskite oxides continue to attract huge interest due to their fascinating and wide-ranging properties for diverse applications. The tunability of these properties may be further enhanced by increasing their compositional complexity via double perovskite-ordered configurations containing multiple cations. In this work, we focus on an exhaustive chemical space of single and double oxide perovskites and optimally explore this space to identify novel compositions that are likely to form stable compounds. Critically, we examine the relationship between formability, the practical ability to synthesize a compound, and stability, the thermodynamic preference to form the structure. Our formability and stability training data sets were enumerated from the available experimental literature and in-house density functional theory computations and contained 1505 and 3469 examples, respectively, representing state-of-the-art in the current open literature in perovskite and double perovskite compounds. Subsequently, cross-validated and highly accurate machine learning classification models are built using these training data sets and employed to screen for novel stable oxide perovskites. The study identifies (1) atomic features relevant to prediction of formability and stability in perovskite and double perovskite compounds, (2) the importance of including energy contributions due to local structural relaxations going beyond the high symmetry perovskite phase, and (3) 437,828 double perovskite compounds that are likely to be stable and 891,188 compounds that are likely to be formable. From the intersection of this large chemical space of formable and stable oxide perovskites, 414 compositions are identified as the most promising candidates for future experimental synthesis of novel oxide perovskites. The developed models may be generalized and have implications beyond perovskite discovery if applied to other families of compounds.

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

Assessing the accuracy of compound formation energies with quantum Monte Carlo

Accurately predicting the formation energy of a compound, which describes its thermodynamic stability, is a key challenge in materials physics. Here, we employ a many-body quantum Monte Carlo (QMC) with single-reference trial functions to compute the formation energy of two electronically disparate compounds, the intermetallic VPt2 and the semiconductor CuI, for which standard density functional theory (DFT) predictions using both the Perdew-Burke-Ernzerhof (PBE) and the strongly constrained and appropriately normed (SCAN) density functional approximations deviate markedly from available experimental values. For VPt2, we find an agreement between QMC, SCAN, and PBE0 estimates, which therefore remain in disagreement with the much less exothermic experimental value. For CuI, the QMC result agrees with neither SCAN nor PBE pointing towards DFT exchange-correlation biases, likely related to the localized Cu 3d electrons. Compared to the behavior of some density functional approximations within DFT, spin-averaged QMC exhibits a smaller but still appreciable deviation when compared to experiment. The QMC result is slightly improved by incorporating spin-orbit corrections for CuI and solid I 2 , so that experiment and theory are brought into imperfect but reasonable agreement within about 120 meV/atom.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Investigation of Frozen Chemistry for Molten Salt Reactors (MSRs)

Understanding accident progression and the potential/conditions for fission product release from fuel is necessary to evaluate safety for any nuclear reactor system. Molten Salt Reactors (MSRs) under development need such analysis to support safety evaluations. Fission product chemistry specific to MSR concepts is a critical area that introduces distinct considerations relative to the current state-of-knowledge in reactor safety, primarily developed for water-moderated nuclear reactor systems. In Light Water Reactor (LWR) systems, it is necessary to capture the chemical interaction of fission products with the reactor evironment, containment and confinement systems. The overall effects at this point are relatively well understood for the purposes of performing safety evaluations. A key insight from LWR studies is that fission product chemical behavior can be reasonably captured by modeling approaches where the chemistry is "frozen". These modeling approaches assume that radionuclide reaction and speciation can be represented by chemical classes, each with characteristic transport behavior that is invariant under a broad range of thermochemical conditions. However, radionuclides can exhibit a range of behavior in the liquid salt-melt phase of the coolant used in MSRs. Radionuclides, salt, and the metal containment surfaces (i.e. pipes) can co-exist in dynamic equilibrium that could evolve with small system mass changes. A detailed investigation to the degree the equilibrium state can dynamically evolve with changes in the conditions of the molten salt mixture has not been previously conducted. It is currently not well understood where frozen chemistry assumptions are valid. Expanding the state-of-knowledge in this regard is relevant to better assessing the range of chemical effects that should be incorporated as part of MSR safety assessments. This investigation used the Oak Ridge Isotope GENeration (ORIGEN) module of the Standardized Computer-Analysis for Licensing Evaluation (SCALE) code to generate simulated radionuclide inventories for the MSR Experiment (MSRE) and then modeled reactor chemical speciation using the Molten Salt Thermodynamic Database – Thermochemical (MSTDB-TC) coupled with Thermochimica. The effect of composition variation during decay of fission product inventory in a molten salt over a period of 500 days prolonged post- at multiple temperatures was studied. Mass fractions for fluorine and berilium were varied in order to probe the effects of free fluorine control. Finally, speciation of fluoride reactors were showed by comparing MSRE readionuclide inventories with a FLiBe based molten salt breeder reactor (MSBR). The results showed that fission product mass change has little effect on phase mass changes and vapor pressures for fluoride species, but differ with varying carrier and fuel salt compositions. However, iodine species were found to have a vapor pressure not only dependent on temperature, but also the free fluorine potential, releasing iodine when the free fluorine potential is equal to the iodine inventory. This observation, however, arose under free fluorine potentials that are very unlikely to be realized in typical molten salt mixtures. Despite this observation, temperature was found to be the dominant parameter that drove phase change and fission product species vapor pressure. The results indicate that the current frozen chemistry approach is adequate for MSR analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗