Parallel-in-Time Integration of Transient Phenomena in No-Insulation Superconducting Coils Using Parareal

Erik Schnaubelt,Mariusz Wozniak,Julien Dular,Idoia Cortes Garcia,Arjan Verweij,Sebastian Schöps
2024-04-20
Abstract:High-temperature superconductors (HTS) have the potential to enable magnetic fields beyond the current limits of low-temperature superconductors in applications like accelerator magnets. However, the design of HTS-based magnets requires computationally demanding transient multi-physics simulations with highly non-linear material properties. To reduce the solution time, we propose using Parareal (PR) for parallel-in-time magneto-thermal simulation of magnets based on HTS, particularly, no-insulation coils without turn-to-turn insulation. We propose extending the classical PR method to automatically find a time partitioning using a first coarse adaptive propagator. The proposed PR method is shown to reduce the computing time when fine engineering tolerances are required despite the highly nonlinear character of the problem. The full software stack used is open-source.
Computational Engineering, Finance, and Science,Superconductivity,Computational Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: How to reduce the computational time in the simulation of the transient magnetothermal phenomena of high - temperature superconductor (HTS) non - insulated coils through the Parallel - in - Time (PinT) method. Specifically, the authors propose using the Parareal algorithm to accelerate these computationally intensive multi - physics simulations, especially for the simulation of the transient behavior of non - insulated (NI) high - temperature superconductor disk - shaped coils. ### Problem Background 1. **Application Potential of High - Temperature Superconductors** - High - temperature superconductors (HTS) have the potential to make the magnetic field strength in particle accelerator magnets exceed the limit of low - temperature superconductors. - In particular, non - insulated (NI) disk - shaped coils are of great concern in practical applications due to their high electrical and thermal stability. 2. **Complexity of Simulation** - Simulating the transient behavior of HTS coils (such as quench phenomena) involves coupled electromagnetics, thermodynamics, and even solid mechanics problems. - These problems are highly nonlinear because material properties vary greatly with temperature. - Therefore, the traditional single - core time - stepping method is very time - consuming, and an efficient parallelization method is required to accelerate the calculation. ### Proposed Solutions 1. **Parareal Algorithm** - Parareal is a parallel - time - integration method that can divide the time interval into multiple sub - intervals and solve the solution of each sub - interval in parallel, thereby accelerating the overall calculation. - The authors propose an improved Parareal method, which automatically divides the time interval through an adaptive coarse propagator to ensure load balancing and efficient calculation. 2. **Adaptive Time Partitioning** - An initial adaptive coarse propagator is used to determine the boundaries of the time interval, avoiding the problem of uneven load caused by simple equidistant time steps. - This method enables non - expert users to easily apply the Parareal algorithm. ### Experimental Verification - The authors carried out numerical experiments through a specific three - dimensional non - insulated HTS disk - shaped coil model to verify the effectiveness of the proposed method. - The results show that under high - precision requirements, when using 16 time windows, the maximum acceleration ratio can reach 2.74 times. ### Conclusions - The Parareal algorithm can significantly reduce the computational time in the simulation of the transient magnetothermal phenomena of HTS non - insulated coils. - Future work should further optimize the load balancing of fine - grained propagators and develop more efficient coarse propagators to improve the acceleration ratio. ### Formula Summary The main formulas involved in the paper include: - Weak - form magnetothermal equations: \[ (C_V \partial_t T, T')_\Omega + (\kappa \nabla T, \nabla T')_\Omega = (\rho \vec{J} \cdot \vec{J}, T')_\Omega \] \[ (\partial_t (\mu \vec{H}), \vec{H}')_\Omega + (\rho \nabla \times \vec{H}, \nabla \times \vec{H}')_\Omega = 0 \] - Power - law expression of HTS resistivity: \[ \rho_{\text{HTS}} = E_c \left( \frac{\|\vec{J}_{\text{HTS}}\|}{J_c} \right)^{n - 1} \] These formulas describe the complex magnetothermal coupling behavior in HTS coils and provide a theoretical basis for the application of the Parareal algorithm.