  We report the first simulations of the Faraday instability using the full
three-dimensional Navier-Stokes equations in domains much larger than the
characteristic wavelength of the pattern. We use a massively parallel code
based on a hybrid Front-Tracking/Level-set algorithm for Lagrangian tracking of
arbitrarily deformable phase interfaces. Simulations performed in rectangular
and cylindrical domains yield complex patterns. In particular, a
superlattice-like pattern similar to those of [Douady & Fauve, Europhys. Lett.
6, 221-226 (1988); Douady, J. Fluid Mech. 221, 383-409 (1990)] is observed. The
pattern consists of the superposition of two square superlattices. We
conjecture that such patterns are widespread if the square container is large
compared to the critical wavelength. In the cylinder, pentagonal cells near the
outer wall allow a square-wave pattern to be accommodated in the center.
