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,585 records · Page 88

Development of lean, efficient, and fast physics-framed deep-learning-based proxy models for subsurface carbon storage

In this work, we present deep-learning-based surrogate models for CCUS developed with four different algorithms and a physics-framed two-phase flow problem involving displacement of water by CO 2 . The deep-learning models were trained using 3D datasets describing the pressure plume, CO 2 saturation plume, and water extraction rate generated by numerical simulation. The hyperparameters defining the architecture of the neural networks were optimized to determine the slimmest network size and training parameters that give the most efficient performance at the least training cost. To develop a robust model that closely mimics the governing physical laws, the discretized form of the two-phase fluid transport equation was used to formulate the supervised deep-learning task. The algorithms investigated in this study predicted the data to above 95% accuracy, with the multi-layer perceptron model demonstrating the best performance by balancing training speed, prediction time, and prediction accuracy with lean network capacity. Furthermore, the surrogate models simultaneously predict reservoir pressure and CO 2 saturation in every grid block, including the surface well extraction rate and bottomhole pressure, at all simulation times for a given static model realization in just a few seconds on a standard desktop computer. A key outcome of this study is that limits can be placed on network design parameters to avoid over designing neural networks, with associated efficiencies in training and prediction times. This is very useful because large volumes of data may be generated in CCUS projects and over-design of neural network architectures imposes penalties that are antithetical to the goal of near-real time forecasting.

58 GEOSCIENCES↗

Multifidelity deep operator networks for data-driven and physics-informed problems

Operator learning for complex nonlinear systems is increasingly common in modeling multi-physics and multi-scale systems. However, training such high-dimensional operators requires a large amount of expensive, high-fidelity data, either from experiments or simulations. In this work, we present a composite Deep Operator Network (DeepONet) for learning using two datasets with different levels of fidelity to accurately learn complex operators when sufficient high-fidelity data is not available. Additionally, we demonstrate that the presence of low-fidelity data can improve the predictions of physics-informed learning with DeepONets. We demonstrate the new multi-fidelity training in diverse examples, including modeling of the ice-sheet dynamics of the Humboldt glacier, Greenland, using two different fidelity models and also using the same physical model at two different resolutions.

97 MATHEMATICS AND COMPUTING↗

A hybrid approach to estimate the complex motions of clouds in sky images

Tracking the motion of clouds is essential to forecasting the weather and to predicting the short-term solar energy generation. Existing techniques mainly fall into two categories: variational optical flow, and block matching. In this article, we summarize recent advances in estimating cloud motion using ground-based sky imagers and quantitatively evaluate state-of-the-art approaches. Then we propose a hybrid tracking framework to incorporate the strength of both block matching and optical flow models. To validate the accuracy of the proposed approach, we introduce a series of synthetic images to simulate the cloud movement and deformation, and thereafter comprehensively compare our hybrid approach with several representative tracking algorithms over both simulated and real images collected from various sites/imagers. The results show that our hybrid approach outperforms state-of-the-art models by reducing at least 30% motion estimation errors compared with the ground-truth motions in most of simulated image sequences. Furthermore, our hybrid model demonstrates its superior efficiency in several real cloud image datasets by lowering at least 15% Mean Absolute Error (MAE) between predicted images and ground-truth images.

54 ENVIRONMENTAL SCIENCES↗

Systematic softening in universal machine learning interatomic potentials

Machine learning interatomic potentials (MLIPs) have introduced a new paradigm for atomic simulations. Recent advancements have led to universal MLIPs (uMLIPs) that are pre-trained on diverse datasets, providing opportunities for universal force fields and foundational machine learning models. However, their performance in extrapolating to out-of-distribution complex atomic environments remains unclear. In this study, we highlight a consistent potential energy surface (PES) softening effect in three uMLIPs: M3GNet, CHGNet, and MACE-MP-0, which is characterized by energy and force underprediction in atomic-modeling benchmarks including surfaces, defects, solid-solution energetics, ion migration barriers, phonon vibration modes, and general high-energy states. The PES softening behavior originates primarily from the systematically underpredicted PES curvature, which derives from the biased sampling of near-equilibrium atomic arrangements in uMLIP pre-training datasets. Our findings suggest that a considerable fraction of uMLIP errors are highly systematic, and can therefore be efficiently corrected. We argue for the importance of a comprehensive materials dataset with improved PES sampling for next-generation foundational MLIPs.

36 MATERIALS SCIENCE↗

Enabling accurate chemical modeling of shocked energetic materials using a machine learning interatomic potential

Understanding the complex chemistry of organic materials under dynamic compression is important for many applications, but it is challenging due to the large number of reactions occurring at various time scales. Here, in this study, we develop a machine learning potential based on Chebyshev polynomials to study the insensitive energetic material 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) under detonation. We discuss a strategy for constructing diverse training data needed to capture the complex chemistry of TATB. Our potential demonstrates strong transferability across a wide range of thermodynamic conditions and other explosives, enabling accurate and reliable chemical modeling of organic materials under extreme conditions. The efficiency of our approach allows for simulations over several nanoseconds and for large system sizes, providing detailed insights into the chemistry of shocked TATB. The model accurately reproduces experimental Hugoniot equation of state data, and our simulations reveal the rapid formation of nitrogen-rich carbon clusters following shock. The methods and datasets developed here offer a robust framework for accurate chemical modeling of other shocked organic energetic materials.

Chemistry↗

A Novel Method for Preoperative Positioning of Total Ankle Replacement Using 3D Digital Model

Abstract Objective To establish a digital model of the ankle joint through 3D imaging technology and explore the preoperative placement of ankle replacement prostheses. Methods Computed tomography images of intact ankle joints from 54 cases in the outpatient and inpatient departments of our hospital were collected; according to the INBONE® total ankle system surgery process, the surgery model and surgical osteotomy were finished using MIMICS based on computer simulation method. The shortest distance was measured between the center point and the anterior, posterior, medial, and lateral, respectively, to ensure the precise position of the ankle replacement prosthesis by digital simulation surgery. The relationship between the two variables was analyzed by bivariate correlation analysis. Results The dataset of this study included 48 cases of the sub‐data set (26 males and 22 females) and included 27 cases of left ankle and 21 cases of right ankle. The average medial malleolar angle was 18.67°± 2.87°, the average amount of bone resection was 12.13 ± 1.86 cm 3 , the mid‐anterior distance was 1.72 ± 0.19 cm, the mid‐posterior distance was 2.00 ± 0.19 cm, the ratio of mid‐anterior to mid‐posterior was 0.87, the mid‐medial distance was 1.26 ± 0.17 cm, the mid‐lateral distance was 1.19 ± 0.16 cm, and the ratio of mid‐medial to mid‐lateral was 1.06. After osteotomy, the anteroposterior diameter was 3.73 ± 0.32 cm, the transverse diameter was 2.46 ± 0.27 cm, and the ratio of anteroposterior diameter to transverse diameter was 1.53. In the bottom view, the shape after osteotomy is rectangular. The mid‐anterior distance was strongly negatively correlated with age, the mid‐anterior distance and the amount of bone resection, the mid‐medial distance and the amount of bone resection, the mid‐lateral distance and the amount of bone resection, the mid‐lateral distance and the anteroposterior diameter, the anteroposterior diameter and the transverse diameter were all strongly positively correlated. Conclusion The projection point of the lower tibia centerline on the tibial horizontal osteotomy surface is located at a position slightly anterior to the midpoint of the transverse diameter after ankle arthroplasty. The rational positioning of the total ankle replacement is located at both a position slightly anterior to the midpoint of the transverse diameter and midpoint of the anteroposterior diameter, which can be used as a reference method before total ankle arthroplasty surgery.

Wu, Shi‐xun↗

SU-E-J-164: Estimation of DVH Variation for PTV Due to Interfraction Organ Motion in Prostate VMAT Using Gaussian Error Function

Purpose: We developed a method to predict the change of DVH for PTV due to interfraction organ motion in prostate VMAT without repeating the CT scan and treatment planning. The method is based on a pre-calculated patient database with DVH curves of PTV modelled by the Gaussian error function (GEF). Methods: For a group of 30 patients with different prostate sizes, their VMAT plans were recalculated by shifting their PTVs 1 cm with 10 increments in the anterior-posterior, left-right and superior-inferior directions. The DVH curve of PTV in each replan was then fitted by the GEF to determine parameters describing the shape of curve. Information of parameters, varying with the DVH change due to prostate motion for different prostate sizes, was analyzed and stored in a database of a program written by MATLAB. Results: To predict a new DVH for PTV due to prostate interfraction motion, prostate size and shift distance with direction were input to the program. Parameters modelling the DVH for PTV were determined based on the pre-calculated patient dataset. From the new parameters, DVH curves of PTVs with and without considering the prostate motion were plotted for comparison. The program was verified with different prostate cases involving interfraction prostate shifts and replans. Conclusion: Variation of DVH for PTV in prostate VMAT can be predicted using a pre-calculated patient database with DVH curve fitting. The computing time is fast because CT rescan and replan are not required. This quick DVH estimation can help radiation staff to determine if the changed PTV coverage due to prostate shift is tolerable in the treatment. However, it should be noted that the program can only consider prostate interfraction motions along three axes, and is restricted to prostate VMAT plan using the same plan script in the treatment planning system.

60 APPLIED LIFE SCIENCES↗

SU-E-J-253: The Radiomics Toolbox in the Computational Environment for Radiological Research (CERR)

Purpose: To present an open source and free platform to facilitate radiomics research — The “Radiomics toolbox” in CERR. Method: There is scarcity of open source tools that support end-to-end modeling of image features to predict patient outcomes. The “Radiomics toolbox” strives to fill the need for such a software platform. The platform supports (1) import of various kinds of image modalities like CT, PET, MR, SPECT, US. (2) Contouring tools to delineate structures of interest. (3) Extraction and storage of image based features like 1st order statistics, gray-scale co-occurrence and zonesize matrix based texture features and shape features and (4) Statistical Analysis. Statistical analysis of the extracted features is supported with basic functionality that includes univariate correlations, Kaplan-Meir curves and advanced functionality that includes feature reduction and multivariate modeling. The graphical user interface and the data management are performed with Matlab for the ease of development and readability of code and features for wide audience. Open-source software developed with other programming languages is integrated to enhance various components of this toolbox. For example: Java-based DCM4CHE for import of DICOM, R for statistical analysis. Results: The Radiomics toolbox will be distributed as an open source, GNU copyrighted software. The toolbox was prototyped for modeling Oropharyngeal PET dataset at MSKCC. The analysis will be presented in a separate paper. Conclusion: The Radiomics Toolbox provides an extensible platform for extracting and modeling image features. To emphasize new uses of CERR for radiomics and image-based research, we have changed the name from the “Computational Environment for Radiotherapy Research” to the “Computational Environment for Radiological Research”.

60 APPLIED LIFE SCIENCES↗

SU-E-T-630: Predictive Modeling of Mortality, Tumor Control, and Normal Tissue Complications After Stereotactic Body Radiotherapy for Stage I Non-Small Cell Lung Cancer

Purpose: While rates of local control have been well characterized after stereotactic body radiotherapy (SBRT) for stage I non-small cell lung cancer (NSCLC), less data are available characterizing survival and normal tissue toxicities, and no validated models exist assessing these parameters after SBRT. We evaluate the reliability of various machine learning techniques when applied to radiation oncology datasets to create predictive models of mortality, tumor control, and normal tissue complications. Methods: A dataset of 204 consecutive patients with stage I non-small cell lung cancer (NSCLC) treated with stereotactic body radiotherapy (SBRT) at the University of Pennsylvania between 2009 and 2013 was used to create predictive models of tumor control, normal tissue complications, and mortality in this IRB-approved study. Nearly 200 data fields of detailed patient- and tumor-specific information, radiotherapy dosimetric measurements, and clinical outcomes data were collected. Predictive models were created for local tumor control, 1- and 3-year overall survival, and nodal failure using 60% of the data (leaving the remainder as a test set). After applying feature selection and dimensionality reduction, nonlinear support vector classification was applied to the resulting features. Models were evaluated for accuracy and area under ROC curve on the 81-patient test set. Results: Models for common events in the dataset (such as mortality at one year) had the highest predictive power (AUC = .67, p < 0.05). For rare occurrences such as radiation pneumonitis and local failure (each occurring in less than 10% of patients), too few events were present to create reliable models. Conclusion: Although this study demonstrates the validity of predictive analytics using information extracted from patient medical records and can most reliably predict for survival after SBRT, larger sample sizes are needed to develop predictive models for normal tissue toxicities and more advanced machine learning methodologies need be consider in the future.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

TU-AB-BRC-03: Accurate Tissue Characterization for Monte Carlo Dose Calculation Using Dual-and Multi-Energy CT Data

Purpose: To develop a general method for human tissue characterization with dual-and multi-energy CT and evaluate its performance in determining elemental compositions and the associated proton stopping power relative to water (SPR) and photon mass absorption coefficients (EAC). Methods: Principal component analysis is used to extract an optimal basis of virtual materials from a reference dataset of tissues. These principal components (PC) are used to perform two-material decomposition using simulated DECT data. The elemental mass fraction and the electron density in each tissue is retrieved by measuring the fraction of each PC. A stoichiometric calibration method is adapted to the technique to make it suitable for clinical use. The present approach is compared with two others: parametrization and three-material decomposition using the water-lipid-protein (WLP) triplet. Results: Monte Carlo simulations using TOPAS for four reference tissues shows that characterizing them with only two PC is enough to get a submillimetric precision on proton range prediction. Based on the simulated DECT data of 43 references tissues, the proposed method is in agreement with theoretical values of protons SPR and low-kV EAC with a RMS error of 0.11% and 0.35%, respectively. In comparison, parametrization and WLP respectively yield RMS errors of 0.13% and 0.29% on SPR, and 2.72% and 2.19% on EAC. Furthermore, the proposed approach shows potential applications for spectral CT. Using five PC and five energy bins reduces the SPR RMS error to 0.03%. Conclusion: The proposed method shows good performance in determining elemental compositions from DECT data and physical quantities relevant to radiotherapy dose calculation and generally shows better accuracy and unbiased results compared to reference methods. The proposed method is particularly suitable for Monte Carlo calculations and shows promise in using more than two energies to characterize human tissue with CT.

60 APPLIED LIFE SCIENCES↗

Forcing data (CESM2/CMIP6) for projection of drought impacts (2015-2100) at the K34 site in Manaus, Brazil

Historical and projected output data variables extracted and derived from Community Earth System Model 2 (CESM2) runs from the Coupled Model Intercomparison Project Phase 6 (CMIP6) archive. CESM2 is a fully coupled Earth system model used in simulations of Earth's past, present, and future climates (Danabasoglu et al., 2020). Variables in this dataset are in six-hourly resolution, and include air temperature (both in K and ℃), specific humidity (kg kg-1), air pressure (Pa), relative humidity (%), and vapor pressure deficit (kPa). CESM2 was the only CMIP6 model that provided VPD at the high temporal resolution required for this analysis. Data are included in .csv files, and the text file CESM2-CMIP6_forcing_K34-Manaus_headers.txt provides descriptions of data file headers.

54 ENVIRONMENTAL SCIENCES↗

Subsurface Characterization of Hydraulic Fracture Test Site-2 (HFTS-2), Delaware Basin

Hydraulic Fracturing Test Site-2 (HFTS-2) is a field-based research experiment performed in the Wolfcamp Formation of the Permian (Delaware) Basin. This paper focuses on integration, advanced geological characterization, and 3D subsurface modeling of the comprehensive HFTS-2 dataset. The study showcases a multidisciplinary reservoir characterization approach that incorporates geology, petrophysics, geochemistry, geomechanics, microseismic, and subsurface engineering analysis. Subsurface characterization of organic-rich mudstone formations requires understanding complex hydraulic fracture network growth in relation to inherent lithology, geomechanical properties, and interaction with pre-existing natural fractures. This paper presents a characterization workflow incorporating pre- and post-stimulation subsurface data, unique to the HFTS-2 dataset. The study integrated: (1) rock properties from logs, cores, and thin sections; (2) natural and hydraulic fracture descriptions from cores and image logs; (3) local and regional stresses; (4) geomechanics; (5) microseismic; (6) fiber optic (FO) and bottomhole pressure gauge (BHPG) response; and (7) produced fluids analysis. During a stimulation treatment, creation of the stimulated rock volume (SRV) is influenced by several subsurface factors. Key contributing factors include structural context, stress conditions, lithology, facies architecture, pre-existing natural fractures, and geomechanical properties. The HFTS-2 subsurface data integration indicates that the SRV is comprised of a complex juxtaposition of hydraulic fracture swarms, as evidenced by image logs analysis, core description, and microseismic monitoring. The HFTS-2 microseismic event density was used to generate 3D heat maps that serve as a representative SRV footprint, corroborated by secondary datasets. These maps were further integrated with petrophysical and geomechanical characteristics, as well as responses from FO and BHPG, to estimate the lateral and vertical dimensions of the effective fractures. The geological characterization for the HFTS-2 dataset combined with 3D modeling for petrophysical and geomechanical properties provides a strong foundation for subsurface simulation and optimization studies. Downloaded from http://onepetro.org/URTECONF/proceedings-pdf/21URTC/1-21URTC/D011S005R001/2477501/urtec-2021-5243-ms.pdf/1 by Carol Worster on 28 February 2022 URTeC 5243 The workflow improved our understanding of HFTS-2 hydraulic fracture propagation and characteristics in relation to offset pressure depletion and interaction with pre-existing natural fractures. Analysis showed that fracture geometry varies by stage and by well, and a complex fracture network is generated with varying fracture density. The multidisciplinary workflow presented herein for integration and characterization serves as a foundation to evaluate completion efficiency and estimate areal and vertical stimulation and depletion extent for the project. Furthermore, the workflow and learnings can also be transferred to other unconventional plays.

58 GEOSCIENCES↗

Assessing the Performance of the South American Land Data Assimilation System Version 2 (SALDAS-2) Energy Balance across Diverse Biomes

Understanding the exchange of energy between the surface and the atmosphere is important in view of the climate scenario. However, it becomes a challenging task due to a sparse network of observations. This study aims to improve the energy balance estimates for the Amazon, Cerrado, and Pampa biomes located in South America using the radiation and precipitation forcing obtained from the Clouds and the Earth’s Radiant Energy System (CERES) and the precipitation CPTEC/MERGE datasets. We employed three surface models—Noah-MP, Community Land Model (CLSM), and Integrated Biosphere Simulator (IBIS)—and conducted modeling experiments, termed South America Land Data Assimilation System (SALDAS-2). The results showed that SALDAS-2 radiation estimates had the smallest errors. Moreover, SALDAS-2 precipitation estimates were better than the Global Land Data Assimilation System (GLDAS) in the Cerrado (MBE = −0.16) and Pampa (MBE = −0.19). Noah-MP presented improvements compared with CLSM and IBIS in 100% of towers located in the Amazon. CLSM tends to overestimate the latent heat flux and underestimate the sensible heat flux in the Amazon. Noah-MP and Ensemble outperformed GLDAS in terms latent and sensible heat fluxes. The potential of SALDAS-2 should be emphasized to provide more accurate estimates of surface energy balance.

modeling↗

An Improved Wake Vortex Tracking Algorithm for Multiple Aircraft

The accurate tracking of vortex evolution from Large Eddy Simulation (LES) data is a complex and computationally intensive problem. The vortex tracking requires the analysis of very large three-dimensional and time-varying datasets. The complexity of the problem is further compounded by the fact that these vortices are embedded in a background turbulence field, and they may interact with the ground surface. Another level of complication can arise, if vortices from multiple aircrafts are simulated. This paper presents a new technique for post-processing LES data to obtain wake vortex tracks and wake intensities. The new approach isolates vortices by defining "regions of interest" (ROI) around each vortex and has the ability to identify vortex pairs from multiple aircraft. The paper describes the new methodology for tracking wake vortices and presents application of the technique for single and multiple aircraft.

Switzer, George F.↗

Simulation of pseudo-CT images based on deformable image registration of ultrasound images: A proof of concept for transabdominal ultrasound imaging of the prostate during radiotherapy

Purpose: Imaging of patient anatomy during treatment is a necessity for position verification and for adaptive radiotherapy based on daily dose recalculation. Ultrasound (US) image guided radiotherapy systems are currently available to collect US images at the simulation stage (US{sub sim}), coregistered with the simulation computed tomography (CT), and during all treatment fractions. The authors hypothesize that a deformation field derived from US-based deformable image registration can be used to create a daily pseudo-CT (CT{sub ps}) image that is more representative of the patients’ geometry during treatment than the CT acquired at simulation stage (CT{sub sim}). Methods: The three prostate patients, considered to evaluate this hypothesis, had coregistered CT and US scans on various days. In particular, two patients had two US–CT datasets each and the third one had five US–CT datasets. Deformation fields were computed between pairs of US images of the same patient and then applied to the corresponding US{sub sim} scan to yield a new deformed CT{sub ps} scan. The original treatment plans were used to recalculate dose distributions in the simulation, deformed and ground truth CT (CT{sub gt}) images to compare dice similarity coefficients, maximum absolute distance, and mean absolute distance on CT delineations and gamma index (γ) evaluations on both the Hounsfield units (HUs) and the dose. Results: In the majority, deformation did improve the results for all three evaluation methods. The change in gamma failure for dose (γ{sub Dose}, 3%, 3 mm) ranged from an improvement of 11.2% in the prostate volume to a deterioration of 1.3% in the prostate and bladder. The change in gamma failure for the CT images (γ{sub CT}, 50 HU, 3 mm) ranged from an improvement of 20.5% in the anus and rectum to a deterioration of 3.2% in the prostate. Conclusions: This new technique may generate CT{sub ps} images that are more representative of the actual patient anatomy than the CT{sub sim} scan.

60 APPLIED LIFE SCIENCES↗

Realistic Simulation of Property Graphs

Building on the existing Multiplicative Attribute Graph (MAG) model, the software relaxes the assumptions of the MAG model allowing dataset expansion while replicating statistical properties.

Source record↗

WRF Output from 270m domain running simulation with no 3D morphology, NUDAPT 3D morphology, 100m resolution 3D morphology and 10m resolution 3D morphology

This is a group of four datasets that were run for an experiment testing the effect of the resolution and the coverage of 3D urban morphological inputs on meteorological output. This data is associated with the paper "Sensitivity of Urban Morphology to the Resolution of Urban Morphological Feature Inputs: Implications for Characterizing Urban Sustainability."

integrated multisector multiscale modeling↗

WRF Output from 270m domain running simulation with no 3D morphology, NUDAPT 3D morphology, 100m resolution 3D morphology and 10m resolution 3D morphology

This is a group of four datasets that were run for an experiment testing the effect of the resolution and the coverage of 3D urban morphological inputs on meteorological output. This data is associated with the paper "Sensitivity of Urban Morphology to the Resolution of Urban Morphological Feature Inputs: Implications for Characterizing Urban Sustainability."

integrated multisector multiscale modeling↗