Engineering Papers⌕ Search

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.

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

At least 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

The collapse and three-dimensional explosion of three-dimensional massive-star supernova progenitor models

ABSTRACT The explosion outcome and diagnostics of core-collapse supernovae depend sensitively on the nature of the stellar progenitor, but most studies to date have focused exclusively on one-dimensional, spherically symmetric massive star progenitors. We present some of the first core-collapse supernovae simulations of three-dimensional massive star supernovae progenitors, a 12.5- and a 15-M⊙ model, evolved in three dimensions from collapse to bounce through explosion with the radiation-hydrodynamic code fornax. We compare the results using those starting from three-dimensional progenitors to three-dimensional simulations of spherically symmetric, one-dimensional progenitors of the same mass. We find that the models evolved in three dimensions during the final stages of massive star evolution are more prone to explosion. The turbulence arising in these multidimensional initial models serves as seed turbulence that promotes shock revival. Detection of gravitational waves and neutrinos signals could reveal signatures of pre-bounce turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Gray two-moment neutrino transport: Comprehensive tests and improvements for supernova simulations

Aims. In this work we extended an energy-integrated neutrino transport method to facilitate efficient, yet precise, modeling of compact astrophysical objects. We particularly focus on core-collapse supernovae. Methods. We implemented a gray neutrino-transport framework from the literature into FLASH and performed a detailed evaluation of its accuracy in core-collapse supernova simulations. Based on comparisons with results from simulations using energy-dependent neutrino transport, we incorporated several improvements to the original scheme. Results. Our analysis shows that our gray neutrino transport method successfully reproduces key aspects from more complex energy-dependent transport across a variety of progenitors and equations of state. We find both qualitative and reasonable quantitative agreement with multi-group M1 transport simulations. However, the gray scheme tends to slightly favor shock revival. In terms of gravitational wave and neutrino signals, there is a good alignment with the energy-dependent transport, although we find 15–30% discrepancies in the average energy and luminosity of heavy-lepton neutrinos. Simulations using the gray transport are around four times faster than those using energy-dependent transport.

Astronomy & Astrophysics↗
DOE OSTI2022

Energy conserving and well-balanced discontinuous Galerkin methods for the Euler–Poisson equations in spherical symmetry

ABSTRACT This paper presents high-order Runge–Kutta (RK) discontinuous Galerkin methods for the Euler–Poisson equations in spherical symmetry. The scheme can preserve a general polytropic equilibrium state and achieve total energy conservation up to machine precision with carefully designed spatial and temporal discretizations. To achieve the well-balanced property, the numerical solutions are decomposed into equilibrium and fluctuation components that are treated differently in the source term approximation. One non-trivial challenge encountered in the procedure is the complexity of the equilibrium state, which is governed by the Lane–Emden equation. For total energy conservation, we present second- and third-order RK time discretization, where different source term approximations are introduced in each stage of the RK method to ensure the conservation of total energy. A carefully designed slope limiter for spherical symmetry is also introduced to eliminate oscillations near discontinuities while maintaining the well-balanced and total-energy-conserving properties. Extensive numerical examples – including a toy model of stellar core collapse with a phenomenological equation of state that results in core bounce and shock formation – are provided to demonstrate the desired properties of the proposed methods, including the well-balanced property, high-order accuracy, shock-capturing capability, and total energy conservation.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1982

Galaxy clusters with multiple components. I - The dynamics of Abell 98

The nature and evolution of rich galaxy clusters containing multiple condensations is studied. Velocities for 13 members of A98 are presented and the determination of dynamical properties of the individual subclusters is substantially improved. The application of techniques for quantifying gravitational scale lengths are stressed without invocation of the symmetry arguments demanded by the usual parameter fitting procedures. The entire A98 cluster is modeled as a two-body system. The probability that A98 is bound is found to be 98 percent. From this model and the results of previous N-body simulations, it is shown that the most likely description of the system is that it reached maximum expansion about 3.5 billion years ago and is presently in a state of collapse. The subclusters will merge in another three billion years.

Beers, T. C.↗
NASA NTRS1994

Nonaxisymmetric evolution in protostellar disks

We present a two-dimensional, multigridded hydrodynamical simulation of the collapse of an axisymmetric, rotating, 1 solar mass protostellar cloud, which forms a resolved, hydrotastic disk. The code includes the effects of physical viscosity, radiative transfer and radiative acceleration but not magnetic fields. We examine how the disk is affected by the inclusion of turbulent viscosity by comparing a viscous simulation with an inviscid model evolved from the same initial conditions, and we derive a disk evolutionary timescale on the order of 300,000 years if alpha = 0.01. Effects arising from non-axisymmetric gravitational instabilities in the protostellar disk are followed with a three-dimensional SPH code, starting from the two-dimensional structure. We find that the disk is prone to a series of spiral instabilities with primary azimulthal mode number m = 1 and m = 2. The torques induced by these nonaxisymmetric structures elicit material transport of angular momentum and mass through the disk, readjusting the surface density profile toward more stable configurations. We present a series of analyses which characterize both the development and the likely source of the instabilities. We speculate that an evolving disk which maintains a minimum Toomre Q-value approximately 1.4 will have a total evolutionary span of several times 10(exp 5) years, comparable to, but somewhat shorter than the evolutionary timescale resulting from viscous turbulence alone. We compare the evolution resulting from nonaxisymmetric instabilities with solutions of a one-dimensional viscous diffusion equation applied to the initial surface density and temperature profile. We find that an effective alpha-value of 0.03 is a good fit to the results of the simulation. However, the effective alpha will depend on the minimum Q in the disk at the time the instability is activated. We argue that the major fraction of the transport characterized by the value of alpha is due to the action of gravitational torques, and does not arise from inherent viscosity within the smoothed particle hydrodynamics method.

Laughlin, Gregory↗
NASA NTRS1993

A cosmological hydrodynamic code based on the total variation diminishing scheme

An explicitly second-order finite difference code based on a total variation diminishing scheme for self-gravitating cosmological hydrodynamic systems is described. The code was developed to follow correctly the adiabatic changes of extremely supersonic preshock flows with a Mach number larger than 100 as well as very strong shocks. In calculations of the formation of a 1D Zel'dovich pancake, an energy accuracy of 1 percent is obtained for 32 cells per unit wavelength, and the accuracy reaches 0.01 percent as the number of cells approach 1024. It is shown that shocks are well resolved and separate cleanly the hot, dense, collapsed peaks from the cold, low-density, expanding voids.

Ryu, Dongsu↗
NASA NTRS2006

The Evolution of the Cluster X-ray Scaling Relations in the Wide Angle ROSAT Pointed Survey Sample at 0.6 < z < 1.0

The X-ray properties of a sample of 11 high-redshift (0.6 < z < 1 .O) clusters observed with Chardm and/or XMM-Newton are used to investigate the evolution of the cluster scaling relations. The observed evolution in the normalization of the L-T, M-T, M(sub 2)-T and M-L relations is consistent with simple self-similar predictions, in which the properties of clusters reflect the properties of the Universe at their redshift of observation. Under the assumption that the model of self-similar evolution is correct and that the local systems formed via a single spherical collapse, the high-redshift L-T relation is consistent with the high-z clusters having virialized at a significantly higher redshift than the local systems. The data are also consistent with the more realistic scenario of clusters forming via the continuous accretion of material. The slope of the L-T relation at high redshift (B = 3.32 +/- 0.37) is consistent with the local relation, and significantly steeper than the self-similar prediction of B = 2. This suggests that the same non-gravitational processes are responsible for steepening the local and high-z relations, possibly occurring universally at z is approximately greater than 1 or in the early stages of the cluster formation, prior to their observation. The properties of the intracluster medium at high redshift are found to be similar to those in the local Universe. The mean surface-brightness profile slope for the sample is Beta = 0.66 +/- 0.05, the mean gas mass fractions within R(sub 2500(z)) and R(200(z)) are 0.069 +/- 0.012 and 0.11 +/- 0.02, respectively, and the mean metallicity of the sample is 0.28 +/- 0.11 Z(sub solar).

Maughan, B. J.↗
NASA NTRS1992

Gravitational influences on the behavior of confined diffusion flames

Closed form, mathematical representations for the shapes of laminar diffusion flames stabilized on two-dimensional slot burners are derived for arbitrary gravitational levels. The solution technique employs the Schvab-Zeldovich conserved species variable. Effects of gravity on the velocity field are incorporated through the momentum equation; these effects, in turn, influence the species field. The results reduce smoothly to the classical Burke-Schumann solution as gravity approaches zero. It is shown that both flame height and width increase with decreasing gravity levels. Several interesting cases are considered explicitly; these include effects of negative and extremely high values of gravity. In the latter case, it is shown that the flame collapses onto the burner surface. The analysis finds application in the modeling of steady and unsteady gas and solid propellant diffusion flames.

Hegde, U.↗
DOE OSTI2024

Robustness of neutron star merger simulations to changes in neutrino transport and neutrino-matter interactions

Binary neutron star mergers play an important role in nuclear astrophysics: Their gravitational wave and electromagnetic signals carry information about the equation of state of cold matter above nuclear saturation density, and they may be one of the main sources of r-process elements in the Universe. Neutrino-matter interactions during and after merger impact the properties of these electromagnetic signals and the relative abundances of the produced r-process elements. Existing merger simulations are, however, limited in their ability to realistically model neutrino transport and neutrino-matter interactions. Here, we perform a comparison of the impact of the use of state-of-the art two-moment or Monte Carlo transport schemes on the outcome of merger simulations, for a single binary neutron star system with a short-lived neutron star remnant [(5–10) ms]. We also investigate the use of different reaction rates in the simulations. While the best transport schemes generally agree well on the qualitative impact of neutrinos on the system, differences in the behavior of the high-density regions can significantly impact the collapse time and the properties of the hot tidal arms in this metastable merger remnant. The chosen interaction rates and transport algorithm as well as recent improvements by Radice et al. to the two-moment algorithms can all contribute to changes at the (10–30)% level in the global properties of the merger remnant and outflows. Furthermore, the limitations of previous moment schemes fixed by Radice et al. also appear sufficient to explain the large difference that we observed in the production of heavy-lepton neutrinos in a previous comparison of Monte Carlo and moment schemes in the context of a low-mass binary neutron star system.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

A simulation of the intracluster medium with feedback from cluster galaxies

We detail method and report first results from a three-dimensional hydrodynamical and N-body simulation of the formation and evolution of a Coma-sized cluster of galaxies, with the intent of studying the history of the hot, X-ray emitting intracluster medium. Cluster gas, galaxies, and dark matter are included in the model. The galaxies and dark matter fell gravitational forces; the cluster gas also undergoes hydrodynamical effects such as shock heating and PdV work. For the first time in three dimensions, we include modeling of ejection of processed gas from the simulated galaxies by winds, including heating and heavy element enrichment. For comparison, we employ a `pure infall' simulation using the same initial conditions but with no galaxies or winds. We employ an extreme ejection history for galactic feedback in order to define the boundary of likely models. As expected, feedback raises the entropy of the intracluster gas, preventing it from collapsing to densities as high as those attained in the infall model. The effect is more pronounced in subclusters formed at high redshift. The cluster with feedback is always less X-ray luminous, but experiences more rapid luminosity evolution, than the pure infall cluster. Even employing an extreme ejection model, the final gas temperature is only approximately 15% larger than in the infall model. The radial temperature profile is very nearly isothermal within 1.5 Mpc. The cluster galaxies in the feedback model have a velocity dispersion approximately 15% lower than the dark matter. This results in the true ratio of specific energies in galaxies to gas being less than one, beta(sub spec) approximately 0.7. The infall model predicts beta(sub spec) approximately 1.2. Large excursions in these values occur over time, following the complex dynamical history of the cluster. The morphology of the X-ray emission is little affected by feedback. The emission profiles of both clusters are well described by the standard beta-model with beta(sub fit) approximately equal to 0.7 - 0.9. X-ray mass estimates based on the assumptions of hydrostatic equilibrium and the applicability of the beta-model are quite accurate in both cases. A strong, radial iron abundance gradient is present, which develops as a consequence of the steepening of the galaxy density profile over time. Spectroscopic observations using nonimaging detectors with wide (approximately 45 min) fields of view dramatically smear the gradient. Observations with arcminute resolution, made available with the ASCA satellite, would readily resolve the gradient.

Metzler, Christopher A.↗
DOE OSTI2024

Classical field approximation of ultralight dark matter: Quantum break times, corrections, and decoherence

The classical field approximation is widely used to better understand the predictions of ultralight dark matter. Here, in this work, we use the truncated Wigner approximation method to test the classical field approximation of ultralight dark matter. This method approximates a quantum state as an ensemble of independently evolving realizations drawn from its Wigner function. The method is highly parallelizable and allows the direct simulation of quantum corrections and decoherence times in systems many times larger than have been previously studied in reference to ultralight dark matter. Our study involves simulation of systems in 1, 2, and 3 spatial dimensions. We simulate three systems, the condensation of a Gaussian random field in three spatial dimensions, a stable collapsed object in three spatial dimensions, and the merging of two stable objects in two spatial dimensions. We study the quantum corrections to the classical field theory in each case. We find that quantum corrections grow exponentially during nonlinear growth with the timescale being approximately equal to the system dynamical time. In stable systems the corrections grow quadratically. We also find that the primary effect of quantum corrections is to reduce the amplitude of fluctuations on the de Broglie scale in the spatial density. Finally, we find that the timescale associated with decoherence due to gravitational coupling to baryonic matter is at least as fast as the quantum corrections due to gravitational interactions. These results are consistent with the predictions of the classical field theory being accurate.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1979

On the origin of asteroids

A general scenario is described for the early history of the solar system. The primitive solar nebula is formed from the infall of gas from a collapsing interstellar cloud fragment. It becomes repeatedly unstable against collapse to form giant gaseous protoplanets. In the course of protoplanet evolution the center of the protoplanet enters a thermodynamic regime in which common rocky minerals become liquids; convection brings solids to the central region where a substantial fraction of them rain out to form a protoplanetary core. In the inner solar system protoplanetary envelopes are tidally stripped away, thus injecting into the solar nebula large equantities of chondrules and inclusions. Late in the development of the solar nebula, after most of the gas has disappeared, turbulence dies out and the small solids settle into a thin layer at midplane of the nebula. Gravitational instabilities in this layer form asteroidal and cometary bodies. Some further consequences of this scenario are discussed.

Cameron, A. G. W.↗
NASA NTRS2003

Formation of Giant Planets and Brown Dwarves

According to the prevailing core instability model, giant planets begin their growth by the accumulation of small solid bodies, as do terrestrial planets. However, unlike terrestrial planets, the growing giant planet cores become massive enough that they are able to accumulate substantial amounts of gas before the protoplanetary disk dissipates. Models predict that rocky planets should form in orbit about most stars. It is uncertain whether or not gas giant planet formation is common, because most protoplanetary disks may dissipate before solid planetary cores can grow large enough to gravitationally trap substantial quantities of gas. Ongoing theoretical modeling of accretion of giant planet atmospheres, as well as observations of protoplanetary disks, will help decide this issue. Observations of extrasolar planets around main sequence stars can only provide a lower limit on giant planet formation frequency . This is because after giant planets form, gravitational interactions with material within the protoplanetary disk may cause them to migrat~ inwards and be lost to the central star. The core instability model can only produce planets greater than a few jovian masses within protoplanetary disks that are more viscous than most such disks are believed to be. Thus, few brown dwarves (objects massive enough to undergo substantial deuterium fusion, estimated to occur above approximately 13 jovian masses) are likely to be formed in this manner. Most brown dwarves, as well as an unknown number of free-floating objects of planetary mass, are probably formed as are stars, by the collapse of extended gas/dust clouds into more compact objects.

Lissauer, Jack J.↗
NASA NTRS1992

Spin-up of a rapidly rotating star by angular momentum loss - Effects of general relativity

It has recently been shown that a rapidly rotating Newtonian star can spin up by radiating angular momentum. Extremely fast pulsars losing energy and angular momentum by magnetic dipole radiation or gravitational radiation may exhibit this behavior. Here, we show that this phenomenon is more widespread for rapidly rotating stars in general relativity. We construct and tabulate polytropic sequences of fully relativistic rotating stars of constant rest mass and entropy. We find that the range of adiabatic indices allowing spin-up extends somewhat above 4/3 because of the nonlinear effects of relativistic gravity. In addition, there is a new class of 'supramassive' stars which will inevitably spin up by losing angular momentum regardless of their equation of state. A supramassive star, spinning up via angular momentum loss, will ultimately evolve until it becomes unstable to catastrophic collapse to a black hole. Spin-up in a rapidly rotating star may thus be an observational precursor to such collapse.

Cook, Gregory B.↗
DOE OSTI2022

Numerical relativity simulations of prompt collapse mergers: Threshold mass and phenomenological constraints on neutron star properties after GW170817

We determine the threshold mass for prompt (no bounce) black hole formation in equal-mass neutron star (NS) mergers using a new set of 227 numerical relativity simulations. We consider 23 phenomenological and microphysical finite-temperature equations of state (EOS), including models with hyperons and first-order phase transitions to deconfined quarks. We confirm the existence of EOS-insensitive relations between the threshold mass, binary tidal parameter at the threshold (Λth), maximum mass of nonrotating NSs, and radii of reference mass NSs. We combine the EOS-insensitive relations, phenomenological constraints on NS properties, and observational data from GW170817 to derive an improved lower limit on radii of maximum mass and a 1.6 M ⊙ NS of 9.81 and 10.90 km, respectively. We also constrain the radius and quadrupolar tidal deformability (Λ) of a 1.4 M ⊙ NS to be larger than 10.74 km and 172, respectively. We consider uncertainties in all independent parameters—fitting coefficients as well as GW170817 masses while reporting the range of radii constraints. We discuss an approach to constrain the upper as well as lower limit of NS maximum mass using future binary NS detections and their identification as prompt or delayed collapse. With future observations, it will be possible to derive even tighter constraints on the properties of matter at and above nuclear density using the method proposed in this work.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1980

Criteria for fragmentation in a collapsing rotating cloud

The isothermal collapse of a rotating protostar with pressure and self-gravity is calculated by a hydrodynamic computer code that treats three space dimensions. A number of initial conditions are tested to determine the conditions under which a cloud is unstable to fragmentation. It is shown that fragmentation does not proceed to a significant extent during the first free-fall time, but that in most cases fragmentation has occurred by 1.5 to 2 times the initial free-fall time. The extent to which the density in the fragments is enhanced over that in their surroundings on this time scale depends on the initial ratio of thermal to gravitational energy, but is only weakly dependent on the type and amplitude of the perturbation and on the initial rotational energy.

Bodenheimer, P.↗
DOE OSTI2024

Crescendo beyond the horizon: More gravitational waves from domain walls bounded by inflated cosmic strings

Gravitational-wave (GW) signals offer a unique window into the dynamics of the early universe. GWs may be generated by the topological defects produced in the early universe, which contain information on the symmetry of UV physics. We consider the case in which a two-step phase transition produces a network of domain walls bounded by cosmic strings. Specifically, we focus on the case in which there is a hierarchy in the symmetry-breaking scales, and a period of inflation pushes the cosmic string generated in the first phase transition outside the horizon before the second phase transition. We show that the GW signal from the evolution and collapse of this string-wall network has a unique spectrum, and the resulting signal strength can be sizeable. In particular, depending on the model parameters, the resulting signal can show up in a broad range of frequencies and can be discovered by a multitude of future probes, including the pulsar timing arrays and space- and ground-based GW observatories. As an example that naturally gives rise to this scenario, we present a model with the first phase transition followed by a brief period of thermal inflation driven by the field responsible for the second stage of symmetry breaking. Here, the model can be embedded into a supersymmetric setup, which provides a natural realization of this scenario. In this case, the successful detection of the peak of the GW spectrum probes the soft supersymmetry breaking scale and the wall tension.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗