Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulated 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 1,405 records · Page 78

Generative models for simulation of KamLAND-Zen

Abstract The next generation of searches for neutrinoless double beta decay ($$0 \nu \beta \beta $$ 0 ν β β ) are poised to answer deep questions on the nature of neutrinos and the source of the Universe’s matter–antimatter asymmetry. They will be looking for event rates of less than one event per ton of instrumented isotope per year. To claim discovery, accurate and efficient simulations of detector events that mimic$$0 \nu \beta \beta $$ 0 ν β β is critical. Traditional Monte Carlo (MC) simulations can be supplemented by machine-learning-based generative models. This work describes the performance of generative models that we designed for monolithic liquid scintillator detectors like KamLAND to produce accurate simulation data without a predefined physics model. We present their current ability to recover low-level features and perform interpolation. In the future, the results of these generative models can be used to improve event classification and background rejection by providing high-quality abundant generated data.

Physics↗

Citizen science for IceCube: Name that Neutrino

Name that Neutrino is a citizen science project where volunteers aid in classification of events for the IceCube Neutrino Observatory, an immense particle detector at the geographic South Pole. From March 2023 to September 2023, volunteers did classifications of videos produced from simulated data of both neutrino signal and background interactions. Name that Neutrino obtained more than 128,000 classifications by over 1800 registered volunteers that were compared to results obtained by a deep neural network machine-learning algorithm. Possible improvements for both Name that Neutrino and the deep neural network are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Characterization and identification of HPC applications at leadership computing facility

High Performance Computing (HPC) is an important method for scientific discovery via large-scale simulation, data analysis, or artificial intelligence. Leadership-class supercomputers are expensive, but essential to run large HPC applications. The Petascale era of supercomputers began in 2008, with the first machines achieving performance in excess of one petaflops, and with the advent of new supercomputers in 2021 (e.g., Aurora, Frontier), the Exascale era will soon begin. However, the high theoretical computing capability (i.e., peak FLOPS) of a machine is not the only meaningful target when designing a supercomputer, as the resources demand of applications varies. A deep understanding of the characterization of applications that run on a leadership supercomputer is one of the most important ways for planning its design, development and operation. In order to improve our understanding of HPC applications, user demands and resource usage characteristics, we perform correlative analysis of various logs for different subsystems of a leadership supercomputer. This analysis reveals surprising, sometimes counter-intuitive patterns, which, in some cases, conflicts with existing assumptions, and have important implications for future system designs as well as supercomputer operations. For example, our analysis shows that while the applications spend significant time on MPI, most applications spend very little time on file I/O. Combined analysis of hardware event logs and task failure logs show that the probability of a hardware FATAL event causing task failure is low. Combined analysis of control system logs and file I/O logs reveals that pure POSIX I/O is used more widely than higher level parallel I/O. Based on holistic insights of the application gained through combined and co-analysis of multiple logs from different perspectives and general intuition, we engineer features to "fingerprint" HPC applications. We use t-SNE (a machine learning technique for dimensionality reduction) to validate the explainability of our features and finally train machine learning models to identify HPC applications or group those with similar characteristic. To the best of our knowledge, this is the first work that combines logs on file I/O, computing, and inter-node communication for insightful analysis of HPC applications in production.

Liu, Zhengchun↗

Axom

The Axom library consists of n collection of software components that provide core computer science infrastructure capabilities that can be shared by diverse high performance computing (HPC) applications. The current set of capabilities that Axom provides includes: Customizable support for error/warning and diagnostic message reporting and coordination among components of integrated applications (e.g., physics packages, libraries, etc.); Scalable parallel aggregation and filtering of message logging; In memory datastore for hierarchical, mesh-aware simulation data. The datastore supports data description, allocation, deallocation. parallel l/0, etc; Mesh data model abstractions that enable the development of mesh-based computational algorithms that work with many diverse, application· specific mesh types; and Computational geometry primitives for developing high performance parallel algorithms based on spatial geometry and meshes, e.g., spatial acceleration data structures, sets, maps, relations, etc.

Hornung, RichardD.↗

geothermal_osr (Energy Predictor for Geothermal Open Source Reservoir) [SWR-22-23]

This is a repository with the tools developed for predicting energy produced at Open-Source Reservoir (OSR). It includes both simulation data for OSR, as well as Jupyter notebooks for training and evaluating prediction models. OSR was constructed based on the data from Brady Hot Springs reservoir (Nevada, USA) but has a number of sufficiently modified characteristics and does not disclose any sensitive data.

Duplyakin, Dmitry↗

eagles-project/asediag

Aerosol process analysis in model-native Spectral Element (SE) grid. Aerosol Diagnostics on Model Native Grid is a Python-based tool designed for diagnosing aerosol processes in E3SM. It is particularly designed to analyze simulation data on the model-native spectral element (SE) grid, soit is also known as Aerosol SE Diagnostics or “asediag”

Hassan, Taufiq↗

NOODLES Grid [SWR-25-93]

NOODLES Grid is a real-time visualization server for power system simulation data. It loads precomputed datasets, builds optimized instance renderings, and streams live interactive scenes to connected clients. It is built for high scalability, flexible visualization, and fast interaction. NOODLES is a cross-platform/device/tool protocol for collaborative visualization. NOODLES was Developed at the National Renewable Energy Laboratory (NREL) as a capability of the Insight Center https://www.nrel.gov/computational-science/insight-center.html

Brunhart-Lupo, Nicholas [National Renewable Energy↗

Comparing Regional Energy Consumption for Direct Drone and Truck Deliveries

Drone delivery, once thought of as fictitious, is becoming a reality with the efforts of both forward-looking enterprises and supportive government policies. This emerging mode of e-commerce delivery raises many concerns. One important concern is the energy efficiency of direct delivery drones compared with conventional delivery trucks at a regional systems level. Here, in this study, we develop and apply methods to quantify the regional energy impacts of drone delivery, then we assess these impacts and compare them with the impacts of truck delivery. To study this problem, we develop an optimization model that determines an optimal set of fulfillment centers (FCs) with variable service capacities that allow drones to make direct e-commerce deliveries. We adopt two drone delivery energy estimation models from the literature and use them as inputs to demonstrate the potential range of energy needs. We also develop another optimization model to account for the energy consumption of diesel trucks (DTs) and battery electric vehicles (BEVs). We test the models using validated simulation data for the Chicago metropolitan area in the U.S. to quantify the energy implications of these three delivery modes. For drone delivery, we further extend our analyses by considering the impact of wind speed and flight patterns. Our results show that direct delivery drones require 15.8% more energy than BEVs on an average windy day, and they need 15% more energy than DTs on a very windy day. We provide essential parameter values for reproducibility and list relevant open problems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Quantifying shifts in natural selection on codon usage between protein regions: a population genetics approach

Codon usage bias (CUB), the non-uniform usage of synonymous codons, occurs across all domains of life. Adaptive CUB is hypothesized to result from various selective pressures, including selection for efficient ribosome elongation, accurate translation, mRNA secondary structure, and/or protein folding. Given the critical link between protein folding and protein function, numerous studies have analyzed the relationship between codon usage and protein structure. The results from these studies have often been contradictory, likely reflecting the differing methods used for measuring codon usage and the failure to appropriately control for confounding factors, such as differences in amino acid usage between protein structures and changes in the frequency of different structures with gene expression. Here we take an explicit population genetics approach to quantify codon-specific shifts in natural selection related to protein structure in S. cerevisiae and E. coli. Unlike other metrics of codon usage, our approach explicitly separates the effects of natural selection, scaled by gene expression, and mutation bias while naturally accounting for a region’s amino acid usage. Bayesian model comparisons suggest selection on codon usage varies only slightly between helix, sheet, and coil secondary structures and, similarly, between structured and intrinsically-disordered regions. Similarly, in contrast to previous findings, we find selection on codon usage only varies slightly at the termini of helices in E. coli. Using simulated data, we show this previous work indicating “non-optimal” codons are enriched at the beginning of helices in S. cerevisiae was due to failure to control for various confounding factors (e.g. amino acid biases, gene expression, etc.), and rather than selection to modulate cotranslational folding. Our results reveal a weak relationship between codon usage and protein structure, indicating that differences in selection on codon usage between structures are slight. In addition to the magnitude of differences in selection between protein structures being slight, the observed shifts appear to be idiosyncratic and largely codon-specific rather than systematic reversals in the nature of selection. Overall, our work demonstrates the statistical power and benefits of studying selective shifts on codon usage or other genomic features from an explicitly evolutionary approach. Limitations of this approach and future potential research avenues are discussed.

59 BASIC BIOLOGICAL SCIENCES↗

Estimating Sparse Direct Effects in Multivariate Regression With the Spike-and-Slab LASSO

The multivariate regression interpretation of the Gaussian chain graph model simultaneously parametrizes (i) the direct effects of p predictors on q outcomes and (ii) the residual partial covariances between pairs of outcomes. We introduce a new method for fitting sparse versions of these models with spike-and-slab LASSO (SSL) priors. We develop an Expectation Conditional Maximization algorithm to obtain sparse estimates of the p × q matrix of direct effects and the q × q residual precision matrix. Our algorithm iteratively solves a sequence of penalized maximum likelihood problems with self-adaptive penalties that gradually filter out negligible regression coefficients and partial covariances. Because it adaptively penalizes individual model parameters, our method is seen to outperform fixed-penalty competitors on simulated data. We establish the posterior contraction rate for our model, buttressing our method’s excellent empirical performance with strong theoretical guarantees. Using our method, we estimated the direct effects of diet and residence type on the composition of the gut microbiome of elderly adults.

EM algorithm↗

Fluorogenic aptamers resolve the flexibility of RNA junctions using orientation-dependent FRET

To further understand the transcriptome, new tools capable of measuring folding, interactions, and localization of RNA are needed. Although Förster resonance energy transfer (FRET) is an angle- and distance-dependent phenomenon, the majority of FRET measurements have been used to report distances, by assuming rotationally averaged donor–acceptor pairs. Angle-dependent FRET measurements have proven challenging for nucleic acids due to the difficulties in incorporating fluorophores rigidly into local substructures in a biocompatible manner. Fluorescence turn-on RNA aptamers are genetically encodable tags that appear to rigidly confine their cognate fluorophores, and thus have the potential to report angular-resolved FRET. In this work, we use the fluorescent aptamers Broccoli and Mango-III as donor and acceptor, respectively, to measure the angular dependence of FRET. Joining the two fluorescent aptamers by a helix of variable length allowed systematic rotation of the acceptor fluorophore relative to the donor. FRET oscillated in a sinusoidal manner as a function of helix length, consistent with simulated data generated from models of oriented fluorophores separated by an inflexible helix. Analysis of the orientation dependence of FRET allowed us to demonstrate structural rigidification of the NiCo riboswitch upon transition metal-ion binding. This application of fluorescence turn-on aptamers opens the way to improved structural interpretation of ensemble and single-molecule FRET measurements of RNA.

59 BASIC BIOLOGICAL SCIENCES↗

Investigations of nonlinear polarization transfer between obliquely intersecting beams

Electromagnetic waves propagating through plasma can interact nonlinearly through a variety of different mechanisms. The excitation of a plasma beat wave (ions or electrons) can create a refractive index modulation that changes the dispersion of the interacting beams. Alternatively, high-intensity beams can enter the regime where relativistic nonlinearities influence the propagation dynamics. In recent studies [ Opt. Express 29 , 1162 ( 2021 ) OPEXFF 1094-4087 10.1364/OE.413064 ], it was proposed that two beams propagating along the same axis can exchange their polarization state due to nonlinear interaction. Here we present a numerical analysis of two laser beams intersecting in a nonlinear medium at varying angles. Polarization transfer is observed as predicted by analytical theory for a range of angles. For small angles, it is found that filamentation of the interacting beams becomes important. Analytical estimates of the filamentation threshold are presented, and good agreement is found with the simulation data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Extracting information from X-ray diffraction patterns containing Laue oscillations

The presence of Laue oscillations in a film grown on a solid surface is broadly taken as indicating a high quality, crystallographically aligned film of the targeted compound. In this paper we briefly review the origins of both Laue oscillations and Kiessig fringes and show how they can be used together to determine if extra thickness exists above or below the coherently diffracting domains. The differences between experimental and “ideal” films are discussed and the effect of structural features (roughness, different thickness coherently diffracting domains and thickness in addition to the coherently diffracting domains) are illustrated with experimental and simulated data for metal and mixed-metal chalcogenide films of titanium, bismuth, vanadium/iron, and bismuth/molybdenum. Examples are given showing how quantitative information can be extracted from experimental diffraction patterns.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Detectability of Varied Hybridization Scenarios Using Genome-Scale Hybrid Detection Methods

Hybridization events complicate the accurate reconstruction of phylogenies, as they lead to patterns of genetic heritability that are unexpected under traditional, bifurcating models of species trees. This phenomenon has led to the development of methods to infer these varied hybridization events, both methods that reconstruct networks directly, as well as summary methods that predict individual hybridization events from a subset of taxa. However, a lack of empirical comparisons between methods – especially those pertaining to large networks with varied hybridization scenarios – hinders their practical use. Here, we provide a comprehensive review of popular summary methods: TICR, MSCquartets, HyDe, Patterson’s D-Statistic (ABBA-BABA), D3, and Dp. TICR and MSCquartets are based on quartet concordance factors gathered from gene tree topologies and HyDe, Patterson’s D-Statistic, D3, and Dp use site pattern frequencies to identify hybridization events between sets of three taxa. We then use simulated data to address questions of method accuracy and ideal use scenarios by testing methods against complex networks which depict gene flow events that differ in depth (timing), quantity (single vs. multiple, overlapping hybridizations), and rate of gene flow (γ). We find that deeper or multiple hybridization events may introduce noise and weaken the signal of hybridization, leading to higher relative false negative rates across all methods. Despite some forms of hybridization eluding quartet-based detection methods, MSCquartets displays high precision in most scenarios. While HyDe results in high false negative rates when tested on hybridizations involving extinct or unsampled ghost lineages, HyDe is the only method able to identify the direction of hybridization, distinguishing the source parental lineages from recipient hybrid lineages. Lastly, we test the methods on a dataset of ultraconserved elements from the bee subfamily Nomiinae, finding possible hybridization events between clades which correspond to regions of poor support in the species tree estimated in a previous study.

Bjorner, Marianne B.↗

ESTIMATION OF THE WAKEFIELD RESONANT FREQUENCY USING DIFFERENT SIMULATION TOOLS

Reliable and fast wakefield calculations are important for beam dynamics and THz generation in accelerators. We compare approaches and limitations with different available simulation tools and compare results. As an example, we analyse a cylindrical corrugated waveguide with narrower and wider radii of 5 and 6 mm, and aperture width and periodicity of 1 and 2 mm using ABCI by KEK and ANSYS HFSS simulation software. First, ABCI was used to determine the resonant frequency of the cylindri-cal corrugated waveguide with different lengths. The corresponding results were taken as a reference to exam-ine the simulation method in the ANSYS HFSS environ-ment. The eigenmode solver of ANSYS HFSS was used to determine different resonant frequencies. This was fol-lowed by examining a cylindrical corrugated waveguide that was implemented having a length of 17 mm for dif-ferent conductor settings. The corresponding waveguide was fed by a plane wave having a resonant frequency that satisfies the intersection of the light line and dispersion curve of the corresponding mode. The analysis showed that the maximum loss factor was achieved at 39.20 and 40.07 GHz using ANSYS HFSS simulation data for differ-ent conductors, whereas ABCI resulted in 38.50 GHz. The talk will present the results gathered by different simula-tion setups implemented in ANSYS HFSS.

Unutmaz, Muhammed [SLAC National Accelerator Labor↗

ATOM - Accelerating therapeutics through opportunities in medicine

Create a new paradigm of drug discovery that would reduce the time from an identified drug target to clinical candidate from the current ~6 years to just 12 months. ATOM will develop, test, and validate a multidisciplinary approach to drug discovery in which modern science, technology and engineering, supercomputing, simulations, data science, and artificial intelligence are highly integrated into a single drug-discovery platform that can ultimately be shared with the drug development community at-large.

60 APPLIED LIFE SCIENCES↗

The MELCOR Plot File Format

This work provides a detailed explanation of the MELCOR Plot File format. MELCOR is a fully integrated, engineering-level computer code that models the progression of severe accidents in light water reactor nuclear power plants, and it primarily exports its simulation data to a binary file in the Plot File format. This work documents the Plot File's overall structure and the data types of exported information. The process of interpreting and associating time series information with MELCOR's Plot Variables for complete time histories is also presented. The format and meaning of time-independent information, called Special information, is also discussed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗