Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computing Methodologies”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Application of high-speed, species-specific chemiluminescence imaging for laminar flame speed and Markstein length measurements in spherically expanding flames

A high-speed chemiluminescence imaging diagnostic system has been developed and characterized for measuring species-specific properties of spherically expanding flames, with the objective of obtaining synchronized measurements of multiple flame properties from a single flame experiment. The imaging system comprises a high-speed camera coupled to a high-speed image intensifier and optical filters to isolate emission from desired excited species. Here this study validates the diagnostic scheme specifically for measuring laminar flame speed and Markstein length from chemiluminescence at three wavelength bands corresponding to OH*, CH*, and total chemiluminescence, respectively. Mixtures of methane-air and methane/ethane-air were tested at initial conditions of atmospheric pressure and room temperature. An under-utilized methodology for computing unburned Markstein length from burned Markstein length was implemented. Flame speed modeling of each mixture was conducted in Chemkin using AramcoMech 2.0. No significant differences were observed between the results derived from the three emission wavelength bands for the flame speeds nor the Markstein lengths. For both of these properties, excellent agreement was observed between the experimental results and the values reported in the literature. Because laminar flame speed and Markstein length show no dependency on the imaging wavelength, the present results imply that multiple flame properties, such as laminar flame speed, Markstein length, flame thickness, and excited-species intensity profiles, can all be acquired simultaneously for the same flame using this method.

42 ENGINEERING↗
DOE OSTI2015

Generalized gas-solid adsorption modeling: Single-component equilibria

Over the last several decades, modeling of gas–solid adsorption at equilibrium has generally been accomplished through the use of isotherms such as the Freundlich, Langmuir, Tóth, and other similar models. While these models are relatively easy to adapt for describing experimental data, their simplicity limits their generality to be used with many different sets of data. This limitation forces engineers and scientists to test each different model in order to evaluate which one can best describe their data. Additionally, the parameters of these models all have a different physical interpretation, which may have an effect on how they can be further extended into kinetic, thermodynamic, and/or mass transfer models for engineering applications. Therefore, it is paramount to adopt not only a more general isotherm model, but also a concise methodology to reliably optimize for and obtain the parameters of that model. A model of particular interest is the Generalized Statistical Thermodynamic Adsorption (GSTA) isotherm. The GSTA isotherm has enormous flexibility, which could potentially be used to describe a variety of different adsorption systems, but utilizing this model can be fairly difficult due to that flexibility. To circumvent this complication, a comprehensive methodology and computer code has been developed that can perform a full equilibrium analysis of adsorption data for any gas-solid system using the GSTA model. The code has been developed in C/C++ and utilizes a Levenberg–Marquardt’s algorithm to handle the non-linear optimization of the model parameters. Since the GSTA model has an adjustable number of parameters, the code iteratively goes through all number of plausible parameters for each data set and then returns the best solution based on a set of scrutiny criteria. Data sets at different temperatures are analyzed serially and then linear correlations with temperature are made for the parameters of the model. The end result is a full set of optimal GSTA parameters, both dimensional and non-dimensional, as well as the corresponding thermodynamic parameters necessary to predict the behavior of the system at temperatures for which data were not available. It will be shown that this code, utilizing the GSTA model, was able to describe a wide variety of gas-solid adsorption systems at equilibrium.In addition, a physical interpretation of these results will be provided, as well as an alternate derivation of the GSTA model, which intends to reaffirm the physical meaning.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Incorporating Coverage-Dependent Reaction Barriers into First-Principles-Based Microkinetic Models: Approaches and Challenges

Mean-field microkinetic models (MKMs) are appealing for their relatively facile construction, computational tractability, and high-throughput catalyst screening capabilities. As such, they will continue to be a valuable tool for materials design in heterogeneous catalysis even as the field aims to describe more complex systems. Numerous prior reports have provided the groundwork for constructing first-principles-based MKMs, including the analysis of strategies for incorporating lateral interactions into thermodynamic parameters (e.g., adsorption energies). Yet, there remains a need for concerted dialogue on methods for calculating and incorporating coverage-dependent kinetic parameters into MKMs. In this Perspective, we assess strategies for doing so, including the corresponding key physical implications and computational challenges. Here, we emphasize that decoupling thermodynamic and kinetic parameters within MKMs can violate thermodynamic consistency and risk unphysical solutions. For some reactions and catalyst materials, scaling relationships can predict coverage-dependent activation energies, but there are several exceptions evident in the literature, indicating that this approach is not universally applicable and that the field could benefit from research aimed at elucidating the limitations. Conducting high-coverage transition state searches is a rigorous but computationally costly strategy, and the effects of various methods for mitigating this cost on resulting energetics have yet to be broadly explored and validated. The goal of this Perspective is to generate discussion on and inspire focused research into the physical relevance of approaches for describing coverage-dependent reaction barriers in MKMs, including the development of computationally tractable methodologies, to advance the applicability of MKMs across diverse reaction chemistries and conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Quantification of flexoelectricity in PbTiO 3 /SrTiO 3 superlattice polar vortices using machine learning and phase-field modeling

Flexoelectricity refers to electric polarization generated by heterogeneous mechanical strains, namely strain gradients, in materials of arbitrary crystal symmetries. Despite more than 50 years of work on this effect, an accurate identification of its coupling strength remains an experimental challenge for most materials, which impedes its wide recognition. Here, we show the presence of flexoelectricity in the recently discovered polar vortices in PbTiO 3 /SrTiO 3 superlattices based on a combination of machine-learning analysis of the atomic-scale electron microscopy imaging data and phenomenological phase-field modeling. By scrutinizing the influence of flexocoupling on the global vortex structure, we match theory and experiment using computer vision methodologies to determine the flexoelectric coefficients for PbTiO 3 and SrTiO 3 . Here, our findings highlight the inherent, nontrivial role of flexoelectricity in the generation of emergent complex polarization morphologies and demonstrate a viable approach to delineating this effect, conducive to the deeper exploration of both topics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ab initio-based metric for predicting the protectiveness of surface films in aqueous media

Abstract Materials can passivate by forming surface films when placed in aqueous media. However, these films may or may not be stable, and their stability can be predicted by a metric called the Pilling-Bedworth Ratio (PBR). In this article, we extend PBR to predict passivation protectiveness of multi-component materials. We then evaluate this PBR (ePBR)’s effectiveness by comparing its predictions against experimental studies of 21 multi-element materials of diverse chemistries, with agreement for 17 of the materials. Finally, we encode the methodology to compute ePBR in a web-application to predict the protectiveness of 140,000+ materials in the Materials Project database.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Stress in titania nanoparticles: An atomistic study

Stress engineering is becoming an increasingly important method for controlling electronic, optical, and magnetic properties of nanostructures, although the concept of stress is poorly defined at the nanoscale. We outline a methodology for computing bulk and surface stress in nanoparticles using atomistic simulation. The method is applicable to ionic and non- ionic materials alike and may be extended to other nanostructures. We apply it to spherical anatase nanoparticles ranging from 2 to 6 nm in diameter and obtain a surface stress of 0.89 N/m, in agreement with experimental measurements. Based on the extent that stress inhomogeneities at the surface are transmitted into the bulk, two characteristic length-scales are identified: below 3 nm bulk and surface regions cannot be defined and the available analytic theories for stress are not applicable, and above about 5 nm the stress becomes well-described by the theoretical Young-Laplace equation. The effect of a net surface charge on the bulk stress is also investigated. It is found that moderate surface charges can induce significant bulk stresses, on the order of 100 MPa, in nanoparticles within this size range.

molecular dynamics simulations↗
DOE OSTI2013

Improving the accuracy and efficiency of time-resolved electronic spectra calculations: Cellular dephasing representation with a prefactor

Time-resolved electronic spectra can be obtained as the Fourier transform of a special type of time correlation function known as fidelity amplitude, which, in turn, can be evaluated approximately and efficiently with the dephasing representation. Here we improve both the accuracy of this approximation—with an amplitude correction derived from the phase-space propagator—and its efficiency—with an improved cellular scheme employing inverse Weierstrass transform and optimal scaling of the cell size. We demonstrate the advantages of the new methodology by computing dispersed time-resolved stimulated emission spectra in the harmonic potential, pyrazine, and the NCO molecule. In contrast, we show that in strongly chaotic systems such as the quartic oscillator the original dephasing representation is more appropriate than either the cellular or prefactor-corrected methods.

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

Quadratic partial eigenvalue assignment in large-scale stochastic dynamic systems for resilient and economic design

Failure of structural systems under dynamic loading can be prevented via active vibration control which shifts the damped natural frequencies of the systems away from the dominant range of loading spectrum. The damped natural frequencies and the dynamic load typically show significant variations in practice. A computationally efficient methodology based on quadratic partial eigenvalue assignment technique and optimization under uncertainty has been formulated in the present work that will rigorously account for these variations and result in an economic and resilient design of structures. A novel scheme based on hierarchical clustering and importance sampling is also developed in this work for accurate and efficient estimation of probability of failure to guarantee the desired resilience level of the designed system. Numerical examples are presented to illustrate the proposed methodology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Bubble nucleation in simple and molecular liquids via the largest spherical cavity method

In this work, we propose a methodology to compute bubble nucleation free energy barriers using trajectories generated via molecular dynamics simulations. We follow the bubble nucleation process by means of a local order parameter, defined by the volume of the largest spherical cavity (LSC) formed in the nucleating trajectories. This order parameter simplifies considerably the monitoring of the nucleation events, as compared with the previous approaches which require ad hoc criteria to classify the atoms and molecules as liquid or vapor. The combination of the LSC and the mean first passage time technique can then be used to obtain the free energy curves. Upon computation of the cavity distribution function the nucleation rate and free-energy barrier can then be computed. We test our method against recent computations of bubble nucleation in simple liquids and water at negative pressures. We obtain free-energy barriers in good agreement with the previous works. The LSC method provides a versatile and computationally efficient route to estimate the volume of critical bubbles the nucleation rate and to compute bubble nucleation free-energies in both simple and molecular liquids.

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

Hadron structure from lattice QCD

Recent progress in lattice QCD calculations of nucleon structure will be presented. Calculations of nucleon matrix elements and form factors have long been difficult to reconcile with experiment, but with advances in both methodology and computing resources, this situation is improving. Some calculations have produced agreement with experiment for key observables such as the axial charge and electromagnetic form factors, and the improved understanding of systematic errors will help to increase confidence in predictions of unmeasured quantities. The long-omitted disconnected contributions are now seeing considerable attention and some recent calculations of them will be discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Calculating shock Hugoniot and isentropes using multiphase equation of state tables and application to shock and release of diamond ablators in inertial confinement fusion implosions

Advances in shock and ramp compression techniques now allow experimental access to unprecedented extreme conditions of pressure and temperature, providing a means to test theoretical models. Here, we describe a simple methodology to compute multi-phase shock Hugoniot and isentropes using multiphase equation of state tables. We treat explicitly the phase coexistence along the phase boundary to reveal the evolution of the sample as it undergoes the phase transformation in adiabatic conditions. We illustrate the method by calculating the predicted shock and shock-and-release behavior of diamond at conditions relevant for the initial stage of inertial confinement fusion implosions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Uncovering universal rules governing the selectivity of the archetypal DNA glycosylase TDG

Thymine DNA glycosylase (TDG) is a pivotal enzyme with dual roles in both genome maintenance and epigenetic regulation. TDG is involved in cytosine demethylation at CpG sites in DNA. Here we have used molecular modeling to delineate the lesion search and DNA base interrogation mechanisms of TDG. First, we examined the capacity of TDG to interrogate not only DNA substrates with 5-carboxyl cytosine modifications but also G:T mismatches and nonmismatched (A:T) base pairs using classical and accelerated molecular dynamics. To determine the kinetics, we constructed Markov state models. Base interrogation was found to be highly stochastic and proceeded through insertion of an arginine-containing loop into the DNA minor groove to transiently disrupt Watson–Crick pairing. Next, we employed chain-of-replicas path-sampling methodologies to compute minimum free energy paths for TDG base extrusion. We identified the key intermediates imparting selectivity and determined effective free energy profiles for the lesion search and base extrusion into the TDG active site. Here, our results show that DNA sculpting, dynamic glycosylase interactions, and stabilizing contacts collectively provide a powerful mechanism for the detection and discrimination of modified bases and epigenetic marks in DNA.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Organ and effective dose coefficients for cranial and caudal irradiation geometries: Neutrons

Dose coefficients based on the recommendations of International Commission on Radiological Protection (ICRP) Publication 103 were reported in ICRP Publication 116, the revision of ICRP Publication 74 and ICRU Publication 57 for the six reference irradiation geometries: anterior–posterior, posterior–anterior, right and left lateral, rotational and isotropic. In this work, dose coefficients for neutron irradiation of the body with parallel beams directed upward from below the feet (caudal) and downward from above the head (cranial) using the ICRP 103 methodology were computed using the MCNP 6.1 radiation transport code. The dose coefficients were determined for neutrons ranging in energy from 10 –9 MeV to 10 GeV. Here, at energies below about 500 MeV, the cranial and caudal dose coefficients are less than those for the six reference geometries reported in ICRP Publication 116.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2019

Impact of phonons and spin-orbit coupling on Auger recombination in InAs

Here we extend the methodology for calculating Auger recombination from first principles to include spin-orbit interaction, which allows us to study materials in which the spin-orbit splitting and the band gap are of similar magnitude. We use this methodology to compute the direct and indirect Auger coefficients in InAs and related alloys. The direct process involves only Coulomb interaction, while the indirect process is mediated by phonons. We show how the spin-orbit split-off states allow for efficient direct Auger recombination, greatly enhancing the excitation rate of Auger holes due to a near resonance between the spin-orbit splitting and the band gap of the material. We find that the direct hole-hole-electron Auger recombination coefficient decays exponentially with increasing band gap, while the electron-electron-hole Auger recombination coefficient decays faster than exponentially. Additionally, we modeled the phonon-assisted process in InAs and show that the direct Auger process is dominant.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Geometrically necessary dislocation fingerprints of dislocation loop absorption at grain boundaries

Here we present a numerical methodology to compute the Nye-tensor fingerprints of dislocation loop absorption at grain boundaries (GBs) for comparison with TEM observations of irradiated polycrystals. Our approach links atomistic simulations of self-interstitial atom (SIA) prismatic loops gliding toward and interacting with GBs in body-centered cubic iron with experimentally extracted geometrically necessary dislocation (GND) maps to facilitate the interpretation of damage processes. The Nye-tensor analysis is strongly mesh-size dependent—corresponding to resolution-dependent TEM observations. The method computes GND fingerprints from discretized dislocation line segments extracted from molecular dynamics simulations of dislocation loops being absorbed at a GB. Specifically, we perform MD simulation of prismatic loops of two diameters and monitor the three stages of the absorption process: loop glide, the partial, and full absorption of the loops at a [1 0 0] symmetric tilt GB. These methods provide a framework for future investigations of the nature of defect absorption by grain boundaries under irradiation conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

First-principles determination of the phonon-point defect scattering and thermal transport due to fission products in ThO 2

This work presents the first-principles calculations of the lattice thermal conductivity degradation due to point defects in thorium dioxide using an iterative solution of the Pierels-Boltzmann transport equation. We have used the non perturbative Green’s function methodology to compute the phonon- point defect scattering rates that consider the local distortion around the point defect, including the mass difference changes, interatomic force constants, and structural relaxation near the point defects. The point defects considered in this work include the vacancy of thorium ($V$ $TH$ ) and oxygen ($V$ $O$ ), substitutions of helium ($He$ $Th$ ), krypton ($Kr$ $Th$ ), zirconium ($Zr$ $Th$ ), iodine ($I$ $Th$ ) and xenon ($Xe$ $Th$ ) in the thorium site, and the three different configurations of the Schottky defects. The results of the phonon-defect scattering rate reveal that, among all the considered intrinsic defects, the thorium vacancy and helium substitution in the thorium site scatter the phonon most due to the substantial changes in the force constant and structural distortions. The scattering of phonons due to the substitutional defects unveils that the zirconium atom scatters phonons the least, followed by xenon, iodine, krypton, and helium. This is contrary to the intuition that the scattering strength follows $He$ $Th$ > $Kr$ $Th$ > $Zr$ $Th$ > $I$ $Th$ > $Xe$ $Th$ based on the mass difference. This striking difference in the zirconium phonon scattering is due to the local chemical environment changes. Zirconium is an electropositive element with valency similar to thorium and, therefore, can bond with the oxygen atoms, thus creating less force constant variance compared to iodine, an electronegative element noble gas helium, xenon, and krypton. In conclusion, these results can serve as the benchmark for the analytical models and help the engineering-scale modeling effort for nuclear design.

36 MATERIALS SCIENCE↗