Engineering topics

Narayan, Ramesh

Publications and source records attributed to Narayan, Ramesh.

At least 19 records

DOE OSTI2023

Bridging Scales in Black Hole Accretion and Feedback: Magnetized Bondi Accretion in 3D GRMHD

Fueling and feedback couple supermassive black holes (SMBHs) to their host galaxies across many orders of magnitude in spatial and temporal scales, making this problem notoriously challenging to simulate. We use a multi-zone computational method based on the general relativistic magnetohydrodynamic (GRMHD) code KHARMA that allows us to span 7 orders of magnitude in spatial scale, to simulate accretion onto a non-spinning SMBH from an external medium with a Bondi radius of R B ≈ 2 × 10 5 GM•/c 2 , where M• is the SMBH mass. For the classic idealized Bondi problem, spherical gas accretion without magnetic fields, our simulation results agree very well with the general relativistic analytic solution. Meanwhile, when the accreting gas is magnetized, the SMBH magnetosphere becomes saturated with a strong magnetic field. The density profile varies as ~r -1 rather than r -3/2 and the accretion rate $\dot{M}$ is consequently suppressed by over 2 orders of magnitude below the Bondi rate $\dot{M}$ B . Here, we find continuous energy feedback from the accretion flow to the external medium at a level of $\thicksim$ 10 -2 $\dot{M}$c 2 $\thicksim$ 5 X 10 -5 $\dot{M}$ B c 2 . Energy transport across these widely disparate scales occurs via turbulent convection triggered by magnetic field reconnection near the SMBH. Thus, strong magnetic fields that accumulate on horizon scales transform the flow dynamics far from the SMBH and naturally explain observed extremely low accretion rates compared to the Bondi rate, as well as at least part of the energy feedback.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Pathways to Discovery in Astronomy and Astrophysics for the 2020s

We live in an extraordinary period of discovery in astronomy and astrophysics. Six Nobel Prizes have been awarded over the past decade alone for discoveries based on astronomical data (dark energy, gravitational waves, neutrino oscillations, the discovery of exoplanets, cosmology, and supermassive black holes). Many of the ambitious scientific visions of the 2010 New Worlds, New Horizons (Astro2010)1 decadal survey are being fulfilled, but momentum has only grown. We stand on the threshold of new endeavors that will transform not only our understanding of the universe and the processes and physical paradigms that govern it but also humanity’s place in it. This report of the Committee for a Decadal Survey on Astronomy and Astrophysics 2020 (Astro2020) proposes a broad, integrated plan for space- and ground-based astronomy and astrophysics for the decade 2023–2032.2 It also lays the foundations for further advances in the following decade. This is the seventh in a sequence of decadal survey studies in this field from the National Academies of Sciences, Engineering, and Medicine. This survey examines the program of record, providing advice on the major projects from prior surveys that are yet to be completed. It also lays out priorities for future investments driven by scientific opportunities. The recommendations in this report advance foundational activities that support the people who drive innovation and discovery, and that promote the technologies and tools needed to carry out the science. The report also recommends sustaining activities on a broad range of cost and time scales, as well as activities that enable future visionary projects by maturing them scientifically and technically. Last, the recommendations set in motion the construction of frontier facilities that will change the view and understanding of the cosmos. The survey is bounded by plausible budget scenarios based on briefings from the sponsoring agencies— the Department of Energy (DOE), the National Aeronautics and Space Administration (NASA), and the National Science Foundation (NSF). Within these bounds, the survey aims high, reflecting this time of great scientific promise and progress, with opportunities to pursue some of the most compelling scientific quests of our times.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Linear analysis of the Kelvin–Helmholtz instability in relativistic magnetized symmetric flows

ABSTRACT We study the linear stability of a planar interface separating two fluids in relative motion, focusing on the symmetric configuration where the two fluids have the same properties (density, temperature, magnetic field strength, and direction). We consider the most general case with arbitrary sound speed cs, Alfvén speed vA, and magnetic field orientation. For the instability associated with the fast mode, we find that the lower bound of unstable shear velocities is set by the requirement that the projection of the velocity on to the fluid-frame wavevector is larger than the projection of the Alfvén speed on to the same direction, i.e. shear should overcome the effect of magnetic tension. In the frame where the two fluids move in opposite directions with equal speed v, the upper bound of unstable velocities corresponds to an effective relativistic Mach number $M_{\rm re}\equiv v/v_{\rm {f}\perp }\sqrt{(1-v_{\rm {f}\perp }^2)/(1-v^2)} \cos \theta =\sqrt{2}$, where $v_{\rm {f}\perp }=[v_{\rm {A}}^2+c_{\rm s}^2(1-v_{\rm {A}}^2)]^{1/2}$ is the fast speed assuming a magnetic field perpendicular to the wavevector (here, all velocities are in units of the speed of light), and θ is the laboratory-frame angle between the flow velocity and the wavevector projection on to the shear interface. Our results have implications for shear flows in the magnetospheres of neutron stars and black holes – both for single objects and for merging binaries – where the Alfvén speed may approach the speed of light.

Astronomy & Astrophysics↗

DOE OSTI2022

Probing plasma physics with spectral index maps of accreting black holes on event horizon scales

The Event Horizon Telescope (EHT) collaboration has produced the first resolved images of the supermassive black holes at the centre of our galaxy and at the centre of the elliptical galaxy M87. As both technology and analysis pipelines improve, it will soon become possible to produce spectral index maps of black hole accretion flows on event horizon scales. Here, in this work, we predict spectral index maps of both M87* and Sgr A* by applying the general relativistic radiative transfer (GRRT) code IPOLE to a suite of general relativistic magnetohydrodynamic (GRMHD) simulations. We analytically show that the spectral index increases with increasing magnetic field strength, electron temperature, and optical depth. Consequently, spectral index maps grow more negative with increasing radius in almost all models, since all of these quantities tend to be maximized near the event horizon. Additionally, photon ring geodesics exhibit more positive spectral indices, since they sample the innermost regions of the accretion flow with the most extreme plasma conditions. Spectral index maps are sensitive to highly uncertain plasma heating prescriptions (the electron temperature and distribution function). However, if our understanding of these aspects of plasma physics can be tightened, even the spatially unresolved spectral index around 230 GHz can be used to discriminate between models. In particular, Standard and Normal Evolution (SANE) flows tend to exhibit more negative spectral indices than Magnetically Arrested Disc (MAD) flows due to differences in the characteristic magnetic field strength and temperature of emitting plasma.

(magnetohydrodynamics) MHD↗

DOE OSTI2022

Selective Dynamical Imaging of Interferometric Data

Recent developments in very long baseline interferometry (VLBI) have made it possible for the Event Horizon Telescope (EHT) to resolve the innermost accretion flows of the largest supermassive black holes on the sky. The sparse nature of the EHT’s (u, v)-coverage presents a challenge when attempting to resolve highly time-variable sources. We demonstrate that the changing (u, v)-coverage of the EHT can contain regions of time over the course of a single observation that facilitate dynamical imaging. These optimal time regions typically have projected baseline distributions that are approximately angularly isotropic and radially homogeneous. We derive a metric of coverage quality based on baseline isotropy and density that is capable of ranking array configurations by their ability to produce accurate dynamical reconstructions. We compare this metric to existing metrics in the literature and investigate their utility by performing dynamical reconstructions on synthetic data from simulated EHT observations of sources with simple orbital variability. We then use these results to make recommendations for imaging the 2017 EHT Sgr A* data set.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Reconnection-driven Particle Acceleration in Relativistic Shear Flows

Particle energization in shear flows is invoked to explain nonthermal emission from the boundaries of relativistic astrophysical jets. Yet the physics of particle injection, i.e., the mechanism that allows thermal particles to participate in shear-driven acceleration, remains unknown. With particle-in-cell simulations, we study the development of Kelvin–Helmholtz (KH) instabilities seeded by the velocity shear between a relativistic magnetically dominated electron–positron jet and a weakly magnetized electron–ion ambient plasma. Here, we show that, in their nonlinear stages, KH vortices generate kinetic-scale reconnection layers, which efficiently energize the jet particles, thus providing a first-principles mechanism for particle injection into shear-driven acceleration. Our work lends support to spine-sheath models of jet emission—with a fast core/spine surrounded by a slower sheath —and can explain the origin of radio-emitting electrons at the boundaries of relativistic jets.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2015

X-Ray Polarimetry with the Polarization Spectroscopic Telescope Array (PolSTAR)

This paper describes the Polarization Spectroscopic Telescope Array (PolSTAR), a mission proposed to NASA's 2014 Small Explorer (SMEX) announcement of opportunity. PolSTAR measures the linear polarization of 3-50 keV (requirement; goal: 2.5-70 keV) X-rays probing the behavior of matter, radiation and the very fabric of spacetime under the extreme conditions close to the event horizons of black holes, as well as in and around magnetars and neutron stars. The PolSTAR design is based on the technology developed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission launched in June 2012. In particular, it uses the same X-ray optics, extendable telescope boom, optical bench, and CdZnTe detectors as NuSTAR. The mission has the sensitivity to measure approx. 1% linear polarization fractions for X-ray sources with fluxes down to approx. 5 mCrab. This paper describes the PolSTAR design as well as the science drivers and the potential science return.

Krawczynski, Henric S.↗

NASA NTRS2011

The Extreme Spin of the Black Hole in Cygnus X-1

The compact primary in the X-ray binary Cygnus X-1 was the first black hole to be established via dynamical observations. We have recently determined accurate values for its mass and distance, and for the orbital inclination angle of the binary. Building on these results, which are based on our favored (asynchronous) dynamical model, we have measured the radius of the inner edge of the black hole s accretion disk by fitting its thermal continuum spectrum to a fully relativistic model of a thin accretion disk. Assuming that the spin axis of the black hole is aligned with the orbital angular momentum vector, we have determined that Cygnus X-1 contains a near-extreme Kerr black hole with a spin parameter a* > 0.95 (3(sigma)). For a less probable (synchronous) dynamical model, we find a. > 0.92 (3 ). In our analysis, we include the uncertainties in black hole mass, orbital inclination angle, and distance, and we also include the uncertainty in the calibration of the absolute flux via the Crab. These four sources of uncertainty totally dominate the error budget. The uncertainties introduced by the thin-disk model we employ are particularly small in this case given the extreme spin of the black hole and the disk s low luminosity.

Gou, Lijun↗

NASA NTRS2011

The Extreme Spin of the Black Hole Cygnus X-1

Remarkably, an astronomical black hole is completely described by the two numbers that specify its mass and its spin. Knowledge of spin is crucial for understanding how, for example, black holes produce relativistic jets. Recently, it has become possible to measure the spins of black holes by focusing on the very inner region of an accreting disk of hot gas orbiting the black hole. According to General Relativity (GR), this disk is truncated at an inner radius 1 that depends only on the mass and spin of the black hole. We measure the radius of the inner edge of this disk by fitting its continuum X-ray spectrum to a fully relativistic model. Using our measurement of this radius, we deduce that the spin of Cygnus X-1 exceeds 97% of the maximum value allowed by GR.

Gou, Lijun↗

NASA NTRS2005

The Extreme Spin of the Black Hole in Cygnus X-1

The compact primary in the X-ray binary Cygnus X-1 was the first black hole to be established via dynamical observatIOns. We have recently determined accurate values for its mass and distance, and for the orbital inclination angle of the binary. Building on these.results, which are based on our favored (asynchronous) dynamical model, we have measured the radius of the inner edge of the black hole's accretion disk by fitting its thermal continuum.spectrum to a fully relativistic model of a thin accretion disk. Assuming that the spin axis of the black hole is aligned with the orbital angular momentum vector, we have determined that Cygnus X-I contains a near-extreme Kerr black hole with a spin parameter a* > 0.95 (3(sigma)). For a less probable (synchronous) dynamIcal model, we find a* > 0.92 (3(sigma)). In our analysis, we include the uncertainties in black hole mass orbital inclination angle and distance, and we also include the uncertainty in the calibration of the absolute flux via the Crab. These four sources of uncertainty totally dominate the error budget. The uncertainties introduced by the thin-disk model we employ are particularly small in this case given the extreme spin of the black hole and the disk's low luminosity.

Gou, Lijun↗

NASA NTRS2003

Bondi Accretion and the Problem of the Missing Isolated Neutron Stars

A large number of neutron stars (NSs), approximately 10(exp 9), populate the Galaxy, but only a tiny fraction of them is observable during the short radio pulsar lifetime. The majority of these isolated NSs, too cold to be detectable by their own thermal emission, should be visible in X-rays as a result of accretion from the interstellar medium. The ROSAT All-Sky Survey has, however, shown that such accreting isolated NSs are very elusive: only a few tentative candidates have been identified, contrary to theoretical predictions that up to several thousand should be seen. We suggest that the fundamental reason for this discrepancy lies in the use of the standard Bondi formula to estimate the accretion rates. We compute the expected source counts using updated estimates of the pulsar velocity distribution, realistic hydrogen atmosphere spectra, and a modified expression for the Bondi accretion rate, as suggested by recent MHD simulations and supported by direct observations in the case of accretion around supermassive black holes in nearby galaxies and in our own. We find that, whereas the inclusion of atmospheric spectra partly compensates for the reduction in the counts due to the higher mean velocities of the new distribution, the modified Bondi formula dramatically suppresses the source counts. The new predictions are consistent with a null detection at the ROSAT sensitivity.

Perna, Rosalba↗

NASA NTRS2000

Self-Consistent Models of Accretion Disks

Research was carried out on several topics in the theory of astrophysical accretion flows around black holes, neutron stars and white dwarfs. The focus of our effort was the advection-dominated accretion flow (ADAF) model which the PI and his collaborators proposed and developed over the last several years. Our group completed a total of 46 papers, of which 36 are in refereed journals and 12 are in conference proceedings. All the papers have either already appeared in print or are in press.

Narayan, Ramesh↗

NASA NTRS1997

Self-Consistent Models of Accretion Disks

The investigations of advection-dominated accretion flows (ADAFs), with emphasis on applications to X-ray binaries containing black holes and neutron stars is presented. This work is now being recognized as the standard paradigm for understanding the various spectral states of black hole X-ray Binaries (BHXBs). Topics discussed include: (1) Problem in BHXBS, namely that several of these binaries have unusually large concentrations of lithium in their companion stars; (2) A novel test to show that black holes have event horizons; (3) Application of the ADAF model to the puzzling X-ray delay in the recent outburst of the BHXB, GRO J1655-40; (4) Description of the various spectral states in BHXBS; (5) Application of the ADAF model to the famous supermassive black hole at the center of our Galaxy, Sgr A(*); (6) Writing down and solving equations describing steady-state, optically thin, advection-dominated accretion onto a Kerr black hole; (7) The effect of "photon bubble" instability on radiation dominated accretion disks; and (8) Dwarf nova disks in quiescence that have rather low magnetic Reynolds number, of order 10(exp 3).

Narayan, Ramesh↗

NASA NTRS1995

Unified description of accretion flows around black holes

We provide a unified description of thermal equilibria of black hole accretion disks, including the newly discovered advection-dominated solutions. We classify the solutions on the basis of optical depth and impotance of advection cooling. We demonstrate that only four physically distinct topological types of equilibria exist. Two of the types correspond to optically thin and optically thick equilibria, while the other two types are distinguished by whether advection is negligible or dominant. A stable Shakura-Sunyaev disk exists only for accretion rates dot-M below a certain maximum. However, there is a critical viscosity parameter alpha(sub crit), which is a function of radius, such that for alpha greater than alpha(sub crit) advection-dominated solutions exist for all dot-M. Even when alpha less than alpha(sub crit), the advection-dominated solutions are available for a wide range of dot-M except for a gap around the Eddington rate. We therefore suggest that advection-dominated flows may be more common than standard thin disks in black hole systems. For certain ranges of radii and dot-M, no stable steady state solution is possible. In these cases, we suggest that limit cycle behavior may occur, leading to variability.

Chen, Xingming↗

NASA NTRS1995

Accretion disk boundary layers in cataclysmic variables. 1: Optically thick boundary layers

We develop numerical models of accretions disks in cataclysmic variables (CVs), including and emphasizing the boundary layer region where the accretion disk meets the accreting white dwarf. We confine ourselves to solutions where the boundary layer region is vertically optically thick, and find that these solutions share several common features. The angular and radial velocities of the accreting material drop rapidly in a dynamical boundary layer, which has a radial width approximately 1%-3% of the white dwarf radius. The energy dissipated in this region diffuses through the inner part of the disk and is radiated from the disk surface in a thermal boundary layer, which has a radial width comparable to the disk thickness, approximately 5%-15% of the white dwarf radius. We examine the dependence of the boundary layer structure on the mass accretion rate, the white dwarf mass and rotation rate, and the viscosity parameter alpha. We delineate the boundary between optically thick and optically thin boundary layer solutions as a function of these parameters and suggest that by means of a careful comparison with observations it may be possible to estimate alpha in CVs. We derive an expression for the total boundary layer luminosities as a function of the parameters and show that it agrees well with the luminosites of our numerical solutions. Finally, we calcuate simple blackbody continuum spectra of the boundary layer and disk emission for our solutions and compare these to soft X-ray, EUV, and He II emission-line observations of CVs. We show that, through such comparisons, it may be possible to determine the rotation rates of the accreting stars in CVs, and perhaps also the white dwarf masses and the accretion rates. The spectra are quite insensitive to alpha, so the uncertainty in this parameter does not affect such comparisons.

Popham, Robert↗

NASA NTRS1994

On the bimodal distribution of gamma-ray bursts

Kouveliotou et al. recently confirmed that gamma-ray bursts are bimodal in duration. In this paper we compute the statistical properties of the short (less than or = 2 s) and long (greater than 2 s) bursts using a method of analysis that makes no assumption regarding the location of the bursts, whether in the Galaxy or at a cosmological distance. We find the 64 ms channel on Burst and Transient Source Experiment (BATSE) to be more sensitive to short bursts and the 1024 ms channel to be more sensitive to long bursts. We show that all the currently available data are consistent with the simple hypothesis that both short and long bursts have the same spatial distribution and that within each population the sources are standard candles. The rate of short bursts per unit volume is about 40% of the rate of long bursts. Although the durations of short and long gamma-ray bursts span several orders of magnitude and the total energy of a typical short burst is smaller than that of a typical long burst by a factor of about 20, surprisingly the peak luminosities of the two kinds of bursts are equal to within a factor of about 2.

Mao, Shude↗

NASA NTRS1993

Steady flow on to a conveyor belt - Causal viscosity and shear shocks

Some hydrodynamical consequences of the adoption of a causal theory of viscosity are explored. Causality is introduced into the theory by letting the coefficient of viscosity go to zero as the flow velocity approaches a designated propagation speed for viscous signals. Consideration is given to a model of viscosity which has a finite propagation speed of shear information, and it is shown that it produces two kinds of shear shock. A 'pure shear shock' corresponds to a transition from a superviscous to a subviscous state with no discontinuity in the velocity. A 'mixed shear shock' has a shear transition occurring at the same location as a normal adiabatic or radiative shock. A generalized version of the Rankine-Hugoniot conditions for mixed shear shocks is derived, and self-consistent numerical solutions to a model 2D problem in which an axisymmetric radially infalling stream encounters a spinning star are presented.

Syer, D.↗

NASA NTRS1993

A theory of nonlocal mixing-length convection. I - The moment formalism

A flexible and potentially powerful theory of convection, based on the mixing length picture, is developed to make unbiased self-consistent predictions about overshooting and other complicated phenomena in convection. The basic formalism is set up, and the method's power is demonstrated by showing that a simplified version of the theory reproduces all the standard results of local convection. The second-order equations of the theory are considered in the limit of a steady state and vanishing third moments, and it is shown that they reproduce all the standard results of local mixing-length convection. There is a particular value of the superadiabatic gradient, below which the only possible steady state of a fluid is nonconvecting. Above this critical value, a fluid is convectively unstable. Two distinct regimes of convection, which are identified as efficient and inefficient convection, are determined.

Grossman, Scott A.↗