In this project, students explored black hole accretion physics through numerical simulations using the PLUTO hydrodynamics code. The work began with standard benchmark tests to validate the numerical setup, followed by the implementation of the Bondi accretion model, which describes spherically symmetric, adiabatic accretion onto a compact object.

By configuring the simulation geometry, boundary conditions, and physical parameters, students modeled the evolution of density, velocity, and accretion rate around the accretor. The results reproduced the expected Bondi accretion behavior, including smooth radial density gradients and transonic flow.
This project provided hands-on experience in computational astrophysics and numerical modeling of accretion processes.