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.

The rotation curve (RC) of our Galaxy, the Milky Way, is constructed starting from its very inner regions (few hundred pc) out to a large Galactocentric distance of $\sim 200$ kpc using kinematical data on a variety of tracer objects moving in the gravitational potential of the Galaxy, without assuming any theoretical models of the visible and dark matter components of the Galaxy. We study the effect on the RC due to the uncertainties in the values of the Galactic Constants (GCs) $R_\odot$ and $V_\odot$ (these being the sun's distance from and circular rotation speed around the Galactic center, respectively) and the velocity anisotropy parameter $\beta$ of the halo tracer objects used for deriving the RC at large Galactocentric distances. The resulting RC in the disk region is found to depend significantly on the choice of the GCs, while the dominant uncertainty in the RC at large distances beyond the stellar disk comes from the uncertainty in the value of $\beta$. In general we find that the mean RC steadily declines at distances beyond $\sim 60$ kpc, independently of the value of $\beta$. Also, at a given radius, the circular speed is lower for larger values of $\beta$ (i.e., for more radially biased velocity anisotropy). Considering that the largest possible value of $\beta$ is unity, which corresponds to stellar orbits being purely radial, our results for the case of $\beta=1$ give a lower limit to the total mass of the Galaxy within $\sim 200$ kpc, $M(200 {\rm kpc}) \gsim (6.8\pm4.1) \times 10^{11} M_\odot$, independently of any model of the dark matter halo of the Galaxy.

The velocity distribution function (VDF) of the hypothetical Weakly Interacting Massive Particles (WIMPs), currently the most favored candidate for the Dark Matter (DM) in the Galaxy, is determined directly from the circular speed ("rotation") curve data of the Galaxy assuming isotropic VDF. This is done by "inverting" --- using Eddington's method --- the Navarro-Frenk-White universal density profile of the DM halo of the Galaxy, the parameters of which are determined, by using Markov Chain Monte Carlo (MCMC) technique, from a recently compiled set of observational data on the Galaxy's rotation curve extended to distances well beyond the visible edge of the disk of the Galaxy. The derived most-likely local isotropic VDF strongly differs from the Maxwellian form assumed in the "Standard Halo Model" (SHM) customarily used in the analysis of the results of WIMP direct-detection experiments. A parametrized (non-Maxwellian) form of the derived most-likely local VDF is given. The astrophysical "g-factor" that determines the effect of the WIMP VDF on the expected event rate in a direct-detection experiment can be lower for the derived most-likely VDF than that for the best Maxwellian fit to it by as much two orders of magnitude at the lowest WIMP mass threshold of a typical experiment.

Direct detection of Weakly Interacting Massive Particle (WIMP) candidates of Dark Matter (DM) is studied within the context of a self-consistent truncated isothermal model of the finite-size dark halo of the Galaxy based on the "King model" of the phase space distribution function of collisionless DM particles. Our halo model takes into account the modifications of the phase-space structure of the halo due to the gravitational influence of the observed visible matter in a self-consistent manner. The parameters of the halo model are determined by a fit to a recently determined circular rotation curve of the Galaxy that extends up to $\sim$ 60 kpc. Unlike in the Standard Halo Model (SHM) customarily used in the analysis of the results of WIMP direct detection experiments, the velocity distribution of the WIMPs in our model is non-Maxwellian with a cut-off at a maximum velocity that is self-consistently determined by the model itself. For our halo model that provides the best fit to the rotation curve data, the 90% C.L. upper limit on the WIMP-nucleon spin-independent cross section from the recent results of the CDMS-II experiment, for example, is $\sim 5.3\times10^{-8}\pb$ at a WIMP mass of $\sim$ 71 GeV. We also find, using the original 2-bin annual modulation amplitude data of the DAMA experiment, that there exists a range of small WIMP masses, typically $\sim$ 2 -- 16 GeV, within which DAMA collaboration's claimed annual modulation signal purportedly due to WIMPs is compatible with the null results of other experiments. These results strengthen the possibility of low-mass ($\lsim 10\gev$) WIMPs as a candidate for dark matter as indicated by several earlier studies performed within the context of the SHM. A more rigorous analysis using DAMA bins over smaller intervals should be able to better constrain the "DAMA regions" in the WIMP parameter space within the context of our model.

For inviscid, rotational accretion flows, both isothermal and polytropic, a simple dynamical systems analysis of the critical points has given a very accurate mathematical scheme to understand the nature of these points, for {\em any} pseudo-potential by which the flow may be driven on to a Schwarzschild black hole. This allows for a complete classification of the critical points for a wide range of flow parameters, and shows that the only possible critical points for this kind of flow are saddle points and centre-type points. A restrictive upper bound on the angular momentum of critical solutions has been established. A time-dependent perturbative study reveals that the form of the perturbation equation, for both isothermal and polytropic flows, is invariant under the choice of any particular pseudo-potential. Under generically true outer boundary conditions, the inviscid flow has been shown to be stable under an adiabatic and radially propagating perturbtion. The perturbation equation has also served the dual purpose of enabling an understanding of the acoustic geometry for inviscid and rotational flows.