More Like this

1. Accretion–ablation mechanics

   https://doi.org/10.1098/rsta.2022.0373  

   Pradhan, Satya Prakash; Yavari, Arash (December 2023, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences)

   In this paper, we formulate a geometric nonlinear theory of the mechanics of accreting–ablating bodies. This is a generalization of the theory of accretion mechanics of Sozio & Yavari (Sozio & Yavari 2019J. Nonlinear Sci.29, 1813–1863 (doi:10.1007/s00332-019-09531-w)). More specifically, we are interested in large deformation analysis of bodies that undergo a continuous and simultaneous accretion and ablation on their boundaries while under external loads. In this formulation, the natural configuration of an accreting–ablating body is a time-dependent Riemannian $3$-manifold with a metric that is an unknowna prioriand is determined after solving the accretion–ablation initial-boundary-value problem. In addition to the time of attachment map, we introduce a time of detachment map that along with the time of attachment map, and the accretion and ablation velocities, describes the time-dependent reference configuration of the body. The kinematics, material manifold, material metric, constitutive equations and the balance laws are discussed in detail. As a concrete example and application of the geometric theory, we analyse a thick hollow circular cylinder made of an arbitrary incompressible isotropic material that is under a finite time-dependent extension while undergoing continuous ablation on its inner cylinder boundary and accretion on its outer cylinder boundary. The state of deformation and stress during the accretion–ablation process, and the residual stretch and stress after the completion of the accretion–ablation process, are computed. This article is part of the theme issue ‘Foundational issues, analysis and geometry in continuum mechanics’. 

   more » « less
2. Bridging Scales in Black Hole Accretion and Feedback: Magnetized Bondi Accretion in 3D GRMHD

   https://doi.org/10.3847/2041-8213/ad1048  

   Cho 조, Hyerin 혜린; Prather, Ben S.; Narayan, Ramesh; Natarajan, Priyamvada; Su, Kung-Yi; Ricarte, Angelo; Chatterjee, Koushik (December 2023, The Astrophysical Journal Letters)

   Abstract 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 ofR_B≈ 2 × 10⁵GM_•/c², whereM_•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⁻¹rather thanr^−3/2and the accretion rate $\dot{M}$is consequently suppressed by over 2 orders of magnitude below the Bondi rate ${\dot{M}}_{\mathrm{B}}$. We find continuous energy feedback from the accretion flow to the external medium at a level of $\sim {10}^{-2}\dot{M}{c}^2\sim 5\times {10}^{-5}{\dot{M}}_{\mathrm{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. 

   more » « less
3. Core segregation during pebble accretion

   https://doi.org/10.1016/j.epsl.2022.117537  

   Olson, Peter; Sharp, Zachary; Garai, Susmita (June 2022, Earth and Planetary Science Letters)
4. GRRMHD simulations of MAD accretion discs declining from super-Eddington to sub-Eddington accretion rates

   https://doi.org/10.1093/mnras/stac3330  

   Curd, Brandon; Narayan, Ramesh (November 2022, Monthly Notices of the Royal Astronomical Society)

   ABSTRACT We present two general relativistic radiation magnetohydrodynamics (GRRMHD) simulations of magnetically arrested discs (MADs) around non-spinning (a* = 0) and spinning (a* = 0.9) supermassive black holes (BHs). In each simulation, the mass accretion rate is decreased with time such that we sample Eddington-scaled rates over the range $$3 \gtrsim \dot{M}/\dot{M}_{\rm {Edd}}\gtrsim 0.3$$. For the non-spinning BH model, the total and radiative efficiencies increase as the accretion rate decreases, varying over the range $$\eta _{\rm {tot}}\sim 9\!-\!16{{\ \rm per\ cent}}$$ and $$\eta _{\rm {rad}}\sim 6{-}12{{\ \rm per\ cent}}$$, respectively. This model shows very little jet activity. In contrast, the spinning BH model has a strong relativistic jet powered by spin energy extracted from the BH. The jet power declines with accretion rate such that $$\eta _{\rm {jet}}\sim 18{-}39{{\ \rm per\ cent}}$$ while the total and radiative efficiencies are $$\eta _{\rm {tot}}\sim 64{-}100{{\ \rm per\ cent}}$$ and $$\eta _{\rm {rad}}\sim 45{-}79{{\ \rm per\ cent}}$$, respectively. We confirm that mildly sub-Eddington discs can extract substantial power from a spinning BH, provided they are in the MAD state. The jet profile out to $$100\, GM/c^2$$ is roughly parabolic with a power-law index of k ≈ 0.43−0.53 during the sub-Eddington evolution. Both models show significant variability in the outgoing radiation which is likely associated with episodes of magnetic flux eruptions. The a* = 0.9 model shows semiregular variations with a period of $$\sim 2000\, GM/c^3$$ over the final $$\sim 10\, 000\, GM/c^3$$ of the simulation, which suggests that magnetic flux eruptions may be an important source of quasi-periodic variability. For the simulated accretion rates, the a* = 0 model is spinning up while the a* = 0.9 model is spinning down. Spinup–spindown equilibrium of the BH will likely be achieved at 0.5 < a*, eq < 0.6, assuming continuous accretion in the MAD state. 

   more » « less
5. Nuclear burning in collapsar accretion discs

   https://doi.org/10.1093/mnras/staa3002  

   Zenati, Yossef; Siegel, Daniel M; Metzger, Brian D; Perets, Hagai B (October 2020, Monthly Notices of the Royal Astronomical Society)

   ABSTRACT The core collapse of massive, rapidly-rotating stars are thought to be the progenitors of long-duration gamma-ray bursts (GRB) and their associated hyperenergetic supernovae (SNe). At early times after the collapse, relatively low angular momentum material from the infalling stellar envelope will circularize into an accretion disc located just outside the black hole horizon, resulting in high accretion rates necessary to power a GRB jet. Temperatures in the disc mid-plane at these small radii are sufficiently high to dissociate nuclei, while outflows from the disc can be neutron-rich and may synthesize r-process nuclei. However, at later times, and for high progenitor angular momentum, the outer layers of the stellar envelope can circularize at larger radii ≳ 107 cm, where nuclear reactions can take place in the disc mid-plane (e.g. 4He + 16O → 20Ne + γ). Here we explore the effects of nuclear burning on collapsar accretion discs and their outflows by means of hydrodynamical α-viscosity torus simulations coupled to a 19-isotope nuclear reaction network, which are designed to mimic the late infall epochs in collapsar evolution when the viscous time of the torus has become comparable to the envelope fall-back time. Our results address several key questions, such as the conditions for quiescent burning and accretion versus detonation and the generation of 56Ni in disc outflows, which we show could contribute significantly to powering GRB SNe. Being located in the slowest, innermost layers of the ejecta, the latter could provide the radioactive heating source necessary to make the spectral signatures of r-process elements visible in late-time GRB-SNe spectra. 

   more » « less