Technical Progress in HTS Magnetic Bulk Application Development

Similar documents
Development of a cryogenic induction motor for use with a superconducting magnetic bearing

Physica C 468 (2008) Contents lists available at ScienceDirect. Physica C. journal homepage:

Effect of Pre-magnetization on Quasistatic Force Characteristics in a Space Superconducting Interface Structure Adopting High T c Superconductors

Effects of the field dependent J/sub c/ on the vertical levitation force between a superconductor and a magnet

Members Subjected to Torsional Loads

Mitigation of Demagnetization of Bulk Superconductors by Time-Varying External Magnetic Fields

Reduction of Mechanical Loss of Flywheel Energy Storage System with Spherical Spiral Groove Bearing

Maglev by a passive HT C superconducting short runner above the magnetic guideway

Concept Design and Performance Analysis of HTS Synchronous Motor for Ship Propulsion. Jin Zou, Di Hu, Mark Ainslie

Study on Trapped Field Characteristics of HTS Bulk Annuli With Iron Rings for Ferromagnetic Shimming of a Compact NMR Magnet

Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science Electric Machines

Study of the Pulsed Field Magnetization Strategy for the Superconducting Rotor

Comparison of Trapped Field Characteristic of Bulk Magnet System Using Various Type Refrigerators

Advanced Higher Physics. Rotational motion

AS mentioned in [1], the drift of a levitated/suspended body

Equal Pitch and Unequal Pitch:

Analysis of Hybrid Magnetic Bearing for High Speed. Spindle

World Record, High Magnetic Fields from Bulk Superconductors

Design and analysis of Axial Flux Permanent Magnet Generator for Direct-Driven Wind Turbines

Theory & Practice of Rotor Dynamics Prof. Rajiv Tiwari Department of Mechanical Engineering Indian Institute of Technology Guwahati

Physics. Student Materials Advanced Higher. Tutorial Problems Mechanics HIGHER STILL. Spring 2000

Conception of a YBCO Superconducting ZFC-Magnetic Bearing Virtual Prototype

Research Article Trial Application of Pulse-Field Magnetization to Magnetically Levitated Conveyor System

Measurement Techniques for Engineers. Motion and Vibration Measurement

Prototype ring spinning tester with superconducting magnetic bearing system for high productivity

Experimental Assessment of Unbalanced Magnetic Force according to Rotor Eccentricity in Permanent Magnet Machine

Chapter 14. Optical and Magnetic Materials. 경상대학교 Ceramic Design Lab.

SAMCEF For ROTORS. Chapter 1 : Physical Aspects of rotor dynamics. This document is the property of SAMTECH S.A. MEF A, Page 1

Stress Analysis Lecture 3 ME 276 Spring Dr./ Ahmed Mohamed Nagib Elmekawy

6) Motors and Encoders

Nonlinear effects on the rotor driven by a motor with limited power

Loss analysis of a 1 MW class HTS synchronous motor

New Electric Reluctance Motor with Bulk Superconducting Materials on the Rotor

1 MR SAMPLE EXAM 3 FALL 2013

The first man-loading high temperature superconducting Maglev test vehicle in the world

System design of 60K Stirling-type co-axial pulse tube coolers for HTS RF filters

Physics 212 Question Bank III 2010

Pulsed field magnetization of 0-0 and bridge-seeded Y-Ba-Cu-O bulk superconductors

Design and Experimental Evaluation of the Flywheel Energy Storage System

Doubly salient reluctance machine or, as it is also called, switched reluctance machine. [Pyrhönen et al 2008]

STATIC AND DYNAMIC ANALYSIS OF A PUMP IMPELLER WITH A BALANCING DEVICE PART I: STATIC ANALYSIS

Nonlinear Rolling Element Bearings in MADYN 2000 Version 4.3

DESIGN AND ANALYSIS OF LIGHT WEIGHT MOTOR VEHICLE FLYWHEEL M.LAVAKUMAR #1, R.PRASANNA SRINIVAS* 2

MAGNETIC PROBLEMS. (d) Sketch B as a function of d clearly showing the value for maximum value of B.

Superconductivity. The Discovery of Superconductivity. Basic Properties

Rotational Dynamics Smart Pulley

3D Finite Element Modeling and Vibration Analysis of Gas Turbine Structural Elements

Tests on Superconductor Gravitational Effects

Physics 12. Unit 5 Circular Motion and Gravitation Part 1

Physics 212 Question Bank III 2006

High-Performance Y-based Superconducting Wire and Their Applications

TOPIC D: ROTATION EXAMPLES SPRING 2018

LECTURE 12. STEADY-STATE RESPONSE DUE TO ROTATING IMBALANCE

A Preliminary Investigation of an Electrodynamic Wheel for Simultaneously Creating Levitation and Propulsion

3 d Calculate the product of the motor constant and the pole flux KΦ in this operating point. 2 e Calculate the torque.

Stepping Motors. Chapter 11 L E L F L D

Chapter 9 [ Edit ] Ladybugs on a Rotating Disk. v = ωr, where r is the distance between the object and the axis of rotation. Chapter 9. Part A.

Switching of magnetic domains reveals spatially inhomogeneous superconductivity

Development of a Compressor for a Miniature Pulse Tube Cryocooler of 2.5 W at 65 K

Postprint. This is the accepted version of a paper presented at ACTUATOR 2014, Bremen, Germany, June 2014.

Development of axial flux HTS induction motors

A Linear Motor Driver with HTS Levitation Techniques

EE 410/510: Electromechanical Systems Chapter 4

CHAPTER 3 INFLUENCE OF STATOR SLOT-SHAPE ON THE ENERGY CONSERVATION ASSOCIATED WITH THE SUBMERSIBLE INDUCTION MOTORS

SUPERCONDUCTOR JOULE LOSSES IN THE ZERO-FIELD- COOLED (ZFC) MAGLEV VEHICLE

Suggested Problems. Chapter 1

A turnout switch for the superconductively levitated transport system SupraTrans

Keywords: Electric Machines, Rotating Machinery, Stator faults, Fault tolerant control, Field Weakening, Anisotropy, Dual rotor, 3D modeling

Influence of the ramp angle on levitation characteristics of HTS maglev

Simultaneous measurement of critical current, stress, strain and lattice distortions in high temperature superconductors

PhysicsAndMathsTutor.com 1

Field homogeneity improvement of maglev NdFeB magnetic rails from joints

Name: Fall 2014 CLOSED BOOK

LAB REPORT: THREE-PHASE INDUCTION MACHINE

A Cryogenic Magnetostrictive Actuator using a Persistent High Temperature Superconducting Magnet, Part 1: Concept and Design

FACTORS AFFECTING CHARACTERIZATION OF' BULK HIGHTEMPERATURE SUPERCONDUCTORS*

A trapped field of 17.7 T in a stack of high temperature superconducting tape

Development of the Floating Centrifugal Pump by Use of Non Contact Magnetic Drive and Its Performance

IEEE Transactions on Applied Superconductivity. Copyright IEEE.

Third harmonic current injection into highly saturated multi-phase machines

Dynamic levitation performance of Gd Ba Cu O and Y Ba Cu O bulk superconductors under a varying external magnetic field

8.012 Physics I: Classical Mechanics Fall 2008

Title use of Bi-2223/Ag squirrel-cage rot IEEE TRANSACTIONS ON APPLIED SUPERCONDUCTIVITY (2006), 16(2): 14.

Chapter 7: Stepper Motors. (Revision 6.0, 27/10/2014)

ET-105(A) : PHYSICS. Show that only an infinitesimal rotation can be regarded as a vector.

Dynamics of Rotor Systems with Clearance and Weak Pedestals in Full Contact

Superconductor Joule Losses in the Zero-Field-Cooled Maglev Vehicle

High-temperature superconducting magnet for use in Saturated core FCL

Vibration characteristics of a multi-block high-temperature superconducting maglev system

Generators for wind power conversion

Jitter and Basic Requirements of the Reaction Wheel Assembly in the Attitude Control System

Lecture Module 5: Introduction to Attitude Stabilization and Control

Piedmont Chapter Vibration Institute Training Symposium 10 May, 2012 FIELD BALANCING OF ROTATING MACHINERY.

Circular motion minutes. 62 marks. theonlinephysicstutor.com. facebook.com/theonlinephysicstutor Page 1 of 22. Name: Class: Date: Time: Marks:

Simulation and Analysis of Hybrid Magnetic Bearing for High Speed Spindles

This equation of motion may be solved either by differential equation method or by graphical method as discussed below:

Cycloidal Magnetic Gear Speed Reducer

DESIGN AND APPLICATION

Instrumental technique. Manju C K

Strong High-Temperature Superconductor Trapped Field Magnets

Transcription:

p dia / µ H c 2 /2 ID 1953993 1 Technical Progress in HTS Magnetic Bulk Application Development F. N. Werfel, U. Floegel-Delor, T. Riedel, R. Rothfeld, P. Schirrmeister, D. Wippich, R. Koenig Abstract Experimental studies and theoretical modeling of High Temperature Superconductor (HTS) magnetic bulk applications are conducted and reported. The magnetic force density as a function of the magnet superconductor separation in the mirror image FEM 2D picture shows a strong dependency on the distance and PM mass. The behavior contrasts all attempts of generating high magnetic levitation forces with a minimum of material effort. An optimum configuration PM HTS coupler has been designed and built for a 3 mm wafer spin test apparatus. With rotational speeds up to 15 rpm a safe operation and handling the imbalance of the wafer spin clutch by passing the resonance speed at 37 rpm is investigated. The smart rotor rotation is studied to drive a half-wave plate of a polarization apparatus to detect cosmic microwave radiation safely through resonances and imbalance effects. Index Terms HTS bulk magnet application, FEM, magnetic wafer spin clutch, coupling phenomena, resonance effects, CMB polarization, and soft rotation. T I. INTRODUCTION he basic parameters of magnetic levitation are well undunderstand and are described in a number of excellent reviews and reports [1]-[4]. Technological skills as efficient cryogenics, minimal rotational loss or damping properties to utilize the stable levitation in high sensitive apparatus has hitherto been less carefully investigated. We investigate and report about the construction of HTS thrust bearings to support a 3 mm spin wafer magnetic clutch and rotate a half-wave plate of a polarimeter for measuring cosmic microwave background (CMB) radiation. A HTS magnetic clutch in semiconductor wafer processing is a new and unique technique. Equally effective is the soft continuous rotation of a half - wave plate at liquid Helium temperature for modulating the polarization of incident cosmic background radiation. In both cases a number of requirements and advantages favor the use of HTS bearing as low rotational friction, no wear and the absence of any motion control. In both applications reported here, the trapped magnetic flux capability in bulk HTS is basic for the function. Very recently, an extremely high trapped magnetic field of 17.6 T between a stack of two GdBCO samples has been successfully demonstrated by the Cambridge group [5]. Following their conclusions, the limitations of maximum fields are not determined by electric parameters of the HTS samples it seems rather a question of the mechanical stability Manuscript received August 12, 214 F. N. Werfel, U. Floegel - Delor, T. Riedel, R. Rothfeld, D. Wippich, P. Schirrmeister and R. Koenig are with the Adelwitz Technologiezentrum GmbH (ATZ), Adelwitz, Germany, Phone +49 34222 452, e-mail: werfel@t-online.de; of the material. Similar specific are the technical requirements of magnetic levitation to support and to spin a 12-inch silicon wafer in a hermetically closed semiconductor apparatus. The shaft-less rotor equipment is designed, built and tested to reduce particle generation during high sensitive wet and dry wafer spin processes. As a key feature the rotational control of the spinning wafer in its rotor dynamics and high precision is analyzed. Loading and de-loading the wafer uses an external robot system. In a first step a bearing of comparable size was modified to a magnetic coupling or clutch system. As a great advantage the magnetic clutch is capable to support and spin the wafer simultaneously and allows an exact control of the motion in all operational phases (processing, cleaning, drying). 1,1,1 II. MATERIALS A Optimum magnetic biasing From the viewpoint of force functionality a bulk superconductor basically acts as a passive element. An external magnetic field whether it is generated from permanent magnet (PM) or coil will interact with the HTS. Either the magnetic field is screened perfectly by surface currents (diamagnetic screening) or the magnetic flux is trapped by pinning forces within the bulk. Meissner effect is usually ignored by PM - HTS force consideration because of its weakness (a few tenths of a percent). Curiously, levitation experiments and pictures showing a flying PM above a cooled HTS block are still insistently described as Meissner effect. However, the Meissner part is very weak, it s close to Mirror principle FEM (2D) relative force density B,,5 1, 1,5 2, C D C B gap / diameter h = 2 x d = 1 x d =.5 x d g D d PM HTS mirror Fig. 1. Superconductor magnet interaction: 2D FEM mirror principle calculation for designing PM configuration

axial force F (mn) ID 1953993 2 one hundred percent pinning forces that levitates a PM several centimeter stably above a bulk HTS. For efficient application the magnitude of magnetic forces determines the necessary amount of HTS material and design. Thereby magnetic circuits obey the basic laws of electricity and magnetism, in detail the rules for conservation of magneto-motoric force. The PM property depends on the energy product (BH) max with an operating point dependent on the application. In most levitation cases the magnetic flux density being about half of the remanence force B PM ~ B r /2 determines the performance. It follows the basic experience that higher magnetic forces require larger magnetic flux density excitations equivalent to a larger volume and mass input of magnet material. To investigate this relation detailed either an analytical [6]- [7] or at least semi-quantitative picture is helpful. A first and useful step is the so called mirror image concept in figure 1. In this picture one assumes an ideal diamagnetic superconductor response to an external magnetic excitation. The magnetic interaction is thereby considered as a one-to-one reaction. In this model both the field magnitude of the response excitation as well as the position of the diamagnetic image follow perfectly the origin. The horizontal stability shows some itinerancy due to edge effects of the superconductor dimension. In figure 1 FEM calculations of permanent magnet interaction with a HTS in a principle 2D mirror image are shown. The gap distance is a simple ratio of the diameter of the PM giving three curve parameters for the obtained relative diamagnetic pressure. The resultant force density is referenced to the magnet performance/ filling factor by increasing the PM height as half, equal and double the gap. The calculated curves give a first approach of the expected levitation performance dependent on magnet size and distance of the magnet from the superconductor. B Design of HTS bulk The second part of the progress in HTS bearing and levitation is the superconductor. Figure 2 displays the concept and geometry of seeded REBCO melt textured growth. Typical bulk stator geometries displayed in the lower part of figure 2 represents a 3 mm YBCO ring (left) with orthogonal c axes and a 8 mm journal-type cylinder with radial c axes (right). Both components consist of assembled, glued and machined single grain bulks. Fig. 2. HTS bulk melt texturing technique (top); well oriented thrust type and journal- type YBCO stator components. A geometrical condition for generating high currents and high trapped field values is an optimum geometry of the magnetic field component which should be orthogonal to the crystallographic a, b plane of YBCO single grain bulk. Figure 3 displays measurements of different YBCO cylinders with a c axis orientation parallel and perpendicular to the radial magnetic flux distribution. Surprisingly, the axial force seems to be less sensitive to the crystal orientation of the YBCO stator for displacement values up to 4 mm. The polycrystalline stator exhibits lower magnetic forces as expected. While stable external levitation is a prominent performance of bulk superconductors the same forces act on the inner structure of the material too. Magnetic energy and pressure behave much like mechanical work. The force density exerted on a unit superconductor volume depend on the induced current density j and the magnetic induction B or the field inside the superconductor, The first term is the magnetic pressure whereby the second term corresponds to a tension of B/µ along the lines of force. 3 2 1 radial bearing, 4 x PM 2 x 1 x 5 mm³, 1 mm gap, Fe collector 1 mm, YBCO hollow cylinder 45x22x3 B --> versus bulk orientation HTS 1 2 3 4 Fig. 3. Measurement of the YBCO orientation relative to magnetic flux direction together with polycrystalline YBCO stator F displacement (mm) polycrystalline single grain B c B ab +Fe B c +Fe B ab III. HTS MAGNETIC WAFER CLUTCH (1) A Design Magnetic couplings and clutches are well known in modern machinery to transfer forces or momentum forces from one element to another one without mechanical connection. However, these types of coupling are unstable in at least one geometrical direction. Either the magnets attract each other or show repulsive forces. Therefore, magnetic bearings made of permanent magnets need an additional stabilization either in radial or in axial direction. In contrast, the interaction between a superconductor and magnet result in self-stabilization without additional effort. Based on the concept of trapped flux we designed and built a wafer clutch and used magnetic forces between permanent magnet and HTS bulk elements. Cooling of the HTS is performed using a developed rotatable vacuum cryostat with LN 2 cryogenics. To maintain a homogeneous

ID 1953993 3 moment of inertia of the cryostat the dynamics of LN 2 liquid in the cryostat during the rotation has been studied. The uncontrolled LN 2 distribution and the resulting momentum forces during acceleration and slowing-down the cryostat were especially critical for a save wafer spin operation. With the results of experiments and calculations we solved the challenge robust cryogenics by utilizing principles of ship construction. The cryostat was subdivided into a shell and chamber structure. The stator consists of a superconductor ring axially spaced to the rotor with a corresponding axial gap. Figure 4 shows the magnetic wafer clutch consisting of a PM rotor part and the ring-shaped YBCO in the rotatable cryostat. Both are operating magnetically through the bottom of a semiconductor processing chamber. The silicon wafer is fixed on a magnetic wafer carrier and spin without any mechanical contact in the sealed chamber. The processing chamber is fabricated in the shape of a Kloepper bottom to withstand vacuum or pressure forces at lowest wall thickness. The material of the chamber is high-grade non-magnetic CuAl - alloy known from submarine construction. B Function By consideration the force density dependence of Fig. 1 we use the mirror model PM configuration with height-to-lateral dimension of 1:1. The 12 cuboid-shape PM pairs share the outer magnetic clutch part that comprises of an YBCO annular ring. The ring is mounted of individual YBCO blocks machined to fit the ring shape. After field cooling (fc) the cryostat is driven by a motor on a shaft and the wafer turn table rotates together with the semiconductor wafer. With such a magnetic wafer spin apparatus it is possible to process a semiconductor at elevated temperatures up to 2 C, higher pressure and without possible induction of dangerous currents in electric circuits. According to figure 1 the transferable torque force can be increased by larger size, smaller distance to the HTS bulk, and higher pole number of the designed magnet circuit. For smaller twist angles of 5 8 at fixed stator the rotor has a tendency to switch back into original position. This retractable force is the limit for the applicable guidance force to accelerate or break the rotor. Even, if the driven rotor is removed from the process chamber the pinning forces in the superconductor allow an exact re-positioning of the rotor with reference to the driving rotor position. The technical data are collected in Table I. A rotor exchange with similar magnet circuit can be performed and repeated as long as the superconductor is under cryogenic conditions. C Wafer carrier in the critical rpm The critical speed of the wafer carrier occurs at about 37 rpm. With this value the calculated radial stiffness is k r = 6 N/mm. For a safe operation it is necessary to pass the critical rpm region in a short time to prevent large elongating amplitudes of the carrier. To obtain a first approximation how fast one should pass the critical rpm and to limit the magnitude of the amplitude we consider the rotor dynamic equation (2) Solving this equation delivers a homogeneous and a particular R res Im Fig. 4. 3 mm wafer spin rotor and unbalance effects Table I Wafer carrier parameters Excitation B 12 x PM PM HTS distance 12 mm Axial stiffness S ax 16 N/mm Radial stiffness S rad 7 N/mm Rotational stiffness (arc) 4 N/mm Moment of inertia.673 kgm²(carrier) t.42 kgm² (total) Rotational moment 2.1 Nm (carrier) Dynamic moment M 13.2 Nm Spin-up 15 rpm time 5 s Motor power x M 2 kw Maximum speed 15 rpm Re Fig. 5. Orbital radius vector under resonance imbalance excitation followed archimedic helix geometry solution for the imbalance excitation. According to the solutions the orbit radius R increases with the rotational speed, at the beginning rather fast followed by a gradual deceleration to the maximum value of the amplitude in the resonance. The orbital radius in the resonance is R res = /2 D with the eccentricity and the damping factor D. In the limit case of no damping (D = ) the center of the wafer ring as virtual axis follows an archimedic helix. The (virtual) axis point moves at the resonance frequency on an

magnetic flux variation [mt] m T ID 1953993 4 orbit of the shape of an Archimedean spiral. Geometrically the resultant curvature is described by a vector rotating with a constant angular velocity and a continuous increase of the vector value (length) as it is shown in figure 5. In the absence of any damping the magnitude of the orbital radius of the wafer carrier shown in figure 4 with an assumed eccentricity of.2 mm can be estimated within one second (1 s) to R res = ½.2 mm 2π (37/6s) 1 s = 3.9 mm (3) In this case the wafer ring would touch the chamber housing and needing a distinct damping behavior. The quantity U = m r is known as unbalance and the resulting centrifugal force is the product of unbalance and angular velocity squared. While unbalance forces (F centrifugal ) increase rapidly with speed, the unbalance quantity itself (m r) does not change at all. F centrifugal = m r ω 2 = U ω 2 (4) In fast- rotating machines the unbalance is a combination of static and dynamic unbalance outlined in figure 4. IV. POLARIZATION EXPERIMENT WITH HTS BEARING Polarization experiments of the Cosmic Microwave Background (CMB) uses HTS bearings because of advantages over mechanical bearings having stick-slip friction effect that induce vibrations. The vibrations deposit mechanical energy in the detectors and cause wires to vibrate thereby inducing excess noise. The best technique currently available for modulating the polarization of the incident radiation is by We have developed thrust magnetic bearings for use with polarization experiments that is based on magnetic levitation HWP above HTS rings in sizes 2 4 mm. The HWP is mounted on corresponding assembled PM ring as rotor. The apparatus were designed to fit both in balloon borne bolometric experiments as well as stationary in European Southern Observatory ESO operating in the Atacama Desert region of Chile. The airborne experiment is displayed in figure 6 with the Minnesota group of S. Hanany [8]. The stator parts of the magnetic thrust bearings consist of single grain YBCO bulks fitted in ring shape and supported by annular holders of G-1 or brass. Both materials are nonmagnetic. The magnetic stiffnesses of a 3 mm HTS thrust bearing operated at a height of 38 km in the airborne experiment of figure 6 were determined to 8 N/mm both in axial and radial direction at 4 mm gap distance. At a larger distance of 1 mm between stator and rotor we measured 43 N/mm axially but only 12 N/mm radially. Similarly to the wafer spin apparatus rotation specific parameters play a dominant role: solid state critical rpm, stiffness, damping and losses. The rotational drag or friction depends on the magnetic distance and the homogeneity of the circular magnetization of the assembled PM ring with the HWP integrated. In figure 7 the recorded magnetic flux variation is measured above the rim center of the assembled PM ring. Direct on the surface the maximum deviation is about 6 mt over one rotation. At 2 mm distance the inhomogeneity is reduced to 4 mt which corresponds to about 1% of the total magnetic flux. 1 15 45 PM segment alignment 9 12 6 3-1 18 21 33 height above PM mm 2 mm 1 24 27 3 zero trace line Fig. 7. Recorded inhomogeneity of the assembled magnet ring carrying a half- wave plate for cosmic microwave polarization; Hall probe measurement (insert). V. CONCLUSIONS Fig. 6. Helium balloon borne EBEX polarization experiment for detection of cosmic microwave radiation (U. Minnesota). rotation a half wave plate (HWP). The radiation field encodes new and valuable information about the earlier development phases of galaxies. When a HWP is placed in the path of the incident light and turned at a rotation frequency f, the polarization vector of the output radiation rotates at a frequency of 4f. The signal recorded by a Bolometer detector is modulated at a frequency of 4f and standard lock-in techniques to extract its magnitude. Large-diameter HTS magnetic bearing and magnetic clutch configurations have been investigated for optimum magnet biasing and smart rotational behavior. 2D FEM calculations give a first approach for design and construction. The levitation experiments with a new HTS magnetic wafer spin clutch for semiconductor processing exhibit critical issues. The combination of unbalance and rotor dynamics requires damping and sufficient fast acceleration. Polarization experiments of cosmic microwave background radiation utilize the soft rotation performance of HTS bearing as low noise, low vibration and ultra-low friction.

ID 1953993 5 REFERENCES [1] E. H. Brandt, Rigid levitation and suspension of high-temperature superconductors by magnets, Am. J. Phys. 58, 43-49, 199. [2] F. C. Moon, Superconducting levitation- Applications to Bearings and Magnetic Transportation J. Wiley and Sons, New York, 1994 [3] J. R. Hull, Superconducting bearings Supercond. Sci. Technol. 13 R1 - R15, 2. [4] F. N. Werfel, U. Floegel Delor, R. Rothfeld, T. Riedel, G. Goebel, D. Wippich, and P. Schirrmeister, Superconducting Bearings, Flywheels and Transportation Supercond. Sci. Technol. 25, 147, 212. [5] J. H. Durrell, A.R. Dennis, J. Jaroszynski, M. D. Ainslie, K.G.B. Palmer, Y.- H. Shi, A. M. Campbell, J. Hull, M. Strasik, E. E. Hellström and D. A. Cardwell, A trapped field of 17.6 T in melt processed bulk Gd-Ba-Cu-O reinforced with shrink-fit steel, Supercond. Sci.Technol. 27, 821, 21. [6] C. Navau, N. De-Valle, A. Sanchez, Macroscopic Modeling of Magnetization and Levitation of Hard Type-II Superconductors: The Critical-State, Model, IEEE Trans. Appl. Supercond. vol.32, no.1, 82123, 213. [7] G. T. Ma, Consideration on the Finite-Element Simulation of High Temperature Superconductors for Magnetic Levitation Purpose, IEEE Trans. Appl. Supercond. vol. 23 no 5, p. 36169, 213. [8] S. Hanany, T. Matsumura, B. Johnson, T. Jones, J.R. Hull, and K.B. Ma, A cosmic microwave background polarimeter using superconducting bearings, IEEE Trans. Appl. Supercond.13,2128-2133, 23.