2.

5 CONTRACCIÓN Y MULTIPLICACIÓN
Cuando los dos índices de un tensor se igualan, y se realiza una suma sobre esto índice repetido, el
proceso se llama contracción. Un ejemplo es:

Que es la suma de los términos diagonales de Aij. La suma Ajj

Es un escalar y es independiente del sistema coordinado. En otras palabras, Ajjes un invariante. (Hay
tres invariantes independientes de un tensor de segundo orden, y Ajjes uno de ellos; vea el ejercicio
2.11.) Los tensores de orden superior pueden formarse multiplicando los tensores de orden inferior. Si
A y B son dos tensores de segundo orden, luego los 81 números definidos por Pijkl≡ AijBkl transformar
de acuerdo a (2.13), y por lo tanto forman un tensor de cuarto orden.
Los tensores de orden inferior se pueden obtener realizando una contracción dentro de un multiplicado
formar. Las cuatro posibles contracciones de AijBklson:

Donde todas las definiciones de igualdad se desprende de la convención de suma. Los cuatro productos
en (2.14) son tensores de segundo orden. Nota también en (2.14) cómo se han reorganizado los
términos hasta que el índice sumado sea adyacente; en este punto se pueden escribir como un producto
de matrices el producto contratado de un tensor de segundo orden. A y un vector u es un vector Los dos
posibilidades son:

Donde de nuevo la definición de igualdad se desprende de la convención de suma. El doble producto

contraído de dos tensores de segundo orden. A y B es un escalar Usando las tres notaciones, la dos
posibilidades son:

Donde el colon en negrita (:) implica una doble producto de contracción o punto doble.
Ejemplo 2.5
𝜕
La fuerza superficial Fj por unidad de volumen en un fluido es la derivada vectorial, del tensor de
𝜕𝑥𝑖
esfuerzo Tij. Determinar los tres componentes del vector Fj.
Solución
𝜕
Aquí el vector derivado, , debe aparecer a la izquierda del tensor de tensión Tij. Así, para Fj ser
𝜕𝑥𝑖
expresado en términos de multiplicación de matrices, debe escribirse como un vector de fila
multiplicado por el matriz de esfuerzo.

Desde el segundo índice de Tijespecifica la dirección del componente de tensión, los tres componentes
de fuerza. F1, F2 y F3 Solo dependen de los componentes de esfuerzo en las
‘1’, ‘2’ y ‘3’ direcciones, respectivamente, inclusoaunque se toman derivados en las tres direcciones.

2.6 FUERZA EN UNA SUPERFICIE

Un elemento de área de superficie tiene un tamaño (o magnitud) y una orientación, por lo que puede
tratarse como un vector dA.SidAes el tamaño del elemento de superficie, y n es su vector unidad
normal, entoncesdA=ndA.

FIGURA 2.5 Fuerza F por unidad de área en un elemento de superficie cuya normal exterior es n, las
áreas de las caras del tetraedro son perpendiculares al i-ésimo eje de coordenadas son dAi, El área de la
cara más grande del tetraedro es dA.Como en
Figura 2.4, se muestran las direcciones de tensiones normales y cortantes.

Supongamos que los nueve componentes, Tij, del tensor de esfuerzo con respecto a un conjunto dado
de coordenadas cartesianas O123, y queremos encontrar la fuerza por unidad de área, f(n) con
componentes fi, en un elemento superficial arbitrariamente orientado con normal n (ver Figura 2.5).
Una forma de completar esta tarea es cambiar a un sistema de coordenadas rotado y usar(2.12) para
encontrar las tensiones normales y cortantes en el elemento de superficie. Se describe un método
alternativo aquí. Considere el elemento tetraédrico mostrado en
Figura 2.5. La fuerza neta f1en el elemento en la primera dirección producida por las tensiones Tijes:
La geometría del tetraedro requiere , donde ni son los componentes del vector normal de
superficie n. Así, la ecuación de fuerza neta se puede reescribir:

Dividiendo por dAentonces produce (con la suma implícita), o por algún componente de f.

Donde la versión solo en negrita de (2.15) es inequívoco cuando, Tij=Tji, una afirmación que es
demostrado en el Capítulo 4. Por lo tanto, el producto contraído o interno del tensor de esfuerzo T y el
vector normal de la unidad n da la fuerza por unidad de área en una superficie perpendicular a n. Este
resultado es análogo a 𝑢𝑛 = 𝑢 ∗ 𝑛, dónde unes el componente del vector u a lo largo n; sin embargo,
mientras un es un escalar, F en (2.15) es un vector.

2.7 DELTA DE KRONECKER Y TENSOR ALTERNO

El delta de Kronecker se define como:

En tres dimensiones espaciales es la matriz identidad 3x3.

En las operaciones de multiplicación de matrices que involucran el delta de Kronecker, simplemente

reemplaza su
Índice sumado por su otro índice. Considerar:

En la derecha está u1 cuando i=1, u2 cuando i=2 y u3 cuando i=3; entonces:

Se puede demostrar que δij es un tensor isotrópico en el sentido de que sus componentes no han
cambiado por una rotación del marco de referencia, es decir . Tensores isotrópicos puede ser
de varios órdenes. No hay tensor isotrópico de primer orden, y δij es el único tensor isotrópico de
segundo orden. También hay un solo tensor isotrópico de tercer orden. Se llama el tensor alterno o
símbolo de permutación, y es definido como:

De esta definición, está claro que un índice enεijk puede ser movido dos lugares (ya sea a la derecha o
para la izquierda) sin cambiar su valor. Por ejemplo, εijk=εjki dónde i se ha movido dos lugares para la
derecha, yεijk=εkji dónde k se ha movido dos lugares a la izquierda. Para un movimiento de un lugar, sin
embargo, el signo se invierte. Por ejemplo, εijk=-εjki dónde j se ha movido un lugar a la derecha.
Una relación muy utilizada es la relación delta epsilon:

La validez de esta relación puede ser verificada eligiendo algunos valores para los índices ijlm.
Esta relación puede recordarse observando los siguientes dos puntos: 1) El índice adyacente k se suma
y 2) Los primeros dos índices en el lado derecho, a saber, i y l son los primeros índices de εijk y el
primer índice libre de εklm. Los índices restantes en el lado derecho siguen inmediatamente.

2.8 PROUCTO PUNTO Y PRODUCTO CRUZ.

El producto punto de dos vectores u y v está definido en (2.2). También se puede escribir u*v=uvcos θ,
dónde u y v son las magnitudes de los vectores y θ es el ángulo entre los vectores (ver ejercicios 2.14 y
2.15). El producto punto es, por lo tanto, la magnitud de un vector multiplicado por el componente del
otro en la dirección de la primera. El producto punto u*v es igual a la suma de los términos diagonales
del tensor uivjy es cero cuando u y v son ortogonales
El producto cruzado entre dos vectores u y v es definido como el vector w cuya magnitud esuv sin θ
dónde θ es el ángulo entre u y v, y cuya dirección es perpendicular al plano de u y v tal que u, v y w
forman un sistema derecho. En este caso, uxv=-vxu. Además, los vectores unitarios en sistemas de
coordenadas diestros obedecen la regla cíclica e1xe2=e3. Estos requisitos son suficientepara
determinar:

(Ver ejercicio 2.16). Ecuación (2.20) Se puede escribir como el determinante de una matriz:

En notación indicial, el componente k de uxvSe puede escribir como:

Como revisión, para k=1 los términos distintos de cero en la suma doble en (2.21) es el resultado de
i=2, j=3 y de i=3, j=2. Esto se desprende de la definición (2.18) que el símbolo de permutación es cero
si alguno de los dos índices son iguales. Así, (2.21) da:

Que está de acuerdo con(2.20). Tenga en cuenta que la cuarta forma de (2.21) se obtiene a partir del
segundo por moviendo el índice k dos lugares a la izquierda; ver el comentario siguiente (2.18).

2.9 GRADIENTE, DIVERGENCIA Y CURVA

El operador de diferenciación vectorial “delta” es definida simbólicamente por:

Cuando se opera en una función escalar de posición ϕ genera el vector:

 𝛿𝜙
Cuyo componente i es . El vector se llama el gradiente ϕ y es perpendicular a superficies
𝛿𝑥𝑖
definida por ϕ=constante. Además, se especifica es la magnitud y dirección de la máxima tasa espacial
de cambio de ϕ (Figura 2.7).

FIGURA 2.7 Una ilustración del gradiente , de una función escalar ϕ. Las curvas de constante ϕ y
son perpendiculares, y la derivada espacial de ϕ en la dirección n es dado por n* . El cambio
más rápido en Φ es encontrado cuando n y son paralelos

La tasa espacial de cambio de ϕ en cualquier otra dirección n es dado por:

En coordenadas cartesianas, la divergencia de un vector de campo u es definida como el escalar.

Hasta ahora, tenemos definidas las operaciones de la pendiente de un escalar y la divergencia de una
vector. Podemos, sin embargo, generalizar estas operaciones. Por ejemplo, la divergencia , de un
tensor de segundo orden T puede ser definida como el vector cuyo componente j es:
La operación de divergencia disminuye el orden del tensor por uno. Por el contrario, el gradiente
operación aumenta el orden de un tensor por uno, cambiando un tensor de orden cero a un tensor de
𝛿𝑢
primer orden, y un tensor de primer orden a un tensor de segundo orden, es decir, i.e. 𝑖 . La curva de
𝛿𝑥_𝑗
un vector de campo u está definida como el vector cuya componente i se puede escribir como:

Utilizando (2.21) y (2.22). Los tres componentes del vector se pueden encontrar fácilmente
desde el lado derecho de (2.24). Para el componente i=1, los términos distintos de cero en el resultado
de doble suma desde j=2, k=3 y desde j=3, k=2. Los tres componentes de son finalmente
encontrados como:

Un vector de campo u se llama solenoide o libre de divergencias si e irrotacional o sin curva

si La palabra solenoidal se refiere al hecho de que la divergencia de la inducción
magnética siempre es cero debido a la ausencia de monopolosmagnéticos. La razón de la palabra.
La irrotación se aclara en el capítulo 3.

2.10 TENSORES SIMÉTRICOS Y ANTISIMÉTRICOS

Un tensor B se llama simétrico en los índices i y j si los componentes no cambian cuando i y j se
intercambian, es decir, si Bij = Bji. Por lo tanto, la matriz de un tensor de segundo orden simétrico se
compone de solo seis componentes distintos (los tres en la diagonal donde i = j, y los tres arriba o
debajo de la diagonal donde i s j). Por otro lado, un tensor se llama antisimétrico si Bij = -Bji. Los
componentes diagonales de un tensor antisimétrico son cada uno cero, y solo tiene tres componentes
distintos (los tres arriba o abajo de la diagonal). Cualquier tensor puede representarse como la suma de
tensores simétricos y antisimétricos. Porque si escribimos:

entonces la operación de intercambiar i y j no cambia el primer término, pero cambia el signo del
segundo término. Por lo tanto, (Bij + Bji) / 2 = Sij se llama la parte simétrica de Bij, y (Bij - Bji)/2 = Aij
se llama la parte antisimétrica de Bij.
Cada vector puede asociarse con un tensor antisimétrico, y viceversa. Por ejemplo, podemos asociar el
vector 𝜔 que tiene componentes 𝜔𝑖 , con un tensor antisimétrico:

Los dos están relacionados a través de:

Como comprobación, (2.27) da R11 = 0 y R12= -ε123𝜔3 = −𝜔3 , de acuerdo con (2.26). (En el Capítulo 3,
R se reconoce como el tensor de rotación correspondiente al vector de vorticidad 𝜔). Una operación
que ocurre comúnmente es el producto doblemente contraído, P, de un tensor T simétrico y otro tensor
B:

donde S y A son las partes simétricas y antisimétricas de B. La igualdad final se desprende de la

convención de suma de índices; las sumas se completan tanto en k como en l, por lo que estos índices
pueden reemplazarse por dos índices distintos. Intercambiar los índices de A en el término final de
(2.28) produce
P = TijSij - TijAji, pero esto también puede escribirse P = TjiSji - TjiAji porque Sij y Tij son simétricos.
Ahora, reemplace el índice j por k y el índice i por l para encontrar:

Esta relación y la cuarta parte de la igualdad extendida en (2.28) requieren que T ijAij =TklAkl = 0, y:

Por lo tanto, el producto doblemente contraído de un tensor T simétrico con cualquier tensor B es igual
a T contraído doblemente con la parte simétrica de B, y el producto doblemente contraído de un tensor
simétrico y un tensor antisimétrico es cero. El último resultado es análogo al hecho de que la integral
definida en un intervalo uniforme (simétrico) del producto de una función simétrica y antisimétrica es
cero.

2.11 VALOR PROPIO Y VECTORES PROPIOS DE UN TENSOR SIMÉTRICO

Se supone que el lector está familiarizado con los conceptos de valores propios y vectores propios de
una matriz, por lo que aquí solo se proporciona una breve revisión de los principales resultados.
Supongamos que 𝜏 es un tensor simétrico con elementos reales. Entonces se pueden probar los
siguientes hechos:

1. Hay tres valores propios reales 𝜆𝑘 = (𝑘 = 1,2,3), que pueden o no ser todos distintos. (Aquí, el
superíndice k no es un exponente, y 𝜆𝑘 no denota el componente de un vector.) Estos valores propios
(𝜆1 , 𝜆2 , y 𝜆3 ) son las raíces o soluciones del polinomio de tercer grado:

2. Los tres vectores propios 𝑏 𝑘 que corresponden a valores propios distintos son ortogonales entre sí.
Estos vectores propios definen las direcciones de los ejes principales de 𝜏. Cada b se encuentra
resolviendo tres ecuaciones algebraicas:

(𝑖 = 1, 2 o 3), donde el superíndice k en 𝜆 y b se ha omitido para mayor claridad porque no hay suma
sobre k.

3. Si el sistema de coordenadas se gira de modo que sus vectores unitarios coincidan con los vectores
propios, entonces 𝜏 es diagonal con los elementos 𝜆𝑘 en este sistema de coordenadas rotado:
4. Aunque los elementos 𝜏𝑖𝑗 cambian a medida que se gira el sistema de coordenadas, no pueden ser
más grandes que el mayor 𝜆 o más pequeño que el menor 𝜆; Los 𝜆𝑘 representan los valores extremos de
𝜏𝑖𝑗 .

EJEMPLO 2.11
El tensor de velocidad de deformación S está relacionado con el vector de velocidad u mediante:

Para un flujo bidimensional paralelo a la dirección , 𝑢 = (𝑢1 (𝑥2 ), 0), muestra cómo S está diagonalizada
en un marco de referencia girado para coincidir con los ejes principales.

Solución
Para el perfil de velocidad dado u1 (x2), es evidente que

. El tensor de velocidad de deformación en el sistema

de coordenadas original es por lo tanto:

Los valores propios se determinan a partir de:

que tiene soluciones . El primer vector propio b1 viene dado por:

que tiene solución , cuando b1 se normaliza para tener unidad de magnitud. El

segundo eigenvector se encuentra de manera similar para que:

Estos vectores propios se muestran en la Figura 2.8. La dirección de la matriz del coseno del original y
el sistema de coordenadas girado es, por lo tanto:

que representa la rotación del sistema de coordenadas por 45°. Usando la regla de transformación
(2.12), los componentes de S en el sistema rotado se encuentran de la siguiente manera:
(En lugar de usar (2.12), todos los componentes de S en el sistema girado se pueden encontrar
realizando el producto de matriz 𝐶 𝑇 . 𝑆. 𝐶.). La matriz de S en el marco girado es, por lo tanto:

FIGURA 2.8 Sistema de coordenadas original Ox1x2 y el sistema de coordenadas girado Ox’1x’2 que
tienen vectores unitarios que coinciden con los vectores propios del tensor de velocidad de
deformación en el Ejemplo 2.4. Aquí, la velocidad de deformación se determina a partir de un flujo
unidireccional que tiene solo una variación de la corriente transversal, y se determina que el ángulo de
rotación es de 45 °.

La matriz anterior contiene solamente términos diagonales. Para Γ positivo, se mostrará en el siguiente
capítulo que representa un estiramiento lineal a una velocidad Γ a lo largo de un eje principal, y una
compresión lineal a una tasa -Γ a lo largo del otro; las deformaciones de corte son cero en el sistema de
coordenadas del eje principal del tensor de velocidad de deformación.

2.12 TEOREMA DE GAUSS

Este teorema muy útil relaciona el volumen y las integrales de superficie. Sea V un volumen delimitado
por una superficie cerrada A. Considere un elemento de superficie infinitesimal dA que tiene una
unidad externa normal n con componentes ni (Figura 2.9), y sea Q(x) un campo escalar, vectorial o
tensorial de cualquier orden. El teorema de Gauss establece que:
La forma más común del teorema de Gauss es cuando Q es un vector, en cuyo caso el teorema es:

lo que comúnmente se llama el teorema de la divergencia. En palabras, el teorema establece que la

integral de volumen de la divergencia de Q es igual a la integral de superficie del flujo de salida de Q.

FIGURA 2.9 Ilustración del teorema de Gauss para un volumen V encerrado por el área de superficie A.
Se muestra un elemento de volumen pequeño, dV, y un elemento de área pequeña dA con n normal
hacia afuera.

Alternativamente, (2.30) define una derivada de campo generalizada, denotada por D, de Q cuando se
considera en su forma limitante para un volumen extremadamente pequeño:

Curiosamente, esta forma se especializa fácilmente en el gradiente, divergencia y curvatura de

cualquier escalar, vector o tensor Q. Además, considerando (2.31) como una definición, las recetas para
el cálculo de las derivadas de campos vectoriales se pueden obtener de cualquier sistema coordinado.
Como se indica, (2.31) define el gradiente de un tensor Q de cualquier orden. Para un tensor de orden
uno o superior, la divergencia y la curvatura se definen al incluir un producto puntual (escalar) o un
producto cruzado (vector), respectivamente, bajo la integral:
EJEMPLO 2.12 Obtenga la receta para la divergencia de un vector Q (x) en coordenadas cartesianas a
partir de la definición integral (2.32).

Solución

Considere un volumen rectangular elemental centrado en x con caras perpendiculares a los ejes de
coordenadas (consulte la Figura 2.4). Indique las longitudes de los lados paralelos a cada eje de
coordenadas mediante Δx1, Δx2 y Δx3, respectivamente. Hay seis caras en este volumen rectangular.
Primero considera los dos que son perpendiculares al eje x1, EADH con n = e1 y FBCG con n=-e1. Una
expansión de Taylor de Q (x) desde el centro del volumen al centro de cada uno de estos lados produce:

de modo que la contribución de la dirección x a la integral de superficie en (2.32) es:

Del mismo modo para las otras dos direcciones:

Reunir las contribuciones de las seis caras (o las tres direcciones) para evaluar (2.32) produce:

y cuando se toma el límite, surge la forma de la divergencia en la coordenada cartesiana esperada:

2.13 TEOREMA DE STOKES

El teorema de Stokes relaciona la integral sobre una superficie abierta A con la integral alrededor de la
curva de delimitación de la superficie C. Aquí, a diferencia del teorema de Gauss, el interior y el exterior
de A no están bien definidos, por lo que se debe hacer una elección arbitraria para la dirección de la
normalidad externa n (aquí siempre se origina en el exterior de A). Una vez hecha esta elección, el
vector tangente unitario a C, t, apunta en sentido contrario a las agujas del reloj cuando se mira al
exterior de A. El vector unitario final, p, es perpendicular a la curva C y tangente a la superficie, por lo
que es perpendicular a n y t. Juntos, los tres vectores unitarios forman un sistema diestro: t× n = p (ver
Figura 2.10). Para esta geometría, el teorema de Stokes establece:

FIGURA 2.10 Ilustración del teorema de Stokes para la superficie A delimitada por la curva cerrada C. A
los efectos de definir vectores unitarios, se debe elegir el interior y el exterior de A, y aquí se ilustra una
de estas opciones. El vector unitario n es perpendicular a A y se origina desde el exterior de A. El vector
unitario, t, es localmente tangente a la curva C. El vector unitario p es perpendicular a C pero
localmente es tangente a la superficie A, por lo que es perpendicular tanto a n como a t. Los vectores
unitarios n, t y p definen una tríada de direcciones para diestros, t × n = p.

donde s es la longitud de arco de la curva cerrada C. Este teorema significa que la integral de superficie
de la curvatura de un campo vectorial u es igual a la integral de línea de u a lo largo de la curva de límite
de la superficie. En mecánica de fluidos, el lado derecho de (2.34) se denomina circulación de u
alrededor de C. Además, (2.34) se puede usar para definir el rizo de un vector a través del límite de la
circulación sobre una superficie infinitesimal como:

La ventaja de las definiciones integrales de las derivadas de campo es que tales definiciones no
dependen del sistema de coordenadas.

EJEMPLO 2.13
Obtenga la receta para el rizo de un vector u (x) en coordenadas cartesianas a partir de la definición
integral dada por (2.35).
Solución
Esta receta se obtiene considerando contornos rectangulares en tres planos perpendiculares que se
intersecan en el punto (x, y, z) como se muestra en la Figura 2.11. Primero, considere el rectángulo
elemental que se encuentra en un plano definido por x = const. El punto central en este plano es (x, y,
z) y el área del rectángulo es ∆y∆z. Se puede demostrar mediante una integración cuidadosa de una
expansión de Taylor del integrando que la integral a lo largo de cada segmento de línea puede ser
representada por el producto del integrando en el centro del segmento y la longitud del segmento con
atención a la dirección de integración ds. Así obtenemos:

FIGURA 2.11 Dibujo para el Ejemplo 2.11. El punto (x, y, z) está ubicado en el centro del cruciforme. Las
flechas muestran la dirección de ds en el borde del contorno rectangular que se encuentra en la
constante x=constante. El área dentro de este contorno es ∆y∆z.

Tomando los límites produce:

De manera similar, la integración alrededor de los rectángulos elementales en los otros dos planos
conduce a: