Teoría Electromagnetismo
Teoría Electromagnetismo
Teoría Electromagnetismo
1
Ilustre ilustrador.
2
Los demás autores han hecho contribuciones menores y aparecen mencionados tras el índice.
Índice general
1. El campo electroestático 5
1.1. Principio de superposición . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.2. Ley de Coulomb . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.2.1. ¿Por qué sabemos que el denominador es r2 ? . . . . . . . . . . . . . . . . . . . . . . . 7
1.2.2. Comentarios sobre los límites de ley de Coulomb . . . . . . . . . . . . . . . . . . . . . 7
1.3. Campo eléctrico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.3.1. Caso de una carga puntual . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.3.2. Caso de varias cargas puntuales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.3.3. Distribuciones de carga . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.3.3.1. Distribución volumétrica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.3.3.2. Distribución superficial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.3.3.3. Distribuciones lineales de carga . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.3.3.4. Comentario final . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
1.4. Teorema de Gauss . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
1.4.1. Diferenciales en coordenadas cilíndricas . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.4.2. Repaso de matemáticas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.4.3. Flujo eléctrico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
1.4.4. Laplaciano de 1r . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
1.4.5. Delta de Dirac . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
1.4.6. Forma diferencial del Teorema de Gauss . . . . . . . . . . . . . . . . . . . . . . . . . . 22
1.4.7. Obtención de la forma integral a partir de la forma diferencial . . . . . . . . . . . . . . 25
2. El potencial electroestático 28
2.1. Para una carga puntual . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
2.2. Para una distribución de carga . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
2.3. El significado físico del potencial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
2.4. Ecuaciones de Poisson y Laplace . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
2.4.1. *Resolución en cartesianas (No entra): . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
2.5. Teorema de Green . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
2.5.1. Teorema de Green en física . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
2.5.2. Teorema de unicidad de la solución de la ecuación de Poisson . . . . . . . . . . . . . . 40
3. Distribuciones de dipolos 44
3.1. El dipolo eléctrico. Momento dipolar: . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
3.2. Potencial generado por un dipolo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
3.3. Campo que genera un dipolo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
~ y φ generado por una distribución de dipolos
3.4. E . . . . . . . . . . . . . . . . . . . . . . . . . . 50
1
Laín-Calvo ÍNDICE GENERAL
5. Medios dieléctricos 84
5.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
5.2. Polarización dieléctrica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
5.3. Campo en el interior de un dieléctrico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
5.4. Vector desplazamiento eléctrico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88
5.5. Constante dieléctrica y susceptibilidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
5.6. L.H.I. (Materiales lineales homogéneos e isótropos) . . . . . . . . . . . . . . . . . . . . . . . . 90
5.7. Condiciones frontera entre dos medios dieléctricos . . . . . . . . . . . . . . . . . . . . . . . . . 92
5.7.1. Componente normal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
5.7.2. Componente tangencial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 94
6. Sistemas de conductores 96
6.1. Curiosidades introductorias . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
6.2. Equilibrio electroestático (propiedades de conductores) . . . . . . . . . . . . . . . . . . . . . . 96
6.3. Teorema de reciprocidad e identidad de Gauss . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
6.3.1. Casos concretos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 104
6.4. Relación entre los potenciales y las cargas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 104
6.5. Influencia total: condensadores . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 106
Capítulo 1
El campo electroestático
F~23
F~13
q3
q1
q2
Entonces:
1 q1 q3
F~13 = 2 r̂
4πε0 r13
1 q2 q3
F~23 = 2 r̂
4πε0 r23
Además:
Proposición 1. El principio de superposición afirma que cuando las ecuaciones de comportamiento que rigen
un problema físico son lineales, entonces el resultado de una medida o la solución O de un problema práctico
relacionado con una magnitud extensiva asociada al fenómeno, cuando están presentes los conjuntos de factores
causantes I1 , . . . , Ik con k ∈ N, puede obtenerse como la suma de los efectos de cada uno.
Demostración. Sea f : Rn −→ Rm una aplicación lineal que relaciona dos magnitudes físicas I (input) y
O (output). Supongamos, además, que contamos con k ∈ N subsistemas tales que cada uno produce una
magnitud física de entrada Iu . Queremos calcular cuál es el valor de O. Para ello:
k
!
X
O = f (I1 + · · · + Ik ) = f Iu
u=1
Y, como f es lineal por hipótesis:
k
X
O= f (Iu )
u=1
Q.E.D.
q0
~0
R ~r
q
~
R
1 q0q 1 q0q
F~q0 →q = r̂ = ~r
4πε0 r2 4πε0 r3
donde r̂ = |~~rr| y ε0 ' 8, 85 m
F
es una constante llamada permitividad eléctrica del vacío y r es la distancia
que separa ambas cargas. Nótese que esta fuerza será atractiva si las cargas son del mismo signo y será repulsiva
en caso contrario.
F~q→q0 = −F~q0 →q
1. Cavendish: |ε| ≤ 0, 02
~
~ := lı́m F
E
q0 →0 q0
q2
F~30
q1
q0 F~10
F~20
q5
q3
q4
El concepto de límite está para indicar que la nueva carga añadida es tan pequeña que su interacción con
el resto de cargas del sistema es despreciable. Nótese que la propia carga q0 es necesaria para el cálculo de F~ ,
0
luego el límite es, en realidad, una indeterminación .
0
1 q 0 q0
F~ = r̂
4πε0 r2
Aplicando la definición:
1 q 0 q0
~ R~ = lı́m 4πε0 r2 r̂ 1 q0
E = r̂
q0 →0 q0 4πε0 r2
Q.E.D.
Notación 1. Usaremos ’ para denotar las cargas generadoras del campo. Utilizaremos R ~ 0 para referirnos a
las posiciones de las cargas respecto a nuestro sistema de coordenadas, mientras que usaremos ~r para las
distancias relativas.
~ −R
donde ~ri = R ~ 0.
i
Ejemplo 1. Imaginemos:
q3
~r20
X
q0
(x, y, z) r̂10
En este caso:
~ R~ = 1 q1 q2 q3
E 2 r̂10 + r 2 r̂20 + r 2 r̂30
4πε0 r10 20 30
donde:
2 2 2
2
r10 = x0 − x01 + y0 − y10 + z0 − z10
En consecuencia:
q
|~r10 | = (x0 − x01 )2 + (y0 − y10 )2 + (z0 − z10 )2
Por ende:
Nota 1. En general se evitará el uso del subíndice 0. Es decir, se usará r12 en vez de r10
2 ; x, en vez de x y así,
0
sucesivamente.
V0
dV 0
~0
R ~r
~
R
Y
X
~ =R
Nótese que R ~ 0 + ~r ⇔ ~r = R
~ −R
~ 0 . La densidad de carga del sólido viene dada por una función:
~ 0 = dQ ⇔
ρ R
dV 0
⇔ dQ = ρ R~ 0 dV 0
El diferencial de volumen que tomo va a ser lo suficientemente pequeño para poder considerarlo puntual.
De esta forma:
ρ R~ 0 dV 0
dE~ R ~ = 1 r̂
4πε0 r2
A continuación, para hallar el campo total integramos al volumen:
˚ ρ R ~0
~ ~ 1
E R = 2
r̂dV 0
4πε0 V 0 r
Q.E.D.
S0
dS 0
~0
R ~r
~
R
Y
X
~ es:
El campo que genera el diferencial de carga en el punto R
σ ~0
R
~ R~ = 1
dE r̂dS 0
4πε0 r2
Y, de forma similar al caso volumétrico:
¨ ~0
σ R
~ R~ = 1
E r̂dS 0
4πε0 S0 r2
Q.E.D.
dl0
Todo es semejante a los casos anteriores:
dQ
~0 =
λ R
dl0
λ ~0
R
~ ~ 1
dE R = m r̂dl0
4πε0 r2
Integrando a ambos lados llegamos a:
ˆ λ R
~ 0 r̂
~ R~ = 1
E dl0
4πε0 l0 r2
Q.E.D.
dz 0 ~
dE
~
R
dz 0 ~
dE
Sabemos:
dQ = λdz 0
El campo creado por el diferencial de carga es:
~ = 1 λdz 0
dE r̂
4πε0 r2
Por simetría sabemos que se va a anular el campo en el eje z, pues por cada dz 0 que existe por encima del
punto, existe otro dz 0 por debajo; de manera que los campos que se generan en la componente z se anulan
entre sí.
Por tanto, definamos:
~ cos α
dER := dE
R
Como cos α = :
r
1 λdz 0 R
dER (R) =
4πε0 r2 r
El resultado final es:
ˆ ∞
1 λR 0
ER (R) = 2 dz
0 4πε0 r3
Como r2 = R2 + z 02 :
ˆ ∞ ∞
2λR dz 0 λR z0
ER (R) = 3 = √ =
4πε0 0 (R2 + z 02 ) 2 2πε0 R2 R2 + z 02 z 0 =0
z 0 →∞
−−z}|{
−−→1
z 0
λR z 0 λ
= r − 0
= 2πε0 R
2πε0 R 2
R2 +1
z 02
| {z }
−−0−−→R2
z →∞
S
S
φ
φ
~ siempre va hacia fuera de la superficie.
En una superficie cerrada, el vector S
Teorema 1 (Teorema de Gauss (forma integral)). La integral del campo eléctrico a lo largo de una superficie
cerrada es igual a la carga total encerrada dentro de dicha superficie partido por la permitividad eléctrica del
vacío.
‹
E ~ = QT
~ · dS
SC ε0
donde SC es una superficie cerrada, QT es la carga encerrada dentro de la superficie y ε0 es la permitividad
eléctrica del vacío.
‹ ¨ ¨ ¨
~ · dS
E ~= ~ · dS
E ~ + ~ · dS
E ~ + ~ · dS
E ~ =
| {z } | {z } | {z }
SC tapa superior tapa inferior superficie lateral
~
E⊥d ~
S⇒=0 ~
E⊥d ~
S⇒=0 ~ S⇒=E(r)dS
Ekd ~
¨
= E (r) dS
superficie lateral
Como el campo es constante a lo largo de la superficie y la superficie lateral del cilindro es S = 2πrh:
¨
E (r) dS = E (r) 2πrh
superficie lateral
dlcírculo = rdϕ
dlaltura = dz
dradio = dr
De superficie:
dSlateral = rdϕdz
dStotal = rdrdϕ
De volumen:
dV = rdrdϕdz
î ĵ k̂
~≡∇
~ ×A
~ := ∂ ∂ ∂
rotA , , × (Ax , Ay , Az ) = ∂
∂x
∂
∂y
∂
∂z
∂x ∂y ∂z
Ax Ay Az
Definición 9. Llamamos ángulo sólido al ángulo espacial tal que su expresión diferencial es:
r̂ · dS~ ~r · dS~
dΩ := 2
= 3
r r
Gráficamente:
~
dS
Z
~r
dΩ
La idea del ángulo sólido1 es el ángulo 2D (la «anchura» del cono) que se forma tras proyectar la superficie
sobre una esfera de radio unidad. Es decir, es una medida del agujero que tendría que hacer en la esfera de
~ Nótese que efectivamente es un ángulo, pues es un parámetro adimensional.
radio unidad para poder ver el dS.
Teorema 3 (Teorema de Stokes). Sea Ω un abierto en R3 y sea A ~ : Ω −→ R3 una función vectorial tal que
∂Ai ~ a lo largo de la
∃ ∀i, j = x, y, z y sea S una superficie abierta. Entonces la integral del rotacional de A
∂j
1
Más información en el siguiente enlace: https://en.wikipedia.org/w/index.php?title=Solid_angle&oldid=880143402.
superficie S tiene el mismo valor que la integral de A~ a lo largo de la curva C que delimita la superficie abierta
S. ¨ ˛
~ ~
∇ × A · dS = ~ ~ · d~l
A
S C
Observación 5. Siempre debemos recorrer la curva siguiendo el sentido dado por la regla de la mano derecha,
~
donde el pulgar apunta en la dirección de dS.
~ = 1 q
E r̂
4πε0 r2
De manera que:
‹ ¨ ~ ¨
~ · dS
~= q r̂ · dS q
Φ= E 2
= dΩ
SC 4πε0 | r{z } 4πε0
=dΩ
Q.E.D.
1
Proposición 7. Siempre que sea r 6= 0, el gradiente de respecto a las coordenadas sin primar puede
r
expresarse como:
~ 1 = −~r = − r̂
∇
r r3 r2
Demostración. Expresamos r en función de las coordenadas x, y, z, x0 , y 0 , z 0 :
q
r = |~r| = (x − x0 )2 + (y − y 0 )2 + (z − z 0 )2
∂ 1 1 1 0
=− h i 3 2 x − x =
∂x r 2 0 2 0 2 0 2 2
(x − x ) + (y − y ) + (z − z )
x0 − x
=h i3
x0 )2 y 0 )2 z 0 )2
2
(x − + (y − + (z −
Análogamente:
∂ 1 y0 − y
=h i3
∂y r
x0 )2 y 0 )2 z 0 )2
2
(x − + (y − + (z −
∂ 1 z0 − z
=h i3
∂z r
(x − x0 )2 + (y − y 0 )2 + (z − z 0 )2
2
Luego:
=−~
r
z }| {
0 0 0
~ 1 x − x, y − y, z − z
∇ = 3 =
r 2
0
2 0
2 0
2
x−x + y−y + z−z
| {z }
=r2
−~r r̂
= 3
=− 2
r r
Q.E.D.
1
Proposición 8. Siempre que sea r 6= 0, el gradiente de respecto a las coordenadas primadas puede expresarse
r
como:
~ 0 1 ~r r̂ ~ 1
∇ = 3 = 2 = −∇
r r r r
Demostración. La demostración es casi idéntica a la anterior, sólo que cambian los signos de las derivadas
parciales. Q.E.D.
1
1.4.4. Laplaciano de r
1
Proposición 9. Siempre que sea r 6= 0, el laplaciano de respecto a las coordenadas sin primar puede
r
expresarse como:
2 1
∇ =0
r
Demostración. Por la proposición 7 en la página anterior:
~ 1 (x0 − x, y 0 − y, z 0 − z)
∇ =h i3
r 0 2 0 2 0 2 2
(x − x ) + (y − y ) + (z − z )
Ahora, sólo queda hallar la divergencia. Para ello, hallemos las derivadas parciales:
∂ x − x0
∂x h i3 =
(x − x0 )2 + (y − y 0 )2 + (z − z 0 )2
2
h i3 h i1
x0 )2 y 0 )2 z 0 )2 2 2 2 2
2
(x − + (y − + (z − − (x − x0 ) 32
0 0 0
(x − x ) + (y − y ) + (z − z ) 2 (x − x0 )
= h i3 =
0 2 0 2 0 2
(x − x ) + (y − y ) + (z − z )
(x − x0 )2 + (y − y 0 )2 + (z − z 0 )2 − 3 (x − x0 )2
= h i5 =
0 2 0 2 0 2 2
(x − x ) + (y − y ) + (z − z )
r2 − 3 (x − x0 )2
=
r5
Análogamente:
∂ y− y0
∂y
h i3 =
(x − x0 )2 + (y − y 0 )2 + (z − z 0 )2
2
r2 − 3 (y − y 0 )2
=
r5
∂ z− z0
∂z h i3 =
(x − x0 )2 + (y − y 0 )2 + (z − z 0 )2
2
r2 − 3 (z − z 0 )2
=
r5
Podemos ver que:
1 3 h 2 2 2 i 3
∇ = 5 r 2 − x − x0 + y − y 0 + z − z 0
2
= 5 r2 − r2 = 0
r r r
Q.E.D.
Y
1
a
x0 X
δ (x − x0 ) := lı́m f (x − x0 )
a→0
Demostración. ˆ ∞ ˆ ∞
def
δ (x − x0 ) g (x) dx = lı́m f (x − x0 ) g (x) dx =
−∞ −∞ a→0
ˆ x0 + a2 ˆ x0 + a2 ˆ x0 + a2
1 1
= lı́m f (x − x0 ) g (x) dx = lı́m g (x) dx = lı́m g (x) dx =
a→0 x − a
0 2
a→0 x − a
0 2
a a→0 a x0 − a2
1
= lı́m ag (x0 ) = g (x0 )
a→0 a
Q.E.D.
~δ (~r − ~r0 ) := δ (x − x0 ) δ (y − y0 ) δ (z − z0 )
pues la integral a todos los posibles ángulos sólidos de dΩ es 4π. A continuación, nótese que por la proposi-
ción 10 en la página anterior:
˚
−4π = −4πδ (~r) dV 0
V0
Por tanto:
˚ ˚
1
∇2 dV 0 = −4πδ (~r) dV 0
V0 r V 0
Como lo anterior debe cumplirse para cualquier volumen V 0 , llegamos a que debe ser:
1
∇2 = −4πδ (~r)
r
Q.E.D.
Proposición 12. La forma integral del teorema de Gauss implica su forma diferencial. Es decir, si el flujo
del campo eléctrico a lo largo de una superficie cerrada SC es la carga encerrada en dicha superficie partida
~ dentro de dicha superficie SC debe
por ε0 , entonces la divergencia del campo eléctrico en cualquier punto R
~
ser la función densidad de carga ρ evaluada en dicho punto R partido por ε0 .
‹ ρ ~
R
E ~ = QT ⇒ ∇
~ · dS ~ ·E~ R ~ =
SC ε0 ε0
Demostración. Primero,
Recordemos, ahora, que, por la proposición 2 en la página 8, el campo generado por una distribución de
carga es:
˚ ρ R ~ 0 · r̂
~ ~ 1
E R = dV 0
4πε0 V0 r2
Q0
dV 0
~r
~0
R
~
R
Nótese que como la divergencia y la integral actúan respecto a coordenadas diferentes, el orden de los opera-
dores no importa, es decir, conmutan. De forma que puedo cambiar el orden:
˚ ρ R ~ 0 · r̂ ˚
~ ~ ~ 1 ~ 0 1 ~ 0 ~ r̂
∇·E R = ∇· 2
dV = ρ R ∇· dV 0
4πε0 V0 r 4πε0 V0 r2
~ −R
Recordando ~r = R ~ 0: ˚
1
~ ~ ~
∇·E R = ~0 δ R
ρ R ~ −R
~ 0 dV
ε0 V0
Por último, por la propiedad de traslación de la delta (ver proposición 10 en la página 21):
ρ R ~
~ ~
∇·E R = ~
ε0
Q.E.D.
Teorema 4 (Forma diferencial del teorema de Gauss). En cualquier punto R ~ ∈ R3 se cumple que la divergencia
~ coincide con el valor de la función densidad de carga evaluada
del campo eléctrico evaluada en dicho punto R
~ partido por ε0 .
en el punto R
~
ρ R
~ ~ ~
∇·E R =
ε0
Observación 8. Una divergencia positiva indica que las líneas de campo salen de la carga, mientras que una
divergencia negativa indica lo contrario.
˚ ρ R ~0 R ~ −R ~0
E~ R~ = 1 d3 ~0
R
4πε0 3 |{z}
V R~ −R~0 =dV
donde ρ R~ 0 es la función densidad de carga. Si hacemos la divergencia a ambos lados de la ecuación y
usamos la proposición 11 en la página 22, llegamos a:
˚
~ ~ ~ 1 ~0 δ R ~ −R ~ 0 d3 R~0
∇·E R = ρ R
ε0 V
Y lo anterior, por la proposición 10 en la página 21, es equivalente a:
ρ R ~
~ ~
∇·E R = ~
ε0
Q.E.D.
2
https://es.wikipedia.org/w/index.php?title=Ley_de_Gauss&oldid=107400507#Forma_diferencial_de_la_ley_de_
Gauss
Corolario 2. Gracias a las proposiciones 12 en la página 22 y 13, sabemos que las formulaciones del teorema
de Gauss en su forma integral (ver teorema 1 en la página 14) y en su forma diferencial (ver teorema 4 en la
página anterior) son equivalentes.
~ generado por una distribución volumé-
Proposición 14. En electroestática, la integral del campo eléctrico E
trica de cargas a través de una curva cerrada C es siempre nula.
˛
~ R
E ~ · d~r = 0
C
Demostración. Por la proposición 3 en la página 10, sabemos que el campo eléctrico generado por una distri-
bución volumétrica de cargas viene dado por la expresión:
˚ ρ R ~0
~ ~ 1
E R = 2
r̂dV 0
4πε0 V 0 r
Como las integrales se realizan a coordenadas diferentes, la integral triple conmuta con la realizada a lo largo
de la línea:
˛ ˚ ˛
~ R~ · d~r = 1 ~0 r̂ 0
E ρ R 2
·d~
r dV
C 4πε0 V0 r
C |{z}
~
=∇( r )
1
Q.E.D.
~ generado por una distribución volumétrica de carga es conservativo. Es
Corolario 3. El campo eléctrico E
decir, se da:
~ ×E
∇ ~ R~ = ~0
C
S2 Y
S1
X
Como lo anterior se cumple para cualquier superficie abierta S que se apoye en la curva C, necesariamente
el núcleo de la integral debe ser nulo:
~
∇×E R~ ~ = ~0
Como, además, el campo eléctrico generado por una distribución volumétrica de cargas (ver proposición 3
~ el campo eléctrico generado por una distribución de
en la página 10) únicamente depende de la posición R,
cargas es conservativo. Q.E.D.
Observación 9. El corolario 3 en la página anterior, nos va a permitir replantear todo lo visto hasta la fecha
utilizando funciones escalares (la función potencial). De eso es de lo que va a tratar el siguiente capítulo.
Capítulo 2
El potencial electroestático
Definición 12. Llamamos potencial eléctrico φ R ~ a la función escalar que cumple −∇φ ~ = E.
~ Esta
~ cualquiera, siempre que sea conservativo.
definición tiene carácter general, para un campo E
1 q0
~ =
φ R
4πε0 r
Demostración. Tenemos:
q0 ~r
Z
~0 ~
R
R
X
Por la proposición 2 en la página 8, el campo generado por una carga puntual lo podemos expresar como:
1 0 r̂
~ R
E ~ = q
4πε0 r 2
|{z}
~ (1)
=−∇ r
Por analogía con la definición de potencial eléctrico (ver definición 12 en la página anterior):
E~ R ~ = −∇φ ~ R~
~ es la función escalar:
Donde φ R
1 q0
~ =
φ R
4πε0 r
Q.E.D.
˚ ~0
ρ R
~ = 1
φ R dV 0 (2.2.0.1)
4πε0 V0 r
Demostración. Fijémonos en el dibujo:
V0
~r
Z
~0 ~
R
R
X
~ viene dado por:
Por la proposición 3 en la página 10, sabemos que el campo en un punto R
˚ ρ R ~ 0 r̂dV 0
~ ~ 1
E R =
4πε0 V0 r2
Por la proposición 7 en la página 18, tenemos:
˚ 1
~ R~ =− 1 ~
~0 ∇
E ρ R dV 0
4πε0 V 0 r
Como la integral se hace respecto a las coordenadas primadas y el gradiente se hace respecto a las coordenadas
sin primar, ambos operadores conmutan:
˚ ρ R ~0
~ R ~ = −∇ ~ 1
E dV 0
4πε0 V 0 r
| {z }
~
=:φ(R)
Es decir:
˚ ~0
ρ R
~ = 1
φ R dV 0
4πε0 V0 r
~ R
E ~ = −∇φ
~ ~
R
Q.E.D.
Observación 11. Nótese que diferenciando a ambos lados de la expresión para el potencial dada por la propo-
sición 16 en la página 29, tenemos:
ρ ~0
R
dφ R ~ = 1 dV 0
4πε0 r
Usaremos esto para calcular el potencial generado por distribuciones de carga volumétricas en la práctica.
~ 2)
(R
q
~ 1)
(R
Demostración. Por el teorema del trabajo en su forma diferencial, tenemos:
ˆ ~2
R
dW
¯ R~ 1 →R~ 2 = F~ · d~r = q E
~ · d~r ⇔ W1→2 = ~ · d~r
qE
~1
R
Por la definición de potencial eléctrico (ver definición 12 en la página 28), podemos reescribir lo anterior como:
ˆ ~2
R
¯ R~ 1 →R~ 2 = −q
dW ~ · d~r
∇φ
~1
R
ˆ ~2
R ˆ ~2
R h i
dφ ~1
~2 − φ R
¯ R~ 1 →R~ 2 = −q
dW · d~r = −q dφ = −q φ R
~1
R d~r ~1
R
Q.E.D.
Demostración. Por el teorema de Gauss en forma diferencial (ver teorema 4 en la página 24) sabemos:
ρ R ~
~ ~
∇·E R = ~
ε0
~ R
Por la definición de potencial eléctrico (ver definición 12 en la página 28), tenemos E ~ = −∇φ;
~ sustituyendo
llegamos a:
ρ R ~ ρ R ~
∇~ · −∇φ~ = ⇔∇ ~ · ∇φ
~ =−
ε0 ε0
ρ R ~
⇔ ∇2 φ = −
ε0
Q.E.D.
d2 φ d2 φ d2 φ ρ (x, y, z)
2
+ 2 + 2 =−
dx dy dz ε0
~ entonces
Corolario 4 (La ecuación de Laplace). Si la densidad de carga volumétrica es nula en un punto R,
el laplaciano del potencial se anula en ese punto.
~ =0
∇2 φ R
~ por 0.
Demostración. Se sigue trivialmente del teorema 5 sustituyendo ρ R Q.E.D.
d2 φ d2 φ d2 φ
∇2 φ = + 2 + 2 =0
dx2 dy dz
Tenemos garantizado que la solución de la ecuación diferencial (lo probamos más adelante) es única y dicha
solución debe ser de la forma:
d2 X d2 Y d2 Z
YZ + XZ + XY =0
dx2 dy 2 dz 2
Descartando la solución trivial (X (x) = 0, Y (y) = 0, Z (z) = 0), obtenemos:
1 d2 X 1 d2 Y 1 d2 Z 2
+ + =0⇔
X dx2 Y dy 2 Z dz 2
1 d2 Y 1 d2 Z 1 d2 X
⇔ + = −
Y dy 2 Z dz 2 X dx2
Como
1 d2 X
k12 := − 6= F (y, z)
X dx2
obtenemos la ecuación diferencial:
1 d2 X 2 d2 X
− = k 1 ⇔ = −k12 X
X dx2 dx2
Suponiendo k12 6= 0, sabemos que la solución de esta ecuación diferencial es:
X (x) = Ae±ik1 x
Si k1 = 0, la solución es de la forma:
X (x) = Bx + C
En consecuencia, obtenemos:
1 d2 Y 1 d2 Z 2 1 d2 Y 2 1 d2 Z
+ = k1 ⇔ = k1 −
Y dy 2 Z dz 2 Y dy 2 Z dz 2
Ahora, como
1 d2 Z
−k22 := k12 −6= F (y)
Z dz 2
podemos hacer lo mismo que antes, obteniendo la solución:
0 ±ik y
Ae 2 si k2 6= 0
Y (y) 0
B y+C 0 si k2 = 0
Ahora, finalmente, tendríamos:
1 d2 Z
= k12 + k22 =: k32
Z dz 2
cuya solución es:
A00 e±jk2 z si 6 0
k3 =
Z (z)
B 00 z + C 00 si k3 = 0
En conclusión, la solución general es una combinación lineal de senos y consenos trigonométricos e hiper-
bólicos.
Ejemplo 4 (*Ejemplo de resolución de la ecuación diferencial (No entra)). Tenemos la siguiente situación:
φ(b) = 0
b φ(0) = 0
X
Z
Conocemos las siguientes condiciones de contorno:
1. φ = 0 cuando y = 0 o y = b.
3. Cuando x → ∞, φ (x, y, z) = 0.
Si α, β o γ fuesen alguno de ellos nulos, sustituiríamos el término correspondiente por un término del estilo
c01 z + c02 .
Por simetría traslacional de nuestro problema, el potencial no puede depender de z, eso nos implica que
la parte de φ que depende de z debe ser obligatoriamente de la forma:
c01 = 0
Y, en consecuencia:
Z (z) = c02
X h i
φ (x, y, z) = a1 (α) eαx + a2 (α) e−αx b1 (β) eβy + b2 (β) e−βy c02
α,β
Sabemos que debe cumplirse (por las condiciones de las k que hallamos antes en la deducción de la solución
general):
α2 + β 2 = 0
Suponiendo α ∈ R:
−α2 = β 2 ⇔ iα = β
Y, por tanto:
X
φ (x, y, z) = a01 (α) eαx + a02 (α) e−αx b01 (iα) eiαy + b02 (iα) e−iαy
α
Suponiendo α > 0; por la condición de que el potencial en el infinito debe ser cero, para que la función no
«explote», obligatoriamente debe ser:
a01 (α) = 0
Luego la solución es necesariamente de la forma:
X
φ (x, y, z) = a02 (α) e−αx b01 (iα) eiαy + b02 (iα) e−iαy
α
Ahora, de nuevo, podemos hacer el mismo truco de antes y eliminarnos la constante a02 (α), metiéndola por
la propiedad distributiva en el corchete, cambiando las constantes b01 y b02 . Nos queda:
X
φ (x, y, z) = e−αx b001 (iα) eiαy + b002 (iα) e−iαy
α
φ (x, 0) = 0
Como la exponencial nunca se anula debe ser:
Sabiendo:
eiθ − e−iθ
sen θ =
2i
Llegamos a:
X
φ (x, y, z) = b001 (iα) e−αx 2i sen αy
α
Ahora aplicamos otra condición de contorno: φ (x, b) = 0. Para ello debe ser:
nπ
sen αb = 0 ⇔ αb = nπ ⇔ α =
b
Así, nuestra solución queda:
∞
X nπ
− nπ x
φ (x, y, z) = Bn e β sen y
b
n=0
donde:
ˆ b nπ
2
Bn = φ (0, y) sen y dy
b 0 b
4V0
Bn =
nπ
A continuación, definimos V0 := φ (0, y). Con esto, podemos escribir:
∞
4V0 X 1 (2n − 1) π (2n−1)π
φ (x, y) = sen y e− b x
π 2n − 1 b
n=1
Ez = 0
1. ‹ ˚
∂ψ ~ · ∇φ
~
φ dS 0 = φ · ∇2 ψ + ∇ψ dV 0
S 0 ∂ n̂ V0
2. ‹ ˚
∂φ ∂ψ
ψ −φ dS 0 = ψ∇2 φ − φ∇2 ψ dV 0
S0 ∂ n̂ ∂ n̂ V0
˚ ‹ ‹
~ φ∇ψ
~ ~ · dS
~0 = ∂ψ
∇ dV 0 = φ∇ψ φ dS (2.5.0.2)
V 0 S 0 S 0 ∂ n̂
~
donde n es el vector unitario que tiene como dirección la dirección de dS. 1
~n
(0, dS, 0) ~i2
~i1
Y
el potencial en
la superficie
(capa dipolar)
=−E·n̂
~
z}|{ z }| {
˚ ρ ~
R 0 dV 0 ‹ 1
1 1 1 ∂φ ∂ r 0
~ 0
φ R = + − φ dS
4πε0 V0 r 4π S 0
r ∂ n̂ ∂ n̂
| {z }
potencial de las cargas
dentro del volumen
| {z }
lo que sucede fuera del volumen
S0
~
ρ(R)
ρ=0
V
Tomamos:
1
ψ :=
r
y, ahora, por la segunda identidad de Green (ver teorema 6 en la página 37), tenemos:
˚ ~0 ‹ !
1ρ R
1
~ ∇
0 1 0 1 ∂φ ∂
− −φ R 2
dV = −φ r
dS 0
V 0 r ε 0 r S 0 r ∂ n̂ ∂ n̂
| {z }
=−4πδ(~
r)
~ −R
Por la proposición 11 en la página 22 y recordando que ~r = R ~ 0 , tenemos que lo anterior es equivalente a:
˚ ρ R ~ 0 dV 0 ˚ ‹ 1
!
1 ~ 4πδ R
0 ~ −R ~ dV =
0 0 1 ∂φ ∂
− + φ R −φ r
dS 0
ε0 r V 0 S 0 r ∂ n̂ ∂ n̂
Por la propiedad de traslación de la delta (ver proposición 10 en la página 21), lo anterior es equivalente a:
˚ ρ R ~ 0 dV 0 ‹ 1
!
1 ~ = 1 ∂φ ∂
− + 4πφ R −φ r
dS 0 ⇔
ε0 r S 0 r ∂ n̂ ∂ n̂
˚ ~ 0 dV 0 ‹ !
1 ρ R 1 1 ∂φ ∂ 1
~ =
φ R + −φ r dS 0
4πε0 V0 r 4π S0 r ∂ n̂ ∂ n̂
Q.E.D.
Lema 1. Sean Ω un abierto en R3 y ψ : R3 −→ R 3 ψ ∈ C 2 R3 , R , entonces se cumple:
2
~ · ψ ∇ψ
∇ ~ ~ · ∇ψ
− ψ∇2 ψ = ∇ψ ~ = ∇ψ
~
*Demostración (No entra). Por la definición de divergencia (ver definición 6 en la página 16) y por la definición
de gradiente (ver definición 5 en la página 16), tenemos:
∂ ∂ ∂ ∂ψ ∂ψ ∂ψ ∂
∂ψ
∂
∂ψ
∂
∂ψ
~ ~
∇ · ψ ∇ψ = , , · ψ ,ψ ,ψ = ψ + ψ + ψ =
∂x ∂y ∂z ∂x ∂y ∂z ∂x ∂x ∂y ∂y ∂z ∂z
∂ψ ∂ψ ∂ 2 ψ ∂ψ ∂ψ ∂ 2 ψ ∂ψ ∂ψ ∂2ψ
= +ψ 2 + +ψ 2 + +ψ 2 =
∂x ∂x ∂x ∂y ∂y ∂y ∂z ∂z ∂z
2 2 2
∂ψ ∂2ψ ∂ψ ∂2ψ ∂ψ ∂2ψ
= +ψ 2 + +ψ 2 + +ψ 2 =
∂x ∂x ∂y ∂y ∂z ∂z
2 2 2 2
∂ψ ∂ψ ∂ψ ∂ ψ ∂2ψ ∂2ψ
= + + +ψ + + =
∂x ∂y ∂z ∂x2 ∂y 2 ∂z 2
∂ψ ∂ψ ∂ψ ∂ψ ∂ψ ∂ψ ~ · ∇ψ
~ + ψ∇2 ψ
= , , · , , + ψ∇2 ψ = ∇ψ
∂x ∂y ∂z ∂x ∂y ∂z
Entonces, claramente:
2
~ ~ ~ · ∇ψ
∇ · ψ ∇ψ − ψ∇2 ψ = ∇ψ ~ + ψ∇2 ψ − ψ∇2 ψ = ∇ψ
~ · ∇ψ
~ = ∇ψ
~
Q.E.D.
1. (Condición de Dirichlet) es única si conocemos el potencial en la superficie que limita el volumen estu-
diado y dicho potencial es constante a lo largo de ella.
2. (Condición de von Neumann) no es única, pero todas las soluciones difieren únicamente en una constante
si conocemos la componente perpendicular a la superficie que limita el volumen estudiado en dicha
superficie al completo.
Demostración. Tenemos:
~
ρ(R)
ρ=0
S0
c0
∇2 ψ = ∇2 φ1 − ∇ 2 φ2 = 0
Justificación de Dirichlet:
Supongamos que conocemos el potencial en la superficie S y que es constante en ella, entonces:
‹
φ1 (~x) = φ2 (~x) ∀~x ∈ S ⇒ ψ (~x) = φ1 (~x) − φ2 (~x) = 0 ∀~x ∈ S ⇒ ~ · dS
ψ ∇ψ ~=0
S
Por tanto, por el teorema de la divergencia (ver teorema 2 en la página 17):
‹ ˚
0= ~ ~
ψ ∇ψ · dS = ∇~ · ψ ∇ψ
~ dV
S V
2
~
Como era ψ = 0, debe ser necesariamente ∇ψ = 0. Ahora bien:
2
~
∇ψ ~ = ~0 en todo V ⇒ ∇φ
= 0 en todo V ⇒ ∇ψ ~ 1 = ∇φ
~ 2 en todo V ⇒ E
~1 = E
~ 2 en todo V
Por otra parte, también por el teorema de la divergencia (ver teorema 2 en la página 17):
‹ ˚
~
∇ψ = 0 ⇒ 0 = ~ ~
ψ ∇ψ · dS = ∇~ · ψ ∇ψ
~ dV
S V
2
~
Como era ∇2 ψ = 0, debe ser necesariamente ∇ψ = 0:
2
~
∇ψ ~ = ~0 en todo V ⇒ ∇φ
= 0 en todo V ⇒ ∇ψ ~ 1 = ∇φ
~ 2 en todo V ⇒ E
~1 = E
~ 2 en todo V
Q.E.D.
Ejemplo 5. Imaginemos:
−
r r
q −q q
−
φ=0 φ=0
A B
Ambos problemas A y B son equivalentes. La función φ es la misma en ambos problemas.
Capítulo 3
Distribuciones de dipolos
Un dipolo es lo siguiente:
+q
−q ~δ
Dipolo
Un cuadrupolo es lo siguiente:
−q
+q
+q
−q
Cuadrupolo
Nótese que ambos se caracterizan porque la carga total es nula. Análogamente, un octupolo es:
+q
−q
−q −q
+q
+q +q
−q
Octupolo
Y, así, sucesivamente.
+q
~δ
Z (x0 , y 0 , z 0 )
−q ~r
(x, y, z)
~0
R
~
R
Y
X
p~ := q~δ
donde ~δ es el vector que va de la carga negativa a la positiva.
1 p~ · r̂
~ ≈
φdipolo R
4πε0 r2
Demostración. El potencial generado por el dipolo (dos cargas puntuales de sentido contrario) en un punto
~ es:
R
~ = φR~ 0 +~δ + φR~ 0 = φR~ 0 +~δ − φR~ 0
φdipolo R +q −q +q +q
donde con los superíndices denotamos la posición y con los subíndices indicamos la carga. Es decir, el
potencial que genera el dipolo en el punto R ~ es el que que genera una carga −q en el punto R ~ 0 más el que
~ ~
genera una carga q en el punto R + δ, pero esto último podemos reescribirlo como el potencial que genera
0
p~
Z (x0 , y 0 , z 0 )
~r
(x, y, z)
~0
R
~
R
Y
X
En consecuencia, nuestro problema inicial es equivalente a estudiar lo que cambiaría el potencial si moviéramos
una carga +q del punto R~ 0 al punto R ~ 0 + ~δ.
Por otra parte, ~δ |~r| ⇒ R~ 0 + ~δ ≈ R ~ 0 + lı́m ~δ. Si ~δ → ~0, entonces ∆φ R
~ → ~0. Es decir, como ~δ es muy
~
δ→~0
~ 0 , podemos aproximar la variación de potencial por un diferencial de potencial:
pequeño respecto a R
∆φ R ~ ≈ dφ R ~
q R ~ − − ~δ ~0
R ~−
R q −~δ
~0
R
2 − 2 ≈ 2
4πε0 ~ 0 − ~δ
~ −R ~ −R
~0 4πε0 ~ ~ 0
R R R−R
Esto nos va a permitir usar conceptos ya conocidos del cálculo. Ahora bien, si recordamos lo que era el
diferencial de una función1 , vemos que claramente:
~ = ∂φ dx0 + ∂φ dy 0 + ∂φ dz 0 = ∇
dφ R ~ · ~δ
~ 0φ R (3.2.0.1)
∂x0 ∂y 0 ∂z 0
donde ~δ = (dx0 , dy 0 , dz 0 ).
Ahora, por la proposición 8 en la página 19:
~ 1 ~ 0 1 q ~ 1 q ~0 1 ~ q ~ 0 q
∇ = −∇ ⇒ ∇ =− ∇ ⇔∇ = −∇
r r 4πε0 r 4πε0 r 4πε0 r 4πε0 r
p · ~r) ~r − r2 p~
1 3 (~
~ R
E ~ ≈
4πε0 r5
n
1 f: Ω −→ R X ∂f
Sean Ω abierto en Rn y . Entonces se define df como df := dxi .
(x1 , . . . xn ) f (x1 , . . . , xn ) ∂xi
i=1
Demostración. Simplemente tenemos que aplicar la definición de potencial eléctrico (ver definición 12 en la
página 28):
~ R
E ~ = −∇φ ~ ~
R
Por la proposición 18 en la página 46, sabemos:
X
pj j − j 0
1 p~ · ~r 1 j=x,y,z
~ ≈
φdipolo R =
4πε0 r 3 4πε0 3
2
0 2
0 2
0 2
x−x + y−y + z−z
| {z }
=r2
Derivando, obtenemos:
X X
r 3 px − pj (j − j 0 ) 32 r · 2 (x − x0 ) · 1 px r 2 − 3 pj (j − j 0 ) (x − x0 )
∂φ 1 j=x,y,z
1 j=x,y,z
=
= 6 = 5
∂x 4πε0 r 4πε0 r
p · ~r (x − x0 )
1 px r2 − 3~
=
4πε0 r5
Análogamente:
∂φ 1 p · ~r (y − y 0 )
py r2 − 3~
=
∂y 4πε0 r5
∂φ 1 p · ~r (z − z 0 )
pz r2 − 3~
=
∂z 4πε0 r5
Ahora:
~ R
E ~ = −∇φ
~ ~ =
R
1 1
= 5
p · ~r x − x0 , −py r2 + 3~
−px r2 + 3~ p · ~r y − y 0 , −pz r2 + 3~
p · ~r z − z 0 =
4πε0 r
1 1 2
= 5
p · ~r x − x0 , y − y 0 , z − z 0 =
−r (px , py , pz ) +3~
4πε0 r | {z } | {z }
=~
p =~
r
p · ~r) ~r −
1 3 (~ r2 p~
=
4πε0 r5
Q.E.D.
~ = α~r + β~
E p
donde α, β ∈ R. Es decir, el vector campo eléctrico está contenido en el plano formado por los vectores ~r
y p~.
Z
Er
θ
Eθ
~r
1 pr cos θ 1 p cos θ
φ (r, θ) = 3
=
4πε0 r 4πε0 r2
Si recordamos la expresión del gradiente en coordenadas esféricas, obtenemos:
=0
z}|{
~ ∂φ 1 ∂φ 1 ∂φ
E=− r̂ − θ̂ − ϕ̂
∂r} | r{z∂θ}
| {z r sen θ ∂ϕ
=:Er =:Eθ
∂φ 1 p cos θ
Er = − =
∂r 2πε0 r3
1 ∂φ 1 1 p sen θ 1 p sen θ
Eθ = − = 2
=
r ∂θ r 4πε0 r 4πε0 r3
Por consiguiente:
~ = 1 2p cos θ p sen θ p
E 3
r̂ + θ̂ = 2 cos θr̂ + sen θ θ̂
4πε0 r r3 4πε0 r3
Ahora, podemos hallar las llamadas posiciones de Gauss:
1 p
θ = 0 ⇒ Er = y Eθ = 0
2πε0 r3
π 1 p
θ= ⇒ Er = 0 y Eθ =
2 4πε0 r3
V0
n̂
p~ S0
dV 0
Z ~r
~0
R
~
R
Y
X
Definición 14. Definimos el momento dipolar total como la suma de los momentos dipolares de todos los
dipolos involucrados:
n
X
p~ := p~i
i=1
Observación 15. Como es usual en física, si el número de partículas es muy alto, podemos aproximar la suma
finita por una integral, de manera que:
˚
p~ ≈ d~
p
V0
Definición 16. Llamamos densidad de carga equivalente volumétrica o densidad de Poisson volu-
métrica a:
~ 0 · P~ R
ρb := −∇ ~0
~ = (A1 , . . . , An ):
*Demostración (No entra). Sea A
∂ ∂
n
X ∂
~ 0 ~
∇ · ϕA = , . . . , · (ϕA 1 , . . . , ϕA n ) = (ϕAi ) =
∂x10 ∂xn0 ∂x0i
i=1
n
X X n X n
∂ϕ ∂Ai ∂ ∂ϕ
= 0 A i + ϕ 0 = ϕ 0 (A i ) + Ai · 0 =
∂xi ∂xi ∂xi ∂xi
i=1 i=1 i=1
∂ ∂ ∂ϕ ∂ϕ
=ϕ , . . . , 0 · (A1 , . . . , An ) + (A1 , . . . , An ) · ,..., 0 =
∂x01 ∂xn | {z } | {z } ∂x01 ∂xn
| {z } =A~ =A ~ | {z }
~0
=∇ ~ 0ϕ
=∇
~0·A
= ϕ∇ ~+A
~·∇
~ 0ϕ
Q.E.D.
Proposición 20. Sea una distribución de dipolos que ocupa un volumen V 0 y sea S 0 la superficie cerrada que
engloba dicho volumen, entonces el potencial y el campo eléctrico generado por dicha distribución de dipolos
pueden expresarse como:
1. ‹ ˚
1 σb 0 ρb 0
~ =
φ R dS + dV
4πε0 S0 r V0 r
2. ‹ ˚
1 σb~r 0 ρb~r 0
~ R
E ~ = dS + dV
4πε0 S0 r3 V0 r3
1 ~
Por otra parte, haciendo uso del lema 2 en la página anterior y tomando φ = y A = P~ R
~ , llegamos a:
r
P~ R~0 ~ 0 P~ R
∇ ~0
~ 0 ~ ~ 0 ~ 0 1
∇ = +P R ·∇ ⇔
r r r
P~ R~0 ~ 0 P~ R
∇ ~0
~ 0 − ~ ~ 0 ~ 0 1
⇔∇ =P R ·∇
r r r
pues la divergencia y la integral conmutan al ser la integral a las coordenadas primadas y la divergencia a las
coordenadas sin primar. Por la proposición 7 en la página 18, la expresión anterior queda:
‹ ˚
~ ~ 1 σb~r 0 ρb~r 0
E R = 3
dS + 3
dV
4πε0 S0 r V0 r
Q.E.D.
Ejemplo 6. Tenemos:
p~
− + − +− +
− +
− +
− +
−x̂ − + x̂
− +
En el interior los dipolos se anulan y no hay ρb , pero lo habría si el volumen fuese irregular (con, por
ejemplo, huecos («burbujas») en su interior).
La carga ligada sería:
˚ ‹
Qb = ~ 0 ~ ~0
−∇ P R dV + 0
P~ dS
~0
V0 S0
Por el teorema de la divergencia (ver teorema 2 en la página 17), lo anterior es equivalente a:
‹ ‹
Qb = − ~ ~
P dS +0
P~ dS
~0 = 0
S0 S0
n̂ = ẑ
σb = pb
σb = −pb
P~ = P0 ẑ
~ · P~ = 0 + 0 + ∂Pz
ρb = −∇ =0
∂z
|{z}
=0
tapa superior: P~ · n̂ = P0
σb = superficie lateral: P0 ẑ · r̂ = 0
tapa inferior: P0 ẑ · (−ẑ) = −P0
Hemos reducido al vacío; es decir: resulta que tener un cilindro polarizado como éste es equivalente a tener
dos discos.
Capítulo 4
~r
~0
R
~
R
1 1
Nuestro objetivo va a ser reemplazar por , donde R representa la distancia a nuestro origen de
r R
coordenadas esféricas. Es decir conseguir:
1 1
u + ...
r R
φ R~ = 1 q
4πε0 r
donde:
r a 2
r= (x − 0)2 + (y − 0)2 + z −
2
Z
(x, y, z)
~r
(0, 0, a2 )
~
R
Y
X
1 dn h 2 2 i
Pn (x) := x − 1
2n n! dxn
Ejemplo 8 (La carga puntual). Imaginemos:
~r
q
θ ~
R r
R0
Y
X
Claramente:
~ 0 = a · ẑ = (0, 0, a)
R
1 q
~ =
φ R
4πε0 r
Aplicando el teorema del coseno, obtenemos:
r2 = a2 + R2 − 2aR cos θ
En consecuencia:
1 q
~ =
φ R √
4πε0 R + a − 2aR cos θ
2 2
a 2 a − 1
1 2
1+ −2 cos θ =
r r r
" 2 #
1 1 a 2 a 3 a 2 a
= 1− − 2 cos θ + − cos θ + . . . =
r 2 r r 8 r r
a 2 a 3
1 1 a 2 a 3 a 4 2
= 1− − 2 cos θ + +4 cos θ − 2 cos θ + . . . =
r 2 r r 8 r r r
1 a 3 2 1 a 2
= 1 + cos θ + cos θ − + ...
r r 2 2 r
De esta forma:
q a 3 cos2 θ − 1 a 2
φ (r) = 1 + cos θ + + ...
4πε0 r r 2 r
i
Podemos ver que justo las funciones que acompañan a ar son justo Pi (cos θ) donde con Pi denotamos los
polinomios de Legendre (ver la definición 17 en la página 56). Consecuentemente:
a 2 a 3
q a
φ (r) = P0 (cos θ) + P1 (cos θ) + P2 (cos θ) + P3 (cos θ) + . . . =
4πε0 r r r
1 q X a l
∞
= Pl (cos θ)
4πε0 r r
l=0
Observación 16. Por la expresión anterior, podemos ver que tener una carga a una distancia a es equivalente
a tener una carga puntual en el centro de coordenadas más un dipolo, más un cuadrupolo, etc.
Observación 17. Podemos expresar la densidad de carga cuando sólo tenemos una carga puntual en el origen
como:
ρ R~ = δ (x) δ (y) δ (z) q
| {z }
=δ(x,y,z)
Como ya sabemos:
˚ ˚
QT = ~
ρ R dV = qδ (x, y, z) dV
V V
por la propiedad de traslación de la delta (proposición 10 en la página 21) la integral anterior valdrá q si
la carga está contenida en V y 0 en caso contrario.
qn zn
q3 z3
q2 z2
q1 z1
Z
(r, θ, φ)
dz
θ ~r
Y
X
Anteriormente, hemos obtenido:
1 q X a l
∞
φ (r, θ) = Pl (cos θ)
4πε0 r r
l=0
Si la carga está en el eje Z, podemos escribir lo anterior como:
∞
1 q X z0 l
φ (r, θ) = Pl (cos θ)
4πε0 r r
l=0
Si ahora tenemos varias cargas alineadas en el eje Z, debe ser:
∞ 0 l
1 λ (z 0 ) X z
dφ (r, θ) = Pl (cos θ) dz 0
4πε0 r r
l=0
En consecuencia:
ˆ ∞ 0 l
∞
1 λ (z 0 ) X z
φ (r, θ) = Pl (cos θ) dz 0 =
−∞ 4πε0 r r
l=0
∞ ˆ ∞
X
1 λ (z 0 ) · (z 0 )l Pl (cos θ) 0
= dz =
4πε0 −∞ rl+1
l=0
∞ ˆ
1 X Pl (cos θ) ∞ 0
0
0 l
= λ z · z dz =
4πε0 rl+1
l=0 | −∞ {z }
momento dipolar
axial de orden l
=: Ml
∞
1 X Ml Pl (cos θ)
=
4πε0 rl+1
l=0
Por tanto:
∞
1 X Ml Pl (cos θ)
φ (r, θ) = (4.4.0.1)
4πε0 rl+1
l=0
q1
a
2
q2
a Y
2
X
q3
Ejemplo 12 (el dipolo). Usando el ejemplo anterior, tomemos q = q1 = −q3 y q2 = 0. En este caso:
a 3
M0 = 0 M1 = q · a M2 = 0 M3 = 2a 2 M4 = 0
Para un l cualquiera, Ml únicamente será distinto de cero para los l impares. Obtenemos, al sustituir
en 4.4.0.1, el siguiente potencial:
" #
1 q · a cos θ q a3 5 − cos3 θ − 3 cos θ
φ (r, θ) = 0+ +0+ 4 + ...
4πε0 r2 r 8
Ejemplo 13 (el cuadrupolo). Usando el ejemplo 11 en la página anterior con los valores q = q3 = q y
q2 = −2q, llegamos a:
2
M0 = 0 M1 = 0 M2 = 2q a4 M3 = 0
En general:
(
Ml = 0 si l es impar
a l
Ml = Ml = 2q 2 si l es par
De esta forma, sustituyendo en 4.4.0.1 en la página anterior, llegamos a:
1 qa2 P2 (cos θ) qa4 P4 (cos θ)
φ (r, θ) = 0+0+ +0+ + ...
4πε0 2 r3 8 r5
Demostración. Tenemos:
dV 0
r
R0
R
Sabemos:
˚ ~0
ρ R
~ = 1
φ R dV 0
4πε0 V 0 r
~r = x − x0 , y − y 0 , z − z 0
q
r= (x − x0 )2 + (y − y 0 )2 + (z − z 0 )2
Vamos a hacer el desarrollo en serie de Taylor en torno al punto que cumple R0 = 0:
3
" # 3 X 3
" #
1 1 X ∂ 1r 1 X ∂ 2 1
= + x0i + r
x0i x0j + · · · =
r r R0 =0 ∂x0i 2 ∂x 0 ∂x0
i j
i=1 R0 =0 i=1 j=1 R0 =0
1 1 ~0 1 ~0 + 1R ~0 H 1 ~0 + . . .
= = + ∇ ·R R
r r R0 =0 r R0 =0 2 r R0 =0
donde H es la matriz hessiana.
∂ 1r ∂ 1
Recordando = − r0 ∀i = 1, 2, 3, podemos expresar lo anterior como:
∂xi ∂xi
3
" # 3 X
3
" #
1 1 X ∂ 1r 1 X ∂ 2 1
= − x0i + r
x0i x0j + . . .
r r R0 =0 ∂xi 0
2 ∂x i ∂x j 0
i=1 R =0 i=1 j=1 R =0
Ahora:
!
1
∂ 1
∂ 1r [(x−x0 )2 +(y−y0 )2 +(z−z 0 )2 ] 2 − (xi − x0i )
= =
∂xi ∂xi r3
∂ 2 1r −r3 + (xi − x0i ) 3 21 r (xi − x0i ) 2 −r2 + 3 (xi − x0i )2
= =
∂x2i r6 r5
0 ) 3r 2 1 1 2 x − x0 0 ) x − x0
2
∂ r 1 (x i − x i 2r j j 3 (xi − x i j j
= 6
= 5
∂xi ∂xj r r
En consecuencia:
" #
∂ 1r − (xi − x0i ) −xi
= 3
= 3
∂xi r R0 =0 R
R0 =0
" # " #
∂ 2 1r −r2 + 3 (xi − x0i )2 −R2 + 3x2i
= =
∂x2i r5 R5
R0 =0 R0 =0
"
# 0 ) x − x0
2
∂ r 1 3 (xi − x i j j 3xi xj
= 5
=
∂xi ∂xj 0 r R5
R =0 R0 =0
Podemos compactar las dos expresiones anteriores en:
" #
∂ 2 1r 3xi xj − R2 δij
=
∂xi ∂xj 0 R5
R =0
~0 1 X
~ ·R 3 X
3
1 R 3xi xj − δij R2 0 0
= + + xi xj + . . .
R R3 2 R5
i=1 j=1
~
= p~·R
3
R
˚ z }| {
~ ˚ X3 X 3 2 ˚
1 1 R 1 3x x − δ R
= ρ R~ 0 dV 0 + · ~ 0 ρ R0 dV 0 +
R
i j ij
x 0 0
i xj ρ ~ 0 dV 0 + . . .
R
4πε0 R 0 R3 2 R5 0
| V {z } | {z } i=1 j=1 | V {z }
=q =~
p =:Q0ij
˝ 0 0 0
donde p~ = ~ R
ρ R ~ dV .
Es decir:
1 q ~
p~ · R
3 3
1 X X 3xi xj − δij R2 0
~ =
φ R + 3 + Qij + . . .
4πε0 R R 2 R5
i=1 j=1
Q.E.D.
4.7. Ejemplos
Ejemplo 14 (Problema 27). Tenemos un anillo:
~
R
R0 Y
X
dl
El valor del potencial en el eje Z es muy fácil de calcular, de hecho, lo podemos hallar de forma analítica:
q
λ=
2πR0
1 q
φ (0, 0, z) = √
4πε0 R + z 2
2
Sin embargo, fuera de el eje es difícil, por lo que usamos el desarrollo multipolar. Ya tenemos q; ahora,
calculemos p~:
ˆ ˆ 2π
px = λx0 dl0 = λ R0 cos ϕ0 R0 dϕ0 = 0
C0 0
ˆ
py = λy 0 dl0 = 0
C0
ˆ
pz = λz 0 dl0 = 0
C0
Es decir:
p~ = ~0
Podemos ver esto gráficamente mediante:
~
dl
d~ ~ 0 dl0
p = λR
~
dl
d~ ~ 0 )dl0
p = λ(−R
Debido a la simetría respecto al eje Z no hay momento dipolar. Ahora, toca calcular los Q0xy .
ˆ ˆ 2π
Q0xy = Q0yx = 0 0
x y λdl = λ 0
R0 cos ϕ0 R0 sen ϕ0 R0 dϕ0 =
C0 0
ˆ 2π
λR03 2 2π
= λR03 sen ϕ0 cos0 dϕ0 = sen ϕ 0 = 0
0 2
Q0zz = 0
ˆ ˆ ˆ 2π
Q0xx = 02
x λdl = λ 0
R02 cos2 ϕ0 R0 dϕ0 = λR03 cos2 ϕ0 dϕ0
C0 0
´ ϕ sen 2ϕ
Sabiendo cos2 ϕdϕ = 2 + 4 , obtenemos:
2π
ϕ sen 2ϕ 2π
= λR03 + = λR03 + 0 − 0 − 0 = λπR03
2 4 0 2
Nótese que:
Q0yy = Q0xx
por simetría. Por último:
Sustituyendo:
~ = q R2
φ R 1 + 04 3 x2 + y 2 − 2R2 + . . .
4πε0 R 4R
Teniendo en cuenta R2 = x2 + y 2 + z 2 :
~ = q R2 x2 + y 2 − 2z 2
φ R 1 + 04 + ...
4πε0 R R 4
Ejemplo 15 (R = 10R0 en el eje z). En las mismas condiciones que en el ejemplo anterior, vamos a calcular
φ (0, 0, 10R0 ). Para este caso, también tenemos una solución exacta; de manera que podremos comparar.
!
1 Q R02 −2 (10R0 )2 1 Q 1
φ≈ (0, 0, 10R0 ) ≈ 1+ = 1−
4πε0 10R0 (10R0 )4 4 4πε0 10R0 200
Podemos ver que:
φcuadrupolar 1
= = 0, 5 %
φmonopolar 200
Por otra parte, el potencial exacto es:
1 Q
φ= (z = 10R0 ) = √
4πε0 101R0
√
φmonopolar 101
= = 1, 00498 ⇒ Error relativo del 0, 5 %
φ= 10
φ≈
= 0, 9999626 ⇒ Error relativo de 3, 74 · 10−5
φ=
Ejemplo 16 (z = 0, R = 10R0 ). Nos encontramos en las mismas condiciones que en el ejemplo 14 en la
página 64 y queremos calcular el potencial para las condiciones z = 0, R = 10R0 . Entonces:
R2 = x2 + y 2 = (10R0 )2
El término aproximado es:
!
1 Q R02 (10R0 )2 1 Q 1
φ≈ (x, y, 0) ≈ 1+ = 1+
4πε0 10R0 (10R0 )4 4 4πε0 10R0 400
Podemos calcular el cociente:
φcuadrupolar 1
= = 0, 25 %
φmonopolar 400
Ejemplo 17 (Problema 26). Tenemos la siguiente distribución de cargas puntuales (un octupolo).
−q +q
+q
−q
+q
−q
−q +q
qT = 0
Ahora, obtengamos la densidad de carga expresada mediante la delta de Dirac.
ρ R0 = qδ x0 , y 0 , z 0 − qδ x0 − L, y 0 , z 0 + qδ x0 − L, y 0 − L, z 0 − q x0 , y 0 − L, z 0 +
qδ x0 , y 0 − L, z 0 − L − qδ x0 , y 0 , z 0 − L + qδ x0 − L, y 0 , z 0 − L − qδ x0 − L, y 0 − L, z 0 − L
Lo siguiente es calcular el momento dipolar.
˚
p~ = ρ R ~ 0 dV 0 = q (0, 0, 0) − q (L, 0, 0) + q (L, L, 0) − q (0, L, 0) + q (0, L, L) − q (0, 0, L) +
~0 · R
V0
Qyy = 0 = Qzz
φ≈ ≈ 0
˚ ˚ ˚
3x2 − R2 ~ 3x22 − R2 ~ 3x23 − R2 ~ 0 dV 0 =
= 1 5 x0 21 ρ 0
R dV +0
x 02
2 ρ R 0
dV 0
+ x02
3 ρ R
R V0 R5 V0 R5 V0
˚
3x21 − R2 x0 21 + 3x22 − R2 x0 22 + 3x23 − R2 x0 23 ~ 0 0
= ρ R dV =
V0 R5
˚
3x21 x0 21 + 3x22 x0 22 + 3x23 x0 23 − R2 x0 21 + x0 22 + x0 23 ~ 0 0
= ρ R dV =
V0 R5
=x21 +x22 +x23
˚ z}|{
3x21 x0 21 + 3x22 x0 22 + 3x23 x0 23 − R2 ·R0 2 ~ 0 0
= ρ R dV =
V0 R5
˚
3x0 21 − R0 2 x21 + 3x0 22 − R0 2 x22 + 3x0 23 − R0 2 x23 ~ 0 0
= ρ R dV =
V0 R5
˚ ˚ ˚
x2 0 0 x22 0 x23 0 0
= 15 x0 21 02~
− R ρ R dV + 5 02 02 ~
3x 2 − R ρ R dV + 5 ~ dV
3x0 23 − R0 2 ρ R
R V0 R V0 R V0
Ahora, si definimos:
˚
Qij := ~ 0 3x0 x0 − R02 δij dV 0
ρ R i j
V0
La expresión anterior quedaría:
3
X X x2 3
3xi xi − δii R2 x21 x22 x23
Q0ii = Q11 + Q22 + Q33 = i
Qii
R5 R5 R5 R5 R5
i=1 i=1
˚
3xi xj − δij R2 0 3xi xj 0 3xi xj 0 0 ~ 0 dV 0 =
Qij = Q ij = x x
i j ρ R
R5 i6=j R 5 R 5
V 0
˚
xi xj ~ 0 dV 0 = xi xj [Qij ]
= 5 3x0i x0j ρ R i6=j
R V0 R5
Es decir, podemos reescribir la suma:
3 X
X 3 3
X 3
X
3xi xj − δij R2 0 3x2i − R2 0 3xi xj 0
Qij = Qii + Qij =
R5 R5 R5
i=1 j=1 i=1
i, j = 1
i 6= j
3
X 3
X 3
X
x2i xi xj xi xj
= Qii + Qij = Qij
R5 R 5 R5
i=1 i,j=1
i, j = 1
i 6= j
En consecuencia, obtenemos una expresión equivalente para el polinomio de Taylor:
~ X3 X 3
φ R~ = 1 q + p~ · R + 1 xi xj
Qij + . . .
4πε0 R R 3 2 R5
i=1 j=1
Q.E.D.
2. Si por el origen pasa un plano de simetría, entonces p~ estará contenido en dicho plano.
3. Si por el origen pasa un eje de simetría, entonces p~ estará contenido en dicho eje.
Demostración.
~ = ρ −R
1. Si el origen de cargas es un centro de simetría, entonces ρ R ~ ∀R~ ∈ V . Recordemos que,
por la proposición 22 en la página 69 es:
˚
p~ = ρ R ~
~ RdV
V
Ahora, podemos dividir la integral al volumen en dos integrales, lo que hay del plano z = 0 hacia arriba
(a lo que llamaremos V +) y lo que hay del plano z = 0 hacia abajo (lo que llamaremos V −). Así:
˚ ˚ ˚ ˚ h i
p~ = ~ ~
ρ R RdV = ~ ~
ρ R RdV + ~ ~
ρ R RdV = ρ R~ R ~ + ρ −R ~ −R ~ dV
V V+ V− V+
~ =ρ R
Como ρ −R ~ , podemos expresar la integral anterior como:
˚
p~ = ~ R
ρ R ~ −R
~ dV = ~0
V+
2. Podemos suponer, sin pérdida de generalidad, que el plano de simetría es el XY . Por tanto, sabemos
que: ρ (x, y, z) = ρ (x, y, −z) ∀ (x, y, z) ∈ V . Si la componente pz es nula, entonces p~ necesariamente
estará en el plano XY . Veámoslo, por la proposición 22 en la página 69, tenemos:
˚
pz = ρ (x, y, z) zdV
V
De nuevo, podemos dividir la integral al volumen en dos integrales, lo que hay del plano z = 0 hacia
arriba (a lo que llamaremos V +) y lo que hay del plano z = 0 hacia abajo (lo que llamaremos V −). Así:
˚ ˚ ˚
pz = ρ (x, y, z) zdV + ρ (x, y, z) dV = [ρ (x, y, z) z + ρ (x, y, −z) (−z)] dV
V+ V− V+
3. Podemos suponer, sin pérdida de generalidad, que el eje de simetría es el eje z. Si z es eje de simetría,
entonces los planos XZ e Y Z también serán de simetría. En consecuencia, por (2) será px = 0 = py y,
por ende, p~ únicamente tiene componente z.
Q.E.D.
porque, si el origen es el centro de simetría, entonces, por cada d~ p que ejerza un dV de la distribución
debe existir otro dV que ejerza menos −d~ p. Si la distribución no es uniforme, de manera que el origen no es
centro de simetría, podemos ver claramente que p~ 6= ~0 con el siguiente ejemplo:
ρ0 cos θ
(r, θ, φ0 )
θ0 θ
(r, θ + π, φ0 )
1
Qxx = Qyy = − Qzz
2
Qzz 3 cos2 θ − 1
φcuadrupolar =
4πε0 4R3
Qxz = 0 = Qyz
(Qij ) = (0)
˚ p
= ρ (r, z) 3r2 cos2 θ − r2 + z 2 dV =
V
ˆ ∞ ˆ ∞ ˆ 2π p
= ρ (r, z) 3r2 cos2 θ − r2 + z 2 rdrdθdz
r=0 z=−∞ θ=0
Ahora, como la densidad de carga ρ no depende de θ, sale fuerza de la integral. En consecuencia, llegamos
a: ˆ ∞ˆ ∞ ˆ 2π
p
Qxx = ρ (r, z) 3r 2 2
cos θdθ − 2π r2 + z 2 rdrdz (4.9.0.1)
r=0 z=−∞ θ=0
Resolvamos la integral del coseno aparte:
ˆ 2π
cos2 θdθ (4.9.0.2)
0
1 1
sen π − 0 + sen 0 = π π−
2 2
En consecuencia, podemos sustituir en la ecuación 4.9.0.1 la integral del coseno por la integral del seno,
pues dan el mismo resultado:
ˆ ∞ˆ ∞ ˆ 2π p
2 2 2 2
Qxx = ρ (r, z) 3r sen θdθ − 2π r + z rdrdz =
r=0 z=−∞ θ=0
ˆ ∞ ˆ ∞ ˆ 2π h p i ˚ h p i
2 2
= 2 2
ρ (r, z) 3r sen θ − r + z rdrdzdθ = ρ (r, z) 3r2 sen2 θ − r2 + z 2 dV
r=0 z=−∞ θ=0 V
Deshaciendo el cambio a cartesianas, recordando y = r sen θ, obtenemos:
˚ h p i
Qxx = ρ (x, y, z) 3y 2 − x2 + y 2 + z 2 dV = Qyy
V
3 Si la distribución de carga es antisimétrica, entonces es ρ (x, y, z) = −ρ (−x, −y, −z) ∀ (x, y, z). Por la
proposición 22 en la página 69 tenemos que:
˚
Qij = ρ R~ 3xi xj − R2 δij dV
V
Podemos dividir el volumen V en la parte que se encuentra sobre el plano z = 0 (V +) y la parte que se
encuentra por debajo del plano z = 0 (V −). Así:
˚ ˚
Qij = ρ (x1 , x2 , x3 ) 3xi xj − R2 δij dV + ρ (x1 , x2 , x3 ) 3xi xj − R2 δij dV =
V+ V−
˚
= ρ (x1 , x2 , x3 ) 3xi xj − x2 + y 2 + z 2 δij +
V+
i
+ρ (−x1 , −x2 , −x3 ) 3 (−xi ) (−xj ) − (−x)2 + (−y)2 + (−z)2 δij dV
Como ρ (x1 , x2 , x3 ) = −ρ (x1 , x2 , x3 ), tenemos que:
˚
Qij = ρ (x1 , x2 , x3 ) 3xi xj − R2 δij − 3xi xj − R2 δij dV = 0
V+
donde el primer sumando se anula al ser Q = 0 por hipótesis. En consecuencia, el momento dipolar es
el mismo desde cualquier sistema de referencia; dicho de otra forma, no depende del origen.
En este caso, como tenemos un número finito de cargas, la integral que usamos para calcular el momento
dipolar p~ se reduce a un sumatorio.
˚ X n
p~ = ~ ~
ρ R RdV = ~ i = ~0
qi R
V i=1
Como, por hipótesis, es p~ = ~0 y Q = 0, obtenemos que Qij = Q0ij . En otras palabras, el momento
cuadrupolar no depende del sistema de referencia.
Q.E.D.
σdS
+ + ~
dS
- + +
+ +
-
+ +
~0
R −σdS - ~r
-
e
~
R
Es decir, una capa dipolar son dos placas planas paralelas conductoras cargadas separadas una distancia e
tales que la densidad de carga evaluada en la intersección de cualquier recta t perpendicular a las superficies
con la primera superficie tiene justo el valor opuesto de la densidad de carga evaluada en la intersección de
esa misma recta t con la segunda superficie. En otras palabras, se trata de dos superficies planas paralelas S1
y S2 separadas una distancia e tales que:
σ1 R ~ = −σ2 R ~ + eŜ ∀R ~ ∈ S1
Proposición 26. El potencial generado por una capa dipolar de densidad de carga σ y espesor e viene dado
por la expresión:
¨
φ R~ = −τ dΩ = −
τ
Ω
4πε0 4πε0
~ y τ = σe.
donde Ω es el ángulo sólido con el que se ve la capa dipolar desde el punto R
Demostración. Usando la definición 13 en la página 45, obtenemos el diferencial de momento dipolar es:
+
p~ = q · d~
-
p = dq · ~e = σdS · ~e
d~
Recordando la proposición 18 en la página 46, obtenemos:
=:τ
z}|{ ~ ~
~ = 1 d~p · ~r 1 σe ~r · dS τ ~r · dS
dφ R 3
= 3
= 3
4πε0 r 4πε0 r 4πε0 r
dp
donde τ = = σe. Nótese que el dS va de la carga negativa a la positiva:
dS
dS
+ + + + + + +
- - - - - - -
~
−~r · dS
Por otra parte, por la definición 9 en la página 17, es justo el ángulo sólido visto desde el punto
r3
~ Por tanto:
R.
~
~ = − τ −~r · dS = −τ dΩ
dφ R
4πε0 r 3 4πε0
Integrando, obtenemos:
¨
~ −τ τ
φ R = dΩ = − Ω
4πε0 4πε0
Q.E.D.
~ en otro lugar, la representación gráfica habría sido:
Observación 19. Si hubiéramos cogido el punto R
σdS
+ ~
dS
+
−σdS -
+
-
+
-
~r -
e
Observación 20. Observemos una representación del potencial frente a la distancia a la capa dipolar en el caso
particular de σ = cte:
φ
τ
2ε0
r
−e2 e
2
−2τε
0
El eje horizontal representa la distancia con signo al centro de ambas placas. Es decir, la capa dipolar se
e e e
sitúa entre r = − y r = . En r = −∞ el ángulo sólido es 0 y justo al lado de la distribución (r = − ) el
2 2 2
τ
ángulo sólido es 2π; luego, justo en el borde de la capa, el potencial vale − . Al otro lado de la capa dipolar
2ε0
e τ
(r = ) el ángulo sólido vale −2π y el potencial . En r = ∞ el ángulo sólido se hace cero.
2 2ε0
Observación 21. Nótese que el ángulo sólido cambia de signo al pasar de un lado al otro de la capa dipolar.
− +~r
~r ~
dS
~r − +
− +
− +
~ = ~0 es:
Proposición 27. El campo eléctrico generado por una capa dipolar con τ tal que ∇τ
˛
~ = τ d~l × ~r
E
4πε0 C r3
donde la integral se realiza a través de la curva C que delimita la superficie de la capa dipolar.
~ = −∇φ
E ~
~ = τ ~
E ∇Ω
4πε0
Para ver cuánto cambia el ángulo sólido, vamos a hacer el siguiente dibujo:
P 0 ~δ
P
Ω
-δ dl
Así, podemos ver que lo que cambia el ángulo sólido cuando el punto de referencia se mueve ~δ es equivalente
a que el punto de referencia no se mueva, pero que la capa dipolar se desplace −~δ. En el plano, podemos ver
que la variación del ángulo sólido es sumar la «parte nueva» y restar «la parte que desaparece»; gráficamente
tendríamos:
dl
De manera que:
¨ ~
(−~r) · dS
dΩ = 3
S r
~ = −~δ × d~l:
Como dS
!
˛ ˛ ~r · d~l × ~δ ˛ d~l × ~r · ~δ ˛
~r d~l × ~r
dΩ = − 3 · −~δ × d~l = = = · ~δ
C r C r3 C r3 C r3
Por otra parte, por la regla de la cadena:
~
dΩ = ∇Ω · ~δ
Hemos obtenido dos expresiones diferentes para dΩ, así que ambas deben ser iguales:
˛ ~ !
dl × ~r
· ~δ = ∇Ω
~ · ~δ
C r3
~ ~ ~
~ B,
Ahora, por regla general, dados cuatro vectores A, ~ C ~ que cumplen A = B · C no es necesario
~ y D
~=D
A ~ ·C~
que B ~ Sin embargo, en nuestro caso dado que la igualdad debe cumplirse para cualquier vector ~δ,
~ = D.
forzosamente debe ser:
˛
~ = d~l × ~r
∇Ω
C r3
Y en consecuencia:
˛
~ = τ d~l × ~r
E
4πε0 C r3
Q.E.D.
Capítulo 5
Medios dieléctricos
5.1. Introducción
En la naturaleza aparecen varias estructuras polarizadas como por ejemplo una red de átomos de carbono
en el que uno de dichos átomos ha sido sustituido por uno de nitrógeno:
C C C C C
C C C N+ C
C C C C C
C C C C
También, al introducir un átomo neutro (como un gas noble) en un campo eléctrico los electrones cambian
sus disposición en torno al núcleo:
~
E
− Centro cargas − −
−
− −
− ++
− +++
++ +
− − + +
+
+ −
−
− −−
−
p~ 6= ~0
Definición 20. Llamamos sustancias polares a aquellas que, en ausencia de campos eléctricos externos, pre-
sentan momento dipolar. Las sustancias que no cumplen la condición anterior se llaman sustancias apolares
o sustancias no polares.
p~ 6= ~0
O
H
H 105o
De esta forma, en un vaso de agua, en general, al estar los dipolos distribuidos aleatoriamente, el momento
dipolar resultante es nulo (o casi nulo). Sin embargo, si sometemos las moléculas de agua a un campo eléctrico,
todas se orientan en la dirección de dicho campo eléctroestático:
~
E
Definición 21. Llamamos polarización por orientación al cambio de orientación que sufre una sustancia
polar al verse afectada por un campo eléctrico externo.
~
E
F
+
− +
F −
p~ = q · d~
donde d~ es el vector que une los centros de carga negativa y positiva.
Un ejemplo de polarización por deformación es:
~
E
− + − − + −
O C O O C O
p~T = ~0 p~T 6= ~0
Observación 22. La polarización por deformación sólo juega un papel importante en las sustancias no polares,
pues si bien existe la polarización por deformación en las sustancias polares, ésta es despreciable frente a la
polarización por orientación.
Observación 23. Nótese que la temperatura influye en el grado de polarización.
Conclusión 1.
Sustancia ¿Presenta momento dipolar en asuencia de campo eléctrico externo? Polarización por
polar Sí orientación
no polar No deformación
P~ := αE
~
donde α ∈ R es la polarizabilidad.
Proposición 28. Tener un medio material relleno de dipolos es equivalente a tener unas σb y ρb equivalentes,
donde:
(
σb = P~ · n̂
ρb = −∇~ · P~
σb
P~ (R)
~
σb = P~ · n̂
~r ρb = −∇~ · P~ ρb
~
R
~0
R
Corolario 6. Recordando la proposición 20 en la página 51, el potencial y el campo generados por un dieléctrico
(en ausencia de cargas reales (libres [no emparejadas]) deben ser:
‹ ˚
~ 1 σb 0 ρb 0
φ R = dS + dV
4πε0 S0 r V0 r
‹ ˚
~ R~ = 1 σb r̂ 0 ρb r̂ 0
E dS + dV
4πε0 S0 r2 V0 r2
En el caso en que además de carga equivalente hubiera cargas libres, por el principio de superposición,
tenemos:
‹ ˚
~ 1 σb + σf 0 ρb + ρf 0
φ R = dS + dV
4πε0 S0 r V0 r
‹ ˚
~ R
~ = 1 (σ b + σf ) r̂ 0 (ρ b + ρf ) r̂ 0
E dS + dV
4πε0 S0 r2 V0 r2
donde con f denotamos free (libre, en inglés).
~ := ε0 E
D ~ + P~
Teorema 9 (Teorema de Gauss generalizado (forma diferencial)). Las fuentes del desplazamiento eléctrico
~ evaluada
son las cargas reales (las cargas libres). Es decir, la divergencia del vector desplazamiento eléctrico D
~
en un punto R coincide con el valor de la función densidad de carga libre evaluada en dicho punto R. ~
~ ·D
∇ ~ R ~ = ρf R ~
Demostración. Por el teorema de Gauss en forma diferencial (ver teorema 4 en la página 24) tenemos:
ρtotal R ~ ρb R ~ + ρf R ~
~ ·E
∇ ~ R ~ = =
ε0 ε0
Por otro lado, habíamos dicho que:
~ ·E
ε0 ∇ ~ = −∇
~ · P~ + ρf ⇔ ε0 ∇
~ ·E
~ +∇
~ · P~ = ρf
~ ·
∇ ~ + P~
ε|0 E = ρf
{z }
~
=D
Q.E.D.
Teorema 10 (Teorema de Gauss generalizado (forma integral)). La integral del vector desplazamiento eléctrico
a lo largo de una superficie cerrada es igual a la carga libre encerrada en dicha superficie.
‹
D~ · dS
~ = Qf
S
Y, por el teorema de Gauss generalizado en forma diferencial (ver teorema 9 en la página anterior), llegamos
a: ˚ ˚
~ ~
∇ · DdV = ρf dV = Qf
V V
Q.E.D.
Proposición 29. Sea D ~ el vector desplazamiento eléctrico y P~ el vector polarización. Entonces, sus rotacio-
nales son iguales, siempre que el campo E~ sea conservativo.
~ ×D
∇ ~ =∇
~ × P~
Demostración* (No entra). Por la definición del vector desplazamiento eléctrico (ver definición 25 en la página
anterior), tenemos:
~ ×D
∇ ~ =∇ ~ × ε0 E ~ + P~
~ ×D
∇ ~ = ε0 ∇
~ ×E
~ +∇
~ × P~
~ ×E
Como el campo eléctrico es conservativo por hipótesis, por el corolario 3 en la página 26, es ∇ ~ = ~0, de
forma que obtenemos:
~ ×D
∇ ~ =∇
~ × P~
Q.E.D.
Observación 24. En general, el rotacional del desplazamiento eléctrico será cero, pues para conseguir que
el rotacional de la densidad de polarización sea no nulo hace falta recurrir a materiales artificiales como
metamateriales o materiales zurdos.
P~ = ε0 χE
~
La susceptibilidad eléctrica puede no ser constante. También puede ser un escalar o un tensor 0, 2 (una matriz),
como veremos más adelante:
ε := ε0 (1 + χ)
Proposición 30. En ausencia de polarización permanente, el vector desplazamiento eléctrico puede expresarse
como:
~ = ε0 E
D ~ + ε0 χ E
~ = εE
~
~ = ε0 E
D ~ + P~
~ = ε0 E
D ~ + ε0 χE
~ = ε0 (1 + χ) E
~
~ = εE
D ~
Q.E.D.
χ≥0
Sólo en materiales artificiales se puede conseguir χ < 0.
Material εr
vacío 1
aire 1, 006
teflón 2, 1
parafina 2, 1
poliestireno 2, 7
nylon 3, 5
cuarzo 5
amoníaco 22
metanol 34
glicerina 50
agua 81
titanato de bario (BaTiO3 ) (ferroeléctrico) 1200 − 1500
Definición 29. Llamamos medio isótropo a aquel en el que el vector polarización y el vector campo eléctrico
son paralelos en todo momento:
P~ k E
~
P~
~
E
χ11 χ12 χ13
La susceptibilidad sería en este caso un tensor (una matriz 3 × 3): χ = χ21 χ22 χ23
χ31 χ32 χ33
Definición 30. Llamamos medios lineales a aquellos en los que la susceptibilidad no depende del campo
eléctrico:
χ 6= F E~
~ .
Análogamente, llamamos medios no lineales a aquellos en los que χ = F E
Definición 31. Llamamos medios homogéneos a aquellos en los que la susceptibilidad no depende de la
posición:
~
χ 6= F R
~ .
Análogamente, llamamos medios no homogéneos a aquellos en los que χ = F R
A B C
Demostración. Nuestra situación es la siguiente:
(2) S2 ~2
dS ~2
D
~2 (2)
D
(1) (1) ~1
dS ~1
D
Es decir, Tenemos una pill-box («caja de pastillas»).
(2)
σf
(1)
Lo anterior se debe a que el valor de la integral a lo largo de la superficie lateral tiende a cero cuando la altura
transversal tiende a cero.
A continuación, tomamos un n̂ tal que dS ~2 = dS2 n̂ y dS~1 = dS1 (−n̂), de forma que:
¨ ¨ ¨ ¨
~ ~
D2 · dS2 + ~ ~
D1 · dS1 = D2n dS2 − D1n dS1
S2 S1 S2 S1
~ 1 · dS
D ~1 = D~ 1 (−n̂) dS1 = −D1n dS1
˜
Por otra parte, como podemos expresar Qreal = Qf = S σf dS, llegamos a:
¨ ¨
~ ~
D · dS = σf dS
S S
En consecuencia, haciendo tender la altura transversal a cero, S1 = S2 = S, por lo que obtenemos:
¨ ¨ ¨
D2n dS − D1n dS = σf dS
S S S
De esta forma, al ser la suma de las integrales la integral de la suma:
¨
(D2n − D1n − σf ) dS = 0
S
Como este argumento debe cumplirse para cualquier superficie sobre la que apliquemos el teorema de Gauss,
necesariamente:
D2n − D1n = [σf ]S
Q.E.D.
Corolario 7. Si en la frontera entre dos medios dieléctricos no hay carga libre, entonces las componentes
normales a la superficie del vector desplazamiento eléctrico serán iguales en ambos medios dieléctricos; es
decir, si [σf ]S = 0:
D2n = D1n
ε2 E2n − ε1 E1n = σf
y el caso particular σf = 0:
ε2 E2n = ε1 E1n
E1t = E2t
Demostración. Como el campo eléctrico es conservativo por hipótesis, por el corolario 3 en la página 26, se
cumple:
~ ×E
∇ ~ = ~0
h→0
~2
E ε2
(2) 2 t̂
4 3
t
(1) 1
ε1
~1
E
De nuevo, si hacemos tender la altura del rectángulo a cero, las integrales en los tramos (3) y (4) van a
tender a cero. Luego, obtenemos:
~ 2 · d~l2 = −E
E ~ 2 · t̂dl2 = −E2t dl
En consecuencia:
ˆ ˆ ˆ
E1t dl − E2t dl = 0 ⇔ (E1t − E2t ) dl = 0
Como este argumento funciona para cualquier espira que escojamos, debe ser:
Q.E.D.
D1t D2t
=
ε1 ε2
Demostración. Se obtiene trivialmente al aplicar la definición de vector desplazamiento eléctrico (ver defini-
ción 25 en la página 88) a la proposición 32 en la página anterior. Q.E.D.
Capítulo 6
Sistemas de conductores
Gas ionizado
Cátodo(−) Ánodo(+)
+ −
Como la masa de cualquier ion es mucho mayor que la masa de un electrón, el cátodo es el que recibe el
mayor aporte de energía. Por eso, en la práctica es preciso refrigerar el cátodo y reforzarlo para que no se
deforme.
~ = ~0
E
Demostración. Recordemos que nos encontramos en condiciones de electroestática (por hipótesis). Usaremos
reducción al absurdo. Supongamos que E ~ 6= ~0 en el conductor, entonces, una carga que se encuentre dentro
~ ~
del conductor sufrirá una fuerza F = q E ⇒ ~a 6= ~0, luego la carga sufrirá una aceleración y se moverá. En
consecuencia no nos encontramos en condiciones de equilibrio electroestático, pero habíamos supuesto que sí
que estábamos en dicha situación; hemos llegado a un absurdo. Q.E.D.
φ = cte
Demostración. Para cualquier trayectoria entre dos puntos A y B dentro de un conductor se da:
B
E=0
A E=0
ˆ A
φA − φB = − ~ · d~r
E
B
Como el campo eléctrico es nulo en un conductor:
ˆ A
~ · d~r = 0 ⇒ φA − φB = 0 ⇔ φA = φB
E
B
Ahora bien, en la superficie del conductor, como el campo eléctrico existe, su potencial asociado debe ser
una función diferenciable y, en particular, continua. Por tanto, el valor del potencial en la superficie de un
conductor debe ser su valor dentro del conductor. Q.E.D.
~ = ρ
~ ·E
∇
ε0
Como en un conductor (exceptuando la superficie) es E ~ = ~0, debe ser ρ = 0. Ahora, como no puede haber
cargas en un conductor (exceptuando la superficie), si éste está cargado, su carga debe situarse necesariamente
en su superficie. Q.E.D.
Corolario 12. La componente normal del vector desplazamiento eléctrico en el exterior del conductor coincide
con la densidad de carga libre que hay en su superficie. Además, la componente del campo eléctrico tangente
a la superficie del conductor es nula tanto en el interior como en el exterior del conductor.
D2n = σf
E1t = E2t = 0
D2n − D1n = σf
~ 1 = ~0, debe ser, por la definición 25
Ahora, como el campo eléctrico en el interior de un conductor es nulo E
en la página 88, D1n = 0. De forma que llegamos a:
D2n = σf
~ 1 = ~0, debe ser E1t = 0. Ahora, por la
Por otra parte, como el campo en el interior de un conductor es nulo E
proposición 32 en la página 94:
0 = E1t = E2t
Q.E.D.
~ es siempre perpendicular a la superficie en un conductor.
Corolario 13. El campo eléctrico E
Demostración. Se sigue trivialmente del corolario anterior. Si la componente tangencial del campo eléctrico en
la superficie de un conductor es siempre nula, entonces, únicamente puede tener componente normal. Q.E.D.
Proposición 34. Un conductor con un agujero en su interior apantalla el interior del campo eléctrico que halla
~ = ~0 en el agujero, independientemente
en su exterior. Es decir, si no hay cargas en dicho agujero, deber ser E
del campo que haya fuera.
Demostración. Por el teorema de Gauss (ver teorema 1 en la página 14) aplicado a la superficie interior del
conductor Si , sabemos:
‹
~ · dS
E ~ = 0 ⇒ QS = 0
i
Si
− +
− +B
A
− +
Como el campo eléctrico es conservativo, por el corolario 3 en la página 26, debe ser:
˛ ˆ B ˆ A
~ · d~l = 0 ⇔
E E~ · d~l + ~ · d~l = 0
E
C
|A {z } B
=0
~ = ~0
E
Q.E.D.
Observación 25. Nótese que el recíproco no es cierto. Si tenemos cargas en el interior de un conductor, fuera
del conductor sí que habrá campo eléctrico. Se puede comprobar fácilmente mediante la forma integral del
teorema de Gauss (ver teorema 1 en la página 14).
+ q
+
+
− −−
+ − q −
−− +
E=0
+ +
Observación 26. Si tenemos un conductor con cargas en su interior y lo conectamos a tierra, el conductor
apantalla el interior del exterior y el exterior del interior. Podemos verlo con el ejemplo del cable coaxial. Esto
se debe a que al conectar la superficie a tierra, todo el exceso de carga se va a tierra y desaparece.
E=0
−
−
q
−
−
−
φ=0
(2) (2)
φ1 , . . . , φn los valores del potencial en la superficie de los conductores en el estado (2). Por último, sean
(1) (2)
ρf y ρf las funciones densidades de carga libre del dieléctrico en los estados (1) y (2), respectivamente.
Asimismo, sean φ(1) y φ(2) las funciones potencial en el dieléctrico en los estados (1) y (2), respectivamente.
Supongamos, también, que el potencial en la superficie exterior que delimita nuestro dieléctrico vale cero.
Entonces, se cumple:
˚ n
X
(1) (2) (2) (1) (1) (2)
Pn φ(2) ρf − φ(1) ρf dV = φ j qj − φ j qj
T− j=1 Vj j=1
(1) (j)
Demostración. Como la carga en un conductor está distribuida a lo largo de la superficie, debe ser:
‹
(1) (1) ~j
qj = σj · dS (6.3.0.1)
Sj
donde n̂ es el vector unitario perpendicular a la superficie del conductor j-ésimo que se dirige hacia fuera de
dicho conductor.
Recordemos que por el corolario 13 en la página 98, en la superficie de un conductor no puede haber
componente tangencial. De aquí, obtendríamos una función potencial φ(1) .
Por ejemplo, en el caso de un condensador, tendríamos:
+Q −Q
(1) (2)
(1) (1)
∆V = φ1 − φ2
(1) (1)
q1 = Q q2 = −Q
(2) (2)
Volviendo a nuestro caso original, en una situación (2), las cargas ahora son q1 , . . . , qn . De nuevo:
‹
(2) (2)
qj = σj dSj (6.3.0.3)
Sj
Al igual que antes, de aquí obtendríamos una función potencial φ(2) . Si ahora aplicamos la segunda
identidad de Green (ver teorema 6 en la página 37):
˚ ‹
2 2
∂φ ∂ψ
ψ∇ φ − φ∇ ψ dV = ψ −φ dS
V S ∂ n̂ ∂ n̂
tomando ψ = φ(1) y φ = φ(2) , obtenemos:
˚ h ‹ !
i (2) (1)
(1) ∂φ (2) ∂φ
φ (1) 2 (2)
∇ φ −φ (2) 2 (1)
∇ φ dV = φ −φ dS (6.3.0.5)
V S ∂ n̂ ∂ n̂
Por otra parte, por la definición 12 en la página 28, por el teorema de Gauss Generalizado en forma
diferencial (ver teorema 9 en la página 88) y por la proposición 30 en la página 90, tenemos:
~ ·D
∇ ~ = ρf ⇔ ∇
~ · εE
~ = ρf ⇔ ε∇
~ ·E
~ = ρf ⇔ −ε∇ ~ = ρf ⇔ ∇2 φ = − ρf
~ · ∇φ
ε
En consecuencia, 6.3.0.5 es equivalente a:
!
˚ (2) (1) ‹
ρ ρ (1) ∂φ
(2)
(2) ∂φ
(1)
φ(1) − f − φ(2) − f dV = φ −φ dS ⇔
V ε ε S ∂ n̂ ∂ n̂
!
˚ (1) (2) ‹
ρ
(2) f
ρ
(1) f (1) ∂φ
(2)
(2) ∂φ
(1)
⇔ φ −φ dV = φ −φ dS
V ε ε S ∂ n̂ ∂ n̂
Haciendo uso de las ecuaciones 6.3.0.2 en la página 101 y 6.3.0.4 en la página anterior, la expresión anterior
se transforma en:
˚ (1) (2) ‹ (2) ! (1) !!
ρf ρf σj σj
φ (2)
−φ (1) dV = φ (1)
− −φ (2)
− dS ⇔
V ε ε S ε ε
˚ (1) (2) ‹ (1) (2) !
ρ
(2) f
ρ
(1) f
σ
(2) j
σ
(1) j
⇔ φ −φ dV = φ −φ dS (6.3.0.6)
V ε ε S ε ε
Como en el interior de los conductores no hay cargas libres, podemos restringir la integral al volumen del
dieléctrico menos el volumen de los conductores. Ahora, recordemos que T es el volumen de todo el espacio
considerado, entonces el volumen del dieléctrico excluyendo los conductores es:
n
X
V =T− Vj
j=1
Asimismo, podemos descomponer la integral en superficie en una suma de integrales a lo largo de las
superficies de todos los conductores más la integral a lo largo de la superficie exterior. Pero el potencial se
anula por hipótesis sobre dicha superficie exterior, luego la integral también.
h i h i ‹ (1) (2)
σ σ
φ(1) = φ(2) =0⇒ φ(2) j − φ(1) j dS = 0
Sexterior Sexterior Sexterior |{z} ε |{z} ε
=0 =0
Por último, multiplicando por ε a ambos lados y haciendo uso de las ecuaciones 6.3.0.1 en la página 101
y 6.3.0.3 en la página anterior, obtenemos:
˚ n
X
(1) (2) (2) (1) (1) (2)
Pn φ(2) ρf − φ(1) ρf dV = φ j qj − φ j qj
T− j=1 Vj j=1
Q.E.D.
Corolario 14 (Identidad de Gauss). Sean (1) y (2) dos estados de un sistema que contiene n conductores
(1) (1) (2) (2)
cargados con carga q1 , . . . , qn respectivamente en el estado (1) y con carga q1 , . . . , qn en el estado (2).
(1) (2)
Dichos conductores se encuentran dentro de un dieléctrico l.h.i. Además, sean φ1 , . . . , φn los valores del
(2) (2)
potencial sobre la superficie de los conductores en el estado (1) y sean φ1 , . . . , φn los valores del potencial
(1) (2)
en la superficie de los conductores en el estado (2). Por último, sean ρf = 0 = ρf las funciones densidades
de carga libre del dieléctrico en los estados (1) y (2), respectivamente. Supongamos, también, que el potencial
en la superficie exterior que delimita nuestro dieléctrico vale cero. Entonces, se cumple:
n
X
(2) (1) (1) (2)
φ j qj − φ j qj =0
j=1
(1) (2)
Demostración. Se sigue trivialmente al sustituir ρf = 0 = ρf en el teorema 11 en la página 100 y se usa
que la integral definida de cero es cero. Q.E.D.
q (1) q (2)
q (1) φ(2) = q (2) φ(1) ⇔ = = cte := capacidad del conductor
φ(1) φ(2)
C := P−1
φ1
Proposición 35. Supongamos que tenemos n conductores en un dieléctrico homogéneo. Sea Φ = ... el
φn
vector que contiene el valor del potencial sobre la superficie de cada uno de los conductores. Análogamente,
q1
..
sea Q = . el vector que contiene la carga de cada uno de los conductores. Además, sea P la matriz de
qn
coeficientes del potencial asociada a nuestro sistema. Entonces:
Φ=P·Q
1, q1
i, qi rij
dSi j
n
1 σi dSi
dφji =
4πε rij
En consecuencia, el potencial generado en el conductor j debido al conductor i es:
‹
1 σi dSi
φji =
Si 4πε rij
Como el dieléctrico es homogéneo, entonces:
‹
1 σi dSi
φji =
4πε Si rij
Por consiguiente el potencial en el conductor j es:
Xn ‹
1 σi dSi
φj =
4πε Si rij
i=1
‚
Si multiplicamos y dividimos por qj = Sj σj dSj y, además, también multiplicamos y dividimos por qi ,
llegamos a:
n
X ‹ ‹ X n
1 σi σj
φj = dSi dSj qi = pji qi
4πεqi qj Si Sj rij
i=1 | {z } i=1
=pji
⇔Φ=P·Q
Q.E.D.
−
−
+ −
−
−
+ −
+
−
+ −
−
− + −
− + −
Ejemplo 20. Para el caso de dos conductores, tendríamos:
φ1 = p11 q1 + p12 q2
φ2 = p21 q1 + p22 q2
Ahora, si ambas cargas son iguales (pero de signo contrario) q := q1 = −q2 , lo anterior se reduce a:
φ1 = p11 q − p12 q
φ2 = p21 q − p22 q
Si hallamos φ1 − φ2 :
φ1 − φ2
q=
p11 + p22 − 2p12
q = C (φ1 − φ2 )
Utilizando las cij (las inversas de pij ), llegamos a:
q1 = c11 φ1 + c12 φ2
q2 = c21 φ1 + c22 φ2
Si suponemos q1 = −q2 y sumamos obtendríamos:
c = −c21
0 = (c11 + c21 ) φ1 + (c12 + c22 ) φ2 ⇔ 11
c12 = −c22
Capítulo 7
Energía electrostática
Definición 35. Llamaremos energía electroestática de un sistema de cargas puntuales al trabajo necesario
para traer todas ellas desde el infinito una a una.
Para el caso de distribuciones volumétricas de carga en un medio dieléctrico, la energía electroestática será
el trabajo necesario para aumentar la densidad volumétrica de carga libre ρf desde la función nula en todos
los puntos hasta la distribución final.
Demostración. Podemos suponer, sin pérdida de generalidad, que nuestro sistema de cargas puntuales se
encuentra en el interior de un conductor, ya que por la proposición 34 en la página 98, un conductor aísla el
interior del exterior.
P
∞
q1 q2
q3
P
∞
Veamos qué trabajo tenemos que realizar para llevar la q1 desde la superficie del ∞ hasta el punto 1:
w1 = 0
El trabajo es nulo, pues no tengo ninguna carga en mi sistema todavía.
Ahora, quiero llevar la carga q2 al punto 2:
w2 = q2 φ2←1 ≡ q2 φ21
De manera que el trabajo total para dos cargas es:
1 q1 1 1 q1 q2 1 q1 q2
W2 = w1 + w2 = q2 φ21 = q2 = +
4πε0 r12 2 4πε0 r12 4πε0 r12
Como r12 = r21 , podemos expresar lo anterior como:
1 1
W2 = (φ21 q2 + φ12 q1 ) = (q1 φ12 + q2 φ21 )
2 2
Ahora, traigo la carga 3:
w3 = q3 φ3 = q3 (φ31 + φ32 )
Haciendo lo mismo que antes, obtenemos:
1 q1 1 q2 1
w3 = q3 + = (q3 φ31 + q1 φ13 + q2 φ23 + q3 φ32 )
4πε0 r13 4πε0 r23 2
1
w3 = (q3 φ31 + q1 φ13 + q3 φ32 + q2 φ23 )
2
En consecuencia:
1
W3 = (q1 φ12 + q1 φ13 + q2 φ21 + q2 φ23 + q3 φ31 + q3 φ32 )
2
De forma, que, en general:
n n n
1X X 1X
Wn = qi φij = qi φi = Ue
2 2
i=1 i=1
j=1
i 6= j
| {z }
=φi
Q.E.D.
P
∞
S0
ρ(R0 )
dV 0
C0
Como ahora estamos en un medio continuo, el número de partículas tiende a infinito, por lo que podemos
aproximar la suma por una integral:
˚
1
Ue u dq φq
2
Como es dq = ρ R~ 0 dV 0 , podemos expresar lo anterior como:
˚
1 ~ 0 dV 0
~0 φ R
Ue = ρ R
2 V0
Por el teorema de Gauss generalizado (ver teorema 9 en la página 88) ρ R~0 = ∇ ~ · D,
~ de forma que:
˚ ˚
1 ~0 φ R ~ 0 dV 0 = 1 ~ ·D
~ φ R ~ 0 dV
Ue = ρ R ∇ (7.2.0.1)
2 V 0 2 T.E
Por otra parte, por la regla del producto:
~ · φD
∇ ~ = φ∇
~ ·D
~ +D
~ · ∇φ
~ ⇔
~ φD
⇔∇ ~ −D
~ · ∇φ
~ = φ∇
~ ·D
~
1 Q2
Ue =
2 4πε0 a
También podemos hallar la solución de esta otra forma:
˚ ˆ 2π ˆ π ˆ ∞ 2
1 ~ · DdV
~ 1 1 Q
Ue = E =0+ ε0 r2 sen θdrdθdϕ =
2 T.E 2 ϕ=0 θ=0 r=a 4πε0 r2
←ϕ ←θ ˆ ∞
1 Q2 z}|{ z}|{ ∞ 1 1 Q2 1 1 Q2
= 2 2π (2) 2
dr = − =
2 (4π) ε0 a r 2 4πε0 r a 2 4πε0 a
Ejemplo 22. Tenemos una esfera dieléctrica:
Q
ε a
E=0
con densidad:
Q
ρ0 = 4 3
3 πa
Para r > a:
~ = 1 Q
E r̂
4πε0 r2
Y en consecuencia:
1 Q
φ=
4πε0 r
Y para r < a, por el teorema de Gauss generalizado (ver teorema 10 en la página 88):
4
D · 4πr2 = ρ0 πr3 ⇔
3
~ (r) = ρ0 r r̂
⇔D
3
En consecuencia:
~ (r) = ρ0 r r̂
E
3ε
De manera que:
ρ0 r2
φ (r) = − + k0
3ε 2
Claramente para r = a, por continuidad, debe darse:
ρ0 a2 1 Q Q ρ0 a2
− + k0 = ⇔ k0 = +
6ε 4πε0 a 4πε0 a 6ε
De forma que:
ρ0 2 Q
φ (r) = a − r2 +
6ε 4πε0 a
Bien, ahora para calcular la energía electrostática, hacemos:
˚ "ˆ ˆ ∞ 2 #
a 2
~ · DdV
~ 1 ρ0 2 2 1 Q
Ue = E = 2π · 2 r r dr + ε0 2
r2 dr
T.E 2 r=0 9ε a 4πε0 r
O, alternativamente, podríamos hacer:
˚ ˆ a
1 1 ρ0 2 2
Q
Ue = ρ0 φ (r) dV = 4πρ0 a −r + r2 dr
2 2 r=0 6ε 4πε0 a
P
∞
Ahora, vamos a estudiar cómo varía la energía potencial eléctrica U al variar la densidad voluḿetrica de carga
libre ρ. Para ello, aplicamos nuestro operador de variación virtual δ a ambos lados, obteniendo, por la regla
del producto:
h i h i
δU R ~ = δ ρf R ~ φ R ~ dV = δ ρf R ~ φ R ~ +ρ R ~ δ φ R ~ dV 0
¯ e = δU = δρf φdV 0
dδ̄U
Integrando, tendríamos:
ˆ ρfinal ˚
Ue = ~ 0
δρ φ R dV 0
(7.3.0.1)
ρ=0 V0
Físicamente, podemos ver la integral anterior como el hecho de ir introduciendo en cada dV un δρ poco a
poco hasta llegar a la ρfinal deseada en cada punto.
Por el teorema de Gauss generalizado (ver teorema 9 en la página 88) se da:
~ ·D
ρ=∇ ~
~ · φ δD
∇ ~ = φ∇
~ · δD
~ + δD
~ · ∇φ
~ ⇔ φ∇
~ · δD
~ =∇
~ · φ δD
~ − δD
~ · ∇φ
~ =∇
~ · φ δD
~ +E~ · δD
~
Corolario 15. La energía electroestática almacenada en una distribución de cargas en un medio dieléctrico
l.h.i. de permitividad ε puede expresarse como:
˚
1 2
Ue = εE dV
T.E. 2
Demostración. En un medio l.h.i tenemos ε = cte y al ser la diferencial un operador lineal, por la proposición 30
en la página 90:
~ = εE
D ~ ⇒ δD
~ = εδ E
~
Q.E.D.
(i)
(1)
(j) (n)
P
∞
Demostración. Por la proposición 37 en la página 110, restringiéndonos a la superficie de los conductores,
debe ser:
n ‹
1X
Ue = φi σi dSi
2 Si
i=1
Como los conductores son equipotenciales (ver corolario 10 en la página 97), φi es constante a lo largo de Si .
En consecuencia: ‹
n n n
1X 1X 1X
Ue = φi σi dSi = φ i qi = qi φ i
2 Si 2 2
i=1 i=1 i=1
Q.E.D.
Ejemplo 23. En el caso de dos conductores con cargas opuestas (influencia total), se obtiene:
1
Ue = C (∆V )2
2
Capítulo 8
Fuerzas electroestáticas
F~ = −∇U
~ e
(i)
(1)
ε = cte
(j) (n)
P
∞
¯ = −dUe
dW
Por otra parte, por definición de trabajo:
¯ = F~ · d~r
dW
Por tanto:
¯ = −dUe = F~ · d~r ⇔
dW
dUe
⇔ F~ = − ~ e
≡ −∇U
d~r
Q.E.D.
Observación 29. En general no podemos asegurar que los coeficientes del potencial se mantengan constan-
tes ante algún cambio en la geometría. Así que puede que no podamos calcular la energía electroestática
almacenada mediante la proposición 39 en la página 116.
− +
(1) (3)
− +
+Q
+Q +
ε0
−
ε ε0
ε
−Q −Q (2)
F~ = +∇U
~ e
P
∞
V1 Vj
1
ε j
i n
Vi Vn
Demostración. En este caso, como estamos manteniendo el potencial constante de forma externa, estamos
introduciendo energía en el sistema. Podemos suponer, sin pérdida de generalidad, que mantenemos el potencial
constante en los conductores por medio de baterías. Llamaremos Ub a la energía electrostática de las baterías.
Por el teorema de la energía mecánica se tiene que cumplir:
dUb = dW
¯ + dUe (8.2.0.1)
¯ = F~ · d~r. Por otra parte, por la proposición 37 en la página 110, restringiéndonos a la superficie de
donde dW
los conductores, debe ser:
n ‹
1X
Ue = σi φi dSi
2 Si
i=1
Como los conductores son equipotenciales (ver corolario 10 en la página 97), el potencial φi es constante para
la superficie Si . Así:
n ‹ n
1X 1X
Ue = φi σi dSi = φ i qi
2 2
i=1 i=1
donde el segundo sumando se anula, porque el potencial es constante por hipótesis. En consecuencia:
n
1X
dUe = dqi φi (8.2.0.2)
2
i=1
Por otra parte, como se cumple Ub,i = qi φi por definición de potencial, tenemos:
n
X
Ub = qi φ i
i=1
Tomando la diferencial a ambos lados, por el mismo argumento que para dUe , llegamos a:
n
X
dUb = dqi φi (8.2.0.3)
i=1
dUb = 2dUe
En consecuencia, sustituyendo en la ecuación diferencial del trabajo (ver 8.2.0.1 en la página anterior),
obtenemos:
¯ + dUe ⇔ dW
2dUe = dW ¯ = dUe
Por la definición de trabajo, tenemos:
dUe
F~ · d~r = dUe ⇔ F~ = ~ e
= ∇U
d~r
Q.E.D.
Observación 30. Nótese la diferencia entre la situación a carga constante (signo negativo del gradiente) y la
situación a potencial constante (signo positivo del gradiente).
8.3. Ejemplos
Ejemplo 24 (Problema 48). Tenemos la siguiente situación:
y
+Q
F
D1 ε D 2 ε0
x
−Q
L
x
~ = E ŷ
E
~ = (0, D, 0)
D
Nos piden determinar la fuerza total que actúa sobre la superficie del dieléctrico en dos situaciones.
σ σ0
ε E ε0
V0
Obtendríamos:
~ = V0 (−ŷ)
E
d
2
1 V0 V02
Ue = ε 2 xhd + ε0 2 (L − x) hd =
2 d d
1 V02
h [εx + ε0 (L − x)]
=
2 d
Por la proposición 41 en la página 119, tenemos:
dUe 1 V02
Fx = = h (ε − ε0 )
dx 2 d
dx
Demostración. Tomando la diferencial a ambos lados en la proposición 37 en la página 110 y considerando
dV = dSdx, tenemos:
1~ ~
dUe = D · EdSdx
2
Como las placas están sometidas a potencial constante, por la proposición 41 en la página 119, tenemos:
~ e ⇔ dF dUe
¯ = ∇U
dF ¯ =
d~r
donde la notación de diferencial inexacta proviene de que estamos restrigiéndonos a la fuerza que actúa sobre
un dS.
Como en este caso la fuerza actúa exclusivamente en el eje X:
1~ ~
dUe D · EdSdx 1~ ~
dF
¯ = = 2 = D · E dS
dx dx 2
Despejando, llegamos a la definición de presión:
dF
¯ 1~ ~
P= = D ·E
dS 2
Por la proposición 30 en la página 90, podemos escribir lo anterior como:
1 D2
P=
2 ε
Por último, por el corolario 12 en la página 98, tenemos D = σ en las placas del condensador. Así:
σ2
P=
2ε
Q.E.D.
Capítulo 9
9.1. Introducción
~ t a la variación de la carga con respecto al
Definición 37. Llamamos intensidad de corriente eléctrica I R,
~ en un instante t.
tiempo que tiene lugar en un punto R
~ t := dQ R,
I R, ~ t
dt
Observación 31. A partir de este momento el campo eléctrico ya no tendrá por que ser conservativo. Es decir,
en general, será:
~ ×E
∇ ~ 6= ~0
D E
donde ρ R, ~ t
~ t es la densidad volumétrica de carga libre y ~v R, es la velocidad promedio de las cargas.
Sus unidades son:
h i C m C
J~ = 3 = 2
m s m s
En caso de régimen estacionario de corrientes, podemos prescindir de la dependencia del tiempo, obteniendo
una función que depende únicamente de la posición:
D E
J~ R~ =ρ R ~ v R~
Proposición 43. La intensidad de corriente eléctrica total I que circula por un conductor cilíndrico es igual
a la integral a lo largo de una sección del conductor de área S ~ del producto escalar de la densidad volumétrica
de corriente por el diferencial de superficie.
¨
I= J~ · dS
~
S
Demostración. Imaginemos:
dS
dl
~ ~l
dQ = ρdSd
Por definición de intensidad de corriente eléctrica (ver definición 37 en la página 125), tenemos:
dQ ~
dI
¯ = = ρ h~v i · dS
dt
donde la diferencia inexacta viene de que estamos hallando la corriente eléctrica dI¯ que circula por un dS. Por
definición de densidad de corriente eléctrica (ver definición 39 en la página 125), podemos expresar lo anterior
como:
¯ = J~ · dS
dI ~
Para hallar la que circula por todo nuestro cilindro, debemos integrar a su superficie transversal:
¨
I= J~ · dS
~
S
Q.E.D.
Observación 33. Para ver la interpretación física del producto escalar que aparece en el enunciado de la
proposición 43 en la página anterior, podemos recurrir al siguiente dibujo de dos cañerías que trasportan
agua:
~
dS
En ambas situaciones, el agua que sale de las tuberías es la misma, aunque la superficie del segundo caso
sea mayor.
Corolario 17. Si la densidad de corriente eléctrica que circula a lo largo de un conductor cilíndrico es
constante, entonces, la intensidad de corriente eléctrica que circula por una sección del cilindro de área S
viene dada por:
I = J~ · S
~
−→
Demostración. Se sigue trivialmente de la proposición 43 en la página anterior. Al ser J~ = cte, J~ sale de la
integral y se llega al resultado. Q.E.D.
Corolario 18. La intensidad de corriente eléctrica total I que circula por un conductor plano es igual a la
integral a lo largo de una sección (unidimensional) del conductor de longitud L del producto de la densidad
superficial de corriente por el diferencial de superficie.
ˆ
I= JS dl
L
dl
Sea
Demostración. SCPla superficie
del conductor plano. Partimos de la
proposición 43 en la página 126
~ t = J~S R,
tomando J~ R, ~ t = ~0 en el resto de casos. Como
~ t en la superficie del conductor plano y J~ R,
las superficies SCP y S son perpendiculares, el núcleo de la integral únicamente será no nulo en la intersección
de ambas superficies. De esta forma: ˆ
I= JS dl
L
Q.E.D.
dS
~ t)
J(R,
~ t)
ρ(R,
Por una parte, la carga que «se escapa» a lo largo de un dS de la superficie de nuestro volumen V es:
J~ · dS
~
~ va hacia fuera del volumen. Por tanto, la variación de carga en dicho punto de la superficie es:
ya que dS
dQ
d¯ = −J~ · dS
~
dt
Para hallar cuál es la variación de carga a lo largo de toda la superficie, tenemos que integrar:
‹
dQ
=− J~ · dS
~
dt S
Sustituyendo, obtenemos: ˚ ‹
∂
ρdV =− J~ · dS
~
∂t V S
Como la integral es a lo largo del espacio y la derivada es con respecto al tiempo, ambos operadores conmutan.
Por otra parte, mediante el teorema de la divergencia (ver teorema 2 en la página 17), podemos reescribir la
parte derecha de la ecuación; así, obtenemos:
˚ ˚ ˚
∂ρ ~ ~ ∂ρ ~ ~
dV = − ∇ · JdV ⇔ + ∇ · J dV = 0
V ∂t V V ∂t
Como lo anterior debe ocurrir también para cualquier subvolumen de V , necesariamente será:
∂ρ ~ ~
+∇·J =0
∂t
Q.E.D.
Corolario 19. En régimen estacionario de corrientes, la densidad de corriente eléctrica no tiene fuentes
escalares; dicho de otra forma, la divergencia de la densidad volumétrica de corriente es nula en cualquier
punto.
~ · J~ R
∇ ~ =0
∂ρ
Demostración. Por el corolario 16 en la página 125, en régimen estacionario de corrientes, se cumple = 0.
∂t
∂ρ
El resultado se obtiene trivialmente al sustituir = 0 en el enunciado del teorema 12 en la página anterior.
∂t
Q.E.D.
Demostración. Por lo visto en la sección 2.4.1 en la página 32, sabemos que la solución es de la forma:
con α ∈ C si α 6= 0 y:
X (x) = Ax + B
~
si α = 0. Fijémonos en la primera componente de ∇ϕ:
∂ϕ ∂X ∂ ∂X
= Y (y) Z (z) + X (x) [Y (y) Z (z)] = Y (y) Z (z)
∂x ∂x |∂x {z } ∂x
=0
∂ϕ
Si evaluamos en el punto (x, 0, 0), podemos escribir lo anterior como:
∂x
∂ϕ ∂X
(x, 0, 0) = (x) C
∂x ∂x
donde C es una constante, de hecho C = Y (0) Z (0).
Antes de proseguir, nótese que:
p
~ =
(x, 0, 0) → (∞, 0, 0) ⇒ R x2 + y 2 + z 2 → ∞
∂ϕ
(x, 0, 0) = Aαeαx − Bαe−αx C
∂x
Claramente, en la expresión anterior, si hacemos tender x → ∞ o x → −∞, únicamente puede suceder que
~ ∈ R3 y, por tanto,
la expresión se anule en ambos casos si A = 0 = B. Entonces, necesariamente es ϕ = 0 ∀R
~ = ~0 ∈ R3 .
∇ϕ
Por otra parte, si X (x) es de la forma X (x) = Ax + B, entonces:
∂ϕ
(x, 0, 0) = AC
∂x
Y, evidentemente, si tomamos el límite de la expresión anterior cuando x → ∞ o x → −∞, la función se
anula únicamente si A = 0. En ese caso, podríamos repetir el razonamiento anterior con las funciones Y (y)
y Z (z) y llegaríamos a que necesariamente ϕ = K ∀R ~ ∈ R3 donde K es una constante. Por tanto sería,
~ = ~0 R
∇ϕ ~ ∈R . 3 Q.E.D.
~ ∈ R3 .
para al menos un R
~ ∈ R3 .
para al menos un R
Demostración.
1. Sabemos que si es ∇ ~ × F~ = ~0 ∀R ~ ∈ R3 , entonces existe una función ϕ : R3 −→ R tal que F~ = −∇ϕ. ~
~ ~ ~
Vamos a probar el enunciado por reducción al absurdo. Supongamos que se cumple ∇ · F = 0 ∀R ∈ R3 . ~
Entonces, debe cumplirse ∇2 ϕ = 0 ∀R ~ ∈ R3 . Por el lema 3 en la página 129, sabemos que toda solución
de la ecuación anterior que satisface que lı́m ∇ϕ~ = ~0, cumple también ∇ϕ ~ R~ = ~0 ∀R ~ ∈ R3 . En otras
|R|→∞
~
palabras, en ese caso sería F = 0 y eso es absurdo, porque por hipótesis era F~ 6= ~0. Así, necesariamente,
~ ~
~ · F~ 6= 0 en al menos un R
∇ ~ ∈ R3 .
~ × F~ = ~0 ∀R
2. De nuevo, vamos a aplicar reducción al absurdo. Si fuese ∇ ~ ∈ R3 , entonces se cumpliría (1)
~ ~ ~
y tendría que ser ∇ · F 6= 0 en al menos algún R ∈ R ; pero, por hipótesis, era ∇
3 ~ · F~ = ~0 ∀R
~ ∈ R3 .
Hemos llegado a contradicción. En consecuencia, deber ser ∇~ × F~ 6= ~0 en al menos un R ~ ∈R .
3
Q.E.D.
~
E
e-
Este modelo se llama modelo de Drüde. La idea es que un electrón se ve acelerado por el campo eléctrico
y luego choca con un núcleo, transfiriéndole toda su energía, deteniéndose. De nuevo, se ve acelerado por el
campo eléctrico y choca contra otro núcleo y así, sucesivamente. Por la segunda ley de Newton, tenemos:
~ = m~a ⇔ a = qE
F~ = q E
m
Como tenemos tres grados de libertad, por termodinámica, tenemos:
3 1
kT = mv 2
2 2
También se da:
qE
vf = τ
m
donde τ es el tiempo de vuelo libre del electrón (el tiempo entre choque y choque).
Con todo esto, podemos suponer que la velocidad de arrastre es proporcional al campo eléctrico:
~
hvi = µE
donde µ es la movilidad electrónica.
En consecuencia, macroscópicamente, no es descabellado pensar que:
J~ ∝ E
~
J~ = ρ hvi = ρµ E~ =σ·E
~
|{z}
=:σ
¡No confundir con la densidad de carga superficial! Las unidades de la conductividad son:
A
m2 A 1 S
[σ] = = = = Ω−1 m−1 =
V
m
V·m Ω·m m
donde S denota «Siemens». Dependiendo del medio σ puede ser un escalar o un tensor (0, 2).
Definición 42. Llamamos resistividad ρ de un material a la inversa de su conductividad. ¡No confundir
con la densidad de carga volumétrica!
Análogamente a lo visto en electrostática, tenemos las siguientes definiciones:
Definición 43. Llamamos medio isótropo a aquel en el que el vector densidad volumétrica de corriente y el
vector campo eléctrico son paralelos en todo momento:
J~ k E
~
Análogamente, llamamos medios anisótropos a aquellos en los que J~ ~
∦ E.
σ11 σ12 σ13
La conductividad sería en este caso un tensor (una matriz 3 × 3): σ = σ21 σ22 σ23
σ31 σ32 σ33
Definición 44. Llamamos medios lineales a aquellos en los que la conductividad no depende del campo
eléctrico:
σ 6= F E~
~ .
Análogamente, llamamos medios no lineales a aquellos en los que σ = F E
Definición 45. Llamamos medios homogéneos a aquellos en los que la conductividad no depende de la
posición:
~
σ 6= F R
~ .
Análogamente, llamamos medios no homogéneos a aquellos en los que σ = F R
~ ×E
∇ ~ = ~0
J~
B
A
~ ×E
∇ ~ 6= 0
Demostración. Como estamos en régimen estacionario de corrientes, por el corolario 19 en la página 129, es
~ · J~ = 0. Ahora, como el campo eléctrico se anula por el infinito y es J~ 6= ~0 ya que por el cable circula
∇
una corriente no nula, por el teorema 13 en la página 130, en al menos algún punto debe ser ∇ ~ × J~ 6= ~0.
Por la definición 41 en la página anterior, tenemos J~ = σ E.
~ De esta forma, como el cable es un medio lineal,
homogéneo e isótropo por hipótesis, necesariamente en al menos algún punto R ~ se cumplirá:
~ × J~ 6= ~0 ⇔ ∇
∇ ~ × σE ~ 6= ~0 ⇔ σ ∇
~ ×E~ 6= ~0 ⇒ ∇
~ ×E~ 6= ~0
de forma que el campo eléctrico no será conservativo en dichos puntos. En consecuencia, nos conviene des-
componer el campo eléctrico en dos sumandos:
~ =E
E ~c + E
~ nc
donde E~ c es conservativo y E
~ nc no lo es.
A continuación, examinemos qué ocurre en el «punto de contacto» de ambos extremos del cable. Por
hipótesis, por ese punto no circula corriente, luego será J~ = ~0. Por la definición 41 en la página 132, tenemos:
~0 = J~ = σ E
~ ⇔E
~ = ~0 ⇔ E
~c + E
~ nc = ~0 ⇔ E
~ nc = −E
~c
en dicho punto.
Por otra parte, en el resto del cable, tenemos I = cte y, como (por hipótesis) la densidad volumétrica de
corriente J~ es constante a lo largo de cada sección transversal del cable y, además, es paralela a las paredes
del cable (en consecuencia, J~ tiene la misma dirección que S,
~ pues S es una sección transversal), se tiene por
el corolario 17 en la página 127:
−
→
cte = I = JS ⇒ J~ = cte ⇒ ∇ ~ × J~ = ~0 ⇒ ∇~ ×E ~ = ~0
en todos los puntos del cable, salvo en el punto de «contacto» de ambos extremos. El último paso se debe a
la definición 41 en la página 132. Dicho de otra forma, únicamente hay un campo eléctrico no conservativo en
el punto de contacto.
Por el teorema de Stokes (ver teorema 3 en la página 17), tenemos:
¨ ¨ ˛
0 6= ~ ~ ~
∇ × E · dS = ~ ~ ~
∇ × Enc · dS = E~ nc · d~l (9.5.2.1)
SA SA C
donde SA es cualquier superficie abierta que se apoya sobra la curva que forma el cable y C es dicha curva.
~ nc = ~0 en todo punto salvo en el punto de contacto; en consecuencia:
Por lo mencionado antes, E
˛ ˆ A
~ nc · d~l =
E ~ nc · d~l
E
C B
donde B y A son los extremos del cable (suponemos que la corriente circula de A a B). Pero, habíamos dicho
que en el punto de contacto era:
~ nc = −E
E ~ c = ∇φ
~
donde la última igualdad se debe a que, como E ~ c es conservativo, sabemos que existe una función escalar tal
~ c = −∇φ.
que E ~ Por tanto, tenemos:
ˆ A ˆ A ˆ A ˆ A
~ ~ ~ ~ dφ ~
Enc · dl = ∇φ · dl = · dl = dφ = φ (A) − φ (B) = ε (9.5.2.2)
B B B d~ l B
Por otro lado, por el teorema de Stokes (ver teorema 3 en la página 17):
¨ ˛ ˛
∇~ ×E ~ · dS~= ~ ×E
∇ ~ · dS
~= ~ · d~l
E (9.5.2.3)
SA C C
Por último, podemos calcular la integral del lado izquierdo de esta expresión de otra forma, a través de la
definición 41 en la página 132:
˛ ˛ ~
~ ~ J ~
ε= E · dl = · dl
C C σ
Como, por hipótesis J~ era paralelo a las paredes del cable, será J~ k d~l y, por consiguiente, podemos reescribir
lo anterior como: ˛
J
ε= dl
C σ
Para concluir, podemos expresar lo anterior de otra manera a partir de la definición 42 en la página 132, al
1
ser ρ = :
σ ˆ B ˆ B
dl ρ
ε=I =I dl
A σS A S
Q.E.D.
Definición 46. En régimen estacionario de corriente, tenemos un cable (un medio l.h.i.) por el que circula una
corriente I 6= 0. Supongamos que en cada sección transversal del cable se cumple que la densidad volumétrica
de corriente J~ es constante, J~ es paralela a la pared del cable en todo momento, los extremos de dicho cable
están «conectados», pero que entre ellos no circula corriente y el campo eléctrico se anula en el infinito.
Entonces, llamaremos fuerza electromotriz ε a la diferencia de potencial ε entre los dos extremos A y B
del cable que aparece en la proposición 14 en la página 133.
Definición 47. En régimen estacionario de corriente, tenemos un cable (un medio l.h.i.) por el que circula una
corriente I 6= 0. Supongamos que en cada sección transversal del cable se cumple que la densidad volumétrica
de corriente J~ es constante, J~ es paralela a la pared del cable en todo momento, los extremos de dicho cable
están «conectados», pero que entre ellos no circula corriente y el campo eléctrico se anula en el infinito.
Entonces llamaremos resistencia entre dos puntos A y B de un material al término R que aparece en el
enunciado del teorema 14 en la página 133.
I
R
ε
Demostración. El trabajo realizado por la fuerza eléctrica sobre un dV de carga dq viene dado por:
¯ = F~e · d~r = dq E
dW ~ · d~r = E
~ · (dq d~r) (9.6.0.1)
Por la definición 39 en la página 125, tenemos que la densidad volumétrica de corriente eléctrica en dicho dV
es :
J~ = ρ h~v i
dq d~r
donde ρ = y podemos expresar h~v i como . Así, obtenemos:
dV dt
dq d~r
J~ = ~ dt = dq d~r
⇔ JdV
dV dt
Sustituyendo dq d~r en la ecuación 9.6.0.1, obtenemos:
dW ~ · JdV
¯ =E ~ dt
Q.E.D.
2.
E2t = E1t
Demostración.
1.
Por el corolario 19 en la página 129, en régimen estacionario de corrientes, se da:
~ · J~ = 0
∇
Ahora, escojamos un cilindro con base de área S y de altura h tal que una de sus bases está en un
medio conductor, la otra en el otro y ambas bases son paralelas a la superficie de contacto entre ambos
conductores. Es decir, estamos cogiendo una «caja de pastillas» o «pillbox». Entonces:
˚
~ ~
∇·J =0⇒ ∇~ · JdV
~ =0
V
donde SC es la superficie del cilindro. Ahora bien, podemos descomponer dicha integral en:
‹ ¨ ¨ ¨
0= J~ · dS
~= J~ · dS
~+ J~ · dS
~+ J~ · dS
~
SC tapa superior tapa inferior superficie lateral
Lo anterior se debe a que el valor de la integral a lo largo de la superficie lateral tiende a cero cuando
la altura transversal h tiende a cero.
A continuación, tomamos un n̂ tal que dS ~2 = dS2 n̂ y dS
~1 = dS1 (−n̂), de forma que:
¨ ¨ ¨ ¨
0= J~2 · dS
~2 + J~1 · dS
~1 = J2n dS2 − J1n dS1
S2 S1 S2 S1
J~1 · dS
~1 = J~1 (−n̂) dS1 = −J1n dS1
Q.E.D.
Dieléctrico Conductores
~ = εE
D ~ J~ = σ E
~
~ ~
∇·D =0 ~ ~
∇·J =0
D2n = D1n J2n = J1n
E2t = E1t E2t = E1t
+q
~ · J~ + ∂ρ = 0
∇ (9.8.0.1)
∂t
Como el medio es l.h.i., tenemos:
J~ = σ E
~ ⇒∇
~ · J~ = σ ∇
~ ·E
~
~ = εE
D ~ = 1D
~ ⇔E ~
ε
~ · J~ como sigue:
Como el medio es l.h.i., estamos en disposición de expresar ∇
~ · J~ = σ ∇
∇ ~ ·D
~
ε
Por el teorema de Gauss generalizado en forma diferencial (ver teorema 9 en la página 88):
~ ·D
∇ ~ =ρ
σ
σ
~ t = ρ0 R
⇔ ρ = ρ0 e− ε t ⇔ ρ R, ~ e− ε t
Q.E.D.
Observación 35. Lo que nos da el carácter dieléctrico o conductor de un material es, precisamente, el cociente
σ
que aparece en la proposición 45 en la página anterior. Si este cociente es muy alto, estamos ante un
ε
conductor, mientras que si es muy bajo, estamos ante un dieléctrico.
Definición 48. Llamamos tiempo de relajación τ al tiempo que transcurre hasta que la densidad de carga
de un punto de un medio l.h.i. se ha reducido en un factor e. Es decir:
ε
τ := = ερ
σ
donde ε es la permitividad eléctrica del medio, σ es su conductividad y ρ es su resistividad.
Observación 36. En la prática, podremos considerar que el conductor alcanza el equilibrio electroestático en
un tiempo 5τ . Si cogemos como ejemplo el cobre, obtenemos que su tiempo de relajación es:
9.9. Ejemplos
Ejemplo 25 (Problema 52).
2d σ2 = σ0
σ1 (x) = σ0 2d+x
(1) (2)
S0
x=0 x = 2d x = 3d
2d
σ1 (x) = σ0
2d + x
σ2 (x) = σ0
φ (x = 0) = V0
φ (x = 3d) = 0
Nos piden determinar J~ (x), E
~ (x), P . Para ello:
~ i = Ei x̂ ∀i = 1, 2
E
J~i = Ji x̂
Como estamos en régimen estacionario de corrientes, por el corolario 19 en la página 129:
~ · J~1 = 0 ⇔ J1 = cte
∇
~ · J~2 ⇔ 0 ⇔ J2 = cte
∇
En x = 2d:
J2 = J1 =: J
ˆ
~ J~ J (2d + x) J x2
E1 = = x̂ ⇒ φ1 (x) = − E1 dx = − 2xd + + k1
σ1 2dσ0 2dσ0 2
~
~ 2 = J = J x̂ ⇔ σ2 (x) = − J x + k2
E
σ2 σ0 σ0
φ1 (x = 0) = V0 = k1
J 3Jd
φ2 (x = 3d) = 0 = − 3d + k2 ⇔ k2 =
σ0 σ0
J 3Jd J2d
φ1 (x = 2d) = φ2 (x = 2d) ⇔ − 4d2 + 2d2 + V0 = − ⇔
2dσ0 V0 σ0
J3d Jd 4Jd σ0 V 0
⇔− + V0 = ⇔ V0 = ⇔J =
σ0 σ0 σ0 4d
De forma que:
σ0 V 0
J~ = x̂
4d
σ0 V0
~1 = 4d (2d + x) V0
E x̂ = 2 (2d + x) x̂
2dσ0 8d
~ 2 = V0 x̂
E
4d
Para calcular la potencia, aplicamos la ley de Joule (ver teorema 15 en la página 135):
˚ ˚
P = u1 dV + u2 dV
V1 V2
(1) (2)
a X
Por la proposición 44 en la página 136, se cumple J1 = J2 y se crea una densidad de carga superficial en la
l1 l2
superficie de contacto entre ambos dieléctricos para que se cumpla D2n − D1n = [σf ]S . R1 = . R2 = .
σ1 S σ2 S
Capítulo 10
El campo magnetoestático
~ · J~ = 0
∇ I = JS ε = IR
donde R es la resistencia del circuito, S es el área de la sección transversal del cilindro, I es la intensidad de
corriente eléctrica que circula por él y ε es la diferencia de potencial entre sus extremos.
dl
I
J~ R
~ ×E
∇ ~ = ~0 ε
Demostración. La primera propiedad viene dada directamente en la definición 49; la segunda se deduce del
corolario 17 en la página 127, siendo por hipótesis J~ k dS.
~ Por último, la tercera se cumple porque se dan los
requisitos para poder aplicar la ley de Ohm (ver teorema 14 en la página 133). Q.E.D.
Observación 37. Los circuitos lineales estacionarios nos serán realmente útiles porque nos permitirán descom-
poner cualquier corriente que circula en un volumen en un conjunto de corrientes en circuitos elementales.
dl
J~ k dS
~
J~
10.2. ~
Ley de Biot-Savart y la inducción magnética B
Observación 38 (¡Ambigüedad de la nomenclatura!). Existen dos magnitudes físicas a las que usualmente
se les da el nombre de «campo magnético» o «intensidad de campo magnético»; dichas magnitudes son la
inducción magnética B~ y la excitación magnética H. ~ En los capítulos siguientes se desarrollarán y definirán
ambos conceptos. Con el fin de evitar confusiones, en este texto se usarán los nombres «inducción magnética» y
«excitación magnética»; pero el lector debe estar al tanto de que en muchos textos se usa «campo magnético»
o «intensidad de campo magnético» indistintamente para ambas.
dl
~r
~0
R
~
R
Proposición 46. La inducción magnética B ~ en un punto R ~ generada por un conductor lineal por el que
circula una corriente eléctrica I viene dada por la expresión:
µ I ˛ d~l × ~r
~
B R~ = 0
4π L r3
donde la integral se realiza a lo largo del conductor lineal, d~l lleva la dirección de la corriente I y ~r = R
~ −R
~ 0,
~ 0
siendo R la posición de cada dl. ~
Ejemplo 27. Queremos obtener el campo magnético generado por un anillo de radio a a lo largo de su eje
de simetría que es perpendicular al plano que lo contiene.
α θ
Z
α θ
µ0 I dl r
dBz = cos θ
4π r3 | {z }
=sen α
Proposición 47. La inducción magnética B ~ generada por un conductor rectilíneo indefinido por el que circula
una corriente I en coordenadas cilíndricas viene dada por la expresión:
~ (r) = µ0 I ϕ̂
B
2πr
donde r indica la distancia al conductor.
I ~
B
I
dl
θ
θ
r
dl
µ0 I d~l × ~r µ0 I dl r π µ I dl
~ = 0
dB = sen + θ = cos θϕ̂
4π r3 4π r3 2 4π r2
Por otra parte:
dl = rdθ
Así, sustituyendo, tenemos:
~ = µ0 I 1
dB cos θdθϕ̂
4π r
Integrando a ambos lados:
ˆ π
~ = µ0 I µ0 I µ0 I µ0 I
2 π
B cos θdθϕ̂ = [sen θ]−2 π ϕ̂ = [1 − (−1)] ϕ̂ = ϕ̂
4πr − π2 4πr 2 4πr 2πr
Q.E.D.
dF~ = Id~l × B
~
Teorema 16 (Ley de Fuerza de Ampère). La fuerza que un circuito elemental por el que circula una corriente
I1 ejerce sobre un circuito elemental por el que circula una corriente I2 obedece la expresión:
˛ ˛ d~l2 × d~l1 × ~r12
µ0
F~2←1 = I1 I2 3
4π L2 L1 r12
donde ~r12 es el vector que une los diferenciales d~l1 y d~l2 , llevando d~l1 la dirección de I1 y d~l2 la dirección de
I2 . Las integrales se realizan a lo largo de la longitud de cada circuito lineal.
C1
I1
dl1
C2
I2
~r12
dl2
Demostración. Por la proposición 46 en la página 144, sabemos que la inducción magnética generada por el
circuito elemental 1 en un punto del circuito elemental 2 viene dada por:
µ I ˛ d~l × ~r
~
B R~ = 0 1 1 12
(10.3.0.1)
4π L1 3
r12
Por otra parte, por el axioma 3 en la página anterior, sabemos que la fuerza dF~ producida en un dl2 a causa
de la inducción magnética viene dada por:
dF~ = I2 d~l2 × B
~
Q.E.D.
Explicaremos el factor µ0
4π más adelante. Por la tercera ley de Newton debe ser F~21 = −F~12
~ B,
Lema 4. ∀A, ~ C
~ ∈ R3 se cumple:
~× B
A ~ ×C
~ =B~· A~·C
~ −C~ A~·B
~
Ejercicio 1. Comprobar que las fuerzas dadas por la proposición 16 en la página anterior satisfacen:
F~1←2 = −F~2←1
Podemos ver cómo el segundo sumando en ambas expresiones tiene mismo módulo y sentido contrario, luego
esa parte de las expresiones ya cumple la igualdad buscada. En consecuencia, únicamente tenemos que probar:
˛ ˛ d~l1 d~l2 · ~r12 ˛ ˛ d~l2 d~l1 · ~r12 ˛ ˛ d~l2 d~l1 · ~r21
?
3 = − 3 = 3
L2 L1 r12 L2 L1 r12 L2 L1 r21
Desarrollando el término de la izquierda, obtenemos:
˛ ˛
~ ~ ~r12
dl1 dl2 · 3
L2 L1 r12
Ahora, notamos:
d 1
r12
~ 1 ~r12
− = −∇ =
d~l2 r12 3
r12
Sustituyendo, obtenemos:
˛ ˛ ˛ ˛ d 1 ˛ ˛
~r12 r12 1
d~l1 d~l2 · 3 = −d~l1 d~l2 · =− d~l1 d =0
L2 L1 r12 L2 L1 d~l2 L1 L2 r12
Por otra parte, desarrollando el término de la derecha llegamos al mismo resultado aplicando el mismo pro-
cedimiento. En consecuencia, se cumple F~1←2 = −F~2←1 .
Proposición 48. La fuerza por unidad de longitud que se ejercen entre sí dos conductores rectilíneos paralelos
indefinidos por los que circulan intensidades de corriente eléctrica I1 y I2 , respectivamente, viene dada por la
expresión:
dF µ0 I1 I2
=
dl 2πr
donde r es la distancia que separa ambos conductores.
I1 I2
F~21
B
Demostración. Por la proposición 47 en la página 145, sabemos que la inducción magnética que genera un
conductor rectilíneo indefinido (2) en un punto a distancia r de él viene dada por la expresión:
~ (r) = µ0 I2 ϕ̂
B
2πr
Por otra parte, por el axioma 3 en la página 146, sabemos que la fuerza que sufre un dl del otro conductor
~ es:
rectilíneo indefinido (1) debido a la inducción magnética B
~ = µ0 I1 I2 d~l × ϕ̂
dF~ = I1 d~l × B
2πr
Observación 40. Es del enunciado de la proposición 48 en la página anterior de donde viene la definición de
amperio.
V0
~ 0)
J(R
dS dl
~r
~0
R
~
R
Demostración. Vamos a descomponer nuestro volumen de conductor en infinitos circuitos elementales. Por
la proposición 46 en la página 144, sabemos que el campo magnético generado por cada de dichos circuitos
lineales puede expresarse como:
˛ µ ˛ dId
~ ~ = 0 ¯ ~l0 × ~r
dB R (10.4.0.1)
C0 4π C 0 r3
donde la notación de dI
¯ viene de que nos referimos a la corriente que circula por dicho conductor elemental y
no por todo el conductor
volumétrico. Esa es también la razón por la que en la parte izquierda de la ecuación,
no hemos escrito B~ R ~ ; ya que en este caso dB
~ es un diferencial de orden 3 (un diferencial de volumen) y
sólo hemos integrado a lo largo de una longitud. Tomando la diferencial a ambos lados en la proposición 43
en la página 126, obtenemos:
¯ = J~ · dS
dI ~0
~ 0 k d~l0 , de manera
Por la definición 49 en la página 142, sabemos que en un circuito elemental se cumple J~ k dS
que podemos simplificar la expresión anterior a:
˛ ˛ ˛ ~
~ = µ0 dS 0 dl0 J~ × ~r µ0 J × ~r 0
dB 3
= dV
C0 4π C0 r 4π C 0 r3
Q.E.D.
~ : R3 −→ R3 una función vectorial, ambas diferenciables,
Lema 5. Sean ϕ : R3 −→ R una función escalar y A
entonces:
~ × ϕA
∇ ~ = ϕ∇ ~ ×A ~ − A × ∇ϕ
~
~
siendo A:
˚ J~ R
~0
~ = µ0
A dV 0
4π V0 r
Demostración. Por la proposición 49 en la página 149, tenemos que la inducción magnética generada por un
conductor volumétrico viene dada por:
µ ˚ J~ R ~ 0 × ~r
µ0
˚
~ ~
B R =
0
dV =0 ~
J R~ 0 × ~r dV 0
4π 0 r 3 4π 0 r3
V V
~ 1 ~r
Por la proposición 7 en la página 18, tenemos que ∇ = − 3 , de forma que podemos reescribir la expresión
r r
anterior como: µ ˚
~ ~ 0 ~ ~ 0 ~ 1
B R = −J R × ∇ dV 0 (10.4.0.2)
4π V0 r
1
Por otra parte, por el lema 5 tomando ϕ = , tenemos:
r
J~ R~0
~ 1~ ~ ~ 0 ~ ~ 0 ~ 1
∇× = ∇ × J R −J R × ∇
r r| {z } r
=0
donde el término marcado se anula porque la densidad volumétrica de corriente depende de las coordenadas
pri-
1
~ ~ 0
madas, mientras que el rotacional es respecto a las coordenadas sin primar. Así, despejamos −J R × ∇ ~ ,
r
obteniendo:
J~ R~0
~ ~ 0 ~ 1 ~
−J R × ∇ =∇×
r r
Como la integral es a las coordenadas primadas, mientras que el rotacional según las coordenadas no primadas,
la integral y el rotacional conmutan, por lo que se obtiene:
˚ J~ R ~0
~ R
B ~ =∇ ~ × µ0 dV 0
4π V 0 r
| {z }
~
=:A
Lema 6. Sea F~ : R3 −→ R3 una función vectorial de clase C (2) , entonces la divergencia de su rotacional es
nula.
~ · ∇
∇ ~ × F~ = 0
Corolario 21 (2ª ecuación de Maxwell). La inducción magnética B ~ no tiene fuentes escalares; dicho de otra
~
forma, la divergencia de la inducción magnética B es siempre nula.
~ ·B
∇ ~ =0
para cualquier volumen V . Por el teorema de la divergencia (ver teorema 2 en la página 17), tenemos que:
˚ ‹
0= ~ ~
∇ · BdV = ~ · dS
B ~
V S
y, en consecuencia, se cumple que «las gallinas que entran por las que salen» y el flujo de la inducción
magnética a lo largo de cualquier superficie cerrada es siempre nulo. Por ende, las líneas de campo de la
inducción magnética siempre son cerradas. Q.E.D.
Lema 7. Sea F~ : R3 −→ R3 una función vectorial de clase C (2) . Entonces el rotacional de su rotacional puede
expresarse como:
~ ×∇
∇ ~ × F~ = ∇~ ∇ ~ · F~ − ∇2 F~
Teorema 17 (Teorema de Ampère en forma diferencial). En régimen estacionario de corrientes, siempre que
no haya densidades volumétricas de corriente que entran o salen del volumen estudiado, el rotacional de la
inducción magnética B~ en un punto R~ coincide con la densidad volumétrica de corriente eléctrica J~ en dicho
~
punto R multiplicada por la permeabilidad magnética del vacío µ0 .
~ ×B
∇ ~ R~ = µ0 J~ R ~
Como la integral es a las coordenadas primadas, mientras que la divergencia y la laplaciana son a las coorde-
nadas sin primar, ambos operadores conmutan.
˚ ~ R
J ~0 ˚ J~ R~0
~ ~ ~ ~ µ0 ~ 0 µ0 2 dV 0
∇×B R =∇ ∇· dV − ∇ (10.4.1.1)
4π V 0 r 4π V 0 r
1
Por el lema 2 en la página 51 tomando ϕ = , tenemos:
r
J~ R~0
~ · 1
= ∇ ~ +J~ R
~ · J~ R 0 ~ ·∇
0 ~ 1
∇
r r | {z } r
=0
donde el término marcado se anula porque la divergencia se realiza respecto a las coordenadas sin primar,
pero la densidad volumétrica de corriente depende de las coordenadas primadas. Por tanto:
J~ R~0
~ ~ ~ 0 ~ 1
∇· =J R ·∇
r r
1
De nuevo, por el lema 2 en la página 51 tomando ϕ = , tenemos:
r
J~ R~0
∇~0· = 1∇ ~ 0 · J~ R~ 0 +J~ R ~ 0 1 = J~ R
~0 · ∇ ~0 1
~0 · ∇
r r| {z } r r
=0
donde el término marcado se anula como consecuencia del corolario 19 en la página 129, pues estamos en
régimen estacionario de corrientes. Sustituyendo lo anterior en 10.4.1.2 en la página anterior, obtenemos:
J~ R~0 J~ R~0
~ ·
∇ = −∇ ~0·
r r
Ahora, como la divergencia y la integral son respecto a las mismas coordenadas, podemos aplicar el teorema
de la divergencia (ver teorema 2 en la página 17), obteniendo:
=0
z }| {
˚ ~
J R~ 0 ‹ ~ ~
~ · J R · dS 0
0
∇ dV 0 = − =0
V0 r S0 r
donde el término marcado se anula porque, por hipótesis no hay densidades de corriente eléctrica que entren o
salgan del volumen estudiado. En consecuencia, todo el primer sumando de la expresión 10.4.1.1 en la página
anterior se anula.
Por consiguiente:
˚ J~ R~0
~ ~ ~ µ0 2 dV 0
∇×B R =− ∇ (10.4.1.3)
4π V0 r
J~ R~0
2 . Por la definición 51 en la página 152, tenemos:
Concentrémonos en el término ∇
r
J~ R~0 Jx R ~0 Jy R ~0 Jz R~0
2 2 2 2
∇ = ∇ ,∇ ,∇
r r r r
Como J~ depende de las coordenadas primadas, mientras que el laplaciano actúa con respecto a las coordenadas
sin primar, Jx , Jy y Jz son constantes para el laplaciano. En consecuencia:
J~ R~0
2 ~ 0 2 1 ~ 0 2 1 ~ 0 2 1
∇ = Jx R ∇ , Jy R ∇ , Jz R ∇ =
r r r r
12 ~ 0 ~ 0 ~ 0 2 1 ~ ~0
=∇ Jx R , Jy R , Jz R =∇ J R
r r
Por la proposición 11 en la página 22, podemos expresar lo anterior como:
J~ R~0
∇ 2 = −4πδ (~r) J~ R
~ 0 = −4πδ R ~ −R~ 0 J~ R
~0
r
Q.E.D.
Teorema 18 (Teorema de Ampère en forma integral). En régimen estacionario de corrientes, siempre que no
haya densidades volumétricas de corriente que entran o salen del volumen estudiado, la integral de la inducción
magnética B ~ a lo largo de cualquier curva C es igual a la intensidad de corriente eléctrica IT que atraviesa la
superficie interior de dicha curva multiplicada por la permeabilidad magnética del vacío µ0 .
˛
~ R
B ~ · d~l = µ0 IT
C
IT se considera positiva si circula en la dirección dada como positiva por la regla de la mano derecha al recorrer
la curva.
nI
I
H
~ · d~l = 0
B
C
I I d~l
d~l
Demostración. Por el teorema de Ampère en forma diferencial (ver teorema 18), tenemos:
~ ×B
∇ ~ R~ = µ0 J~ R ~
Por el teorema de Stokes (ver teorema 3 en la página 17 y por la proposición 43 en la página 126, podemos
expresar lo anterior como: ˛
~ R
B ~ · d~l = µ0 IT
C
donde C es la curva que delimita la superficie S. Q.E.D.
Capítulo 11
Proposición 51.
Demostración. Los resultados se siguen trivialmente al sustituir la expresión dada para la inducción magnética
~ en la proposición 50 en la página 151 en el corolario 21 en la página 152 y en el teorema 17 en la página 153.
B
Q.E.D.
Lema 8. Sea ψ : R3 −→ R una función escalar de clase C (2) , entonces la el rotacional de su gradiente es
nulo.
~ × ∇ψ
∇ ~ = ~0
~ como:
Supongamos que podemos expresar A
~=A
A ~ 0 + ∇ψ
~
donde A~ 0 es una función vectorial y ψ es una función escalar. Como el rotacional de un gradiente es siempre
nulo por el lema 8 en la página anterior, tenemos:
~ ×A
∇ ~=∇
~ ×A
~0 + ∇
~ × ∇ψ
~ =∇
~ ×A
~0
| {z }
=0
~=
En consecuencia, existen A ~ 0 tales que:
6 A
~ =∇
B ~ ×A
~=∇
~ ×A
~0
Demostración.
Partimos de la proposición 49 en la página 149, que nos da la inducción magnética generada por cualquier
distribución volumétrica de corrientes.
µ ˚ J~ R ~ 0 × ~r
~ R
B ~ = 0 dV 0
4π 0 r 3
V
Por último, si contamos con una corriente que es lineal, la densidad de corriente J~C es paralela a d~l en
todo punto y su módulo es la intensidad de corriente I. Además, J~C se anula en todo punto que no
pertenece a la curva cerrada C. En consecuencia, podemos simplificar la expresión a:
µ ˛ Id~l × ~r
~ R
B ~ = 0
4π C 0 r3
Para concluir, por el principio de superposición (ver proposición 1 en la página 6), tenemos que la
inducción magnética generada en cualquier punto R ~ podrá expresarse como la suma de las inducciones
magnéticas generadas por cada tipo de distribuciones de corrientes, con lo que llegamos al enunciado.
Q.E.D.
(0, 0, z)
~r
z=0
a
B~ (0, 0, z) = µ0 I sen3 α
2a
ˆ ˛
µ0 I adϕ µ0 I
~
A (0, 0, z) = ϕ̂ = d~l = ~0
4π r 4πr C
Ejemplo 29 (Problema 61). Tenemos un anillo:
ω
dl = R0 · dθ
θ
R0
dQ dl
I= = λ = λR0 ω
dt dt
dl = R0 dθ = R0 ωdt
Por el axioma 2 en la página 143, tenemos:
~ = µ0 I d~l × ~r
dB
4π r3
Sustituyendo, obtenemos:
Integrando:
N, I
L
α2 α α1 Z
dz 0
z=0 z=z
z0
z − z0
N
n :=
L
N 0
dz = ndz 0
L
~ = µ0
dB ndz 0 I sen3 α ẑ
2a
µ0 I na2
dBz (0, 0, z) = h i dz 0 ⇔
2 a2 + (z − z 0 )2
ˆ L
µ0 nI a2 0
⇔ Bz = i 3 dz
z 0 =0 2 h 2
a + (z − z 0 )2
2
adα
⇒ −dz 0 = −
sen2 α
En consecuencia, queda:
ˆ ˆ
µ0 nI 3 0 µ0 nI α2 sen3 α a µ0 nI
Bz = sen αdz = dα = [cos α1 − cos α2 ] =
2a 2a α1 sen α 2
µ0 nI z z−L
= √ −q
2 2
a +z 2
2 2
a + (z − L)
¿Qué pasa si el solenoide es infinito? Entonces α1 → 0 y α2 → π. En consecuencia:
µ0 nI
Bz = [1 − (−1)] = µ0 nI
2
Demostración. Por el teorema de Ampère en forma diferencial (ver teorema 17 en la página 153), tenemos
que es:
~ ×B
∇ ~ R~ = µ0 J~ R ~
~ ~ ~ ~
Como en G es ∇ × B = 0, entonces existe una función escalar φm : R −→ R de clase C tal que B R
3 (1) ~ =
~ m R
−∇φ ~ ∀R ~ ∈ G. Q.E.D.
Proposición 55. El potencial magnético escalar φm generado por un circuito elemental por el que circula
una corriente I viene dado por:
φm R ~ = − µ0 I Ω
4π
donde Ω es el ángulo sólido con el que se ve el área delimitada por el circuito lineal estacionario desde el punto
~
R.
dl
I ~r
−d~r
Demostración. En todo punto ajeno al circuito elemental es J~ = ~0, luego, por la proposición, existe una
~ = −∇φ
función escalar φm : R3 −→ R de clase C (1) tal que B ~ m en todo punto del espacio salvo en el circuito
lineal estacionario.
Por la proposición 46 en la página 144, sabemos que la inducción magnética generada por un circuito
elemental viene dada por la expresión:
µ I ˛ d~l × ~r
~
B R~ = 0
4π L r3
Por otra parte, por el teorema del gradiente, tenemos:
ˆ
~ = −∇φ
B ~ m ⇔ φm = − ~ · d~r + C
B
donde C es una constante, que podemos suponer nula sin pérdida de generalidad. Así, sustituyendo, obtenemos
ˆ " ˛ ~ # ˆ ˛ ~
µ0 I dl × ~r µ0 I dl × ~r
φm = − · d~r = − · d~r
4π L r3 4π L r3
Como en el producto triple un cambio circular de los factores no altera el resultado1 , obtenemos:
ˆ ˛ ˆ ˛
µ0 I d~r × d~l µ0 I −d~r × d~l
φm =− · ~
r = · ~r
4π L r3 4π L r3
~ := −d~r × d~l, llegamos a:
Llamando dS
ˆ ˛ ~ · ~r
µ0 I dS
φm =
4π L r3
1
Es una de las propiedades del producto triple. Puede verse en el siguiente enlace: https://en.wikipedia.org/w/index.php?
title=Triple_product&oldid=877202524#Properties
Observación 41. Nótese que por la proposición 55 en la página anterior, la diferencia de potencial magnético
entre dos puntos viene dada por:
φm R~ 2 − φm R ~ 1 = − µ0 I Ω ~ − Ω ~
4π R2 R1
µ I ~r · S
~
~ ≈ 0
φm R
4π r 3
~ −R
donde ~r = R ~ 0 e I es la intensidad que circula por la espira.
~
dS
I
a
Demostración. Por la proposición 55 en la página anterior, sabemos que el potencial magnético generado por
un circuito elemental viene dado por:
φm R ~ = − µ0 I Ω
4π
Diferenciando a ambos lados, obtenemos:
µ0 I
dφm = − dΩ
4π
Como S ~ |~r|, podemos considerar que ~r no varía a lo largo de la superficie y, por tanto, ~r sale fuera de la
integral.
¨ ~
µ0 I ~r ~ µ0 I ~r · S
φm ≈ 3
· dS = 3
4π r S 4π r
Q.E.D.
Definición 53. Llamamos momento dipolar magnético m ~ de una espira al vector superficie de dicha
espira multiplicado por la intensidad que circula por ella.
m ~
~ := I S
Proposición 57. El potencial magnético escalar generado por una espira por la que circula una corriente
I tal que su área es muy pequeña en comparación con la distancia al punto en el que se está calculando el
potencial, siendo m
~ el momento dipolar magnético dipolar, puede aproximarse como:
µ m
φm R ~ ≈ 0 ~ · ~r
4π r3
~ −R
donde ~r = R ~ 0.
donde:
~ ~ ~
~ ·∇
C ~ D ~ = Cx ∂ D + Cy ∂ D + Cz ∂ D
∂x ∂y ∂z
Proposición 58. La inducción magnética generada por una espira de momento dipolar m
~ y de área muy
~
pequeña en comparación con la distancia al punto R puede aproximarse por:
µ 3 (m 2~
~ R
B ~ ≈ 0 ~ · ~r) ~r − r m
4π r5
Demostración. Por la proposición 57, tenemos que el potencial magnético escalar generado por una espira que
cumple las condiciones del enunciado puede aproximarse como:
µ m
φm R ~ ≈ 0 ~ · ~r
4π r3
Por la proposición 54 en la página 160, sabemos:
B ~ m ≈ − µ0 ∇
~ = −∇φ ~ m ~r
~ · 3 (11.3.0.1)
4π r
Análogamente:
~r !
∂ r3 −3ry rx r2 − 3ry2 −3ry rz
= , ,
∂y r5 r5 r5
∂ ~r
r3 −3rz rx −3rz ry r2 − 3rz2
= , ,
∂z r5 r5 r5
Así:
~ m ~
r mx r2 − 3mx rx2 − 3my rx ry − 3mz rx rz my r2 − 3mx ry rx − 3my ry2 − 3mz ry rz
∇ ~ · 3 = , ,
r r5 r5
mz r2 − 3mx rz rx − 3my rz ry − 3mz rz2
,
r5
r2 r2 m
~ − 3m
~ · ~r
= (m x , my , mz ) − 3 m
~ · ~
r (rx , ry , rz ) =
r5 r5
Sustituyendo en la ecuación 11.3.0.1 en la página anterior, obtenemos:
~ · ~r) ~r − r2 m
~ ≈ µ0 3 (m ~
B 5
4π r
Q.E.D.
Proposición 59. El potencial magnético vector A~ generado por una espira de momento dipolar m
~ y de área
~ puede aproximarse por:
muy pequeña en comparación con la distancia al punto R
µ m
A~ R~ ≈ 0 ~ × ~r
4π r3
~ −R
donde ~r = R ~ 0.
~ · ~r
µ0 m µ0 ~r
φm ≈ 3
= ~ · 3
m
4π r 4π r
Por la proposición 7 en la página 18, tenemos:
µ0 ~ 1
φm ≈− m ~ ·∇ (11.3.0.2)
4π r
1
Por otra parte, por el lema 2 en la página 51 tomando ϕ = , tenemos:
r
~ m
~ 1~ ~ 1
∇· = ∇ | ·m
{z ~ ·∇
~} +m
r r r
=~0
−
→
donde el término marcado se anula porque m ~ = cte.
Así, sustituyendo en la ecuación 11.3.0.2, obtenemos:
µ0 ~ m
~
φm ≈− ∇ ·
4π r
En consecuencia, por medio de la ecuación 11.3.0.4, podemos expresar 11.3.0.3 como sigue:
~ ≈ µ0 ∇
B ~ × ∇ ~ ×m ~
=∇~ × µ0 ∇~ ×m ~
4π r 4π r
~ ≈ µ0 ∇
A ~ ×m
~
(11.3.0.5)
4π r
1
Por último, por el lema 5 en la página 151 tomando ϕ = , obtenemos:
r
~ m~ 1~ ~ 1
∇× = ∇ | ×m
{z ~} −m
~ ×∇
r r r
=~0
−→
donde el término marcado se anula porque m
~ = cte. Por la proposición 7 en la página 18, tenemos:
~ m
~ ~r m~ × ~r
∇× =m ~ × 3 =
r r r3
~ × ~r
~ ≈ µ0 m
A
4π r3
Q.E.D.
Proposición 60. Una distribución de dipolos de volumen V 0 puede ser modelada mediante una serie de
densidades de corriente equivalentes (llamadas de imanación) que vienen dadas por:
J~m R~0 = ∇~0×M ~ R ~0
J~S,m R
~0 = M
~ R~ 0 × n̂0
V0
~r
~0
R
~
R
Demostración. Haciendo la diferencial a ambos lados de la proposición 59 en la página 164, obtenemos que el
potencial magnético vector generado por un dV 0 de distribución de dipolos viene dado por:
µ dm
~ R
dA ~ = 0 ~ × ~r (11.4.0.1)
4π r3
Por otra parte, por la definición 54 en la página anterior, tenemos:
~ 0 = dm
M R
~
⇔ d m
~ = M ~ dV 0
R
dV 0
Sustituyendo en 11.4.0.1, obtenemos:
µ M~ R~ 0 × ~r
~ ~ 0
dA R = dV 0 ⇔
4π r3
µ ˚ M ~ R
~ 0 × ~r
⇔A~ R~ = 0 dV (11.4.0.2)
4π V0 r3
1
Por otra parte, por el lema 5 en la página 151, tomando ϕ = , obtenemos:
r
!
~
M
~ 0 ~ 0 1 ~ + 1∇~0×M~
∇ × =∇ ×M
r r r
~ ~ ~ ~
~ 0 × M = ~r × M + 1 ∇
∇ ~ ⇔ M × ~r = −∇
~0×M ~ 0 × M + 1∇~0×M
~
r r 3 r r 3 r r
JS,m
θ
µ n̂
dm ~ dV 0 = M0 ẑdV 0 ⇔
~ =M
˚
4
⇔m
~ = M0 ẑdV 0 = πR03 M0 ẑ
V0 3
Obtenemos las corrientes equivalentes:
~ ×M
J~m = ∇ ~ = ~0
J~S,m = M
~ × n̂ = M0 sen θ ϕ̂
B B
π
θ− 2
θ
~ ~r
dB
~ = µ0 J~S,m × ~r
dB dS
4π r3
µ0 JS,m ~ =
dB
dS
4π r2
Por el dibujo podemos ver que las las componentes que no están en el eje z se anulan entre sí; de manera
que:
=JS,m
z }| {
µ0 M0 sen θ π
dBz = cos − θ R02 sen θdϕ ⇔
4π R02 2
ˆ θ=π ˆ
µ0 M 0 3 µ0 M 0 π 2
⇔ Bz = 2π sen θdθ = sen θ sen
| {z θ} dθ =
4π θ=0 2 0
=1−cos2 θ
ˆ π ˆ π 3 π
µ 0 M0 2 µ0 M 0 π cos θ
sen θdθ − cos θ sen θdθ = [− cos θ]0 − =
2 0 0 2 3 0
µ0 M 0 2 µ0 M 0 4 2
= 2− = = µ 0 M0
2 3 2 3 3
Capítulo 12
V0
~ = cte
M
Consideramos que la inducción magnética B ~ R ~ que hay en un punto R ~ es la suma de la inducción
~ ext y la inducción magnética generada por los átomos más próximos B
magnética externa B ~ 0 . Esto es sólo un
modelo (los hay mejores).
~ =B
B ~ ext + B
~0
Para hallar el campo B ~ 0 , vamos a suponer que tenemos un agujero esférico en nuestro medio natural. La idea
es que la imanación del medio va a generar unas corrientes equivalentes que, a su vez, generarán inducción
magnética. Supongamos que el mencionado agujero tiene radio R0 y es tal que J~m = ~0. Es decir, la densidad
volumétrica equivalente de corriente se anula en el agujero esférico. Sin embargo, sí que tendremos una densidad
de corriente superficial equivalente J~S,m sobre la superficie interior del agujero esférico.
θ
~
M
µ0
n̂
J~S.m = M
~ × n̂ = M · 1 · sen (π − θ) (−ϕ̂) = −M sen θϕ̂
~ θ
M
π
θ− 2
µ0 θ n̂
R
~
dB
µ0 M sen θ π µ0
dBz = cos − θ R02 sen θdθdϕ = M sen3 θdθdϕ ⇔
4π R02 | 2
{z } 4π
=sen θ
2 ~ 0 = − 2 µ0 M
~ ẑ
⇔ Bz = µ0 M ⇒ B
3 3
En consecuencia, la inducción magnética total que hay en un punto R ~ viene dada como:
B ~ ext + − 2 µ0 M
~ =B ~
3
El resultado obtenido no es más que un caso particular de la llamada ecuación de Weiss, que sirve para
cualquier geometría (cualquier tipo de agujero):
~ =B
B ~ ext + µ0 (γ − 1) M
~
12.2. ~
La excitación magnética H
Proposición 61. En un medio con imanación magnética M ~ 6= ~0, el teorema de Ampère para la inducción
~
magnética B puede escribirse como sigue:
h i
~ ×B
∇ ~ R ~ = µ0 J~C R
~ +∇~ ×M ~ R ~
donde J~C es la densidad de corriente volumétrica de las llamadas corrientes de conducción o corrientes reales
y∇~ ×M ~ = J~m es la densidad de corriente volumétrica de las corrientes equivalentes de imanación.
~
J~sm
M
J~ J~m
J~
V0
Demostración. Partimos del teorema de Ampère en forma diferencial (ver teorema 17 en la página 153):
~ ×B
∇ ~ = µ0 J~ (12.2.0.1)
~ 6= ~0, la densidad de corriente total puede expresarse como suma de dos
En un medio con imanación, M
términos:
J~ = J~C + J~m (12.2.0.2)
donde J~C es la densidad de corriente volumétrica asociada a las corrientes reales o libres y J~m es la densidad
volumétrica de corriente equivalente de imanación. Por la proposición 60 en la página 166, tenemos:
J~m = ∇
~ ×M
~
~
~ := B − M
H ~
µ0
~ en función de la excitación
Observación 43. Nótese que, según la definición 55, la inducción magnética B
~
magnética Hviene dada como:
B~ = µ0 H ~ +M ~
Teorema 19 (Teorema de Ampère para H). ~ En régimen estacionario de corrientes, siempre que no haya
densidades volumétricas de corriente que entran o salen del volumen estudiado, el rotacional de la excitación
magnética H~ en un punto R
~ coincide con la densidad volumétrica de corriente eléctrica real J~C en dicho punto
~
R.
∇~ ×H ~ R~ = J~C R ~
~ ×B
∇ ~
⇔ ~ ×M
−∇ ~ = J~C
µ0
Como el rotacional es un operador lineal, obtenemos:
!
~
B
~ ×
∇ ~
−M = J~C
µ0
Ahora, integrando la igualdad anterior a una superficie abierta S cualquiera a ambos lados, tenemos:
¨ ¨
~ ~ ~
∇ × H R · dS = ~ J~C R~ · dS~
S S
Por el teorema de Stokes (ver teorema 3 en la página 17 y por la proposición 43 en la página 126, podemos
expresar lo anterior como: ˛
~ R
H ~ · d~l = IC
C
donde C es la curva que delimita la superficie S. Q.E.D.
~ = ~0 en un volumen V ⊆ R3 , entonces existe una función escalar ψm : R3 −→ R
Proposición 62. Si es J~C R
de clase C (2) tal que:
~ =∇
∇2 ψm R ~ ·M
~ R~ = −∇
~ ·H
~ R~ ~ ∈V
∀R
∇ ~ = 1∇
~ ·H ~ ·B
~ −∇
~ ·M
~
µ0
~ ·B
pues la divergencia es un operador lineal. Por el corolario 21 en la página 152, es ∇ ~ = ~0 ∀R
~ ∈ R3 y, en
~ ·B
particular, es ∇ ~ = ~0 ∀R
~ ∈ V . Por tanto, obtenemos:
~ ·H
∇ ~ = −∇
~ ·M
~ (12.2.0.3)
~ = ~0 ∀R
Por otra parte, por el teorema 19 en la página anterior, si es J~C R ~ ∈ V , será:
~ ×H
∇ ~ R~ = ~0 ∀R
~ ∈V
~ ×H
y, en consecuencia, sabemos que existe una función escalar ψm : R3 −→ R de clase C (2) (para que ∇ ~ sea
continua) tal que:
~ R
H ~ = −∇ψ ~ m R ~ ∀R ~ ∈V
Q.E.D.
Definición 56. Si se dan las condiciones para aplicar la proposición 62, llamaremos densidad volumétrica
de masa magnética ρm a la función escalar:
ρm R ~ = −∇ ~ ·M ~ R~
Bien, ahora supongamos que en nuestro volumen tenemos únicamente una masa magnética mm puntual (algo
que no tiene sentido físico estrictamente) e integremos la expresión anterior al volumen V :
˚ ˚
~ ~ ~
∇ · H R dV = ρm R ~ dV = mm
V V
La fórmula anterior debe valer para cualquier volumen; luego, en particular, podemos escoger una esfera cuyo
centro esté en la posición de mm . Como la excitación magnética generada por una masa magnética puntual
no puede depender de la dirección (una masa magnética puntual tiene simetría esférica), sabemos que H ~ R ~
~ y además, será H
será constante a lo largo de dS ~ R~ k dS ~ en todo momento. Así, obtenemos:
~ = mm
~ 4πr2 = mm ⇔ H R
H R
4πr2
Como habíamos dicho que la excitación magnética generada una masa magnética puntual debe tener simetría
esférica, debe ser:
~ R
H ~ = 1 mm r̂
4π r2
A continuación, vemos que si tomamos:
1 mm
ψm =
4π r
entonces:
~ m=−d
−∇ψ
1 mm
=−
∂ 1 mm
r̂ =
1 mm
r̂ = ~ R
H ~
d~r 4π r ∂r 4π r 4π r2
Luego, ya hemos encontrado el potencial magnético generado por una «masa magnética puntual». Nótese que
la expresión es análoga a la obtenida en electroestática.
Siguiendo los mismos pasos que en la proposición 16 en la página 29, llegamos a que el potencial magnético
generado por una distribución volumétrica de masas magnéticas es:
˚ ρm R ~0
~ 1
ψm R = dV 0
4π V 0 r
Si nos restringimos a una superficie, entonces podemos expresar el potencial anterior como:
¨ σm R ~0
~ = 1
ψm R dS 0
4π S 0 r
Así, el potencial magnético generado por cualquier distribución volumétrica y superficial de masas magnéticas
puede expresarse como:
˚ ρm R ~0 ¨ σm R ~0
ψm R ~ = 1 dV + dS 0
4π V0 r S0 r
~ = −∇ψ
Para obtener las expresiones para la excitación magnética, simplemente aplicamos H ~ m , obteniendo:
˚ ¨
~ ~ 1 ρm~r 0 σm~r 0
H R = 3
dV + 3
dS
4π V0 r S0 r
Q.E.D.
~ = M0 ẑ.
Ejemplo 32. Tenemos un imán de forma esférica con una imanación uniforme con M
µ n̂
~ ·M
ρm = −∇ ~ =0
~ · n̂ = M0 cos θ
σm = M
Con esto, podemos obtener el diferencial de excitación magnética:
~ = 1 σm~r 0
dH dS
4π r3
~ = 1 M0 cos θ 2
dH R0 sen θdθdϕ
4π R02
1 M0 cos θ
dHz = cos θR02 sen θdθdϕ
4π R02
~ = Hz (−ẑ). Integrando, obtenemos:
siendo H
ˆ π π
1 2 1 cos3 θ 1 M0
Hz = M0 cos θ sen θdθ = M0 − = M0 (1 − (−1)) =
2 0 2 3 0 6 3
Así:
~ = − M0 ẑ
H
3
nótese que [H] = A
m = [M ].
Observación 44. Por la definición de excitación magnética (ver definición 55 en la página 174) al aplicarle la
definición 57, obtenemos:
~ ~ ~
M B ~ ~ 1 B ~ = 1 χm B ~
:= −M ⇔M +1 = ⇔M
χm µ0 χm µ0 µ0 1 + χm
Observación 45. Por la observación 43 en la página 174, al aplicar la definición 57, llegamos a:
~ = µ0 H
B ~ + χm H~ = µ0 (1 + χm ) H ~
~ definimos:
Análogo a lo visto en electroestática y para la densidad volumétrica de corriente E,
~ y la excitación magnética H
Definición 58. Llamamos medio isótropo a aquel en el que la imanación M ~
son paralelos en todo momento:
~ kH
M ~
Definición 60. Llamamos medios homogéneos a aquellos en los que la susceptibilidad magnética no depende
de la posición:
χm 6= F R~
Análogamente, llamamos medios no homogéneos a aquellos en los que χm = F R ~ .
~ = µH
B ~
µ = µ0 (1 + χm )
Definición 62. Llamamos permeabilidad magnética relativa µr de un medio al cociente entre la per-
meabilidad magnética de dicho medio µ entre la permeabilidad magnética del vacío µ0 .
µ
µr := = (1 + χm )
µ0
Definición 63. Decimos que un medio magnético es
Ferromagnéticos: El cobalto: 2, 5 · 102 , el níquel: 2 · 102 , el acero: 2 · 103 , el permalloy78 : 105 , el hierro
(< 0, 05 % de impurezas): 2 · 105 o el sumpermalloy: 106 .
Ms
Mr
Hc Hc H
Mr
−Ms
La idea es la siguiente, yo tengo un material no imanado M = 0 al que le aplico una excitación magnética
~ hasta que el material no es capaz de aumentar más su imanación (1), llegando a la imanación de saturación
H
Ms . Si ahora, disminuyo la excitación magnética hasta que H ~ = ~0 (2), entonces la magnetización no es nula,
sino que tiene cierto valor Mr llamado magnetización remanente. A continuación, aplicamos el campo H en
sentido contrario, de manera que la magnetización disminuye hasta que cambia de signo (3) y alcanza el
valor −Ms . El corte con el eje de abscisas se llama campo coercitivo o −HC . Ahora, disminuimos el campo
magnético a cero (4), nuevamente y la magnetización cuando H = 0 es −Mr . Por último, volvemos aumentar
el campo magnético (5) y, en consecuencia, la magnetización aumenta hasta que cambia de signo cuando
H = HC y, de nuevo, alcanzamos la magnetización de saturación Ms . Esto se llama el ciclo de histéresis.
En estos diagramas, puede obtenerse la susceptibilidad magnética como:
dM
χm (H) =
dH
Cada material tiene un ciclo de histéresis característica. Veamos unos ejemplos:
B B
Aproximadamente lineal
H H
B2n = B1n
Demostración. Tomamos un cilindro tal que sus tapas son paralelas a la superficie de contacto entre ambos
medios magnéticos de forma que una de sus tapas esté en un medio y la otra en el otro; es decir, una «pillbox».
Gráficamente:
(2) S2
~2
B
µ2
µ1
(1) ~1
B
Por otra parte, por el corolario 21 en la página 152, sabemos que la divergencia de la inducción magnética
~
B es nula en todo punto:
~ ·B
∇ ~ =0
donde S es la superficie cerrada que delimita el cilindro. Ahora bien, podemos descomponer dicha integral en:
‹ ¨ ¨ ¨
0= ~ · dS
B ~= B~ · dS
~+ ~ · dS
B ~+ B~ · dS
~
S tapa superior tapa inferior superficie lateral
~ 1 · dS
B ~1 = B
~ 1 (−n̂) dS1 = −B1n dS1
Q.E.D.
Corolario 22. En la frontera entre dos medios magnéticos de permeabilidades magnéticas µ1 y µ2 se da:
µ2 H2n = µ1 H1n
Proposición 65. En régimen estacionario de corrientes, siempre que no haya densidades volumétricas de
corriente que entran o salen del volumen estudiado, en la frontera entre dos medios magnéticos la resta de las
componentes tangenciales de la excitación magnética es igual a la densidad superficial de corriente eléctrica
libre (de conducción) que atraviesa la superficie de contacto entre ambos medios en dirección perpendicular a
la dirección tangencial tomada.
Demostración. Tomamos un rectángulo tal que dos de sus lados son paralelos a la superficie de contacto entre
ambos medios magnéticos. Gráficamente:
h→0
~2
B µ2
(2) t̂
t
(1) µ1
~1
B
~ en forma integral (ver teorema 19 en la página 174), tenemos:
Por el teorema de Ampère para el campo H
˛
H~ · d~l = IC
C
donde n̂ es el vector unitario normal a la superficie de contacto entre ambos medios que apunta hacia el medio
dos. Si hacemos tender la altura del rectángulo a cero (h → 0), entonces nuestra expresión queda:
˛ ˆ ˆ
~ ~ ~
H · dl = H1 · t̂dl − H ~ 2 · t̂dl = IC,x̂
C 1 2
donde IC,x̂ es la corriente de conducción que circula por la frontera entre ambos medios y x̂ = t̂ × n̂. Como en
un rectángulo dos lados paralelos miden lo mismo (llamemos a dicha longitud L), obtenemos:
˛ ˆ L ˆ L
~ · d~l =
H H1t dl − H2t dl = IC,x̂ ⇔
C 0 0
IC,x̂
(H1t − H2t ) L = Ix̂ ⇔ H1t − H2t =
L
Por el corolario 18 en la página 127, podemos expresar el término de la derecha como JS,C,x̂ . Así, llegamos al
enunciado. Q.E.D.
Corolario 23. En régimen estacionario de corrientes, siempre que no haya densidades volumétricas de co-
rriente que entran o salen del volumen estudiado, en la frontera entre dos medios magnéticos se cumple:
B1t B2
− t = JS,C,x̂
µ1 µ2
donde JS,x̂ es la densidad de corriente superficial que circula por la frontera entre ambos medios magnéticos
en la dirección x̂, siendo x̂ = t̂ × n̂ donde n̂ es el vector normal a la superficie de contacto entre ambos medios
que apunta hacia el medio dos.
Corolario 24. En régimen estacionario de corrientes, siempre que no haya densidades volumétricas de co-
rriente que entran o salen del volumen estudiado, si no circula corriente eléctrica en la frontera entre dos
medios magnéticos, entonces la componente tangencial de la excitación magnética H ~ se conserva.
B2t B1t
H2t = H1t ⇔ =
µ2 µ1
Demostración. La demostración se sigue trivialmente de la proposición 65 en la página 183 y del corolario 23
al sustituir JS,C,x̂ = 0. Q.E.D.
Proposición 66 (Ley de Ohm para el circuito magnético). Supongamos que nos encontramos en régimen
estacionario de corrientes y, además, que no hay densidades volumétricas de corriente que entren o salgan
del volumen estudiado. Por añadidura, tenemos un medio magnético con forma de cable, por el que circula un
flujo magnético Φm . Si en cada sección transversal del medio (transversal con respecto a la dirección de Φm )
se cumple que la inducción magnética B ~ es constante y, al tomar una curva cerrada C que atraviesa el centro
del medio magnético en todo momento, B ~ es siempre paralela a dicha curva C; entonces la intensidad libre
IC que atraviesa la superficie abierta definida por dicha curva es igual al producto del flujo magnético Φ que
circula por el medio y una cierta magnitud R.
m = Φm R
donde: ˛
m= ~ · d~l = IC
H
C
y: ˛
dl
R=
C µS
siendo µ la permeabilidad magnética del medio y S la superficie transversal del cable.
I
I
n n
Demostración. Por el teorema de Ampère para la excitación magnética en forma integral (ver teorema 19 en
la página 174) aplicado a la curva cerrada C que recorre el centro del medio magnético, sabemos:
˛
H~ · d~l = IC = m
C
Por otra parte, por la definición 61 en la página 179, tenemos:
˛ ˛ ~
~ · d~l = B ~
H · dl
C C µ
siendo µ la permeabilidad magnética del medio. Si la inducción magnética B ~ es paralela a la curva C en todo
momento, entonces: ˛ ˛
~ · d~l = B dl
H (12.6.0.1)
C C µ
Por otra parte, siendo S una superficie transversal cualquiera del medio, por la definición 64 en la página
anterior, tenemos: ¨
Φm = ~ · dS
B ~
S
~ k d~l. Como B
donde dS ~ es constante a lo largo de toda superficie transversal del medio y B
~ es paralelo a la
curva C en todo momento, lo anterior se simplifica a:
Φm
Φm = BS ⇔ B =
S
Sustituyendo en la ecuación 12.6.0.1, obtenemos:
˛ ˛ ˛
~ ~ Φm dl dl
m= H · dl = = Φm = Φm R
C S C µ C µS
Q.E.D.
Definición 65. Llamaremos fuerza magnetomotriz al factor m que aparece en la proposición 66 en la
página anterior.
Definición 66. Llamaremos reluctancia a la magnitud R que aparece en la proposición 66 en la página
anterior. ˛
dl
R :=
C µS
siempre que se cumplan las condiciones de la mencionada proposición.
12.7. Ejemplos
Ejemplo 33 (Problema 76). Tenemos:
R1 Rx
µ S1 S2
I
Φ
l1 l2 NI = m R2
N
Rx R1
x
Tenemos, por la definición 66 en la página anterior:
l1 x l2
R1 = Rx = R2 =
2µS1 µ0 S1 µS2
Y, por la proposición 66 en la página 185, obtenemos:
NI
N I = m = Φ (R1 + Rx + R2 + Rx + R1 ) ⇔ Φ = l1 2x l2
µS1 + µ0 S1 + µS2
Este sistema es un relé que, en la práctica, sirve para trasladar un interruptor de un sitio donde podría
ser peligroso a otro donde no lo es tanto.
Ejemplo 34 (Problema 75). Tenemos un imán esférico con una imanación M ~ = (M0 + M1 z) ẑ.
~
M n̂
θ R0
Para resolver este problema o bien aplicamos corrientes equivalentes, o bien aplicamos las masas magnéti-
cas. Lo vamos a hacer de ambas formas. Empezamos con corrientes equivalentes. Por la proposición 60 en la
página 166, tenemos:
J~S,m = M
~ × n̂ = M0 + M1 , z · sen θϕ̂
|{z}
=R0 cos θ
J~m = ∇
~ ×M
~ = ~0
~
Calculemos la inducción magnética B.
B B
π
θ− 2
θ
~ ~r
dB
~ = µ0 JS
dB dS
4π r2
µ0 JS π
dBz = cos − θ R02 sen θdθdϕ ⇔
4π R02 2
ˆ ˆ ˆ
µ0 π µ0
π π
=
⇔ Bz = 2π (M0 + M1 R0 cos θ) sen3 θdθ = M sen 3
θdθ + M R cos θ sen 3
θdθ
2 0
0 1 0
4π θ=0 0
| {z }
=0
2 ~ = 2 µ0 M0 ẑ
= µ 0 M0 ⇒ B
3 3
Ahora, resolvamos el problema por masas magnéticas. Por la definición 56 en la página 175, tenemos:
~ ·M
ρm = −∇ ~ = −M1
σm = M cos θ
A continuación, por la proposición 63 en la página 175, sabemos:
¨ ˚
H~ = 1 σm~r
dS +
ρm~r
dV
4π 3 3
S r V r
Vamos con la componente volumétrica:
~V = 1 ρm~r
dH dV
4π r3
1 −M1
dHz = cos θr2 sen θdrdθdϕ
4π r2
ˆ 2π ˆ π ˆ R0 π
−M1 M1 sen2 θ
−Hz = cos θ sen θdrdθdϕ = − 2πR0 =0
4π ϕ=0 0 0 4 2 0
| {z }
=0
Así:
~ = 1 M cos θ 2
dH R0 sen θdθdϕ
r2
4π |{z}
=R02
ˆ
1 1 π
−dHz = M cos θ sen θdθdϕ ⇒ −Hz = (M0 + M1 R0 cos θ) cos2 θ sen θdθ =
4π 2 θ=0
π π
1 cos3 θ cos4 θ 1 2 M0
= M0 − + M1 R0 − = M0 + 0 = ẑ
2 3 0 4 0 2 3 3
En consecuencia:
H~ = − M0 ẑ
3
Por la observación 43 en la página 174, tenemos:
~ = µ0 − M 0 2
B ẑ + M0 ẑ = µ0 M0 ẑ
3 3
Imaginemos que ahora hacemos un agujero cilíndrico en la esfera.
θ0
Teníamos:
~ = (M0 + M1 z) ẑ
M
J~m = ~0
Por el método de corrientes equivalentes, por la proposición 60 en la página 166, tendríamos en cilíndricas:
J~S,m = M
~ × n̂ = (M0 + M1 z) (−ϕ̂)
No obstante, vamos a hacerlo mediante masas magnéticas. De nuevo, por la definición 56 en la página 175:
ρm = −M1
σm = (M0 + M1 z) cos θ
Nótese que:
1 M0 M0
−2 cos3 θ0 ẑ = −
= cos3 θ0 ẑ
2 3 3
Por la observación 43 en la página 174, tenemos:
Capítulo 13
Inducción electromagnética
∂Φm
ε=−
∂t
Ejemplo 35 (Dinamo). Tenemos una espira que rota con velocidad constante ω dentro de un campo magné-
tico.
~
B
N S
~ ·S
B ~ = BS cos ωt
¨
Φm = B cos ωtdS = BS cos ωt
S
Así, por la ley de inducción de Faraday-Lenz (ver axioma 4), tenemos:
∂Φm
V (t) := ε (t) − = BS sen ωt
∂t
Observación 46. Podemos entender la ley de Faraday-Lenz de forma intuitiva de la siguiente forma: a una
bobina «no le gusta» que cambie el flujo magnético que la atraviesa y se opone a cualquier cambio, generando
para ello un flujo. Pero para generar un flujo, es necesario que circule una intensidad.
~ = µnI (t)
B
I1
N1 N2
I2
N2
V2 = V1 = nV1
N1
Ejemplo 37. Tenemos dos masas una de aluminio y otra de neodimio (imanado) y las dejamos caer por el
interior de un tubo de latón.
F~m
p~
La masa de aluminio caerá a través del tubo de latón como si el tubo no estuviera. Sin embargo, mientras
cae la masa de neodimio imanada, la inducción magnética en el interior del cilindro variará. De hecho, si
vemos el cilindro como un montón de espiras circulares apiladas una encima de otra, conforme caiga la masa
imanada, el flujo magnético a través de ellas cambiará. Al cambiar, por el axioma 4 en la página 192, se creará
una corriente inducida en ellas; dicha corriente generará una inducción magnética en el interior del cilindro
que se opondrá a la inducción magnética generada por el imán. Dicha inducción generará una fuerza sobre la
masa imanada que se opone a su caída. Por tanto, la masa imanada caerá más lentamente que la masa sin
imanar.
V(t) (V)
300
200
100
t (s)
0.0100 0.0200 0.0300 0.0400 0.0500 0.0600 0.0700 0.0800
−100
−200
−300
Si conectamos un aparato de resistencia R, entonces, por la ley de Ohm, la intensidad que circula vendrá
dada por:
√
V 220 2
I= = cos (2π50t)
R R
Estudiemos la potencia:
V (t)2
P (t) = V (t) I (t) = = 2 · 2202 cos2 (2π50t)
R
P(t) (W)
8.00 · 10 4
6.00 · 10 4
4.00 · 10 4
2.00 · 10 4
t (s)
0.0100 0.0200 0.0300 0.0400 0.0500 0.0600 0.0700 0.0800
√
Por otra parte, la potencia media transmitida por la señal en un periodo T (llamando V0 = 220 2 y
ω0 = 50) es:
ˆ T ˆ T
1 1 V02 cos2 ωt V2
P (t) = V (t) I (t) dt = dt = 0
T 0 T 0 R 2R
Como la integral es al espacio, mientras que la derivada parcial es con respecto al tiempo, ambos operadores
conmutan. ¨ ¨ ¨ !
∂ ~
B ∂ ~
B
~ ×E
∇ ~ · dS
~= − ~⇔
· dS ~ ×E
∇ ~+ ~=0
· dS
S S ∂t S ∂t
Como lo anterior debe ser cierto para cualquier superficie S que se apoye sobre la curva C, el integrando debe
ser nulo:
~ ~
∇~ ×E ~ + ∂ B = ~0 ⇔ ∇ ~ ×E~ = − ∂B
∂t ∂t
Q.E.D.
Corolario 25. El campo eléctrico E ~ siempre puede expresarse como la suma del opuesto del gradiente de un
~
potencial φ y el opuesto de la derivada con respecto al tiempo del potencial magnético vector A.
~
~ R,
E ~ t = −∇φ
~ ~ t − ∂ A R,
R, ~ t
∂t
Demostración. Partimos del teorema 20 en la página anterior:
~
~ = − ∂B
~ ×E
∇
∂t
~ como el
Por la proposición 50 en la página 151, sabemos que podemos expresar la inducción magnética B
rotacional del potencial magnético vector. Así:
∇ ~ =−∂ ∇
~ ×E ~ ×A
~
∂t
Como el rotacional es un operador espacial y la derivada con respecto al tiempo es un operador temporal,
conmutan y obtenemos:
~ ~
∇~ ×E ~ = −∇ ~ × ∂A ⇔ ∇ ~ ×E ~ +∇ ~ × ∂ A = ~0
∂t ∂t
Como el rotacional es un operador lineal, obtenemos:
!
∂ ~
A
~ × E
∇ ~+ = ~0
∂t
~ ~
~ + ∂ A = −∇φ
E ~ ⇔E ~ − ∂A
~ = −∇φ
∂t ∂t
Q.E.D.
[Φm ] Wb
donde [Mij ] = = = H = Henrio.
[I] A
Además, llamamos coeficiente de autoinducción Mii al opuesto del factor que relaciona la fuerza
electromotriz presente en el conductor i-ésimo con la variación temporal de la corriente eléctrica que circula
por él.
∂Ii
εii = −Mii
∂t
Proposición 67. El coeficiente de inducción mutua Mij entre dos conductores lineales cerrados Ci y Cj puede
calcularse como: ¨ ˛
µ0 d~lj × ~r ~
Mij = · dSi
4π Si Cj r3
siendo Si la superficie definida por la curva Ci .
(j)
(i)
Ii (t) Ij (t)
Demostración. Por la proposición 46 en la página 144, sabemos que la inducción magnética generada por un
~ es:
conductor lineal cerrado Cj en un punto R
µ I ˛ d~l × ~r
~
B R~ = 0 j j
4π Cj r3
Por la definición 64 en la página 185, sabemos que el flujo magnético generado por el conductor j-ésimo a
través de la superficie Si viene dado por:
¨ ¨ ˛ ¨ ˛
~ · dS
~i = µ0 I j d~lj × ~r ~ µ0 Ij d~lj × ~r ~
Φij = B · dSi = · dSi
Si Si 4π Cj r3 4π Si Cj r3
Por el axioma 4 en la página 192, tenemos que la fuerza electromotriz presente en el conductor i-ésimo debido
al conductor j-ésimo puede expresarse como:
¨ ˛
∂Φij µ0 d~lj × ~r ~ ∂Ij
εij = − =− · dSi
∂t 4π Si Cj r3 ∂t
Q.E.D.
Proposición 68. Sea un sistema de n conductores lineales cerrados por los que circulan corrientes eléctricas
I1 , . . . , In . Sean, Φ1 , . . . , Φn los flujos magnéticos que circulan a través de las superficies S1 , . . . , Sn delimitadas
por las curvas C1 , . . . , Cn . Entonces los flujos magnéticos están relacionados con las intensidades de la siguiente
forma:
Φ=M·I
donde:
Φ1 I1 M11 ··· M1n
.. .. ..
Φ = ... I=. M= . .
Φn In Mn1 · · · Mnn
siendo Mij el coeficiente de inducción mutua entre los circuitos i-ésimo y j-ésimo.
(j)
(1)
I1
Ij
(n)
(i)
Ii In
Φij = Mij Ij + C
donde C es una constante. Dado que si Ij = 0, debe ser Φij = 0, deducimos que C = 0. Así, tenemos:
Φij = Mij Ij
Para obtener el flujo magnético total que fluye a través de la superficie Si únicamente tenemos que sumar a
todos los conductores j. Así:
Xn Xn
Φi = Φij = Mij Ij
j=1 j=1
Φ=M·I
Q.E.D.
Demostración. Por la definición 64 en la página 185, sabemos que el flujo magnético generado por el circuito
j-ésimo en el conductor i-ésimo viene dado por:
¨
Φij = B~ j · dS
~i
Si
siendo Si la superficie delimitada por la curva Ci . Por la proposición 54 en la página 160, podemos expresar
~j = ∇
B ~ ×A ~ j , obteniendo:
¨
Φij = ~ ×A
∇ ~ j · dS
~i
Si
Por último, por la proposición 53 en la página 157, sabemos que el potencial magnético vector generado por
el conductor j-ésimo viene dado por la expresión:
˛
~ j = µ0 Ij d~lj
A
4π Cj r
Q.E.D.
Demostración. Se sigue trivialmente del enunciado de la proposición 69 en la página anterior, pues el producto
escalar es conmutativo. Q.E.D.
Ejemplo 38. Queremos obtener el coeficiente de inducción mutua entre un circuito lineal indefinido y una
espira rectangular.
c
(1) (2)
Como, por el corolario 26 en la página anterior es M12 = M21 , podemos calcular este último que es más
sencillo. Por el teorema de Ampère en forma integral (ver teorema 18 en la página 155), obtenemos fácilmente
que la inducción magnética generada por un conductor rectilíneo indefinido es:
~ = µ0 I ϕ̂
B
2πr
Por la definición 64 en la página 185, tenemos:
¨ ˆ z0 +a ˆ c+b ˆ c+b
~ · dS
~= µ0 I µ0 I dr µ0 Ia c+b
Φ21 = B dzdr = a = ln
S z=z0 r=c 2πr 2π c r 2π c
Así:
µ0 a b+c
M12 = M21 = ln
2π c
F~O = q E
~ O + q~v × B
~O
Sin embargo, para el observador O0 , la carga está quieta; luego para él sobre la carga q únicamente actúa
~ O0 . Es decir:
el campo eléctrico E
F~O0 = q E
~ O0
Dicho de otra forma, para O0 se genera un campo eléctrico inducido por el movimiento.
Capítulo 14
Energía magnética
siendo Mij el coeficiente de inducción mutua entre los conductores i-ésimo y j-ésimo.
Demostración. Podemos suponer sin pérdida de generalidad que la corriente eléctrica Ii que circula por el
conductor i-ésimo es debida a una batería. De esta forma, la fuerza electromotriz presente en el conductor
i-ésimo tiene dos términos, uno de ellos viene dado por la batería, mientras que el otro se corresponde con la
fuerza electromotriz inducida. Por la ley de Ohm (ver teorema 14 en la página 133) y por la ley de inducción
de Faraday-Lenz (ver axioma 4 en la página 192), tendríamos:
∂Φi
εi − = Ii Ri
∂t
Multiplicando a ambos lados por Ii (t), tenemos:
∂Φi ∂Φi
εi Ii − Ii = Ii2 Ri ⇔ εi Ii = Ii + Ii2 Ri
∂t ∂t
El término de la izquierda se corresponde con la potencia suministrada por la batería, mientras que el segundo
sumando de la derecha se corresponde con la potencia disipada en el circuito por efecto Joule. De esta forma,
la variación con respecto al tiempo de energía potencial magnética en el circuito debe ser necesariamente el
primer sumando del lado derecho de la ecuación:
∂Ui ∂Φi
=P = Ii
∂t ∂t
Multiplicando ambos lados por ∂t, obtenemos:
∂Ui = ∂Φi Ii
Vamos a ir trayendo del infinito cada uno de los n conductores lineales cerrados. De esta forma, el trabajo
que tenemos que hacer para aumentar el flujo que circula por el primer conductor de cero a su flujo final
Φ1,final cuando sólo existe el primer conductor es:
ˆ Φ1,final
w1 = I1 dΦ1
0
Ahora, hacemos aparecer el conductor 2 y, de nuevo, calculamos el trabajo que tenemos que hacer para
aumentar el flujo que pasa por él entre 0 y su flujo final Φ2,final . De nuevo, por la proposición 68 en la
página 198, tenemos:
Φ2 = M21 I1 + M22 I2 ⇒ dΦ2 = M21 dI1 + M22 dI2
En consecuencia: ˆ ˆ
I1,final I2,final
1
w2 = I2 M21 dI1 + I2 M22 dI2 = M21 I1 I2 + M22 I22
0 0 2
Y, en general, para el conductor i-ésimo:
i−1
X i−1
X
Φi = Mij Ij + Mii Ii ⇒ dΦi = Mij dIj + Mii dIi
j=1 j=1
i−1 ˆ
X Ij,final ˆ Ii,final i−1
X 1
wi = Ii Mij dIj + Ii Mii dIi = Mij Ii Ij + Mii Ii2
0 0 2
j=1 j=1
Como, por el corolario 26 en la página 199, la matriz de coeficientes de inducción mutua es simétrica, podemos
escribir lo anterior como:
n
X X i−1
n X n
X Xn X n
1 1 1 1 1
Um = Mii Ii2 + Mij Ii Ij + Mji Ij Ii = Mii Ii2 + Mij Ii Ij =
2 2 2 2 2
i=1 i=1 j=1 i=1 i=1
j=1
j 6= i
n n
1 XX
= Mij Ii Ij
2
i=1 j=1
Q.E.D.
Proposición 71. La energía magnética almacenada Um en una distribución (con volumen V ) de conductores
viene dada por:
˚ ˆ A~
Um = J~C · δ AdV
~
V ~ ~0
A=
~ es la excitación magnética y B
donde T.E. indica todo el espacio, H ~ es la inducción magnética.
Demostración. Vamos a dividir nuestra distribución volumétrica en infinitos circuitos elementales por los que
circulan corrientes reales (de conducción).
dS
~ R)
J( ~
siendo U la energía potencial magnética. Diferenciando a ambos lados la expresión anterior, obtenemos:
d (dU ) = dI
¯ C dΦm (14.2.0.1)
Por la proposición 50 en la página 151, podemos expresar lo anterior en función del potencial magnético vector:
¨
Φm = ~ ×A
∇ ~ · dS~C
SC
donde C es la curva delimitada por la superficie SC . Tomando la diferencial a ambos lados, obtenemos:
~ · d~l
dΦm = A
Sustituyendo en la ecuación Por otra parte, por la proposición 43 en la página 126, tenemos:
¯ C = J~C · dS
dI ~
Nótese que S 6= SC .
Sustituyendo en la ecuación 14.2.0.1 en la página anterior, obtenemos:
d (dU ) = J~C · dS
~A~ · d~l
S1
S2
d~l ~
dS
~ k d~l (esto es trivialmente cierto en un circuito lineal), podemos reescribir lo anterior como:
Como J~C k dS
d (dU ) = J~C · A
~ dSdl = J~C · A
~ dV
Bien, para obtener la energía magnética almacenada, vamos a ver cuánto varía la energía potencial mag-
~ De esta forma, aplicamos el operador de diferencia virtual δ a ambos
nética al variar el potencial vector A.
lados, obteniendo:
δ [d (dU )] = J~C · δ A
~ dV
pues no estamos variando la densidad volumétrica de corriente real. Ahora, el término izquierdo de nuestra
ecuación se corresponde con la variación del diferencial inexacto de la energía magnética almacenada. Así:
¯ m = J~C · δ A
dδ̄U ~ dV
~ en su versión diferencial
En régimen estacionario de corrientes, por el teorema de Ampère para el campo H
(ver teorema 19 en la página 174), podemos reescribimos lo anterior como:
˚ ˆ ~
A
Um = ~ ×H
∇ ~ · δ AdV
~ (14.2.0.2)
V ~ ~0
A=
Nótese que, por la proposición 50 en la página 151 y dado que los operadores virtuales y espaciales conmutan,
se tiene:
δB~ =δ ∇ ~ ×A~ =∇ ~ × δA ~
Sustituyendo esto y por el teorema de la divergencia (ver teorema 2 en la página 17), llegamos a:
‹ ˆ ~
A ˚ ˆ ~
B
Um = ~ ~ ~
H × δ A · dS + ~ · δ BdV
H ~ (14.2.0.3)
S ~ ~0
A= V ~ ~0
B=
Si extendemos la integral a todo el espacio (no sólo a nuestro volumen V ), la superficie de nuestra integral
de la izquierda es la superficie del infinito y allí:
˚ ~
~ = 1 J × ~r 0 ~ ∼ 1
H 3
dV ⇒ H
4π V0 r r2
˚
~= 1 J~ 0 ~ = 1 1 ⇒ δA
~ ∼1
A dV ⇒ δ A
4π V0 r 4π r r
S ∼ r2
Por tanto:
‹ ˆ ~
A ˆ ˆ
~ ~ ~ 1 1 1 1
H × δ A · dS ∼ lı́m 2
rdr = lı́m 2
dr = − −−−→ 0
Σ∞ ~ ~0
A= r→∞ r r r→∞ r r r→∞
En consecuencia, al extender las integrales de la ecuación 14.2.0.3 a todo el espacio, obtenemos:
˚ ˆ ~
B
Um = ~ · δ BdV
H ~
T.E. ~ ~0
B=
Q.E.D.
Corolario 27. En régimen estacionario de corrientes, la energía magnética almancenada Um en una distri-
bución l.h.i. (con volumen V ) de conductores viene dada por:
˚ ˚ ˚
1 1 ~ ·H~ dV = 1
Um = 2
µH dV = B J~C · A
~ dV
2 T.E. 2 T.E. 2 V
donde T.E. indica todo el espacio, J~ es la densidad volumétrica de corriente eléctrica, A ~ es el potencial
~ es la excitación magnética, B
magnético vector, H ~ es la inducción magnética y µ es la permeabilidad magnética
del medio.
S ∼ r2
Por tanto: ‹ ˆ ˆ
~ ~ ~ 1 1 1 1
A × H · dS ∼ lı́m 2
rdr = lı́m 2
dr = − −−−→ 0
Σ∞ r→∞ rr r→∞ r r r→∞
Así, la expresión 14.2.0.5 se simplifica a:
˚
1
Um = J~C · A
~ dV
2 T.E.
Finalmente, como J~C = ~0 fuera del volumen del conductor, la integral anterior se reduce a:
˚
1
Um = J~C · A
~ dV
2 V
Q.E.D.
Ejemplo 39 (Problema 67). Tenemos un hilo que atraviesa una corriente lineal cerrada de radio a según
descrito en el dibujo:
(0, 0, −z)
d~l
α
z
~r
θ
=sen α
z }| {
µ0 I dy z µ0 I dy r sen θ µ0 I dy z µ0 Iz dy
dBx = = = =
4π r 3 4π r 3 4π r r2 4π (y 2 + z 2 ) 23
ˆ ∞ " #∞
µ0 I dy µ0 I y µ0 I
Bx = z 3 = p =
4π y=0 (y 2 + z 2 ) 2 4π z 2 y 2 + z 2 4πz 2
y=0
Ejemplo 40 (Problema 68). Tenemos una superficie plana situada en plano XY por el que circula una
corriente superficial:
J~ (x, y, 0) = J0 x̂
Nos piden determinar el potencial magnético vector. Para ello, podríamos hallar primer la inducción magnética.
¨
~ = µ0 J~S × ~r 0
B dS
4π S0 r3
µ0 J0 0
~ =
dB dS
4π r2
Sin embargo, en vez de hacerlo como hemos puesto antes, podemos aplicar la ley de Ampère para el campo
~
H en su forma integral (ver teorema 19 en la página 174) en un rectángulo:
(2)
(1)
(3)
(4)
˛
~ · d~l = I
H
C
ˆ ˆ ˆ ˆ
~ · d~l +
H ~ · d~l +
H ~ · d~l +
H ~ · d~l = I ⇔
H
| {z } | {z }
(1) (2) (3) (4)
~
H⊥d~l⇒=0 ~
H⊥d~l⇒=0
a a a
⇔H z= b+H z =− b = I ⇔ 2H z = b = I = JS b ⇔
2 | {z 2 } 2
=H (z= a2 )
a JS
H z= =
2 2
Vemos, que efectivamente se cumple la proposición 65 en la página 183:
H2t − H1t = JS
Por otra parte, como estamos en el vacío, por la definición 61 en la página 179, tenemos:
~ = µ0 H
B ~
En consecuencia:
~ z = µ0 J0 (−ŷ) si z > 0
B
2
~ z = µ0 J0 ŷ si z < 0
B
2
Por la proposición 50 en la página 151, sabemos que es:
~ =∇
B ~ ×A
~
~ = (0, B, 0) = ∇
B ~ ×A
~
y
Por otra parte, por la proposición 50 en la página 151, sabemos que podemos calcular el potencial magnético
vector de la siguiente forma:
¨ ~
~ = µ0
A
JS
dS
4π S r
Así, vemos que, necesariamente será Ay = 0 = Az . Por tanto:
~ ~ ∂Ax
∇ × A = 0, ,0
∂z
De esta forma, para z > 0:
∂Ax µ0 J 0 µ0 J0
=− ⇔ Ax = − z + k1
∂z 2 2
y para z < 0:
∂Ax µ0 J 0 µ0 J 0
= ⇔ Ax = z + k2
∂z 2 2
Como el potencial es una función continua, debe ser k1 = k2 . Podemos tomar el origen de potencial en el
~ = ~0.
plano de forma que allí sería A
Demostración. Por el axioma 3 en la página 146, tenemos que la fuerza que actúa sobre un diferencial de
corriente lineal por la que circula una corriente dI
¯ viene dada por:
¯ ~l × B
dF~ = dId ~
d~l
~ k d~l × J,
Como en una corriente lineal es dS ~ podemos reescribir lo anterior como:
dF~ = J~ × B
~ dV
Definición 68. Llamamos densidad de fuerza magnética f~ a la fuerza por unidad de volumen que sufre
una distribución de corrientes eléctricas.
~
~ := dF R
f~ R ~
dV
Por la proposición 72 en la página anterior, tenemos:
f~ = J~ × B
~
suministrada por la batería está fija, la única forma en la que se puede influir en el sistema es a través de la
variación del flujo magnético. De esta forma, la potencia P que un agente hace sobre el circuito viene dada
por:
dW
¯ ∂Φi
=P =
dt ∂t
Por ende, si el trabajo que hacemos sobre el sistema es nulo, la variación de flujo también debe ser nula.
Q.E.D.
F~ = −∇U
~ m
¯ = −dUm
dW
Por otra parte, por definición de trabajo:
¯ = F~ · d~r
dW
Por tanto:
¯ = −dUm = F~ · d~r ⇔
dW
dUm
⇔ F~ = − ~ m
≡ −∇U
d~r
Q.E.D.
F~ = +∇U
~ m
Demostración. Podemos suponer sin pérdida de generalidad que la corriente eléctrica Ii que circula por el
conductor i-ésimo es debida a una batería. De esta forma, la fuerza electromotriz presente en el conductor
i-ésimo tiene dos términos, uno de ellos viene dado por la batería, mientras que el otro se corresponde con la
fuerza electromotriz inducida. Por la ley de Ohm (ver teorema 14 en la página 133) y por la ley de inducción
de Faraday-Lenz (ver axioma 4 en la página 192), para cada uno de los circuitos elementales tenemos:
∂Φi
εi − = Ii Ri
∂t
Multiplicando a ambos lados por Ii (t), tenemos:
∂Φi ∂Φi
εi Ii − Ii = Ii2 Ri ⇔ εi Ii = Ii + Ii2 Ri
∂t ∂t
El término de la izquierda se corresponde con la potencia suministrada por la batería, mientras que el segundo
sumando de la derecha se corresponde con la potencia disipada en el circuito por efecto Joule. Multiplicando
la ecuación anterior por ∂t, obtenemos:
εi Ii ∂t = ∂Φi Ii + Ii2 Ri ∂t
Llamaremos Ub a la energía potencial que mantiene constantes las intensidades I1 , . . . , In . De esta forma, debe
ser:
dUb,i = dΦi Ii
Sumando a todos los conductores, obtenemos:
n
X
dUb = Ii dΦi (14.4.2.1)
i=1
dUb = dW
¯ + dUm
¯ = F~ · d~r, obtenemos:
Usando 14.4.2.2 y aplicando la definición de trabajo dW
dUm
2dUm = F~ · d~r + dUm ⇔ dUm = F~ · d~r ⇔ F~ = ~ m
= ∇U
d~r
Q.E.D.
Capítulo 15
Ecuaciones de Maxwell
∇ ~ = ρtotal
~ ·E ~ ·B
∇ ~ =0
ε0
~
∇ ~ = − ∂B
~ ×E ~ ×H
∇ ~ = JC
∂t
~ es el campo eléctrico, B
donde E ~ es la inducción magnética, H ~ es la excitación magnética, ρtotal es la densidad
volumétrica de carga total y ε0 es la permitividad eléctrica del vacío.
Demostración. Lo anterior es simplemente una recopilación de los teoremas 1 en la página 14, 19 en la
página 174 y 20 en la página 195 y del corolario 21 en la página 152. Q.E.D.
~ del espacio se cumple que:
Teorema 22 (Ley de Ampère generalizada). En todo punto R
~
~ ×B
∇ ~ R ~ + µ0 ∂ D R
~ = µ0 J~ R ~
∂t
donde B ~ es la inducción magnética, J~ es la densidad de corriente volumétrica y D
~ es el vector desplazamiento
eléctrico.
Demostración. Partimos de la ecuación de continuidad (ver teorema 12 en la página 128):
~ · J~ + ∂ρf = 0
∇
∂t
Por el teorema de Gauss generalizado en forma diferencial (ver teorema 9 en la página 88), tenemos:
~ ·D
ρf = ∇ ~
Sustituyendo, tenemos:
~ · J~ + ∂ ∇
∇ ~ ·D
~ =0
∂t
Como la divergencia es un operador espacial mientras que la parcial con respecto al tiempo es un operador
temporal, ambos conmutan:
! !
∂ ~
D ∂ ~
D
∇~ · J~ + ∇
~ · =0⇔∇ ~ · J~ = −∇
~ · (15.1.0.1)
∂t ∂t
Por el lema 6 en la página 152, la divergencia del término izquierdo de la ecuación es nula. Sin embargo, por
el teorema 12 en la página 128, sabemos que la divergencia del término derecho no es nulo a menos que nos
encontremos en régimen estacionario de corrientes. Por tanto, tiene sentido pensar que hace falta añadir un
término en el lado derecho de la ecuación que permita que se cumpla la ecuación de continuidad. Llamemos
a dicho término X~ : R3 −→ R3 . Así, tendríamos:
∇ ~ = µ0 J~ + X
~ ×B ~
∂D~
~ ~ ~ ~ ~
∇ × H R = JC R + ~
R
∂t
donde H ~ es la excitación magnética, J~C es la densidad de corriente volumétrica real (de conducción) y D
~ es
el vector desplazamiento eléctrico.
Demostración. Partimos del teorema 22 en la página anterior:
~
~ = µ0 J~ + µ0 ∂ D
~ ×B
∇
∂t
Por la observación 43 en la página 174, tenemos:
h i ~
∂D ~
~ × µ0 H
∇ ~ +M
~ = µ0 J~ + µ0 ⇔ µ0 ∇~ ×H
~ +∇ ~ = µ0 J~ + µ0 ∂ D ⇔
~ ×M
∂t ∂t
~
∂D
~ ×H
⇔∇ ~ = J~ − ∇
~ ×M
~ +
∂t
Por la proposición 60 en la página 166, tenemos:
~ ~
~ = J~ − J~m + ∂ D = J~C + ∂ D
~ ×H
∇
∂t ∂t
Teorema 24 (Ecuaciones de Maxwell). El electromagnetismo clásico queda resumido con las ecuaciones para
~ y la inducción magnética B:
el campo eléctrico E ~
∇ ~ = ρtotal
~ ·E ~ ·B
∇ ~ = ~0
ε0
!
~ ~
∇ ~ = − ∂B
~ ×E ~ ×B
∇ ~ = µ0 J~C + J~m +
∂D
∂t ∂t
Demostración. Lo anterior no es más que una recopilación de los teoremas 1 en la página 14, 22 en la página 214
y 20 en la página 195 y del corolario 21 en la página 152. Q.E.D.
Teorema 25 (Resumen del electromagnetismo). Todo el electromagnetismo clásico queda descrito con las
ecuaciones:
~ ·D
∇ ~ = ρf ~ ·B
∇ ~ = ~0
~ ~
∇ ~ = − ∂B
~ ×E ~ = J~C + ∂ D
~ ×H
∇
∂t ∂t
~ = µ0 H
B ~ +M
~ ~ = ε0 E
D ~ + P~
~ · J~ + ∂ρf = 0
∇
∂t
donde E ~ es el campo eléctrico, D~ es el vector desplazamiento eléctrico, P~ es el vector polarización, ρf es la
densidad volumétrica de carga libre, B ~ es la inducción magnética, H ~ es la excitación magnética, J~C es la
~
densidad volumétrica de corriente real, J es la densidad voumétrica de corriente total, ε0 es la permitividad
eléctrica del vacío y µ0 es la permeabilidad magnética del vacío.
Demostración. Lo anterior no es más que una recopilación de los teoremas 9 en la página 88, 23 en la página
anterior, 20 en la página 195, 12 en la página 128, del corolario 21 en la página 152, de la definición 25 en la
página 88 y de la observación 43 en la página 174. Q.E.D.
~ R,
B ~ t =B
~ ~r R
~ cos (ωt + ϕB )
donde E~r R~ y B~r R ~ son funciones cualquiera que incluyen la dependencia del campo eléctrico y de la
inducción magnética. Veamos que esto es consistente con las ecuaciones de Maxwell:
~ ·E
∇ ~ (R, t) = ∇~ ·E ~ cos (ωt + ϕE )
~ ~r R
pues el término del coseno es constante para la divergencia. Por el teorema 24 en la página anterior, tenemos:
ρ R,~ t
~ ·E
∇ ~ ~r R
~ cos (ωt + ϕE ) = ∇~ ·E
~ (R, t) = ⇔ ρ R,~ t = ε0 ∇~ ·E
~ ~r R
~ cos (ωt + ϕE )
ε0
De esta forma, la densidad volumétrica de carga ρ tiene forma armónica. Prosigamos:
~ ·B
∇ ~ (R, t) = ∇
~ ·B
~ ~r R
~ cos (ωt + ϕB )
Prosigamos:
~ ×E
∇ ~ R,~ t =∇×E
~ ~r R
~ cos (ωt + ϕE )
~
∂D
Luego, necesariamente tanto J~ como deben tener forma armónica para que se cumpla la ley de Ampère
∂t
generalizada.
Es decir, hemos probado que existen soluciones de la ecuación de Maxwell tales que todas las funciones
que aparecen tiene forma armónica. Q.E.D.
Observación 47. La ecuación de ondas anterior también es válida para funciones escalares, tomando m = 1.
1.
~ = ρf
~ ·E
∇
ε
2.
~ ·B
∇ ~ =0
3.
~
~ = − ∂B
~ ×E
∇
∂t
4.
~
~ = µJ~C + εµ ∂ E
~ ×B
∇
∂t
y se hizo la luz.
La velocidad de propagación de la luz v en un medio de permitividad eléctrica ε y permeabilidad magnética
µ viene dada por:
1
v=√
εµ
Además, las ondas electromagnéticas vienen descritas por:
∂2E~ ~
~ − εµ
∇2 E = ~ ρ + µσ ∂ E
∇
∂t2 ε ∂t
∂2B~ ~
∂B
~ − εµ
∇2 B = µσ
∂t2 ∂t
Demostración. A lo largo de la demostración usaremos ρ ≡ ρf y J~ ≡ J~C pero el lector debe recordar que en
realidad estamos hablando de la densidad de carga libre y de la densidad volumétrica de corriente eléctrica.
Partimos de la tercera ecuación de Maxwell (la que contiene ∇ ~ × E)
~ y hacemos el rotacional a ambos
lados. !
∂ ~
B
~ × ∇
∇ ~ ×E~ =∇ ~ × −
∂t
Por el lema 7 en la página 152 y como el rotacional es un operador espacial y la derivada con respecto al
tiempo es un operador temporal, ambos conmutan, obtenemos:
~ ∇
∇ ~ ·E ~ =−∂ ∇
~ − ∇2 E ~ ×B
~
∂t
Por el lema 7 en la página 152 y como el rotacional es un operador espacial y la derivada con respecto al
tiempo es un operador temporal, ambos conmutan, obtenemos:
~ ∇
∇ ~ ·B
~ − ∇2 B~ = µ∇~ × J~ + εµ ∂ ∇ ~ ×E~
∂t
Por la segunda y la tercera ecuación de Maxwell, lo anterior es equivalente a:
2~
~ × J~ − εµ ∂ B ⇔
~ = µ∇
−∇2 B
∂t2
~
∂2B
~ − εµ
⇔ ∇2 B ~ × J~
= −µ∇
∂t2
Por la definición 41 en la página 132, podemos expresar la densidad volumétrica de corriente eléctrica J~ en
~ obteniendo:
función del campo eléctrico E,
~
∂2B
~ − εµ
∇2 B ~ ×E
= −µσ ∇ ~
∂t2
Y, de nuevo, por la tercera ecuación de Maxwell, llegamos a:
∂2B~ ~
∂B
~ − εµ
∇2 B = µσ (15.3.0.2)
∂t2 ∂t
Como podemos ver, las ecuaciones 15.3.0.1 y 15.3.0.2 satisfacen la ecuación de ondas dada en la definición 69
en la página anterior. En consecuencia, E~ y B~ se propagan en el espacio como ondas. Por analogía con la
ecuación de ondas (ver definición 69 en la página anterior), tenemos la velocidad de propagación en ambas
ondas es la misma y viene dada por:
1 1
v2 = ⇔v= √
εµ εµ
De hecho, en el vacío, obtenemos:
1
v=√ =c
ε0 µ 0
donde c es la velocidad de la luz en el vacío. Q.E.D.
Esto último es una coincidencia casual, es decir, nada de la teoría que hemos visto hasta ahora justifica que
justo el valor del cociente anterior sea la velocidad de la luz. Sin embargo; este hecho, unido a que las ecuaciones
de Maxwell admiten una solución ondulatoria, apunta a que la luz es, precisamente, una onda electromagnética.
En otras palabras, hemos visto que las ecuaciones de Maxwell son compatibles con la existencia de la luz.
Observación 48. Advertimos al lector de que si se encuentra en cuarto mal iluminado y enuncia las ecuaciones
de Maxwell, no se creará luz que le permita ver mejor.
Observación 49. De hecho, el índice de refracción de un medio con permitividad eléctrica ε y permeabilidad
magnética µ viene dado por: r
c εµ √
n= = = εr µ r
v ε0 µ 0
Corolario 28 (Ecuaciones de la telegrafía). En un medio sin cargas libres ρf = 0, las ondas electromagnéticas
vienen descritas por las ecuaciones:
2~ ~
∇2 E~ − εµ ∂ E − µσ ∂ E = 0
∂t 2 ∂t
2 ~ ~
∇2 B~ − εµ ∂ B − µσ ∂ B = 0
∂t 2 ∂t
Demostración. La demostración se sigue trivialmente del teorema 26 en la página 218 al sustituir ρf =
0. Q.E.D.
~
Proposición 77. Recordemos que, por la proposición 52 en la página 156, el potencial magnético vector A
no es único. Si imponemos la condición (llamada de Lorenz):
~ + εµ ∂φ = 0
~ ·A
∇
∂t
siendo φ el potencial escalar que aparece en el corolario 25 en la página 196; entonces, tanto el potencial
~ como el potencial escalar φ satisfacen la ecuación de ondas dada por la definición 69 en
magnético vector A
la página 218.
Demostración. Por el corolario 25 en la página 196, podemos expresar el campo eléctrico como sigue:
~
~ − ∂A
~ = −∇φ
E
∂t
Haciendo la divergencia a ambos lados, obtenemos:
!
~
∂A
~ ·E
∇ ~ = −∇
~ · ∇φ
~ ~ ·
−∇
∂t
Por la primera ecuación del 24 en la página 216 y, dado que la divergencia es un operador espacial mientras que
la parcial con respecto al tiempo es un operador temporal, ambos operadores conmutan, por lo que obtenemos:
ρ (R, t) ∂ ~ ~
= −∇2 φ − ∇·A
ε ∂t
Por hipótesis, tenemos:
ρ (R, t) ∂ ∂φ
= −∇2 φ − −εµ ⇔
ε ∂t ∂t
ρ (R, t) ∂2φ
⇔ = −∇2 φ + εµ 2 ⇔
ε ∂t
∂2φ ρ (R, t)
⇔ ∇2 φ − εµ 2
=−
∂t ε
Y esta última ecuación satisface la ecuación de ondas dada en la definición 69 en la página 218.
Por otra parte, por la cuarta ecuación del teorema 24 en la página 216:
~
~ = µJ~ + εµ ∂ E
~ ×B
∇
∂t
Por la proposición 50 en la página 151, podemos escribir lo anterior como:
~
~ ×∇
∇ ~ = µJ~ + εµ ∂ E
~ ×A
∂t
Por el lema 7 en la página 152 y el corolario 25 en la página 196, tenemos:
" #
∂ ∂ ~
A
~ ∇
∇ ~ ·A~ − ∇2 A ~ = µJ~ + εµ ~ −
−∇φ
∂t ∂t
Como el gradiente es un operador espacial, mientras que la parcial con respecto al tiempo es un operador
temporal, ambos operadores conmutan. En consecuencia, llegamos a:
~
~ ~ ~ 2~ ~ ~ ∂φ ∂2A
∇ ∇ · A − ∇ A = µJ − ∇ εµ − εµ 2 ⇔
∂t ∂t
~
2~ ∂2A ~ ∂φ ~ ∇ ~ ·A
~ = −µJ~
⇔ ∇ A − εµ 2 − ∇ εµ −∇
∂t ∂t
Como el gradiente es un operador lineal, obtenemos:
~
2~ ∂2A ~ ∂φ ~ ~
∇ A − εµ 2 − ∇ εµ + ∇ · A = −µJ~
∂t ∂t
Pero, por hipótesis, lo que hay dentro del gradiente se anula. Así, obtenemos:
~
∂2A
~ − εµ
∇2 A = −µJ~
∂t2
Y esta última ecuación satisface la ecuación de ondas dada en la definición 69 en la página 218. Q.E.D.
Como el gradiente y la parcial con respecto al tiempo son operadores lineales y el orden del rotacional conmuta
con la derivada parcial:
∂A~0 ~ ~
~ 0−
−∇φ ~ − ∂ ∇ψ
= −∇φ ~ − ∂ A + ∂ ∇ψ ~ − ∂A
~ = −∇φ
∂t ∂t ∂t ∂t ∂t
Y, por el corolario 7.3.0.3 en la página 115, obtenemos:
∂A~0
~ 0−
−∇φ ~
=E
∂t
Por otra parte, calculemos:
~ ×A
∇ ~0 = ∇
~ × A~ − ∇ψ
~
Y, por el lema 8 en la página 156, el segundo sumando es nulo. Así, por la proposición 50 en la página 151,
obtenemos:
~ ×A
∇ ~0 = ∇
~ ×A~=B ~
∂A~0
~ 0−
En consecuencia, como −∇φ yB ~ =∇ ~ ×A ~ 0 , por la proposición 77 en la página 220, como A
~ 0 y φ0
∂t
satisfacen la condición de Lorenz, entonces A ~ 0 y φ0 satisfacen la ecuación de ondas.
Por último, por hipótesis, tenemos:
∂ψ
0 ∂ φ +
0=∇ ~ 0 + εµ ∂φ = ∇
~ ·A ~ · A ~ − ∇ψ
~ + εµ
∂t
∂t ∂t
Como la parcial y la divergencia son operadores lineales, obtenemos:
2
⇔∇~ ·A~ + εµ ∂φ − ∇2 ψ + εµ ∂ ψ = 0 ⇔
∂t ∂t2
~ ~ ∂φ 2 ∂2ψ
⇔ ∇ · A + εµ − ∇ ψ − εµ 2 = 0
∂t ∂t
Así, si es:
∂2ψ
∇2 ψ − εµ =0
∂t2
~ y φ satisfacen la condición de
ψ satisface la ecuación de ondas dada en la definición 69 en la página 218 y A
~
Lorenz. Por la proposición 77 en la página 220, A y φ satisfacen la ecuación de ondas. Q.E.D.
Capítulo 16
Energía electromagnética
~ R,
B ~ 0 cos ~k · R
~ t =B ~ − ωt
~
∂2E
~ − εµ
∇2 E =0
∂t2
~
∂2B
~ − εµ
∇2 B =0
∂t2
~ la habremos resuelto para B.
Como vemos, por simetría de las ecuaciones, si resolvemos la ecuación para E, ~
Vamos a trabajar en el campo de los complejos. Como E ~ tiene dependencia armónica con respecto al tiempo,
tenemos: h i
~ R,
E ~ t = Re E ~ ~r R
~ e−iωt (16.1.0.1)
~ C − εµ ∂2 ~
∇2 E EC = ~0
∂t2
Como, en el campo de los complejos derivar con respecto al tiempo es equivalente a multiplicar por −iω,
obtenemos:
~ ~r e−iωt + εµω 2 E
∇2 E ~ ~r e−iωt = ~0 ⇔
⇔ ∇2 E ~ ~r + εµω 2 E
~ ~r e−iωt = ~0
De nuevo, como las ecuaciones son simétricas, se resolvemos una, habremos resuelto el resto. Vamos con la
ecuación con respecto a la coordenada x. Tenemos:
∂ 2 Ex ∂ 2 Ex ∂ 2 Ex
+ + + k 2 Ex = 0 (16.1.0.2)
∂x2 ∂y 2 ∂z 2
Hagamos separación de variables, es decir, supondremos que la solución es de la forma:
d2 X d2 Y d2 Z
Y Z + XZ + XY + k 2 XY Z = 0
dx2 dy 2 dz 2
Como no nos interesa la solución trivial Ex = 0, podemos dividir por XY Z, con lo que llegamos a:
d2 X 1 d2 Y 1 d2 Z 1 2 d2 X 1 2 d2 X
+ + + k = 0 ⇔ + C x = 0 ⇔ + Cx2 X = 0
dx2 X dy 2 Y dz 2 Z dx2 X dx2
| {z }
=:Cx2
Nótese que todo el término Cx2 es constante para x y que Cx2 ∈ R. Así, hemos llegado a una ecuación diferencial
de una variable. Siempre que sea Cx 6= 0, la solución de la ecuación anterior viene dada por:
d2 Y 1 d2 Z 1 2 2 d2 Y 1 d2 Z 1
+ + k = C x ⇔ + + k 2 − Cx2 = 0 ⇔
dy 2 Y dz 2 Z dy 2 Y |dz 2 Z
{z }
=:Cy2
d2 Y 1 2 d2 Y
⇔ + C y = 0 ⇔ + Cy2 Y = 0
dy 2 Y dy 2
De nuevo, nótese que todo el término Cy es constante para y y que Cy2 ∈ R. Así, hemos llegado a una ecuación
diferencial de una variable. Siempre que sea Cy 6= 0, la solución de la ecuación anterior viene dada por:
d2 Z 1 2 2 2 d2 Z 1
+ k − C x = C y ⇔ + k 2 − Cx2 − Cy2 = 0 ⇔
dz 2 Z dz 2 Z | {z }
=:Cz2
d2 Z 1 2 d2 Z
⇔ + C z = 0 ⇔ + Cz2 Z = 0
dz 2 Z dz 2
De nuevo, nótese que todo el término Cz es constante para z y que Cz2 ∈ R. Así, hemos llegado a una ecuación
diferencial de una variable. Siempre que sea Cz 6= 0, la solución de la ecuación anterior viene dada por:
Como Cx2 , Cy2 , Cz2 ∈ R, las soluciones obtenidas son las de un oscilador armónico. Recordemos que debe
cumplirse:
Cx2 + Cy2 + Cz2 = k 2
Como k 2 ∈ R, podemos imaginar que existe un vector ~k = (Cx , Cy , Cz ) que es solución de la ecuación anterior.
En consecuencia, sabemos que existen Dx , Dy , Dz , φx , φy , φz tales que podemos expresar la solución como:
h i
Ex (x, y, z) = Re Dx eikx x+iφx Dy eiky y+iφy Dz eikz z+iφz
Como, sólo nos interesan soluciones que sean isótropas en la dependencia espacio-temporal, por analogía,
tenemos que: h i
~ ~
Ey (x, y, z) = Ny Re eik·R
h i
~ ~
Ez (x, y, z) = Nz Re eik·R
~ 0 = (Nx , Ny , Nz ), obtenemos:
En consecuencia, llamando E
h i
~ ~r R
E ~ =E ~ 0 Re ei~k·R~
~ obtendríamos:
Actuando análogamente para B,
~ R,
B ~ 0 cos ~k · R
~ t =B ~ − ωt
Q.E.D.
Demostración. Siguiendo los mismos pasos que en la demostración de la proposición 79 en la página 223, pero
sin suponer σ = 0, obtenemos que el valor de k 2 viene dado por:
k 2 = ω 2 εµ + iωµσ
Definición 70. Denominamos factor de calidad Q al cociente entre la permitividad eléctrica ε de un medio
y su conductividad σ multiplicado por la frecuencia de la onda.
εω
Q :=
σ
Observación 50. Según la definición 70, un medio con un alto factor de calidad Q 1 se comportará como
un dieléctrico puro, mientras que un medio con un bajo factor de calidad Q 1, se comportará como un
conductor puro.
V A J W
[S] = = = 2
mm 2 s·m 2 m
Definición 73. Llamamos energía electromagnética almacenada de una onda a la suma de su energía
electroestática almacenada y su energía magnética almacenada.
Uem := Ue + Um
δ̄¯ ~ · δD
dUe = E ~
δ̄¯ ~ · δB
dUm = H ~
δ̄¯
dUe ~
~ · δD
=E
δt δt
δ̄¯
dUm ~
~ · δB
=H
δt δt
Sumando ambas ecuaciones, obtenemos:
δ̄¯
dUe δ̄¯
dUm ~ ~
+ ~ · δD + H
=E ~ · δB
δt δt δt δt
Como, la diferencial es un operador lineal, obtenemos:
δ̄¯
d (Ue + Um ) ~ ~
~ · δD + H
=E ~ · δB
δt δt δt
δ̄¯
dUem ~ ~
~ · δD + H
=E ~ · δB (16.4.0.1)
δt δt δt
Por otra parte, por la tercera ecuación de Maxwell (ver teorema 24 en la página 216), tenemos:
~
~ = − ∂B
~ ×E
∇
∂t
~ obtenemos:
Multiplicando escalarmente a ambos lados por H,
~
~ · ∇
H ~ ×E ~ · ∂B
~ = −H (16.4.0.2)
∂t
~ (ver teorema 25 en la página 216), tenemos:
Además, por la cuarta ecuación de Maxwell para H
~
~ = J~C + ∂ D
~ ×H
∇
∂t
~ a ambos lados, llegamos a:
Multiplicando escalarmente por E
~
~· ∇
E ~ ×H
~ = J~C · E ~ ∂D
~ +E (16.4.0.3)
∂t
Restándole la ecuación 16.4.0.3 a la ecuación 16.4.0.2, obtenemos:
~ ~
~ · ∇
H ~ ×E
~ −E~· ∇~ ×H
~ = −J~C · E ~ · ∂D − H
~ −E ~ · ∂B
∂t ∂t
Por el lema 11 en la página 204, podemos expresar el lado izquierdo de la ecuación anterior como:
~ ~
~ · E
∇ ~ ×H
~ = −J~C · E ~ · ∂D − H
~ −E ~ · ∂B
∂t ∂t
Usando la ecuación 16.4.0.1, podemos reescribir el término derecho de la ecuación como sigue:
~ · E
∇ ~ ×H ~ = −J~C · E ~ − δ̄¯
dUem
δt
Por la definición 72 en la página anterior, podemos expresar la ecuación anterior en función del vector de
Poynting.
~ ·S
∇ ~ − δ̄¯
~ = −J~C · E dUem
δt
Ahora, integramos al volumen V a ambos lados:
˚ ˚ ˚
~ ~ ~ ~ δ̄¯
dUem
∇ · SdV = − JC · EdV − dV
V V V δt
Por el teorema de la divergencia (ver teorema 2 en la página 17), lo anterior es equivalente a:
‹ ˚
~ · dS
~0 = − δ̄Uem
S J~C · EdV
~ −
S 0 V δt
siendo S 0 la superficie que delimita el volumen V . Pasando términos al otro lado, obtenemos:
‹ ˚
δ̄Uem ~ ~
⇔ =− 0
S · dS − J~C · EdV
~
δt S0 V
Q.E.D.
Apéndice A
Registro de cambios
por:
~ · H=
~ 1 ~ ~ ~ ~
∇ ∇·B−∇·M
µ0
~ es la excitación magnética y B
donde T.E. indica todo el espacio, H ~ es la inducción magnética.
εi Ii ∂t = ∂Φi Ii + Ii2 Ri
por:
εi Ii ∂t = ∂Φi Ii + Ii2 Ri ∂t
por: D E
J~S R,
~ t := σ R,
~ t ~v R,
~ t
~ es la excitación magnética y B
donde T.E. indica todo el espacio, H ~ es la inducción magnética.
por:
La energía magnética almacenada Um en una distribución (con volumen V ) de conductores viene dada
por:
˚ ˆ A~
Um = J~C · δ AdV
~
V ~ ~0
A=
~ es la excitación magnética y B
donde T.E. indica todo el espacio, H ~ es la inducción magnética.
En la demostración del teorema 17 en la página 149, justo encima de la ecuación 10.4.1.1, se ha cambiado
«Como la integral es a las coordenadas primadas, mientras que la divergencia y la laplacian son a las
coordenadas sin primar, ambos operadores conmutan.» por «Como la integral es a las coordenadas
primadas, mientras que la divergencia y la laplaciana son a las coordenadas sin primar, ambos operadores
conmutan.»
por:
NI
N I = m = Φ (R1 + Rx + R2 + Rx + R1 ) ⇔ Φ = l1 2x l2
µS1 + µ0 S1 + µS2
por:
˚ ρm R ~ 0 ~r ¨ σm R ~ 0 ~r
~ R
H ~ = 1 dV 0 + dS 0
4π 0 r 3 0 r 3
V S
por: ˚ ¨
~ R
H ~ = 1 ρm~r 0
dV +
σm~r 0
dS
4π 3 3
V0 r S0 r
En la proposición 30 en la página se ha cambiado «El vector» por «En ausencia de polarización perma-
nente, el vector».
En la definición 57 en la página 174 se ha cambiado «Llamamos» por «En un medio sin imanación
permanente, llamamos».
En la demostración del teorema 26 en la página 219 se ha cambiado el último párrafo por «Esto último
es una coincidencia casual, es decir, nada de la teoría que hemos visto hasta ahora justifica que justo
el valor del cociente anterior sea la velocidad de la luz. Sin embargo; este hecho, unido a que las ecua-
ciones de Maxwell admiten una solución ondulatoria, apunta a que la luz es, precisamente, una onda
electromagnética. En otras palabras, hemos visto que las ecuaciones de Maxwell son compatibles con la
existencia de la luz.».
[1] Clases impartidas por Juan Pablo Martínez Jiménez en la asignatura de Electromagnetismo del curso
2018-2019 de la carrera de Física en la Universidad de Zaragoza.
[2] Del. [En línea]. En: Wikipedia, the free encyclopedia, 2019. Disponible en: https://en.wikipedia.org/
w/index.php?title=Del&oldid=877898711 [Consultado 03-02-2019].
235