# Low-Frequency Analytical Model of Superconducting Magnet Impedance

Authors: Michele Martino
Journal: Modelling and Simulation in Engineering
 (2022)

---

## Abstract

A superconducting magnet for particle accelerators is often modeled as an ideal inductor, as it indeed exhibits a completely negligible resistance; this is fully satisfactory, as an example, for control purposes, as the time constant formed by the magnet inductance and the resistance of normal conducting cables connecting it to the power converter accurately describe the essentially dominant dynamics of the circuit. Such a model would however fail to correctly represent the noise attenuation mechanism at play in practical superconducting magnets, which also include a vacuum pipe or a beam screen in the inner part of the aperture, an iron yoke on the outer part, and, potentially, a stainless steel or aluminum collar in between. Even at relatively low frequencies, a more accurate model is therefore needed. A sufficiently general one is proposed and illustrated.


---

## Body

## 1. Introduction

This work is aimed at improving the impedance modeling of superconducting magnets for particle accelerators in the low-frequency range (i.e., from DC to a few kHz). Its scope is also limited to individual magnet modelling; in practical magnet circuits, several magnets might be connected in series so the circuit can have significant physical length; its impedance might then need to be modelled as a transmission line which outgoes the scope of this work.Availability of a low-frequency model is relevant as it allows for the proper specification of noise performances, within the above-mentioned range of frequencies, in magnet circuits of particle accelerators, as an example for the HL-LHC (High-Luminosity Large Hadron Collider) project [1–3]. Clearly, noise in higher frequency can also have a significant impact on the beam, but such an analysis far outgoes the scope of this contribution and would need to be carried out decade by decade [3] as different phenomena come to play in different frequency ranges and from different causes (not necessarily related to the noise produced by the power converters). The exact type of noise analysis carried out in this work will be clarified in the following; its focus is, however, limited to the impact of voltage noise (inevitably produced by power converters) on the magnetic field to be experienced by the beam. The proposed model should nevertheless be considered instrumental to further analysis on beam dynamics figures of merit which however outgoes the scope of this work. It is fundamentally a generalization of the content presented in the seminal paper [4] carried out with a more rigorous and coherent notation. Indeed, in [4], the frequency domain, Laplace domain, and time domain notations are often used all at once, which does not allow for a concise and rigorous analytical formulation of the overall magnet impedance as seen by the circuit terminals (i.e., by the power converter terminals). Such a drawback is overcome here, and several generalizations are also presented and discussed together with an illustrative case study.The content is organized as follows: in Section2, the basic assumptions are presented together with the relevant results and their complete derivation. In Section 3 the presented results are translated into a simple equivalent circuit, generalizing the one presented in [4]. Particular emphasis is given to the noise analysis which can be thought of as the main goal of this work. A dominant effect is identified and described in detail in Section 4.At the same time, an exact formulation is presented for an ideal dipole magnet, and a simple case study, from the Large Hadron Collider, is discussed to better illustrate the soundness of the approximations introduced. With regard to important practical features, only nominal operating conditions are considered, as they are the only relevant ones for the quality of circulating beams. The impact of the different structural elements (collars, yoke, etc.) on quench events or on faults (electromechanical, thermal, etc.), although quite critical for practical operation, is therefore not investigated herein. The presence of an iron yoke and the impact of magnetoresistance, together with a more accurate calculation of the inductance (to the best of the author’s knowledge, such a formula has not been presented elsewhere), are briefly addressed in the Appendix to help readability. Summarizing remarks are finally given in Section 5.


## 2. Low-Frequency Field Modeling

### 2.1. Quasistationary Magnetic Modelling

The fundamental assumption herein is that the displacement current densityJD=∂D/∂t is negligible with respect to the intensity of the magnetic field internal stray capacitive effects (in the nF range), such as, for example, the ones between coils and stainless steel or aluminum collars, are herein neglected as they normally kick in at higher frequencies (tens of kHz) (D being the electric displacement field). Furthermore, the physical dimensions of superconducting magnets are such that no propagation phenomena need to be considered. Maximum physical dimensions considered are indeed negligibly small compared to the shortest reasonable wavelength; as an example, one can consider a maximum frequency of 1 MHz (well beyond the range of frequencies of interest in this work) for which the corresponding wavelength would be of about 300 m, so a 10 m-long magnet would still be accurately approximated.With both these conditions, the modeling presented here falls within the so-calledQuasistationary Magnetics.


### 2.2. 2D Modelling

In addition to the above hypotheses, to further simplify the modeling, the lengths of magnets (along the conventional axisz) are considered large enough, so a 2D approximation is deemed satisfactory (edge effects are neglected).


### 2.3. Absence of Nonlinear Materials

All materials considered here are assumed to be linear. The more realistic case of superconducting magnets with a concentric iron yoke is discussed in the Appendix although a full investigation of its nonlinear and frequency-dependent impact on the presented results is beyond the scope of this work. Furthermore, conductive materials experience magnetoresistance; this is also briefly addressed in the Appendix; although the impact of this effect is likely to be smaller than the one caused by the saturation of the iron yoke, its full investigation is beyond the scope of the work presented.


### 2.4. Thin Layer Approximation

Following closely the analysis conducted in [4], and using the same notation, it can be shown that a current Is (with s being the variable of the Laplace domain, the initial conditions are all assumed to be zero) flowing perpendicularly through an annulus having an average radius r=b and a thickness Δb, as shown in Figure 1, produces, for a magnet of order n (n=1 for a dipole, n=2 for a quadrupole, etc.), a magnetic vector potential:
(1)Azr,θ,s=rbnsgnb−rAzb,θ,s,Azb,θ,s=μ0bΔb2nJzb,θ,s.Figure 1
Simplified cross-section (not to scale) of an ideal superconducting magnet with superimposed adopted cylindrical coordinates): coils are located atr=b, an inner conductive shell (beam screen or vacuum pipe for example) is located at r=a, and an outer conductive shell (stainless steel or Al collar as an example) is located at r=c.The quasistatic current density is(2)Jzr,θ,s=Ntcosnθ2Is1bΔbδr−bΔb,where Nt is the number of turns, and therefore, IsNt/2cosnθ represents the azimuthal current density distribution within the annulus (it is assumed that Jzb,θ,s=Ntcosnθ/2Is1/bΔb as equation (2) is to be correctly interpreted as a generalized function).It must be noted that in equation (1), the factor bΔb (proportional to the area of the annulus) is present both at the numerator and at the denominator; therefore, the intensity of the potential vector depends only on the intensity of the current I and its azimuthal distribution. This expedient allows factoring the expression of the potential vector, which turns out to be useful in the following; furthermore, this formulation does not have the drawback of the one proposed in [4] where the current density was expressed in Am−1 instead of the correct Am−2.It is probably worth highlighting the notation introduced in equation (1), and adopted throughout the work, which exploits the signum function sgn as an exponent. This is done to make the expression of the magnetic vector potential more compact, although less conventional than in [4–6], by avoiding the explicit description of two cases at each occurrence.It is fundamental to note that the complete mathematical derivation has nevertheless been carried out assuming that the thicknessesΔa, Δb, and Δc are (negligibly) small compared to the respective radii a, b, and c such that all fields in the annular sections can be considered radially uniform. This clearly represents an idealization which might not always accurately represent practical superconducting magnets for particle accelerators; these aspects will be addressed, when relevant, in the next sections and in the Appendix.


### 2.5. Inner Conductive Layer (Beam Screen)

Atr=a (a<b), there is a thin sheet of conductor with thickness Δa representing a so-called beam screen (which could also represent a simple vacuum pipe where an actual beam screen is absent). A current is induced which can be written as
(3)Jza,θ,s=σiEz=−sσiAza,θ,s+Azia,θ,s,where σi is the conductivity of the layer and Az+Azi represents the total vector potential Aztot.The dependence onθ and s will be omitted in the following. Azir can be expressed by considering the new source Jzia as
(4)Azir=μ0aΔa2nransgna−rJzia=μ0aΔa2nransgna−r−sσiAztota=−sσiμ0aΔa2nransgna−rAza+Azia=−sτiransgna−rμ0bΔb2nabnsgnb−aJzb+ransgna−rAzia,where(5)τi=σiμ0aΔa2nhas the dimensions of time.Azir can therefore be rewritten as follows:
(6)Azir=−sτiabnransgna−rbrnsgnb−rμ0bΔb2nrbnsgnb−rJzb+Azir=−sτikirAzr+Azir,where
(7)kir=1,r≤a,ar2n,a<r<b,ab2n=ki,r≥b.Therefore,(8)Azir=−sτikir1+sτiAzr.The total vector potential in this case isAztotr=Azir+Azr which can be finally calculated as
(9)Aztotr=1+sτi1−kir1+sτiAzr.Therefore, the inner conductive layer introduces a rational transfer function (between thesource potential vector and the total one) which has a pole with time constant τi and a zero with time constant τi1−kir. The position of the pole is determined only by the conductivity of the annulus and its area (actually, by the product of the two) together with the order of the magnet, while the position of the zero also depends on the radial coordinate r. In particular, in the interior of the conductive layer r≤a, kir=1, the zero disappears. This implies that within the range of frequencies considered herein, the field in the interior of the conductive layer keeps being attenuated as the frequency increases.


### 2.6. Outer Conductive Layer (Stainless Steel or Al Collar)

Atr=c (c>b), there is a thin sheet of conductor with thickness Δc representing a stainless steel or Al collar often present in the magnet design. A current is induced, which can be written as
(10)Jzc,θ,s=σoEz=−sσoAzc,θ,s+Azoc,θ,s,where σo is the conductivity of the sheet and Az+Azo represents the total vector potential Aztot.Azor can be expressed by considering the new source Jzoc as
(11)Azor=μ0cΔc2nrcnsgnc−rJzoc=μ0cΔc2nrcnsgnc−r−sσoAztotc=−sσoμ0cΔc2nrcnsgnc−rAzc+Azoc=−sτorcnsgnc−rμ0bΔb2ncbnsgnb−cJzb+rcnsgnc−rAzoc,where
(12)τo=σoμ0cΔc2n.Azor can therefore be rewritten as follows:
(13)Azor=−sτobcnrcnsgnc−rbrnsgnb−rμ0bΔb2nrbnsgnb−rJzb+Azor=−sτokorAzr+Azor,where
(14)kor=bc2n=ko,r≤b,br2n,b<r<c,1,r≥c.Hence,(15)Azor=−sτokor1+sτoAzr,from which the total vector potential can be calculated as(16)Aztotr=1+sτo1−kor1+sτoAzr.The outer conductive layer also introduces a rational transfer function which has a pole with time constantτo and a zero with time constant τo1−kor. The position of the pole, as for the inner layer, is determined only by the conductivity of the annulus and its area (actually by the product of the two) together with the order of the magnet, while the position of the zero also depends on the radial coordinate r. In particular, in the interior of the coils r≤b, kor=ko; hence, the zero is at constant position and its time constant is smaller than the one of the pole (zero is at higher frequency w.r.t. the pole); in this case, no further attenuation of the source field is experienced beyond the frequency of the zero.


### 2.7. Inner and Outer Conductive Layers

In this case, the total potential vector must account for both contributions:(17)Aztot=Azr+Azir+Azor,which will be discussed first individually and then combined in the following subsections.


#### 2.7.1. Inner Layer Contribution

(18)Azir=μ0aΔa2nransgna−rJzia=μ0aΔa2nransgna−r−sσiAztota=−sσiμ0aΔa2nransgna−rAza+Azia+Azoa=−sτiransgna−rμ0bΔb2nabnsgnb−aJzb+Azia+μ0cΔc2nacnsgnc−aJzoc.Rearranging and factoring terms, it yields(19)Azir=−sτikirAzr+Azir+koirAzor,where
(20)koir=abnransgna−rcrnsgnc−r=1,r≤a,ar2n,a<r<c,ac2n,r≥c.


#### 2.7.2. Outer Layer Contribution

(21)Azor=μ0cΔc2nrcnsgnc−rJzoc=μ0cΔc2nrcnsgnc−r−sσoAztotc=−sτorcnsgnc−rAzc+Azic+Azoc=−sτorcnsgnc−rμ0bΔb2ncbnsgnb−cJzb+μ0aΔa2ncansgna−cJzia+Azoc.Rearranging and factoring terms, it yields(22)Azor=−sτokorAzr+kiorAzir+Azor,where
(23)kior=acnrcnsgnc−rarnsgna−r=ac2n,r≤a,rc2n,a<r<c,1,r≥c.


#### 2.7.3. Combining Contributions

In the following, the dependence onr will be omitted unless expressly needed.
(24)Azi=−sτikiAz+Azi+koiAzo,Azo=−sτokoAz+kioAzi+Azo,(25)Azi=−sτi1+sτikiAz+koiAzo,Azo=−sτo1+sτokoAz+kioAzi.Two transfer functions are defined:(26)His=sτi1+sτi,(27)Hos=sτo1+sτo,so that the system can be rewritten as(28)Azi+HiskoiAzo=−HiskiAz,HoskioAzi+Azo=−HoskoAz.Its solution is as follows.(29)Azi=−ki−koiHosko1−HiHokiokoiHisAz,Azo=−ko−kiHiskio1−HiHokiokoiHosAz,which can be expanded as follows:(30)Azi=−sτi1+sτiki−sτo/1+sτokoiko1−sτi/1+sτisτo/1+sτokoikioAz,Azo=−sτo1+sτoko−sτi/1+sτikikio1−sτi/1+sτisτo/1+sτokoikioAz.The final expression, in canonical form, is(31)Azi=−sτiki−sτokoiko−ki1+sτi+τo+s21−koikioAz,Azo=−sτoko−sτikikio−ko1+sτi+τo+s21−koikioAz.It is now possible to write the total potential vectorAztotr=Azr+Azir+Azor as a function of Azr, produced by the source current, as follows:
(32)Aztotr=Azr1+sτi1−ki+τo1−ko+s2τiτoγr1+sτi+τo+s2τiτo1−koikio,where γr=1−ki−ko+kikio+kokoi−koikio.Aski,ko,koi,kio all depend on r, the above expression is fully general and allows for the calculation of the total potential vector anywhere. From this general expression, it is easy to highlight the presence of a rational transfer function between Azr and Aztotr that has two poles and two zeros. Only the position of the two zeros depends on r as
(33)koirkior=ac2n=kiko≤1,therefore, the coefficients of the denominator are fixed, hence the position of the poles. It is also straightforward to verify that such a transfer function is a minimum phase: the real part of both poles and zeros is negative (second-order polynomials with all positive coefficients; it can also be easily shown that the two poles are both real).


## 2.1. Quasistationary Magnetic Modelling

The fundamental assumption herein is that the displacement current densityJD=∂D/∂t is negligible with respect to the intensity of the magnetic field internal stray capacitive effects (in the nF range), such as, for example, the ones between coils and stainless steel or aluminum collars, are herein neglected as they normally kick in at higher frequencies (tens of kHz) (D being the electric displacement field). Furthermore, the physical dimensions of superconducting magnets are such that no propagation phenomena need to be considered. Maximum physical dimensions considered are indeed negligibly small compared to the shortest reasonable wavelength; as an example, one can consider a maximum frequency of 1 MHz (well beyond the range of frequencies of interest in this work) for which the corresponding wavelength would be of about 300 m, so a 10 m-long magnet would still be accurately approximated.With both these conditions, the modeling presented here falls within the so-calledQuasistationary Magnetics.


## 2.2. 2D Modelling

In addition to the above hypotheses, to further simplify the modeling, the lengths of magnets (along the conventional axisz) are considered large enough, so a 2D approximation is deemed satisfactory (edge effects are neglected).


## 2.3. Absence of Nonlinear Materials

All materials considered here are assumed to be linear. The more realistic case of superconducting magnets with a concentric iron yoke is discussed in the Appendix although a full investigation of its nonlinear and frequency-dependent impact on the presented results is beyond the scope of this work. Furthermore, conductive materials experience magnetoresistance; this is also briefly addressed in the Appendix; although the impact of this effect is likely to be smaller than the one caused by the saturation of the iron yoke, its full investigation is beyond the scope of the work presented.


## 2.4. Thin Layer Approximation

Following closely the analysis conducted in [4], and using the same notation, it can be shown that a current Is (with s being the variable of the Laplace domain, the initial conditions are all assumed to be zero) flowing perpendicularly through an annulus having an average radius r=b and a thickness Δb, as shown in Figure 1, produces, for a magnet of order n (n=1 for a dipole, n=2 for a quadrupole, etc.), a magnetic vector potential:
(1)Azr,θ,s=rbnsgnb−rAzb,θ,s,Azb,θ,s=μ0bΔb2nJzb,θ,s.Figure 1
Simplified cross-section (not to scale) of an ideal superconducting magnet with superimposed adopted cylindrical coordinates): coils are located atr=b, an inner conductive shell (beam screen or vacuum pipe for example) is located at r=a, and an outer conductive shell (stainless steel or Al collar as an example) is located at r=c.The quasistatic current density is(2)Jzr,θ,s=Ntcosnθ2Is1bΔbδr−bΔb,where Nt is the number of turns, and therefore, IsNt/2cosnθ represents the azimuthal current density distribution within the annulus (it is assumed that Jzb,θ,s=Ntcosnθ/2Is1/bΔb as equation (2) is to be correctly interpreted as a generalized function).It must be noted that in equation (1), the factor bΔb (proportional to the area of the annulus) is present both at the numerator and at the denominator; therefore, the intensity of the potential vector depends only on the intensity of the current I and its azimuthal distribution. This expedient allows factoring the expression of the potential vector, which turns out to be useful in the following; furthermore, this formulation does not have the drawback of the one proposed in [4] where the current density was expressed in Am−1 instead of the correct Am−2.It is probably worth highlighting the notation introduced in equation (1), and adopted throughout the work, which exploits the signum function sgn as an exponent. This is done to make the expression of the magnetic vector potential more compact, although less conventional than in [4–6], by avoiding the explicit description of two cases at each occurrence.It is fundamental to note that the complete mathematical derivation has nevertheless been carried out assuming that the thicknessesΔa, Δb, and Δc are (negligibly) small compared to the respective radii a, b, and c such that all fields in the annular sections can be considered radially uniform. This clearly represents an idealization which might not always accurately represent practical superconducting magnets for particle accelerators; these aspects will be addressed, when relevant, in the next sections and in the Appendix.


## 2.5. Inner Conductive Layer (Beam Screen)

Atr=a (a<b), there is a thin sheet of conductor with thickness Δa representing a so-called beam screen (which could also represent a simple vacuum pipe where an actual beam screen is absent). A current is induced which can be written as
(3)Jza,θ,s=σiEz=−sσiAza,θ,s+Azia,θ,s,where σi is the conductivity of the layer and Az+Azi represents the total vector potential Aztot.The dependence onθ and s will be omitted in the following. Azir can be expressed by considering the new source Jzia as
(4)Azir=μ0aΔa2nransgna−rJzia=μ0aΔa2nransgna−r−sσiAztota=−sσiμ0aΔa2nransgna−rAza+Azia=−sτiransgna−rμ0bΔb2nabnsgnb−aJzb+ransgna−rAzia,where(5)τi=σiμ0aΔa2nhas the dimensions of time.Azir can therefore be rewritten as follows:
(6)Azir=−sτiabnransgna−rbrnsgnb−rμ0bΔb2nrbnsgnb−rJzb+Azir=−sτikirAzr+Azir,where
(7)kir=1,r≤a,ar2n,a<r<b,ab2n=ki,r≥b.Therefore,(8)Azir=−sτikir1+sτiAzr.The total vector potential in this case isAztotr=Azir+Azr which can be finally calculated as
(9)Aztotr=1+sτi1−kir1+sτiAzr.Therefore, the inner conductive layer introduces a rational transfer function (between thesource potential vector and the total one) which has a pole with time constant τi and a zero with time constant τi1−kir. The position of the pole is determined only by the conductivity of the annulus and its area (actually, by the product of the two) together with the order of the magnet, while the position of the zero also depends on the radial coordinate r. In particular, in the interior of the conductive layer r≤a, kir=1, the zero disappears. This implies that within the range of frequencies considered herein, the field in the interior of the conductive layer keeps being attenuated as the frequency increases.


## 2.6. Outer Conductive Layer (Stainless Steel or Al Collar)

Atr=c (c>b), there is a thin sheet of conductor with thickness Δc representing a stainless steel or Al collar often present in the magnet design. A current is induced, which can be written as
(10)Jzc,θ,s=σoEz=−sσoAzc,θ,s+Azoc,θ,s,where σo is the conductivity of the sheet and Az+Azo represents the total vector potential Aztot.Azor can be expressed by considering the new source Jzoc as
(11)Azor=μ0cΔc2nrcnsgnc−rJzoc=μ0cΔc2nrcnsgnc−r−sσoAztotc=−sσoμ0cΔc2nrcnsgnc−rAzc+Azoc=−sτorcnsgnc−rμ0bΔb2ncbnsgnb−cJzb+rcnsgnc−rAzoc,where
(12)τo=σoμ0cΔc2n.Azor can therefore be rewritten as follows:
(13)Azor=−sτobcnrcnsgnc−rbrnsgnb−rμ0bΔb2nrbnsgnb−rJzb+Azor=−sτokorAzr+Azor,where
(14)kor=bc2n=ko,r≤b,br2n,b<r<c,1,r≥c.Hence,(15)Azor=−sτokor1+sτoAzr,from which the total vector potential can be calculated as(16)Aztotr=1+sτo1−kor1+sτoAzr.The outer conductive layer also introduces a rational transfer function which has a pole with time constantτo and a zero with time constant τo1−kor. The position of the pole, as for the inner layer, is determined only by the conductivity of the annulus and its area (actually by the product of the two) together with the order of the magnet, while the position of the zero also depends on the radial coordinate r. In particular, in the interior of the coils r≤b, kor=ko; hence, the zero is at constant position and its time constant is smaller than the one of the pole (zero is at higher frequency w.r.t. the pole); in this case, no further attenuation of the source field is experienced beyond the frequency of the zero.


## 2.7. Inner and Outer Conductive Layers

In this case, the total potential vector must account for both contributions:(17)Aztot=Azr+Azir+Azor,which will be discussed first individually and then combined in the following subsections.


### 2.7.1. Inner Layer Contribution

(18)Azir=μ0aΔa2nransgna−rJzia=μ0aΔa2nransgna−r−sσiAztota=−sσiμ0aΔa2nransgna−rAza+Azia+Azoa=−sτiransgna−rμ0bΔb2nabnsgnb−aJzb+Azia+μ0cΔc2nacnsgnc−aJzoc.Rearranging and factoring terms, it yields(19)Azir=−sτikirAzr+Azir+koirAzor,where
(20)koir=abnransgna−rcrnsgnc−r=1,r≤a,ar2n,a<r<c,ac2n,r≥c.


### 2.7.2. Outer Layer Contribution

(21)Azor=μ0cΔc2nrcnsgnc−rJzoc=μ0cΔc2nrcnsgnc−r−sσoAztotc=−sτorcnsgnc−rAzc+Azic+Azoc=−sτorcnsgnc−rμ0bΔb2ncbnsgnb−cJzb+μ0aΔa2ncansgna−cJzia+Azoc.Rearranging and factoring terms, it yields(22)Azor=−sτokorAzr+kiorAzir+Azor,where
(23)kior=acnrcnsgnc−rarnsgna−r=ac2n,r≤a,rc2n,a<r<c,1,r≥c.


### 2.7.3. Combining Contributions

In the following, the dependence onr will be omitted unless expressly needed.
(24)Azi=−sτikiAz+Azi+koiAzo,Azo=−sτokoAz+kioAzi+Azo,(25)Azi=−sτi1+sτikiAz+koiAzo,Azo=−sτo1+sτokoAz+kioAzi.Two transfer functions are defined:(26)His=sτi1+sτi,(27)Hos=sτo1+sτo,so that the system can be rewritten as(28)Azi+HiskoiAzo=−HiskiAz,HoskioAzi+Azo=−HoskoAz.Its solution is as follows.(29)Azi=−ki−koiHosko1−HiHokiokoiHisAz,Azo=−ko−kiHiskio1−HiHokiokoiHosAz,which can be expanded as follows:(30)Azi=−sτi1+sτiki−sτo/1+sτokoiko1−sτi/1+sτisτo/1+sτokoikioAz,Azo=−sτo1+sτoko−sτi/1+sτikikio1−sτi/1+sτisτo/1+sτokoikioAz.The final expression, in canonical form, is(31)Azi=−sτiki−sτokoiko−ki1+sτi+τo+s21−koikioAz,Azo=−sτoko−sτikikio−ko1+sτi+τo+s21−koikioAz.It is now possible to write the total potential vectorAztotr=Azr+Azir+Azor as a function of Azr, produced by the source current, as follows:
(32)Aztotr=Azr1+sτi1−ki+τo1−ko+s2τiτoγr1+sτi+τo+s2τiτo1−koikio,where γr=1−ki−ko+kikio+kokoi−koikio.Aski,ko,koi,kio all depend on r, the above expression is fully general and allows for the calculation of the total potential vector anywhere. From this general expression, it is easy to highlight the presence of a rational transfer function between Azr and Aztotr that has two poles and two zeros. Only the position of the two zeros depends on r as
(33)koirkior=ac2n=kiko≤1,therefore, the coefficients of the denominator are fixed, hence the position of the poles. It is also straightforward to verify that such a transfer function is a minimum phase: the real part of both poles and zeros is negative (second-order polynomials with all positive coefficients; it can also be easily shown that the two poles are both real).


## 2.7.1. Inner Layer Contribution

(18)Azir=μ0aΔa2nransgna−rJzia=μ0aΔa2nransgna−r−sσiAztota=−sσiμ0aΔa2nransgna−rAza+Azia+Azoa=−sτiransgna−rμ0bΔb2nabnsgnb−aJzb+Azia+μ0cΔc2nacnsgnc−aJzoc.Rearranging and factoring terms, it yields(19)Azir=−sτikirAzr+Azir+koirAzor,where
(20)koir=abnransgna−rcrnsgnc−r=1,r≤a,ar2n,a<r<c,ac2n,r≥c.


## 2.7.2. Outer Layer Contribution

(21)Azor=μ0cΔc2nrcnsgnc−rJzoc=μ0cΔc2nrcnsgnc−r−sσoAztotc=−sτorcnsgnc−rAzc+Azic+Azoc=−sτorcnsgnc−rμ0bΔb2ncbnsgnb−cJzb+μ0aΔa2ncansgna−cJzia+Azoc.Rearranging and factoring terms, it yields(22)Azor=−sτokorAzr+kiorAzir+Azor,where
(23)kior=acnrcnsgnc−rarnsgna−r=ac2n,r≤a,rc2n,a<r<c,1,r≥c.


## 2.7.3. Combining Contributions

In the following, the dependence onr will be omitted unless expressly needed.
(24)Azi=−sτikiAz+Azi+koiAzo,Azo=−sτokoAz+kioAzi+Azo,(25)Azi=−sτi1+sτikiAz+koiAzo,Azo=−sτo1+sτokoAz+kioAzi.Two transfer functions are defined:(26)His=sτi1+sτi,(27)Hos=sτo1+sτo,so that the system can be rewritten as(28)Azi+HiskoiAzo=−HiskiAz,HoskioAzi+Azo=−HoskoAz.Its solution is as follows.(29)Azi=−ki−koiHosko1−HiHokiokoiHisAz,Azo=−ko−kiHiskio1−HiHokiokoiHosAz,which can be expanded as follows:(30)Azi=−sτi1+sτiki−sτo/1+sτokoiko1−sτi/1+sτisτo/1+sτokoikioAz,Azo=−sτo1+sτoko−sτi/1+sτikikio1−sτi/1+sτisτo/1+sτokoikioAz.The final expression, in canonical form, is(31)Azi=−sτiki−sτokoiko−ki1+sτi+τo+s21−koikioAz,Azo=−sτoko−sτikikio−ko1+sτi+τo+s21−koikioAz.It is now possible to write the total potential vectorAztotr=Azr+Azir+Azor as a function of Azr, produced by the source current, as follows:
(32)Aztotr=Azr1+sτi1−ki+τo1−ko+s2τiτoγr1+sτi+τo+s2τiτo1−koikio,where γr=1−ki−ko+kikio+kokoi−koikio.Aski,ko,koi,kio all depend on r, the above expression is fully general and allows for the calculation of the total potential vector anywhere. From this general expression, it is easy to highlight the presence of a rational transfer function between Azr and Aztotr that has two poles and two zeros. Only the position of the two zeros depends on r as
(33)koirkior=ac2n=kiko≤1,therefore, the coefficients of the denominator are fixed, hence the position of the poles. It is also straightforward to verify that such a transfer function is a minimum phase: the real part of both poles and zeros is negative (second-order polynomials with all positive coefficients; it can also be easily shown that the two poles are both real).


## 3. Circuital Model

### 3.1. Inductance

Equation (32) allows the calculation of the total magnetic potential vector everywhere; in particular, it allows to calculate it in the coils (r=b), and this is enough to calculate the inductance of the magnet (electrical terminals of the magnet are indeed located at r=b).


#### 3.1.1. Potential Vector in the Coils

Note that(34)kib=ki,kob=ko,koib=ki,kiob=ko,and introducing the following
(35)k′i=1−ki,k′o=1−ko,the expression of the potential vector in the coils can be simplified as follows
(36)Aztotb=Azb1+sτiki′+τoko′+s2τiτo1−ki−ko+kiko1+sτi+τo+s2τiτo1−kiko,which can be rewritten to highlight its dynamics by means of the factor Λs as
(37)Aztotb,θ,s=Azb,θ,sΛs.


#### 3.1.2. From Potential Vector to Inductance

In the quasistationary magnetic conditions assumed throughout this work, the stored magnetic energy can be calculated as(38)U=12∭VJ·Adv.For the 2D approximation used herein, the stored energy per unit length is therefore(39)u=dUdz=12∬SJ·Ads,where the surface S is the cross-section of the magnet; however, J is nonzero only for r=b.Since the total potential vector atr=b is simply the source potential vector times Λs (the transfer function in equation (37)), the energy per unit length can be easily calculated.
(40)u=12∬SJzr,θ,sΛsAzr,θ,srdrdθ=12ΛsI2sμ0Nt28n1bΔbbΔb∫02πcos2nθdθ=Λsμ0Nt2I2sπ16n,where δr−b/Δbrdr has been integrated over all possible values of the radius whose result is equal to bΔb.The inductance per unit lengthℓ can be derived from the energy per unit length as follows:
(41)ℓ=dLdz=1I∂u∂I,which gives(42)ℓ=Λsμ0Nt2π8n=ΛsℓDC.Thestatic or DC inductance can be written as
(43)LDC=ℓDCl=μ0πNt28nl,where l is the magnet length over the z axis (this expression coincides with the one reported in [4]).Finally, the soughtdynamic inductance can be expressed as
(44)Ls=LDCΛs.


### 3.2. Equivalent Circuit

It can be verified that(45)Λs=1+sτi1−ki1+sτo1−ko1+sτi+τo+s2τiτo1−kiko,to which corresponds the equivalent circuit in Figure2.Figure 2
Full circuit model for a lossy superconducting magnet: inner and outer conducting shells.The proposed equivalent circuit represents a generalization of the one presented in [4]. It allows several considerations to be made:
(i)
Not all the current supplied by the power convertericircuit is producing magnetic field in the theoretical aperture of the magnet (r≤b); indeed, a fraction io is shunted by the outer conductive layer(ii)
The magnetic field within the magnet’s aperture can be thought as produced by the currentiBm(iii)
Not all of the field produced byiBm is actually seen by the beam in the interior of the inner conductive shell (such as the interior of a beam screen, r≤a)(iv)
There is indeed a shielding effect, and the difference between the field outside the inner layer (a<r<b) and the field inside the inner layer, (r≤a) can be thought as being produced by the difference current ii=iBm−iBb(v)
In other terms, only the currentiBb is producing the magnetic field seen by the beamThe equivalent circuit allows for an easy calculation of many interesting features. In particular, it allows studying the effect of power converter noise (both voltage and current) on the magnetic field seen by the beam inside the beam screen or vacuum pipe.As an example, from equation (32), the TF (transfer function) between the vector potential in the region r≤a with and without the inner and outer layers is the following:
(46)Aztotr≤aAzr≤a=1+sτo1−ko1+sτi+τo+s2τiτo1−kiko.This TF can be deduced from the equivalent circuit, as in Figure2, as
(47)iBbicircuit=iBbiBmiBmicircuit=11+sτi1+sτi1+sτo1−ko1+sτi+τo+s2τiτo1−kiko.The analogy between the equivalent circuit and the dynamics of the magnetic vector potential cannot be further generalized though as, for example, currentio would be flowing at r=c whereas the proposed circuit is limited to the electrical terminals of the magnet, i.e., r=b. Nevertheless, the equivalent circuit represents still a very useful model, and indeed, one can derive many other figures of interest; as an example, the losses during the energy ramp up and ramp down phases can be calculated by the knowledge of ii and io along with their respective resistances (easily derived by inspection of the equivalent circuit).


### 3.3. Noise Analysis

Equation (47) answers fully the question of how much noise passes from the icircuit to the magnetic field experienced by the beam (r≤a). The overall effect is an attenuation that at high frequency (higher than the frequencies of all poles and zeros) gets stronger by 20 per frequency decade. It is nevertheless clear from the equivalent circuit that the overall magnet impedance is smaller than that of a superconducting magnet that has no inner or outer layer because of the parallel branches (circuit admittance is larger). As such, for a given voltage noise contribution, the overall circuit current noise is larger for a magnet with inner/outer conductive layers with respect to an ideal one. Therefore, the relevant question is how a real magnet compares to an ideal one in terms of noise transfer from the power converter voltage to the magnetic induction field experienced by the beam.The admittance of the circuit shown in Figure3 is expressed by the following equation:
(48)Acircuit=icircuitvcircuit=1Rc1+sτi+τo+s2τiτo1−kiko1+sτΣ+s2τcτi′+τo′+τiτoγio+s3τcτi′τo′,where
(49)τc=LDCRc,τΣ=τc+τi+τo,τi′=1−kiτi,τo′=1−koτo,γio=1−kiko.Figure 3
Full circuit model as seen from a power converter of a lossy superconducting magnet: inner and outer conducting shells.Rc represents the resistance of the normal conducting cables connecting the power converter to the superconducting magnet.The TF of interest is the one fromvcircuit to iBb which is expressed in the following equation.
(50)iBbvcircuit=iBbicircuiticircuitvcircuit=1Rc1+sτo1−ko1+sτΣ+s2τcτi′+τo′+τiτoγio+s3τcτi′τo′.Equation (50) shows that there are one zero and three poles, indicating a stronger attenuation of noise (in high frequency) w.r.t. the no layers case as expressed in the following equation
(51)iBbvcircuitnolayers=1Rc11+sτc.However, in such a form, i.e., equation (50), there is not much insight about the extra filtering; such expression can be considerably simplified noticing that τc>>τi,τo,τi′,τo′; hence, the following approximations hold:
(52)τΣ=τc+τi+τo≈τc+τi′+τΣ′=τi′,τcτi′+τo′+τiτoγio≈τcτi′+τo′+τi′τo′.Combining equation (50) and the approximations in equation (52) gives
(53)iBbvcircuit≈1Rc1+sτo1−ko1+sτΣ′+s2τcτi′+τo′+τi′τo′+s3τcτi′τo′=1Rc1+sτo′1+sτc1+sτi′1+sτo′=1Rc11+sτc11+sτi′.From equation (53), it is straightforward to deduce that
(54)iBbvcircuit≈iBbvcircuitnolayers11+sτi′=iBbvcircuitnolayers11+sτi1−ki.Equation (54) shows that
(i)
there is a first-order additional filtering effect(ii)
only the internal conductive layer (such as a beam screen) contributes to this additional filtering from voltage noise (to the magnetic field experienced by the beam) (in this respect, the usefulness of the proposed modelling holds even for practical magnets when the assumption of athin layer for the stainless steel or Al collars does not)(iii)
the filtering occurs at higher frequencies compared to the filtering effect of the currents (or equivalently of the fields themselves)(iv)
the ratioiBb/iBm indeed depends on the time constant τi, whereas the ratio iBb/iBbnolayers depends on τi′ and τi is larger than τi′=τi1−ki as 0<ki<1


## 3.1. Inductance

Equation (32) allows the calculation of the total magnetic potential vector everywhere; in particular, it allows to calculate it in the coils (r=b), and this is enough to calculate the inductance of the magnet (electrical terminals of the magnet are indeed located at r=b).


### 3.1.1. Potential Vector in the Coils

Note that(34)kib=ki,kob=ko,koib=ki,kiob=ko,and introducing the following
(35)k′i=1−ki,k′o=1−ko,the expression of the potential vector in the coils can be simplified as follows
(36)Aztotb=Azb1+sτiki′+τoko′+s2τiτo1−ki−ko+kiko1+sτi+τo+s2τiτo1−kiko,which can be rewritten to highlight its dynamics by means of the factor Λs as
(37)Aztotb,θ,s=Azb,θ,sΛs.


### 3.1.2. From Potential Vector to Inductance

In the quasistationary magnetic conditions assumed throughout this work, the stored magnetic energy can be calculated as(38)U=12∭VJ·Adv.For the 2D approximation used herein, the stored energy per unit length is therefore(39)u=dUdz=12∬SJ·Ads,where the surface S is the cross-section of the magnet; however, J is nonzero only for r=b.Since the total potential vector atr=b is simply the source potential vector times Λs (the transfer function in equation (37)), the energy per unit length can be easily calculated.
(40)u=12∬SJzr,θ,sΛsAzr,θ,srdrdθ=12ΛsI2sμ0Nt28n1bΔbbΔb∫02πcos2nθdθ=Λsμ0Nt2I2sπ16n,where δr−b/Δbrdr has been integrated over all possible values of the radius whose result is equal to bΔb.The inductance per unit lengthℓ can be derived from the energy per unit length as follows:
(41)ℓ=dLdz=1I∂u∂I,which gives(42)ℓ=Λsμ0Nt2π8n=ΛsℓDC.Thestatic or DC inductance can be written as
(43)LDC=ℓDCl=μ0πNt28nl,where l is the magnet length over the z axis (this expression coincides with the one reported in [4]).Finally, the soughtdynamic inductance can be expressed as
(44)Ls=LDCΛs.


## 3.1.1. Potential Vector in the Coils

Note that(34)kib=ki,kob=ko,koib=ki,kiob=ko,and introducing the following
(35)k′i=1−ki,k′o=1−ko,the expression of the potential vector in the coils can be simplified as follows
(36)Aztotb=Azb1+sτiki′+τoko′+s2τiτo1−ki−ko+kiko1+sτi+τo+s2τiτo1−kiko,which can be rewritten to highlight its dynamics by means of the factor Λs as
(37)Aztotb,θ,s=Azb,θ,sΛs.


## 3.1.2. From Potential Vector to Inductance

In the quasistationary magnetic conditions assumed throughout this work, the stored magnetic energy can be calculated as(38)U=12∭VJ·Adv.For the 2D approximation used herein, the stored energy per unit length is therefore(39)u=dUdz=12∬SJ·Ads,where the surface S is the cross-section of the magnet; however, J is nonzero only for r=b.Since the total potential vector atr=b is simply the source potential vector times Λs (the transfer function in equation (37)), the energy per unit length can be easily calculated.
(40)u=12∬SJzr,θ,sΛsAzr,θ,srdrdθ=12ΛsI2sμ0Nt28n1bΔbbΔb∫02πcos2nθdθ=Λsμ0Nt2I2sπ16n,where δr−b/Δbrdr has been integrated over all possible values of the radius whose result is equal to bΔb.The inductance per unit lengthℓ can be derived from the energy per unit length as follows:
(41)ℓ=dLdz=1I∂u∂I,which gives(42)ℓ=Λsμ0Nt2π8n=ΛsℓDC.Thestatic or DC inductance can be written as
(43)LDC=ℓDCl=μ0πNt28nl,where l is the magnet length over the z axis (this expression coincides with the one reported in [4]).Finally, the soughtdynamic inductance can be expressed as
(44)Ls=LDCΛs.


## 3.2. Equivalent Circuit

It can be verified that(45)Λs=1+sτi1−ki1+sτo1−ko1+sτi+τo+s2τiτo1−kiko,to which corresponds the equivalent circuit in Figure2.Figure 2
Full circuit model for a lossy superconducting magnet: inner and outer conducting shells.The proposed equivalent circuit represents a generalization of the one presented in [4]. It allows several considerations to be made:
(i)
Not all the current supplied by the power convertericircuit is producing magnetic field in the theoretical aperture of the magnet (r≤b); indeed, a fraction io is shunted by the outer conductive layer(ii)
The magnetic field within the magnet’s aperture can be thought as produced by the currentiBm(iii)
Not all of the field produced byiBm is actually seen by the beam in the interior of the inner conductive shell (such as the interior of a beam screen, r≤a)(iv)
There is indeed a shielding effect, and the difference between the field outside the inner layer (a<r<b) and the field inside the inner layer, (r≤a) can be thought as being produced by the difference current ii=iBm−iBb(v)
In other terms, only the currentiBb is producing the magnetic field seen by the beamThe equivalent circuit allows for an easy calculation of many interesting features. In particular, it allows studying the effect of power converter noise (both voltage and current) on the magnetic field seen by the beam inside the beam screen or vacuum pipe.As an example, from equation (32), the TF (transfer function) between the vector potential in the region r≤a with and without the inner and outer layers is the following:
(46)Aztotr≤aAzr≤a=1+sτo1−ko1+sτi+τo+s2τiτo1−kiko.This TF can be deduced from the equivalent circuit, as in Figure2, as
(47)iBbicircuit=iBbiBmiBmicircuit=11+sτi1+sτi1+sτo1−ko1+sτi+τo+s2τiτo1−kiko.The analogy between the equivalent circuit and the dynamics of the magnetic vector potential cannot be further generalized though as, for example, currentio would be flowing at r=c whereas the proposed circuit is limited to the electrical terminals of the magnet, i.e., r=b. Nevertheless, the equivalent circuit represents still a very useful model, and indeed, one can derive many other figures of interest; as an example, the losses during the energy ramp up and ramp down phases can be calculated by the knowledge of ii and io along with their respective resistances (easily derived by inspection of the equivalent circuit).


## 3.3. Noise Analysis

Equation (47) answers fully the question of how much noise passes from the icircuit to the magnetic field experienced by the beam (r≤a). The overall effect is an attenuation that at high frequency (higher than the frequencies of all poles and zeros) gets stronger by 20 per frequency decade. It is nevertheless clear from the equivalent circuit that the overall magnet impedance is smaller than that of a superconducting magnet that has no inner or outer layer because of the parallel branches (circuit admittance is larger). As such, for a given voltage noise contribution, the overall circuit current noise is larger for a magnet with inner/outer conductive layers with respect to an ideal one. Therefore, the relevant question is how a real magnet compares to an ideal one in terms of noise transfer from the power converter voltage to the magnetic induction field experienced by the beam.The admittance of the circuit shown in Figure3 is expressed by the following equation:
(48)Acircuit=icircuitvcircuit=1Rc1+sτi+τo+s2τiτo1−kiko1+sτΣ+s2τcτi′+τo′+τiτoγio+s3τcτi′τo′,where
(49)τc=LDCRc,τΣ=τc+τi+τo,τi′=1−kiτi,τo′=1−koτo,γio=1−kiko.Figure 3
Full circuit model as seen from a power converter of a lossy superconducting magnet: inner and outer conducting shells.Rc represents the resistance of the normal conducting cables connecting the power converter to the superconducting magnet.The TF of interest is the one fromvcircuit to iBb which is expressed in the following equation.
(50)iBbvcircuit=iBbicircuiticircuitvcircuit=1Rc1+sτo1−ko1+sτΣ+s2τcτi′+τo′+τiτoγio+s3τcτi′τo′.Equation (50) shows that there are one zero and three poles, indicating a stronger attenuation of noise (in high frequency) w.r.t. the no layers case as expressed in the following equation
(51)iBbvcircuitnolayers=1Rc11+sτc.However, in such a form, i.e., equation (50), there is not much insight about the extra filtering; such expression can be considerably simplified noticing that τc>>τi,τo,τi′,τo′; hence, the following approximations hold:
(52)τΣ=τc+τi+τo≈τc+τi′+τΣ′=τi′,τcτi′+τo′+τiτoγio≈τcτi′+τo′+τi′τo′.Combining equation (50) and the approximations in equation (52) gives
(53)iBbvcircuit≈1Rc1+sτo1−ko1+sτΣ′+s2τcτi′+τo′+τi′τo′+s3τcτi′τo′=1Rc1+sτo′1+sτc1+sτi′1+sτo′=1Rc11+sτc11+sτi′.From equation (53), it is straightforward to deduce that
(54)iBbvcircuit≈iBbvcircuitnolayers11+sτi′=iBbvcircuitnolayers11+sτi1−ki.Equation (54) shows that
(i)
there is a first-order additional filtering effect(ii)
only the internal conductive layer (such as a beam screen) contributes to this additional filtering from voltage noise (to the magnetic field experienced by the beam) (in this respect, the usefulness of the proposed modelling holds even for practical magnets when the assumption of athin layer for the stainless steel or Al collars does not)(iii)
the filtering occurs at higher frequencies compared to the filtering effect of the currents (or equivalently of the fields themselves)(iv)
the ratioiBb/iBm indeed depends on the time constant τi, whereas the ratio iBb/iBbnolayers depends on τi′ and τi is larger than τi′=τi1−ki as 0<ki<1


## 4. General Model for Inner Layer

In the previous section, it was concluded that the dominant effect is the one due to the inner layer; in the Appendix, it will be shown that, considering only the inner layer, the proposed equivalent circuit is well defined even in the presence of an outer iron yoke. A generalization is now presented by means of the symbolic equivalent circuit shown in Figure4. For such a circuit, where
(55)Hs=AztotsAzs,Figure 4
Circuit model seen from a power converter of a lossy superconducting magnet: most general inner conducting shell.the TF of interest (i.e., between the power converter voltage and the currentiBb) can be written as
(56)iBbvcircuit=iBbvinner·vinnervmagnet·vmagnetvcircuit=1skiHLDC·skiHLDCHssLDC1−kiH+kiHHs·sLDC1−kiH+kiHHsRc+sLDC1−kiH+kiHHs=1RcHs1+sτc1−kiH+kiHHs,where an equivalent geometrical factorkiH is replacing the one introduced in equation (7).Equation (56) simplifies into equation (53) in the thin layer approximation whereby, from equation (9),
(57)Hs=AztotsAzs=11+sτi.All conclusions drawn from equation (54) concerning the noise attenuation are hence valid irrespective to the validity of the thin layer approximation; in particular, the relevant fact that the additional noise attenuation (from voltage-to-current or voltage-to-field) becomes dominant at higher frequencies w.r.t. the current-to-current (or field-to-field) attenuation as kiH≤1.


### 4.1. Full Analytical Formulation for a Dipole

For a dipole magnet, i.e.,n=1, an exact expression is available for the case of the inner conductive layer in the region r<a.In [7], a general formula, equation (58), is reported where the thickness of the layer is not bounded to be much smaller than the radius  a ; such a formula is also a generalization of the one presented in [8] (both of them only apply to dipole magnets, i.e., n=1),
(58)AztotrAzr=2μra−Δ/2/a+Δ/2μrK1γa−−γa−K1′γa−μrI1γa++γa+I1′γa+−μrI1γa−−γa−I1′γa−μrK1γa++γa+K1′γa+,The quantities involved are as follows:(59)γ≈j2πfμ0μrσi=1+jδ,γa−=γ·a−Δ2,γa+=γ·a+Δ2. where δ is the penetration depth of the conductive layer having an electrical conductivity σi. The functions involved are the first-order modified Bessel functions of first I1 and second kind K1 and their first derivatives (w.r.t. their argument). Assuming that μr=1 (which is an excellent approximation) and noting that
(60)zK1′z=−zK0z−K1z,zI1′z=zI0z−I1z, Equation (58) can be considerably simplified to become equation (61),
(61)AztotrAzr=2a−Δ/2/a+Δ/22K1γa−+γa−K0γa−γa+I0γa++2I1γa−−γa−I0γa−γa+K0γa+.


#### 4.1.1. Case of LHC Dipole

A numerical example is illustrated considering the case of the LHC dipoles. The cross-section of the LHC dipole beam screen is not perfectly circular; an equivalent first-order approximated TF is presented in [9]; however, such a cross-section can be assumed to be circular for comparison between the exact analytical model and the first-order one considered within the scope of this work. The whole beam screen is further approximated only with the thin layer of copper, for which, however, the magnetoresistance effect (briefly discussed in the Appendix) is taken into account. The parameters used for the beam screen are the following: a−=a−Δa/2=18.45mm, Δa=75μm, and σi=2.09×109Sm−1 (at 20 K and 8 T). The comparison between the exact expression in equation (61) and the first-order approximation, equation (57), is illustrated in Figure 5 in the frequency domain and Figure 6 in the time domain. It can be observed that the impulse response (impulse response hδt=F−1Hs=j2πf) of the exact TF is zero at t=0 (Figure 6, top) which is a consequence of the fact that the TF is of an order larger than 1; indeed, its phase response (Figure 5, bottom) goes beyond the asymptote of -90° of the 1st order, thin layer, approximation. The overall agreement is excellent up to about 10 kHz whereas, looking only at the amplitude (which is the relevant one for noise considerations), a very good agreement is maintained even beyond about 20 kHz when the penetration depth δ becomes equal to the shell thickness; this is an important result which confirms the validity of the thin layer approximation presented and its relevance for existing and next-to-come particle accelerators.Figure 5
Bode plot comparison: equation (61) in blue vs. equation (57) in red. Top: amplitude; bottom: phase.Figure 6
Impulse (top) and step (bottom) response comparison: equation (61) in blue vs. equation (57) in red.


## 4.1. Full Analytical Formulation for a Dipole

For a dipole magnet, i.e.,n=1, an exact expression is available for the case of the inner conductive layer in the region r<a.In [7], a general formula, equation (58), is reported where the thickness of the layer is not bounded to be much smaller than the radius  a ; such a formula is also a generalization of the one presented in [8] (both of them only apply to dipole magnets, i.e., n=1),
(58)AztotrAzr=2μra−Δ/2/a+Δ/2μrK1γa−−γa−K1′γa−μrI1γa++γa+I1′γa+−μrI1γa−−γa−I1′γa−μrK1γa++γa+K1′γa+,The quantities involved are as follows:(59)γ≈j2πfμ0μrσi=1+jδ,γa−=γ·a−Δ2,γa+=γ·a+Δ2. where δ is the penetration depth of the conductive layer having an electrical conductivity σi. The functions involved are the first-order modified Bessel functions of first I1 and second kind K1 and their first derivatives (w.r.t. their argument). Assuming that μr=1 (which is an excellent approximation) and noting that
(60)zK1′z=−zK0z−K1z,zI1′z=zI0z−I1z, Equation (58) can be considerably simplified to become equation (61),
(61)AztotrAzr=2a−Δ/2/a+Δ/22K1γa−+γa−K0γa−γa+I0γa++2I1γa−−γa−I0γa−γa+K0γa+.


### 4.1.1. Case of LHC Dipole

A numerical example is illustrated considering the case of the LHC dipoles. The cross-section of the LHC dipole beam screen is not perfectly circular; an equivalent first-order approximated TF is presented in [9]; however, such a cross-section can be assumed to be circular for comparison between the exact analytical model and the first-order one considered within the scope of this work. The whole beam screen is further approximated only with the thin layer of copper, for which, however, the magnetoresistance effect (briefly discussed in the Appendix) is taken into account. The parameters used for the beam screen are the following: a−=a−Δa/2=18.45mm, Δa=75μm, and σi=2.09×109Sm−1 (at 20 K and 8 T). The comparison between the exact expression in equation (61) and the first-order approximation, equation (57), is illustrated in Figure 5 in the frequency domain and Figure 6 in the time domain. It can be observed that the impulse response (impulse response hδt=F−1Hs=j2πf) of the exact TF is zero at t=0 (Figure 6, top) which is a consequence of the fact that the TF is of an order larger than 1; indeed, its phase response (Figure 5, bottom) goes beyond the asymptote of -90° of the 1st order, thin layer, approximation. The overall agreement is excellent up to about 10 kHz whereas, looking only at the amplitude (which is the relevant one for noise considerations), a very good agreement is maintained even beyond about 20 kHz when the penetration depth δ becomes equal to the shell thickness; this is an important result which confirms the validity of the thin layer approximation presented and its relevance for existing and next-to-come particle accelerators.Figure 5
Bode plot comparison: equation (61) in blue vs. equation (57) in red. Top: amplitude; bottom: phase.Figure 6
Impulse (top) and step (bottom) response comparison: equation (61) in blue vs. equation (57) in red.


## 4.1.1. Case of LHC Dipole

A numerical example is illustrated considering the case of the LHC dipoles. The cross-section of the LHC dipole beam screen is not perfectly circular; an equivalent first-order approximated TF is presented in [9]; however, such a cross-section can be assumed to be circular for comparison between the exact analytical model and the first-order one considered within the scope of this work. The whole beam screen is further approximated only with the thin layer of copper, for which, however, the magnetoresistance effect (briefly discussed in the Appendix) is taken into account. The parameters used for the beam screen are the following: a−=a−Δa/2=18.45mm, Δa=75μm, and σi=2.09×109Sm−1 (at 20 K and 8 T). The comparison between the exact expression in equation (61) and the first-order approximation, equation (57), is illustrated in Figure 5 in the frequency domain and Figure 6 in the time domain. It can be observed that the impulse response (impulse response hδt=F−1Hs=j2πf) of the exact TF is zero at t=0 (Figure 6, top) which is a consequence of the fact that the TF is of an order larger than 1; indeed, its phase response (Figure 5, bottom) goes beyond the asymptote of -90° of the 1st order, thin layer, approximation. The overall agreement is excellent up to about 10 kHz whereas, looking only at the amplitude (which is the relevant one for noise considerations), a very good agreement is maintained even beyond about 20 kHz when the penetration depth δ becomes equal to the shell thickness; this is an important result which confirms the validity of the thin layer approximation presented and its relevance for existing and next-to-come particle accelerators.Figure 5
Bode plot comparison: equation (61) in blue vs. equation (57) in red. Top: amplitude; bottom: phase.Figure 6
Impulse (top) and step (bottom) response comparison: equation (61) in blue vs. equation (57) in red.


## 5. Conclusion

A quasistationary magnetic model of an ideal superconducting magnet has been presented. Although idealized, the model includes all the important constituents of practical particle accelerators magnets: the beam screen, the iron yoke (addressed in the Appendix), and the collar. The model has been translated into an equivalent circuit generalizing the one presented in [4]. Such an equivalent circuit has been exploited to carry out a noise analysis focusing on the final impact of power converter voltage noise on the magnetic field to be experienced by the beam. A quantitative, although approximated, estimation of additional noise attenuation (w.r.t. a pure inductance model of the magnet) has also been presented and represents the main result of this work. Non-thin inner and outer layers have been briefly addressed in this context, whereas the coils’ layer is discussed with more details in the Appendix where the validity of the presented model is shown to hold.



---

*Source: 2105847-2022-11-07.xml*