Scheil equation


In metallurgy, the Scheil-Gulliver equation describes solute redistribution during solidification of an alloy.

Assumptions

Four key assumptions in Scheil analysis enable determination of phases present in a cast part. These assumptions are:
  1. No diffusion occurs in solid phases once they are formed
  2. Infinitely fast diffusion occurs in the liquid at all temperatures by virtue of a high diffusion coefficient, thermal convection, Marangoni convection, etc.
  3. Equilibrium exists at the solid-liquid interface, and so compositions from the phase diagram are valid
  4. Solidus and liquidus are straight segments
The fourth condition may be relaxed when numerical techniques are used, such as those used in CALPHAD software packages, though these calculations rely on calculated equilibrium phase diagrams. Calculated diagrams may include odd artifacts that influence Scheil calculations.

Derivation

The hatched areas in the figure represent the amount of solute in the solid and liquid. Considering that the total amount of solute in the system must be conserved, the areas are set equal as follows:
Since the partition coefficient is
and mass must be conserved
the mass balance may be rewritten as
Using the boundary condition
the following integration may be performed:
Integrating results in the Scheil-Gulliver equation for composition of the liquid during solidification:
or for the composition of the solid:

Applications of the Scheil equation: Calphad Tools for the Metallurgy of Solidification

Nowadays, several Calphad softwares are available - in a framework of computational thermodynamics - to simulate solidification in systems with more than two components; these have recently been defined as Calphad Tools for the Metallurgy of Solidification. In recent years, Calphad-based methodologies have reached maturity in several important fields of metallurgy, and especially in solidification-related processes such as semi-solid casting, 3d printing, and welding, to name a few. While there are important studies devoted to the progress of Calphad methodology, there is still space for a systematization of the field, which proceeds from the ability of most Calphad-based software to simulate solidification curves and includes both fundamental and applied studies on solidification, to be substantially appreciated by a wider community than today. The three applied fields mentioned above could be widened by specific successful examples of simple modeling related to the topic of this issue, with the aim of widening the application of simple and effective tools related to Calphad and Metallurgy. See also "Calphad Tools for the Metallurgy of Solidification" in an ongoing issue of an Open Journal. https://www.mdpi.com/journal/metals/special_issues/Calphad_Solidification
Given a specific chemical composition, using a software for computational thermodynamics - which might be open or commercial - the calculation of the Scheil curve is possible if a thermodynamic database is available. A good point in favour of some specific commercial softwares is that the install is easy indeed and you can use it on a windows based system - for instance with students or for self training.
One should get some open, chiefly binary, databases, one could find - after registering - at Computational Phase Diagram Database of the National Institute for Materials Science of Japan, NIMS https://cpddb.nims.go.jp/index_en.html. They are available - for free - and the collection is rather complete; in fact currently 507 binary systems are available in the thermodynamic data base format.
Some wider and more specific alloy systems partly open - with tdb compatible format - are available with minor corrections for Pandat use at Matcalc https://www.matcalc.at/index.php/databases/open-databases.

Numerical expression and numerical derivative of the Scheil curve: application to grain size on solidification and semi-solid processing

A key concept that might be used for applications is the derivative of the solid fraction fs with temperature. A numerical example using a copper zinc alloy at composition Zn 30% in weight is proposed as an example here using the opposite sign for using both temperature and its derivative in the same graph.
Kozlov and Schmid-Fetzer have calculated numerically the derivative of the Scheil curve in an open paper https://iopscience.iop.org/article/10.1088/1757-899X/27/1/012001 and applied it to the growth restriction factor Q in Al-Si-Mg-Cu alloys.
Q=fs→0
This - Calphad calculated value of numerical derivative - Q has some interesting applications in the field of metal solidification. In fact, Q reflects the phase diagram of the alloy system and its reciprocal has been found to have a relationship with grain size d on solidification, which empirically has been found in some cases to be linear:
d=a+b/Q
where a and b are constants, as illustrated with some examples from the literature for Mg and Al alloys. Before Calphad use, Q values were calculated from the conventional relationship:
Q=mc0
where m is the slope of the liquidus, c0 is the solute concentration, and k is the equilibrium distribution coefficient.
More recently some other possible correlation of Q with grain size d have been found, for instance:
d=B/1/3
where B is a constant independent of alloy composition

External Links