Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computing methods”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Efficient high-fidelity TRISO statistical failure analysis using Bison: Applications to AGR-2 irradiation testing

The ability of tri-structural isotropic (TRISO) fuel to contain fission products is largely dictated by the quality of the manufacturing process, since most of the fission product release is expected to occur due to coating layer failure in a small number of particles containing defects. The Bison fuel performance code has capabilities to predict failure in individual particles, accounting for the presence of defects, and to apply statistical analysis methods to compute the probability of failure in a set of fuel particles. Bison has recently undergone significant development both to improve its physical representations of fuel particle behavior and to improve the efficiency of its statistical failure calculations. Physical model improvements include new capabilities to account for the pressure generated by fission gases on inner pyrolytic carbon (IPyC) crack surfaces and to use local material coordinate orientation to accurately incorporate the anisotropy in the material properties in aspherical particles. To improve statistical modeling efficiency, a direct integration approach which involves directly integrating the failure probability function associated with statistically varying parameters has been developed. The direct integration approach is much more efficient than the Monte Carlo (MC) schemes commonly employed, and allows Bison to directly run high-dimensional fuel performance models, which improves the accuracy of failure probability calculations. Finally, a set of benchmark problems is considered here to compare the MC and direct integration approaches, and a statistical failure analysis of compacts in the Advanced Gas Reactor (AGR)-2 experiments is performed using the direct integration approach.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Affinity propagation clustering of full-field, high-spatial-dimensional measurements for robust output-only modal identification: A proof-of-concept study

Determination of the model order is a challenging problem in system identification, especially in output-only or operational modal identification where some modes are weakly excited. While existing methods such as the stabilization diagram method (spectral information) are effective, they do not scale to high-dimensional data, which is usually needed for high-fidelity characterization of structural dynamics and has been made available in the emerging full-field measurement techniques using optical methods such as photogrammetry and laser vibrometers. In this proof-of-concept study we present a new non-parametric, data-driven approach for robust output-only identification of high-spatial-dimensional modal parameters of basic structures by efficiently processing and interactively exploiting the full-field measurement (i.e., very dense spatial measurement points). Specifically, we first over-estimate the system model once, producing a pool of candidate modes associated with their modal frequencies and full-field, high-spatial-dimensional mode shapes. This is accomplished by a data-driven method termed affinity propagation clustering (APC), where the active clusters, which are the active modes in our formulations, emerge from the "message-passing" procedure and does not require a pre-determination of the cluster number (mode or model order). Next, rather than using the spectral information to distinguish the physical and spurious modes in the stabilization diagram method, we exploit and visualize the spatial, full-field mode shape associated with each candidate mode to do so. We conduct extensive experiments on basic structural models with comparisons to a few existing methods. The results indicate that the new method is computationally efficient for identifying high-spatial-dimensional modal parameters, and robust to identify weak modes by exploiting the full-field measurement. We also discuss its applicability and limitations for structures with complex geometry (shapes).

42 ENGINEERING↗
DOE OSTI2024

The carbon challenge: Design, synthesis, and chemisorption behavior of solid sorbents in direct air capture of carbon dioxide

Direct air capture (DAC) of CO 2 is a promising solution for reducing the carbon footprint through "negative emission" technology. However, the low CO 2 concentration (~400 ppm) and the dynamic nature of DAC processes present challenges in designing effective sorbent systems. Recent advancements in material design and structural engineering have led to the development of high-performance solid sorbents, offering a more stable, safe, and energy-efficient alternative to traditional liquid CO 2 capture methods. This review highlights progress in solid sorbent-based DAC, focusing on amine-modified materials, hydrogen-bonded frameworks, and ionic liquid-engineered scaffolds. The discussion covers design principles, synthesis methodologies, and their impact on CO 2 chemisorption, comparing the advantages and limitations of each approach. Characterization techniques, especially operando methods and computational tools, are reviewed to understand sorbent behavior during CO 2 integration and release. The reaction pathways and interaction mechanisms of these sorbents with CO 2 are analyzed to guide future design. Additionally, the CO 2 chemisorption behaviors, including capacity, sorption kinetics, recyclability, and durability in the presence of gaseous impurities and under humid conditions will be evaluated and compared. Further, the review offers unique insights into the physical properties, chemical structures, and surface engineering effects of these sorbents, based on comprehensive characterization and evaluation techniques.

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

Tutorial: Extracting entanglement signatures from neutron spectroscopy

This tutorial is a pedagogical introduction to recent methods of computing quantum spin entanglement witnesses from spectroscopy, with a special focus on neutron scattering on quantum spin systems. We offer a brief introduction to the concepts and equations, define a data analysis protocol, and discuss the interpretation of three entanglement witnesses: one-tangle, two-tangle, and Quantum Fisher Information. We also discuss practical experimental considerations, and give three examples of extracting entanglement witnesses from experimental data: Copper Nitrate, KCuF 3 , and NiPS 3 .

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Monte Carlo simulations for understanding the transport properties of metallic nanowires

We report the conductance calculation of metallic nanowires of different metals under the application of a tensile force. The elongation induced by the tensile force was simulated by means of the canonical Monte Carlo method combined with the embedded atom method to compute the energy of the system. We studied nanowires composed of Ag, Au, Pt and Cu. The conductance was calculated by means of the Landauer formula and the non-equilibrium Green function formalism. As previously verified experimentally and using other simulation techniques, the conductance is quantized. In addition to this, we have found that there is a strong correlation between the magnitude of the conductance and the width of the nanowire at its thinnest point. Conductance histograms were analyzed for all the different metals.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Analysis of spin frustration in an Fe III 7 cluster using a combination of computational, experimental, and magnetostructural correlation methods

The synthesis, structure, and magnetic properties are reported for [Fe 7 O 3 (O 2 C t Bu) 9 (mda) 3 (H 2 O) 3 ] ( 1 ), where mdaH 2 is N -methyldiethanolamine. 1 was prepared from the reaction of [Fe 3 O(O 2 C t Bu) 6 (H 2 O) 3 ](NO 3 ) with mdaH 2 in a 1:~3 ratio in MeCN. The core of 1 consists of a central octahedral Fe III ion held within a non-planar Fe 6 loop by three μ 3 -O 2- and three μ 2 -RO - arms from the three mda 2- chelates. Variable-temperature dc and ac magnetic susceptibility studies revealed dominant antiferromagnetic coupling, leading to a ground state spin of S = 5 / 2 . The ground state was confirmed by a fit of magnetization data collected in the 0.1–7.0 T and 1.8–10.0 K ranges. The four Fe 2 pairwise exchange parameters ( J 1 - J 4 ) were estimated by independent methods: theoretical calculations using either broken symmetry energy differences (-46.3, -16.2, -3.9, and - 28.1 cm -1 , respectively) or Green’s function approximation methods (-41.4, -14.8, -13.2, and - 24.7 cm -1 ), and a magnetostructural correlation (MSC) previously developed for high nuclearity Fe III /O complexes (-39.5, -13.8, -6.7, and - 23.5 cm -1 ). Additionally, the J 1 - J 4 obtained from the MSC and theoretical methods were used with the program PHI to both simulate χ M T vs T as well as to serve as reasonable input values to fit the experimental data (-41.0, -11.4, -5.0, and - 27.3 cm -1 ). Analysis of the J ij led to identification of the spin frustration effects operative and the resultant spin vector alignments at each Fe III ion, thus allowing for the rationalization of the experimental ground state.

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

ab initio Sub-Mechanism Development for Cyclopentene Oxidation

To accurately predict low-temperature oxidation behavior, chemical kinetics mechanisms must contain complete reaction networks that include detailed consumption reactions of intermediates produced directly from hydroperoxyalkyl radicals, Q̇OOH, which undergo competing unimolecular reactions and bimolecular reactions with O2. Rates of chain-branching are governed by the flux between the two competing pathways, and inherently depend on temperature, pressure, and oxygen concentration. Neglect of consumption pathways for major oxidation intermediates leads to mechanism truncation error that is ameliorated by expanding the level of detail included in sub-mechanisms and employing ab initio methods for computing rates of elementary reactions and thermochemical properties of species involved. In the present work, an ab initio-derived sub-mechanism is developed using AutoMech to model the chemical kinetics of cyclopentene, a major product of cyclopentane oxidation. The ab initio sub-mechanism builds on a detailed mechanism developed using Reaction Mechanism Generator (RMG) for the specific purpose of determining the extent to which replacing cyclopentene-specific reactions and species with quantum chemical computations reduces model inaccuracies resulting from mechanism truncation error. In an effort to minimize interference from other reactions present during the formation of cyclopentene from cyclopentyl + O2, providing a narrower experimental scope, the model is compared against speciation measurements from jet-stirred reactor (JSR) experiments on cyclopentene oxidation. The experiments utilize vacuum ultraviolet-absorption spectroscopy and mass spectrometry for isomer-resolved speciation of intermediates at 835 Torr from 700 – 950 K. [O2]-dependent experiments were also conducted from 0.057 – 2.01 · 1018 molecules cm–3 at 825 K to examine the influence of oxygen on species profiles. Model predictions using the ab initio-revised mechanism yielded significant improvements in species profiles for both the temperature- and [O2]-dependent measurements, owing in part to increased rates of HOȮ and H2O2 production, which underscores the influence of theoretical calculations of reaction rates involving species produced from Ṙ + O2 such as cyclopentene.

AutoMech↗
DOE OSTI2020

Atomic resolution convergent beam electron diffraction analysis using convolutional neural networks

Two types of convolutional neural network (CNN) models, a discrete classification network and a continuous regression network, were trained to determine local sample thickness from convergent beam diffraction (CBED) patterns of SrTiO 3 collected in a scanning transmission electron microscope (STEM) at atomic column resolution. Acquisition of atomic resolution CBED patterns for this purpose requires careful balancing of CBED feature size in pixels, acquisition speed, and detector dynamic range. The training datasets were derived from multislice simulations, which must be convolved with incoherent source broadening. Sample thicknesses were also determined using quantitative high-angle annular dark-field (HAADF) STEM images acquired simultaneously. The regression CNN performed well on sample thinner than 35 nm, with 70% of the CNN results within 1 nm of HAADF thickness, and 1.0 nm overall root mean square error between the two measurements. The classification CNN was trained for a thicknesses up to 100 nm and yielded 66% of CNN results within one classification increment of 2 nm of HAADF thickness. Our approach depends on methods from computer vision including transfer learning and image augmentation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Computing the viscous effect in early-time drop impact dynamics

The impact of a liquid drop on a solid surface involves many intertwined physical effects, and is influenced by drop velocity, surface tension, ambient pressure and liquid viscosity, among others. Experiments by Kolinski et al. ( Phys. Rev. Lett. , vol. 112, no. 13, 2014 b , p. 134501) show that the liquid–air interface begins to deviate away from the solid surface even before contact. They found that the lift-off of the interface starts at a critical time that scales with the square root of the kinematic viscosity of the liquid. To understand this, we study the approach of a liquid drop towards a solid surface in the presence of an intervening gas layer. We take a numerical approach to solve the Navier–Stokes equations for the liquid, coupled to the compressible lubrication equations for the gas, in two dimensions. With this approach, we recover the experimentally captured early time effect of liquid viscosity on the drop impact, but our results show that lift-off time and liquid kinematic viscosity have a more complex dependence than the square-root scaling relationship. We also predict the effect of interfacial tension at the liquid–gas interface on the drop impact, showing that it mediates the lift-off behaviour.

42 ENGINEERING↗
DOE OSTI2021

Identification of Known and Novel Monomers for Poly(hydroxyurethanes) from Biobased Materials

Sourcing polymers from biobased materials is desirable for a transition to a more sustainable world. However, finding new monomers sourced from biobased molecules requires extensive experimental effort, due to the fact that most biomolecules lack the functional groups or properties needed to serve as a monomer in a particular polymeric system without chemical modification. In this work, we demonstrate a computational screening method to obtain monomers of a non-isocyanate polyurethane (NIPU), polyhydroxyurethane (PHU), from biomolecules. Our method generates candidate monomers from a list of 15 selected biobased starting molecules by performing chemical transformations that are standard to chemical catalysis and organic synthesis on these starting molecules and then screening the results for molecules with moieties that can be used as monomers for a PHU. Product molecules from two generations of reaction, a total number of 20,057, were analyzed automatically to identify molecules that have two or more epoxy or cyclic carbonate functional groups as candidates for further derivatization. These 21 candidate molecules were then subjected to another two generations of reaction, producing a set of 17,913 additional molecules. Further scrutiny of this set for desired functional groups characterizing PHUs and querying against the PubChem database was carried out, and 15 known molecules that are commercially available were identified, with the desirable features that they are predicted to be derived from biobased starting molecules in four or fewer steps. Moreover, molecule space generated from this process was primarily comprised of molecules that were not known to PubChem, indicating the framework also has the ability to generate novel candidate monomers that originate from biobased starting points.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Discovery of Novel Rhizoctonia solani DHFR Inhibitors as Fungicides Using Virtual Screening

Dihydrofolate reductase (DHFR) is an essential enzyme in the folate pathway and has been recognized as a well-known target for antibacterial and antifungal drugs. We discovered eight compounds from the ZINC database using virtual screening to inhibit Rhizoctonia solani (R. solani), a fungal pathogen in crops. These compounds were evaluated with in vitro assays for enzymatic and antifungal activity. Among these, compound Hit8 is the most active R. solani DHFR inhibitor, with the IC 50 of 10.2 μM. The selectivity of inhibition is 22.3 against human DHFR with the IC 50 of 227.7 μM. Moreover, Hit8 has higher antifungal activity against R. solani (EC 50 of 38.2 mg L –1 ) compared with validamycin A (EC 50 of 67.6 mg L –1 ), a well-documented fungicide. These results suggest that Hit8 may be a potential fungicide. Finally, our study exemplifies a computer-aided method to discover novel inhibitors that could target plant pathogenic fungi.

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

Toward DMC Accuracy Across Chemical Space with Scalable Δ-QML

In the past decade, quantum diffusion Monte Carlo (DMC) has been demonstrated to successfully predict the energetics and properties of a wide range of molecules and solids by numerically solving the electronic many-body Schrödinger equation. With O(N 3 ) scaling with the number of electrons N, DMC has the potential to be a reference method for larger systems that are not accessible to more traditional methods such as CCSD(T). Assessing the accuracy of DMC for smaller molecules becomes the stepping stone in making the method a reference for larger systems. We show that when coupled with quantum machine learning (QML)-based surrogate methods, the computational burden can be alleviated such that quantum Monte Carlo (QMC) shows clear potential to undergird the formation of high-quality descriptions across chemical space. We discuss three crucial approximations necessary to accomplish this: the fixed-node approximation, universal and accurate references for chemical bond dissociation energies, and scalable minimal amons-set-based QML (AQML) models. Numerical evidence presented includes converged DMC results for over 1000 small organic molecules with up to five heavy atoms used as amons and 50 medium-sized organic molecules with nine heavy atoms to validate the AQML predictions. Finally, numerical evidence collected for Δ-AQML models suggests that already modestly sized QMC training data sets of amons suffice to predict total energies with near chemical accuracy throughout chemical space.

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

Stochastic Real-Time Second-Order Green’s Function Theory for Neutral Excitations in Molecules and Nanostructures

We present a real-time second-order Green’s function (GF) method for computing excited states in molecules and nanostructures, with a computational scaling of O(N e 3 ), where N e is the number of electrons. The cubic scaling is achieved by adopting the stochastic resolution of the identity to decouple the 4-index electron repulsion integrals. To improve the time propagation and the spectral resolution, we adopt the dynamic mode decomposition technique and assess the accuracy and efficiency of the combined approach for a chain of hydrogen dimer molecules of different lengths. We find that the stochastic implementation accurately reproduces the deterministic results for the electronic dynamics and excitation energies. Furthermore, we provide a detailed analysis of the statistical errors, bias, and long-time extrapolation. Overall, the approach offers an efficient route to investigate excited states in extended systems with open or closed boundary conditions.

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

Mean-Field Ring Polymer Rates Using a Population Dividing Surface

Mean-field Ring Polymer Molecular Dynamics (MF-RPMD) offers a computationally efficient method for the simulation of reaction rates in multi-level systems. Previous work has established that, to model a nonadiabatic state-to-state reaction accurately, the dividing surface must be chosen to explicitly sample kinked ring polymer configurations where at least one bead is in a different electronic state than the others. Building on this, we introduce a population difference coordinate and a kink-constrained dividing surface, and we test the accuracy of the resulting mean-field rate theory on a series of bilinear coupling model systems. Here, we demonstrate that this new MF-RPMD rate approach is efficient to implement and quantitatively accurate for models over a wide range of driving forces, coupling strengths, and temperatures.

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

Characterizing Surface Ice-Philicity Using Molecular Simulations and Enhanced Sampling

The formation of ice, which plays an important role in diverse contexts ranging from cryopreservation to atmospheric science, is often mediated by solid surfaces. Although surfaces that interact favorably with ice (relative to liquid water) can facilitate ice formation by lowering nucleation barriers, the molecular characteristics that confer ice-philicity to a surface are complex and incompletely understood. To address this challenge, here we introduce a robust and computationally efficient method for characterizing surface ice-philicity that combines molecular simulations and enhanced sampling techniques to quantify the free energetic cost of increasing surface–ice contact at the expense of surface–water contact. Using this method to characterize the ice-philicity of a family of model surfaces that are lattice matched with ice but vary in their polarity, we find that the nonpolar surfaces are moderately ice-phobic, whereas the polar surfaces are highly ice-philic. In contrast, for surfaces that display no complementarity to the ice lattice, we find that ice-philicity is independent of surface polarity and that both nonpolar and polar surfaces are moderately ice-phobic. Furthermore, our work thus provides a prescription for quantitatively characterizing surface ice-philicity and sheds light on how ice-philicity is influenced by lattice matching and polarity.

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

Adsorbate Partition Functions via Phase Space Integration: Quantifying the Effect of Translational Anharmonicity on Thermodynamic Properties

A new method for computing anharmonic thermophysical properties for adsorbates on metal surfaces is presented. Classical Monte Carlo phase space integration is performed to calculate the partition function for the motion of a hydrogen atom on Cu(111). Here, a minima-preserving neural network potential energy surface is used within the integration routine. Two different sampling schema for generating the training data are presented, and two different density functionals are used. The results are benchmarked against direct state counting results by using discrete variable representation. The phase space integration results are in excellent quantitative agreement with the benchmark results. Additionally, both the discrete variable representation and the phase space integration results confirm that the motion of H on Cu(111) is highly anharmonic. The results were applied to calculate the free energy of dissociative adsorption of H 2 and the resulting Langmuir isotherms at 400, 800, and 1200 K in a partial pressure range of 0–1 bar. It shows that the anharmonic effects lead to significantly higher predicted surface site fractions of hydrogen.

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

Accurate and Efficient Prediction of p K w in Aqueous Electrolytes Using Local Electrostatic Potentials

The K w =10 -pK w in aqueous electrolytes exhibit significant variations as a function of ion concentration. This behavior is not well-described by the Pitzer model, particularly at high concentrations. Here, we provide a molecular interpretation of the concentration dependent trends in p K w and develop a new method that accurately predicts the dissociation constant based upon the local electrostatic potential of the water oxygen that is determined by its nearest neighbors. The method is computationally efficient, relying only upon sampling of the local configurations (via molecular dynamics) and a small set of experimentally measured p K w . It is extensible to a range of salt compositions and the molecular understanding of the local potentials provides new routes toward the development of electrolytes with tailored physicochemical properties like p K w .

anions↗