Simulation on helium expansion and quench in CICC for KSTAR with moving mesh methods

A numerical code, MFEM1D, has been developed to study the helium expansion and quench propagation in CICC. Since the heat induced flow gives cause to a high heat transfer coefficient between supercritical helium and the strands, the temperature difference between the strands and helium is assumed to be very small. The strong coupling of heat transfer at the front of the normal zone generates a contact discontinuity in temperature and density. In order to obtain convergence of numerical solutions, a moving mesh method is used to capture the contact discontinuity in the short front region of the normal zone. The coupled equation is solved by the moving mesh finite element method with artificial viscosity. Details of the numerical implementation are discussed and the quench study for the KSTAR coils is performed.