Abstract:
We numerically simulate the dynamical relaxation to equilibrium of a sunspot. Treating the sunspot as a thick axisymmetric flux tube in cylindrical geometry, we solve the time dependent MHD equation to examine the evolution of a sunspot towards equilibrium, starting from an arbitrary initial state.