Search NASA⌕ Search

SEARCH · Search NASA

Results for “astrophysical simulations”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Sugars and Sugar Derivatives in Residues Produced from the UV Irradiation of Astrophysical Ice Analogs

Carbonaceous chondrites contain a large variety of organic compounds of prebiotic interest, which include amino acids, amphiphiles, nucleobases, and sugar derivatives. The presence of these compounds strongly suggests that molecules essential to life can form abiotically under astrophysical conditions. Among the sugar derivatives reported in the Murchison and Murray meteorites, only one sugar (dihydroxyacetone) was found, together with a variety of sugar alcohols and sugar acids containing up to 6 carbon atoms, including sugar acid derivatives of the biological sugars ribose and glucose. On the other hand, laboratory studies on the formation of complex organic molecules from the ultraviolet (UV) irradiation of simulated astrophysical ice mixtures consisting of H2O, CO, CO2, CH3OH, CH4, NH3, etc., at low temperature have been routinely carried out in the past 15 years. These studies have shown that the organic residues recovered at room temperature contain amino acids, amphiphiles, nucleobases, as well as other complex organics, supporting a scenario in which molecules of prebiotic interest can form in extra-terrestrial environments.

Astrochemistry↗
NASA NTRS2019

Enabling Real-time Multi-messenger Astrophysics Discoveries with Deep Learning

Multi-messenger astrophysics is a fast-growing, interdisciplinary field that combines data, which vary in volume and speed of data processing, from many different instruments that probe the Universe using different cosmic messengers: electromagnetic waves, cosmic rays, gravitational waves and neutrinos. In this Expert Recommendation, we review the key challenges of real-time observations of gravitational wave sources and their electromagnetic and astroparticle counterparts, and make a number of recommendations to maximize their potential for scientific discovery. These recommendations refer to the design of scalable and computationally efficient machine learning algorithms; the cyber-infrastructure to numerically simulate astrophysical sources, and to process and interpret multi-messenger astrophysics data; the management of gravitational wave detections to trigger real-time alerts for electromagnetic and astroparticle follow-ups; a vision to harness future developments of machine learning and cyber-infrastructure resources to cope with the big-data requirements; and the need to build a community of experts to realize the goals of multi-messenger astrophysics.

E A Huerta↗
DOE OSTI2023

A Fully Explicit Integrator for Modeling Astrophysical Reactive Flows

Simulating complex astrophysical reacting flows is computationally expensive—reactions are stiff and typically require implicit integration methods. The reaction update is often the most expensive part of a simulation, which motivates the exploration of more economical methods. In this research note, we investigate how the explicit Runge–Kutta–Chebyshev (RKC) method performs compared to an implicit method when applied to astrophysical reactive flows. These integrators are applied to simulations of X-ray bursts arising from unstable thermonuclear burning of accreted fuel on the surface of neutron stars. We show that the RKC method performs with similar accuracy to our traditional implicit integrator, but is more computationally efficient when run on CPUs.

GPU computing↗
DOE OSTI2025

Toward Higher-order Accuracy in Self-gravitating Hydrodynamics

High-order algorithms have emerged in numerical astrophysics as a promising avenue to reduce truncation error (proportional to a power of the linear resolution Δ x ) with only a moderate increase to computational expense. Significant effort has been placed in the development of finite-volume algorithms for (magneto)hydrodynamics; however, state-of-the-art astrophysical simulations tightly couple a plenitude of physics, additionally including gravity, photon transport, cosmic-ray transport, chemistry, and/or diffusion, to name a few. Algorithms frequently operator-split this additional physics (often a first-order error in time) and/or adopt a model wherein their evaluation is limited to second-order accuracy in space. In this work, we present a fourth-order-accurate finite-volume scheme for self-gravitating hydrodynamics on a uniform Cartesian grid. The method supplies source terms for the gravitational acceleration ( ρ g ) and gravitational energy release ( ρ v · g ) associated with fourth-order-accurate solutions to the Poisson equation. Our scheme (1) guarantees the conservation of total linear momentum while (2) decreasing (in proportion to Δ x 4 ) the effects of spurious heating and/or cooling associated with truncation error in the gravity. We demonstrate expected convergence rates for the algorithm by measuring errors in test problems evolving self-gravity modified linear waves and 3D polytropic equilibria. We test robustness of the algorithm by integrating an induced “inside-out” adiabatic collapse. We also discuss a method to smoothly downgrade the solution to second-order spatial accuracy to avoid spurious overshoots near steep density and/or pressure gradients.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Direct cross section measurement of 102 Pd ⁢(𝛾,𝑝) and 102 Pd ⁢(𝛾,𝛼) for the astrophysical 𝑝 process

Background: A handful of neutron-deficient stable nuclei, known as the “p nuclei,” cannot be produced through astrophysical neutron capture processes. Instead, some of these nuclei are proposed to be produced by 𝛾-induced reactions on existing r- and s-process seeds. The specific astrophysical site or sites are not yet identified, however, with uncertainties in the cross sections of these 𝛾-induced reactions playing a role. Databases of reaction rates for astrophysical simulations often rely on theoretical statistical model calculations, such as Hauser-Feshbach, for rates where no experimental information is known. However, reasonable variations in the choice of parametrizations of various nuclear properties can create order-of-magnitude variations in the final predicted cross sections and reaction rates, which are then propagated through the models to the predicted final abundances. Purpose: To better constrain these statistical model calculations and ultimately reduce the uncertainties from the nuclear physics on our understanding of the p nuclei, a measurement of the cross sections of 𝛾-induced reactions on the p-nucleus 102 Pd was undertaken. This work represents the first measurement of its kind, using segmented silicon detectors to measure prompt charged particle emission from 𝛾-induced reactions. Methods: Quasimonoenergetic gamma beams from the High Intensity 𝛾 Source facility bombarded an enriched 102 Pd target. A segmented silicon array was arranged to detect the particles emitted from (𝛾,𝑝) and (𝛾,𝛼) reactions. Results: Reaction cross sections were deduced at multiple 𝛾-beam energies between 10 and 19 MeV, and compared to statistical model calculations using talys-1.96. The 102 Pd ⁢(𝛾,𝑝)⁢ 101 Rh reaction cross section was reasonably well reproduced by a subset of photon strength functions and level densities, though the strength to the ground state of 101 Rh was underestimated at higher incident gamma energies. The 102 Pd ⁢(𝛾,𝛼)⁢ 98 Ru was in general overpredicted by the various alpha-nucleus optical model potentials. Conclusions: While the theoretical cross sections used to model the (𝛾,𝑝) reactions for the p process may be reasonable, a more careful approach is needed in the case of (𝛾,𝛼). Further work to probe gamma-induced reaction cross sections at and near the p nuclei is warranted.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

OpenACC Unified Programming Environment for Multi-hybrid Acceleration with GPU and FPGA

Accelerated computing in HPC such as with GPU, plays a central role in HPC nowadays. However, in some complicated applications with partially different performance behavior is hard to solve with a single type of accelerator where GPU is not the perfect solution in these cases. We are developing a framework and transpiler allowing the users to program the codes with a single notation of OpenACC to be compiled for multi-hybrid accelerators, named MHOAT (Multi-Hybrid OpenACC Translator) for HPC applications. MHOAT parses the original code with directives to identify the target accelerating devices, currently supporting NVIDIA GPU and Intel FPGA, dispatching these specific partial codes to background compilers such as NVIDIA HPC SDK for GPU and OpenARC research compiler for FPGA, then assembles binaries for the final object with FPGA bitstream file. In this paper, we present the concept, design, implementation, and performance evaluation of a practical astrophysics simulation code where we successfully enhanced the performance up to 10 times faster than the GPU-only solution.

Boku, Taisuke↗
DOE OSTI2026

MAGIC: M arching Cubes Isosurface Uncertainty Visualization for G auss i an Uncertain Data With Spatial C orrelation

Here, in this paper, we study the propagation of data uncertainty through the marching cubes algorithm for isosurface visualization for correlated uncertain data. Consideration of correlation has been shown paramount for avoiding errors in uncertainty quantification and visualization in multiple prior studies. Although the problem of isosurface uncertainty with spatial data correlation has been previously addressed, there are two major limitations to prior treatments. First, there are no analytical formulations for uncertainty quantification of isosurfaces when the data uncertainty is characterized by a Gaussian distribution with spatial correlation. Second, as a consequence of the lack of analytical formulations,existing techniques resort to a Monte Carlo sampling approach, which is expensive and difficult to integrate into visualization tools. To address these limitations, we present a closed-form framework to efficiently derive uncertainty in marching cubes level-sets for Gaussian uncertain data with spatial correlation (MAGIC). To derive closed-form solutions, we leverage the Hinkley's derivation on the ratio of Gaussian distributions. With our analytical framework, we achieve a significant speed-up and enhanced accuracy of uncertainty quantification over classical Monte Carlo methods. We further accelerate our analytical solutions using many-core processors to achieve speed-ups up to 585× and integrability with production visualization tools for broader impact. We demonstrate the effectiveness of our correlation-aware uncertainty framework through experiments on meteorology, urban flow, and astrophysics simulation datasets.

Gaussian↗
NASA NTRS1995

Probing the structure of cometary ice

Computer simulations of bulk and vapor deposited amorphous ices are presented. The structure of the bulk low density amorphous ice is in good agreement with experiments on pressure disordered amorphous ice. Both the low density bulk ice and the vapor deposited ices exhibit strong ordering. Vapor deposition of hot (300 K) water molecules onto a cold (77 K) substrate yields less porous ices than deposition of cold (77 K) water molecules onto a cold substrate. Both vapor deposited ices are more porous than the bulk amorphous ice. The structure of bulk high density amorphous ice is only in fair agreement with experimental results. Attempts to simulate high density amorphous ice via vapor deposition were not successful. Electron diffraction results on vapor deposited amorphous ice indicate that the temperature of the nucleation of the cubic phase depends upon the amount of time between the deposition and the onset of crystallization, suggesting that freshly deposited ice layers reconstruct on time of the order of hours. The temperature dependence of the microporosity of the vapor deposited amorphous ices might affect laboratory experiments that are aimed at simulating astrophysical ices in the context of the origin of prebiotic organic material and its transport to the Earth.

Wilson, Michael A.↗
NASA NTRS2000

Recent Advances in Interstellar Astrochemistry

This presentation will address the possibility that molecules created in the interstellar medium may play a role in the origin and evolution of life on planetary surfaces. The argument will be presented in three main parts. First, I will review the various organic and volatile compounds that are now known or suspected to exist in a variety of circumstellar and interstellar environments. This information is largely the result of the combined applications of observational infrared and radio spectroscopy, laboratory astrophysical simulations, and theoretical astrochemistry. Second, I will devote additional attention to those materials that exhibit chemical properties that are of biological interest. Finally, I will discuss the evidence, largely gathered from the laboratory study of extraterrestrial materials (meteorites and cosmic dust), that interstellar materials, including organics, can survive the transition from the dense cloud phase into forming stellar systems. Once there, some of this material can be delivered largely unaltered to planetary surfaces where it can play key roles in the origin and subsequent evolution of life. Thus, the guiding principle of this presentation will be to attempt to 'put the astro into astrobiology.'

Sandford, Scott A.↗
NASA NTRS2002

Organic Synthesis in Space

This talk will review our current understanding of the synthesis of organic molecules in space, with particular emphasis on the synthesis of those compounds that may be of prebiotic interest. The talk will address the possibility that molecules created in the interstellar medium may play a role in the origin and evolution of life on planetary surfaces. The various organic and volatile compounds that are now known or suspected to exist in a variety of space environments (stellar outflows, the diffuse interstellar medium, dense molecular clouds, protostellar nebulae, and planetesimal parent bodies in planetary systems) will be reviewed. This information comes largely from the combined applications of observational infrared and radio spectroscopy, laboratory astrophysical simulations, and theoretical astrochemistry. This will be followed by a discussion of the evidence, largely gathered from the laboratory isotopic study of extraterrestrial materials (meteorites and cosmic dust), that interstellar materials, including organics, can and do survive the transition from the interstellar space into forming stellar systems. Once there, some of this material can be delivered largely unaltered to planetary surfaces where it can play key roles in the origin and subsequent evolution of life.

Sandford, Scott A.↗
NASA NTRS2003

Laboratory Measurements of Optical Properties of Micron Size Individual Dust Grains

A laboratory program is being developed at NASA Marshall Space Flight Center for experimental determination of the optical and physical properties individual dust grains in simulated astrophysical environments. The experimental setup is based on an electrodynamic balance that permits levitation of single 0.1 - 10 micron radii dust grains in a cavity evacuated to pressures of approx. 10(exp -6) torr. The experimental apparatus is equipped with observational ports for measurements in the UV, visible, and infrared spectral regions. A cryogenic facility for cooling the particles to temperature of approx. 10-50K is being installed. The current and the planned measurements include: dust charging processes, photoelectric emissions and yields with UV irradiation, radiation pressure measurements, infrared absorption and scattering properties, and condensation processes, involving the analogs of cosmic dust grains. Selected results based on photoemissions, radiation pressure, and other laboratory measurements will be presented.

Abbas, M. M.↗
NASA NTRS1995

Experimental Evaluation and Workload Characterization for High-Performance Computer Architectures

This research is conducted in the context of the Joint NSF/NASA Initiative on Evaluation (JNNIE). JNNIE is an inter-agency research program that goes beyond typical.bencbking to provide and in-depth evaluations and understanding of the factors that limit the scalability of high-performance computing systems. Many NSF and NASA centers have participated in the effort. Our research effort was an integral part of implementing JNNIE in the NASA ESS grand challenge applications context. Our research work under this program was composed of three distinct, but related activities. They include the evaluation of NASA ESS high- performance computing testbeds using the wavelet decomposition application; evaluation of NASA ESS testbeds using astrophysical simulation applications; and developing an experimental model for workload characterization for understanding workload requirements. In this report, we provide a summary of findings that covers all three parts, a list of the publications that resulted from this effort, and three appendices with the details of each of the studies using a key publication developed under the respective work.

El-Ghazawi, Tarek A.↗
NASA NTRS2004

Cloud Detection from Satellite Imagery: A Comparison of Expert-Generated and Automatically-Generated Decision Trees

Automated cloud detection and tracking is an important step in assessing global climate change via remote sensing. Cloud masks, which indicate whether individual pixels depict clouds, are included in many of the data products that are based on data acquired on- board earth satellites. Many cloud-mask algorithms have the form of decision trees, which employ sequential tests that scientists designed based on empirical astrophysics studies and astrophysics simulations. Limitations of existing cloud masks restrict our ability to accurately track changes in cloud patterns over time. In this study we explored the potential benefits of automatically-learned decision trees for detecting clouds from images acquired using the Advanced Very High Resolution Radiometer (AVHRR) instrument on board the NOAA-14 weather satellite of the National Oceanic and Atmospheric Administration. We constructed three decision trees for a sample of 8km-daily AVHRR data from 2000 using a decision-tree learning procedure provided within MATLAB(R), and compared the accuracy of the decision trees to the accuracy of the cloud mask. We used ground observations collected by the National Aeronautics and Space Administration Clouds and the Earth s Radiant Energy Systems S COOL project as the gold standard. For the sample data, the accuracy of automatically learned decision trees was greater than the accuracy of the cloud masks included in the AVHRR data product.

Shiffman, Smadar↗
NASA NTRS2011

Atomic Data for Neutron-capture Elements I. Photoionization and Recombination Properties of Low-charge Selenium Ions

We present multi-configuration Breit-Pauli AUTOSTRUCTURE calculations of distorted-wave photoionization (PI) cross sections. and total and partial final-state resolved radiative recombination (RR) and dielectronic recombination (DR) rate coefficients for the first six ions of the trans-iron element Se. These calculations were motivated by the recent detection of Se emission lines in a large number of planetary nebulae. Se is a potentially useful tracer of neutron-capture nucleosynthesis. but accurate determinations of its abundance in photoionized nebulae have been hindered by the lack of atomic data governing its ionization balance. Our calculations were carried out in intermediate coupling with semi re1ativistic radial wavefunctions. PI and recombination data were determined for levels within the ground configuration of each ion, and experimental PI cross-section measurements were used to benchmark our results. For DR, we allowed (Delta)n = 0 core excitations, which are important at photoionized plasma temperatures. We find that DR is the dominant recombination process for each of these Se ions at temperatures representative of photoionized nebulae (approx.10(exp 4) K). In order to estimate the uncertainties of these data, we compared results from three different configuration-interaction expansions for each ion, and also tested the sensitivity of the results to the radial scaling factors in the structure calculations. We find that the internal uncertainties are typically 30-50% for the direct PI cross sections and approx.10% for the computed RR rate coefficients, while those for low-temperature DR can be considerably larger (from 15-30% up to two orders of magnitude) due to the unknown energies of near-threshold autoionization resonances. These data are available at the CDS, and fitting coefficients to the total RR and DR rate coefficients are presented. The results are suitable for incorporation into photoionization codes used to numerically simulate astrophysical nebulae, and will enable robust determinations of nebular Se abundances.

Sterling, N. C.↗
NASA NTRS2012

Imposing a Lagrangian Particle Framework on an Eulerian Hydrodynamics Infrastructure in Flash

In many astrophysical simulations, both Eulerian and Lagrangian quantities are of interest. For example, in a galaxy cluster merger simulation, the intracluster gas can have Eulerian discretization, while dark matter can be modeled using particles. FLASH, a component-based scientific simulation code, superimposes a Lagrangian framework atop an adaptive mesh refinement Eulerian framework to enable such simulations. The discretization of the field variables is Eulerian, while the Lagrangian entities occur in many different forms including tracer particles, massive particles, charged particles in particle-in-cell mode, and Lagrangian markers to model fluid structure interactions. These widely varying roles for Lagrangian entities are possible because of the highly modular, flexible, and extensible architecture of the Lagrangian framework. In this paper, we describe the Lagrangian framework in FLASH in the context of two very different applications, Type Ia supernovae and galaxy cluster mergers, which use the Lagrangian entities in fundamentally different ways.

Dubey, A.↗
NASA NTRS1977

Laboratory simulation of infrared astrophysical features

Intermediate resolution emission spectroscopy was used to study a group of 9 terrestrial silicates, 1 synthetic silicate, 6 meteorites and 2 lunar soils; comparisons were made with the intermediate resolution spectra of Comet Kohoutek in order to determine which materials best simulate the 10um astrophysical feature. Mixtures of silicates which would yield spectra matching the spectrum of the comet in the 10um region include: (1) A hydrous layer lattice silicate in combination with a high temperature condensate; (2) an amorphous magnesium silicate in combination with a high temperature condensate and (3) glassy olivine and glassy anorthite in approximately equal proportions.

Rose, L. A.↗
NASA NTRS1979

Laboratory simulation of infrared astrophysical features

Laboratory infrared emission and absorption spectra have been taken of terrestrial silicates, meteorites, and lunar soils in the form of micrometer and submicrometer grains. The emission spectra were taken in a way that imitates telescopic observations. The purpose was to see which materials best simulate the 10-micron astrophysical feature. The emission spectra of dunite, fayalite, and Allende give a good fit to the 10-micron broadband emission feature of comets Bennett and Kohoutek. A study of the effect of grain size on the presence of the 10-micron emission feature of dunite shows that for particles larger than 37 microns no feature is seen. The emission spectrum of the Murray meteorite, a Type 2 carbonaceous chrondrite, is quite similar to the intermediate-resolution spectrum of comet Kohoutek in the 10-micron region. Hydrous silicates or amorphous magnesium silicates in combination with high-temperature condensates, such as olivine or anorthite, would yield spectra that match the intermediate-resolution spectrum of comet Kohoutek in the 10-micron region. Glassy olivine and glassy anorthite in approximately equal proportions would also give a spectrum that is a good fit to the cometary 10-micron feature.

Rose, L. A.↗
DOE OSTI2021

Well-balanced treatment of gravity in astrophysical fluid dynamics simulations at low Mach numbers

We report accurate simulations of flows in stellar interiors are crucial to improving our understanding of stellar structure and evolution. Because the typically slow flows are merely tiny perturbations on top of a close balance between gravity and the pressure gradient, such simulations place heavy demands on numerical hydrodynamics schemes. We demonstrate how discretization errors on grids of reasonable size can lead to spurious flows orders of magnitude faster than the physical flow. Well-balanced numerical schemes can deal with this problem. Three such schemes were applied in the implicit, finite-volume SEVEN-LEAGUE HYDRO code in combination with a low-Mach-number numerical flux function. We compare how the schemes perform in four numerical experiments addressing some of the challenges imposed by typical problems in stellar hydrodynamics. We find that the α-β and deviation well-balancing methods can accurately maintain hydrostatic solutions provided that gravitational potential energy is included in the total energy balance. They accurately conserve minuscule entropy fluctuations advected in an isentropic stratification, which enables the methods to reproduce the expected scaling of convective flow speed with the heating rate. The deviation method also substantially increases accuracy of maintaining stationary orbital motions in a Keplerian disk on long timescales. The Cargo–LeRoux method fares substantially worse in our tests, although its simplicity may still offer some merits in certain situations. Overall, we find the well-balanced treatment of gravity in combination with low Mach number flux functions essential to reproducing correct physical solutions to challenging stellar slow-flow problems on affordable collocated grids.

79 ASTRONOMY AND ASTROPHYSICS↗