The Lagrangian perturbation theory on Friedmann-Lemaitre cosmologies is compared with numerical simulations (tree-, adaptive P$^3$M- and PM codes). In previous work we have probed the large-scale performance of the Lagrangian perturbation solutions up to the third order by studying their cross- correlations with N-body simulations for various power spectra (Buchert etal 1994, Melott etal 1995, Weiss etal 1996). Thereby, spatial optimization techniques were applied by (high-frequency-)filtering of the initial power spectra. In this work the novel method of temporal optimization [Shifted-Time- Approximation (STA) and Frozen-Time-Approximation (FTA)] is investigated and used. The method is designed to compensate the native property of Lagrangian perturbation solutions to delay the collapse of structures. The method can be treated analytically. Applying the STA and FTA prescriptions a significant improvement of the performance of Lagrangian perturbation schemes (as measured by cross-correlation statistics) is observed. Using this tool we investigate a local study of special clustering models of dark matter as candidates for typical elements of the large-scale structure in the Universe, and so also focus on the performance of the perturbation solutions on smaller scales at high-spatial resolution. The models analyzed were presented in (Buchert etal 1996) and allow studying typical features of the clustering process in the weakly non-linear regime. The spatial and temporal limits of applicability of the solutions at second and third order are determined and compared with the first-order solution, which is equivalent to the ``Zel'dovich approximation'' (Zel'dovich 1970, 1973) for the type of initial data analyzed.

We present high--spatial resolution studies of the density field as predicted by Lagrangian perturbation approximations up to the third order. The first--order approximation is equivalent to the ``Zel'dovich approximation'' for the type of initial data analyzed. The study is performed for two simple models which allow studying of typical features of the clustering process in the early non--linear regime. We calculate the initial perturbation potentials as solutions of Poisson equations algebraically, and automate this calculation for a given initial random density field. The presented models may also be useful for other questions addressed to Lagrangian perturbation solutions and for the comparison of different approximation schemes. In an accompanying paper we investigate a detailed comparison with various N--body integrators using these models (Karakatsanis \& Buchert 1995). Results of the present paper include the following: 1. The collapse is accelerated significantly by the higher--order corrections confirming previous results by Moutarde \etal (1991); 2. the spatial structure of the density patterns predicted by the ``Zel'dovich approximation'' differs much from those predicted by the second-- and third--order Lagrangian approximations; 3. Second--order effects amount to internal substructures such as ``second generation'' --pancakes, --filaments and --clusters, as are also observed in N--body simulations; 4. The third--order effect gives rise to substructuring of the secondary mass--shells. The hierarchy of shell--crossing singularities that form features small high--density clumps at the intersections of caustics which we interprete as gravitational fragmentation.