Physical damping, regarding the nonlinear Navier-Stokes viscous flow dynamics, refers to a tensorial turbulent dissipation term, attributed to adjacent moving macroscopic flow components. Mutual dissipation among these parts of fluid is described by a braking term in the momentum equation together with a heating term in the energy equation, both responsible of the damping of the momentum variation and of the viscous conversion of mechanical energy into heat. A macroscopic mixing scale length is currently the only characteristic length needed in the nonlinear modelling of viscous fluid dynamics describing the nonlinear eddy viscosity through the kinematic viscosity coefficient in the viscous stress tensor, without any reference to the chemical composition and to the atomic dimensions. Therefore, in this paper, we write a new formulation for the kinematic viscosity coefficient to the turbulent viscous physical dissipation in the Navier-Stokes equations, where molecular parameters are also included. Results of 2D tests are shown, where comparisons among flow structures are made on 2D shockless radial viscous transport and on 2D damping of collisional chaotic turbulence. An application to the 3D accretion disc modelling in low mass cataclysmic variables is also discussed. Consequences of the kinematic viscosity coefficient reformulation in a more strictly physical terms on the thermal conductivity coefficient for dilute gases are also discussed. The physical nature of the discussion here reported excludes any dependence by the pure mathematical aspect of the numerical modelling.

We show how to improve constraints on \Omega_m, \sigma_8, and the dark-energy equation-of-state parameter, w, obtained by Mantz et al. (2008) from measurements of the X-ray luminosity function of galaxy clusters, namely MACS, the local BCS and the REFLEX galaxy cluster samples with luminosities L> 3 \times 10^{44} erg/s in the 0.1--2.4 keV band. To this aim, we use Tinker et al. (2008) mass function instead of Jenkins et al. (2001) and the M-L relationship obtained from Del Popolo (2002) and Del Popolo et al. (2005). Using the same methods and priors of Mantz et al. (2008), we find, for a \Lambda$CDM universe, \Omega_m=0.28^{+0.05}_{-0.04} and \sigma_8=0.78^{+0.04}_{-0.05}$ while the result of Mantz et al. (2008) gives less tight constraints $\Omega_m=0.28^{+0.11}_{-0.07}$ and \sigma_8=0.78^{+0.11}_{-0.13}. In the case of a wCDM model, we find \Omega_m=0.27^{+0.07}_{-0.06}, $\sigma_8=0.81^{+0.05}_{-0.06}$ and $w=-1.3^{+0.3}_{-0.4}$, while in Mantz et al. (2008) they are again less tight \Omega_m=0.24^{+0.15}_{-0.07}, \sigma_8=0.85^{+0.13}_{-0.20} and w=-1.4^{+0.4}_{-0.7}. Combining the XLF analysis with the f_{gas}+CMB+SNIa data set results in the constraint \Omega_m=0.269 \pm 0.012, \sigma_8=0.81 \pm 0.021 and w=-1.02 \pm 0.04, to be compared with Mantz et al. (2008), \Omega_m=0.269 \pm 0.016, \sigma_8=0.82 \pm 0.03 and w=-1.02 \pm 0.06. The tightness of the last constraints obtained by Mantz et al. (2008), are fundamentally due to the tightness of the $f_{gas}$+CMB+SNIa constraints and not to their XLF analysis. Our findings, consistent with w=-1, lend additional support to the cosmological-constant model.

Current theories on planetary formation establish that giant planet formation should be contextual to their quick migration towards the central star due to the protoplanets-disc interactions on a timescale of the order of $10^5$ years, for objects of nearly 10 terrestrial masses. Such a timescale should be smaller by an order of magnitude than that of gas accretion onto the protoplanet during the hierarchical growing-up of protoplanets by collisions with other minor objects. These arguments have recently been analysed using N-body and/or fluid-dynamics codes or a mixing of them. In this work, inviscid 2D simulations are performed, using the SPH method, to study the migration of one protoplanet, to evaluate the effectiveness of the accretion disc in the protoplanet dragging towards the central star, as a function of the mass of the planet itself, of disc tangential kinematics. To this purpose, the SPH scheme is considered suitable to study the roles of turbulence, kinematic and boundary conditions, due to its intrinsic advective turbulence, especially in 2D and in 3D codes. Simulations are performed both in disc sub-Keplerian and in Keplerian kinematic conditions as a parameter study of protoplanetary migration if moderate and consistent deviations from Keplerian Kinematics occur. Our results show migration times of a few orbital periods for Earth-like planets in sub-Keplerian conditions, while for Jupiter-like planets estimates give that about $10^4$ orbital periods are needed to half the orbital size. Timescales of planet migration are strongly dependent on the relative position of the planet with respect to the shock region near the centrifugal barrier of the disc flow.

Viscous Keplerian discs become sub-Keplerian close to a black hole since they pass through sonic points before entering into it. We study the time evolution of polytropic viscous accretion discs (both in one and two dimensional flows) using Smoothed Particle Hydrodynamics. We discover that for a large region of the parameter space, when the flow viscosity parameter is less than a critical value, standing shock waves are formed. If the viscosity is very high then the shock disappears. In the intermediate viscosity the disc oscillates very significantly in viscous time-scale. Our simulations indicate that these centrifugally supported high density region close to a black hole plays an active role in the flow dynamics, and consequently, the radiation dynamics.

We present fully time-dependent solutions of accretion processes on black holes and compare them with analytical solutions. We use Smoothed Particle Hydrodynamics and Total Variation Diminishing mathods as the numerical tenchniques. Apart from steady state solutions, we also obtain time dependent results where the shock oscillates with the characteristic period of quasi-periodic oscillations (QPOs) observed in galactic and extragalactic black hole candidates. This oscillation produces 10-15% variation in the hard X-rays in galactic candidates and in soft X-rays in Extragalactic candidates.

We present results of numerical simulation of inviscid thick accretion disks and wind flows around black holes. We use Smoothed Particle Hydrodynamics (SPH) technique for this purpose. Formation of thick disks are found to be preceded by shock waves travelling away from the centrifugal barrier. For a large range of the parameter space, the travelling shock settles at a distance close to the location obtained by a one-and-a-half dimensional model of inviscid accretion disks. Occasionally, it is observed that accretion processes are aided by the formation of oblique shock waves, particularly in the initial transient phase. The post-shock region (where infall velocity suddenly becomes very small) resembles that of the usual model of thick accretion disk discussed in the literature, though they have considerable turbulence. The flow subsequently becomes supersonic before falling into the black hole. In a large number of cases which we simulate, we find the formation of strong winds which are hot and subsonic when originated from the disk surface very close to the black hole but become supersonic within a few tens of the Schwarzschild radius of the blackhole. In the case of accretion of high angular momentum flow, very little amount of matter is accreted directly onto the black hole. Most of the matter is, however, first squeezed to a small volume close to the black hole, and subsequently expands and is expelled as a strong wind. It is quite possible that this expulsion of matter and the formation of cosmic radio jets is aided by the shock heating in the inner parts of the accretion disks.