Free boundary and moving boundary problems, that can be used to model crystal growth or the solidification and melting of pure materials, receive growing attention in science and technology. The optimal control of these problems appear even more interesting since certain desired shapes of the boundaries improve, e.g., the material quality in the case of crystal growth. We consider the so called two-phase Stefan problem that models a solid and a liquid phase separated by a moving interface.In the work presented, we take a sharp interface model approach and define a quadratic tracking-type cost functional that penalizes the deviation of the interface from the desired state at a final time as well as the control costs. Following the "optimize-then-discretize" paradigm, we formulate a first order optimality system using the formal Lagrange approach and derive the adjoint PDE system that provides the needed gradient of the cost functional.By means of an example setup of a container with an in-and outflow of water and a cooling unit at the bottom, we illustrate how the derived formulations can be used to achieve a desired interface between the solid and the fluid phase by controlling the flow at the inlet.
Two-Phase Stefan ProblemWe consider a two-dimensional two-phase Stefan problem as in [2,3] where a domain Ω ∈ R 2 is split into the solid phase Ω s (t) and the liquid phase Ω l (t) as in figure 1. The phases are separated by the moving interface Γ I (t) which we model by the Stefan condition and which we represent as a graph, h : [0, T ] × Γ cool → R, over the lower boundary Γ cool . Further, Γ cool is a cooling boundary, the inlet Γ in is a heat source and Γ out is the outlet. The temperature distribution is characterized by the heat equation, and the system is fully coupled with Navier-Stokes equations which describe the fluid flow in Ω l (t). A detailed description of the system can be found in [1]. We solve the system numerically with finite elements and mesh movement methods with the software FEniCS. The code is available at https://gitlab.mpi-magdeburg.mpg.de/baran/Stefan_Problem_in_FEniCS.git.