Effect of magnetization process on levitation force between a superconducting. disk and a permanent magnet

Similar documents
Heat propagation and stability in a small high T superconductor. coil

Modeling and simulation on the magnetization in field-cooling and zero-field-cooling processes

MultiPhysics Analysis of Trapped Field in Multi-Layer YBCO Plates

ADHESION MEASURES OF ELASTO-PLASTIC THIN FILM VIA BUCKLE-DRIVEN DELAMINATION

physica status solidi current topics in solid state physics

RESEARCH ON RANDOM FOURIER WAVE-NUMBER SPECTRUM OF FLUCTUATING WIND SPEED

Effect of Rotation, Magnetic Field and Initial Stresses on Propagation of Plane Waves in Transversely Isotropic Dissipative Half Space

EXACT TRAVELLING WAVE SOLUTIONS FOR THE GENERALIZED KURAMOTO-SIVASHINSKY EQUATION

MATHEMATICAL AND NUMERICAL BASIS OF BINARY ALLOY SOLIDIFICATION MODELS WITH SUBSTITUTE THERMAL CAPACITY. PART II

Millennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion

Electromagnetic radiation of the travelling spin wave propagating in an antiferromagnetic plate. Exact solution.

arxiv:gr-qc/ v2 6 Feb 2004

Ph1c Analytic Quiz 2 Solution

DETERMINATION OF CARRIER TEMPERATURE FROM JUNCTION I(V) MEASUREMENTS

Non-Markovian study of the relativistic magnetic-dipole spontaneous emission process of hydrogen-like atoms

ELECTROMAGNETIC WAVES WITH NONLINEAR DISPERSION LAW. P. М. Меdnis

COMPUTER METHODS FOR THE DETERMINATION OF THE CRITICAL PARAMETERS OF POLLUTED INSULATORS

The homopolar generator: an analytical example

The simulation analysis of the bridge rectifier continuous operation in AC circuit

Where as discussed previously we interpret solutions to this partial differential equation in the weak sense: b

Part G-4: Sample Exams

"Research Note" ANALYSIS AND OPTIMIZATION OF A FISSION CHAMBER DETECTOR USING MCNP4C AND SRIM MONTE CARLO CODES *

Towards an Absolute Cosmic Distance Gauge by using Redshift Spectra from Light Fatigue.

Thermal Mechanisms of Stable Macroscopic Penetration of Applied Currents in High Temperature Superconductors and their Instability Conditions

Measuring & Inducing Neural Activity Using Extracellular Fields I: Inverse systems approach

ENERGY AND MOMENTUM IN ELECTROMAGNETIC WAVES

Control Theory association of mathematics and engineering

A model for measurement of the states in a coupled-dot qubit

Relativity in Classical Physics

19 Lecture 19: Cosmic Microwave Background Radiation

Developing Excel Macros for Solving Heat Diffusion Problems

Nuclear Shell Structure Evolution Theory

Dr G. I. Ogilvie Lent Term 2005

INTRO VIDEOS. LESSON 9.5: The Doppler Effect

Maximum Entropy and Exponential Families

Aharonov-Bohm effect. Dan Solomon.

Stress triaxiality to evaluate the effective distance in the volumetric approach in fracture mechanics

CAN SLOW ROTATING NEUTRON STAR BE RADIO PULSAR?

Chapter 9. The excitation process

The gravitational phenomena without the curved spacetime

Fiber Optic Cable Transmission Losses with Perturbation Effects

An Integrated Architecture of Adaptive Neural Network Control for Dynamic Systems

Four-dimensional equation of motion for viscous compressible substance with regard to the acceleration field, pressure field and dissipation field

arxiv:cond-mat/ v1 [cond-mat.str-el] 3 Aug 2006

Finite Element Analysis on Pulsed Laser Forming of Sheet Metal

Relativistic effects in earth-orbiting Doppler lidar return signals

arxiv:gr-qc/ v7 14 Dec 2003

STUDY OF INTERFACIAL BEHAVIOR OF CNT/POLYMER COMPOSITE BY CFE METHOD

University of Bristol - Explore Bristol Research

RULE OF MIXTURE FOR COMPOSITE THERMOELECTRICS

Possibility of Refuse Derived Fuel Fire Inception by Spontaneous Ignition

1 Josephson Effect. dx + f f 3 = 0 (1)

Fracture analysis of a functionally graded interfacial zone between two dissimilar homogeneous materials

Class XII - Physics Electromagnetic Waves Chapter-wise Problems

General Closed-form Analytical Expressions of Air-gap Inductances for Surfacemounted Permanent Magnet and Induction Machines

Calculation of Desorption Parameters for Mg/Si(111) System

A Spatiotemporal Approach to Passive Sound Source Localization

Physical Laws, Absolutes, Relative Absolutes and Relativistic Time Phenomena

materials ISSN

Review of classical thermodynamics

Final Review. A Puzzle... Special Relativity. Direction of the Force. Moving at the Speed of Light

Speed-feedback Direct-drive Control of a Low-speed Transverse Flux-type Motor with Large Number of Poles for Ship Propulsion

the following action R of T on T n+1 : for each θ T, R θ : T n+1 T n+1 is defined by stated, we assume that all the curves in this paper are defined

Determination of the Aerodynamic Characteristics of Flying Vehicles Using Method Large Eddy Simulation with Software ANSYS

A Time-Dependent Model For Predicting The Response Of A Horizontally Loaded Pile Embedded In A Layered Transversely Isotropic Saturated Soil

TEMPERATURE EFFECTS ON MECHANICAL PROPERTIES OF ZINC DITHIOPHOSPHATE TRIBOFILMS

ABSTRACT. and since the phase transition of asymmetric one-component fluid or mixture

Bäcklund Transformations: Some Old and New Perspectives

IMPACT MODELLING OF THE COEFFICIENT OF RESTITUTION OF POTATOES BASED ON THE KELVIN- VOIGHT PAIR

Nonreversibility of Multiple Unicast Networks

Simple Considerations on the Cosmological Redshift

1 sin 2 r = 1 n 2 sin 2 i

Natural Convection Experiment Measurements from a Vertical Surface

Nd-Fe-B permanent magnets for passive and hybrid magnetic bearings with extended application temperature range

Models for the simulation of electronic circuits with hysteretic inductors

MODELLING THE POSTPEAK STRESS DISPLACEMENT RELATIONSHIP OF CONCRETE IN UNIAXIAL COMPRESSION

non Sinusoidal Conditions

Modeling of Threading Dislocation Density Reduction in Heteroepitaxial Layers

What are the locations of excess energy in open channels?

Complexity of Regularization RBF Networks

Generalized Dimensional Analysis

CONVECTION AT A.MODEL ICE EDGE

Doppler effect of the rupture process of the great M W 7.9 Wenchuan earthquake

MODE I FATIGUE DELAMINATION GROWTH ONSET IN FIBRE REINFORCED COMPOSITES: EXPERIMENTAL AND NUMERICAL ANALYSIS

Is classical energy equation adequate for convective heat transfer in nanofluids? Citation Advances In Mechanical Engineering, 2010, v.

Accelerator Physics Particle Acceleration. G. A. Krafft Old Dominion University Jefferson Lab Lecture 4

A model of crack based on dislocations in smectic A liquid crystals. School of Physics, Beijing Institute of Technology

Directional Coupler. 4-port Network

Transient wave propagation analysis of a pantograph- catenary system

UTC. Engineering 329. Proportional Controller Design. Speed System. John Beverly. Green Team. John Beverly Keith Skiles John Barker.

The transition between quasi-static and fully dynamic for interfaces

The Concept of Mass as Interfering Photons, and the Originating Mechanism of Gravitation D.T. Froedge

ELECTROMAGNETIC NORMAL MODES AND DISPERSION FORCES.

Einstein s Three Mistakes in Special Relativity Revealed. Copyright Joseph A. Rybczyk

Three-dimensional morphological modelling in Delft3D-FLOW

Singular Event Detection

Stabilization of the Precision Positioning Stage Working in the Vacuum Environment by Using the Disturbance Observer

Calculation of Magnetic Field of Two Coaxially Located Circular Permanent Magnets Using a Method of Physical Analogies

Gravitomagnetic Effects in the Kerr-Newman Spacetime

Advances in Radio Science

Transcription:

Effet of magnetization proess on levitation fore between a superonduting disk and a permanent magnet L. Liu, Y. Hou, C.Y. He, Z.X. Gao Department of Physis, State Key Laboratory for Artifiial Mirostruture and Mesosopi Physis, Peking University, Beijing 1871, P. R. China L. Xiao, H.T. Ren, Y.L. iao, M.H. Zheng General Researh Institute for Non-ferrous Metals, Beijing 188, China Abstrat The levitation fores between a permanent magnet and a oaxial superonduting disk after different magnetization proesses are measured. Signifiant effet of the magnetization proess on the levitation fore is observed. Theoretial alulations of levitation fore based on the ritial state model with temperature and field-dependent ritial urrent density, and the heat dissipation due to the flux motion are in perfet agreement with the experimental data. Keywords: superonduting disk, levitation fore, field ooling, zero-field-ooling PACS number(s): 74.6.2w, 74.25.Ha, 74.25.Ld Introdution It is well known that the very high trapped magneti field in bulk high temperature superonduting disk (SD) and the levitation fore between SD and a permanent magnet (PM) have possible industrial appliations [1-7]. Several groups reported their theoretial and experimental results in reent years. Qin et al. reported the theoretial alulation of the hystereti levitation fore between

a PM and a non-magnetized SD after zero-field-ooling (ZFC) from first priniples [8]. Wang et al. onsidered the magnetization of the SD, and alulated the levitation fore between a PM and a SD after field-ooling (FC) [9]. In the above alulations, however, the ritial urrent density as a funtion of the magneti field and temperature, as well as the heat dissipation in SD was not taken into aount. Ohsaki et al. onsidered both fators and alulated the distributions of the urrent density and temperature in SD after a magneti field pulse, but their alulations did not distinguish the normal state region from the superonduting state region in the pulsed field magnetization proess [1]. Hou et al. proposed a systemati theoreti framework to disuss the basi equation of the urrent motion in SD by taking into aount the heat dissipation, the temperature and field-dependent ritial urrent density, and the different urrent dependent eletri fields E( j ) in the normal state region and the superonduting state region [11]. Aording to this framework they simulated the trapped magneti field in SD, whih is in good agreement with experimental results of several groups [12-13]. However, no experiment results on levitation fore are reported so far. In this paper, the levitation fores between a PM and a oaxial SD were measured after different magnetization proesses. The rate β at whih the applied magneti field is swept down in FC proess was hanged and signifiant effet on the levitation fore was observed. The alulated levitation fores aording to the systemati theoreti framework [11] oinide with the experimental data very well. Experimental The SD used in experiments is a melt-textured-grown (MTG) YBCO sample of 3 mm in diameter, 2 mm in thikness [14]. The permanent magnet (Nd-Fe-B) is the same size and the magneti field on its surfae is.53t. The SD sample is set in a ontainer first. In ZFC proess the SD is ooled by liquid nitrogen to 77 K below its ritial temperature (92K) without applied magneti field. In FC ase, a

hosen magneti field (.5T,.375T,.25T generated by the eletromagnets) is applied on the SD along the axial diretion of the disk at room temperature, then the SD is ooled below its ritial temperature in liquid nitrogen and the applied magneti field is swept down at a given rate β (.1T/s or.25t/s) to zero. After above proess, we wait for 1 minute and set the SD on the sample holder in the levitation fore measurement system, whih is shown in Fig. 1. With a jak the SD an approah the PM the levitation fores in this proess are measured by a fore sensor, and the distane between the SD and the PM is measured preisely by a displaement sensor. These eletrial signals are transmitted to the omputer and proessed by LABVIEW program. Theoretial Basi equations for urrent motion Aording to the systemati theoreti framework [11] we onsider the SD with radius a and thikness b. The desired equation of motion for the urrent density (, t) r in the superonduting disk as following [9]: ( r, ) 1 a b = µ 1 ρ' ' (, ') ( ) ( ', ') r r yl + φ ρ t d dzq E A z (1) where A ϕ is the vetor potential of applied magneti field, yl (, ') = ( ρ, ρ', ') Q rr f z z (2) and 1 Q is the reiproal kernel, whih is defined by a b ' ' 1 (, ') ( ', '') = δ ( '') dρ dzq rr Q r r r r (3) with ( ρ, ρ', ') f z z = π dφ ρ'osφ 2π 2 2 2 ( z z' ) + ρ + ρ' 2 ρρ'osφ 12

1 ρ ' 1 2 2 2 = 1 k K 2 ( k ) E( k ) πk ρ and k 2 = 4 ρρ ' ( ρ + ρ' ) + ( z z' ) 2 2 Eq. (1) an be easily time integrated by starting with ( ) zt ρ,, = ( is the initial urrent density distribution in the SD disk) and then by putting (,, = + ) = (,, ) + (,, ) soon as the indued urrent density (, zt, ) ρ zt t dt ρ zt ρ zt dt. As ρ is obtained, the vetor potential generated by the indued urrent density A an be derived as, (, ) a b ρ µ ρ' ' (, ') ( ') A z d dz Q = rr r. (4) And the radial and axial trapped magneti field an be written in the form of, B ρ A =, z B z ( ρ A ) 1 = ρ ρ (5) respetively. These equations should be supplemented by relationships between and the magneti field B and the eletri field E, whih depends on the material. In the superonduting state region, where <, is the ritial urrent density, the power-n model is used to desribe the nonlinear harateristis of the superondutor [9]: = E E n (6a) n = σ + 1 and σ is the flux reep exponent, and in the normal state region, where, we use E( ) = E (6b) Generally the ritial urrent density depends on both the loal field B and the temperature T. The Kim model is used to desribe the flux density dependene of [1]: = B B + B (7) where is when B = and B is a parameter. We inlude the temperature dependene of as the following equation [1]:

2 2 T α = 1 T (8) where T is the ritial temperature at B = T and α is onstant. Heat dissipation When a superondutor is subjeted to a non-stationary external magneti field, the heat generation rate per unit volume is W = E (9) The temperature hange due to the heat generation is desribed by the heat diffusion equation. T κ = t 2 C T W (1) here κ is the thermal ondutivity, and C is the heat apaity per unit volume. The levitation fore Aording to the similar disussion as before [9], the top surfae enter of the SD approahes the bottom surfae enter of PM as, s = z vt (11) where veloity v represents the speed at whih the PM approahes the SD, z is the initial distane between the top surfae of the SD and the bottom surfae of the PM. As the urrent density (, zt, ) ρ and the radial magneti field PM B ρ inside the SD have been derived, the vertial levitation fore along the z-axis an be readily obtained as [9], Seletion of Parameters z a b PM = 2 π ρ ρ ( ρ, ) ( ρ, ) (12) F d dz z B z b ρ The parameters used in our alulation are taken as follows: T = 92 K [9]; µ 4π 1 7 = ; 6 3 C =.88 1 /m K [8]; κ 1 1 = 6 W m K [8]; 4 E = 1 1 V m [1]; 9 2 α = 1.9 1 A / m ; v = 2 mm s ; B =.5Τ and the flux reep exponent σ = 4. All these parameters remain unhanged

in the following disussions. Results and Disussions The experimental results of the levitation fore measured in both the FC and ZFC ases are shown in Fig.2. It an be observed that the levitation fore in the FC ase ( β =.25 T/s) has a muh larger value than that in the ZFC ase. Taking the distane of 1 mm for example, the levitation fore is 45 N and 1 N for the FC and ZFC ases respetively. In the FC magnetization proess, after the applied magneti field is swept down to zero, large indued urrent remains in the SD, whih dominates the levitation fore, while the indued urrent in the SD by PM s approahing is muh smaller than that after FC [9]. But no indued urrent exists in the SD after the ZFC. There is only a small indued urrent in the SD when the SD approahes the PM. So it is reasonable that the high levitation fore is obtained after the FC. The solid urves in Fig.2 represent our alulated levitation fores using Equation (12). The alulation results are in perfet agreement with the experimental data. How the rate β affets the levitation fore is also shown in Fig.2. At lower rate the gradient of magneti field in the edge is smaller, driving fore is weaker, and less magneti flux esape from the SD. Thus the trapped magneti field at the enter of the SD an be higher after the applied magneti field is swept down to zero and the indued urrent density in the SD is larger, whih is orresponding to higher levitation fore. When the rate β keeps onstant.25 T/s, the magnitude B of the applied magneti field beomes a ruial fator. Fig.3 shows the levitation fore as a funtion of the distane s with different B. We measure larger levitation fores with a larger applied magneti field, whih is reasonable beause a larger urrent is indued in the SD with a larger applied magneti field and thus generates a larger levitation fore when the PM approahes at a onstant veloity v.

Conlusion We have set up a levitation fore measurement system and measured the levitation fores between a permanent magnet and a oaxial superonduting disk after different magnetization proesses. The effets of the magnitude of applied magneti field and the rate on the levitation fore are studied in both experiment and theoretial alulation. The alulated results oinide with the experimental data very well, whih strongly supports the systemati theoreti framework proposed by Hou et al.. Aknowledgments This work was supported by the National Siene Foundation of China (NSFC 11744), the Ministry of Siene and Tehnology of China (Projet No. NKBRSF-G19996462) and the TaiZhao Foundation of Peking University. Referenes [1] B.R. Weinberger, L. Lynds,.R. Hull, and U. Balahandran, Appl. Phys. Lett. 59, 1132 (1991). [2] M. Matsunaga, et al. Superond. Si. Tehnol. 15, 842 (22). [3] H. Minami, and. Yuyama, pn. Appl. Phys. 34, 346 (1995). [4] L.K. Kovalev, et al. Superond. Si. Tehnol. 15, 817 (22). [5] T. Oka, et al. Physia C 335, 11 (2). [6] S.I. Yoo, T. Higuhi, N. Sakai, H. Fujimoto and M. Murakami, Mater. Si. Eng. B 52, 23 (1998). [7]. Wang, et al. IEEE Trans. Appl. Superond. 9, 94 (1999). [8] M.. Qin, et al. Phys. Rev. B. 66, 24516(22). [9].. Wang, et al, Superond. Si. Tehnol. 16, 527 (23). [1] H. Ohsaki, T. Shimosaki and N. Nozawa, Superond. Si. Tehnol. 15, 754 (22). [11] Y. Hou, et al. http://arxiv.org/abs/ond-mat/311499.

[12] G. Fuhs, et al. Appl. Phys. Lett. 76, 217 (2). [13] S. Gruss, et al. Appl. Phys. Lett. 79, 3131 (21). [14] H.T. Ren, et al. Physia C,282-287,485 (1997) Figure aptions: Fig. 1: Experimental setup for the levitation fore measurement. Fig. 2: The measured levitation fores F as a funtion of distane s with different β. The solid urves z are alulated results. Fig. 3: The measured levitation fores F as a funtion of distane s with different magnitudes of the z applied magneti field. Fig. 1

Fig. 2

Fig. 3