Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational 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 1,243 records · Page 69

Ab initio ground states of strongly-correlated materials on quantum computers

The accurate first-principles description of strongly-correlated materials is an important and challenging problem in condensed matter physics. Ab initio downfolding has emerged as a way of deriving accurate many-body Hamiltonians including strong correlations, representing a subspace of interest of a material, using density functional theory calculations as a starting point. However, the solution of these material-specific models can scale exponentially on classical computers, constituting a challenge. Here we propose that utilizing quantum computers for obtaining the properties of downfolded Hamiltonians yields an accurate description of the ground state properties of strongly-correlated systems, while circumventing the exponential scaling problem. We benchmark the solution of Hubbard-like models obtained through downfolding by utilizing a classical tensor network implementation of variational quantum eigensolvers (VQE), and we reveal a strategy for driving the optimization through a hybrid minimization of the energy and maximization of the overlap with an approximate solution obtained through low-cost computational methods. This results in a reduction of the energy error by orders of magnitude compared to conventional VQE approaches, and allows us to reproduce long-range correlations for the first time. We demonstrate our first-principles approach for diverse strongly-correlated materials, correctly predicting the antiferromagnetic state of one-dimensional cuprate Ca 2 CuO 3 , the excitonic ground state of monolayer WTe2, and the charge-ordered state of correlated metal SrVO 3 . Our efficient computational implementation allows us to simulate large systems with up to 54 qubits and encompassing up to four correlated bands, which is indicative of the complexity that our framework can address.

Antonios M Alvertis↗

Simulated transonic flows for aircraft with nacelles, pylons, and winglets

A computational method which simulates transonic flow about wing-fuselage configurations has been extended to include the treatment of multiple body and non-planar wing surfaces. The finite difference relaxation scheme is characterized by a modified small disturbance flow equation and multiple embedded grid system. Wing-body combinations with as many as four nacelles/pods, four pylons, and wing-tip-mounted winglets can be analyzed. A scheme for modeling inlet spillage and engine exhaust interference effects has been included. Computed results are correlated with experimental data for three transport configurations.

Boppe, C. W.↗

Computational Analysis in Support of the SSTO Flowpath Test

A synergistic approach of combining computational methods and experimental measurements is used in the analysis of a hypersonic inlet. There are four major focal points within this study which examine the boundary layer growth on a compression ramp upstream of the cowl lip of a scramjet inlet. Initially, the boundary layer growth on the NASP Concept Demonstrator Engine (CDE) is examined. The follow-up study determines the optimum diverter height required by the SSTO Flowpath test to best duplicate the CDE results. These flow field computations are then compared to the experimental measurements and the mass average Mach number is determined for this inlet.

Duncan, Beverly S.↗

Accelerated basis adaptation in homogeneous chaos spaces

Polynomial chaos expansions (PCE) provide an efficient approach to uncertainty quantification (UQ) and have been adapted to diverse applications across the spectrum of science and engineering. For situations involving large stochastic parameterizations, the curse of dimensionality renders PCE-based methods computationally prohibitive. A basis adaptation approach for PCE was proposed by Tipireddy and Ghanem (2014) which transforms the input random variables through an isometry such that sufficient probabilistic characterization of specific quantities of interest (QoI) is concentrated in an algebraic manifold embedded in the linear span of the dominant transformed variables. While quite versatile, that original version of the PCE basis adaptation exhibited slow convergence for a number of problems of practical significance. In the present paper, we propose two novel methods to accelerate the convergence of the original basis adaptation approach, thus expanding its reach while also providing insight into its performance. In the first method, information gained from a pilot PCE representation is used to correct the mean and Gaussian coefficients in the adapted space. By taking advantages of probabilistic information in higher dimensional adaptation gleaned from an initial adaptation, the second method updates the rotation matrix used to identify the dominant transformed variables. In this manner, the new rotation matrix concentrates even more probabilistic information in its first few dimensions. These two method can be combined to achieve even better performance, the combined method is referred to as sequentially optimized adaptation method. Here, the methods are demonstrated on an analytical test function and a model of a space structure with several sub-components and a non-smooth quantity of interest representing the maximum acceleration over time. Both methods achieve accelerated convergence of the basis adaptation approach with negligible additional costs.

Accelerated convergence↗

Nuclear–Electronic Orbital General Rate Theory: Predicting Hydrogen Kinetic Isotope Effects in the Deep Tunneling Regime

Hydrogen transfer is a critical component of many chemical and biological processes. The ratio of rate constants for hydrogen and deuterium transfer defines the H/D kinetic isotope effect (KIE), which is a powerful tool for elucidating hydrogen transfer mechanisms. Interpretation of experimental H/D KIEs relies on accurate and affordable computational methods. However, due to their light mass, hydrogen and deuterium can undergo tunneling, which is challenging to describe in multidimensional molecular systems. Herein, we introduce the nuclear–electronic orbital general rate theory (NEO-GRT), which enables the efficient prediction of H/D KIEs based on full-dimensional molecular quantum chemistry calculations. The NEO-GRT approach describes the hydrogen transfer rate constant with a general expression that spans the vibrationally adiabatic and nonadiabatic hydrogen tunneling regimes. The input quantities are computed using NEO density functional theory, which treats the transferring hydrogen or deuterium nucleus quantum mechanically on the same level as the electrons. We investigate two intramolecular proton transfer reactions in organic molecules at temperatures down to 50 K to evaluate the performance of NEO-GRT by comparison to transition state theory and ring-polymer instanton theory. The KIEs computed with NEO-GRT agree with those calculated using ring-polymer instanton theory for the full-dimensional molecular systems at the same level of electronic structure theory. This agreement indicates that NEO-GRT captures the deep hydrogen tunneling effects, in contrast to transition state theory, which neglects such effects. Given its relatively low computational cost, NEO-GRT is a promising approach for predicting H/D KIEs in large organic and organometallic systems.

Hydrogen↗

The effect of subpixel clouds on remote sensing

A method for estimating the cloud effect on remote sensing is described, and it is applied to cloudiness in several climatological conditions. The algorithm is based on the Haurwitz (1948) measurements of the cloud layer transmission of solar radiation for an overcast sky and on an empirical interpolation of data for broken cloudiness by Pochop et al. (1968). Radiances for a sunny area observed directly from space and through a cloud, and for a shady area observed from space and through a cloud are computed. Methods for detecting the cloud effect from satellite images are discussed. The relation between cloud reflectance and cloud size is studied. It is observed that the subpixel clouds affect the detected radiance and vegetation index, and the effect depends on the cloud types and the dependence of the cloud transmissivity on cloud fraction. Procedures for decreasing or eliminating cloud effect are examined.

Kaufman, Yoram J.↗

An M-step preconditioned conjugate gradient method for parallel computation

This paper describes a preconditioned conjugate gradient method that can be effectively implemented on both vector machines and parallel arrays to solve sparse symmetric and positive definite systems of linear equations. The implementation on the CYBER 203/205 and on the Finite Element Machine is discussed and results obtained using the method on these machines are given.

Adams, L.↗

Effects of Mistuning on the Forced Response of Turbomachinery Rotors

The basic purpose of this research has been to develop a fundamental understanding of the effects of blade-to-blade dissimilarities, or mistuning, on the dynamics of nearly cyclic bladed-disk assemblies. This topic is of importance as mistuning has been shown to increase the forced response amplitudes of some blades significantly, and even to lead to blade failure. Furthermore, the current trend toward high performance propulsion turbomachinery designed for finite service life demands an accurate prediction of system performance and dynamics at the design stage. This objective has been achieved by carrying out the following tasks. First, the investigation of the free and forced responses of representative, yet sufficiently simple blade assembly models that capture all the important characteristics of typical turbomachinery rotors. Second, the development of computational methods that predict the effects of mistuning in a systematic and reliable way, along with the development of a systematic reduced-order modeling procedure for mistuned bladed disks. Third, the application of these findings and the tools developed to an industrial rotor, namely the first stage of turbine blades of the oxidizer turbopump in the space shuttle main rocket engine (SSME). The research supported by NASA has led to the development of a coherent theory for mistuned blade assemblies. The further implementation of these computational tools into the forced response prediction system currently under development in the Structural Dynamics Branch-FREPS-would enable the designer and the analyst to: (1) identify types of blade assemblies highly sensitive to mistuning in various frequency and other parameter ranges and (2) characterize mistuning effects on their forced response by predicting true response amplitudes and fatigue life estimates.

Pierre, Christophe↗

Effects of Mistuning on the Forced Response of Turbomachinery Rotors

The basic purpose of this research has been to develop a fundamental understanding of the effects of blade-to-blade dissimilarities, or mistuning, on the dynamics of nearly cyclic bladed-disk assemblies. This topic is of importance as mistuning has been shown to increase the forced response amplitudes of some blades significantly, and even to lead to blade failure. Furthermore, the current trend toward high performance propulsion turbomachinery designed for finite service life demands an accurate prediction of system performance and dynamics at the design stage. This objective has been achieved by carrying out the following tasks. First, the investigation of the free and forced responses of representative, yet sufficiently simple blade assembly models that capture all the important characteristics of typical turbomachinery rotors. Second, the development of computational methods that predict the effects of mistuning in a systematic and reliable way, along with the development, of a systematic reduced-order modeling procedure for mistuned bladed disks. Third, the application of these findings and the tools developed to an industrial rotor, namely the first stage of turbine blades of the oxidizer turbopump in the space shuttle main rocket engine (SSME). The research supported by NASA has led to the development of a coherent theory for mistuned blade assemblies. The further implementation of these computational tools into the forced response prediction system currently under development in the Structural Dynamics Branch Forced Response Prediction System (FREPS) would enable the designer and the analyst (1) to identify types of blade assemblies highly sensitive to mistuning in various frequency and other parameter ranges and (2), to characterize mistuning effects on their forced response by predicting true response amplitudes and fatigue life estimates.

Stefko, George L.↗

Angular-momentum projection in coupled-cluster theory: Structure of 34 Mg

Single- reference coupled-cluster theory is an accurate and affordable computational method for the nuclear many-body problem. For open-shell nuclei, the reference state typically breaks rotational invariance and angular momentum must be restored as a good quantum number. We perform angular-momentum projection after variation and employ the disentangled coupled-cluster formalism and a Hermitian approach. We compare our results with benchmarks for 8 Be and 20 Ne using a two-nucleon interaction from chiral effective field theory and for pf-shell nuclei within the traditional shell model. We compute the rotational band in the exotic nucleus 34 Mg and find agreement with data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Synchrotron-based techniques for characterizing STCH water-splitting materials

Understanding the role of oxygen vacancy–induced atomic and electronic structural changes to complex metal oxides during water-splitting processes is paramount to advancing the field of solar thermochemical hydrogen production (STCH). The formulation and confirmation of a mechanism for these types of chemical reactions necessitate a multifaceted experimental approach, featuring advanced structural characterization methods. Synchrotron X-ray techniques are essential to the rapidly advancing field of STCH in part due to properties such as high brilliance, high coherence, and variable energy that provide sensitivity, resolution, and rapid data acquisition times required for the characterization of complex metal oxides during water-splitting cycles. X-ray diffraction (XRD) is commonly used for determining the structures and phase purity of new materials synthesized by solid-state techniques and monitoring the structural integrity of oxides during water-splitting processes (e.g., oxygen vacancy–induced lattice expansion). X-ray absorption spectroscopy (XAS) is an element-specific technique and is sensitive to local atomic and electronic changes encountered around metal coordination centers during redox. While in operando measurements are desirable, the experimental conditions required for such measurements (high temperatures, controlled oxygen partial pressures, and H 2 O) practically necessitate in situ measurements that do not meet all operating conditions or ex situ measurements. Here, we highlight the application of synchrotron X-ray scattering and spectroscopic techniques using both in situ and ex situ measurements, emphasizing the advantages and limitations of each method as they relate to water-splitting processes. The best practices are discussed for preparing quenched states of reduction and performing synchrotron measurements, which focus on XRD and XAS at soft (e.g., oxygen K-edge, transition metal L-edges, and lanthanide M-edges) and hard (e.g., transition metal K-edges and lanthanide L-edges) X-ray energies. The X-ray absorption spectra of these complex oxides are a convolution of multiple contributions with accurate interpretation being contingent on computational methods. The state-of-the-art methods are discussed that enable peak positions and intensities to be related to material electronic and structural properties. Through careful experimental design, these studies can elucidate complex structure–property relationships as they pertain to nonstoichiometric water splitting. A survey of modern approaches for the evaluation of water-splitting materials at synchrotron sources under various experimental conditions is provided, and available software for data analysis is discussed.

08 HYDROGEN↗

Combining protein sequences and structures with transformers and equivariant graph neural networks to predict protein function

Abstract Motivation Millions of protein sequences have been generated by numerous genome and transcriptome sequencing projects. However, experimentally determining the function of the proteins is still a time consuming, low-throughput, and expensive process, leading to a large protein sequence-function gap. Therefore, it is important to develop computational methods to accurately predict protein function to fill the gap. Even though many methods have been developed to use protein sequences as input to predict function, much fewer methods leverage protein structures in protein function prediction because there was lack of accurate protein structures for most proteins until recently. Results We developed TransFun—a method using a transformer-based protein language model and 3D-equivariant graph neural networks to distill information from both protein sequences and structures to predict protein function. It extracts feature embeddings from protein sequences using a pre-trained protein language model (ESM) via transfer learning and combines them with 3D structures of proteins predicted by AlphaFold2 through equivariant graph neural networks. Benchmarked on the CAFA3 test dataset and a new test dataset, TransFun outperforms several state-of-the-art methods, indicating that the language model and 3D-equivariant graph neural networks are effective methods to leverage protein sequences and structures to improve protein function prediction. Combining TransFun predictions and sequence similarity-based predictions can further increase prediction accuracy. Availability and implementation The source code of TransFun is available at https://github.com/jianlin-cheng/TransFun.

59 BASIC BIOLOGICAL SCIENCES↗

Overview of DOE-NE NEAMS Program

The NEAMS (Nuclear Energy Advanced Modeling and Simulation) program is a multinational lab team effort aiming to develop and deploy predictive computer methods for the analysis and design of advanced nuclear reactors. US Department of Energy has made significant investments in advanced mod-sim during last decade between CASL and NEAMS. Next year, aim is for all NE mod sim to be centralized in a single program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Review on modeling of the anode solid electrolyte interphase (SEI) for lithium-ion batteries

A passivation layer called the solid electrolyte interphase (SEI) is formed on electrode surfaces from decomposition products of electrolytes. The SEI allows Li + transport and blocks electrons in order to prevent further electrolyte decomposition and ensure continued electrochemical reactions. The formation and growth mechanism of the nanometer thick SEI films are yet to be completely understood owing to their complex structure and lack of reliable in situ experimental techniques. Significant advances in computational methods have made it possible to predictively model the fundamentals of SEI. This review aims to give an overview of state-of-the-art modeling progress in the investigation of SEI films on the anodes, ranging from electronic structure calculations to mesoscale modeling, covering the thermodynamics and kinetics of electrolyte reduction reactions, SEI formation, modification through electrolyte design, correlation of SEI properties with battery performance, and the artificial SEI design. Multi-scale simulations have been summarized and compared with each other as well as with experiments. Computational details of the fundamental properties of SEI, such as electron tunneling, Li-ion transport, chemical/mechanical stability of the bulk SEI and electrode/(SEI/) electrolyte interfaces have been discussed. This review shows the potential of computational approaches in the deconvolution of SEI properties and design of artificial SEI. We believe that computational modeling can be integrated with experiments to complement each other and lead to a better understanding of the complex SEI for the development of a highly efficient battery in the future.

25 ENERGY STORAGE↗

Analysis of axisymmetric body effects on rotor aerodynamics using modified slender body theory

A modified slender body theory was used in modelling the flowfield around an axisymmetric body at zero angle of attack. The model was defined in an attempt to characterize body-induced effects on helicopter rotor performance and loading using efficient computation methods. The streamfunction for an axial source distribution was applied in the body model, which was then adapted to the helicopter analysis code (CAMRAD). The body-induced velocities were calculated for separate locations on the rotor disk. Three different rotor models in nonuniform inflow with a free wake geometry were considered. Comparison of predictions with data on ellipsoidal and axisymmetric bodies showed good agreement for moderately thick ellipsoidal bodies. Oscillatory blade bending moments were increased by the presence of the body, but overall performance impact was negligible. Computational efficiency was 10-20 percent better than with a potential flow panel code.

Yamauchi, G. K.↗

Comparison of Methods for Determining Boundary Layer Edge Conditions for Transition Correlations

Data previously obtained for the X-33 in the NASA Langley Research Center 20-Inch Mach 6 Air Tunnel have been reanalyzed to compare methods for determining boundary layer edge conditions for use in transition correlations. The experimental results were previously obtained utilizing the phosphor thermography technique to monitor the status of the boundary layer downstream of discrete roughness elements via global heat transfer images of the X-33 windward surface. A boundary layer transition correlation was previously developed for this data set using boundary layer edge conditions calculated using an inviscid/integral boundary layer approach. An algorithm was written in the present study to extract boundary layer edge quantities from higher fidelity viscous computational fluid dynamic solutions to develop transition correlations that account for viscous effects on vehicles of arbitrary complexity. The boundary layer transition correlation developed for the X-33 from the viscous solutions are compared to the previous boundary layer transition correlations. It is shown that the boundary layer edge conditions calculated using an inviscid/integral boundary layer approach are significantly different than those extracted from viscous computational fluid dynamic solutions. The present results demonstrate the differences obtained in correlating transition data using different computational methods.

Liechty, Derek S.↗

Composite sandwich thermostructural behavior - Computational simulation

Computational methods have been developed for simulating the thermomechanical behavior of composite sandwiches, in which the analyses with several levels of progressive sophistication were used in conjunction with composite hygrothermomechanical theory. The sophistication levels include: (1) three-dimensional detailed finite element modeling of the honeycomb, the adhesive, and the composite faces; (2) three-dimensional finite element modeling assuming a homogeneous core; (3) laminate theory simulation; and (4) simple equations for predicting the equivalent properties of the honeycomb core. These levels have been packaged into a procedure embedded in a computer code streamlined for the simulation of the composite sandwich hygrothermal and structural behavior. It is shown that in order to properly simulate the thermomechanical response of the composite sandwich, all the honeycomb thermal and mechanical properties must be used.

Chamis, C. C.↗