Abstract. We consider a fully discrete approximation of the H gradient flow of an energy integral where the energy density is given by the sum of a nonnegative multi-well potential term and a bending energy term. The spatial discretization is based on a Fourier spectral method, which is combined with an implicit Euler time discretization. The numerical method is shown to be stable and to exhibi...