JOSÉ MILCIADEZ DÍAZ, RAFAEL CASTILLO, FERNANDO VEGA

PONTIFICIA UNIVERSIDAD JAVERIANA, DEPARTAMENTO DE FÍSICA

INTRODUCCION AL ANALISIS VECTORIAL

(Parte 3)

Donde encuentre el símbolo …..

! conduce a una animación ejemplo.
" le lleva a una solución detallada.

Introducción | Parte 1 | Parte 2 | Parte 4

13.- Derivada de un vector con respecto a una variable escalar.

Consideremos que t es una variable escalar y asumamos que hay una función escalar u = u (t ) y dos
r r r r
funciones vectoriales A = A(t ) y B = B (t ) entonces, de acuerdo premisas del cálculo diferencial,
tenemos que:
r r

dt
(
A+ B = )
d r r dA dB
+
dt dt
r
d r
dt
( )
uA = u
dA du r
+
dt dt
A

r r

dt
(
d r r
)r dB dA r
A• B = A• +
dt dt
•B (13-1)

r r

dt
(
d r r
)r dB dA r
A× B = A× +
dt dt
×B

Sí las direcciones de los ejes coordenados no cambian y las componentes ortogonales de un vector
son funciones de la variable escalar t, digamos que Ax = Ax (t ) , Ay = Ay (t ) y Az = Az (t ) , entonces

r
dA dAx ˆ dAy ˆ dAz ˆ
= i+ j+ k (13-2)
dt dt dt dt

El diferencial de una función vectorial con respecto a una variable escalar es igual a la derivada de la
función por el diferencial de la variable escalar:
r
r  dA 
dA =  dt = iˆdAx + ˆjdAy + kˆdAz (13-3)
 dt 
INTRODUCCION AL ANALISIS VECTORIAL
(JM Díaz, R. Castillo, F. Vega).

r
Ejemplo Ilustrativo 13.1. a) El vector A = Ax (t )iˆ + Ay (t ) ˆj + Az (t )kˆ es función del parámetro t. La
dependencia de las componentes escalares es: Ax (t ) = sen(2t ) ; Ay (t ) = cos(2t ) y Az (t ) = 3t 2 − 5t .
Considerando los ángulos en radianes, se requiere: a) las dos primeras derivadas del vector; b)
evaluar el vector y sus dos primeras derivadas cuando t =1; c) hallar el ángulo entre el vector y su
primera derivada para cualquier valor de t, para t =1 y para t =2.

Solución.
a) Las dos primeras derivadas con respecto al parámetro t son:
r r
d2A
= iˆ(2 cos 2t ) − ˆj (2 sen 2t ) + kˆ(6t − 5) ; = −iˆ(4sen 2t ) − ˆj (4 cos 2t ) + kˆ(6)
dA
dt dt 2

b) Cuando t =1, los valores obtenidos

r serán: r
r dA(1) d2A
A(1) ≅ 0.91iˆ − 0.42 ˆj − 2kˆ ; ≅ −0.83iˆ − 1.82 ˆj + kˆ ; ≅ −3.64iˆ + 1.66 ˆj + 6kˆ
dt dt 2

c) Utilizando el producto escalar, el ángulo entre el vector y su primera derivada para cualquier valor
de t, será

cosθ =
(3t 2
− 5t )(6t − 5)

[1 + (3t ] [4 + (6t − 5) ]
1

− 5t )
1
2 2 2 2 2

para t = 1:
 − 2
θ = cos − 1   ≅ 1.982 radianes.
 5 
Para t =2:
 14 
θ = cos −1   = 0.54 radianes.
 265 

35
INTRODUCCION AL ANALISIS VECTORIAL
(JM Díaz, R. Castillo, F. Vega).

Ejemplo Ilustrativo 13.2. Desplazamiento y velocidad de una partícula.

La figura muestra la trayectoria L de una partícula que se mueve en el espacio; los puntos A y B son
r r
los ocupados por la partícula en los instantes tA y tB y cuyos vectores de posición son rA y rB ,
respectivamente. Durante el intervalo de tiempo ∆t = t B − t A la partícula siguió la trayectoria L y al
r r r
final del intervalo tuvo un cambio de posición ∆r = rB − rA que se denomina desplazamiento de la
partícula y que ocurrió en el intervalo de tiempo ∆t ; el desplazamiento es un vector cuya dirección
r
r r ∆r
es la de la secante AB y cuya magnitud es rB − rA . El cociente representa el desplazamiento
∆t
r
r ∆r
por unidad de tiempo de la partícula y se denomina velocidad media, v m = . La velocidad media
∆t
es un vector que tiene la dirección del desplazamiento y su magnitud (en metros/segundo) es igual a
la distancia en línea recta entre los puntos A y B dividida por el tiempo gastado por la partícula en ir
de A a B.

A D

C
→
→ ∆r
rA ∆t ∆r
z B
→ L
rB Tr a ia
O yector

O y
x

Sí el intervalo de tiempo considerado es menor, digamos la mitad, la partícula probablemente solo

llega hasta el punto C; el desplazamiento cambia en dirección y magnitud; la velocidad media
tendría la dirección de la secante AC. Al disminuir mucho más el intervalo de tiempo el
desplazamiento de la partícula se reduce y solo alcanzaría una posición muy cercana al punto A,
digamos la posición D en la figura; en esta situación observamos que el desplazamiento tiene la
dirección de la cuerda AD y que esa dirección es más cercana a la dirección de la tangente de la
trayectoria L en el punto A. En el límite, cuando ∆t → 0 , la dirección de la cuerda tiende a la
dirección de la tangente de la trayectoria en el punto A y la velocidad media tiende a la velocidad
r r
r r  ∆r  dr
instantánea de la partícula (v ) . Así, la velocidad instantánea será: v = lim ∆t →0   ≡ y es un
 ∆t  dt
r
dr
vector tangente a la trayectoria de la partícula cuya magnitud es v = .
dt

Como ilustración numérica, consideremos el movimiento del pasador móvil de una manivela en
rotación y cuya posición está definida por x = 2 cos πt metros; y = senπt metros. Hallar a) el vector
de posición y la velocidad para t = 1 y t =3/2 segundos; b) el desplazamiento y la velocidad media en
ese intervalo de tiempo.

36
INTRODUCCION AL ANALISIS VECTORIAL
(JM Díaz, R. Castillo, F. Vega).

Solución.
El vector de posición y la velocidad de la partícula para cualquier instante de tiempo son:
r
r (t ) = (2 cos πt )iˆ + (2senπt ) ˆj ; v (t ) = = (− 2π senπt )iˆ + (2π cos πt ) ˆj
r r dr
dt
r r
para t = 1 segundo r (1) = 2iˆ (m); v (1) = −2πˆj (m/s)
r r
para t = 3/2 segundos r (3 / 2) = −2 ˆj (m) v (3 / 2) = 2πiˆ (m/s)

r r r
b)- ∆r = r (3 / 2 ) − r (1) = 2iˆ − 2 ˆj (metros);
r
vm =
2iˆ − 2 ˆj
1/ 2
( )
= 4 iˆ − ˆj (m/s)

( )
Ejemplo Ilustrativo 13.3. Los vectores unitarios ρˆ , φˆ, ẑ del sistema de coordenadas cilíndricas
están relacionados (ecuación 9.3) con los vectores unitarios del sistema cartesiano iˆ, ˆj , kˆ por las ( )
dρˆ & ˆ
expresiones: ρˆ = iˆ cos φ + ˆjsenφ ; φˆ = −iˆsenφ + ˆj cos φ ; y zˆ = kˆ . Demostrar que = φφ y que
dt
dφˆ
= −φ&ρˆ ; donde el punto sobre la variable indica la derivada de la variable con respecto al tiempo.
dt

Solución. Por derivación directa:

dρˆ dφ dφ dφ & ˆ
(a) = iˆ(− senφ ) + ˆj (cos φ ) = (− iˆsenφ + ˆj cos φ ) = φφ
dt dt dt dt
dφˆ dφ dφ dφ
(b) = −iˆ(cos φ ) + ˆj (− senφ ) = −(iˆ cos φ + ˆjsenφ ) = −φ&ρˆ
dt dt dt dt

37
INTRODUCCION AL ANALISIS VECTORIAL
(JM Díaz, R. Castillo, F. Vega).

Ejemplo Ilustrativo 13.4. Dado el vector de posición de un punto en coordenadas cartesianas como
r
función del tiempo, r (t ) = iˆx(t ) + ˆjy (t ) + kˆz (t ) , Hallar la velocidad y la aceleración a) en coordenadas
cartesianas, b) en coordenadas cilíndricas. c) considere que x(t ) = t 2 ; y (t ) = 2t ; z (t ) = t 3 , determinar el
vector de posición, la velocidad y la aceleración en coordenadas cilíndricas cuando t=1.

Solución.
La longitud la tomamos en metros y el tiempo en segundos.

a) Como la velocidad es la variación de la posición con respecto al tiempo, entonces:

r
r
v (t ) =
dr
= x&iˆ + y&ˆj + z&kˆ = v x iˆ + v y ˆj + v z kˆ
dt
r
r
a (t ) =
dv
= &x&iˆ + &y&ˆj + &z&kˆ = a x iˆ + a y ˆj + a z kˆ
dt
r
b) El vector de posición en coordenadas cilíndricas es r (t ) = ρρˆ + zzˆ , hallado por inspección directa
o utilizando las ecuaciones (9.1). Las expresiones para la velocidad y la aceleración serán:

r
v (t ) = (ρρˆ + zzˆ ) =
d
dt
dρ
dt
ρˆ + ρ
dρˆ dz
+ zˆ = ρ&ρˆ + ρ φ&φˆ + z&zˆ = ρ&ρˆ + ρφ&φˆ + z&zˆ
dt dt
( )
r
a (t ) =
d
dt
( )
ρ&ρˆ + ρφ&φˆ + z&zˆ = ρ&&ρˆ + ρ&
dρˆ
dt
+ ρ&φ&φˆ + ρφ&&φˆ + ρφ&
dφˆ
dt
+ &z&zˆ

r
( ) ( )
a (t ) = ρ&&ρˆ + ρ& φ&φˆ + ρ&φ&φˆ + ρφ&&φˆ + ρφ& − φ&ρˆ + &z&zˆ ( ) (
= ρ&& − ρφ& 2 ρˆ + ρφ&& + 2 ρ&φ& φˆ + &z&zˆ )

(c ) Para cualquier tiempo t, el valor de ρ (t ) y sus primeras derivadas son:

ρ (t ) = t t 2 + 4 ; ρ& =
2t 2 + 4
; ρ&& =
(
2t t 2 + 6 );
t +42
(t 2
+4 )
3/ 2

el valor de φ (t ) y sus primeras derivadas son:

φ (t ) = tan −1 (2t −1 );
2 4t
φ& = − ; φ&& =
t +4
2
(t 2
+ 4)
2

y el valor de z (t ) y sus derivadas son:

z (t ) = t 3 ; z& (t ) = 3t 2 ; &z&(t ) = 6t .

Reemplazando los valores correspondientes en las expresiones obtenidas en la parte b), tenemos:
r
r (1) = ρ (1)ρˆ + z (1)zˆ = 5 ρˆ + zˆ

v(1) = ρ& (1)ρˆ + ρ (1)φ&(1)φˆ + z& (1)zˆ =

6 2 5 ˆ
ρˆ − φ + 3 zˆ
5 5
 14 4 5  4 5 24  ˆ
r
[ ] [ ]
a (1) = ρ&&(1) − ρ (1)φ& 2 (1) ρˆ + ρ (1)φ&&(1) + 2 ρ& (1)φ&(1) φˆ + &z&zˆ =  3 / 2 −  ρˆ + 
25 
− φ + 6 zˆ
5  5 5 5

38
INTRODUCCION AL ANALISIS VECTORIAL
(JM Díaz, R. Castillo, F. Vega).

14.- Derivada direccional y Gradiente de una Función Escalar..

Como una manera de ver el campo escalar consiste en estudiar las superficies que se forman al
hacer Φ ( x, y, z ) = C (una constante); Consideremos que Φ ( x, y, z ) representa un campo escalar y
que Φ es una función de posición de valor único y diferenciable. Una superficie equipotencial
entonces, tiene la ecuación

Φ ( x, y , z ) = C (14-1)
donde C es una constante. Al asignar a C una sucesión de valores se obtiene una serie de
superficies equipotenciales. No habrá intersección de dos de dichas superficies debido a que Φ solo
puede tomar un solo valor. La figura 12 ilustra tres contornos en cada uno de los cuales la función
tiene un valor constante.

Φ3
z
Φ2

Φ1 p2
y
→
dl

x C1 p1 θ
C2 →
G

C3

Tomemos dos puntos cercanos, P1 sobre C1 y P2 sobre ella o fuera de ella; sí P1 = P1 ( x, y, z )

r
entonces P2 = P2 ( x + dx, y + dy, z + dz ) . El desplazamiento infinitesimal dl desde P1 a P2 puede
r
expresarse en términos de coordenadas rectangulares: dl = iˆdx + ˆjdy + kˆdz . Queremos evaluar la
diferencia de potencial entre P1 y P2 , la cual podemos hallar por los métodos del cálculo :

∂Φ ∂Φ ∂Φ
dΦ = dx + dy + dz (14-2)
∂x ∂y ∂z
r
Es útil definir el vector G , como

r ∂Φ ˆ ∂Φ ˆ ∂Φ
G = iˆ +j +k (14-3)
∂x ∂y ∂z
r r
efectuando ahora el producto punto G • dl tenemos que dΦ también puede expresarse como
r r
dΦ = G.dl (14-4)
r
Veamos ahora algunas propiedades de G :
r
Sí P2 está sobre la superficie C1 , dl es tangente a la superficie en P1 y dΦ = 0 porque Φ es
r
constante en la superficie; esto implica que G debe ser normal a la superficie en P1 .

39
INTRODUCCION AL ANALISIS VECTORIAL
(JM Díaz, R. Castillo, F. Vega).

r r
Para determinar la magnitud de G , escojamos P2 fuera de C1 tal forma que dl haga un ángulo θ
r r r
con la normal a C, como G es en la dirección normal, entonces dΦ = G • dl = G (cosθ )dl y
consecuentemente,

dΦ
= G cosθ
dl
r r
La componente de G en la dirección dl es la tasa de crecimiento de Φ en esa dirección. Esta
r
expresión es la derivada direccional de Φ en la dirección de dl .
r
Sí θ = 0 , dl viene a ser un elemento normal a la superficie equipotencial, escrito dn , y la derivada
r
direccional es máxima e igual a la magnitud de G , esto es

dΦ
=G (14-5)
dn

y es la máxima tasa de crecimiento de Φ .

Debido a que G coincide con la máxima tasa espacial de crecimiento de Φ , en dirección y magnitud,
r ∂Φ ˆ ∂Φ ˆ ∂Φ
recibe el nombre de Gradiente de Φ . Por esta razón, el vector G = iˆ +j +k se escribe
∂x ∂y ∂z
como
r
G = gradΦ = ∇Φ (14-6)

∂ ˆ ∂ ˆ ∂
donde ∇ ≡ iˆ + j + k se conoce como el operador NABLA (expresado en coordenadas
∂x ∂y ∂z
cartesianas).

De lo anterior, concluimos y enfatizamos nuevamente que el gradiente de un campo escalar Φ es

un vector que representa la máxima tasa de cambio, tanto en magnitud como en dirección, de Φ .

El gradiente es muy utilizado en Física, por ejemplo, la conducción de calor en una dimensión se
expresa como

dQ dT
= −κA (14-7)
dt dx

dQ dT
donde es la tasa de flujo de calor a través del área A, es el gradiente de temperatura y κ
dt dx
r
es la conductividad térmica. Definiendo la corriente de calor J Q como la tasa de calor que fluye a
través del área unitaria, un vector normal al área, la ley de conducción de calor puede expresarse
r
como J Q = −κ∇T .

40
INTRODUCCION AL ANALISIS VECTORIAL
(JM Díaz, R. Castillo, F. Vega).

r
Similarmente, la difusión de materia en un gas puede escribirse como J = − D∇n , donde n es la
r
densidad de partículas y J es la corriente de partículas.

Ejemplo Ilustrativo 14.1. Hallar el gradiente de la función Φ ( x, y, z ) = x 2 yz y la derivada direccional

de Φ en la dirección del vector unitario eˆ =

1
14
(iˆ + 2 ˆj + 3kˆ) en el punto (2,3,1) .
Solución.
r
Llamamos G al gradiente de Φ ; en cualquier punto, el gradiente de la función y la derivada
direccional en la dirección de ê , respectivamente, son:

r  ∂ ∂ ∂
G = ∇Φ =  iˆ + ˆj + kˆ (x 2 yz ) = iˆ(2 xyz ) + ˆj (x 2 z ) + kˆ(x 2 y )
 ∂x ∂y ∂z 

dΦ r
[ ]  1 ˆ
= G • eˆ = iˆ(2 xyz ) + ˆj (x 2 z ) + kˆ(x 2 y ) •  ( )
i + 2 ˆj + 3kˆ  =
1
(2 xyz + 2 x 2 z + 3 x 2 y )
dl  14  14
Evaluados en el punto (2,3,1) :
r
∇Φ = G (2,3,1) = 12iˆ + 4 ˆj + 12kˆ

dΦ 56
= ≅ 14.97
dl 14

41
INTRODUCCION AL ANALISIS VECTORIAL
(JM Díaz, R. Castillo, F. Vega).

Ejemplo Ilustrativo 14.2. Gradiente de la función de potencial de un sistema de N-cargas puntuales.

La función de potencial de un sistema de N cargas puntuales está dado por :
r N kq j
Φ ( x, y, z ) = Φ (r ) = ∑ r r (14-9)
j =1 r − r j
r
donde k es la constante dieléctrica, r es el vector de posición del punto donde se evalúa el potencial,
r
r j es el vector de posición de la j-ésima carga q j . ( )
Se puede demostrar que el campo electrostático es el negativo del gradiente de la función de
r
potencial: E = −∇Φ . Para el sistema de las N cargas puntuales tenemos que:

r  N q j  N  1 
E = −∇Φ = −∇ k ∑ r r  = − k ∑ q j ∇  r r 
 j =1 r − r j  j =1  r − r j 

Así, vemos que se requiere evaluar la expresión:

 1   ∂
[
 =  iˆ + ˆj ∂ + kˆ ∂  (x − x )2 + ( y − y )2 + (z − z )2 ]
1
∇ r r
−
2
 r −r   ∂x ∂y ∂z 
j j j
 j 

 = − (x − x j )i + ( y − y j ) j + (z − z j )k = − (rr − rr )
 1  ˆ ˆ ˆ
∇ r r
j
 r −r
 j

 [ ]
(x − x j )2 + (y − y j )2 + (z − z j )2 2
3 r r
r − rj
3

y el campo electrostático del sistema de N cargas puntuales será entonces:

N q (r − r )
r r
r
E = −∇Φ = k ∑ r r 3
j j
(14-9)
j =1 r − r
j

Expresión que concuerda con la hallada en el ejemplo ilustrativo 6.3.

42
INTRODUCCION AL ANALISIS VECTORIAL
(JM Díaz, R. Castillo, F. Vega).

Ejemplo Ilustrativo 14.3. Potencial y campo electrostático de un Dipolo Eléctrico. Un dipolo eléctrico
está constituido por una carga positiva y otra negativa separadas una distancia muy pequeña (2a, en
nuestro caso); el momento del dipolo es un vector dirigido de la carga negativa a la carga positiva y
r
cuya magnitud es el producto de la distancia que se para las cargas por la carga: p = 2aq .
La figura muestra un dipolo centrado en origen y paralelo al eje-z. Se quiere hallar la función de
r
potencial del dipolo Φ (r ) en un punto lejano y a partir de ella encontrar el campo electrostático del
dipolo.
P

r1
z
r r2
+q

θ
2a
y

-q
x
r 1 1
El potencial Φ (r ) en el punto P es Φ =
kq kq
− = kq −  .
r1 r2  r1 r2 
Usando el teorema de los cosenos podemos expresar r1 y r2 en términos de r y θ :

1 1
 a 
[ ]  
2 2
1 a a 2
r1 = a + r − 2ar cosθ
2 2 2
= r 1 − 2  cosθ +    ≅ r 1 − 2  cosθ 
 r  r    r 

1
−
1 1 a  2 1 a 
≅ 1 − 2  cosθ  ≅ 1 + cosθ 
r1 r  r  r r 
Similarmente,
1 1 a 
≅ 1 − cosθ 
r2 r  r 
 2a  kp
y el potencial lo podemos escribir como Φ = kq 2  cos θ = 2 cos θ
r  r
donde p ≡ 2aq .

r  ∂ 1 ∂  kp   2 cos θ   senθ  ˆ
El campo electrostático es E = −∇Φ = − rˆ + θˆ  cos θ  = kp  rˆ +  θ
 ∂r r ∂θ  r 2   r 3   r 3  

r  2kp cos θ   kpsenθ 

E = rˆ +  θˆ (14-10)
 r3   r 
3

43
INTRODUCCION AL ANALISIS VECTORIAL
(JM Díaz, R. Castillo, F. Vega).

Ejemplo Ilustrativo 14.4. Considere que la temperatura en una cierta región del espacio es
independiente del tiempo y está dada en grados centígrados por la ecuación T ( x, y ) = 30 x + 20 y ,
donde x, y están dados en metros; a) hallar las superficies isotérmicas; b) determinar la máxima
variación espacial de la temperatura y c) la dirección de máximo flujo de calor.

Solución.
a)- las superficies isotérmicas se obtienen al asignar un valor específico y constante para la
temperatura; en la siguiente tabla se da la ecuación de los planos isotérmicos para 20, 40 y 60
grados centígrados y la figura se muestra un corte de dichos planos.
T (oC) Superficie isotérmica
20 2y + 3x – 2 = 0
40 2y + 3x – 4 = 0
60 2y + 3x – 6 = 0
y

^
j 20
T = ^i 30 + ^j 20
∆

^i 30

O x

20ºC 40ºC 60ºC

b)- La máxima variación espacial la da el gradiente de temperatura:

∂T ∂T
∇T = iˆ +j = 30iˆ + 20 ˆj ; ∇T = 36 o C / m
∂x ∂y
que es un vector normal a las superficies isotérmicas; las componentes del gradiente dan la
variación espacial de la temperatura en la dirección de los ejes-coordenados.

c)- El calor fluye de mayor a menor temperatura y la dirección de máximo flujo de calor corresponde
a la dirección lo opuesta al gradiente.

15.- El operador nabla [∇ ] y algunas de sus propiedades.

Al final de la sección anterior definimos el operador nabla en coordenadas cartesianas como un

operador vectorial que cuando opera sobre una función escalar nos da un vector que se denomina
gradiente de la función. Este operador es de gran utilidad para expresar matemáticamente
conceptos como:

El gradiente de una función escalar, escrito como ∇Φ .

r
La divergencia de una función vectorial, escrita como ∇ • A .
r
El rotacional de una función vectorial, escrito como ∇ × A .
El Laplaciano de una función escalar, escrito como ∇ 2 Φ .

44
INTRODUCCION AL ANALISIS VECTORIAL
(JM Díaz, R. Castillo, F. Vega).

r
El Laplaciano de una función vectorial, escrito como ∇ 2 A .

Cada uno de estos temas serán considerados en las siguientes secciones. Por ahora es conveniente
ver que el operador nabla puede expresarse en otros sistemas de coordenadas. Aquí estamos
interesados en la forma que toma el operador nabla en coordenadas esféricas y en coordenadas
cilíndricas.

Para obtener el operador nabla en coordenadas cilíndricas partimos de su expresión en

coordenadas cartesianas. Los vectores iˆ, ˆj , kˆ en términos de los vectores unitarios ρˆ , φˆ, ẑ están
∂ ∂ ∂ ∂ ∂ ∂
dados por las relaciones (9-4). Los operadores , , expresados en términos de , ,
∂x ∂y ∂z ∂ρ ∂φ ∂z
se obtienen por los métodos del cálculo y el uso de las relaciones (9-1), de esta manera:

∂ ∂ senφ ∂ ∂ ∂ cos φ ∂ ∂ ∂
= cos φ − ; = senφ + ; =
∂x ∂ρ ρ ∂φ ∂y ∂ρ ρ ∂φ ∂z ∂z

∂ ∂ ∂
efectuando los productos y sumas indicados en ∇ ≡ iˆ + ˆj + kˆ tenemos que la expresión
∂x ∂y ∂z
del operador ∇ coordenadas cilíndricas es:

∂ 1 ∂ ∂
∇ = ρˆ + φˆ + zˆ (15-1)
∂ρ ρ ∂φ ∂z

En forma análoga se obtiene el operador ∇ en coordenadas esféricas. La relación (10 − 6 ) da los

vectores iˆ, ˆj , kˆ en términos de los vectores unitarios r̂ ,θˆ, φˆ y para expresar los operadores
∂ ∂ ∂ ∂ ∂ ∂
, , en términos de , , usamos las relaciones (10-1), así:
∂x ∂y ∂z ∂r ∂θ ∂φ

∂ ∂ cos θ ⋅ cos φ ∂ senφ ∂

= senθ ⋅ cos φ + −
∂x ∂r r ∂θ ρ ∂φ
∂ ∂ cos θ ⋅ senφ ∂ cos φ ∂
= senθ ⋅ senφ + +
∂y ∂r r ∂θ ρ ∂φ
∂ ∂ senθ ∂
= cos θ −
∂z ∂r r ∂θ

efectuando los productos y las sumas respectivos, tenemos que el operador ∇ en coordenadas
esféricas puede escribirse como
∂ ˆ1 ∂ 1 ∂
∇ = rˆ +θ + φˆ (15-2)
∂r r ∂θ rsenθ ∂φ

Entre las principales propiedades del operador ∇ podemos anotar:

Sí f y g son dos campos escalares; entonces ∇( f + g ) = ∇f + ∇g , y ∇( fg ) = f∇g + g∇f .

Sí f = f (u ) , donde u es un campo escalar; entonces ∇f =

df
∇u .
du

45
INTRODUCCION AL ANALISIS VECTORIAL
(JM Díaz, R. Castillo, F. Vega).

Sí f = f (u1 , u2 ,....un ) , donde los ui n son campos escalares; entonces

∂f ∂f ∂f n
∂f
∇f = ∇u1 + ∇u2 + ...... + ∇u n = ∑ ∇u j .
∂u1 ∂u2 ∂un j =1 ∂u j

Ejemplo Ilustrativo 15.1. Los vértices de un prisma están localizados en los puntos
(0,0,0); (1,0,0); (0,2,0); (0,0,1) ; (a) hallar el vector superficie de cada cara del prisma y demostrar que el
vector superficie total es cero; (b) ¿Cuánto vale la superficie de cada cara y cuál es el valor de la
superficie total?

Solución.
r r r
Sea a , b ,... f los vectores paralelos a los lados del prisma, como se indica en la figura.

(0,0,1)

→
c→ e→
f
→ (0,2,0)
→ b
a y
(1,0,0) →
d
x

a)-Los vectores de superficie de las caras del prisma serán:

2
(
r 1 r r 1
) ( ) (
2
)
S1 = d × f = 2 ˆj − iˆ × − iˆ + k = iˆ − 0.5 ˆj + kˆ
r
(
1 r r
) () ( )
S 2 = c × b = 0.5 kˆ × 2 ˆj = −iˆ
2
() ( )
r 1 r r
S 3 = (a × c ) = 0.5 iˆ × kˆ = 0.5 ˆj
2
r
(
1 r r 1
) ( ) ()
S 4 = − a × b = 2 ˆj × iˆ = −kˆ
2 2
La suma vectorial de estos cuatro vectores es nula.
r r r r
b)- La magnitud de cada vector superficie da el área de es cara: S1 = 1.5; S 2 = 1; S 3 = 0.5; S 4 = 1 ; la
suma aritmética representa el área total.

16.- Diferenciales de longitud, de área y de volumen.

Antes de continuar con el análisis y operaciones involucradas en la divergencia y el rotacional de un

r
campo vectorial, es conveniente revisar los conceptos diferencial de longitud, dl , diferencial de
r
área, dS , y diferencial de volumen, dV , en los sistemas de coordenadas cartesiano, cilíndrico y
esférico.

46
INTRODUCCION AL ANALISIS VECTORIAL
(JM Díaz, R. Castillo, F. Vega).

a) Coordenadas cartesianas.

En la figura 13 se muestra un elemento diferencial que tiene forma de un paralelepípedo de lados

r
dx, dy, dz con los vértices marcados de 1 a 8. Un desplazamiento de P1 a P2 que llamamos dl puede
escribirse como

4 3

5 8
dz
dr

1
2
dx
r + dr
6 7
dy
z
r
^
k

^j y
Fig. 13. Paralelipipedo diferencial.
^i

r r
dl = dr = iˆdx + ˆjdy + kˆdz (16-1)
El paralelepípedo tiene seis caras y el área de cada una de ellas se puede definir como el producto
del vector unitario normal a dicha cara por el valor de la superficie de la cara; sin embargo, esta
definición es ambigua debido a que si llamamos n̂ al vector unitario perpendicular a la cara, − n̂
también será perpendicular a dicha cara. Para evitar esta ambigüedad debemos especificar el
sentido de n̂ ; para la superficie de un cuerpo sólido la convención es tomar a n̂ en la dirección que
sale del cuerpo. Para otras superficies es necesario definir explícitamente la dirección de n̂ ; una vez
se haya encontrado el vector unitario normal entonces podemos expresar el vector elemento de
r
área como dS = nˆ dS . Así, para el paralelipípedo infinitesimal de la figura tenemos que:
r
dS x = iˆdydz para la cara 5,6,7, 8.
r
dS y = ˆjdydz para la cara 2,3,8,7 (16-2)
r
dS z = kˆdxdy para la cara 3,4,5,8

El volumen de este paralelipípedo será dV = dxdydz (16-3)

47
INTRODUCCION AL ANALISIS VECTORIAL
(JM Díaz, R. Castillo, F. Vega).

b) Coordenadas cilíndricas.
La figura 14 ilustra la porción diferencial cuyos vértices se han numerado del 1 al 8. Un
desplazamiento del punto 6 al punto 3 puede escribirse como
z

ρ
dρ
4
3
dz 5
8

1
2 r + dr
ρd φ 6
7
z z^
r
y
φ

φ^ O
x
ρ^
Fig. 14. Elemento diferenc ial en coordenadas cilíndricas.

r r
dl = dr = ρˆ (dρ ) + φˆ( ρdφ ) + zˆdz (16-4)

El área diferencial de una de las caras del elemento es el producto de la normal unitaria a esa cara
por el valor de la superficie de dicha cara. Así, tomando las caras cuyas normales son los vectores
unitarios ρˆ , φˆ, ẑ tenemos:
r
dS ρ = ρˆ ( ρdφdz )
r
dS φ = φˆ(dρdz ) (16-5)
r
dS z = zˆ (ρdρdφ )

El volumen diferencial será dV = (dρ )( ρdφ )(dz ) = ρdρdφdz (16-6)

c) Coordenadas esféricas.

La figura 15 ilustra la porción diferencial cuyos vértices se han numerado del 1 al 8. Un despla-
zamiento del punto 6 al punto 3 puede escribirse como

48
INTRODUCCION AL ANALISIS VECTORIAL
(JM Díaz, R. Castillo, F. Vega).

3
z

ρd φ=r·senθ·d φ 4
8 rdθ

ρd φ
5 2
dr rdθ
1 r·sen θ·d φ
θ

dr 7
y 6
φ dφ r + dr

r
x

Fig. 15. Elemento diferencial en coordenadas es féricas. O

r r
dl = dr = rˆdr + θˆ(rdθ ) + φˆ(rsenθdφ ) (16-7)

El área diferencial de una de las caras del elemento es el producto de la normal unitaria a esa cara
por el valor de la superficie de dicha cara. Así, tomando las caras cuyas normales son los vectores
r
unitarios r ,θˆ, φˆ , tenemos:
r
dS r = rˆ(rdθ )(rsenθdφ ) = rˆ(r 2 senθdθdφ )
r
dSθ = θˆ(dr )(rsenθdφ ) = θˆ(rsenθdrdφ ) (16-8)
r
dS φ = φˆ(dr )(rdθ ) = φˆ(rdrdθ )

El volumen diferencial será dV = (dr )(rsenθdφ )(rdθ ) = r 2 senθdrdθdφ (16-9)

49
INTRODUCCION AL ANALISIS VECTORIAL
(JM Díaz, R. Castillo, F. Vega).

17.- Integrales de Línea, de Superficie y de Volumen.

En los problemas físicos son de gran utilidad las curvas uniformes o regulares. Una curva uniforme
está constituida por uno o más arcos uniformes unidos unos a otros y que no se corta a así misma.
Usaremos los términos línea, trayectoria, curva, contorno intercambiablemente con el mismo
significado.
→→
A( r )
P
s u→
P0 b
→
a r
→
r0

O
Fig. 16. Porción de curva ab

r
La figura 16 ilustra una porción de curva ab con un punto fijo Po (ro ) y un punto variable
r
P (r ) = P (x, y, z ) . Llamamos s a la suma de las longitudes de los arcos que corresponden a las
diferentes posiciones de P en la curva y tomaremos a s como positiva o negativa según sea la
dirección Po P ; sí la trayectoria es abierta cada punto corresponde a un valor único de s y sí la curva
es cerrada el mismo punto puede diferir en un múltiplo entero del perímetro de la curva. El vector
r
r = OP es una función continua de s y a lo largo de cada arco uniforme esta función tiene una
r
r dr
derivada continua u = que es un vector unitario tangente a la curva en la dirección que s
dsr r
aumenta. El diferencial dr = u ds representa el vector que es tangente a la curva en el punto P y
cuya magnitud es ds.

a) Integral de línea.
r r r r
Consideremos ahora al vector A(r ) = A( x, y, z ) definido en cada punto P (r ) de la curva ab , la
r r r
componente tangencial de A a lo largo de la trayectoria ( A • u ) es una función de la longitud del
arco (s) y la integral curvilínea

b r r b r r
∫a
A • u ds = ∫ A • dr
a
(17-1)

se denomina integral de línea del vector A, a lo largo de la curva que se extiende desde a hasta
b.
r r
Debemos observar que el integrando, A • dr , es un escalar y la integral de línea, en general,
depende de la trayectoria o curva a lo largo de la cual se toma; cuando entre los puntos a y b hay
varias trayectorias, las integrales de línea a lo largo de las diferentes trayectorias no necesariamente
son iguales.

50
INTRODUCCION AL ANALISIS VECTORIAL
(JM Díaz, R. Castillo, F. Vega).

Cuando la trayectoria es cerrada, la integral de línea de la componente tangencial de un campo

r
vectorial A a lo largo de una curva cerrada se denomina CIRCULACIÓN DEL CAMPO y se
escribe como
r r r r
∫ A • dl ≡ Circulación de A(r ) (17-2)
C

r r
en esta expresión dl = dr y el símbolo de la integral utilizado significa que la trayectoria alrededor de
la cual se realiza la integración es cerrada.
r
Si el vector A es el gradiente de una función uniforme, la integral de línea es independiente de la
r r r
trayectoria y depende solamente de los puntos inicial y final debido a que A • dr = dr • ∇Φ = dΦ ,y

r r
A • dr = ∫ dΦ = Φ (b ) − Φ(a )
b b
∫
a a
(17-3)

La circulación en este caso es nula y la expresamos con la notación:

r r
∫ • dl = 0
C
A (17-4)

Consideremos a S como una superficie encerrada por el contorno C (figura 17a), uno de cuyos lados
se considera positivo y el otro negativo. Se define como dirección positiva a lo largo del contorno C
aquella en que el observador, en el lado positivo de la superficie recorre el contorno dejando siempre
a la izquierda la superficie S.

C C

a) b)
Fig. 17. Circulación.

Si dividimos ahora la superficie S en regiones más pequeñas por medio de arcos (figura 17b) la
r
suma de las circulaciones del campo vectorial A en la dirección positiva de estas regiones será
r r r r
∑∫ A • d r
Ci
= ∫ • dr
A
C

debido a que en esta suma la integral a lo largo de cada trayectoria común a dos áreas se presenta
dos veces en direcciones opuestas y no contribuye al resultado; las integrales que se presentan solo
una vez corresponden al contorno C. Esto es, cuando se divide una superficie en una serie de áreas
elementales, la circulación de un campo vectorial alrededor de su contorno es igual a la suma de las
circulaciones de este campo alrededor de dichas áreas.

b) Integrales de superficie.
Con el propósito de la integración, es conveniente dividir la superficie en porciones o elementos
regulares de superficie, cada uno de estos elementos regulares de superficie para una cierta

51
INTRODUCCION AL ANALISIS VECTORIAL
(JM Díaz, R. Castillo, F. Vega).

orientación de los ejes coordenados tiene por proyección sobre el plano x-y el interior de una curva
uniforme cerrada que se puede representar por una ecuación de la forma z = f ( x, y ) en la cual
f ( x, y ) y sus derivas con respecto a y a x son continuas.

Consideremos ahora la integración sobre un elemento regular de superficie S, para ello

Observamos que en cada punto P ( x, y, z ) de la superficie S existe un plano tangente definido y una
r
normal n̂ ; un elemento de área dS = nˆ dS en el punto P puede definirse como el área en el plano
tangente que tiene por proyección sobre el plano xy el elemento de área dxdy . Esto significa que
r
dS • nˆ = (nˆ • zˆ )dS = (cos θ )dS = dxdy

donde θ es el ángulo formado por la normal y el eje-z en el punto P. El área total será

S = ∫ dS = ∫∫ (sec θ )dxdy
S

en esta expresión la integral doble se extiende sobre la proyección de S sobre el plano xy. Sí,
f ( x, y, z ) es una función escalar definida en todos los puntos de la superficie S, entonces la
integral de f sobre la superficie S se escribe como

∫S
fdS (17-5)

y puede definirse como la integral doble ∫∫ f (secθ )dxdy extendida sobre la proyección de la
superficie S sobre el plano xy.
r
Análogamente, sí A( x, y, z ) es un campo vectorial definido en cada punto P de la superficie S y n̂
r
es el vector unitario sobre el lado positivo de la superficie, la componente normal de A será
r r
An = A • nˆ y la Integral de superficie de A sobre la superficie S se puede escribir como
r r r
∫S
An dS = ∫ A • nˆ dS = ∫ A • dS
S S
(17-6)
Otra manera de definir la integral de superficie consiste en un proceso de suma extendida sobre la
superficie directamente. Dividimos el área total en un gran número de elementos pequeños de
r
superficie, cada uno de los cuales tendrá un vector área dS j = nˆ j dS j y luego los sumamos. En esta
forma, la integral de superficie de un campo vectorial se define como
r r r r
lim dS j →0 ∑ A • nˆ j dS j = ∫ A • nˆ dS = ∫ A • dS
S
S

Como n̂ y − n̂ son vectores perpendiculares a la superficie y ambos satisfacen la condición

nˆ • nˆ = 1 , existe una ambigüedad en la definición. Para evitarla debemos especificar la dirección del
vector unitario; en los sólidos por convención se toma como vector unitario al que sale del sólido,
para otras superficies es necesario definir explícitamente su dirección. Una vez conocido n̂ ,
r
podemos definir el elemento vectorial de área dS = nˆ dS y entonces,

52
INTRODUCCION AL ANALISIS VECTORIAL
(JM Díaz, R. Castillo, F. Vega).

r r r
∫
S
A • ˆ
n dS = ∫ A
S
• d S

c) Integrales de Volumen.
En la Sección 15 estudiamos los diferenciales de volumen y se halló sus expresiones en cada uno
de los principales sistemas de coordenadas (cartesianas, cilíndricas y esféricas). Aquí simplemente
r r s
definiremos las integrales de volumen de un campo escalar ρ (r ) y de un campo vectorial A(r ) ,
respectivamente, como:
r r
∫ ρ (r )dV = ∫∫∫ ρ (r )dV (17-7)
V
V

r r r r
∫ A(r )dV = ∫∫∫ A(r )dV (17-8)
V
V

53