The dynamical mechanisms controlling the rheology of dense suspensions close to jamming are investigated numerically, using simplified models for the relevant dissipative forces. We show that the velocity fluctuations control the dissipation rate and therefore the effective viscosity of the suspension. These fluctuations are similar in quasi-static simulations and for finite strain rate calculations with various damping schemes. We conclude that the statistical properties of grain trajectories -in particular the critical exponent of velocity fluctuations with respect to volume fraction φ -only weakly depend on the dissipation mechanism. Rather they are determined by steric effects, which are the main driving forces in the quasistatic simulations. The critical exponent of the suspension viscosity with respect to φ can then be deduced, and is consistent with experimental data. It is generally assumed that QS accurately describe the dynamics of the true system in the limit of asymptotically small shear rateγ. However, the existence of a proper quasistatic limit remains controversial, and there is growing evidence that quasistatic flows actually correspond to a finite-size dominated regime, with a correlation length that saturates at the system size [8,14].Symmetrically, for φ < φ c , amorphous materials can flow under an infinitesimal shear stress σ and present a viscosity η diverging at φ c like η ∝ (φ c − φ) −α . Scaling laws are expected to be different in thermal (glassy) systems and in athermal sytems [16]. In the case of a suspension of non-Brownian particles, the best fit of recent experimental results give a critical exponent of α = 2.4 for volume-controlled experiments [15] and of α = 1.9 for pressure-controlled experiments [3]. The explanation of 1 (labelled ν = 1) or by the lubrication like mechanism of Eq. 3 (labelled LF). In MD, the viscosity is measured in the low shear rate regime, forγ = 10 −6 and ζ = 10 −1 . In this regime, η/ζ does not depend on the precise value of these parameters, as shown in Fig. 2a. The quantitative agreement between MD(ν = 1) and MD(LF) is coincidental. Inset: compilation of experimental data available in the literature at imposed pressure P (with φc = 0.587) and imposed volume fraction φ (with φc = 0.615) for the ratio of suspension to solvent viscosity.