The nonlinear dynamics of gravitational clustering in the expanding universe has been an active area of numerical research for many years. During past few years the formation, motion, and internal dynamics of galaxies and galaxy clusters became the most important and the hottest issues in cosmology and in extragalactic astronomy.
High resolution is essential for the problem.
Only the largest supercomputers are able to handle this kind of simulations. This will be done by setting initial conditions using large-scale constraints obtained from observational data. The focus will be on confronting the model with extensive observational data available for the Local Supercluster and on the history of formation of the Local Supercluster. This produces a hierarchy of irregular meshes at different levels of resolution on which the Poisson equation is solved efficiently. In order to increase the efficiency, the code has adaptive spatial, time, and mass resolution.