The evolution of a planar perturbation in a Einstein-de Sitter Universe is studied using a previously introduced Lagrangian scheme. An approximate discrete dynamical system is derived, which describes the mass agglomeration process qualitatively. Quantitative predictions for the late density profile are obtained therefrom, and validated by numerical simulations. The main result is a scaling regime for the density profile of a collapsing object of mass $M$ around cosmological coordinate $r^*$, $\rho(r)\sim \frac{M}{d}(\frac{|r-r^*|}{d})^{-{1/4}}$. The characteristic scale of the agglomeration, $d\sim (t/t_0)^{{4/9}}$, is an increasing function of cosmological time $t$. The major part of the mass hence always lies in a region with decreasing mass density. This shows that one-dimensional self-gravitating motion is not sufficient to effectively drive structure formation in an Einstein-de Sitter Universe. These results are compared with analogous investigations for the adhesion model (Burgers equation with positive viscosity), where the agglomeration is faster, and one-dimensional dynamics is effective. We further study the mutual motion of two mass agglomerations, and show that they oscillate around each other for long times, like two ``heavy particles''. Individual particles in the two agglomerations do not mix effectively on the time scale of the interagglomeration motion.

The non-linear evolution of a stratified perturbation in a three dimensional expanding Universe is considered. A general Lagrangian scheme (Q model) is introduced and numerical investigations are performed. The asymptotic contraction of the core of the agglomeration is studied. A power-law scaling is detected and an heuristic interpretation of the numerical findings is provided. An asymptotic equation for the multi-stream velocity flow is derived and it is shown to agree quantitatively with the dynamics of the Q model. The relation to the adhesion model is discussed.

The non-linear evolution of one-dimensional perturbations in a three-dimensional expanding Universe is considered. A general Lagrangian scheme is derived, and compared to two previously introduced approximate models. These models are simulated with heap-based event-driven numerical procedure, that allows for the study of large systems, averaged over many realizations of random initial conditions. One of the models is shown to be qualitatively, and, in some respects, concerning mass aggregation, quantitatively similar to the adhesion model.

A fast algorithm to study one-dimensional self-gravitating systems, and, more generally, systems that are Lagrangian integrable between collisions, is presented. The algorithm is event-driven, and uses a heap-ordered set of predicted future events. In the limit of large number of particles $N$, the operation count is dominated by the cost of reordering the heap after each event, which goes asymptotically as $\log N$. Some applications are discussed in detail.

The dynamics of a one-dimensional self-gravitating medium, with initial density almost uniform is studied. Numerical experiments are performed with ordered and with Gaussian random initial conditions. The phase space portraits are shown to be qualitatively similar to shock waves, in particular with initial conditions of Brownian type. The PDF of the mass distribution is investigated.