ANTONIO TAPIA SILVA VECTORES 1
CAPÍTULO A
INTRODUCCIÓN A LOS VECTORES.
Este capítulo se inicia con los conceptos más elementales de vectores en el plano
cartesiano y en el espacio y son constantes en magnitud y en dirección; pero esto,
tan sólo es el inicio de los conceptos generales. Es nuestro punto de partida. Los
estudiantes del primer curso de Física solamente abarcan este acápite.
Se continúa con una introducción general del estudio sobre el gradiente,
divergencia y rotacional, resaltando unas tantas algunas interpretaciones de su
significado físico al igual que el teorema de Stokes.
Por la simetría de los problemas es necesario el uso de las coordenadas cilíndricas
y esféricas, por lo que se realiza la transformación de coordenadas rectangulares a
curvilíneas, definiendo sus vectores unitarios, factores de escala y haciendo
algunas aplicaciones en coordenadas esféricas y cilíndricas, para hallar su
diferencial de línea, su elemento de área y de volumen, al igual que su gradiente
divergencia y su rotacional. Es muy necesario saber utilizar la transformación de
las diferentes operaciones en los principales sistemas de coordenadas. Esta es la
parte más importante del capítulo de “introducción a los vectores”
El estudiante tendrá a su alcance una herramienta muy útil para el tratamiento
vectorial en coordenadas de cualquier sistema y podrá despertar en él, el deseo de
investigar en este campo que tiene aplicaciones prácticas muy importantes.
A.1 MAGNITUDES ESCALARES Y VECTORIALES
A.1.1. MAGNITUD ESCALAR. - Es aquella magnitud que queda bien definida
cuando se indica tan sólo un número (cantidad con respecto a la unidad de
medida) y la especie. Ejemplo: la masa, la temperatura, el volumen, etc.
Si decimos que la masa de un objeto es de 60 kg., intuitivamente entendemos
comparando con la unidad de medida que es la unidad de masa y queda bien
comprendida esta magnitud.
A.1.2 MAGNITUD VECTORIAL. - Es aquella que además de tener una
magnitud escalar, se necesita conocer su dirección. Es decir, saber hacia dónde se
dirige, para que quede bien definida.
Ejemplo la velocidad, la fuerza, etc.
A.1.3 CONCEPTO DE DIRECCIÓN
Una línea recta orientada define una dirección. Podemos movernos a lo largo de
una recta. La recta puede tener dos sentidos y acordamos definir el sentido
positivo y negativo a uno y otro extremo respectivamente. Necesitamos tener
algún punto o sistema de referencia para indicar la dirección.
ANTONIO TAPIA SILVA VECTORES 2
En el plano nos basta tener una línea recta para indicar la dirección. Usamos la
recta de las abscisas o eje de las “x”, que es paralela a la línea tangente a las aguas
en reposo. Hacia la derecha consideramos positivo y hacia la izquierda es
negativo.
Para definir la dirección que hace una recta en un plano, es suficiente indicar el
ángulo que forma esa recta con el eje X positivo.
En el espacio, necesitamos como mínimo dos ángulos, para definir una dirección.
Para esto tomaremos el eje “y”, que es paralelo a la dirección de la recta de la
plomada y perpendicular al eje X.
Hacia arriba es positivo y hacia abajo es negativo. El punto donde se intersecan el
eje X y el eje Y es el origen de coordenadas. El eje Z es perpendicular al plano
que forman X e Y.
1.4 VECTOR. - Es una cantidad física que requiere el conocimiento de una
magnitud escalar (o módulo) y su dirección. Se representa por un segmento
orientado.
A.1.5. Elementos de un vector. -
El módulo o parte escalar: es lo que representa físicamente la longitud del
segmento.
Ejemplo si su longitud es 10 cm. Y cada centímetro nos representa 2 N, entonces
su módulo es 20 N.
Fig. A.3
- El origen es el punto “O” (inicio del vector).
- El extremo. Es el punto “P”.
ANTONIO TAPIA SILVA VECTORES 3
- Dirección. Es la de la recta que la contiene.
- Sentido. Hacia donde apunta la flecha.
A.1.5. REPRESENTACIÓN DEL VECTOR.
Gráficamente se representa por el segmento orientado.
Analíticamente se representa por dos letras (primero se indica el origen y
luego el extremo con una flechita encima de izquierda a derecha; y se lee: “Vector
OP”; o con una sola letra con una flechita encima de izquierda a derecha y se
lee:
= “Vector A”.
En algunos libros se usa el símbolo a negrilla; es decir, los bordes de la letra algo
más ancho y más oscuro que lo normal.
A= Vector A
El módulo es
Se lee módulo de A; o simplemente A= A
A.1.7 VECTOR UNITARIO
Es un vector de módulo la unidad y adimensional, que solamente indica dirección.
Se define como la relación entre el vector y su módulo
UA = A/A … (A.1)
Entonces A = A UA
En cada sistema de coordenadas se define un sistema de vectores unitarios que
indican solamente dirección. Así en cartesiana tenemos los vectores i, j, k que nos
indican direcciones en los ejes X, Y, Z y son constantes en su módulo (=1) y en su
dirección.
ANTONIO TAPIA SILVA VECTORES 4
Vectores paralelos. Tienen la misma dirección y por lo tanto el mismo vector
unitario
Negativo de un vector. El negativo de un vector A es otro vector -A, que tiene
igual módulo, pero dirección opuesta
Fig.A.6
Los vectores pueden trasladarse sobre su misma recta que los contiene o
paralelamente sin que altere sus propiedades.
A.1.8. ADICIÓN DE VECTORES
Dados dos vectores A y B, entonces definimos LA SUMA VECTORIAL o
RESULTANTE de los vectores A Y B como un nuevo vector que se obtiene de la
siguiente manera:
Fig.A.7
1º Trasladamos el origen del vector B a continuación del extremo del vector A,
2º Unimos el origen del vector A con el extremo de B en su nueva ubicación. El
vector que resulta de esta unión es el vector suma o resultante R =A+B
O bien podemos trasladar el vector A paralelamente a continuación del vector B y
unimos el origen de B con el extremo de A y el vector cuyo origen es el origen
de B y el extremo de A en su nueva posición, es el vector resultante R=B+A
Fig.A.8
Como podemos ver, la suma vectorial de dos vectores A y B es conmutativa:
A+B = B+A … (A.2)
Módulo de la resultante. Se aplica la ley del coseno.
Fig.A.9
En la figura
ANTONIO TAPIA SILVA VECTORES 5
… (A.3)
Siendo el ángulo que forman A y B
A.1.9. DIFERENCIA DE VECTORES
Dados los vectores A y B, entonces definimos la diferencia D = A - B como la
suma de A con el negativo de B.
D = A + (- B)
Fig.A.10
En la gráfica es la otra diagonal que une los extremos. En el extremo de B, va el
origen de D y en el extremo de A va el extremo de D.
En donde su módulo es
… (A.4)
SUMA DE VARIOS VECTORES
Para sumar más de dos vectores A + B + C+ ... + N, se trasladan paralelamente
cada vector a continuación del anterior y se une el origen de A con el extremo del
último vector trasladado.
A B C
B D
C R
D Fig.A.11
La resultante es el vector cuyo origen coincide con el origen del vector A y el
extremo del extremo del vector enésimo.
A.1.10 COMPONENTES DE UN VECTOR
Cualquier conjunto de vectores que al sumarse dan un nuevo vector R, entonces
esos vectores del conjunto, se dice que son componentes de R.
Las posibilidades son infinitas, pero nos importan en especial los vectores
componentes que se orientan según los ejes del sistema cartesiano: dos en el plano
ANTONIO TAPIA SILVA VECTORES 6
y tres en el espacio.
Así el vector A, puede considerarse como la suma de dos vectores componentes
perpendiculares.
A=Ax+Ay … (A.6)
Ax en la dirección del eje X y Ay en la dirección del eje Y.
Fig.A.12
Y en función de sus vectores unitarios
Por consiguiente.
Y su módulo
La dirección
También se puede ver que:
Ax = A Cos
Ay = A Sen
A = A Cos i + A Sen j … (A.7)
El vector unitario UA es:
UA= A/A = Cos i + Sen j … (A.8)
Y su módulo es: UA= (Cos2 + Sen2 ) ½
=1
En el espacio, consideramos una coordenada z, en cuyo caso.
Ax=ASen Cos
Ay=A Sen Sen
Az= A cos
ANTONIO TAPIA SILVA VECTORES 7
Fig.A.13
Y por consiguiente
A = Ax i + Ayj+ Azk
A= (Ax2+Ay2+Az2) ½
Si , y son los ángulos que hace el vector A con los ejes X, Y y Z
respectivamente y se denominan los ángulos directores, entonces.
Ax= ACos
Ay= ACos
Az= ACos
se denominan los cosenos directores.
Como se puede notar en la figura, A x, Ay y Az son las coordenadas del punto
P(x,y,z) es decir, si consideramos.
A=
x= ASen Cos
y= A Sen Sen
z= A Cos
Entonces tendremos el vector posición de origen en O y extremo en P, dado por:
… (A.9)
Si tenemos dos puntos P1 y P2 con coordenadas P1 (X1, Y1, Z1) y P2 (X2,Y2, Z2)
ANTONIO TAPIA SILVA VECTORES 8
Fig.A.14
En la figura vemos que
ó
… (A.10)
Para sumar (o restar) vectores se suman (o restan) componente a componente
Ejemplo: Sean los vectores
A = 2 i + 5 j +3 k N
B = 3 i + 6 j –2 k N
C = -2 i + 3 j – 5 k N
Hallar A + B – C (su módulo y dirección).
A+B-C = 7 i + 8 j + 6 k
Su módulo es r= (49+64+36) ½ r =12,2 N
La dirección nos da los cosenos directores
Cos = 7/r =54,98º
Cos = 8/r =49,02º
Cos = 6/r =60,54º
A.1.11 PRODUCTO DE UN ESCALAR POR UN VECTOR:
El producto de un escalar m por un vector A es otro vector paralelo (o de sentido
contrario) al vector A. Si la multiplicación es por un número puro m, entonces si
m > 1 lo amplía, si m=1 es el mismo vector, si m < 1 lo reduce y si m = 0 lo anula.
Pero si m tiene unidades entonces cambia completamente la unidad del vector y
tan sólo es un vector paralelo.
Ejemplo: Si A = 20 j m/seg2
m = 2 kg.
Entonces:
mA = 40 j N
Que es completamente diferente al vector A (tan sólo es paralelo)
A.1.12 LEYES DEL ÁLGEBRA VECTORIAL:
1) A + B = B + A Conmutatividad.
2) A + (B+C) = (A +B) + C Asociatividad.
3) mA = A m Conmutatividad con respecto a un coeficiente adimensional.
4) m(nA) = (mn) A = m (n A) m y n son adimensionales.
5) m(A+B) = m A + m B Distributiva con respecto a un coeficiente adimensional
ANTONIO TAPIA SILVA VECTORES 9
EJEMPLO DE APLICACIÓN
Dos ciudades A y B están situadas una frente a la otra en las dos orillas de un río,
de 8 Km. de ancho, siendo la velocidad del agua de 4 Km/h. La ciudad C se
encuentra en la misma rivera de B y a 6 Km aguas arriba de B. Si un hombre
quiere ir de la ciudad A a la ciudad C y utiliza una embarcación que desarrolla una
velocidad máxima de 10 km/h. ¿Cuál será la dirección que debe seguir para llegar
en el menor tiempo posible? ¿Y cuánto tarda?
Solución:
Fig.A.15
La resultante debe ir en la dirección AC.
En la figura:
Es decir:
Igualando en i y en j
Y resolviendo, se tiene que
Sen = 0,821
= 55,18º
Es decir, haciendo un ángulo de 34, 82o con la ribera a la izquierda del río aguas
arriba.
El tiempo 1 h, 25 min.
Otro ejemplo. - En el paralelogramo de la figura los segmentos AQ y AP dividen
a los vectores a y b en dos partes iguales, respectivamente. Halle la relación entre
los módulos de los vectores generados por los puntos DE y DB; FD y DE.
ANTONIO TAPIA SILVA VECTORES 10
Fig.A.16
En efecto:
...(*)
Entonces: en (*):
Factorizando en función de los vectores y tenemos:
De donde:
;
Entonces
n= 1/3
Si tomamos en los puntos F, P y B, se obtiene el mismo resultado. Es decir que lo
triseca en partes iguales.
A.1.13 PRODUCTO ESCALAR (o producto INTERNO, o producto punto)
Dados dos vectores A y B, entonces definimos que el producto escalar o producto
interno, o producto punto A.B, como el producto de sus módulos por el coseno del
ángulo que forman.
(0<<) … (A.11)
(Conmutatividad)
El producto es un ESCALAR carece de dirección; (es un número con su
unidad). Con esto:
i.i = (1) (1)Cos 0 = 1
j.j= (1) (1)Cos 0 = 1
k.k= (1) (1)Cos 0 = 1
i.j = (1) (1) cos 90 = 0
i.k= (1) (1) cos 90 = 0
j.k= (1) (1) cos 90 = 0
Es decir si Un y Um son vectores unitarios (o versores) de un sistema
ANTONIO TAPIA SILVA VECTORES 11
ortogonal,
1 si n = m
Un. Um = nm = … (A.12)
0 si n m
nm se denomina DELTA DE KRONECKER.
De la definición vemos que: A.A=A2
Ya que A.A= AA cos=A2
Si
Entonces
=
+
+
De donde:
Es decir que el producto escalar del vector A por sí mismo es la suma de los
cuadrados de sus componentes
De igual forma se puede demostrar que si
… (A.13)
Es decir que el producto escalar de A y B es la suma del producto de sus
componentes.
Ejemplo: Si A= 2i + 3j + 6k N B= 2i + j- k m
A.B = 4+3-6 = 1 Joule
Como A.B = AB cos
ANTONIO TAPIA SILVA VECTORES 12
Fig.A.17
BCos , representa gráficamente, la proyección del vector B sobre el vector A
Ejemplo: Hallemos los ángulos que forman las diagonales de un cubo
Los vectores son:
A= B=
A. B= = = ABcos
Cos =
Entonces el ángulo es:
=70,53º
A.1.14 PRODUCTO VECTORIAL (O PRODUCTO EXTERNO, O
PRODUCTO CRUZ)
Dados dos vectores A y B, entonces se define el producto vectorial de A x B como
otro vector C talque C = A x B, cuyo módulo es el producto de los módulos de los
vectores A y B , por el seno del ángulo que forman y su dirección es
perpendicular al plano que forman los vectores A y B, en sentido de un tornillo de
roscas, como si girara de A hacia B.
… (A.14)
es un vector unitario perpendicular al plano que forman y . Se sigue la
dirección de la regla de la mano derecha: Se apuña de A hacia B y a donde señala
el pulgar derecho, es la dirección. Sigue la regla del tirabuzón. En espiral.
ANTONIO TAPIA SILVA VECTORES 13
Fig.A.18
El producto vectorial no es conmutativo.
A x B = -B x A … (A.15)
Es antisimétrico.
El módulo del producto vectorial representa el área del paralelogramo que generan
los vectores A y B
Fig.A.19
… (A.16)
Es el área del paralelogramo
Con esto A x A = 0
ya que AxA= AA Sen0= 0 BxB = 0
Para un sistema cartesiano trirrectangular:
ANTONIO TAPIA SILVA VECTORES 14
Fig.A.20
ixi=0; j x j = 0, k x k =0
(A.17) ixj=k jxk=i kxi=j
j x i = -k kxj=-i ixk=-j
Será fácil acordarse si seguimos el orden de giro antihorario como positivo; en el
sentido horario, es negativo.
Fig.A.21
Así:
AxB=
= +
+ … (A.18)
Entonces:
= +
Acomodando:
= … (A.19)
Esta última expresión es el desarrollo del determinante
… (A.20)
A.1.15 PRODUCTOS TRIPLES:
Sean A, B, C tres vectores no coplanares, entonces el producto triple A.(BxC)
representa el volumen del paralelepípedo que generan.
ANTONIO TAPIA SILVA VECTORES 15
Es fácil ver que:
i j k
A.(BxC) = A. Bx By Bz
Cx Cy Cz
Y desarrollando es fácil obtener que
Ax Ay Az
[Link] = Bx By Bz … (A.21)
Cx Cy Cz
Se demuestra que
EJEMPLO. Dado los vectores:
A =2 i - j +3k m; B = 3i +5 j – k m; C = i +2j – 3 k m
Hallar: a) (AxB)xC b) Ax (BxC) c) A· (BXC) d) (AxB)·C
Resp: (a) 76.4m3 (b) 25.67m3 (c)-31 m3 d) -31 m3
A.1.16 ALGUNAS PROPIEDADES DEL PRODIUCTO VECTORIAL
Es fácil demostrar que
AxB = -BxA
(A.B) C (B.C) A
Ax (BxC) (AxB)xC
Ax (BxC) = (A.C) B - (A.B) C
(AxB) xC = (A.C) B-(B.C) A
A.1.17 VECTORES RECÍPROCOS
Dados los conjuntos de vectores a,b , c y a´,b´,c´ y si [Link]0 y sí
ANTONIO TAPIA SILVA VECTORES 16
; … (A.22)
Además, vemos que
; ... (A.23)
Entonces estos vectores se dicen que son recíprocos.
A.1.18 CAMPOS ESCALARES Y VECTORIALES
.
Lo tomamos en cuenta desde un punto de vista físico. Si en cada punto del
espacio , se le puede asociar un escalar , entonces habremos
definido un campo escalar en la región. Por lo tanto, la función es una
función de escalar de posición, o bien, función del punto escalar. Ejemplo, si
tenemos una función de la temperatura .
Si en cada punto de una región del espacio se puede asociar un vector
, entonces hemos definido un campo vectorial de en esa región. La
función vectorial depende del punto, por lo cual se le llama función vectorial de
posición o función vectorial de punto. Por ejemplo, si tenemos una función de
velocidades dada por
Si los campos son independientes del tiempo, se denominan estacionarios.
CAPÍTULO B
ANTONIO TAPIA SILVA VECTORES 17
GRADIENTE, DIVERGENCIA Y ROTACIONAL
1. GRADIENTE
Acordémonos de la expresión:
... (B.1)
Y si está en el espacio tridimensional
f = f(x, y, z)
Siendo su derivada total:
... (B.2)
Está expresión representa el incremento de magnitud de la función f, si es que va
del punto (x, y, z) al punto (x + x, y + y, z + z).
Si tenemos que:
Y:
Entonces podemos escribir (B.2) como:
.... (B.3)
Donde se ha puesto la expresión:
... (B.4)
f, se denomina GRADIENTE de la función escalar f.
El operador (NABLA) es un PSEUDO VECTOR y es:
... (B.5)
Por sí solo no tiene significado físico. Al multiplicarse por una escalar queda
definida la operación gradiente de la función (B.4).
La ecuación (B.3) también puede ser escrita como:
df = f. d = ... (B.6)
ANTONIO TAPIA SILVA VECTORES 18
Donde es el ángulo entre f y el vector desplazamiento de posición . Es
paralelo (colineal) con f cuando el ángulo = 0. El gradiente representa la
dirección del máximo cambio de f.
B.2 DERIVADA DIRECCIONAL
Es la componente de f en la dirección de un vector unitario y es igual a:
… (B.7)
También podemos decir que la derivada direccional es la razón de cambio de la
función f por unidad de longitud.
Y se puede partir de (B.3) y dividir por dr.
Su significado físico de f depende de los que signifiquen f o represente
físicamente. Pero podemos decir que el gradiente es un vector cuyo módulo es la
máxima derivada direccional en el punto donde lo estamos considerando y su
dirección es la que tiene la máxima derivada en ese punto.
Ejemplo 1. Supongamos que tenemos la función f = 2x³yz – 4x voltios que nos
representa un potencial eléctrico. Hallemos la derivada direccional en el punto
P(1, 2, 3) y en la dirección y sentido del vector posición que pasa por el punto P(3,
4, 0).
En este caso el gradiente f nos representa el negativo de un campo eléctrico
El vector unitario en el punto (3, 4, 0) es:
Entonces:
y en el punto P(1, 2, 3)
ANTONIO TAPIA SILVA VECTORES 19
Podemos dar un concepto físico de a partir del concepto de derivada
direccional y de vector unitario.
Así:
(versor tangente a la curva)
=
Donde es el ángulo entre y la normal a la curva.
El máximo valor de esta última expresión es cuando el coseno del ángulo es igual
a la unidad. Es decir, cuando . Entonces:
… (B.8)
B.3 DIVERGENCIA.
La divergencia de un campo vectorial F, podemos definirla como el flujo por
elemento de volumen
... (B.8)
La expresión:
Representa el flujo de un campo vectorial F, a través de una superficie .
Donde es el vector diferencial de área y su dirección es perpendicular a la
superficie que envuelve al campo. Puede ser saliente (manantial) o entrante
(sumidero) a la superficie. El elemento de área es:
=
Donde n es un vector unitario perpendicular a la superficie, conocido mayormente
como “NORMAL”
Para una superficie cerrada:
... (B.9)
La convención es en el sentido de , saliente del
volumen limitado por la superficie S y a la vez es
perpendicular a S.
Usemos la ecuación (B.9) a un volumen de tamaño diferencial, al cual
consideramos como un prisma recto. Usemos coordenadas rectangulares como en
la figura 2.
ANTONIO TAPIA SILVA VECTORES 20
Fig.2
Para hallar el flujo total, tomemos el flujo en cada superficie y ordenando
tendremos en los tres pares de superficies:
... (B.10)
Las superficies con signo “prima” están a distancias diferenciales de las
superficies sin prima.
Como las superficies son de tamaño diferencial, las componentes de F son
aproximadamente constantes a lo largo de la superficie, por lo cual, las integrales
de superficie en (B.10) vienen a ser simples multiplicaciones de la componente de
F perpendiculares a la superficie.
El flujo será, entonces:
(B.11)
es el elemento de volumen encerrado por la superficie; y en el
límite cuando el elemento es un diferencial de volumen:
... (B.12)
Si
Haciendo:
... (B.13)
ANTONIO TAPIA SILVA VECTORES 21
Vemos que esta expresión es un escalar y la llamamos DIVERGENCIA del
vector F y, por consiguiente:
… (B.14)
Es decir:
; …
O lo que es lo mismo
Con lo cual demostramos la ecuación (B.8).
En coordenadas curvilíneas es preferible tomar la ecuación (B.8).
Para demostrar la ecuación (B.13), tenemos que:
Como podemos notar, es un escalar.
Significado físico
Es obvio que depende de lo que F represente físicamente. Supongamos que
Donde
Masa por unidad de volumen (densidad de la masa) en
movimiento
Vemos que la divergencia físicamente representa la variación de la masa del
fluido en la unidad de volumen, por unidad de tiempo con que se mueve con
referencia al punto P.
EJEMPLO 2. Evaluar F en el punto P (1, 0, 2), si N
En efecto:
Derivando y reemplazando en el punto, tenemos que
ANTONIO TAPIA SILVA VECTORES 22
F= 4 N/m
B.4 TEOREMA DE LA DIVERGENCIA:
Dada una función vectorial F de una superficie cerrada S, entonces se cumple que:
Podemos partir de la definición de divergencia, entonces:
y tenemos
. … (B.14)
Es el teorema de la divergencia de Gauss.
EJEMPLO. Dado el campo vectorial F = r, donde r es un vector posición,
compruebe el teorema de la divergencia para un cubo de lado c de la figura
r= x i + y j + z k
Fig.B.3
=
Por otro lado:
= 3V
=
Veamos al teorema de la divergencia (O DE GAUSS) de otra forma
ANTONIO TAPIA SILVA VECTORES 23
Este importante teorema, se puede enunciar de la siguiente manera:
«La integral de superficie de la componente normal a un vector A, a través de
una superficie cerrada, es igual a la integral de la divergencia de A en todo el
volumen encerrado por dicha superficie».
Si A = Axi + Ayj + Azk
Veámoslo de otra manera
Sea S una superficie cerrada en la que toda recta paralela a los ejes coordenados la
corta a lo más en dos puntos.
Supongamos que las ecuaciones de las superficies límites inferior S 1 y superior S2
son:
Z1 = f1(x,y) y Z 2= f2(x,y) respectivamente
Y si llamamos R a la proyección de la superficie sobre el plano XY, tendremos
entonces que
Fig.B.4
En la cara inferior S2
ANTONIO TAPIA SILVA VECTORES 24
Por lo tanto:
Entonces
Es decir:
Y en forma análoga:
Sumando:
B.5 ROTACIONAL.
La circulación de un campo vectorial F a lo largo de una trayectoria cerrada, la
obtenemos usando la integral de línea. Es decir:
... (B.15)
En donde, por ejemplo, si F es una fuerza, C representa el trabajo realizado a lo
largo de la trayectoria
Consideremos el entorno infinitesimal de la figura
ANTONIO TAPIA SILVA VECTORES 25
Fig.5
...
(B.16)
Las componentes de F son casi constantes en cada tramo de tamaño diferencial y
acomodando, la ecuación (B.16) queda aproximadamente como:
... (B.17)
Y en el límite, cuando . Además, el producto de los incrementos
es el área en la dirección del eje Z y tendremos que:
Además, notamos que
; y
Representan las derivadas parciales; por lo que la circulación es
... (B.18)
Donde el elemento de área tiene la dirección perpendicular al plano XY (en
)
Y por analogía:
En el plano XZ ... (B.19)
Y En el plano YZ ... (B.20)
Las ecuaciones (B.18), (B.19) y (B.20) son las componentes rectangulares del
producto vectorial del operador (nabla) con la función F. Esto es EL ROTOR
DE F.
ANTONIO TAPIA SILVA VECTORES 26
... (B.21)
Podemos expresar la circulación como
... (B.22)
Podemos expresar el elemento de área en su forma vectorial y tenemos que en
coordenadas cartesianas viene dado por:
O bien:
... (B.23)
Siendo n un vector unitario perpendicular a la superficie siguiendo la regla de
la mano derecha. Es decir, recorriendo en el sentido de la circulación; apoyando
los dedos en ese sentido y hacia donde apunta el dedo pulgar (como se indicó) en
el capítulo de producto vectorial.
Y reemplazando (B.21) en (B.15), tendremos que:
Por lo que el rotacional es:
... (B.24)
es la componente del diferencial de área en dirección de la normal.
B.6. TEOREMA DE STOKES
Si reunimos muchos contornos de líneas diferenciales, como el de la figura, hasta
formar una superficie macroscópica S, entonces cada pequeño contorno genera
una contribución a la circulación
Y la circulación total será
Y teniendo en cuenta la ecuación (B.15), tenemos la ecuación
... (B.25)
Esta última expresión se denomina LEY DE STOKES
.
ANTONIO TAPIA SILVA VECTORES 27
Ejemplo 4.- Un vector E es irrotacional si xE=0: determinar las constantes A,
B y C si:
E= (x + 2y + Az)i + (Bx – 3y –z)j + (4x+cy+2z)k N/Coul y
demuestre también que puede expresarse como un escalar.
Solución:
De donde:
c= -1; A= 4; B= 2
Supongamos que E=V, entonces tendremos que
Integrando en x
E integrando en forma similar en y e z, se tiene que:
De donde el potencial es +C
B.7 ALGUNAS FÓRMULAS IMPORTANTES:
1.
2.
3.
4.
5.
6. es el LAPLACIANO de la función. En coordenadas
cartesianas:
7.
8.
9. y
ALGUNOS CASOS PARTICULARES DEL TEOREMA DE LA
DIVERGENCIA. - FÓRMULAS DE GREEN
ANTONIO TAPIA SILVA VECTORES 28
FORMA INTEGRAL DEL OPERADOR NABLA
Se da por la fórmula
: … (33)
Donde el símbolo o del subíndice significa un producto simple (ej: gradiente), un
producto escalar o un producto vectorial).
Si o del subíndice significa un producto escalar:
Si es producto cruz:
Si es gradiente:
LAS IDENTIDADES DE GREEN
En el teorema de la divergencia hagamos:
A=
Pero:
.( ) = (.) + ().()
= 2 + ().()
ANTONIO TAPIA SILVA VECTORES 29
Entonces:
Es la primera identidad de la Green
Si cambiamos por un en (*)
Restando de (*)
Este es el teorema recíproco o segunda identidad de Green.
TEOREMA DE GREEN EN EL PLANO
Supongamos Que tenemos una región R cerrada del plano XY limitada por una
curva simple y cerrada C. Sean M y N dos funciones continuas de X e Y
derivables en y sus derivadas también son continuas entonces tenemos que:
(1)
Cuando se recorre en dirección contraria a las agujas del reloj (sentido horario).
Este es un caso particular del teorema de la divergencia.
El teorema de Green establece que:
(2)
Llamado también PRIMERA IDENTIDAD DE GREEN.
El teorema de la reciprocidad o simétrico, llamado también «SEGUNDA
IDENTIDAD DE GREEN», establece que:
(3)
Trataremos de demostrar la ecuación (1)
ANTONIO TAPIA SILVA VECTORES 30
Sean y1(x) é y2(x) las ecuaciones las curva AEB y AFB respectivamente donde R
es la región limitada por la curva C.
Es decir:
=-
En forma análoga x = X1(x) y X2(x) son las ecuaciones de las curvas EAF y EBF
respectivamente entonces deducimos que:
Es decir, sumando
Se puede extender la demostración del teorema de Green en el plano, al caso en
que la curva se corta en más de dos puntos por restas paralelas a los ejes
coordenados.
ANTONIO TAPIA SILVA VECTORES 31
O si la región es múltiplemente convexa.
Ó bien si la forma es:
Ejemplo: Comprobemos el teorema de Green en el plano para la integral:
, donde C es la curva cerrada que limita la región definida
por las curvas y = 2x e y = .
Estas curvas se cortan en los puntos P1(0,0) y P2(4,8)
.A lo largo de
= =
ANTONIO TAPIA SILVA VECTORES 32
La integral a lo largo x = 2x desde P2(0,4) a P1(0,0)
Es decir que la integral curva línea pedida es:
141,22 – 110,60 = 23,46
Por otro lado 1 = xy + y2
N = x2
Entonces:
EL TEOREMA DE GREEN EN FORMA VECTORIAL (TEOREMA DE
STOKES)
Si hacemos Mdx + Ndy = (Mi + Nj).(dxi + dyj)
Y si A = Mi + Nj
ANTONIO TAPIA SILVA VECTORES 33
r = xi + yj
Entonces Mdx + Ndy = [Link]
Y si tomamos el A, tendremos que:
y si multiplicamos por la normal al plano igual este caso es n = k.
Y podemos escribir el teorema de Green, en el plano como:
Con da = dxdy
Esta es la ley de Stokes (dada en la página 31).
En la expresión
Si A es una fuerza o un campo de fuerzas en el que se halla una partícula,
entonces la integral representa el trabajo para desplazar la partícula a lo largo de la
curva cerrada C.
CAPÍTULO C
COORDENADAS CURVILÍNEAS
ANTONIO TAPIA SILVA VECTORES 34
TRANSFORMACIÓN DE COORDENADAS.
Consideremos las coordenadas rectangulares de un punto P(x,y,z), expresadas en
función de las variables , en la forma:
... (C.1)
Y si despajamos
... (C.2)
Suponemos que las funciones de (C.1) y de (C.2) son uniformes y derivables, con
derivadas continuas, de tal manera que la correspondencia en las ternas (x,y,z) y
es recíproca (biunívoca)
Dado un conjunto P de coordenadas , se le puede asociar, según (C.2),
un conjunto único de números , a los que llamamos coordenadas
curvilíneas en P. Los sistemas de las ecuaciones (C.1) y (C.2) definen las
transformaciones de coordenadas.
C.1 COORDENADAS CURVILÍNEAS ORTOGONALES
Sean las superficies cuyas ecuaciones son , , , donde
son constantes, como se muestran en la figura. A estas se les denominan
“superficies coordenadas”.
La intersección de cada par de superficies define las líneas coordenadas. Si se
cortan en ángulo recto, entonces el sistema curvilíneo es ortogonal. Existe una
analogía entre las líneas con los ejes .
C.2 VECTORES UNITARIOS Y FACTORES DE ESCALA
ANTONIO TAPIA SILVA VECTORES 35
Sea un vector de posición en el punto P, lo podemos expresar en
la forma:
… (Según (C.1))
El vector tangente en P a la línea (para la cual son constantes) es
.
Por consiguiente, el vector unitario en dirección de la tangente es
… (C.3)
Hagamos
Donde es el módulo del vector tangente en el punto P.
En forma análoga
Las magnitudes se denominan “FACTORES DE ESCALA” donde
; ; … (C.4)
El sentido de los vectores es el del crecimiento de .
Los vectores unitarios normales a las superficies. Como , es un vector
normal en P a la superficie , el vector unitario es esta dirección y sentido
viene dado por:
, … (C.5)
Y en forma análoga en las otras dos superficies sus versores normales son:
Se puede expresar un vector un vector en función de los vectores unitarios
o de , como:
ANTONIO TAPIA SILVA VECTORES 36
Todo vector también se puede
expresar en función de los vectores
; o de
: A los que también
se les denomina vectores unitarios
en la base. Entonces
;o
bien
A las componentes , que son los módulos de las componentes del
vector en dirección de los vectores unitarios tangentes a la curva, (juegan el papel
de en el sistema cartesiano). Se denominan “COMPONENTES
COVARIANTES” y a las componentes , coeficientes de los gradientes y
que señalan en dirección perpendicular, se denominan “CONTRAVARIANTES”
del vector. Existe una similitud con los vectores unitarios cartesianos i, j, k. En
las coordenadas curvas, la dirección de los vectores cambia de un punto a otro.
Fig.C.2. Se representan los vectores tangentes a las líneas y normales a las superficies.
El conjunto de vectores COVARIANTES y los CONTRAVARIANTES, forman
dos sistemas de vectores recíprocos.
C.3 EL ELEMENTO DE LÍNEA
Partiremos de la derivada total de un vector
ANTONIO TAPIA SILVA VECTORES 37
Por lo que la expresión
… (C.6)
es el elemento de línea.
C.4 LOS DIFERENCIALES DE ÁREA
Sabemos que el área que forman dos vectores A y B se da por y su
dirección es perpendicular al plano que forman A y B. Los vectores unitarios
, determinan los productos escalar y vectorial respectivamente; usamos
sus propiedades:
1, si
=
0 si
En sentido antihorario es positivo y en sentido horario es negativo. Por
consiguiente, los diferenciales de área son:
=
Por consiguiente, el diferencial de área en coordenadas curvilíneas es:
… (C.7)
C.5 El elemento de volumen en coordenadas curvilíneas:
El elemento de volumen viene dado por
V= A.(BxC)
Por consiguiente:
ANTONIO TAPIA SILVA VECTORES 38
… (C.8)
Apliquemos de manera inmediata a coordenadas esféricas y cilíndricas.
C.6 COORDENADAS ESFÉRICAS.
Es aplicable en los problemas que tienen simetría esférica.
Un punto P del espacio tridimensional se representa en coordenadas esféricas por
medio de tres números reales (r, , ); entonces:
La correspondencia es biunívoca, por consiguiente,
El vector r en cartesianas es
Y su elemento de línea es
El vector posición en función de las coordenadas esféricas y los vectores
unitarios es
… (C.9)
ANTONIO TAPIA SILVA VECTORES 39
Determinemos sus factores de escala:
=
=
Resumiendo:
Los vectores unitarios son entonces:
Pero
Entonces
En forma análoga con y resumiendo, los vectores unitarios están dados
por
… (C.10)
ANTONIO TAPIA SILVA VECTORES 40
Despejemos i. j y k de este conjunto de ecuaciones para hallar los vectores
unitarios recíprocos.
Hallemos el determinante
=1
er Sen Sen Cos
i= e Cos Sen -Sen
e Cos 0
i= er Sen Cos - e (-cos cos )+ e (
Acomodando, y en forma similar para los demás, tendremos que los recíprocos
son:
… (C.11)
El elemento de línea. O “desplazamiento elemental”; se da por la derivada total
por la Ec. (C.6)
Y reemplazando sus valores
… (C.12)
Los diferenciales de área.
Usamos la ecuación (C.7)
Por lo que el vector diferencial de área es:
+ +
ANTONIO TAPIA SILVA VECTORES 41
y reemplazando los valores de los factores de escala
… (C.13)
Y las componentes de los diferenciales son:
… (C.14)
El elemento de volumen
Se da por
(C.15) …
Ubiquémoslos en el sistema cartesiano:
Fig.C.5
C.7 CORDENADAS CILÍNDRICAS
Veamos la figura:
ANTONIO TAPIA SILVA VECTORES 42
Fig. C.6
El vector posición r
Entonces
(u1, u2, u3)
Hallemos los factores de escala
Resumiendo:
Determinemos los vectores unitarios:
ANTONIO TAPIA SILVA VECTORES 43
(C.16)
Acomodando:
(C.17)
Su determinante
Los vectores unitarios i, j, k son
Sen 0
Cos 0 =
0 1
Hallando de manera similar para los vectores unitarios j y k y resumiendo:
... (C.18)
El elemento de desplazamiento:
dr = (1) de + d e + dz ez
… (C.19)
El elemento de área (o los elementos de área)
Tendremos que
=
ANTONIO TAPIA SILVA VECTORES 44
Entonces
… (C.20)
El diferencial de volumen es:
dV = h1h2h3 dV1 dV2 dV3
= (1) (1) d d dz
… (C.21)
En la figura se nota con claridad los elementos de área y de volumen.
Fig.C.7
C.8 Reciprocidad en la transformación de vectores en los tres sistemas de
coordenadas.
Entre coordenadas cartesianas y esféricas. Tengamos el vector
ANTONIO TAPIA SILVA VECTORES 45
Usando las relaciones de (C.11) y efectuando los reemplazos respectivos en los
versores de las coordenadas rectangulares del vector y
factorizando los versores esféricos
Podemos expresarlo por la matriz
… (C.22)
Pero el determinante es la UNIDAD, por lo que podemos expresar su matriz
transpuesta.
… (C.23)
Entre coordenadas cilíndricas y rectangulares. De manera similar, tengamos el
vector
Aplicando a las fórmulas (C.16) y (C.18).
Como matriz
ANTONIO TAPIA SILVA VECTORES 46
… (C.24)
Y su recíproco
… (C.25)
Entre cilíndricas y esféricas.
Sea el vector
Veamos los vectores unitarios del vector
Factorizando y y simplificando, encontramos que
Para el vector
De donde
De donde
Coinciden los versores en ambas bases de vectores.
Para
Resumiendo, los versores de ambos sistemas se relacionan por
ANTONIO TAPIA SILVA VECTORES 47
De forma matricial
… (C.26)
Su transpuesta
… (C.27)
Para la reciprocidad entre los dos vectores, se induce que
… (C.28)
Y su inversa
… (C.29)
C.9. GRADIENTE EN COORDENADAS CURVILINEAS
Consideremos el gradiente de la función f : f = f1 e1+ f2 e2 +f3 e3 ...
(C.30)
Siendo f1, f2 y f3 coeficientes por determinar.
Por el concepto de la derivada direccional, sabemos que
ANTONIO TAPIA SILVA VECTORES 48
Y también ya conocemos que el vector desplazamiento en curvilíneas es:
Multiplicando (C.30) por el vector desplazamiento y reemplazando en df
… (C.31)
Además, sabemos que la derivada total de la función f es
Igualando df
Igualemos los coeficientes con
Y reemplazando las 3 ecuaciones en (C.30) se tiene el operador NABLA en el
gradiente de coordenadas curvilíneas y se expresa por
… (C.32)
Esta expresión viene a ser El Gradiente En Coordenadas Curvilíneas.
Reemplacemos en las coordenadas cilíndricas y esféricas:
En coordenadas cilíndricas
Es decir:
… (C.33)
En coordenadas esféricas. El gradiente se da por
ANTONIO TAPIA SILVA VECTORES 49
O sea
… (C.34)
C.10. DIVERGENCIA EN COORDENADAS CURVILÍNEAS
La operación de divergencia no es la simple multiplicación escalar entre el
operador nabla y el vector, debido a que los versores son función de las
coordenadas, por lo que las derivadas de los vectores unitarios tienen valores
diferentes de cero. Es preferible usar la definición generalizada de la divergencia,
independiente del sistema de coordenadas, partimos de la definición de
divergencia que es el flujo vectorial por unidad de volumen y est dada por
ecuación (B.8).
Entonces
Tomemos un elemento V de volumen
Fig.C.8
El vector
Y el vector diferencial de área
Multipliquemos escalarmente en cada cara. El flujo es entonces
+ + + +
+ … (C.35)
El sentido de la normal en las integrales con prima es negativo. Teniendo en
cuenta la ubicación de las superficies constantes ui y teniendo en cuenta que el
vector en cada superficie prima es función de su posición. Podemos considerar:
ANTONIO TAPIA SILVA VECTORES 50
Por consiguiente, tendremos para el flujo
= +
+ … (C.36)
Como las dimensiones del elemento de volumen son infinitésimas, entonces la
variación en la integral es muy pequeña; por lo que la integral es casi constante.
Entonces podemos escribirlo como una buena aproximación:
+
+
+ … (C.37)
Dividiendo por el elemento de volumen para hallar la divergencia tenemos
= +
+ … (C.38)
Pero
De manera similar para y para .
Reemplazando (C.38)
… (C.39)
ANTONIO TAPIA SILVA VECTORES 51
En coordenadas esféricas
= … (C.40)
En coordenadas cilíndricas
… (C.41)
C.11 ROTACIONAL Y TEOREMA DE STOKES
Tomemos el rotacional a partir de la circulación, como la integral de línea de un
vector alrededor de una trayectoria cerrada.
La circulación a lo largo de una línea se da por
C= A. dl … (C.42)
Donde, si C es el trabajo, entonces A es la fuerza. Supongamos que la figura C.9
nos representa un contorno muy pequeño. Calculemos (C.32) para el contorno
infinitesimal de la figura. El vector A es:
Y el vector elemento de línea
Entonces:
… (C.43)
Podemos hallar la circulación en la superficie S 1 (dirección perpendicular al
elemento de volumen) y le llamaremos C1.
=
ANTONIO TAPIA SILVA VECTORES 52
Agrupando
Y con los argumentos del tema anterior, tendremos en el límite y acomodando:
Reemplazando en la circulación:
C1 =
Dividiendo entre el elemento de área , tendremos el teorema de
Stockes. La circulación por unidad de área es:
Que es la componente del rotacional en la dirección de e1. Es decir
= … (C.44)
En forma análoga para las demás componentes tendremos que la circulación por
unidad de elemento de área será, entonces:
… (C.45)
Y lo podemos expresar como un determinante, mediante:
ANTONIO TAPIA SILVA VECTORES 53
= … (C.46)
A la vez estamos comprobando que se cumple para coordenadas curvilíneas la
definición del rotacional dada por
Determinemos el rotacional en coordenadas esféricas y cilíndricas.
En coordenadas esféricas:
Ya sabemos que sus coordenadas curvilíneas y sus factores de escala,
respectivamente son
Por consiguiente:
… (C.47).
En coordenadas cilíndricas
Remplazando en la fórmula del rotacional en coordenadas curvilíneas tendremos
que:
… (C.48)
Es decir
… (C.49)
ANTONIO TAPIA SILVA VECTORES 54
C.12. EL LAPLACIANO EN COORDENADAS CURVILÍNEAS
Partimos de la definición de Laplaciano
… (C.50)
y Sabemos que el gradiente y la divergencia se dan, respectivamente por
Donde se ha puesto
C.8
Reemplazando, tenemos que el laplaciano es
… (C.51)
En coordenadas cilíndricas
= … (C.52)
En coordenadas esféricas
Tendremos que
= … (C.53).
ANTONIO TAPIA SILVA VECTORES 55
UN EJEMPLO DE APLICACIÓN
Supongamos que se desea hallar la circulación del vector fuerza dado por el vector
, en la trayectoria indicada
Permitámonos convertir a coordenadas cilíndricas dado que la circulación es en
Z=0.
De donde
El vector
La componente en z no es necesaria. El elemento de desplazamiento
La circulación es
Podemos utilizar el rotacional en coordenadas cilíndricas:
Y el teorema de Stokes
Es obvio que debemos de tomar solamente la componente en Z ya que el elemento
de área tiene la dirección perpendicular al plano donde circula.
ANTONIO TAPIA SILVA VECTORES 56
COMPONENTES CONTRAVARIANTES DE UN VECTOR
Veamos para el caso de dos sistemas de coordenadas curvilíneas
Sea el vector dado en dos dimensiones de coordenadas curvilíneas y
, hallemos la relación contravariante del vector en ambos sistemas.
ANTONIO TAPIA SILVA VECTORES 57
En cartesiana, tenemos que = + +
La transformación de coordenadas se da por
1.a
Y en el otro sistema:
.
1.b
z
Busquemos una forma de transformación directa de a
Entonces tendremos
(2)
Por la regla de la cadena
=
es decir que
= (3)
De (2)
Reemplazando en (3)
ANTONIO TAPIA SILVA VECTORES 58
+ +
+ =
= + +
Igualando coeficientes:
(4)
Podemos expresar el vector en ambos sistemas de coordenadas, como:
(5)
Los coeficientes y son las componentes contravariantes del
vector en ambos sistemas.
Entonces
=
Entonces
+ +
+ =
+ +
Por consiguiente:
ANTONIO TAPIA SILVA VECTORES 59
(6)
Podemos simplificar
(8) , con p=1,2,3
Ó más aún
p=1,2,3 (9)
COMPONENTES COVARIANTES EN LOS DOS SISTEMAS DE
COORDENADAS
Representemos las componentes covariantes en los sistemas y
, por y , respectivamente.
Hagamos
(1)
Y si
(2)
Con p = 1, 2, 3.
Además:
= +
+ (3)
ANTONIO TAPIA SILVA VECTORES 60
+ (4)
Igualando los coeficientes de los vectores unitarios , en (3) y (4), tenemos:
= (5)
En analogía, reemplazando como en el caso de los contravariantes:
(6)
Definamos:
“Si tres magnitudes , en un sistema de coordenadas , están
relacionadas con otras tres en otro sistema de coordenadas , mediante las
transformaciones de la ecuaciones (6), (7), (8) o (9), entonces dichas magnitudes
se denominan componentes de un vector contravariante, o tensor contravariante
de primer orden”.
ANTONIO TAPIA SILVA VECTORES 61
ANÁLISIS TENSORIAL
Espacios en N dimensiones. Ampliemos el caso de 3 dimensiones (x,y,z), (
a espacio de N dimensiones.
Un punto en N dimensiones se caracteriza por un conjunto de números a los que
nombramos por , en donde m1, 2, 3,…N no son exponentes y por
supuesto que no se pueden ser representados gráficamente.
TRANSFORMACIÓN DE COORDENADAS
Si ,y son las coordenadas de un mismo punto
en dos sistemas de referencia distintos.
Supongamos que existen N relaciones independientes entre las coordenadas
anteriores, de la forma.
= ,
= , (1)
ANTONIO TAPIA SILVA VECTORES 62
= .
.
. .
.
= .
Ó
(2)
Con k=1,2, 3, …, N
Se supone que todas las funciones son uniformes, continuas y con derivadas
también continuas.
Bajo estas condiciones a cada conjunto de coordenadas , le
corresponde un conjunto único , de manera que
Con k=1,2, 3, …, N
Las relaciones (2) ó (3) definen las fórmulas de transformación de un sistemas a
otro.
Convenio de los índices repetidos. Se deben a Einstein.
Podemos expresar:
como
, como
Todo índice que se repita en los términos que implica una sumación respecto de
él, ej: . A veces se le denomina SEUDO ÍNDICE O ÍNDICE UMBRAL.
Veamos el ejemplo:
Los términos de la suma
= =
Los términos de la suma
ANTONIO TAPIA SILVA VECTORES 63
= =
Supongamos que
es el elemento de línea
Y el cuadrado del elemento de línea se puede expresar por
= ,
donde
Esta es la fórmula cuadrática fundamental o forma simétrica. Los coeficientes
se denominan coeficientes métricos y tienen simetría.
En coordenadas ortogonales
con p=q
si
Es fácil darse cuenta que
, ,
Escribamos todos los términos de
, con N=1, 2, 3.
= + +
= + + +
+ + +
+ + + .
ANTONIO TAPIA SILVA VECTORES 64
TENSORES CONTRAVARIANTES, COVARIANTES Y MIXTOS
Supongamos que magnitudes en un sistema de coordenadas
se relacionan con algún otro sistema , mediante las leyes de
transformación:
con p,r = 1, 2, 3, …,N
Y según el convenio de la sumación de Einstein (de los índices repetidos)
denominamos las componentes contravariantes de un tensor de segundo orden
Se denominan componentes covariantes de un tensor de segundo orden a las
magnitudes si se verifica que
Y Se denominan componentes de un tensor mixto de segundo orden a las
magnitudes si se verifica que
Como ejemplo escribamos la ley de transformación del tensor .
En efecto
La Ley de transformación del tensor .
La Ley de transformación del tensor =
EJEMPLO: las componentes de un tensor covariante en coordenadas cartesianas
son . Hallemos sus componentes covariantes en coordenadas
esféricas.
En efecto, tendremos que las componentes covariantes en rectangulares son:
ANTONIO TAPIA SILVA VECTORES 65
(Ojo que son súper índices, no exponentes)
Además
En coordenadas esféricas las componentes covariantes son:
y la transformación
con las fórmulas de transformación
En asterisco:
+ +
+ +
= + +
= + +
= + +
= + + :
+ +
= +
= +
ANTONIO TAPIA SILVA VECTORES 66
= +
+ +
+ +
= +
= + -
Son 9 términos.
DELTA DE KRONECKER
Se representa por
cómo se tiene en el producto escalar de dos vectores
OPERACIONES CON TENSORES
Suma (o diferencia) de dos tensores del mismo orden y el mismo tipo; es decir
que tienen el mismo número de índices covariantes y contravariantes, su suma es
otro tensor del mismo orden y tipo.
Multiplicación externa (producto externo). El producto externo de dos tensores
es otro tensor cuyo orden es la suma de los órdenes de los tensores dados.
Es decir
=
no es conmutativo.
Contracción. Si en un tensor se igualan un índice contravariante a otro
covariante, entonces debe de sumarse respecto de dicho índice. La suma es un
tensor de orden inferior en dos unidades respecto al tensor de origen. Por ejemplo,
si en el tensor , que es de orden 5, si r=m, entonces se obtiene que = ,
que es otro tensor de orden 3. Si en este último tensor se tiene que n=q , entonces
se obtiene = , que es un tensor de orden 1. Es decir, un vector.
Multiplicación interna. Resulta del proceso de multiplicación externa de dos
tensores y luego de una contracción, igualando un índice del primer factor a un
índice del segundo factor.
ANTONIO TAPIA SILVA VECTORES 67
Ejemplo, sean los tensores y , entonces su producto externo es ,
igualando q=r, p=s se obtiene .
Ejemplos de aplicación
-expresemos la velocidad de una partícula en forma tensorial.
Supongamos que la ley de desplazamiento sobre una trayectoria curvilínea es
entonces su tensor es
-Expresemos la aceleración de una partícula en forma tensorial.
Solamente lo tomamos como una derivada de la velocidad:
La fuerza
FÓRMULAS DEL GRADIENTE ROTACIONAL, DIVERGENCIA
ROTACIONAL Y LAPLACIANO COMO TENSORES
GRADIENTE. Si es un escalar, entones el gradiente de este escalar se expresa
como
p =
p representa la derivada covariante respecto de xp
DIVERGENCIA. La divergencia de es la contracción de su derivada
covariante con respecto de . Es decir, la contracción de .
Entonces
= =
ROTACIONAL. Se define como el tensor de segundo orden, dado por
LAPLACIANO. Lo definimos como
= =
=
Es preferible tomar el valor absoluto de
ANTONIO TAPIA SILVA VECTORES 68
Apliquemos a la expansión multipolar del potencial electrostático
La expresión para el potencial electrostático para una distribución continua de
cargas se da mediante la expresión:
Donde puede tener una distribución lineal, superficial o volumétrica.
Si la evaluación de la integral es difícil, esa preferible tomar una aproximación de
la integral, expandiendo el denominador según la serie del binomio. En este caso
nos estamos refiriendo a que podemos expandir el potencial de la forma:
Donde lo asociamos con el potencial debido a la carga neta de la
distribución y lo denominamos potencial monopolar. lo asociamos con el
potencial debido al momento dipolar de la distribución de carga y por supuesto
que se denomina potencial dipolar; al con el momento cuadrupolar; el
siguiente el octopolar; las subsiguientes las denominamos de ORDEN
SUPERIOR.
Si la distribución de carga cuyo potencial se desea evaluar se extiende en una
región finita (para lo cual podemos elegir el sistema de coordenadas dentro o fuera
de la distribución de la carga), evaluaremos el potencial en puntos alejados de la
distribución; es decir en puntos >>
Expandiendo el denominador:
Y si
ANTONIO TAPIA SILVA VECTORES 69
Con
Tomemos la serie en segunda aproximación con
Tendremos que:
Agrupando términos semejantes que multipliquen a , etc. y
multiplicando por , se tiene:
Reemplazando esta última expresión remarcada con el doble asterisco en la
expresión de asterisco
Es decir, la suma de los potenciales
Donde
Es el potencial debido al momento monopolar (carga neta) de la distribución de
carga:
La cual puede considerarse como una carga puntual en el origen de coordenadas.
Desde luego que el diferencial de carga prima se puede expresar en función de las
densidades lineal, o superficial o volumétrica, según la circunstancia. Esto es:
, o bien ; ó
El término
ANTONIO TAPIA SILVA VECTORES 70
=
Es el potencial debido a un momento dipolar de la distribución de carga. Es
decir
, ó , ó
Al cual se considera como un dipolo puntual situado en el origen de coordenadas
El término
(*)
es el MOMENTO CUADRUPOLAR.
Antes de interpretarlo físicamente, tengamos en cuenta que:
=
Es decir:
En total son nueve términos, que se pueden escribir en forma matricial como:
Por otro lado:
Que puede ser expresado en forma matricial de la siguiente forma:
y
z
ANTONIO TAPIA SILVA VECTORES 71
Y la expresión:
Reemplazando esta expresión matricial en (*) y como la integración es sobre la
distribución de carga; es decir sobre los puntos , en los términos
que dependen de pueden salir de la integral como un factor común.
Entonces tendremos:
Siendo:
es una expresión tensorial llamada MOMENTO CUADRUPOLAR
ELÉCTRICO.
Si observamos que el momento cuadrupolar
Donde se ha puesto:
ANTONIO TAPIA SILVA VECTORES 72
Nos basta con evaluar 5 integrales, ya que
Entonces nos es suficiente evaluar dos términos de la diagonal y el siguiente se
puede obtener en forma algebraica
Si la distribución de carga fuera discreta, para evaluar debemos de sustituir la
integral por la sumatoria por .
La ecuación del potencial se puede escribir como:
Apliquemos al anillo de densidad uniforme:
Sea el anillo de densidad uniforme y radio . Calculemos el potencial eléctrico
en cualquier punto del espacio, en una aproximación hasta el orden cuadrupolar.
Hallemos también para puntos sobre el eje Z perpendicular al eje del anillo
>>
ANTONIO TAPIA SILVA VECTORES 73
Hallemos
Entonces el primer término
Calculemos
Tenemos que:
Con
=
=
entonces
Y para :
Con
Y en coordenadas cilíndricas:
Entonces hallemos :
=0
Hallemos
Con
ANTONIO TAPIA SILVA VECTORES 74
Entonces:
Pero tenemos que:
Por lo que
=
=
Entonces:
Donde
Por lo tanto
Es decir:
Calculemos :
ANTONIO TAPIA SILVA VECTORES 75
Y con
Reemplacemos:
Donde:
Por lo que
Entonces
Podemos calcular ; para lo cual tenemos que:
Calculemos :
ANTONIO TAPIA SILVA VECTORES 76
Calculemos
Pero
Entonces
En forma similar para
Entonces
Entonces estamos en condición de hallar el TENSOR MOMENTO
CUADRUPOLAR
Entonces
ANTONIO TAPIA SILVA VECTORES 77
Entonces
En coordenadas esféricas:
En el paréntesis:
Por consiguiente:
Por lo tanto, el potencial en coordenadas esféricas es
Si el punto está sobre el eje Z:
ANTONIO TAPIA SILVA VECTORES 78
En este caso
Por consiguiente
Este mismo resultado se pudo haber obtenido de
Expandiendo:
ANTONIO TAPIA SILVA VECTORES 79