Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation techniques”

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

Machine learning-based surrogate modeling for data-driven optimization: a comparison of subset selection for regression techniques

Optimization of simulation-based or data-driven systems is a challenging task, which has attracted significant attention in the recent literature. A very efficient approach for optimizing systems without analytical expressions is through fitting surrogate models. Due to their increased flexibility, nonlinear interpolating functions, such as radial basis functions and Kriging, have been predominantly used as surrogates for data-driven optimization; however, these methods lead to complex nonconvex formulations. Alternatively, commonly used regression-based surrogates lead to simpler formulations, but they are less flexible and inaccurate if the form is not known a priori. In this work, we investigate the efficiency of subset selection regression techniques for developing surrogate functions that balance both accuracy and complexity. Subset selection creates sparse regression models by selecting only a subset of original features, which are linearly combined to generate a diverse set of surrogate models. Five different subset selection techniques are compared with commonly used nonlinear interpolating surrogate functions with respect to optimization solution accuracy, computation time, sampling requirements, and model sparsity. Furthermore, our results indicate that subset selection-based regression functions exhibit promising performance when the dimensionality is low, while interpolation performs better for higher dimensional problems.

42 ENGINEERING↗
NASA NTRS2008

Reducing Spaceborne-Doppler-Radar Rainfall-Velocity Error

A combined frequency-time (CFT) spectral moment estimation technique has been devised for calculating rainfall velocity from measurement data acquired by a nadir-looking spaceborne Doppler weather radar system. Prior spectral moment estimation techniques used for this purpose are based partly on the assumption that the radar resolution volume is uniformly filled with rainfall. The assumption is unrealistic in general but introduces negligible error in application to airborne radar systems. However, for spaceborne systems, the combination of this assumption and inhomogeneities in rainfall [denoted non-uniform beam filling (NUBF)] can result in velocity measurement errors of several meters per second. The present CFT spectral moment estimation technique includes coherent processing of a series of Doppler spectra generated in a standard manner from data over measurement volumes that are partially overlapping in the along-track direction. Performance simulation of this technique using high-resolution data from an airborne rain-mapping radar shows that a spaceborne Ku-band Doppler radar operating at signal-to-noise ratios greater than 10 dB can achieve root-mean-square accuracy between 0.5 and 0.6 m/s in vertical-velocity estimates.

Tanelli, Simone↗
DOE OSTI2022

Stable Element Doping of Sol-gel Toward Simulating Environmental Matrix in Surrogate Explosive Nuclear Debris

Training nuclear forensic analysis personnel in post-detonation scenarios is of critical importance to nuclear threat response capabilities. Thus, realistic nuclear debris simulants which resemble the size, color, elemental composition, and radionuclide content of actual nuclear fallout from a recent detonation would be valuable for training nuclear first responders in realistic scenarios. As nuclear fallout types vary significantly in each of these parameters based on detonation environment (rural, urban, maritime etc.) and collection location, the ability to tailor each of these parameters accurately in simulated debris would be of immense benefit to the post-detonation analysis community for training both in-field collections and triage and the validation of laboratory level nuclear forensic techniques. Sol-gel synthesis techniques can provide the tunability of size, shape and composition required for producing surrogate nuclear debris of a wide variety. The sol-gel process consists of forming a metal oxide material, often silica, through polymerization of a metal-alkoxy precursor. In this work, we characterize the ability to load the sol-gel particles with secondary elemental components such as iron, aluminum, and calcium toward approximating the elemental composition of debris from various detonation environments and demonstrate the ability to produce particles with controllable size, shape, and color. We also demonstrate quantitative radionuclide encapsulation toward reproducing the radionuclide content of actual fallout from a recent detonation. Finally, we then employ these techniques in producing simulated aerodynamic debris samples with realistic elemental matrix composition and radionuclide content simulating a recent uranium-fueled detonation taking place in a rural desert environment and compare it to historic fallout from the Nevada Nuclear Security Site.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Stable Element Doping of Sol-gel Toward Simulating Environmental Matrix in Surrogate Explosive Nuclear Debris

Training nuclear forensic analysis personnel in post-detonation scenarios is of critical importance to nuclear threat response capabilities. Thus, realistic nuclear debris simulants which resemble the size, color, elemental composition, and radionuclide content of actual nuclear fallout from a recent detonation would be valuable for training nuclear first responders in realistic scenarios. As nuclear fallout types vary significantly in each of these parameters based on detonation environment (rural, urban, maritime etc.) and collection location, the ability to tailor each of these parameters accurately in simulated debris would be of immense benefit to the post-detonation analysis community for training both in-field collections and triage and the validation of laboratory level nuclear forensic techniques. Sol-gel synthesis techniques can provide the tunability of size, shape and composition required for producing surrogate nuclear debris of a wide variety. The sol-gel process consists of forming a metal oxide material, often silica, through polymerization of a metal-alkoxy precursor. In this work, we characterize the ability to load the sol-gel particles with secondary elemental components such as iron, aluminum, and calcium toward approximating the elemental composition of debris from various detonation environments and demonstrate the ability to produce particles with controllable size, shape, and color. We also demonstrate quantitative radionuclide encapsulation toward reproducing the radionuclide content of actual fallout from a recent detonation. Finally, we then employ these techniques in producing simulated aerodynamic debris samples with realistic elemental matrix composition and radionuclide content simulating a recent uranium-fueled detonation taking place in a rural desert environment and compare it to historic fallout from the Nevada Nuclear Security Site.

36 MATERIALS SCIENCE↗
NASA NTRS1992

A brief overview of computational structures technology related activities at NASA Lewis Research Center

The presentation gives a partial overview of research and development underway in the Structures Division of LeRC, which collectively is referred to as the Computational Structures Technology Program. The activities in the program are diverse and encompass four major categories: (1) composite materials and structures; (2) probabilistic analysis and reliability; (3) design optimization and expert systems; and (4) computational methods and simulation. The approach of the program is comprehensive and entails exploration of fundamental theories of structural mechanics to accurately represent the complex physics governing engine structural performance, formulation, and implementation of computational techniques and integrated simulation strategies to provide accurate and efficient solutions of the governing theoretical models by exploiting the emerging advances in computer technology, and validation and verification through numerical and experimental tests to establish confidence and define the qualities and limitations of the resulting theoretical models and computational solutions. The program comprises both in-house and sponsored research activities. The remainder of the presentation provides a sample of activities to illustrate the breadth and depth of the program and to demonstrate the accomplishments and benefits that have resulted.

Hopkins, Dale A.↗
NASA NTRS2015

Check-Cases for Verification of 6-Degree-of-Freedom Flight Vehicle Simulations

The rise of innovative unmanned aeronautical systems and the emergence of commercial space activities have resulted in a number of relatively new aerospace organizations that are designing innovative systems and solutions. These organizations use a variety of commercial off-the-shelf and in-house-developed simulation and analysis tools including 6-degree-of-freedom (6-DOF) flight simulation tools. The increased affordability of computing capability has made highfidelity flight simulation practical for all participants. Verification of the tools' equations-of-motion and environment models (e.g., atmosphere, gravitation, and geodesy) is desirable to assure accuracy of results. However, aside from simple textbook examples, minimal verification data exists in open literature for 6-DOF flight simulation problems. This assessment compared multiple solution trajectories to a set of verification check-cases that covered atmospheric and exo-atmospheric (i.e., orbital) flight. Each scenario consisted of predefined flight vehicles, initial conditions, and maneuvers. These scenarios were implemented and executed in a variety of analytical and real-time simulation tools. This tool-set included simulation tools in a variety of programming languages based on modified flat-Earth, round- Earth, and rotating oblate spheroidal Earth geodesy and gravitation models, and independently derived equations-of-motion and propagation techniques. The resulting simulated parameter trajectories were compared by over-plotting and difference-plotting to yield a family of solutions. In total, seven simulation tools were exercised.

Murri, Daniel G.↗
NASA NTRS2015

Challenges of NDE Simulation Tool Challenges of NDE Simulation Tool

Realistic nondestructive evaluation (NDE) simulation tools enable inspection optimization and predictions of inspectability for new aerospace materials and designs. NDE simulation tools may someday aid in the design and certification of advanced aerospace components; potentially shortening the time from material development to implementation by industry and government. Furthermore, modeling and simulation are expected to play a significant future role in validating the capabilities and limitations of guided wave based structural health monitoring (SHM) systems. The current state-of-the-art in ultrasonic NDE/SHM simulation cannot rapidly simulate damage detection techniques for large scale, complex geometry composite components/vehicles with realistic damage types. This paper discusses some of the challenges of model development and validation for composites, such as the level of realism and scale of simulation needed for NASA' applications. Ongoing model development work is described along with examples of model validation studies. The paper will also discuss examples of the use of simulation tools at NASA to develop new damage characterization methods, and associated challenges of validating those methods.

Leckey, Cara A. C.↗
NASA NTRS2019

Quadruple Simulations of Thermographic Inspections of Impacted Composites

Thermography has been shown to be a viable technique for inspection of composites. Impact damage in composites typically contains multiple overlapping delaminations at different depths. Understanding the limitations of the thermographic inspection is enhanced by performing simulations of the technique. Most simulations of composite thermographic inspections have focused on simulations of a single delamination at a fixed depth. The quadrupole method has been shown as a viable technique for rapid three-dimensional thermographic simulations of a delamination. This method is expanded to enable rapid simulation of multiple overlapping delaminations at different depths. Quadrupole simulations are compared to finite element simulations of multiple delaminations at different depths. The simulations are also compared to the thermographic measurements on impacted composites where shape and depth of the delaminations are known from x-ray computed tomography data.

Winfree, William P.↗
DOE OSTI2023

Modeling the 4D discharge of lithium-ion batteries with a multiscale time-dependent deep learning framework

The lithium-ion battery (LIB) field is moving towards the direction of investigating spatially resolved physical phenomena in the 3D porous microstructure of electrodes. These pore-scale simulations give new insights into the local dynamics of lithiation/de-lithiation and charge transport, Nevertheless, the computational time of these simulations limits the integration of these models in optimization workflows of cycling conditions or electrode manufacturing processes. Machine learning models present a way of assessing in real-time the performance of materials. While several successful techniques for replicating simulations with machine learning have been proposed, this case study presents a more demanding problem, due to the necessity of understanding the behavior of heterogeneous 3D local data, as it evolves in time: this poses both a scientific and a technical challenge. To this end, we propose an autoregressive multiscale convolutional neural network model to predict relevant quantities at the pore-scale in the solid phase: the lithium concentration (in the active material) and potential (in the active material and carbon binder). Here, these are ultimately used to reconstruct the battery discharge curve. 3D images of the electrode microstructures are the input to the network, trained with a dataset of finite element method simulations to predict the discharge behavior of the cathode side in lithium ion batteries. We propose this machine learning model as a proof-of-concept of the applicability of multiscale networks for time-dependent physics problems. The trained model exhibits very high accuracy (with errors lower than 2 %) in forecasting the discharge behavior of new unseen cathodes.

25 ENERGY STORAGE↗
NASA NTRS1987

Evaluation of constraint stabilization procedures for multibody dynamical systems

Comparative numerical studies of four constraint treatment techniques for the simulation of general multibody dynamic systems are presented, and results are presented for the example of a classical crank mechanism and for a simplified version of the seven-link manipulator deployment problem. The staggered stabilization technique (Park, 1986) is found to yield improved accuracy and robustness over Baumgarte's (1972) technique, the singular decomposition technique (Walton and Steeves, 1969), and the penalty technique (Lotstedt, 1979). Furthermore, the staggered stabilization technique offers software modularity, and the only data each solution module needs to exchange with the other is a set of vectors plus a common module to generate the gradient matrix of the constraints, B.

Park, K. C.↗
DOE OSTI2026

Observable optimization for precision theory: machine learning energy correlators

The practice of collider physics typically involves the marginalization of multi-dimensional collider data to uni-dimensional observables relevant for some physics task. In many cases, such as classification or anomaly detection, the observable can be arbitrarily complicated, such as the output of a neural network. However, for precision measurements, the observable must correspond to something computable systematically beyond the level of current simulation tools. In this work, we demonstrate that precision-theory-compatible observable space exploration can be systematized by using neural simulation-based inference techniques from machine learning. We illustrate this approach by exploring the space of marginalizations of the energy 3-point correlator to optimize sensitivity to the top quark mass. We first learn the energy-weighted probability density from simulation, then search in the space of marginalizations for an optimal triangle shape. Although simulations and machine learning are used in the process of observable optimization, the output is an observable definition which can be then computed to high precision and compared directly to data without any memory of the computations which produced it. We find that the optimal marginalization is isosceles triangles on the sphere with a side ratio approximately $1 : 1 : \sqrt{2}$ (i.e. right triangles) within the set of marginalizations we consider.

Jets and Jet Substructure↗
DOE OSTI2024

Phase-space methods for neutrino oscillations: Extension to multibeams

The phase-space approach (PSA), which was originally introduced in Lacroix [] to describe neutrino flavor oscillations for interacting neutrinos emitted from stellar objects is extended to describe arbitrary numbers of neutrino beams. The PSA is based on mapping the quantum fluctuations into a statistical treatment by sampling initial conditions followed by independent mean-field evolution. A new method is proposed to perform this sampling that allows treating an arbitrary number of neutrinos in each neutrino beams. We validate the technique successfully and confirm its predictive power on several examples where a reference exact calculation is possible. We show that it can describe many-body effects, such as entanglement and dissipation induced by the interaction between neutrinos. Due to the complexity of the problem, exact solutions can only be calculated for rather limited cases, with a limited number of beams and/or neutrinos in each beam. The PSA approach considerably reduces the numerical cost and provides an efficient technique to accurately simulate arbitrary numbers of beams. Examples of PSA results are given here, including up to 200 beams with time-independent or time-dependent Hamiltonians. We anticipate that this approach will be useful to bridge exact microscopic techniques with more traditional transport theories used in neutrino oscillations. It will also provide important reference calculations for future quantum computer applications where other techniques are not applicable to classical computers. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1978

Thematic mapper design parameter investigation

Simulated multispectral data sets used in the Thematic Mapper design parameter investigation are analyzed. The sets were produced by processing high-resolution digital multispectral data through models of the various Thematic Mapper configurations and then through three resampling techniques as a simulation of ground processing. Emphasis is placed on classification processing and analysis of results.

Colby, C. P., Jr.↗
NASA NTRS1989

Three-dimensional simulation of underexpanded plumes using upwind algorithms

Different numerical techniques have been examined for simulating the underexpanded plumes in supersonic flows. The results revealed that the monotone differentiable limiter is superior to the TVD (total variation diminishing) minmod limiter when the MUSCL (monotonic upstream schemes for conservation laws) approach is used for Roe's upwind scheme. Roe's upwind scheme is found to be superior to the central difference scheme not only for capturing shock waves, but also for capturing shear flows.

Obayashi, Shigeru↗
DOE OSTI2021

Ab Initio Molecular Dynamics Simulations of Amorphous Calcium Carbonate: Interpretation of Pair Distribution Function and X-Ray Absorption Spectroscopy Data

The structure and transformation of hydrous amorphous calcium carbonate (ACC) are key to understanding biomineralization pathways and their relationship with the properties of the resulting material. Quantitative interpretation of scattering experiments aimed at elucidating the structure of ACC is challenging, due to the amorphous nature of this material, and, therefore, requires models for the structure and the scattering physics. Here, we generate physically realistic ensembles of hydrated ACC structures and their vibrational disorder from ab initio molecular dynamics (AIMD) simulations with an emphasis on enabling the consistent interpretation of the finer details of three complementary structural probes: neutron and x-ray pair distribution experiments and x-ray absorption spectroscopy (XAS). In each case, we simulate the signal directly in reciprocal space and then manipulate it into the real-space pair distribution function (PDF) or spectrum using the same procedures for the experimental and theoretical data. Good agreement with experiment was obtained across the three techniques with the simulations accounting well for all features in the spectra. Remaining small discrepancies pointed to differences between real samples and the idealized simulated systems such as deviations from the nominal CaCO3·nH2O stoichiometry. Additionally, the simulations offered a more accurate description of the local coordination environment of calcium than previous shell-by-shell fits to spectra of synthetic ACC and classical molecular dynamics simulations. This work demonstrates that AIMD is a powerful approach for extracting detailed structural information from neutron PDF, x-ray PDF, and XAS of amorphous carbonate phases.

Prange, Micah P.↗
DOE OSTI2018

The F-t-Pj-RG method: An adjacent-rolling-windows based steady-state detection technique for application to kinetic Monte Carlo simulations

A window-based steady-state detection algorithm has been developed for application to kinetic Monte Carlo simulation data. The algorithm, termed F-t-Pj-RG sequentially applies an F-test, a t-test, and a projection test on adjacent windows of the data while rolling (or shifting) and growing the windows when any of the tests fail. In aggregate, the algorithm is able to (a) automatically reject the warm-up period as not being at steady-state, as well as (b) determine an appropriate window size for converged statistics when sampling the data, which is necessary for detection of steady-state, and (c) detect steady-state within a particular tolerance. The last step, the projection test, is actually an oscillating-slope projection test, and is performed on j sequential data windows (i.e., more than two adjacent windows). It requires more than simply being within the user defined tolerance: the oscillating-slope projection test includes a condition that the slope must oscillate around zero when 2, which is an additional indication of steady-state. When all three tests are passed, the F-t-Pj test is passed, indicating that the prerequisites of steady-state detection have been met and also that conditions consistent with the definition of steady-state have been realized. This algorithm is applied to a variety of data sets that correspond to the diverse type of data trends that can be produced by kinetic Monte Carlo simulations. The algorithm is shown to be robust in its ability to handle differing functional forms, and is able to detect steady-state with low computational cost. Finally, the low computational cost of this method and its robustness towards varied data trends make it suitable for on-the-fly use in kinetic Monte Carlo simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Discrete ordinates analysis of the forced-flight variance reduction technique in Monte Carlo neutral particle transport simulations

This paper presents mathematical formulations and methods to predict the effect of forced-flight variance reduction on Monte Carlo tally variance and calculation time. This includes deducing biasing operators that are then used to construct a history-score probability density function (HSPDF), which represents all possible Monte Carlo random walks and gives the probability of a Monte Carlo history scoring in a tally from a particular phase-space position. The history-score moment equations (HSMEs), the statistical moments of the HSPDF, are then derived to calculate the statistical behavior of the Monte Carlo tally when forced-flight variance reduction is applied. In addition, the future-time equation (FTE) is derived to predict the Monte Carlo computational time as a result of applying forced-flight variance reduction. The solutions of the HSMEs and FTE can be used to predict Monte Carlo computational cost. Furthermore, this work also describes a discrete ordinates method to solve the forced-flight HSMEs and FTE. Several 1-D and 2-D test problems verify that the derivations are performed and implemented correctly.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Comparison of Conservation Behavior in Several Variationally Derived Implementations for Simulating Electromagnetic Plasmas

Variational techniques based on Low's Lagrangian formulation have demonstrated significant improvement in conservation properties as compared to traditional particle-in-cell (PIC) algorithms. Previous work showed that equations derived from a discretized Lagrangian preserve a discretized version of the connection between Lagrangian symmetries and conservation laws, i.e., Noether's theorem. Specifically, a thorough analysis of energy conservation within the new system explained the absence of “grid heating” phenomena. In the present work, variationally derived implementations in a gridded periodic domain were compared to analogous implementations using Fourier bases. Since the Fourier basis allows for exact momentum conservation, this approach enables us to compare the relative merits of exact and approximate momentum conservation. Furthermore, to establish a tangible comparison for numerical methods, we investigate two simple physical examples; (1) the coupling of an electrostatic and electromagnetic plasma wave and (2) an equilibrium supporting Weibel instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗