Settling of metal droplets in a terrestrial magma ocean: On the correction of the Stokes velocity

The formation of magma oceans on at least the major terrestrial planets is widely assumed even if the full accretion history and early evolution of terrestrial planets is discussed rather controversial. Various processes occur within these magma oceans, among them the settling of small iron drops and their chemical equilibration with the silicate environment. Different models were proposed to explain the differentiation within a magma ocean. In order to model a magma ocean and give constraints about settling time of iron droplets, a constant velocity is used in most models, which is calculated using Stokes’ famous formula. According to the other model parameters this assumption is invalid, since Stokes’ assumption of a creeping flow is violated. However, former models investigated the behaviour of a large number of droplets, which makes it is impossible to solve the flow around each single iron drop. We introduce a model, in which the flow around iron drops in a molten silicate environment is computed using the appropriate hydrodynamic equations. We investigate the terminal velocity of a single drop descending through a magma ocean and additionally study the effects of the presence of other obstacles in the neighbourhood. We determine a new mean velocity, which may serve as an input parameter for the existent models of magma oceans. Although we used a full fluid-dynamical approach, the velocity computed with our model is not very different from the simple Stokes case and thus proves that the former assumptions were legitimate and can be used in future as well as the new velocity presented here.