  We compare N-body simulations performed in MOND with analogs in Newtonian
gravity with dark matter (DM). We have developed a code which solves the
Poisson equation in both gravity models. It is a grid solver using adaptive
mesh refinement techniques, allowing us to study isolated galaxies as well as
interacting galaxies. Galaxies in MOND are found to form bars faster and
stronger than in the DM model. In Newton dynamics, it is difficult to reproduce
the observed high frequency of strong bars, while MOND appears to fit better
the observations. Galaxy interactions and mergers, such as the Antennae, are
also simulated with Newton and MOND dynamics. In the latter, dynamical friction
is much weaker, and merging time-scales are longer. The formation of tidal
dwarf galaxies in tidal tails are also compared in MOND and Newton+DM models.
