Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational astrophysics”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Computational Astrophysics in the Era of Technological Heterogeneity [Slides]

The HPC landscape is changing and we’re headed toward an era where compute specialization will be prevalent. There are opportunities for co-design that can influence this future. Technological heterogeneity will be a major challenge unless we shift our approach to developing the computational tools for astrophysics. Parthenon provides convenient functionality and a bright future for block structured AMR applications. Phoebus is a new (soon-to-be) open source code for relativistic astro that promises excellent performance, portability, and unique physics capabilities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Modeling Solids in Nuclear Astrophysics with Smoothed Particle Hydrodynamics

Abstract Smoothed particle hydrodynamics (SPH) is a frequently applied tool in computational astrophysics to solve the fluid dynamics equations governing the systems under study. For some problems, for example when involving asteroids and asteroid impacts, the additional inclusion of material strength is necessary in order to accurately describe the dynamics. In compact stars, that is white dwarfs and neutron stars, solid components are also present. Neutron stars have a solid crust, which is the strongest material known in nature. However, their dynamical evolution, when modeled via SPH or other computational fluid dynamics codes, is usually described as a purely fluid dynamics problem. Here, we present the first 3D simulations of neutron star crustal toroidal oscillations including material strength with the Los Alamos National Laboratory SPH code FleCSPH. In the first half of the paper, we present the numerical implementation of solid material modeling together with standard tests. The second half is on the simulation of crustal oscillations in the fundamental toroidal mode. Here, we dedicate a large fraction of the paper to approaches that can suppress numerical noise in the solid. If not minimized, the latter can dominate the crustal motion in the simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Towards Exascale Astrophysics of Mergers and Supernovae (TEAMS)

The TEAMS project brought together cutting-edge simulations, theoretical insights, and collaborative efforts to deepen our understanding of some of the universe’s most extreme phenomena—supernovae, neutron star mergers, and the powerful signals they emit. Using one of the largest suites of 3D supernova simulations ever conducted, researchers uncovered new insights into how massive stars explode, how those explosions vary by stellar mass, and what conditions lead to the birth of neutron stars or black holes. They also studied the radiation and gravitational wave signals emitted during these events, revealing how future observations can be used to uncover what happens deep inside collapsing stars. The team developed improved tools for modeling how light and neutrinos behave in such explosive environments, enabling more accurate predictions of what astronomers might observe. Work also explored how the chemical composition and geometry of kilonovae—the visible explosions that follow neutron star mergers—influence their signals and can reveal the origins of heavy elements like gold. These efforts not only advanced scientific knowledge, but also trained a new generation of researchers at the intersection of astrophysics, computational science, and nuclear theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Resolution requirements for numerical modeling of neutrino quantum kinetics

Neutrino quantum kinetics is a rapidly evolving field in computational astrophysics, with a primary focus on collective neutrino oscillations in core-collapse supernovae and postmerger phases of binary neutron star mergers. In recent years, there has been considerable debate concerning resolution dependence in numerical simulations. In this paper, we conduct a comprehensive resolution study in both angular- and spatial directions by using two independent schemes of quantum kinetic neutrino transport: finite volume and pseudospectral methods. We complement our discussion by linear stability analysis including inhomogeneous modes. Our result suggests that decreasing spatial resolutions underestimates the growth of flavor instability, and then leads to wrong asymptotic states of flavor conversions, which potentially has a critical impact on astrophysical consequences. We further delve into numerical results of low resolution simulations, that reveals the underlying mechanism responsible for numerical artifacts caused by insufficient resolutions. Finally, this study settles the debate on requirements of resolutions and serves as a guideline for numerical modeling of quantum kinetic neutrino transport.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Not-quite-transcendental Functions for Logarithmic Interpolation of Tabulated Data

From tabulated nuclear and degenerate equations of state to photon and neutrino opacities and nuclear reaction rates, tabulated data is ubiquitous in computational astrophysics. The dynamic range that must be covered by these tables typically spans many orders of magnitude. Here we present a novel strategy for accurately and performantly interpolating tabulated data that spans these large dynamic ranges. We demonstrate the efficacy of this strategy in tabulated lookups for nuclear and terrestrial equations of state. We show that this strategy is a faster drop-in replacement for linear interpolation of logarithmic grids.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

High-frequency heating of the solar wind triggered by low-frequency turbulence

The fast solar wind’s high speeds and non-thermal features require that considerable heating occurs well above the Sun’s surface. Two leading theories seem incompatible: low-frequency ‘Alfvénic’ turbulence, which transports energy outwards and is observed ubiquitously by spacecraft but seems insufficient to explain the observed dominance of ion over electron heating; and high-frequency ion-cyclotron waves, which explain the non-thermal heating of ions but lack an obvious source. Here we argue that the recently proposed ‘helicity barrier’ effect, which limits electron heating by inhibiting the turbulent cascade of energy to the smallest scales, can unify these two paradigms. Our six-dimensional simulations show how the helicity barrier causes the large-scale energy to grow through time, generating small parallel scales and high-frequency ion-cyclotron-wave heating from low-frequency turbulence, while simultaneously explaining various other long-standing observational puzzles. Furthermore, the predicted causal link between plasma expansion and the ion-to-electron heating ratio suggests that the helicity barrier could contribute to key observed differences between fast and slow wind streams.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

On Using Linux Kernel Huge Pages with FLASH, an Astrophysical Simulation Code

We present efforts at improving the performance of FLASH, a multi-scale, multi-physics simulation code principally for astrophysical applications, by using huge pages on Ookami, an HPE Apollo 80 A64FX platform. FLASH is written principally in modern Fortran and makes use of the PARAMESH library to manage a block-structured adaptive mesh. We explored options for enabling the use of huge pages with several compilers, but we were only able to successfully use huge pages when compiling with the Fujitsu compiler. As a result, the use of huge pages substantially reduced the number of translation lookaside buffer misses, but overall performance gains were marginal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Performance of an Astrophysical Radiation Hydrodynamics Code under Scalable Vector Extension Optimization

We present results of a performance study of an astrophysical radiation hydrodynamics code, V2D, on the Arm-based A64FX processor developed by Fujitsu. The code solves sparse linear systems, a task for which the A64FX architecture should be well suited. Here, we performed the performance analysis study on Ookami, an Apollo 80 platform utilizing the A64FX processor. We explored several compilers and performance anal-ysis packages and found the code did not perform as expected under scalable vector extension optimization, suggesting that a “deeper dive” into analyzing the code is worthwhile. However, a simple driver program that exercised basic sparse linear algebra routines used by V2D did show significant speedup with the use of the scalable vector extension optimization. We present the initial results from the study which used V2D on a relatively simple test problem that emphasized the repeated solution of sparse linear systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Nuclear Reaction Network WinNet

Abstract We present the state-of-the-art single-zone nuclear reaction network WinNet , which is capable of calculating the nucleosynthetic yields of a large variety of astrophysical environments and conditions. This ranges from the calculation of the primordial nucleosynthesis, where only a few nuclei are considered, to the ejecta of neutron star mergers with several thousands of involved nuclei. Here we describe the underlying physics and implementation details of the reaction network. We additionally present the numerical implementation of two different integration methods, the implicit Euler method and Gears method, along with their advantages and disadvantages. We furthermore describe basic example cases of thermodynamic conditions that we provide together with the network and demonstrate the reliability of the code by using simple test cases. With this publication, WinNet will be publicly available and open source at GitHub and Zenodo.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

XXVIIth International Conference on Supersymmetry and Unification of Fundamental Interactions (SUSY 2019) (Final Report)

Supersymmetry (SUSY) is one of most elegant extensions of the Standard Model (SM) and explains the puzzles of the SM by providing a candidate to explain the dark matter content of the universe, allowing scientists to understand the origin of the electroweak scale requiring the top mass to be around 170 GeV and leading to the unification of forces at a grand unified scale. Further the minimal supersymmetric standard model (MSSM) predicts the Higgs boson mass to be less than 135 GeV. The discovery of the Higgs Boson with mass around 125 GeV at the LHC has provided a major support to SUSY ideas. Searches for SUSY are ongoing at the Large Hadron Collider (LHC). Direct and indirect dark matter experiments are searching for a particle dark matter candidate which arises most naturally in SUSY models. Proton decay predicted by SUSY grand unified theories is being searched for at deep underground experiments. In addition, recent advances in neutrino and dark matter physics, observational astrophysics, precision cosmology and the promising new window into the cosmos opened by the direct detection of gravitational waves, have brought new ideas on the potential connections between new fundamental particles and our understanding of their impact on the early universe and its evolution. At present, the major questions include: Is SUSY still the best candidate for models beyond the SM? Do we have any well motivated alternative to SUSY? Have we exhausted all possibilities to search for new physics at high and low energy scales? XXVIIth International Conference on Supersymmetry and Unification of Fundamental Interactions (SUSY 2019), hosted by Texas A&M University – Corpus Christi during May 20-24, 2019, provided a unique venue to discus and understand the status of SUSY, connection between particle physics and cosmology, supersymmetry and its alternative, Higgs sector, neutrino sector, flavor sector, dark matter, electroweak phase transition, astroparticle physics, gravitational waves and string theory. Discussion of results from the LHC, recent neutrino experiments and observations, direct and indirect dark matter detection experiments, detection of gravitational waves, data from particle colliders, as well as measurements of the CMB and Large Scale Structure were an integral part of SUSY 2019. To ensure the younger participants will benefit from the conference the most, the conference was preceded by the 4 day long pre-SUSY summer school for graduate students and postdocs. The invited speakers were leading scientists in the fields of SUSY interest. The school took place on Texas A&M University – Corpus Christi campus during the week prior the SUSY 2019 conference (May 15 – 18, 2019). Since its inception in 1993, SUSY has become one of the most important and widely attended international meetings in high energy physics, devoted to new ideas in fundamental particle physics. SUSY 2019 brought together approximately 250 scientists, theorists, phenomenologists, experimentalists and cosmologists, (including over 60 graduate students and 70 postdocs) representing 22 nations: Australia, Belgium, Canada, Chile, China, Colombia, France, Germany, India, Italy, Japan, Mexico, Peru, Portugal, Romania, South Korea, Spain, Sweden, Switzerland, Taiwan, United Kingdom and United States. SUSY 2019 provided a stimulating venue for the exchange of scientific ideas among experts in dark matter, neutrino physics, particle physics, astrophysics and cosmology. The following scientific topics were delivered during SUSY 2019 in form of 44 plenary talks and over 200 parallel talks: Unification of Forces; Electroweak, Top and Higgs Physics; Precision Calculations and MC tools; BSM in Flavor Physics; Neutrino Masses: Models and Phenomenology; Cosmology and Gravitational Waves; Dark Matter, Astroparticle Physics; Formal Field Theory and Strings; Alternatives to Supersymmetry; Quantum Information: Machine Learning/Big Data. 28 talks were given during the pre-SUSY program related to the following topics: Neutrino Physics; Big Data; Collider Physics & SUSY; String Phenomenology; Cosmology; Dark Matter; SUSY Models and Phenomenology

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

From Interstellar Cloud to Star to Laboratory: Frontier HEDP Studies of Magnetized Colliding Plasma Flows with Strong Radiative Cooling (Final Report)

High-speed flows of gas are the norm in astrophysics where exploding stars, winds from the regions around black holes and streams of galactic material create “supersonic flows”. As the flows collide, strong shock waves are formed where gas is compressed and heated. The regions behind these shock waves are of great interest to scientists. New stars can form in these shock regions in the context of galaxies. In the context of planets like our Earth, the shocks waves come from the interaction of the solar wind with our planets’ magnetic field causing Aurora and placing astronauts and satellites in jeopardy. Thus, understanding supersonic flows of magnetized plasmas is essential for progress in many fields of astronomy and space physics (relevant to planets and the Earth). Remarkably, we can also create these kinds of flows in the laboratory using Z-pinches and laser driven experiments that include magnetic fields. The exploration of these forms of “High Energy Density Physics” (HEDP) flows is also critical in order to achieve Inertial Confinement Fusion, which is a step towards developing a sustainable source of Fusion energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Dynamics of Laterally Propagating Flames in X-Ray Bursts. I. Burning Front Structure

We investigate the structure of laterally propagating flames through the highly stratified burning layer in an X-ray burst. Two-dimensional hydrodynamics simulations of flame propagation are performed through a rotating plane-parallel atmosphere, exploring the structure of the flame. We discuss the approximations needed to capture the length and timescales at play in an X-ray burst and describe the flame acceleration observed. Finally, our studies complement other multidimensional studies of burning in X-ray bursts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Insights From the Software Design of a Multiphysics Multicomponent Scientific Code

Using simulations for scientific discovery requires that the software used in the simulations undergoes a rigorous design and development process similar to that of the lab instruments in the experimental sciences. To devise a good design methodology, it is critical to understand the requirements, constraints, and challenges. Furthermore, this article describes insights from the long-term stewardship of a multiphysics multicomponent software, FLASH, that was designed more than 20 years ago for astrophysics, now serves multiple communities, and has been successful in adapting to the changing world of high-performance computing.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Can Cooling and Heating Functions Be Modeled with Homogeneous Radiation Fields?

Abstract Cooling and heating functions describe how radiative processes impact the thermal state of a gas as a function of its temperature and other physical properties. In a most general case the functions depend on the detailed distributions of ionic species and on the radiation spectrum. Hence, these functions may vary on a very wide range of spatial and temporal scales. In this paper, we explore cooling and heating functions between 5 ≤ z ≤ 10 in simulated galaxies from the Cosmic Reionization On Computers project. We compare three functions: (1) the actual cooling and heating rates of hydrodynamic cells as a function of cell temperature, (2) the median cooling and heating functions computed using median interstellar medium (ISM) properties (median ISM), and (3) the median of the cooling and heating functions of all gas cells (instantaneous). We find that the median ISM and instantaneous approaches to finding a median cooling and heating function give identical results within the spread due to cell-to-cell variation. However, the actual cooling (heating) rates experienced by the gas at different temperatures in the simulations do not correspond to either summarized cooling (heating) functions. In other words, the thermodynamics of the gas in the simulations cannot be described by a single set of a cooling plus a heating function with a spatially constant radiation field that could be computed with common tools, such as CLOUDY.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

In-flight performance of S PIDER'S 280-GHz receivers

S PIDER is a balloon-borne instrument designed to map the cosmic microwave background at degree-angular scales in the presence of Galactic foregrounds. S PIDER has mapped a large sky area in the Southern Hemisphere using more than 2000 transition-edge sensors (TESs) during two NASA Long Duration Balloon flights above the Antarctic continent. During its first flight in January 2015, S PIDER observed in the 95 GHz and 150 GHz frequency bands, setting constraints on the B-mode signature of primordial gravitational waves. Its second flight in the 2022-23 season added new receivers at 280 GHz, each using an array of TESs coupled to the sky through feedhorns formed from stacks of silicon wafers. Here, these receivers are optimized to produce deep maps of polarized Galactic dust emission over a large sky area, providing a unique data set with lasting value to the field. In this work, we describe the instrument’s performance during S PIDER'S second flight.

280 GHz cosmology↗
DOE OSTI2024

Light scalars at the cosmological collider

We study the self-energies of weakly interacting scalar fields in de Sitter space with one field much lighter than the Hubble scale. We argue that self-energies drastically simplify in this light limit. We illustrate this in theories with two scalar fields, one heavy and one light, interacting with one another through either cubic or quartic interactions. To regulate infrared divergences, we compute these self-energies in Euclidean de Sitter space and then carefully analytically continue to Lorentzian signature. In particular, we do this for the most general renormalizable theory of two scalar fields with even interactions to leading order in the coupling and the mass of the light field. These self-energies are determined by de Sitter sunset diagrams, whose analytic structure and UV divergences we derive. Even at very weak couplings, the light field can substantially change how the heavy field propagates over long distances. The light field’s existence may then be inferred from how it modifies the heavy field’s oscillatory contribution to the primordial bispectrum in the squeezed limit, i.e. its cosmological collider signal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A reusable neural network pipeline for unidirectional fiber segmentation

Abstract Fiber-reinforced ceramic-matrix composites are advanced, temperature resistant materials with applications in aerospace engineering. Their analysis involves the detection and separation of fibers, embedded in a fiber bed, from an imaged sample. Currently, this is mostly done using semi-supervised techniques. Here, we present an open, automated computational pipeline to detect fibers from a tomographically reconstructed X-ray volume. We apply our pipeline to a non-trivial dataset by Larson et al . To separate the fibers in these samples, we tested four different architectures of convolutional neural networks. When comparing our neural network approach to a semi-supervised one, we obtained Dice and Matthews coefficients reaching up to 98%, showing that these automated approaches can match human-supervised methods, in some cases separating fibers that human-curated algorithms could not find. The software written for this project is open source, released under a permissive license, and can be freely adapted and re-used in other domains.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

C HIMERA : A Massively Parallel Code for Core-collapse Supernova Simulations

Herein we provide a detailed description of the C HIMERA code, a code developed to model core collapse supernovae (CCSNe) in multiple spatial dimensions. The CCSN explosion mechanism remains the subject of intense research. Progress to date demonstrates that it involves a complex interplay of neutrino production, transport, and interaction in the stellar core, three-dimensional stellar core fluid dynamics and its associated instabilities, nuclear burning, and the fundamental physics of the neutrino–stellar core weak interactions and the equations of state of all stellar core constituents—particularly, the nuclear equation of state associated with core nucleons, both free and bound in nuclei. C HIMERA , by incorporating detailed neutrino transport, realistic neutrino–matter interactions, three-dimensional hydrodynamics, realistic nuclear, leptonic, and photonic equations of state, and a nuclear reaction network, along with other refinements, can be used to study the role of neutrino radiation, hydrodynamic instabilities, and a variety of input physics in the explosion mechanism itself. It can also be used to compute observables such as neutrino signatures, gravitational radiation, and the products of nucleosynthesis associated with CCSNe. The code contains modules for neutrino transport, multidimensional compressible hydrodynamics, nuclear reactions, a variety of neutrino interactions, equations of state, and modules to provide data for post-processing observables such as the products of nucleosynthesis, and gravitational radiation. C HIMERA is an evolving code, being updated periodically with improved input physics and numerical refinements. We detail here the current version of the code, from which future improvements will stem, which can in turn be described as needed in future publications.

79 ASTRONOMY AND ASTROPHYSICS↗