For low angular momentum axially symmetric accretion flow maintained in hydrostatic equilibrium along the vertical direction, the value of the Mach number at the critical points deviates from unity, resulting in the non-isomorphism of the critical and the sonic points. This introduces several undesirable complexities while analytically dealing with the stationary integral accretion solutions and the corresponding phase portraits. We propose that introduction of an effective dynamical sound speed may resolve the issue in an elegant way. We linear perturb the full spacetime-dependent general relativistic Euler and the continuity equations governing the structure and the dynamics of accretion disc in vertical equilibrium around Schwarzschild black holes and identify the sonic metric embedded within the stationary background flow. Such metric describes the propagation of the linear acoustic perturbation inside the accretion flow. We construct the wave equation corresponding to that acoustic perturbation and find the speed of propagation of such perturbation. We finally show that the ordinary thermodynamic sound speed should be substituted by the speed of propagation of the linear acoustic wave which has been obtained through the dynamical perturbation. Such substitution will make the value of Mach number at the critical point to be equal to unity. Use of the aforementioned effective sound speed will lead to a modified stationary disc structure where the critical and the sonic points will be identical.

This work illustrates how the formation of energy-preserving shocks for polytropic accretion and temperature-preserving shocks for isothermal accretion are influenced by various geometrical configurations of general relativistic, axisymmetric, low-angular-momentum flow in the Kerr metric. Relevant pre- and post-shock states of the accreting fluid, both dynamical and thermodynamic, have been studied comprehensively. Self-gravitational back-reaction on the metric has not been considered in the present context. An elegant eigenvalue-based analytical method has been introduced to provide qualitative descriptions of the phase-orbits corresponding to stationary transonic accretion solutions, without resorting to involved numerical schemes. Effort has been made to understand how the weakly-rotating flow behaves in close proximity of the event horizon and how such `quasi-terminal' quantities are influenced by the black hole spin for different matter geometries. Our main purpose is thus to mathematically demonstrate that for non-self-gravitating accretion, separate matter geometries, in addition to the corresponding space-time geometry, control various shock-induced phenomena observed within black hole accretion discs. This work is expected to reveal how such phenomena observed near the horizon depend on physical environment of the source harbouring a supermassive black hole at its centre. It is also expected to unfold correspondences between the dependence of accretion-related parameters on flow geometries and on black hole spin. Temperature-preserving shocks in isothermal accretion may appear bright as substantial amount of rest-mass energy of the infalling matter gets dissipated at the shock surface, and the prompt removal of such energy to maintain isothermality may power the X-ray/IR flares emitted from our Galactic centre.

In this work, we have studied the effect of higher order perturbations, particularly the second order in details, on the sonic horizon. We have considered two different schemes of perturbations which are velocity potential perturbation and mass acceleration rate perturbation. These two schemes give us qualitatively similar behaviour. We have found that the analogue gravity formalism also holds for the higher order perturbations.

For three different geometric configurations of matter and two different thermodynamic equations of state, low angular momentum, multi-transonic, axially symmetric accretion flow of matter having fractional dimension of mass distribution onto a rotating black hole has been studied by employing certain post-Newtonian pseudo-Kerr black hole potential. Such task has been accomplished mathematically by mapping the fractal nature of accreted medium onto its continuum counterpart. The difference between spin dependence of accretion dynamics of the fractal medium and the continuous medium has been highlighted. By employing a time dependent linear perturbation scheme, it has been demonstrated that accretion of matter with fractional dimension of density distribution can be considered as a natural example of classical analogue model. The corresponding acoustic surface gravity has been estimated in terms of accretion variables. The value of the surface gravity changes as the accreted matter makes a transition from its fractal nature to the corresponding continuum distribution.

We investigate the role of the spin angular momentum of astrophysical black holes in controlling the special relativistic chaotic dynamics of test particles moving under the influence of a post-Newtonian pseudo-Kerr black hole potential, along with a perturbative potential created by a asymmetrically placed (dipolar) halo. Proposing a Lyapunov-like exponent to be the effective measure of the degree of chaos observed in the system under consideration, it has been found that black hole spin anti-correlates with the degree of chaos for the aforementioned dynamics. Our findings have been explained applying the general principles of dynamical systems analysis.

Analytical studies of black hole accretion usually presumes the stability of the stationary transonic configuration. Various authors in the past several decades demonstrated the validity of such an assumption for inviscid hydrodynamic flow. Inviscid approximation is a reasonable approach for low angular momentum advection dominated flow in connection to certain supermassive black holes at the centres of the galaxies (including our own) fed from a number of stellar donors. Introduction of a weak viscosity, as a first order linear correction involving the viscosity parameter, however, may sometimes provide a more detail understanding of the observed black hole spectra. The transonic behaviour of the stationary solutions have been studied for the aforementioned quasi-viscous accretion for all possible geometric configurations of axisymmetric flow. For a sufficiently low range of the viscosity parameter, transonic solutions containing one or three critical points have been found for allowed ranges in the astrophysical parameters under the post-Newtonian pseudo-Schwarzschild scheme. With the introduction of such viscosity parameter, the only feasible critical points are of saddle and spiral types in contrary to the inviscid case where centre type points were formed instead of spiral ones. Introduction of linear perturbations on stationary flow solutions and their time evolution in both standing and radially propagating wave forms have been examined (completely analytically) in detail. Our analysis shows that similar kind of secular instability exists in all the considered disk models at large distance in the asymptotic limit, however, the model itself is valid only within a certain length scale and the disks sustain within that length scale only for at least a considerable time scale.

Multi-transonic accretion for a spinning black hole has been compared among different disc geometries within post Newtonian pseudo potential framework. The variation of stationary shock characteristics with black hole spin has been studied in details for all the disc models and compared for adiabatic as well as for isothermal scenario. The variations of surface gravity with spin for all these cases have also been investigated.

We introduce a novel formalism to investigate the role of the spin angular momentum of astrophysical black holes in influencing the behaviour of low angular momentum general relativistic accretion. We propose a metric independent analysis of axisymmetric general relativistic flow, and consequently formulate the space and time dependent equations describing the general relativistic hydrodynamic accretion flow in the Kerr metric. The associated stationary critical solutions for such flow equations are provided and the stability of the stationary transonic configuration is examined using an elegant linear perturbation technique. We examine the properties of infalling material for both prograde and retrograde accretion as a function of the Kerr parameter at extremely close proximity to the event horizon. Our formalism can be used to identify a new spectral signature of black hole spin, and has the potential of performing the black hole shadow imaging corresponding to the low angular momentum accretion flow.

Analytical treatment of black hole accretion generally presumes the stability of the stationary configuration. Various authors in the past several decades demonstrated the validity of such an assumption for inviscid hydrodynamic flow. Inviscid assumption is a reasonable approximation for low angular^M momentum advection dominated flow in connection to certain^M supermassive black holes at the^M centres of the galaxies (including our own) fed from a number of stellar donors.^M Introduction of a weak viscosity, however, may sometimes provide a more detail understanding of the observed spectrum. Recently it has been demonstrated that introduction of small amount of viscosity in the form of quasi-viscous flow makes a stationary accretion disc -- where the geometric configuration of matter is described by axisymmetric flow in hydrostatic equilibrium -- unstable. We perform similar analysis for other disc models (for all three possible geometric configurations of matter) for quasi-viscous models under the post-Newtonian scheme. We introduced perturbations on the stationary flow solution particularly in standing wave form and studied their time evolution to observe whether they grow with time. Our analysis shows that same sort of secular instability exists in other disc models too. We further argued that with sufficiently low value of viscosity in the realistic astrophysical perspective, the instability does not effectively jeopardize the stationary condition.

For axially symmetric flow of dissipationless inhomogeneous fluid onto a non rotating astrophysical black hole under the influence of a generalized pseudo-Schwarzschild gravitational potential, we investigate the influence of the background flow configuration on determining the salient features of the corresponding acoustic geometry. The acoustic horizon for the aforementioned flow structure has been located and the corresponding acoustic surface gravity $\kappa$ as well as the associated analogue Hawking temperature $T_{\rm AH}$ has been calculated {\it analytically}. The dependence of $\kappa$ on the flow geometry as well as on the nature of the back ground black hole space time (manifested through the nature of the pseudo-Schwarzschild potential used) has been discussed. Dependence of the value of $\kappa$ on various initial boundary conditions governing the dynamic and the thermodynamic properties of the background fluid flow has also been studied.

We apply the theory of algebraic polynomials to analytically study the transonic properties of general relativistic hydrodynamic axisymmetric accretion onto non-rotating astrophysical black holes. For such accretion phenomena, the conserved specific energy of the flow, which turns out to be one of the two first integrals of motion in the system studied, can be expressed as a 8$^{th}$ degree polynomial of the critical point of the flow configuration. We then construct the corresponding Sturm's chain algorithm to calculate the number of real roots lying within the astrophysically relevant domain of $\mathbb{R}$. This allows, for the first time in literature, to {\it analytically} find out the maximum number of physically acceptable solution an accretion flow with certain geometric configuration, space-time metric, and equation of state can have, and thus to investigate its multi-critical properties {\it completely analytically}, for accretion flow in which the location of the critical points can not be computed without taking recourse to the numerical scheme. This work can further be generalized to analytically calculate the maximal number of equilibrium points certain autonomous dynamical system can have in general. We also demonstrate how the transition from a mono-critical to multi-critical (or vice versa) flow configuration can be realized through the saddle-centre bifurcation phenomena using certain techniques of the catastrophe theory.

Using mathematical formalism borrowed from dynamical systems theory, a complete analytical investigation of the critical behaviour of the stationary flow configuration for the low angular momentum axisymmetric black hole accretion provides valuable insights about the nature of the phase trajectories corresponding to the transonic accretion in the steady state, without taking recourse to the explicit numerical solution commonly performed in the literature to study the multi-transonic black hole accretion disc and related astrophysical phenomena. Investigation of the accretion flow around a non rotating black hole under the influence of various pseudo-Schwarzschild potentials and forming different geometric configurations of the flow structure manifests that the general profile of the parameter space divisions describing the multi-critical accretion is roughly equivalent for various flow geometries. However, a mere variation of the polytropic index of the flow cannot map a critical solution from one flow geometry to the another, since the numerical domain of the parameter space responsible to produce multi-critical accretion does not undergo a continuous transformation in multi-dimensional parameter space. The stationary configuration used to demonstrate the aforementioned findings is shown to be stable under linear perturbation for all kind of flow geometries, black hole potentials, and the corresponding equations of state used to obtain the critical transonic solutions. Finally, the structure of the acoustic metric corresponding to the propagation of the linear perturbation studied are discussed for various flow geometries used.