Optimization of a Superconducting Magnetic Energy Storage Device via a CPU-Efficient Semi-Analytical Simulation

TL;DR

This paper presents a CPU-efficient semi-analytical simulation method for optimizing superconducting magnetic energy storage devices, significantly reducing computation time while accurately modeling critical current densities and ac losses.

Contribution

The study introduces a novel Radia-based semi-analytical algorithm that outperforms finite element methods in simulating and optimizing superconducting energy storage devices.

Findings

Radia-based algorithm achieves twenty-fold CPU time reduction.

Simulations compare MgB2 and YBCO devices at 4.2 K.

Accurate ac loss calculations inform device efficiency.

Abstract

Recent advances in second generation (YBCO) high temperature superconducting wire could potentially enable the design of super high performance energy storage devices that combine the high energy density of chemical storage with the high power of superconducting magnetic storage. However, the high aspect ratio and considerable filament size of these wires requires the concomitant development of dedicated optimization methods that account for both the critical current density and ac losses in type II superconductors. Here, we report on the novel application and results of a CPU-efficient semi-analytical computer code based on the Radia 3D magnetostatics software package. Our algorithm is used to simulate and optimize the energy density of a superconducting magnetic energy storage device model, based on design constraints, such as overall size and number of coils. The rapid performance of the code is pivoted on analytical calculations of the magnetic field based on an efficient implementation of the Biot-Savart law for a large variety of 3D "base" geometries in the Radia package. The significantly-reduced CPU time and simple data input in conjunction with the consideration of realistic input variables, such as material-specific, temperature and magnetic field-dependent critical current densities have enabled the Radia-based algorithm to outperform finite element approaches by a twenty fold reduction in CPU time at the same accuracy levels. Comparative simulations of MgB2 and YBCO-based devices are performed at 4.2 K, and finally, calculations of the ac losses are computed in order to ascertain the realistic efficiency of the design configurations.