1. Derivation of Nabarro-Herring Creep Rate^[1]

The Nabarro-Herring creep rate, [math]\displaystyle{ \dot\epsilon_{NH} }[/math], can be derived by considering an individual rectangular grain (in a single or polycrystal) as outlined in Reference.^[1] Two opposing sides have a compressive stress applied and the other two have a tensile stress applied. The atomic volume is decreased by compression and increased by tension. Under this change, the activation energy to form a vacancy is altered by [math]\displaystyle{ \pm \sigma\Omega }[/math]. The atomic volume is [math]\displaystyle{ \Omega }[/math] and the stress is [math]\displaystyle{ \sigma }[/math]. The plus and minus indication is an increase or decrease in the activation energy due to the tensile and compressive stresses, respectively. The fraction of vacancy concentrations in the compressive ([math]\displaystyle{ N_{\nu}^C }[/math]) and tensile ([math]\displaystyle{ N_{\nu}^T }[/math]) regions are given as:

[math]\displaystyle{ N_{\nu}^C \approx \exp \left(-\frac{Q_f}{kT} \right) \exp \left(-\frac{\sigma \Omega}{kT} \right) }[/math], [math]\displaystyle{ N_{\nu}^T \approx \exp \left(-\frac{Q_f}{kT} \right) \exp \left(\frac{\sigma \Omega}{kT} \right) }[/math]

In these equations [math]\displaystyle{ Q_f }[/math] is the vacancy formation energy, [math]\displaystyle{ k }[/math] is the Boltzmann constant, and [math]\displaystyle{ T }[/math] is the absolute temperature. These vacancy concentrations are maintained at the lateral and horizontal surfaces in the grain. These net concentrations drive vacancies to the compressive regions from the tensile ones which causes grain elongation in one dimension and grain compression in the other. This is creep deformation caused by a flux of vacancy motion.

The vacancy flux, [math]\displaystyle{ J_{\nu} }[/math], associated with this motion is given by:

[math]\displaystyle{ J_{\nu} = -D_{\nu} \left(\frac{\delta N_{\nu}}{\delta x}\right) }[/math]

where [math]\displaystyle{ D_{\nu} }[/math] is the vacancy diffusivity. This is given as:

[math]\displaystyle{ D_{\nu} = D_{0\nu} \exp \left(\frac{-Q_m}{kT} \right) }[/math]

where [math]\displaystyle{ D_{0\nu} }[/math] is the diffusivity when 0 vacancies are present and [math]\displaystyle{ Q_m }[/math] is the vacancy motion energy. The term [math]\displaystyle{ \frac{\delta N_{\nu}}{\delta x} }[/math] is the vacancy concentration gradient. The term [math]\displaystyle{ \delta x }[/math] is proportional to the grain size [math]\displaystyle{ d }[/math] and [math]\displaystyle{ \delta N_{\nu}= N_{\nu}^T - N_{\nu}^C }[/math]. If we multiply [math]\displaystyle{ J_{\nu} }[/math] by [math]\displaystyle{ d^2 }[/math] we obtain:

[math]\displaystyle{ \frac{\delta V}{\delta t} \approx D_{0\nu} d \exp \left[- \frac{Q_f+Q_m}{kT} \right] \left[\exp\left(\frac{\sigma \Omega}{kT} \right)- \exp\left( -\frac{\sigma \Omega}{kT}\right) \right] }[/math]

where [math]\displaystyle{ \frac{\delta V}{\delta t} }[/math] is the volume changed per unit time during creep deformation. The change in volume can be related to the change in length along the tensile axis as [math]\displaystyle{ \delta V \approx d^2 \delta d }[/math]. Using the relationship between [math]\displaystyle{ \delta V }[/math] and [math]\displaystyle{ \delta d }[/math] the NH creep rate is given by:

[math]\displaystyle{ \dot\epsilon_{NH} = \frac{1}{d}\frac{\delta d}{\delta t} }[/math]

[math]\displaystyle{ \dot\epsilon_{NH} = \left(\frac{D_{0\nu}}{d^2}\right) \exp \left[-\frac{Q_f+Q_m}{kT} \right] \left[\exp \left(\frac{\sigma \Omega}{kT}\right) -\exp\left(-\frac{\sigma \Omega}{kT} \right)\right] }[/math]

This equation can be greatly simplified. The lattice self-diffusion coefficient is given by:

[math]\displaystyle{ D_L = D_{0 \nu} \exp \left[-\frac{Q_f+Q_m}{kT} \right] }[/math]

As previously stated, NH creep occurs at low stresses and high temperatures. In this range [math]\displaystyle{ \sigma \Omega \lt \lt kT }[/math]. For small [math]\displaystyle{ x }[/math], [math]\displaystyle{ \exp(\pm x) \approx 1 \pm x }[/math]. Thus we can re-write [math]\displaystyle{ \dot\epsilon_{NH} }[/math] as:

[math]\displaystyle{ \dot\epsilon_{NH} = A_{NH} \left(\frac{D_L}{d^2}\right)\left(\frac{\sigma \Omega}{kT} \right) }[/math]

where [math]\displaystyle{ A_{NH} }[/math] is a constant that absorbs the approximations in the derivation.

Alternatively, this can be derived in a different method where the constant [math]\displaystyle{ A_n }[/math] has different dimensions. In this case, the NH creep rate [math]\displaystyle{ \dot\epsilon }[/math] is given by:^[2]

[math]\displaystyle{ \dot\epsilon = {A_n} \frac{Db^3\sigma}{d^2kT} }[/math]

2. Comparison to Coble Creep

Coble creep is closely related to Nabarro-herring creep and is controlled by diffusion as well. Unlike Nabarro-Herring creep, mass transport occurs by diffusion along the surface of single crystals or the grain boundaries in a polycrystal.^[1] For a general expression of creep rate, the comparison between Nabarro–Herring and Coble creep can be presented as follows:^[3]

[math]\displaystyle{ \dot\epsilon = \frac{ADGb}{kT}{(\frac{\sigma}{G})}^n{(\frac{b}{d})}^p }[/math]

G is the shear modulus. The diffusivity is obtained form the tracer diffusivity, [math]\displaystyle{ D^* }[/math]. The dimensionless constant [math]\displaystyle{ A_n }[/math] depends intensively on the geometry of grains. The parameters [math]\displaystyle{ A }[/math], [math]\displaystyle{ n }[/math] and [math]\displaystyle{ p }[/math] are dependent on creep mechanisms. Nabbaro-Herring creep does not involve the motion of dislocations. It predominates over high-temperature dislocation-dependent mechanisms only at low stresses, and then only for fine-grained materials. Nabarro–Herring creep is characterized by creep rates that increase linearly with the stress and inversely with the square of grain diameter. In contrast, in Coble creep atoms diffuse along grain boundaries and the creep rate varies inversely with the cube of the grain size.^[4] Lower temperatures favor Coble creep and higher temperatures favor Herring–Nabarro creep because the activation energy for vacancy diffusion within the lattice is typically larger than that along the grain boundaries, thus lattice diffusion slows down relative to grain boundary diffusion with decreasing temperature.^[4]

3. Experimental and Theoretical Examples of Nabarro-Herring Creep

• Creep in dense, polycrystalline magnesium oxide^[5] and iron-doped polycrystalline magnesia^[6]
• Compressive creep in polycrystalline beryllium oxide^[7]
• Creep in polycrystalline Al₂O₃ that has been doped with Cr, Fe, or Ti^[8]
• Creep in dry synthetic dunite^[9] which results in trace melt and some grain growth
• Reproduced for nanopolycrystalline systems in Phase Field Crystal simulations (theory matched in terms of creep stress and grain size exponents)^[10]

The content is sourced from: https://handwiki.org/wiki/Engineering:Herring%E2%80%93Nabarro_creep