Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation datasets”

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

Reactive Transport Simulation of Fracture Channelization and Transmissivity Evolution

Underground fractures serve as flow conduits, and they may produce unwanted migration of water and other fluids in the subsurface. An example is the migration and leakage of greenhouse gases in the context of geologic carbon sequestration. This study has generated new understanding about how acids erode carbonate fracture surfaces and the positive feedback between reaction and flow. A two-dimensional reactive transport model was developed and used to investigate the extent to which geochemical factors influence fracture permeability and transmissivity evolution in carbonate rocks. The only mineral modeled as reactive is calcite, a fast-reacting mineral that is abundant in subsurface formations. The X-ray computed tomography dataset from a previous experimental study of fractured cores exposed to carbonic acid served as a testbed to benchmark the model simulation results. The model was able to capture not only erosion of fracture surfaces but also the specific phenomenon of channelization, which produces accelerating transmissivity increase. Results corroborated experimental findings that higher reactivity of the influent solution leads to strong channelization without substantial mineral dissolution. Simulations using mineral maps of calcite in a specimen of Amherstburg limestone demonstrated that mineral heterogeneity can either facilitate or suppress the development of flow channels depending on the spatial patterns of reactive mineral. In these cases, fracture transmissivity may increase rapidly, increase slowly, or stay constant, and for all these possibilities, the calcite mineral continues to dissolve. Collectively, these results illustrate that fluid chemistry and mineral spatial patterns need to be considered in predictions of reaction-induced fracture alteration and risks of fluid migration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2016

An Efficient Ray-Tracing Method for Determining Terrain Intercepts in EDL Simulations

The calculation of a ray's intercept from an arbitrary point in space to a prescribed surface is a common task in computer simulations. The arbitrary point often represents an object that is moving according to the simulation, while the prescribed surface is fixed in a defined frame. For detailed simulations, this surface becomes complex, taking the form of real-world objects such as mountains, craters or valleys which require more advanced methods to accurately calculate a ray's intercept location. Incorporation of these complex surfaces has commonly been implemented in graphics systems that utilize highly optimized graphics processing units to analyze such features. This paper proposes a simplified method that does not require computationally intensive graphics solutions, but rather an optimized ray-tracing method for an assumed terrain dataset. This approach was developed for the Mars Science Laboratory mission which landed on the complex terrain of Gale Crater. First, this paper begins with a discussion of the simulation used to implement the model and the applicability of finding surface intercepts with respect to atmosphere modeling, altitude determination, radar modeling, and contact forces influencing vehicle dynamics. Next, the derivation and assumptions of the intercept finding method are presented. Key assumptions are noted making the routines specific to only certain types of surface data sets that are equidistantly spaced in longitude and latitude. The derivation of the method relies on ray-tracing, requiring discussion on the formulation of the ray with respect to the terrain datasets. Further discussion includes techniques for ray initialization in order to optimize the intercept search. Then, the model implementation for various new applications in the simulation are demonstrated. Finally, a validation of the accuracy is presented along with the corresponding data sets used in the validation. A performance summary of the method will be shown using the analysis from the Mars Science Laboratory's terminal descent sensing model. Alternate uses will also be shown for determining horizon maps and orbiter set times.

Shidner, Jeremy D.↗
DOE OSTI2012

Sensitivity Analysis of FEAST-Metal Fuel Performance Code: Initial Results

This memo documents the completion of the LANL milestone, M3FT-12LA0202041, describing methodologies and initial results using FEAST-Metal. The FEAST-Metal code calculations for this work are being conducted at LANL in support of on-going activities related to sensitivity analysis of fuel performance codes. The objective is to identify important macroscopic parameters of interest to modeling and simulation of metallic fuel performance. This report summarizes our preliminary results for the sensitivity analysis using 6 calibration datasets for metallic fuel developed at ANL for EBR-II experiments. Sensitivity ranking methodology was deployed to narrow down the selected parameters for the current study. There are approximately 84 calibration parameters in the FEAST-Metal code, of which 32 were ultimately used in Phase II of this study. Preliminary results of this sensitivity analysis led to the following ranking of FEAST models for future calibration and improvements: fuel conductivity, fission gas transport/release, fuel creep, and precipitation kinetics. More validation data is needed to validate calibrated parameter distributions for future uncertainty quantification studies with FEAST-Metal. Results of this study also served to point out some code deficiencies and possible errors, and these are being investigated in order to determine root causes and to improve upon the existing code models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRSDate not supplied

The SCIFLI Airborne Observation of Artemis 1 Ascent

SCIFLI, The Scientifically Calibrated In-Flight Imagery (https://scifli.larc.nasa.gov/), team at NASA Langley Research Center specializes in the collection of multispectral data for space vehicles during Ascent and Entry, Descent, and Landing (EDL) observations. Multispectral datasets are invaluable across the NASA Agency and to commercial stakeholders for evaluating the performance of launch and re-entry space vehicles and ensuring the safety of the scientific research community. The launch of Artemis 1 marked the initiation of NASA returning to lunar exploration. SCIFLI’s Airborne Multispectral Imager (SAMI) was deployed by the SCIFLI team to scientifically document the historic event and provide the NASA research community with aerial footage of the event in wavelength channels ranging from the ultraviolet (UV) to the visible (VIS) to mid-wave infrared (MWIR). SAMI was designed to spectrally image unique aerothermal phenomena during the Artemis 1 launch. The use of these datasets captured in-flight during the launch will provide insight to research organizations across the Agency; aiding in validation efforts for simulations and modeling that contributed to the kickoff of the Agency’s resurgence to lunar exploration. This presentation will focus on the imagery captured on SAMI by the WB-57 team on November 16th, 2022, during the observation. Imaging objectives were identified and considered beforehand to determine the configuration of the SAMI instrument for the imaging mission. SCIFLI and Opto-Knowledge Systems, Inc. (OKSI) performed a thorough review of the datasets collected to identify and characterize aerothermal phenomena occurring during the observation. Additional post-processing was completed to provide quantitatively calibrated temperature images of the rocket during the observation. Some of the candidates were not optimal for quantitative temperature extraction due to common degradation factors, but after various image enhancements they proved useful for qualitatively characterizing different phenomena during the observation.

Artemis1↗
DOE OSTI2022

Dataset from ORNL Flexible Research Platform (FRP)

The data comes from a two-story light commercial building which can be used to physically simulate light commercial buildings common in the nation's existing building stock. Based on a data collection plan, these measurements were performed for five different building operations. In the data set for each period, 98 data variables were selected and collected, in which 7 variables are weather data and the rest are all building and system operation data.

Cui, Borui↗
DOE OSTI2024

a priori uncertainty quantification of reacting turbulence closure models using Bayesian neural networks

While many physics-based closure model forms have been posited for the sub-filter scale (SFS) in large eddy simulation (LES), vast amounts of data available from direct numerical simulations (DNS) create opportunities to leverage data-driven modeling techniques. Albeit flexible, data-driven models still depend on the dataset and the functional form of the model chosen. Increased adoption of such models requires reliable uncertainty estimates both in the data-informed and out-of-distribution regimes. Here, in this work, we employ Bayesian neural networks (BNNs) to capture both epistemic and aleatoric uncertainties in a reacting flow model. In particular, we model the filtered progress variable scalar dissipation rate which plays a key role in the dynamics of turbulent premixed flames. We demonstrate that BNN models can provide unique insights about the structure of uncertainty of the data-driven closure models. We also propose a method for the incorporation of out-of-distribution information in a BNN, which can be used for out-of-distribution query detection. The efficacy of the model is demonstrated by a priori evaluation on a dataset consisting of a variety of flame conditions and fuels.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Advantages of a variable‐resolution global climate model in reproducing the seasonal evolution of East Asian summer monsoon

Abstract The East Asian summer monsoon (EASM) is unique among monsoon systems that it features meridional evolution of the summer monsoon. In this study, we evaluate the performances of a Variable‐Resolution Community Earth System Model (VR‐CESM) regionally refined over eastern China (14 km) in reproducing the seasonal evolution of EASM precipitation over China. Compared with reference datasets, VR‐CESM shows better performance than the corresponding globally uniform coarse‐resolution model CESM (quasi‐uniform 1°), especially over western China where complex local topography exists. The northward monsoon migration is closely related to low‐level southerly flows and vertical moisture advection, which are more reasonably simulated in VR‐CESM. The four critical timings of the EASM (monsoon onset, withdrawal, peak, and duration) are also better captured in VR‐CESM than in CESM. The corresponding spatial Pearson correlation coefficients of the four critical timings with respect to reference datasets are about 0.1 higher in VR‐CESM than those in CESM. Both models are most accurate in simulating monsoon onset and least accurate at simulating the monsoon peak. The overestimated zonal thermal contrast in CESM is responsible for the earlier monsoon onset and excessive precipitation in September over the Yangtze River valley. Finer resolution in VR‐CESM, especially over the Tibetan Plateau (TP), appears to be a main factor in simulating better zonal thermal contrast and seasonal evolution of the EASM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Micropolar Elastoplasticity Using a Fast Fourier Transform‐Based Solver

ABSTRACT This work presents a micromechanical spectral formulation for obtaining the full‐field and homogenized response of elastoplastic micropolar composites. A closed‐form radial‐return mapping is derived from thermodynamics‐based micropolar elastoplastic constitutive equations to determine the increment of plastic strain necessary to return the generalized stress state to the yield surface, and the algorithm implementation is verified using the method of numerically manufactured solutions. Then, size‐dependent material response and micro‐plasticity are shown as features that may be efficiently simulated in this micropolar elastoplastic framework. The computational efficiency of the formulation enables the generation of large datasets in reasonable computing times.

42 ENGINEERING↗
DOE OSTI2025

Development of an interatomic potential for the Ta–Li system

A new interatomic potential for the Ta–Li system is introduced to facilitate the study of phase stability, mechanical properties and non-equilibrium dynamics after Li implantation in Ta. Here, this potential is based on a generalization of the embedded atom method (GEAM) and includes contributions from embedding energy, explicit two- and three-body interactions, and nonlocal many-body interaction terms. The parameters of the potential are optimized using energies and atomic forces for a wide range of configurations obtained from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, including elastic constants, equations of state, phonon dispersion curves, point defect properties, and melting temperatures for different compositions. Although our potential is trained on a small dataset, its accuracy is comparable to that of available machine learning potentials for Li and Ta. Our simulations show that at temperatures below 500 K, Li atoms in Ta–Li alloys form clusters separated by Ta-rich domains, and we find no evidence of ordered phase formation. For Li concentrations below a few percent, Li atoms preferentially segregate to surfaces and grain boundaries. However, in alloys containing more than ~10% Li, the accumulation of Li in symmetric-tilt grain boundaries can lead to one of the following effects: formation of amorphous-like regions, changes in grain boundary structural units, or lateral movement of the grain boundary.

GEAM potential↗
DOE OSTI2021

Recurrent neural network-based multiaxial plasticity model with regularization for physics-informed constraints

We report a recurrent neural network (RNN) based model is developed as a surrogate to predict nonlinear plastic response under multiaxial loading. The RNN-based model is trained and tested on stress versus strain curves generated using a numerical solution based on the classical radial return method. Besides simply learning the basic constitutive relationship, a novel approach is taken to enforce certain physical conditions. Specifically, regularization is employed to maintain non-negative plastic power density throughout the loading history thereby ensuring monotonically increasing plastic work and thermodynamic consistency. Enforcing physics in this manner permits coupling of the data-driven RNN approach with physics-based knowledge and laws. This has the effect of reducing the necessary amount of data and ensuring known physical laws are not violated. Since, once trained, the model need not perform the expensive task of solving nonlinear equations, its efficiency is orders of magnitude greater than its numerical counterpart. The RNN-based model has been trained on varied sets of data and the accuracy on test datasets validated. The developed model is general and robust and has widespread application such as in the simulation of metal forming, large scale plasticity, and part life prediction.

42 ENGINEERING↗
DOE OSTI2012

Quantifying the Interfractional Displacement of the Gastroesophageal Junction During Radiation Therapy for Esophageal Cancer

Purpose: Accounting for interfractional changes in tumor location improves the accuracy of radiation treatment delivery. The purpose of this study was to quantify the interfractional displacement of the gastroesophageal junction (GEJ) based on standard treatment setup in patients with esophageal cancer undergoing radiation therapy. Methods and Materials: Free-breathing four-dimensional computed tomography (4D-CT) datasets were acquired weekly from 22 patients during treatment for esophageal adenocarcinoma. Scans were registered to baseline (simulation) 4D-CT scans by using bony landmarks. The distance between the center of the GEJ contour on the simulation scan and the mean location of GEJ centers on subsequent scans was used to assess changes in GEJ location between fractions; displacement was also correlated with clinical and respiratory variables. Results: The mean absolute random error was 1.69 mm (range, 0.11-4.11 mm) in the lateral direction, 1.87 mm (range, 0.51-4.09 mm) in the anterior-posterior (AP) direction, and 3.09 mm (range, 0.99-6.16 mm) in the superior-inferior (SI) direction. The mean absolute systemic GEJ displacement between fractions was 2.88 mm lateral ({>=}5 mm in 14%), mostly leftward; 2.90 mm ({>=}5 mm in 14%) AP, mostly anterior; and 6.77 mm ({>=}1 cm in 18%) SI, mostly inferior. Variations in tidal volume and diaphragmatic excursion during treatment correlated strongly with systematic SI GEJ displacement (r = 0.964, p < 0.0001; and r = 0.944, p < 0.0001, respectively) and moderately with systematic AP GEJ displacement (r = 0.678, p = 0.0005; r = 0.758, p < 0.0001, respectively). Systematic displacement in the inferior direction resulted in higher-than-intended doses ({>=}60 Gy) to the GEJ, with increased hot-spot to the adjacent stomach and lung base. Conclusion: We found large (>1-cm) interfractional displacements in the GEJ in the SI (especially inferior) direction that was not accounted for when skeletal alignment alone was used for patient positioning. Because systematic displacement in the SI direction had dosimetric impact and correlated with tidal volume, better accounting for depth of breathing is needed to reduce interfractional variability.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2025

3D strain field reconstruction by inversion of dynamical scattering

Strain governs not only the mechanical response of materials but also their electronic, optical, and catalytic properties. For this reason, the measurement of the 3D strain field is crucial for a detailed understanding and for further development of material properties through strain engineering. However, measuring strain variations along the electron beam direction has remained a major challenge for (scanning-) transmission electron microscopy (S/TEM). In this article, we present a method for 3D strain field determination using 4D-STEM. The method is based on the inversion of dynamical diffraction effects, which occur at strain field variations along the beam direction. We test the method against simulated data with a known ground truth and demonstrate its application to an experimental 4D-STEM dataset from an inclined pseudomorphically grown Al0.47Ga0.53N layer.

Niermann, Laura↗
DOE OSTI2019

Hierarchical top-down bottom-up calibration with consideration for uncertainty and inter-scale discrepancy of Peierls stress of bcc Fe

A hierarchical multiscale model of plasticity in single crystal bcc Fe, framed at the mesoscale pertaining to a statistically representative set of dislocations for each slip system, relies on specification of a set of parameters which are informed using both bottom-up (BU) (atomistic) and top-down (TD) (experimental) information. We term this specification process a 'connection' between BU and TD pathways to inform the mesoscale model. The connection is considered in the presence of error, uncertainty, and discrepancy between the models. Furthermore, we expand upon a previously developed reconciled TD and BU calibration method to account for anticipated discrepancy between information from different length scales. The results of a previously formulated likelihood-based connection test of the multiscale model suggest a 'missing link' may be responsible for part of the inter-scale discrepancy. In this case, this link is assumed to be a relation between the unit-process (single dislocation line) of coordinated kink-pair nucleation on a screw dislocation segment and the many-body dislocation process which manifests in the onset of experimentally observed plastic deformation in a single crystal. A physics-informed discrepancy layer is formulated to improve the connection in the presence of additional persistent uncertainties. This physics-informed hypothesis testing is demonstrated as an alternative to fully data-driven search methods, particularly applicable to datasets of limited size that are often encountered in such multiscale material modeling applications. Discrete dislocation dynamics simulations are used to investigate the effect of the Peierls stress on straight screw dislocations and in the operation of Frank–Read sources. Discussion concerns the critical role of uncertainty quantification in providing a basis for this form of incremental hypothesis testing.

42 ENGINEERING↗
DOE OSTI2011

S-Preconditioner for Multi-fold Data Reduction with Guaranteed User-Controlled Accuracy

The growing gap between the massive amounts of data generated by petascale scientific simulation codes and the capability of system hardware and software to effectively analyze this data necessitates data reduction. Yet, the increasing data complexity challenges most, if not all, of the existing data compression methods. In fact, lossless compression techniques offer no more than 10% reduction on scientific data that we have experience with, which is widely regarded as effectively incompressible. To bridge this gap, in this paper, we advocate a transformative strategy that enables fast, accurate, and multi-fold reduction of double-precision floating-point scientific data. The intuition behind our method is inspired by an effective use of preconditioners for linear algebra solvers optimized for a particular class of computational dwarfs (e.g., dense or sparse matrices). Focusing on a commonly used multi-resolution wavelet compression technique as the underlying solver for data reduction we propose the S-preconditioner, which transforms scientific data into a form with high global regularity to ensure a significant decrease in the number of wavelet coefficients stored for a segment of data. Combined with the subsequent EQ-calibrator, our resultant method (called S-Preconditioned EQ-Calibrated Wavelets (SPEQC-WAVELETS)), robustly achieved a 4- to 5- fold data reduction while guaranteeing user-defined accuracy of reconstructed data to be within 1% point-by-point relative error, lower than 0:01 Normalized RMSE, and higher than 0:99 Pearson Correlation. In this paper, we show the results we obtained by testing our method on six petascale simulation codes including fusion, combustion, climate, astrophysics, and subsurface groundwater in addition to 13 publicly available scientific datasets. We also demonstrate that application-driven data mining tasks performed on decompressed variables or their derived quantities produce results of comparable quality with the ones for the original data.

Jin, Ye↗
DOE OSTI2020

eFESTA: Ensemble Feature Exploration with Surface Density Estimates

We propose surface density estimate (SDE) to model the spatial distribution of surface features-isosurfaces, ridge surfaces, and streamsurfaces-in 3D ensemble simulation data. The inputs of SDE computation are surface features represented as polygon meshes, and no field datasets are required (e.g., scalar fields or vector fields). The SDE is defined as the kernel density estimate of the infinite set of points on the input surfaces and is approximated by accumulating the surface densities of triangular patches. We also propose an algorithm to guide the selection of a proper kernel bandwidth for SDE computation. An ensemble Feature Exploration method based on Surface densiTy EstimAtes (eFESTA) is then proposed to extract and visualize the major trends of ensemble surface features. For an ensemble of surface features, each surface is first transformed into a density field based on its contribution to the SDE, and the resulting density fields are organized into a hierarchical representation based on the pairwise distances between them. The hierarchical representation is then used to guide visual exploration of the density fields as well as the underlying surface features. We demonstrate the application of our method using isosurface in ensemble scalar fields, Lagrangian coherent structures in uncertain unsteady flows, and streamsurfaces in ensemble fluid flows.

Computational modeling↗
DOE OSTI2025

Quantifying and Zoning Urban Heat Island Effects Using Unsupervised Machine Learning

This work explores the Urban Heat Island (UHI) effects in Maricopa County, Arizona, employing a simulation-based approach that combines large-scale building energy modeling with advanced spatial analysis. Utilizing the Automatic Building Energy Modeling (AutoBEM) software suite, we simulated the energy consumption for approximately 1.35 million buildings based on the Model America version 1.0 (MAv1) dataset. Our methodology incorporated spatial analysis at multiple scales, including individual buildings, clusters of zones determined by K-means clustering, and geographical level evaluation based on Zip codes. The results revealed significant variations in energy consumption and heat emissions across different building types and urban zones. High-emission hotspots identified through clustering pointed to areas most contributing to the UHI effects. Zip code-based area analysis further contextualized these findings, offering an urban context-based perspective on emission distribution and informing potential urban energy policies for mitigating UHI effects.

Chowdhury, Shovan [ORNL]↗
DOE OSTI2020

SYSTEMATIC MULTI-INDEX VALIDATIONS OF SIDI ENGINE FLOW FIELD LES COMPUTATIONS USING CRANK ANGLE-RESOLVED PIV MEASUREMENTS

In-cylinder flow fields make significant impacts on the fuel atomization, fuel mixture formation, and combustion process in spark ignition direct injection (SIDI) engines. In recent years, model-based simulation approaches are preferred in regard to investigating the transient in-cylinder flow field characteristics. Most commonly, the simulation models are validated using single representative flow field at a typical crank angle measured by particle image velocimetry (PIV). However, it provides only limited knowledge about the flow field which is highly three-dimensional and of transient nature. In this study, crank angle-resolved Ply measurements are conducted on three distinct planes inside the cylinder to capture the transient process of flow field characteristics which vary with the crank angle. These three planes consist of one tumble plane through the spark plug tip, one tumble plane along two intake ports, and one swirl plane at 30 mm below the cylinder head. Large eddy simulation (LES) is employed for the numerical computations using the CONVERGE codes. On the basis of large datasets for both temporal and spatial domains, a multi-index systematic validation approach is conducted. Crank angle-resolved calculations of global indices and local indices are implemented using the flow fields velocity data obtained from both PIV and LES on select planes. Global indices reveal the trends in similarities of different crank angle degrees and locations, while local indices give the detail comparison results. In summary, with the systematic multi-index validation approach, the crucial crank angle degrees and locations for model verification will be detected. Furthermore, the corresponding critical flow features are analyzed. Practical guideline of flow field validation is proposed.

Crank angle-resolved PIV↗
DOE OSTI2015

WE-G-207-06: 3D Fluoroscopic Image Generation From Patient-Specific 4DCBCT-Based Motion Models Derived From Physical Phantom and Clinical Patient Images

Purpose: Respiratory-correlated cone-beam CT (4DCBCT) images acquired immediately prior to treatment have the potential to represent patient motion patterns and anatomy during treatment, including both intra- and inter-fractional changes. We develop a method to generate patient-specific motion models based on 4DCBCT images acquired with existing clinical equipment and used to generate time varying volumetric images (3D fluoroscopic images) representing motion during treatment delivery. Methods: Motion models are derived by deformably registering each 4DCBCT phase to a reference phase, and performing principal component analysis (PCA) on the resulting displacement vector fields. 3D fluoroscopic images are estimated by optimizing the resulting PCA coefficients iteratively through comparison of the cone-beam projections simulating kV treatment imaging and digitally reconstructed radiographs generated from the motion model. Patient and physical phantom datasets are used to evaluate the method in terms of tumor localization error compared to manually defined ground truth positions. Results: 4DCBCT-based motion models were derived and used to generate 3D fluoroscopic images at treatment time. For the patient datasets, the average tumor localization error and the 95th percentile were 1.57 and 3.13 respectively in subsets of four patient datasets. For the physical phantom datasets, the average tumor localization error and the 95th percentile were 1.14 and 2.78 respectively in two datasets. 4DCBCT motion models are shown to perform well in the context of generating 3D fluoroscopic images due to their ability to reproduce anatomical changes at treatment time. Conclusion: This study showed the feasibility of deriving 4DCBCT-based motion models and using them to generate 3D fluoroscopic images at treatment time in real clinical settings. 4DCBCT-based motion models were found to account for the 3D non-rigid motion of the patient anatomy during treatment and have the potential to localize tumor and other patient anatomical structures at treatment time even when inter-fractional changes occur. This project was supported, in part, through a Master Research Agreement with Varian Medical Systems, Inc., Palo Alto, CA. The project was also supported, in part, by Award Number R21CA156068 from the National Cancer Institute.

60 APPLIED LIFE SCIENCES↗