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,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Weak gravitational lensing by Kerr-MOG black hole and Gauss–Bonnet theorem

Highlights: • The deflection angle of the Kerr-MOG black hole is significantly greater than the original Kerr black hole. • In the future astrophysical observation might shed light on the effects of modified gravity theory. • MOG effect could be taken into account in the gravitational lensing experiments. • Our results put forward that many features of weak gravitational lensing arise from the topological properties of an underlying optical manifold. • For weak gravitational lensing of (at least) spherically symmetric distributions, the Gauss–Bonnet approach holds several advantages over common methods. - Abstract: The deflection angle of Kerr-MOG black holes is studied for different values of the parameter in modified gravity (MOG). To this end, we employ the Gauss–Bonnet theorem, which was first studied by Gibbons and Werner and then extended by Ono, Ishihara and Asada, who use a generalized optical metric where the deflection of light for an observer and source at finite distance. By using this method, we study the weak gravitational lensing by Kerr-MOG black hole. Our computations show that with an increase in the MOG parameter (α), the deflection angle becomes significantly larger than that of Kerr black hole. The results obtained show that MOG effect could be taken into account in the gravitational lensing experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Improved Method for Evaluating and Specifying the Chromaticity of Light Sources

This article recommends a method for calculating and specifying light source chromaticity using the CIE 2015 10° color matching functions (CMFs), which, according to analysis of existing psychophysical experiment data, can reduce visual mismatch compared to specifications based on the traditional CIE 1931 2° CMFs. Specifically, this work evaluates, documents, and recommends for adoption a supporting system of measures to be used with the CIE 2015 10° CMFs: a new uniform chromaticity scale (UCS) diagram with coordinates (s, t), a measure of correlated color temperature (CCTst), and a measure of distance from the Planckian locus (Dst). It also presents options for updating nominal classification quadrangles. A complete method of this nature has not yet been standardized, which may be contributing to the slow uptake of the CIE 2015 CMFs. The proposed tools are analogous to u, v, CCT, Duv, and the ANSI C78.377 nominal classification quadrangles that are currently specified in the CIE 1960 UCS using the CIE 1931 2° CMFs. While conceptually equivalent, the differences between the current standard method and the new st system are important for reducing unintended mismatch in the chromaticity of light. Furthermore, the implications of changing chromaticity specification methods are identified by a comparison over a diverse set of real light source SPDs. Finally, as a step toward improving the design and specification of light sources for illumination purposes, we propose a process for implementing the recommendations.

42 ENGINEERING↗
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 OSTI2017

Universality of scanning tunneling microscopy in cuprate superconductors

We consider the problem of local tunneling into cuprate superconductors, combining model-based calculations for the superconducting order parameter with wave function information obtained from first-principles electronic structure. For some time it has been proposed that scanning tunneling microscopy (STM) spectra do not reflect the properties of the superconducting layer in the CuO 2 plane directly beneath the STM tip, but rather a weighted sum of spatially proximate states determined by the details of the tunneling process. These ``filter'' ideas have been countered with the argument that similar conductance patterns have been seen around impurities and charge ordered states in systems with atomically quite different barrier layers. Here we use a recently developed Wannier function-based method to calculate topographies, spectra, conductance maps, and normalized conductance maps close to impurities. We find that it is the local planar Cu d x 2 –y 2 Wannier function, qualitatively similar for many systems, that controls the form of the tunneling spectrum and the spatial patterns near perturbations. Furthermore we explain how, despite the fact that STM observables depend on the materials-specific details of the tunneling process and setup parameters, there is an overall universality in the qualitative features of conductance spectra. In particular, we discuss why STM results on Bi 2 Sr 2 CaCu 2 O 8 (BSCCO) and Ca 2–x Na x CuO 2 Cl 2 (NaCCOC) are essentially identical.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2017

Getting super-excited with modified dispersion relations

We demonstrate that in some regions of parameter space, modified dispersion relations can lead to highly populated excited states, which we dub as 'super-excited' states. In order to prepare such super-excited states, we invoke dispersion relations that have negative slope in an interim sub-horizon phase at high momenta. This behaviour of quantum fluctuations can lead to large corrections relative to the Bunch-Davies power spectrum, which mimics highly excited initial conditions. We identify the Bogolyubov coefficients that can yield these power spectra. In the course of this computation, we also point out the shortcomings of the gluing method for evaluating the power spectrum and the Bogolyubov coefficients. As we discuss, there are other regions of parameter space, where the power spectrum does not get modified. Therefore, modified dispersion relations can also lead to so-called 'calm excited states'. We conclude by commenting on the possibility of obtaining these modified dispersion relations within the Effective Field Theory of Inflation.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2018

Supersymmetric symplectic quantum mechanics

Symplectic Quantum Mechanics SQM considers a non-commutative algebra of functions on a phase space Γ and an associated Hilbert space H{sub Γ} to construct a unitary representation for the Galilei group. From this unitary representation the Schrödinger equation is rewritten in phase space variables and the Wigner function can be derived without the use of the Liouville–von Neumann equation. In this article we extend the methods of supersymmetric quantum mechanics SUSYQM to SQM. With the purpose of applications in quantum systems, the factorization method of the quantum mechanical formalism is then set within supersymmetric SQM. A hierarchy of simpler hamiltonians is generated leading to new computation tools for solving the eigenvalue problem in SQM. We illustrate the results by computing the states and spectra of the problem of a charged particle in a homogeneous magnetic field as well as the corresponding Wigner function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗