In a novel approach to studying viscous accretion flows, viscosity has been introduced as a perturbative effect, involving a first-order correction in the $\alpha$-viscosity parameter. This method reduces the problem of solving a second-order nonlinear differential equation (Navier-Stokes equation) to that of an effective first-order equation. Viscosity breaks down the invariance of the equilibrium conditions for stationary inflow and outflow solutions, and distinguishes accretion from wind. Under a dynamical systems classification, the only feasible critical points of this "quasi-viscous" flow are saddle points and spirals. A linearised and radially propagating time-dependent perturbation gives rise to secular instability on large spatial scales of the disc. Further, on these same length scales, the velocity evolution equation of the quasi-viscous flow has been transformed to bear a formal closeness with Schr\"odinger's equation with a repulsive potential. Compatible with the transport of angular momentum to the outer regions of the disc, a viscosity-limited length scale has been defined for the full spatial extent over which the accretion process would be viable.

A first-order correction in the $\alpha$-viscosity parameter of Shakura and Sunyaev has been introduced in the standard inviscid and thin accretion disc. A linearised time-dependent perturbative study of the stationary solutions of this "quasi-viscous" disc leads to the development of a secular instability on large spatial scales. This qualitative feature is equally manifest for two different types of perturbative treatment -- a standing wave on subsonic scales, as well as a radially propagating wave. Stability of the flow is restored when viscosity disappears.

The stationary background flow in the spherically symmetric infall of a compressible fluid, coupled to the space-time defined by the static Schwarzschild metric, has been subjected to linearized perturbations. The perturbative procedure is based on the continuity condition and it shows that the coupling of the flow with the geometry of space-time brings about greater stability for the flow, to the extent that the amplitude of the perturbation, treated as a standing wave, decays in time, as opposed to the amplitude remaining constant in the Newtonian limit. In qualitative terms this situation simulates the effect of a dissipative mechanism in the classical Bondi accretion flow, defined in the Newtonian construct of space and time. As a result of this approach it becomes impossible to define an acoustic metric for a conserved spherically symmetric flow, described within the framework of Schwarzschild geometry. In keeping with this view, the perturbation, considered separately as a high-frequency travelling wave, also has its amplitude reduced.

Transonicity in a spherically symmetric accreting system has been considered in both the stationary and the dynamic regimes. The stationary flow, set up as a dynamical system, has been shown to be greatly unstable to even the minutest possible deviation in the boundary condition for transonicity. With the help of a simple analytical model, and some numerical modelling, it has then been argued that the flow indeed becomes transonic and stable, when the evolution of the flow is followed through time. The time-dependent approach also shows that there is a remarkable closeness between an equation of motion for a perturbation in the flow, and the metric of an analog acoustic black hole.

For inviscid, rotational accretion flows driven by a general pseudo-Newtonian potential on to a Schwarzschild black hole, the only possible fixed points are saddle points and centre-type points. For the specific choice of the Newtonian potential, the flow has only two critical points, of which the outer one is a saddle point while the inner one is a centre-type point. A restrictive upper bound is imposed on the admissible range of values of the angular momentum of sub-Keplerian flows through a saddle point. These flows are very unstable to any deviation from a necessarily precise boundary condition. The difficulties against the physical realisability of a solution passing through the saddle point have been addressed through a temporal evolution of the flow, which gives a non-perturbative mechanism for selecting a transonic solution passing through the saddle point. An equation of motion for a real-time perturbation about the stationary flows reveals a very close correspondence with the metric of an acoustic black hole, which is also an indication of the primacy of transonicity.

Modelling the flow in a thin accretion disc like a dynamical system, we analyse the nature of the critical points of the steady solutions of the flow. For the simple inviscid disc there are two critical points, with the outer one being a saddle point and the inner one a centre type point. For the weakly viscous disc, there are four possible critical points, of which the outermost is a saddle point, while the next inner one is a spiral. Coupling the nature of the critical points with the outer boundary condition of the flow, gives us a complete understanding of all the important physical features of the flow solutions in the subsonic regions of the disc. In the inviscid disc, the physical realisability of the transonic solution passing through the saddle point is addressed by considering a temporal evolution of the flow, which is a very likely non-perturbative mechanism for selecting the transonic inflow solution from among a host of other possible stable solutions. For the weakly viscous disc, while a linearised time-dependent perturbation imposed on the steady mass inflow rate causes instability, the same perturbative analysis reveals that for the inviscid disc, there is a very close correspondence between the equation for the propagation of the perturbation and the metric of an acoustic black hole. Compatible with the transport of angular momentum to the outer regions of the disc, a viscosity-limited length scale is defined for the spatial extent of the inward rotational drift of matter.

The inviscid and thin accretion disc is a simple and well understood model system in accretion studies. In this work, modelling such a disc like a dynamical system, we analyse the nature of the fixed points of the stationary solutions of the flow. We show that of the two fixed points, one is a saddle and the other is a centre type point. We then demonstrate, using a simple but analogous mathematical model, that a temporal evolution of the flow is a very likely non-perturbative mechanism for the selection of an inflow solution that passes through the saddle type critical point.

The role of turbulence in a spherically symmetric accreting system has been studied on very large spatial scales of the system. This is also a highly subsonic flow region and here the accreting fluid has been treated as nearly incompressible. It has been shown here that the coupling of the mean flow and the turbulent fluctuations, gives rise to a scaling relation for an effective "turbulent viscosity". This in turn leads to a dynamic scaling for sound propagation in the accretion process. As a consequence of this scaling, the sonic horizon of the transonic inflow solution is shifted inwards, in comparison with the inviscid flow.

The spherically symmetric stationary transonic (Bondi) flow is considered a classic example of an accretion flow. This flow, however, is along a separatrix, which is usually not physically realizable. We demonstrate, using a pedagogical example, that it is the dynamics which selects the transonic flow.