Computational stability of an initially radial solution of a growth/dissolution problem in a nonradial implementation

We consider a free boundary problem modelling the growth/dissolution of a crystal. The aim is to investigate the following question: Does the solution to the crystal growth problem posed in two dimensions with radially symmetric initial and boundary condition evolve as a radially symmetric solution?