SEARCH · Engineering Papers

Results for “gravitational collapse”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1978

General-relativistic astrophysics

The overall relevance of general relativity to astrophysics is considered, and some of the knowledge about the ways in which general relativity should influence astrophysical systems is reviewed. Attention is focused primarily on finite-sized astrophysical systems, such as stars, globular clusters, galactic nuclei, and primordial black holes. Stages in the evolution of such systems and tools for studying the effects of relativistic gravity in these systems are examined. Gravitational-wave astronomy is discussed in detail, with emphasis placed on estimates of the strongest gravitational waves that bathe earth, present obstacles and future prospects for detection of the predicted waves, the theory of small perturbations of relativistic stars and black holes, and the gravitational waves such objects generate. Characteristics of waves produced by black-hole events in general, pregalactic black-hole events, black-hole events in galactic nuclei and quasars, black-hole events in globular clusters, the collapse of normal stars to form black holes or neutron stars, and corequakes in neutron stars are analyzed. The state of the art in gravitational-wave detection and characteristics of various types of detector are described.

Thorne, K. S.↗

NASA NTRS1994

The r-process and neutrino-heated supernova ejecta

As a neutron star is formed by the collapse of the iron core of a massive star, its Kelvin-Helmholtz evolution is characterized by the release of gravitational binding energy as neutrinos. The interaction of these neutrinos with heated material above the neutron star generates a hot bubble in an atmosphere that is nearly in hydrostatic equilibrium and heated, after approximately 10 s, to an entropy of S/N(sub AS)k greater than or approximately = 400. The neutron-to-proton ratio for material moving outward through this bubble is set by the balance between neutrino and antineutrino capture on nucleons. Because the electron antineutrino spectrum at this time is hotter than the electron neutrino spectrum, the bubble is neutron-rich (0.38 less than or approximately = Y(sub e) less than or approximately = 0.47). Previous work using a schematic model has shown that these conditions are well suited to the production of heavy elements by the r-process. In this paper we have advanced the numerical modeling of a 20 solar mass 'delayed' supernova explosion to the point that we can follow the detailed evolution of material moving through the bubble at the late times appropiate to r-process nucleosynthesis. The supernova model predicts a final kinetic energy for the ejecta of 1.5 x 10(exp 51) ergs and leaves behind a remnant with a baryon mass of 1.50 solar mass (and a gravitational mass of 1.445 solar mass). We follow the thermodynamic and compositional evolution of 40 trajectories in rho(t), T(t), Y(sub e)(t) for a logarithmic grid of mass elements for the last approximately = 0.03 solar mass to be ejected by the proto-neutron star down to the last less than 10(exp -6) solar mass of material expelled at up to approximately = 18 s after core collapse. We find that an excellent fit to the solar r-process abundance distribution is obtained with no adjustable parameters in the nucleosynthesis calculations. Moreover, the abundances are produced in the quantities required to account for the present Galactic abundances. However, at earlier times, this one-dimensional model ejects too much material with entropies S/N(sub A)k approximately 50 and Y(sub e) approximately 0.46. This leads to an acceptable over production of N = 50 nuclei, particularly Sr-88, Y-89, and Zr-90, relative to their solar abundances. We speculate on various means to avoid the early overproduction and/or ejection of N = 50 isotonic nuclei while still producing and ejecting the correct amount of r-process material.

Woosley, S. E.↗

DOE OSTI2024

Numerical challenges in modeling gravothermal collapse in Self-Interacting Dark Matter halos

When dark matter has a large cross section for self scattering, halos can undergo a process known as gravothermal core collapse, where the inner core rapidly increases in density and temperature. To date, several methods have been used to implement Self-Interacting Dark Matter (SIDM) in N-body codes, but there has been no systematic study of these different methods or their accuracy in the core-collapse phase. In this paper, we compare three different numerical implementations of SIDM, including the standard methods from the GIZMO and Arepo codes, by simulating idealized dwarf halos undergoing significant dark matter self interactions (σ/m = 50 cm 2 /g). When simulating these halos, we also vary the mass resolution, time-stepping criteria, and gravitational force-softening scheme. The various SIDM methods lead to distinct differences in a halo's evolution during the core-collapse phase, as each results in spurious scattering rate differences and energy gains/losses. The use of adaptive force softening for gravity can lead to numerical heating that artificially accelerates core collapse, while an insufficiently small simulation time step can cause core evolution to stall or completely reverse. Additionally, particle numbers must be large enough to ensure that the simulated halos are not sensitive to noise in the initial conditions. Even for the highest-resolution simulations tested in this study (10 6 particles per halo), we find that variations of order 10% in collapse time are still present. Here, the results of this work underscore the sensitivity of SIDM modeling on the choice of numerical implementation and motivate a careful study of how these results generalize to halos in a cosmological context.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Signature of hadron-quark crossover in binary-neutron-star mergers

Here, we study observational signatures of the hadron-quark crossover in binary-neutron-star mergers by numerical-relativity simulations with various mass configurations. We employ two equations of state (EoSs) for matter consistent with inference from the observational data. In the crossover scenario the EoS is softened in a density realized in binary-neutron-star mergers and is smoothly continued to quark matter. In the phase transition scenario without crossover, the EoS remains stiff and a first-order phase transition takes place in a density out of reach of mergers. A GW170817-like system forms a remnant massive neutron star in both scenarios, and it collapses into a black hole only in the crossover scenario due to the softening while gravitational-wave emission is strong. This difference is clearly reflected in the sudden shutdown of gravitational waves. For a given EoS, the lifetime of the merger remnant is determined primarily by the total mass of the system. Identifying these features in a variety of future events with the next generation of ground-based gravitational-wave detectors will enable us to clarify details of hadron-quark transition. The mass of the accretion disk surrounding the remnant black hole is affected not only by the lifetime of the remnant but also by the mass ratio of the system. Electromagnetic emission associated with the disk outflow will also be useful for detailed investigation of the hadron-quark transition.

QCD phase transitions↗

NASA NTRS1984

Formation of chondrules by drag heating in dust-enriched environments

The solar system which had a beginning, 4.6 billion years ago is discussed. The atoms that make up the solar system are much older than that, however: they were also part of the universe that existed before the solar system formed. The presence and location of the atoms of silicon, iron, calcium, that make up the Earth and other planets before there was a solar system is questioned. Astronomical it is shown that they existed in tiny dust grains which, along with molecules of gas, wandered about in the vast reaches of space between the stars of our galaxy. The solar system was formed when one particular cloud of interstellar dust and gas became dense enough to be gravitationally unstable: the gravitational pull that dust grains and gas molecules exerted on one another caused the cloud to collapse.

Wood, J. A.↗

DOE OSTI2023

Gravitational Waves from Rapid Structure Formation on Microscopic Scales before Matter-Radiation Equality

We propose a new scenario for the formation of asteroid-mass primordial black holes (PBHs). Our mechanism is based on the annihilation of the string-wall network associated with thebreaking of a U(1) global symmetry into a discrete ZN symmetry. If the potential hasmultiple local minima (N > 1) the network is stable, and the annihilation is guaranteed by a biasamong the different vacua. The collapse of the string-wall network is accompanied bycatastrogenesis, a large production of pseudo-Goldstone bosons (pGBs) — e.g. axions,ALPs, or majorons — gravitational waves, and PBHs. If pGBs rapidly decay into products thatthermalize, as predicted e.g. in the high-quality QCD axion and heavy majoron models, they do notcontribute to the dark matter population, but we show that PBHs can constitute 100% of the darkmatter. The gravitational wave background produced by catastrogenesis with heavy unstable axions,ALPs, or majorons could be visible in future interferometers.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Massive black holes seeded by dark matter: Implications for little red dots and gravitational wave signatures

Observations of supermassive black holes (SMBHs) at high redshifts challenge standard seeding scenarios. We examine a dissipative self-interacting dark matter (dSIDM) model in which gravothermal collapse leads to the formation of massive BH seeds ab initio. We utilize a semianalytical framework to predict properties of the dSIDM-seeded SMBH population. Billion solar mass quasars are reproduced along with low-mass faint active galactic nuclei (known as little red dots) with SMBH-to-galaxy stellar mass ratios consistent with recent James Webb Space Telescope observations. To match the abundance of the observed bright quasars, a subpercent-level duty cycle is suggested, implying a large population of dormant SMBHs. The gravitational wave (GW) signals from mergers of these massive SMBHs can be detected by LISA while remaining within the NANOGrav constraints on the GW background. These results provide testable signatures of DM-driven SMBH formation, offering a pathway to probe hidden-sector physics through SMBH and GW observables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1990

Fueling active galactic nuclei by magnetic braking

Recent detections of massive concentrations of molecular gas near the centers of galaxies hosting active nuclei suggest that these concentrations may be the source of accretion fuel for the nucleus. However, for that to be true, an angular momentum barrier must be overcome before the material in such a cloud can reach the nucleus. It is suggested that magnetic braking of the cloud may remove sufficient angular momentum to permit its material to draw considerably closer to the central object. The mechanism is particularly effective in the limit that the gas becomes self-gravitating because removal of a fraction of the initial angular momentum can lead to dynamical instability and collapse. Any small misalignment between the initial rotation axis of the cloud and the rotation axis of the galaxy can be substantially amplified as a result of the braking. It is argued that mass accretion onto the central object may occur in episodes, in some cases with a constant mass accretion rate during each episode.

Krolik, Julian H.↗

NASA NTRS2003

How Can The SN-GRB Time Delay Be Measured?

The connection between SNe and GRBs, launched by SN 1998bw / GRB 980425 and clinched by SN 2003dh / GRB 030329-with the two GRBs differing by a factor of approximately 50000 in luminosity-so far suggests a rough upper limit of approximately 1-2 days for the delay between SN and GRB. Only four SNe have had nonnegligible coverage in close coincidence with the initial explosion, near the W shock breakout: two Qpe II, and two Type IC, SN 1999ex and SN 1998bw. For the latter, only a hint of the minimum between the UV maximum and the radioactivity bump served to help constrain the interval between SN and GRB. Swift GRB alerts may provide the opportunity to study many SNe through the UV breakout phase: GRB 980425 look dikes -apparently nearby, low- luminosity, soft-spectrum, long-lag GRBs-accounted for half of BATSE bursts near threshold, and may dominate the Swift yield near threshold, since it has sensitivity to lower energies than did BATSE. The SN to GRB delay timescale should be better constrained by prompt UV/optical observations alerted by these bursts. Definitive delay measurements may be obtained if long-lag bursters are truly nearby: The SNe/GRBs could emit gravitational radiation detectable by LIGO-II if robust non-axisymmetric bar instabilities develop during core collapse, and/or neutrino emission may be detectable as suggested by Meszaros et al.

Norris, J. P.↗

NASA NTRS1994

Theory of protostellar accretion disks

I will present an overview of the current paradigm for the theory of gaseous accretion disks around young stars. Protostellar disks form from the collapse of rotating molecular cloud cores. The disks evolve via outward angular momentum transport provided by several mechanisms: gravitational instabilities, thermal convective turbulence, and magnetic stresses. I will review the conditions under which these mechanisms are efficient and consistent with the observed disk evolutionary timescales of several million years. Time permitting, I will discuss outbursts in protostellar disks (FU Orionis variables), the effect of planet formation on disk structure, and the dispersal of remnant gas.

Ruden, S.↗

DOE OSTI2025

Accurately simulating core-collapse self-interacting dark matter halos

The properties of satellite halos provide a promising probe for dark matter (DM) physics. Observations have motivated current efforts to explain surprisingly compact DM halos. If DM is not collisionless, but has strong self-interactions, halos can undergo gravothermal collapse, leading to higher densities in the central region of the halo. However, it is challenging to model this collapse phase from first principles. To improve on this, we sought to better understand the numerical challenges and convergence properties of self-interacting dark matter (SIDM) N-body simulations in the collapse phase. Especially, our aim was to better understand the evolution of satellite halos. To do so, we ran SIDM N-body simulations of a low-mass halo in isolation and within an external gravitational potential. The simulation set-up was motivated by the perturber of the stellar stream GD-1. We find that the halo evolution is very sensitive to energy conservation errors, and a SIDM kernel size that is too large can artificially speed up the collapse. Moreover, we demonstrate that the King model can describe the density profile at small radii for the late stages that we have simulated. Furthermore, for our most highly resolved simulation (N = 5 × 10 7 ) we have made the data public. It can serve as a benchmark. Overall, we find that the current numerical methods do not suffer from convergence problems in the late collapse phase and provide guidance on how to choose numerical parameters, for example that the energy conservation error is better kept well below 1%. This allows simulations to be run of halos that become concentrated enough to explain observations of GD-1-like stellar streams or strong gravitational lensing systems.

dark matter↗

DOE OSTI2024

Cosmological dynamics of string theory axion strings

The quantum chromodynamics (QCD) axion may solve the strong CP problem and explain the dark matter (DM) abundance of our Universe. The axion was originally proposed to arise as the pseudo-Nambu-Goldstone boson of global U⁢(1) PQ Peccei-Quinn (PQ) symmetry breaking, but axions also arise generically in string theory as zero modes of higher-dimensional gauge fields. In this work we show that string theory axions behave fundamentally differently from field theory axions in the early Universe. Field theory axions may form axion strings if the PQ phase transition takes place after inflation. In contrast, we show that string theory axions do not generically form axion strings. In special inflationary paradigms, such as D-brane inflation, string theory axion strings may form; however, their tension is parametrically larger than that of field theory axion strings. We then show that such QCD axion strings overproduce the DM abundance for all allowed QCD axion masses and are thus ruled out, except in scenarios with large warping. A loop-hole to this conclusion arises in the axiverse, where an axion string could be composed of multiple different axion mass eigenstates; a heavier eigenstate could collapse the network earlier, allowing for the QCD axion to produce the correct DM abundance and also generating observable gravitational wave signals.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1991

The biasing scheme in N-body simulations

A popular prescription for the galaxy formation, i.e., identifying peaks in the linear density field as sites of galaxy formation, is numerically justified under the Cold Dark Matter cosmogony. In particular, the peak-background scheme, originally developed to simplify analytic calculations of peak properties in a Gaussian random density field, is proven to be an excellent prescription for assigning peak tracers in low-resolution N-body simulations through many statistical comparisons. It is found that statistical properties of peak tracers allocated by the peak-background scheme closely resemble those of true peaks when the background smoothing length is about three times the galaxy scale and the correspondence improves as they evolve gravitationally. It is also shown that these peak particles are indeed strongly associated with nonlinear objects that are presently collapsed.

Park, Changbom↗

NASA NTRS2003

Mixing and Transport in the Solar Nebula

Boss & Vanhala (2000, 2001) prepared reviews of triggered collapse and injection models, using Prudence Foster's finite differences code at very high spatial resolution (440 x 1440 cells) to demonstrate the convergence of the R-T fingers in triggered injection models. A two dimensional hydrodynamical calculation with unprecedentedly high spatial resolution (960 x 2880 zones, or almost 3 million grid points) demonstrated that it suitable shock front can both trigger the collapse of an otherwise stable presolar cloud, and inject shock front particles into the collapsing cloud through the formation of what become Rayleigh-Taylor fingers of compressed fluid layers falling into the gravitational potential well of the growing protostar. These calculations suggest that heterogeneity derived from these R-T fingers will persist down to the scale of their injection onto the surface of the solar nebula. Haghighipour developed a numerical code capable of calculating the orbital evolution of dust grains of varied sizes in a gaseous nebula, subject to Epstein and Stokes drag as well as the self-gravity of the disk. In collaboration with the PI and George W. Wetherill, Haghighipour has been involved in development of a new idea on the possibility of rapid formation of ice giant planets via the disk instability mechanism. Haghighipour studied the stability of a five-body system consisting of the Sun and four protoplanets by numerically integrating their equations of motions. Using Levison and Duncan s SWIFT integrator, Haghighipour showed that, depending on the orbital parameters of the bodies, such a system can be stable for 0.1-10 Myr. Time periods of 1 Myr or more are long enough to be consistent with the time scale proposed for the formation of giant planets by the disk instability mechanism and the photoevaporation of the gaseous envelopes of the outermost protoplanets by a nearby OB star, resulting in the formation of ice giant planets. The PI has used his three dimensional models of marginally gravitationally unstable disks to study the preservation of isotopic heterogeneity in evolving protoplanetary disks. Such heterogeneity might arise from the infall onto the disk s surface of solids processed in the X-wind region of the disk, or derived from stellar nucleosynthesis and injected by R-T fingers. The technique used consists of solving a color equation, identical to the gas continuity equation, which follows the time evolution in three space dimensions of an arbitrarily placed initial color field, i.e., a dye inserted the disk. The models show that significant concentrations of color could persist for time periods of about a thousand years or more, even in the most dynamically active region of such a disk. Such a time period might be long enough for solids to coagulate and grow to significant sizes while retaining the isotopic signature of their birth region in the nebula.

Boss, Alan P.↗

DOE OSTI2023

Ab-initio General-relativistic Neutrino-radiation Hydrodynamics Simulations of Long-lived Neutron Star Merger Remnants to Neutrino Cooling Timescales

Abstract We perform the first 3D ab-initio general-relativistic neutrino-radiation hydrodynamics of a long-lived neutron star merger remnant spanning a fraction of its cooling timescale. We find that neutrino cooling becomes the dominant energy loss mechanism after the gravitational-wave dominated phase (∼20 ms postmerger). Electron flavor antineutrino luminosity dominates over electron flavor neutrino luminosity at early times, resulting in a secular increase of the electron fraction in the outer layers of the remnant. However, the two luminosities become comparable ∼20–40 ms postmerger. A dense gas of electron antineutrinos is formed in the outer core of the remnant at densities ∼10 14.5 g cm −3 , corresponding to temperature hot spots. The neutrinos account for ∼10% of the lepton number in this region. Despite the negative radial temperature gradient, the radial entropy gradient remains positive, and the remnant is stably stratified according to the Ledoux criterion for convection. A massive accretion disk is formed from the material squeezed out of the collisional interface between the stars. The disk carries a large fraction of the angular momentum of the system, allowing the remnant massive neutron star to settle to a quasi-steady equilibrium within the region of possible, stable, rigidly rotating configurations. The remnant is differentially rotating, but it is stable against the magnetorotational instability. Other MHD mechanisms operating on longer timescales are likely responsible for the removal of the differential rotation. Our results indicate the remnant massive neutron star is thus qualitatively different from a protoneutron stars formed in core-collapse supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Dynamics of a relativistic jet through magnetized media

The merger of two neutron stars (NSs) produces the emission of gravitational waves, the formation of a compact object surrounded by a dense and magnetized environment. If the binary undergoes delayed collapse a collimated and relativistic jet, which will eventually produce a short gamma-ray burst (SGRB), may be launched. The interaction of the jet with the environment has been shown to play a major role in shaping the structure of the outflow that eventually powers the gamma-ray emission. In this paper, we present a set of 2.5D RMHD simulations that follow the evolution of a relativistic non-magnetized jet through a medium with different magnetization levels, as produced after the merger of two NSs. We find that the predominant consequence of a magnetized ambient medium is that of suppressing instabilities within the jet and preventing the formation of a series of collimation shocks. One implication of this is that internal shocks lose efficiency, causing bursts with low-luminosity prompt emission. On the other hand, the jet-head velocity and the induced magnetization within the jet are fairly independent of the magnetization of the ambient medium. Future numerical studies with a larger domain are necessary to obtain light curves and spectra in order to better understand the role of magnetized media.

(magnetohydrodynamics) MHD↗

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↗

NASA NTRS1982

Axisymmetric collapse of rotating, isothermal clouds

The results of over 50 new models of the axisymmetric collapse of rotating, isothermal clouds are presented, with the following objectives: (1) to fully explore the initial conditions necessary for collapse from uniform density and uniform rotation, subject to constant volume and constant pressure boundary conditions; (2) to catalog the possible end states for cloud collapse from these initial conditions; and (3) to determine if there is a critical value of rotational energy/gravitational energy associated with ring formation, as appears to be the case for adiabatic clouds. Three end states are obtained: Bonnor-Ebert spheroids, rings and collapsing disks. The rings are formed with values of the ratio of rotational energy to the absolute value of the gravitational energy typically less than the Maclaurin spheroid value for dynamic instability to ring formation.

Boss, A. P.↗