improvements in modelling of hts - eventos fct/unl...disk hts, 2 – inductor coil diameter of hts...

21
Kurbatova Ekaterina National Research University “Moscow Power Engineering Institute” [email protected] O4-14

Upload: others

Post on 21-Jul-2020

7 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Improvements in modelling of HTS - Eventos FCT/UNL...disk HTS, 2 – inductor coil Diameter of HTS 30 mm Thickness of HTS 10 mm Inner diameter of coil 80 mm Outer diameter of coil

Kurbatova Ekaterina

National Research University

“Moscow Power Engineering Institute”

[email protected]

O4-14

Page 2: Improvements in modelling of HTS - Eventos FCT/UNL...disk HTS, 2 – inductor coil Diameter of HTS 30 mm Thickness of HTS 10 mm Inner diameter of coil 80 mm Outer diameter of coil

Introduction

Model of HTS properties with transport currents and magnetization

Distribution of electromagnetic field sources in HTS

Comparison of calculation results with experiments

2

Page 3: Improvements in modelling of HTS - Eventos FCT/UNL...disk HTS, 2 – inductor coil Diameter of HTS 30 mm Thickness of HTS 10 mm Inner diameter of coil 80 mm Outer diameter of coil

Calculations of HTS bearings with small gaps show, that there are difference between experiments and calculations using critical-state models.

The critical state models do not describe the partial Meissner effect at FC modes.

It is proposed to expand existing models of HTS properties with the addition type of electromagnetic field source.

3

Page 4: Improvements in modelling of HTS - Eventos FCT/UNL...disk HTS, 2 – inductor coil Diameter of HTS 30 mm Thickness of HTS 10 mm Inner diameter of coil 80 mm Outer diameter of coil

m

iM=i

For ideal superconductor:

B=B0=µ0(M+H)=const

M=B0 /µ0 -H

4

Page 5: Improvements in modelling of HTS - Eventos FCT/UNL...disk HTS, 2 – inductor coil Diameter of HTS 30 mm Thickness of HTS 10 mm Inner diameter of coil 80 mm Outer diameter of coil

( ) ( ) ( ) maxρ( , , ) 0.5ρ 1 th 1 / 1 / ( ) 1 / ( , ) /C C C JT T T H T J T = + − − − − H J H J H

Mathematical model of the transport current is approximation of the critical state model. It connects the nonlinear resistivity of HTS material with current density, magnetic field strength and temperature.

0

0.2

0.4

0.6

0.8

1

1.2

0 100 200 300 400

ρ/ρ 0

(pu

)

Current density (106 А/м2)

T=85 К

T=80 К

T=75 К

5

Page 6: Improvements in modelling of HTS - Eventos FCT/UNL...disk HTS, 2 – inductor coil Diameter of HTS 30 mm Thickness of HTS 10 mm Inner diameter of coil 80 mm Outer diameter of coil

,max( , ) ( ) [1 ( / ( )) ] ,C C CJ T J T H T = −H H

( , ) 0,CJ T =H cHH

cHH

( )( ) ( )( )2 2

0 ,max 0( ) 1 / , ( ) 1 / ,c c c c c cH T H T T J T J T T= − = −

Nonlinear dependence of the critical current density on the magnetic fieldis one of the main characteristic of HTS material. To obtain accuratecalculation results this dependence should be determined.

6

Page 7: Improvements in modelling of HTS - Eventos FCT/UNL...disk HTS, 2 – inductor coil Diameter of HTS 30 mm Thickness of HTS 10 mm Inner diameter of coil 80 mm Outer diameter of coil

B=B0 =const, dM/dH=-1, if H1<H<H2,

else M=Mc

Linear sections are limited by thecritical magnetization line Mc. In thiscase the magnetic field is higher,than superconducting cylinder cancompensate.

At the linear part of the model magneticinduction is constant. Superconductingcylinder saves the trapped magnetic field.

For blue line in the picture:

Magnetization model

7

Page 8: Improvements in modelling of HTS - Eventos FCT/UNL...disk HTS, 2 – inductor coil Diameter of HTS 30 mm Thickness of HTS 10 mm Inner diameter of coil 80 mm Outer diameter of coil

C,max CM( , ) ( ) [1 ( / ( )) ] ,CM T H M T H T = − H

( , ) 0,CM T H =

CM ( )H TH

CM ( )H TH

( )( ) ( )( )2 2

CM CM0 C,max C0( ) 1 / , ( ) 1 / ,c cH T H T T M T M T T= − = −

1 – α=1.5, β = 0.2

2 – α=1.5, β = 1.0

3 – α=1.5, β = 2.0

4 – α=1.5, β = 5.0

8

Page 9: Improvements in modelling of HTS - Eventos FCT/UNL...disk HTS, 2 – inductor coil Diameter of HTS 30 mm Thickness of HTS 10 mm Inner diameter of coil 80 mm Outer diameter of coil

Elementary cylinders are identical

Elementary cylinders with

different Mс and Нс

Smooth transition

9

Page 10: Improvements in modelling of HTS - Eventos FCT/UNL...disk HTS, 2 – inductor coil Diameter of HTS 30 mm Thickness of HTS 10 mm Inner diameter of coil 80 mm Outer diameter of coil

α

m

φ

θ

m

β

γ

Elementary cylinders with different Мс, Нс, φ. Anisotropy

10

Page 11: Improvements in modelling of HTS - Eventos FCT/UNL...disk HTS, 2 – inductor coil Diameter of HTS 30 mm Thickness of HTS 10 mm Inner diameter of coil 80 mm Outer diameter of coil

1

2 5

4

Model parameters may be defined from the experimental measurementsor chosen from the comparison of calculation with data from somesimple experiments with HTS sample (measurements of magnetic fieldnear the surface of superconductor, force interaction betweenpermanent magnet and superconductor).

11

Page 12: Improvements in modelling of HTS - Eventos FCT/UNL...disk HTS, 2 – inductor coil Diameter of HTS 30 mm Thickness of HTS 10 mm Inner diameter of coil 80 mm Outer diameter of coil

( ) ( ) ( ) ( )

( ) ( )

1 2

3

0 0

1 1 1 1

10

1;

4 2 4

1 1;

4

j

j j

lN N ljj ст

k j k j j k kj jS S

N

e k j kj S

t t dS t dS tr r

t t dSr

= = = =

=

= − +

=

j

n rA J n M A

The calculation of J is performed using

equations of a nonstationary magnetic field

B=∇xA

J = E/ρ = -dA/dt - ∇φ

Calculation M is made similarly to

the stationary magnetic field

B = ∇xA, B = µ0(M+H)

M = f(H)

The model was realized in home-made software EasyMag3D, which is based on the

spatial integral equations method.

12

Page 13: Improvements in modelling of HTS - Eventos FCT/UNL...disk HTS, 2 – inductor coil Diameter of HTS 30 mm Thickness of HTS 10 mm Inner diameter of coil 80 mm Outer diameter of coil

1 2

1 – disk HTS, 2 – inductor coil

Diameter of HTS 30 mm

Thickness of HTS 10 mm

Inner diameter of coil 80 mm

Outer diameter of coil 60 mm

Length of coil 40 mm

0

5000

10000

15000

20000

25000

0 0.5 1 1.5 2 2.5 3

МДС

t, c

1

2A

B

MMF

13

Page 14: Improvements in modelling of HTS - Eventos FCT/UNL...disk HTS, 2 – inductor coil Diameter of HTS 30 mm Thickness of HTS 10 mm Inner diameter of coil 80 mm Outer diameter of coil

Point A

Point B

Model for current densityModel for magnetization

14

Page 15: Improvements in modelling of HTS - Eventos FCT/UNL...disk HTS, 2 – inductor coil Diameter of HTS 30 mm Thickness of HTS 10 mm Inner diameter of coil 80 mm Outer diameter of coil

Magnetization Transport current

Point A

Point B

Combined model

15

Page 16: Improvements in modelling of HTS - Eventos FCT/UNL...disk HTS, 2 – inductor coil Diameter of HTS 30 mm Thickness of HTS 10 mm Inner diameter of coil 80 mm Outer diameter of coil

Point A

Point B

Model for magnetization Model for current density

16

Page 17: Improvements in modelling of HTS - Eventos FCT/UNL...disk HTS, 2 – inductor coil Diameter of HTS 30 mm Thickness of HTS 10 mm Inner diameter of coil 80 mm Outer diameter of coil

Magnetization Transport current

Point A

Point B

Combined model

17

Page 18: Improvements in modelling of HTS - Eventos FCT/UNL...disk HTS, 2 – inductor coil Diameter of HTS 30 mm Thickness of HTS 10 mm Inner diameter of coil 80 mm Outer diameter of coil

0

20

40

60

80

100

0 10 20 30 40 50

Forc

e,

N

Displacement, mm

18

Page 19: Improvements in modelling of HTS - Eventos FCT/UNL...disk HTS, 2 – inductor coil Diameter of HTS 30 mm Thickness of HTS 10 mm Inner diameter of coil 80 mm Outer diameter of coil

1

2

Ø40

Ø20

Ø31

Ø60

30

10

N

S

ρ

z 0

3

5

4

1

2

4

2

3

1

5

-25

-20

-15

-10

-5

0

0 1 2 3 4

F, H

z, мм

19

Page 20: Improvements in modelling of HTS - Eventos FCT/UNL...disk HTS, 2 – inductor coil Diameter of HTS 30 mm Thickness of HTS 10 mm Inner diameter of coil 80 mm Outer diameter of coil

It is proposed to expand the model of the critical state by additional sources of a magnetic field-related currents. These sources are modeled by the distribution of magnetization.

The performed calculations and experiments showed that the developed models can modeling the partial Meissner effect at FC mode and have better results in the analysis of magnetic HTS bearing under cyclic loading.

It is necessary to perform more computational and experimental studies to justify the effectiveness of such models.

20