Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulator data”

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 127 records · Page 7

Comparative Visualization of Vector Field Ensembles Based on Longest Common Subsequence

We propose a longest common subsequence (LCS) based approach to compute the distance among vector field ensembles. By measuring how many common blocks the ensemble pathlines passing through, the LCS distance defines the similarity among vector field ensembles by counting the number of sharing domain data blocks. Compared to the traditional methods (e.g. point-wise Euclidean distance or dynamic time warping distance), the proposed approach is robust to outlier, data missing, and sampling rate of pathline timestep. Taking the advantages of smaller and reusable intermediate output, visualization based on the proposed LCS approach revealing temporal trends in the data at low storage cost, and avoiding tracing pathlines repeatedly. Finally, we evaluate our method on both synthetic data and simulation data, which demonstrate the robustness of the proposed approach.

Liu, Richen↗

Implications of gnomonic distortion on electron backscatter diffraction and transmission Kikuchi diffraction

The effect of gnomonic distortion on orientation indexing of electron backscatter diffraction patterns is explored through simulation of electron diffraction patterns for sample-to-detector geometries associated with transmission Kikuchi diffraction (TKD) and electron backscatter diffraction (EBSD). Simulated data were analysed by computing a similarity index for both Hough transformed data and simulated patterns to determine the sensitivity of each method for detecting subtle differences in the effect of gnomonic distortions on electron diffraction patterns. These results indicate that the increased gnomonic distortions in electron diffraction patterns for a TKD geometry enhance the sensitivity for detecting subtle differences in interband angles. Additionally, the utilisation of a Hough transform-based indexing approach further enhances the sensitivity.

36 MATERIALS SCIENCE↗

POD-based constrained sensor placement and field reconstruction from noisy wind measurements: A perturbation study

Sensor placement at the extrema of Proper Orthogonal Decomposition (POD) is efficient and leads to accurate reconstruction of the wind field from a limited number of measure- ments. In this paper we extend this approach of sensor placement and take into account measurement errors and detect possible malfunctioning sensors. We use the 48 hourly spa- tial wind field simulation data sets simulated using the Weather Research an Forecasting (WRF) model applied to the Maine Bay to evaluate the performances of our methods. Specifically, we use an exclusion disk strategy to distribute sensors when the extrema of POD modes are close. It turns out that this strategy can also reduce the error of recon- struction from noise measurements. Also, by a cross-validation technique, we successfully locate the malfunctioning sensors.

54 ENVIRONMENTAL SCIENCES↗

Methane Emissions Produced by Pneumatic Devices and Produced Water Tanks on Natural Gas Wellsites

In recent years, more scrutiny has been placed on the release of greenhouse gasses (GHGs) by the oil and gas industry in North America as climate concerns increase. In the coming years, stricter regulations proposed by the Environmental Protection Agency (EPA) concerning the release of GHGs will further direct attention to GHG emissions quantification, modeling, and overall reduction. Active natural gas wellsites employ production equipment that routinely emits methane during normal operation. Production equipment often includes, but is not limited to, gas production units (GPUs), pneumatic controllers (PCs) and actuators powered by produced natural gas, and produced water storage tanks. This work details direct methane emission measurements taken at active natural gas wellsites originating from intermittent pneumatic controllers and produced water storage tanks. Efforts to model methane flows originating from these sources are also detailed, and those results are compared to direct measurement. The direct methane emissions measurement campaign collected flow data from two active dry wellsites in the Marcellus shale region. These measurement campaigns captured, in total, several weeks of continuous data from these sites. Data included multiple pneumatic device vent flowrate channels, a single channel monitoring total emissions vented from the produced water storage tanks, and corresponding site weather data. Periodic composition measurements were taken from flows vented from produced water tanks. Computational modeling efforts attempted to mimic a given site’s methane emissions behavior given limited input parameters. Direct measurement data and simulation data yielded, through reduction and analysis, apparent methane emissions factors, statistics on pneumatic device actuation events, important relationships between GPU parameters and pneumatic device actuation events, produced water storage tank vent flows and flow compositions, and parameters most contributing to the magnitude and frequency of these fugitive methane emissions. Measurements, considering data from a 72hr period, produced PC emissions factors from 0.008 to 0.03 SCFH per device with an average device emitting 0.026 SCFH at the first site. The second site produced PC emissions factors spanning from 0.05 to 25.47 SCFH with an average device emitting 1.93 SCFH not including the 25.47 SCFH device, which was assumed to be malfunctioning. Two unique water tank emission measurements were made, 618 g CH4 per hour and 32.6 g CH4 per hour. On a water basis, the water tank emissions factors were 0.66 and 0.042 kg CH4 per bbl of water, respectively.

03 NATURAL GAS↗

Simulation process and data flow for a large system dynamics model

This paper documents the workflow and supporting technologies that a large system dynamics model, the biomass scenario model, employs to streamline the data preparation, simulation, quality control, and analysis process at the National Renewable Energy Laboratory. The workflow centers on automation of routine aspects of the flow of data between data stores, simulations, and visualizations. It enforces quality checks on data, reproducibility of computations, and traceability of results, while maintaining complete archives of modeling and analysis artifacts. The resulting frictionless simulation/analysis environment supports large-scale sensitivity analysis, interactive creation of ensembles of simulations, and rapid visualization-based exploration of simulation results.

09 BIOMASS FUELS↗

Efficient scatter distribution estimation and correction in CBCT using concurrent Monte Carlo fitting

Purpose: X-ray scatter is a significant impediment to image quality improvements in cone-beam CT (CBCT). The authors present and demonstrate a novel scatter correction algorithm using a scatter estimation method that simultaneously combines multiple Monte Carlo (MC) CBCT simulations through the use of a concurrently evaluated fitting function, referred to as concurrent MC fitting (CMCF). Methods: The CMCF method uses concurrently run MC CBCT scatter projection simulations that are a subset of the projection angles used in the projection set, P, to be corrected. The scattered photons reaching the detector in each MC simulation are simultaneously aggregated by an algorithm which computes the scatter detector response, S{sub MC}. S{sub MC} is fit to a function, S{sub F}, and if the fit of S{sub F} is within a specified goodness of fit (GOF), the simulations are terminated. The fit, S{sub F}, is then used to interpolate the scatter distribution over all pixel locations for every projection angle in the set P. The CMCF algorithm was tested using a frequency limited sum of sines and cosines as the fitting function on both simulated and measured data. The simulated data consisted of an anthropomorphic head and a pelvis phantom created from CT data, simulated with and without the use of a compensator. The measured data were a pelvis scan of a phantom and patient taken on an Elekta Synergy platform. The simulated data were used to evaluate various GOF metrics as well as determine a suitable fitness value. The simulated data were also used to quantitatively evaluate the image quality improvements provided by the CMCF method. A qualitative analysis was performed on the measured data by comparing the CMCF scatter corrected reconstruction to the original uncorrected and corrected by a constant scatter correction reconstruction, as well as a reconstruction created using a set of projections taken with a small cone angle. Results: Pearson’s correlation, r, proved to be a suitable GOF metric with strong correlation with the actual error of the scatter fit, S{sub F}. Fitting the scatter distribution to a limited sum of sine and cosine functions using a low-pass filtered fast Fourier transform provided a computationally efficient and accurate fit. The CMCF algorithm reduces the number of photon histories required by over four orders of magnitude. The simulated experiments showed that using a compensator reduced the computational time by a factor between 1.5 and 1.75. The scatter estimates for the simulated and measured data were computed between 35–93 s and 114–122 s, respectively, using 16 Intel Xeon cores (3.0 GHz). The CMCF scatter correction improved the contrast-to-noise ratio by 10%–50% and reduced the reconstruction error to under 3% for the simulated phantoms. Conclusions: The novel CMCF algorithm significantly reduces the computation time required to estimate the scatter distribution by reducing the statistical noise in the MC scatter estimate and limiting the number of projection angles that must be simulated. Using the scatter estimate provided by the CMCF algorithm to correct both simulated and real projection data showed improved reconstruction image quality.

60 APPLIED LIFE SCIENCES↗

Managing uncertainty in data-driven simulation-based optimization

Optimization using data from complex simulations has become an attractive decision-making option, due to ability to embed high-fidelity, non-linear understanding of processes within the search for optimal values. Due to lack of tractable algebraic equations, the link between simulations and optimization is oftentimes a surrogate metamodel. However, several forms of uncertainty exist within the cycle that links simulation data, to metamodels, to optimization. Uncertainty may originate from parameters of the simulation, or the form and fitted parameters of the metamodel. This paper reviews different literatures that are relevant to surrogate-based optimization and proposes different strategies for handling uncertainty, by combining machine learning with stochastic programming, robust optimization, and discrepancy modeling. We show that incorporating uncertainty management within simulation-based optimization leads to more robust solutions, which protect the decision-maker from infeasible solutions. Here, we present the results of our proposed approaches through a case study for direct-air capture through temperature swing adsorption.

42 ENGINEERING↗

Constraining the physical conditions in the jets of γ-ray flaring blazars using centimeter-band polarimetry and radiative transfer simulations. I. Data and models for 0420–014, OJ 287, and 1156+295

To investigate parsec-scale jet flow conditions during GeV γ-ray flares detected by the Fermi Large Angle Telescope, we obtained centimeter-band total flux density and linear polarization monitoring observations from 2009.5 through 2012.5 with the 26 m Michigan radio telescope for a sample of core-dominated blazars. We use these data to constrain radiative transfer simulations incorporating propagating shocks oriented at an arbitrary angle to the flow direction in order to set limits on the jet flow and shock parameters during flares temporally associated with γ-ray flares in 0420–014, OJ 287, and 1156+295; these active galactic nuclei exhibited the expected signature of shocks in the linear polarization data. Both the number of shocks comprising an individual radio outburst (3 and 4) and the range of the compression ratios of the individual shocks (0.5-0.8) are similar in all three sources; the shocks are found to be forward-moving with respect to the flow. While simulations incorporating transverse shocks provide good fits for 0420–014 and 1156+295, oblique shocks are required for modeling the OJ 287 outburst, and an unusually low value of the low-energy cutoff of the radiating particles' energy distribution is also identified. Our derived viewing angles and shock speeds are consistent with independent Very Long Baseline Array results. While a random component dominates the jet magnetic field, as evidenced by the low fractional linear polarization, to reproduce the observed spectral character requires that a significant fraction of the magnetic field energy is in an ordered axial component. Both straight and low pitch angle helical field lines are viable scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗

Improving and Assessing the Quality of Uncertainty Quantification in Deep Learning

Deep learning (DL) models have enjoyed increased attention in recent years because of their powerful predictive capabilities. While many successes have been achieved, standard deep learning methods suffer from a lack of uncertainty quantification (UQ). While the development of methods for producing UQ from DL models is an active area of current research, little attention has been given to the quality of the UQ produced by such methods. In order to deploy DL models to high-consequence applications, high-quality UQ is necessary. This report details the research and development conducted as part of a Laboratory Directed Research and Development (LDRD) project at Sandia National Laboratories. The focus of this project is to develop a framework of methods and metrics for the principled assessment of UQ quality in DL models. This report presents an overview of UQ quality assessment in traditional statistical modeling and describes why this approach is difficult to apply in DL contexts. An assessment on relatively simple simulated data is presented to demonstrate that UQ quality can differ greatly between DL models trained on the same data. A method for simulating image data that can then be used for UQ quality assessment is described. A general method for simulating realistic data for the purpose of assessing a model’s UQ quality is also presented. A Bayesian uncertainty framework for understanding uncertainty and existing metrics is described. Research that came out of collaborations with two university partners are discussed along with a software toolkit that is currently being developed to implement the UQ quality assessment framework as well as serve as a general guide to incorporating UQ into DL applications.

97 MATHEMATICS AND COMPUTING↗

Particle Acceleration and Fractional Transport in Turbulent Reconnection

We consider a large-scale environment of turbulent reconnection that is fragmented into a number of randomly distributed unstable current sheets (UCSs), and we statistically analyze the acceleration of particles within this environment. We address two important cases of acceleration mechanisms when particles interact with the UCS: (a) electric field acceleration and (b) acceleration by reflection at contracting islands. Electrons and ions are accelerated very efficiently, attaining an energy distribution of power-law shape with an index 1–2, depending on the acceleration mechanism. The transport coefficients in energy space are estimated from test-particle simulation data, and we show that the classical Fokker–Planck (FP) equation fails to reproduce the simulation results when the transport coefficients are inserted into it and it is solved numerically. The cause for this failure is that the particles perform Levy flights in energy space, while the distributions of the energy increments exhibit power-law tails. We then use the fractional transport equation (FTE) derived by Isliker et al., whose parameters and the order of the fractional derivatives are inferred from the simulation data, and solving the FTE numerically, we show that the FTE successfully reproduces the kinetic energy distribution of the test particles. We discuss in detail the analysis of the simulation data and the criteria that allow one to judge the appropriateness of either an FTE or a classical FP equation as a transport model.

79 ASTRONOMY AND ASTROPHYSICS↗

Macrovoid resolved simulations of transport through HPRO relevant membrane geometries

Modeling the transport properties such as diffusivity and permeability of high pressure reverse osmosis (HPRO) membranes is critical for the selection and manufacture of membranes suitable for operation under high pressure. These properties can be significantly affected by the changes in heterogeneous pore structures due to compaction. Here, the modeling platform presented in this work resolves the two scale porosity in HPRO relevant membranes. A synthetic membrane geometry is constructed based on available experimental visualizations of pore structures. The simulations directly capture the flow channeling that results from a combination of material properties and geometric features of the macrovoids. A parametric study is presented to account for the transition of flow characteristics from a material governed regime to a macrovoid governed regime. Permeability of the membrane is evaluated using the simulation data and compared with an existing model that scales with the square of tortuosity over a range of material properties. The model is found to perform well under a narrow range of tortuosity while deviating from the calculated permeabilities at extreme conditions. The effective permeability of the membrane is found to vary by at least two orders of magnitude between the two flow regimes. It is also observed that tortuosity is a bounded property with its upper limit determined by the macrovoid geometry. Consequently, the tortuosity based correlations fail near a flow regime that is mainly governed by the macrovoids. The modeled permeability can be more than an order of magnitude smaller than the simulation result. A new model based on flux-weighted porosity of a membrane is introduced and its correlation with tortuosity is studied. The model agrees with the simulated data as, in addition to tortuosity, it also accounts for the flux partition within and outside the flow channels. Such correlations enable extending the existing understanding of flow characteristics to enhance predictability of porous media models.

42 ENGINEERING↗

Realized ecological forecast through an interactive Ecological Platform for Assimilating Data (EcoPAD, v1.0) into models

Abstract. Predicting future changes in ecosystem services is not only highly desirable but is also becoming feasible as several forces (e.g., available big data, developed data assimilation (DA) techniques, and advanced cyber-infrastructure) are converging to transform ecological research into quantitative forecasting. To realize ecological forecasting, we have developed an Ecological Platform for Assimilating Data (EcoPAD, v1.0) into models. EcoPAD (v1.0) is a web-based software system that automates data transfer and processing from sensor networks to ecological forecasting through data management, model simulation, data assimilation, forecasting, and visualization. It facilitates interactive data–model integration from which the model is recursively improved through updated data while data are systematically refined under the guidance of model. EcoPAD (v1.0) relies on data from observations, process-oriented models, DA techniques, and the web-based workflow. We applied EcoPAD (v1.0) to the Spruce and Peatland Responses Under Climatic and Environmental change (SPRUCE) experiment in northern Minnesota. The EcoPAD-SPRUCE realizes fully automated data transfer, feeds meteorological data to drive model simulations, assimilates both manually measured and automated sensor data into the Terrestrial ECOsystem (TECO) model, and recursively forecasts the responses of various biophysical and biogeochemical processes to five temperature and two CO 2 treatments in near-real time (weekly). Forecasting with EcoPAD-SPRUCE has revealed that mismatches in forecasting carbon pool dynamics are more related to model (e.g., model structure, parameter, and initial value) than forcing variables, opposite to forecasting flux variables. EcoPAD-SPRUCE quantified acclimations of methane production in response to warming treatments through shifted posterior distributions of the CH 4 :CO 2 ratio and the temperature sensitivity ( Q 10 ) of methane production towards lower values. Different case studies indicated that realistic forecasting of carbon dynamics relies on appropriate model structure, correct parameterization, and accurate external forcing. Moreover, EcoPAD-SPRUCE stimulated active feedbacks between experimenters and modelers to identify model components to be improved and additional measurements to be taken. It has become an interactive model–experiment (ModEx) system and opens a novel avenue for interactive dialogue between modelers and experimenters. Altogether, EcoPAD (v1.0) acts to integrate multiple sources of information and knowledge to best inform ecological forecasting.

54 ENVIRONMENTAL SCIENCES↗

Development of Estimation Method for Radioactivity Concentration Distribution of Radioactive Wastes in a Box-shaped Container - 17034

We are developing a new estimation method suitable for determining radioactivity concentration distribution of radioactive wastes stored in a square box-shaped container. This method uses the Maximum A Posteriori (MAP) method which is a way to solve the inverse problem, and it estimates the radioactivity concentration distribution by comparing simulated data with measured data. In the measurement, gamma rays emitted from the radioactive wastes are measured by several gamma ray spectrometers which are set at different positions in relation to each other. On the other hand, the gamma ray spectrum data of each spectrometer are simulated by a Monte Carlo simulation code. The method uses an iterative calculation to solve the concentration distribution with measured spectrum data and simulated spectrum data. We conducted a simulation test to check the algorithm of our estimation method. And we also conducted an experiment using activated objects. We found that the concentration of the wastes in a box-shaped container could be estimated within a 20 % error for the simulation, and within a 40 % error for the experiment and we confirmed that feasibility of our method was favorable. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Elastic Bayesian Model Calibration

Functional data are ubiquitous in scientific modeling. For instance, quantities of interest are modeled as functions of time, space, energy, density, etc. Uncertainty quantification methods for computer models with functional response have resulted in tools for emulation, sensitivity analysis, and calibration that are widely used. However, many of these tools do not perform well when the computer model’s parameters control both the amplitude variation of the functional output and its alignment (or phase variation). This paper introduces a framework for Bayesian model calibration when the model responses are misaligned functional data. The approach generates two types of data out of the misaligned functional responses: (1) aligned functions so that the amplitude variation is isolated and (2) warping functions that isolate the phase variation. These two types of data are created for the computer simulation data (both of which may be emulated) and the experimental data. The calibration approach uses both types so that it seeks to match both the amplitude and phase of the experimental data. The framework is careful to respect constraints that arise, especially when modeling phase variation, and is framed in a way that it can be done with readily available calibration software. In conclusion, we demonstrate the techniques on two simulated data examples and on two dynamic material science problems: a strength model calibration using flyer plate experiments and an equation of state model calibration using experiments performed on the Sandia National Laboratories’ Z-machine.

97 MATHEMATICS AND COMPUTING↗

Equation-based object-oriented modeling and simulation for data center cooling: A case study

© 2019 Elsevier B.V. Data center cooling accounts for about 1% of electricity usage in the United States. Computer models are pivotal in designing and operating energy-efficient cooling systems. Compared to conventional building performance simulation programs, the equation-based object-oriented modeling language Modelica is an emerging approach that can enable fast prototyping and dynamic simulation of cooling systems. In this case study, we first modeled the cooling and control systems of an actual data center located in Massachusetts using the open-source Modelica Buildings library, and then calibrated a baseline model based on measurement data. The simulation of the baseline model identified several operation-related issues in the cooling and control systems, such as degraded cooling coils, improper dead band in control settings, and simultaneous cooling and heating in air handlers. Afterwards, we used a sequential search technique as well as an optimization scheme to investigate the energy saving potentials for different energy efficiency measures aiming to address the abovementioned issues. Simulation results show potential energy savings up to 24% by resolving identified control-related issues and optimizing the supply air temperature.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NuMI Beam Monitoring Simulation and Data Analysis Status and Progress

With the Main Injector Neutrino Oscillation Search (MINOS) experiment decommissioned, muon and hadron monitors became an important diagnostic tool for the NuMI O↵- axis vμ Appearance (NOvA) experiment at Fermilab to monitor the Neutrinos at the Main Injector (NuMI) beam. The goal of this study is to maintain the quality of the monitor signals and to establish correlations with the neutrino beam profile. And we carry out a systematic study of the response of the muon monitors to the changes in the parameters of the proton beam and lattice parameters. We report here on the progress of the beam data analysis and comparison with the simulation results.

43 PARTICLE ACCELERATORS↗

Machine Learning Potentials with the Iterative Boltzmann Inversion: Training to Experiment

Methodologies for training machine learning potentials (MLPs) with quantum-mechanical simulation data have recently seen tremendous progress. Experimental data have a very different character than simulated data, and most MLP training procedures cannot be easily adapted to incorporate both types of data into the training process. Here, we investigate a training procedure based on iterative Boltzmann inversion that produces a pair potential correction to an existing MLP using equilibrium radial distribution function data. By applying these corrections to an MLP for pure aluminum based on density functional theory, we observe that the resulting model largely addresses previous overstructuring in the melt phase. Interestingly, the corrected MLP also exhibits improved performance in predicting experimental diffusion constants, which are not included in the training procedure. The presented method does not require autodifferentiating through a molecular dynamics solver and does not make assumptions about the MLP architecture. Our results suggest a practical framework for incorporating experimental data into machine learning models to improve the accuracy of molecular dynamics simulations.

36 MATERIALS SCIENCE↗

Combining simulations and data with deep learning and uncertainty quantification for advanced energy modeling

A novel and modern framework for energy modeling is developed in this paper with a focus on nuclear energy modeling and simulation. The framework combines multiphysics simulations and real data, with validation by uncertainty quantification tasks and facilitation by machine and deep learning methods. The hybrid framework is built on the basis of a wide range of physical models, real data, mathematical and statistical methods, and artificial intelligence techniques. The framework is demonstrated in different applications, including quantifying uncertainties in computer simulations, multiphysics coupling, analysis of variance using machine learning surrogate models, deep learning of time series phenomena, and propagating parametric uncertainties of nuclear data. The applications demonstrated are oriented to nuclear engineering simulations, even though majority of the methods are applicable to other energy sources (eg, renewable). Efficient utilization of this framework is expected to yield a much better understanding of the physical phenomena analyzed as well as an improvement in the performance of the energy design/model under construction.

42 ENGINEERING↗