Scheil equation

Last updated

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

Contents

Solidification of a binary Cu Zn alloy, lever rule, with composition of 30% of Zinc in weight, using open version of Computherm Pandat. Cu70Zn30 solidif lever.png
Solidification of a binary Cu Zn alloy, lever rule, with composition of 30% of Zinc in weight, using open version of Computherm Pandat.
Solidification -using the Scheil rule - of a binary Cu Zn alloy, at composition Zn in weight 30% using free version of Computherm Pandat software Cu70Zn30 solidif scheil.png
Solidification -using the Scheil rule - of a binary Cu Zn alloy, at composition Zn in weight 30% using free version of Computherm Pandat software
Solidification of a binary Cu Zn alloy, with composition of 30% of Zinc in weight, using open version of Computherm Pandat. Red line is following lever rule, while Scheil model applies to the blue one Cu70Zn30 solidification.png
Solidification of a binary Cu Zn alloy, with composition of 30% of Zinc in weight, using open version of Computherm Pandat. Red line is following lever rule, while Scheil model applies to the blue one

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 (straight solidus/liquidus segments) 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 (i.e. retrograde solubility) that influence Scheil calculations.

Derivation

Scheil solidification.svg

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 (related to solute distribution) is

(determined from the phase diagram)

and mass must be conserved

the mass balance may be rewritten as

.

Using the boundary condition

at

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 (extension *.tdb), one could find - after registering - at Computational Phase Diagram Database (CPDDB) 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 (tdb) 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.

Different levels of solid fractions ( in red ) in the phase diagram of Copper and Zinc. Levels are from solid fraction fs=0.8 in steps down to 0.2 Solid Fraction Cu-Zn.png
Different levels of solid fractions ( in red ) in the phase diagram of Copper and Zinc. Levels are from solid fraction fs=0.8 in steps down to 0.2

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 (numerical) 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.

Scheil solidification of a copper zinc alloy, temperature in blue, numerical derivative of temperature with the opposite of solid fraction is red Scheilnumderivative.jpg
Scheil solidification of a copper zinc alloy, temperature in blue, numerical derivative of temperature with the opposite of solid fraction is red
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.

Application to grain size on solidification

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:

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=m*c0(k−1)

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:

where B is a constant independent of alloy composition.

Application to solidification cracking

In recent publications, prof. Sindo Kou has proposed an approach to evaluate susceptibility to solidification cracking; this approach is based on a similar approach where a quantity, , which has the dimensions of a temperature is proposed as an index of the cracking susceptibility. Again one could exploit Scheil based solidification curves to link this index to the slope of the (Scheil) solidification curve: ∂T/(∂(fS)^1/2)= ∂T/(∂(fS)*(∂(fS)^1/2)/∂(fS))= (1/2)∂T/∂(fS)*(fS)^1/2=

Application to semi-solid processing

Last but not least prof. E.J.Zoqui has summarized in his work the approach proposed by several researchers in the criteria for semi-solid processing, which involves the stability of the solid phase fs with the temperature; to process semisolid alloys the sensitivity to variation of solid fraction with temperature should be minimal: in one direction it could evolve to a difficult to deform solid, on the other to a liquid which may be difficult to shape without proper moulding. It turns out that we can express this criterion again by evaluating the slope of the solidification curve, in fact ∂(fS)/∂T should be less than a certain threshold, which is commonly accepted in the scientific and technical literature to be below 0.03 1/K. Mathematically this may be expressed by an inequation, ∂(fS)/∂T < 0.03 (1/K) - where K stands for Kelvin degrees - could be equally assumed for a rough estimate of the two main semi-solid casting processing: both rheocasting ( 0.3<fs<0.4 ) and thixoforming (0.6<fs<0.7). If one would go back just to the (numerical) and functional approaches above, one should consider the reciprocal value i.e. ∂T/∂(fS)> 33 (K)

Related Research Articles

<span class="mw-page-title-main">Solubility</span> Capacity of a substance to dissolve in a solvent in a homogeneous way

In chemistry, solubility is the ability of a substance, the solute, to form a solution with another substance, the solvent. Insolubility is the opposite property, the inability of the solute to form such a solution.

Solubility equilibrium is a type of dynamic equilibrium that exists when a chemical compound in the solid state is in chemical equilibrium with a solution of that compound. The solid may dissolve unchanged, with dissociation, or with chemical reaction with another constituent of the solution, such as acid or alkali. Each solubility equilibrium is characterized by a temperature-dependent solubility product which functions like an equilibrium constant. Solubility equilibria are important in pharmaceutical, environmental and many other scenarios.

<span class="mw-page-title-main">Eutectic system</span> Mixture with a lower melting point than its constituents

A eutectic system or eutectic mixture is a homogeneous mixture that has a melting point lower than those of the constituents. The lowest possible melting point over all of the mixing ratios of the constituents is called the eutectic temperature. On a phase diagram, the eutectic temperature is seen as the eutectic point.

<span class="mw-page-title-main">Zone melting</span> Purification process by moving a molten zone along a metal bar

Zone melting is a group of similar methods of purifying crystals, in which a narrow region of a crystal is melted, and this molten zone is moved along the crystal. The molten region melts impure solid at its forward edge and leaves a wake of purer material solidified behind it as it moves through the ingot. The impurities concentrate in the melt, and are moved to one end of the ingot. Zone refining was invented by John Desmond Bernal and further developed by William G. Pfann in Bell Labs as a method to prepare high-purity materials, mainly semiconductors, for manufacturing transistors. Its first commercial use was in germanium, refined to one atom of impurity per ten billion, but the process can be extended to virtually any solute–solvent system having an appreciable concentration difference between solid and liquid phases at equilibrium. This process is also known as the float zone process, particularly in semiconductor materials processing.

<span class="mw-page-title-main">Supercooling</span> Lowering the temperature of a liquid below its freezing point without it becoming a solid

Supercooling, also known as undercooling, is the process of lowering the temperature of a liquid below its freezing point without it becoming a solid. It is achieved in the absence of a seed crystal or nucleus around which a crystal structure can form. The supercooling of water can be achieved without any special techniques other than chemical demineralization, down to −48.3 °C (−54.9 °F). Supercooled water can occur naturally, for example in the atmosphere, animals or plants.

<span class="mw-page-title-main">Creep (deformation)</span> Tendency of a solid material to move slowly or deform permanently under mechanical stress

In materials science, creep is the tendency of a solid material to undergo slow deformation while subject to persistent mechanical stresses. It can occur as a result of long-term exposure to high levels of stress that are still below the yield strength of the material. Creep is more severe in materials that are subjected to heat for long periods and generally increase as they near their melting point.

In the physical sciences, a partition coefficient (P) or distribution coefficient (D) is the ratio of concentrations of a compound in a mixture of two immiscible solvents at equilibrium. This ratio is therefore a comparison of the solubilities of the solute in these two liquids. The partition coefficient generally refers to the concentration ratio of un-ionized species of compound, whereas the distribution coefficient refers to the concentration ratio of all species of the compound.

Precipitation hardening, also called age hardening or particle hardening, is a heat treatment technique used to increase the yield strength of malleable materials, including most structural alloys of aluminium, magnesium, nickel, titanium, and some steels, stainless steels, and duplex stainless steel. In superalloys, it is known to cause yield strength anomaly providing excellent high-temperature strength.

<span class="mw-page-title-main">Cottrell atmosphere</span> Concept in materials science

In materials science, the concept of the Cottrell atmosphere was introduced by A. H. Cottrell and B. A. Bilby in 1949 to explain how dislocations are pinned in some metals by boron, carbon, or nitrogen interstitials.

A solid solution, a term popularly used for metals, is a homogeneous mixture of two different kinds of atoms in solid state and having a single crystal structure. Many examples can be found in metallurgy, geology, and solid-state chemistry. The word "solution" is used to describe the intimate mixing of components at the atomic level and distinguishes these homogeneous materials from physical mixtures of components. Two terms are mainly associated with solid solutions – solvents and solutes, depending on the relative abundance of the atomic species.

Critical radius is the minimum particle size from which an aggregate is thermodynamically stable. In other words, it is the lowest radius formed by atoms or molecules clustering together before a new phase inclusion is viable and begins to grow. Formation of such stable nuclei is called nucleation.

<span class="mw-page-title-main">Dendrite (metal)</span>

A dendrite in metallurgy is a characteristic tree-like structure of crystals growing as molten metal solidifies, the shape produced by faster growth along energetically favourable crystallographic directions. This dendritic growth has large consequences in regard to material properties.

<span class="mw-page-title-main">Spinodal decomposition</span> Mechanism of spontaneous phase separation

Spinodal decomposition is a mechanism by which a single thermodynamic phase spontaneously separates into two phases. Decomposition occurs when there is no thermodynamic barrier to phase separation. As a result, phase separation via decomposition does not require the nucleation events resulting from thermodynamic fluctuations, which normally trigger phase separation.

Hume-Rothery rules, named after William Hume-Rothery, are a set of basic rules that describe the conditions under which an element could dissolve in a metal, forming a solid solution. There are two sets of rules; one refers to substitutional solid solutions, and the other refers to interstitial solid solutions.

In materials science, segregation is the enrichment of atoms, ions, or molecules at a microscopic region in a materials system. While the terms segregation and adsorption are essentially synonymous, in practice, segregation is often used to describe the partitioning of molecular constituents to defects from solid solutions, whereas adsorption is generally used to describe such partitioning from liquids and gases to surfaces. The molecular-level segregation discussed in this article is distinct from other types of materials phenomena that are often called segregation, such as particle segregation in granular materials, and phase separation or precipitation, wherein molecules are segregated in to macroscopic regions of different compositions. Segregation has many practical consequences, ranging from the formation of soap bubbles, to microstructural engineering in materials science, to the stabilization of colloidal suspensions.

A phase-field model is a mathematical model for solving interfacial problems. It has mainly been applied to solidification dynamics, but it has also been applied to other situations such as viscous fingering, fracture mechanics, hydrogen embrittlement, and vesicle dynamics.

<span class="mw-page-title-main">Langmuir adsorption model</span> Model describing the adsorption of a mono-layer of gas molecules on an ideal flat surface

The Langmuir adsorption model explains adsorption by assuming an adsorbate behaves as an ideal gas at isothermal conditions. According to the model, adsorption and desorption are reversible processes. This model even explains the effect of pressure i.e. at these conditions the adsorbate's partial pressure, , is related to the volume of it, V, adsorbed onto a solid adsorbent. The adsorbent, as indicated in the figure, is assumed to be an ideal solid surface composed of a series of distinct sites capable of binding the adsorbate. The adsorbate binding is treated as a chemical reaction between the adsorbate gaseous molecule and an empty sorption site, S. This reaction yields an adsorbed species with an associated equilibrium constant :

Dislocation creep is a deformation mechanism in crystalline materials. Dislocation creep involves the movement of dislocations through the crystal lattice of the material, in contrast to diffusion creep, in which diffusion is the dominant creep mechanism. It causes plastic deformation of the individual crystals, and thus the material itself.

<span class="mw-page-title-main">Enthalpy of fusion</span> Enthalpy change when a substance melts

In thermodynamics, the enthalpy of fusion of a substance, also known as (latent) heat of fusion, is the change in its enthalpy resulting from providing energy, typically heat, to a specific quantity of the substance to change its state from a solid to a liquid, at constant pressure.

Miedema's model is a semi-empirical approach for estimating the heat of formation of solid or liquid metal alloys and compounds in the framework of thermodynamic calculations for metals and minerals. It was developed by the Dutch scientist Andries Rinse Miedema while working at Philips Natuurkundig Laboratorium. It may provide or confirm basic enthalpy data needed for the calculation of phase diagrams of metals, via CALPHAD or ab initio quantum chemistry methods. For a binary system composed by elements A and B, a generic Miedema Formula could be cast as where terms Phi and nwS are explained and reported below. For a binary system the physical picture could be simplified by considering a relatively simpler function of the difference of these three physical parameters resulting in a more complex form

References