Download Ecuaciones y Energía del campo magnético

1
Temas teóricos
Electromagnetismo
Ecuaciones y Energía del campo magnético. Lino Spagnolo.
A semejanza de la Electrostática, en Magnetostática se tienen las siguientes
ecuaciones básicas:
Electrostática
 S D  dS  Q
 E  dl  0
  D  v
Magnetostática
 S B  dS  0
 B  0
 B  dl  o I  o S J  dS
 B  dl  S   B  dS  o I
 E  0
Ley de Ampère   B  o J
E  V ; V( r ) 
V( r ) 
V( r ) 
1
4 o
1
4 o
( r ')
1
dr '
4 o l r  r '
 s ( r ')
S r  r ' dS '
v ( r ')
v r  r '
  J  0 B   A   A  0
A( r ) 
 J  0
A( r ) 
dv '
A( r ) 
o
4
o
4
v
S
o
4
l
I dl '
r r'
K ( r ') dS '
r r '
J ( r ') dv '
r r '
2
Campo magnético creado por un
elemento de corriente en un punto P.
Ley de Biot-Savart.
dH ( r ) 
dB( r ) 
I ( r ') dl ' (r  r ')
4 r  r '
3
o I ( r ') dl ' (r  r ')
3
4
r r '
Si en lugar de un elemento de corriente es
una densidad volumétrica de corriente, el
campo es:
B( r ) 
o
4
(r  r ')
v J (r ')  r  r ' 3 dv '
Fuerza sobre un elemento de
corriente debido a la presencia de un
B.
dF  I dl  B
campo magnético
Fuerza entre elementos de corriente
de dos circuitos eléctricos.
Ley de Ampère para las fuerzas.
dF2 
o
dl  (r  r )
I1I 2 dl2  1 2 3 1
4
r2  r1
Efectuando el doble producto vectorial, queda la integral:
3
F2  
o
I1I 2
4
(r2  r1 )(dl2  dl1 )

r2  r1
3
La Ley de Ampère puede obtenerse también a partir de la Ley de Biot-Savart, como
comprobación de que ambas leyes son complementarias y que pueden usarse una o la
otra según mejor convenga en la solución de cada problema.
A partir de la ecuación general del campo creado por una densidad volumétrica
de corriente:
o
4
B( r ) 
Dado que
v
J ( r ') 
(r  r ')
r r '
3
dv '
 1 



 y llamando R  r  r ' podemos escribir:
3
r

r
'
r r '



1
B( r )   o  J ( r ') ( )dv '
4 v
R
(r  r ')
Teniendo en cuenta la propiedad vectorial del Rotor del producto de una función escalar
por otra vectorial:
 (G)   ( G)    G
1
queda:
R
 J ( r ')  1
1
( )  J ( r ')    
   J ( r ')
 R  R
R


Además, para corrientes estacionarias el rotor es nulo, o sea:   J ( r ')  0
Asimilando que
G  J ( r ') y  
Por lo tanto el campo será:
B( r ) 
o
4
v

J ( r ')
R
dv '   
o
4
v
J ( r ')
R
Expresión que es precisamente el potencial vectorial magnético
A( r ) 
o
4
v
J ( r ') dv '
r r '
dv '
(**)
A.
 B   A
Debe observarse que la expresión del campo vectorial A está definida a menos de una
función gradiente, puesto que si definimos un nuevo campo vectorial:
4
A( r ) 
o
4
v
J ( r ') dv '
r r '
Seguiremos teniendo el mismo campo magnético:
Dado que el rotor del gradiente es nulo.
  ( r )
B   A
Si ahora hallamos el rotor de la expresión (**):
  B( r )     A 
o
4
v
  
J(r )
R
dv ' 
o
4
v    G dv '
Considerando nuevamente la propiedad vectorial:
 A  (  A)  2 A
(***)
Como se vio, este potencial vectorial está definido a menos de un gradiente.
Además, por el teorema de Helmholtz, un campo vectorial como el
se conocen su rotor y su divergencia.
B
estará definido si
Como su rotor está definido, su divergencia podemos fijarla con libertad a partir
de la ecuación (***). Con ello se fijará una condición a cumplir por el potencial
vectorial A , originalmente llamada condición de Coulomb, pero que posteriormente se
vio que este tipo de condiciones aparecería en muchas partes del Electromagnetismo, ya
sea clásico como cuántico.
En la teoría de los campos, hay un importante concepto llamado simetría de
Gauge, que define la propiedad de las ecuaciones que los caracterizan, de no cambiar la
naturaleza de los campos si se les aplica una determinada operación a todos los
componentes del mismo y en todos los puntos del espacio. Los campos más conocidos
con tal propiedad son el gravitatorio y el electromagnético.
Así por ejemplo, si al campo eléctrico definido por
campo magnético definido por
H
1

E
A
  y al
t
  A , redefinimos al potencial escalar 
mediante la operación:
(r ; t )   '(r ; t )  (r ; t ) 
Y redefinimos al potencial:
A(r ; t )  A '(r ; t )  F (r ; t )
Se puede verificar que los respectivos campos
F (r ; t )
t
E y H en nada cambiarán, bastando
para ello reemplazar los nuevos valores  ' y A '. . Esta invariancia se define como
Invariancia de Gauge ante una operación o transformación Gauge.
5
Cuando se está tratando con campos Gauge como el Electromagnético, y por
alguna razón de simplificación matemática de las ecuaciones de campo, resulte
conveniente modificar su estructura, siempre se podrá recurrir a una operación Gauge, o
condición Gauge, que permita su reacomodamiento, o como ahora se denomina el
término: recalibración de escala.
Tal condición fija un nuevo dimensionamiento del campo vectorial, siempre con
la tendencia de simplificar su expresión. A esto se le denominó condición de
calibración o de escala, más conocido por condición de Coulomb.
(Más adelante se hallará un nueva condición de calibración o de escala, conocida
como condición de L. V. Lorenz.)
En la expresión (***) se fija, como condición de Gauge, que la divergencia del
campo potencial vector
A sea nula:    A  0 A fin de simplificar la expresión.
Por lo cual sólo nos quedará la ecuación:
  B( r )     A   2 A
O sea, en la forma:
  B( r ) 
o
4
v
   (
J(r )
Puede demostrarse que el Laplaciano
Y la expresión de
en:

o
4
v
2 (
Con lo cual:
R
)dv '  
o
4
v
2 (
J(r )
R
)dv '
1
2 ( )  4  R  4  (r  r ')
R
 (r  r ') corresponde a la delta de Dirac, cuya integral se convierte

1
J ( r ') )dv   o
R
4
2 A 
o
4
v
v
2 (
J ( r ') (4  (r  r '))dv 
1
J ( r ') )dv   o J ( r ')
R
Quedando finalmente la integral del rotor como:
  B  o J ( r ')
Que es la Ley de Ampère en formato de ecuación diferencial.
o
4 J ( r ')  o J ( r ')
4
6
Energía del campo magnético.
Analizaremos ahora la energía del campo magnético la cual ofrece también
varios enfoques de análisis pero sus resultados son más inmediatos y directos que para
el campo eléctrico.
Comenzaremos el estudio para el caso simple de un circuito eléctrico con un
elemento magnético como el solenoide para analizar el almacenamiento de la energía
magnética en ese elemento.
La ecuación del circuito eléctrico para el caso de tensión aplicada constante y
dI
dt
T
T 2
T dI

W

E
I
dt

R
I
dt

L
0 o
0
0 I dt dt

Y la energía eléctrica está dada por: 
W  1 LI 2  R T I 2 dt
(01)
T
0

2
una corriente que a partir del instante
t  0 de la conexión, es: Eo  RI  L
Fig. 1
T
En la cual el término
el otro término Wm
R  I 2 dt es la energía consumida por efecto Joule, mientras que

0
1 2
LIT corresponde a la energía del campo magnético.
2
Resolviendo la ecuación diferencial de la carga del circuito será:
RI  L
I  0
dI
 E0 y en t  0 
dt
 Eo  Cte.
R
 t
dI
R t
Operando: 
   dt  I  IT (1  e L )
0
E
L 0
(I  o )
R
En la cual IT  E0 / R
1 2
La energía magnética almacenada en el solenoide es Wm  L IT
2
I
(02)
(03)
7
Se le puede dar otro formato ya que para el caso especial de una bobina o inductor con
N vueltas, se conocen sus parámetros tales como su inductancia L y el campo
magnético B generado.
En efecto, para solenoides o bobinas muy largas (longitud l ) y de pequeño radio
inductancia viene dada por la fórmula:
r
, su
NI
N 2S
L  o
Donde S   r 2 El campo magnético está dado por: B  o
l
l
Por lo cual la energía magnética resultará:
1 2 1 B2
1 B2
Wm  L I 
Sl 
vol. ( J )
2
2 o
2 o
(04)
En consecuencia la densidad de energía magnética resulta:
B2
J
m 
( 3)
2 o m
(05)
Energía Magnética del campo.
Plantearemos ahora el análisis de la energía magnética desde el aspecto de un
campo magnético generado por un corriente constante o estacionaria.
Como el elemento magnético más común es la bobina (inductor o solenoide),
analicemos el flujo generado por una espira en la cual circula una corriente I .
i   B  dSi 
S
La
l A dl
con B   A
fem inducida, según la ecuación de Faraday, es: fem  V  
La energía suministrada por la batería para que circule una corriente
dW  V I dt  I
d
dt
I por la espira, es:
d
dt  dW  I d  dWm
dt
(06)
Esta expresión de la energía será equivalente a la magnética en el caso en que el
medio no tenga histéresis ni la espira tenga un movimiento mecánico ya que en tal caso
se gastará parte en producir ese movimiento.
Si imaginamos que existen otras espiras adyacentes formando un solenoide con
las características antes detalladas, los campos magnéticos generados por cada espira
forman un conjunto sumable de valor:
N
dW   I i di  dWm
(07)
i 1
Recurriendo al concepto de inductancia mutua de las espiras, la misma se define como:
8
k  Lk I k
(08)
Para cada una de las espiras diferente de la primera.
Con lo cual puede integrarse la energía pero teniendo presente que el cálculo debe
efectuarse para la variación de la corriente desde i  0 hasta i  I , su valor final.
Si en el medio de las espiras no hay histéresis, por lo cual puede asumirse como lineal,
la variación de la corriente puede tomarse como lineal (la energía no depende de la
forma como se llegue al valor final), es decir:
I' I
Y como el flujo es proporcional a la intensidad:
d '   Lk dI k'   Lk I k d    d 
k
(09)
k
'
'
 d W    Ii di    Ii  i d    Ii i   d 
1
i
i
dW 
i
0
1
 Ii i  Wm
2 i
(10)
Fórmula similar a la obtenida para la energía electrostática:
1
 Qi Vi  We
2 i
Con este resultado se puede hallar la energía magnética en función del
campo:
Por las ecuaciones de Maxwell: m
  B  nˆ dS   A  nˆ dS   A  dl
S
S
l
Se puede poner:
Wm 
Reemplazando
1
1
1
I i i   I i   A  dl  

 2
2 i
2 i  l
I dl por J dv
Wm 
1
2
l I A  dl
(11)
, elemento de volumen:
1
l A  dl I  2 v J  Adv
 H  J
1
1
Wm   J  Adv   (  H )  Adv
2 v
2 v
Y con la propiedad vectorial:   ( H  A)  A  H  H  A
Utilizando la ecuación de Maxwell-Ampère:
(12)
9
1
1
  ( H  A) dv   H   Adv

2 v
2 v
1
1
Wm   ( H  A)  nˆ dS   H  B dv
2 S
2 v
Wm 
Reemplazando:
(13)
Como la integral es sobre todo el volumen del espacio, la integral de superficie debe
extenderse a una esfera con radio tendiente a infinito. Pero en tal caso como los campos
AyH
son proporcionales al producto
proporcional a
1 1

r r2
mientras que la superficie sólo es
r 2 , la integral de superficie tiende a cero.
La última integral se reduce a:
Wm 
1
H  B dv ( J )
2 v
(14)
Y la densidad de energía magnética será:
1
B2 1
J
m  H  B 
 o H 2 ( 3 )
2
2 o 2
m
(15)
Otra forma sencilla de expresar la Energía Magnética del campo.
 C al flujo magnético total o concatenado a través de una

superficie S , y al valor  m  C como el flujo magnético correspondiente a una
N
espira, N  1, por la ecuación de Faraday-Maxwell se tiene que la fem inducida es:
d C
(16)
fem    
dt
Además el flujo se relaciona con el vector “densidad de flujo magnético B ” por
Si denominamos
la fórmula:
 m   B  dS   B  nˆ dS
S
S
(17)
Y si el campo B es uniforme y constante, como en el interior de una bobina
suficientemente larga, entonces:
C  NBS    
d C
dB
  NS
(V ).
dt
dt
(18)
10
Por definición, la inductancia mutua de una bobina, o solenoide, recorrida por
una corriente uniforme I , vale:
L
C
d C
Wb
dI
(en Hy o
)   
 L
I
A
dt
dt
(19)
Por lo cual, la energía definida como:
dW   I dt  L
dI
I dt  L I dI
dt
Integrando:
W   L I dI 
1 2
LI
2
(20)
Valor de la energía magnética en el interior del solenoide. Como vimos antes, esta
energía puede atribuirse al campo magnético existente en su interior.
La fórmula anterior es útil para hallar el valor de la inductancia mutua del
solenoide. El valor del campo magnético en su interior se halló como:
N
dB
I y como    NS
l
dt
Con la relación constitutiva B   H , reemplazamos:
H
dH
 N 2 S dI
dI
   NS 

 L
dt
l dt
dt
De lo cual:
L
 N 2S
l
(21)