We explore the MRI driven dynamo in a radiatively inefficient accretion flow (RIAF) using the mean field dynamo paradigm. Using singular value decomposition (SVD) we obtain the least squares fitting dynamo coefficients $\alpha$ and $\gamma$ by comparing the time series of the turbulent electromotive force and the mean magnetic field. Our study is the first one to show the poloidal distribution of these dynamo coefficients in global accretion flow simulations. Surprisingly, we obtain a high value of the turbulent pumping coefficient $\gamma$ which transports the mean magnetic flux radially outward. This would have implications for the launching of magnetised jets which are produced efficiently in presence a large-scale poloidal magnetic field close to the compact object. We present a scenario of a truncated disc beyond the RIAF where a large scale dynamo-generated poloidal magnetic field can aid jet-launching close to the black hole. Magnitude of all the calculated coefficients decreases with radius. Meridional variations of $\alpha_{\phi \phi}$, responsible for toroidal to poloidal field conversion, is very similar to that found in shearing box simulations using the `test field' (TF) method. By estimating the relative importance of $\alpha$-effect and shear, we conclude that the MRI driven large-scale dynamo, which operates at high latitudes beyond a disc scale height, is essentially of the $\alpha-\Omega$ type.

We study the long-term evolution of non-radiative geometrically thick ($H/R\approx 0.5$) accretion flows using 3D global ideal MHD simulations and a pseudo-Newtonian gravity. We find that resolving the scale height with 42 grid points is adequate to obtain convergence with the product of quality factors $\langle \langle Q_{\theta} \rangle \rangle \langle \langle Q_{\phi} \rangle \rangle \geq 300$ and magnetic tilt angle $\theta_B \sim 13^{\circ}-14^{\circ}$. Like previous global isothermal thin disk simulations, we find stronger mean magnetic fields for the restricted azimuthal domains. Imposing periodic boundary conditions with the azimuthal extent smaller than $2\pi$ make the turbulent field at low $m$ appear as a mean field in the runs with smaller azimuthal extent. But unlike previous works, we do not find a monotonic trend in turbulence with the azimuthal extent. We conclude that the minimum azimuthal extent should be $\geq \pi/2$ to capture the flow structure, but a full $2 \pi$ extent is necessary to study the dynamo. We find an intermittent dynamo cycle, with $\alpha$-quenching playing an important role in the nonlinear saturated state. Unlike previous local studies, we find almost similar values of kinetic and magnetic $\alpha$-s, giving rise to an irregular distribution of dynamo-$\alpha$. The effects of dynamical quenching are shown explicitly for the first time in global simulations of accretion flows.

The presence of a surface at the inner boundary, such as in a neutron star or a white dwarf, allows the existence of a standing shock in steady spherical accretion.The standing shock can become unstable in 2D or 3D; this is called the {\em standing accretion shock instability} (SASI).Two mechanisms -- advective-acoustic and purely acoustic -- have been proposed to explain SASI. Using axisymmetric hydrodynamic (HD) and magnetohydrodynamic (MHD) simulations, we find that the advective-acoustic mechanism better matches the observed oscillation timescales in our simulations. The global shock oscillations present in the accretion flow can explain many observed high frequency ($\gtrsim 100$ Hz) quasi-periodic oscillations (QPOs) in X-ray binaries (XRBs). The presence of a moderately strong magnetic field adds more features to the shock oscillation pattern, giving rise to low frequency modulation in the computed light curve. This low frequency modulation can be responsible for $\sim 100$ Hz QPOs (known as hHz QPOs). We propose that the appearance of hHz QPO determines the separation of twin peak QPOs of higher frequencies.

By analyzing the time series of RXTE/PCA data, the nonlinear variabilities of compact sources have been repeatedly established. Depending on the variation in temporal classes, compact sources exhibit different nonlinear features. Sometimes they show low correlation/fractal dimension, but in other classes or intervals of time they exhibit stochastic nature. This could be because the accretion flow around a compact object is a nonlinear general relativistic system involving magnetohydrodynamics. However, the more conventional way of addressing a compact source is the analysis of its spectral state. Therefore, the question arises: What is the connection of nonlinearity to the underlying spectral properties of the flow when the nonlinear properties are related to the associated transport mechanisms describing the geometry of the flow? The present work is aimed at addressing this question. Based on the connection between observed spectral and nonlinear (time series) properties of two X-ray binaries: GRS 1915+105 and Sco X-1, we attempt to diagnose the underlying accretion modes of the sources in terms of known accretion classes, namely, Keplerian disc, slim disc, advection dominated accretion flow (ADAF) and general advective accretion flow (GAAF). We explore the possible transition of the sources from one accretion mode to others with time. We further argue that the accretion rate must play an important role in transition between these modes.

Bondi accretion assumes that there is a sink of mass at the center -- which in case of a black hole (BH) corresponds to the advection of matter across the event horizon. Other stars, such as a neutron star (NS), have surfaces and hence the infalling matter has to slow down at the surface. We study the initial value problem in which the matter distribution is uniform and at rest at $t=0$. We consider different inner boundary conditions for BHs and NSs: outflow boundary condition (mimicking mass sink at the center) valid for BHs; and {\em reflective} and steady-shock (allowing gas to cross the inner boundary at subsonic speeds) boundary conditions for NSs. We also obtain a similarity solution for cold accretion on to BHs and NSs. 1-D simulations show the formation of an outward propagating and a standing shock in NSs for reflective and steady-shock boundary conditions, respectively. Entropy is the highest at the bottom of the subsonic region for reflective boundary conditions. In 2-D this profile is convectively unstable. Using steady-shock inner boundary conditions, the flow is unstable to the standing accretion shock instability (SASI) in 2-D, which leads to global shock oscillations and may be responsible for quasi-periodic oscillations (QPOs) seen in the lightcurves of accreting systems. For steady accretion in the quiescent state, spherical accretion rate on to a NS can be suppressed by orders of magnitude compared to that on to a BH.