Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coordinate 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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

SAS4A/SASSYS-1 Validation with the FFTF LOFWOS Tests #10-12

The FFTF liquid sodium cooled test reactor was built to assist the development and testing of advanced fuels and material for fast reactors. The goals of the FFTF Passive Safety Testing program included confirming the safety margins of FFTF as a liquid metal reactor, providing data for computer code validation, and demonstrating the inherent and passive safety benefits of its specific design features. The program included thirteen unprotected loss of flow without scram tests. The results of modeling and simulation of LOFWOS Tests #10, #11, and #12 are presented in this report to support validation of Argonne’s SAS4A/SASSYS-1 fast reactor safety analysis code. These three tests did not use the primary loop pony motors and thus transitioned to natural circulation flow rates from varying initial power levels. This validation activity leverages modeling and simulation efforts that originated under an International Atomic Energy Agency Coordinated Research Project (CRP) for benchmark analysis of the FFTF LOFWOS Test #13. The SAS4A/SASSYS-1 model originally developed for the CRP has been extended to include LOFWOS Tests #10-12 and modernized to utilize new code features and modeling practices. Predicted results from SAS4A/SASSYS-1 were compared to measured test data and differences were investigated. Overall, it was concluded that, for LOFWOS Tests #10-12, there was good agreement between the flow, power, and temperature predictions and the measured data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Dangling Bonds as Possible Contributors to Charge Noise in Silicon and Silicon–Germanium Quantum Dot Qubits

Spin qubits based on Si and Si 1–x Ge x quantum dot architectures exhibit among the best coherence times of competing quantum computing technologies, yet they still suffer from charge noise that limit their qubit gate fidelities. Identifying the origins of these charge fluctuations is therefore a critical step toward improving Si quantum-dot-based qubits. Here, we use hybrid functional calculations to investigate possible atomistic sources of charge noise, focusing on charge trapping at Si and Ge dangling bonds (DBs). We evaluate the role of global and local environment in the defect levels associated with DBs in Si, Ge, and Si 1–x Ge x alloys, and consider their trapping and excitation energies within the framework of configuration coordinate diagrams. We additionally consider the influence of strain and oxidation in charge-trapping energetics by analyzing Si and GeSi DBs in SiO 2 and strained Si layers in typical Si 1–x Ge x quantum dot heterostructures. Our results identify that Ge dangling bonds are more problematic charge-trapping centers both in typical Si 1–x Ge x alloys and associated oxidation layers, and they may be exacerbated by compositional inhomogeneities. In conclusion, these results suggest the importance of alloy homogeneity and possible passivation schemes for DBs in Si-based quantum dot qubits and are of general relevance to mitigating possible trap levels in other Si, Ge, and Si 1–x Ge x -based metal-oxide-semiconductor stacks and related devices.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Predicting metal–metal interactions. II. Accelerating generalized schemes through physical insights

Operando-computational frameworks that integrate descriptors for catalyst stability within catalyst screening paradigms enable predictions of rates and selectivity on chemically faithful representations of nanoparticles under reaction conditions. These catalyst stability descriptors can be efficiently predicted by density functional theory (DFT)-based models. The alloy stability model, for example, predicts the stability of metal atoms in nanoparticles with site-by-site resolution. In this work, we use physical insights to present accelerated approaches of parameterizing this recently introduced alloy-stability model. These accelerated approaches meld quadratic functions for the energy of metal atoms in terms of the coordination number with linear correlations between model parameters and the cohesive energies of bulk metals. By interpolating across both the coordination number and chemical space, these accelerated approaches shrink the training set size for 12 fcc p- and d-block metals from 204 to as few as 24 DFT calculated total energies without sacrificing the accuracy of our model. We validate the accelerated approaches by predicting adsorption energies of metal atoms on extended surfaces and 147 atom cuboctahedral nanoparticles with mean absolute errors of 0.10 eV and 0.24 eV, respectively. This efficiency boost will enable a rapid and exhaustive exploration of the vast material space of transition metal alloys for catalytic applications.

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

Axisymmetric hydrodynamics in numerical relativity using a multipatch method

In this paper, we describe a method of implementing the axisymmetric evolution of general-relativistic hydrodynamics and magnetohydrodynamics through modification of a multipatch grid scheme. In order to ease the computational requirements required to evolve the post-merger phase of systems involving binary compact massive objects in numerical relativity, it is often beneficial to take advantage of these system's tendency to rapidly settle into states that are nearly axisymmetric, allowing for 2D evolution of secular timescales. We implement this scheme in the spectral Einstein code and show the results of application of this method to four test systems including viscosity, magnetic fields, and neutrino radiation transport. Our results show that this method can be used to quickly allow already existing 3D infrastructure that makes use of local coordinate system transformations to be made to run in axisymmetric 2D with the flexible grid creation capabilities of multipatch methods. Our code tests include a simple model of a binary neutron star postmerger remnant, for which we confirm the formation of a massive torus which is a promising source of post-merger ejecta.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A redox model for NO oxidation, NH 3 oxidation and high temperature standard SCR over Cu-SSZ-13

A kinetic model is developed to predict the influence of temperature and hydrothermal aging on the redox of active Cu sites under standard SCR, NO oxidation and NH 3 oxidation conditions over a practically relevant fully-formulated Cu-SSZ-13 catalyst. NO 2 /N 2 production during NO/NH 3 titration of Cu II sites is utilized to identify rate parameters associated with NO-only RHC (reduction half cycle) and NH 3 -only RHC respectively. Integral N 2 formation during subsequent NO + NH 3 titration is consistent with the production of one NO 2 per two CuII sites reduced during NO-only RHC and one N 2 per six Cu II sites reduced during NH 3 -only RHC. Decreased reduction of Cu II sites by NO/NH 3 upon hydrothermal aging, along with the production of one NO 2 per two CuII sites during NO-only RHC, is accordant with the involvement of proximal ZCuOH and oxygen-bridged dimeric Cu II sites. Oxidation of partially solvated and framework coordinated ZCu (Cu I ) sites occurs in presence of O 2 , does not produce N 2 and can lead to the consumption of Brønsted acid sites. A global OHC kinetic model is developed to predict transient and integral N 2 formation during exposure of CuI sites to a mixture of NO and O 2 . The resulting redox kinetic model quantitatively predicts NO and NH 3 consumption during isothermal transient response Cu redox (TRCR) protocols, along with temperature and age dependent steady-state standard SCR and oxidation conditions. The redox model presented in this work synthesizes recent kinetic, spectroscopic and computational findings to provide a foundational description of active site redox during standard SCR, NO oxidation and NH3 oxidation over Cu-SSZ-13.

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

Leveraging Intramolecular π-Stacking to Access an Exceptionally Long-Lived 3 MC Excited State in an Fe(II) Carbene Complex

The ability to manipulate excited-state decay cascades using molecular structure is essential to the application of abundant-metal photosensitizers and chromophores. Ligand design has yielded some spectacular results elongating charge-transfer excited state lifetimes of Fe(II) coordination complexes, but triplet metal-centered ( 3 MC) excited states - recently demonstrated to be critical to the photoactivity of isoelectronic Co(III) polypyridyls - have to date remained elusive, with temporally isolable examples limited to the picosecond regime. Here, with this report, we show how strong-field donors and intramolecular π-stacking can conspire to stabilize a long-lived 3 MC excited state for a remarkable 4.1 ± 0.3 ns in fluid solution at ambient temperature. Analysis of variable-temperature time-resolved absorption data using theoretical models ranging from Arrhenius to semiclassical Marcus theory, combined with computational modeling and X-ray crystallography, reveal a Jahn−Teller stabilized excited state with a high activation barrier for ground-state recovery. The net result is a chromophore with a 3 MC excited-state lifetime that is orders of magnitude longer than anything yet observed for an Fe(II) complex.

carbene compounds↗
DOE OSTI2025

High-Pressure Effects on an Octa-Hydrated Curium Complex: An Experimental and Theoretical Investigation

An octa-hydrated curium compound [Cm(H 2 O) 8 ](Hdtp)(dtp)·H 2 O (Cm1,H 2 dtp = 2,3-di(tetrazol-5-yl)pyrazine) along with its lanthanide analogues [Ln(H 2 O) 8 ](Hdtp)-(dtp)·H 2 O (Ln1, Ln 3+ = La 3+ –Nd 3+ , Sm 3+ –Lu 3+ ) were synthesized and characterized using single crystal X-ray diffraction and spectroscopic methods. Bond length analysis of VIII Cm(III)–OH 2 (where VIII refers to the coordination number) was compared to VIII Ln(III)–OH 2 (Ln 3+ = Nd 3+ and Sm 3+ ), indicating similar VIII M(III)–OH 2 bond lengths owing to the similar eight-coordinate ionic radii of these VIII M(III) cations. Owing to the reduced coordination number, the VIII Cm(III)–OH 2 bond lengths were shorter than previously reported IX Cm(III)–OH 2 bonds in [Cm(H 2 O) 9 ](CF 3 SO 3 ) 3 . The octa-aquo complexes were also characterized by solid-state UV–vis–NIR spectroscopy in addition to variable-temperature and variable-pressure photoluminescence. Variable-pressure absorption spectra of Cm1 were compared with Ln1 and show that the Cm(III) f → f transitions have a stronger dependence on pressure than that observed in Ln1 (Ln 3+ = Nd 3+ and Sm 3+ ). The experimental and computational analyses reveal that the monotonic decrease in the computed energy difference between the ground state and the first excited state corresponds to the observed red shift of the photoluminescence peak. This is accompanied by a gradual reduction in the average Cm(III)–OH 2 bond length and a delocalization of spin densities, alongside an intensified interaction involving the 5f orbitals under increasing pressure. Furthermore, these changes accommodate the new geometry and collectively modify the energy landscape, resulting in peak broadening and quenching.

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

Non-Equilibrium Actinide Radiation Chemistry and the Nuclear Fuel Cycle

Actinides are inherently unstable elements that frequently coexist with other radioisotopes, generating intense ionizing radiation fields that drive the formation of non equilibrium oxidation states. These transient species exert a profound mechanistic influence on the radiation response of actinide containing systems due to their unique redox chemistry. Despite their importance, they remain poorly understood, yet such insight is essential for advancing actinide science and accurately predicting radiation driven behavior. Actinide separations—critical for nuclear energy technologies, strategic deterrence, space exploration, and nuclear medicine—depend on precise control of actinide oxidation states to recover targeted elements from complex matrices such as used nuclear fuel. However, during these processes, actinides, their coordination complexes, and the separation media are all exposed to intense, multicomponent (alpha, beta, gamma, etc.) radiation fields that can alter process efficiency, selectivity, and chemical stability. Understanding, controlling, and mitigating radiation induced reactions is therefore key to innovating and optimizing next generation separation technologies. This seminar will provide an overview of the nuclear fuel cycle and non equilibrium actinide radiation chemistry in the context of recovering actinides from used nuclear fuel, with a particular emphasis on direct dissolution–based reprocessing strategies. We will explore time resolved electron pulse radiolysis and alpha and gamma dose accumulation studies, integrated with multiscale computational modeling, to elucidate the molecular level roles of radiation driven, non equilibrium actinide species in process performance and in the radiolytic stability of organic ligands used for actinide recovery. These insights offer new pathways for designing advanced separation methods and next generation solvent systems, with broad implications for the future of the nuclear fuel cycle.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2021

On the eigenvector bias of Fourier feature networks: From regression to solving multi-scale PDEs with physics-informed neural networks

Physics-informed neural networks (PINNs) are demonstrating remarkable promise in integrating physical models with gappy and noisy observational data, but they still struggle in cases where the target functions to be approximated exhibit high-frequency or multi-scale features. Here in this work we investigate this limitation through the lens of Neural Tangent Kernel (NTK) theory and elucidate how PINNs are biased towards learning functions along the dominant eigen-directions of their limiting NTK. Using this observation, we construct novel architectures that employ spatio-temporal and multi-scale random Fourier features, and justify how such coordinate embedding layers can lead to robust and accurate PINN models. Numerical examples are presented for several challenging cases where conventional PINN models fail, including wave propagation and reaction–diffusion dynamics, illustrating how the proposed methods can be used to effectively tackle both forward and inverse problems involving partial differential equations with multi-scale behavior.

42 ENGINEERING↗
DOE OSTI2024

Smooth Dispersion Is Physically Appropriate: Assessing and Amending the D4 Dispersion Model

The addition of dispersion corrections to density functionals is essential for accurate energy and geometry predictions. Among them, the D4 scheme is popular due to its low computational cost and high accuracy. However, due to its design, the D4 correction can occasionally lead to anomalies, such as unphysical curvature and bumps in the potential energy surface. We find these anomalies are common in the D4 model, although observable consequences are rarer than in the D3 model for reasons we explain. Nevertheless, we uncover instances of unphysical local minima and stationary points with the D4 scheme and propose two solutions that yield smoother dispersion energy as a function of nuclear position. One is trivial to implement, based on a smoother reparametrization of Gaussian weighting (D4S) to find the effective coordination number. The other replaces Gaussian weighting with soft linear interpolation (D4SL). These new approaches usually remove artificial extremum points, while maintaining accuracy.

Tkachenko, Nikolay V↗
DOE OSTI2022

Reinforcement Learning for Volt- Var Control: A Novel Two-stage Progressive Training Strategy

This paper develops a reinforcement learning (RL) approach to solve a cooperative, multi-agent Volt-Var Control (VVC) problem for high solar penetration distribution systems. The ingenuity of our RL method lies in a novel two-stage progressive training strategy that can effectively improve training speed and convergence of the machine learning algorithm. In Stage 1 (individual training), while holding all the other agents inactive, we separately train each agent to obtain its own optimal VVC actions in the action space: fconsume, generate, do-nothingg. In Stage 2 (cooperative training), all agents are trained again coordinatively to share VVC responsibility. Rewards and costs in our RL scheme include (i) a system-level reward (for taking an action), (ii) an agent-level reward (for doing-nothing), and (iii) an agent-level action cost function. This new framework allows rewards to be dynamically allocated to each agent based on their contribution while accounting for the trade-off between control effectiveness and action cost. The proposed methodology is tested and validated in a modified IEEE 123-bus system using realistic PV and load profiles. Simulation results confirm that the proposed approach is robust and computationally efficient; and it achieves desirable volt-var control performance under a wide range of operation conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Automated Development of Molten Salt Machine Learning Potentials: Application to LiCl

The in silico modeling of molten salts is critical for emerging "carbon-free" energy applications but is inhibited by the cost of quantum mechanically treating the high polarizabilities of molten salts. Here, we integrate configurational sampling using classical force fields with active learning to automate and accelerate the generation of Gaussian approximation potentials (GAP) for molten salts. This methodology reduces the number of expensive ab initio evaluations required for training set generation to O(100), enabling the facile parametrization of a molten LiCl GAP model that exhibits a 19 000-fold speedup relative to AIMD. The developed molten LiCl GAP model is applied to sample extended spatiotemporal scales, permitting new physical insights into molten LiCl's coordination structure as well as experimentally validated predictions of structures, densities, self-diffusion constants, and ionic conductivities. The developed methodology significantly lowers the barrier to the in silico understanding and design of molten salts across the periodic table.

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

User-Oriented Improvements in the MOOSE framework in support of Multiphysics Simulation

The MOOSE Framework is a foundational capability used by NEAMS to create over 15 different simulation tools for advanced nuclear reactors. Due to this ubiquity, improvements to the framework in support of modeling and simulation goals are critical to the program. These improvements can take many forms including optimization, improved user experience, streamlined APIs, parallelism, and new capability. The work transcribed in the report was in direct support of the simulation tools and is already deployed or will be deployed in the coming months. The capabilities implemented were, in the same order as this report, chainable execution objects or executors, support for transfers between applications at the same level in a coupling scheme, support for boundary/subdomain restricted transfers, support for transfers between applications with different coordinate or unit systems, support for MOOSE applications in the NEAMS workbench, deployment of MOOSE application of the INL HPC OnDemand platform, addition of a triangular meshing library in libMesh and increased support of face variables.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Steric and Lewis Basicity Influence of the Second Coordination Sphere on Electrocatalytic CO 2 Reduction by Manganese Bipyridyl Complexes

This study aims to provide a greater insight into the balance between steric (bpy vs (Ph) 2 bpy vs mes 2 bpy ligands) and Lewis basic ((Ph) 2 bpy vs (MeOPh) 2 bpy vs (MeSPh) 2 bpy ligands) influence on the efficiencies of the protonation-first vs reduction-first CO 2 reduction mechanisms with [Mn I (R 2 bpy)(CO) 3 (CH 3 CN)] + precatalysts, and on their respective transition-state geometries/energies for rate-determining C–OH bond cleavage toward CO evolution. The presence of only modest steric bulk at the 6,6'-diphenyl-2,2'-bipyridyl ((Ph) 2 bpy) ligand has here allowed unique insight into the mechanism of catalyst activation and CO 2 binding by navigating a perfect medium between the nonsterically encumbered bpy-based and the highly sterically encumbered mes2bpy-based precatalysts. Cyclic voltammetry conducted in CO 2 -saturated electrolyte for the (Ph) 2 bpy-based precatalyst [2-CH 3 CN] + confirms that CO 2 binding occurs at the two-electron-reduced activated catalyst [2] – in the absence of an excess proton source, in contrast to prior assumptions that all manganese catalysts require a strong acid for CO 2 binding. This observation is supported by computed free energies of the parent–child reaction for [Mn–Mn] o dimer formation, where increased steric hindrance relative to the bpy-based precatalyst correlates with favorable CO 2 binding. A critical balance must be adhered to, however, as the absence of steric bulk in the bpy-based precatalyst [1-CH 3 CN] + maintains a lower overpotential than [2-CH 3 CN] + at the protonation-first pathway with comparable kinetic performance, whereas an ~2-fold greater TOF max is observed at its reduction-first pathway with an almost identical overpotential as [2-CH 3 CN] + . Notably, excessive steric bulk in the mes2bpy-based precatalyst [3-CH 3 CN] + results in increased activation free energies of the C–OH bond cleavage transition states for both the protonation-first and the reduction-first pathways relative to both [1-CH 3 CN] + and [2-CH 3 CN] + . In fact, [3-CH 3 CN] + requires a 1 V window beyond its onset potential to reach its peak catalytic current, which is in contrast to the narrower (<0.30 V) potential response window of the remaining catalysts here studied. Furthermore, voltammetry recorded under 1 atm of CO 2 with 2.8 M (5%) H 2 O establishes [2-CH 3 CN] + to have the lowest overpotential (η = 0.75 V) in the series here studied, attributed to its ability to lie “on the fence” when providing sufficient steric bulk to hinder (but not prevent) [Mn–Mn] o dimerization, while simultaneously having a limited steric impact on the free energy of activation for the rate-determining C–OH bond cleavage transition state. While the methoxyphenyl bpy-based precatalyst [4-CH 3 CN] + possesses an increased steric presence relative to [2-CH 3 CN] + , this is offset by its capacity to stabilize the C–OH bond cleavage transition states of both the protonation-first and the reduction-first pathways by facilitating second coordination sphere H-bonding stabilization.

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

Real-time equation-of-motion CC cumulant and CC Green’s function simulations of photoemission spectra of water and water dimer

Newly developed coupled-cluster (CC) methods enable simulations of ionization potentials and spectral functions of molecular systems in a wide range of energy scales ranging from core-binding to valence. This paper discusses the results obtained with the real-time equation-of-motion CC cumulant (RT-EOM-CC) approach and CC Green’s function (CCGF) approaches in applications to the water and water dimer molecules. Here, we compare the ionization potentials obtained with these methods for the valence region with the results obtained with the coupled-cluster with singles, doubles, and perturbative triples formulation as a difference of energies for N and N – 1 electron systems. All methods show good agreement with each other. They also agree well with the experiment with errors usually below 0.1 eV for the ionization potentials. We also analyze unique features of the spectral functions, associated with the position of satellite peaks, obtained with the RT-EOM-CC and CCGF methods employing single and double excitations, as a function of the monomer OH bond length and the proton transfer coordinate in the dimer. Finally, we analyze the impact of the basis set effects on the quality of calculated ionization potentials and find that the basis set effects are less pronounced for the augmented-type sets.

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

Behavior of solitary waves of coupled nonlinear Schrödinger equations subjected to complex external periodic potentials with odd-$\mathcal{PT}$ symmetry

In this work, we discuss the response of both moving and trapped solitary wave solutions of a two-component nonlinear Schrödinger system in 1 + 1 dimensions to an odd-$\mathcal{PT}$ external periodic complex potential. The dynamical behavior of perturbed solitary waves is explored by conducting numerical simulations of the nonlinear system and using a collective coordinate variational approximation. We present case examples corresponding to choices of parameter values and initial conditions involved therein. The results of the collective coordinate approximation are compared against numerical simulations where we observe qualitatively good agreement between the two. Unlike the case for a single-component solitary wave in a complex periodic $\mathcal{PT}$-symmetric potential, the collective coordinate equations do not have a small oscillation regime, and initially the height of the two components changes in opposite directions often causing instability. We find that the dynamic stability criteria we have used in the one-component case are a good indicator for the onset of dynamic instabilities in the present setup.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Classification of Orbits in Poincare Maps Using Machine Learning

Poincare plots, also called Poincare maps, are used by plasma physicists to understand the behavior of magnetically confined plasma in numerical simulations of a tokamak. These plots are created by the intersection of field lines with a two-dimensional poloidal plane that is perpendicular to the axis of the torus representing the tokamak. A plot is composed of multiple orbits, each created by a different field line as it goes around the torus. Each orbit can have one of four distinct shapes, or classes, that indicate changes in the topology of the magnetic fields confining the plasma. Given the (x, y) coordinates of the points that form an orbit, the analysis task is to assign a class to the orbit, a task that appears ideally suited for a machine learning approach. In this paper, we describe how we overcame two major challenges in solving this problem - creating a high-quality training set, with few mislabeled orbits, and converting the coordinates of the points into features that are discriminating, despite the variation within the orbits of a class and the apparent similarities between orbits of different classes. Our automated approach is not only more objective and accurate than visual classification, but is also less tedious, making it easier for plasma physicists to analyze the topology of magnetic fields from numerical simulations of the tokamak.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Efficient Failure Probability Calculations and Modeling Interface Debonding in TRISO Particles with BISON

This document constitutes completion of the NEAMS milestone, which is titled Demonstrate efficiency improvements in TRISO failure probability calculations and the effect of debonding on fission product transport in TRISO fuel particles in Bison . In this report we present the development of the following: (1) the anisotropic elasticity model of pyrolytic carbon and extended material models using the local coordinate system for aspherical particle geometry; (2) a capability of modeling interface debonding in TRISO particles; and (3) an efficient high-fidelity approach to calculate the failure probability of TRISO particles. Benchmark problems and AGR-2 tests are used to verify the models and demonstrate these new capabilities.

36 MATERIALS SCIENCE↗