Numerical quenchback in thermofluid simulations of superconducting magnets

摘要

One of the most important thermofluid processes encountered in internally cooled superconducting magnets is that of quenching. Numerical simulation of the quench propagation involves accurately modelling a moving boundary layer at the quench front. Due to the highly non-linear nature of the quench process, slightest numerical errors can rapidly grow to unacceptable limits. The quench propagation in such a non-converged solution exhibits a very rapid propagation velocity which resembles a 'quenchback' effect. Hence, the term 'Numerical Quenchback' is used to characterize a numerically unstable solution of the governing quench model. This paper presents the underlying physical phenomena that causes a numerical discretization scheme to have error terms that increase exponentially with time, causing the numerical quenchback effect. Specifically, by analytically solving the equivalent differential equation of the numerical scheme, we are able to obtain closed-form relations for the error terms associated with the propagation velocity. This allows us to define error criteria on the space and time steps used in the simulation. The reliability of the error criteria is proven by detailed convergence studies of the quench process.