Analysis of a mathematical model related to Czochralski crystal growth

This paper is devoted to the study of a stationary problem consisting of the Boussinesq approximation of the Navier–Stokes equations and two convection–diffusion equations for the temperature and concentration, respectively. The equations are considered in 3D and a velocity–pressure formulation of the Navier–Stokes equations is used. The problem is complicated by nonstandard boundary conditions for velocity on the liquid–gas interface where tangential surface forces proportional to surface gradients of temperature and concentration (Marangoni effect) and zero normal component of the velocity are assumed. The velocity field is coupled through this boundary condition and through the buoyancy term in the Navier–Stokes equations with both the temperature and concentration fields. In this paper a weak formulation of the problem is stated and the existence of a weak solution is proved. For small data, the uniqueness of the solution is established.