Search NASA⌕ Search

SEARCH · Search NASA

Results for “General Relativity”

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 1,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Verification and validation of detonation-shock-dynamics relations for explosives described by general equation of state and chemical reaction models

Detonation shock dynamics is a powerful method to model the behaviour of High Explosives (HE). However in order to use this method, the underlying relationship between the local radius of curvature and the detonation speed must be known. Previous work has developed methods to calculate this effect using simple, single-step Arrhenius and polytropic gas, models for the chemical reaction and the equation of state, respectively. In recent years, more complex models for both reaction rates and equations of state have been developed which show better agreement with experimental data than these simple models, especially when considering condensed phase explosives.. This work presents the governing equations for solving these problems in a way that is generalised to use arbitrary equations of state as well as reaction models which may have more than a single step and multiple product species. This implementation is verified against exact solutions, demonstrating that the equations were implemented properly. The verified algorithm is then validated against experimental data and high fidelity simulations, showing that it is able to make accurate predictions in a regime where the underlying assumptions of the governing equations are valid. Importantly, this approach has many applications: from creating equivalent detonation shock dynamics models for existing reactive burn calibrations for HE; to developing new functional forms and calibrations of reactive burn models for condensed phase high explosives.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Effective Q criterion for disk stability in an external potential

The standard Q criterion (with Q > 1) describes the stability against local, axisymmetric perturbations in a disk supported by rotation and random motion. Most astrophysical disks, however, are under the influence of an external gravitational potential, which can significantly affect their stability. A typical example is a galactic disk embedded in a dark matter halo. Here, we do a linear perturbation analysis for a disk in an external potential and obtain a generalized dispersion relation and the effective stability criterion. An external potential, such as that due to the dark matter halo concentric with the disk, contributes to the unperturbed rotational field and significantly increases its stability. We obtain the values for the effective Q parameter for the Milky Way and for a low surface brightness galaxy, UGC 7321. We find that in each case the stellar disk by itself is barely stable and it is the dark matter halo that stabilizes the disk against local, axisymmetric gravitational instabilities. Thus, the dark matter halo is necessary to ensure local disk stability. This result has been largely missed so far because in practice the Q parameter for a galactic disk is obtained using the observed rotational field that already includes the effect of the halo.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Wide field near-infrared photometry of 12 galactic globular clusters: Observations versus models on the red giant branch

We present wide field near-infrared (near-IR) photometry of 12 Galactic globular clusters, typically extending from the tip of the cluster red giant branch (RGB) to the main sequence turnoff. Using recent homogenous values of cluster distance, reddening and metallicity, the resulting photometry is directly compared to the predictions of several recent libraries of stellar evolutionary models. Of the sets of models investigated, Dartmouth and Victoria-Regina models best reproduce the observed RGB morphology, albeit with offsets in J−K{sub S} color which vary in their significance in light of all sources of observational uncertainty. Therefore, we also present newly recalibrated relations between near-IR photometric indices describing the upper RGB versus cluster iron abundance as well as global metallicity. The influence of enhancements in alpha elements and helium are analyzed, and we find that the former affect the morphology of the upper RGB in accord with model predictions. Meanwhile, the empirical relations we derive are in good agreement with previous results, and minor discrepancies can likely be attributed to differences in the assumed cluster distances and reddenings. In addition, we present measurements of the horizontal branch (HB) and RGB bump magnitudes, finding a non-negligible dependence of the near-IR HB magnitude on cluster metallicity. Lastly, we discuss the influence of assumed cluster distances, reddenings and metallicities on our results, finding that our empirical relations are generally insensitive to these factors to within their uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

GENERAL-RELATIVISTIC SIMULATIONS OF THREE-DIMENSIONAL CORE-COLLAPSE SUPERNOVAE

We study the three-dimensional (3D) hydrodynamics of the post-core-bounce phase of the collapse of a 27 M{sub Sun} star and pay special attention to the development of the standing accretion shock instability (SASI) and neutrino-driven convection. To this end, we perform 3D general-relativistic simulations with a three-species neutrino leakage scheme. The leakage scheme captures the essential aspects of neutrino cooling, heating, and lepton number exchange as predicted by radiation-hydrodynamics simulations. The 27 M{sub Sun} progenitor was studied in 2D by Mueller et al., who observed strong growth of the SASI while neutrino-driven convection was suppressed. In our 3D simulations, neutrino-driven convection grows from numerical perturbations imposed by our Cartesian grid. It becomes the dominant instability and leads to large-scale non-oscillatory deformations of the shock front. These will result in strongly aspherical explosions without the need for large-scale SASI shock oscillations. Low-l-mode SASI oscillations are present in our models, but saturate at small amplitudes that decrease with increasing neutrino heating and vigor of convection. Our results, in agreement with simpler 3D Newtonian simulations, suggest that once neutrino-driven convection is started, it is likely to become the dominant instability in 3D. Whether it is the primary instability after bounce will ultimately depend on the physical seed perturbations present in the cores of massive stars. The gravitational wave signal, which we extract and analyze for the first time from 3D general-relativistic models, will serve as an observational probe of the postbounce dynamics and, in combination with neutrinos, may allow us to determine the primary hydrodynamic instability.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Fast rise of 'Neptune-size' planets (4-8 R {sub ⊕}) from P ∼ 10 to ∼250 days—statistics of Kepler planet candidates up to ∼0.75 au

We infer the period (P) and size (R{sub p} ) distribution of Kepler transiting planet candidates with R{sub p} ≥ 1 R {sub ⊕} and P < 250 days hosted by solar-type stars. The planet detection efficiency is computed by using measured noise and the observed time spans of the light curves for ∼120,000 Kepler target stars. We focus on deriving the shape of planet periods and radius distribution functions. We find that for orbital periods P > 10 days, the planet frequency dN{sub p} /dlog P for 'Neptune-size' planets (R{sub p} = 4-8 R {sub ⊕}) increases with period as ∝P {sup 0.7±0.1}. In contrast, dN{sub p} /dlog P for 'super-Earth-size' (2-4 R {sub ⊕}) as well as 'Earth-size' (1-2 R {sub ⊕}) planets are consistent with a nearly flat distribution as a function of period (∝P {sup 0.11±0.05} and ∝P {sup –0.10±0.12}, respectively), and the normalizations are remarkably similar (within a factor of ∼1.5 at 50 days). Planet size distribution evolves with period, and generally the relative fractions for big planets (∼3-10 R {sub ⊕}) increase with period. The shape of the distribution function is not sensitive to changes in the selection criteria of the sample. The implied nearly flat or rising planet frequency at long periods appears to be in disagreement with the sharp decline at ∼100 days in planet frequency for low-mass planets (planet mass m{sub p} < 30 M {sub ⊕}) recently suggested by the HARPS survey. Within 250 days, the cumulative frequencies for Earth-size and super-Earth-size planets are remarkably similar (∼28% and 25%), while Neptune-size and Jupiter-size planets are ∼7% and ∼3%, respectively. A major potential uncertainty arises from the unphysical impact parameter distribution of the candidates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

RADIAL STABILITY IN STRATIFIED STARS

We formulate within a generalized distributional approach the treatment of the stability against radial perturbations for both neutral and charged stratified stars in Newtonian and Einstein's gravity. We obtain from this approach the boundary conditions connecting any two phases within a star and underline its relevance for realistic models of compact stars with phase transitions, owing to the modification of the star's set of eigenmodes with respect to the continuous case.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

bhlight: GENERAL RELATIVISTIC RADIATION MAGNETOHYDRODYNAMICS WITH MONTE CARLO TRANSPORT

We present bhlight, a numerical scheme for solving the equations of general relativistic radiation magnetohydrodynamics using a direct Monte Carlo solution of the frequency-dependent radiative transport equation. bhlight is designed to evolve black hole accretion flows at intermediate accretion rate, in the regime between the classical radiatively efficient disk and the radiatively inefficient accretion flow (RIAF), in which global radiative effects play a sub-dominant but non-negligible role in disk dynamics. We describe the governing equations, numerical method, idiosyncrasies of our implementation, and a suite of test and convergence results. We also describe example applications to radiative Bondi accretion and to a slowly accreting Kerr black hole in axisymmetry.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

NEUTRINO-DRIVEN TURBULENT CONVECTION AND STANDING ACCRETION SHOCK INSTABILITY IN THREE-DIMENSIONAL CORE-COLLAPSE SUPERNOVAE

We conduct a series of numerical experiments into the nature of three-dimensional (3D) hydrodynamics in the postbounce stalled-shock phase of core-collapse supernovae using 3D general-relativistic hydrodynamic simulations of a 27 M{sub ⊙} progenitor star with a neutrino leakage/heating scheme. We vary the strength of neutrino heating and find three cases of 3D dynamics: (1) neutrino-driven convection, (2) initially neutrino-driven convection and subsequent development of the standing accretion shock instability (SASI), and (3) SASI-dominated evolution. This confirms previous 3D results of Hanke et al. and Couch and Connor. We carry out simulations with resolutions differing by up to a factor of ∼4 and demonstrate that low resolution is artificially favorable for explosion in the 3D convection-dominated case since it decreases the efficiency of energy transport to small scales. Low resolution results in higher radial convective fluxes of energy and enthalpy, more fully buoyant mass, and stronger neutrino heating. In the SASI-dominated case, lower resolution damps SASI oscillations. In the convection-dominated case, a quasi-stationary angular kinetic energy spectrum E(ℓ) develops in the heating layer. Like other 3D studies, we find E(ℓ) ∝ℓ{sup −1} in the “inertial range,” while theory and local simulations argue for E(ℓ) ∝ ℓ{sup −5/3}. We argue that current 3D simulations do not resolve the inertial range of turbulence and are affected by numerical viscosity up to the energy-containing scale, creating a “bottleneck” that prevents an efficient turbulent cascade.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

PLASMA INSTABILITIES IN THE CONTEXT OF CURRENT HELIUM SEDIMENTATION MODELS: DYNAMICAL IMPLICATIONS FOR THE ICM IN GALAXY CLUSTERS

Understanding whether Helium can sediment to the core of galaxy clusters is important for a number of problems in cosmology and astrophysics. All current models addressing this question are one-dimensional and do not account for the fact that magnetic fields can effectively channel ions and electrons, leading to anisotropic transport of momentum, heat, and particle diffusion in the weakly collisional intracluster medium (ICM). This anisotropy can lead to a wide variety of instabilities, which could be relevant for understanding the dynamics of heterogeneous media. In this paper, we consider the radial temperature and composition profiles as obtained from a state-of-the-art Helium sedimentation model and analyze its stability properties. We find that the associated radial profiles are unstable to different kinds of instabilities depending on the magnetic field orientation at all radii. The fastest growing modes are usually related to generalizations of the magnetothermal instability (MTI) and the heat-flux-driven buoyancy instability which operate in heterogeneous media. We find that the effect of sedimentation is to increase (decrease) the predicted growth rates in the inner (outer) cluster region. The unstable modes grow quickly compared to the sedimentation timescale. This suggests that the composition gradients as inferred from sedimentation models, which do not fully account for the anisotropic character of the weakly collisional environment, might not be very robust. Our results emphasize the subtleties involved in understanding the gas dynamics of the ICM and argue for the need of a comprehensive approach to address the issue of Helium sedimentation beyond current models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Gauge invariants of linearized gravity with a general background metric

In linearized gravity with distributed matter, the background metric has no generic symmetries, and decomposition of the metric perturbation into global normal modes is generally impractical. This complicates the identification of the gauge-invariant part of the perturbation, which is a concern, for example, in the theory of dispersive gravitational waves (GWs) whose energy–momentum must be gauge-invariant. Here, we propose how to identify the gauge-invariant part of the metric perturbation and the six independent gauge invariants per se for an arbitrary background metric. For the Minkowski background, the operator that projects the metric perturbation on the invariant subspace is proportional to the well-known dispersion operator of linear GWs in vacuum. For a general background, this operator is expressed in terms of the Green’s operator of the vacuum wave equation. Further, if the background is smooth, it can be found asymptotically using the inverse scale of the background metric as a small parameter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Gauge-invariant gravitational waves in matter beyond linearized gravity

Modeling the propagation of gravitational waves (GWs) in media other than vacuum is complicated by the gauge freedom of linearized gravity in that, once nonlinearities are taken into consideration, gauge artifacts can cause spurious acceleration of the matter. To eliminate these artifacts, we propose how to keep the theory of dispersive GWs gauge-invariant beyond the linear approximation and, in particular, obtain an unambiguous gauge-invariant expression for the energy–momentum of a GW in a dispersive medium. Using analytic tools from plasma physics, we propose an exactly gauge-invariant 'quasilinear' theory, in which GWs are governed by linear equations and also affect the background metric on scales large compared to their wavelength. As a corollary, the gauge-invariant geometrical optics of linear dispersive GWs in a general background is formulated. As an example, we show how the well-known properties of vacuum GWs are naturally and concisely yielded by our theory in a manifestly gauge-invariant form. We also show how the gauge invariance can be maintained within a given accuracy to an arbitrary order in the GW amplitude. These results are intended to form a physically meaningful framework for studying dispersive GWs in matter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Reaction–drift–diffusion models from master equations: application to material defects

We present a general method to produce well-conditioned continuum reaction–drift–diffusion equations directly from master equations on a discrete, periodic state space. We assume the underlying data to be kinetic Monte Carlo models (i.e. continuous-time Markov chains) produced from atomic sampling of point defects in locally periodic environments, such as perfect lattices, ordered surface structures or dislocation cores, possibly under the influence of a slowly varying external field. Our approach also applies to any discrete, periodic Markov chain. Here, the analysis identifies a previously omitted non-equilibrium drift term, present even in the absence of external forces, which can compete in magnitude with the reaction rates, thus being essential to correctly capture the kinetics. To remove fast modes which hinder time integration, we use a generalized Bloch relation to efficiently calculate the eigenspectrum of the master equation. A well conditioned continuum equation then emerges by searching for spectral gaps in the long wavelength limit, using an established kinetic clustering algorithm to define a proper reduced, Markovian state space.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Semiclassical geons as solitonic black hole remnants

We find that the end state of black hole evaporation could be represented by non-singular and without event horizon stable solitonic remnants with masses of the order the Planck scale and up to ∼ 16 units of charge. Though these objects are locally indistinguishable from spherically symmetric, massive electric (or magnetic) charges, they turn out to be sourceless geons containing a wormhole generated by the electromagnetic field. Our results are obtained by interpreting semiclassical corrections to Einstein's theory in the first-order (Palatini) formalism, which yields second-order equations and avoids the instabilities of the usual (metric) formulation of quadratic gravity. We also discuss the potential relevance of these solutions for primordial black holes and the dark matter problem.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Cosmological post-Newtonian equations from nonlinear perturbation theory

We derive the basic equations of the cosmological first-order post-Newtonian approximation from the recently formulated fully nonlinear and exact cosmological perturbation theory in Einstein's gravity. Apparently the latter, being exact, should include the former, and here we use this fact as a new derivation of the former. The complete sets of equations in both approaches are presented without fixing the temporal gauge conditions so that we can use the gauge choice as an advantage. Comparisons between the two approaches are made. Both are potentially important in handling relativistic aspects of nonlinear processes occurring in cosmological structure formation. We consider an ideal fluid and include the cosmological constant.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

GRMHD formulation of highly super-Chandrasekhar magnetized white dwarfs: stable configurations of non-spherical white dwarfs

The topic of magnetized super-Chandrasekhar white dwarfs is in the limelight, particularly in the last few years, since our proposal of their existence. By full-scale general relativistic magnetohydrodynamic (GRMHD) numerical analysis, we confirm in this work the existence of stable, highly magnetized, significantly super-Chandrasekhar white dwarfs with mass more than 3 solar mass. While a poloidal field geometry renders the white dwarfs oblate, a toroidal field makes them prolate retaining an overall quasi-spherical shape, as speculated in our earlier work. These white dwarfs are expected to serve as the progenitors of over-luminous type Ia supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

ADM Analysis of gravity models within the framework of bimetric variational formalism

Bimetric variational formalism was recently employed to construct novel bimetric gravity models. In these models an affine connection is generated by an additional tensor field which is independent of the physical metric. In this work we demonstrate how the ADM decomposition can be applied to study such models and provide some technical intermediate details. Using ADM decomposition we are able to prove that a linear model is unstable as has previously been indicated by perturbative analysis. Moreover, we show that it is also very difficult if not impossible to construct a non-linear model which is ghost-free within the framework of bimetric variational formalism. However, we demonstrate that viable models are possible along similar lines of thought. To this end, we consider a set up in which the affine connection is a variation of the Levi-Civita one. As a proof of principle we construct a gravity model with a massless scalar field obtained this way.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Is there scale-dependent bias in single-field inflation?

Scale-dependent halo bias due to local primordial non-Gaussianity provides a strong test of single-field inflation. While it is universally understood that single-field inflation predicts negligible scale-dependent bias compared to current observational uncertainties, there is still disagreement on the exact level of scale-dependent bias at a level that could strongly impact inferences made from future surveys. In this paper, we clarify this confusion and derive in various ways that there is exactly zero scale-dependent bias in single-field inflation. Much of the current confusion follows from the fact that single-field inflation does predict a mode coupling of matter perturbations at the level of f{sub NL}{sup local}; ≈ −5/3, which naively would lead to scale-dependent bias. However, we show explicitly that this mode coupling cancels out when perturbations are evaluated at a fixed physical scale rather than fixed coordinate scale. Furthermore, we show how the absence of scale-dependent bias can be derived easily in any gauge. This result can then be incorporated into a complete description of the observed galaxy clustering, including the previously studied general relativistic terms, which are important at the same level as scale-dependent bias of order f{sub NL}{sup local} ∼ 1. This description will allow us to draw unbiased conclusions about inflation from future galaxy clustering data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Cosmological singularity theorems for f ( R ) gravity theories

In the present work some generalizations of the Hawking singularity theorems in the context of f ( R ) theories are presented. The main assumptions are: the matter fields stress energy tensor satisfies the condition ( T {sub ij} −( g {sub ij} /2) T ) k {sup i} k {sup j} ≥ 0 for any generic unit time like field k {sup i} ; the scalaron takes bounded positive values during its evolution and the resulting space time is globally hyperbolic. Then, if there exist a Cauchy hyper-surface Σ for which the expansion parameter θ of the geodesic congruence emanating orthogonally from Σ satisfies some specific bounds, then the resulting space time is geodesically incomplete. Some mathematical results of reference [92] are very important for proving this. The generalized theorems presented here apply directly for some specific models such as the Hu-Sawicki or Starobinsky ones [27,38]. For other scenarios, some extra assumptions should be implemented in order to have a geodesically incomplete space time. The hypothesis considered in this text are sufficient, but not necessary. In other words, their negation does not imply that a singularity is absent.

79 ASTRONOMY AND ASTROPHYSICS↗