3D relativistic accretion disk simulation around a 10 solar mass black hole, built with Python and Vispy (GPU-accelerated OpenGL).

• Paczyński-Wiita pseudo-Newtonian potential — approximates general relativistic effects (correct ISCO, perihelion precession, unstable orbits)
• Keplerian orbital initialization with radial density profile and vertical scale height
• 800 particles integrated with scipy.integrate.odeint

• Custom colormap (dark red → orange → yellow → white → blue-white)
• Doppler beaming — approaching side of the disk appears brighter
• Variable particle size — inner particles are larger/more intense
• Fading trails — transparent at the tail, semi-opaque at the head
• Solid black hole shadow with photon ring at ~2.6 r_s
• Positions normalized to Schwarzschild radii (r_s)
• Interactive camera — drag to rotate, scroll to zoom, middle-click to pan

python3 -m venv .venv
source .venv/bin/activate
pip install -r requirements.txt

python3 accretion_disk.py

The animation window will open with a slowly rotating camera. Drag to rotate, scroll to zoom, middle-click to pan. Close the window to exit.