An Electromagnetic–Thermal Coupling Numerical Study of the Synchronous Generator with Second-Generation High-Temperature Superconducting Armatures

: Generators with high-temperature superconducting armatures have an advantage in the fact that they can carry high currents. However, the AC loss of high-temperature superconducting (HTS) armatures is di ﬃ cult to calculate precisely because HTS coils exist in a complex and time-varying electromagnetic environment. In addition, when the HTS coil is carrying a short circuit fault overcurrent, an electromagnetic–thermal simulation study of this process is required to ensure that the HTS coil is not damaged. In this paper, ﬁrst, a fully coupled T-A formulation model is used to calculate the AC loss of HTS armatures. Then, the current and temperature distributions are simulated, considering the intrinsic characteristic of superconducting coated conductors, when the generator su ﬀ ers the worst short circuit fault accidently. It is found that the turn with the lowest critical current quenches after 0.01 s, but the temperature rise cannot damage the coil if the circuit breaker can clear the fault quickly. The e ﬀ ects of the copper stabilizer thickness on the thermal stability of the HTS coil during the worst short circuit fault are also investigated. A thicker copper stabilizer improves the thermal stability of the HTS coil in the event of a short circuit fault, but the use of a simulation model is needed to make trade-o ﬀ s between the engineering current density and the thermal stability of the HTS tapes. The work in this paper is necessary and can provide an important reference for manufacturing superconducting generators. the simulation results show that a thicker copper stabilizer improves the thermal stability of the HTS coil in the event of a short circuit fault, but for the HTS generator designer, this requires the use of a simulation model to make trade-o ﬀ s between the engineering current density and the thermal stability of the HTS tapes.


Introduction
Second-generation (2G) high-temperature superconducting (HTS) generators are promising because they can achieve a higher power density compared to conventional generators [1]. For synchronous generator applications, HTS coils placed in the stator have an advantage in terms of static cooling, but HTS coils will suffer AC loss because they are under alternating magnetic fields and currents.
The AC loss of an isolated HTS coil has been studied by many groups [2][3][4][5][6]. To calculate the AC loss of HTS coils in the actual generator environment, a two-stage segregated model consisting of a machine model and an AC loss model has been proposed [7] and used [8,9]. The rotating machine is first simulated by the A-formulation, and then the HTS coil is simulated by the H-formulation [10], with the boundary condition calculated from the A-formulation model.
The T-A formulation model [11,12] can fully couple the A-formulation model of the generator and the HTS coil model. The efficiency of the T-A formulation model is higher than the H-formulation model due to the decrease in the degrees of freedom [11]. Moreover, the T-A formulation model can directly achieve the fully coupled simulation of the HTS coils under the generator's electromagnetic environment because rotating generators are usually simulated by the A-formulation [13].
When a short circuit fault happens, the magnetic field on the surface of the HTS coils and the current flowing in the coils will change drastically and cause a greater AC loss. The ohmic loss also occurs when the current exceeds the critical current of the coils. If the cooling capacity of the cryocooler is not enough to remove the heat immediately, the HTS coils may be permanently damaged [14]. To understand whether the short circuit fault will cause permanent damage to HTS coils, conducting an electromagnetic thermal coupling simulation is of great importance. It should be considered that these factors will have an effect on HTS materials, including the magnetic field, the current and the temperature.
To date, some numerical simulations of the short circuit faults of HTS generators have been done [15][16][17][18][19]. However, none of these simulations achieves fully coupled simulations, considering the dependence of the critical current density on the temperature and the magnetic field.
The main novelty of this paper is to conduct an electromagnetic-thermal coupling numerical study of a synchronous generator with 2G HTS armatures. In Section 2, the field-circuit thermal coupled model is introduced. In Section 3, the parameters of the HTS generator and the HTS coils are introduced. In Section 4, the model is first used to calculate the AC loss of the stator HTS coils under the rated condition of the generator. Then, the transient electromagnetic-thermal behavior of the HTS coils under a sudden three-phase short circuit fault at the no-load operation will be studied. This situation causes a maximum peak short circuit current, which is the worst case for the HTS coil. Moreover, the effects of the copper stabilizer thickness on the thermal stability of the HTS coil at the worst short circuit fault will also be investigated. Finally, our conclusions are drawn in Section 5.

Formulation and Model
To consider the complex anisotropy of YBCO tapes in our numerical model, the dependence of the critical current density on the magnetic flux density B and the temperature T is expressed as [4,[20][21][22]: where B par and B per represent the magnetic field component parallel and perpendicular to the tape surface, respectively. J c0 is the self-field critical current density at the operating temperature. b, k and Bc are the curve-fitting parameters. T c = 92 K is the critical temperature. T 0 = 77 K is the operating temperature. It is important to note that Equations (2) and (3) are only an approximate expression of the dependence of the critical current density on the magnetic flux density and the temperature. The actual properties of the HTS tapes will not follow these expressions completely.
The superconducting layer resistivity is calculated by the E-J power law relation [23]:

The Field Model
The current density vector potential T is defined in the superconducting layer to solve the current density distribution in HTS-coated conductors: As can be seen in Figure 1, set the current density vector potential of the endpoint on one side of each HTS tape as zero, and the boundary condition for the 2-D T-formulation is: where I su is the superconducting layer current of the HTS coil.

The Field Model
The current density vector potential T is defined in the superconducting layer to solve the current density distribution in HTS-coated conductors: As can be seen in Figure 1, set the current density vector potential of the endpoint on one side of each HTS tape as zero, and the boundary condition for the 2-D T-formulation is: where Isu is the superconducting layer current of the HTS coil. In all domains of the generator, the magnetic vector potential A is defined and the A-formulation is used: superconucting layers 1 () copper layers where Jcu is the copper layer current density of the HTS coil, calculated by coupling with the circuit part (introduced in Section 2.2). μ is the permeability of the material. The relative permeability for all the materials in the model is one. The boundary condition for the A-formulation is set as: In all domains of the generator, the magnetic vector potential A is defined and the A-formulation is used: where J cu is the copper layer current density of the HTS coil, calculated by coupling with the circuit part (introduced in Section 2.2). µ is the permeability of the material. The relative permeability for all the materials in the model is one.
The boundary condition for the A-formulation is set as:

The Circuit Model
The equivalent circuit model of the whole superconducting synchronous generator is shown in Figure 2. For the stator windings, since the three-phase stator windings are symmetrical, only the variables related to the phase-A winding are explained.

The Circuit Model
The equivalent circuit model of the whole superconducting synchronous generator is shown in Figure 2. For the stator windings, since the three-phase stator windings are symmetrical, only the variables related to the phase-A winding are explained. The circuit equation of the phase-A is: where UA, LLA, RLA, iA are the output voltage, the load inductance, the load resistance and the current of the phase-A winding, respectively. For the HTS winding, it can be seen as a series of turns (right part of the Figure 2), and the circuit equations are [24]: where the subscript k represents the k-th turn of the HTS coil. N is the total number of turns of the coil. ek is the generated electromotive force. Ez is the electric field strength in the coil domain. Sk is the cross-sectional area of the conductor. L is the core length. ik is the total current. i, R and T represent the current, the resistance and the temperature, respectively. The subscript su and cu represent the superconducting layer and copper layer, respectively. Bsu,k is the magnetic flux density in the superconducting layer. lk is the tape length. ρcu is the copper layer resistivity which is from [25]. Scu,k is the cross-section area of the copper layer.

The Thermal Model
The governing equation of the heat transfer model in the HTS coil is [25,26]: The circuit equation of the phase-A is: where U A , L LA , R LA , i A are the output voltage, the load inductance, the load resistance and the current of the phase-A winding, respectively. For the HTS winding, it can be seen as a series of turns (right part of the Figure 2), and the circuit equations are [24]: where the subscript k represents the k-th turn of the HTS coil. N is the total number of turns of the coil. e k is the generated electromotive force. E z is the electric field strength in the coil domain. S k is the cross-sectional area of the conductor. L is the core length. i k is the total current. i, R and T represent the current, the resistance and the temperature, respectively. The subscript su and cu represent the superconducting layer and copper layer, respectively. B su,k is the magnetic flux density in the superconducting layer. l k is the tape length. ρ cu is the copper layer resistivity which is from [25]. S cu,k is the cross-section area of the copper layer.

The Thermal Model
The governing equation of the heat transfer model in the HTS coil is [25,26]: where d, C and k are the mass density, the heat capacity and the thermal conductivity of the materials of each layer of the HTS coil, respectively. The material properties of the Kapton layer are from [26] and other layers are from [25]. The magnetization loss density Q su is generated in the superconducting layers. The ohmic loss density Q cu is generated in the copper layer, and it is caused by the shunting of the copper layer, which will occur after the total current exceeds the critical current of the HTS tapes.
The distribution of the same kind of loss in each turn is non-uniform because, in the current distribution, the magnetic field distribution and the temperature distribution of each turn are different. The values of these heat sources for each turn need to be calculated by coupling with the field model and the circuit model introduced previously: During the operation of the HTS machine, the cryocooler will take away a certain amount of heat in coils. However, in order to estimate the most serious consequences of short circuit faults conservatively, the adiabatic boundary condition is set on the surface of the coils.
The coupling relationship among the field, circuit and thermal model of the HTS coils is summarized in Figure 3: where d, C and k are the mass density, the heat capacity and the thermal conductivity of the materials of each layer of the HTS coil, respectively. The material properties of the Kapton layer are from [26] and other layers are from [25]. The magnetization loss density Qsu is generated in the superconducting layers. The ohmic loss density Qcu is generated in the copper layer, and it is caused by the shunting of the copper layer, which will occur after the total current exceeds the critical current of the HTS tapes.
The distribution of the same kind of loss in each turn is non-uniform because, in the current distribution, the magnetic field distribution and the temperature distribution of each turn are different. The values of these heat sources for each turn need to be calculated by coupling with the field model and the circuit model introduced previously: (14) During the operation of the HTS machine, the cryocooler will take away a certain amount of heat in coils. However, in order to estimate the most serious consequences of short circuit faults conservatively, the adiabatic boundary condition is set on the surface of the coils.
The coupling relationship among the field, circuit and thermal model of the HTS coils is summarized in Figure 3: The field model uses the temperature distribution calculated from the thermal model, and the current values of the superconducting layers and copper layers calculated from the circuit model, to calculate the magnetization loss of the superconducting layers, the magnetic field distribution, and the generated electromotive force of the HTS coils. The circuit model uses the temperature distribution calculated from the thermal model, and the magnetic field distribution calculated from the field model, to calculate the current values of the superconducting layers and copper layers and the copper layer loss.
All these models are implemented in the COMSOL Multiphysics software.

Parameters
The topology of the studied HTS generator is shown in Figure 1, each phase is composed of two coils in parallel. For example, the phase A winding consists of two coils: A1 and A2. A1+ and A1− are the positive current domain and the negative current domain of the A1 coil, respectively.
The studied generator has four poles and six slots. Each slot has one double racetrack coil with 84 turns. The rated current of these HTS field coils is about 37 A and the rated output power is 30 Kw. The generator is designed to operate at 1500 rpm, corresponding to a stator armature frequency of 50 Hz. Other design parameters for the generator can be found in a previous paper [27]. The parameters of the HTS coils are summarized in Table 1. b, k and Bc are calculated by fitting the The field model uses the temperature distribution calculated from the thermal model, and the current values of the superconducting layers and copper layers calculated from the circuit model, to calculate the magnetization loss of the superconducting layers, the magnetic field distribution, and the generated electromotive force of the HTS coils. The circuit model uses the temperature distribution calculated from the thermal model, and the magnetic field distribution calculated from the field model, to calculate the current values of the superconducting layers and copper layers and the copper layer loss.
All these models are implemented in the COMSOL Multiphysics software.

Parameters
The topology of the studied HTS generator is shown in Figure 1, each phase is composed of two coils in parallel. For example, the phase A winding consists of two coils: A1 and A2. A1+ and A1− are the positive current domain and the negative current domain of the A1 coil, respectively.
The studied generator has four poles and six slots. Each slot has one double racetrack coil with 84 turns. The rated current of these HTS field coils is about 37 A and the rated output power is 30 kW. The generator is designed to operate at 1500 rpm, corresponding to a stator armature frequency of 50 Hz.
Other design parameters for the generator can be found in a previous paper [27]. The parameters of the HTS coils are summarized in Table 1. b, k and B c are calculated by fitting the experimental data of a sample HTS tape. The surface plot of the normalized critical current density for the HTS tape under a parallel and perpendicular magnetic field is shown in Figure 4. experimental data of a sample HTS tape. The surface plot of the normalized critical current density for the HTS tape under a parallel and perpendicular magnetic field is shown in Figure 4.  Since the six coils are symmetrical, only the results of the A1 coil are shown in this paper. The numbering sequence of each turn of the studied coil is shown in Figure 1, and this arrangement will be used throughout Section 4. Although this paper is only a numerical simulation study, it is worth mentioning that our group has now manufactured the stator of the generator, which is shown in Figure 5. Since the six coils are symmetrical, only the results of the A1 coil are shown in this paper. The numbering sequence of each turn of the studied coil is shown in Figure 1, and this arrangement will be used throughout Section 4. Although this paper is only a numerical simulation study, it is worth mentioning that our group has now manufactured the stator of the generator, which is shown in Figure 5.

AC Loss of HTS Coils under the Rated Operation
The induced AC loss on HTS armatures will directly affect the efficiency of the generator and thereby increase the cost of cooling. Hence, accurately and quickly estimating the AC loss of the HTS generator is urgent and beneficial. The T-A formulation model has offered a convenient tool to solve this problem. In this section, the AC loss of HTS coils in the generator during the rated operation will be calculated and discussed.
The calculated AC loss values of several representative turns, which are placed at different locations in the HTS coil, are summarized in Table 2. It can be seen that different turns of the coil have significantly different AC loss values. The calculated current waveform of the A1 coil in one cycle at the rated condition is shown in Figure 6, and four typical time moments, A, B, C and D, are chosen to show the calculation results.

AC Loss of HTS Coils under the Rated Operation
The induced AC loss on HTS armatures will directly affect the efficiency of the generator and thereby increase the cost of cooling. Hence, accurately and quickly estimating the AC loss of the HTS generator is urgent and beneficial. The T-A formulation model has offered a convenient tool to solve this problem. In this section, the AC loss of HTS coils in the generator during the rated operation will be calculated and discussed.
The calculated AC loss values of several representative turns, which are placed at different locations in the HTS coil, are summarized in Table 2. It can be seen that different turns of the coil have significantly different AC loss values. The calculated current waveform of the A1 coil in one cycle at the rated condition is shown in Figure 6, and four typical time moments, A, B, C and D, are chosen to show the calculation results.

AC Loss of HTS Coils under the Rated Operation
The induced AC loss on HTS armatures will directly affect the efficiency of the generator and thereby increase the cost of cooling. Hence, accurately and quickly estimating the AC loss of the HTS generator is urgent and beneficial. The T-A formulation model has offered a convenient tool to solve this problem. In this section, the AC loss of HTS coils in the generator during the rated operation will be calculated and discussed.
The calculated AC loss values of several representative turns, which are placed at different locations in the HTS coil, are summarized in Table 2. It can be seen that different turns of the coil have significantly different AC loss values. The calculated current waveform of the A1 coil in one cycle at the rated condition is shown in Figure 6, and four typical time moments, A, B, C and D, are chosen to show the calculation results.  6. Current waveform of the studied A1 coil in one cycle. Figure 6. Current waveform of the studied A1 coil in one cycle.
The critical current of the HTS tapes is mainly influenced by the magnetic field perpendicular to the tape surface. Figure 7 shows the perpendicular magnetic flux density magnitude distributions of the area around the HTS tapes. First, due to the lower coil being closer to the rotor than the upper coil, its surface magnetic field is significantly larger than that of the upper coil. The turn in the lower coil has a greater AC loss value than the corresponding turn with the same distance to the stator iron teeth in the upper coil (e.g., turn 42 and turn 84). In addition, the perpendicular magnetic field in region A1+ is much greater than that in the A1− region. Although the average perpendicular magnetic field of turn 42 in the A1− region is less than turn 21, the average perpendicular magnetic field of turn 42 in the A1+ region is much greater than turn 21, which makes the AC loss of turn 42 greater than turn 21.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 8 of 15 The critical current of the HTS tapes is mainly influenced by the magnetic field perpendicular to the tape surface. Figure 7 shows the perpendicular magnetic flux density magnitude distributions of the area around the HTS tapes. First, due to the lower coil being closer to the rotor than the upper coil, its surface magnetic field is significantly larger than that of the upper coil. The turn in the lower coil has a greater AC loss value than the corresponding turn with the same distance to the stator iron teeth in the upper coil (e.g., turn 42 and turn 84). In addition, the perpendicular magnetic field in region A1+ is much greater than that in the A1− region. Although the average perpendicular magnetic field of turn 42 in the A1− region is less than turn 21, the average perpendicular magnetic field of turn 42 in the A1+ region is much greater than turn 21, which makes the AC loss of turn 42 greater than turn 21.  Figure 8 shows the normalized current density distribution of the HTS tapes. The normalized current density distribution is similar to the perpendicular magnetic field distribution. The critical current degradation is more severe in each turn of the lower coil due to the greater perpendicular magnetic field, so the current penetration into the interior of the HTS tapes is more severe, making the AC loss greater than the upper coil. Turn 42 in the lower coil suffers the maximum surface perpendicular magnetic field and has the minimum critical current because it is nearest to the iron teeth, which attracts the flux lines. Its current penetration is also the most severe and it has the greatest AC loss value.  Figure 8 shows the normalized current density distribution of the HTS tapes. The normalized current density distribution is similar to the perpendicular magnetic field distribution. The critical current degradation is more severe in each turn of the lower coil due to the greater perpendicular magnetic field, so the current penetration into the interior of the HTS tapes is more severe, making the AC loss greater than the upper coil. Turn 42 in the lower coil suffers the maximum surface perpendicular magnetic field and has the minimum critical current because it is nearest to the iron teeth, which attracts the flux lines. Its current penetration is also the most severe and it has the greatest AC loss value.
The generator has three phases and is symmetrical. Since the current and magnetic field of each coil have three phases and are symmetrical, with only the phase angle being different, the AC loss is the same for all six coils. The calculated total AC loss of the 6 HTS coils is 229.38 W. Considering that the cooling penalty is about 10 [28], the loss ratio is about 5.29%.

Transient Electromagnetic-Thermal Behavior of HTS Coils under the L-L-L Fault
In this section, the worst short circuit fault in the generator output terminal is studied: a sudden three-phase short circuit fault at the no-load operation. When a short circuit fault occurs at a time when the electromotive force of the phase-A is passing through zero, the maximum fault current occurs in phase A [16], which is the worst scenario for the A1 HTS coil studied in this paper. The generator has three phases and is symmetrical. Since the current and magnetic field of each coil have three phases and are symmetrical, with only the phase angle being different, the AC loss is the same for all six coils. The calculated total AC loss of the 6 HTS coils is 229.38 W. Considering that the cooling penalty is about 10 [28], the loss ratio is about 5.29%.

Transient Electromagnetic-Thermal Behavior of HTS Coils under the L-L-L Fault
In this section, the worst short circuit fault in the generator output terminal is studied: a sudden three-phase short circuit fault at the no-load operation. When a short circuit fault occurs at a time when the electromotive force of the phase-A is passing through zero, the maximum fault current occurs in phase A [16], which is the worst scenario for the A1 HTS coil studied in this paper.
When the fault happens, the generator short circuit fault protection system will act to cut off the fault current to remove the fault. Therefore, it is only necessary to pay attention to whether the short circuit fault will cause permanent harm to the superconducting coil within a short period of time after the fault occurs.
The short circuit fault occurs at t = 0.04 s of the simulation and the circuit breaker can clear the fault in 0.08 s [17], which means that longer simulation is unnecessary. Figure 9 depicts the variations in the total current, the critical current (CC), the superconducting layer current (SLC) and the copper layer current (CLC) of turn 42 and turn 43 of the A1 coil before and after the occurrence of the worst short circuit fault. The scale of the x-axis has been changed to make it easier for the reader to observe the waveform at 0.04-0.06 s. Figure 10   When the fault happens, the generator short circuit fault protection system will act to cut off the fault current to remove the fault. Therefore, it is only necessary to pay attention to whether the short circuit fault will cause permanent harm to the superconducting coil within a short period of time after the fault occurs.
The short circuit fault occurs at t = 0.04 s of the simulation and the circuit breaker can clear the fault in 0.08 s [17], which means that longer simulation is unnecessary. Figure 9 depicts the variations in the total current, the critical current (CC), the superconducting layer current (SLC) and the copper layer current (CLC) of turn 42 and turn 43 of the A1 coil before and after the occurrence of the worst short circuit fault. The scale of the x-axis has been changed to make it easier for the reader to observe the waveform at 0.04-0.06 s. Figure 10 depicts the variations in the ohmic loss values and temperature values of turn 42 and turn 43 after the occurrence of the worst short circuit fault.
Turn 42 and turn 43 are chosen because the former has the lowest critical current and the most obvious temperature rise in each turn of the coil, while the latter has the highest critical current. As will be shown below, although the total current is the same for both them, during a short circuit fault, due to their different magnetic field environments, their critical current values are different. This difference in the critical current values at the beginning results in their different subsequent copper layer shunting and temperature rise phenomena in each HTS tape.
It can be clearly found that the total current increases rapidly after the occurrence of the short circuit fault. The critical current of turn 42 is much lower than that of turn 43 because of the different magnetic field environments on their surfaces. When the total current first exceeds the CC of turn 42 (0.0428 s), the excess current flows through the copper layer, producing the ohmic loss and raising the temperature. However, turn 43 is still in the superconducting state. As the short circuit current increases further, the total current exceeds the CC of turn 43 (0.0436 s); due to its high critical current, the copper layer current of turn 43 is less than turn 42 s, resulting in less ohmic loss and a lower temperature rise.     For turn 43, due to the non-periodic current component in the short circuit transient process, a peak value (0.0471 s) appears in the first half-cycle after the short circuit fault happens. With this component gradually decaying to zero, the waveform will eventually stabilize into a steady state. Since there is not much of a temperature rise for turn 43 after the first and also the highest peak current (0.0471 s), the subsequent smaller current peaks (for example, the 0.0671 s moment) also cannot induce much of a temperature rise in turn 43.
As for turn 42, after the 0.0471 s moment, although the total current is decreasing, the CC is still decreasing as well and the CLC can still increase for a short time. This phenomenon is caused by the temperature rise due to the ohmic loss. The higher temperature means lower CC, making the copper layer shunt phenomenon more severe and accelerating the temperature rise. After the 0.0495 s moment, the temperature of turn 42 rises to the critical temperature. Turn 42 loses its superconducting characteristic, and all the short circuit currents go through the copper layer.
A large amount of ohmic loss in a short time may cause damage to the HTS coil and this phenomenon is worthy of vigilance. One solution is to design a HTS coil with a greater margin to get through the short time overcurrent process. However, this method is not economical. Another method is to allow the current to exceed the critical current value in a short time as long as the temperature rise is limited and the HTS coil can return to the superconducting state after the fault disappears [29].
For the whole HTS coil, although the temperature rise in the turn with the lowest critical current is obvious, the highest temperature at t = 0.12 s is not high enough to exceed the melting points of the materials of the HTS coil. As long as the machine short circuit protection device can remove the fault, the HTS coil will not be damaged.

Effects of the Copper Stabilizer Thickness on the Thermal Stability of the HTS Coil
Copper stabilizer can improve the thermal stability of the 2G HTS tapes [30]. In this section, the effects of the copper stabilizer thickness on the thermal stability of the HTS coil during the worst short circuit fault are numerically studied. Figure 11 depicts the variations in the total current, the critical current and the copper layer current of turn 42 after the occurrence of the worst short circuit fault, with copper stabilizer thicknesses of 10 µm and 20 µm. Figure 12 depicts the ohmic loss values and the temperature values of turn 42 after the occurrence of the worst short circuit fault with copper stabilizer thicknesses of 10 µm and 20 µm.
It is found that, as the copper stabilizer thickens, the temperature rise phenomenon becomes weaker. The temperature rise in the copper layer in a short period can be approximately considered to be proportional to the ohmic loss volume density (ignore the heat transfer from the copper layer to other layers to simplify the explanation). The ohmic loss volume density is proportional to the square of the current density of the copper layer, so the temperature rise in the copper layer can be approximately regarded as inversely proportional to the square of the thickness of the copper stabilizer. Thus, although the time when the total current exceeds the critical current is almost the same for the two cases, the final temperature will be lower for the thicker copper stabilizer.
However, though the thicker copper stabilizer improves the thermal stability of the HTS coil in the event of a short circuit fault, the engineering current density also decreases. For the HTS generator designer, this requires the use of an electromagnetic-thermal coupling numerical simulation model to make trade-offs between the engineering current density and the thermal stability of the HTS tapes.     One noteworthy phenomenon is that the total current eventually becomes inconsistent after a short circuit fault occurs in both cases. The total current in both cases is the same at the beginning of the short circuit fault because the coil is initially in the superconducting state and the coil resistance is almost zero. The short circuit impedance is mainly the inductive resistance of the coil, which is the same for the two cases. Then, as the temperature rises, the coil resistance is gradually dominated by the copper layer resistance and the short circuit impedance becomes greater. For the thicker copper stabilizer case, the copper layer resistance is smaller, so the short circuit impedance is smaller, causing a greater short circuit current. Moreover, since the copper layer loss is proportional to the square of the current, finally, at t = 0.12 s, the 20-µm copper stabilizer case's ohmic loss even exceeds that of the 10-µm copper stabilizer case by a little.

Conclusions
In this paper, an electromagnetic-thermal coupling numerical study of a synchronous generator with 2G HTS armatures is conducted. The main results in this paper are summarized as follows: First, by coupling the field part model with the circuit part model, the AC loss of HTS coils during rated operation is calculated. Second, the transient electromagnetic-thermal behavior of the HTS coil under the worst short circuit fault is studied. The results show that although the total current is the same for each turn of the HTS coil during the short circuit fault, the critical current value of each turn is different due to their diverse surrounding magnetic fields. This difference in the critical current value at the beginning leads to their subsequently different copper layer shunting and temperature rise phenomena. The turn with the lowest critical current will eventually quench completely, and the temperature rise is obvious, but it is not high enough to exceed the melting points of the materials of the HTS coil. It is concluded from this work that, if the generator short circuit protection device can remove the fault as soon as possible, the HTS coil in this high-speed HTS synchronous generator will not suffer unrecoverable quench and be damaged. Last but not least, the simulation results show that a thicker copper stabilizer improves the thermal stability of the HTS coil in the event of a short circuit fault, but for the HTS generator designer, this requires the use of a simulation model to make trade-offs between the engineering current density and the thermal stability of the HTS tapes.