Motion of a continuous fluid can be decomposed into an "incompressible" rearrangement, which preserves the volume of each infinitesimal fluid element, and a gradient map that transfers fluid elements in a way unaffected by any pressure or elasticity (the polar decomposition of Y. Brenier). The Euler equation describes a system whose kinematics is dominated by the incompressible rearrangement. The opposite limit, in which the incompressible component is negligible, corresponds to the Zel'dovich approximation, a model of motion of self-gravitating fluid in cosmology. We present a method of approximate reconstruction of the large-scale proper motions of matter in the Universe from the present-day mass density field. The method is based on recovering the corresponding gradient transfer map. We discuss its algorithmics, tests of the method against mock cosmological catalogues, and its application to observational data, which result in tight constraints on the mean mass density Omega_m and age of the Universe.

We consider exact solutions to the problem of ballistic aggregation in a flow of adhesive particles, providing a model for large-scale structure formation in cosmology within the framework of the Zel'dovich approximation. Two different explicit variational constructions, suggested for this problem by (A. Shnirelman, 1986) and (E Weinan et al, 1996) in the case of flat symmetry, are shown to be equivalent. Explicit solutions are constructed for flows of cylindric and spherical symmetry. Two explicit counterexamples, showing that natural generalizations of both variational constructions fail to provide solutions for nonsymmetric multidimensional flows, are presented.

We show that the deterministic past history of the Universe can be uniquely reconstructed from the knowledge of the present mass density field, the latter being inferred from the 3D distribution of luminous matter, assumed to be tracing the distribution of dark matter up to a known bias. Reconstruction ceases to be unique below those scales -- a few Mpc -- where multi-streaming becomes significant. Above 6 Mpc/h we propose and implement an effective Monge-Ampere-Kantorovich method of unique reconstruction. At such scales the Zel'dovich approximation is well satisfied and reconstruction becomes an instance of optimal mass transportation, a problem which goes back to Monge (1781). After discretization into N point masses one obtains an assignment problem that can be handled by effective algorithms with not more than cubic time complexity in N and reasonable CPU time requirements. Testing against N-body cosmological simulations gives over 60% of exactly reconstructed points. We apply several interrelated tools from optimization theory that were not used in cosmological reconstruction before, such as the Monge-Ampere equation, its relation to the mass transportation problem, the Kantorovich duality and the auction algorithm for optimal assignment. Self-contained discussion of relevant notions and techniques is provided.

A method for the reconstruction of the primordial density fluctuation field is presented. Various previous approaches to this problem rendered {\it non-unique} solutions. Here, it is demonstrated that the initial positions of dark matter fluid elements, under the hypothesis that their displacement is the gradient of a convex potential, can be reconstructed uniquely. In our approach, the cosmological reconstruction problem is reformulated as an assignment problem in optimisation theory. When tested against numerical simulations, our scheme yields excellent reconstruction on scales larger than a few megaparsecs.