Calculating the Inrush Current of Superconducting Transformers

Under certain circumstances, after connecting a superconducting transformer to the power network, a high value current may flow through its windings. This current can exceed the critical value of the superconductor many times and cause the windings to lose their superconductive state. Loss of superconductive state of the windings may result in thermal interruption of their continuity as a result of conduction of a current of very high density. The mathematical relationships used to calculate the inrush current of conventional transformers do not work well for the calculation of superconducting transformers. This is due to the properties of superconducting materials used in the windings, first of all to the stepwise changes of the windings’ resistance when exiting the superconducting state and when returning to this state. The article presents the mathematical dependencies allowing to calculate the pulse waveforms of the inrush current of these transformers are derived. Basic electrical circuit sizes are used in the calculations, making the calculations quick and easy. Using the formulas, calculations of the inrush current of 8.5 kVA and 13.5 kVA superconducting transformers. The results were verified with the results of the inrush current measurements, achieving good compliance.


Introduction
A high value current may flow through the transformer windings not only at short circuit or overload, but under certain circumstances even after the transformer has been connected to the mains. The current at the switching on of an unloaded transformer can reach a value of up to 20 times higher than the rated current [1,2]. The time of the inrush current fade can, in the case of large power transformers, reach several thousand periods of supply voltage.
Conventional transformers, i.e., with copper or aluminum windings, are elements of the power system with a relatively low failure rate. The maintenance of superconducting transformers (HTS transformer) in a faultless state is hampered by the tendency of their windings to easily emerge from the superconductive state. The inrush current of a high value and long decay time may not only make it difficult to connect the superconducting transformer to the mains, but also leads to thermal damage of its windings. Therefore, it is important to counteract the effects of this current [1,3,4]. This requires learning about the mechanism for the development of the inrush current in superconducting transformers, which are different from those in conventional transformers, as well as for the capturing of this phenomenon by mathematical relationships. This will allow the methods of limiting the inrush current of conventional transformers, which are still being developed [5][6][7][8], to be transferred to superconducting transformers.
The inrush current of conventional transformers is a relatively well-recognised phenomenon, both through theoretical analysis and experimental research on real transformers [9][10][11][12]. Specht and Holcomb [13][14][15] in 1951 and 1961 gave relationships for the inrush current of conventional transformers that are still used in calculations today [16]. Formulas to approximate the shape of the inrush current are proposed in [17,18]. Transformer models during inrush current are presented in [19][20][21][22][23][24][25]. Some of these models are based on electromagnetic-circuit analysis [21,25]. These models need to have transformer detail data which may be difficult to obtain. The B-H curve of the transformer core plays an important role in the turn-on current analysis. Analytical equations using the B-H curve to model the turn-on current are given in papers [17,18,20,22]. The magnitude of the first cycle inrush current is important for transformer design and protection. Therefore, analytical formulas to estimate the peak of the inrush current during the first cycle as a function of switching angle are presented in [26][27][28][29]. Due to the importance of the inrush current envelope for transformer and system protection, analytical formulas to estimate the inrush current envelope are presented in [13,14,30]. These formulas require some detailed data like detailed core dimensions, in addition to the electrical parameters and magnetic core characteristics.
The literature on the inrush current of superconducting transformers is relatively scarce and concentrates mainly on its waveforms obtained experimentally for existing transformers [31][32][33]. Meanwhile, calculation methods used in calculating the inrush current of conventional transformers do not work well in calculating this current for superconducting transformers, mainly because of the different parameters and characteristics of the superconducting wires with respect to conventional winding wires. This paper proposes a solution to this problem.
The mathematical descriptions proposed in the presented literature does not take into account dynamic changes of the windings' resistance, which leads to incorrect calculations for superconducting transformers. This article proposes a solution to this problem and derives mathematical relationships describing the course of the inrush current of superconducting transformers. Due to the limited amount of data concerning the inrush current phenomenon of superconducting transformers, the given dependencies were verified on the designed and built by the author of the paper transformers of this type.

Superconducting Transformer
Superconducting transformers differ from conventional ones mainly in the material from which their superconducting windings are made [34][35][36]. The use of superconducting winding materials determines the remaining differences in the construction of both types of transformers [37][38][39].
The diagram of an HTS transformer in the idle state ( Figure 1) does not differ significantly from that of a conventional transformer, with the difference that resistance R representing power losses in the high-temperature superconducting primary winding (HTS winding) of the HTS transformer is nonlinear as a function of changes in current I, magnetic field strength H and temperature T. The nonlinearity of this resistance has a significant impact on the inrush current of the transformer. In the superconductive state of the winding, resistance R is equal to zero. This resistance increases rapidly, in approximately 10 −12 s, when the critical parameters of the superconductor are exceeded, when the winding goes to the resistive state. Critical current Ic, critical temperature Tc and critical magnetic field strength Hc determine the critical area of the superconducting winding ( Figure 2) [40,41]. The real critical parameters Tcr, Hcr, Icr of the superconducting winding have lower values and determine the location of the superconductor's working point (Pcr) on the critical surface. Maintaining the superconductivity of the winding is only possible when the temperature,  Critical current I c , critical temperature T c and critical magnetic field strength H c determine the critical area of the superconducting winding ( Figure 2) [40,41]. The real critical parameters T cr , H cr , I cr of the superconducting winding have lower values and determine the location of the superconductor's working point (P cr ) on the critical surface. Maintaining the superconductivity of the winding is only possible when the temperature, magnetic field strength and current strength are simultaneously lower than the real critical values, i.e., when the superconductor's working point (P w ) is located under a critical surface and is determined by T w , H w , I w . Critical current Ic, critical temperature Tc and critical magnetic field strength Hc determine the critical area of the superconducting winding ( Figure 2) [40,41]. The real critical parameters Tcr, Hcr, Icr of the superconducting winding have lower values and determine the location of the superconductor's working point (Pcr) on the critical surface. Maintaining the superconductivity of the winding is only possible when the temperature, magnetic field strength and current strength are simultaneously lower than the real critical values, i.e., when the superconductor's working point (Pw) is located under a critical surface and is determined by Tw, Hw, Iw. Due to the high critical values of the magnetic field strength of modern superconducting materials, its influence on the course of the inrush current was omitted from the considerations. The impact of temperature changes is considered indirectly by assuming lower critical current values of the superconductor with its temperature increase. This approach allowed to simplify the transition characteristics of the HTS transformer's winding to the function of two variables R = f(I), as shown in Figure 3. Due to the high critical values of the magnetic field strength of modern superconducting materials, its influence on the course of the inrush current was omitted from the considerations. The impact of temperature changes is considered indirectly by assuming lower critical current values of the superconductor with its temperature increase. This approach allowed to simplify the transition characteristics of the HTS transformer's winding to the function of two variables R = f (I), as shown in Figure 3.

Stream in the Transformer Core after Switching on
Knowing the time course of the flux, after switching on the transformer, it is possible to determine the time course of the inrush current from the magnetisation characteristic Φ = f(I) [42,43].
Initially, the connection of an unloaded single-phase transformer ( Figure 1) to a sinusoidal source ( Figure 4) will be considered. The constant resistance R is assumed for the entire duration of the inrush current. The impedance of the power grid and the impedance of the power source are ignored.

Stream in the Transformer Core after Switching on
Knowing the time course of the flux, after switching on the transformer, it is possible to determine the time course of the inrush current from the magnetisation characteristic Φ = f (I) [42,43].
Initially, the connection of an unloaded single-phase transformer ( Figure 1) to a sinusoidal source ( Figure 4) will be considered. The constant resistance R is assumed for the entire duration of the inrush current. The impedance of the power grid and the impedance of the power source are ignored.

Stream in the Transformer Core after Switching on
Knowing the time course of the flux, after switching on the transformer, it is possible to determine the time course of the inrush current from the magnetisation characteristic Φ = f(I) [42,43].
Initially, the connection of an unloaded single-phase transformer ( Figure 1) to a sinusoidal source ( Figure 4) will be considered. The constant resistance R is assumed for the entire duration of the inrush current. The impedance of the power grid and the impedance of the power source are ignored. The course of the inrush current is influenced by the shape of the magnetisation curve in the saturation area of the transformer core [44,45]. Figure 5 shows a simplified shape of the magnetisation curve adopted for the calculation [15,46,47].
The circuitry in Figure 1 is described by the equations [18]:  The course of the inrush current is influenced by the shape of the magnetisation curve in the saturation area of the transformer core [44,45]. Figure 5 shows a simplified shape of the magnetisation curve adopted for the calculation [15,46,47].  Assuming that the transformer is switched on at the angle −π ( Figure 4) the basic equation of the transformer operation can be written as: After the integration of Equation (3), the expression for flux is obtained: where A is the integration constant determined from the initial conditions. After the transformer has been connected to the power network, the value of the excitation current at the first moment must be equal to zero. However, there is a residual magnetism flux Φr (φ = Φr dla ωt = −π) and therefore the constant A is: Ultimately, the flux equation has a form: The circuitry in Figure 1 is described by the equations [18]: Assuming that the transformer is switched on at the angle −π ( Figure 4) the basic equation of the transformer operation can be written as: Energies 2021, 14, 6714

of 19
After the integration of Equation (3), the expression for flux is obtained: where A is the integration constant determined from the initial conditions. After the transformer has been connected to the power network, the value of the excitation current at the first moment must be equal to zero. However, there is a residual magnetism flux Φ r (ϕ = Φ r dla ωt = −π) and therefore the constant A is: Ultimately, the flux equation has a form: where the angle of the phase shift δ is: Equation (6) contains two components-fixed: and disturbing: From relation (6) it follows that the magnetic flux ϕ reaches its highest peak moment after switching on the transformer. This value depends on the phase of the supply voltage e and the value of the residual magnetism flux Φ r at the moment of switching on the transformer and the value of the ratio R/X of the electrical circuit [48].
Let us assume that we are omitting the resistance of the primary winding and the active component to cover iron losses (R = 0 i R Fe = ∞) and the flux in the transformer core is delayed to the voltage by an angle π/2 ( Figure 4). The inrush current will not appear when the transformer is switched on when the voltage e reaches its maximum value and there is no residual magnetism flux in the transformer core (Φ r = 0).
The least favourable conditions for switching on the transformer occur when the supply voltage e passes through zero. Because at the moment of switching on the transformer the flux in the core must be equal to zero (if the residual magnetism flux is omitted, i.e., Φ r = 0), a one-way flux component must occur. This component, when added to the course of a fixed flux, causes the peak value of the flux to double: As can be seen from relation (6), the residual magnetism flux Φ r has a significant influence on the course of the flux ϕ. Expression (6) reaches its highest value: Such a high flux value leads to saturation of the transformer core, at which the peak value of the inrush current reaches significant values.
From Equation (6) it appears that the reaction X is a measure of circuit inertia. Resistance R is responsible for damping the unidirectional component of the disturbing flux, and thus reducing the peak value of the inrush current and shortening its duration [16]. The value of the flux is reduced by the value expressed in the equation: If the resistance is omitted in the circuit from Figure 1 (R = 0), then the attenuation of the one-way flux component does not occur and the flux has an undamped course described by the relation: Therefore, the inrush current is also undamped.

Inrush Current of Superconducting Transformer
The inrush current, similar to the magnetic flux described by Equation (6), has two components: fixed and transient. The fixed component is the periodically alternating current drawn by the transformer to excite the flux in the core and cover the power losses in the transformer in its steady state. This current is between 1 and 10% of the rated current, so it will be omitted from consideration. The disruptive component of the inrush current is a unidirectional current impulse that appears when the transformer core is saturated and this component will be taken into account.
The resistance R of the equivalent circuit diagram of the HTS transformer ( Figure 1) depends on the state of the primary superconductor winding. We assume that the transition characteristics of the superconducting windings are as shown in Figure 3. For the calculations it has been assumed that the resistance R in the superconductive state has a value greater than zero. Such an approach allows for taking into account the residual resistance of the superconductor, connection resistances, current bushings resistances and the power network resistance. It is assumed that when the transformer is switched on, its windings are in a superconductive state. It is assumed that the temperature of the windings does not exceed the critical temperature and the effect of the magnetic field strength is ignored for the duration of the inrush current.
In the saturation state of the core, the scheme of an HTS transformer in the idle state ( Figure 1) is simplified to a serial connection of the resistance of the primary winding R and the inductance of that winding L, as shown in Figure 6 [14,15]. The unidirectional current appears when the momentary value of the flux ϕ exceeds the value of the saturation flux of the core Φn, which takes place for ωt = −θ (Figure 7). By converting Equation (13), including ωt = −θ and ϕ = Φn, a formula is obtained for the angle at which the unidirectional current pulse appears: In Figure 7, the transformer is switched on when the voltage e passes through the zero value, i.e., for angle −π. The primary winding of the HTS transformer is in the superconductive state and its resistance is close to zero (R = R0). The flux reaches the saturation value of the transformer core Φs for the angle −θ. At this point, a unidirectional current pulse i appears. The increasing current pulse reaches the value of critical current of the primary winding i = Ic for the angle −θc. At this point, the primary winding of the transformer passes to the resistive state and its resistance increases to R = Rm. The falling When the W switch is closed at t = 0, the circuit from Figure 6 is described by the equation: The unidirectional current appears when the momentary value of the flux Φ exceeds the value of the saturation flux of the core Φ n , which takes place for ωt = −θ (Figure 7). By converting Equation (13), including ωt = −θ and ϕ = Φ n , a formula is obtained for the angle at which the unidirectional current pulse appears: In Figure 7, the transformer is switched on when the voltage e passes through the zero value, i.e., for angle −π. The primary winding of the HTS transformer is in the superconductive state and its resistance is close to zero (R = R0). The flux reaches the saturation value of the transformer core Φs for the angle −θ. At this point, a unidirectional current pulse i appears. The increasing current pulse reaches the value of critical current of the primary winding i = Ic for the angle −θc. At this point, the primary winding of the transformer passes to the resistive state and its resistance increases to R = Rm. The falling unidirectional current pulse reaches the critical current at which the primary winding of the transformer returns to the superconductive state i = Icw for angle δcw. For the δ angle, the transformer core comes out of saturation and the unidirectional current i disappears. As a result of the two changes in the state of the circuit from Figure 6, the unidirectional current pulse i is described by three different equations: (1) for the value of i from −θ to −θc, (2) for the value of ic from −θc to δcw and (3) for the value of icw from δcw to δ. Under initial conditions, for ωt = −θ, the unidirectional current must have a value equal to zero. By substituting i = 0 and ωt = −θ for Equation (14), the formula for the unidirectional component of the inrush current is obtained: In Figure 7, the transformer is switched on when the voltage e passes through the zero value, i.e., for angle −π. The primary winding of the HTS transformer is in the superconductive state and its resistance is close to zero (R = R 0 ). The flux reaches the saturation value of the transformer core Φ s for the angle −θ. At this point, a unidirectional current pulse i appears. The increasing current pulse reaches the value of critical current of the primary winding i = I c for the angle −θ c . At this point, the primary winding of the transformer passes to the resistive state and its resistance increases to R = R m . The falling unidirectional current pulse reaches the critical current at which the primary winding of the transformer returns to the superconductive state i = I cw for angle δ cw . For the δ angle, the transformer core comes out of saturation and the unidirectional current i disappears. As a result of the two changes in the state of the circuit from Figure 6, the unidirectional current pulse i is described by three different equations: (1) for the value of i from −θ to −θ c , (2) for the value of i c from −θ c to δ cw and (3) for the value of i cw from δ cw to δ.
Under initial conditions, for ωt = −θ, the unidirectional current must have a value equal to zero. By substituting i = 0 and ωt = −θ for Equation (14), the formula for the unidirectional component of the inrush current is obtained: where: and If the HTS transformer winding does not come out of the superconduction state, the unidirectional current pulse over its entire duration can be described by Equation (16).
By converting Equation (16) under boundary conditions ωt = δ and i = 0, the formula is obtained for the angle δ at which the unidirectional current disappears: Now let us assume that the primary winding of the HTS transformer comes out of the superconducting state for angle −θ c , as shown in Figure 7. At that point, Equation (16) is no longer correct. Assuming the boundary conditions i = I c and ωt = −θ c and introducing them into Equation (16), the formula for angle θ c is obtained: where parameter A is a given dependence: and B is: By solving the general Equation (14) of the circuit from Figure 6, for the initial conditions i = I c and ωt = −θ c , the equation for the unidirectional current is obtained after the HTS transformer winding has passed to the resistive state: Equation (23) is valid as long as the primary winding of the HTS transformer is in the resistive state and for the longest time until the current pulse disappears.
Let us now assume that the primary winding of the HTS transformer returns to the superconductive state for angle δ cw , as shown in Figure 7. The angle δ cw is determined by substituting i c = I cw and ωt = δ cw for Equation (23) and after transformations, the formula is obtained: where parameter G is a given dependence: and H is: By solving the general Equation (14) of the circuit from Figure 6, for the initial conditions i = I cw and ωt = δ cw , the formula for the unidirectional current is obtained after recovering the state of superconductivity of the transformer winding: For the δ angle ( Figure 7) the unidirectional current pulse fades. The angle δ is determined by substituting i cw = 0 and ωt = δ for formula (27) and after transformations, the relation is obtained: where parameter K is a given dependence: and M is: A block diagram of the algorithm for calculating the turn-on current using the given formulas is shown in Figure 8.

Experimental Verification of Calculations
The verification of the obtained formulas was carried out experimentally on two singlephase 8.5 kVA and 13.8 kVA superconducting transformers, the construction of which is shown in Figures 9 and 10. The rating data of the transformers are given in Table 1. For the core of the 8.5 kVA transformer, a wound and cut core with a rectangular cross-section of 0.0049 m 2 , made of PN ET52-27 sheet metal with the trade symbol RZC-70/230/70, was used. The core of the 13.8 kVA transformer was made as stepped and diagonally braided, with an area of 0.008 m 2 , from M150-30S electrotechnical sheet metal. diagonally braided, with an area of 0.008 m 2 , from M150-30S electrotechnical sheet metal.
The current bushings of the transformers were made of bar-shaped copper elements for the 8.5 kVA transformer and cylindrical for the 13.8 kVA transformer. The current bushings were connected to the superconducting windings through tin solder joints.    The current bushings of the transformers were made of bar-shaped copper elements for the 8.5 kVA transformer and cylindrical for the 13.8 kVA transformer. The current bushings were connected to the superconducting windings through tin solder joints.
The superconducting windings were wound on glass composite carcasses and placed in cryostats made of the same material. The cooling liquid was nitrogen.
All electrical connections inside the cryostat were made with superconducting wires by tin soldering. The superconducting windings were wound on glass composite carcasses and placed in cryostats made of the same material. The cooling liquid was nitrogen.
All electrical connections inside the cryostat were made with superconducting wires by tin soldering.
The primary winding of the 8.5 kVA transformer uses the SCS4050 superconducting tape (SuperPower Inc., Glenville, WV, USA) with a minimum critical current of 115 A, at 77 K and in its own field. For the primary winding of the 13.8 kVA transformer, the SCS4050-AP type tape with a minimum critical current of 87 A. The secondary winding of the 8.5 kVA transformer is made of SCS4050 tape and the 13.8 kVA transformer is made of SCS12050-AP tape with a minimum critical current of 333 A, at 77 K and in its own field. The parameters of the transformer windings used in the measurements are listed in Table 2. The SCS tapes have a layered structure consisting of a 1.6 μm thick (RE)BCO type superconductor and several layers of resistive materials ( Figure 11). The total thickness of the SCS4050 tape is 0.1 mm and the width is 4 mm. The SCS12050-AP differs in width of 12 mm. The winding resistances of the two transformers for different states of operation are given in Table 2.  The primary winding of the 8.5 kVA transformer uses the SCS4050 superconducting tape (SuperPower Inc., Glenville, WV, USA) with a minimum critical current of 115 A, at 77 K and in its own field. For the primary winding of the 13.8 kVA transformer, the SCS4050-AP type tape with a minimum critical current of 87 A. The secondary winding of the 8.5 kVA transformer is made of SCS4050 tape and the 13.8 kVA transformer is made of SCS12050-AP tape with a minimum critical current of 333 A, at 77 K and in its own field. The parameters of the transformer windings used in the measurements are listed in Table 2. The SCS tapes have a layered structure consisting of a 1.6 µm thick (RE)BCO type superconductor and several layers of resistive materials ( Figure 11). The total thickness of the SCS4050 tape is 0.1 mm and the width is 4 mm. The SCS12050-AP differs in width of 12 mm. The winding resistances of the two transformers for different states of operation are given in Table 2. made of SCS12050-AP tape with a minimum critical current of 333 A, at 77 K and in its own field. The parameters of the transformer windings used in the measurements are listed in Table 2. The SCS tapes have a layered structure consisting of a 1.6 μm thick (RE)BCO type superconductor and several layers of resistive materials ( Figure 11). The total thickness of the SCS4050 tape is 0.1 mm and the width is 4 mm. The SCS12050-AP differs in width of 12 mm. The winding resistances of the two transformers for different states of operation are given in Table 2.  The transformers are described in more detail in the articles [49,50]. The data needed for the calculations were obtained from the inrush current measurements described in the article [51]. The measurement system in which the inrush current was measured is shown in Figures 12 and 13 shows the electrical circuit diagram. The current waveform was recorded indirectly by measuring the voltage drop across the shunt. The power circuit parameters are given in Table 3. The resistance of the shunt was 1 mΩ a class of accuracy 0.5. The electrical connections were made with copper wires with a conductor diameter of 6 mm and a total length of 3 m. The resistance of the power cables was 11 mΩ, while the inductance is negligibly small. The power source was a 230 V/0 ÷ 230 V autotransformer with a power output of 23.4 kVA. The resistance of the power source was 1.5 Ω and its impedance was 128 mH. Power on the HTS transformer was performed when the automatic switching system detected the grid voltage passing through zero. The switching circuit was an electronic circuit of our own design. Simultaneously, the data acquisition system using the NI DAQ PCI-6280 card and a PC with an application written in the LabVIEW environment was started.  The transformers are described in more detail in the articles [49,50]. The data needed for the calculations were obtained from the inrush current measurements described in the article [51]. The measurement system in which the inrush current was measured is shown in Figures 12 and 13 shows the electrical circuit diagram. The current waveform was recorded indirectly by measuring the voltage drop across the shunt. The power circuit parameters are given in Table 3. The resistance of the shunt was 1 mΩ a class of accuracy 0.5. The electrical connections were made with copper wires with a conductor diameter of 6 mm and a total length of 3 m. The resistance of the power cables was 11 mΩ, while the inductance is negligibly small. The power source was a 230 V/0 ÷ 230 V autotransformer with a power output of 23.4 kVA. The resistance of the power source was 1.5 Ω and its impedance was 128 mH. Power on the HTS transformer was performed when the automatic switching system detected the grid voltage passing through zero. The switching circuit was an electronic circuit of our own design. Simultaneously, the data acquisition system using the NI DAQ PCI-6280 card and a PC with an application written in the LabVIEW environment was started.    The transformers are described in more detail in the articles [49,50]. The data needed for the calculations were obtained from the inrush current measurements described in the article [51]. The measurement system in which the inrush current was measured is shown in Figures 12 and 13 shows the electrical circuit diagram. The current waveform was recorded indirectly by measuring the voltage drop across the shunt. The power circuit parameters are given in Table 3. The resistance of the shunt was 1 mΩ a class of accuracy 0.5. The electrical connections were made with copper wires with a conductor diameter of 6 mm and a total length of 3 m. The resistance of the power cables was 11 mΩ, while the inductance is negligibly small. The power source was a 230 V/0 ÷ 230 V autotransformer with a power output of 23.4 kVA. The resistance of the power source was 1.5 Ω and its impedance was 128 mH. Power on the HTS transformer was performed when the automatic switching system detected the grid voltage passing through zero. The switching circuit was an electronic circuit of our own design. Simultaneously, the data acquisition system using the NI DAQ PCI-6280 card and a PC with an application written in the LabVIEW environment was started.    Comparisons of the first inrush current pulses shown in Figures 14 and 15 show the discrepancies in the waveforms of the measured and calculated pulses.

Power Source
Supply Line Shunt  Figures 14 and 15 show the discrepancies in the waveforms of the measured and calculated pulses.

Comparisons of the first inrush current pulses shown in
For the 8.5 kVA HTS transformer, the calculated unidirectional current impulse starts approximately 14 degrees earlier and ends approximately 17 degrees after the current pulse obtained from the measurements. The offset between the pulse peaks is approximately 9 degrees. The maximum measured current is 178 A and the calculated current is 178.2 A.
With the HTS 13.8 kVA transformer, the length of both pulses is similar. However, there is a delay, an impulse calculated relative to the measured one, of approximately 18 degrees at the base and 30 degrees at the top. The maximum measured current is 257.2 A and the calculated current is 260.7 A.     Figures 14 and 15 show the discrepancies in the waveforms of the measured and calculated pulses.
For the 8.5 kVA HTS transformer, the calculated unidirectional current impulse starts approximately 14 degrees earlier and ends approximately 17 degrees after the current pulse obtained from the measurements. The offset between the pulse peaks is approximately 9 degrees. The maximum measured current is 178 A and the calculated current is 178.2 A.
With the HTS 13.8 kVA transformer, the length of both pulses is similar. However, there is a delay, an impulse calculated relative to the measured one, of approximately 18 degrees at the base and 30 degrees at the top. The maximum measured current is 257.2 A and the calculated current is 260.7 A.   For the 8.5 kVA HTS transformer, the calculated unidirectional current impulse starts approximately 14 degrees earlier and ends approximately 17 degrees after the current pulse obtained from the measurements. The offset between the pulse peaks is approximately 9 degrees. The maximum measured current is 178 A and the calculated current is 178.2 A.
With the HTS 13.8 kVA transformer, the length of both pulses is similar. However, there is a delay, an impulse calculated relative to the measured one, of approximately 18 degrees at the base and 30 degrees at the top. The maximum measured current is 257.2 A and the calculated current is 260.7 A.
The time for which the 8.5 kVA transformer winding is in the resistive state is 3.6 ms for the results obtained from the measurements and 4.4 ms for the results obtained from the calculations. In the case of a 13.8 kVA transformer, the time for which the primary winding is in the resistive state is 5.5 ms for the results obtained from the measurements and 6.3 ms for the results obtained from the calculations. Figures 16 and 17 compare the first five inrush current pulses of both transformers. Table 4 shows a comparison of the measured and calculated maximum values of the first four inrush current pulses of the discussed transformers. Table 5 shows the duration of the first four current pulses for the transformers and Table 6 shows the durations of the interruptions without current between the individual inrush current pulses. The time for which the 8.5 kVA transformer winding is in the resistive state is 3.6 ms for the results obtained from the measurements and 4.4 ms for the results obtained from the calculations. In the case of a 13.8 kVA transformer, the time for which the primary winding is in the resistive state is 5.5 ms for the results obtained from the measurements and 6.3 ms for the results obtained from the calculations. Figures 16 and 17 compare the first five inrush current pulses of both transformers. Table 4 shows a comparison of the measured and calculated maximum values of the first four inrush current pulses of the discussed transformers. Table 5 shows the duration of the first four current pulses for the transformers and Table 6 shows the durations of the interruptions without current between the individual inrush current pulses.   The time for which the 8.5 kVA transformer winding is in the resistive state is 3.6 ms for the results obtained from the measurements and 4.4 ms for the results obtained from the calculations. In the case of a 13.8 kVA transformer, the time for which the primary winding is in the resistive state is 5.5 ms for the results obtained from the measurements and 6.3 ms for the results obtained from the calculations. Figures 16 and 17 compare the first five inrush current pulses of both transformers. Table 4 shows a comparison of the measured and calculated maximum values of the first four inrush current pulses of the discussed transformers. Table 5 shows the duration of the first four current pulses for the transformers and Table 6 shows the durations of the interruptions without current between the individual inrush current pulses.

Conclusions
The proposed method allows the calculation of the inrush current parameters for superconducting transformers with the knowledge of basic electrical quantities (resistance, inductance, voltage) using relatively simple mathematical formulas. These formulas can be easily used in one's own computer programs for calculating the inrush current of both newly designed and existing transformers. Based on the formulas given, not only the maximum value of the inrush current can be estimated, but also the maximum value of all subsequent impulses of this current. It can be estimated by how much the inrush current of the individual pulses exceeds the critical value for the superconductor from which the transformer windings are made. It can also be estimated for how long this individual pulses overshoot will last, and thus how long the superconducting windings will be in a resistive state. The resistive state of superconducting windings should be treated as a state threatening to break the continuity of windings as a result of their thermal damage.
Verification of the results obtained from the calculations with the results obtained from the measurements was carried out only for two superconducting transformers of 8.5 kVA and 13.8 kVA. Although good compliance of the inrush current waveforms has been achieved, there is no confirmation of equally good compliance for more transformers, especially for medium and high-power transformers.
Differences in measured and calculated current pulses cannot be overlooked. The simplified shape of the magnetisation curve used in the calculations ( Figure 5) translates into differences in the basis of the calculated and actual current pulses. These differences are exacerbated by the omission in the calculations of the steady-state component of the turn-on current, which is in fact the idle current of the transformer. The disadvantage of this method is the lack of consideration of changes in resistance of the transformer winding as a function of temperature changes of the superconductor resulting from the Joule-Lenz phenomenon. This is visible in the symmetry of the calculated impulse of the inrush current, which has the same course of ascending and descending slope, differing from the real course, in which the time of current loss is longer than the time of ascending (Figures 14 and 15). The differences between the turn-on current obtained from the measurements and simulations are also caused by the simplified transition characteristics of the superconductor adopted for the calculations.
The present considerations should make us aware of the differences in the inrush current of conventional and superconducting transformers. The calculation of the inrush current of a superconducting transformer from the formulae valid for conventional transformers must lead to large differences between the measured and calculated values.