15.7. Species Equations

For solidification and melting of a pure substance, phase change occurs at a distinct melting temperature, . For a multicomponent mixture, however, a mushy freeze/melt zone exists between a lower solidus and an upper liquidus temperature. When a multicomponent liquid solidifies, solutes diffuse from the solid phase into the liquid phase. This effect is quantified by the partition coefficient of solute , denoted , which is the ratio of the mass fraction in the solid to that in the liquid at the interface.

Ansys Fluent computes the solidus and liquidus temperatures in a species mixture as,

(15–8)

(15–9)

where is the partition coefficient of solute , is the mass fraction of solute , and is the slope of the liquidus surface with respect to . If the value of the mass fraction exceeds the value of the eutectic mass fraction , then is clipped to when calculating the liquidus and solidus temperatures. It is assumed that the last species material of the mixture is the solvent and that the other species are the solutes.

Ansys Fluent expects that you will input a negative value for the liquidus slope of species (). If you input a positive slope for , Ansys Fluent will disregard your input and instead calculate it using the Eutectic temperature and the Eutectic mass fraction :

(15–10)

Updating the liquid fraction via Equation 15–3 can cause numerical errors and convergence difficulties in multicomponent mixtures. Instead, the liquid fraction is updated as,

(15–11)

where the superscript indicates the iteration number, is a relaxation factor with a default value of 0.9, is the cell matrix coefficient, is the time step, is the current density, is the cell volume, is the current cell temperature and is the interface temperature.

Ansys Fluent offers two models for species segregation at the micro-scale, namely the Lever rule and the Scheil rule. The former assumes infinite diffusion of the solute species in the solid, and the latter assumes zero diffusion. For the Lever rule, the interface temperature , is calculated as:

(15–12)

where is the number of species.

The Scheil rule evaluates as:

(15–13)

For information about how back diffusion (that is, a finite amount of diffusion of the solute species) can be incorporated into the formulation, see the section that follows.

For the Lever rule, species transport equations are solved for the total mass fraction of species , :

(15–14)

where is the reaction rate and is given by

(15–15)

is the velocity of the liquid and is the solid (pull) velocity. is set to zero if pull velocities are not included in the solution. The liquid velocity can be found from the average velocity (as determined by the flow equation) as

(15–16)

The liquid () and solid () mass fractions are related to each other by the partition coefficient :

(15–17)

When the Scheil model is selected, Ansys Fluent solves for as the dependent variable  [677]:

(15–18)