Discrete symmetries are commonplace in field theoretical models but pose a severe problem for cosmology since they lead to the formation of domain walls during spontaneous symmetry breaking in the early universe. However if one of the vacuua is favoured over the others, either energetically, or because of initial conditions, it will eventually come to dominate the universe. Using numerical methods, we study the evolution of the domain wall network for a variety of field configurations in two and three dimensions and quantify the rate at which the walls disappear. Good agreement is found with a recent analytic estimate of the termination of the scaling regime of the wall network.

We study the cosmology of the Supersymmetric Standard Model augmented by a gauge singlet to solve the $\mu$-problem and describe the evolution of the domain walls which are created during electroweak symmetry breaking due to the discrete $Z_{3}$ symmetry in this model. The usual assumption, that non-renormalizable terms induced by gravity (which explicitly break this symmetry) may cause the walls to collapse on a cosmologically safe timescale, is reconsidered. Such terms are constrained by considerations of primordial nucleosynthesis, and also by the fact that by not respecting the $Z_{3}$ symmetry they induce divergences which destabilise the hierarchy and reintroduce the $\mu$--problem. We find that, even when the K\"ahler potential is `non-minimal' (i.e. when the hidden sector couples directly to the visible) the model is either ruled out cosmologically or suffers from a naturalness problem.