Thermal shock

Prevention

Methods to prevent thermal shock include:

• Minimizing the thermal gradient by changing the temperature gradually
• Increasing the thermal conductivity of the material
• Reducing the coefficient of thermal expansion of the material
• Increasing the strength of the material
• Introducing compressive stress in the material, such as in tempered glass
• Decreasing the Young's modulus of the material
• Increasing the toughness of the material through crack tip blunting or crack deflection, utilizing the process of plastic deformation, and phase transformation

Thermal shock resistance

The thermal shock resistance, ${\displaystyle \mathrm{\Delta }{T}_s}$, is the maximal temperature difference at which a material can be quenched without sustaining damage.

Strength-controlled thermal shock resistance

Thermal shock resistance is used for material selection in applications subject to rapid temperature changes. The maximum temperature jump, ${\displaystyle \mathrm{\Delta }T}$, sustainable by a material can be approximated for strength-controlled models by: $${\displaystyle B\mathrm{\Delta }T=\frac{{\sigma }_f}{\alpha E}}$$ where ${\displaystyle {\sigma }_f}$ is the failure stress (which can be yield or fracture stress), ${\displaystyle \alpha }$ is the coefficient of thermal expansion, ${\displaystyle E}$ is the Young's modulus, and ${\displaystyle B}$ is a constant depending upon the part constraint, material properties, and thickness.

$${\displaystyle B=\frac{A}{C}}$$ where ${\displaystyle C}$ is a system constrain constant dependent upon the Poisson's ratio, ${\displaystyle \nu }$, and ${\displaystyle A}$ is a non-dimensional parameter dependent upon the Biot number, ${\displaystyle \mathrm{B}\mathrm{i}}$.

$${\displaystyle C=\{\begin{array}{ll}1& \text{axial stress}\\ (1-\nu )& \text{biaxial constraint}\\ (1-2\nu )& \text{triaxial constraint}\end{array}}$$

${\displaystyle A}$ may be approximated by: $${\displaystyle A=\frac{Hh/k}{1+Hh/k}=\frac{\mathrm{B}\mathrm{i}}{1+\mathrm{B}\mathrm{i}}}$$ where ${\displaystyle H}$ is the thickness, ${\displaystyle h}$ is the heat transfer coefficient, and ${\displaystyle k}$ is the thermal conductivity.

Perfect heat transfer

If perfect heat transfer (${\displaystyle \mathrm{B}\mathrm{i}=\mathrm{\infty }}$) is assumed, the maximum heat transfer supported by the material is:$${\displaystyle \mathrm{\Delta }T=A_1\frac{{\sigma }_f}{E\alpha }}$$where,

• ${\displaystyle A_1\approx 1}$ for cold shock in plates
• ${\displaystyle A_1\approx 3.2}$ for hot shock in plates

A material index for material selection according to thermal shock resistance in the fracture stress derived perfect heat transfer case is therefore: $${\displaystyle \frac{{\sigma }_f}{E\alpha }}$$

Poor heat transfer

For cases with poor heat transfer (${\displaystyle \mathrm{B}\mathrm{i}<1}$), the maximum heat differential supported by the material is: $${\displaystyle \mathrm{\Delta }T=A_2\frac{{\sigma }_f}{E\alpha }\frac{1}{\mathrm{B}\mathrm{i}}=A_2\frac{{\sigma }_f}{E\alpha }\frac{k}{hH}}$$where,

• ${\displaystyle A_2\approx 3.2}$ for cold shock
• ${\displaystyle A_2\approx 6.5}$ for hot shock

In the poor heat transfer case, a higher thermal conductivity is beneficial for thermal shock resistance. The material index for the poor heat transfer case is often taken as: $${\displaystyle \frac{k{\sigma }_f}{E\alpha }}$$According to both the perfect and poor heat transfer models, larger temperature differentials can be tolerated for hot shock than for cold shock.

Fracture toughness controlled thermal shock resistance

In addition to thermal shock resistance defined by material fracture strength, models have also been defined within the fracture mechanics framework. Lu and Fleck produced criteria for thermal shock cracking based on fracture toughness controlled cracking. The models were based on thermal shock in ceramics (generally brittle materials). Assuming an infinite plate, and mode I cracking, the crack was predicted to start from the edge for cold shock, but the center of the plate for hot shock. Cases were divided into perfect, and poor heat transfer to further simplify the models.

Perfect heat transfer

The sustainable temperature jump decreases, with increasing convective heat transfer (and therefore larger Biot number). This is represented in the model shown below for perfect heat transfer (${\displaystyle \mathrm{B}\mathrm{i}=\mathrm{\infty }}$).

$${\displaystyle \mathrm{\Delta }T=A_3\frac{K_{Ic}}{E\alpha \sqrt{\pi H}}}$$ where ${\displaystyle K_{Ic}}$ is the mode I fracture toughness, ${\displaystyle E}$ is the Young's modulus, ${\displaystyle \alpha }$ is the thermal expansion coefficient, and ${\displaystyle H}$ is half the thickness of the plate.

• ${\displaystyle A_3\approx 4.5}$ for cold shock
• ${\displaystyle A_4\approx 5.6}$ for hot shock

A material index for material selection in the fracture mechanics derived perfect heat transfer case is therefore: $${\displaystyle \frac{K_{Ic}}{E\alpha }}$$

Poor heat transfer

For cases with poor heat transfer, the Biot number is an important factor in the sustainable temperature jump.

$${\displaystyle \mathrm{\Delta }T=A_4\frac{K_{Ic}}{E\alpha \sqrt{\pi H}}\frac{k}{hH}}$$

Critically, for poor heat transfer cases, materials with higher thermal conductivity, k, have higher thermal shock resistance. As a result, a commonly chosen material index for thermal shock resistance in the poor heat transfer case is: $${\displaystyle \frac{k{K}_{Ic}}{E\alpha }}$$

Kingery thermal shock methods

The temperature difference to initiate fracture has been described by William David Kingery to be: $${\displaystyle \mathrm{\Delta }{T}_c=S\frac{k{\sigma }^{\ast }(1-\nu )}{E\alpha }\frac{1}{h}=\frac{S}{h{R}^{{}^{\prime }}}}$$ where ${\displaystyle S}$ is a shape factor, ${\displaystyle {\sigma }^{\ast }}$ is the fracture stress, ${\displaystyle k}$ is the thermal conductivity, ${\displaystyle E}$ is the Young's modulus, ${\displaystyle \alpha }$ is the coefficient of thermal expansion, ${\displaystyle h}$ is the heat transfer coefficient, and ${\displaystyle R^{\prime }}$ is a fracture resistance parameter. The fracture resistance parameter is a common metric used to define the thermal shock tolerance of materials.

$${\displaystyle R^{\prime }=\frac{k{\sigma }^{\ast }(1-v)}{E\alpha }}$$

The formulas were derived for ceramic materials, and make the assumptions of a homogeneous body with material properties independent of temperature, but can be well applied to other brittle materials.

Testing

Thermal shock testing exposes products to alternating low and high temperatures to accelerate failures caused by temperature cycles or thermal shocks during normal use. The transition between temperature extremes occurs very rapidly, greater than 15 °C per minute.

Equipment with single or multiple chambers is typically used to perform thermal shock testing. When using single chamber thermal shock equipment, the products remain in one chamber and the chamber air temperature is rapidly cooled and heated. Some equipment uses separate hot and cold chambers with an elevator mechanism that transports the products between two or more chambers.

Glass containers can be sensitive to sudden changes in temperature. One method of testing involves rapid movement from cold to hot water baths, and back.

Examples of thermal shock failure

• Hard rocks containing ore veins such as quartzite were formerly broken down using fire-setting, which involved heating the rock face with a wood fire, then quenching with water to induce crack growth. It is described by Diodorus Siculus in Egyptian gold mines, Pliny the Elder, and Georg Agricola.
• Ice cubes placed in a glass of warm water crack by thermal shock as the exterior surface increases in temperature much faster than the interior. The outer layer expands as it warms, while the interior remains largely unchanged. This rapid change in volume between different layers creates stresses in the ice that build until the force exceeds the strength of the ice, and a crack forms, sometimes with enough force to shoot ice shards out of the container.
• Incandescent bulbs that have been running for a while have a very hot surface. Splashing cold water on them can cause the glass to shatter due to thermal shock, and the bulb to implode.
• An antique cast iron cookstove is a simple iron box on legs, with a cast iron top. A wood or coal fire is built inside the box and food is cooked on the top outer surface of the box, like a griddle. If a fire is built too hot, and then the stove is cooled by pouring water on the top surface, it will crack due to thermal shock.
• The strong gradient of temperature (due to the dousing of a fire with water) is believed to have caused the breakage of the third Tsar Bell.
• Thermal shock is a primary contributor to head gasket failure in internal combustion engines.