Oscillating shocks in the transonic, viscous, variable Gamma, accretion flows around black holes

We investigate the time evolution of the transonic-viscous accretion flow around a non-rotating black hole. The input parameters used for the simulation are obtained from semi-analytical solutions. This code is based on the TVD routine and correctly handles the angular momentum transport due to viscosity. The thermodynamic properties of the flow are described by a variable adiabatic index equation of state. We regenerate the inviscid and viscous steady-state solutions, including shocks, using the simulation code and compare them with the semi-analytical solutions. The angular momentum piles up across a shock due to shock-jump conditions and viscous transport of angular momentum. This will push the shock-front outward and can result in shock oscillation or a complete destabilization of shock. We study how shocks behave in the presence of viscosity. As the viscosity parameter ($\alpha$) crosses a critical value, the previously steady shock becomes time-dependent, eventually leading to oscillations. The value of this critical viscosity depends on the injection angular momentum ($\\lambda_{ou}$) and the specific energy ($\epsilon$). We estimated the posteriori bremsstrahlung and synchrotron cooling, and the net radiative output also oscillates with the frequency of the shock. We also study the variation of frequency, amplitude, and mean position of oscillation with $\alpha$. Considering a black hole with a mass of $10M_{\odot}$, we observed that the power spectrum exhibits a prominent peak at the fundamental frequency of a few to about tens of Hz, accompanied by multiple harmonics. This characteristic is frequently observed in numerous accreting black hole candidates.