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 379 records · Page 21

Facilitating Staging-based Unstructured Mesh Processing to Support Hybrid In-Situ Workflows

In-situ and in-transit processing alleviate the gap between the computing and I/O capabilities by scheduling data analytics close to the data source. Hybrid in-situ processing splits data analytics into two stages: the data processing that runs in-situ aims to extract regions of interest, which are then transferred to staging services for further in-transit analytics. To facilitate this type of hybrid in-situ processing, the data staging service needs to support complex intermediate data representations generated/consumed by the in-situ tasks. Unstructured (or irregular) mesh is one such derived data representation that is typically used and bridges simulation data and analytics. However, how staging services efficiently support unstructured mesh transfer and processing remains to be explored. This paper investigates design options for transferring and processing unstructured mesh data using staging services. Using polygonal mesh data as an example, we show that hybrid in-situ workflows with staging-based unstructured mesh processing can effectively support hybrid in-situ workflows, and can significantly decrease data movement overheads.

data-driven↗

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

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

97 MATHEMATICS AND COMPUTING↗

HERO WEC V1.0 - WEC-Sim Detailed Simulation Runs and Summary Data

Data sets from simulation runs include 144 wave cases that were run based on the WEC-Sim model (https://mhkdr.openei.org/submissions/483). The 144 wave cases represent waves with the following wave height and wave period ranges: - Significant Wave Height: 0.25 - 4.0m in 0.25m increments - Wave Period: 5 - 13 sec in 1 sec increments Each run was simulated using a Pierson-Moskowitz irregular wave spectrum with a 200 second ramp time and a total simulation time of 2,200 seconds. Summary data set includes a spreadsheet and image files with matrices that are associated with data from simulation runs. All matrices cover the same significant wave height and wave periods from the sim runs, in the same increments. The following matrices are included: - Mechanical Power Input: Absorbed power calculated using linear input velocity at PTO (power take-off) multiplied by winch tension. - Mechanical Power Output: The hydraulic power calculated at the pump output prior to the pressure relief valve dump. - Feed Pressure: Average pressure at RO (reverse osmosis) input prior to pre-filters - RO Power Input: The hydraulic power calculated at the RO system inlet, accounts for flow loss due to pressure relief valve - Water Production: Average clean water flow at RO output

16 TIDAL AND WAVE POWER↗

Expanded analysis of machine learning models for nuclear transient identification using TPOT

Industries around the world are becoming more and more data driven. The nuclear field is no exception with several different applications being proposed. One popular area of research is the use of machine learning in transient detection. This paper seeks to build upon a previous study which made use of the AutoML package TPOT to train traditional machine learning models to classify transient events occurring with a reactor. Synthetic data was once again collected using a GPWR reactor simulator. Data on 12 different events was collected using 15 different initial conditions. Here, a dataset consisting of over 100,000 data points was compiled and used to train 7 different machine learning models using a pre-defined TPOT dictionary with 12 different preprocessing techniques. Three of the trained models were able to produce validation results in the 90s with the expanded dataset. Once the models were trained, it was possible to look into where during the simulation, misclassifications occurred. Using these three models, analysis was done to determine if TPOT could be used to train models that were effective if important features were missing. The results from this were positive with the newly trained models scoring close to the original models. Finally, to conclude this study, the three high performing models were retrained using different random states to see if there was any major variation when different states were used.

42 ENGINEERING↗

Detecting vortices in superconductors: Extracting one-dimensional topological singularities from a discretized complex scalar field

In type II superconductors, the dynamics of superconducting vortices determine their transport properties. In the Ginzburg-Landau theory, vortices correspond to topological defects in the complex order parameter. Extracting their precise positions and motion from discretized numerical simulation data is an important, but challenging, task. In the past, vortices have mostly been detected by analyzing the magnitude of the complex scalar field representing the order parameter and visualized by corresponding contour plots and isosurfaces. However, these methods, primarily used for small-scale simulations, blur the fine details of the vortices, scale poorly to large-scale simulations, and do not easily enable isolating and tracking individual vortices. In this paper, we present a method for exactly finding the vortex core lines from a complex order parameter field. With this method, vortices can be easily described at a resolution even finer than the mesh itself. The precise determination of the vortex cores allows the interplay of the vortices inside a model superconductor to be visualized in higher resolution than has previously been possible. Finally, by representing the field as the set of vortices, this method also massively reduces the data footprint of the simulations and provides the data structures for further analysis and feature tracking.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Recasting a model atomistic glassformer as a system of icosahedra

We consider a binary Lennard-Jones glassformer whose super-Arrhenius dynamics are correlated with the formation of icosahedral structures. Upon cooling, these icosahedra organize into mesoclusters. We recast this glassformer as an effective system of icosahedra which we describe with a population dynamics model. This model we parameterize with data from the temperature regime accessible to molecular dynamics simulations. We then use the model to determine the population of icosahedra in mesoclusters at arbitrary temperature. Using simulation data to incorporate dynamics into the model, we predict relaxation behavior at temperatures inaccessible to conventional approaches. Our model predicts super-Arrhenius dynamics whose relaxation time remains finite for non-zero temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Liquid-state polaron theory of the hydrated electron revisited

The quantum path integral/classical liquid-state theory of Chandler and co-workers, created to describe an excess electron in solvent, is re-examined for the hydrated electron. The portion that models electron-water density correlations is replaced by two equations: the range optimized random phase approximation (RO-RPA), and the Donley, Rajasekaran, and Liu (DRL) approximation to the “two-chain” equation, both shown previously to describe accurately the static structure and thermodynamics of strongly charged polyelectrolyte solutions. The static equilibrium properties of the hydrated electron are analyzed using five different electron-water pseudopotentials. The theory is then compared with data from mixed quantum/classical Monte Carlo and molecular dynamics simulations using these same pseudopotentials. It is found that the predictions of the RO-RPA and DRL-based polaron theories are similar and improve upon previous theory, with values for almost all properties analyzed in reasonable quantitative agreement with the available simulation data. Also, it is found using the Larsen, Glover, and Schwartz pseudopotential that the theories give values for the solvation free energy that are at least three times larger than that from experiment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Validation of a Global Geospace Model With a Systems Science Approach Based on Canonical Correlation Analysis

A systems science approach based on canonical correlation analysis (CCA) is applied as a new, behavioral way to validate global geospace models. The biggest novelty of the technique is that it validates models at a system level, whereby a side‐by‐side comparison is performed of CCA applied to a 30‐day observational and the corresponding simulation data sets comprising quiet, moderate and active times. The simulation used the Multiscale Atmosphere‐Geospace Environment (MAGE) model. It is shown that (a) CCA must be combined with sensitivity analysis to be effective, (b) the MAGE model generally reproduces the observed behavior (more so for quieter time intervals), quantified by the intercorrelations between different variables and (c) the technique identifies the SuperMAG SML index as a quantity for which refinements of the model are needed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact and mitigation of angular uncertainties in Bragg coherent x-ray diffraction imaging

Bragg coherent diffraction imaging (BCDI) is a powerful technique to explore the local strain state and morphology of microscale crystals. The method can potentially reach nanometer-scale spatial resolution thanks to the advances in synchrotron design that dramatically increase coherent flux. However, there are experimental bottlenecks that may limit the image reconstruction quality from future high signal-to-noise ratio measurements. In this work we show that angular uncertainty of the sample orientation with respect to a fixed incoming beam is one example of such a factor, and we present a method to mitigate the resulting artifacts. On the basis of an alternative formulation of the forward problem, we design a phase retrieval algorithm which enables the simultaneous reconstruction of the object and determination of the exact angular position corresponding to each diffraction pattern in the data set. We have tested the algorithm performance on simulated data for different degrees of angular uncertainty and signal-to-noise ratio.

36 MATERIALS SCIENCE↗

Modification of wind turbine wakes by large-scale, convective atmospheric boundary layer structures

In this study, we consider the impact of large-scale, convective structures in an unstable atmospheric boundary layer on wind turbine wakes. Simulation data from a high-fidelity large-eddy simulation (LES) of the AWAKEN wind farm site matching unstable atmospheric conditions were analyzed, and both turbine performance and wake behavior were affected based on their location relative to the convective structures. Turbines located in updraft regions of the flow experienced lower inflow velocity and generated less power, but their wakes were observed to recover faster and saw greater turbulent kinetic energy mixing higher in the boundary layer. The opposite effect was found for turbines in the downdraft regions of the convective structures. A simplified model of this wake behavior was also developed based on a two-dimensional k–ε Reynolds-Averaged Navier–Stokes formulation. This simplified model included the effects of vertical transport, but could be efficiently solved as a parabolic system, and was found to capture similar wake modifications observed in the high-fidelity LES computations.

16 TIDAL AND WAVE POWER↗

Square Kilometre Array Science Data Challenge 1: analysis and results

ABSTRACT As the largest radio telescope in the world, the Square Kilometre Array (SKA) will lead the next generation of radio astronomy. The feats of engineering required to construct the telescope array will be matched only by the techniques developed to exploit the rich scientific value of the data. To drive forward the development of efficient and accurate analysis methods, we are designing a series of data challenges that will provide the scientific community with high-quality data sets for testing and evaluating new techniques. In this paper, we present a description and results from the first such Science Data Challenge 1 (SDC1). Based on SKA MID continuum simulated observations and covering three frequencies (560, 1400, and 9200 MHz) at three depths (8, 100, and 1000 h), SDC1 asked participants to apply source detection, characterization, and classification methods to simulated data. The challenge opened in 2018 November, with nine teams submitting results by the deadline of 2019 April. In this work, we analyse the results for eight of those teams, showcasing the variety of approaches that can be successfully used to find, characterize, and classify sources in a deep, crowded field. The results also demonstrate the importance of building domain knowledge and expertise on this kind of analysis to obtain the best performance. As high-resolution observations begin revealing the true complexity of the sky, one of the outstanding challenges emerging from this analysis is the ability to deal with highly resolved and complex sources as effectively as the unresolved source population.

Bonaldi, A.↗

Online learning of quadratic manifolds from streaming data for nonlinear dimensionality reduction and nonlinear model reduction

Here, this work introduces an online greedy method for constructing quadratic manifolds from streaming data, designed to enable in situ analysis of numerical simulation data on the Petabyte scale. Unlike traditional batch methods, which require all data to be available upfront and take multiple passes over the data, the proposed online greedy method incrementally updates quadratic manifolds in one pass as data points are received, eliminating the need for expensive disk input/output operations as well as storing and loading data points once they have been processed. A range of numerical examples demonstrate that the online greedy method learns accurate quadratic manifold embeddings while being capable of processing data that far exceed common disk input/output capabilities and volumes as well as main-memory sizes.

97 MATHEMATICS AND COMPUTING↗

Fast radio burst dispersion measure distribution as a probe of helium reionization

Fast radio burst (FRB) discoveries are occurring rapidly, with thousands expected from upcoming surveys. The dispersion measures (DM) observed for FRB include important information on cosmological distances and the ionization state of the universe from the redshift of emission until today. Rather than considering the DM–redshift relation, we investigate the statistical ensemble of the distribution of dispersion measures. Here, we explore the use of this abundance information, with and without redshift information, to probe helium reionization through simulated data to redshift z = 6. Carrying out Monte Carlo simulations of FRB survey samples, we examine the effect of different source redshift distributions, host galaxy models, sudden vs gradual reionization, and covariance with cosmological parameters on determination of helium reionization properties. We find that a fluence limited survey with 10 4 FRBs can discriminate different helium reionization histories at ~6σ using the DM-distribution of bursts, without redshift information (and ~10σ with redshifts).

79 ASTRONOMY AND ASTROPHYSICS↗

A Framework for Compressing Unstructured Scientific Data via Serialization

We present a general framework for compressing unstructured scientific data with known local connectivity. A common application is simulation data defined on arbitrary finite element meshes. The framework employs a greedy topology preserving reordering of original nodes which allows for seamless integration into existing data processing pipelines. This reordering process depends solely on mesh connectivity and can be performed offline for optimal efficiency. However, the algorithm’s greedy nature also supports on-the-fly implementation. The proposed method is compatible with any compression algorithm that leverages spatial correlations within the data. The effectiveness of this approach is demonstrated on a large-scale real dataset using several compression methods, including MGARD, SZ, and ZFP.

Reshniak, Viktor [ORNL] (ORCID:0000000315454462)↗

Scalable Hybrid Learning Techniques for Scientific Data Compression

Data compression is becoming critical for storing scientific data because many scientific applications need to store large amounts of data and post process this data for scientific discovery. Unlike image and video compression algorithms that limit errors to primary data (PD), scientists require compression techniques that accurately preserve derived quantities of interest (QoIs). Here, this article presents a physics-informed compression technique implemented as an end-to-end, scalable, GPU-based pipeline for data compression that addresses this requirement. Our hybrid compression technique combines machine learning techniques and standard compression methods. Specifically, we combine an autoencoder, an error-bounded lossy compressor to provide guarantees on raw data error, and a constraint satisfaction post-processing step to preserve the QoIs within a minimal error (generally less than floating point error). The effectiveness of the data compression pipeline is demonstrated by compressing nuclear fusion simulation data generated by a large-scale fusion code, XGC, which produces hundreds of terabytes of data in a single day. Our approach works within the ADIOS framework and results in compression by a factor of more than 150 while requiring only a few percent of the computational resources necessary for generating the data, making the overall approach highly effective for practical scenarios.

ITER↗

Monitoring an Electro-Refining Process (Rev. 1)

In this report we examine measurements from two different sensors for detecting a diversion of plutonium in an electro-refining (ER) process. One method is based on counts from a high dose neutron detector (HDND). The other is based on spectra from a micro-calorimeter sensor. The spectra come from ER-salt, fission product waste, and from a source term. The goal is to differentiate between normal and off normal operating conditions using neutron counts or changes in counts in the spectra. Unfortunately, there is no electro-refining process that we can actually monitor to obtain real data from various operating conditions so we will use simulated data to illustrate our approach.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗