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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Photoionization of Benzene and Small Polycyclic Aromatic Hydrocarbons in Ultraviolet-Processed Astrophysical Ices: A Computational Study

We employed density functional theory (DFT) calculations to model the photoionization behavior of benzene and small polycyclic aromatic hydrocarbons when they are embedded in a matrix of water ice in order to investigate issues raised by recent experimental work by Gudipati and Allamandola. The ionization energies of benzene, naphthalene, anthracene, and pyrene were found to be lowered by 1.5-2.1 eV in water ice. Low-lying vertical electronic excitation energies were computed with time-dependent DFT for both neutral and ionized species and are found in both cases to be remarkably unaffected by the ice matrix. Chemical behavior in ultraviolet-photoprocessed ices is also discussed, with a focus on electron recombination and pathways leading to phenol and analogous products.

Woon, D. E.↗
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↗
NASA NTRS1991

Dynamo magnetic field modes in thin astrophysical disks - An adiabatic computational approximation

An adiabatic approximation is applied to the calculation of turbulent MHD dynamo magnetic fields in thin disks. The adiabatic method is employed to investigate conditions under which magnetic fields generated by disk dynamos permeate the entire disk or are localized to restricted regions of a disk. Two specific cases of Keplerian disks are considered. In the first, magnetic field diffusion is assumed to be dominated by turbulent mixing leading to a dynamo number independent of distance from the center of the disk. In the second, the dynamo number is allowed to vary with distance from the disk's center. Localization of dynamo magnetic field structures is found to be a general feature of disk dynamos, except in the special case of stationary modes in dynamos with constant dynamo number. The implications for the dynamical behavior of dynamo magnetized accretion disks are discussed and the results of these exploratory calculations are examined in the context of the protosolar nebula and accretion disks around compact objects.

Stepinski, T. F.↗
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↗
DOE OSTI2022

Inferring the Thermal History of the Intergalactic Medium from the Properties of the Hydrogen and Helium Lyα Forest

The filamentary network of intergalactic medium (IGM) gas that gives origin to the Lyα forest in the spectra of distant quasars encodes information on the physics of structure formation and the early thermodynamics of diffuse baryonic material. Here we use a massive suite of more than 400 high-resolution cosmological hydrodynamical simulations run with the Graphics Processing Unit–accelerated code Cholla to study the IGM at high spatial resolution maintained over the entire computational volume. The simulations capture a wide range of possible IGM thermal histories by varying the photoheating and photoionizing background produced by star-forming galaxies and active galactic nuclei. A statistical comparison of synthetic spectra with the observed 1D flux power spectra of hydrogen at redshifts 2.2 ≤ z ≤ 5.0 and with the helium Lyα opacity at redshifts 2.4 < z < 2.9 tightly constrains the photoionization and photoheating history of the IGM. By leveraging the constraining power of the available Lyα forest data to break model degeneracies, we find that the IGM experienced two main reheating events over 1.2 Gyr of cosmic time. For our best-fit model, hydrogen reionization completes by z R ≈ 6.0 with a first IGM temperature peak of T 0 ≃ 1.3 × 10 4 K and is followed by the reionization of He ii that completes by z R ≈ 3.0 and yields a second temperature peak of T 0 ≃ 1.4 × 10 4 K. We discuss how our results can be used to obtain information on the timing and the sources of hydrogen and helium reionization.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Connecting the Light Curves of Type IIP Supernovae to the Properties of Their Progenitors

Abstract Observations of core-collapse supernovae (CCSNe) reveal a wealth of information about the dynamics of the supernova ejecta and its composition but very little direct information about the progenitor. Constraining properties of the progenitor and the explosion requires coupling the observations with a theoretical model of the explosion. Here we begin with the CCSN simulations of Couch et al., which use a nonparametric treatment of the neutrino transport while also accounting for turbulence and convection. In this work we use the SuperNova Explosion Code to evolve the CCSN hydrodynamics to later times and compute bolometric light curves. Focusing on Type IIP SNe (SNe IIP), we then (1) directly compare the theoretical STIR explosions to observations and (2) assess how properties of the progenitor’s core can be estimated from optical photometry in the plateau phase alone. First, the distribution of plateau luminosities ( L 50 ) and ejecta velocities achieved by our simulations is similar to the observed distributions. Second, we fit our models to the light curves and velocity evolution of some well-observed SNe. Third, we recover well-known correlations, as well as the difficulty of connecting any one SN property to zero-age main-sequence mass. Finally, we show that there is a usable, linear correlation between iron core mass and L 50 such that optical photometry alone of SNe IIP can give us insights into the cores of massive stars. Illustrating this by application to a few SNe, we find iron core masses of 1.3–1.5 M ⊙ with typical errors of 0.05 M ⊙ . Data are publicly available online on Zenodo: doi: 10.5281/zenodo.6631964 .

79 ASTRONOMY AND ASTROPHYSICS↗