Disentangling redshift-space distortions and non-linear bias using the 2D power spectrum
Explore the source record for details and available documents.
Jennings, Elise↗
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.
Bowtie filters for dedicated breast CT: Theory and computational implementation
Purpose: To design bowtie filters with improved properties for dedicated breast CT to improve image quality and reduce dose to the patient. Methods: The authors present three different bowtie filters designed for a cylindrical 14-cm diameter phantom with a uniform composition of 40/60 breast tissue, which vary in their design objectives and performance improvements. Bowtie design #1 is based on single material spectral matching and produces nearly uniform spectral shape for radiation incident upon the detector. Bowtie design #2 uses the idea of basis material decomposition to produce the same spectral shape and intensity at the detector, using two different materials. Bowtie design #3 eliminates the beam hardening effect in the reconstructed image by adjusting the bowtie filter thickness so that the effective attenuation coefficient for every ray is the same. All three designs are obtained using analytical computational methods and linear attenuation coefficients. Thus, the designs do not take into account the effects of scatter. The authors considered this to be a reasonable approach to the filter design problem since the use of Monte Carlo methods would have been computationally intensive. The filter profiles for a cone-angle of 0° were used for the entire length of each filter because the differences between those profiles and the correct cone-beam profiles for the cone angles in our system are very small, and the constant profiles allowed construction of the filters with the facilities available to us. For evaluation of the filters, we used Monte Carlo simulation techniques and the full cone-beam geometry. Images were generated with and without each bowtie filter to analyze the effect on dose distribution, noise uniformity, and contrast-to-noise ratio (CNR) homogeneity. Line profiles through the reconstructed images generated from the simulated projection images were also used as validation for the filter designs. Results: Examples of the three designs are presented. Initial verification of performance of the designs was done using analytical computations of HVL, intensity, and effective attenuation coefficient behind the phantom as a function of fan-angle with a cone-angle of 0°. The performance of the designs depends only weakly on incident spectrum and tissue composition. For all designs, the dynamic range requirement on the detector was reduced compared to the no-bowtie-filter case. Further verification of the filter designs was achieved through analysis of reconstructed images from simulations. Simulation data also showed that the use of our bowtie filters can reduce peripheral dose to the breast by 61% and provide uniform noise and CNR distributions. The bowtie filter design concepts validated in this work were then used to create a computational realization of a 3D anthropomorphic bowtie filter capable of achieving a constant effective attenuation coefficient behind the entire field-of-view of an anthropomorphic breast phantom. Conclusions: Three different bowtie filter designs that vary in performance improvements were described and evaluated using computational and simulation techniques. Results indicate that the designs are robust against variations in breast diameter, breast composition, and tube voltage, and that the use of these filters can reduce patient dose and improve image quality compared to the no-bowtie-filter case.
60 APPLIED LIFE SCIENCES↗
Errors introduced in fission neutron spectrum measurements using a single reference
Measurements of physical quantities are often made relative to a reference quantity, simplifying the interpretation of the measurement and its uncertainties. However, if applied without careful considerations, this technique can carry with it a significant and complex systematic error that is usually ignored. For example, measurements of continuous neutron spectra relative to a reference spectrum yield a result that does not properly include the important effects of neutron scattering in the experimental environment. These effects act to distort the measured spectrum shape and create notable errors in spectrum-integrated quantities. While this effect is one of a series of known potential sources of systematic error including, but not limited to target impurities, angular distributions, incident beam flux, and others, the error introduced from differences in environmental scattering effects between the reference and desired spectra are frequently overlooked. In this work, we demonstrate the origin of this error using simulated data corresponding to measurements of the neutron-induced prompt fission neutron spectra of 233 U, 235 U, 238 U, and 239 Pu as well as the spontaneous fission neutron spectrum of 252 Cf. Errors in the spectrum shape can be as high as 10%–15% depending on the combination of incident and outgoing neutron energy. The spectrum integral and average spectrum energy also contain errors that change with incident neutron energy and can be as high as ~2%. Although this study is based on simulated data, we also show a method for proving that this error exists in an experimental data set and we show a correction that can be applied to results of this type to obtain a more accurate result, thereby largely removing systematic errors from environmental scattering effects. The potential for these kinds of errors should be considered in measurements producing continuous neutron spectra.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Domain Adaptation for Measurements of Strong Gravitational Lenses
Upcoming surveys are predicted to discover galaxy-scale strong lenses on the magnitude of 105, making deep learning methods necessary in lensing data analysis. Currently, there is insufficient real lensing data to train deep learning algorithms, but training only on simulated data results in poor performance on real data. Domain adaptation can bridge the gap between simulated and real datasets. We adopt domain adaptation on the estimation of Einstein radius in simulated galaxy-scale gravitational lensing images. We evaluate two domain adaptation techniques - domain adversarial neural networks (DANN) and maximum mean discrepancy (MMD). We train on a source domain of simulated lenses and apply it to a target domain with emulation of DES survey conditions. We show that both domain adaptation techniques can significantly improve the model performance on the more complex target domain datasets. Our results show the potential of using domain adaptation to perform analysis on future survey data with a deep neural network trained on simulated data.
79 ASTRONOMY AND ASTROPHYSICS↗
Discovering faint and high apparent motion rate near-Earth asteroids using a deep learning program
ABSTRACT Although many near-Earth objects have been found by ground-based telescopes, some fast-moving ones, especially those near detection limits, have been missed by observatories. We developed a convolutional neural network for detecting faint fast-moving near-Earth objects. It was trained with artificial streaks generated from simulations and was able to find these asteroid streaks with an accuracy of 98.7 per cent and a false positive rate of 0.02 per cent on simulated data. This program was used to search image data from the Zwicky Transient Facility (ZTF) in four nights in 2019, and it identified six previously undiscovered asteroids. The visual magnitudes of our detections range from ∼19.0 to 20.3 and motion rates range from ∼6.8 to 24 deg d−1, which is very faint compared to other ZTF detections moving at similar motion rates. Our asteroids are also ∼1–51 m diameter in size and ∼5–60 lunar distances away at close approach, assuming their albedo values follow the albedo distribution function of known asteroids. The use of a purely simulated data set to train our model enables the program to gain sensitivity in detecting faint and fast-moving objects while still being able to recover nearly all discoveries made by previously designed neural networks which used real detections to train neural networks. Our approach can be adopted by any observatory for detecting fast-moving asteroid streaks.
Wang, Franklin (ORCID:0000000344887355)↗
High-Accuracy Module Emulators from Physically-Constrained AI Algorithms
Focal Area(s): How do we use AI tools to integrate observations, simulated data and physical and chemical fundamentals (Focal Area 3) into model components (Focal Area 2) that have high accuracy and stability and low computational burden to improve Earth System Predictability? Science Challenge: Earth system modeling of the hydrological cycle involves compute-intensive modules representing complex chemical and physical process. Recently, AI tools that are far less compute intensive have been developed that emulate these modules, but many of these efforts are not yet sufficiently accurate or even stable. We know a lot about the physics and chemistry of earth system processes. The Science Challenge is developing AI tools that not only incorporate observations and simulated data, but also incorporate the physics and chemistry of the process, while still maintaining the compute efficiency.
54 ENVIRONMENTAL SCIENCES↗
Scalability Testing Approach for Internet of Things for Manufacturing SQL and NoSQL Database Latency and Throughput
The proliferation of low-cost sensors and industrial data solutions has continued to push the frontier of manufacturing technology. Machine learning and other advanced statistical techniques stand to provide tremendous advantages in production capabilities, optimization, monitoring, and efficiency. The tremendous volume of data gathered continues to grow, and the methods for storing the data are critical underpinnings for advancing manufacturing technology. This work aims to investigate the ramifications and design tradeoffs within a decoupled architecture of two prominent database management systems (DBMS): sql and NoSQL. A representative comparison is carried out with Amazon Web Services (AWS) DynamoDB and AWS Aurora MySQL. The technologies and accompanying design constraints are investigated, and a side-by-side comparison is carried out through high-fidelity industrial data simulated load tests using metrics from a major US manufacturer. The results support the use of simulated client load testing for comparing the latency of database management systems as a system scales up from the prototype stage into production. As a result of complex query support, MySQL is favored for higher-order insights, while NoSQL can reduce system latency for known access patterns at the expense of integrated query flexibility. Here, by reviewing this work, a manufacturer can observe that the use of high-fidelity load testing can reveal tradeoffs in IoTfM write/ingestion performance in terms of latency that are not observable through prototype-scale testing of commercially available cloud DB solutions.
AWS↗
Race conditions and data partitioning: risks posed by common errors to reproducible parallel simulations
When parallel algorithms for simulation were introduced in the 1970s, their development and use interested only experts in parallel computation. This circumstance changed as multi-core processors became commonplace, putting a parallel computer into the hands of every modeler. A natural outcome is growing interest in parallel simulation among persons not intimately familiar with parallel computing. At the same time, parallel simulation tools continue to be developed with the implicit assumption that the modeler is knowledgeable about parallel programming. The unintended consequence is a rapidly growing number of users of parallel simulation tools that are unlikely to recognize when the interaction of race conditions, partitioning strategies, and simultaneous action in their simulation models make results non-reproducible, thereby calling into question the validity of conclusions drawn from the simulation data. Here, we illustrate the potential dangers of exposing parallel algorithms to users who are not experts in parallel computation with example models constructed using existing parallel simulation tools. By doing so, we hope to refocus tool developers on usability, even if this new focus incurs loss of some performance.
97 MATHEMATICS AND COMPUTING↗
Curated Modeled Fault Data Set
The curated fault model simulation data set consists of tagged and fully-described time series representing simulated faults for the AFDD test building (ORNLs Flexible Research Platform (FRP)), including baseline performance, faulty performance, and corresponding energy impact. A total of 26 fault models are considered for 99 simulation scenarios with various fault intensity levels. Additional Contacts: Principal investigator: Matt Leach Matt.Leach@nrel.gov Simulation performer: Janghyun Kim Janghyun.Kim@nrel.gov
AFDD↗
Structure sensitivity in oxide catalysis: First-principles kinetic Monte Carlo simulations for CO oxidation at RuO{sub 2}(111)
We present a density-functional theory based kinetic Monte Carlo study of CO oxidation at the (111) facet of RuO{sub 2}. We compare the detailed insight into elementary processes, steady-state surface coverages, and catalytic activity to equivalent published simulation data for the frequently studied RuO{sub 2}(110) facet. Qualitative differences are identified in virtually every aspect ranging from binding energetics over lateral interactions to the interplay of elementary processes at the different active sites. Nevertheless, particularly at technologically relevant elevated temperatures, near-ambient pressures and near-stoichiometric feeds both facets exhibit almost identical catalytic activity. These findings challenge the traditional definition of structure sensitivity based on macroscopically observable turnover frequencies and prompt scrutiny of the applicability of structure sensitivity classifications developed for metals to oxide catalysis.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Understanding Complex Magnetic Spin Textures with Simulation-Assisted Lorentz Transmission Electron Microscopy
There is an increased interest in topologically nontrivial magnetic spin textures such as skyrmions and chiral domain-wall solitons, both from a point of fundamental physics understanding as well as potential technological interest in low-power memory applications. In order to control their behavior, it is necessary to understand their complex spin texture at the nanoscale. Lorentz transmission electron microscopy (LTEM) is a suitable technique for studying these systems due to its high spatial resolution and capability to simultaneously characterize magnetic texture and microstructure. In this work, we present the application of PyLorentz, an open-source software suite that we have developed, for quantitative image analysis of Neel-type skyrmions in thin-film heterostructures. PyLorentz enhances LTEM capabilities by enabling reconstruction of magnetic induction maps from experimental images, as well as simulating LTEM images using micromagnetic simulation data. We demonstrate this for simulated Neel skyrmions as well as experimental data from [Pt/Co/W] multilayer heterostructures. Finally, we also show how simulation-assisted LTEM analysis is crucial for understanding these complex magnetic spin textures, in which the reconstructed magnetic induction map (seen in the LTEM images) differs significantly from the magnetization configuration.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Online Static Load Model Estimation in Distribution Systems: Preprint
The paper focuses on parameter identification for time varying voltage dependent polynomial load models, namely the constant impedance (Z), constant current (I), and constant power (P), in the distribution systems. Typically, parameter identification has been done for aggregated loads at substations for the transmission planning studies. Given the higher voltage volatility faced by active distribution systems, this paper takes the load model estimation to the distribution networks. The paper assumes the presence of measurements at the distribution nodes from various sources, each at varying temporal resolutions. Least squares estimation (LSE) method with bounded variables is used to convert the measurements into ZIP parameter estimations for loads at various locations of the feeder. The results and discussions focus on some of the estimation issues faced in a distribution system due to varying voltage sensitivities of loads along a feeder, measurement resolution, impact of window size and sampling rate of measurement data, and the use of dynamic simulations to create synthetic measurement data in the absence of real ones.
24 POWER TRANSMISSION AND DISTRIBUTION↗
Development of New g-function Data for Simulating a Novel Shallow Bore Ground Heat Exchanger
Geothermal heat pump (GHP), which is also referred as ground source heat pump, is the most energy-efficient technology for space heating and cooling. However, the application of GHP is hindered by its high initial cost, of which approximately 30% is for the ground heat exchanger (GHE). In recent years, researchers have developed different types of shallow bore ground heat exchangers (SBGHE) intending to reduce the cost of drilling. The ability to predict the thermal response of an SBGHE is critical for sizing SBGHE. While g-functions have been commonly used for predicting the thermal response of conventional vertical bore ground heat exchangers (VBGHE), they cannot be directly used for predicting the performance of SBGHE because they did not account for the impacts of the seasonal variation of the soil temperature along with the depth of an SBGHE. In addition, an SBGHE has a larger thermal mass within the borehole than the VBGHE due to the larger borehole diameter of SBGHE. This study develops new g-function data pairs for a new design of SBGHE, which is named Underground Thermal Battery (UTB). Impacts of the seasonal variation of soil temperature on the thermal response of UTB were accounted for by superposing a time-dependent soil temperature onto the g-functions calculated with a numerical model that assumes constant undisturbed soil temperature. The TOUGH program was used to predict the thermal response of several configurations involving multiple UTBs. The results indicate that the proposed methodology is appropriate to generate g-functions for the UTB, and the g-function value of UTB is much lower than that of VBGHE in the time range of 15 min to 1 year due to the large thermal mass and convection heat transfer within the UTB.
Shi, Liang↗
Design Considerations of a Coordinative Demand Charge Mitigation Strategy
This paper presents a coordinative demand charge mitigation (DCM) strategy for reducing electricity consumption during system peak periods. Available DCM resources include batteries, diesel generators, controllable loads, and conservation voltage reduction. All resources are directly controlled by load serving entities. A mixed integer linear programming based energy management algorithm is developed to optimally coordinate of DCM resources considering the load payback effect. To better capture system peak periods, two different kinds of load forecast are used: the day-ahead load forecast and the peak-hour probability forecast. Five DCM strategies are compared for reconciling the discrepancy between the two forecasting results. The DCM strategies are tested using actual utility data. Simulation results show that the proposed algorithm can effectively mitigate the demand charge while preventing the system peak from being shifted to the payback hours. We also identify the diminishing return effect, which can help load serving entities optimize the size of their DCM resources.
Hu, Rongxing↗
Detecting electron neutrinos from solar dark matter annihilation by JUNO
We explore the electron neutrino signals from light dark matter (DM) annihilation in the Sun for the large liquid scintillator detector JUNO. In terms of the spectrum features of three typical DM annihilation channels χχ→νν-bar,τ{sup +}τ{sup −},bb-bar, we take two sets of selection conditions to calculate the expected signals and atmospheric neutrino backgrounds based on the Monte Carlo simulation data. Then the JUNO sensitivities to the spin independent DM-nucleon and spin dependent DM-proton cross sections are presented. It is found that the JUNO projected sensitivities are much better than the current spin dependent direct detection experimental limits for the νν-bar and τ{sup +}τ{sup −} channels. In the spin independent case, the JUNO will give the better sensitivity to the DM-nucleon cross section than the LUX and CDMSlite limits for the νν-bar channel with the DM mass lighter than 6.5 GeV. If the νν-bar or τ{sup +}τ{sup −} channel is dominant, the future JUNO results are very helpful for us to understand the tension between the DAMA annual modulation signal and other direct detection exclusions.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Production of alternate realizations of DESI fiber assignment for unbiased clustering measurement in data and simulations
A critical requirement of spectroscopic large scale structure analyses is correcting for selection of which galaxies to observe from an isotropic target list. This selection is often limited by the hardware used to perform the survey which will impose angular constraints of simultaneously observable targets, requiring multiple passes to observe all of them. In SDSS this manifested solely as the collision of physical fibers and plugs placed in plates. In DESI, there is the additional constraint of the robotic positioner which controls each fiber being limited to a finite patrol radius. A number of approximate methods have previously been proposed to correct the galaxy clustering statistics for these effects, but these generally fail on small scales. To accurately correct the clustering we need to upweight pairs of galaxies based on the inverse probability that those pairs would be observed (Bianchi & Percival 2017). This paper details an implementation of that method to correct the Dark Energy Spectroscopic Instrument (DESI) survey for incompleteness. To calculate the required probabilities, we need a set of alternate realizations of DESI where we vary the relative priority of otherwise identical targets. These realizations take the form of alternate Merged Target Ledgers (AMTL), the files that link DESI observations and targets. We present the method used to generate these alternate realizations and how they are tracked forward in time using the real observational record and hardware status, propagating the survey as though the alternate orderings had been adopted. We detail the first applications of this method to the DESI One-Percent Survey (SV3) and the DESI year 1 data. We include evaluations of the pipeline outputs, estimation of survey completeness from this and other methods, and validation of the method using mock galaxy catalogs.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Towards a unified description of the hydrogen bond network of liquid water: A dynamics based approach
The definition of a hydrogen bond (H-bond) is intimately related to the topological and dynamic properties of the hydrogen bond network within liquid water. The development of a universal H-bond definition for water is an active area of research as it would remove many ambiguities in the network properties that derive from the fixed definition employed to assign whether a water dimer is hydrogen bonded. This work investigates the impact that an electronic-structure based definition, an energetic, and a geometric definition of the H-bond has upon both topological and dynamic network behavior of simulated water. In each definition, the use of a cutoff (either geometric or energetic) to assign the presence of a H-bond leads to the formation of transiently bonded or broken dimers, which have been quantified within the simulation data. The relative concentration of transient species, and their duration, results in two of the three definitions sharing similarities in either topological or dynamic features (H-bond distribution, H-bond lifetime, etc.), however no two definitions exhibit similar behavior for both classes of network properties. In fact, two networks with similar local network topology (as indicated by similar average H-bonds) can have dramatically different global network topology (as indicated by the defect state distributions) and altered H-bond lifetimes. A dynamics based correction scheme is then used to remove artificially transient H-bonds and to repair artificially broken bonds within the network such that the corrected network exhibits the same structural and dynamic properties for two H-bond definitions (the properties of the third definition being significantly improved). The algorithm described represents a significant step forward in the development of a unified hydrogen bond network whose properties are independent of the original hydrogen bond definition that is employed.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Learning effective physical laws for generating cosmological hydrodynamics with Lagrangian deep learning
The goal of generative models is to learn the intricate relations between the data to create new simulated data, but current approaches fail in very high dimensions. When the true data-generating process is based on physical processes, these impose symmetries and constraints, and the generative model can be created by learning an effective description of the underlying physics, which enables scaling of the generative model to very high dimensions. In this work, we propose Lagrangian deep learning (LDL) for this purpose, applying it to learn outputs of cosmological hydrodynamical simulations. The model uses layers of Lagrangian displacements of particles describing the observables to learn the effective physical laws. The displacements are modeled as the gradient of an effective potential, which explicitly satisfies the translational and rotational invariance. The total number of learned parameters is only of order 10, and they can be viewed as effective theory parameters. Further, we combine N-body solver fast particle mesh (FastPM) with LDL and apply it to a wide range of cosmological outputs, from the dark matter to the stellar maps, gas density, and temperature. The computational cost of LDL is nearly four orders of magnitude lower than that of the full hydrodynamical simulations, yet it outperforms them at the same resolution. We achieve this with only of order 10 layers from the initial conditions to the final output, in contrast to typical cosmological simulations with thousands of time steps. This opens up the possibility of analyzing cosmological observations entirely within this framework, without the need for large dark-matter simulations.
79 ASTRONOMY AND ASTROPHYSICS↗