Issn: Rmf@ciencias - Unam.mx: Revista Mexicana de Física

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 9

Machine Translated by Google

Revista Mexicana de Física

ISSN: 0035­001X

rmf@ciencias.unam.mx
Sociedad Mexicana de Física A.C.
México

Sosa, V.
Levitation forces between a finite rectangular superconductor and a spherical magnet Revista Mexicana de
Física, vol. 63, núm. 6, noviembre­diciembre, 2017, pp. 497­504 Sociedad Mexicana de Física A.C.

Distrito Federal, México

Disponible en: http://www.redalyc.org/articulo.oa?id=57053213001

como citar

Problema completo
Sistema de información científica
Más información sobre este artículo
Red de Revistas Científicas de América Latina, el Caribe, España y Portugal
Página de inicio de la revista en redalyc.org Proyecto académico sin fines de lucro, desarrollado bajo la iniciativa de acceso abierto
Machine Translated by Google

INVESTIGACIÓN Revista Mexicana de F´ısica 63 (2017) 497­504 NOVIEMBRE­DICIEMBRE 2017

Fuerzas de levitación entre un rectángulo finito.


superconductor y un imán esférico
V. Sosa
´ ´
Universidad Autonoma de Yucat an, Facultad de Ingenier ´ıa,
´ ´ ´
Av. Industrias no Contaminantes por Periferico Norte S/N, Apartado Postal 150 Cordemex, M erida, Yuc., M exico.
´
Permanent address: Cinvestav­IPN Unidad Merida, Departamento de F ´ısica Aplicada,
´ ´
Apartado Postal 73 Cordemex, 97310, Merida, Yuc., M exico. e­mail:
victor.sosa@cinvestav.mx
Recibido el 17 de octubre de 2016; aceptado el 25 de julio de 2017

Se elaboró un modelo para describir la interacción entre un imán permanente esférico y un superconductor cuadrado en estado mixto.
La propuesta se basa en ideas de dipolos de imagen como es el caso del modelo de dipolo Frozen que se ha utilizado en la literatura. Más precisamente,
se colocó un dipolo congelado en el centro de la muestra superconductora; este dipolo “nació” cuando se alcanzó el primer campo crítico en ese
posición. Esto condujo a la fuerza de atracción del imán­superconductor en el estado de campo penetrado. En la descripción del estado de Meissner, el
El campo desmagnetizador no se incluyó en el cálculo. Se analizan tanto las fuerzas verticales como las horizontales.

Palabras clave: Levitación superconductora; analítico; cálculo de fuerza.

PACOS: 74.20­z; 74.25.Ha

1. Introducción
Varias revisiones sobre las fuerzas de levitación superconductoras.
han sido publicados, ver por ejemplo el artículo de Hull [6]
La levitación superconductora es un fenómeno fascinante que o más recientemente, el excelente trabajo de Navau et al. [7].
ofrece el desarrollo potencial de volantes superconductores rápidos para Como se mencionó anteriormente, en el modelo de dipolos de la
aplicaciones de almacenamiento de energía [1]. para hacer estos imagen se supone que el superconductor tiene un espesor infinito. Bajo esto
dispositivos posibles, es necesario tener una descripción confiable Según la hipótesis, las imágenes se encuentran dentro del superconductor,
de las complejas interacciones del superconductor con un campo magnético independientemente de la distancia a la que se encuentre el imán. Pero, para un superconductor más delgado
axialmente simétrico externo. Avances relevantes han que la altura del imán, esta idea no es factible. A
Se ha hecho desde el modelo inicial de Hellman et al. [2] Según ese modelo, La descripción basada en las ideas del modelo dipolo congelado es
un dipolo puntual (el imán) induce una propuesto en este trabajo. Para justificar la propuesta,
dipolo de imagen dentro del superconductor, que se considera realizó algunas mediciones en un sistema real y comparó
para llenar el espacio debajo de un plano infinito, es decir, el superconductor ellos con los cálculos. A continuación, los datos utilizados para alimentar al
tiene un espesor infinito. De esta manera, el magnético Se presenta el modelo teórico.
El campo en la interfaz superconductor/vacío permanece paralelo a la
superficie de la muestra. Esta condición describe
2. Experimental
el efecto Meissner en el superconductor, es decir, la expulsión del campo.
Describir las interacciones entre el imán y
El superconductor en estado mixto representa una tarea importante. Se utilizó una muestra de YBa2Cu3O7 (YBCO) dopada con CeO2.

A continuación, se analizan brevemente dos métodos comunes: el congelado (0,5% en peso) Este superconductor fue preparado y caracterizado en el

modelo dipolo [3] y el modelo Bean [4]. En el primer modelo Laboratorie CRISMAT (Francia). El procedimiento

se induce un segundo dipolo puntual dentro del superconductor, de fabricación y texturizado se describe en detalle en [8]. El

orientado de tal manera que siempre atraiga el imán. Como El superconductor se presionó isostáticamente antes de texturizarlo en un

el imán desciende sobre el superconductor y se acerca Disco de 20 mm de diámetro. La muestra se extrajo del pellet puliéndolo

la superficie, el dipolo se acerca también desde el otro lado. hasta formar un prisma rectangular de ancho

Una vez que el imán llega al punto más cercano al superconductor, el dipolo 2R = 10 mm, longitud 2S = 10 mm y espesor L = 3,5 mm.

se "congela" y permanece en una posición fija. El superconductor de textura fundida tiene una temperatura crítica

a pesar de futuros desplazamientos del imán. Es necesario para de 88,4 K y una densidad de corriente crítica de 6,5 × 104 A/cm2 en
77 K y campo cero.
añadir más dipolos de puntos congelados para reproducir la histéresis
observada en el comportamiento de la fuerza vertical [5]. En La configuración experimental se muestra en la Fig. 1. La muestra
En el modelo Bean la idea es estudiar la distribución local de superconductora se colocó en un baño de nitrógeno líquido (campo cero).
tanto la densidad de corriente ­ya sea dependiente o independiente enfriamiento). Entonces, un imán esférico con su magnetización dirigida
del campo magnético interno y el campo magnético interno verticalmente, que tiene una masa de 0,712 gy un radio de
el superconductor sometido a una fuerza externa (generalmente uniforme) 2 mm, se bajó con la mano hacia el centro del superconductor siguiendo un
campo magnético. Las fuerzas totales se calculan sumando camino vertical. Por lo repulsivo
las contribuciones de toda la muestra. fuerza ejercida por el superconductor, fue necesario empujar
Machine Translated by Google

498 V.SOSA

TABLA I. Lista de imanes utilizados para magnetizar el superconductor. Nd significa un imán de NdFeB y Sm significa un imán de SmCo. Bmagnet es el
campo promedio aplicado al superconductor en un área de 1 cm × 1 cm y Bsc es el campo magnético atrapado en la muestra.

Imán Tamaño (cm) Bimán(T) Licenciatura(T)

Nd­disco 5 × 5 × 0,63 0,366 ± 0,011 0,290 ± 0,015

Sm­Cuadrado 2,5 × 2,5 × 0,63 0,475 ± 0,013 0,310 ± 0,008

Nd­disco 3 × 3 × 1,25 0,54 ± 0,015 0,320 ± 0,010

Sm­Rectángulo 2,5 × 1,25 × 0,63 0,640 ± 0,013 0,330 ± 0,007

Disco pequeño 2,5 × 2,5 × 0,63 0,66 ± 0,008 0,330 ± 0,005

Nd­disco 2,5 × 2,5 × 0,63 0,75 ± 0,019 0,330 ± 0,007

Disco pequeño 2,5 × 2,5 × 0,95 0,80 ± 0,011 0,340 ± 0,014

Disco pequeño 1,25 × 1,25 × 0,63 0,85 ± 0,013 0,340 ± 0,008

Nd­rectángulo 5 × 2,5 × 1 0,88 ± 0,025 0,350 ± 0,011

firmemente la esfera. El imán fue liberado después de lograr un 3.1. Expulsión de campo
levitación estable, una vez que las líneas magnéticas fueron fijadas al
superconductor. Se midió una altura de a = 3,5 mm en Primero, la fuerza repulsiva entre el imán y un objeto perfecto.

esa posición. diamagnet (es decir, un superconductor en el estado de Meissner)


ser calculado. Para realizar este cálculo es necesario
El imán se caracterizó previamente midiendo
conocer la magnetización de la muestra superconductora,
el campo magnético en diferentes puntos a lo largo del eje de simetría,
M(Ha), donde Ha es el campo aplicado. En el caso de un
y se concluyó que tiene un momento dipolar magnético
expulsión completa del campo del superconductor, tenemos
µ = 2,533 × 10−2 A m2 . Se caracterizó la capacidad de la
B = 0, y por lo tanto M = −H = −H a − SM. Aquí, SM.
muestra superconductora para atrapar campos magnéticos
es el campo magnético producido por el superconductor con
de la siguiente manera. Se acercó a un imán permanente.
magnetización M, es decir, el campo desmagnetizante. Se obtiene la
a la superficie del superconductor sumergido en nitrógeno líquido hasta
siguiente ecuación integral:
que sus superficies hicieron contacto; entonces, el imán estaba
recuperado. Para estimar el campo atrapado en el superconductor,
se midió el campo en el centro de su superficie superior (Bsc) M (r) = −H a(r) − dHM. (1)
con sonda Hall. Esto se realizó con diferentes imanes. Carolina del Sur

más amplio que la muestra de YBCO. En todos los casos, el campo


En principio, esta ecuación podría resolverse mediante una iterativa.
producido por el imán (Bmagnet) se midió en su superficie
proceso. Generalmente, la magnetización del superconductor es
y escaneado sobre un área de 1 × 1 cm (el área del superconductor).
logrado mediante cálculos numéricos, principalmente utilizando el
Para tener un campo representativo, los 30
método de elementos finitos [11­13].
Se promediaron los valores medidos de Bmagnet . Bsc se midió en
La intención de este trabajo fue obtener análisis completos.
ejecuciones sucesivas al menos tres veces para cada imán.
expresiones de las fuerzas, pero para calcular la desmagnetización
y los valores también se promediaron. Los resultados se muestran en campo bajo un campo externo no uniforme H a es muy compli­
Tabla I.
tarea indicada. En una primera aproximación, el superconductor en
A medida que aumenta Bmagnet , Bsc aumenta hasta un valor de El estado de Meissner tiene una magnetización.
0,35 T. Está claro que la muestra podría atrapar campos más altos si
Se aplicaron campos más grandes. La capacidad de los superconductores M (r) = −H a(r). (2)
de alta Tc para atrapar campos magnéticos se ha estudiado
exhaustivamente en otros lugares [9,10]. Ahora bien, la intención de incluir Se demostrará que ya, en este nivel, el modelo describe
esta caracterización en este estudio fue comparar los valores medidos adecuadamente las características generales del sistema.
del campo atrapado con una estimación obtenida Es interesante recordar que según el modelo Bean
con el modelo desarrollado en la siguiente sección. del estado mixto, la magnetización inicial de una losa superconductora
infinita de densidad de corriente crítica Jc y espesor L, cuando se le
aplica un campo magnético externo, puede
3. Teoría escribirse como M = −Ha + H2 a/Hp, donde Hp = JcL es el
el llamado campo de penetración. Esta expresión es válida siempre que
La fuerza entre un imán y un superconductor en el Ha no supera Hp; Se supone que Ha es uniforme y paralelo a uno de los
El estado mixto se describirá utilizando dos modelos: magnético lados del superconductor. Se puede ver
expulsión total del campo (estado de Meissner), art. 3.1, y penetración que el supuesto de descuidar HM en este trabajo conduce a una
del campo magnético, Sec. 3.2. conclusión similar cuando se toma un valor infinito de Jc en la

Rev. Méx. Fis. 63 (2017) 497­504


Machine Translated by Google

FUERZAS DE LEVITACIÓN ENTRE UN SUPERCONDUCTOR RECTANGULAR FINITO Y UN IMÁN ESFÉRICO 499

Modelo Bean, en el sentido de que se obtiene la misma M(H a) en Como se muestra en la Fig. 1, el imán esférico con momento
ambos casos. dipolar µ está ubicado sobre un superconductor rectangular.
La magnetización de un superconductor de forma ortorrómbica El origen del sistema de coordenadas está en el centro del imán. La
en estado de Meissner ha sido previamente estudiada teórica y distancia entre este punto y la superficie superior de la losa es a. El
experimentalmente [14]. En ese trabajo, se sometieron pequeñas centro de la losa está ubicado en (0, 0, −(a + L/2)). El superconductor
muestras de niobio (típicamente 1,5 mm × 1,5 mm × 1 a 6 mm) a un se encuentra horizontal entre los planos z = −a y z = −(a + L) con
campo externo uniforme Ha de 0 a 7 × 104 A/m. sus lados paralelos a los ejes X e Y.
Para valores bajos de Ha, se concluyó que M = χHa (χ es la
susceptibilidad inicial). Para un superconductor con proporciones Es bien sabido que el campo magnético debido a un imán
geométricas similares a las utilizadas en este trabajo (ancho/largo esférico es el mismo que el producido por un dipolo puntual con el
0,35), su resultado es χ −2,5. A medida que crece la relación ancho/ mismo momento magnético ubicado en el centro de la esfera. Este
largo del superconductor, χ tiende a −1, es decir, el campo campo es:
desmagnetizador tiende a desaparecer. 1 3r(r µ) −
La fuerza sobre el superconductor se puede calcular sumando H a(r) = (4)
4p 5
rublos metro 3 r

todas las interacciones infinitesimales dF entre cada dipolo dμ = M


dV y el campo aplicado B a: Sustituyendo la ecuación. (4) en la ecuación. (3) y cambiando el signo
se obtiene la fuerza repulsiva eyectada sobre el imán. A continuación,
F= (dµ ∙ B a) = µ0 (dV M ∙ H a) se calcula la componente vertical de esta fuerza. Se analizaron las
Carolina del Sur Carolina del Sur
configuraciones vertical (µ = µˆk) y horizontal (µ = µˆi).

2
= −μ0 dV (Ha) (3)
Carolina del Sur
3.1.1. Configuración vertical

En este caso, la integración se realiza en todo el volumen del Cuando µ = µˆk, la fuerza vertical sobre el imán está dada por
superconductor. Según la ecuación. (3), se puede concluir que la 2

fuerza sobre el superconductor proviene de una densidad de energía µ0 metros

Fz = [g(a) − g(a + L)] 4π 16π (5)


potencial (dU /dV ) = µ0H2 . Este .
a
modelo dipolo­dipolo se ha utilizado previamente para calcular donde se define g en función de la variable independiente ξ
fuerzas de levitación en diferentes sistemas imán­superconductor, por como

ejemplo, cilindro dipolo [15], cilindro de barra [16], anillo dipolo [17], RS
cilindro cuadrado [18], esfera de alambre [19] y cilindro de anillo [20]. gramo(ξ)
= ξ 2(ξ 2 + R2) 2(ξ 2 + S2) 2(ξ 2 + R2 + S2) 2

6 ×R 2 2 4 4 4
(9x +6S ) +R (35x + 41ξ 2S 2 + 12S )

2 6 6
+R (48x + 80x 4S 2 + 41x 2S 4 +6S )

8
+ 22x + 48x 6S 2 + 35x 4S 4 + 9x 2S 6

4 S
3R(7ξ + 10x 2R2 + 4R4 )
+ arctán
ξ 4(ξ 2 + R2) 5/2 X 2 + R2
4 4 R
3S(7ξ + 10ξ 2S 2 +4S
+ ) arctán (6)
ξ 4(ξ 2 + S2) 5/2 X2 + S2

Tomando S → ∞, se obtiene la expresión para una losa


infinitamente larga:

3πR 7ξ 4
+ 10x 2R2 + 4R4
g(ξ) = (7)
2 ξ 4(ξ 2 + R2) 5/2

Además, cuando R → ∞, el resultado para una hoja infinita se


obtiene:

FIGURA 1. Configuración experimental del sistema. Se coloca un g(ξ) = ξ (8)
4
imán permanente esférico encima de un superconductor rectangular,
que se sumergió en un baño de nitrógeno líquido (no mostrado). Los Finalmente, si L → ∞, la fuerza vertical entre el imán
centros de la esfera y del superconductor permanecen en la misma y un plano superconductor semiinfinito resulta ser
línea vertical. El centro del imán está ubicado a una distancia a por
2
encima de la superficie superior del rectángulo. La flecha vacía indica µ µ03
Fz = (9)
el momento dipolar magnético µ de la esfera. 8a 4 4π

Rev. Méx. Fis. 63 (2017) 497­504


Machine Translated by Google

500 V.SOSA

la fuerza vertical es prácticamente la misma que la obtenida para una


muestra infinita. Esto explica el éxito del modelo de Hellman incluso
para estas muestras pequeñas.

3.1.2. Configuración horizontal

El resultado de Fz para µ orientado a lo largo de la dirección x en la


misma expresión que la ecuación. (5), pero en esta configuración la
función g(ξ) viene dada por

RS
gramo(ξ)
= ξ 2(ξ 2 + R2) 2(ξ 2 + S2)( ξ2 + R2 + S2) 2

×3R 6 + 6R 4 (2x 2 +S 2 ) +x 4 (X 2 +S 2 )

+R 2 (10x 4 + 10ξ 2S 2 +3S 4 )

FIGURA 2. Fuerza vertical (ecuaciones (5) y (6)) sobre un imán


3R(2ξ 4 + 5ξ 2R2 + 2R4 ) S
esférico con µ = 2,533 × 10−2 Am2 en función del ancho del + arctán
superconductor rectangular en estado de Meissner. Tomamos R=S ξ 4(ξ 2 + R2) 5/2 X2 + R2
por simplicidad. La esfera se coloca a = 3,5 mm por encima del
3S(3ξ ) 2arctán
+ 2S 2 R
superconductor, con su magnetización a lo largo de la dirección z. La + (10)
fuerza crece monótonamente y se satura en diferentes niveles, ξ 4(ξ 2 + S2) 3/2 X 2 + S2
dependiendo del valor de L. Para un espesor de 10 mm y superior, la
fuerza de saturación prácticamente coincide con el resultado dado Para S → ∞, es decir, una losa infinitamente larga:
por la ecuación (9) para una muestra semiinfinita (horizontal discontinua). línea).
3πR(2ξ 4 + 5ξ 2R2 + 2R4 )
g(ξ) = (11)
2x 4(x 2 + R2) 5/2

Ahora, tomando R → ∞, la expresión para una hoja infinita es la


siguiente:

g(ξ) = ξ (12)
4

Nuevamente, en el límite de un plano semiinfinito (L → ∞) ),


la fuerza vertical tiende al valor (µ0/4π)(3µ 2/16a que coincide con 4
el valor obtenido con el método de imágenes [2].

3.2. Penetración de campo

Para describir la penetración del campo magnético en el


superconductor se parte de una suposición basada en el modelo de
dipolo congelado. Un dipolo de punto magnético representa el
comportamiento de la región penetrada de la muestra. Esto equivale
a realizar una expansión multipolar del campo magnético y retener
FIGURA 3. Fuerza vertical (ecuaciones (5) y (10)) sobre un imán
sólo el término dipolar. En el presente modelo se asumirá que el
esférico con µ = 2,533 × 10−2 Am2 en función del espesor del
dipolo aparece y permanece en el centro de la muestra. Más
superconductor en estado de Meissner. Nuevamente, R = S, a = 3,5
precisamente, el dipolo aparece en este punto cuando se alcanza el
mm y µ en la dirección vertical. La fuerza crece monótonamente con
primer campo crítico Hc . Después de eso, a medida que el imán se
L y alcanza un valor de saturación que depende de R. La línea
discontinua horizontal representa el valor dado por el modelo de va acercando al superconductor, la magnitud del momento dipolar
Hellman para un superconductor semiinfinito. inducido µs aumenta hasta que el imán alcanza la posición donde se
alcanza el equilibrio mecánico. Se supone que la magnetización del
que es la expresión temprana obtenida por Hellman et al. [2] con el imán permanece orientada verticalmente en todo su recorrido. Por
método de imágenes. tanto, el dipolo del punto magnético inducido estará orientado siempre
La variación de Fz con las dimensiones del superconductor se en la misma dirección.
muestra en las Figs. 2 y 3. Se puede observar que esta fuerza crece
de forma monótona y alcanza un valor de saturación. Es interesante El primer campo magnético crítico se puede obtener de ]. Para
2
observar que para una muestra con 2R = 2S = 2 cm (aproximadamente Hc(T) = H0[1−(T /Tc) en una muestra de YBCO sumergida
seis veces la distancia al imán) y L = 1 cm, nitrógeno líquido se tomará H0 = 100 Oe, Tc = 90 K

Rev. Méx. Fis. 63 (2017) 497­504


Machine Translated by Google

FUERZAS DE LEVITACIÓN ENTRE UN SUPERCONDUCTOR RECTANGULAR FINITO Y UN IMÁN ESFÉRICO 501

debe ser igual al peso del imán más la fuerza entre el imán y µs.
De esta forma se obtiene el valor µs = 6,95 × 10−3 A m2 .

Para evaluar qué tan razonable es este valor, se calcularon


dos campos en el centro de la superficie superior del
superconductor: el debido al imán esférico que levita en su
posición de equilibrio y el producido por el dipolo inducido µs .
Estos campos son 0,12 T y 0,26 T respectivamente. Estos dos
números se compararon con los valores medidos del campo
atrapado en el superconductor (Tabla I), como se muestra en la
Fig. 4. Es fácil ver que la estimación teórica encaja muy bien con
la tendencia experimental.
Con base en las observaciones experimentales, se supone
que la magnitud µs varía de 0 a 6,95 × 10−3 A m2 a medida que
el imán desciende hasta el superconductor. Se supone que la
dependencia precisa de µs con la altura a dada en mm (3,5 < a <
FIGURA 4. Campo magnético atrapado en el superconductor
10,6) tiene un comportamiento exponencial:
rectangular en función del campo magnético de los imanes permanentes
descritos en la Tabla I. Los datos con barras de error corresponden a 10.6 − un norte

los valores experimentales allí reportados. El punto señalado con una µs(a) = 6,95 × 10−3Am2 (13)
7.1
flecha corresponde a los campos calculados producidos por el imán
esférico y el dipolo inducido µs en la superficie superior del donde n es un parámetro de ajuste. Siguiendo la idea del modelo
superconductor. Este punto teórico concuerda bien con el de dipolo congelado, se supone que cuando el imán se separa de
su posición de equilibrio, µs permanece congelado en el centro
comportamiento experimental observado. La línea de puntos es una guía para
ojo. del superconductor. Este comportamiento se esboza en la Fig. 5.

4. Fuerza total
En esta sección se analizan las contribuciones combinadas de las
fuerzas de atracción y repulsión sobre el imán. La fuerza eyectada
por un dipolo magnético µ situado en el origen sobre un segundo
dipolo µs situado en el punto r viene dada en general por: [21]

3μ0
Fs = (µ ∙ r) µs + ( µs ∙ r)µ +
4pr5

5(µ ∙ r)( µs ∙ r) +
(µ ∙ µs)r − r (14)
2r

En la configuración vertical de nuestra geometría, µ=µˆk,


FIGURA 5. Diagrama que ilustra el mecanismo propuesto de cómo la
µs=µs ˆk, y r = (x, 0, z). Sustituyendo en la ecuación. (14) y
muestra superconductora es penetrada por el campo magnético del cambiando de signo, la fuerza de atracción eyectada sobre el
imán es:
imán. (a) No se produce penetración mientras a>10,6 mm (Ha < Hc).
(b) Una vez que se alcanza el campo magnético crítico, aparece un
3 µ0 µs 2 5z
dipolo magnético µs que permanece en el centro del superconductor; Fm = 2 2 + zx
− 1 xˆi
2
µs aumenta cuando a disminuye, siguiendo la ecuación. (13); 3,5 mm 4π(x + de 2) 5/2
< a < 10,6 mm. (c) µs alcanza su valor máximo cuando el imán alcanza
la posición donde se alcanza el equilibrio mecánico. (d) 5z 2
+ 2 2 + zx
− 3 de ˆk (15)
El dipolo inducido permanece congelado (mismo valor y misma
posición) cuando se recupera el imán.
4.1. Vertical
y T = 77 K, y se obtiene Hc(77K) = 27 Oe. Este valor se alcanza
en el centro del superconductor cuando el centro del imán se La fuerza vertical total proyectada por el superconductor sobre el
encuentra a una distancia de 12,3 mm, es decir, cuando a = 10,6 imán es la suma de la calculada a partir de las ecuaciones. (5)­(6)
mm. y la componente z de la ecuación. (15). En particular, el análisis
Cuando se alcanza la levitación estable de la esfera en a = se restringirá al caso x = 0. Como se mencionó anteriormente, µs
3,5 mm, la fuerza repulsiva dada por las ecuaciones. (5) y (6) aumenta desde 0 hasta su valor máximo según la ecuación. (13)

Rev. Méx. Fis. 63 (2017) 497­504


Machine Translated by Google

502 V.SOSA

4.2. Horizontal

Cuando x = 0, no actúa ninguna fuerza horizontal sobre el imán.


Sin embargo, cuando se desplaza lateralmente a lo largo del eje
X, la fuerza horizontal resultante es la combinación del
componente x en la ecuación. (15) y la fuerza repulsiva FP debida
al estado de Meissner del superconductor. Este último se calculó
C
a partir de la ecuación. (3) usando Mathematica, y el resultado es

2 µ0 µ
PF = [P(R, S, −a, x) − P(R, S, −(a + L), x) 128π 2

+ P(−R, S, −a, x) − P(−R, S, −(a + L), x)] (16)

dónde

2 2 2
FIGURA 6. Fuerza vertical total sobre el imán en función de la Con 3S 2 2 z + [z + (α − η) + ](a­h)
distancia a, calculada con µ = 2,533 × 10−2 A m2 = 5 mm, , R = S P(a, S, z, n) = 2 2
T3 [Con (α − η) 2] 2(α − η)
L = 3,5 mm y n = 0,8. No se incluyó el peso del imán (= 7×10−3 N).
Las flechas indican si el imán se acerca (rama superior) o se aleja 2 2 3z 2
2z +S 2+ + (α − η) +
(rama inferior) del superconductor. Una curva FC histerética típica − S2+(α−η) 2
2
de la Ref. 6 se muestra en el recuadro y se observa una buena [S2 +z (α − η) 2] 2
concordancia cualitativa con el modelo. 2 2
3S[2S + 3(a­h) ]
+ 4
[S2 + (α − η) 2] 3/2(α − η)
Con

× arctán
2
S2 + (a −h)

2a 4 +2z 2
− 8α 3η + 5z 2η 2 + 2η
4

+
(α − η) 4[z 2 + (α − η) 2] 5/2

2 2 2 2 2
un
(5z + 12 ) − 2αn(5z + 4to )
+
(α − η) 4[z 2 + (α − η) 2] 5/2

S
x 3z arctán (17)
2
2z + (a­h)

Para una muestra infinita (R, S, L → ∞), la fuerza repulsiva


tiende a cero para cualquier valor finito de x, como debería ser.
La fuerza horizontal total Fx se muestra en la Fig. 8. Se puede
ver que para desplazamientos laterales inferiores a 4,4 mm, la
FIGURA 7. Fuerza vertical total sobre el imán en función de la fuerza de atracción domina sobre la repulsiva y se tira del imán
distancia a, calculada con µ = 2,533×10−2 A m2 L , R = S = 5 mm, para regresar a su posición de equilibrio. Por encima de 4,4 mm,
= 3,5 mm y n = 5. Nuevamente, aquí no se consideró el peso del la fuerza de repulsión es mayor que la fuerza de atracción y el
imán. Las flechas indican que el imán se acerca y se aleja del imán es expulsado. Este resultado es particularmente importante
superconductor (rama superior e inferior, respectivamente). Se para describir la estabilidad del imán alrededor de su posición de
observa una curva ZFC histerética típica. levitación. El cálculo está cualitativamente de acuerdo con los
resultados experimentales [23,24]. Para xz (| x |≤ 1 mm), la fuerza
durante el descenso del imán, y permanece con el mismo valor se acerca a la línea recta discontinua que se muestra en la Fig. 8.
durante el ascenso. Como Wu et al. [5], utilizamos el valor n = 0,8. La pendiente de esta recta es k = 50 N/m. Este valor es
Como se muestra en la Fig. 6, se observa una curva histerética notablemente similar al obtenido por Johansen et al [22] en su
típica para un proceso de enfriamiento de campo (FC). Los bucles estudio teórico de las oscilaciones laterales de un superconductor
menores en puntos intermedios del camino no fueron considerados rectangular en estado mixto que levita sobre un imán rectangular
en el modelo. Se ha demostrado que estos bucles pueden (k = 56 N/m para una separación de 3,5 mm). ). Si se desprecia la
explicarse añadiendo más dipolos congelados en el superconductor [5].amortiguación causada por la fricción, se espera que el imán de
Por otro lado, si n = 5, obtenemos un comportamiento típico masa 0,712 g oscile armónicamente en esta región con una
para un proceso de enfriamiento de campo cero (ZFC), como se muestra en la Fig. 7. frecuencia de 42 Hz. Para verificar esto, se le dio al imán una ligera

Rev. Méx. Fis. 63 (2017) 497­504


Machine Translated by Google

FUERZAS DE LEVITACIÓN ENTRE UN SUPERCONDUCTOR RECTANGULAR FINITO Y UN IMÁN ESFÉRICO 503

La frecuencia cal cambió a 39 Hz, todavía lejos de la medida.


valor. Es evidente que para mejorar el acuerdo es necesario
tener en cuenta la amortiguación del movimiento y calcular la
campo desmagnetizante del superconductor bajo la específica
campo magnético no uniforme del imán utilizado en este estudio.

5. Conclusiones

El modelo propuesto en el presente trabajo ignora el campo


desmagnetizante pero proporciona expresiones analíticas que describen
varios aspectos de la fuerza entre un imán permanente y un
superconductor finito en estado mixto. Más precisamente, estas
cuestiones son: i) Las dimensiones del superconductor se incorporan
en el cálculo. en el limite de
un superconductor infinito, la fuerza repulsiva tiende a la
FIGURA 8. Fuerza horizontal total (línea continua) sobre el imán valor adecuado. De esta manera, el modelo nos permite cuantificar cómo
en función de su posición a lo largo del eje X. Los datos utilizados en el bueno/malo es la aproximación de tomar la muestra como infinita. ii) El
cálculo son: µ = 2.533×10−2 A m2 , µs = 6,95×10−3 A m2 , modelo de dipolo congelado se adaptó para que sea aplicable a un
R = S = 5 mm, L = 3,5 mm y a = 3,5 mm. La curva completa es superconductor con espesor finito. El momento magnético de este
antisimétrica en x. La recta discontinua corresponde a un armónico dipolo se obtuvo de la condición
comportamiento del oscilador. El recuadro superior muestra datos experimentales tomados del equilibrio mecánico observado experimentalmente. Este
de la referencia. 23, que concuerdan bien con el modelo. El
El valor fue consistente con nuestras mediciones del campo magnético.
El recuadro inferior muestra la variación del campo magnético medido en
campo atrapado en el superconductor. iii) Las diferentes formas
un punto fijo, cuando el imán oscila horizontalmente alrededor de su
en los que se puede magnetizar la muestra.
posición de equilibrio.
por un exponente n. Variando este parámetro fue posible reproducir el
empuje y su movimiento posterior fue monitoreado midiendo la típico acercamiento histerético/alejamiento.
campo magnético en un punto fijo sobre el imán, cerca de su Curvas de la fuerza vertical total en procesos FC y ZFC.
posición de equilibrio. La variación del campo con el tiempo es iv) También se calculó la fuerza horizontal total y se
se muestra en el recuadro inferior de la Fig. 8. Se observó que el Se observó una buena concordancia con resultados experimentales
anteriores.
El imán oscilaba con una frecuencia de sólo 7 Hz, es decir, uno
sexta parte del valor calculado. Esta medida fue confirmada con un
estroboscopio. Para considerar el efecto de Agradecimientos
El campo desmagnetizador sobre la magnetización en Meissner.
En este estado, se utilizó el resultado para un prisma rectangular Quiero agradecer el apoyo económico brindado por Cona­cyt (beca No.
superconductor bajo un campo externo uniforme: M = −2,5Ha [14]. 245453) y la ayuda técnica del Dr. Fidel
en lugar de M = −Ha, Ec. (2). Como resultado, la teoría Gamboa y Osvaldo Gómez.

1. FN Werfel et al., Supercond. Ciencia. Tecnología. 25 (2012) 014007. 9. Z. Deng, K. Tsuzuki, M. Miki, B. Felder, S. Hara y M.
2. F. Hellman, EM Gyorgy, DW Johnson Jr., HM O'Bryan, Izumi, Física C 471 (2011) 1459.
y RC Sherwood, J. Appl. Física. 63 (1988) 447.
10. M. Tomita y M. Murakami, Naturaleza 421 (2003) 334.
3. AA Kordyuk, J. Appl. Física. 83 (1998) 610.
11. E. Diez­Jimenez, E. Valiente­Blanco, y J.­L. Perez­Diaz, J.
4. CP Bean, Phys. Rev. Lett. 8 (1962) 250.
Supercond. Magn. de noviembre. 26 (2013) 71.
5. Xing­da Wu, Ke­Xi Xu, Yue Cao, Shun­bo Hu, Peng­xiang
Zuo y Guan­dong Li, Physica C 486 (2013) 17. 12. G. Iannone, S. Farinon, G. De Marzi, P. Fabbricatore y
U. Gambardella, IEEE Trans. Aplica. Supercond. 25 (2015)
6. JR Hull, Aplicaciones de levitación de superconductores de alta
8200107.
temperatura. En: AV Narlikar (ed.), Superconductividad de alta
temperatura 2: aplicaciones de ingeniería (Springer Verlag, 13. IA Rudnev, AI Podlivaev, IEEE Trans. Aplica. Supercond. 26
Berlín­Heidelberg 2004), págs. 91­142. (2016) 8200104.
7. C. Navau, N. Del Valle y A. Sánchez, IEEE Trans. Aplica.
14. C. Navau, CA Cardoso, OF de Lima y FM Araujo­Moreira, J. Appl.
Supercond. 23 (2013) 8201023.
Física. 96 (2004) 486.
8. C. Leblond, I. Monot, J. Provost y G. Desgardin, Physica C
311 (1999) 211. 15. J. Lugo y V. Sosa, Física C 324 (1999) 9.

Rev. Méx. Fis. 63 (2017) 497­504


Machine Translated by Google

504 V.SOSA

16. MK Alqadi, FY Alzoubi, HM Al­Khateeb y NY Ayoub, 20. HM Al­Khateeb, FY Alzoubi, MK Alqadi y NY Ayoub, Turk. J. Física. 31
Modificación. Física. Letón. B 20 (2006) 1549. (2007) 271.

17. FY Alzoubi, MK Alqadi, HM Al­Khateeb y NY Ayoub, IEEE Trans. 21. KW Yung, PB Landecker y DD Villani, Magn. eléctrico.
Aplica. Supercond. 17 (2007) 3814. Separar. 9 (1998) 39.

18. SM Al­Khateeb, MK Alqadi, FY Alzoubi y NY Ayoub, 22. TH Johansen y H. Bratsberg, J. Appl. Física. 74 (1993) 4060.
Física china. Letón. 24 (2007) 2700. 23. T. Hikihara y G. Isozumi, Physica C 270 (1996) 68.

19. HM Al­Khateeb, MK Alqadi, FY Alzoubi y NY y Ay­oub, J. Supercond. 24. Yong Yang y Xiaojing Zheng, Supercond. Sci. Technol. 21 (2008)
Nov. Magn. 21 (2008) 93. 015021.

Rev. Méx. Fis. 63 (2017) 497­504

También podría gustarte