Search NASA⌕ Search

SEARCH · Search NASA

Results for “calculation method”

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 1,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Maximum von Mises Stress in the Loading Environment of Mass Acceleration Curve

Method for calculating stress due to acceleration loading: 1) Part has been designed by FEA and hand calculation in one critical loading direction judged by the analyst; 2) Maximum stress can be due to loading in another direction; 3) Analysis procedure to be presented determines: a) The maximum Mises stress at any point; and b) The direction of maximum loading associated with the "stress". Concept of Mass Acceleration Curves (MAC): 1) Developed by JPL to perform preliminary structural sizing (i.e. Mariners, Voyager, Galileo, Pathfinder, MER,...MSL); 2) Acceleration of physical masses are bounded by a curve; 3) G-levels of vibro-acoustic and transient environments; 4) Convergent process before the couple loads cycle; and 5) Semi-empirical method to effectively bound the loads, not a simulation of the actual response.

von Mises stress↗
DOE OSTI2018

Dissection of the multichannel reaction of acetylene with atomic oxygen: from the global potential energy surface to rate coefficients and branching dynamics

The O( 3 P) + C 2 H 2 reaction is the first step in acetylene oxidation. The accurate kinetic data and the understanding of the reaction dynamics is of great importance. To this end, in this study a full-dimensional global potential energy surface (PES) for the ground triplet state of the O( 3 P) + C 2 H 2 reaction is constructed based on approximately 85,000 ab initio points calculated at the level of explicitly correlated unrestricted coupled cluster single, double, and perturbative triple excitations with the explicitly correlated polarized valence triple zeta basis set (UCCSD(T)-F12b/VTZ-F12). The PES is fit using the permutation invariant polynomial-neural network (PIP-NN) approach with a total root mean square error of 0.21 kcal mol –1 . The key topographic features of the PES, including multiple potential wells and saddle points along different reaction pathways, are well represented by this fit PES. The kinetics and dynamics of the O( 3 P) + C 2 H 2 reaction are investigated using the quasi-classical trajectory (QCT) method. The calculated rate coefficients are in good agreement with experimental data over a wide temperature range, especially when the temperature is lower than 1500 K. The product branch ratio has also been determined, which indicates the H + HCCO channel as the dominant reaction pathway at 298–3000 K, accounting for 80–90% of the overall rate coefficient, in agreement with experimental observations. The dynamics of the reaction is analyzed in detail.

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

Probabilistic Nanomagnetic Memories for Uncertain and Robust Machine Learning

This project evaluated the use of emerging spintronic memory devices for robust and efficient variational inference schemes. Variational inference (VI) schemes, which constrain the distribution for each weight to be a Gaussian distribution with a mean and standard deviation, are a tractable method for calculating posterior distributions of weights in a Bayesian neural network such that this neural network can also be trained using the powerful backpropagation algorithm. Our project focuses on domain-wall magnetic tunnel junctions (DW-MTJs), a powerful multi-functional spintronic synapse design that can achieve low power switching while also opening the pathway towards repeatable, analog operation using fabricated notches. Our initial efforts to employ DW-MTJs as an all-in-one stochastic synapse with both a mean and standard deviation didn’t end up meeting the quality metrics for hardware-friendly VI. In the future, new device stacks and methods for expressive anisotropy modification may make this idea still possible. However, as a fall back that immediately satisfies our requirements, we invented and detailed how the combination of a DW-MTJ synapse encoding the mean and a probabilistic Bayes-MTJ device, programmed via a ferroelectric or ionically modifiable layer, can robustly and expressively implement VI. This design includes a physics-informed small circuit model, that was scaled up to perform and demonstrate rigorous uncertainty quantification applications, up to and including small convolutional networks on a grayscale image classification task, and larger (Residual) networks implementing multi-channel image classification. Lastly, as these results and ideas all depend upon the idea of an inference application where weights (spintronic memory states) remain non-volatile, the retention of these synapses for the notched case was further interrogated. These investigations revealed and emphasized the importance of both notch geometry and anisotropy modification in order to further enhance the endurance of written spintronic states. In the near future, these results will be mapped to effective predictions for room temperature and elevated operation DW-MTJ memory retention, and experimentally verified when devices become available.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS2012

Causal Correlation Functions and Fourier Transforms: Application in Calculating Pressure Induced Shifts

By adopting a concept from signal processing, instead of starting from the correlation functions which are even, one considers the causal correlation functions whose Fourier transforms become complex. Their real and imaginary parts multiplied by 2 are the Fourier transforms of the original correlations and the subsequent Hilbert transforms, respectively. Thus, by taking this step one can complete the two previously needed transforms. However, to obviate performing the Cauchy principal integrations required in the Hilbert transforms is the greatest advantage. Meanwhile, because the causal correlations are well-bounded within the time domain and band limited in the frequency domain, one can replace their Fourier transforms by the discrete Fourier transforms and the latter can be carried out with the FFT algorithm. This replacement is justified by sampling theory because the Fourier transforms can be derived from the discrete Fourier transforms with the Nyquis rate without any distortions. We apply this method in calculating pressure induced shifts of H2O lines and obtain more reliable values. By comparing the calculated shifts with those in HITRAN 2008 and by screening both of them with the pair identity and the smooth variation rules, one can conclude many of shift values in HITRAN are not correct.

correlation function↗
DOE OSTI2022

Sensitivity Coefficients Calculated for the Prompt Neutron Decay Constant at or Near Delayed Critical [Abstract]

Experimenters at Los Alamos National Laboratory (LANL) measure the prompt neutron decay constant for many experiments at the National Criticality Experiments Research Center (NCERC) to infer reactivity and the effective neutron multiplication factor. These quantities are very important for nuclear criticality safety and validating nuclear data. Uncertainty in measures of criticality of an experimental configuration can be determined prior to physically performing the experiment by applying first order perturbation theory to Monte Carlo codes, such as MCNP®. The first order perturbation theory produces first derivatives of some nuclear parameter to nuclear data (e.g., cross section data). This first derivative is commonly referred to as a sensitivity coefficient. Currently, the MCNP® software has the capability of computing effective neutron multiplication factor sensitivity coefficients to cross section data. This work builds off of this MCNP® capability and the first order perturbation theory to provide a method of calculating sensitivity coefficients for the prompt neutron decay constant at or near delayed critical to cross section data. The prompt neutron decay constant sensitivity coefficient calculated in this work does not depend on any modification of the MCNP® source code. Prompt neutron decay constant sensitivity coefficient calculations can be used to infer reactivity and effective neutron multiplication factor sensitivity coefficient values as well. By investigating the trends of prompt neutron decay constant sensitivity coefficients for nuclide-reaction pairs across energy spectra, experiments can be designed to maximize or minimize the uncertainty in the prompt neutron decay constant in a particular energy region, which can lead to further optimization studies. Prompt neutron decay constant sensitivity coefficients will be calculated for the Jezebel benchmark. Subsequently, these sensitivity coefficients will be used in a data assimilation process to determine if there is or are optimal experiments that can be performed to provide insight into adjustments of uncertain/inaccurate cross section data. Specifically, the effect of the prompt neutron decay constant sensitivity coefficients on the nuclear data-induced uncertainty in the effective neutron multiplication factor will be examined

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 OSTI2019

Modeling and Computation of Casting Process by Particle Method

Since particle-based computational methods, such as smoothed particle hydrodynamics (SPH) method, are suitable to treat free surfaces and moving boundaries, it is expected that the simulation with particle methods can represent more accurate molten metal flow. Hence, we apply the SPH method to some casting processes and validate the numerical results by comparing with experiments. While fluid behaviors of the transportation and the mold filling processes are expressed sufficiently, ladle pouring process of aluminum alloy is not reproduced by the existing scheme. So we developed a new numerical model of surface oxide film, and our numerical results show good agreement with the experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Neutron Capture Cross Sections of Radioactive Nuclei

Alternative methods to calculate neutron capture cross sections on radioactive nuclei are reported using the theory of Inclusive Non-Elastic Breakup (INEB) developed by Hussein and McVoy. The statistical coupled-channels theory proposed in Ref. [2] is further extended in the realm of random matrices. The case of reactions with the projectile and the target being twocluster nuclei is also analyzed and applications are made for scattering from a deuteron target. An extension of the theory to a three-cluster projectile incident on a two-cluster target is also discussed. The theoretical developments described here should open new possibilities to obtain information on the neutron capture cross sections of radioactive nuclei using indirect methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2017

Kepler Planet Detection Metrics: Per-Target Detection Contours for Data Release 25

A necessary input to planet occurrence calculations is an accurate model for the pipeline completeness (Burke et al., 2015). This document describes the use of the Kepler planet occurrence rate products in order to calculate a per-target detection contour for the measured Data Release 25 (DR25) pipeline performance. A per-target detection contour measures for a given combination of orbital period, Porb, and planet radius, Rp, what fraction of transit signals are recoverable by the Kepler pipeline (Twicken et al., 2016; Jenkins et al., 2017). The steps for calculating a detection contour follow the procedure outlined in Burke et al. (2015), but have been updated to provide improved accuracy enabled by the substantially larger database of transit injection and recovery tests that were performed on the final version (i.e., SOC 9.3) of the Kepler pipeline (Christiansen, 2017; Burke Catanzarite, 2017a). In the following sections, we describe the main inputs to the per-target detection contour and provide a worked example of the python software released with this document (Kepler Planet Occurrence Rate Tools KeplerPORTs)1 that illustrates the generation of a detection contour in practice. As background material for this document and its nomenclature, we recommend the reader be familiar with the previous method of calculating a detection contour (Section 2 of Burke et al.,2015), input parameters relevant for describing the data quantity and quality of Kepler targets (Burke Catanzarite, 2017b), and the extensive new transit injection and recovery tests of the Kepler pipeline (Christiansen et al., 2016; Burke Catanzarite, 2017a; Christiansen, 2017).

Planet Detection Metrics↗
NASA NTRS2000

Effect of Negative Ions on the Conductivity of the Titan Atmosphere

In an earlier paper, Borucki et al (1987) calculated the electrical conductivity and electrical charge on aerosols in Titan's atmosphere due to the ionization by galactic cosmic rays and electron precipitation from Saturn's magnetosphere. The lower atmosphere was predicted to be substantially more conducting than the atmospheres of Earth and Venus because of the high concentration of free electrons. The prediction of a high conductivity is based on the lack of electrophillic species which form negative ions with low mobility and which reduce the number of free electrons. At that time, no molecular species capable of forming negative ions in concentrations sufficient to perturb the atmospheric conductivity were identified. Recently, E. Bakes and her colleagues have been investigating the formation of nitrogenous macromolecules using quantum mechanical methods. Their calculations indicate that the molecules will be highly electrophillic and are likely to be present in the atmosphere at mixing ratios of order 10(exp -7). This mixing ratio is sufficiently large that a substantial reduction in the conductivity is expected at altitudes below 100 km. Revision of the atmospheric model to accommodate the presence of negative ions and to increase the fidelity of the modeling will be described.

Borucki, W. J.↗
DOE OSTI2022

Formulation of the density eigenvalue problem in neutron transport for relevant engineering applications

A new formulation of the density eigenvalue problem for the neutron transport equation is presented. This new eigenvalue, named ζ eigenvalue can be introduced freely in the transport model, acting on a selected portion of the phase space. Despite its broader applications and its connection with the nature of the transport operator, the ζ eigenvalue has been presented here mainly as a design-oriented technique for the efficient evaluation of the critical concentration for a specific nuclide (or mixture of nuclides). This new eigenvalue is particularly adequate to study the definition of the material composition in the criticality design process of a multiplying system. The method is then applied for the study of classical problems such as the critical moderation ratio and the poison concentration to control the reactor. Some results are presented in one dimensional configuration using the multi- group spherical harmonics approach. This eigenvalue formulation proves to be a convenient and useful way to attain criticality, also for complex, realistic systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Extension of the PINN diffusion model to k-eigenvalue problems

This paper extends our recent work on the Physics-Informed Neural Networks (PINN) approach for the fixed source diffusion models and applies it to the diffusion theory based k-eigenvalue problems. To make the PINN equitable for the eigenvalue problems, we introduce a novel integral regularization term to the loss function in the framework, and allow the direct inference of the principal eigenvalue and the associated eigenfunction. The regularization term enforces a pre-defined value on the integration of the model predictions, and this value can be directly related to a physical property of the system. We also introduce an additional learnable parameter to approximate the principal eigenvalue. As a proof of principle, we solve the one-group two-dimensional k-eigenvalue neutron diffusion equation in this work. We then provide two numerical examples to demonstrate the applicability of the PINN approach. In each example, we solve the k-eigenvalue diffusion equation in a multi-region configuration constrained with a set of Robin boundary conditions for generality. We use a FEM solution based on the power-iteration method to verify the results of the PINN solution. The results showed relative percentage error in the predicted eigenvalue of about 0.77% and about 1.2% for example 1 and example 2, respectively. The mean absolute error in the predicted flux for example 1 is ∼ 0.002 and for example 2 is ∼ 0.0024. These results indicate some preliminary successes of the PINN application to k-eigenvalue problems. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1989

Estimation of velocity perturbations in satellite fragmentation events

The magnitude, variance and directionality of the velocity perturbations of the fragments of a satellite can shed valuable information regarding the nature and intensity of the fragmentation. Up until now, the only method used to calculate the three orthogonal components of the velocity change consisted of inverting the process of evaluating the changes in the orbital elements of the fragments due to velocity perturbing forces. But the traditional method failed in five different cases: (1) when the parent satellite's orbit was circular; (2) and (3) when the true anomaly of the parent was either 0 deg or 180 deg; and (4) and (5) when the argument of latitude of the parent was 90 deg or 270 deg. Described here is a new method of calculating the velocity perturbations which is free from the shortcomings of the traditional method and could be used in all occasions, provided the fragmentation data and the orbital elements data are consistent with one another.

Tan, Arjun↗
NASA NTRS1972

A theoretical and experimental investigation of the hypersonic turbulent boundary layer subject to normal and longitudinal pressure gradients and cross flow along a windward plane of symmetry.

A theoretical and experimental investigation for the hypersonic turbulent boundary layer undergoing both normal and longitudinal pressure gradients and cross flow along the plane of symmetry on the windward side of a body of revolution is presented. The purpose of this investigation is to provide a method for calculating boundary layers that occur on the centerline of symmetry of an inlet. The three-dimensional compressible integral equations are written along the symmetry plane and integrated numerically, and the effect of the cross flow on the behavior of the boundary layer is also studied. The experiments were performed at Mach number 6 and Reynolds number of 44,000,000/ft for the axisymmetric case and for the case with cross flow. The velocity profiles and the boundary layer parameters agree relatively well with theoretical values with and without cross flow effects. The presence of the cross flow alters significantly the value of the profiles. Total temperature-velocity correlations deviate from the linear Crocco distribution. Heat transfer rates agree with theoretical values and values obtained by the flat plate reference enthalpy method. Cross flow influence on the heat transfer rates is of minor relevance.

Zakkay, V.↗
DOE OSTI2025

Machine-learning-informed scattering correlation analysis of sheared colloids

We have carried out theoretical analysis, Monte Carlo simulations and machine-learning analysis to quantify microscopic rearrangements of dilute dispersions of spherical colloidal particles from coherent scattering intensity. Both monodisperse and polydisperse dispersions of colloids were created and underwent a rearrangement consisting of an affine simple shear and non-affine rearrangement using the Monte Carlo method. We calculated the coherent scattering intensity of the dispersions and the correlation function of intensity before and after the rearrangement and generated a large data set of angular correlation functions for varying system parameters, including number density, polydispersity, shear strain and non-affine rearrangement. Singular value decomposition of the data set shows the feasibility of machine-learning inversion from the correlation function for the polydispersity, shear strain and non-affine rearrangement using only three parameters. A Gaussian process regressor is then trained on the data set and can retrieve the affine shear strain, non-affine rearrangement and polydispersity with relative errors of 3%, 1% and 6%, respectively. Altogether, our model provides a framework for quantitative studies of both steady and non-steady microscopic dynamics of colloidal dispersions using coherent scattering methods.

Gaussian process regression↗
DOE OSTI2021

Complexity and near extremal charged black branes

Highlights: • Trying to explore a remedy to overcome the puzzle in holographic complexity by assuming a cutoff behind the horizon and the corresponding boundary term. • Holographic complexity growth obtained by the WDW patch on shell action is insensitive to the UV cutoff at the late time. Therefore one may conclude that either on shell action does not compute complexity, or the late time is not given by the physical charges, or in our computations we have missed something! • We find that one needs to carefully compute on shell action taking into account the behind the horizon and the corresponding boundary term. Doing so, one arrives at a consistent picture. • We check this consistency by evaluating the late time complexity growth for near horizon limit of the near extremal black hole using two different methods. We would have obtained non-consistent picture, if we had not considered the behind the horizon cutoff as well as the corresponding boundary term. We compute holographic complexity of charged black brane solutions in arbitrary dimensions for the near horizon limit of near extremal case using two different methods. The corresponding complexity may be obtained either by taking the limit from the complexity of the charged black brane, or by computing the complexity for near horizon limit of near extremal solution. One observes that these results coincide if one assumes to have a cutoff behind horizon whose value is fixed by UV cutoff and also taking into account a proper counterterm evaluated on this cutoff. We also consider the situation for Vaidya charged black branes too.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Application of chord length distributions and principal component analysis for quantification and representation of diverse polycrystalline microstructures

Highlights: • Microstructure quantification based on directional chord length distribution and principal component analysis is presented • A boundary-based algorithm of calculations of chord length distributions is introduced • Utility of the quantification framework is shown in three case studies on synthetic, experimental, and simulation datasets • Angularly resolved chord length distribution captures essential details of grain size and shape in diverse microstructures • Low-dimensional representations of directional chord length distributions are obtained using principal component analysis - Abstract: Quantification of mesoscale microstructures of polycrystalline materials is important for a range of practical tasks of materials design and development. The current protocols of quantifying grain size and morphology often rely on microstructure metrics (e.g., mean grain diameter) that overlook important details of the mesostructure. In this work, we present a quantification framework based on directionally resolved chord length distribution and principal component analysis as a means of extracting additional information from 2-D microstructural maps. Towards this end, we first present in detail a method for calculating chord length distribution based on boundary segments available in modern digital datasets (e.g., from microscopy post-processing) and their low-rank representations by principal component analysis. The utility of the proposed framework for capturing grain size, morphology, and their anisotropy for efficient visualization, representation, and specification of polycrystalline microstructures is then demonstrated in case studies on datasets from synthetic generation, experiments (on Ni-base superalloys), and simulations (on steel during recrystallization).

36 MATERIALS SCIENCE↗
DOE OSTI2025

Understanding Pore Filling Processes and Adsorption/Desorption Hysteresis in Nanoporous Metal–Organic Frameworks: Insights from Grand Canonical Monte Carlo Simulations and Free Energy Calculations

Grand canonical Monte Carlo (GCMC) simulations were used to investigate pore filling and hysteresis in nanoporous metal-organic frameworks (MOFs). Adsorption and desorption isotherms were calculated for argon at 87 K in 1866 MOFs from the CoRE MOF database and for short n-alkanes in selected MOFs, keeping the adsorbent structure rigid. Analysis of the molecular configurations showed two different mechanisms and origins of hysteresis: one involving a transition of the adsorbate arrangement in the pores similar to a gas-to-liquid transition associated with a large change in the loading and one more similar to a liquid-to-solid transition associated with a relatively small change in the loading. Our GCMC simulations in MOFs with diverse pore topologies indicate exceptions to an empirical relationship for the minimum diameter of a cylindical pore required for hysteresis as a function of the adsorbate diameter and reduced temperature. The simulations reveal some structures where isotherms exhibit two steps in the adsorption branch and only one step in the desorption branch. Hysteresis loops with a different number of adsorption and desorption steps are not common. Here, to better understand why hysteresis is observed in the GCMC simulations, the concept of the transition probability for observing a step in the adsorption isotherm at a given pressure in a GCMC simulation is introduced. We used two different methods to calculate the transition probabilities and find that these yield comparable results. Furthermore, the transition probability provides a measure for the length of GCMC simulations to yield reliable results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗