The field of mathematical and numerical analysis of systems of nonlinear partial differential equations involving interfaces and free boundaries is a flourishing area of research. Many such systems arise from mathematical models in material science, fluid dynamics and biology, for example phase separation in alloys, epitaxial growth, dynamics of multiphase fluids, evolution of cell membranes and in industrial processes such as crystal growth. The governing equations for the dynamics of the interfaces in many of these applications involve surface tension expressed in terms of the mean curvature and a driving force. Here the forcing terms depend on variables that are solutions of additional partial differential equations which hold either on the interface itself or in the surrounding bulk regions. Often in applications of these mathematical models, suitable performance indices and appropriate control actions have to be specified. Mathematically this leads to optimization problems with partial differential equation constraints including free boundaries. Because of the maturity of the field of computational free boundary problems it is now timely to consider such control problems.