Modeling and Dynamics of HTS Motors for Aircraft Electric Propulsion

In this paper, the methodology of how a dynamic model of a conventional permanent magnet synchronous motor (PMSM) may be modified to model the dynamics of a high-temperature superconductor (HTS) machine is illustrated. Simulations of a typical PMSM operating under room temperature conditions and also at temperatures when the stator windings are superconducting are compared. Given a matching set of values for the stator resistance at superconducting temperature and flux-trapped rotor field, it is shown that the performance of the HTS PMSM is quite comparable to a PMSM under normal room temperature operating conditions, provided the parameters of the motor are appropriately related to each other. From these simulations, a number of strategies for operating the motor so as to get the propeller to deliver thrust with maximum propulsive efficiency are discussed. It is concluded that the motor–propeller system must be operated so as to deliver thrust at the maximum propulsive efficiency point. This, in turn, necessitates continuous tracking of the maximum propulsive efficiency point and consequently it is essential that the controller requires a maximum propulsive efficiency point tracking (MPEPT) outer loop.


Introduction
High-temperature superconductors were discovered in 1986 by two physicists, Johannes Georg Bednorz and Karl Müller, working at the International Business Machines (IBM) Labs.Bednorz and Müller [1] identified superconductivity in a sample of LaBaCuO material at the much higher temperature of 36 K. Subsequently, high-temperature superconductivity was identified in structured materials containing planes of copper oxide, such as YBa 2 Cu 3 O 7 , by a group led by Chu [2].These materials had a transition temperature beginning at 93 K and had zero resistance at 80 K.The phenomenal growth in the research and development of devices and applications based on this effect have been documented by Ford and Saunders [3].Yet, there are indeed a number of key issues that must be resolved before useful application products can be designed and mass produced.These issues include the problem of flux creep, weak links, and poor mechanical properties, which have also been identified by Goyal [4] and are being resolved.
An electrical motor can be termed as high temperature superconducting if its design contains elements that are high-temperature superconductors and these elements are working in the superconducting state at much higher temperatures than elements exhibiting the normal feature of electrical conduction.The motor usually contains a high-temperature superconducting element, such as an armature or a field winding or both.It is common practice to use a high current density and low loss high-temperature superconductor (HTS) wire for winding the armature or field coils or both.Although conventional superconducting materials with characteristically low resistance to electron flow were first discovered by the Dutch physicist Heike Kammerlingh Onnes in 1911, it was thought to have been completely explained in 1957 by the Bardeen, Cooper, and Schrieffer (BCS) theory [5,6] based on quantum mechanics.The BCS theory completely failed to predict the existence of HTS, although it did provide a framework for characterizing HTS.A superconducting material is in the superconducting state when a point determined by the operating temperature, current density, and magnetic field characteristics lies below a surface defined by the material's critical characteristics.Critical surface characteristics are defined in terms of three variables: temperature, current density, and magnetic field intensity.Thus, superconductivity is generally observed below a critical temperature at a given critical (magnetic) field strength and a critical current density.A superconductor needs to be kept below a critical temperature T c and the applied field must be less than the critical field H c .If an applied magnetic field is stronger than the critical field, it can penetrate the superconductor, causing a quenching of the superconducting state.Even though the temperature may be below the critical temperature, it will no longer be superconducting.Materials exhibiting this behavior are known as type I superconductors.However, HTS materials are known as type II, which exhibit two critical magnetic field states before losing the property of superconductivity.While the most well-known characteristic of superconductivity is zero resistance, the magnetic field inside a superconductor, which is proportional to the sum of the magnetization and the applied field, adds up to zero.The exclusion of magnetic fields within a superconductor is the defining characteristic of a superconductor and creates a phenomenon known as the Meissner effect.Thus, the disappearance of the electrical resistance and the Meissner effect when the magnetic field inside the superconductor is shielded by a lossless current flowing in a thin layer on the surface of the superconducting material (the perfect case of diamagnetism, which occurs when an external magnetic field penetrates only a finite, outer region of the material and does not influence the remaining inner domain of the material) are the fundamental properties of superconductors.Furthermore, in a cooled superconductor after the magnetization and the applied field add up to zero, if the applied field is removed, the magnetization remains and the internal field is now acting in a direction opposite to the originally applied field.Thus, the flux is essentially "trapped" within the superconductor.These superconductors are referred to as flux pinned or trapped superconductors and the trapped field could be as high as 15+ Tesla.The trapped field can be maintained, in principle, by a quantum process known as flux pumping whereby small amounts of magnetic flux are pumped in by the repeated application of magnetic transients, which require relatively low power.
Yttrium barium copper oxide (YBCO) is a family of crystalline chemical compounds known for exhibiting high-temperature superconductivity.This family of superconductors includes the first material that was discovered to become superconducting above the boiling point of liquid nitrogen (77 K) at about 93 K. Typically, YBCO has a layered structure consisting of copper-oxygen planes with yttrium and barium atoms in the crystal structure as well.The resulting crystal structure is similar to a perovskite with a unit cell consisting of stacked cubes of BaCuO 3 and YCuO 3 .YBCO conductors are manufactured as wire assemblies and as tapes that are now available from several industrial suppliers and are very promising conductors for the design of HTS electric motors and their applications.Many YBCO compounds have the general formula YBa 2 Cu 3 O 7−x (also known as Y123), although materials with other Y:Ba:Cu ratios exist, such as YBa 2 Cu 4 O y (Y124) and Y 2 Ba 4 Cu 7 O y (Y247).The bismuth strontium calcium copper oxide superconductor family (BSCCO tapes)-with the following compounds: Bi-2201 or Bi 2 Sr 2 CuO 6 , Bi-2212 or Bi 2 Sr 2 CaCu 2 O 8 , Bi-2223 or Bi 2 Sr 2 Ca 2 Cu 3 O 10 -have also been used in the construction of HTS motors.The first HTS motors began appearing just before the end of the last millennium and in the following year [7][8][9][10][11][12].Kalsi [13] has described in some detail the different types of superconducting machines and has provided models that could be used for simulation.

Architectures of HTS Machines for Hybrid and All-Electric Propulsion
Aircraft propulsion systems can benefit from electric HTS machines in two distinct ways.First, power can be generated by a conventional gas turbine and converted into electrical energy in an HTS generator, which in turn is used to power a controlled HTS motor driving a propeller.Due to the onboard availability of electric power, an electric motor may be used to drive the compressor in a high bypass ratio gas "turbine" engine, thus making the turbine redundant.Thus, such a hybrid propulsion system uses both HTS generators and motors.This turbo-generator motor drive has several advantages over a direct drive gas turbine propulsion system.The turbine main shaft speed is not limited by the fan and, therefore, can be operated at higher efficiencies.The use of a turbo-generator allows for greater redundancy, which in turn improves the overall reliability of the propulsion system.Moreover, a turbo-generator motor drive facilitates greater control of thrust and consequently makes propulsion-based flight control feasible, particularly when the gas turbine can be distributed along the lifting surfaces.
On the other hand, all-electric would use only HTS motors without the need for an onboard generator and would use only batteries or fuel cells to drive the propulsion HTS motors.A wide variety of HTS reluctance motors with diamagnetic rotors, HTS permanent magnets, and wound field synchronous motors and HTS brushless Direct Current (DC) motors are being developed for aircraft and other applications [14][15][16][17][18][19][20].Thus, several superconducting field synchronous machines with room temperature or cryogenic high purity metal armatures are under development.In these applications, it is advantageous to reduce the number of moving parts in the HTS machines.For this reason, for most aircraft drive systems, rotary machines with HTS stator windings and permanent magnets or wound and trapped-field rotors are preferred.The stator is generally constructed by winding HTS tape, about 4-6 mm wide, flat on spools of the same width, over each of the poles, which are quite large.Such windings are referred to as pancake coils.Alternately, the windings are inserted into circular tracks on a disc and these types of windings are referred to as racetrack windings.
Most superconducting machines have employed liquid neon or cold gaseous helium for cooling HTS coils.Liquid nitrogen-cooled HTS motors have also been developed for ship propulsion by Okazaki, Sugimoto, and Takeda [21].The coolant is generally transferred from an external refrigerator to the rotor through a rotating coupling.The field winding consists of several HTS coils that are conduction-cooled through the support structure or exposed to the tubes carrying the coolant to enhance the cooling rate.The torque tube is mainly responsible for transferring torque from the "cold" (cryogenically-cooled) environment to the "warm" shaft end.The pole pairs and the support structure are enclosed in the vacuum-sealed cryostat to minimize heat input due to radiation and to provide an insulated operating environment for the HTS field coils.A refrigeration system or cylinder, which uses the cold circulating gas in a closed loop, continuously maintains the HTS field winding at the desired low, cryogenic temperature.
One approach to keeping the field coils magnetized is to use flux-trapped superconductors to maintain the flux in the field coils.Such machines are modeled as high-field and low-resistance machines, but the torque constant is relatively high.They need minimal energy for super-cooling over an entire flight and continuous magnetization is provided by periodic flux pumping.
In this paper, the methodology of how a dynamic model of a conventional permanent magnet synchronous motor (PMSM) may be modified to model the dynamics of an HTS machine is illustrated.Simulations of a typical PMSM operating under room temperature conditions and at temperatures when the stator windings are superconducting are compared.Given a matching set of values for the stator resistance at superconducting temperature and flux-trapped rotor field, it is shown that the performance of the HTS PMSM is quite comparable to a PMSM under normal room temperature operating conditions, provided the parameters of the motors are appropriately related to each other.From these simulations, several strategies for operating the motor to get the propeller to deliver thrust with maximum propulsive efficiency are discussed.Thus, it is concluded that the motor-propeller system must be operated to deliver thrust at the maximum propulsive efficiency point.This, in turn, necessitates the continuous tracking of the maximum propulsive efficiency point and, consequently, it is essential that the controller requires a maximum propulsive efficiency point tracking (MPEPT) outer loop.

Dynamics of HTS Motors
HTS motors-that is, motors with armature and field windings made of HTS wires or tapes that are encased within a cryostat so that the temperature of the stator and rotor are in the superconductivity domain-may be modeled just as conventional motors.However, the values of the resistances, torque, and back Electro-Motive Force (EMF) constants are substantially different from the corresponding values in a conventional motor.The models of DC, brushless DC (BLDC), and synchronous motors are briefly reviewed in the following subsections.

Dynamics of an HTS DC Motor
The mechanical shaft/armature dynamics of a DC motor are: where ω m is the mechanical shaft speed, θ m is the mechanical shaft angular position, J is the shaft polar moment of inertia, B is the viscous friction coefficient in the bearings, K T is the driving electrical torque constant, and T L is the load torque.Additionally, i drive = i a is the rotor (armature) current, L a is the armature circuit equivalent inductance, and R a is the armature circuit equivalent resistance.The armature equivalent circuit (current) dynamics are: where V in is the input control voltage and K B is the back EMF constant.Ideally, K T = K B .In the case of an HTS motor, R a is considerably reduced while K T (and consequently K B ) is considerably increased.

Dynamics of an HTS Brushless DC Motor
A typical star connected three-phase BLDC motor is considered.The BLDC three-phase motor is driven by electronically switching on two of the three phases at a time, either in the forward or in the reverse direction while the remaining third phase is switched off.
The mechanical shaft/armature dynamics are: where θ m is the mechanical shaft angular position, θ e is the electrical shaft position, P is the number of pole pairs, ω m is the mechanical shaft speed, J is the shaft polar moment of inertia, B is the viscous friction coefficient in the bearings, T drive is the driving electrical torque, and T L is the load torque.Additionally, i ab is the stator phase differential current, v ab is the stator phase differential voltage, e ab is the stator phase differential back EMF, L s is the stator phase circuit equivalent inductance, and R s is the stator phase circuit equivalent resistance.Thus, for a three-phase machine with three stator pole pairs, The back EMF in each phase is modeled as (Rambabu [22]) The drive torque is where each of the functions F j (θ e ) is a trapezoidal function of θ e given by F 1 (θ e ) = 6θ e /π when 0 ≤ θ e ≤ π/6, F 1 (θ e ) = 1 when π/6 ≤ θ e ≤ 5π/6, F 1 (θ e ) = 6 − 6θ e /π when 5π/6 ≤ θ e ≤ 7π/6, and For a star connected motor, i a + i b + i c = i 0 .Using Equations ( 4) and ( 5), the complete three-phase BLDC motor state-space equations are obtained as follows: or In the case of an HTS motor, R s is considerably reduced while K T (and consequently K B ) is considerably increased.

Dynamics of an HTS Synchronous Motor
In a typical three-phase, two-pole permanent magnet synchronous motor (PMSM), the dynamic modeling begins with the choice of two reference frames.The three phases are labeled 'a', 'b' and 'c' respectively.The stator reference axis for the 'a' phase is generally chosen to the direction of maximum magneto-motive force or MMF, which is analogous to the EMF in electric circuits when a positive 'a'-phase current is supplied at its maximum level.Reference axes for 'b' and 'c' stator phases are chosen 120 • and 240 • (electrical angle) ahead of the 'a'-axis, respectively.Following the universal convention for choosing the rotor reference frame, the direction of permanent magnet flux is chosen as the d-axis, while the q-axis is 90 • ahead of the d-axis.Thus, one may also define a stator d-axis along the direction of the resultant stator field, while the stator q-axis is 90 • ahead of the stator d-axis.The angle of the rotor d-axis with respect to the stator d-axis is defined as θ r .The strength of the permanent magnetic field is defined by ψ f .
From Vepa [23], in the stator dand q-axes, Differentiating the latter with time, d dt Hence, it follows that In the rotor α-and β-axes, (stationary axes aligned with the stator dand q-axes) assuming the rotor position θ m , speed ω m , and electrical frequency ω r dynamic equations are as follows: The equation for the stator currents is as follows: To relate the stator currents to the three phase currents, Clarke's transformation from the stator three phase to the rotor α-and β-axes may be applied and this gives the following: Similar equations hold for the phase voltages and voltages along the α-and β-axes.Applying Parks transformation to transform from the rotor α-and β-axes to the rotor dand q-axes gives the following: Similar equations hold for the voltages along the α-and β-axes and the voltages along the dand q-axes.Using Equations ( 19) and (20) in the rotor's dand q-axes with L α = L β and R α = R β , and with one has for the rotor and current dynamics In general, for any of the DC, BLDC, or synchronous motors, In a typical flux-trapped HTS motor, R is reduced by a significant factor due to the property of superconductivity, while K T and e back are both significantly increased due to the flux trapping.If only R is reduced by a significant factor, the motor would not be very stable in operation.However, because of the introduction of flux trapping, both K T and e back are significantly and proportionately increased and, consequently, the stability of the motor is not degraded.

Propeller Load Torque
The blade element momentum (BEM) theory described in Vepa [23] is briefly summarized.It combines two methods of modeling the operation of a propeller where the first method is based on the balance of momentum within a rotating annular stream tube passing through a propeller disc, while the second method, which is known as the blade element theory, employs the airfoil lift and drag forces along various spanwise sections of the blade.Both methods are utilized to obtain the induced flow velocities across the propeller while only the second method is applied to obtain the blade thrust and blade torque.
In the blade element theory applied to a propeller blade, the blade is modeled by many elements along its length.On each element, following Vepa [23], it is assumed that: (i) Aerodynamic interference between any of the blade elements is absent; (ii) Only the two-dimensional section lift and drag forces are utilized to evaluate the forces on the blade elements.
To apply the BEM theory, following Vepa [23], the blade is divided into N elements along its length.A slightly different flow is experienced at each blade element as they have different translational speeds.Furthermore, the chord lengths and inflow angles are different at each element.Following the division of the blade into a large number (usually about twenty or more) of elements, the blade element theory involves calculating the flow about each of the elements.By numerically integrating the forces and moments along the blade span, performance characteristics of the entire blade, such as torque, thrust and power coefficients, are obtained.
From Vepa [21], the load torque T L due to an N bladed propeller may be written in terms of the blade element momentum (BEM) theory as follows: The propeller blade is divided into a finite number of blade sections, N s , with uniform section properties and each with a length ∆r j and writing the integral as a summation across the length of the blades from the hub to the tip, where at a point r along the radius of the propeller blade with blade chord c, blade length R and a total of N blades, The axial and radial induction factors satisfy the following: where Q tip,j is the Prandtl's tip flow correction factor at a section defined as The airfoil section drag coefficient is known to be composed of the profile and induced drag components and is assumed to be given by an expression representing a drag polar of the form: and where α s is the value of α when stall is initiated and is determined by setting dC L /dα = 0. Furthermore, α j = θ j − φ j , and θ j is the local blade pitch angle (including the blade twist).Differentiating φ with respect to time gives the following: .
Assuming that a j and a j are slowly varying functions and that .

Dynamic Model of Motor-Propeller System
According to the Snel model, the lift coefficient is expressed as follows (see also Han, et al. [24]): In Equation (36), ∆C L 1j , ∆C L 2j satisfy where at a point r j along the radius of the blade with blade chord c and blade length R, The coefficients of the model are defined as follows: and α zj is the blade section zero-lift angle-of-attack.In matrix form: Thus, the propeller load torque is given by Equation ( 2) and the motor dynamics defined by Equation ( 24).

Typical Simulation Results
In this section, typical simulation-based responses of an HTS PMSM-driven propeller are illustrated.The typical parameters of a PMSM-propeller system are listed in Table 1.The parameters corresponding to the case when a motor is within a cryogenic cooler at superconducting temperatures were estimated from the properties of a superconducting HTS tape and flux-trapped magnetic material.The typical properties of a superconducting HTS tape used in HTS PMSM, as well as the field due to flux trapping, were obtained from [10][11][12][13][14]. The stator resistance shown in the table includes the resistance of the source, as well as the stator winding.Initially, the PMSM propeller combination was simulated at normal room temperature with no stator dand q-axes voltage inputs.Thus, the u d and u q voltage inputs in the stator dand q-axes were set to zero.Following the open loop simulation, the PMSM was simulated at superconducting temperatures, again with no control inputs.Figure 1 shows the characteristics of a typical loaded synchronous motor.Figure 1 shows the variation of the q-axis current i q as a function of time over the timeframe of 1 s at normal room temperature.Figure 1 also shows the corresponding variation of the d-axis current i d as a function of time, to the same scale.Figure 1 also shows the corresponding variation of the mechanical rotor speed ω m as a function of time over the same time frame.Figure 1 finally also shows the corresponding variation of the mechanical angular rotor position θ m as a function of time.
were set to zero.Following the open loop simulation, the PMSM was simulated at superconducting temperatures, again with no control inputs.Figure 1 shows the characteristics of a typical loaded synchronous motor.Figure 1 shows the variation of the q-axis current q i as a function of time over the timeframe of 1 s at normal room temperature.Figure 1 also shows the corresponding variation of the d-axis current d i as a function of time, to the same scale.Figure 1 also shows the corresponding variation of the mechanical rotor speed m ω as a function of time over the same time frame.Figure 1 finally also shows the corresponding variation of the mechanical angular rotor position m θ as a function of time.Figure 2 shows the variation of the dimensional load torque in Newton-meters as the motor starts and reaches steady state.Figure 3 shows the corresponding non-dimensional characteristics of the propeller, namely the coefficients of thrust and torque, the propulsive efficiency, and the advance ratio as functions of time over the timeframe of 1 s at normal room temperature.Figures 4-6 show the same characteristics as in Figures 1-3 but for the HTS motor.Figure 2 shows the variation of the dimensional load torque in Newton-meters as the motor starts and reaches steady state.Figure 3 shows the corresponding non-dimensional characteristics of the propeller, namely the coefficients of thrust and torque, the propulsive efficiency, and the advance ratio as functions of time over the timeframe of 1 s at normal room temperature.Figures 4-6 show the same characteristics as in Figures 1-3 but for the HTS motor.

Discussion and Conclusions
The above study demonstrates that HTS motors-that is, motors with armature and field windings made of HTS wires or tapes that are encased within a cryostat so that the temperature of the stator and rotor are in the superconductivity domain-may be modeled just as conventional motors.It must be recognized though that the values of the resistances, torque, and back EMF constants are substantially different from the corresponding values in a conventional motor.
The simulation of HTS PMSM was carried out in this paper with the aim of studying the nature of the control strategies that must be applied to a PMSM driving a propeller.The traditional method of controlling a PMSM is to use the strategy of field-oriented control (FOC).FOC is generally implemented by resolving the phase currents in the directions of the direct and quadrature axes, so the current in one axis generates the driving torque while the other generates the flux.One method of implementing such a control law is to use a voltage source inverter.A voltage source inverter (VSI) drive provides a voltage source as compared to a load commutated inverter (LCI) drive, which provides a current source.
The simulations carried out in this paper indicate that it is quite important to drive the propeller at an advance ratio such that it is always operating at the maximum propulsive efficiency point.Moreover, it is also important to operate the propeller while minimizing the overall energy losses.As precision positioning or speed controls are not really required, it was found that the use of a variable frequency drive and a direct current to alternating current power electronic converter may provide the optimum energy solution for driving the propeller.However, open loop variable frequency drive-based voltage source control may not be stable for all frequencies and it may be necessary to include a stabilizing inner loop.The stability and inner loop feedback stabilization of a PMSM with an outer open loop variable frequency drive and a constant load have been considered by many authors [25][26][27][28][29] and will not be reconsidered here.Yet, the methodology is particularly suitable for HTS PMSMs, as the stator resistance of these machines is particularly low.Thus, an MPEPT controller, where an outer loop identifies the operating speed at the maximum propulsive efficiency point while the variable frequency drive adjusts the drive frequency so as to force the motor to track the operating speed at the maximum propulsive efficiency point, is essential.For a PMSM, the parameters of the model of the PMSM and the desired advance ratio will influence the selection of the drive frequency.Thus, not only must the parameters of the PMSM be known accurately, it must be ensured that the HTS PMSM is always operating in identical conditions to ensure that the parameters are not changing continuously.Alternatively, one must determine the motor parameters in real-time and adapt the controller accordingly.
Currently, a generic voltage source controller based on a multi-loop control structure and an inner stabilization loop are being developed to control a typical HTS PMSM with a dynamic propeller load that is not constant, and the results of that study will be reported independently.The dynamic model is altered to include the load angle dynamics and the inner loop control law is synthesized by the method reported in [29].Scaling laws for matching the results of this study to other HTS PMSMs are used to facilitate the design of a distributed propulsion system for an all-electric aircraft.
a − e b , e bc = e b − e c , e ca = e c − e a .

Figure 1 .
Figure 1.Variation of the permanent magnet synchronous motor (PMSM) characteristics as a function of time over the timeframe of 1 s at normal room temperature.

Figure 1 .
Figure 1.Variation of the permanent magnet synchronous motor (PMSM) characteristics as a function of time over the timeframe of 1 s at normal room temperature.

Figure 2 .
Figure 2. Variation of the PMSM load torque as a function of time over the timeframe of 1 s at normal room temperature.

Figure 2 .
Figure 2. Variation of the PMSM load torque as a function of time over the timeframe of 1 s at normal room temperature.

Figure 2 .
Figure 2. Variation of the PMSM load torque as a function of time over the timeframe of 1 s at normal room temperature.

Figure 3 .
Figure 3. Variation of the PMSM-driven propeller non-dimensional coefficients of thrust and torque, propulsive efficiency, and advance ratio as a function of time over the timeframe of 1 s at normal room temperature.

Figure 3 .
Figure 3. Variation of the PMSM-driven propeller non-dimensional coefficients of thrust and torque, propulsive efficiency, and advance ratio as a function of time over the timeframe of 1 s at normal room temperature.Aerospace 2017, 4, x FOR PEER REVIEW 12 of 15

Figure 4 .
Figure 4. Variation of the high-temperature superconductor (HTS) PMSM characteristics as a function of time over the timeframe of 1 s at the superconducting temperature.

Figure 5 .
Figure 5. Variation of the HTS PMSM load torque as a function of time over the timeframe of 1 s at the superconducting temperature.

Figure 4 .
Figure 4. Variation of the high-temperature superconductor (HTS) PMSM characteristics as a function of time over the timeframe of 1 s at the superconducting temperature.

Figure 4 .
Figure 4. Variation of the high-temperature superconductor (HTS) PMSM characteristics as a function of time over the timeframe of 1 s at the superconducting temperature.

Figure 5 .
Figure 5. Variation of the HTS PMSM load torque as a function of time over the timeframe of 1 s at the superconducting temperature.

Figure 6 .
Figure 6.Variation of the HTS PMSM-driven propeller non-dimensional coefficients of thrust and torque, propulsive efficiency, and advance ratio as a function of time over the timeframe of 1 s at the superconducting temperature.

Figure 5 .
Figure 5. Variation of the HTS PMSM load torque as a function of time over the timeframe of 1 s at the superconducting temperature.

Figure 4 .
Figure 4. Variation of the high-temperature superconductor (HTS) PMSM characteristics as a function of time over the timeframe of 1 s at the superconducting temperature.

Figure 5 .
Figure 5. Variation of the HTS PMSM load torque as a function of time over the timeframe of 1 s at the superconducting temperature.

Figure 6 .
Figure 6.Variation of the HTS PMSM-driven propeller non-dimensional coefficients of thrust and torque, propulsive efficiency, and advance ratio as a function of time over the timeframe of 1 s at the superconducting temperature.

Figure 6 .
Figure 6.Variation of the HTS PMSM-driven propeller non-dimensional coefficients of thrust and torque, propulsive efficiency, and advance ratio as a function of time over the timeframe of 1 s at the superconducting temperature.

Table 1 .
Assumed parameters of a permanent magnet synchronous motor (PMSM)-propeller system.(The parameters corresponding to a superconducting temperature are indicated by 'sc').