Modèle intégral axisymétrique pour le calcul des pertes AC dans les bobinages supraconducteurs

RESUME L’objectif de ce memoire est de developper un modele numerique precis et rapide qui permette de determiner la distribution de courant et/ou des champs magnetiques dans les bobinages supraconducteurs, notamment pour aider a la conception de transformateurs supraconducteurs. Ce modele doit donc etre capable de determiner la distribution de la densite de courant et les pertes en regime alternatif (par la suite appelees pertes AC) dans les rubans supraconducteurs constituant la bobine simulee. En effet, il n’existe pas encore de modele numerique efficace et suffisamment rapide qui resolve les pertes AC de bobinages supraconducteurs comportant un grand nombre de tours, de l’ordre de 100 a 1000. Or depuis la decouverte des materiaux supraconducteurs dits a haute temperature critique (HT_c) en 1986, la recherche dans le domaine de la supraconductivite ne cesse d’evoluer, et leur utilisation dans les appareillages electrotechniques devient de plus en plus realisable et meme avantageux. En effet, non seulement les transformateurs supraconducteurs ont un rendement superieur a celui des transformateurs conventionnels, mais ils occupent un espace physique beaucoup plus reduit par rapport aux transformateurs conventionnels (reduction de l’ordre de 30% en taille et en masse). Cependant, de nouveaux outils de conception doivent devenir disponibles afin de diminuer les couts de fabrication des prototypes (en optimisant l’utilisation du supraconducteur) et ainsi assurer le developpement des appareillages supraconducteurs. La resistivite fortement non-lineaire des materiaux supraconducteurs et le grand nombre de tours presents dans un bobinage reel rend cette tâche difficile avec les outils de conception existants. En particulier, les logiciels d’elements finis commerciaux ne sont pas du tout adaptes a cette tâche. Dans le contexte de ce memoire, afin de trouver la distribution de courant et les pertes AC dans les rubans supraconducteurs d’une bobine, nous proposons une variation d’une methode numerique developpee par Frederic Sirois et son equipe. Le developpement du modele est fait sur Matlab et l’equation de base du modele, a partir de laquelle la resolution de la distribution du courant est effectuee, est l’equation des courants de Foucault sous forme integrale en basses frequences (i.e. le terme de courant de deplacement est neglige).---------- ABSTRACT The goal of this Master’s thesis is to develop an accurate and a fast numerical model to solve current distribution and/or magnetic fields in superconductor windings, and especially to improve the conception of superconducting transformers. The model must be able to determine current distribution and alternative current losses (called AC losses) in the superconducting tapes of the simulated winding. Indeed, efficient and fast numerical models that compute AC losses of a superconducting winding with a large number of turns (10 to 1000 turns) do not exist yet. However, since the discovery of high temperature superconductors (HTS) in 1986, research in the superconductivity field never stopped evolving and the use of these materials in power devices becomes more and more realistic and even advantageous. Indeed, superconducting transformers have a greater efficiency than conventional transformers and are much more compact (about 30% smaller in terms of size and weight compared to conventional ones). Yet, new conception tools are required in order to reduce the cost of prototypes (by optimizing the use of superconductors) and thus, ensure the development of superconducting devices. The highly non-linear resistivity of superconductor materials and the large number of turns of real windings make the task challenging with the existing conception tools. In particular, commercial finite element software packages are not adapted to this task. In this Master’s thesis, in order to find the current distribution and the AC losses in superconducting tapes of a winding, we present a variation of a numerical method developed by my research supervisor Frederic Sirois and his team. The model was developped in Matlab environment and the main equation of the model, from which the solution of the current distribution is found consist of the integral form of the eddy current equation at low frequency (i.e. displacement current neglected). All the superconducting tapes are discretized in elements in which the current distribution is uniform. In the problem formulation, the Biot-Savart integral equation is used in order to establish a direct relationship between the current density and the vector potential. Moreover, in order to compute the AC losses of a winding with a large number of turns within an acceptable time, we introduce in this Master’s thesis a new approximation method called “far tape approximation”. The idea behind it is to reduce the size of the problem by using different discretization for each tape of the winding.