Cosmological linear perturbation theory predicts that the peculiar velocity
, V(x), and matter overdensity, delta(x), at the same point x will be stati
stically independent quantities, as long as the initial density fluctuation
s are random Gaussian distributed. However, nonlinear gravitational effects
might change the situation. Using a framework of second-order perturbation
theory and the Edgeworth expansion method, we study the local density depe
ndence of bulk velocity dispersion that is coarse-grained at a weakly nonli
near scale. For a typical cold dark matter (CDM) model, the first nonlinear
correction of this constrained bulk velocity dispersion amounts to similar
to 0.3 delta (Gaussian smoothing), at a weakly nonlinear scale, with a ver
y weak dependence on cosmological parameters. We also compare our analytica
l prediction with published numerical results given at nonlinear regimes.