AddThis Social Bookmark Button


AddThis Social Bookmark Button




One of the fundamental goals in metallurgical engineering is to be able to control the final physical, mechanical, and/or chemical properties of an alloy. To achieve this goal, one must understand the relationship between the alloy's chemical composition, the processing conditions, the resulting microstructure, and the alloy's final properties. During processing, and often also during use, most alloys undergo one or more phase transformations, which are best understood through the use of phase diagrams. Experimentally determined phase diagrams, however, are usually available for binary systems only (e.g., Co-Cr, Fe-C, Fe-Co, Fe-Cr, Fe-Mn, Fe-Mo, Fe-Ni, Fe-Ti, Fe-V, Ni-Cr, Ti-Al, Ti-B, Ti-N, etc), to some extent for ternary systems (e.g., Fe-Cr-C, W-Co-C, W-Fe-C, W-Ni-C, etc), and very rarely for higher-order systems. This is where the CALPHAD (CALculation of PHAse Diagrams) method comes in.

In the early seventies computational thermodynamics started with the pioneering work of Dr. Larry Kaufman, evolving to what is known today as the CALPHAD method. The method is based on the concept of deriving the thermodynamic functions of a system from all available experimental data. The thermodynamic functions are expressed as polynomials of temperature and chemical composition. The numerical values of the polynomial coefficients are obtained using numerical optimization techniques.

The CALPHAD method is based on the fact that a phase diagram is a representation of the thermodynamic properties of a system. Thus, if the thermodynamic properties are known, it would be possible to calculate the multi-component phase diagrams. Thermodynamic descriptions of lower-order systems (e.g., the Gibbs energy of each phase) are combined to extrapolate higher-order systems.

The Gibbs energy of a phase is described by a model that contains a relatively small number of experimentally optimized variable coefficients. Examples of experimental information used include melting and other transformation temperatures, solubilities, as well as thermodynamic properties such as heat capacities, enthalpies of formation, and chemical potentials.

For pure elements and stoichiometric compounds, the following model is most commonly used:

Gm - HmSER = a + b·T +c·T·ln(T) + Σdi·Ti


Gm - HmSER is the Gibbs energy relative to a standard element reference state (SER),

HmSER is the enthalpy of the element in its stable state at the temperature of 298.15 Kelvin and the pressure of 105 Pascal (1 bar), and

a, b, c, and di are the model parameters.

For multi-component solution phases, the following expression for the Gibbs energy is used:

G = G° + idGmix + xsGmix


G° is the Gibbs energy due to the mechanical mixing of the constituents of the phase,

idGmix is the ideal mixing contribution, and

xsGmix is the excess Gibbs energy of mixing (the non-ideal mixing contribution).

If the sub-lattice model is used to describe solution phases, then the G°, idGmix, and xsGmix of an A-B binary system with two sub-lattices, (A,B)p(A,B)q, can be expressed as follows:

G° = yIA·yIIA·G°A:A + yIA·yIIB·G°A:B + yIB·yIIA·G°B:A + yIB·yIIB·G°B:B

idGmix = p·R·T·[yIA·ln(yIA) + yIB·ln(yIB)] + q·R·T·[yIIA·ln(yIIA) + yIIB·ln(yIIB)]

xsGmix = yIA·yIB·[yIIA·Σk=0 LkA,B:A·(yIA - yIB)k + yIIB·Σk=0 LkA,B:B·(yIA - yIB)k] +

    + yIIA·yIIB·[yIA·Σk=0 LkA:A,B·(yIIA - yIIB)k + yIB·Σk=0 LkB:A,B·(yIIA - yIIB)k]


yI and yII are the site fractions of components A and B in the first and second sub-lattices, respectively,

G°A:A, G°A:B, G°B:A, and G°B:B is the Gibbs energy of ApAq, ApBq, BpAq, BpBq compounds, respectively,

LkA,B:A and LkA,B:B are the kth-order interaction parameter between component A and B in the first sub-lattice, and

LkA:A,B and LkB:A,B are the kth-order interaction parameter between component A and B in the second sub-lattice.

In the expressions for interaction parameters, a comma separates interacting components in the same sub-lattice, and a colon separates components that occupy different sub-lattices. The equations for G°, idGmix, and xsGmix can be generalized for multi-component and multi-sub-lattice phases.

If a phase in a multi-component solution is described with a single sub-lattice model, then the G°, idGmix, and xsGmix contributions to the Gibbs energy can be expressed as follows:

G° = Σi ci·G°i

idGmix = R·T·Σi ci·ln(ci)

xsGmix = Σi Σj>i ci·cj·Σk Lki,j·(ci - cj)k


ci and cj is the mole fraction of species i and j, respectively, and

Lki,j is a binary interaction parameter between species i and j.

The binary interaction parameter Lki,j is dependant on the value of k. When the value of k is equal to zero or one, the equation for xsGmix becomes regular or sub-regular, respectively.

Consulting Services

To cover the costs of running this site, we accept consulting assignments to perform customer tailored Thermo-Calc and DICTRA calculations. If we cannot solve your problem, we will help you find at least one organization which has the right human and computational resources to address your specific needs.

We offer a money back guaranty for our consulting services if you are not satisfied. Drop us a line; our e-mail address is: