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,495 records · Page 83

Towards district scale flood simulations using conventional and anisotropic porosity shallow water models with high-resolution topographic information

Current topographic survey technology provides high-resolution (HR) datasets for urban environments. Incorporating this HR information in models aiming to provide flood risk assessment is desirable because the flood wave propagation is depending on the urban topographic features, i.e. buildings, bridges and street networks. Conceptual, numerical and practical challenges arise from the application of shallow water models to HR urban flood modeling. For example, numerical challenges are occurrence of wet-dry fronts, geometric discontinuities in the urban environment and discontinuous solutions, i.e. shock waves. These challenges can be overcome by using a Godunov-type scheme. However, the computational cost of this type of schemes is high, such that HR two-dimensional shallow water simulations with practical relevance have to be run on supercomputers. The porous shallow water model is an alternative approach that aims to reduce computational cost by using a coarse resolution and accounting for unresolved processes by means of the porosity terms. Usually, a speedup between two and three orders of magnitude in comparison to HR simulations can be obtained. This study presents preliminary results of a practical test case concerning pluvial flooding in a district of the city of Nice, France, caused by the intense rainfall event on October 3rd, 2015. HR topography data set on a 1 m resolution is available for the district, whereby street features of infra-metric dimensions have been included. A reference solution is calculated by a HR shallow water model on a 1 m by 1 m structured computational grid. The porous shallow water model is run on a 10 m by 10 m grid and the influence of the drag source term is studied. The model calculations show a large deviation, which is caused by the poor meshing strategy of the porous shallow water (AP) model. The study also summarizes practical challenges that arise during the application of the AP and HR models to a large urban catchment. The main difficulty is to obtain a good mesh. In smaller scale investigations, the mesh is currently constructed by hand such that the cell edges align with buildings. This approach is not feasible for large scale urban catchments with a large number of buildings. Future steps that have to be taken, such as a strategy for automatic mesh generation, are reported on.

42 ENGINEERING↗

Global dust model intercomparison in AeroCom phase I

This study presents the results of a broad intercomparison of a total of 15 global aerosol models within the AeroCom project. Each model is compared to observations related to desert dust aerosols, their direct radiative effect, and their impact on the biogeochemical cycle, i.e., aerosol optical depth (AOD) and dust deposition. Additional comparisons to Angström exponent (AE), coarse mode AOD and dust surface concentrations are included to extend the assessment of model performance and to identify common biases present in models. These data comprise a benchmark dataset that is proposed for model inspection and future dust model development. There are large differences among the global models that simulate the dust cycle and its impact on climate. In general, models simulate the climatology of vertically integrated parameters (AOD and AE) within a factor of two whereas the total deposition and surface concentration are reproduced within a factor of 10. In addition, smaller mean normalized bias and root mean square errors are obtained for the climatology of AOD and AE than for total deposition and surface concentration. Characteristics of the datasets used and their uncertainties may influence these differences. Large uncertainties still exist with respect to the deposition fluxes in the southern oceans. Further measurements and model studies are necessary to assess the general model performance to reproduce dust deposition in ocean regions sensible to iron contributions. Models overestimate the wet deposition in regions dominated by dry deposition. They generally simulate more realistic surface concentration at stations downwind of the main sources than at remote ones. Most models simulate the gradient in AOD and AE between the different dusty regions. However the seasonality and magnitude of both variables is better simulated at African stations than Middle East ones. The models simulate the offshore transport of West Africa throughout the year but they overestimate the AOD and they transport too fine particles. The models also reproduce the dust transport across the Atlantic in the summer in terms of both AOD and AE but not so well in winter-spring nor the southward displacement of the dust cloud that is responsible of the dust transport into South America. Based on the dependency of AOD on aerosol burden and size distribution we use model bias with respect to AOD and AE to infer the bias of the dust emissions in Africa and the Middle East. Finally, according to this analysis we suggest that a range of possible emissions for North Africa is 400 to 2200 Tg yr -1 and in the Middle East 26 to 526 Tg yr -1 .

54 ENVIRONMENTAL SCIENCES↗

Utah FORGE: Temperature-Dependent Fracture Seismicity from Fluid Injection Experiments

This dataset contains experimental data from fluid injection experiments conducted to investigate the influence of temperature on fracture seismicity. The experiments were performed on granite samples from Utah FORGE. The samples were prepared with a 30-degree inclined fracture and subjected to controlled stress and temperature conditions. Data were collected under three distinct temperature settings: 24 C, 78 C, and 137 C. During the experiments, a constant confining pressure of 10 MPa and a constant shear stress at 80% of the shear strength of the sample were maintained. Pore pressure was incrementally increased at a rate of 300 kPa every three minutes to simulate fluid injection. Temperature was raised rapidly and then stabilized for the duration of each test. The dataset includes shear stress and displacement measurements under each temperature condition, along with supplementary figures illustrating the experimental setup and time-series plots of pressures and temperature.

15 GEOTHERMAL ENERGY↗

BuildStockQuery [SWR-23-58]

BuildStockQuery is a python library designed to simplify and streamline the process of querying massive, terabyte-scale datasets generated by ResStock(TM). ResStock (SWR-19-15) is a U.S. DOE-supported, NREL-built, national residential building energy stock model that enables a new approach to large-scale residential energy analysis across the U.S. by combining large public and private data sources, statistical sampling, detailed sub-hourly building simulations, and high-performance computing. BuildStockQuery offers an intuitive Object-Oriented Programming (OOP) interface to the ResStock output dataset allowing users to easily perform common queries and receive results in familiar pandas DataFrame format, abstracting away the need for complex SQL query. By initializing a query object with the pertinent Athena database and table names, users can easily query for various kinds of insights, for example, timeseries electricity for an end use for a given state grouped by building types.

Adhikari, Rajendra↗

Sci-Thur PM - Colourful Interactions: Highlights 04: A Fast Quantitative MRI Acquisition and Processing Pipeline for Radiation Treatment Planning and Simulation

MRI-only Radiation Treatment Planning (RTP) is becoming increasingly popular because of a simplified work-flow, and less inconvenience to the patient who avoids multiple scans. The advantages of MRI-based RTP over traditional CT-based RTP lie in its superior soft-tissue contrast, and absence of ionizing radiation dose. The lack of electron-density information in MRI can be addressed by automatic tissue classification. To distinguish bone from air, which both appear dark in MRI, an ultra-short echo time (UTE) pulse sequence may be used. Quantitative MRI parametric maps can provide improved tissue segmentation/classification and better sensitivity in monitoring disease progression and treatment outcome than standard weighted images. Superior tumor contrast can be achieved on pure T{sub 1} images compared to conventional T{sub 1}-weighted images acquired in the same scan duration and voxel resolution. In this study, we have developed a robust and fast quantitative MRI acquisition and post-processing work-flow that integrates these latest advances into the MRI-based RTP of brain lesions. Using 3D multi-echo FLASH images at two different optimized flip angles (both acquired in under 9 min, and 1mm isotropic resolution), parametric maps of T{sub 1}, proton-density (M{sub 0}), and T{sub 2}{sup *} are obtained with high contrast-to-noise ratio, and negligible geometrical distortions, water-fat shifts and susceptibility effects. An additional 3D UTE MRI dataset is acquired (in under 4 min) and post-processed to classify tissues for dose simulation. The pipeline was tested on four healthy volunteers and a clinical trial on brain cancer patients is underway.

60 APPLIED LIFE SCIENCES↗

Performance Year 1 Technical Report - OPEN COG Grid: Extendable Coherent Models-Datasets for Cognitive Power Grids

The OPEN COG Grid project is a collaborative effort between LLNL, NREL, and Texas A&M University (TAMU) to develop synthetic power system datasets that (i) contain all technical information that would be available in a real system, allowing to conduct studies ranging from dynamic simulation to long term planning studies; ii) are accessible to researchers from the broader data sciences community, as oppossed to power system experts only; and (iii) This report summarizes the work conducted during the first 15 months of execution of the project. These activities encompassed: 1. Conduct a survey of existing open data sets and open source power systems simulators, their supported use cases, and accessibility (Chapter 1). 2. Define a new extensible specification for power system data, covering all parameters necessary for most computational use cases (Chapter 2). 3. Collecting real technical system data to complete missing parameters in existing open source datasets (Chapter 3). 4. Develop models that capture the behavior of emergent actors in power grids, neglected by existing datasets; aggregated residential demand response (Chapter 4) and demand response of cryptocurrency miners (Chapter 5). 5. Collect detailed spatial information on distributed energy resources, particular, solar photovoltaic facilities (Chapter 6). The following chapters provide detailed descriptions of these tasks, the assumptions taken, and their findings. In conducting these tasks, the project team produced: two (accepted) conference papers; one journal paper under submission; one draft journal paper pending submission; released one repository with the developed power system data specification, with documentation and examples; and one extended dataset for the Texas power grid under review for release. The team hopes these contributions will enhance access to power system data and remove barriers to the development of new computational techniques for power systems, particularly, those inspired by cognitive sciences.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Current and Future Performance of the CMS Simulation

The CMS full simulation using Geant4 has delivered billions of simulated events for analysis during Runs 1 and 2 of the LHC. However, the HL-LHC dataset will be an order of magnitude larger, with a similar increase in occupancy per event. In addition, the upgraded CMS detector will be considerably more complex, with an extended silicon tracker and a high granularity calorimeter in the endcap region. Increases in conventional computing resources are subject to both technological and budgetary limitations, so novel approaches are needed to improve software efficiency and to take advantage of new architectures and heterogeneous resources. Several projects are in development to address these needs, including the vectorized geometry library Vec-Geom and the GeantV transport engine, which uses track-level parallelization. The current computing performance of the CMS simulation will be presented as a baseline, along with an overview of the various optimizations already available for Geant4. Finally, the progress and outlook for integrating VecGeom and GeantV in the CMS software framework will be discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Towards Exascale Scientific Metadata Management

Advances in technology and computing hardware are enabling scientists from all areas of science to produce massive amounts of data using large-scale simulations or observational facilities. In this era of data deluge, effective coordination between the data production and the analysis phases hinges on the availability of metadata that describe the scientific datasets. Existing workflow engines have been capturing a limited form of metadata to provide provenance information about the identity and lineage of the data. However, much of the data produced by simulations, experiments, and analyses still need to be annotated manually in an ad hoc manner by domain scientists. Systematic and transparent acquisition of rich metadata becomes a crucial prerequisite to sustain and accelerate the pace of scientific innovation. Yet, ubiquitous and domain-agnostic metadata management infrastructure that can meet the demands of extreme-scale science is notable by its absence. To address this gap in scientific data management research and practice, we present our vision for an integrated approach that (1) automatically captures and manipulates information-rich metadata while the data is being produced or analyzed and (2) stores metadata within each dataset to permeate metadata-oblivious processes and to query metadata through established and standardized data access interfaces. We motivate the need for the proposed integrated approach using applications from plasma physics, climate modeling and neuroscience, and then discuss research challenges and possible solutions.

97 MATHEMATICS AND COMPUTING↗

Susceptibility of the Colorado River basin to megadroughts in a warming climate. Final Report.

During the course of the project, we completed the following: Task 1: Off-line model simulations. We proposed to extend Maurer’s data sets through at least 2005 (to include extreme drought years in the Colorado basin). Task 2: Coupled model implementation. We proposed to implement the “standard” climate version of WRF, as used by collaborator Ruby Leung in NARCCAP simulations (see Section 5.4), and perform tests to assure that model output for runs equivalent to NARCCAP Phase 1 (reanalysis boundary conditions) are consistent. We proposed do test sensitivity to higher spatial resolution. Task 3: Retrospective coupled model simulations. We proposed to compare the climatology of key land surface and atmospheric variables from the retrospective RCM and GCM runs, as well as 1/8-degree as the historical North American Land Data Assimilation System dataset (Maurer et al., 2002) to evaluate the ability of each of the RCMs and GCMs to reproduce characteristics of observed climate. Task 4: Future climate simulations. We proposed to follow the approach from Task 3 in analysis of a 21st Century simulation, which was to be performed using the coupled WRF-VIC simulations over the Colorado River basin. Task 5: Diagnostic studies, including partially coupled model runs. We proposed to perform further diagnostic studies of the nature of land hydrologic feedbacks as they are predicted to affect runoff response to climate change in the Colorado River basin. Task 6: Evaluation of conditions in post-2000 drought years.

99 GENERAL AND MISCELLANEOUS↗

Land Model Testbed: Accelerating Development, Benchmarking and Analysis of Land Surface Models

A Land Model Testbed (LMT), designed to provide a computational framework for systematically assessing model fidelity and supporting rapid development of complex multiscale models, offers a general-purpose workflow for conducting large ensemble simulations of multiple land surface models, post-processing large volumes of model output, and evaluating model results. It leverages existing tools for launching model simulations and the International Land Model Benchmarking (ILAMB) package for assessing model fidelity through comparison with best-available observational datasets. Increased complexity and proliferation of uncertain parameters in process representations in land surface models has driven the need for frequent and intensive testing and evaluating of models to quantify uncertainties and optimize parameters such that results are consistent with observations. The LMT described here meets these needs by providing tools to run thousands of ensemble simulations simultaneously and post-process their output files, by automating execution of an enhanced version of ILAMB with site-specific benchmarks and multivariate functional relationships, and by offering ensemble diagnostics and a customizable dashboard for displaying model performance metrics and associated graphics. We envision the LMT capabilities will serve as a foundational computational resource for a proposed user facility focused on terrestrial multiscale model--data integration.

Sreepathi, Sarat↗

Comparing simulated demand flexibility against actual performance in commercial office buildings

Commercial building energy benchmarking has been used as a mechanism to evaluate energy use of a single building over time, relative to other similar buildings, or to simulations of a reference building conforming to various energy standards. Lack of empirical demand flexibility data and consistent flexibility metrics has limited the ability to compare demand flexibility performance with estimated demand flexibility in buildings. In this study, we collected demand response performance data for a total of 831 demand response events from 192 sites as a first step to build such a building demand flexibility dataset, and propose a standard core data schema to consolidate field data from different sources. We also performed parametric simulations of a control strategy called “global temperature adjustment” using commercial office prototype building models. We then compared the simulated demand flexibility performance against the actual data for offices with global temperature adjustment strategy implemented. During demand response events with an average outside air temperature of 34 °C (range 23 °C–42 °C), the measured demand decrease intensity of the demand flexibility metrics were 6.1 watts per square meter (W/m 2 ), 10.0 W/m 2 , 11.1 W/m 2 , 7.1 W/m 2 , and 4.7 W/m 2 for small, small–medium, medium, medium–large, and large office buildings, respectively. Compared to the measured data in medium- and large-size buildings, the simulated demand decrease intensity was 0.7 W/m 2 (17%) lower on average. The discrepancy between simulated and measured peak demand intensities fell within one standard deviation of the mean measured data. Here, the comparison results validate the credibility of simulations in capturing real building data for assessing the technical potential of building demand flexibility.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

An Assessment of the Laminar Hypersonic Double-Cone Experiments in the LENS-XX Tunnel

This is an investigation on two experimental datasets of laminar hypersonic flows, over a double-cone geometry, acquired in Calspan—University at Buffalo Research Center’s Large Energy National Shock (LENS)-XX expansion tunnel. These datasets have yet to be modeled accurately. A previous paper suggested that this could partly be due to mis-specified inlet conditions. The authors of this paper solved a Bayesian inverse problem to infer the inlet conditions of the LENS-XX test section and found that in one case they lay outside the uncertainty bounds specified in the experimental dataset. However, the inference was performed using approximate surrogate models. Here in this paper, the experimental datasets are revisited and inversions for the tunnel test-section inlet conditions are performed with a Navier–Stokes simulator. The inversion is deterministic and can provide uncertainty bounds on the inlet conditions under a Gaussian assumption. It was found that deterministic inversion yields inlet conditions that do not agree with what was stated in the experiments. An a posteriori method is also presented to check the validity of the Gaussian assumption for the posterior distribution. This paper contributes to ongoing work on the assessment of datasets from challenging experiments conducted in extreme environments, where the experimental apparatus is pushed to the margins of its design and performance envelopes.

42 ENGINEERING↗

A Global Model Simulation of Aerosol Effects of Surface Radiation Budget- Toward Understanding of the "Dimming to Brightening" Transition

We present a global model study on the role aerosols play in the change of solar radiation at Earth's surface that transitioned from a decreasing (dimming) trend to an increasing (brightening) trend. Our primary objective is to understand the relationship between the long-term trends of aerosol emission, atmospheric burden, and surface solar radiation. More specifically, we use the recently compiled comprehensive global emission datasets of aerosols and precursors from fuel combustion, biomass burning, volcanic eruptions and other sources from 1980 to 2006 to simulate long-term variations of aerosol distributions and optical properties, and then calculate the multi-decadal changes of short-wave radiative fluxes at the surface and at the top of the atmosphere by coupling the GOCART model simulated aerosols with the Goddard radiative transfer model. The model results are compared with long-term observational records from ground-based networks and satellite data. We will address the following critical questions: To what extent can the observed surface solar radiation trends, known as the transition from dimming to brightening, be explained by the changes of anthropogenic and natural aerosol loading on global and regional scales? What are the relative contributions of local emission and long-range transport to the surface radiation budget and how do these contributions change with time?

Chin, Mian↗

SU-E-T-528: Robustness Evaluation for Fiducial-Based Accelerated Partial Breast Proton Therapy

Purpose: To investigate the robustness of the proton treatment plans in the presence of rotational setup error when patient is aligned with implanted fiducials. Methods: Five Stage I invasive breast cancer patients treated with the APBP protocol (PCG BRE007-12) were studied. The rotational setup errors were simulated by rotating the original CT images around the body center clockwise and counterclockwise 5 degrees (5CW and 5CCW). Manual translational registration was then performed to match the implanted fiducials on the rotated images to the original dataset. Patient contours were copied to the newly created CT set. The original treatment plan was applied to the new CT dataset with the beam isocenter placed at the geometrical center of PTV. The dose distribution was recalculated for dosimetric parameters comparison. Results: CTV and PTV (D95 and V95) coverages were not significantly different between the two simulated plans (5CW and 5CCW) and the original plan. PTV D95 and CTV D95 absolute difference among the three plans were relatively small, with maximum changes of 0.28 CGE and 0.15 CGE, respectively. PTV V95 and CTV V95 absolute differences were 0.79% and 0.48%. The dosage to the thyroid, heart, contralateral breast and lung remained zero for all three plans. The Dmax and Dmean to the volume of ipsilateral breast excluding CTV were compared, with maximum difference values of 1.02 CGE for Dmax and 3.56 CGE for Dmean. Ipsilateral lung Dmean maintained no significant changes through the three plan comparison, with the largest value 0.32 CGE. Ipsilateral lung Dmax was the most sensitive parameter to this simulation study, with a maximum difference at 20.2 CGE. Conclusion: Our study suggests that fiducial-based Accelerated Partial Breast Proton Therapy is robust with respect to +/− 5 degree patient setup rotational errors, as long as the internal fiducial markers are used for patient alignment.

60 APPLIED LIFE SCIENCES↗

SNAD transient miner: Finding missed transient events in ZTF DR4 using k-D trees

Here we report the automatic detection of 11 transients (7 possible supernovae and 4 active galactic nuclei candidates) within the Zwicky Transient Facility fourth data release (ZTF DR4), all of them observed in 2018 and absent from public catalogs. Among these, three were not part of the ZTF alert stream. Our transient mining strategy employs 41 physically motivated features extracted from both real light curves and four simulated light curve models (SN Ia, SN II, TDE, SLSN-I). These features are input to a k-D tree algorithm, from which we calculate the 15 nearest neighbors. After pre-processing and selection cuts, our dataset contained approximately a million objects among which we visually inspected the 105 closest neighbors from seven of our brightest, most well-sampled simulations, comprising 89 unique ZTF DR4 sources. Our result illustrates the potential of coherently incorporating domain knowledge and automatic learning algorithms, which is one of the guiding principles directing the SNAD team. It also demonstrates that the ZTF DR is a suitable testing ground for data mining algorithms aiming to prepare for the next generation of astronomical data.

79 ASTRONOMY AND ASTROPHYSICS↗

Neural network emulation of flow in heavy-ion collisions at intermediate energies

Applications of new techniques in machine learning are speeding up progress in research in various fields. In this work, we construct and evaluate a deep neural network (DNN) to be used within a Bayesian statistical framework as a faster and more reliable alternative to the Gaussian process (GP) emulator of an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model simulator of heavy-ion reactions at intermediate beam energies. We found strong evidence of the DNN being able to emulate the IBUU simulator's prediction on the strengths of protons' directed and elliptical flow very efficiently even with small training datasets and with accuracy about ten times higher than the GP. Here, limitations of our present work and future improvements are also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Variable Grid Method

NETL’s Variable Grid Method (VGM) is a novel approach that leverages GIS capabilities to simultaneously visualize and quantify spatial data trends and underlying data uncertainty. The flexible VGM approach utilizes a range of spatial datasets and uncertainty quantifications, which can be calculated using data related to sample density, sample variance, interpolation error, multiple simulations, etc., to create an integrated visualization of data and uncertainty. The intuitive manner of the VGM helps communicate the relationship between uncertainty and spatial data to effectively guide research, support advanced computation analyses, and helps inform research, management and policy decisions. The VGM approach is being developed into a standalone tool designed to integrate into a variety of research workflows, providing users the capabilities to utilize the VGM for their datasets, analysis, and models in support of their decision making needs. The VGM approach was developed as part of NETL’s Offshore portfolio. It was further matured and utilized in projects associated with both NETL’s Unconventional Resources and Carbon Storage portfolios, and is now being further developed by Varigrid, LLC.

Carbon Storage↗