Researcher profile

Francesco Grilli

Francesco Grilli contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2022arXiv

2D Modeling of HTS Coils with $T$-$A$ Formulation: How to Handle Different Coupling Scenarios

Numerical models based on the finite-element method (FEM) are popular tools for investigating the macroscopic electromagnetic behavior of high-temperature superconductor (HTS) applications. This article explains how to use the $T$-$A$ formulation for modeling HTS coils in 2D with different coupling scenarios between the turns. First we consider a racetrack coil wound from one piece of superconducting tape. Then we consider a coil obtained by winding a cable composed of different HTS tapes. In the latter case, the tape turns are either electrically connected along their entire length or just at the two ends of the coil: in the model, these two different types of electrical connection are implemented with the help of the electrical circuit module. The current density distributions and the AC losses of the coils in the different coupling scenarios are compared and discussed. The limits of applicability of the presented approach are pointed out. The model is developed for the straight section of racetrack coils, but can be easily adapted to axisymmetric geometries.

preprint2022arXiv

3D homogenization of the T-A formulation for the analysis of coils with complex geometries

The modeling and analysis of superconducting coils is an essential task in the design stage of most devices based on high-temperature superconductors (HTS). These calculations allow verifying basic estimations and assumptions, proposing improvements, and computing quantities that are not easy to calculate with an analytical approach. For instance, the estimation of losses in HTS is fundamental during the design stage since losses can strongly influence the cooling system requirements and operating temperature. Typically, 2D finite element analysis is used to calculate AC losses in HTS, due to the lack of analytical solutions that can accurately represent complex operating conditions such as AC transport current and AC external applied magnetic field in coils. These 2D models are usually a representation of an infinitely long arrangement. Therefore, they cannot be used to analyze end effects and complex 3D configurations. In this publication, we use the homogenization of the T-A formulation in 3D for the analysis of superconducting coils with complex geometries where a 2D approach can not provide accurate analyses and verification of assumptions. The modeling methodology allows an easier implementation in commercial software (COMSOL Multiphysics) in comparison with the currently available 3D H homogenization, despite the complexity of the geometry. This methodology is first validated with a racetrack coil (benchmark case) by comparing the results with the well-established H formulation. Then, the electromagnetic behavior of coils with more complex geometries is analyzed.

preprint2021arXiv

A new benchmark problem for electromagnetic modelling of superconductors: the high-Tc superconducting dynamo

The high-Tc superconducting (HTS) dynamo is a promising device that can inject large DC supercurrents into a closed superconducting circuit. This is particularly attractive to energise HTS coils in NMR/MRI magnets and superconducting rotating machines without the need for connection to a power supply via current leads. It is only very recently that quantitatively accurate, predictive models have been developed which are capable of analysing HTS dynamos and explain their underlying physical mechanism. In this work, we propose to use the HTS dynamo as a new benchmark problem for the HTS modelling community. The benchmark geometry consists of a permanent magnet rotating past a stationary HTS coated-conductor wire in the open-circuit configuration, assuming for simplicity the 2D (infinitely long) case. Despite this geometric simplicity the solution is complex, comprising time-varying spatially-inhomogeneous currents and fields throughout the superconducting volume. In this work, this benchmark problem has been implemented using several different methods, including H-formulation-based methods, coupled H-A and T-A formulations, the Minimum Electromagnetic Entropy Production method, and integral equation and volume integral equation-based equivalent circuit methods. Each of these approaches show excellent qualitative and quantitative agreement for the open-circuit equivalent instantaneous voltage and the cumulative time-averaged equivalent voltage, as well as the current density and electric field distributions within the HTS wire at key positions during the magnet transit. A critical analysis and comparison of each of the modelling frameworks is presented, based on the following key metrics: number of mesh elements in the HTS wire, total number of mesh elements in the model, number of degrees of freedom (DOFs), tolerance settings and the approximate time taken per cycle for each model.

preprint2020arXiv

Electromagnetic modeling of large-scale high-temperature superconductor systems

The development of the high-temperature superconductors (HTS) conductors has allowed the development of diverse superconductor devices. Some of these devises, like the power generators and high-field magnets, are classified as large-scale HTS systems, because they are made of hundreds or thousands of turns. Mathematical models are required to address the analysis of these kind of systems. This task cannot be done by means of analytical models, because they are limited to the analysis of simple assemblies. The finite-element models using the H formulation have been extensively used during the last years. Nevertheless, the use of H formulation models to analyze large-scale systems is hindered by the excessive computational load. The recently proposed T-A formulation models have allowed building more efficient models for systems made of HTS tapes. Additionally, the homogenization and multi-scaling methods have been successfully applied in conjunction with the H and T-A formulations, these simplification methods allows reducing the required computational resources. In this article a new simplification method, called densification, is proposed. The strategies emerging from the combined use of the formulations and the simplification methods already mentioned are extensively explored, and the comprehensive validation and comparison of all the resulting strategies is presented.

preprint2020arXiv

Electromagnetic Modeling of Superconductors with Commercial Software: Possibilities with Two Vector Potential-Based Formulations

In recent years, the $H$ formulation of Maxwell's equation has become the de facto standard for simulating the time-dependent electromagnetic behavior of superconducting applications with commercial software. However, there are cases where other formulations are desirable, for example for modeling superconducting turns in electrical machines or situations where the superconductor is better described by the critical state than by a power-law resistivity. In order to accurately and efficiently handle those situations, here we consider two published approaches based on the magnetic vector potential: the $T$-$A$ formulation of Maxwell's equations (with power-law resistivity) and Campbell's implementation of the critical state model. In this contribution, we extend the $T$-$A$ formulation to thick conductors so that large coils with different coupling scenarios between the turns can be considered. We also revise Campbell's model and discuss it in terms of its ability to calculate AC losses: in particular, we investigate the dependence of the calculated AC losses on the frequency of the AC excitation and the possibility of using quick one-step (instead of full cycle) simulations to calculate the AC losses.

preprint2020arXiv

Implementation of the H-$ϕ$ formulation in COMSOL Multiphysics for simulating the magnetization of bulk superconductors and comparison with the H-formulation

The H-formulation, used abundantly for the simulation of high temperature superconductors, has shown to be a very versatile and easily implementable way of modeling electromagnetic phenomena involving superconducting materials. However, the simulation of a full vector field in current-free domains unnecessarily adds degrees of freedom to the model, thereby increasing computation times. In this contribution, we implement the well known H-$ϕ$ formulation in COMSOL Multiphysics in order to compare the numerical performance of the H and H-$ϕ$ formulations in the context of computing the magnetization of bulk superconductors. We show that the H-$ϕ$ formulation can reduce the number of degrees of freedom and computation times by nearly a factor of two for a given relative error. The accuracy of the magnetic fields obtained with both formulations are demonstrated to be similar. The computational benefits of the H-$ϕ$ formulation are shown to far outweigh the added complexity of its implementation, especially in 3-D. Finally, we identify the ideal element orders for both H and H-$ϕ$ formulations to be quartic in 2-D and cubic in 3-D, corresponding to the highest element orders implementable in COMSOL.

preprint2020arXiv

Review of the AC Loss Computation for HTS using the H-formulation

This article presents a review of the finite element method (FEM) model based on the $H$ formulation of Maxwell's equations used to calculate AC losses in high temperature superconductor (HTS) tapes, cables and windings for different applications. This model, which uses the components of the magnetic field as state variables, has been gaining a great popularity and has been in use in tens of research groups around the world. This contribution first reviews the equations on which the model is based and their implementation in finite element method programs for different cases, such 2D longitudinal and axis-symmetric geometries, 3D geometries. Modeling strategies to tackle large number of HTS tapes, such as multi-scale and homogenization methods, are also introduced. Then, the second part of the article reviews the applications for which the $H$ formulations has been used to calculate AC losses, ranging from individual tapes, to complex cables and large magnet windings. Afterwards, a section is dedicated to the discussion of the $H$ formulation in terms of computational efficiency and easiness of implementation. Its pros and cons are listed. Finally, the last section draws the main conclusions.

preprint2019arXiv

Cross-field demagnetization of stacks of tapes: 3D modelling and measurements

Stacks of superconducting tapes can trap much higher magnetic fields than conventional magnets. This makes them very promising for motors and generators. However, ripple magnetic fields in these machines present a cross-field component that demagnetizes the stacks. At present, there is no quantitative agreement between measurements and modeling of cross-field demagnetization, mainly due to the need of a 3D model that takes the end effects and real micron-thick superconducting layer into account. This article presents 3D modeling and measurements of cross-field demagnetization in stacks of up to 5 tapes and initial magnetization modeling of stacks of up to 15 tapes. 3D modeling of the cross-field demagnetization explicitly shows that the critical current density, $J_c$, in the direction perpendicular to the tape surface does not play a role in cross-field demagnetization. When taking the measured anisotropic magnetic field dependence of $J_c$ into account, 3D calculations agree with measurements with less than 4 % deviation, while the error of 2D modeling is much higher. Then, our 3D numerical methods can realistically predict cross-field demagnetization. Due to the force-free configuration of part of the current density, J, in the stack, better agreement with experiments will probably require measuring the Jc anisotropy for the whole solid angle range, including $J$ parallel to the magnetic field.

preprint2019arXiv

Non-twisted stacks of coated conductors for DC magnets: analysis of inductance and AC losses

In the last 10-15 years, the most common strategy in the development of High Temperature Superconducting (HTS) cable for magnets has been to imitate Low Temperature Superconducting (LTS) cable designs. However, requirements for LTS materials are not valid for HTS materials, which are extremely stable. For example, non-twisted multifilamentary Bi-2223 tapes have been successfully used in several magnets. This paper review stability and analyse inductance and AC losses in non-twisted stack of HTS tapes. Numerical calculations show that twisting has negligible effect on inductance variations in a stack of tapes. Regarding AC losses, any magnet built with coated conductors have larger losses than LTS ones, because of the aspect ratio and large width of the tape. If a wide tape is replaced by a non-twisted stack of narrow tapes, losses and residual magnetisation could be reduced. In contrast with multifilamentary wires, twisting a stack of tapes reduces losses only marginally. Therefore, cables composed of non-twisted stack could be designed to have losses comparable to the one of twisted stack concepts. Few examples of large cables for fusion applications are discussed. Designs based on non-twisted stacks can be simpler, more robust and cost effective than twisted ones.

preprint2019arXiv

Screening Current-Induced Field and Field Drift Study in HTS coils using T-A homogenous model

The emergence of second generation (2G) high-temperature superconductor (HTS) tapes has favored the development of HTS magnets for their applications in areas such as NMR, MRI and high field magnets. The screening current-induced field and the field drift are two major problems hindering the use of HTS tapes in the mentioned areas. Both problems are caused by the screening current, then it is necessary to have a modeling strategy capable to estimate such phenomena. Thus far, the H formulation has been the most successful and used approach to model medium-size systems (hundreds of tapes). However, its application to large-scale systems is still impaired by excessive computation times and memory requirements. Homogenization and multi-scaling strategies have been successfully implemented to increase the computational efficiency. In this contribution, we show that using the homogenization technique with the recently developed T-A formulation allows reducing the computation time and the amount of memory up to the point that real-time simulations of slow ramping cycles of large-scale systems are possible. The T-A homogeneous model also allows systematically investigating the screening current using numerical simulations.

preprint2019arXiv

T-A Formulation to Model Electrical Machines with HTS Coated Conductor Coils

Modelling high temperature superconductor (HTS) motors remains challenging mainly due to the high aspect ratio of these conductors but also because of the properties of the magnetic materials. This paper presents a 2D time dependent model to assess the AC losses of superconducting motors based on the new T-A formulation, which by using Finite Element Methods (FEM), allows its implementation in commercial software. The T-A formulation computes the magnetic flux density with different Maxwell's equations depending on the areas of the motor and makes it possible to use the thin strip approximation i.e. the HTS tapes are modelled as infinitely thin lines. The model is then expected to tackle the high aspect ratio of the HTS as well as decreasing both the mesh complexity and the computing time. The first objective of the paper is to validate the method in 2D by evaluating the AC losses of a specific synchronous motor called SUTOR; the computed results are compared with good agreements to those assessed with the MEMEP method, already validated. In a second part, the same losses are computed, taking into account the anisotropy of Jc with the implementation of a data set based on experimentally measured Ic at 65 K and 77 K.