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.

We introduce a one-dimensional toy model of globular clusters. The model is a version of the well-known gravitational sheets system, where we take additionally into account mass and energy loss by evaporation of stars at the boundaries. Numerical integration by the "exact" event-driven dynamics is performed, for initial uniform density and Gaussian random velocities. Two distinct quasi-stationary asymptotic regimes are attained, depending on the initial energy of the system. We guess the forms of the density and velocity profiles which fit numerical data extremely well and allow to perform an independent calculation of the self-consistent gravitational potential. Some power-laws for the asymptotic number of stars and for the collision times are suggested.

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.