This document was uploaded by user and they confirmed that they have the permission to share it. If you are author or own the copyright of this book, please report to us by using this DMCA report form. Report DMCA
en d es la totalidad de los vectores f(x, y, z) asignados a cada punto (x, y, 2 ) de una región R del espacio. Sea P(x, y , z) un punto no singular de una curva C en el espacio y Q(x + Ax, y + Ay, z + Az), otro punto de C. (Véase figura 3.8). Entonces los vectores de posición r y r+ A rd e P y Q son r = xi + yj + zk, r t Ar = (x + \ x ) i 4- (y + A y)j + (z + A z)k. {x, y, z ) en P en la dirección del vector unitario T tangente a C en P se define como ,es un vector que se define d = dr-ycf> + . Entonces
69 d
F¡ gura 3.7
U n e le m e n to d ife re n c ia l de área de la su p e rficie de una esfera.
PROBLEMA 3.27 Expresar el elemento diferencial de área de la superficie dS para una superficie S representada por (3.102)
z = z (x , y ).
S o lu c ió n :
Por (3.88),
Entonces, usando (3.96), d S = | r x x xy | dxdy
dxdy.
www.FreeLibros.me
(3,103)
r
58
Análisis vectorial
3.6
Z
Derivada direccional y gradiente
Un campo escalar
F ig u ra 3 .8
La d e riv a d a d ir e c c io n a l.
Sea
<>s
A*, y > Ay, z + ¡\z) As
_ ]jm AS-o
(x, y , z)
^
donde As es la longitud de arco de C desde P hasta Q. El gradiente de la función escalar
, , grad ó ;
fió . defc , i <■-r-1- J+ ~~~ k. dy dz
n 1AC^ (3.105)
PROBLEMA 3.28 Mostrar que la derivada/íireccional de
ds
-1 grad c¡> • T,
(3.106)
donde T es el vector unitario tangente a C. S o lu c ió n :
Por (3.68), el vector unitario T tangente a C en cualquier punto es T
d x d x . d y . d z --- j f — k. ds ds ds ds
(3.107)
de/) dx dx ds
(3.108)
------ — i i
.
Por cálculo,
dcj> ds
dtp dy d y ds
d(f> dz d z ds
que es la forma de un producto escalar de grad
r) s
grad 0 • T.
Introduciendo el operador diferencial V (léase nabla o del),
üx el gradiente de
www.FreeLibros.me
™ j+ - k , dy dz
(3.109)
59
Cálculo diferencial vectorial
grad
~
dx
i *
dy
i i-
dz
k.
( 3 .1 10)
Por (3.110), las propiedadas del gradiente son 7 (c<¡6) = c V<5, V (0 + 0) =
(3.111)
' V iA .
(3.112)
V(<W = 0 V 0 +
(3 -il3 )
donde 0 y 0 son funciones escalares difereneiables enalguna región del espacio y c es una constante. PROBLEMA 3.29 Solución:
V(
Verificar (3.113)
Por la definición de gradiente (3.110),
0A)i + o
x
a y
W>
r) rA \ / rí ?// r) rh\ w +f i ± \ i +L d ± + f i ± \ (9y (9y / \ dz <9z /
W
+f i ± \ i +L dx <9x / \
dx
<9y
dz
3— (0^ ) k
ó z
)
\clx
dy
<9z
= 4 rel="nofollow">V i/j + if/S7 cf).
PROBLEMA 3.30 Mostrar que la magnitud y dirección del gradiente V0 es independiente del sistema coordenado. Esto es, mostrar que la magnitud del gradiente IV 0 Ies igual al valor máximo de la derivada direccional de jp(jt, y, z) y su dirección es la de la rapidez máxima de crecimiento de la función 0. Por las definiciones de derivada direccional (3.106) y el producto escalar (1.25), la derivada direccional de 0 en el punto (x, y, z) es Solución:
^
ds
= V 0 - T = |V 0 | | T | eo s 0 = j V 0 | eo s 0,
(3.114)
donde 9 es el ángulo entre V 0 y el vector unitario T tangente a la curva C. Como -1 < eos 9 < 1, 30/ds es máximo cuando 0 = 0, i.e., cuando la dirección de T es la dirección de V0 y dcf> ds
= | V 0 I-
(3.115)
Así que la magnitud del gradiente es 90/9slmax y su dirección es la de la máxima rapidez de crecimiento de la función
PROBLEMA 3.31 Si 0(x, y, z) es una función escalar y V0 =£ 0 en un punto P en el espacio, entonces mostrar que V 0 es perpendicular a la superficie S definida por (f)(x,y, z) - c,
(3.116)
donde c es una constante. Usando (3.94) del problema 3.25, un vector unitario normal n para la superficie S definida por (3.116) es
Solución:
www.FreeLibros.me
Análisis vectorial
60
n =
d(fc . d ó . d
MV +/MV +(Éá
(3.117) IV 0
dx) \dy) \<9z Entonces, concluimos que V0 es perpendicular a S. S o lu c ió n a lte rn a : Si suponemos que una curva C del espacio, situada en una superficie S está representada por r(í) = x ( í ) i + y ( í ) j + z ( t ) k , entonces por (3.116) tenemos, para c constante, <£[x(0. y( t ) , z (0 ] = c. Diferenciando esto con respecto a t,
dx
dy
z'(O = V f r ' ( í ) = 0,
dz
(3.118)
donde por (3.40), el vector tangente a C en P es r '( 0 = x '( í ) i + y '( t) j + z' ( t ) k. Por otra parte, (3.118) implica que V $ l / ( f ) , i.e., es perpendicular a C en P.
Este razonamiento se puede aplicar a cualquier curva suave sobre la superficie S que pase por P. Por consiguiente, V0 es perpendicular a todas esas curvas en P, lo que puede cumplirse sólo si V0 es perpendicular a la superficie S. PROBLEMA 3.32 Hallar (a) la derivada direccional de
Solución:
=
r2 - T;
(2 - l ) i + (1 - l ) i + (1 - 0 ) k
i + k
r, - r,
= [(2 - l ) 2 + (1 - l ) 2 + (1 - O)2]^ ”
V2
V 2 Í+ V 2 k lV 2 ’
° ’ 72J
Por otra parte, por (3.110) con (p(x, y , z) dado, V 0 = 2 x i + 2 y j + 2 zk . Así que en (1, 1, 0) el gradiente es V<£ - 2 i + 2 j + Ok = [2 ,2 , 0], y su magnitud es |V<¿ I = (22 + 22 + 0)y> = 2\¡2. Por consiguiente, por (3.106) la derivada direccional de 0 en (1, 1, 0) es ( 2 ) ( 0) + ( 0 )
www.FreeLibros.me
:V2-
61
Cálculo diferencial vectorial
Esto significa que el valor de 0 se incrementa en %/2 por unidad de distancia mientras pasa de P( 1, 1 ,0 ) 3 0 (2 , 1, 1). (b) El valor máximo de 90/9s en (1, 1, 0) es I A 0 I= 2 \/2 y su dirección es la de A 0 = [2 ,2 ,0 ]. PROBLEMA 3.33 Hallar un vector unitario n normal a la superficie dada por z = x 2 + y 2 en el punto (1, 0, 1). S o lu c ió n : C om oz = x 2 4- y 2 , la superficie se define como 0 (x, y, z) = x 2 + y 2 - z = 0. El gradiente es
V<£ = 2 x i + 2 y j - k.
Así en (1, 0, 1), y> = 2 i - k. Por tanto, según (3 .1 1 7 ),el vector unitario es _ V 2i - k 2i-k n ~ | V 0 | = [22 + 0 + ( - l ) 2]^ = V5
PROBLEMA 3.34
V5 ’ ° ’
>/5
'
Para un vector arbitrario constante a, mostrar que V (a • r) = a,
(3.119)
donde r es el vector de posición. Solu ción :
Sea a = [a ,, a2, a 3]; puesto que r = [x, y , z], a • r = a l x + a 2y + a 3z.
Así que, d S d V (a - r) = — ( a tx + a 2y + a 3z) i + — (a jX + a 2y + a 3 z) j + — (at x + a , y + a 3z ) k ax dy dz = a l i + a 2j + a 3 k = a.
PROBLEMA 3.35
Si
- d t «s d4>. S olu ción :
Si r = xi + y¡ + zk , entonces su diferencial es d r = d x i + dy j + d z k .
Por (3.110), V0 =
k. Por consiguiente, tomando el producto escalar,
V
dx
dx +
dy
dy +
dz
dz = dá.
La identidad (3.120) se puede expresar como (3.121)
(dr-S7)
Si 0 =
dt.
(3.122)
dt S o lu c ió n :
Si 0 =
www.FreeLibros.me
Análisis vectorial
62
d(f>= * ± d x + * ± d y + ? * d z + ^ d t. dx
dy
dz
dt
Por los resultados del problema 3.35, dcf> = (V 0 ) - d r +
= d r ■Sjcf) +
PROBLEMA 3.37
dt
dt.
dt
Hallar V0 si 0 = r = l r l = ( x 2 + y 2 + z2) 1/2.
La superficie definida por 0 = r = (x2 + y 2 + z2) 1/2 = const es una esfera con centro en el origen (0, 0, 0). Por consiguiente, de acuerdo con el resultado del problema 3.31, V r es normal a la esfera y por lo tanto es paralelo al vector de posición r = [x, y, z]. Así, podemos escribir V r = kr. Por (3.120),
S olución:
dr = V r • d r = k r • d r = kr dr.
Así que k ~
r
,
y por tanto V r = — r = er, r
(3.123)
donde er es el vector unitario en la dirección del vector de posición r. S o lu c ió n a lte rn a : Por (3.110) el gradiente es V r = V (x2 + y 1 + z 2)I/j = —
dx
(x2 + y 2 + z 2)1/j i + —
dy
1 2x . 1 2 (x2 + y 2 + z 2)1^ =
x . y . — i +— r r
z r
j+
(x2 + y 2 + z 2)l/l j + — (x2 + y 2 + z 2) /2 k dz
_____ 2y_____ . ^ 1 _____ 2z_____ 2 (x2 + y 2 + z 2)1^2(x2 + y 2 + z 2)1^ —
k
1 = —r r
=erPROBLEMA 3.38
Si <¡>= 0(m), donde u = u(x, y, z ), entonces mostrar que V 0 = V ^ ( u ) = 0 '( « )V u .
S o lu c ió n :
Por (3.110), el gradiente es
www.FreeLibros.me
(3.124)
Cálculo diferencial vectorial
PROBLEMA 3.39
63
Si r = I r I = (x 2 + y 2 + z 2) 112 , hallar Vr" y V (l/r), donde n es cualquier
número real. S o lu c ió n : 3.37,
Sea
~ r ” - Entonces por (3.124) y el resultado (3.123) del problema
V (rn) - — ( r" )V r
n rn ~ 1 V r
dr
- n r" _ 1 e r n r " - 2 r.
(3.125)
Haciendo w = -1 en (3.125), v ( r- ) - - ^ r
3.7
(3.126)
El operador V
El operador vectorial diferencial V » - i - i f A j+ A l t , dx ay oz
[3.1091
no es un vector, sino un operador. Puede considerarse como un vector simbólico. Asi si 0 es un campo escalar,entonces <¡>V es un operador, mientras que V0 da la im portante función vectorial llamada gradiente. De modo similar, si f es una función vectorial diferenciable, entonces fV y f X V son operadores, mientras que V f y V X f dan importantes funciones escalar y vectorial respectivamente. (Vcanse secciones 3.8-9). Por tanto, multiplicando V a la izquierda da operadores, mientras que multiplicándolo a la derecha da importantes funciones escalares y vectoriales. En otras palabras, V sólo opera sobre lo que le sigue. 3.8
Divergencia de una función vectorial
Si una función vectorial es f = [ /,, / : . /, ], donde / , . f2, /, son funciones escalares, entonces su producto punto o escalar con el vector simbólico V es 7
f
f — i + — i + — k ) • (t, ¡ + l i \d x
ay
/,
di2
dx
dy
dz
J
/. k '
dz
Por consiguiente, una función‘vectorial f se transforma en una función escalar cuando se opera con V por la izquierda, lista función escalar se llama la divergencia de la función vectorial f y se escribe como div f. Asi que, . div f = V • f
di.
di
di.
dx
dy
dz
— L + —i + — í .
(3.12?)
Aunque la definición (3.127> no da significado físico o geométrico alguno al concepto de divergencia, da su forma de fácil cálculo. En el capítulo 4 se dan otras definiciones de divergencia con su interpretación física también
PROBLEMA 3.40
Si f y g son dos funciones vectoriales, mostrar que V-(f+g)=V-f+V-g.
www.FreeLibros.me
(3.128)
64
Por (3.127), V • (f + g) = ~ (f> + ¿i) + T " ^ dx dy
+ é j) +
+
dz
<9/1
d /2
d f 3\ ^ /
dx
dy
dz/
\d x
^é_3
dy
dz
= V • f + V • g. Si r es el vector de posición, hallar la divergencia V • r.
PROBL EMA 3.41
S o lu c ió n :
Si r = x\ + y j + zk, entonces p o r (3.127), V • r = — + — + — = 3. dx dy dz
PROBLEMA 3.42
(3.129)
Hallar la divergenciade f = x y z i + x 2y 2z') + y z 3k.
S o lu c ió n :Por (3.127), la divergencia
es
d /, d i, d /3 div f = V - f = — + — + — dx dy dz = —— (xyz) +
[
( y z 3)
(x2y 2z) +
^ = y z + 2 x 2yz
+ 3 y z 2.
(3.130)
Si
S o lu c ió n :
„ n ,v (3.131)
Verificar (3.131).
El gradiente es
dd> . dd) .
dd)
Vc6 = — i + — j +— k. dx dy dz Entonces,
dy \ d y j
dx \ d x / = d^0 dx2
d^> dy2
d z \d
d 2<¿> dz2
La divergencia del gradiente V-Vse escribe como V 2. Entonces V *(V<£) se escribe comoV • V
- V3 = dx2
dy2
+ —— , dz2
. f ü l + * * J ^ . £ ¿ + Ü á + Ü ¿. \d x 2 dyJ d z 2/ dx2 dy2 dz2
www.FreeLibros.me
(3.132)
13.131]
65
Cálculo diferencial vectorial
Una función escalar
(3.133)
PROBLEMA 3.44 Mostrar que la función l/r, donde r= Ir 1= (x 2 + y 2 + z 2) 1' 2 ps una función armónica, siempre que r¥= 0. Claramente l/r es continua, p u esx 2 , y 2, z 2 y r son continuas. Entonces por (3.131), el Laplaciano es Solución:
Por otra parte la primera y la segunda derivadas parciales de l/r con respecto a x son
dx
(x 2 + y 2 + z 2)~ ^ = - x ( x 2 + y 2 + z 2)~
— (x2 + y 2 + z 2T Vl = — [ ~ x ( x 2 + y 2 + z 2)~¡/’] dx2 dx = 3 x 2(x2 + y 2 + z 2)- ^ —(x2 + y 2 + z 2)- ^ = 3 x 2r -5 —r -3.
(3.134)
De modo similar, las segundas derivadas parciales de \/r con respecto a y y a 2 son (x2 + y 2 + z 2r l/3 = 3 y 2r - 5 - r - 3,
(3.135)
dy — (x2 + y 2 + z 2T Vi = 3 z 2r-5 - r " 3. dz2
(3.136)
Como x , y , z y r son continuas, las segundas derivadas parciales son continuas también si r # 0. Agregando (3.134-6), V2
= 3r-5 (x2
+ y 2 + z 2) - 3 r " 3= 3 r“ 5r 2
con tal que r =£ 0. Así que se satisface la ecuación
3.9
- 3 r-3
= 3 r - 3- 3 r " 3 =0,(3.137)
de Laplace y la función l/r es armónica.
Rotacional de una función vectorial
Si una función vectorial es f = [/, / 2/ 3 ], donde / , , f 2 y / 3 son funciones escalares con primeras derivadas continuas, entonces su producto cruz o vectorial con el vector simbólico V es
www.FreeLibros.me
A núlisis vectorial
66
Llamamos a esta función el rotacional (rot f) de la función vectorial f; i.e., . -á * ' .... rotacional f - V x f ,
V (3.139)
Debe observarse que V X f no es necesariamente perpendicular a f aunque tratamos a V como un vector simbólico. En el capítulo 4 se discute el significado físico del rotacional de f. PROBLEMA 3.45
Si f y g son dos funciones vectoriales, mosírar que V / (1 • g) = V y f
Solución:
4
(3.140)
V x g.
Sean f = [ / , , f 2, / 3] y g = [£ , , g2, g3 ]. Entonces, i'-*- g =• (/! + á ,) i 4- (f2 + é2) i + (h +
Por consiguiente, V x (f 4- g) =
¿ (/3 - g3)
d(( 2 4- g2)
dy
dz
d /,
dl2
dy
dz
¿éi I
dz
' ¿y
= V x f + V x g.
PROBLEMA 3.46
S o lu c ió n :
Si r es el vector de posición, hallar V X r.
Como r = xi + > j + zk, por (3.138) el rotacional es i d
V x r =
dx x
PROBLEMA 3.47 Solución:
j
k
d dy
d
y
(3.141)
0.
dz z
Hallar el rotacional de f = xyzi + x 2y 2z] +>>z3k.
Por (3.138), el rotacional es
V x f
i
j
k
<9 dx
d dy
d dz
xyz
x 2y 2z
yz3
**
d
£
dx
(x2y 2z) -
.
.
d
.
3
— (x y z ) - — ( y z 3) dz dx
( y z 3) - - f OZ
dy
(xyz )
( z3 - x 2y 2) i + (xy) j 4- (2 x y 2z - xz) k.
Si 0 es una función escalar con :>n segundas derivadas continuas, entonces el rotacional del gradiente de
Verificar (3.142)
El gradiente de 0 es
www.FreeLibros.me
(3.142)
Cálculo diferencial vectorial
,
dcf) . d(f> . dx dy
V<¿ =
67
dcf> dz
—;- i + —— J + —£-
k.
Por consiguiente, suponiendo que
j
k
A
JL
JL
dx
dy
dz
dcf)
dcf)
dcf>
dx
dy
dz
d 2cf)
d 2cf)
d 2cf)
dy dz
d z dy
d z dx
J+
dxdz
d 2cf)
d 2cf)
dxdy
dy dx
=0 pues las segundas derivadas son continuas y por lo tanto, d 2cf)
d 2cf)
dydz
dzdy
d 2cf) dzdx
d 2cf)
d 2cf)
d 2cf)
dxdz
dxdy
dydx
Ni la función vectorial es f = ff,. f t , f 3]. donde las com ponentes escalares / , , f 2, f 3 tienen segundas derivadas continuas, entonces la divergencia del rotacional de f es cero, V ■(V X Í) = 0.
(3.143)
Verificar (3.143).
PROBLEMA 3.49
S olución : Si f = fi'\ + f 2j + / 3k, entonces por (3.138) su rotacional es V x f
fdf3
d fi 2: 2\
/dfi
\dy
d z )/
; 1+ . \ d\(9z z
d í 3\
/ df2
dfí
dx)*
\<9x
dy
k.
Por lo tanto, suponiendo que f tiene segundas derivadas parciales continuas, por (3.127), la divergencia del rotacional es d
fdf3
d f 2\
V • (V x f) = — ( — i ----- i dx \d y
d 2L dxdy
+
dz )
d 2f. dxdz
+
d
(di,
dy
\dz
d
di.
—i
di j
d z \d x
d 21, dydz
/ dí2
dydx
dy
d 2í 2
d 2(l
d z d>
dz dy
= 0 pues las segundas derivadas son continuas y por lo tanto, d 2i 3
d 2(3
d 2i 2
d 2í 7
d 2f í
d2í l
dxdy
dydx
dxdz
dzdx’
dydz
dzdy
PROBLEMA 3.50
Mostrar que si B = V X A, entonces A no está unívocamente determinada
por B. S olución :
Si \p es una función escalar arbitraria, sea A' = A+V
(3.144)
Entonces por (3.140) el rotacional de A' es V x A' = V x (A + V
(3.145)
www.FreeLibros.me
68
Análisis vectorial
pues V X (Vi//) = O por (3.142). Así que B se puede expresar como B = V X A ', lo que indica que A no está determinada unívocamente a partir de B, pues i// es arbitrario.
3 .1 0 Operaciones con V y algunas identidades vectoriales Usando el operador V, hemos definido gradiente, divergencia y rotacional para obtener cantidades escalares y vectoriales. En esta sección consideraremos varías combinaciones del operador V con funciones escalares y vectoriales. Debe recalcarse que e! orden en el cual aparecen los símbolos en las expresiones que usan V es muy im portante, pues el operador V solamente opera en lo que le sigue. Así por ejemplo, (f-V )g é g (f-V ).
El lado izquierdo de esta expresión es un campo vectorial, mientras que el derecho es un operador. (Véase problema 3.51). PROBLEMA 3.51 Determinar el significado que puede asignarse a (f*V )0 y ( f -V)g. (Véase problema 3.35).
S o lu c ió n : Como V es un operador diferencial con respecto a las coordenadas en el espacio; sólo opera con lo que le sigue.Así que interpretamos (f*V )0 como producto vectorial del vector f con el gradiente V 0; i.e., (f-V)<¿ = f-(V<¿).
(3.146)
Se pueden obtener resultados similares interpretando a V como un vector simbólico; i.e., „
d
.
d
.
d .
V = —— i + —— j + —— k. dx
dy
dz
Entonces tomando su producto escalar a izquierda con f, f . V = /i J L + Í2 ± +Í3 ± . dx dy dz
(3.147)
A sí que,
(f-v)0= /, \
+ f¡ r + f, — 0
dx
dy
dx
dy
dz/
dz
= f-(V 0 ). Como Vg no está definida, usamos el operador diferencial (3.147) que opera sobre g para obtener ( f - V ) g = (i, \
dx
+ f2
dy
= /, ^ + f 2 dx dy
+ f3 7 r ) g dz]
I1 . dz
(3.148)
PROBLEMA 3.52 Hallar (f*V )0 y (f-V )g en (1, 1, 1) si f = - y i + x] + zk, g = 3 xyz2i + 2xy3j —x 2y z k y
S o lu c ió n :
Como
www.FreeLibros.me
69
Cálculo diferencial vectorial
V<£ = — (xyz) i + — (xyz) j + —— (xyz) k = y z i + x z j + x y k, dx dy dz tenemos (f • V) 0 = f • V (f rel="nofollow">= ( - y ) ( y z) + (x) (xz) + (z) (xy) = - y 2z + x 2z + x y z . Por ta n to e n ( l, 1,1), (f • V)0 = -1 + 1 + 1 = 1. Multiplicando f • V por g a la derecha da (f ■V) g = ( - y )
dx
( 3x y z 2i + 2xy3j - x2yzk)
+ (x) — (3xyz2i + 2xy3j - x2yzk) H d y\T _d + (z) — (3xyz2i + 2xy3j - x2yzk) dz = ( ~ y ) ( 3 y z 2i + 2 y 3j - 2 x y z k ) + (x) (3xz2i + 6xy2j - x2zk) + (z) ( 6 x y z i - x2yk) = ( - 3 y 2z 2.+ 3 x 2z 2 + 6xyz2) i + ( - 2 y 4 + 6 x 2y 2) j + (2 x y 2z - x 3z - x2y z)k . Entonces en (1, 1, 1), ( f - V ) g = ( - 3 + 3 + 6 )i + ( - 2 + 6 )j + (2 - 1 - l) k = 4 j = [0 ,4 ,0 ],
PROBLEMA 3.53
Si r es el vector de posición, m ostrar que (f ■V) r = f.
Solución:
(3.149)
Si r = x i + yj + zk, entonces
c ~ + í ci (ref • vV7 '\) r = fj dx
=
r
dy
+ ít j 7
az
+ f2i.+ f3k
= f.
PROBLEMA 3.54
Mostrar que ( d r - V ) f = di .
S o lu c ió n :
(3.150)
P or(3.147), , „ , d á r • V = d x -----h dy dx
d — + dz dy
d — . dz
Por consiguiente, (d r • V) f = c/x —— + dy —— + d z —— dx dy dz di . dt , dt = ---- dx + —— ay + —— az <9x dy <9z = d f.
PROBLEMA 3.55
Si f = f(x, y, z, t), entonces m ostrar que
www.FreeLibros.me
Análisis vectorial
70
d i = (dr.V)f + —
dt
(3.151)
dt.
Como f = f(x, y, z, t), tenemos por la definición de diferenciales,
Solución:
d t = — dx dx ,
dt
dt
,
di
,
dt
,
---- dy + ---- d z + ----dt dy ,
dz
dt
,
dt
dt dt
,
d x ------ h d y ------ h d z ------ h — dt dx
dy
dz
dt
( d r - V ) t + — dt. dt
PROBLEMA 3.56
Determinar el significado que puede darse a (f X V )0 y (f X V)g.
Como el operador V sólo actúa sobre lo que le sigue, interpretamos ( f X V ) i como el producto vectorial de f por el gradiente V0; i.e., Solución:
(f x V ) 0 = f x (V 0 ).
'
(3.152)
Se pueden obtener resultados similares interpretando V como un vector simbólico; i.e., V = Jd x~ 1 + 7d y~ J + Td z“ k> Como f X V es un operador,
f xV =
i
j
k
fi
u
U
d
d
d
dx
dy
dz
(^ \
~
dz
f 3
d dy
3 dy
2 di
k.
(3.153)
Asi que, (fx V )0
= f x (V 0 ).
No se asigna definición alguna a (f X V)g ya que éste es un tipo de operador diferencial con cantidades vectoriales.
Aunque no se asigne significado u cantidades tales como Vg y (f X V ig. ellas son generalizaciones de vectores y se llaman diádicos. Formalmente, un diádico es una generalización de un vector.
PROBLEMA 3.57
Si f y g son dos funciones vectoriales, mostrar que (f x V ) . g = f - ( V x g).
Solución:
Por(3.153),
www.FreeLibros.me
(3.154)
71
Cálculo diferencial vectorial
- (/,
— ) + (f, ^
dz
\
dy ¡
\d y
dz /
-t,
dx
\ dz
V
dz I
dy
dx
ir)'‘Á dx / \t d x1 - ¥ dy
= f • (V x g) pues ¡*g = g i , j ‘g = £2 y k - g = ^ 3.
Ahora consideraremos un grupo de identidades vectoriales bien conocidas que incluyen al operador V. Aun cuando estas identidades se pueden verificar por desarrollo directo usando las componentes de la función vectorial estableceremos estas identidades heurísticamente tratando V como un vector simbólico y también como un operador diferencial. Entonces las expresiones se manipulan de acuerdo con las fórmulas apropiadas del álgebra vectorial. En el resultado final V se da con su significado operacional.
PROBLEMA 3.58
Demostrar que V-(0f) = 0V -f + f-V0.
(3.155)
So l u c i ó n : Como V es un operador diferencial, la regla para diferenciación de un producto puede expresarse como V - ( 0 I ) = V * •(<£*) + V f - ( 0 f ) . El producto punto C 0 -(0f) significa que f es fijo y V opera sobre 0. De modo similar, V f*(0f) significa que es fijo y V opera sobre f. Como V *0 no está definido, V<¿-W>f) = f - V 0 , V f ( 0 f) = 0 V - f . Sumando las dos ecuaciones anteriores, V - ( 0 f ) = f - 7 0 + 0-Vf. Sol uc i ón al t er na: Usando componentes,
v • (0 r) = / - (0/,)+ ~ (0/2) +~ (0/3) dx
dy
dz
d(1 d
df , \
(
+d y + d z ) +[ fl
d
dcf)
dcf)
dx + 2 d y + 3 d z
- 0 V • f + f • V>. PROBLEMA 3.59
Demostrar que V x ( 0 f ) = 0V
Sol uci ón:
X
f + (V 0 ) x f = 0 V x f - f x V 0 .
Como V X (0f) = V0 X (0f) + Vf X (0f) y V X0 no está definido,
www.FreeLibros.me
(3.156)
72
Análisis vectorial
V(£ x (0 f) = (V>) x f = - f x V 0, V f x ( 0 f ) = 0V x f. Entonces, sumando estas dos ecuaciones, V x ( 0 f ) = (V 0) x f + >V x f = 0V x f - f x V 0 . Demostrar que
PROBLEMA 3.60
V • (f x g) = g • (V x f ) - f • (V x g). Solución:
(3.157)
Podemos escribir V - ( f x g) = Vf *( f x g ) + V g - ( f x g),
donde los subíndices tienen el mismo significado de antes; esto es, V f*(f X g) significa que g es fijo y V opera sobre f. Por (1.77), V f • (f x g) = g -V x f, Vg*( f x g) = - f • V x g. Así que, \ 7 * ( f x g ) = g- V x f - f - V x g . Verificar
PROBLEMA 3.61
V
X
(f
X
g) = f ( V - g ) - g ( v - f ) + ( g - V ) f - (f - v ) g .
(3.158)
Como V X ( f X g) = V f X (f X g) + Vg X ( f X g), aplicando la regla del término central para el triple producto vectorial (1.83) da
Solución:
V f x (f x g) = (g* V )f - g(V • f), Vg x(f xg) = f(V-g)-(f-V)g. Por consiguiente, V x (f x g) = f (V • g) - g(V • f ) + (g • V) f - (f • V) g. Las cantidades (g 'V )f y (f'V )g se definen en (3.148). PROBLEMA 3.62
Demostrar que
V (f • g) = f x (V x g) + g x (V x f) + (f- V )g + ( g- V) f .
(3.159)
Aplicamos la regla del término central para el triple producto vectorial (1.83) a f X (V X g), donde la función vectorial f es una constante. Entonces,
Solución:
f x ( V x g ) = V g( f - g ) - ( f - V ) g .
(3.160)
Vg ( f - g ) = f x ( V x g ) + ( f . V ) g ,
(3.161)
Por consiguiente,
e intercambiando f y g, V f ( g - f ) = V f ( f - g ) = g x ( V x f ) + ( g- V) f .
(3.162)
En consecuencia, usando (3.161-2), V ( f - g ) = Vf ( f - g ) + Vg ( f - g ) = g x (V x f ) + f x (V x g) + ( g - V ) f + ( f - V) g .
Si f es una función vectorial con segundas derivadas continuas, entonces el rotacional del rotacional de f es
www.FreeLibros.me
Cálculo diferencial vectorial
rot. (ro t f) = V x ( V x f) = V (V • f) -
73
(3.163)
= grad (div f) - V*f. PROBLEMA 3.63 Solución:
Verificar (3.163).
Por la fórmula para el triple producto vectorial (1.83), A x (B x C) = ( A • C) B - ( A • B) C;
[1.83]
entonces haciendo A = B = V y C = f, Vx(Vxf)=V(V-f)-(V-V)f = V ( V - f ) - V 2r.
En lugar de escribir (V*f)V, debemos escribir V( V*f) de modo que V opera sobre f. Entonces por (3.163), V 2f = V ( V - f ) - V x ( V x f) ,
(3.164)
que es la fórmula para calcular V2f. PROBLEMA 3.64 Solución:
Calcular V* [f(r)r] si r es el vector de posición y r = 1r 1.
Por (3.155) V- [ / ( r ) r ] = / ( r ) V - r + r - V[ / ( r ) ] ,
y por (3.129), V • r = 3. Entonces por (3.123-4), V [/ (r)] = / '( r ) V r = f '( r ) - = f ' ( r ) e r. r
(3.165)
por tanto, puesto que r* r = r2 , V • [f (r) r] = 3 / ( r ) +
PROBLEMA 3.65 Solución:
r
r . r = 3 f (r) + t i '(r).
(3.166)
Calcular V - ( r " _1 r).
Haciendo f ( r ) = r n~1 en (3.166), V . (r n_1 r) = 3 r n ~ l + (n - l ) r n_1 = (n + 2)rn“ 1.
(3.167)
Haciendo n = - 2 en (3.167), obtenem os, para r # 0, (3.168)
v . ( i ) . 0.
Ahora, por (3.126) ten em o sV (l/r) = - l / r } r; por lo tanto (3.168) indica que p a r a r l o , • .........\ .............................................................................................. ..... v . v ( i ) = v ( í ) = °.
PROBLEMA 3.66 Solución:
Calcular V X [f(r)r].
Usando (3.156), V X [ / (r) r] = / (r)V X r + [ y /( r ) ] x r.
www.FreeLibros.me
[3.1371
74
Análisis vectorial
Ahora, por (3.141), V X r = 0 y por (3.165), V [/(r)] = f ' ( r ) r/r. Por tanto, Vx[/(r)r] = ^ >
rx r= 0
(3.169)
pues r X r = 0 por (1.58).
PROBLEMA 3.67 Solución:
Calcular V • V [/(r)] = V 2/(r).
P o r(3 .1 6 5 ), V [ / (r)] = L Q
r.
Así que por (3.155), V • V [ / (r)] = V -
/'( r )
í ' (r)
V •r + r-V
Entonces por (3.129) y (3.165), obtenemos V • x/= 3 y V
1 d
> (0
1 V ( r ) ' |* —_ r
r . ~ r dr . r
-/"(r)--/'(r)
Así, puesto que r • r = r 2, V • V [/ (r)] = - f '( r ) + -
I r
/'( r ) r2
r •r
(3.170)
PROBLEMA 3.68 Solución:
C alcularV -V r" = V 2r ".
H aciendo/(r) = r n en (3.170), V - V r n = V 2r" = - — r n + — r" r dr dr2 = 2nrn " 2 + n(n - l ) r " - 2 = n(n + l ) r n ~ 2.
(3.171)
Haciendo n = —1 en (3.171), tenemos para r =£ 0, V 2 ^ J= 0, que es (3.137).
PROBLEMA 3.69
Mostrar que V- ( 0 Vi A - iA V 0) = 0 V 2v!' - ,/'V 20 .
Solución:
(3.172)
Haciendo f = Vi//en (3.155), V • (0V
(3.173)
entonces intercambiando 0 y i// en (3.173), V - ( ^ V 0 ) =
Sustrayendo (3.174) de (3.173) se obtiene V - ( 0 V 0 - < / f V 0 ) = 0 V 2( / ' - i / ' V 20 pues V0* Vi// = V i// • V0.
www.FreeLibros.me
(3.174)
75
Cálculo diferencial vectorial
3.11 Problemas suplementarios PROBLEMA 3.70 Sean A (í) = íi + í 2j + t 3k y B (í) = i + íj + (1 - f)k. En t = 1, calcular (a) A '(t), (b) (d/dt) [A (0*B (0] Y (c) (d/dt) [A (0 X B (í)], Respuesta: (a) [1, 2, 3], (b) 3, (c) [-1 , 1, 0],
Mostrar que si las terceras derivadas de la función vectorial f( t ) existen,
PROBLEMA 3.71
entonces
dt
(b)
-
dt
[ f ( í) V ( í ) t " X i ) 1 = [ f ( 0 f '( 0 r —(f)].
PROBLEMA 3.72 Si f(í) X f '(í) = 0, entonces mostrar que f(í) tiene una dirección fija. [Sugerencia: Sea f(t) = f ( t ) e (t), donde e(?) es un vector unitario. Deducir que e (t) X e' (t) = 0 y e (í ) * e '( í) = 0], PROBLEMA 3.73 PROBLEMA 3.74
(a)
(b)
Verificar la regla de la cadena para derivadas (3.19). Si f = uvwi + u w 2j - v3k y g = u 3¡ - u v w j + u 2w k, calcular
;— a en el origen, dudv d2£ en el punto (1, 1, 0). d2f x -----dv2
du
Respuesta: (a) 0, (b) [0, -3 6 , 0]. PROBLEMA 3.75
Trazar las curvas representadas por las siguientes funciones vectoriales:
(a)
r (í) = e o s í i + sen í j + sen, ík ,
0 < t < 2 tt-,
(b)
r(í) = a e o s í i + a sen t j + bt k,
6 rel="nofollow">0
(c)
r(f) = Í 3Í + t í , —oo < t < oa.
y
0 <í<<x>;
Respuesta: (a) Una elipse, (b) una hélice cilindrica. PROBLEMA 3.76 Hallar la longitud de las curvas representaáas por las siguientes funciones vectoriales: (a)
r(í) = e o s í i 4- sen í j + sen í k, 0 < í < 2n,
(b)
r(í) = í i + sen 277 í j + e o s 2t7 ík ,
(c)
r (í) = e o s 3 í i + sen 3 t j, 0 < t < 2n.
Respuesta:(a)
377, (b)
\Jl + 4772, (c)
0 < tM
1,
677.
PROBLEMA 3.77 Hallar el vector unitario T tangente a las curvas representadas por las siguientes funciones vectoriales en los puntos especificados: (a)
r(í) - eos í i + sen í j + í k,
en t = 77;
(b) (c)
r(í)
a ( í - sen í ) i t a ( l - eo s í) j, en cualquier punto t.
www.FreeLibros.me
76
Análisis vectorial
PROBLEMA 3.78 Hallar un vector unitario n normal en (0, 0, 0) a la superficie S representada por z = 3x2 + 4y 2 . Respuesta: [0, 0 ,- 1 ] . PROBLEMA 3.79 Hallar un vector unitario n normal en (1, 1, 1) a la superficie S representada p o r * 2 + y 2 - z - 1 = 0 . [2
Respuesta: —,
2
1 —
Hallar un vector unitario n normal a la superficie S representada por la ecuaciones paramétricas x = u eos v, y = u sen v, z = z(u). PROBLEMA 3.80
Respuesta:
1 = = [~z'(u) eos v, - z ' ( u) sen v, 1], V I + [(z O )]2
PROBLEMA 3.81 Hallar la derivada direccional de
Respuesta: 18/VT4. PROBLEMA 3.82
(a)
Si 0(x, y, z) = x y + y z + zx , hallar
V 0 en (1, 1, 3),
(b)
en (1 ,1 , 3)
en la dirección de [1, 1, 1],
ds
(c) la derivada normal t —= V 0 * n en (1, 1, 3), donde n es un vector unitario normal an a la superficie S definida por una constante 0(x, y , z). Respuesta: (a) [4, 4, 2 ], PROBLEMA 3.83
(b) 1 0 /\/3 , (c) 6.
Hallar la divergencia y el rotacional de f = (x - j>)i + (y - z)j + (z - x ) k
Y g = C*2 + y z ) i + (y 2 + z x )j + (z2 + x^)k . Respuesta: V - f = 3, VX f = [1, 1, 1], V-g = 2(x + y + z) y V X g = 0. Si 0 = 3x2 - y z y i = 3 xyz2i + 2 xy3j - x ^ z k , hallar en (1, -1 , 1), V - f , (c) V x f, (d) f - V 0 , (e) V • (
PROBLEMA 3.84
(a)
V 0,
(b)
(g) V 20. Respuesta: (a) [ 6 ,- 1 ,1 ] , (b) 4, (f) [ - 3 , - 4 1 ,- 3 5 1 , y (g) 6. PROBLEMA 3.85
(c)
[ - 1 ,- 8 ,- 5 ] ,
(d)
-1 5 ,
(e)
1,
Si V X f = 0, donde f = (xyz)m (x"i + y n) + z"k), mostrar que m = 0
o bien n = -1 . PROBLEMA 3.86
Para un vector arbitrario constante a, y el vector de posición r, mostrar
que (a)
V • (a x r) = 0,
(b)
V x (a x r) = 2 a,
(C)
v M
(d)
V (a-V - ) + v
. V
PROBLEMA 3.87 r, m ostrar que
x H
X
^ O
,
( a x V - )= 0.
Para cualquier función vectorial f diferenciable y el vector de posición
www.FreeLibros.me
■K
77
C ilc u b diferencial v ecttritl
V - ( f x r) = r - ( V x f). Mostrar que si 0 y \p son armónicas, entonces 0 + \p y c 0, para una constante arbitraria c, son también armónicas. PROBLEMA 3.88
PROBLEMA 3.89
Si las segundas derivadas de las funciones 0 y 0 existen, entonces
demostrar que (a)
V2( 0 ^ ) = 0 W
+ 2V 0 • V
(b)
V / £ ) = ’AW» ~ \ ^
_
PROBLEMA 3.90 Si u es un vector unitario constante y f es una función vectorial diferenciable, m ostrar que u * [V (f *u) - V X ( f X u) = V ‘ f.
PROBLEMA 3.91
Si f = 2zi + x 2j + x k y 0 = 2x2.y2z 2 , hallar ( f X V) 0 en el punto
(1,- 1, 1). Respuesta: [ - 8 ,- 4 , 4]. PROBLEMA 3.92 Se dice que el campo vectorial f es solenoidal si f es diferenciable y V • f = 0; se llama irrotacional si f es diferenciable y V X f = 0. Mostrar que si f y g son irrotacionales, entonces f X g es solenoidal. PROBLEMA 3.93 Si 0 y \p son funciones escalares que tienen segundas derivadas continuas, entonces m ostrar que V0 X Vi// es solenoidal. PROBLEMA 3.94
Si 0 es armónica, entonces mostrar que V0 es solenoidal e irrotacional.
PROBLEMA 3.95
Si cf = V0, donde c es una constante, entonces mostrar que f • V X f =
PROBLEMA 3.96
Mostrar que si f es cualquier función vectorial diferenciable, entonces
0.
(f • V) f = i V (f • f) - f x (V x f).
PROBLEMA 3.97 Mostrar que para funciones escalares difereneiables cualesquiera u (x ,y , z), v(x, y , z) y w(x, y , z),
V v Vw] =
du
d_u
dx
dy
dv dx
dv dy
dv
dw
dw d i
dw dz
dx
du dz dz
v, w y se denota por (u, v ,w ) .(x, y , z ).
d (u, v, w) d (x, y, z)
PROBLEMA 3.98 Demostrar que una condición necesaria y suficiente es que el jacobiano de u, v, w se anule para u, v, w para ser funcionalm ente dependiente: esto es, u, v, w satisfacen la ecuación F(u, v, w) = 0.
Verificar que u = x + y , v = x - y + z , y w = (2x + z) 2 + (2y - z ) 2 son funcionalmente dependientes.
PROBLEMA 3.99
www.FreeLibros.me
4
CAPITULO
CALCULO INTEGRAL VECTORIAL 4.1
Integrales de línea pura a -Z t < b, se representa por
En la sección 3.4, una curva
diferenc ? d i a lo largo de C es dt -
dx i +
dy j + dz
(4.1)
k.
Las integrales que incluyen vectores de desplazamiento diferencial d t se llaman integrales de linea. Consideremos las siguientes integrales de línea a lo largo de una curva C que puede ser abierta o cerrada:
J í
<¡ rel="nofollow"> d t ,
(4.2)
í dt,
(4.3) jppW SÉ É l
llllM lM j
(4.4)
Con un escalar i
la integral de línea (4.2) se reduce a * d ,' l
1 | j j j | p ; j|| | | |
■
i p
| £
ó dy + k
Para reducir (4.2), hemos usado ta relación i <¿ dx - i I 'c
é dx
e tc .,
(4.o)
constantes. Cuando ía curva del espacio C forma una trayectoria cerrada, (4.2) se escribe como
En particular, si C es una curva plana cerrada simple, y la dirección de integración puede describirse, los símbolos que se usan son £ ¡C
www.FreeLibros.me
<¿>d r,
79
Cálculo integral vectorial
La primera integral indica movimiento a lo largo de la curva cerrada C en la dirección positiva o sea contraria a la de las manecillas del reloj; i.e., el movimiento a lo largo de C es tal que la región que encierra está siempre a la izquierda. La segunda integral indica movimiento en la dirección negativa o sea la misma en que se mueve un reloj. PROBLEMA 4.1 Si 0 = xy, calcular
I >dx desde (0, 0, 0) hasta (1, 1, 0) a lo largo de Je (a) la curva C especificada por y = x 2, z =0 y (b) la recta C que une (0, 0, 0) y (1, 1, 0). Solución:
(a) La curva C puede representarse paramétricamente por x = t,
y = t 2,
z = 0,
0 < t < 1.
Por consiguiente, a lo largo de esta trayectoria, dx
dx i + dy j + d z k = dt i + 2t dt j
y >= (t) (t2) = t3 . Entonces,
J. d . 1 . 0 ) (0,0,0)
r1 4>dr = I f3(i + 2 t\)d t “'o = i f ’ f3 d t + j f 1 2 1* dt •Jo Jo =
1 . 2 . i + J. 4 5
(b) Por el resultado del problema 3.15, la línea recta C que conecta (0, 0, 0) y (1, 1 ,0 ) puede representarse paramétricamente por x = t, y = t, z = 0, 0 < Í < 1 . Por lo tanto, a lo largo de esta trayectoria, dr = d x i + d y j + d z k = d t i + d tj y 0 = (f)(0 = t2 . Entonces, /- ( i , i , o) r\ I '
Jo 1 .
Jo 1.
= 3^ 1 + 3^ J-
Como muestran los resultados del problema 4.1, los valores de este tipo de integral de línea dependen generalmente del camino de integración.
PROBLEMA 4.2
Calcular
z = 0. So l u c i ó n : Por el resultado del problema 3.13, el círculo de x 2 + y 2 =a2 , z = 0 puede representarse paramétricamente por x = a eos t,
y = a sen t,
z = 0,
0 < t < 2n.
Por consiguiente,
www.FreeLibros.me
r 80
Análisis vectorial
dr = dx i + dy j + d z k = - a sen t dt i + a eos t dt j. Entonces, 2 77-
( - a sen t i + a eos f j) c/í 2 77
/> 2 77
- a sen í dt + j I
= i I j0
,
a eos t dt ^0
= i a eos íI g77"+ j a sen í IQ77 = i (a - a) +j(0 - 0) = 0.
(4.8)
La integral de 1ínea (4.3) se llama a veces la integral escalar de linea o simplemente integral de linea de un campo vectorial f. La integral de línea de f alrededor de una curva cerrada C se llama la circulación de f alrededor de C y es c i r c f 3<£ f-dr. Je PROBLEMA 4.3
(4.9)
Si T es el vector unitario tangente a lo largo de una curva C, mostrar
que
r
[ ■d r =
•Je
r
f - T ds.
(4.10)
''o
Solución: Si en (3.38), t = s es la longitud de arco de C medida desde algún punto fijo, entonces C puede representarse por r(s) = x ( s ) i + y ( s ) j + z (s)k ,
a < s < b.
(4.11)
Por (3.68) el vector unitario tangente T a lo largo de C es dr ds Entonces, por (4.12), d r = — ds = T ds. ds
(4.13)
Por lo tanto, J~
f 'J'í d s =
f
f -T d s ■
(4 J 4 )
Como el producto vectorial f • T es igual a la componente de f en la dirección de T, i.e., la dirección de C, (4.10) muestra que la integral de línea de f es equivalente a la integración de la componente tangencial de f a lo largo de C con respecto a la longitud de arco. PROBLEMA 4.4
Si f = f •Je
Solución:
(x, y , z)i + f 2(x, y , z)j + f 3(x, y , z)k, mostrar que f-cfr= r
( f í dx + f2 d y + { l d z ).
Je
Como d r = dx i + dy j + d z k, f • d r = f, dx + í2 dy + f 3 dz,
www.FreeLibros.me
(4.15)
81
Cálculo integral vectorial
y por consiguiente,
í
f-dr=í
•Je
dx + (2 dy
(/,
l z dz).
+
-Je
Si f = r i + y'¡ + x y z k , calcular I f ' d r desde (0, 0, 0) a (1, 1, 1) a lo
PROBLEMA 4.5
•Je
largo de (a) una recta que conecta estos dos puntos y (b) un camino C, como se muestra en la figura 4.1, que consiste en tres segmentos de recta C it C2, C3 que ligan estos dos puntos vía ( 1, 0, 0 ) y (1, 1 ,0 ).
(1, l, l)
So l u c i ó n : (a) Por el resultado del problema 3.15,1a rectaC que conecta (0, 0, 0) y ( 1, 1, 1) puede representarse paramétricamente por x = t,
y = t,
(t, 0, 0)
z=t.
Por consiguiente, la curva C está representada por r = t i + t j + t k,
Figura 4.1
( 1, 1, 0) S o lu c ió n al problem a 4 .5 .
0 < t < 1.
'i Entonces, a lo largo de esta trayectoria, dr = dt i + dt j + d t k ,
f = t 2i + t j + í 3k,
asi que,
f
( t 2 + t + t l) d t = - + - + - = — 3 2 4 12
, d" J
(b) En este caso, j
!• d r = í
Jn
f-dr+í
f-dr+l
J r. .
[• dr.
J e . . Je.
Para C v dr = dx i y f = x 2i; por consiguiente, í d r = f 1 x 2 dx = - ,
f 1
í
3
Para C2, dr = dy j y f = ( l ) 2i + yj = i + y ¡ ; por consiguiente,
L ' - d t - J ' y d , = r1 Para C3, dr = dz k y f = (1) i + ( l)j + (1) (l)z k = i + j + zk; por tanto,
f
i,
f- d r =
f
Jo
z dz= - .
2
Así que,
I Como lo muestran los resultados del problema 4.S, los valores de las integrales de línea de un campo vectorial dependen también en general, de la trayectoria de integración. PROBLEMA 4.6
Si f = x i + 2yj + zk, calcular I f ‘dr desde, (0, 0, 0) hasta (1, 1, 1) a •Je
lo largo de (a) una recta que conecte estos dos puntos y (b) una trayectoria C, como
www.FreeLibros.me
82
Análisis vectorial
se muestra en la figura 4.1, que consiste en tres segmentos de recta C j, C2, C3 que ligan estos dos puntos vía (1, 0, 0) y (1, 1, 0). (Cf. prob. 4.5). Solución:
(a) Por los resultados del problema 4.5, d r = dt i + dt j + dt k.
Como la representación paramétrica de f es f = íi + 2/j + /k, la integral es f
f •dr = f
Je
(t + 2t + t ) d t = f
-Jo
4t
dt = 2.
*'o
(b) Para una trayectoria C, consistente en segmentos lineales C x, C2, C3 f
f •dr = f
•'C
f •dr + f Je t
Para Ci, dr = dx i y
f-dr+ f
Jc2
f-dr.
Jc3
f = x i; por consiguiente,
f
f
f • dr = J C iJo
x dx = - , 2•
Para C2, dr = dy j y f = i + 2 j j; por consiguiente, r f •dr = f -'Cj ^0
2y dy = 1.
Para C3, dr = dz k y f = i + 2j + zk; por consiguiente,
Así que, í
1
f • d r = —+ 1 + - = 2 .
2
2
L1 resultado del problema 4.6 muestra que en este caso el valor de la integral de línea es el mismo para dos trayectorias diferentes que unan (0. 0, Oí u
j
j
f-dr
X dx * 2 y dy + z dz .
tenemos
'
^£0,0,0) ... *— 4. y. 2 + -2
|(i .i .n
2 | ( 0 i 0 ,0)
- 1- • 1, f 1 2
2
1 2> lo que muestra que el resultado depende sólo de los puntos (0. 0. 0) v i l . 1, i ) y es independiente de la trayectoria de integración. La discusión anterior muestra que si la integral / , d x + / 3 dv + / 3 dz de
www.FreeLibros.me
Cálculo integral vectorial
83
(4.15) es la diferencial total de alguna función, entonces la integral es la misma para todas la trayectorias que unen los puntos. Sin embargo esto no es cierto siempre y las condiciones suficientes para que la integral de línea de un campo vectorial f sea independiente de la trayectoria se discuten en la sec. 4.10, Hallar (j) r • dx a lo largo del círculo C representado por x 2 + y 2 =a2,
PROBLEMA 4.7
z = 0. Con referencia al problema 4.2, la ecuación paramétrica de la curva es
Solución:
r = a e o s t i + a sen t j,
0
Por consiguiente, dx = - a sen t dti + a eos t dtj. Entonces X- dx = - a 2 eos t sen t dt + a 2 senf eos t dt = 0 dt = 0. Por lo tanto, r - cí r = 0.
(4.16)
i Si f = f xi + / 2j + / 3k, m ostrar que
PROBLEMA 4.8
f
f x d x = i
í
Je
(f2
dz -
f3
dy) + j
í
Je
( (}
d x - f l dz)
Je (f¡ d y - f2 dx).
+k í
(4.17)
Je Solución:
Por (2.27), el producto vectorial de f con dx es
r
f x dr dx
dy
dz
= i ( / 2 dz - f3 dy) + j ( / 3 dx - f í d z) + k
dy - l2 dx).
Por lo tanto, la integral es f
fxrfr=i
Je
í
(f2 d z - / 3 dy) + j í
Je
PROBLEMA 4.9
( f3 dx -
dz) + k í
Je
(f¡
Je
X dx a lo largo del círculo representado p o r * 2 + y 2 = « 2,
Calcular
z = 0. S o lu ció n :
Con referencia al problema 4.2, C se representa paramétricamente por r = a e o s t i + a sen t j ,
y por consiguiente, dx = - a sen t d t i + a eos t d t j. i r x dx =
a eos t -a sen t
a
0 < t < 2tt, Así que j
k
sen t
0
dt a eos t dt
j \ G>gs e
0
slt)
= k (a 2 e o s 2 f + a 2 sen2 t) dt = a 2 d t k. CA
De modo que integrando sobre C,
www.FreeLibros.me
A nálisis vectorial
84
r r2rr r2rr {prxdr- I ka2dt= ka2 I dt = 2 n a 2k, c *'o Jo
(4.18)
lo que muestra que el resultado es un vector en la dirección de z con magnitud 2na2, que es el doble del área de un círculo. (Véanse problemas 4.68 y 4.69).
4.2
Integrales de superficie
Una superficie S se representa por r(u, v) = x( u, v )i + y ( u , v ) j + z ( u, v)k,
[3.79]
donde r es el vector de posición y u, v son parámetros. (Cf., sección 3.5) El elemento diferencial dS del área de la superficie es d S = t ux r v dudv.
[3.95]
Como se muestra en el capítulo 3, (3.95) se puede escribir como d S ~ n dS,
[3.97]
donde n es el vector unitario normal a la superficie en el punto que corresponde a las coordenadas u y v, de modo que n = I ru x r v |
[3.80]
dS = | ru x r v| dudv.
[3.96]
Las integrales que incluyen el elemento diferencia] de área dS se llaman integrales de superficie. Consideremos las integrales de superficie jjó
dS,
(4.19)
jj f ' d S ,
(4.20)
s
:W'::
JJ'
f
X
dS,
(4.21)
s
cada una de las cuales esta sobre una superficie S que puede estar abierta o cerrada. S iS es una superficie cerrada, las integrales de superficie se expresan
P ap superficies cerradas, es común suponer que la dirección positiva cié la normal está . .i. ■ ■■ PROBLEMA 4.10
Calcular
JJ' dS sobre la parte z >0de la esferax 2 +y 2 +z2 = a2. s
So l u c i ó n : Por el resultado del problema 3.26, la esfera*2 + y 2 + z 2 = a 2 puede representarse paramétricamente cambiando d en u y tp en v. Así que, r (u, v) = a sen u eos v i + a sen u sen v j + a eos u k,
(4.22)
donde 0 < « < 7 r , 0 < y < 2n. El elemento diferencial d S del área de la superficie es
www.FreeLibros.me
85
Cálculo integral vectorial
V
d S = ru x rv dudv -
a 2 sen2 u e o s v d udv i + a 2 sen2 u sen v d udv j (4.23)
+ a 2 sen u e o s u d udv k.
Como S es la parte de la esfera donde z > 0, el recorrido u debe cambiarse a 0 < « < ít/2 pues eos u < 0 para 7t/2 < m < n. Por consiguiente, -7 7 /2
ru x rv dudv
í - n
27t í'Tr/2
■t
i 277
+k
n Como
f
¡ ¿o
r2rr rn/2 a 2 sen2 u eos v dudv + j I I a 2 sen2 u sen v dudv do •'o
/* 7 7 /2
¡ Jo
sen u eo s u dudv.
2v /'2-ít e o s v dv = ( ( sen v d v = 0,
í*2 2 7t7t
•'O
J
r
/* 7 T /2
a sen" do ^2 77
“ '0
•'O
{*2tt o
rn/2 eos v dv I sen2 u du = 0, *'o
r2rr f 77/2 = a2 | sen v d v sen2 do •'o
í*7T/2
•'o
r2rr ('■n/2 = a2 I dv | sen a eos ^0 Jo
r 27T r ^ 2 , I a sen u eos u ''o «'n
r TT/2 = 2 na2 I sen u eo s u du -'0 r w/2 = 2 77a 2 I sen u d (sen u) do 77/2 1 = 2 77a 2- — sen2 u 2 1 0 „ , 1 = 2 77a • — 2
Por consiguiente,
f f
dS =
77
(4.24)
a 2 k,
lo que muestra que el vector resultante está dirigido en la dirección de z y tiene magnitud na2 , i.e., el área rodeada por la intersección de la esfera con el plano xy. (Véase figura 4.2). Fig ura 4 .2
La esfera del problem a 4 .1 0 .
PROBLEMA 4.11
Mostrar que sobre la esferax 2 + y 2 + z 2 =a 2 , (4.25)
J~f dS= °
www.FreeLibros.me
86
Análisis vectorial
So l u c i ó n : Si S x y S2 son las partes de la esferax 2 + y 2 + z2 = a2 para las cuales z > 0 y z < 0 respectivamente, entonces
s
s,
Por los resultados del problema 4.10. d S = ;7 a2 k.
(4.24)
II
S,
Para , puesto que z
ru x r v tfuc/v,
77/ 2
y por cálculos similares, r2n -2 7 ? rr nn rr2.n ín rrnn I I a 2 sen2 u eo s v dudv = I I a 2 s in 2 u sen v cíucív = 0, Jo n/ 2 o 77/2 '277
n
77 J’n/ TT 2
a 2 sen u e o s u c/ucfv = 2 77a 2 • — sen2 u = 27ra2 [——] = —n a 2. 2 77/2 '
Por consiguiente, d S - - n a 2k.
(4.26)
I I s,
Sumando (4.24) y (4.26), J J d S = 77a 2 k — 77a2 k = 0. s La relación (4.25) es válida para cualquier superficie cerrada. (Véase problema 4.50).
La integral de superficie (4.20) de un campo vectorial f se llama el flujo de f a través de S y se representa por
JETs PROBLEMA 4.12
(4.27)
Si S se representa por r(u, v), mostrar que el flujo de f a través de S es
JJ f • d S = JJ [f ru rv] dudv,
(4.28)
S donde [f ru r j = f mru X rv es el triple producto escalar y R uv es la región en (u , v) que corresponde a S. Solución:
El elemento diferencial de área de la superficie es d S = ru x
Por consiguiente el flujo de f a través de S es
www.FreeLibros.me
87
Cálculo integral vectorial
JJ f - d S = JJ f • ru x rv dudv =, JJ [f ru rv] dudv. S
R„
PROBLEMA 4.13
Si f = eos u eos v i + eos u sen v j - sen u k, donde 0 < u < tt y
JJ'
0 < v < 2v, calcular | | f 'd S sobre la parte de la esfera x 2 + y 2 + z2 = a2 para la cual z > 0.
So l u c i ó n : son
Por el problema 3.26, las com ponentes del vector que genera la superficie ru = a e o s u e o s v i + a eo s u sen v j - a sen u k, rv = —a senu sen v i + a sen u eo s v j.
De acuerdo con la definición de triple producto escalar (2.53), eo s u eo s v [f ru rv]
eo s u sen v
a eo s u eo s v a eo s u sen v —a sen u sen v
—sen u
= o,
—a sen u
a sen u e o s v
0
Como las filas primera y segunda del determinante son proporcionales, su valor es cero. En consecuencia, por (4.28)
JJf • c/S = 0.
(4.29)
s
PROBLEMA 4.14
•JJ"
r - dS , donde S es la superficie de la esferax
Hallar
+y
+ z -a .
s
Solución: En el punto (x, y , z) de la superficie de la esfe ra l, el vector de posición r = x i + y ¡ + z k y el vector unitario exterior n normal a la superficie S apuntan directamente hacia el lado opuesto del origen. Así que, n = er =
(4.30)
r Entonces, para puntos de la superficie,
¿7 S<-
y como el área de la superficie de una esfera es 4na2 por el problema 3.26,
JJ r • d S = JJ r • n dS = a JJ dS = a (4 n a 2) = 4 na*. PROBLEMA 4.15 Si f = f i ( x , y, z )i + f 2(x, y , z)j + f-¡{x, y , z)k, mostrar que la integral de superficie de f se puede expresar como
JJ f • d S = JJ f • n dS ^
^
donde R yz, R zx y R
= ±
JJ f ¡ d yd z ± JJ f2 dzdx ± JJ f3 dxdy, R yz
^ zx
R xy
son las proyecciones de S sobre los planos y z, zx y xy ,
www.FreeLibros.me
(4.31)
/
88
Análisis vectorial
respectivamente y los signos de las integrales de la derecha de (4.31) están determinados por los signos de (n -i), (n *j) y (n *k ), respectivamente. Si a, y y son los ángulos entre los ejes coordenados rectangulares y el vector unitario n, entonces, por el resultado del problema 2.18, S o lu ció n :
n • i = eos a ,
n - j = cos^3,
n • k = eos y.
[2.50]
Por consiguiente, por (2.65), n puede escribirse como n = (n • i) i + (n • j) j + (n • k) k
= (eos a ) i + (eos /3)j + (eos y)k,
(4.32)
y por consiguiente, d S = n d S = (eos a ) d S i + (eos /3) c/S j + (eos y ) d S k.
(4.33)
Ahora podemos escribir (véase figura 4.3), dS (n • i) = dS (eos
a ) = ± dydz,
dS (n • j) = dS (eos
/3) = + dzdx,
dS (n • k) = dS (eos
y) = ± dxdy,
donde los signos están determinados por los de (n *i) = eos a (n* j) = eos |3, (n • k) = eos y, respectivamente. Entonces, Figura 4 .3
S o lu c ió n al p ro b le m a
d S = n dS = (± dydz) i + (± d z d x )i + (± dxdy) k.
4 .1 5 .
(4.34)
Así, si f = / j i + / 2j + / 3k,
JJ f • d S = + JJ f l dydz ± JJ l2 d zd x ± JJ /3 dxdy, donde R y z , R zx y R xy son las proyecciones de S sobre los planos yz, zx y xy , respectivamente, y los signos están determinados como en (4.31).
PROBLEMA 4.16
Calcular
JJ r ‘dS, dondeS es la superficie del cubo limitado porjc = 0, s
jc = \ , y = 0 ,y = 1, z = 0 , z = 1 como se muestra en la figura 4.4. El vector normal exterior unitario n y el vector de posición r de los puntos de la superficie del cubo están dirigidos hacia el lado opuesto del origen. Solución:
Dividimos el cubo en las áreasS¡ ~ S 6.
Sobre Sj para AOCB, z = 0, n = - k ; por consiguiente, r - n = (x i + y i + z k ) - ( - k ) = - z = 0. Sobre S 2 para AFEO, y = 0, n = ~ j;p o r consiguiente, r -n = ( xi + y j f z k ) - ( - j ) - - y - 0. 4 .1 6 .
Sobre S 3 para O E D C , \ - 0, n = —i; así que, r • n = (x i f y j f z k) • (—i) Sobre S 4 para AFGB, x = 1,
n = i ; por
r • n - (x i + y j f z k) • (i) - x Sobre S 5 para BGDC, y = 1, r -n
n = j ; por (xi f y j i z k ) - ( j ) - y
Sobre S 6 para FEDG, z = 1, n = k; así que,
www.FreeLibros.me
- x - 0. consiguiente, 1. consiguiente, 1.
89
Cálculo integral vectorial
r- n = ( x i + y j + z k) • (k) = z = 1. Entonces, como las integrales sobre S¡ ~ S 3 son cero y sobre S 4 ~ S 6 son uno,
JJr-cíS = JJr- n dS + JJ r - n dS + JJ r - n dS +JJ r - n dS + JJ r • n dS + JJ r - n dS f J j s . f f c s . J J é s S4
=
1
+
S 5
S 6
1+ 1
= 3.
JT<
Si f = zi + y j + 2zk, calcular | J f*c/S, donde S es la parte de la
PROBLEMA 4.17
superficie del paraboloide x 2 + y 2 = 1 —z para el cual z > 0. (Véase figura 4.5.) So l u c i ó n : En el problema 3.31, un vector unitario n normal a la superficie 4>(x, y , z) = 0 es V0
[3.117]
ív ^ T ’ donde V0 es el gradiente de
Fig ura 4 .5
p ro b le m a 4 .1 7 .
Por consiguiente, (3.117) viene a ser V<£
2xi + 2yj + k
|V 0 I
( 4 x 2 + 4 y 2 + 1)V’
Los productos punto de n con i, j, k son n •i ( 4 x 2 + 4 y 2 + 1)V=’ 2y
n-j
(4 x2 + 4 y 2 + i f i ’ n •k
1 (4 x 2 + 4 y 2 +>¿)1/2'
Así que n*i y n -j son positivos o negativos dependiendo de six e y son positivos o negativos, mientras que n • k es positivo para todos los valores de x e y . Por consiguiente, \ si se usa
Sí ds= r nds 1'
S
I f S
^■
= * I f Ryz
U na pa rte de la su p e rficie del p a ra b o lo id e del
fidydz±Sf4
d zd x ■ / / /3 dxdy
[4.31]
^ zx
para calcular la integral, S debe subdividirse en dos superficies. Para calcular el primer término, debe haber una superficie para x > 0 y otra para x < 0.
www.FreeLibros.me
Análisis vectorial
90
De modo similar, para calcular el segundo término, debe haber una superficie para y > 0 y otra para y < 0. El primer término
J í
/ i dydz de (4.31) se calcula como sigue. P arax > 0,
Ry z
x = (1 - y 2 - z)1/2 sobre S; por consiguiente,
JJ íi dydz = JJ x dydz R yz
R yz
1
r l-y
■f
J y = -1
fz —0
1/ (1 - y 2 - z) 2 dzdy
| j f 1 (1 - y 2)V2 dy
2 3
3 8
— • —
77
1 —
77.
4 P arax < 0, x = - (1 - .y 2 - z ) 1/2 sobre 5; escogiendo el signo negativo,
- JJ f¡ d yd z = — JJ x dydz r l~y
'i
v
(1 — y 2 — z) 2 dzdy
y= 1 4
=
El segundo término
—
77.
/ 2 dzdx de (4.31) se calcula como sigue. Para y > 0,
y = (1 - x2 - z ) '2 sobre 5; por consiguiente,
JJ f2 dz dx = JJ y dzdx ' ' 1~*
y
1/ (1 - x2 - z) 2 c/zcfx
1 =
—
77
4 pues la integral es la misma que para el caso/?yz. Para.y < 0 ,y = -(1 - x 2 - z ) 1/2 sobre 5; escogiendo el signo negativo,
- JJ í2 dzdx = — JJ y dzdx ^ZX
^ZX
■r f Jx = -1
1 —
4
www.FreeLibros.me
77.
J z = 0
(1 - x2 - z) 2 dzdx
Cálculo integral vectorial
91
Para el tercer término de (4.31), como n • k > 0,
Jj"f3dxdy = j~J 2 z dxdy = 2 í í '1
R
i\Xy
/V 1 -x 2
x = -l
“ xy
‘
2
'
3
(1 - x2 - y2) dydx
y "
1“ x2
) * dx
°
1
~
8 3 “ 3 ' 8 77 =
77.
Sumando los resultados anteriores para los tres términos de (4.31) se obtiene
ÍJ
f'dS = | —7
7 + i
77 1
+
(—
77 +
—
771+77
=
277.
Si S se representa por z = z(x, y ) mostrar que el flujo de f a través de
PROBLEMA 4.18 S es
JJt.rfS./{ ,.„ < * s
s s e c y dxdy
" ^xy JJ-
j j '1 ; ,
donde se c y
n •k
dxdy,
(4.35)
+ 1
«5y
Sol uc i ón: Si S se representa por z = z(x, y ) , entonces por el prob. 3.24, el vector normal unitario es dz .
dz .
dx
dy [3.89]
dz\2
Idz
dx)
\<9i
Por consiguiente, si y es el ángulo entre n y k , n - k = e os y =
1
(4.36) d z'2
sec y
dz
+ i
dy ,
Luego, por (3.103), el elemento diferencial del area S, dS, es dS
dz Jx
+ 1 dy
dxdy
.= se c y dxdy 1 dxdy. (n • k)
(4.37)
Entonces el flujo de f a través de S es
www.FreeLibros.me
/ Análisis vectorial
92
JJf.dS ^JJf.ndS s
s
sec y dxdy ^xy dxdy. ■JD T'-í
PROBLEMA 4.19
Sí'
Si f = zi + >>j + 2zk, usar (4.35) para calcular l | f •dS, donde 5 es la
parte de la superficie del paraboloide x 2 + y 2 - 1 —z para la cual z > 0. (Cf., problema 4.17 y figura 4.5). So l u c i ó n : C om oS se representa por z = 1 - x 2 - y 2 con z > 0, el vector unitario normal de (3.89) es dz .
dz .
* + * I + 1 dx) \ dy, 2 x i + 2yj + k (4x2 + 4y2 + 1) ^ Ahora, 1
n •k
( 4x2 + 4y2 + l ) ' 2
1
f •n
2x2 + 2y 2 + 2z [(x)(2x) + (y)(2y) + ( 2 z ) l ] =
(4x2 + 4y 2 + l ) 1^
( 4x2 + 4y2 + l ) l/2
Así que, de acuerdo con (4.35),
rr
f- d S =
JJ
rr
2 x 2 + 2y 2 + 2 z
v
I --------------------- í t ( 4 x 2 + 4y + 1) 2 dxdy J J (2x2 + 4y2 + 1 ) /2
JJ [2x2 + 2y2 + 2 (1 - x2 - y 2)] dxdy R xy s>l
/ ~ / l —x 2
2L , Ls-,
dxdy
y= - /
4 |
f
(1 - x 2)^2 dx
= 4 -2 =
PROBLEMA 4.20
2n.
Calcular J J r X dS, donde S es la superficie del cubalim itado s
*
www.FreeLibros.me
Cálculo integral vectorial
93
por x = 0, x = l , y = 0 ,y = 1, z = 0, z = 1 , r e s el vector de posición y n es el vector normal unitario dirigido hacia afuera. (Cf., problema 4.16 y figura 4.4.) So luc ió n: Con referencia a la figura 4.4, tenemos sobre S¡ paraAOCB, n = - k , dS = n dxdy = - dxdy k; por consiguiente, r x n = ( x i + y j + z k ) x ( - k ) = - y i + x j. Así que, í í r x dS = f
í
Si
0
0
(-y i
+ x j )d xd y
f 1 f 1 y dydx + j r f 1 x dxdy Jo -'o ^0 -'o
=-i
1 . 2 1
1 . 2
+
J '
Sobre S 2 para AFEO, n = -j, dS = n dzdx = - d z d x j ; por consiguiente, r x n = ( x i + y j + z k ) x ( - j ) = z i - xk. Entonces,
JJ r x d S = J >i Jr i ( z i - x k ) d z d x
2
í
í
Jo
J o
z dzd x - k
í
x dxdz
í ^o
J q
2
Sobre S 3 para OEDC, n = -i, dS = n dydz = - d y d z i; así, r x n = (x i + y j + z k) x ( - i) = - zj + y k. Por lo tanto,
JJ r
x
dS =
J' J'
( - z j + y k ) dydz
S3
=- j I J z dzdy + k JQ j*/0 Jo o Jo
¡ y
dydz
*/0 Jo
1 • 1 b =~ 2 J+2 Sobre S4 para AFGB, n = i, dS = n dydz = dydz i; por lo tanto, r x n = (x i + y j + z k ) x i = z j - y k . Así que
Sobre S s para BGDC, n = j, dS = n dzdx = dzdz j ; así que, r x n = ( x i + y j + z k ) x j = - z i + xk. Por tanto,
www.FreeLibros.me
Análisis vectorial
JJrxdS=J Se
J
0
( - z i + x k) dzdx = - ^ i + ^ k. 0
Sobre S 6 para FEDG, n = k, dS = n dxdy = dxdy k; por consiguiente, r x n = ( xi + y j + z k ) x k = y i - x j . Entonces, 1 ri
1 1 ( y i - x i ) dxdy = - i -
ISf.I — - Í I
Sumando todos los resultados anteriores se obtiene r x d S = 0.
(4.38)
JJ' s PROBLEMA 4.21
r X d S donde S es la superficie esférica cerrada representada
Calcular s
por x2 + y 2 + z 2 = a 2. Solución:
Por el problema 4.14, el vector unitario exterior n normal en S es r r n = er = - - = Ir | a
Así que, por (1.58) r X n = (1 /a) r X r = 0. Por consiguiente,
j j r x n dS = 0.
J j ~r x d S = s
(4.39)
s
4.3
Integrales de volumen
Como el elemento de volumen d V es un escalar, consideremos dos integrales de volumen sobre una región R :
j f j t d V '
Í Í J ' dy-
(4A0>
V
(4.41)
R
En el sistema coordenado rectangular dV = dxdydz,
(4.42)
y (4.40) puede expresarse también como (x, y , z) dxdydz,
(4.43)
que es la integral triple de 4>(x, y , z) sobre laregión R. Haciendo <¡>(x, y, z)= 1, el volumen V de la región R es V = J J J dV-,
www.FreeLibros.me
j j j dxdydz.
(4.44)
Cálculo integral vectorial
95
Si f = / ] i + / 2j 1 7.ik, entonces(4.41) puede resolverse en sus componentes, i.e.,
I I I ,M' ‘ j f í ' ',v PROBLEMA 4.22
‘
Calcular | | | V • r dV , donde R es cualquier región con volumen V y
'I I I
r es el vector de posición. So l u c i ó n :
Por (3.129), V *r = 3. Por consiguiente, usando (4.44), d V -3
f f f v - ' j r - f f j s
J f f j V . 3V.
(4.46)
PROBLEMA 4.23 Calcular I í IV X f d V si f = y i —xj y R es cualquier región del espacio con volumen V.
Sol u c i ó n:
Por (3.138), el rotacional de f es
V x f=
i
j
k
d dx
d dy
d dz
y
-X
0
= - 2 k.
Por consiguiente, V x f dV = - 2 k
f f f
4 .4
JJJ
d V = - 2Vk.
Definiciones alternas de gradiente, divergencia y rotacional
El gradiente de una función escalar 0, escrito grad 0 ó V
ax
oy
az
[3.105]
La divergencia de f, escrito div f ó V • f, se define como
im
f ’ d S, V • f - lim h v ->o AV j(-Sf
(4.47)
donde A V es el volumen de la región R limitada por una superficie cerradaS. El volumen A V contiene siempre el punto en el cual se va a evaluar la divergencia V *f cuando A V tiende a cero. El rotacional de f, escrito rot f ó A X f, se define como ó V x f = nma„ lim 1
(4.48)
d ondeA S es la superficie limitada por una curva cerrada simple C y n max es el vector normal unitario asociado con A S tal que la or'tntación del plano de AS dé un valor máximo
www.FreeLibros.me
%
Análisis vectorial
fdr.
(4.49)
El rotacional A X f de f se define también por su componente en una dirección particular; i.e., la componente de A X f en la dirección de n es
n.(V*f)=
lim —
(4.50)
donde A S es una superficie limitada por u ña curva cerrada simple C y n es el vector normal unitario asociado con AS, esto es, AS = n AS. Como en la definición (4.47) de V *f , AS contiene el punto en el cual se va a evaluar el rotacional V X f cuandoA S tiende a cero. PROBLEMA 4.24 So l u c i ó n :
Dar la interpretación física de A • f definido en (4.47).
En (4.27), la integral de superficie
JJ f ' d S se define como el flujo de f que s
pasa por la superficie S. Entonces J J f *dS es el flujo total de salida de f a través de la s superficie cerrada S. Por tanto, (4.47) muestra que la divergencia de f en el punto P es el lím ite del flujo de salida en la red por unidad de volumen mientras S se reduce al punto P. Si el pequeño volumen que rodea a P contiene una fu en te o un sumidero de un campo vectorial f, entonces el flujo de f divergirá de P ó convergerá a P, según que V • f sea positivo o negativo. Por tanto V *f puede considerarse como una medida de la intensidad del vector fu en te o sumidero en P.
PROBLEMA 4.25
Usando (4.47), deducir la fórmula para V*f dada en (3.127).
So l u c i ó n : Si el punto (x , y , z ) en el centro de un pequeño paralepípedo rectángulo tiene aristas Ax, A y, Az, como se mustra en la figura 4.6, entonces A V = A x A y A z. Si f = f \ ¡ + f ú + h k, entonces
d S = d yd z i,
sobre S jp ara^ F G S ,
n
i,
sobre S2paraFED G ,
n = -i,
sobre S3paraBGDC,
n = j,
sobre S4para<4FF//,
n = -j,
sobre S5paraFED G ,
n = k,
sobre S6para/4/íGfí,
n = -k ,
dS = - d y d z i , dS = dzdxj,, d S = - dzdx j, d S = dxdy k, d S = - dxdy k.
Por consiguiente, ignorando las contribuciones diferenciales de orden superior, X
Figura 4 .6
S o lu c ió n al pro b le m a 4 .2 5 .
s
d l l A x\ f - d S = I/, +
dx
2
AyAz,
S,
r
.
d
/
S
di. A x \ . — j 4 y A 2,
I I S2
jy www.FreeLibros.me
d i. A y \ 2 — AzAx,
Cálculo integral vectorial
j j f '" 8 - -('> -* £
j j ' - rfS _ ( í, +
J
J
f
.
97
)
7
^
^ | í ) A*Ay,
3
y ) A ,A y .
«6 Sumando todos los resultados anteriores se obtiene Cf
/dfí
df 2
df 3\
¡d i.
df2
df3.
dx
dy
dz
s Por lo tanto, por (4.47),
f PROBLEMA 4.26 So l u c i ó n :
dí. dí 2 dt, f-dS= — + — + — . dx
dy
[3.127]
dz
Dar la interpretación física de V X f como se define en (4.48).
La integral de la derecha de (4.48) o (4.49), i.e.,
f* r es la circulación de
f alrededor de C definida en (4.9). Por consiguiente, (4.48) muestra que la magnitud de V X f en un punto P es el límite de la circulación por unidad de área mientras la curva C se reduce al punto P, i.e., la intensidad de circulación en P. En general, la circulación de f alrededor de C depende de la orientación del plano de C. La dirección de V X f es la dirección en la cual ocurre la circulación máxima. PROBLEMA 4.27 Usando (4.50) deducir la definición de rotacional de f como el producto cruz de V y f como se da en (3.138). No se tengan en cuenta los términos de orden superior. So l uc i ón : Si un rectángulo EFGH con respecto al punto (x, y , z) tiene lados A y y A z, como se muestra en la figura 4.7, entonces AS = i A S = iA>>Az. Si f = f \ i + / 2j + / s k , entonces dL A z para el lado E F , dr = dy j, Ay,
L d i 3 Ay\
para el lado FG,
dr = d z k ,
para el lado GH,
d r = - dy j,
l'- H df 2 A z \ '■ d r = - [ / , + — y j A y , / J GH
4 .2 7 .
para el lado HE,
/
dr=-dzk,
dí± Ay f • d r = - [(
H E
dy
2
Az.
Sumando todos los resultados anteriores se obtiene
i
1"
/ di ,
dL
a/,
dt,
\dy
dz
dy
dz
AS.
Por lo tanto, por (4.50),
www.FreeLibros.me
i
Análisis vectorial
98
a/,
a/.
ay
az
(4.51)
Mediante integración semejante alrededor de rectángulos en los planos zx, x y ó por simple permutación cíclica de x, y, z y 1, 2, 3,
/a/,_ J - ( V x f) =
\az
dJ i dx )
¡df2
di ,
(4.52)
(4.53)
M v * f)^ ~ a 7 Como para un vector arbitrario A , por (2.65), A = ( A * i)i + (A *j)j + (A 'k )k , obtenemos (3.138): di 3 di, ( di , V x f =(— - — + dy dz 1 + U z i
i
a/3\ ax
j +
a/, ax
ay
k
A A A dx
dy
dz
h
i*
f.
(4.54)
Son definiciones alternas de V<¡>, V ’ f y V X f, las siguientes: s¡4> ~
üm A — r f d s < ¿ , 4 v - o AV J j
(4.55)
5
S
V
X
f-
lim i - f f d S x f , 4 f - o At ' J j s
,4A7' (4.56)
donde A V es el volumen de la región R limitada por una superficie cerrada S.
PROBLEMA 4.28 Mostrar que (4.55) es consistente con la definición V0 dada por(3.105). No se tengan en cuenta los términos de orden superior.
So l u c i ó n :
Se sigue el procedimiento del problema 4.25 y se usan las mismas notaciones.
(Cf., figura 4.6). Así,
S,
/ / ■ « * - - i (* --5 7 y ) V A ». S2
J J d s t - j(* + ^ ) a » a « ,
Jj <js0 . - j(í - ^ ) a»a«, S4
www.FreeLibros.me
99
Cálculo integral vectorial
f j d S5
f
j d
dcf> A z \
y]AxAy.
Sumando todos los resultados anteriores se obtiene AxAyAz
AF. Por consiguiente, por (4.55), obtenemos (3.105): rr
1
dcf)
d S ch = i —
V<¿ = lim — O 4 7 ,0 AV JJ
^
dx
dcf)
dcf)
+j — +k dy
PROBLEMA 4.29 Mostrar que 4.56 es consistente con la definición de V X f dada por (3.138). No se tengan en cuenta los términos de orden superior. S ol u c i ó n : Se sigue el procedimiento del problema 4.25 y se usan las mismas notaciones. (Cf., figura 4.6). Tenemos para la superficie S ¡ , i + í2 j + /3 k) dyd z St
Si
= J J (k f2 - j f3) dyd z
d í2 A x\
k 1' • 4 s r
/
<9/3 A x \
) AyAz _ l r
t
1 AyAz-
j í t
De modo similar, para las otras cinco superficies, rr
JJJ « ,
[
d i/,2 A x \
\~
dx
/
dí,
Ax\
2
s
S2
JJcTS X r - -h (/. ^ |íi
ffds
x ' - k ('■ -
J J ,
3 -r ---- ---- 1 Az Ax, dy
2 í
t ) AzAx ■ ‘ ('• ■ 5 7 t ) AzA* d /2 A z \ Ax Ay - i (f2 + — — J AxAy,
d í, A z
II
rfs* 1"
>'■ - 5 7 T Í A l A y + ‘ l ' 1 * « 7 T 1 i , , A y -
Sumando todos los resultados anteriores se obtiene
www.FreeLibros.me
100
Análisis vectorial
.
dSxí f f
.
fdfi_ _ Wy dz (dU_ \d y
dh dz)
\dz
dx I
. d í,
A x A y Az
\dx
dy /
\dx
d y /.
d í.
\d z
AV.
Por tanto, por (4.56) obtenemos (3.138): V x f = lim a v
-»o
f f dS x f = i Jj \dy
PROBLEMA 4.30 So l u c i ó n :
dz /
+j
\dz
- 1 ^ -) + k dx /
d f2
d í!
<9x
dy
Hallar una representación integral equivalente del operador V.
Por (4.55-6), V se representa por lim — f f d S . V = lim A■ Vv-»o A F JJ
4.5
(4.57)
Teorema de divergencia o de Gauss
La definición de divergencia, [4.47]
div f = lim — f f f d S , 4v-o AV j j
da el valor de div f en un punto. El teorema de divergencia o teorema de Gauss se obtiene ampliando (4.47) a una región finita. Si f es una función vectorial continua en una región R con volumen V y Ja por una superficie cerradaS, entonces (4.58) H
S
Este teorema tiene una consecuencia especial, llamada el teorema de Green. Este teorema se llama también teorema de Green para el espacio porque en dos dimensiones es el teorema de Green para un plano. {Véase sección 4.6 y (4.112).] PROBLEMA 4.31
Verificar el teorema de divergencia (4.58).
So l uc i ón : Considérese una superficie finita cerradaS que encierra una región i? con volumen V. Se divide R en N subregiones con volúmenes A , A V 2, • • • , A VN . En el punto (x ¡, y¡, z¡) dentro de A V¡, donde / = 1, 2, • • •, N , la definición de divergencia (4.47) da V -f ~ ~ ^ y '
Jj
f 'd S + e¿, donde e¡ -*■0 cuando A
-> 0. Así que
&s¡
V - f A V¡ = J j í - d S + E¡ AV¡. AS,
La suma del volumen total V es
www.FreeLibros.me
(4.59)
101
Cálculo integral vectorial
Ahora consideremos el límite de esta expresión cuando N -> °°. La frontera de la superficie de cada AV¡ consiste en numerosos segmentos que son o parte de la frontera S o frontera de dos subregiones adyacentes. Las integrales de superficie de las superficies adjacentes con frontera común se cancelan pues las normales externas tienen direcciones opuestas sobre la superficie frontera común. Así que sólo queda la integral de superficie sobre S. Además por la definición de integral múltiple, lim
V • f AVf =
JJJ V - f dV.
Debido a que
1 ff f - d S = JffJ f - d S ,
lim Y 7V->00 Z—J J J tenemos
JJJ V - f d V = J J f - d S + R
£¡ AV¡.
lim
S
(4.60)
1 =1
Para el segundo término de la derecha de (4.60), jr
e¿A V¡
< ^
| 8; | A V¡ < | e m | '¡ T AV, = | £m | V, i =1
donde em = max e¡. Sin embargo, em -*■ 0 cuando 7V- >° ° yA^- *- 0. Por consiguiente, N
lim V N -»oo «
8; AV; —♦ 0. J
i =1
Así que tomando límites,
PROBLEMA 4.32
Dar la interpretación física del teorema de divergencia (4.58).
So l u c i ó n : Como se muestra en el problema 4.24, la divergencia de un campo vectorial f en un punto dado es la densidad del flujo de salida desde ese punto. La divergencia, o teorema de Gauss (4.58) establece que el flujo total hacia afuera desde una superficie cerradaS es igual a la integral de la divergencia a través de la región/? limitada porS. PROBLEMA 4.33
Mostrar que si r es el vector de posición, entonces
s
r- d S = 2>V,
(4.61)
s
donde V es el volumen de la región R limitada por la superficie cerrada S. So l u c i ó n :
La divergencia V • r es dx dy V • r = —— + — + —— = 3dx dy
Así que aplicando el teorema de Gauss (4.58),
www.FreeLibros.me
dz [3.129] dz
Análisis vectorial
§
' d s = í f j * - ' d r - 3 í s í d v =s v -
S
PROBLEMA 4.34
R
R
Usando 4.61, hallar j j r ‘dS para la superficie S de una esfera s
x 2 + y 2 + z 2 = a2 . (Cf. problema 4.14). Solución:
El volumen V de una esfera con radio a es V = (4/3)77 a3. Por consiguiente,
por (4.61),
s
PROBLEMA 4.35
r • <íS = 3 • — n a 3 = 4 77a 3.
3
Usando (4.61), calcular J J r ' d S , donde S es la superficie de un cubo
limitado por x = 0 ,x = 1, .>>= 0, y - 1, z = 0, z = 1, como se muestra en la figura 4.4. El vector exterior unitario normal n y el vector de posición r de los puntos situados sobre la superficie del cubo están dirigidos hacia el lado opuesto del origen. (Cf., problema 4.16). S o lu ció n :
Como el volumen del cubo es 1, entonces por (4.61),
s
r - d S = 3V = 3 -1 = 3.
PROBLEMA 4.36
Mostrar que
( 4 -6 2 )
R
S
donde S encierra la región R , r es el vector de posición y Ir I= r. Solución:
Por(3.167), V - I r " ' 1 r] = (n + 2) r " - 1.
Si n = - l , = (—1 + 2 ) r~2 = —
V-
(4. 63) #<-
Por consiguiente, aplicando el teorema de divergencia (4.58), 7 .dv- $ R
PROBLEMA 4.37
R
S
Mostrar que, para cualquier superficie cerradaS, JJ V
X
f . cf S = 0.
(4.64)
s
Solución:
Por (4.58),
JJ V x f - d S = JJJV - ( V x f ) d V . R
Pero por(3.143),V (V X f) = 0 ; p o r l o tanto,
www.FreeLibros.me
(4.65)
Cálculo integral vectorial
103
V x f • d S = 0. t i
PROBLEMA 4.38 Si f = P(x, y , z)i + Q(x, y , z )j + R (x , y , z)k y R es la región limitada por una superficie cerrada S , mostrar que el teorema de divergencia (4.58) expresado en coordenadas rectangulares es j dxdydz =
jjj" R
Solución:
JJ(P d y d z + Q d zd x + R dxdy).
(4.66)
S
Si escribimos t = P i + Q j + R k, n = e o s a i + e o s /3 j + e o s y k,
entonces, en general, para cualquier superficie S, n • i dS = eos (X dS = dydz, n • j dS = eos ¡3 dS = dzdx, n • k dS = eos y dS = dxdy.
(Cf., problema 4.15). Como la divergencia de f es v .f= ^ + |£ + ^ , dx dy dz
/Jjfv-'dV- Jlf (ff * 57+ff )dxdyi Como antes, J j f - d S = J j (Pi + Qj + Rk ) - n dS s
s (P eos OL + Q eos /3 + R e os y) dS ■
t i
J j (P d yd z + Q d zd x + R dxdy) s Por tanto, (4.58) se reduce a + —— ) d xd yd z = jj (P d yd z + Q d zd x + R dxdy). dy dz fJHff+f? s PROBLEMA 4.39
Usando el teorema de divergencia (4.66) calcular / =
JJ x d y d z + y d zd x + 2 z dxdy, s
donde S es una superficie que consiste en la superficie del paraboloide x 2 + y 2 = 1 ~ z, 0 < z < 1, y el disco x 2 + y 2 < 1, z = 0, como se muestra en la figura 4.8 (Cf., problema 4.17). So l uc i ón :
Por el teorema de divergencia expresado en coordenadas rectangulares (4.66),
www.FreeLibros.me
4 .3 9 .
104
Análisis vectorial
dx
dy dy
■ S i l dx
<9(2z) dxdydz dz
=
/// 4 dxdydz dV
f1 - z) d z
Jo =4
7T- -
= 277, que es el mismo resultado obtenido en el problema 4.17. Obsérvese que la contribución del disco x 2 + y 2 < 1, z = 0 es cero, pues sobre este disco n = - k y f * n l z=0= - 2 z l z=0 =0.
PROBLEMA 4.40
Usando el teorema de divergencia (4.66), calcular ( x 3 d yd z + x 2y d zd x + x 2z dxdy),
z
ÍJ donde S es la superficie cerrada que consiste en el cilindro x 2 + y 2 = a2, 0 < z < b, y los discos circulares x 2 + y 2 < a 2, z = 0 y x 2 + y 2 < a2, z = b , como se muestra en la figura 4.9. So l u c i ó n :
Por el teorema de divergencia (4.66), l =
x 2 dxdydz 2_
La s u p e rfic ie del p ro blem a
=5 I
4 .4 0 .
m
¡
2
x2 dxdydz
| - a J —/ a 2— y'1
•rn fbfaI b
Jo 20 6 —— 20 b 3
pa
Jo
(*•/ a 2 —y 2
x 2 dxdydz
3
ra . 2 2.v , |i " (a 2 -~ y. ) dy Jo 3 16
5 n a 4.b. — 4
Mostrar que el teorema de divergencia (4.58) puede extenderse al volumen infinito exterior a una superficie cerrada con tal que PROBLEMA 4.41
www.FreeLibros.me
105
Cálculo integral vectorial
lim
r 2 |
(4.67)
f | = O,
r-> oc
donde r es la distancia desde el origen a cualquier punto del espacio. So l u c i ó n : Sea S una superficie cerrada limitada por una esfera S r de radio r con su centro en el origen. Si f es un campo vectorial en una región R limitada por S y S r, entonces, aplicando el teorema de divergencia (4.58),
Jjj V - f dV=jj t - d S + jj R
S
(4.68)
f-dS.
Sr
(Véase figura 4.10). Por la desigualdad de Cauchy-Schwartz (1.45), I f #h l< I f I; entonces,
JJ f - d S
=
jJf-nc/S
<
jj | f - n | dS < jj | f |
dS.
Por lo tanto, si If I< e/r2 en todos los puntos de S r, e
dS = 477 8
p u c s j j d S = 4nr2 . Así, si lim r2 I f I = 0, F ig u ra 4 .1 0
Sr
La s u p e rfic ie cerrada del p ro b le m a 4 .4 1.
lim
f f
r->0°
JJ
t - d s = o,
de lo cual se obtiene
JjJv.idv.fr.di. R
S
Obsérvese que en esta integral de superficie, la normal apunta hacia adentro. En la figura 4.11, r es el vector de posición que representa un punto P sobre la superficie S. El ángulo sólido d ü subtendido en el origen O por un elemento dS de la superficies se define como (4.69)
dQ ~rdH. r3
i
El ángulo sólido total O subtendido p o rS es (4.70)
!\ ís f r;ís
Mostrar que el ángulo sólido total subtendido por una superficie cerrada S en el origen O es cero si O está fuera de la región R limitada por S, y que es 47T si O está dentro de i?. PROBLEMA 4.42
So l u c i ó n :
Por el teorema de divergencia(4.58),
J f^ = J jH ^ h S
R
Sin embargo, por (3.168), V *(r/r3) = 0 dentro de i? si r ¥= 0. Por otra parte, si O está fuera de R,
www.FreeLibros.me
Fig ura 4.11
A n g u lo sólido .
106
Análisis vectorial
entonces r ^ O , como se muestra en la figura 4.12 (a). Por tanto, cuando O está fuera de R, el ángulo total subtendido por S es
JÍ1?
n
(4.71)
= 0.
Si O está dentro de R , construim os una pequeña esfera S a de radio ra y consideramos la región R ' limitada por S y S 0, como se muestra en la figura 4.12 (b). Si aplicamos el teorema de divergencia (4.58),
.
r
íí'T
r
S+Sq
pues r
s
.rr
iij
Sq
• r
0 dentro de R '. Así que
En los puntos de S a , r = rQy la normal positiva o exterior está dirigida hacia O, o sea que n = - r/rQ; por consiguiente,
dS
# ^ ■ 1 "
so
so
-Si so
Figura 4 .1 2
dS
- tJ
(b ) S o lu c ió n al p ro blem a
-------- q n ro rTo2
4 .4 2 .
= —477; y cuando O está dentro d e R , el ángulo total subtendido por 5 es
■d S
4.6
= — ( — 4 7 7 ) = 477.
(4.72)
Teoremas de G reen
El primer teorema o identidad de Green establece que si < t>y \p son funciones escalares que tienen segundas derivadas continuas en una región R limitada por una superficie cerrada S, entonces J j J ( < £ V 3<£ +
j
j
j
R
f
^
Vt ¿ - dS.
(4.73)
S
El segundo teorema o identidad de Green establece que si < t>y \l>son funciones escalares que tienen segundas derivadas continuas en una región R limitada por una superficie cerrada S , entonces
JJJ
-
ifiV24>)dV j j
www.FreeLibros.me
(0 V 0 - «A V<^>* c/S.
(4.74)
107
Cálculo integral vectorial
Estos teoremas o identidades se obtienen del teorema de divergencia (4.58) usando unj función vectorial apropiada. funciones vectoriales que tienen segundas derivadas continuas en una región R limitada por una superficie cerrada.?, entonces
■
[f-y:g -
g-v J fl dv
• '. . ^ •
■■ ' (
j j Ir
■ «i - i r
•
• e.i
,
g • (7
.
'
'
' ' ■. '
f ) - g (v • r )]• <js.
(4.75)
donde V2 f - V ( V * f ) ~ V X ( V X í ) por Í3.164). I ste teorenu equivalente al segundo teorema de Green (4.74) que relaciona dos
Verificar el primer teorem a de Green (4.73).
PROBLEMA 4.43
So l u c i ó n :
Usando el vector identidad, V*(0f) = <¡>V"f + f - V0 , d o n d e f = Vi//, V-(0VV2
(4.76)
Integrando sobre la región R ,
JJJ V-(^V^)dV= JJJ W .V2
(4.77)
R
Aplicando el teorema de divergencia (4.58),
JjJ V. ( 0V. / r ) cf K = R
(4.78) S
Así que, J J J * (9S V V + V<¿-V0 d F = R
.
JJí¿V^-c/s. S
Por la definición de derivada direccional (3.98), i.e., 30/ds = grad 0 - T , dt¡/ - Ó V ^ - n d S = <{ rel="nofollow"> ~ - d S , dn
(4.79)
donde d\p/dn es la derivada normal. Por consiguiente, (4.73) puede escribirse como
JJJ ( < ¿ W PROBLEMA 4.44
So l u c i ó n :
+ V 0- V ¿ )d V =
Jjó ^
dS.
(4.80)
Verificar el segundo teorema de Green (4.74).
Intercam biando0 y 0 en (4.73),
JJJ (t/'V2>+ V'/r-V
Jjtf,V
Entonces, sustrayendo (4.81) de (4.73),
www.FreeLibros.me
(4.81)
108
Análisis vectorial
jjj (
jj
i//S72
Q ( 0 V«A - < A \ 7 0 ) - d s .
[4.74]
Por (4.79), dib
dd rel="nofollow">
Vi/f • dS = —— c/S,
\/cf )- dS= - r - d S ;
an
an
Por consiguiente, (4.74) se puede expresar como
J]J(<¿V2
PROBLEMA 4.45
(4-82) S
Si \¡j es armónica en una región R encerrada por S , entonces demostrar
^
- # j j v i { , - d S = j j j ^ d S = 0. s
Solución:
(4.83)
s
Si hacemos 0 = 1 en el primer teorema de Green (4.73),
jj"V'A-dS= jjj S
(4.84)
R
pues V (l) = 0. Como i¡j es armónica, entonces por la ecuación de Laplace (3.133), V2 1// = 0. Por consiguiente, (4.84) se reduce a jjvt/'-dS=
PROBLEMA 4.46
jjv<¿s-ndS=
j j ^ d S = 0.
Si (¡>y \¡¿ son armónicas en una región R encerrada por S, entonces
demostrar que
jj
VAV<£)- dS
jj (cf>
=
s
dS
= 0.
(4.85)
s
Solución: Si <¡>y \¡j son armónicas, entonces por la ecuación de Laplace (3.133) tenemos, V20 = V2 1// = 0. Así que por el segundo teorema de Green (4.74),
jjo w -
PROBLEMA 4.47
Solución: fX(VXg),
R
Deducir el tercer teorema de Green (4.75).
Aplicando el teorema de divergencia (4.58) a los vectores f(V *g) y
jjj
V-[f(V-g)]dK=
R
jjj V • tf
J ^ " f ( V- g ) - c / S ,
(4.86)
S
X
(V x g) ] d F = j j j f x ( v
R
X
g)-c/S.
(4.87)
S
Por las identidades vectoriales (3.155) y (3.157), V-[f ( V- g) ] = ( V - g ) ( V - f ) + f - V( V- g) , V • [f x (V x g)] = (V x g) • (V x f ) - f • V x (V x g).
www.FreeLibros.me
(4.88) (4.89)
Cálculo integral vectorial
109
Por consiguiente, J 'J 'J 'V • [f(V • g) ] dV = J J J í ( V - g ) ( V - f ) + f - V ( V - e ) ] d V , R
JJJ
(4.90)
R
JJJ
V - [ f x ( V x g) ] dV =
R
[ ( V x g M V x f ) - f - V x ( V x g)]c/K.
(4.91)
R
Intercambiando los papeles de f y g y restando, J J J [ r - V ( V - g ) - g - V ( V ‘ f)]dV=
Jf[f(V-g)-g(V-f)]-dS,
R
(4.92)
S
J JJ [f • V x (V x g) -
g • V x (V x f)] dV
R
= - jf j^ tf x (V X g) - g X (V X f ) ] - d s . (4.93) s
Sumando (4.92) y (4.93) y usando V 2f
=
V(V-f) -
V
X ( v
X
f),
[3.164]
obtenemos
JJJ
( f - V 2g - g - V 2f ) d V =
R
f f [f
X
(V
4.7
X
g) + f ( v • g) - g
X (V
X
f)
-
g( V • f )] • d S.
[4.75]
Transformaciones de integrales de volumen a integrales de superficie
El teorema de divergencia (4.58) representa una transformación de integral de volumen a integral de superficie en donde interviene la divergencia de un vector; i.e.,
JJJ
V -f dV *
R
e/S-f.
[4.58]
S
Extendiendo las definiciones de gradiente y rotacional de un vector a volúmenes finitos, obtenemos los siguientes teoremas: El teorema de gradiente expresa que si <¡ rel="nofollow">es una función escalar continua en una región R limitada por una superficie cerradaS, entonces
JJJ
VádV ^
R
JJ
d S 4>.
(4.94)
S
El teorema de rotacional expresa que si f es una función vectorial continua en una región R limitada por una superficie cerradaS, entonces
JJJ V x f dV = R
ll'dS xf. S
Obsérvese que los teoremas (4.58) y (4.94-5) se pueden expresar como
www.FreeLibros.me
(4.95)
110
Análisis vectorial
fv —
(4.96)
R
S
donde a es un escalar o cantidad vectorial y la estrella (*) representa una forma aceptable de multiplicación, i.e., producto punto o cruz, o producto simple. PROBLEMA 4.48
Demostrar el teorema del gradiente (4.94).
S o lu c ió n : Sea f = 0a, donde a es cualquier vector constante. Aplicando el teorema de divergencia (4.58),
JJJ V ‘ (
JJcf>a-dS.
R
(4.97)
S
Pero por (3.155), V • (cf>a) = cf>V • a + a • V<£. Como a es una constante, V *a = 0; por consiguiente, V • (cf, a) = a • V<£. Asi que (4.97) puede escribirse como a-
íJJJ V cf >dV- JJ cf, d s ) = 0. \
R
S
(4.98)
/
Como a es cualquier vector constante, la expresión del paréntesis se anula y queda demostrado (4.94). PROBLEMA 4.49
S o lu c ió n :
Demostrar el teorema de rotacional (4.95).
Si a es cualquier vector constante y sustituimos f en (4.58) por f X a, entonces
JJJ V - ( f x n ) d V = JJ ( í x a)-c/S . R
(4.99)
S
Pero por (3.157), V - ( f x a) = a -(V x f) - f - ( V x a). Como a es constante, V X a = 0; por consiguiente, V • (f x a) = a - (V x f). Por la regla de permutación del triple producto escalar (1.77) o (2.53), f x a- d S = a- d S x f. Así que (4.99) se puede expresar como I
V x F dV -
j j d S x f j = 0.
(4.100)
Otra vez, puesto que a es cualquier vector constante, la expresión del paréntesis se anula y queda demostrado (4.95).
PROBLEMA 4.50
Mostrar que para una superficie cerrada^,
íf d S = 0. www.FreeLibros.me
(4.101)
111
Cálculo integral vectorial
S o lu ció n :
Por el teorema del gradiente (4.94),
J\MS. JJJWdV. S
R
Si 0 = 1, entonces V> = 0; por consiguiente, J J d S = 0.
s Usando la representación integral de V (4.57) i.e.,
PROBLEMA 4.51
V=
lim
AV
0
O ds,
AV
JJ
mostrar que
s
V - ( ^ f ) = <¿>V-f + f-Vs¿>.
[3.155]
So l u c i ó n : Sean 0 O y fo los valores de 0 y f en algún punto F q , y sea A F u ñ a región pequeña que rodea a P0. Sobre la superficie S que limita a A V los valores de 0 y f son 0 = 0 o + A0, f = f0 + Af. Usando (4.57), V- W>f ) =
lim av
=
-> o A k
ffdS.& t)
Jj
s
lim
s
, Í ” o A f [ 0# § L s
d S ' t+
lim^L J J
dS • f + f f d 8 - ( ¿ - ¿ 0) ( t 0 + At)
d S- f
lim ^ \d >o ^ '
+
f0 •
s
jj d S
s
(f> - cf>0 í 0 •
JJ d S
s
s
j^dS-A^Af^.
(4.102)
Pero por el resultado del problema 4.50,
ff Como la última integral
d S = 0.
[4.101]
j j d s ' A 0A f es un término de orden superior, puede descartarse
s en el lím ite. Entonces, en el punto P0 ,
v - w - j f r . h í ' 1* § = 0 O lim av-oAF
d s - < * 1-- f f d s * ]
Q d S • f + f • / lim Q dS JJ L v ^ o A F JJ s \ s
<¿V-f + f - V 0 ,
www.FreeLibros.me
1 t
112
Análisis vectorial
donde todas las expresiones están calculadas en el p u n to /V C o m o />0 es arbitrario, queda demostrado (3.155). PROBLEMA 4.52
Si r es un vector de posición, mostrar que para una superficie cerrada
s
S,
Solución:
r x d S = 0.
(4.103)
Por el teorema del rotacional (4.95), d S x r = I I I V x r dV. S
R
Pero por (3.141), V X r = 0; por consiguiente, f l rr xx Jd SS =. — - pO , d S x r = 0. s
4.8
s
Teorema de Stokes
t i teorema de Stokes establece que si S es una superficie limitada por una curva cerrada simple C y f es una función vectorial que tiene primeras derivadas parciales continuas sobre S y C, entonces J j V x f . c / S = <j) f • dr.
PROBLEMA 4.53
(4.104)
c
s
Demostrar el teorema de Stokes (4.104).
Solución: Considérese una superficie S limitada por una curva cerrada simple C. Se divide 5 en N subregiones tan pequeñas que puedan considerarse planas con áreas A S ,, A S 2 , ' ' ’ , ASyy. En los puntos (x¡, y¡, z¡) de AS¿, la definición (4.50) del rotacional da n • V x f AS, =
(4.105)
donde e .^ 0 , cuando A S¡ -> 0, y n es el vector unitario normal asociado con AS,-. (Véase figura 4.13.) La suma sobre la superficie total S da N
N
-
N
£ n . V x f A S i = £ < h f - d r + 2 ] e¡ AS,, i =1 ci
i =1
Figura 4 .1 3
D e m o stra ció n del teorem a de S tokes.
i=1
Ahora consideremos el límite de esta expresión cuando N ^ - ° ° . La frontera C¡ de cada ASj consiste en pedazos que son o parte de la frontera C o parte de las fronteras de las dos subregiones adyacentes. Las integrales de línea a lo largo de curvas fronteras adyacentes se cancelan, pues los vectores dr tienen direcciones opuestas; así que queda sólo la integral de línea a lo largo de C. Por consiguiente, n • V x f AS, = JJ~ n ■V x feIS - JJ~ V x f • d S,
lim
1 =1
s
lim
y
(f) f - d r = { í )
N ->°c l— l J
www.FreeLibros.me
s
J
í-dr.
113
Cálculo integral vectorial
Así que cuando J J v x f - d S = <j) f - d r + lim s
c
£ ¡A S ,.
(4.106)
>'- 1
Para el término restante,, N
£
£/ AS;
< ^
i= 1
| e , | AS; < I e m I J
AS, = | e m | S,
i =1
donde em = max e¿. Pero em -►0 cuando N
°°,\JSj-* 0. Entonces,
lim V 6 i ASf —♦ 0. 7V->oo L—i i = 1
Así, en el lím ite, j j v
X
f . d s = < £ f.tfr.
s PROBLEMA 4.54
c
Dar la interpretación física del teorema de Stokes.
El rotacional de un campo vectorial f es la intensidad de circulación en un punto para f (véase problema 4.26). El teorema de Stokes (4.104) establece que la circulación total alrededor de una curva cerrada C es igual al flujo del rotacional f a través de una superficie S encerrada por C.
Solución:
PROBLEMA 4.55
Solución:
Mostrar que si r es el vector de posición, entonces
>(j) r-d t = 0.
(4.107)
Por (3.141), V X r = 0. Por consiguiente, por el teorema de Stokes (4.104), r-dr = || c
PROBLEMA 4.56
V x r • c/S = 0.
s
Usando el teorema de Stokes (4.104), mostrar que, para cualquier
superficie cerrada S,
§
V x t - d S = 0.
Sean S¡ y S 2 dos regiones en las cuales una curva cerrada C divide a una superficie cerrada S, como se muestra en la figura 4.14. Aplicando el teorema de Stokes (4.104) a S , y a 5 2,
Solución:
JJ V
x
f - d S = < j ) t-di,
s,
c
J J V x f - d S - (j> f - d r = -< ^) í - d r . s2
c
c
Por consiguiente, para la superficie cerrada S,
www.FreeLibros.me
4 .5 6 .
Análisis vectorial
114
V x f •dS = j j V x f •dS + j j V x f •dS s, s,
s
s.
+ (j> (T) f • d r =<£ f -•dd r + Jc Je = 0. So l u c i ó n al t er na: Considérese una superficie S casi cerrada con una pequeña abertura limitada por una curva cerrada simple C como se muestra, en la figura 4.15. Aplicando el teorema de Stokes (4.104), j j
V x t - d S =<j) t - d r .
s
Figura 4 .1 5
S o lu c ió n a lte rn a al
c
Ahora redúzcase más y más la pequeña abertura de modo que en el límite sea un punto. Entonces la superficie se vuelve cerrada y la integral de línea tiende hacia cero. Por consiguiente,
p ro b le m a 4 .5 6 .
V
X
f - d S = 0.
s
s
PROBLEMA 4.57
Si C es una curva cerrada, mostrar que / { j)
So l u c i ó n :
- d r = 0.
(4.108)
Por el teorema de Stokes, <j) V < f i - d r = j j V x ( V 0 ) - d S , c
[4.104]
s
donde S es la superficie encerrada porC . Pero por (3.142), V X V>=0;por consiguiente, queda demostrado (4.108). PROBLEMA 4.58
i
Demostrar que V X f = 0 es una condición necesaria y suficiente para
i ' d t = 0 alrededor de cualquier curva cerrada C.
So l u c i ó n :
Para la suficiencia, V X f = 0; entonces por el teorema de Stokes (4.104),
(V x f ) - d S = 0.
Para la necesidad, supóngase que que^j) y f 'dr = 0 para cualquier curva cerrada C y que V X f ¥= 0 en algún p u n to P. Entonces si V X f es continua, hay alguna región que rodea a/* donde V X f ¥= 0. Se escoge una pequeña superficie plana S en esta región y un vector unitario normal a S paralelo a V X f, esto es, V X f = an, donde a > 0. Si C es la frontera de S entonces por el teorema de Stokes, (4.104), {j) f- d r = j j e
V x f - d S = j j a n - n d S = a j j dS = aS > 0,
s
s
lo que contradice la hipótesis de que<í)f mdr = 0.
www.FreeLibros.me
s
Cálculo integral vectorial
115
El teorema de Stokes para coordenadas rectangulares establece que si f = Pi + £>j + R k, entonces j> P *
♦ 0
- j j ( íf -
dydz
***
ss
c
(dQ dP\ ( T x ~ ~ J dxdy.
PROBLEMA 4.59
(4 .1 0 9 )
Verificar (4.109).
Sol uc i ón: Si f = Pi + Qj + /? k , entonces f ‘dt = P dx + Q dy + R dz
V
X
f:
i
j
A
A
A
dx
dy
dz
P
Q
R
dR \dy
y
k
dQ\ . tdP 1i + 1 dz
dR
i +
dQ
dP
dx
dy
k.
Ahora, como en el problema 4.38, J J v x ' . d S . J J V x f • n dS s
s
dR
dQ\ . fdP i + dz,
-//[(I W
_
íí—
'dP
eos a +
dR
dP ... dy
• n dS
dQ
dP
dx
dy
eos /3
dz ,
s
■jr s
do;
d R \ . ( dQ a ox ,J +
’d R \dy
dz )
dydz
| P - - d- « \ d Zdx dz dx J
dQ
dP' dy
eos y
dS
dxdy.
Por consiguiente, el teorema de Stokes (4.104) se reduce a Pdx+Q dy +Rdz
i
re ( d R _ do;
JJ
w
dP dR d yd z + —— -
dzdx
dz!
[4.109]
PROBLEMA 4.60
jy
Si f = 4yi + x') + 2zk, calcular / = | | V X f'c/S sobre el hemisferio
x 2 + y 2 + z 2 = a2, z > 0. Sol u c i ó n:
C o m o /e stá en forma de integral de superficie en el teorema de Stokes (4.104),
/ =<j) f • d r
4 y d x + x d y + 2 z dz, Figura 4.16
donde C es el círcu lo x 2 + y 2 = a2, z > 0, dirigido como se muestra en la figura 4.16. La representación paramétrica de C (problema 3.13) es x = a eos t , y = a sen t, z = 0, donde 0 < t < 2 tt. Entonces dx = - a sen t dt, d y = a eos t dt, dz = 0. Por consiguiente,
www.FreeLibros.me
El hemisferio del problema 4.60.
Análisis vectorial
116
r2n 1=1 4 (a sen t) ( - a sen t dt) + (a e o s t) (a e o s t dt) •'o = a2 f •'O
f
a2
( - 4 sen* t + e o s2 t) dt
^
( 1 —5 sen2 t) dt
'o (1 - eos 2 í)
f H = a2
3 5 . — h— eo s 2 í 2 2
f(-
,
dt
.
c/í
= - 3 a 2 ¡7.
El teorema de (¡reen para un plano establece que si P, Q.dP/by, y dQ/dx son continu i una región R del plano x y limitado por una curva cerrada C, entonces (f.110)
<j) P d x + Q d y C
R
La ecuación (4.110) se puede expresar también en forma vectorial; i.e., si f = Pi + entonces V x f h dxdy.
U> f • d r - J j C
PROBLEMA 4.61
(4.111)
R
Verificar el teorema de Green (4.110) para un plano.
Solución: Si f = Pi + Q\ y dS = k dxdy, entonces (4.110) se sigue directamente del teorema de Stokes (4.104). Alternativamente, (4.110) es un caso especial de (4.109) con R = 0. Ahora, si f =Pi + Qj, entonces i
j
k
V x f= A
A
A
dx
dy
dz•
P
Q
0
dQ . dP . fdQ d z 1+ d z , + \<9x
dP \, dy ’
f • d x = P dx + Q dy, V x f - c f S = V x f - k dxdy =
\ox
— ) dxdy. ay/
Por consiguiente, por el teorema de Stokes (4.104), < $P
R
PROBLEMA 4.62 Mostrar que el teorema de Green (4.110) para un plano puede expresarse también como V • f dxdy,
www.FreeLibros.me
(4.112)
Cálculo integral vectorial
117
donde f = Q\ - P j , n es el vector normal unitario hacia afuera de C, como se muestra en la figura 4.17 y s es la longitud del arco. Por (3.107), el vector unitario T tangente a C es
Sol uc i ón:
T = — = — i + — j. ds ds ds Ahora, por la figura 4.17, _ . I dx . dy \ , dy . dx . n = T x k = ( — i + — j x k = — i ------- j. \ds ds / ds ds De acuerdo con esto,
\ds
ds /
ds
Figura 4 .1 7
ds
Así que, {j) f • n d s = ,(j) c
q
\
Él +P ds = {j) P d x + Q dy. ds dsj Jc
Luego V -f =
dQ dx
d{-P) dy
dQ dx
dP dy
Por consiguiente, el teorema de Green (4.110) para un plano puede expresarse como ^ ( • n i / s = JJ" V • f dxdy. c
R
PROBLEMA 4.63 Mostrar que el área A de una región R del plano x y limitada por una curva cerrada simple C es A = i f ) x dy = ’u ) —y d x = — '(h x dy — y Jc “t ^ Jc
dx. (4.113)
Si P(x, y ) = 0 y Q(x, y ) = x en (4.110), entonces
So l uc i ón :
dxdy = A.
Si hacemos P = - y y Q = 0 en (4.110), entonces - y dx - j j
dxdy~= A.
x dy - y dx. Si hacemos P = ~y,
Sumando los resultados anteriores, obtenemos ,4 = ^ Q = x en (4.110),
f
( - y ) d x + x dy = '(D x dy - y dx =
dxdy R
■2
Teorem a de Green en un plano.
JJ d’“'y R
= 2 A.
www.FreeLibros.me
y
118
Análisis vectorial
Por consiguiente, A = —■(t)x dy - y dx.
PROBLEMA 4.64 Determinar el área A de una región R limitada por una elipse C cuyos ejes mayor y menor son 2a y 2b, respectivamente. Solución:
La ecuación de una elipse C con eje mayor 2a y eje menor 2b es i Y2
y 2
— a2
b2
+= 1,
z = 0.
Por consiguiente, C puede representarse paramétricamente por x = a eos t,
y = b sen t,
z = 0,
0 < t ^ 2n.
Entonces dx = - a sen t d t, d y = b eos t d t. Por consiguiente, por (4.112), dy — y dx
A =
277
ü — í 2 Jo
(a e o s í)(6 eo s t dt) - (6 sen t ) (~a sen t d t )
a b ( e o s 2 t + sin 2 í)
'
c /f
•2 7 7
i a 6 f 2 dt 2 J„ ■'o = 7ra6. 4.9
Transformaciones de integrales de superficie a integrales de línea
t i teorema de Stokes (4.104) representa una transformación de integral de superficie a integral de línea en la cual interviene el rotacional de un vector; i.e., |j V x f • t/ S = ( |) f • d t .
s
c
Como
jjv
xf-ds=
s
jj
<ís - y x
f
= J J (d S x
7 ).r ,
s
El teorema de Stokes puede expresarse como JJ(dSxV)-f=<j)dr-!.
(4.114)
s Extendiendo las definiciones del rotacional, el gradiente y el rotacional del rotacional a superficies finitas, obtenemos los siguientes teoremas: Si S es una superficie finita limitada por una curva cerrada simple C y 0 es una función escalar con derivadas continuas, entonces j J d S * VÓ - <j )¿ dt. s
www.FreeLibros.me
c
(4.115)
Cálculo integral vectorial
119
Si S es una superficie finita limitada por una curva cerrada C y f es una función vectorial con derivadas continuas, entonces
J J - *V) x f = ¿
f.
s
(4.110)
c
Obsérvese que ios teoremas (4.104) y (4.115-6) pueden expresarse como j'J '( cfs x V ) * u = ( j) dr * a ,
(4.117)
donde a es cualquier cantidad escalar o vectorial y la estrella <*) representa cualquier forma aceptable de multiplicación, i.e., punto, cruz o producto simple. PROBLEMA 4.65
Verificar (4.115).
Solución: Sea f = 0a, donde a es cualquier vector constante. Aplicando el teorema de Stokes (4.104), J ' J ' V x (4> a) • d S = ^ a ' í f r . s
(4.118)
c
Por (3.156), V x ( 0 a) = < ^ > V x a - a x V<£Como a es un vector constante, V X a = 0; por consiguiente, V x ( 0 a) = SJ(f> x a. Por la regla de permutación del triple producto escalar (1.77), x a)-dS= JJa .d S x V < £ . s
s
Por consiguiente, (4.118) se reduce a X V 0 =
-(j)0 f/r,
a-
s
c
ó a- í J j dS x
V<£-
\ s
=
0.
(4.119)
/
Como a es cualquier vector constante, la expresión del paréntesis se anula y queda demostrado (4.115). PROBLEMA 4.66 Solución:
Verificar (4.116).
Si a es cualquier vector constante y reemplazamos f en (4.104) por f X a,
JJvx(fxa)-dS= s
<j)(f x a)-c /r.
(4.120)
c
Por (3.158), V
x (f x a) = f ( V - a ) - a ( V - f ) + ( a - V ) f - ( f - V ) a .
Como a es un vector constante, V- a = 0 y ( f - V) a = 0; por consiguiente,
www.FreeLibros.me
Análisis vectorial
120
V x (f x a) = (a • V) f - a(V • f).
(4.121)
Así que J J v x ffx a W S .
/ / [(a« V ) f • c/S - a - d S ( V - f ) ]
s
s = a*
JJ [Vf ( f - c / S ) - d S ( V - f ) ] , s
donde Vf indica que V opera solamente sobre f. Aplicando la regla del término medio para el triple producto vectorial (1.98) a (dS X V) X f, ( d S x V) x f = V f(f • dS) - d S(V • f),
(4.122)
y podemos escribir J J V x ( f x a ) - d S = a-
J J ( c / S x V) x f. s
Por la regla de permutación del triple producto escalar (1.77), <£ (f x a) • d r = { p a-(c/ r x f ) = a - < p d r x f. •'c c c Por consiguiente, (4.120) se reduce a drxí,
a- J J 'c d S x V) x f = a-
o sea
JJ(d S x V l x f -
^c/rxfj=0.
(4.123)
Como a es cualquier vector constante, la expresión del paréntesis se anula y queda demostrado (4.116). PROBLEMA 4.67
Demostrar que si r es un vector de posición, entonces
í Solución:
d r = 0.
(4.124)
Por (4.115),
j4.,r, JJd S c
x
V 0.
s
Si 0 = 1, entonces V0 = 0; por consiguiente,
i PROBLEMA 4.68
dr = 0.
Mostrar que, integrando alrededor de una curva cerrada C del plano
xy,
rx rfr
= 2 A,
donde r es el vector de posición y A es el área encerrada por la curva C.
www.FreeLibros.me
(4.125)
121
Cálculo integral vectorial
So l uc i ón :
Si f = r y dS = k dx dy en (4.116), f " " - J J c
(k x V) x r dxdy.
S
Como V t = 3 por (3.129), usando (4.122), (k x V) x r = V(k • r) - k(V • r) = V(z) - 3k = k - 3k = - 2 k . Por consiguiente, j ) r x dr = - ^ d r x r = 2k JJ* dxdy = 2Ak,
(4.126)
c
y
\í
r x dr = 2A.
Si S es una superficie abierta limitada por una curva cerrada simple C,
PROBLEMA 4.69 mostrar que
<^)r x d r = 2 j j dS.
(4.127)
c Solución: entonces
Si en el teorema de Stokes (4.104) f = r X a, donde a es un vector constante,
<^) ( r x a ) - d r = j j V x ( r x a)-<JS.
(4.128)
Como a es un vector constante, tenemos por (4.121), V x (r x a) = (a - V )r - a(V • r). Además por (3.141) y (3.121), como (a - V )r = a, V • r = 3, tenemos V x (r x a) = a - 3 a = - 2 a . Por la regla de permutación del triple producto escalar (1.77), (r X a) • dr = - a • (r X dr). Por consiguiente, (4.128) se puede expresar como a- ^{ pr rx xddrr = 2a- J J " dS. c Como a es cualquier vector constante, ^ r x dr = 2 c
Jj d S. s
Obsérvese que (4.18), i.e.. tp
r x dr = j
k a 'rfr
^
2 ^ a Jk,
y (4.125) son los casos especiales de (4.127). Son importantes para (4.127) porque el doble de la integral de dS sobre la superficie S es igual a la integral de línea de
www.FreeLibros.me
122
Análisis vectorial
r X d r alrededor de la frontera de S, que es independiente del punto escogido como origen del vector de posición r y de la superficie S limitada por C.
4.10 Campos irrotacionales y solenoidales
(a ) Fig ura 4.1 8a Una región sim p le m e n te conexa.
Se dice que una región R es conexa si dos puntos cualesquiera de R se pueden unir por un arco y cada punto del arco pertenece a R. Una región R se dice simplemente conexa si toda curva cerrada de R se puede deformar continuamente hasta un punto de R. La región R de la figura 4.18 (a) es simplemente conexa. Sin embargo, la región R de la figura 4.18 (b) nó es simplemente conexa porque la curva cerrada C que rodea uno de los '‘huecos” no puede ser deformada continuamente hasta un punto sin salirse de R. Una función escalar de potencial 0 es una función uniforme para la cual un campo vectorial continuo f de una región simplemente conexa satisface f -V 0. PROBLEMA 4.70
(4.129)
Mostrar que la condición necesaria y suficiente para que la integral de
r
línea I f • c/r sea independiente de la trayectoria de integración desde el punto P hasta el P punto Q es que el campo vectorial continuo f satisfaga (4.129).
Sol u c i ó n: por (3.120),
Para demostrar la suficiencia, supóngase que f ' dr = S/
y po, consiguiente,
re ( ^
re ^
p
p
^
^ (p >
(4 J3 W
Figura 4 .1 8 b U na región que no es s im p le m e n te conexa.
y
Si 4>es uniforme, el lado derecho de(4.130) tiene un valor definido que depende solamente de los puntos extremos P y Q y no de la trayectoria. rQ Para demostrar la necesidad, supóngase que í ' d r es independiente de la trayectoria de integración. Sea •'p
fQf - c/r,
JP
donde P es un punto fijo y Q es un punto variable d e R . Como la integral de línea es independiente de la trayectoria, Q se mueve a lo largo de una curva que pasa por P sobre la cual es continuo el vector unitario tangente T. (Véase figura 4.19). A lo largo de esta curva, rQ rQ rQ f-dr f - d s = f-T ds 0 JJP P J^ P d“s“ Jj p P es una función de la longitud de arco s. Así que, dd> ds
Figura 4.19
Solución al problema 4.70.
f.
dt
f-T.
(4.131)
ds
La curva PQ podría escogerse de m odo que tenga una dirección dada en Q. Por consiguiente, (4.131) muestra que
Ahora, comparando (4.131) y (3.106),
www.FreeLibros.me
[3.106]
Cálculo integral vectorial
123
0.
(f - V«/>) • T
(4.132)
Como (4.132) se cumple para cualquier dirección de T, f - V 0 = O,
o
[4.1291
f =V0.
PROBLEMA 4.71 Si f = V0 en todas partes de una región simplemente conexa R y C es cualquier curva cerrada de R , mostrar que f •
i Sol uc i ón:
(4.133)
Si f = V $ en todas partes de R , entonces por el problema 4.70, la integral de
línea I f -dr es independiente de la trayectoria de integración. Si la trayectoria de integración *'c es cerrada, entonces P= Q en (4.130). Así que, f-c /r = < f) V>-c/r= f
d
Este problema se resolvió en el problema 4.57 aplicando el teorema de Stokes (4.104). Obsérvese que (4.133) es válido solamente cuando <¡>es una función uniforme y R es simplemente conexo. Si R no es simplemente conexo, entonces es posible tener f = V0 y la circulación de f alrededor de una trayectoria cerrada C en una región R puede no ser cero;i.e., podemos tener < j) f -dx # 0. El problema 4.72 ilustra este punto.
PROBLEMA 4.72
Sea f = -
y
x
.
_|_^2 i + ~2 '^.^2 j- (a) Mostrar que f se puede expresar
como V
f • dr
si C es un círculo de radio a sobre el plano x y cuando su centro está en el origen y cuando está en (a, (3, 0) con a 2 + ¡52 > a2 . Sol uc i ón:
(a) Como
V<¿> =
d —
dx
tan
1(f y- \) w.
v2 \
l +[ ^
d i +
dy
tan
1(í y- S| W
2
U
-y * = x1 — + y i > + x_¡— + y¡ i = f. (b) Cuando el centro de C está en el origen, / =1 =<j)
V 0 - d r = { j ) dcf> ={j) d jten 1 ^ j
La representación paramétrica de C es x = a eo s t,
y = a sen t,
z = 0.
www.FreeLibros.me
n
Análisis vectorial
124
Así que, d t = 2 tt ¿ 0 .
El resultado no es cero porque la región definida para la cual f = V0 no es una región simplemente conexa, pues 0 = tan-1 i y/ x) no está definido en el origen y la región está perforada en el origen. Cuando el centro de C está en (a, (3, 0) su representación paramétrica es x = Ot + a eos t,
y =
+ a sen t,
/3
z = 0.
Por consiguiente, p 2rr
í //S + a sen t
ta r f
ta n -1
2 77
Ot + a eos t
'
= 0.
PROBLEMA 4.73 Si R es una región simplemente conexa, mostrar que la condición necesaria y suficiente para que V X f = 0 es que f = V0 cuando f tiene derivadas continuas. Solución:
Si f = V>, V x f = V x V<¿> = 0.
[3.142]
Por otra parte, si V X f = 0, entonces por la solución del problema 4.58, o por el teorema de Stokes (4.104), f • cfr = 0
para cada curva cerrada simple C de R . Por la solución del problema 4.71, esto indica que la integral de línea de f es independiente de la trayectoria de integración y que f puede expresarse como f = V 0. El problema 4.73 proporciona una manera simple de verificar si un campo vectorial f es el gradiente de un campo escalar 0. Esto se ilustra en el siguiente problema.
J
PROBLEMA 4.74 Dos campos vectoriales están dados por f = 3y2i + z ¡ + 2 y k y g=y z¡ + x z j + x y k . Determinar si estos campos vectoriales son los gradientes de campos escalares. So l uc i ón : Como f = 3 y 2i + z j + 2 y k, i _d_ dx
j
k
d_ d dy dz
3y 2
z
2y
i - 6y j ¿0. Por consiguiente, f no es el gradiente de un campo escalar. Como g = yz¡ + xzj + x y k , i
j
k
d dx
d dy
d dz
yz
xz
xy
= (x - x) i + (y - y ) j + (z - z ) k = 0.
www.FreeLibros.me
125
Cálculo integral vectorial
Por consiguiente, g es el gradiente de un campo escalar. Obsérvese que g = V (xyz + constante). PROBLEMA 4.75
Si en una región simplemente conexa R, V X f = V X g, demostrar que f = g + V<£.
Solución: Como V x f = V x g,
V x ( f - g) = 0. Entonces, por el resultado del problema 4.73, f - g = V<£, y por consiguiente, f = g + V<¿.
(4.134)
npo vectorial f se llama irrotacional en una región R si V xí =0
(4.135)
en todas partes de R. Por losresultados de los problemas 4.70-2,concluimos que para que un campo vectorialf sea irrotacional en una región simplemente conexa,cualquiera de las tres condiciones siguientes es necesaria y suficiente: ' X f = 0; (2) f es el gradiente de un campo escalar; i.e., f = V<¿>; f - d r = 0.
(3) para toda curva cerrada C de Ñ,
PROBLEMA 4.76 entonces
Demostrar que si existe un escalar A ¥= 0 tal que Af sea irrotacional, f-Vxf=0.
Sol uci ón:
(4.136)
Si Af es irrotacional, existe una función escalar 0 tal que Af - V 0 .
Tomando el rotacional en ambos lados, V x (Af) - V x (V
sea AV x f = - V A x f
f x VA.
Si efectuamos el producto punto en ambos lados con f, obtenemos (por el problema 2.21). Af • V x f = f- f x VA = 0. Como A no es cero, f-V x f = 0. Obsérvese que el recíproco es cierto también.
Un campo vectorial f se dice solenoidal si en toda la región R, V - f = 0.
(4.137)
www.FreeLibros.me
Análisis vectorial
126
PROBLEMA 4.77 solenoidal. Solución:
Mostrar que un vector f que sea el rotacional de oiro vector A, es
Si f = V X A, entonces por (3.143), V • f = V - (V x A) = 0;
por consiguiente, f es solenoidal. Una función vectorial de potencial A en una región especifica R * es un campo vectorial para el cual un campo vectorial solenoidal f satisface f = V x A.
(4.138)
No hay un A único para el cual valga (4.138) porque Vx( A~V
f-t/S = 0. s
PROBLEMA 4.78 Si f es un campo vectorial solenoidal, mostrar que existe una función vectorial de potencial A tal que f = V x A.
[4.138]
So l uc i ón : La existencia de A se demuestra calculándolo. Sea f = / i ¡ + / 2j + / 3k y A = ^ 4 1i + ^ 4 2j + ^ 3 k . Ahora necesitamos mostrar que existen funciones escalares A ¡ , A 2 y A 3 tales que f = V x A, o
sea 3A}
dA2
dz
dx
dA2
dAl
dx
dy
(4.139a)
(4.139b)
(4.139c)
Para hallar cualquier A, supongamos que A t = 0; entonces, dAz
dA2
dy
dz
(4.140a)
(4.140b) <M2
www.FreeLibros.me
(4.140c)
127
Cálculo integral vectorial
Integrando se obtiene A2 = ¡
f3 d x + g2{y, z),
- fx0 /1 3 = -
f2 dx + g3(y, z),
i
donde g 2 y £3 son funciones arbitrarias de y y z pero no de jc. La diferencia de las derivadas parciales de A 3 y ,4 2 es d____________ A, c)A r*
/dfj ^ d f A
dy
\<9y
J
dz
J
dz
^
t
dgi dy
dz
Como f es solenoidal, V • f = 0, o sea df¡
( d f2
df3
dx
\dy
dz
Por consiguiente, d A 2 ^ Cx dy
dz
J
dx dx
dz
dz
= fK^ x , y, z)\ - f ^ x o, y, z)\ + —---------dÉ2 -— . dz dz Para satisfacer (4.140a), d g, dg, i 1 = /i(x, y, z) - /j(x 0, y, z) + —-------- ,rf—. dy dz Como por hipótesis, g 2
Y
(4.141)
£3 son funciones arbitrarias á e y y z, (4.141) se satisface si é , - 0,
f
g3 =
(4.142)
fi(x0, y, z) dy,
(4.143)
“Vo donde j o es una constante. Por tanto, sig 2 Y£3 están dadas por (4.142) y (4.143) podemos construir A como f3(x, y, z) c?x ■'*0 + k| T
/^Xo, y, z) dy -
f
f2(x, y, z) dx
(4.144)
ryo En esta demostración, se han hecho varias selecciones arbitrarias; obsérvese también que la región R* en la cual se construyó A se ha supuesto que es un paralelepípedo rectángulo. PROBLEMA 4.79 Si (¡) y \p son escalares con segundas derivadas parciales continuas en una región R, mostrar que f = V 0 x V i/ i
(4.145)
es solenoidal en R. Sol uc i ón:
Por (3.157) y (3.142),
www.FreeLibros.me
/
128
Análisis vectorial
V - f = V-(V<¿ x V
(r) <30 (,r) ^J rt
dr
r b q dr
q
r br ° qdr drqr bdr
I 4jrSr2 = ~ 4nE J^ » •'a
H '
y por consiguiente, ¿(«-0(a) =- M - - - J . 4 77c \ h a/ Si a =
b = r y $(“) = 0,
www.FreeLibros.me
(9.120)
Aplicaciones a la teoría electromagnética
237
(9.121)
4rr£r
La ley de Coulomb para campos electrostáticos establece que la fuerza f entre do cargas q¡ y q 2 es inversamente proporcional ai cuadrado de la distancia r i2 entre ellas: i.e., si er es el vector unitario desde q x a q 2 , entonces,
f=
e
(9.122)
4ff£r PROBLEMA 9.26
Demostrar la ley de Coulomb para campos electrostáticos.
Solución: Según la definición de fuerza de Lorentz (9.11), la fuerza f experimentada por una carga puntual q2 en el campo electrostático es f = q2E .
(9.123)
Si el campo E es producido por una carga puntual q x, entonces sustituyendo (9.119) por E en (9.123), obtenemos la ley de Coulomb para campos electrostáticos (9.122).
Un dipolo eléctrico se forma por dos cargas iguales y opuestas separadas por una distancia arbitrariamente pequeña 1 , donde el vector I comienza en la carga negativa y termina en la carga positiva. El m omento de dipolo eléctrico p es el producto de la carga q y el vector distancia 1 , i.e., P=
1.
(9.124)
Se forma un dipolo puntual cuando I — >0 y q —* - de modo que el m om ento de dipolo eléctrico p permanece constante. PROBLEMA 9.27 Si el dipolo eléctrico está localizado en el origen y el m om ento de dipolo eléctrico está en la dirección positiva de z, mostrar que el potencial y el campo eléctrico del dipolo en un vacío cuya permitividad es € , son q
¿(r) =
P •r 477Enr 3
E(r) =
3 (p • r)r - r 2p
477£„r5
-
47r £ 0r 1
02
4n-£ r
’
1
1
4u S 0 \ r j
r.
Si l < < r, entonces 1 Ti
r2
r ~ ( 1/ 2) eos 0
r + ( 1/ 2) eos
1
_ l eos 6 ______ ________ = ~ 7 2
1
- (l eos d / 2 r ) 2
^ l eos 6
Sustituyendo esto en el valor de
www.FreeLibros.me
Fig ura 9 .3
U n d ip o lo en el orig en y el p u n to de observación en r.
238
Análisis vectorial
al 0(r)=
COS
6
. c 2 • 47780r 2
(9-127)
Como la dirección del momento de dipolo eléctrico p es a lo largo del eje z, se ve en la figura 9.3 que p • r = q l k • r = q l r cos 6. Por consiguiente, según (9.127), >(r) =
•
(9.128)
v ( - ) = - rT r = - J í e r-
[3.126]
4 7t£„
También
Por consiguiente,
, ..
P-r
1
/I
Usando (3.113), el campo eléctrico E es 1 _/PT E - - V 0 = - - ----- V „3 4n- £ 0 t | t V ( p t ) - J ~ 47rE„r 47rE„
v(M • \
(9.129)
Ahora, por (3.119), V (p • r) = p,
(9.130)
y por (3.125),
Sustituyendo (9 .1 3 0 -1 ) en (9.129), E PROBLEMA 9.28
3 (p • r) r - r 2 p 47r£nr 5
Mostrar que en un campo estático producido por una corriente
estacionaria, V - J = 0. Solución:
Como
(9.132)
Las ecuaciones básicas para los campos magnetostáticos son V x H = J,
[9.105]
V • B = 0.
[9.106]
V -(V X H) = 0 según (3.143), tomando la divergencia a ambos lados de (9.105), V - ( V x H) = V - J = 0.
Obsérvese que (9.132) se puede obtener también de la ecuación de continuidad (9.3) haciendo dp/dt = 0.
La ley de circuitos de Ampere en magnetostática establece que en un campo magnetostático, si I es la com ente estacionaria total por una superficie S limitada por una
www.FreeLibros.me
239
Aplicaciones a la teoría electromagnética
curva cerrada C, entonces H dr
(9.133)
l.
i Esta ley corresponde a la ley de Gauss en electrostática. PROBLEMA 9.29
Demonstrar la ley de circuitos de Ampere en magnetostática.
Integrando (9.105) sobre una superficie S limitada por una curva cerrada C, tenemos por (9.2) para la corriente estacionaria total / a través de S,
Solución:
JJ V
x
H- d S -
s
JJ J - d S = ¡.
(9.134)
s
Aplicando el teorema de Stokes (4.104) a (9.134) se obtiene ■
í "
d r = I.
■
Obsérvese que (9.133) puede obtenerse directamente de (9.30) haciendo 9D /91 = 0.
La ley de circuitos de Ampere (9.133) puede usarse para calcular los vectores de campo magnético para casos en que exista un alto grado de simetría. PROBLEMA 9.30 Hallar el campo magnético H de un alambre recto infinito que lleve una corriente estacionaria I. S o lu ció n :
Considérese un alambre recto extendido sobre el eje z de - » a °». Como hay simetría cilindrica, se escoge un camino circular con un punto del eje z como centro con radio a , como se muestra en la figura 9.4. Por la simetría, el vector H es no solamente azimutal sino también tiene la misma dirección de dr y su magnitud es constante alrededor del contorno. Por consiguiente, por (9.133),
*
H • d r = Hcf, (277a) = I.
Así que H = Hcfr e<£ = - — e<¿. 277a PROBLEMA 9.31
(9.135)
Mostrar que en el campo magnetostático, el vector potencial A(r)
satisface
Fig ura 9 .4
V 2A = - fioJ,
(9.136)
U n alam bre in fin ito que lleva una c o rrie n te / y su cam po m a gnético H.
con la condición V • A = 0, donde n o es la permeabilidad del vacío. Solución:
Como B = V X A y B = no H, por (9.105), V x (V x A) = n 03,
(9.137)
que se reduce a V(V • A) - V 2 A =
J.
Así, si imponemos la condición de Lorentz (9.55) bajo estado estacionario, i.e., V • A = 0,
(9.138)
obtenemos
www.FreeLibros.me
Análisis vectorial
240
V 2A = - i i 0 J .
Obsérvese que (9.136) puede obtenerse también de la ecuación de onda no homogénea (9.56) con 9 2 A /9 12 = 0.
9.8
Campos con variación de tiempo armónicos o sinusoidales
Los campos producidos por cargas y corrientes cuya variación con el tiempo es simplemente armónica, o sinusoidal, se llaman campos armónicos o monocromáticos. Supongamos que la fuente varía en el tiempo como J (r, í) = J (r) cos (tü< + a ).
(9.139)
Entonces J(r, /) se representa por J(r, i) - /? e [j(r )e '< " '4(X>] - K e [ J ( r ) e '° e '" '] = R e [ J ( r ) e ^ '} ,
(9.140)
donde R e(z) = parte real de 2 y I(r) = J(r)e/ft es un vector complejo en el espacio y una función de coordenadas del espacio solamente. La fase del vector complejo I(r) es a. Como el factor de tiempo <Jui! es un multiplicador común, todas las derivadas con respecto al tiempo b /d t se pueden reemplazar por jej. Así, las ecuaciones de Maxwell ( 9 .1 6 -7 ) se reducen a
V x É + jc jB - 0 ,
(9.141)
V x H -jwD = I.
(9.142)
La ecuación de continuidad (9.3) se transforma ahora en V • J + jo>p = 0.
(9.143)
Para campos armónicos en medios lineales, isotrópicos y homogéneos, la forma armónica de las ecuaciones de Maxwell (9.141 —2) se reduce a
■
y X E + j(ú pM - 0,
(9.144)
V x H - ; w £ É = J .
(9.145)
PROBLEMA 9.32 Mostrar que las ecuaciones de Maxwell (9.144—5) en regiones sin corriente son invariantes para la transformación
E ' = ± j / j H,
H = T j / - É.
(9.146)
Si J = 0, entonces (9.144—5) se transforma en
Solución:
V x E + jeo/x H = 0,
(9.147)
V x H - ; ü j £ E = 0.
(9.148)
Según (9.146) f - T j/f a - .
h
- ±
^
e
'.
(9.149)
Sustituyendo (9.149) en (9 .1 4 7 -8 ), V x H' - ;
(9.150)
V x E ' + jcDfiW = 0,
(9.151)
www.FreeLibros.me
:
Aplicaciones a la teoría electromagnética
241
lo que expresa que los nuevos vectores E' y H ' satisfacen también las ecuaciones de Maxwell. Obsérvese que la transformación (9.146) es, en esencia, un intercambio de E y H, excepto para factores de escala.
El lema Je Lorentz establece que si E , Hay representan soluciones a las ecuaciones de Maxwell (9.147 8 ) en una región sin fuentes, partiendo de fuentes diferentes que operan a la misma frecuencia fuera de la región considerada, entonces V - ( É a x H b - É „ x H a) - 0. PROBLEMA 9.33
So l u c i ó n :
(9.152)
Demostrar el lema de Lorentz.
Usando la identidad vectorial (3.157), V • ( É a x H b) = H b • V x É a - É a • V x H b.
(9.153)
Como Ea , ETft satisfacen (9.147—8), V x Ea =
V x H b = /cuEEb.
Sustituyendo estas en (9.153), V - ( É a x Hb) = -j(o¡j.Ha • H b - ja>EEa • E b.
(9.154)
Intercambiando subíndices a y b, V • (E fa x Ha) = -;o>/iH a • H b - ;a>EEa • E b
(9.155)
pues el producto punto es conmutativo. Por consiguiente, sustrayendo (9.155) de (9.154), V • (E a x Hb - E b x H a) = 0.
El teorema de reciprocidad de Lorentz establece que si E , y Eft, a las ecuaciones de Maxwell (9.147 - 8 ) en una región sin fuentes, entonces J ^ (É a x
Hb
-
Eb
x
H a) • d S - 0.
[9.152]
son soluciones
(9.156)
Este resultado se obtiene aplicando el teorema de divergencia de Gauss al lema de Lorentz (9.152). PROBLEMA 9.34 a
Mostrar que en el campo armónico, la condición de Lorentz se reduce
V • A + jcúp.E(f) = 0,
(9.157)
y los vectores del campo complejo E(r) y B(r) se pueden expresar como B = V x A,
(9.158)
E = -;
(9.159)
Como en el campo armónico, todas las derivadas con respecto al tiempo 9/91 pueden reemplazarse por /co, (9.157) se deduce de la condición de Lorentz (9.55). De modo similar (9.158) y (9.159) se obtienen respectivamente de (9.42) y (9.43) sustituyendo
Solución:
á = ~ — - -V-Á, JCúfJLE
que se obtiene de (9.157).
www.FreeLibros.me
(9.160)
242
Análisis vectorial
PROBLEMA 9.35 Mostrar que, en una región sin fuentes,todos los vectores de campo armónico, lo mismo que lospotenciales armónicos, satisfacen las ecuaciones deHelmholtz V 2E + K * E - « ,
(9.161)
V 2 í f + K 2H
0,
(9.162)
V2A + K 2A
0,
(9.163)
+ K 2
(9.164)
donde (9.165) So l u c i ó n :
El rotacional de (9.147) es V x (V x E ) =
x H.
(9.166)
Por (9.148), V X H =/cjeE ; por consiguiente, sustituyendo por V XH en (9.166), V x (V x
É ) = -ja>n0'
(9.167)
Según (3.163), V
x (V x E ) = V ( V - E ) - V 2E = - V 2E
pues V • E = 0 en una región sin fuentes. Por tanto, (9.167) puede volver a escribirse como V2E + K 2 E = 0. De modo similar por (9.147) y (9.148), V 2H + K 2E = 0. Obsérvese que (9.163—4) puede deducirse también de (9 .5 6 -7 ) haciendo J = 0,p = 0, y reemplazando 32/9 12 por (yco) ( j u ) = ~ cj 2 .
9.9 Ondas planas El vector de onda o propagación es un vector cuya magnitud es igual a uisjue; i.e., K = K x i + K y i + K z k, donde IK f= K = (K 2 +Kj . + K * ) /2 = <ü onda. PROBLEMA 9.36
(9.168)
La magnitud K se llama un número de
Si r es el vector de posición, E q es un vector complejo constante y É (r) = E 0 e ±;'K’r,
(9.169)
V • É (r) = ± ;K • E (r),
(9.170)
V x E ( r ) = ± ; K x E (r).
(9.171)
mostrar que
Solución:
Si E(r) ='E()±' K ’r, entonces según (3 .1 5 5 -6 ), V ■E (r) = V •(E 0 e t ' K- ') = e ^ K’ rV • E 0 + V ( e ^ ' K-r) • E 0I
(9.172)
V x E ( r ) = V x ( É 0 e ±>K' r ) = e ±iK' rS/ x É 0 + V ( e ±iK,r) x E 0.
(9.173)
Como EQ es un vector complejo constante, se reduce a
www.FreeLibros.me
V • E 0 = 0, V x E 0 = 0. Por lo tanto, (9.172-3)
Aplicaciones a la teoría electromagnética
243
V • E (r) = V ( e±>K’ r) • E o»
(9.174)
V x E (r) = V ( e ±;K’ r) x E 0.
(9.175)
Como K • r = K xx + K yy + K z z, | - ( e ±í K - ' ) = dx
±jKx e ±IK'r
y expresiones similares con respecto a j ^ y z . Por consiguiente, en virtud de (3.102), V ( e ±IK' T) = ± j e ±iK' T(K x i + K y j + K z k) = ± ; K ( e ±>K' r ).
(9.176)
Sustituyendo (9.176) en (9.174—5), V - É ( r ) = + ;e ±,K’ r K • E„ = ± ;K • (É 0 e ^ K-r) = ± ; K - É ( r ) ,
V x E (r) = I je * ! * '* K x E 0 = ± ;K x (E 0
= ± ;K x E (r).
Observando (9.176) y (9 .1 7 0 -1 ), obsérvese que suponiendo (9.169), i.e., si la variación espacial es de la forma e±/K *r, entonces el operador V se puede sustituir por el vector ±/X. PROBLEMA 9.37
Mostrar que el vector de campo (9.169) satisface la ecuación de
Helmholtz (9.161). S o lu ció n :
Por la identidad vectorial (3.163), V 2E = V ( V * E ) - V x (V x É ) .
(9.177)
Si E = E0e±' K ’r, donde E0 es un vector complejo constante, podemos reemplazar Vpor el vector ±/K. Así por (9.177), con el uso de la identidad vectorial (1.83) y / 2 = - 1, V 2E = + ; K ( ± ;K • E ) - ( ± ; K ) x (±;'K x E ) = - ( K - E ) K + K x (K x E )
= - (K • E ) K + (K ■E )K - ( K - K ) É = ~ K 2É.
(9.178)
Por consiguiente, V 2É + K2 E = ( - K 2 + K 2) É = 0.
[9.161]
Una onda plarn (monocromática) que se propaga en la dirección del vector de onda K es un campo armónico representado por E (r, 0 = Re [E (r)
f] = Re [E 0
- K •'>].
(9.179)
La superficie de fase constante se define «¿tí - K • r = c o n st,
(9.180)
La velocidad de propagación v se define como la velocidad a la cual se mueven los planos de fase constante. Una onda transversa es una onda para la cual los vectores de campo magnético y eléctrico E y H son perpendiculares a la direción de propagación K. La impedancia característica V de un medio con constantes e y ¿u es la razón de la magnitud del vector de campo eléctrico E a la del vector de campo magnético H; i.e.,
www.FreeLibros.me
T 244
Análisis vectorial
PROBLEMA 9.38 Mostrar que (a) una superficie de fase constante es un plano normal a K, (b) la velocidad de propagación de una ondaa]plana vp es 1
(9.181)
(c) el vector del campo eléctrico E es normal a la dirección de propagación K, (d) el vector de campo magnético H es perpendicular a K y a E y es „ K x E H = -------- ,
(9.182)
Cúfl
y (e) la impedancia característica r¡ es (9.183)
'- m - /
(a) Si se halla una superficie de fase constante haciendo t = const en (9.180), obtenemos la condición de fase constante
So lu ció n :
K • r = co n st,
(9.184)
que es la ecuación de un plano normal a K. (Cf., problema 6.16.) (b) Si £ es el componente de r en la dirección K (véase figura 9.5), entonces por (1.30) podemos escribir K • r = K | y (9.180) se transforma en cot ~K% = const. Diferenciando con respecto al tiempo t se obtiene d £
cú
dt
K
(o
1
(c) por (9.148), en una región sin fuentes, -
1
jcúE
V x H.
(9.185)
Así que
Figura 9.5
U na onda plana propagándose a lo largo de K.
Reemplazando
V
V -E = — V - ( V x H) = 0. jtú e
(9.186)
- ; K • E - 0,
(9.187)
p o r - /K .
lo que implica que E(r) es normal a la dirección de propagación K. (d) En una región sin fuentes, por (9.147), H =Reemplazando
V
1
JCÚfl
V x E.
(9.188)
p o r - /K , H = - Ti - ( - ; K x E ) = — , ;
[9.182]
lo que muestra que H es perpendicular a K y a E. Por el requisito V • H = 0 podemos obtener también la relación (9.189)
—; K . H = 0, lo que muestra que H es también normal a K. Reemplazando V por -y'K en (9.185), E = --
K x H
www.FreeLibros.me
cú £
íí x K cú
£
(9.190)
Aplicaciones a la teoría electromagnética
245
Así que (9.182) y (9.190) muestran que (E, H, K) forman un conjunto derecho de vectores ortogonales. (Véase figura 9.5.) (e) Por (9.182) y como E, K son ortogonales, | H | = — | K x E | = — | E |. ojfi con Así que la impedancia es |E|
cup.
17 “ ] ¥ [ “ T
9 .10
/ p
V ~e'
” c uype
Problemas suplementarios
PROBLEMA 9.39 Mostrar que en medios isotrópicos homogéneos, E y H satisfacen las ecuaciones de onda no homogéneas,
r , 2 r. ¿ 2)E = p
dt2
V 2 H —p 8 ^ P
dt2
8
^ dt
P
= -VxJ.
PROBLEMA 9.40 Mostrar que en medios isotrópicos homogéneos y una región sin fuentes, E y H satisfacen d 2E
n2n
dE
V E - p 8 --------- aa ------ = 0, '
dt2
n !ll
F
dt
2H
<9H
„
V H —u E --------- u.a ------ = 0. F
dt2
1
dt
PROBLEMA 9.41 Un indicador útil para el campo electromagnético en el caso en que no hay cargas es el indicador de Coulomb, donde V • A = 0. Mostrar que en este caso los potenciales A y <¡>satisfacen d 2A
V A - . s —
dd>
- , J +, e V - .
V 20 = - - . 8 Mostrar también que para que los nuevos potenciales satisfagan la condición del indicador de Coulomb, la función indicadora ip debe satisfacer la ecuación de Laplace V 2 \¡j = 0. Mostrar que los potenciales en el punto definido por el vector de posición r en campos eléctricos y magnéticos uniformes son PROBLEMA 9.42
0 = - E • r,
A = — (B x r). 2
PROBLEMA 9.43
Mostrar que el flujo magnético
y el vector de potencial A están
relacionados por <5 =<|) A - dr, y por consiguiente, que la fem en un circuito fijo C es
www.FreeLibros.me
Análisis vectorial
246
fem = — 4 ) A - dr. dt Mostrar que la fuerza por unidad de volumen fv. llamada a veces la densidad de fuerza de Lorentz, sobre una región de espacio libre (vacío) que contiene cargas y corrientes debidas a un campo electromagnético se puede expresar como PROBLEMA 9.44
fv = p E + J x B . PROBLEMA 9.45
Usando las ecuaciones de Maxwell, mostrar que
t v = - e „ - f ( E x B) + 8 0 E ( V . E ) — dt
I £ 0V ( E 2) + £ 0(E • V )E 2
+ — B ( V - B ) - — V ( B 2) + — ( B - V ) B . Po 2 fio fiD La cantidad e(J E X B = D X B se llama a veces densidad de momentum del campo electromagnético. [Sugerencia: Usar la identidad vectorial g x (V x g) = y ^ ( g 2) - (g • V) g, que puede deducirse de (3.159) tomando f = g] PROBLEMA 9.46 Hallar el campo eléctrico E y el potencial (¡>debido a una línea de carga infinita con densidad de carga por unidad de longitud pr
Respuesta: E = -------- e p,
P¡ 0 ( Pl) - 0 ( p 2) = — -
¿ n Ep
/ P2 ' ln -
2vE
\p¡
PROBLEMA 9.47 Un volumen esférico de radio a centrado en el origen contiene una carga eléctrica de densidad uniforme p Q. Hallar el campo eléctrico E(r) y el potencial (r) debido a esta distribución de carga.
Respuesta: Para 0 < r < a, E (r)
Para r > a ,
= — re „ 3s0
a 3Po E (r) = -------- er, 3 e 0r 2
0 (r) =
^ (r) =
3e0
2
a + —
6e
( r 2 - a 2).
a 3 pn 3e f
Un cable coaxial consiste en un conductor sólido interno de sección transversal circular de radio a rodeado por un cilindro hueco conductor de radio interior b y radio exterior c. Una corriente total / fluye en el conductor interno y regresa por el conductor externo con distribución uniforme. Suponiendo que el eje del cable coincide con el eje z y su longitud es infinita, hallar el campo magnético H dentro de y entre los dos conductores. Respuesta: Fot p < a , H = ------- p e<¿; a < p < b, H = ------ e<¿. 2na2 2 np PROBLEMA 9.48
For b < p < c ,
PROBLEMA 9.49
I
c2- p 2
H = —-------- ------- e ¿ . 2 np c - b
Mostrar que la fuerza sobre un dipolo eléctrico en un campo eléctrico
E es
f = (p •V) E. PROBLEMA 9.50 Mostrar que el momento sobre un dipolo eléctrico colocado en un campo eléctrico uniforme E es
T e = p x E.
www.FreeLibros.me
Aplicaciones a la teoría electromagnética
247
PROBLEMA 9.51 Hallar la expresión para la fuerza ejercida sobre un alambre que lleva una corriente / cuando el alambre se sumerge en un campo magnético B.
Respuesta: f ^ l
ffl d r
x
B.
c PROBLEMA 9.52
En un campo armónico, si el vector complejo de Poynting se define
como P = -
2
(E x H*),
mostrar que V - P + 2 j o > ( ~ ¡i | H
|2 -
i
e | . E | 2j = - |
(E .J *),
donde H* y J* son vectores complejos conjugados de H y J, respectivamente y |H|2= H H * ,
|E |2= E-É*;
dar una interpretación física de esta ecuación. PROBLEMA 9.53 Hallar el campo eléctrico a una distancia a de una línea de carga de longitud infinita y fuerza p¡ coulombs/m. [Sugerencia: Usar la ley de Gauss (9.32).]
Respuesta:E = E p e p = P l/ 2 n Z 0a e p . PROBLEMA 9.54 Mostrar que el campo eléctrico del dipolo eléctrico de la figura 9.3 puede expresarse como
E = — - — (2 co s 6 er + sen 6 efí). 4n-e0r 3 PROBLEMA 9.55
Mostrar que la energía almacenada en un campo electrostático se
puede expresar como
íIUe-d',h
í S h dv’
donde p es la densidad de carga y 0 es el potencial escalar del campo electrostático.
PROBLEMA 9.56
Mostrar que la energía almacenada en un campo magnetostático se
puede expresar como
2JjJB' HdV=2JIC1AdV‘ donde J es la densidad de corriente y A es el vector potencial del campo magnetostático.
www.FreeLibros.me
10 CAPITULO
FORMAS DIFERENCIALES 10.1
Formas diferenciales
Una forma diferencial exterior del espacio 3-dimensional con coordenadas x, y , z es una expresión obtenida sumando y multiplicando funciones de valor real y los diferenciales dx, dy, dz de las coordenadas. Estas operaciones de adición y multiplicación obedecen a las leyes usuales asociativa y distributiva-, sin embargo, la multiplicación no es conmutativa, sino que obedece a la ley anticonmutativa: si reemplazamosx, y , z por jc„. x 2, x 3, respectivamente, entonces d x¡dxj = - d x j d x j para 1 < i, / < 3.
(10.1)
Por la regla anticonmutativa, denotamos esta multiplicación de formas por medio de una cuña A. Si cada sumando de una forma diferencial contiene expresiones p dx¡ donde p = 0, 1, 2 ,3 , la forma se llama forma diferencial de grado p o simplemente, una forma p. Así, una forma una forma una forma una forma
0 I 2 3
es es es es
simplemente una función diferenciable/Or, y, z); una expresión f d x + g d y + h dz\ una expresión f d x A d y + g dy A d z + h d z A dx; una expresión f d x A d y A dz.
Los coeficientes f g. h se suponen funciones escalares de las coordenadas, infinitamente diferenciables. Una forma exterior diferencial es idénticamente cero si y sólo si todos los coeficientes de su definición son idénticamente cero. Como extensión natural de las definiciones anteriores, sobre un espacio «-dimensional con coordenadas.V,, . . . , x „ , una forma diferencial de gradop (o form ap ) es una expresión de la forma »• ■ Ip( x t , 4 , ■ , xn) d x ¡^ A . . . a dx ¡ p , donde la suma se toma sobre todas las combinaciones posibles de ios índices p y los coeficientes an . . . ip ( x t x „ ) se suponen funciones infinitamente diferenciables de las coordenadas (o n variables).
10.2
Suma y producto externo de formas
Las formas diferenciales de la misma clase se suman combinando coeficientes de formas semejantes. Así, en notación de Índices para formas 1, Y ' j i d x¡ +
dxj = £ ( / , + g ¡)d x¡.
(1 0 .2 )
La regla correspondiente vale para formas 2 ó 3. Sin embargo, cuando dos términos
www.FreeLibros.me
248
Formas diferenciales
249
contienen los mismos diferenciales, pero en orden diferente, deben ordenarse sus coeficientes antes de sumar, usando la ley anticonmutativa ( 1 0 . 1 ). La multiplicación se llama a veces producto exterior o producto cuña de formas diferenciales de la misma clase, además de obedecer las leyes as<x-iativa y distributiva, también obedece la ley anticonmutativa ( 1 0 . 1 ) que en notación de índices es dx¡ a dx¡ = - d x t a dXj,
(10.3)
Por esta regla anticonmutativa, debe preservarse el orden de los factores cuando se multiplican las formas. En general,el producto exterior de una forma p y una forma q esuna forma (p + q). Como una forma 0 essimplemente una función, la multiplicación externa por una forma 0 no afecta el grado de una forma. PROBLEMA 10.1 Mostrar que en una forma diferencial las “ repeticiones de los diferenciales son cero” ; i.e., d x j a dx¡
= 0.
(10-4)
Este resultado es una consecuencia natural de la ley anticonmutativa (10.3). Así que haciendo i = j en (10.3),
Solución:
dx¡ A dx¡ = - d x ¡ a dx¡ .
De modo que d x ¡ a dx¡ = 0 .
PROBLEMA 10.2
M ostrar que en el espacio 3-dimensional, todas las formas p con
p > 3 son cero. Esta es también una consecuencia de la ley anticonmutativa (10.3). Un producto de más de dos expresiones dx. debe contener dx. dos veces. Pero de acuerdo con (10.4), las repeticiones dan cero. Por lo tanto todas las formas p con p > 3 son cero. Por ejemplo, como dx A dx = 0, Solución:
dx A dy A dx A d z
= — d x A d x a d y a d z = 0.
Así. en general, el producto exterior de una forma p y una forma q es cero en un espacio «-dimensional si p + q es mayor que n ya que habrá repeticiones. En este capítulo las formas diferenciales están, en general representadas por letras griegas minúsculas. PROBLEMA 10.3
Para las formas diferenciales
(X = x d x - y d y , cú
¡3 = x d x - z d y + y 2 d z ,
= y d x A d z +x d y A d z ,
calcularlos productos exteriores: a a 6.
y
(a) CX a ¡3,
y = z dy,
6 = x y z dx A dy A dz,
(b) a a /3 a y,
(c ) a a co,
y
( d)
Solución: (a) Como d x A d x = d y a d y = 0,
dy a dx
= - d x a dy,
y
dx a dz
- d z A dx, Ct A (8 = ( x d x - y d y ) A ( x d x
- z dy + y 2 dz)
= x 2 dx A dx - xy dy A dx - x z dx a dy + y z dy A dy + x y 2 dx A dz - y 1 dy A dz = x (y - z) d x a d y - y 3 d y a d z - x y 2 d z A d x c
www.FreeLibros.me
=
Análisis vectorial
250
(b ) Como dx A dy A d y y d y A dz A dy contienen repeticiones, ambos son cero y dz A dx A dy = —d x A dz A dy = dx A dy A dz. Por consiguiente, OL A f i A
y = [ x ( y - z ) d x A d y - y 2d y A d z - x y 2d z A d x ] a z d y =
x z ( y - z ) d x A d y a d y - y 2d y A d z a d y - x y 2z d z
= - x y 2z d x A d y
dx
A
dz.
a
(c) Como d x A d x A d z = d y A d y A dz = 0 y a a üj =
a
dy A dx A d z = - dx
Ady
A dz,
(x d x - y d y ) a (y d x A d z + x d y a d z )
= xy dx
a
dx
= ( x 2 + y 2) d x
d z - y 2d y A
a a
dx
dz +
a
x 2 dx Ad y a d z - x y
dy
a
d y A dz,
Como d x A dx A d y A dz y d y A d x A d y A dz contienen repeticiones, ambos
(d)
son cero. Por tanto, a A
0
= (x d x - y d y ) a (xyz dx A
= x 2y z d x
dy A dz)
d x A d y A d z - x y 2z d y A
a
dx
a
dy
A dz
= 0.
PROBLEMA 10.4 m ostrar que
Establecer la ley anticonmutativa para las formas 1 a y (3; i.e., a a /3 = ~/3 a a .
Solución:
(10.5)
En la notación de índices, si las formas diferenciales son a = Y^f¡dXi,
p = Y ^ é j dx j,
i
i
entonces por la ley anticonmutativa (10.3), a A P = (^2
=
U
A 1 ^ 2 é i dXj'j
{‘&1 dx¡ A d x i
^ ‘
’
// f¡ g j dXj a dx ¡
i
i
dx^ j
A
li dx¡
= -jS a a .
PROBLEMA 10.5
Paralas formas 1 a = /, dx
+ / 2 dy + / 3 dz,
( 1 0 .6)
¿3 = g, dx
+ g2 d y + g, d z ,
(10.7)
y
=
dy
A dz
+h2 d z a d x + h 3 d x a d y ,
(10.8)
calcular los productos externos a A ¡3 y a A y y mostrar que estos resultados corresponden a los productos escalar y vectorial de dos vectores ordinarios.
S olución:
Usando (10.3—4),
www.FreeLibros.me
Formas diferenciales
251
O. A /S = ( f l d x + f 2 d y + í 3 d z ) A ( g x d x + g 2 d y + g 3 d z ) = í l g 1 d x a d x + f2é i d y a d x + f3 g l dz. a dx + fl g2 dx
d y + f 2 g 2 d y a d y + t 3g 2 d z a d y
A
+ f l g 3 d x a d z + l 2 g 3 d y a d z + f 3g 3 d z a d z = ( {2 é j -
4 á 2) d y a d z + ( f 3 g í - í ^ g ^ d z a d x
(10.9)
+ (^ ¿2 - 4 á i) d x a dy,
a f\ y = ( f l d x + f 2 d y + f 3 d z ) A (/it d y a d z + h 2 d z A d x + ft3 d x a d y ) = / j / i i d x a d y a d z + / 2/ j t d y A d y +
/i2 d x
a
dz
a
d x + í 2h 2 d y
a
a
dz
d z + f 3h l d z A d y a
d x + f 3h 2 d z
a
a
dz
dz a
dx
+ ( í h 3 d x a d x a d y + t 2h 3 d y A d x A d y + í 3h 3 d z A d x a d y
=
+ í 2 h 2 + í 3h 3) d x a d y A d z .
(flh l
(10.10)
Ahora, se hacen las siguientes correspondencias entre formas diferenciales y funciones vectoriales: forma l a = lv dx + í2 dy + f3 dz <—►vector f = [f1( t2, fs] ,
( 10 .1 1 )
forma 1 ¡3 = g, d x + g 2 d y + g 3 d z <—>vector g = [g1; g 2, é 3],
( 10 . 12 )
forma 2y = hl d y a d z + h2 d z a d x + h3 d x a d y <—►vector h =
[hl , h 2, h 3], (10.13)
Entonces, comparando definiciones de productos exterior y vectorial (10.9) y (2.25), vemos que el producto exterior a A (3 de dos formas 1 (formas 2) corresponde al producto vectorial de f X g. Comparando también las definiciones de productos externo y escalar (10.10) y (2.24), vemos que el producto externo a A 7 de una forma 1 y una forma 2 (forma 3) corresponde al producto escalar de f -h.
Asi', en general a cualquier forma 1 ó 2 corresponderá una función vectorial y a cualquier forma 0 ó 3 corresponderá una función escalar. Las correspondencias entre las operaciones de formas diferenciales y las que incluyen funciones escalares y vectoriales se ilustran en las secciones finales.
10.3
PROBLEMA 10.6 que
Cambio de variables y jacobiano de una transformación
Si se describe una transformación p o rx = x ( u , v ) ,y = y (u , v), mostrar
dx
a
d(x,y) d y = ------------d u d(u,v)
a
(10.14)
dv
donde d(*>y) <9(1!, v)
dx
dx
du dy
dv dy
du
dv
—
*U
*v
(10.15)
yu yv
es el jacobiano de la transformación. Sol uci ón: dx =
A
Como las diferenciales son
du + — d v - x u du + x v d v , du dv
dy dy d y =. — d u + — d v = y u d u + y v d v ,
du
dv
www.FreeLibros.me
252
Análisis vectorial
el producto externo es dx a dy = ( x u d u + x v dv) a ( y u d u + y v dv) = x uy u d u A d u + x u y v d u A d v +
x vy u d v A du + x v y v d v A dv
= (*uyv - x vYu) d u a dv X, ,
X<,
du a dv
yu yv ¿ (x ,y ) du A dv. d (u ,v ) PROBLEMA 10.7 Si jc =x(u, v, w), y =y{u, v, w), z = z(u, v, w) describe una transformación, mostrar que d ( x , y , z) d x A d y a d z = — ------------ d u d (u,v,w )
a
dv
a
dw,
(10.16)
donde dx
dx dv
dx
du d ( * ,y , z )
dy
dy
dv
d ( u , v , w)
du
dv
dw
dz
dz
dz
du
dv
dw
dw Xu
X y
X w
yu
Yv
Yw
Zy
Z w
Zu
(10.17)
es el jacobiano de la transformación. So l u c i ó n :
Los diferenciales son dx , dx , dx , , , , d x = — d u + — d v + —— d w = x u d u + x v d v + x m a w , du dv dw dy dy dy = — du + — dv + du dv
dy — dw = y u d u + y v d v + y w dw, dw
d z = — d u + — d v + —— d w = z u d u + z v d v + z w d w , du dv dw
Entonces, procediendo de manera similar a la del problema 10.6, el producto externo es dx
a
d y a d z = ( x u d u + x v d v + x w dw) A (y u d u + y v d v + y w dw) a (z u du + z v d v + z w dw) =
(
y
v
— y w z u ) ~~ y u ( j ^ v ^ w
+ z u ( x v y w - x wy v) ] c f u
XU Xy
a
dv
—
a
z v)
dw
XW
yu
y v y™
Z u
Zy
du A d v A dw
Z
d (x,y,z) du
a
dv
a
dw.
d ( u , v , w)
PROBLEMA 10.8 Mostrar que (10.14) es equivalente a la ley anticonmutativa (10.3) de los productos externos. Solución:
Intercambiando * e y en (10.14) y por la propiedad del determinante,
www.FreeLibros.me
Formas diferenciales
253
esto es, si se intercambian dos filas de un determinante el valor del determinante cambia de signo, dy
d ( y , x) dx = ---------du a dv
a
d(u,
V)
d (x,y) = - —— —— d u a d v d ( u , V)
= - d x A dy„
(10.18)
Se ve también que (10.14) y (10.16) aparecen en el cálculo elemental, i.e., el teorema sobre transformación de integral múltiple requiere que d x d y se reemplace por 9 (x, y)/d(u,v) dudv cuando se hace la sustitución x - x ( u , t>), y ~ y (u , r). Así que las formas diferenciales se pueden usar para calcular el jacobiano de una transformación. El sistema coordenado cilindrico se define por la transformación x = u eos v, y = u sen v, z = w. [Véase (5.36).J Usando formas diferenciales, hallar el jacobiano 3(x, y , z)/d(u, v, w). PROBLEMA 10.9
Solución:
Los diferenciales son dx = eo s v du - u sen v dv, dy = sen v du + u eo s v dv, dz = dw.
Entonces por (10.16), el producto externo es dx a dy a dz — [(eos v) du - (u sen v) dv] a [(sen v)du + (u eos v) dv] a dw = ( u e o s 2 v du a dv - u s e n 2 v dv a du) a dw = u ( c o s 2 v + s e n 2 v) du a d v a dw = u du a dv A dw d(x,y,z) = ---------------- d u a d v a d w . d ( u , V, w )
Así que el jacobiano es d ( x , y, z) d(u,v,w)
lo que está de acuerdo con el resultado obtenido en (5.37).
10.4
Diferenciación exterior
La derivada exterior de una forma p u e s una forma (p + 1) doj que se obtiene aplicando un operador d para transformar a> en du). Las definiciones de du> para formas del espacio 3 son como sigue. Si la función / diferenciable es una forma 0, entonces d f es la forma 1 (10.19)
d i = ~ dx + ~ dy + f í - dz. dx ay dz
Si lo es una forma l , f , d x + f 2d y entonces du> es la forma 2 doj = d i l
a
+
f \ d z cuyos coeficientes f. son funciones diferenciables,
dx + dí7
a
d y + d(3
a
dz.
www.FreeLibros.me
(10.20)
254
Análisis vectorial
Si co es una forma 2, / , d y A dz + / 2dz A d x + f i d x A d y cuyos coeficientes /,• son funciones diferenciables, entonces du> es la forma 3
da> * df¡ A dy
a
d z + dít
dz
a
dx + dft
a
dx
a
a
dy.
(
10.2 1)
En resumen, ¿{forma 0 ) - forma 1 ¿(forma 1) = forma 2 ¿(forma 2) - forma 3
.
( 10 .2 2 ) (10.23) {10.24)
PROBLEMA 10.10 Mostrar que s i/ es una función diferenciable de las coordenadas, entonces el diferencial d f de / e s una forma 1 . S o lu ció n : Si f es una función diferenciable del espacio 3-dimensional entonces por cálculo elemental, el diferencial d f de / se puede escribir como en (10.19); i.e., ,,
df
df
df
dx
dy
,
.«i
dt .
= — dx + — d y
+— dz,
[10.19]
dz
que es exactamente la expresión de una forma 1 . PROBLEMA 10.11 Usando (10.20), hallar dco para co = y z dx + x 2 dz. Solución: Por la definición de la derivada (10.20), dco
= d(yz)
+ d (x 2)
dx
A
= (y d z + z dy) = y dz
a
a
dz
dx + (2 x dx)
dx + z dy
(y - 2 x) d z
f
a
A
A
a
d x +2 x d x
cfx - z d x
A
dz
dz
dy.
a
PROBLEMA 10.12 Hallar dw cuando co
= f l d x + f 2 d y + f3 d z ,
(10.25)
donde f¡ (i = 1, 2, 3) son funciones diferenciables. Solución:
Por (10.20) y usando (10.19), la derivada exterior es dco
= dti
dx
A
+ dt2
A
dt,
df,
(9/,
—— d x + —
ax
d/2
/d + 1
az
d /2
/3
dt
df ,
d t 2\ dz j
f dd ff 2
ddf jf \
dx
/
d /2
3
dy
dy
\
d
dz) /
/3
a
dy
\
+ t d -z d z / A (dfx
^ r ) dy A dz + [—
+ ( — - — ] dx dx
\
d x + -— d y + — dy dz
\d x
dz
A
d y + - — d z 1a
ay
— dx
dy + d / 3
\dz
A
d t\ - —
dx )
idz
a
dx
dy.
(10.26)
dy
PROBLEMA 10.13 Hallar dco para co dado por co = /j d y A d z + f 2 d z a
dx
+ /3
dx
a
dy,
donde f¡ (i = 1, 2, 3) son funciones diferenciables de x, y, z. Solución:
Usando la definición de derivada externa (10.21) y (10.19),
www.FreeLibros.me
(10.27)
Formas diferenciales
= dfí
dco
a
dy A d z
í dfl
+ df2
dfí
a
dz
dll
255
+ df3
dx
A
a
dy
a
dz
dl2 \
df2
+ l — dx + — dy+ — \dx dy
dz] dz /
A
- — dx + — dy + — d z ) dx dy dz
A
dL
c/x
\
= I — d x + — d y + —— d z I A d y \ox ay dz /
/ df2
A
d z A dx
dL
dU
dfl d(2 d(3\ - — + — + — ) dx dx dy dz!
A
dy
dii
A
dy
(10.28)
dz.
A
PROBLEMA 10.14 Si co = x y dy A dz + x dz A dx + 3z x dx A dy, hallar dco. La derivada exterior es
Solución:
da>
=d (xy)
a
dy
a
dz
+ dx
a
dz
dx
a
+ d (3 zx ) a d x A d y
=
(x dy + y d x ) a d y a d z + (3z dx + 3x dz) A d x a d y
=
(y + 3x) d x A d y a d z .
f
.d f.
fu n c ió n
fo rm a
fo rm a
0
grad
escalar
La operación sola de diferencistción en el sistema de l en tumo a las operaciones de tomar el gradiente de una fur rotacional y la divergencia de una función vectorial. Esto i la figura 10.1.
alar y tomar el
_ — fo rm a
en
Aplicando un operador d a la forma 0 f( x , y, z) que corresponde a una función escalar, obtenemos la forma 1 [10.19]
+ f'-d y M d z dy dz
que a su vez corresponde a la función vectorial V f = [d f/d x , df/dy, 9 //9 z ] , el gradiente de f tal como se define en (3.110). Luego, aplicando el operador d a la forma 1 a> = (l dx + t2 dy + f 3 dz
[10. 25]
que corresponde a una función vectorial f = \ f x, f 2, / 3 ], obtenemos, la forma 2 d(3 dy
dí2 dz
dy A dz +
( d ll d{3\ ( — - — j dz
A
( d l 2 d iA dx + ( — - — ) dx dx
dy
A
dy.
[10.26]
Así que (10.26) corresponde a la función vectorial V X f, el rotacional de f definido en (3.138). Finalmente, si hacemos corresponder la forma 2 a> = f l d y A d z + í 2 d z a d x + f 3 d x a d y
[10,27]
a la función vectorial f = [ / i , f 2, f^ ] , entonces, aplicando el operador d , obtenemos la forma 3 dfl da>
df2
di-,
— + — +— dx d y d z,
dx
A
dy
A
dz,
[10.28]
que corresponde a su vez a la función escalar V * f, la divergencia de f. [Cf., 10.28) y (3.127).]
www.FreeLibros.me
1
fo rm a
2
f
ro t f) fu nción
vecto rial
fu n c ió n
vectorial-
Solución:
da> -
vecto rial
____ dcoí__ - V x f (
PROBLEMA 10.15 Obtener una relación entre las derivadas exteriores de una forma 0 y gradiente, una forma 1 y rotacional y una forma 2 y divergencia.
dx
f)
fu n c ió n
1
“2---------dco. -------fo rm a 2
fo rm a 3
'V • f (div f) fu n c ió n escalar
Figura 10.1
Análisis vectorial
256
Según las definiciones y ejemplos anteriores sobre derivadas exteriores de I vemos que, en general las derivadas exteriores de formas satisfacen la propiedad
(10.2?)
donde a , 0 son formas arbitrarias y a, b son números. PROBLEMA 10.16 Si a y (3 son formas 1, mostrar que
d ( a A/ 3) = d a A / 3 - ® A d / 3 . So l u c i ó n :
(10.30)
En la notación de índices, sea a = Y 2 f i d x i, i
P = J ^ g j d x j, j
(10.31)
donde f¡, gj (1 < i, j < 3) son funciones diferenciabas de las coordenadas. Entonces por la definción de una derviada externa, do. =di¡
a
dx¡,
d /3 = E
dg¡
/
a
dXj,
(10.32)
;
Ahora,
aA^=(2]fi dXi)A(? g¡ dx¡
)'? ?
fjg j dXj a dx¡.
Usando la ley anticonmutativa (10.3) y (10.31—2) la derivada exterior es d (a
a
/3) =
=
E
E
d ( f j g j ) a dXj A dXj
i
i
i
] L,
**
E E fe
adx¡adx>
■ E E Í E l ; á’ H f j \ k + E E Í Z ! i i \ k
A * , A ‘'’‘'
A d x , A d X,
-EE(E £ + E E ( E
E E (d/j i ; ^
A
^ * » J a í ,<íx, a * ,
dxy)
(df, A dx,) A ^
A gy
dx j - E E f i dXj i i
á , dx¡ - E
f¡ dxj A
A
(cfg;
A <ÍX; )
(dgj a dx y)
= d a a /3 - oc a d/3. PROBLEMA 10.17 Verificar (10.30) con a = / d x y 0 = g dy, d o n d e / y g son funciones
diferenciables de las coordenadas x, y, z.
www.FreeLibros.me
Formas diferenciales
257
Como a A (3= f d x A g d y = f g d x A dy, la derivada exterior es
Solución:
d(0 C a jS) = d(fg) a dx a dy d(íg) d(fg)) dd(fg) (fé) — ---- dx + — ----- dy + — ----- dz . dx dy dz d ({g)
1
a
dx A dy
dx a dx A dy
dz df dg\ — g + f — dz dz
dx a dy a dz.
Ahora, d<X = d f A dx = [ ~^~ dx + — dy + — d z ) A d x \dx dy * dz 1 df
— dy
df
dx + — dz a dx,
A
dy
dz
y por lo tanto, = ( ~ - dy a dx + dz a d x ) a g dy \dy dz
da a
dt
g dz A dx
dy
A
— g dx A dy A dz.
dz
De modo similar, ( dg
dg
dg
dfi = dg A dy = I — dx + \dx
— dy +— dz )a dy
dy
dz
dg
dg
dx
dz
/
= ——dx A dy + — dz A dy,
y por consiguiente Ot A
d/3
= f
¡dg
dx
A l —
dg
\
dz
/
dg
dx A dy +——dz A dy 1 = - f — dx
\dx
A
dy A dz.
dz
Así que (d f
d(X
a
/3 -
a
A
dfi = ( — \d z
dg\ g + f — )
dz/
dx
a
dy
a
dz
=
d((X
A
¡3).
[10.30]
Si a es una form ap y ¡3 es una forma q sobre un espacio «-dimensional, la generalización de (10.30) es d (a
a
0) = da
p +. ( - i ) p a
a
a
d p.
(10.33)
PROBLEMA 10.18 Demostrar (10.33). S o lu ció n :
Si a y (3son monomios dados por a - f dx ¡
'1
a
P = g dxh A
...
a
dx ¡ ,
’p ’
(10.34)
...
A
d x jq,
(10.35)
entonces sus derivadas exteriores son
www.FreeLibros.me
Análisis vectorial
258
d a = df
a
dx¡i
d/3 = dg
a
dx;i
...
a
a
dx¡p)
(10.36)
A c?Xy9 ,
(10.37)
a
...
Ahora, con uso repetido de la ley anticonmutativa (10.5), d (0C A /S) = d (/g )
A
d x ¡i A
...
= ( f dé + é df)
A
= Ce//
...
A G?X¿1 A
a
dx,
1
A
c/xip)
A
+ ( - l ) p (/ dx¡t A . . . = d a A /S + ( - l ) p
dx¡p
A
dx¡i
A ...
A
dx;
A
A
’p
(g c/x;i
dx i p)
A
A
dx 1 A 11
A ...
A (dg
dxj(j
A
A
c/x;i
...
A
• dx/q)
A
...
dx;
'q
Ac f Xj
[10.33]
a A d/8.
Por consiguiente, queda demostrado (10.33) para monomios a y j3. El caso general se sigue por linealidad.
10.5
Lema de Poincaré
El lema de Poincaré establece que si cj es cualquier forma de funciones continuamente diferenciables, entonces (10.38)
d(du>) = 0 . PROBLEMA 10.19 Demostrar que para cualquier forma U co, II
O
•i
Solución:
(10.39)
Como una forma 0 es una función escalar, sea co = f( x , y, z),
donde f(x , y, z ) es una función continua dos veces diferenciable. Si denotamos la diferenciación parcial por una notación con subíndices(p.e. df/dx = f x , d2f/ d x d y = f x y , etc.), entonces por (10.19), (10.40)
da> = df = f x dx + f y d y + fz dz. Así que d(da>) = d (d f) = d f x
A
dx + d fy
A
dy + d f z
= (f xx dx + f yx dy + f zx dz)
({ y z
~ iz y )
dy
A
dz
dx + (f xy dx + f yy dy + f zx dz)
a
+ ( f x z dx + f yz dy + f zz dz) =
a
a
a
dy
dz
d z + ( f zx - f x z) d z
a
dx + (f xy —f y x) dx
= 0,
A
dy (10.41)
pues todos los coeficientes son cero porque f xy = f y x , etc. De modo similar, se puede mostrar que (10.38) vale también para cualquier función continuamente diferenciable de n variables. PROBLEMA 10.20 Demostrar que para cualquier forma lw , d(d<ü) Solución:
=
(10.42)
0.
Usando notación de índices, una forma 1 co se puede expresar como (Q =
f i dx¡,
(10.43)
donde f.(i = 1, 2, 3) son funciones diferenciables de x .{i= 1, 2 ,3 ). Ahora, por la definición de derivada exterior, da> =
www.FreeLibros.me
df¡
a
dx j .
(10.44)
Formas diferenciales
259
Así, usando (10.30), d(dúi) = y
' d (df j A dxj)
= ^ ^ [ d ( d / ;) A d x ¡ —d f ¡ A d (d x ¡)].
(10.45)
d (d /¡) = 0
(10.46)
Pero por (10.39), y
d ( d x ;) = 0.
P or consiguiente, d ( d o j ) = 0. PROBLEMA 10.21
Si co es una form a 2, m ostrar que d (d w ) = 0.
S o lu c ió n :
(10.47)
Si co es una form a 2, entonces el grado de d(dco) es 4 , que excede a 3; por
tan to , según el resultado del problem a 10.2, d(da>) = 0 .
PROBLEMA 10.22 D em ostrar (10.38); i.e., m ostrar que para cualquier form a diferencial co, d(d
[10.38]
Supongam os que una form a p co sobre el espacio «-dimensional es un
m onom io dado por dxj^ a . . .
a> = f
a d x ip,
(10.48)
donde f e s una función de n variables continua y doblem ente diferenciable. Ahora, usando (10.33) d ( d c o ) = d ( d f a d x /i A . . . A d x ¡ p ) = d ( d í ) a d x ¿i a . . . A d x ¡p - d i a d ( d x ¡i a . . . a cfx¡p)
= -d /
a
d ( d x ¡(
a
...
a
d x íp)
(10.49)
pues d (d /) = 0 por (10.39). Así, es suficiente m ostrar que d ( d x ¡í a . . . a d x i p ) = 0.
(1 0 .5 0 )
D em ostram os (1 0 .5 0 ) po r inducción sobre p. Supongam os que d(dx¡ia ...
A d x j p _ 1) =
0.
(10.51)
Entonces, usando o tra vez (10.33), d ( d x fi
a
. . . a d x íp) = d ( d x ¡ ¿ )
pues d( dxn ) = 0 por (1 0 .4 0 ) y
a
d x ¡2a . . .
a
dx¡
- d x ¡ i A d (d x ¡2 a . . . A d x jp) = 0
d(x¡2 A . . . dx¡p) = 0 por la hipótesis (1 0.51).
Esto com pleta la dem ostración de (10.47) cuando co es un m onom io. El caso general se sigue por linealidad.
PROBLEMA 10.23 M ostrar que el lem a de Poincaré (1 0 .3 9 ) para una form a 0 co corresponde a la fórm ula vectorial rot (grad /) = V x (V O = 0. S o lu c ió n : C om parándolos resultados del problem a 10.19 y el problem a 10.15, o usando la figura 10.1, vemos que la relación d(d
[10.39]
vale para una form a 0 co y por ta n to , corresponde a la identidad vectorial V X ( y f ) = 0. [Véase (3.142).]
www.FreeLibros.me
260
Análisis vectorial
PROBLEMA 10.24 Mostrar que el lema de Poincaré (1 0 4 2 ) para una forma 1 lo corresponde a la fórmula vectorial div ( ro tf) = V - ( V x f ) = 0 . So l u c i ó n : La figura 10.1 muestra que si lo es una forma 1, dLo corresponde al rotacional de una función vectorial y d(dcj) a su divergencia. Así, la relación d(dcü) = 0
[10.42]
vale para una forma 1 co y por tanto, corresponde a la identidad vectorialV (V X [Véase (3.143).]
10.6
f)
= 0.
Invariancia de derivadas externas bajo transformaciones n T es una regla que asigna a V , un punto (_v, , . . . , x „ ) en un
Un cada p espacio «-dimensional Vf ; i.e..
Tt vm — ¡ | IJna t
s coordenadas como
«»m). * = 1,2,
(10.52)
■ .iii'tformai.iún T es dijen ní table si las funciones coordenadas ( 1 0 * :) son continuam ente diferenciabas. ■de uno y i m to en Vm a a T 1 Je T existe si 7 'es y se d
o
i ■ BM M i
: V
—* V sobre co ti u) i .? , ...........y i i c i ............... dxn siempre que aparezcan por las funciones coordenadas ( 10.5 2) y luego sitri o el resaltado por medio de algebra de formas. s, i.e., suma, producto cufia y (a + /
+ p*,
(10.53)
(dw)* = d(a>*).
(10.55)
Obsérvese que en (10.53) las formas diferenciales a y (3 deben tener los mismos grados, pero en (10.54) a y 0 pueden tener grados diferentes. La ecuación (10.55) muestra que el operador d es invariante bajo una transformación diferencial. PROBLEMA 10.25 Si la transformación T es x = u + v,
y = u-v,
y a = xy, dx, j3 = y dy, mostrar que (a + ¡3)* = a * + ¡3*, Solución:
( a a ¡3)* = a * a /3*,
(d a)* = d (a * ).
Para a, |3, T dados,
0C* = (u + v )(u - v ) d ( u + v) = (u2 - v 2)(du + dv) = ( u 2 - v 2)du + ( u 2 - v 2)dv,
¡3* = (u - v ) d ( u - v) = (u - v)(du - dv) = (u - v )d u - (u - v)dv.
www.FreeLibros.me
Formas diferenciales
261
Entonces, (CX + /3)* = (x y dx + y dy)*
= (u + v ) ( u - v ) d ( u + v) + (u - v ) d ( u - v) = a * + /3*,
(Ot a /3)* = (x y dx a y dy)* = (u + v) (u - v) d ( u + v) a (u - v) d ( u - v) = a*
a
/S*.
Como c/Ot = d (x y dx) = d (x y ) A dx = (x dy + y dx) A dx = - x dx
a
dy,
tenem os
d(0.*) = d [ ( u 2 - v 2) du + ( u2 - v 2)d v ] = d ( u 2 - v 2) a du + d ( u 2 - v 2) a dv = (2 u du - 2v dv) a du + (2u du - 2v dv) a dv = - 2v dv a du + 2u du a dv
= 2 (u + v) du a dv, y por consiguiente,
(d(X)* = ( - x dx A dy)* = - ( u + v ) d ( u + v) a d ( u - v) = - ( u + v )(d u + dv)
a
(d u - dv)
= - ( u + v ) ( - d u a dv + dv a du) = 2 (u + v) du A dv = d(OL*).
PROBLEMA 10.26 V erificar la relación de invariancia (1 0.55), esto es, (dw )* = d(co*) cuando cuando co = x dy A dz y T es la transform ación x = u + v - w, y = u 2 - v , z = v + w 2 S o l u c i ó n : Com o dco = d ( x dy a dz) = dx a d y a dz,
(da>)* =
d(u +
v - w) a d ( u 2 - v) a d ( v + w 2)
= (du + dv - dw) a (2u du - dv) a (dv + 2 w dw) = (2 u dv a du - 2u dw a du - du a dv + dw a dv) a (dv + 2 w dw) = - 2 u dw a du a d v + 4uw dv a du a dw - 2 w du a d v a dw
= - 2 ( u + w + 2uw)du a d v a dw. Ahora, o>* = (x dy a dz)* = (u + v - w ) d ( u 2 - v) a d ( v + w2)
= (u + v - w)(2u du - dv) a (d v + 2 w dw) = (u + v - w)(2u du a d v - 4uw dw a du - 2 w d v A dw) = (2 u 2 + 2uv - 2uw) du a dv + ( - 4 u 2w - 4 u v w + 4 u w 2) d w a du
+ ( - 2 uw - 2vw + 2w2) d v a dw,
www.FreeLibros.me
262
Análisis vectorial
y por consiguiente la derivada exterior es d ( a )*) = d ( 2 u 2 + 2 u v - 2 u w )
du
a
a
dv
+ d ( - 4 u 2w - 4 u v w + 4 u w 2) a d w A d u + d ( - 2 u w - 2 v w + 2 w 2) A d v A d w = ( 4 u d u + 2 u d v + 2 v d u - 2 u d w - 2 w du) a du a d v
+ (-
8 uw
du - 4 u2 d w - 4 v w d u - 4 u w d v - 4 uv d w + 4 w 2 d u
+ 8 u w d w) a d w a d u + ( - 2 w du - 2 u d w - 2 w d v - 2 v d w + 4 w dw) = - 2 u dw
a
du
a
d v - 4uw dv A dw
= —2 ( u + w + 2 u w ) d u
a
dv
a
a
a
dv
dw
a
du - 2w du A dv
a
dw
dw.
Así, hemos verificado que (dw)* = d(co*).
PROBLEMA 10.27 Si T : Vm —> V es una transformación diferenciable y la función f ( x i , . . . , x ) es una forma 0 sobre V mostrar que (di)* = d({*).
( 1 0 .5 6 )
Si T : Vm —►Vn está descrita por ( 1 0 .5 2 ) , i.e.,
Solución:
= x / ( t / 1( • • ■, u m ),
i = 1 , 2 , ■■ ■, n,
entonces (^1 f
* ’ t ^m) =
Um ),
'
, X n ( U j, • • • , Um )].
Como / es diferenciable, la derivada es V dl=
d f ( Xl, ■ ■ ■ , x n) ----------— ------ — d x , .
dx,.
Por consiguiente, d f [ X l ( ul t - ■• , u m ), ■■■ , x „ ( u „ . • • , um)] dx¡
--------------------------------------- dui dx¡ du¡i
(di)* = £
du¡ d((*).
PROBLEMA 10.28 Si gj es una forma 1p sobre Vn Jy T . V m —■* Vn es diferenciable,9
mostrar que [10.55]
( d o j ) * = d(a>*).
Demostramos (1 0 .5 5 ) por inducción. Se ha mostrado en el problema 10.27 que (1 0 .5 5 ) es verdadero cuando cj es una forma 0. Ahora supongamos que (1 0 .5 5 ) es verdadero cuando co es una forma (p - 1). Entonces es suficiente demostrar (1 0 .5 5 ) para una forma p del tipo Solución:
co
a
dxh
Por ( 10.33), la derivada exterior es d(o
a
d x ¡ ) = dco
a
d x ¡ + (—l ) 13- 1 co
= dc o a d x ¡
pues d (d x ¡) = 0 por (10.38). Así
www.FreeLibros.me
a
d(dx¡)
263
Formas diferenciales
[ d ( ( ú a cfx,-)]* = ( d( ú a d x ¡ ) * = (da>)* A (d x ¡)* = d(a>*) A d ( x f ) = d[&>*
a
(c/x,-)*]
= d[(a> A d x ; ) * ] .
La relación (10.55) muestra también que la derivada exterior de una forma diferencial es independiente del sistema coordenado en el cual se calcula.
10.7
Integración de formas
Una curva suave C del espacio en el espacio 3-dimensionai V3 con coordenadas x, y , z se representa por las ecuaciones paramétricas x = x( t ) ,
z = z (t ),
y = y(f),
(10.57) (a)
con a < t < b . Usando el concepto de aplicación (o transformación) discutido en la sección 10.6, (10.57) puede describirse por una transformación diferencial T como T: V, — VV
Así la curva C en V3 es una aplicación de un intervalo cerrado [a, b ] sobre la recta real en V3. (Véase figura 10.2). Si oj es una forma 1 sobre V3 , la integral de w sobre una curva C se define por |
CU »
I
1: 1* ,
(10.58)
Je donde oj * es la forma I transformada de co por T. Así, si co se expresa como ío
» (/, dx + /j d y + /, dz),
donde las/,, (i - 1, 2 ,3 ) son funciones diferenciables de x, y, z, entonces (10.58 ) se expresa como /, dx • /, d y *
L " L
dz - I
<■>*
(b )
JM
r
F ig * 1 0 .2
U * ( t) x '( o + / ,* « / ( < ) * n ( t ) * '( t )J d t, (10.59)
donde/? { t) ~ f ( [x {t),y (t), z ( t ) \ ,x ' ( t) ~ d x ( t) /d t, etc. E¡ lado derecho de (10.59) es sólo una integral ordinaria. PROBLEMA 10.29 Si co = x d x + y d y + x y z dz y C es una curva representada por la transformación x = t, y = t, z = t con 0 < t < 1 , calcular I a>. Je So lu ció n : Como oj* = t2 dt + t d t + t 3 dt = (f 2 + t + t 3)dt, entonces por (10.59),
, , 3N , 1 1 1 13 r „ . f „ . . r (í 2 + t + f3) dt = - + - + - = — . 3 2 4 12 Je •'•'[[00,1 o . 1 ]]
f
PROBLEMA 10.30 Si / es una forma 0 sobre V3, i.e., una función de x, y, z, y C es una curva en V3 representada por (10.57), mostrar que
www.FreeLibros.me
M a p p in g o f th e in t e r v a l on V l in to V 3 .
[a, t
264
Análisis vectorial
f d f= f* (b )-f* (a ),
(10.60)
•Je
donde f* ( t) = f [ x ( t ) , y ( t ) , z ( t ) ] . So l u c i ó n :
Por la definición (10.58) de una integral y usando (10.56), f
cf/ = f
J e
c /(/* ) = f b —
(c/ / ) * = f
* ^ [ a ,b ]
^ [a ,b ]
*/a
( / * ) efí = / * ( 6 ) - f * ( a )
^
por el teorema fundamental del cálculo. PROBLEMA 10.31 Considérese una curva C en un plano representado por la transformación x - t, y - t 2 con - 1 < t < 1. Entonces (a) si co = y 2 d x + 2xy dy, calcular‘ f oj, Je y (b) hallar una función / tal que d f = co y verificar (10.60). Solución: (a) Como co* = t* dt + 2 í 3 d (t2) = 5 í 4 dt, entonces, usando (10.59),
i f co = f
J '[r-i.ii - i.i]
• 'c
co* = f 1 5 t 4 dt = f 5
J-i
= 2.
(b) Para u>=y2 dx + 2x y dy, sea f(x , y ) una función tal que d f = co. Entonces, df Sf d f = —— dx + - — dy = y 2 d x + 2 xy dy. dx dy Así que, dí 2 df — = y \ — = 2 xy. dx dy
(10.61)
Por la primera ecuación de (10.61), í( x , y) = xy 2 + á (y ), y por consiguiente, df 9 ,, , — = 2 xy + g (y). dy Comparando esto con la segunda ecuación de (10.61), tenemos g ' (y) = 0, estos es, £(y) = K = constante. Así, f(x , y) = xy 2 + K, f* ( t) = /[ x ( í) , y (í)] = t5 + K. Por (10.60), f
Je
co= [ d f = f * ( l ) - f * ( - l ) = ( l + K ) - ( - l + K ) = 2, Je
lo que está de acuerdo con el resultado de la parte (a).
Una superficie suave S en el espacio 3-dimensiónal V3, representa por las ecuaciones paramétricas x = x (u , v),
y = y (u, v),
* = *<«,
co n a < u < b , c < v < d . Así,análogamente,(1 0 .6 2 )se transformación diferenciable
www.FreeLibros.me
x , y , z , se
Formas diferenciales
265
T Vt - + V}i ■ por consiguiente, una superficie S en es una aplicación de un rectángulo cerrado D: a < u < b, c < v < d en el espacio 2-dimensional V2 con cordenadas u, v, en V3 (véase figura 10.3). Si co es una forma 2 sobre V3 , entonces la integral de co sobre una superficie S se define por . . < '0 63,
que co* es una form a 2 transform ada de co po r T y el lado derecho de (1 0 .6 3 )
■
es una integral m últiple ordinaria.
PROBLEMA 10.32 Si co = x d y A d z + y d x A d y , y S es la superficie representada por la transformación x = u + v, y = u - v, z = uv, con 0 < m < 1 ,0 < v < 1, calcular
JJ'
(a )
Como la forma 2 transformada de co por T es
Solución:
0
entonces, por (10.63),
JJ
co
= JJ
co* =
J
J
(u2 - 2 u + v 2 + 2 v + 2 u v ) du d v
v 2 + 3 v ------d v
r
Jo 7
6
Fig ura 10 .3
A p lic a c ió n de un re ctáng ulo cerrado de V2 en una
*
su p e rfic ie d e K 3 .
co donde co = x y dy A dz + x dz A dx + 3z x dx A dy,
PROBLEMA 10.33 Calcular
y S es la superficie dada por z = x 1 + y 2 con 0 < x < l , 0 < j c < 1 , 0 < j / < 1 . Como S se puede describir por las transformaciones
So l u c i ó n :
x = u,
y = v,
z = u2 + v 2
con 0 < u < l , 0 < v < l , co* = u v d v a d ( a 2 + v 2) + u d i u 2 + v 2) a d u + 3 u ( u 2 + v 2) d u a d v = uv dv
a
( 2 u d u + 2 v d v ) + u (2 u d u + 2 v d v )
= ( 3 u 3 + 3 u v 2 - 2 u 2v - 2 u v ) d u
a
a
d u + ( 3 u 3 + 3 u v 2) d u
dv,
y por ( 1 0 .6 3 ) ,
J J
co
~JJ -J J a)*
( 3 1/3 + 3 u v 2 — 2 u 2v - 2 u v ) d u d v
■i"(H _5_ 12'
www.FreeLibros.me
a
dv
266
Análisis vectorial
La frontera dS de S está formada por las curvas del borde de S que se describen por / i
3D 3
d
dS¡:
x = x ( u , c),
y = y (u, c), z = z ( u ,c ) ,^
dS{.
x =x (b , v),
y « y ( 6 , v),
dS,:
x = x (u ,d \
y =y (u ,d ),
x =x (a , v),
y = y (a , v), z = z ( a , v ) ,J
z z(b,
V ),
I
(10.64)
z = z(u ,d ),
con a < u < b , c < v < d. Entonces la frontera dS de S se expresa por íD ,¿
i
i¡ 1 3 D ,
D
<9S = d S t + ¿ S , - d S , -
c
,
ao,
(o )
(10.65)
Las cuatro curvas de borde (10.64) son las aplicaciones de los cuatro segmentos de recta que forman la frontera del rectángulo D de V2 en K ,. (Véase figura 10.4.) Los signes menos antes de dS3 y dS* en (10.65) se deben al hecho que un viaje consistente alrededor del borde de D, y por tanto de S, se asegura cambiando el sentido de ios segmentos de recta dD3 y 9Z rel="nofollow">4 . (Véase figura 10.5.) Si <jú es una forma 1 sobre V3 , entonces la integral de co sobre la frontera 35 de S se define como I
Jas
(O =
I
CO +
Jas t
I
(0+ 1
= |
O) +
Jas,
Jast
ÚJ + I
j
J-ds% co -
j
co -
co
J -is, I
( 10 .6 6 )
co
Jast
**3Sj
pues la dirección en la cual se toma la ruta de integración se determina por la relación (10.67)
(O.
f • dSf f dSt W = ~ ^dSjf JdS El teorema de Stokes expresa que si co es una forma 1 sobre V3 y 5 es la superficie S 2-dimensional sobre V3, entonces (b )
( 10 .6 8 ) Figura 10.4
A p lic a c ió n de la fro n te ra de
9
un re c tá n g u lo de V 2 en V 3 .
*
i r
PROBLEMA 10.34 Verificar el teorema de Stokes (10.68) para formas 1.
Se expresa la forma 1 co de V3 como
So l u c i ó n :
= f l d x + f2 d y + f 3 d z ,
(10.69)
donde las f¡ (i = 1,2, 3) son funciones diferenciables de x, y , z. Entonces laforma 1 co*, transformada de co por T : V 2 —* V3 se puede expresar como co*
- f ( u , v) d u + g ( u , v ) d v ,
(10.70)
donde / y g son funciones diferenciables de u, v, i.e., de las coordenadas de V2. definición ( 1 0 .2 0 ) de la derivada exterior de una forma 1 , d(co*)
ídg
df
\d u
dvl
Por la
(10.71)
du A dv.
Ahora, por definición (10.63) de la integral y usando (10.55) y (10.71),
// da -JJ w")* D S
- ff D
D
(ff - f )
' / / lí
www.FreeLibros.me
M
“ /
/
dUch'-
<10 72)
Formas diferenciales
267
A a * iu < :b , c < f < d, entonces, integrando primero con respecto a u, j y
dadv , jf
| í
|á
/(v ) dv,
rfpdv.jf'
d donde / (v) =
í
b dg(u,v)
- dD3
-dDs
du.
du
Como v es una constante en la integral parcial I(v), el integrando es la derivada ordinaria con respecto a u. Así, por el teorema fundamental del cálculo, g ( b, v ) - g
/ (v ) =
D
¡>dD2
c
(a , v ).
u
O
Por tanto, J J ^
dudv
=
J
g (b ,v)d v~ J c
d
g(a,v)dv.
(a)
(10.73)
c
Ahora sobre la curva 9D 2 en 9D, como du = 0, (10.70) se reduce a co* =g(b, v)dv. (Véase figura 10.4.) Así, por definición (10.58) de la integral í
g(b,v)dv=f
cú* =
co.
f
c
JdS,
Por un argumento similar g (a, v ) d v =
f
•Je
co* =
í
f
JdD.
Cú.
•JdS.
Por consiguiente, sustituyendo (1 0 .7 4 -5 ) en (10.73),
lili * * - /dS2
(10.76) Figura 10.5
dS4
que se muestra en la figura
De modo similar, integrando primero con respecto a v,
10.4.
(10.77) dS j
ccS j
Así, sustituyendo (1 0 .7 6 -7 ) en (10.72) y usando (10.66), obtenemos el resultado requerido
JJ
J J J J
dcú =
s
Cú
ds2
dS.
-
J
—
Cú ~
as,
dS7
Cam bio del sentido de los segmentos de la aplicación
Cú +
as3
dS*
(ú
as,
dSa
cú.
es
En el problema 10JO, si los dos extremos de la curva C son la frontera 9C de C, entonces el teorema de Stokes para forma 0 <10.60) se escribe como (10.78)
www.FreeLibros.me
268
A nálisis vectorial
La integral de una forma 0 en un punto/? es &> = « ( p)
I
(10.79)
en P
Como está fuera del propósito de este libro definir la integral de una forma p sobre una región p-dimensional de un espacio «-dimensional V , es suficiente establecer el teorema general de Stokes: Sea u? una forma diferencial de grado (p - 1) definida sobre alguna región acotada p-dimensional M en V (p < n) con una superficie (p - 1 )-dimensional como frontera suave bM. Entonces, [d co M
- f w J¿M
(10.80)
El teorema general de Stokes muestra la relación entre la integral de una forma diferencial w y la de su derivada exterior du>. El teorema general de Stokes está relacionado con los teoremas de Gauss, clásico de Stokes y de Green por el teorema: Sea f - \J\ , f 2, f i ] una función vectorial continua y diferenciable en una región R del espacio 3-dimensional V3 y sea co, = / , dx + f 2 d v + / 3 dz la forma 1 correspondiente. Entonces, comparando = í f, dx + f2 d y + fj d z = I Je J r. J°
í
U * ( t ) x ( t ) + í2* (t)y '(t) + i * ( t ) z ( t ) \ dt
j f - c f r = f ( / , dx + /j dy + f, dz), Je Je
[10.59]
(4.151
obtenemos
/ - /-Jec I d r.
(10.81)
Je
Así, la integral de la forma 1 sobre una curva se llama integral de línea. Si a >2 = / , d y A dz + j \ dz A dx + d x A £/>■ es la forma 2 correspondiente a f, cuya integral de superficie se puede expresar comoJ"J"f * dS por (4.31), entonces
s
Js J - J Ts
í-d S .
(10.82)
Así la integral de la forma 2 sobre una superficie se llama también integral de superficie. Si w 3 = f d x A d y A dz es una forma 3 correspondiente a la función escalar /, cuyo integral de volumen es
f d V por (4.43), entonces JJJ R
ÍÍÍ^ IIJ í d V . R
(10.83)
R
Por consiguiente, la integral de una forma 3 sobre una región/? corresponde a una integral de volumen.
www.FreeLibros.me
Formas diferenciales
269
PROBLEMA 10.35 Una forma 2 sobre V 3, con coordenadas x, y, z dados por cj 2 = f, dy a dz + f2 d z a dx + / 3 dx a dy
f
corresponde a una función vectorial diferenciable = [ / i , f 2, f 3 ] . Mostrar que (10.80) corresponde al teorema de Gauss o teorema de divergencia dado por jj/v .
[4.58]
R
So l u c i ó n :
S
Por el resultado del problema 10.13, la derivada exterior de una forma 2 es fd{1
d f 3\
áfj
d u 2 = 1-T— + - — + \o x dy
- — I dx Ady a azi
dz,
que corresponde a V • f (véase figura 10.1). Así, siM es la región i? en V3 , entonces su frontera 9M = 9R es la superficie cerrada S que encierra R y (10.80) se reduce a
///"“• ■S fa R
(10
dR
84)
Así, con referencia a (10.82) y (10.83), vemos que (10.84) corresponde al teorema de divergencia. Escribiendo (10.84) explícitamente,
JIf §7 + R
+S) * AdyA* ' 5RJí f‘dyA* + f2 d z A dx + f3 dx a dy,
(10.85)
c/zc/x + R dxdy).
[4.66]
que corresponde al teorema de divergencia J 'J 'J ' l~^~~ +
dxdydz =
(P
c/ydz +
PROBLEMA 10.36 Una forma 1 oj] de V3 dada por co¡ = /¡ d x + / 2 d y + / 3 d z ,
, f 2, f 3 ]. Mostrar que (10.80)
corresponde a una función vectorial diferenciable f = corresponde al teorema clásico de Stokes dado por
J J V x t - d S = {j) f - d r . s
Sol uc i ón: es
[4.104]
c
Según el resultado del problema 10.12, la derivada exterior de una forma 1 tdfz
\oy
/d fj
d f 2\
d o Jj = ( ------- —
) dy
dz)
a
dz +
d f 3\
--------- —
\o z
fdf2
) dz
a
dx/
dx +
df^
--------- —— d x
\d x
a
dyj
dy,
que corresponde a V X f (véase figura 10.1). Así, si M es la superficie 2-dimensional finita S en V3, entonces su frontera 9M = 9S es la curva cerrada simple C que limita S y (1 0 .80) se reduce a
J 'j'
d o j, =
J'J"
W j.
( 1 0 .8 6 )
as Así, con referencia a (10.81) y (10.82), vemos que (10.86) corresponde al teorema de Stokes dado por (4.104). Escribiendo (10.86) explícitamente,
www.FreeLibros.me
270
Análisis vectorial
JJ s
d í 2\
dy
dzj
fdíl
dy
dfA
(dí1
d z + I —------- — 1 d z
a
a
dx /
\dz
d{\
dx + I —— I dx -
dí3
\dx
-1
a
dy,
dy
fí dx + í2 dy + f3 dz, '
ds
(10.87)
que corresponde también al teorema de Stokes dado por fj) P dx + Q dy + R d z =
dy dz
c
s
+
dP
dR\
IdQ
dP\
e l - 3 ^ ) d2dx *
dxdy
l4 - 1091
PROBLEMA 10.37 Considérese una forma 1 oj sobre V2 con coordenadas x, y dadas por cú =
P (x, y ) dx + Q (x, y ) dy,
donde P y Q son funciones diferenciables d e x e y . Mostrar que (10.80) corresponde al teorema de Green en el plano dado por P dx + Q dy = J J ( | £ _
[4.110]
R
C
S o lu ció n :
dxdy.
Por la definición (10.20), la derivada exterior de la forma 1 dada es .
(dP
dP
\
fdQ
\d x
dy
j
\d x
d c o = I—— d x + —— d y ) a d x +
,
dQ
\
dy
)
—— d x H-------- a y
a
ay
dQ dP\ -------- — \ d x A d y . dx
( 10.88 )
dx /
Así, si M es la región finita 2-dimensional S sobre V2, i.e., el plano x y , entonces su frontera 9M = 9S es una curva cerrada simple C que lim itaS y (10.80) se reduce a
J
(10.89)
(o = I dco,
dS
Js
o sea que, escribiendo (10.89), explícitamente, J
P dx + Q d y = J
0 ^ -
- |
y
j
dx
a
dy,
(10.90)
que corresponde al teorema de Green en el plano dado por (4.110).
10.8
Formas diferenciales de las ecuaciones de campo de Maxwell
En teoría de campo electromagnético, las ecuaciones básicas de Maxwell en la forma vectorial son
www.FreeLibros.me
Formas diferenciales
271
V • B = O,
[9.18]
V •D -
[9.19]
p.
donde la intensidad de campo eléctrico E, la intensidad de campo magnético H, el desplazamiento eléctrico D, la inducción magnética B, la densidad de com ente eléctrica J y la densidad de carga p son todas funciones de las variables de espacio x, y, z y del tiem p o /. (Cf.,sección 9.3.) Ahora se introducen las formas diferenciales utE - E, dx + E2 d y + E t d z, c¿ h
=
flf,dx + Hj dy 4-
dz,
D,d y a d z + D, d z * dx + D3 dx a
u jb =
B ¡dy a d z + B2 d z a dx + B3 dx a
* Ji dy
dz + J2 dz
a
a
dx + J , dx
dy,
a
PROBLEMA 10.38 Mostrar que (9.16) y (9.18) se pueden expresar como
d a = 0,
(10.91)
y (9.17), (9.19) se pueden expresar como cf/3 + y = 0,
(10.92)
donde 0C —cúg a dt +
/3 = —&>h y = Cúj Sol uc i ón: Ot —
A
= E l dx
(10.93)
dt + cú£) ,
a A
,
(10.94)
dt - u>p.
(10.95)
Sustituyendo los valores de oje y ojb en (10.93), dt + COjg a
dt + E 2 dy
A
dt + E 3 d z
dt + B 1 dy
a
a
dz + B2 dz
dx + B 3 dx
a
A
dy.
Como dx A d x = d y A dy = dz A dz = d t A d t = 0 y d{dt) = 0 por (10.39), entonces por (10.30) la derivada exterior del primer término de lo anterior es d (E l dx
A
dt) = d (E 1 dx) = dEl
a
dt — E¡ dx
A
dx
a
a
d (dt)
dt
(dE l SE, dEl dEl \ - I—— dx h— —— dy + —— d z h — d t) \d x dy dz dt / dEl = —— dy dy ^
A
dx
A
dEl dt H— -— d z dz
A
dx
A
A
dx
a
dt.
De modo similar, las derivadas exteriores de los otros términos son d (E2 dy
A
dE2 dt) - —— dx dx
d (E 3 d z
a
dt) = —— dx dx
A
dy
a
d E2 dt h— -— d z dz
a
dz
A
dt + ------ dy dy
A
dy
a
dz
dE 3
d (B l dy
a
dz)
dB, -
------
dx
a
dy
a
dt,
A
dz
a
dt,
a
dy
a
dz,
dE 3
dx
dBl + ------
dt
dt
www.FreeLibros.me
dt
272
Análisis vectorial
dB2 d (B 2 d z
a
d (S 3 dx
a
dB2
d x ) = ------- d y
dz
A
dy
A
d x + ------- d t
a
dz
a
dx,
d y i------— d t
A
dx
a
dy.
dt
di?, d z /\ d x
dy)
dz
a
dt
Sumando todos estos términos,
da=f d E 3
dE2
dB, dt
dz
dE2
dE,
dx
dy
(dB ,
dx
dE,
dE3
dz
dx
dt +
a
a
dy
A
+ — — -) d x a
dy
a dz,
dt
dB. dz
dt
dx
A
dt
A
dt
(10.96)
dB,
+ ——
\ox
dz
A
dB,
dB2
——
dy
dy
dz
y d a = 0 implica que todos los coeficientes que son escalares son cero y los primeros tres coeficientes son las componentes vectoriales de V X E + (9B/3f), y el último coeficiente es V • B. Así hemos establecido que
da=0
representa (9.16) y (9.18). De modo similar, la derivada exterior de j3 es di3
dH 3
dH2
dD ,
dy
dz
dt
dH 2
dH,
dD3
dy dD , - —
dx
dH , dz
a
A
dy
A
dt
) dx a
dy
a
dz,
dt
dD 2 + ——
a
dx
dt —
dH3
dD.
dx
dt
dz
dz
dD, + ——
dy
dy
dz
a
dx
a
dt
(10.97)
y escribiendo (10.95) explícitamente, y =
J i dy — p
dz
a
dx
a
dy
a
dt a
+]2 dz
a
dx
a
dt
+ J3 dx
a
dy
a
dt
(10.98)
dz.
Entonces sumando (10.97) y (10.98),
d/3+ y = -
'dH 3
dH2
dD,
\d y
dz
dt
'dH,
dH3
dD 2
,dz
dx
dt
fdH2
dH ,
dD3
\d x
dy
'dD, \d x
~
~JT
— JiJ dy
A
dz
A
dt
—J^j dz
a
dx
A
dt
—
A
dy
A
dt
dx
dD2 dD3 -f ------ + —p ) dx dz dy
0.
A
dy
A
dz
(10.99)
La ecuación (10.99) implica que todos los coeficientes que sean escalares son cero y los primeros tres coeficientes son las componentes vectoriales de V X H - (9D /9í) - J, y el último coeficiente es V • D —p. Así hemos establecido también que d/3 + y = 0 representa (9.17) y (9.18).
www.FreeLibros.me
Formas diferenciales
273
PROBLEMA 10.39 Si se da 7 por (10.95), mostrar que
d y = 0,
( 1 0 . 10 0 )
y que esto corresponde a la ecuación de continuidad V-J+^-=0.
[9.3]
dt
Por (10.92),
Solución:
y = - d ¡ 3.
( 10 .1 0 1 )
Así, tom ando la derivada exterior de (10.101) y usando (10.38), d y = —d(df3) = 0. Tomando la derivada exterior de (10.98), dy = d] j
dy
a
—d p
dx
a
/,
= —— dx
a
dx
----- — dt
dz
A a
dy
dy
A
dx
a
dt + d j 3 a dx
a
dy-
a
a
dt
a
dx
dz
a
dz
a
dx
dt + d j 2 a d z
a
a
dy 9
<9/2
dt ^— -— d y
a
dz
a
dy
dx
a
d /3
a
dt A------- d z dz
dy
a
a
dt
dz
a
dt dj<
d/,
(
dx
dy
dz
i9p\
—— H— ------1— -— + ——) dx
a
dtl
dy
a
dz
dt,
a
( 10 .10 2 )
y ( 1 0 . 10 0 ) implica claramente que dj.
d j,
dJ2
dp
1 7 + 1 7 + i r + 5 7 - 0i
0 0 ' 103)
esto es, V • J + dp/dt = 0.
10.9
Problemas suplementarios
PROBLEMA 10.40 Si a = x dx - z dy + y 2 dz, calcular a A /3. Respuesta: ( x 3 - y 3 - 2 z ) dx a dy a dz. PROBLEMA 10.41 Si (X = dx a dy + dy calcular a a ¡3. Respuesta: ( y - x ) dx a dy a d z a dw.
A
dz
¡3 =
—
dz
x 2 dy A dz
a
dw,
dr
a
+ 2
/3 = x
dz A dx
dx
a
dy
y dx A dy,
-
+
y dz
PROBLEMA 10.42 Verificar que
dx
a
dy
A
dz
a
dw
<9(x, v, z , w) = --------------------------
d(r, s, t, u)
ds
a
dt
a
du,
cuandox = x(r, s , t, u), y = y (r, s , t, u), z = z ( r , s , t ,u ) , and w = w (r, s , t, u). PROBLEMA 10.43 Calcular dco donde
(a)
cú =
x 2y dy
(b)
cú — 2
adz
xy dx +
— x z dx
a
dy,
x 2 dy,
(c) cú = 2 y z dy a dz + xydz Respuesta: (a) (2 x y - x ) dx a dy a dz,
a
(b) 0,
dx — x z dx a dy. (c) 0.
www.FreeLibros.me
A
dw,
274
Análisis vectorial
PROBLEMA 10.44 Demostrar las fórmulas de Leibnitz
d(fg) = i dg + g di,
d(ta>) = d f
a
co
+ i da>,
donde f y g son formas 0 y co es una forma 1 . PROBLEMA 10.45 Verificar la relación de invariancia (doj)* = c/(co*), cuando co =xy dx y T es la transformación x = u 2 + v, y = v. PROBLEMA 10.46 Una forma diferencial co se llama cerrada si dco = 0. Se llama exacta si co = da para alguna forma a. Mostrar que toda forma exacta es cerrada. PROBLEMA 10.47 Hallar una forma 1 co para la cual c/co = (x2 + y 2) d x A dy. Respuesta: cú = —x 2 y dx + x y 2 dy. PROBLEMA 10.48 Mostrar que una forma 1 co = 2 x y dx + x 2 dy + 2z dz es exacta. Respuesta: co = do., a = d ( x 2y + z 2). PROBLEMA 10.49 La forma 1
se define en el plano x y sin el origen (0,0). Mostrar que (a) co es cerrada, pero no exacta y (b) sobre cualquier curva cerrada C que no encierra (0, 0)
PROBLEMA 70.5*} Mostrar que el volumen de una región R en V3 es
ar donde &> = x dy
dz + y dz
A
a
dx + z dx
dy.
a
[Sugerencia: da> = 3 dx A dy A dx. ] Si M es una región del espacio 4-dimensional con coordenadas x, y, z, w, y dM es su frontera tridimensional, mostrar que PROBLEMA 10.51
f I J
/<9/t d /j <9f3 <9/4\ [ —— H------- + ----- + —— dx \o x dy d z dw / (/, dy
A
dz.
A
A
dw + f2 d z
dy
A
A
dw
dz
A
A
dw
dx + / 3 dw
A
dx
A
dy + / 4 dx
A
dy
A
dz),
uM donde las f . (x, y , z, w), / = 1 , __ , 4 son funciones diferenciables de x, y, z, w. PROBLEMA 10.52 Sea X una forma diferencial expresada como
A = A 1 dx + A 2 d y + A 3 d z — V dt, donde el vector potencial A = [A j , A 2, A 3 ] y el potencial escalar V son funciones de x, y, z, t. Si a está dada por (10.93), mostrar que d \ = a corresponde a las ecuaciones vectoriales B = V x A,
dt
www.FreeLibros.me
A
APENDICE
NOTACION MATRICIAL ' Y DETERMINANTES A l.
Notación matricial
Una matriz A es un arreglo de m X n números; i.e., a 11
a 12
a 2i
a2 2
ain 2n
9ml Un elemento general a., de A es el elemento de la z'-ésima fila y la /-ésima columna. La dimensión de una matriz de m filas y n columnas es m X n , que se lee “m por n”. Una matriz .4 m X n se escribe también com o/1 = [a ~ ] y la dimensión puede indicarse escribiendo A = [a.. ] m X n. Si n = m, la matriz se llama cuadrada. Si A = [a l..] ] y B = [b I..J ] son dos matrices m X n, entonces las matrices A y B son iguales si suscomponentes correspondientes son iguales, esto es, A = B •*—> a .. = b . $>\AyB son matrices m X n , entonces la suma de las matrices A y B es la matriz C = [ c .. ] cuyos elementos son la suma de los elementos de A y B\ esto es,C = A + B <—* c .. = a .. + b ... Obsérvese que la matriz suma C es también una matriz mXn. Si A es una matriz m X n y B es una matriz n X p , entonces el producto de las matrices A y i? es la matriz m X p n
C = [cij], donde c,¡ = ^
a ikbk¡.
k= 1
Obsérvese que el número de columnas de A debe ser igual al número de filas de B. Así, los productos matriciales A B y BA pueden definirse ambos solamente cuando A y B son matrices cuadradas y en general, A B i ^ B A . El producto XA del escalar \ y la matriz A es una matriz cuyo elemento general es X I]a ; esto es, ~KA = [Xa..]. IJ A 2.
Determinantes
Un número real conocido como el determinante de la matriz se asocia con toda matriz cuadrada que tiene números reales como elementos. La notación \A I o det A denota el determinante de la matriz A ; i.e., determinante de la matriz A = M I = det A. Para una matriz n X n A = [af.], \A \ es la suma de todos los términos posibles de la forma ( - I ) 9-3 ! / , a 2/2 • • • a rtjn >
donde no hay dos índices de columna que sean iguales y q es el número de transposiciones necesario para restaurar los índices
275
www.FreeLibros.me
Análisis vectorial
de columna a su orden natural. Así, si^4 =
au
12
}
a 21
entonces | A | = a i i a
22
a n
a 12
a 13
A = a 21
a 22
a 23
a 31
a 32
a 33
(—1
12
2 1
= a iia
22
—a i 2 a 2 ii
y ^
entonces ¡ A | = a n a 22a33 + ( ~ l ) a l l a 23a 32 + ( ~ l ) a 12a 2ia 33 +
^ ) 2 a 12a 23a 31 + ( ~ l ) 2 a 13a 21a 32 +
l ) a 13a 22a 31
= a n a 22a 33 ~ a l l a 23a 32 — a 12a 21a 33 + a 12a 23a 31 + a 13a 21a 32 ~ a 13a 22a 31-
Si A = [a..] es una matriz n X n, y A .. es la matriz (n - 1) X (n
—
1 ) formada quitando la /-ésima fila y la /-ésima columna
de^4, los cofactores a ., de A son a.. = (—1)!+/ \A..\. Por consiguiente, el desarrollo por la /-ésima fila da
A I= ^
a ikO-¡k>
k= 1
y el desarrollo por la /-ésima columna da
a k j a-k¡-
Así, si la matriz es
/c= 1 y a 12
a [2
A = a 21
a 22
a 23
_a 3 1
a 32
a 33
an
entonces su determinante desarrollado por la primera fila es = an
a 22
a 23
a 32
a 33
“ a 12
a 21
a 23
a 31
a 33
+ a 13
a 21
a 22
a 31
a 32
— a l l ( a 2 2 a 33 — a 23 a 32 ) ~ a 12( a 21 a 33 “ a 2 3 a 3 1 ) + a 13 ( a 2 1 a 32 ~ a 2 2 a 3 l ) = a l l a 2 2 a 33 — a l 1 a 23 a 32 - a 12 a 2 1 a 3 3 + a l 2 a 23 a 3 1 + a 13a 2 1 a 32 - a 1 3 a 2 2 a 31)
y el determinante desarrollado por la primera columna es \A \ = a n
a 22
a 23
a 32
a 33
-
a 2i
a i2
a 13
a 32
a 33
+ a 31
al2
a 13
a 22
a 23
= a n ( a 22a 33 - a 23a 32) ~ a 2 i ( a i 2 a 33 - a 1 3 a 32) + a 31( a 12a 23 - a i 3a 22) = a n a 22a 33 - a n a 23 a 32 - a 2 i a i 2 a 33 + a 2 i a i 3 a 32 + a 3 i a 1 2a2 3 - a 3 i a i3a22-
A3.
Propiedades de los determinantes
Propiedad 1: Si dos filas o columnas de una matriz cuadrada A se intercambian, el signo de I^4 I se cambia. Propiedad 2: Si una fila o columna de una matriz cuadrada ,4 se multiplica por una constante c, el valor del determinante se multiplica por c. Propiedad 3: Si se suma un múltiplo de una fila o columna a otra fila o columna, el valor del determinante no cambia. Propiedad 4: Si una fila o columna de una matriz cuad rad a^ es múltiplo de otra fila o columna, \A I es igual a cero.
276
www.FreeLibros.me
OPERACIONES DIFERENCIALES VECTORIALES EN COORDENADAS ORTOGONALES B1. diP _
=
7 —»+ dx
dip _ dip — J + — k, dy dz
B2.
dip dp
1
+
p
d\p
i
j
k
d
d
d
dx
dy
h
h
d i, d(2 d t3 V-f = — + T- + — dx dy dz
d 2xp
*
+
dz
+—
dx2
d 2xp
+ —
dy2
,
dz2 '
( ? ! ± _ d_í¿ lk \dz
dz )
dy
d 2xp
V 2i ¿ = —
Idx
dx )
dy
U
Coordenadas cilindricas
dip
v- f =
^ e* + ¿7k'
V X f=
B3.
APENDICE
Coordenadas rectangulares
V x f
V>A =
D D
dp
d
dz
íp
ftf>
/3
d
i
(p fp
p
1 )+ -
dftf,
p
dcp
dp
1 fd U
dfcp
p\dcf>
dz
1
d {3
v 2
6n +
dz
p
d
d 2ip
d 2xP
dip
1
dp
V2 d ^ 2 + J z 2 '
—
dp
d in p
d i,
1 / dítj}
d ít
dz
dp
p\dp
dcp
Coordenadas esféricas dtp
1
di p
dip
1
V ^ = — e r + - ~ - e 5 + --------- — (*£, dr r dO r sen 0 d
V-f =
sen 9
d
+T
dr
dd
1
w
r 2 sen 0 1
d
d(f>
sen 6
d
di p\
dr
dr I 1
L
V2 dr V
I d I d - — (r 2 i T) + -------- —— ( sen r dr r sen 0 dd
dicb
( sen 6 t e ) + t ——
dr)
4-
d
+ — dd d
f
dip\
sen 6 —
de dip
sen t) — i +
r
d Q
211
www.FreeLibros.me
I
d2ip
1
+ sen
0
dcf>2 _ d 2xp
1
A. 2 r2 sen2 6 ad
„0 „ 2
ñ
0 i e)
1 dU + -------- — r sen 6 dcp
Análisis vectorial
www.FreeLibros.me
RESUMEN DE A ^ CE RELACIONES VECTORIALES
C1.
Ecuaciones de álgebra vectorial
Ecuación No.
Ecuación
A B = B- A (1.31) (1.32) A • (B + C) = A • B + A • C A x B = -B x A (1.55) (1.56) Ax(B + C )= A x B + AxC (1.58}' Ax A=0 (1.72) A-BxC = AxB-C (1.76) A B x C = [ABC] [ABC] = [BCA] = [CAB] = - [ACB] = - [BAC] = -[CBA] . d.77) (1.83) A x (B x C) = (A ■C)B - (A • B)C (1.98) (A x B) x C = ( A- C) B - ( B- C) A (3.112) V (0 + 0 ) = v> + w (3.113) (3.124) V<£ = V0(u) = 0'(u)Vu (3.128) v-(f + 8) = v - f + v - g (3.131) 1 div (grad 0 ) = v • (V<¿) = (3.140) Vx(f+g) = V x f + Vxg (3.142) V X (V0) = 0 (3.143) V- ( V X f) = 0 ( f xV) - g = f-(Vxg) (3.154) (3.155) V-(^»f) = 0 V - f + f • (V<¿>) (3.156) V X (0f) = 0V X f + (v<¿) X f = 0 v X f - f X V<¿ (3.157) V-(fxg) = g. (V xf ) -f -( V x g) (3.158) V x ( f X g) = f ( V - g ) - g ( V - f ) + ( g - V ) f - ( f - V ) « (3.159) v(f • g) = r X (V x g) + g X (V X f) + (f • V)g + (g • V)f (3.163) rot (rot f) = V x (V x f) = V(V • f) - V2f (3.164) VJf = V ( V - f ) - V x ( V x f ) /0\
279
www.FreeLibros.me
Análisis vectorial
C2.
Ecuación No.
Ecuaciones de cálculo vectorial
Ecuación
(4 .5 8 ) R (4 .9 4 )
S
JjJ
V
0d V = £ J d S c / >
rt (4 .9 5 )
s
ffjvx ÍTv t-ds f-dT jj 9
idV=f f dSxt
R
(4 .1 0 4 )
S
X
s
(4 .1 1 5 )
c
dS x V
s
(4 .1 1 6 )
JJ
=^
0 Jr
c
WSxV)xf=^c/rxf
s
c
280
www.FreeLibros.me
INDICE Acción y reacción, 184 Aceleración, vector, 178 angular, 182-183 centrípeta, 181-183 componente normal, 181-182 componente tangencial, 181-182 en coordenadas cilindricas, 178-179 en coordenadas esféricas, 200 en coordenadas rectangulares, 178 tangencial, 183 Adición de vectores, 2, 19, 20 ley asociativa de la, 2, 3, 20 ley conmutativa de la, 2, 20 ley distributiva de la, 2, 20 Algebra de vectores, 1, 19 Ampere, ley de circuitos de, 238-239 Angular, aceleración, 182 momentum, 185, 186, 191-193, 199 rapidez, 182 velocidad, 182 Angulo entre dos vectores, 4, 26 sólido, 105 Anticonmutativa, ley del producto exterior (externo), 249, 250 del producto vectorial, 8, 22 Aplicación, 260 Apolonio, teorema de, 163 Arco, elemento de, 51 Arco, longitud de, 51, 138 en sistema coordenado curvilíneo ortogonal, 138 Area de la superficie, 84, 86 de un paralelogramo, 9 de un triángulo, 9 limitada por una curva cerrada simple, 117, 118, 120, 121 Asociativa, ley, 2-3, 20, 248
Barotrópido, flúido, 212 Base, 17, 30 autorrecíproca, 33, 38 derecha, 30 izquierda, 31 ortogonal, 32 ortonormal, 32 recíproca, 17, 30, 31
Base, vectores, 21, 134-135, 140-142 del sistema coordenado rectangular, 21 en un sistema de coordenadas cilindricas, 140 en un sistema de coordenadas curvilíneas, 135 en un sistema de coordenadas esféricas, 142 transformación de, 150, 155, 157 Bemoulli, ecuación de, 208, 210-211 Binormal, vector, 169 C
Cadena, regla de la, 44, 75 Cálculo diferencial vectorial, 41, 77 Cálculo integral vectorial, 78-133 Campo armónico, 240 conservativo, 187, 188, 235 eléctrico, 221 electrostático, 234 escalar, 58 gravitacional, 207 irrotacional, 122, 125, 235 magnético, 222 magnetostático, 235 solenoidal, 77, 122, 125, 125 Campo electromagnético, 220, 221 condiciones de frontera de un, 232-234 densidad de energía del, 230 densidad de momentum del, 246 energía en el, 230 funciones potenciales del, 246 Campos electrostáticos, 234 ecuación de Laplace para, 235 ecuación de Poisson para, 235 función de potencial escalar de, 235 ley de Coulomb para, 237 Campos monocromáticos, 240 (véase también campos armónicos) Campos sinusoidales con variación de tiempo, 240 (véase también campos armónicos) Campo vectorial, 58 conservativo, 187, 188, 235 irrotacional, 122, 125, 235 solenoidal, 77, 122, 125, 126 Carga densidad de, 220 distribución de, 221
281
www.FreeLibros.me
fuerza sobre, 201, 221 movimiento de la, 220, 221 principio de conservación de, 220 Cartesiano, sistema coordenado, 134 Cauchy-Schwarz, desigualdad de, 7 Central, fuerza, 184, 193 Centrífuga, fuerza, 198 Centrípeta, aceleración, 181, 183 Centro de masa, 190 Cerrada, curva, 78 Cerrada, forma diferencial, 274 Cerrada, superficie, 84 Circulación, 80, 97, 211 Coeficientes métricos (véase factores de escala) Combinación lineal, 4, 14, 17, 30 Complejo conjugado, vector, 247 Complejo de Poynting, vector, 247 Complejo en le espacio, vector, 240 Componente normal de aceleración, 181 Componentes de un vector, 5, 19, 47 transformación de, 151, 156 Conductividad, 222 térmica, 205 Conexa, región simplemente, 122, 126 Conjuntos recíprocos de vectores, 15, 136 Conmutativa, ley de la adición vectorial, 2, 20 del producto escalar, 5, 22 Conservación de carga, 220 de circulación, 212 de la energía, 189, 190 del momentum angular, 186, 193 del momentum lineal, 186, 191 Conservativa, fuerza, 187 Conservativos, campos, 187, 188, 235 condición necesaria y suficiente para, 114,188 Constituyentes, relaciones, 222 Continuidad ecuación de, 203, 220 de una función vectorial, 41 Coordenadas cartesianas, 148 cilindricas, 139, 149 cilindricas circulares, 149 cilindricas elípticas, 159 cilindricas parabólicas, 159 curvas, 134 curvilíneas. 134
Análisis vectorial
esféricas, 141, 152 esferoidales achatadas, 159 esferoidales alargadas, 159 neutras, 139, 141 ortogonales curvilíneas, 137 parabólicas, 160 rectangulares, 19, 148 superficies, 134, 148, 149, 152 transformación, 134 Coordenadas cilindricas elemento de volumen en, 149 factores de escala en, 140, 149 jacobiano de una transformación en, 139 longitud de arco en, 149 parabólicas, 159 rotacional, divergencia, gradiente, laplaciano en, 149 superficies coordenadas, 149 Coordenadas curvilíneas, 54, 134-137 aceleración en, 178-179, 200 elemento de volumen en, 138 factores de escala en, 138 longitud de arco en, 138 ortogonales, 137 rotacional, divergencia, gradiente, laplaciano en, 143-144 Coordenadas curvilíneas ortogonales, 134-160 rotacional, divergencia, gradiente, laplaciano en, 143-148 Coordenadas esféricas campos estáticos, 234-240 elemento de volumen en, 153 factores de escala en, 141, 153 jacobiano de una transformación en, 141 longitud de arco en, 153 rotacional, divergencia, gradiente, laplaciano en, 153 superficies coordenadas, 152 Coordenadas ortogonales cilindricas, 139, 149 cilindricas elípticas, 159 cilindricas parabólicas, 159 esféricas, 141, 152 esferoidales achatadas, 159 esferoidales alargadas, 159 especiales, 148-160 parabólicas, 160 rectangulares, 19, 148 Coplanares, vectores, 10 condición para, 10, 18, 28 no, 4 Coriolis, fuerza de, 198 Correspondencia uno a uno, 19 Corriente, 220 densidad de, 220 desplazamiento de, 225 estacionaria, 238 Cosenos directores de un vector, 27 ley de los, 7 Coulomb, indicador de, 245 Coulomb, ley de, 237 Cruz, producto, 8, 22 (véase también vectorial, producto) como un caso especial de product externo de formas, 251 en forma de determinante, 22
ley anticonmutativa del, 8, 22 ley distributiva del, 8 ley seudo-simétrica del, 8 Cuerpo, coordenadas del, 195 Cuerpos rígidos, 194-200 Curva borde de una, 266 curvatura de una, 169 ecuaciones paramétricas de una, 49, 51, 169 ecuaciones vectoriales para una, 48, 52 en el espacio, 48, 168 fórmulas de Frenet-Serret para una, 170 integral lineal a lo largo de una, 178 longitud de arco de una, 51, 138 paramétrica, 54 plana, 171 principal, 169 punto no singular de una, 49 punto singular de una, 49 tangente a una, 49 vector tangente a una, 49, 169 torsión de una, 169 vector binormal de una, 169 vector de Darboux de una, 170 vector normal de una binormal, 169 una principal, 169 Curvatura, 181 radio de, 181 Ch
Charles, teorema de, 198-199
O Darboux, vector de, 170 V (véase del) Va (véase laplaciano, operador) del (v), 58, 63 (véase también gradiente, divergencia y rotacional) fórmulas en las cuales interviene, 68-74 representación integral de, 100, 158 Delta de Kronecker, 15, 24, 31, 202 Densidad de carga, 220 de carga de superficie, 233 de corriente, 220 de corriente en la superficie, 234 de flujo de masa, 204 de flujo eléctrico, 222, 225 de flujo magnético, 221 de la fuerza de Lorentz, 246 del flúido, 203 de mornentum, 246 Densidad de flujo de masa, 204 eléctrico, 222 magnético, 221 Dependencia funcional, 77 lineal, 4, 18, 35 Derivada de una función vectorial, 43 parcial, 47 direccional, 58 exterior, 253
282
www.FreeLibros.me
normal, 107 Desigualdad del triángulo, 8 Desplazamiento corriente de, 25 eléctrico, 222, 225 infinitesimal, 195 vector de, 178 Determinante cofactor de un, 276 expresado como rotacional, 65, 147, 149, 153 expresado como triple producto escalar, 24 jacobiano (véase jacobiano) propiedades de un, 276 Diádico, 70 Diferencia de dos vectores, 2, 20 Diferenciable función escalar, 44, 253 función vectorial, 43, 47 transformación, 260 Diferenciación de formas, 253 de vectores, 43 parcial, 47 reglas para, 44 Diferencial elemento (diferencial) del área de una superficie, 56, 84 vector de desplazamiento, 78 Dinámica de los fluidos, 203 de una partícula, 178 de un sistema de partículas, 190 Dipolo campo eléctrico de, 237-238 eléctrico, 237 momento de, 237 potencial de, 237-238 puntual, 237 Dirección de integración, 79 de una curva, 49 de un vector, 1 en sentido contrario al de la manecillas del reloj, 79 en el sentido de las manecillas del reloj, 79 negativa, 79 positiva, 79 Direccional, derivada, 58 Distancia de un punto a una recta, 165 de un punto a un plano, 166, 167 entre dos rectas, 165 Distancia más corta entre dos rectas, 165 Distributiva, ley (véase ley distributiva) Divergencia de la función vectorial, 63, 255 definición alterna, 95, 98 interpretación física de, 96 del gradiente, 64 del rotacional, 67 en coordenadas cilindricas, 149 en coordenadas curvilíneas ortogonales, 143-144 en coordenadas esféricas, 153 en coordenadas rectangulares, 63, 148 teorema de (véase teorema de divergencia)
Indice
Ecua«ión de Bernoulli, 208, 210 de continuidad de dinámica de los fluidos, 203 en teoría electromagnética, 220 de equilibrio, 207 de Euler, 205-206 de fuerza de Lorentz, 201 de Helmholtz, 242 de Laplace, 65, 235, 245 de Maxwell, 222 de movimiento para un flúido, 205, 208 de onda, 228 de una recta, 164 de un plano, 167 Ecuaciones no homogéneas de onda, 228 Ecuación no paramétrica de una recta, 50 Ecuación paramétrica, 49 de una curva, 49 de una esfera, 56, 57 de una recta, 50 de una superficie, 53 de un círculo, 49 Ecuación vectorial de una recta, 164 de un plano, 166-168 Eje fijo, 182 instantáneo de rotación, 195 Eléctrica carga, 220 corriente, 220 densidad de carga, 220 densidad de corriente, 220 Eléctrico campo (eléctrico) de dipolo, 237-238 campo (eléctrico) de una carga puntual, 236 densidad de flujo, 222 desplazamiento, 222, 225 dipolo, 237 flujo, 224, 225 intensidad de campo, 221 vector de campo, 201-221 Electromagnético, campo, 220, 221 Electrostática, 234 Elemento de arco, 51 de una matriz, 275 diferencial del área de una superficie, 84 Elemento de volumen, 94 en coordenadas curvilíneas ortogonales, 138 Energía conservación de la, 189 densidad de, 230 en el campo electromagnético, 230-232 mecánica, 190 potencial, 187 teorema de energía mínima de Kelvin, 217 Energía cinética de una partícula, 187 de un cuerpo rígido, 198 de un flúido, 216 de un sistema, 194 rotacional, 198 translacional, 198
Escalar, 1 campo, 58 función, 41 integral de línea, 80 potencial, 122, 206 producto, 4, 22, 250 (véase también producto punto) triple producto, 10, 27 Espacio coordenadas del, 195 curvas en el, 48, 168 (véase también curva) generan el, 17 m-dimensional, 260-264 n-dimensional, 268 p-dimensional, 268 (p-l)-dimensional, 268 tridimensional, 14, 248, 254 Estática de fluidos, 207 Euler, ecuación de, 205-206 Exterior derivada, 253 diferenciación, 253 forma diferencial, 248 producto, 249 ley anticonmutativa del, 249, 250 F
Factores de escalar, 138 (véase también coeficientes métricos) en coordenadas cilindricas, 140, 143, 149 en coordenadas cilindricas elípticas, 159 en coordenadas cilindricas parabólicas, 159 en coordenadas curvilíneas ortogonales, 138 en coordenadas esféricas, 141-143, 153 en coordenadas esféricas alargadas, 159 en coordenadas parabólicas, 160 en coordenadas rectangulares, 148 en esferoidales achatadas, 159 Faraday, ley de inducción de, 223 Fase del vector complejo, 240 Flúido barotrópido, 212 ideal, 205 incompresible, 204 perfecto, 205 Fluidos dinámica de los, 203 estática de, 207 mecánica de los 203-249 Flujo, 86, % de energía, 232 de masa fluida, 203 eléctrico, 222 ¡■rotacional, 209 magnético, 221 potencial, 209 rotacional, 209 uniforme, 208 vórtice, 209. 211 Forma p (véase formas diferenciales) Formas diferenciales, 248-274 adición y multiplicación de, 248-249 cerradas, 274
283
www.FreeLibros.me
de las ecuaciones de campo de Maxwell, 270-273 diferenciación de, 253 exactas, 274 invariancia de, 260 Frenet-Serret, fórmulas de, 170 Frontera condiciones de, 207, 232 curva, 122 de superficie, 266 superficies de, 101 Fuente, 96 Fuerza central, 184, 193 centrífuga, 198 centrípeta, 181 conservativa, 187, 206 de Coriolis, 198 de Lorentz, 201, 221 efectiva, 198 electromotriz, 224 momento de, 185 resultante, 205 sobre una partícula móvil, 201 sobre un dipolo eléctrico, 246 Fuerzas externas, 190, 193, 205 internas, 190, 193 Función armónica, 65 Funciones escalares, 41, 47 indicadoras, 228 uniformes, 122-123 vectoriales, 41, 47 Funciones vectoriales continuas, 41 continuidad de, 41 curva de, 48, 52 de más de una variable, 47 derivada de, 43 derivada parcial de, 47 de una variable, 41 diferenciables, 43 diferenciación de, 43 divergencia de, 63 lím ite de, 41 rotacional de, 65 suaves, 49 valor de, 41
Gauss, ley de para el campo eléctrico, 225 para el campo magnético, 225 Gauss, teorema de divergencia (véase teorema de divergencia) Geometría, 161-177 analítica, 164-168 del espacio, 161-164 diferencial, 168-175 plana, 161-164 Gradiente de una función escalar, 58, 255 definición alterna de, 59, 58 de una función vectorial, 70 (véase diádico) en coordenadas cilindricas, 149
Análisis vectorial
en coordenadas curvilíneas ortogonales, 143-144 en coordenadas esféricas, 153 en coordenadas rectangulares, 148 teorema de (véase teorema de gradiente) Gram-Schmidt, proceso de ortogonalización, 34-35 Gravitacional, campo, 207 Green, primer teorema o identidad de, 106 segundo teorema o identidad de, 106 tercer teorema o identidad de, 107 Green, teorema de, para el espacio, 100 (véase teorema de divergencia) Green, teorema de, para un plano, 116-117 como case especial del teorema de divergencia, 116-117 como caso especial del teorema de Stokes, 116
H
Hélice, 174 Helmholtz, ecuaciones de, 219, 242 Helmholtz, primer teorema sobre vorticidad de, 214 segundo teorema sobre vorticidad de, 214 Helmholtz, teorema de, 133 Hertz, vector, 229 Hooke, ley de, 189
Ideal, flúido, 205 Igualdad de matrices, 275 de vectores, 1 , 19 Impedancia característica de un medio, 243 Incompresible, flúido, 204 Independencia de la trayectoria, 122-124 Indicadora función, 228 invariancia, 227 transformación, 227 Inercia, momento de, 202 producto de, 202 Infinitesimal, desplazamiento, 195 rotación, 195 Integración de formas, 263-270 de una función vectorial, 78-133 Integral de formas diferenciales, 266 de línea, 78, 268 de superficie, 84, 268 de volumen, 94, 268 Integral de línea, 78, 80, 268 circulación en términos de, 80 independencia de la trayectoria, 122-124 trabajo expresado en términos de, 186 Intensidad de campo eléctrico, 221 de campo magnético, 222 Intensidad de un tubo de vórtice, 211 Intervalo cerrado, 263 Invariancia de derivadas externas, 260
del operador d, 260 Irrotacional .flujo, 209 vector, 77 lrrotacionales campos, 122, 125, 135 condiciones para, 125 J
Jacobi, identidad de, 15 Jacobiano, 77, 132, 134, 139, 141, 158, 251-253
K
Kelvin, teorema de energía mínima de, 217 Kelvin, teorema de, 212, 217 Kepler, segunda ley de, 179-180, 184-185 Kronecker, delta de, 15, 24, 31, 130, 202 L
Lagrange, identidad de, 9 ampliada, 14 Laplace, ecuación de, 65, 210, 235, 245 Laplaciano, operador V 2, 64 en coordenadas cilindricas, 149 en coordenadas curvilíneas ortogonales, 144, en coordenadas esféricas, 153 en coordenadas rectangulares, 149 Leibniz, fórmulas de, 274 Ley cinemática de un cuerpo rígido, 196 Ley conmutativa (véase conmutativa) Ley de conservación de circulación, 212 (véase teorema de Kelvin) Ley de inducción, 223 Ley de cuadrado inverso para fuerza conservativa, 188
M Magnético campo (magnético) de alambre infinito que lleva corriente, 239 densidad de flujo, 221 flujo, 223, 224 inducción, 201, 221 intensidad de campo, 222 Magnetomotriz, fuerza, 224 Magnetostático, campo, 235 Magnitud de un vector, 1, 20 Masa, 184, 203 densidad del flujo de, 204 Matrices (véase también matriz) igualdad de, 275 producto de, 275 suma de, 275 Matriz (véase también matrices) cuadrada, 275 determinante de una, 275 de transformación de vectores base y componentes de un vector, 155-157
Límite de una función, 42 de una función vectorial, 41 vector, 41 Línea de flujo, 208 Linealmente dependientes, vectores, 4, 14, 35 independientes, vectores, 4 Longitud de un vector, 20 (véase también magnitud de un vector)
dimensión de una, 275 elemento de una, 275 notación de, 155-157 Maxwell, ecuaciones de, 222 formas diferenciales de las, 270-271 para campo armónico, 240 Mecánica, 178-202 Medio conductividad del, 222 homogéneo, 222, 240 impedancia característica de un, 243 isotrópico, 222, 240 lineal, 222, 240 permeabilidad del, 222 permitividad del, 222 Momento, 185 de dipolo eléctrico, 237 de fuerza, 185 de inercia, 202 sobre un dipolo eléctrico, 246 Momentum angular, 185 densidad de, 246 lineal, 184 Monocromáticos, campos, 240 (véase también campos armónicos) Monomios, 257 Movimiento de una partícula, 184 de un cuerpo rígido, 194 de un sistema, 190 ecuación del (para un flúido), 205 ley de Newton para el, 184 planetario, 179-180
Lorentz lema de, 241 teorema de reciprocidad, 241 Lorentz, condición de en el campo armónico, 241 para potenciales, 228, 230 Lorentz, fuerza de, 201, 221 densidad de la, 246 ecuación de, 201
Nabla (véase del) Negativo de un vector, 2, 19 Newton, ley del movimiento de, 184, 205 Normal al vector tangente de una curva, 53 a una superficie, 54 principal, 169, 181
Ley do los cosenos, 7 Ley do los senos, 9 Ley distributiva de la adición de vectores, 2, 20 del producto escalar, 5, 23 del producto vectorial, 8, 23 para la adición y la multiplicación de formas diferenciales, 248, 249
284
www.FreeLibros.me
N
Indice
vector, 54, 84 hacia afuera, 220
Onda ecuaciones de, 228 ecuaciones no homogéneas de, 228 monocromática, 243 (véase también onda plana) número de, 242 plana, 243 transversa, 243 vector de, 242 Operador, d, 260 del, 58, 63, 68 (véase también del) diferencial 7 (véase también del) divergencia, 63 gradiente, 59 rotacional, 65 rotacional, 65 Ortogonales bases, 32 vectores, 6 Ortonormales bases, 32 vectores, 33 P
Paralelismo de vectores, condición para, 1. 2, 9 Partícula aceleración angular de, 182 aceleración de, 178, 181 desplazamiento de, 178 energía cinética de una, 187 energía mecánica de una, 190 leyes del movimiento de Newton para, 184 momentum angular de una, 185 momentum de una, 184 potencial (energía) de una, 187 rapidez de, 180 rotación de, 195 velocidad angular de una, 182 velocidad de una, 178, 181 Permeabilidad, 222 Permitividad, 222 Perpendicularidad de vectores, condición para, 6 Pitágoras, teorema de, 17 Pontearé, lema de, 258, 260 Poisson, ecuación de, 235 Potencial, 187 de una carga puntual, 236 de un campo electrostático, 235-236 de un campo magnetostático, 239 de un dipolo, 237-238 de velocidad, 209-210 energía, 187 escalar, 122, 206 flujo (véase flujo irrotacional) funciones del campo electromagnético, 226, 245 escalares, 122, 206 vectoriales, 126
gravitacional, 189 vectorial, 126 Poynting, teorema de, 231 vector complejo de, 247 vector de, 230 Presión de un fluido, 203 Principio de conservación de carga, 220 de circulación, 212 de la energía, 189-190 del momentum angular, 186, 193 del momentum lineal, 191 Producto cruz (véase cruz, producto) cuña, 249, 260 de formas, 248-249 de inercia, 202 de matrices, 275 de un escalar y una matriz, 275 de un vector por un escalar, 1 escalar (véase escalar, producto) exterior (externo), 8, 22, 249 interno, 4, 22 vectorial (véase vectorial, producto) Producto punto, 4, 22 (véase también producto escalar) como caso especial del producto exterior de formas, 251 ley conmutativa para el, 5, 22 Propiedad de linealidad de las derivadas exteriores, 256 Proyección de un vector, 4 Punto de estancamiento, 211 final de un vector, 1 inicial de un vector, 1 Punto no singular de una curva, 49 de una superficie, 54 Punto singular de una curva, 49 de una superficie, 54
Radio de curvatura, 169, 181, 200 de torsión, 169 Rapidez, 180 angular, 182 Rectangulares, coordenadas, 19, 148 Rectángulo cerrado, 265 Región conexa, 122 simplemente conexa, 122, 125 sin corriente, 240 Regla del término medio para el triple producto vectorial, 14 Relaciones entre vectores base rectangulares y cilindricos, 150 entre vectores base rectangulares y esféricos, 155 Relaciones constituyentes, 222 Representación de vectores algebraica, 19 analítica, 19
285
www.FreeLibros.me
geométrica, 1 Representación integral de V , 100, 111 Resultante de vectores, 2, 19 Rotación eje instantáneo de, 195 infinitesimal, 195 Rotacional definición alterna de, 95, 98 del gradiente, 66 del rotacional de una función vectorial, 72 de una función vectorial, 65, 66, 255 en coordenadas cilindricas, 149 en coordenadas curvilíneas ortogonales, 143, 144 en coordenadas esféricas, 153 en coordenadas rectangulares, 149 interpretación física de, 197 teorema del, 109
S Schwarz, Cauchy, desigualdad de, 7 Segmento dirigido, 1, 19 Sinusoidales, campos con variación de tiempo, 240 (véase también campos armónicos) Sistema coordenado rectangular, 19 Sistemas coordenados especiales, 148-160 Sistemas coordenados fijo y del cuerpo, 195 Sistemas de partículas, 190-194 Solenoidal campo, 122, 125, 126 condiciones para, 126 vector, 77 Stokes, teorema de, 112, 266, 268 clásico, 268 demostración del, 112-113 general, 268 en el espacio n-dimensional, 268 interpretación física del, 113 para coordenadas rectangulares, 115 teorema de Green para un plano como caso especial del, 116 especial del, 116 Suave curva, 49, 263 función vectorial, 49 superficie, 264 Suma de formas diferenciales, 248 de matrices, 275 de vectores, 2, 19 Sumidero, 96 Superficie, 53 área de una, 56 elemento diferencial del, 56 cerrada, 84, 113 curvas paramétricas de una, 54 de fase constante, 243 densidad de carga de una, 232 densidad de corriente en una, 234 ecuación paramétrica de una, 53 frontera, 101 frontera de una, 266 integral de, 84, 268
punto no singular de una, 54 punto singular de una, 54 representación vectorial de una suave, 264 vector normal a una, 54 Sustracción de vectores, 2
Teorema de divergencia, 100, 109, 269 como caso especial del teorema general de Stokes en formas diferenciales, 268-269 en formas diferenciales, 268-269 demostración del, 100, 101 expresado en coordenadas rectangulares, 103 interpretación física del, 101 teorema de Green para un plano como caso especial del, 116-117 Teorema de gradiente, 109, 205 demostración del, 110 Teorema de Helmholtz, 133 Teorema de reciprocidad de Lorentz, 241 Teorema general de Stokes, 268 Teoremas de Kelvin, 212, 217 Teoría electromagnética, 220-247 Tema ordenada, 19 Tema positiva de vectores, 10 Tiempo de relajamiento, 226 Torsión, 169 radio de, 169 Totaí energía mecánica, 190 momentum angula;, 199 Trabajo, 186, 193 - ■ Transformación, 134, 2 6 0 ' de integral de superficie » integral de linca, 118 de integra] de volumen a integral de superficie, 109 ■ diferenciable, 260 entre componentes vectoriales rectangulares y cilindricas, 151 entre componentes vectoriales rectangulares y esféricas, 156 entre coordenadas, 134 entre vectores base cilindricos y esféricos, 157 entre vectores base rectangulares y cilindricos, 150 entre vectores base rectangulares y esféricos, 150 integral, 109, 118 inversa, 260 jacobiano de una, 134, 251 uniforme, 134 Trayectoria, 208 Triedro, 169 Triple producto escalar, 10 vectorial, 12 U
128\
Unicidad, teorema efe* Uniforme función, 122-1^3 transform acióií, 134
ción de, 2
vector (unitario) tangente, 52 Uno a uno correspondencia, 19 transformación, 260
Variables, cambio de, 251 Vector aceleración, 178 binormal, 169 cero, 1, 19 complejo, 240 conjugado, 247 de Poynting, 247 en el espacio, 240 componente de un, 5, 20 cosenos directores de un, 27 de Darboux, 170 de potencial, 126 de Poynting, 230 de propagación, 242 desplazamiento, 178 fijo, 1 Hertz, 229 irrotacional, 77 libre, 1 límite, 41 localizado, 1 longitud (magnitud) de un, 20 no localizado, 1 normal, 54, 84 normal principal, 169 operador V (véase del) proyección de un, 5 representación algebraica de un, representación geométrica de u r 39-169 solenoidal, 77 tangente, 49, 135 unitario, 2, 20 velocidad, 178 Vector de posición, 21, 48, 78 identidades en las cuales interviene, 48, 52, 53, 54, 58, 6 1 ,6 2 , 64, 6 6 ,6 9 , 73, 76, 77, 8 3 ,9 5 ,1 0 1 ,1 1 2 -1 1 3 120-121
en coordenadas cilindricas, 139, 178 en coordenadas esféricas, 141, 200 Vectores adición de, 2 álgebra de, 1 ángulo entre, 4, 26 base, 21, 134-135, 140-142 combinación lineal de, 4, 14, 17 conjuntos recíprocos de, 15, 136 coplanares, 10 diferencia de, 2, 20 igualdad de, 1, 19 * linealmente independientes, 4 ortogonales, 6 resultante de, 2, 19 suma de, 2, 19
www.FreeLibros.me
, 2, 19, 20 , 19 campo, 58 producto, 8, 22 triple producto, 12 Velocidad, 178 angular, 182 de propagación, 243 potencial, 209 relativa, 201 vector, 178 en coordenadas cilindricas, 178 en coordenadas esféricas, 200 1- en coordenadas rectangulares, 178 Viscosidad, 205 Volumen, integrales de, 94, 268 Vórtice flujo, 209, 211 (véase también flujo rotacional) línea de, 211 tubo de, 211 vector, 209 Vorticidad, 209
mm
ADDISON-W ESLEY IBEROAM ERICANA Apartado Postal 22-456. M éxico 14060, D.F., M éxico Apartado Aéreo 29696. Bogotá, Colom bia Casilla 70060. Santiago 7, Chile Apartado Postal 29853. Rio Piedras, Puerto Rico 00929 M artín ae los Heros 59, Bis, 4. 28008, M adrid, España 7 Jacob W sy. F «ad ing , M as s a c h u s e tts , E.U.A.
39-169 *
3 9 -
1 6 9 *
978020102943790000
www.FreeLibros.me