Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation 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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Simulation of convective transport during frequency chirping of a TAE using the MEGA code

Here, we present a procedure to examine energetic particle phase-space during long range frequency chirping phenomena in tokamak plasmas. To apply the proposed method, we have performed self-consistent simulations using the MEGA code and analyzed the simulation data. We demonstrate a traveling wave in phase-space and that there exist specific slices of phase-space on which the resonant particles lie throughout the wave evolution. For non-linear evolution of an n = 6 toroidicity-induced Alfvén eigenmode (TAE), our results reveal the formation of coherent phase-space structures (holes/clumps) after coarse-graining of the distribution function. These structures cause a convective transport in phase-space which implies a radial drift of the resonant particles. We also demonstrate that the rate of frequency chirping increases with the TAE damping rate. Our observations of the TAE behavior and the corresponding phase-space dynamics are consistent with the Berk–Breizman theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Tailoring resonant magnetic perturbation to optimize fast-ion confinement during ELM control in KSTAR

Abstract 3D resonant magnetic perturbation (RMP) is one promising way to control edge localized modes that can cause excessive material erosion of tokamak first walls. However, RMP can lead to undesired degradation of plasma confinement, including fast-particle losses, which can impact the performance and safety of the reactor. This work investigates the optimization of the poloidal spectrum of the 3D field to optimize fast ion confinement during edge localized mode (ELM) suppression. In the initial step, the validity of the modeling framework is tested against experimental data. Simulations successfully replicate an increase in poloidal limiter temperature with different poloidal spectra. Then, the simulation shows improvement of fast ion confinement with a reduction of core resonant response, while edge resonant magnetic fields are maintained above the threshold to sustain the ELM suppression. Reduction of the core resonant fields keeps the Kolmogorov–Arnold–Moser surface and reduces the fast particle losses due to the stochastic magnetic field lines. The results highlight the potential of edge localization of the resonant fields to enhance the performance of fusion reactors, but further investigation is needed to improve the validation of this approach.

3D field↗
DOE OSTI2020

Comparison of Rotor Wake Identification and Characterization Methods for the Analysis of Wake Dynamics and Evolution

Optimal wind power plant design requires understanding of wind turbine wake physics and validation of engineering wake models under wake-controlled operating conditions. In this work, we have developed and investigated several different wake identification and characterization methods for analyzing wake evolution and dynamics. The accuracy and robustness of these methods, based on Gaussian function fitting and adaptive contour identification, have been assessed by application to a large-eddy simulation data set. A new contour-based method based on downstream momentum deficit has been considered. Uncertainties arising from wake-identification errors result in characterizations of the wake expansion, recovery, and meandering motion that differ by 19% of the rotor area, 4% of the freestream, and 15% rotor diameter, respectively.

17 WIND ENERGY↗
DOE OSTI2023

Evaluation of Portable Programming Models to Accelerate LArTPC Detector Simulations

The Liquid Argon Time Projection Chamber (LArTPC) technology is widely used in high energy physics experiments, including the upcoming Deep Underground Neutrino Experiment (DUNE). Accurately simulating LArTPC detector responses is essential for analysis algorithm development and physics model interpretations. Accurate LArTPC detector response simulations are computationally demanding, and can become a bottleneck in the analysis workflow. Compute devices such as General-Purpose Graphics Processing Units (GPGPUs) have the potential to substantially accelerate simulations compared to traditional CPU-only processing. The software development for these compute accelerators often carries the cost of specialized code refactorization and porting to match the target hardware architecture. With the rapid evolution and increased diversity of the computer architecture landscape, it is highly desirable to have a portable solution that also maintains reasonable performance. We report our ongoing effort in evaluating Kokkos as a basis for this portable programming model using LArTPC simulations in the context of the Wire-Cell Toolkit, a C++ library for LArTPC simulations, data analysis, reconstruction and visualization.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Flow annealed importance sampling bootstrap meets differentiable particle physics

High-energy physics requires the generation of large numbers of simulated data samples from complex but analytically tractable distributions called matrix elements. Surrogate models, such as normalizing flows, are gaining popularity for this task due to their computational efficiency. We adopt an approach based on flow annealed importance sampling bootstrap (FAB) that evaluates the differentiable target density during training and helps avoid the costly generation of training data in advance. We show that FAB reaches higher sampling efficiency with fewer target evaluations in high dimensions in comparison to other methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Simulation driven adaptive sampling for neutron-diffraction based strain mapping of additively manufactured parts

Neutron diffraction based strain mapping is a useful technique for measuring residual strains in additively manufactured (AM) metal parts. The measurement is traditionally done by scanning the sample in a point-wise raster pattern to extract the strain at each position. Since the overall scan can span several hours, adaptive sampling approaches using Bayesian optimization based on Gaussian process (BO-GP) regression have been introduced—demonstrating that even with a fraction of the typically made measurements the dominant strain patterns in the sample can be reconstructed. However, the parameters of the BO-GP algorithm have to be carefully chosen for best performance, and the movement time between arbitrary points can offset the time savings from a reduced number of measurement locations. In this paper, we propose algorithms to refine the BO-GP based methods by using simulations of strain patterns in AM parts based on the materials and the process used to print them. We demonstrate that the simulated strain patterns can be used to help choose better parameters for the BO-GP based framework—leading to low reconstruction error for the final strain pattern. Furthermore, we show that the strain mapping experiment can be initialized with a sampling pattern learnt from the simulation data and ordered to reduce movement time, dramatically enabling reduction in the overall time required to run the baseline BO-GP method.

Gaussian process regression↗
DOE OSTI2022

Estimating the local dark matter density in a non-axisymmetric wobbling disc

The density of dark matter near the Sun, ρDM, ⊙, is important for experiments hunting for dark matter particles in the laboratory, and for constraining the local shape of the Milky Way’s dark matter halo. Estimates to date have typically assumed that the Milky Way’s stellar disc is axisymmetric and in a steady-state. Yet the Milky Way disc is neither, exhibiting prominent spiral arms and a bar, and vertical and radial oscillations. Here, we assess the impact of these assumptions on determinations of ρDM, ⊙ by applying a free-form, steady-state, Jeans method to two different N-body simulations of Milky Way-like galaxies. In one, the galaxy has experienced an ancient major merger, similar to the hypothesized Gaia–Sausage–Enceladus; in the other, the galaxy is perturbed more recently by the repeated passage and slow merger of a Sagittarius-like dwarf galaxy. We assess the impact of each of the terms in the Jeans–Poisson equations on our ability to correctly extract ρDM, ⊙ from the simulated data. We find that common approximations employed in the literature – axisymmetry and a locally flat rotation curve – can lead to significant systematic errors of up to a factor ~1.5 in the recovered surface mass density ~2 kpc above the disc plane, implying a fractional error on ρDM, ⊙ of the order of unity. However, once we add in the tilt term and the rotation curve term in our models, we obtain an unbiased estimate of ρDM, ⊙, consistent with the true value within our 95 percent confidence intervals for realistic 20 percent uncertainties on the baryonic surface density of the disc. Other terms – the axial tilt, 2nd Poisson and time-dependent terms – contribute less than 10 percent to ρDM, ⊙ (given current data) and can be safely neglected for now. In the future, as more data become available, these terms will need to be included in the analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A novel approach for calculating galaxy rotation curves using spaxel cross-correlation and iterative smoothing

ABSTRACT Precise measurements of the internal dynamics of galaxies have proven of great importance for understanding the internal dark matter distribution of galaxies. We present a novel method for measuring the line-of-sight (LOS) velocities across the face of galaxies by cross-correlation of spectral pixels (spaxels) and an iterative method of smoothing. On simulated data the method can accurately recover the input LOS velocities for different types of spectra (absorption-line dominated, emission-line dominated, and differing shapes of the continuum), and can handle stellar population radial gradients. Most important of all, it continues to provide reliable measurements of LOS velocities with reasonable uncertainties even when the spectra are very low signal-to-noise ratio (approaching ∼1), which is a challenge for traditional template-fitting approaches. We apply our method to data from a real MaNGA galaxy as a demonstration and find promising results with good precision. This novel approach can be complementary to existing methods primarily based on template fitting.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

KilonovaNet : Surrogate models of kilonova spectra with conditional variational autoencoders

ABSTRACT Detailed radiative transfer simulations of kilonova spectra play an essential role in multimessenger astrophysics. Using the simulation results in parameter inference studies requires building a surrogate model from the simulation outputs to use in algorithms requiring sampling. In this work, we present kilonovanet, an implementation of conditional variational autoencoders (cVAEs) for the construction of surrogate models of kilonova spectra. This method can be trained on spectra directly, removing overhead time of pre-processing spectra, and greatly speeds up parameter inference time. We build surrogate models of three state-of-the-art kilonova simulation data sets and present in-depth surrogate error evaluation methods, which can in general be applied to any surrogate construction method. By creating synthetic photometric observations from the spectral surrogate, we perform parameter inference for the observed light-curve data of GW170817 and compare the results with previous analyses. Given the speed with which kilonovanet performs during parameter inference, it will serve as a useful tool in future gravitational wave observing runs to quickly analyse potential kilonova candidates.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

The miniJPAS survey quasar selection – I. Mock catalogues for classification

In this series of papers, we employ several machine learning (ML) methods to classify the point-like sources from the miniJPAS catalogue, and identify quasar candidates. Since no representative sample of spectroscopically confirmed sources exists at present to train these ML algorithms, we rely on mock catalogues. In this first paper, we develop a pipeline to compute synthetic photometry of quasars, galaxies, and stars using spectra of objects targeted as quasars in the Sloan Digital Sky Survey . To match the same depths and signal-to-noise ratio distributions in all bands expected for miniJPAS point sources in the range 17.5 ≤ r < 24, we augment our sample of available spectra by shifting the original r-band magnitude distributions towards the faint end, ensure that the relative incidence rates of the different objects are distributed according to their respective luminosity functions, and perform a thorough modelling of the noise distribution in each filter, by sampling the flux variance either from Gaussian realizations with given widths, or from combinations of Gaussian functions. Finally, we also add in the mocks the patterns of non-detections which are present in all real observations. Although the mock catalogues presented in this work are a first step towards simulated data sets that match the properties of the miniJPAS observations, these mocks can be adapted to serve the purposes of other photometric surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A better way to define dark matter haloes

ABSTRACT Dark matter haloes have long been recognized as one of the fundamental building blocks of large-scale structure formation models. Despite their importance – or perhaps because of it! – halo definitions continue to evolve towards more physically motivated criteria. Here, we propose a new definition that is physically motivated, effectively unique, and parameter-free: ‘A dark matter halo is comprised of the collection of particles orbiting in their own self-generated potential’. This definition is enabled by the fact that, even with as few as ≈300 particles per halo, nearly every particle in the vicinity of a halo can be uniquely classified as either orbiting or infalling based on its dynamical history. For brevity, we refer to haloes selected in this way as physical haloes. We demonstrate that (1) the mass function of physical haloes is Press–Schechter, provided the critical threshold for collapse is allowed to vary slowly with peak height; and (2) the peak-background split prediction of the clustering amplitude of physical haloes is statistically consistent with the simulation data, with accuracy no worse than ≈5 per cent.

Astronomy & Astrophysics↗
DOE OSTI2025

Validating sequential Monte Carlo for gravitational-wave inference

Nested sampling (NS) is the preferred stochastic sampling algorithm for gravitational-wave inference for compact binary coalescences. It can handle the complex nature of the gravitational-wave likelihood surface and provides an estimate of the Bayesian model evidence. However, there is another class of algorithms that meets the same requirements, but has not been used for gravitational-wave analyses: sequential Monte Carlo (SMC), an extension of importance sampling that maps samples from an initial density to a target density via a series of intermediate densities. In this work, we validate a type of SMC algorithm, called persistent sampling (PS), for gravitational-wave inference. We consider a range of different scenarios including binary black holes and binary neutron stars and real and simulated data and show that PS produces results that are consistent with NS whilst being, on average, 2 times more efficient and 2.74 times faster. This demonstrates that PS is a viable alternative to NS that should be considered for future gravitational-wave analyses.

black hole mergers↗
DOE OSTI2021

Quantum pattern recognition algorithms for charged particle tracking

High-energy physics is facing a daunting computing challenge with the large datasets expected from the upcoming High-Luminosity Large Hadron Collider in the next decade and even more so at future colliders. A key challenge in the reconstruction of events of simulated data and collision data is the pattern recognition algorithms used to determine the trajectories of charged particles. The field of quantum computing shows promise for transformative capabilities and is going through a cycle of rapid development and hence might provide a solution to this challenge. This article reviews current studies of quantum computers for charged particle pattern recognition in high-energy physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Nonparametric extensions of nuclear equations of state: Probing the breakdown scale of relativistic mean-field theory

Phenomenological calculations of the properties of dense matter, such as relativistic mean-field theories, represent a pathway to predicting the microscopic and macroscopic properties of neutron stars. However, such theories do not generically have well-controlled uncertainties and may break down within neutron stars. To faithfully represent the uncertainty in this breakdown scale, we develop a hybrid representation of the dense-matter equation of state, which assumes the form of a relativistic mean-field theory at low densities, while remaining agnostic to any nuclear theory at high densities. To achieve this, we use a nonparametric equation of state model to incorporate the correlations of the underlying relativistic mean-field theory equation of state at low pressures and transition to more flexible correlations above some chosen pressure scale. We perform astrophysical inference under various choices of the transition pressure between the theory-informed and theory-agnostic models. Here, we further study whether the chosen relativistic mean-field theory breaks down above some particular pressure and find no such evidence. Using simulated data for future astrophysical observations at about two-to-three times the precision of current constraints, we show that our method can identify the breakdown pressure associated with a potential strong phase transition.

Equations of state of nuclear matter↗
DOE OSTI2025

Optimal observables for the chiral magnetic effect from machine learning

The detection of the chiral magnetic effect (CME) in relativistic heavy-ion collisions remains challenging due to substantial background contributions that obscure the expected signal. In this Letter, we present a novel machine learning approach for constructing optimized observables that significantly enhance CME detection capabilities. By parametrizing generic observables constructed from flow harmonics and optimizing them to maximize the signal-to-background ratio, we systematically develop CME-sensitive measures that outperform conventional methods. Using simulated data from the anomalous viscous fluid dynamics framework, our machine learning observables demonstrate up to 90% higher sensitivity to CME signals compared to traditional 𝛾 and 𝛿 correlators, while maintaining minimal background contamination. The constructed observables provide physical insight into optimal CME detection strategies and offer a promising path forward for experimental searches of the CME at the BNL Relativistic Heavy Ion Collider and the CERN Large Hadron Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Construction of approximate invariants for nonintegrable Hamiltonian systems

We present a method to construct high-order polynomial approximate invariants (AI) for nonintegrable Hamiltonian dynamical systems and apply it to a modern ring-based particle accelerator. Taking advantage of a special property of one-turn transformation maps expressed as square matrices, AIs can be constructed order by order iteratively. Evaluating AI with simulation data, we observe that AI’s fluctuation is actually a measure of chaos. Through minimizing the fluctuations, the stable region of long-term motions, i.e., the dynamic aperture of the accelerator, could be enlarged.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Autotuning of Double-Dot Devices In Situ with Machine Learning

The current practice of manually tuning quantum dots (QDs) for qubit operation is a relatively time-consuming procedure that is inherently impractical for scaling up and applications. In this work, we report on the in situ implementation of a recently proposed autotuning protocol that combines machine learning (ML) with an optimization routine to navigate the parameter space. In particular, we show that a ML algorithm trained using exclusively simulated data to quantitatively classify the state of a double-QD device can be used to replace human heuristics in the tuning of gate voltages in real devices. Here, we demonstrate active feedback of a functional double-dot device operated at millikelvin temperatures and discuss success rates as a function of the initial conditions and the device performance. Modifications to the training network, fitness function, and optimizer are discussed as a path toward further improvement in the success rate when starting both near and far detuned from the target double-dot range.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Improved Vertical Carrier Transport for Green III-Nitride LEDs Using ( In , Ga ) N Alloy Quantum Barriers

We report on experimental and simulation-based results using ( In , Ga ) N alloy quantum barriers in c -plane green light-emitting diode (LED) structures as a means to improve vertical carrier transport and reduce forward voltage ( V F ) . Three-dimensional device simulations that include random alloy fluctuations are used to understand carrier behavior in a disordered potential. The simulated current density–voltage ( J - V ) characteristics and modified electron-hole overlap | F mod | 2 indicate that increasing the indium fraction in the ( In , Ga ) N quantum barriers leads to a reduced polarization discontinuity at the interface between the quantum barrier and quantum well, thereby reducing V F and improving | F mod | 2 . Maps of electron and hole current through the device show a relatively homogenous distribution in the X Y plane for structures using Ga N quantum barriers; in contrast, preferential pathways for vertical transport are identified in structures with ( In , Ga ) N barriers as regions of high and low current. A positive correlation between hole (electron) current in the p -side ( n -side) barrier and indium fraction reveals that preferential pathways exist in regions of high indium content. Furthermore, a negative correlation between the strain ε z z and indium fraction shows that high indium content regions have reduced strain-induced piezoelectric polarization in the Z direction due to the mechanical constraint of the surrounding lower indium content regions. Experimentally, multiple quantum well green LEDs with ( In , Ga ) N quantum barriers exhibit lower V F and blue-shifted wavelengths relative to LEDs with Ga N quantum barriers, consistent with simulation data. These results can be used to inform heterostructure design of low V F , long-wavelength LEDs and provide important insight into the nature of carrier transport in III-nitride alloy materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗