search

bgmeulem / Luminet

This repo provides material for recreating the famous paper by Jean-Pierre Luminet (1979), simulating the first image of a Schwarzschild black hole.
MIT License
86 stars 6 forks source link

Impossible import #2

Closed basz closed 1 year ago

basz commented 1 year ago

from black_hole import *

Nothing leaves a black hole. This is not possible...

bgmeulem commented 1 year ago

Thank you for bringing this oversight to my attention.

Luckily, we are visualizing the accretion disk, rather than the black hole itself. While this mass is not lost forever, it is being consumed ever so slowly by the black hole by accretion effects. Assuming low angular momentum (as is our case with a Swarzschild black hole), and depending on various properties, the accretion rate may vary between $10^{-7} - 10^{-1}\ L_{Edd}/c^2$ (Bu and Yang, 2019). To account for accretion effects, you could always do something like:

acc_rate = 1e-6  # in units of Eddington luminosity per c^2
import numpy as np
if np.random.random() > accretion_rate:
    import black_hole as bh
else:
    raise UserWarning("Your import fell into the black hole :(")

[1] Bu, D.-F. and Yang, X.-H., “What is the real accretion rate on to a black hole for low-angular-momentum accretion?”, Monthly Notices of the Royal Astronomical Society, vol. 484, no. 2, pp. 1724–1734, 2019. doi:10.1093/mnras/stz050.