Mathematical Modeling of the Solid–Liquid Interface Propagation by the Boundary Integral Method with Nonlinear Liquidus Equation and Atomic Kinetics

Abstract

In this paper, we derive the boundary integral equation (BIE), a single integrodifferential equation governing the evolutionary behavior of the interface function, paying special attention to the nonlinear liquidus equation and atomic kinetics. As a result, the BIE is found for a thermodiffusion problem of binary melt crystallization with convection. Analyzing this equation coupled with the selection criterion for a stationary dendritic growth in the form of a parabolic cylinder, we show that nonlinear effects stemming from the liquidus equation and atomic kinetics play a decisive role. Namely, the dendrite tip velocity and diameter, respectively, become greater and lower with the increasing deviation of the liquidus equation from a linear form. In addition, the dendrite tip velocity can substantially change with variations in the power exponent of the atomic kinetics. In general, the theory under consideration describes the evolution of a curvilinear crystallization front, as well as the growth of solid phase perturbations and patterns in undercooled binary melts at local equilibrium conditions (for low and moderate Péclet numbers). In addition, our theory, combined with the unsteady selection criterion, determines the non-stationary growth rate of dendritic crystals and the diameter of their vertices.