We present a numerical study of the relation between the cosmic peculiar velocity field and the gravitational acceleration field. We show that on mildly non-linear scales (4-10 Mpc Gaussian smoothing), the distribution of the Cartesian coordinates of each of these fields is well approximated by a Gaussian. In particular, their kurtoses and negentropies are small compared to those of the velocity divergence and density fields. We find that at these scales the relation between the velocity and gravity field follows linear theory to good accuracy. Specifically, the systematic errors in velocity-velocity comparisons due to assuming the linear model do not exceed 6% in beta. To correct for them, we test various nonlinear estimators of velocity from density. We show that a slight modification of the alpha-formula proposed by Kudlicki et al. yields an estimator which is essentially unbiased and has a small variance.

We use numerical simulations to study the relation between the velocity of the Local Group and its gravitational acceleration. This relation serves as a test for the kinematic origin of the CMB dipole and as a method for estimating beta. We calculate the misalignment angle between the two vectors and compare it to the observed value for the PSCz survey. The latter value is near the upper limit of the 95% confidence interval for the angle; therefore, the nonlinear effects are unlikely to be responsible for the whole observed misalignment. We also study the relation between the amplitudes of the LG velocity and gravity vectors. In an Omega=1 Universe, the smoothed gravity of the LG turns out to be a biased low estimator of the LG (unsmoothed) velocity. In an Omega=0.3 Universe, the estimator is biased high. The discussed biases are, however, only a few per cent, thus the linear theory works to good accuracy. The gravity-based estimator of the LG velocity has also a scatter that limits the precision of the estimate of beta in the LG velocity-gravity comparisons. The random error of beta due to nonlinear effects amounts to several per cent.

We present a numerical study of the cosmic density vs. velocity divergence relation (DVDR) in the mildly non-linear regime. We approximate the dark matter as a non-relativistic pressureless fluid, and solve its equations of motion on a grid fixed in comoving coordinates. Unlike N-body schemes, this method yields directly the volume-averaged velocity field. The results of our simulations are compared with the predictions of the third-order perturbation theory (3PT) for the DVDR. We investigate both the mean `forward' relation (density in terms of velocity divergence) and the mean `inverse' relation (velocity divergence in terms of density), with emphasis on the latter. On scales larger than about 20 megaparsecs, our code recovers the predictions of 3PT remarkably well, significantly better than recent N-body simulations. On scales of a few megaparsecs, the DVDR predicted by 3PT differs slightly from the simulated one. In particular, approximating the inverse DVDR by a third-order polynomial turns out to be a poor fit. We propose a simple analytical description of the inverse relation, which works well for mildly non-linear scales.

We present results of the investigations of the statistical properties of a joint density and velocity divergence probability distribution function (PDF) in the mildly non-linear regime. For that purpose we use both perturbation theory results, extended here for a top-hat filter, and numerical simulations. In particular we derive the quantitative (complete as possible up to third order terms) and qualitative predictions for constrained averages and constrained dispersions -- which describe the nonlinearities and the stochasticity properties beyond the linear regime -- and compare them against numerical simulations. We find overall a good agreement for constrained averages; however, the agreement for constrained dispersions is only qualitative. Scaling relations for the Omega-dependence of these quantities are satisfactory reproduced. Guided by our analytical and numerical results, we finally construct a robust phenomenological description of the joint PDF in a closed analytic form. The good agreement of our formula with results of N-body simulations for a number of cosmological parameters provides a sound validation of the presented approach. Our results provide a basis for a potentially powerful tool with which it is possible to analyze galaxy survey data in order to test the gravitational instability paradigm beyond the linear regime and put useful constraints on cosmological parameters. In particular we show how the nonlinearity in the density--velocity relation can be used to break the so-called Omega-bias degeneracy in cosmic density-velocity comparisons.

We present a new Eulerian code able to follow the evolution of large-scale structures in the Universe in the weakly nonlinear regime. We compare test results with a N-body code and analytical results.