“…Here, χ is the fraction of mechanical energy spent on vacancy generation, Q v f is the activation energy of vacancy formation, ζ is the neutralization effect by vacancy emiting and absorbing jogs, c j is the concentration of jogs, Ω 0 is the atomic volume, D vm is the diffusivity of vacancies andṪ is the time derivative of the temperature field. The coupled solution of the system of ordinary differential equations (ODE) of Equations (9) and (21) is carried out at elementary level with a Newton-Raphson iterative scheme, a similar solution scheme to coupled equations is presented in References [27,31,32]. The coupled solution of Equations Given the new value ofε p and ε p , and the old value X d , X s , X g , ρ i , g and c v .…”