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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Mars Sample Return, Sample Retrieval Lander, Reaction Control System Jet Interaction Supersonic Wind Tunnel Test Overview with CFD Predictions

NASA's Mars Sample Return campaign will be launching several missions over the next decade that will work together to return rock samples from Mars back to Earth. The Sample Retrieval Lander (SRL) will deliver the Mars Ascent Vehicle and fetch rover to the surface of Mars in 2006. Rock samples collected by the Mars 2020 Perseverance rover, landing in early 2021, will be loaded on the the ascent vehicle to be launched into Mars orbit for retrieval by yet another spacecraft. The Sample Retrieval Lander will be a blunt entry capsule similar to past Mars entry vehicles like Mars Science Laboratory and Viking. The vehicle will fly a guided entry, using a small lift vector produced by a non-zero trim angle of attack to eliminate downrange and crossrange position errors at the point of parachute deploy. This energy and heading management is achieved with a reaction control system (RCS) that directs the bank angle of the vehicle and also minimizes unwanted capsule dynamics. The reaction control system and control design is based on the Mars Science Laboratory and Mars 2020 RCS systems. However, due to packaging constraints, the backshell of this new entry vehicle has a different geometry than those earlier designs. To certify the RCS system for flight the project must characterize the jet plume interactions with the capsule backshell that could impair or significantly augment the RCS control authority. This characterization will be done through a combination of computational fluid dynamics (CFD) analysis and wind tunnel test. Two candidate arrangements of the RCS jets have been identified for the SRL vehicle and are currently under evaluation before final selection. The aero/RCS plume interactions of these candidate configurations have been measured in a supersonic wind tunnel test in NASA Langley's Unitary Plan Wind Tunnel. The test was conducted in the fall of 2020 and data is currently being reduced. An overview of the candidate RCS configurations are presented here with an overview of the wind tunnel model design, jet scaling and scaled nozzle design, and the test matrix. Preliminary CFD runs are presented with an assessment of the predicted plumes and their interaction with the wake flow of the vehicle. The predicted effects of the model sting is provided as well. This high fidelity wind tunnel test is being conducted much earlier in the SRL project than would normally be done. The test was funded as part of a CFD evaluation task funded by NASA's Aerosciences Evaluation and Test Capabilities Project. The objective of the evaluation task was to compare the ability of CFD to predict complex flows with data that can be measured in the Langley Unitary Plan Wind Tunnel. RCS Jet interactions were selected as a type of complex flow that is important to NASA missions. In addition to providing useful data to the SRL project, there was added emphasis on quantifying the accuracy of the CFD predictions and wind tunnel test data. An overview of the uncertainty quantification methodologies for computational and experimental portions of this test is presented.

blunt body↗
DOE OSTI2022

Quantitative interpretation of time-lapse seismic data at Farnsworth field unit: Rock physics modeling, and calibration of simulated time-lapse velocity responses

Here, this study investigates the contribution of fluid saturation variation to the time-lapse velocity response by performing fluid substitution modeling. The methodology is exemplified by the time-lapse seismic monitoring of carbon dioxide at Farnsworth field unit (FWU). In order to evaluate the fluid distribution in a matured oil reservoir, the Southwest Regional Partnership (SWP) acquired multiple vertical seismic profile (VSP) surveys at different times during the CO 2 –water alternatinggas (WAG) injection period. In this work, we present a thorough methodology for computing the elastic response of the saturated rock for different fluid saturations using a site-specific petro-elastic model (PEM). The output from the PEM was combined with results from a fluid compositional model to compute the seismic velocities at times corresponding to each VSP survey. To produce a calibrated simulated response, the measured time-lapse seismic velocities were integrated into the numerical simulation model. The mismatches between the predicted and measured time-lapse velocities were minimized through an iterative calibration process using a trained artificial neural network proxy (ANN) coupled with a particle swarm optimizer (PSO). Our study indicates that the hybrid optimization workflow can effectively perform the history matching. With an accurate prediction of the hydrodynamic properties, the migration of CO 2 within the subsurface was modeled by predicting the spatial velocity distribution for a radius of 305 m around the injection well. The technology demonstrated and the expertise gained from this study can guide similar CO 2 -WAG projects.

58 GEOSCIENCES↗
DOE OSTI2014

Progressive construction of a parametric reduced‐order model for PDE‐constrained optimization

Summary An adaptive approach to using reduced‐order models (ROMs) as surrogates in partial differential equations (PDE)‐constrained optimization is introduced that breaks the traditional offline‐online framework of model order reduction. A sequence of optimization problems constrained by a given ROM is defined with the goal of converging to the solution of a given PDE‐constrained optimization problem. For each reduced optimization problem, the constraining ROM is trained from sampling the high‐dimensional model (HDM) at the solution of some of the previous problems in the sequence. The reduced optimization problems are equipped with a nonlinear trust‐region based on a residual error indicator to keep the optimization trajectory in a region of the parameter space where the ROM is accurate. A technique for incorporating sensitivities into a reduced‐order basis is also presented, along with a methodology for computing sensitivities of the ROM that minimizes the distance to the corresponding HDM sensitivity, in a suitable norm. The proposed reduced optimization framework is applied to subsonic aerodynamic shape optimization and shown to reduce the number of queries to the HDM by a factor of 4‐5, compared with the optimization problem solved using only the HDM, with errors in the optimal solution far less than 0.1%. Copyright © 2014 John Wiley & Sons, Ltd.

Zahr, Matthew J.↗
DOE OSTI2014

Systems Life Cycle and Its Relation with the Triple Helix

This chapter examines the life cycle of complex systems in light of the dynamic interconnections among the university, industry and government sectors. Each sector is motivated in its resource allocation by principles discussed elsewhere in this book and yet remains complementary es-tablishing enduring and fundamental relationships. Industry and Government depend upon an educated workforce; universities depend upon industry to spark the R&D which is needed and to sponsor some basic research and much applied research. Government depends upon industry to address operational needs and provide finished products while universities offer government (along with industry) problem solving and problem solving environments. The life cycle of complex systems in this chapter will be examined in this context, providing historical examples. Current examples will then be examined within this multi-dimensional context with respect to the phases of program and project life cycle management from requirements definition through retirement and closeout of systems. During the explanation of these examples, the advances in research techniques to collect, analyze, and process the data will be examined.

and Law↗
DOE OSTI2015

Solving a Large Scale Thermal Radiation Problem Using an Interoperable Executive Library Framework on Petascale Supercomputers

We present a novel methodology to compute the transient thermal condi- tion of a set of objects in an open space environment. The governing energy equa- tion and the convective energy transfer are solved by the sparse iterative solvers. The average radiating energy on a set of surfaces is represented by a linear system of the radiosity equations, which is factorized by an out-of-core parallel Cholesky decomposition solver. The coupling and interplay of the direct radiosity solver us- ing GPUs and the CPU-based sparse solver are handled by a light weight software integrator called Interoperable Executive Library (IEL). IEL manages the distribu- tion of data and memory, coordinates communication among parallel processes, and also directs execution of the set of loosely coupled physics tasks as warranted by the thermal condition of the simulated object and its surrounding environment.

Wong, Kwai↗
DOE OSTI2016

Energy dissipation via acoustic emission in ductile crack initiation

In this article, we present a modeling approach to estimate the energy release due to ductile crack initiation in conjunction to the energy dissipation associated with the formation and propagation of transient stress waves typically referred to as acoustic emission. To achieve this goal, a ductile fracture problem is investigated computationally using the finite element method based on a compact tension geometry under Mode I loading conditions. To quantify the energy dissipation associated with acoustic emission, a crack increment is produced given a pre-determined notch size in a 3D cohesive-based extended finite element model. The computational modeling methodology consists of defining a damage initiation state from static simulations and linking such state to a dynamic formulation used to evaluate wave propagation and related energy redistribution effects. The model relies on a custom traction separation law constructed using full field deformation measurements obtained experimentally using the digital image correlation method. The amount of energy release due to the investigated first crack increment is evaluated through three different approaches both for verification purposes and to produce an estimate of the portion of the energy that radiates away from the crack source in the form of transient waves. The results presented herein propose an upper bound for the energy dissipation associated to acoustic emission, which could assist the interpretation and implementation of relevant nondestructive evaluation methods and the further enrichment of the understanding of effects associated with fracture.

42 ENGINEERING↗
DOE OSTI2022

Ab initio approaches to high-entropy alloys: a comparison of CPA, SQS, and supercell methods

In this work, we present a comparative study of different modeling approaches to the electronic properties of the Hf 0.05 Nb 0.05 Ta 0.8 Ti 0.05 Zr 0.05 high-entropy alloy. Common to our modeling is the methodology to compute the one-particle Green’s function in the framework of density functional theory. We demonstrate that the special quasi-random structures modeling and the supercell, i.e., the locally self-consistent multiple-scattering methods, provide very similar results for the ground state properties such as the spectral function (density of states) and the equilibrium lattice parameter. To reconcile the multiple-scattering single-site coherent potential approximation with the real space supercell methods, we included the effect of screening of the net charges of the alloy components. Based on the analysis of the total energy and spectral functions computed within the density functional theory, we found no signature for the long-range or local magnetic moments formation in the Hf 0.05 Nb 0.05 Ta 0.8 Ti 0.05 Zr 0.05 high-entropy alloy; instead, we find possible superconductivity below ~9K.

36 MATERIALS SCIENCE↗
DOE OSTI2020

CALPHAD Uncertainty Quantification and TDBX

CALPHAD uncertainty quantification (UQ) is the foundation of materials design with quantified confidence. We report a framework and software packages to enable CALPHAD UQ assessment and calculation using commercial CALPHAD software (Thermo-Calc). This Bayesian inference framework is coupled with a Markov chain Monte Carlo algorithm to establish uncertainty traces with a given thermodynamic database file (TDB) and corresponding experimental data points. This general framework is demonstrated with the Ni–Cr binary system. The algorithm is firstly validated on synthetic data with known ground truth. Then it is applied to real experimental data to generate posterior traces. We develop a file format named TDBX, which provides a single source of truth by combining the original TDB content and the traces for each assessed Gibbs energy parameter. CALPHAD UQ calculations are performed based on the TDBX file, from which uncertainties for phase boundaries, enthalpy curves, and solidification range are collected as examples of basic design parameters. This TDBX file with corresponding scripts are made open-source. Finally, the combination of CALPHAD UQ assessments and calculations connected by TDBX supports uncertainty-assisted modeling, enabling the integrated application of modern design with uncertainty methodologies to computational materials design.

36 MATERIALS SCIENCE↗
DOE OSTI2025

SPARTA: A flux adjustment methodology to interpret complex experiments

For the accurate determination of reactivity from a detector count rate, correction of spatial effects is of prime importance. This spatial correction is often provided using simulation methodologies, but this may introduce a bias if the result of the experiment is also used as input data for the simulation. Here, this work presents a flux adjustment methodology able to infer experimental reactivity and correction of spatial effects without the need for a simulation. It can process the signal from a complex experiment such as a heat balance measurement in the TREAT reactor, where control rods are continuously adjusted to maintain a constant power. In the present work, this methodology successfully computed the reactivity and the local spatial variation of the flux of a generated signal. It also proved to be robust against noise and errors on kinetic parameters and provides a credible interpretation of a heat balance experiment in TREAT. Efficiency of flux adjustment methods for complex experiment enable a better experiment interpretation less reliant on nuclear data evaluation.

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