Fenómenos de transporte pa’l que

no sabe

Apuntes de clase, notas de estudio y

ejercicios

1
 Prefacio
Este libro está planeado para que el estudiante se interese en el estudio y aplicación de los
principios de los Fenómenos de Transporte para su beneficio personal y el futuro desarrollo
de su carrera profesional. El libro se presenta de la misma manera como se presenta el curso
de Fenómenos de Transporte para estudiantes de quinto semestre de la Escuela de Ingeniería
Química de la Universidad Industrial de Santander, así: cada unidad corresponde a la
transcripción de una clase del semestre y esta está seguida por una exposición más detallada
de los conceptos discutidos y con material didáctico adicional para ayudar al estudiante a
reflexionar sobre los conceptos presentados y a aplicarlos en la solución de problemas y
ejercicios relacionados con el tema tratado. Por otro lado, el enfoque argumentativo de todo
el libro es la presentación de la traducción de los principios de conservación de la masa,
cantidad de movimiento y energía en las denominadas ecuaciones de cambio. A los autores
nos parece que este libro puede ayudar muchísimo a los estudiantes que por primera vez
toman un curso de Fenómenos de Transporte, pero es un libro complementario a las primeras
dos ediciones del clásico libro de Fenómenos de Transporte de Bird, Stewart y Lightfoot y otros
excelentes libros tales como: Mecánica y Dinámica de Fluidos de Streeter, Mecánica de Fluidos
de White, entre otros.

2
 Unidad I
Preludio matemático
No es posible estudiar ni ser ingeniero sin tener aprecio por las matemáticas, pero tampoco se
necesita saber todo en matemáticas para apreciar la ciencia y la ingeniería.
Las siguientes clases, dos en total, presentan el marco de referencia mínimo para poder traducir
los principios de conservación que rigen los fenómenos de transporte en modelos matemáticos
que nos ayudan a comprender el comportamiento de los procesos naturales y de ingeniería.
En el caso de los procesos naturales, estos modelos nos ayudarán a predecirlos y a tomar las
decisiones que correspondan de acuerdo con tales predicciones. En el caso de los procesos de
ingeniería, estos modelos nos sirven para la creación, control y/o perfeccionamiento de éstos.
Dentro del marco de referencia mínimo están, por un lado, las definiciones conceptuales de los
principios de conservación, sistema y campo. Por otro lado, haremos una presentación de la
traducción de los principios de conservación en términos matemáticos en aras de discutir y
enfatizar el significado y uso de las piezas matemáticas que componen estos modelos. No
seremos ambiciosos. Nos limitaremos a discutir solamente la matemática que necesitamos sin
ahondar, por ahora, en su relación de los procesos físicos que describen. Poco a poco, a medida
que avanzan las clases, iremos conociendo con mayor intimidad esta relación.

3
 Clase I: Principios, sistemas, volúmenes de
control y campos
Lectura inicial
“El aspirante a sacerdote o Jeque desde temprana edad se le aislaba de sus padres y se le
internaba en una “cuca”, lugar donde recibía la formación sacerdotal; se le sometía a ayunos y
duras disciplinas; abstinencia de carne y de todo condimento. Doce años permanecía en el
internado, en perfecta castidad, aprendiendo los secretos curativos de las plantas; se adiestraba
en los conjuros; en el arte de hablar con los demonios y espíritus; se le enseñaba el curso de
los astros; de tal manera que terminaba siendo un erudito en astronomía y en la doctrina y
preceptos de Bachué, en las tradiciones y las glorias pasadas de su pueblo y brindándosele
capacitación en el idioma que tenía que hablarle a los creyentes”.
Eutimio Reyes Manosalva, Leyendas y Mitos Indígenas Cundiboyacenses.
Los cursos de Fenómenos de Transporte tienen muy mala fama tanto entre estudiantes como
entre profesores. En los primeros se genera un temor (apenas entendible) por la naturaleza
abstracta de las matemáticas que se manejan; tensores, operadores diferenciales, ecuaciones
diferenciales en derivadas parciales, soluciones con series y geometría no plana sazonan los
modelos matemáticos que se discuten durante el semestre. En cuanto a los segundos, algunos
colegas observan un respeto distante y a veces un resquemor hacia el curso; es quizá un
sentimiento parecido a no tener ganas de abrazar un erizo. Por experiencia propia, creo que
ambas actitudes pueden en parte entenderse por la frialdad y poco entusiasmo que
tradicionalmente se transmiten a la hora de discutir los fenómenos de transporte y sus modelos
matemáticos. Yendo quizá un poco lejos, me atrevo a decir que, de hecho, una de las causas
de lo anterior son la contemplación, veneración y seguimiento casi ciego del hermoso clásico
de la literatura de ingeniería química llamado “Fenómenos de Transporte” de Robert Byron
Bird, Warren Earl Stewart y Edwin Niblock Lightfoot, Jr; llamado en el ámbito académico
internacional simplemente BSL. Pero, la belleza es un concepto muy subjetivo y transitorio para
los seres humanos. Aunque BSL es sin duda un pilar de la ingeniería química moderna, lo cierto
es que, como muchos clásicos, la mayoría de nosotros no estamos preparados para él cuando
lo conocemos por primera, segunda y hasta quién sabe qué vez. Esto se da en gran medida
por la evolución natural de los modelos educativos. Para la época de la publicación de BSL, la
denominada educación superior seguía un modelo de aprendizaje que podemos llamar de
fuerza bruta amén del hecho que la información no tenía el enorme acceso del que gozamos
hoy. Así, los estudiantes debían romperse el coco durante horas y horas de estudio usando
herramientas que ahora mismo, atrevidamente, podrían ser llamadas precarias. Así, libros como
BSL tratan de manera tosca al estudiante siendo ello notable, por ejemplo, en aquellos pasajes
en que se postulan soluciones infamemente obvias para los escritores del libro. Mas, sin
embargo, todos los esfuerzos intelectuales que desde su publicación han acarreado para
estudiantes y profesores de tratar de seguirle el paso a BSL no han sido, ni son, ni serán en
vano. “La fama cuesta”, dice mi mamá y así como lo expresa la lectura con que abrimos esta
clase. En ella, se hacen patentes las fortísimas demandas intelectuales y físicas que acarreaba
el poder llegar a ocupar un puesto privilegiado en la sociedad indígena cundiboyacense.
Aunque afortunadamente los métodos descritos en la lectura han pasado a la historia, la
esencia del porqué éstos se aplicaban se mantiene. En este sentido, no es mentira que

4
Fenómenos de Transporte es uno de los cursos más difíciles de Ingeniería Química y que como
tal su estudio exige muchísimo tiempo y disciplina. Pero, tal como a menudo sucede en las
románticas historias de búsquedas de tesoros, al final, el héroe (en este caso usted como
estudiante) se verá recompensado por una riqueza intelectual y modo de pensamiento que le
permitirá tanto hacer un ejercicio ingenieril valioso como una serie de reflexiones sobre la
naturaleza misma que antes quizá no habría imaginado.

Objetivo: Presentar la definición de los principios de conservación de la masa, momento y

energía y las definiciones de sistema, volumen de control y campo en el marco de los
fenómenos de transporte.

Temas por tratar

1. Definición de los principios de conservación de la masa, cantidad de movimiento y
energía.
2. Presentación de los principios de conservación en términos de ecuaciones de cambio.
3. Presentación de las definiciones de sistema, volumen de control y campo

Apuntes de clase
Resumen general. Los principios de conservación de la masa, cantidad de movimiento y
energía son el pilar del análisis y modelamiento de los fenómenos de transporte. Estos principios
dicen que las cantidades mencionadas; masa, cantidad de movimiento y energía, se conservan;
esto implica que si éstas se cuentan siempre dan el mismo resultado sin importar ni el método
de conteo ni el marco de referencia usado para contarlas. En ingeniería química, todos estos
principios se expresan como: lo que entra menos lo que sale es igual a lo que se acumula en
un volumen de control dado. El principio de conservación de la masa estipula que no es posible
generar masa de la nada. El de cantidad de movimiento se refiere al hecho que cuando una
cierta masa colisiona con otra la cantidad de movimiento; es decir, el producto de la masa por
la velocidad se conserva; y, el de energía postula que la energía no se crea ni se destruye solo
se transforma. Uno de los objetivos de este curso es traducir estos principios en expresiones
matemáticas que permitan modelar los procesos de transporte tanto en sistemas naturales como
en procesos de ingeniería. Este modelamiento tiene como fin último hacer predicciones
concretas y diseñar estrategias de control y optimización de procesos.

¿Qué son los fenómenos de transporte?

Vamos a definir qué son los fenómenos de transporte a partir de otra rama de la ciencia que
los estudiantes de este curso ya conocen: la termodinámica clásica:
“La termodinámica clásica es una rama de los fenómenos de transporte que estudia las
relaciones entre las diferentes formas de energía cuando éstas han alcanzado un estado de
equilibrio.”
La definición anterior puede sonrojar a más de un colega, pero, una reflexión desligada de
apegos muestra que: (i) los fenómenos de transporte estudian las relaciones entre las diferentes
formas de energía teniendo en cuenta tanto sus mecanismos de transferencia como las
propiedades ligadas a éstas en un estado de equilibrio; (ii) los fenómenos de transporte
estudian la interrelación entre las transferencias de energía, masa y cantidad de movimiento.
Para un ingeniero químico, esto puede resumirse con la siguiente ilustración:

5
Figura I. Representación de las principales ramas que componen los Fenómenos de Transporte.
Los Fenómenos de Transporte lidian esencialmente con la aplicación de los principios de
conservación postulados en física. Particularmente, los tres principios de conservación que se
discuten en este curso de fenómenos de transporte son:
1. La masa no se crea ni se destruye: principio de conservación de la masa.
2. Segunda ley del movimiento de Newton: principio de conservación de la cantidad de
movimiento (momento) lineal.
3. La energía no se crea ni se destruye sólo se transforma, primera ley de la termodinámica:
principio de conservación de la energía.
Estos principios de conservación los traduciremos en términos matemáticos en términos de
“ecuaciones de cambio”. Las ecuaciones de cambio pueden escribirse de manera diferencial o
de manera integral. En una clase futura veremos que las dos maneras de escribir las ecuaciones
de cambio son intercambiables entre sí y que además nos ayudan a comprender y a describir
los fenómenos de transporte con dos filosofías (Euleriana y Lagrangiana) y escalas de análisis
(microscópica y macroscópica) diferentes. Con excepción de la ecuación de cambio que
corresponde al principio de conservación de la energía, las ecuaciones de cambio que
representan los principios de conservación de la masa y la cantidad de movimiento tienen un
nombre particular, a saber: ecuación de continuidad y ecuación de Navier-Cauchy,
respectivamente. La tabla siguiente resume ejemplos de la traducción de los principios de
conservación en ecuaciones de cambio del tipo diferencial.

6
Tabla 1. Ejemplos de la traducción matemática de los principios de conservación en ecuaciones
de cambio del tipo diferencial
Principio de
Ecuación de cambio Expresión matemática diferencial*
conservación

𝜕𝜌𝑖
Masa** Continuidad ⃗ ∙ 𝑛⃗𝑖 ) = 𝑟𝑖
+ (∇
𝜕𝑡

Navier-Cauchy o
Movimiento 𝜕[𝜌𝑣 ]
cantidad de ⃗ ∙ {𝜌𝑣 𝑣 }] = −[∇
+ [∇ ⃗ ∙ 𝜋̿] + [𝜌𝑔]
lineal 𝜕𝑡
movimiento
𝜕 1 1
̌ + ‖𝑣 ‖2 )) + (∇
(𝜌 (𝑈 ⃗ ∙ [𝜌 (𝑈 ̌ + ‖𝑣 ‖2 ) 𝑣 ])
𝜕𝑡 2 2
Energía Energía
= −(∇⃗ ∙ 𝑞 ) − (∇
⃗ ∙ [𝜋̿ ∙ 𝑣 ]) + 𝜌(𝑣 ∙ 𝑔)

*Las expresiones matemáticas mostradas varían de acuerdo con las consideraciones físicas y el método
matemático que se use para escribirlas. Para una revisión de las múltiples traducciones matemáticas de
los principios de conservación, consultaremos (a su debido tiempo) el libro de Fenómenos de Transporte
de Bird, Stewart & Lightfoot y, en particular, su apéndice B.
**El principio de conservación de la masa no se cumple bajo un marco de referencia relativista en el que
es reemplazado por el principio de interconversión de la masa con la energía.

Los principios de conservación se aplican para analizar cómo evolucionan, transportan y

transforman la masa, la cantidad de movimiento y energía en un sistema o en un volumen de
control.

¿Qué es un sistema?
En Fenómenos de Transporte un sistema es una cantidad de materia fija que se desea estudiar.
Todo lo que no está dentro del sistema son los alrededores. Para poder analizar el
comportamiento de un sistema es importante saber cuáles son sus condiciones iniciales, qué
propiedades fisicoquímicas de estado, termodinámicas, y de transporte, dinámicas; i.e.,
constitutivas, tiene y qué pasa en sus fronteras. Más adelante, cuando empecemos a discutir el
significado físico de las ecuaciones de cambio, veremos que en la jerga de Fenómenos de
Transporte que usan los que estudian mecánica de fluidos, la palabra sistema es equivalente a
la definición de “sistema cerrado” que se usa en termodinámica clásica. Es importante hacer
ejemplos para memorizar la definición de sistema puesto que es sobre ellos que postulamos
los principios de conservación y dependiendo de cómo definamos la región del espacio-tiempo
que los constituye plantearemos las ecuaciones de cambio correspondientes.

¿Qué es un volumen de control?

Un volumen de control es una región fija del espacio-tiempo a través de la cual se puede
transportar masa, cantidad de movimiento y energía. Para la termodinámica clásica, un volumen
de control es análogo a un sistema abierto.

¿Qué es un campo?
Desde el punto de vista de la física, un campo es la distribución espacio temporal de una
magnitud física. En otras palabras, un campo es una representación de una propiedad que

7
puede medirse en el entorno de cada punto de una región del espacio para cada instante de
tiempo. Los campos tienen atributos matemáticos de acuerdo con las distribuciones que
representan. Para los alcances de este curso, necesitamos manejar ciertas propiedades de tres
tipos de campos: escalares, vectoriales y tensoriales de orden 2. El siguiente esquema los
describe:

De este esquema, los siguientes mensajes son clave:

(1) Los campos tensoriales engloban todos los campos posibles. Así, un campo escalar es, de
hecho, un tensor de orden 0 y un campo vectorial es un tensor de orden 1. Para los alcances
del curso, nos basta con utilizar tensores de orden 2 cuya representación matemática puede
aproximarse mediante una matriz de tres filas y tres columnas.

(2) Los campos escalares pueden ser números o funciones con dimensiones (compuestos) o sin
dimensionales (adimensionales, simples). Las dimensiones básicas empleadas en este curso
son longitud (L), tiempo (T), temperatura (θ) y cantidad de materia (M). Una dimensión no es
una unidad. Una unidad es una manera de expresar una dimensión. Así, el metro, m, el pie, ft,
y el jeme1 y son maneras de expresar la longitud. Por supuesto, las dimensiones básicas se
combinan para expresar las cantidades físicas que uno requiera cuando estudia un sistema. De
este modo, el producto L×L sirve para expresar un área; que es una longitud en dos
dimensiones. También, el cociente L/T sirve para expresar la rapidez de un movimiento. En
general, cuando se indiquen las dimensiones de una cantidad usaremos el símbolo ⟦=⟧.

Es de cierto modo asombroso el hecho que el estudio de un sistema y sus propiedades puede
abordarse de manera muy sencilla solamente discutiendo sus dimensiones básicas; a esto se
le llama análisis dimensional. Durante el curso, haremos uso del análisis dimensional para
aprender tanto sobre las propiedades de los sistemas como para deducir ciertos aspectos de
su comportamiento. En un sentido amplio, conocer las dimensiones de una propiedad
fisicoquímica -de estado o de transporte- es definirla.

1
Un jeme es la distancia que hay entre los dedos pulgar e índice de una persona cuando pone su mano
de modo que los dos dedos se extienden sobre una superficie tal como una mesa o una pared.

8
Por otro lado, en el contexto de este curso, un campo escalar adimensional es un número Real.
Para los Fenómenos de Transporte, estos números adimensionales sirven para describir y
comparar los tipos de comportamiento que puede presentar el sistema. Por ejemplo, el número
de Reynolds, Re, sirve para comparar las fuerzas de inercia con las fuerzas viscosas cuando un
fluido está en movimiento.
(3) Un campo vectorial es un tensor de orden uno y tiene magnitud, dimensión y dirección.
(4) Un campo tensorial de segundo orden tiene magnitud, dimensión, dirección y planos de
referencia. Estos últimos se construyen sobre la base de vectores unitarios. Existen varias
maneras de representar un tensor de orden dos. Aquí, usaremos o una representación tipo
matricial 3×3 o la notación de Einstein para representarlos, según nos convenga. En todo caso,
los tensores que usaremos tienen nueve componentes escalares. Por razones prácticas, en
adelante diremos los términos tensor y tensor de orden dos se usarán de manera
intercambiable.

Con base en lo anterior, definimos la siguiente notación para representar la naturaleza del
campo que está presente en una ecuación dada y que resulta al realizar operaciones entre
escalares, vectores y tensores.
• Usaremos ( ) si el resultado de la operación es un escalar.
• Usaremos [ ] si el resultado de la operación es un vector.
• Y, usaremos { } si el resultado de la operación es un tensor.

Adoptamos aquí la notación propuesta por Bird, Stewart & Lightfoot [1] en el Apéndice A de
su libro de Fenómenos de Transporte. De hecho, el estudiante debe imprimir este Apéndice y
siempre tenerlo a la mano para consulta durante el curso.

Por otro lado, definiremos que las cantidades vectoriales, vectores, los denotaremos en
adelante usando una flecha encima de ellos; e.g., 𝑣 es el vector velocidad, y que los tensores
los denotaremos con una barra doble encima; e.g., 𝜋̿ es el tensor densidad de flujo (“flux”)
combinado de cantidad de movimiento.

Ejemplo I.1. Características matemáticas de las ecuaciones de cambio.

Diga algunas propiedades matemáticas de las ecuaciones de cambio mostradas en la Tabla 1.
Respuesta. Con base en lo discutido, podemos identificar algunas de las características de las
ecuaciones de cambio que se presentan en la Tabla 1.

(1) Las cantidades: (∇ ̌ + 1 ‖𝑣 ‖2 )), (∇

⃗ ∙ 𝑛⃗𝑖 ) -ecuación de continuidad-, (𝜌 (𝑈 ⃗ ∙ [𝜌 (𝑈
̌+
2
1
‖𝑣 ‖2 ) 𝑣 ]), ⃗ ∙ 𝑞 ), (∇
(∇ ⃗ ∙ [𝜋̿ ∙ 𝑣 ]) y (𝑣 ∙ 𝑔) -ecuación de energía- son escalares.
2

⃗ ∙ {𝜌𝑣 𝑣 }], [∇
(2) Las cantidades [𝜌𝑣 ], [∇ ⃗ ∙ 𝜋̿] y [𝜌𝑔] -ecuación de Navier-Cauchy- son vectores.

(3) La cantidad {𝜌𝑣 𝑣 } es un tensor de segundo orden.

Ejemplo I.2. Consistencia algebraica.

Diga qué tipos de campos describen las ecuaciones de cambio.

9
Respuesta. Uno puede identificar el tipo de campo que una ecuación de cambio describe
después de que haya identificado cualquier término de los que aparecen en ella ya que para
que la ecuación sea consistente todos sus términos deben describir el mismo tipo de campo.
Usando los resultados del ejemplo 1.1 podemos deducir lo siguiente:

(1) Las ecuaciones de continuidad y energía describen campos escalares porque tienen
términos escalares.
(2) La ecuación de Navier-Cauchy (o de cantidad de movimiento) describe un campo vectorial.

Teniendo en cuenta que los términos de las ecuaciones de cambio se suman, la consistencia
aquí descrita se puede comparar con la de cualquier operación algebraica, de manera figurada.
Así, si uno tiene un saco de ñames, papas y yucas solamente puede sumar ñames con ñames,
papas con papas y yucas con yucas. Así mismo, un escalar solamente se suma con otro escalar,
un vector con otro vector y un tensor con otro tensor.

Ejemplo I.3. Consistencia dimensional.

En consonancia con el principio de consistencia algebraica, las ecuaciones de cambio también

son consistentes dimensionalmente. Así, cada término en una de ellas tiene (debe tener) las
mismas dimensiones. Esto nos dice ya bastante sobre el carácter físico de las ecuaciones de
cambio, veamos:

(1) En la ecuación de continuidad, sabiendo que 𝜌𝑖 es la densidad másica -masa por unidad
𝑀
de volumen- de una sustancia i; dimensiones: 𝜌𝑖 ⟦=⟧ 3 , y que 𝑡 ⟦=⟧ 𝑇 es el tiempo,
𝐿
podemos decir que las dimensiones de cada término de la ecuación de continuidad son:
𝑀
𝜕𝜌𝑖 3 𝑀
⟦=⟧ 𝐿 ⟦=⟧ 3
𝜕𝑡 𝑇 𝐿 ×𝑇

⃗ ∙ 𝑛⃗𝑖 ) y 𝑟𝑖 también tienen las dimensiones 3𝑀 .

Por tanto, los términos: (∇ 𝐿 ×𝑇

(2) En la ecuación de cantidad de movimiento, el término [𝜌𝑔] es el producto de la densidad

𝐿
por el vector aceleración por gravedad, 𝑔 ⟦=⟧ 2, por tanto:𝑇
𝑀 𝐿 𝑀
𝜌𝑔 ⟦=⟧ 3
× 2 ⟦=⟧ 2
𝐿 𝑇 𝐿 × 𝑇2
⃗]
𝜕[𝜌𝑣 𝑀
Por tanto, ⃗ ∙ {𝜌𝑣 𝑣 }] ⟦=⟧ [∇
⟦=⟧ [∇ ⃗ ∙ 𝜋̿] ⟦=⟧ .
𝜕𝑡 𝐿2 ×𝑇 2

(3) Con la información que tenemos hasta ahora, podemos también saber las dimensiones de
los términos de la ecuación de energía analizando las dimensiones del término 𝜌(𝑣 ∙ 𝑔):
𝑀 𝐿 𝐿 𝑀
𝜌(𝑣 ∙ 𝑔) ⟦=⟧ × × ⟦=⟧
𝐿3 𝑇 𝑇 2 𝐿 × 𝑇3
Por tanto,

𝜕 1 1 𝑀
̌ + ‖𝑣‖2 )) ⟦=⟧ (∇
(𝜌 (𝑈 ⃗ ∙ [𝜌 (𝑈
̌ + ‖𝑣 ‖2 ) 𝑣 ]) ⟦=⟧ (∇
⃗ ∙ 𝑞 ) ⟦=⟧ (∇
⃗ ∙ [𝜋̿ ∙ 𝑣]) ⟦=⟧
𝜕𝑡 2 2 𝐿 × 𝑇3

10
Ejercicios
I.1. Averigüe a qué tipo de cantidades físicas corresponden las dimensiones de cada ecuación
de cambio. Consejo: comience usando Wikipedia en inglés.

11
Notas de estudio
La Tabla 2 presenta detalles adicionales sobre la traducción de los principios de conservación
en ecuaciones de cambio.
Tabla 2. Descripción detallada de la traducción de los principios de conservación en ecuaciones
de cambio.
Principio de Nombre
Expresión matemática Consideraciones
conservación matemático
Ecuación para un
sistema con reacción
química con múltiples
𝜕𝜌𝑖
⃗ ∙ 𝑛⃗𝑖 ) + 𝑟𝑖
= −(∇ Ecuación de especies químicas.
De la masa 𝜕𝑡 continuidad Su objetivo primario es
i = 1, 2, 3…
determinar el perfil de
concentración de cada
especie en el sistema.
Ecuación para un
sistema con múltiples
especies químicas y
donde solamente se
considera la gravedad
como fuerza externa. La
𝑛
𝜕[𝜌𝑣 ] ecuación es equivalente
⃗ ∙ {𝜌𝑣 𝑣 + 𝜋̿} ] + ∑(𝜌𝑖 𝑔𝑖 ) Ecuación de
De la cantidad de = −[∇ a la segunda ley del
𝜕𝑡 cantidad de
movimiento 𝑖=1 movimiento de Newton:
movimiento 𝑛

∑ 𝐹𝑖 = 𝑚𝑎
𝑖=1
Su objetivo principal es
determinar el perfil de
velocidad de cada
especie en el sistema.
1
̌ + 𝑣2)
𝜕𝜌 (𝑈
2 ⃗ Corresponde a una
= − (∇ expresión de la primera
𝜕𝑡
1 ley de la
∙ (𝜌 (𝑈̌ + 𝑣 2 ) 𝑣) + 𝑞
2 termodinámica:
Ecuación de
De la energía + [𝜋̿ ∙ 𝑣 ])
̌ = 𝛿𝑄 + 𝛿𝑊
𝑑𝑈
energía
𝑛
Su objetivo principal es
determinar el perfil de
+ ∑(𝑛⃗𝑖 ∙ 𝑔𝑖 ) temperatura en el
𝑖=1 sistema.

Definiciones básicas
Sistema: es una región del espacio – tiempo que contiene una cantidad definida de materia.
Esta definición es más fundamental que aquella usada en la termodinámica clásica puesto que
se refiere a las variables fundamentales de la física: posición y tiempo. En este respecto, el
estudiante puede recordar que la termodinámica clásica define los sistemas en términos de
propiedades de estado; específicamente, presión, temperatura, cantidad de sustancia y
volumen. Además, el análisis termodinámico clásico está fundamentado en el concepto de

12
equilibrio para el cual la variable tiempo no es relevante. Así las cosas, uno puede concluir que
las leyes de la termodinámica clásica son un caso particular de las ecuaciones de cambio.
Además, en el sentido estricto, en termodinámica clásica, los principios de conservación se
cumplen a cabalidad en los denominados sistemas aislados; es decir, aquellos no hacen ningún
tipo de intercambio con los alrededores. La Figura 2 muestra un sistema arbitrario, sus
fronteras y los alrededores. Como se puede notar, la forma de un sistema no influye en su
definición.

Figura 2. Representación arbitraria de un sistema, su frontera y los alrededores.

Tres elementos de análisis que significan empezar con pie derecho el análisis del transporte de
una cantidad son: (i) saber definir cuál es el sistema y cuál la frontera. Aunque suene un tanto
inaudito, este primer elemento es en el que más erran los estudiantes de pregrado y posgrado
e incluso los profesionales y hasta profesores de ingeniería (¡!). (ii) Conocer qué ocurre en la
frontera del sistema y cuál es la historia del sistema; condición inicial. Y, (iii) Caracterizar las
propiedades del sistema. En particular, para analizar los procesos de transporte se deben
conocer tanto las propiedades de estado del sistema; e.g., entalpías, entropías, potencial
químico, fugacidad, étc, como las propiedades de transporte del mismo; e.g., difusividad,
viscosidad, conductividad térmica, etc. Estas últimas son propiedades que no se tratan en los
cursos de termodinámica básica y, por tanto, serán objeto de estudio en este curso.
Campo. El otro concepto clave en la comprensión de las ecuaciones de cambio es el concepto
de campo. Este concepto es clave tanto para entender el carácter físico de las ecuaciones como
para su manejo matemático. Una consulta de Wikipedia (ver aquí) muestra la siguiente
definición de campo en el marco de la física: “…es una cantidad física representada por un
escalar, vector o tensor que tiene un valor para cada punto en el espacio-tiempo”. Con base
en esta definición, es claro que las ecuaciones de cambio, Tablas 1 y 2, están expresadas en
términos de campos escalares, vectoriales y tensoriales. La Tabla 3 muestra las propiedades
de cada uno de estos campos y clasifica los términos mostrados en las ecuaciones de cambio
de acuerdo con el tipo de campo representado.

13
Tabla 3. Descripción de las ecuaciones de cambio en términos de campos.
Términos correspondientes
Tipo de campo Propiedades en las ecuaciones de
cambio
Los campos escalares poseen magnitud
y dimensión. Desde el punto de vista
𝜕𝜌𝑖 𝜕∎
Escalar aritmético son números compuestos 𝜌𝑖 , , 𝑟𝑖 , ̌, 𝑣 2 , 𝑡,
𝑈
𝜕𝑡 𝜕𝑡
* e.g., 2 litros (1 L), cinco kilogramos (5
Kg)
Los campos vectoriales poseen
magnitud, dimensión y dirección.

Vectorial 𝜕𝜌𝑣
𝑛
𝑣 , ⃗∇, 𝑔𝑖 , 𝑛⃗𝑖 , 𝑞
⃗ ∙ {𝜌𝑣 𝑣 + 𝜋̿} ] + ∑(𝜌𝑖 𝑔𝑖 )
= −[∇
𝜕𝑡
𝑖=1

Los campos tensoriales poseen

Tensorial magnitud, dimensión, dirección y plano 𝜋̿, 𝑣 𝑣.
de referencia.

*Ver: Aritmética de Baldor.

14
 Clase II: Álgebra y operadores
Lema: “Mathematics is not just a language; mathematics is a language plus reasoning; is like a
language plus logic.”
R.P. Feynman, Feynman’s messenger lectures on The Character of physical law at Cornell
University (1964).
Ver: https://www.youtube.com/watch?v=kEx-gRfuhhk&t=5668.

Objetivo: Ilustrar las operaciones algebraicas entre escalares, vectores y tensores que se usan
en las ecuaciones de cambio, así como presentar los operadores implicados en las mismas;
derivada parcial, nabla, gradiente, divergencia, rotacional, laplaciano y derivada sustancial.

Apuntes de clase
Resumen. Las expresiones matemáticas mostradas en las Tablas 1 y 2; llamadas en general
ecuaciones de cambio, comprenden una serie de términos que representan cantidades de
diverso orden: escalares, vectores y tensores. Estas cantidades se combinan en productos de
diferente índole; escalar, vectorial y diádico, y sobre ellas operan diversos operadores
matemáticos; derivadas parciales, divergencias y gradientes. Por el momento, no ahondaremos
mucho en el carácter físico de las ecuaciones de cambio, sino que estableceremos los elementos
matemáticos necesarios para entenderlas matemáticamente como para aprender a aplicarlas
de manera exitosa en los fenómenos de transporte. Por simplicidad, las operaciones que se
ilustrarán en esta sección están todas referidas a un sistema de coordenadas cartesianas.
Campos escalares, vectoriales y tensoriales
Una nueva mirada a las ecuaciones de cambio revela que cada término de ellas comprende una
multiplicación entre cantidades. Por eso, es bueno ver cuáles son estas operaciones y en qué
resultan.

̌ + 1 ‖𝑣 ‖2 ) -ecuación de energía- es una multiplicación de escalares.

Primero, el término 𝜌 (𝑈
2
Por tanto, resulta en otro escalar.

⃗ ∙ 𝑛⃗𝑖 ) -ecuación de continuidad-, (𝜌 (𝑈

Segundo, como dijimos antes, los términos (∇ ̌+
1
‖𝑣 ‖2 )), (∇ ̌ + 1 ‖𝑣 ‖2 ) 𝑣 ]), (∇
⃗ ∙ [𝜌 (𝑈 ⃗ ∙ 𝑞 ), (∇
⃗ ∙ [𝜋̿ ∙ 𝑣 ]) y (𝑣 ∙ 𝑔) -ecuación de energía- son
2 2
escalares y, por tanto, los productos que están descritos en ellos resultan en escalares. De esto
podemos recordar un poco de álgebra a la vez que deducimos que: (i) la multiplicación de un
escalar y un vector da un nuevo vector, (ii) el producto punto entre dos vectores resulta en un
escalar y por eso también es llamado producto escalar, en álgebra lineal, y (iii) el producto
punto entre un tensor de orden dos y un vector resulta en un vector.

Tercero, el término {𝜌𝑣 𝑣 } es un tensor de orden dos (¡!). Por tanto, la multiplicación (directa)
de dos vectores da un tensor de orden dos, ya que, así como la multiplicación de un escalar
con un vector da otro vector, la multiplicación de un escalar por un tensor da otro tensor del
mismo orden. En general, a los productos directos del tipo 𝑣 𝑣 se les llama productos diádicos.

15
Las observaciones de arriba se pueden sistematizar en cuanto a cuál es el orden de las
cantidades que resultan al hacer diferentes tipos de multiplicaciones entre escalares, vectores
y tensores. Esto está resumido en la siguiente tabla que se adaptó del Apéndice A del libro de
Bird et al.[1]:
Tabla 4. Reglas rápidas para determinar el orden de la cantidad resultante al multiplicar
escalares, vectores y tensores.
Orden de la cantidad del
Nombre del tipo de
Signo de multiplicación resultado de la
multiplicación
multiplicación
El orden es igual a la suma
Ninguno Directa o diádica de los órdenes de cada
término que se multiplica
El orden es igual a la suma
de los órdenes de cada
término que se multiplica
∙ Producto punto
menos dos.
Nota: el producto punto no
está definido para escalares.
El orden es igual a la suma
de los órdenes de cada
término que se multiplica
× Producto cruz* menos uno.
Nota: el producto cruz
solamente lo realizaremos
entre vectores.
El orden es igual a la suma
de los órdenes de cada
término que se multiplica
Producto doble punto o
: menos cuatro.
escalar*
Nota: el producto escalar
solamente lo realizaremos
entre tensores de orden dos.
*¡Calma! Estas operaciones aparecerán en clases futuras.

Tarea II.1. Verifique que la notación adoptada para cada término de las ecuaciones de cambio
( ), [ ] y { } es congruente con el orden de las multiplicaciones que se hacen.

Comenzamos repasando los tipos de cantidades que usaremos durante el curso. Como dijimos
antes, en el curso de fenómenos de transporte se utilizan diferentes cantidades físicas que
sirven para describir tres tipos de campos físicos, a saber:
Campo escalar

Está constituido por números o funciones Reales (o Complejas, en el sentido más amplio de la
mecánica cuántica) que tienen magnitud y pueden tener dimensión, pero no dirección. Ejemplos
de campos escalares son:

16
(i) La función 𝑓(𝑥) = 𝑥 2 definida en todo el dominio de los números Reales. En este caso,
cuando 𝑥 = 1, 𝑓(𝑥) = 1 que es un escalar simple o adimensional. Este campo puede
representarse gráficamente.
(ii) La ecuación de continuidad definida para un sistema con una única sustancia que no
reacciona:
𝜕𝜌
⃗ ∙ 𝜌𝑣 ) = 0
+ (∇
𝜕𝑡

La ecuación describe un campo escalar para el cambio de la densidad, 𝜌, de la sustancia

investigada en función del tiempo y a lo largo de la trayectoria que la sustancia sigue. A
diferencia del ejemplo (i), el campo escalar descrito para la densidad tiene dimensiones que se
pueden determinar a partir del principio de consistencia dimensional tal como lo hicimos en la
clase anterior, así:

𝜕𝜌 𝑀
⟦=⟧ 3
𝜕𝑡 𝐿 ×𝑇
𝜕𝜌
A partir de la expresión anterior, se aprende que el término 𝜕𝑡 sirve para conocer cómo
evoluciona la densidad de la sustancia con el tiempo.
Dimensionalmente, el segundo término de la ecuación de continuidad tiene las mismas
dimensiones:
𝑀
⃗ ∙ 𝜌𝑣 ) ⟦=⟧
(∇ 3
𝐿 ×𝑇

Pero, si hurgamos las dimensiones de cada componente del término podemos aprender otras
cosas sobre su naturaleza física. Así, sabiendo que 𝑣 es el vector velocidad (campo de
⃗?
velocidad) para la sustancia de densidad 𝜌, ¿cuáles son las dimensiones de ∇
𝑀 𝑀 𝐿
⃗ ∙ 𝜌𝑣 ) ⟦=⟧
(∇ ⃗ ⟧×
⟦=⟧ ⟦∇ ×
𝐿3 ×𝑇 𝐿3 𝑇
Simplificando,
1
⃗ ⟦=⟧
∇
𝐿
⃗ es un vector que opera sobre 𝜌𝑣 en términos de la direccionalidad de este último.
Por tanto, ∇
Campo vectorial. A diferencia de un escalar, los vectores representan cantidades en el espacio
que tienen magnitud y dirección y, para este curso, deben tener dimensiones. Un vector se
representa matemáticamente en términos de su magnitud (que es un escalar) y su dirección en
función de vectores unitarios referidos a un determinado sistema de coordenadas. Para este
curso de fenómenos de transporte, los vectores se refieren a coordenadas cartesianas (planas),
cilíndricas y esféricas.

Como ejemplo, para coordenadas cartesianas, el vector velocidad se puede representar de las
siguientes formas equivalentes:

Tabla 5. Representación del vector velocidad en el marco de diferentes campos de estudio.

Notación Vectores unitarios Aplicación Referencias

17
 JW Jewett &
Serway
𝑣 = 𝑣𝑥 𝑖̂ + 𝑣𝑦 𝑗̂+𝑣𝑧 𝑘̂ 𝑖̂, 𝑗̂, 𝑘̂ Mecánica, cinemática Physics for
engineers &
Scientists
Mecánica y dinámica de V.L. Streeter
𝑣 = 𝑢𝑖̂ + 𝑣𝑗̂ + 𝑤𝑘̂ 𝑖̂, 𝑗̂, 𝑘̂ fluidos Fluid dynamics
Bird, Stewart,
Física teórica – Fenómenos Lightfoot
𝑣 = 𝑣𝑥 𝛿𝑥 + 𝑣𝑦 𝛿𝑦 +𝑣𝑧 𝛿𝑧 𝛿𝑥 , 𝛿𝑦 , 𝛿𝑧 de transporte Transport
Phenomena

Durante el curso, vamos a adoptar la representación de vectores unitarios como 𝛿𝑖 en

congruencia con el libro de Fenómenos de Transporte de Bird, Stewart & Lightfoot (BSL) [1].

Campo tensorial. Desde el contexto del curso de Fenómenos de Transporte, un tensor de

segundo orden se refiere a una cantidad que posee magnitud, dirección, dimensión y está
sujeta a un plano de referencia. Los tensores de orden 2 son objetos algebraicos que pueden
representarse como arreglos tipo matricial 3×3. Por tanto, los tensores que usaremos tendrán
9 componentes escalares que estarán referidos a tres vectores unitarios que constituyen la
base vectorial del tensor.

Ejemplo II-1: Considere el envío de un paquete de correspondencia. Para ello, en el sistema de

nomenclatura de Colombia se deben indicar: calle o carrera, número y barrio de la ciudad donde
está el lugar a donde se remite el paquete. En este ejemplo, los dos primeros datos pueden
asociarse con una cantidad vectorial mientras que el barrio (que es una colección de casas que
abarca calles y carreras) puede interpretarse como un tensor. Suponiendo que una persona
vive en la carrera 12 #25-55. Al graficar en un plano cartesiano esa dirección se puede obtener
la siguiente ilustración:

Figura 3. Representación de cantidades escalares, vectoriales y tensorial a través de la

ubicación de una dirección de residencia.

18
En la figura 3, el 55 representa una cantidad escalar: distancia desde la esquina hasta la casa,
mientras que la combinación del número con la carrera representa un vector. Luego, la cuadra
donde se encuentra la casa; que está delimitada por las carreras y calles que van de esquina a
esquina representan un componente del tensor: Barrio.
Actividad. Con base en la información de arriba, determine las componentes del vector posición
de la dirección suministrada. Si cada cuadra es el componente del tensor Barrio, ¿cuántos hay
en la ilustración dada? Asigne una nomenclatura para estos componentes.

Complementando lo dicho en la Clase I, para este curso, denotaremos los tensores con letras
̿ representa el tensor de
griegas usando negrillas y con dos barrar superiores. Por ejemplo: 𝝓
esfuerzos combinados en el marco del modelamiento de los procesos de transporte de cantidad
de movimiento. Dos de los tensores más usados en todo el curso son el tensor de díadas
unitarias, ̿
𝜹, que se representa como:

𝛿11 𝛿12 𝛿13

(𝛿21 𝛿22 𝛿23 )
𝛿31 𝛿32 𝛿33

Y el tensor de esfuerzos cortantes o viscosos, 𝝉̿, que se representa como:

𝜏11 𝜏12 𝜏13
𝜏
( 21 𝜏22 𝜏23 )
𝜏31 𝜏32 𝜏33

Supongamos que se quiere ilustrar geométricamente el tensor de esfuerzos cortantes (𝝉̿) para
un cubo en un sistema de coordenadas cartesianas. La figura 4 ilustra una manera de ubicar
cada uno de los componentes de este tensor.

19
Figura 4. Representación geométrica del tensor de esfuerzos cortantes (𝝉̿) en un sistema de
coordenadas cartesiano.

Orden de las operaciones multiplicativas que están en las ecuaciones de cambio

Aparte de las operaciones de suma y resta que el estudiante sabe realizar entre escalares,
vectores y matrices (que para los alcances del curso representan los tensores de orden 2), en
este curso son esenciales algunas operaciones de multiplicación entre las cantidades antes
citadas. Antes de ello, es necesario recordar que a cada una de las cantidades que se usarán
en el curso se les ha asignado un orden matemático, así: los escalares tienen orden cero, los
vectores orden uno y los tensores son de orden dos. También, adaptaremos la simbología
planteada por BSL para señalar qué tipo de cantidad produce una determinada multiplicación
cuando estas se indiquen en las ecuaciones que se manejarán. Así, para señalar que una
operación resulta en un escalar, ésta se encerrará entre paréntesis, (), para operaciones que
resulten en un vector, éstas se encerrarán dentro de corchetes, [], y para operaciones que
resulten en un tensor de orden dos, éstas se encerrarán entre llaves, {}.

A continuación, se presentan una serie de ejemplos de estas operaciones.

Ejemplo II.2. Use la Tabla 4 para decir cuál es el orden de los resultados de las siguientes
multiplicaciones y use la nomenclatura establecida para indicarlo.

a) 𝒗 ⃗ . En este caso, se trata de una operación entre vectores sin ningún signo de multiplicación,
⃗ 𝒗
por tanto, el resultado tendrá como orden: 1 + 1 = 2; es decir, se obtiene un tensor de orden
dos. La operación se indica entre llaves: {𝒗 ⃗ 𝒗
⃗ }.

b) 𝒗
⃗ × 𝝉̿. Siendo un producto cruz, el resultado tendrá como orden: 1 + 2 – 1 = 2, es decir,
el resultado es un tensor de orden dos y la operación indica así: {𝒗
⃗ × 𝝉̿}.

c) 𝝉̿ ∶ 𝝈
̿ . Este es un producto doble punto entre dos tensores. Por tanto, el resultado tiene
orden: 2 + 2 – 4 = 0. Es un escalar y esto se señala como: (𝝉̿ ∶ 𝝈
̿ ).

Operaciones típicas entre vectores

Durante el curso, será usual hacer productos punto y cruz entre vectores. A manera de repaso,
definimos un vector velocidad: 𝒗
⃗ = 𝒗𝒙 𝛿𝑥 + 𝒗𝒚 𝛿𝑦 + 𝒗𝒛 𝛿𝑧 y hacemos los productos punto y cruz
de este vector consigo mismo, así:

Producto punto

(𝒗 ⃗ ) = (𝒗𝒙 )(𝒗𝒙 )(𝛿𝑥 ∙ 𝛿𝑥 ) + (𝒗𝒚 )(𝒗𝒚 )(𝛿𝑦 ∙ 𝛿𝑦 ) + (𝒗𝒛 )(𝒗𝒛 )(𝛿𝑧 ∙ 𝛿𝑧 ) = 𝒗𝟐𝒙 (𝟏) + 𝒗𝟐𝒚 (𝟏) + 𝒗𝟐𝒛 (𝟏)
⃗ ∙𝒗
=

= 𝒗𝟐𝒙 + 𝒗𝟐𝒚 + 𝒗𝟐𝒛 = ‖𝑣 ‖2 .

Hacemos las siguientes observaciones:

(1) El producto punto de un vector unitario consigo mismo es igual a uno. Esto es
consecuente con la definición trigonométrica del producto punto, que es:

20
 ⃗ ) = ‖𝑎‖‖𝑏⃗‖𝑐𝑜𝑠𝛼. Donde, ‖𝑎‖ y ‖𝑏⃗‖ son las magnitudes de los vectores y 𝛼 es el
⃗ ∙𝒃
(𝒂
ángulo que forman entre ellos. Para un vector unitario, e.g., 𝛿𝑧 , se tiene que su magnitud es
‖𝛿𝑧 ‖ = 1 y que el ángulo forma consigo mismo es 𝛼 = 0°. Así, 𝑐𝑜𝑠𝛼 = 1. Por tanto,
(𝛿𝑧 ∙ 𝛿𝑧 ) = 1.

(2) El producto punto de un vector consigo mismo es el cuadrado de su magnitud, (𝒗 ⃗)=

⃗ ∙𝒗
‖𝑣 ‖2 .

Producto cruz

Usando el mismo vector velocidad de antes, hagamos el producto cruz de éste consigo
mismo recurriendo a la regla del determinante tal como ya se vio en los cursos de álgebra
lineal o física previos a este de fenómenos de transporte.

𝛿𝑥 𝛿𝑦 𝛿𝑧
[𝑣 × 𝑣 ] = |𝑣𝑥 𝑣𝑦 𝑣𝑧 | = (𝑣𝑦 𝑣𝑧 − 𝑣𝑧 𝑣𝑦 )𝛿𝑥 + (𝑣𝑧 𝑣𝑥 − 𝑣𝑥 𝑣𝑧 )𝛿𝑦 + (𝑣𝑥 𝑣𝑦 − 𝑣𝑦 𝑣𝑥 )𝛿𝑧 = 0.
𝑣𝑥 𝑣𝑦 𝑣𝑧

Uno puede llegar rápidamente a la respuesta anterior si se acuerda que geométricamente:

[𝒗
⃗ ×𝒗⃗ ] = ‖𝑣 ‖‖𝑣 ‖𝑠𝑒𝑛𝛼 = 0. De otra parte, se puede notar que el término del medio del
determinante ((𝒗𝒛 𝒗𝒙 − 𝒗𝒙 𝒗𝒛 )𝛿𝑦 ) aparece con los signos invertidos en comparación a como
se acostumbra. Esto se hace para poder expresar la ecuación como una sumatoria.

Una operación novedosa entre vectores, para los estudiantes de este curso, es el producto
diádico entre vectores que resulta en un tensor de orden 2. Usando de nuevo el vector
velocidad, 𝒗
⃗ = 𝒗𝒙 𝛿𝑥 + 𝒗𝒚 𝛿𝑦 + 𝒗𝒛 𝛿𝑧 , veamos cómo se hace un producto diádico.

Producto diádico

Sea 𝑣 = 𝑣𝑥 𝛿𝑥 + 𝑣𝑦 𝛿𝑦 + 𝑣𝑧 𝛿𝑧 , hallamos el producto diádico de 𝒗

⃗ consigo mismo, {𝒗 ⃗ }. Una
⃗ 𝒗
manera de hacer este producto es ordenar el primer vector en forma de columna y el
segundo en forma de fila y hacer la multiplicación respectiva, así:

𝑣𝑥 𝑣𝑥 𝑣𝑥 𝑣𝑥 𝑣𝑦 𝑣𝑥 𝑣𝑧
{𝑣 𝑣} = (𝑣𝑦 ) (𝑣𝑥 𝑣𝑦 𝑣𝑧 ) = ( 𝑣𝑦 𝑣𝑥 𝑣𝑦 𝑣𝑦 𝑣𝑦 𝑣𝑧 )
𝑣𝑧 𝑣𝑧 𝑣𝑥 𝑣𝑧 𝑣𝑦 𝑣𝑧 𝑣𝑧

El procedimiento de la multiplicación consiste en que la primera fila del tensor resulta de

multiplicar el primer componente del vector columna con todos los componentes del vector
fila:

𝑣𝑥 𝑣𝑥 𝑣𝑥 𝑣𝑥 𝑣𝑦 𝑣𝑥 𝑣𝑧
( ) (𝑣𝑥 𝑣𝑦 𝑣𝑧 ) = ( )

La segunda fila del tensor es el producto del segundo componente del vector columna con
los componentes del vector fila:

21
 (𝑣𝑦 ) (𝑣𝑥 𝑣𝑦 𝑣𝑧 ) = ( 𝑣𝑦 𝑣𝑥 𝑣𝑦 𝑣𝑦 𝑣𝑦 𝑣𝑧 )

Y, la tercera fila del tensor es el producto del tercer componente del vector columna con
los componentes del vector fila:

( ) (𝑣𝑥 𝑣𝑦 𝑣𝑧 ) = ( )
𝑣𝑧 𝑣𝑧 𝑣𝑥 𝑣𝑧 𝑣𝑦 𝑣𝑧 𝑣𝑧

Operadores diferenciales involucrados en las ecuaciones de cambio y de transporte o

constitutivas
En lo que sigue, vamos a discutir los principales operadores diferenciales que aparecen en las
ecuaciones de cambio y en las llamadas ecuaciones constitutivas o de transporte.
Puntualmente, discutiremos los operadores: derivada parcial, nabla, divergencia, laplaciano,
rotacional, gradiente, convectivo y derivada material o sustancial
𝝏∎
Operador derivada parcial ( 𝝏𝝌 ). Es un operador escalar que mide la tasa de cambio de un
campo (∎) dado con respecto a una variable independiente. Aquí, estas últimas son las
direcciones en el espacio, e.g., (𝑥, 𝑦, 𝑧), y el tiempo, t. Como ejemplo, supongamos un campo
𝜕𝑓 𝜕𝑓 𝜕𝑓
escalar arbitrario: 𝑓(𝑥, 𝑦, 𝑧), sus derivadas parciales con respecto a (𝑥, 𝑦, 𝑧, 𝑡) son: 𝜕𝑥 , 𝜕𝑦 , 𝜕𝑧 y
0.

Operador nabla (𝛁 ⃗ ). Este operador es un vector que actúa sobre un campo y mide la tasa de
cambio de éste en cada una de las direcciones que se consideran de acuerdo con el marco de
referencia adoptado. En concordancia y como dedujimos en la primera clase, la dimensión de
este operador es 1/L. Para los alcances del curso, el operador nabla actúa sobre escalares y
vectores. Además, al ser un vector, puede aparecer en multiplicaciones del tipo producto punto
y cruz. De hecho, el producto punto de nabla con otro vector define un nuevo operador llamado
divergencia mientras que el producto cruz de nabla con otro vector define el operador
rotacional. En cuanto al producto diádico de nabla con un campo, éste se denomina gradiente.
Para los alcances del curso, nos interesarán solamente los gradientes de campos escalares y
vectoriales. El gradiente de un campo escalar resulta en un campo vectorial y el gradiente de
un campo vectorial en un tensor de orden 2, de acuerdo con la tabla 4.

⃗ ) = 𝜕 𝛿𝑥 + 𝜕 𝛿𝑦 + 𝜕 𝛿𝑧 .
En coordenadas cartesianas, el operador nabla se define como: (∇ 𝜕𝑥 𝜕𝑦 𝜕𝑧

⃗ ∙ ∎). Es el producto punto entre nabla y otro vector o un tensor de

Operador divergencia (𝛁
segundo orden. Sea cual sea el caso, el operador divergencia mide cómo es la tasa de variación
de un campo vectorial o tensorial en cada una de las direcciones del espacio.

En coordenadas cartesianas, la divergencia de un vector velocidad, 𝑣 = 𝑣𝑥 𝛿𝑥 + 𝑣𝑦 𝛿𝑦 + 𝑣𝑧 𝛿𝑧 , es

𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧
⃗ ∙ 𝑣) =
el escalar: (∇ + + .
𝜕𝑥 𝜕𝑦 𝜕𝑧

22
En coordenadas cartesianas, la divergencia de un tensor como el tensor de esfuerzos viscosos
⃗ ∙ 𝜏̿] = (𝜕𝜏𝑥𝑥 + 𝜕𝜏𝑦𝑥 + 𝜕𝜏𝑧𝑥 ) 𝛿𝑥 + (𝜕𝜏𝑥𝑦 + 𝜕𝜏𝑦𝑦 + 𝜕𝜏𝑧𝑦 ) 𝛿𝑦 + (𝜕𝜏𝑥𝑧 + 𝜕𝜏𝑦𝑧 + 𝜕𝜏𝑧𝑧) 𝛿𝑥 .
es el vector: [∇ 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧

Operador laplaciano (𝛁 𝟐 ). El operador laplaciano es un operador escalar que corresponde al

producto punto de nabla consigo mismo. Así, en coordenadas cartesianas:

𝜕2 𝜕2 𝜕2
⃗ ∙∇
∇2 = (∇ ⃗)= + +
𝜕𝑥 2 𝜕𝑦 2 𝜕𝑧 2

El laplaciano puede aplicarse a campos escalares, vectoriales y tensoriales. Aquí, nos interesa
su aplicación a escalares y vectores, 𝛁 𝟐 ∎. Como ilustración, podemos mostrar su aplicación
sobre un vector velocidad, 𝑣 :

𝜕2𝑣 𝜕2𝑣 𝜕2𝑣

⃗ ∙∇
[∇2 𝑣 ] = (∇ ⃗ )𝑣 = + +
𝜕𝑥 2 𝜕𝑦 2 𝜕𝑧 2

En general, el operador laplaciano sirve para medir los cambios de curvatura de un campo
dado.

Operador rotacional (𝛁⃗ × ∎). Es el producto cruz de nabla con un campo vectorial o tensorial.
Para el curso de Fenómenos de Transporte, nos basta con hacer el rotacional de un vector.
Para definir este último, podemos volver a usar la regla del determinante nuevamente para
coordenadas cartesianas:

𝛿𝑥 𝛿𝑦 𝛿𝑧
⃗ × 𝑣 ] = || 𝜕 𝜕 𝜕 | = 𝜕𝑣𝑧 − 𝜕𝑣𝑦 𝛿 + (𝜕𝑣𝑥 − 𝜕𝑣𝑧 ) 𝛿 + 𝜕𝑣𝑦 − 𝜕𝑣𝑥 𝛿
[∇ | ( 𝜕𝑦 ) 𝑥 ( )
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑧 𝜕𝑧 𝜕𝑥 𝑦 𝜕𝑥 𝜕𝑦 𝑧
𝑣𝑥 𝑣𝑦 𝑣𝑧

Como su nombre lo indica, el operador rotacional mide la tendencia de un campo vectorial a

rotar.

Operador gradiente (𝛁⃗ ∎). Es el producto diádico de nabla con un campo dado. Para este curso
son de interés los gradientes de campos escalares y vectoriales. Los primeros resultan en un
vector y los segundos en un tensor de segundo orden. Vamos a sacarle el gradiente a un
campo escalar 𝑓(𝑥, 𝑦, 𝑧, 𝑡) y a un vector velocidad 𝑣 :
𝜕𝑓 𝜕𝑓 𝜕𝑓
⃗ 𝑓] =
[∇ 𝛿𝑥 + 𝛿𝑦 + 𝛿
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝑧
𝜕 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧
𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥
𝜕 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧
⃗ 𝑣} =
{∇ (𝑣𝑥 𝑣𝑦 𝑣𝑧 ) =
𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦
𝜕 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧
( 𝜕𝑧 ) ( 𝜕𝑧 𝜕𝑧 𝜕𝑧 )
Operador convectivo (𝒗 ⃗ ∙ ⃗𝛁)∎. Este operador escalar resulta de hacer el producto punto entre
el vector velocidad y nabla:

23
 𝜕∎ 𝜕∎ 𝜕∎
(𝑣 ∙ ⃗∇)∎ = 𝑣𝑥 + 𝑣𝑦 + 𝑣𝑧
𝜕𝑥 𝜕𝑦 𝜕𝑧
Durante el curso, encontraremos el operador convectivo aplicado tanto a campos escalares
como vectoriales, así:
𝜕𝑓 𝜕𝑓 𝜕𝑓
⃗ )𝑓 = 𝑣𝑥
(𝑣 ∙ ∇ + 𝑣𝑦 + 𝑣𝑧
𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝑣 𝜕𝑣 𝜕𝑣
(𝑣 ∙ ⃗∇)𝑣 = 𝑣𝑥 + 𝑣𝑦 + 𝑣𝑧
𝜕𝑥 𝜕𝑦 𝜕𝑧
El hecho que el operador convectivo sea diferente del operador divergencia aplicado al vector
velocidad contraría el hecho que el producto punto entre dos vectores es conmutativo. Esto es
debido a la naturaleza misma del vector operador nabla.

𝑫∎
Operador derivada sustancial o material ( 𝑫𝒕 ). Este operador escalar es la suma algebraica del
operador derivada parcial con respecto al tiempo (llamada en mecánica de fluidos, derivada
local) con la derivada convectiva:
𝐷∎ 𝜕∎
= ⃗ )∎
+ (𝑣 ∙ ∇
𝐷𝑡 𝜕𝑡
Durante el curso, aplicaremos este operador a campos escalares y vectoriales.
Desde el punto de vista físico, este operador describe la tasa de cambio de una cantidad de
transporte -∎- (e.g., masa, cantidad de movimiento, energía) ligada a un cierto elemento de
transporte (e.g., un área o un volumen) que a su vez está sujeto a un campo de velocidad
dependiente del tiempo y el sistema de coordenadas; e.g., 𝒗 ⃗ = 𝒇(𝒙, 𝒚, 𝒛, 𝒕). En muchos libros
de matemáticas y mecánica de fluidos, la derivada material es simplemente la derivada total,
𝑑
, de un campo dado.
𝑑𝜒

⃗ = 𝑎(𝑥 2 − 𝑦 2 )𝛿𝑥 +
Ejemplo de aplicación de las operaciones descritas. El vector velocidad 𝒗
(−2𝑎𝑥𝑦)𝛿𝑦 describe un flujo ideal y aparecerá varias veces durante el curso. Aquí, vamos a
usarlo para ilustrar las operaciones descritas durante la clase.

(i) Producto punto, (𝒗 ⃗ ):

⃗ ∙𝒗

(𝒗
⃗ ∙𝒗⃗ ) = ‖𝑣 ‖2 = 𝑎2 (𝑥 2 − 𝑦 2 )2 + (−2𝑎𝑥𝑦)2 = 𝑎2 (𝑥 4 − 2𝑥 2 𝑦 2 + 𝑦 4 ) + 4𝑎2 𝑥 2 𝑦 2 =
𝑎2 𝑥 4 + 2𝑎2 𝑥 2 𝑦 2 + 𝑎2 𝑦 4 = 𝑎2 (𝑥 4 + 2𝑥 2 𝑦 2 + 𝑦 4 ) = 𝑎2 (𝑥 2 + 𝑦 2 )2.

(ii) Producto cruz, (𝒗 ⃗ ) = 0.

⃗ ∙𝒗

(iii) Producto diádico, {𝑣 𝑣 }:

𝑣𝑥 𝑎(𝑥 2 − 𝑦 2 ) ∗ 𝑎(𝑥 2 − 𝑦 2 ) 𝑎(𝑥 2 − 𝑦 2 ) ∗ (−2𝑎𝑥𝑦) 𝑎(𝑥 2 − 𝑦 2 ) ∗ 0
{𝑣 𝑣} = (𝑣𝑦 ) (𝑣𝑥 𝑣𝑦 𝑣𝑧 ) = ( (−2𝑎𝑥𝑦) ∗ 𝑎(𝑥 2 − 𝑦 2 ) (−2𝑎𝑥𝑦) ∗ (−2𝑎𝑥𝑦) (−2𝑎𝑥𝑦) ∗ 0 ) =
𝑣𝑧 0 ∗ 𝑎(𝑥 2 − 𝑦 2 ) 0 ∗ (−2𝑎𝑥𝑦) 0∗0

24
 2
(𝑎(𝑥 2 − 𝑦 2 )) −2𝑎2 𝑥𝑦(𝑥 2 − 𝑦 2 ) 0
= ( −2𝑎2 𝑥𝑦(𝑥 2 − 𝑦 2 ) 4𝑎2 𝑥 2 𝑦 2 0)
0 0 0

⃗
𝜕𝑣
(iv) Derivadas parciales, 𝜕𝜒:

𝜕𝑣 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕
= 𝛿𝑥 + 𝛿𝑦 = (𝑎(𝑥 2 − 𝑦 2 )𝛿𝑥 + (−2𝑎𝑥𝑦)𝛿𝑦 )
𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥
𝜕 𝜕
= (𝑎(𝑥 2 − 𝑦 2 )𝛿𝑥 ) + ((−2𝑎𝑥𝑦)𝛿𝑦 ) = 2𝑎𝑥𝛿𝑥 − 2𝑎𝑦𝛿𝑦
𝜕𝑥 𝜕𝑥

𝜕𝑣 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕
= 𝛿𝑥 + 𝛿𝑦 = (𝑎(𝑥 2 − 𝑦 2 )𝛿𝑥 + (−2𝑎𝑥𝑦)𝛿𝑦 )
𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦
𝜕 𝜕
= (𝑎(𝑥 2 − 𝑦 2 )𝛿𝑥 ) + ((−2𝑎𝑥𝑦)𝛿𝑦 ) = −2𝑎𝑦𝛿𝑥 − 2𝑎𝑥𝛿𝑦
𝜕𝑦 𝜕𝑦
𝜕𝑣 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕
= 𝛿𝑥 + 𝛿𝑦 = (𝑎(𝑥 2 − 𝑦 2 )𝛿𝑥 + (−2𝑎𝑥𝑦)𝛿𝑦 ) = 0
𝜕𝑧 𝜕𝑧 𝜕𝑧 𝜕𝑧

𝜕𝑣 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕
= 𝛿𝑥 + 𝛿𝑦 = (𝑎(𝑥 2 − 𝑦 2 )𝛿𝑥 + (−2𝑎𝑥𝑦)𝛿𝑦 ) = 0
𝜕𝑡 𝜕𝑡 𝜕𝑡 𝜕𝑡

⃗ ∙ 𝑣 ) = 𝜕𝑣𝑥 + 𝜕𝑣𝑦 + 𝜕𝑣𝑧:

(v) Divergencia, (∇
𝜕𝑥 𝜕𝑦 𝜕𝑧

En el inciso (iv) ya se hicieron las derivadas parciales indicadas, por tanto:

⃗ ∙ 𝑣 ) = 2𝑎𝑥 − 2𝑎𝑥 + 0 = 0
(∇
𝜕2 𝑣
⃗ 𝜕2 𝑣
⃗ 𝜕2 𝑣
⃗
⃗ ∙∇
(vi) Laplaciano, ∇2 𝑣 = (∇ ⃗ )𝑣 = 2 + 2 + 2
𝜕𝑥 𝜕𝑦 𝜕𝑧

En el inciso (iv) se calcularon las primeras derivadas del vector velocidad y las usaremos para
calcular el laplaciano, así:

𝜕2𝑣 𝜕 𝜕𝑣 𝜕
2
= ( )= (2𝑎𝑥𝛿𝑥 − 2𝑎𝑦𝛿𝑦 ) = 2𝑎𝛿𝑥
𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥

𝜕2𝑣 𝜕 𝜕𝑣 𝜕
2
= ( )= (−2𝑎𝑦𝛿𝑥 − 2𝑎𝑥𝛿𝑦 ) = −2𝑎𝛿𝑥
𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑦

𝜕2𝑣 𝜕 𝜕𝑣 𝜕
2
= ( ) = (0) = 0
𝜕𝑧 𝜕𝑧 𝜕𝑧 𝜕𝑧

Por tanto, ∇2 𝑣 = 2𝑎𝛿𝑥 − 2𝑎𝛿𝑥 = 0.

⃗ × 𝑣 ] = (𝜕𝑣𝑧 − 𝜕𝑣𝑦 ) 𝛿𝑥 + (𝜕𝑣𝑥 − 𝜕𝑣𝑧) 𝛿𝑦 + (𝜕𝑣𝑦 − 𝜕𝑣𝑥 ) 𝛿𝑧

(vii) Rotacional, [∇ 𝜕𝑦 𝜕𝑧 𝜕𝑧 𝜕𝑥 𝜕𝑥 𝜕𝑦

25
De nuevo, las derivadas parciales requeridas para el cálculo ya se hicieron en el inciso (iv) y
podemos usarlas directamente:

⃗ × 𝑣 ] = (0 − 0)𝛿𝑥 + (0 − 0)𝛿𝑦 + (−2𝑎𝑦 − (−2𝑎𝑦))𝛿𝑧 = 0

[∇
𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧
𝜕𝑥 𝜕𝑥 𝜕𝑥
𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧
⃗ 𝑣} =
(viii) Gradiente, {∇ .
𝜕𝑦 𝜕𝑦 𝜕𝑦
𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧
( 𝜕𝑧
𝜕𝑧 𝜕𝑧 )
Conocidas todas las derivadas parciales (inciso (iv)), el gradiente del campo de velocidad
utilizado es:
2𝑎𝑥 −2𝑎𝑦 0
⃗ 𝑣 } = ( −2𝑎𝑦
{∇ −2𝑎𝑥 0)
0 0 0
⃗
𝜕𝑣 ⃗
𝜕𝑣 ⃗
𝜕𝑣
(ix) Derivada convectiva, (𝑣 ∙ ⃗∇)𝑣 = 𝑣𝑥 𝜕𝑥 + 𝑣𝑦 𝜕𝑦 + 𝑣𝑧 𝜕𝑧
Nuevamente, se requieren las derivadas parciales del inciso (iv):

⃗
𝜕𝑣 ⃗
𝜕𝑣 ⃗
𝜕𝑣
𝜕𝑥
= 2𝑎𝑥𝛿𝑥 − 2𝑎𝑦𝛿𝑦 ; 𝜕𝑦
= −2𝑎𝑦𝛿𝑥 − 2𝑎𝑥𝛿𝑦 ; 𝜕𝑧
= 0.

De este modo,

𝜕𝑣
(𝑣 ∙ ⃗∇)𝑣 = 𝑎(𝑥 2 − 𝑦 2 )(−2𝑎𝑦𝛿𝑥 − 2𝑎𝑥𝛿𝑦 ) + (−2𝑎𝑥𝑦)(−2𝑎𝑦𝛿𝑥 − 2𝑎𝑥𝛿𝑦 ) +0=
𝜕𝑦

= (−2𝑎2 (𝑦 3 − 2𝑥𝑦 2 + 𝑥 2 𝑦))𝛿𝑥 + (−2𝑎2 (𝑥 3 − 2𝑥 2 𝑦 − 𝑥𝑦 2 ))𝛿𝑦

⃗
𝐷𝑣 ⃗
𝜕𝑣
(x) Derivada material, = ⃗ )𝑣.
+ (𝑣 ∙ ∇
𝐷𝑡 𝜕𝑡

Para el cálculo, usamos los resultados de incisos (iv) y (ix)., así:

𝐷𝑣 𝜕𝑣
= + (𝑣 ∙ ⃗∇)𝑣 = 0 + (−2𝑎2 (𝑦 3 − 2𝑥𝑦 2 + 𝑥 2 𝑦))𝛿𝑥 + (−2𝑎2 (𝑥 3 − 2𝑥 2 𝑦 − 𝑥𝑦 2 ))𝛿𝑦 =
𝐷𝑡 𝜕𝑡
= (−2𝑎2 (𝑦 3 − 2𝑥𝑦 2 + 𝑥 2 𝑦))𝛿𝑥 + (−2𝑎2 (𝑥 3 − 2𝑥 2 𝑦 − 𝑥𝑦 2 ))𝛿𝑦

26
Notas de estudio
La siguiente tabla muestra los tensores que usaremos durante el curso.

Tabla 6. Descripción de los tensores que se emplearán durante el curso. En la notación

adoptada (Basada en BSL), las filas se representan en el primer índice de cada componente del
tensor mientras que las columnas están representadas en el segundo subíndice de cada
componente.
Nombre y símbolo Componentes* Aplicación Dimensiones
• Base tensorial para representar
cualquier tensor de orden 2 de
manera gráfica o geométrica. El
tensor de diadas unitarias se puede
representar usando el delta de
Tensor de díadas 𝛿11 𝛿12 𝛿13 Kronecker usando la notación de
unitarias (𝛿21 𝛿22 𝛿23 ) Einstein como 𝛿𝑖𝑗 . Adimensional
̿
𝜹 𝛿31 𝛿32 𝛿33 • Caso particular, tensor unitario:
1 0 0
(0 1 0)
0 0 1

Aquí, 𝛿𝑖𝑗 = 0, si j ≠ i y 𝛿𝑖𝑗 = 1, si j = i.

Tensor de • Ley de Newton de la viscosidad

esfuerzos 𝜏11 𝜏12 𝜏13 • Principio de conservación del 𝐹𝑢𝑒𝑟𝑧𝑎 𝑀
cortantes o (𝜏21 𝜏22 𝜏23 ) momento lineal: 2da ley de Newton. [=]
viscosos 𝜏31 𝜏32 𝜏33 Á𝑟𝑒𝑎 𝐿 × 𝑇2
• Principio de conservación de la
𝝉̿ energía: 1ra ley de la termodinámica.

Tensor de 𝜋11 𝜋12 𝜋13 • Principio de conservación del

esfuerzos momento lineal: 2da ley de Newton. 𝐹𝑢𝑒𝑟𝑧𝑎 𝑀
moleculares (𝜋21 𝜋22 𝜋23 ) [=]
𝜋31 𝜋32 𝜋33 • Principio de conservación de la Á𝑟𝑒𝑎 𝐿 × 𝑇2
𝝅
̿ = 𝑝𝜹̿ + 𝝉̿ energía: 1ra ley de la termodinámica
Tensor de flujo
combinado de 𝜙11 𝜙12 𝜙13 • Principio de conservación del
cantidad de momento lineal: 2da ley de Newton. 𝐹𝑢𝑒𝑟𝑧𝑎 𝑀
(𝜙21 𝜙22 𝜙23 ) [=]
movimiento • Principio de conservación de la Á𝑟𝑒𝑎 𝐿 × 𝑇2
𝜙31 𝜙32 𝜙33
̿ =𝝅 energía: 1ra ley de la termodinámica
𝝓 ̿ + 𝜌𝑣𝑣

𝐸𝑛𝑒𝑟𝑔í𝑎 × 𝐷𝑖𝑠𝑡𝑎𝑛𝑐𝑖𝑎
Tensor de 𝜅11 𝜅12 𝜅13 • Conducción de calor en medios
conductividad anisotrópicos Á𝑟𝑒𝑎 × 𝐶𝑎𝑚𝑏𝑖𝑜 𝑑𝑒 𝑡𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎
(𝜅21 𝜅22 𝜅23 )
térmica 𝜅31 𝜅32 𝜅33 • Principio de conservación de la
𝑀×𝐿
̿
𝜿 energía: 1ra ley de la termodinámica [=]
𝑇3 × Θ

Á𝑟𝑒𝑎
∆𝐴𝐵 11 ∆𝐴𝐵 12 ∆𝐴𝐵 13 𝑇𝑖𝑒𝑚𝑝𝑜
Tensor de • Ley de Fick en medios anisotrópicos.
(∆𝐴𝐵 21 ∆𝐴𝐵 22 ∆𝐴𝐵 23 )
difusividad ̿̿̿̿̿
∆𝑨𝑩 • Principio de conservación de la masa
∆𝐴𝐵 31 ∆𝐴𝐵 32 ∆𝐴𝐵 33 𝐿2
[=]
𝑇

La Tabla 7 muestra algunos de los productos en los que participa nabla.

27
Tabla 7. Productos comunes del gradiente que se usan durante el curso. *La notación
utilizada corresponde a la llamada notación de Einstein o de Gibbs.
Nombre y naturaleza del
Operación Resultado*
producto
Gradiente de un campo escalar 3
𝜕𝑝
(e.g., presión: p), producto ⃗ 𝑝]
[∇ ⃗ 𝑝] = ∑ ⃗⃗⃗
[∇ 𝛿𝑖
diádico 𝜕𝑥𝑖
𝑖=1
Gradiente de un campo vectorial 3 3
𝜕𝑣𝑗
(e.g. vector velocidad: 𝑣), ⃗ 𝑣}
{∇ ⃗ 𝑣} = ∑ ∑ ⃗⃗⃗
{∇ 𝛿𝑖 ⃗⃗⃗
𝛿𝑗
producto diádico 𝜕𝑥𝑖
𝑖=1 𝑗=1
Divergencia de un campo 3
𝜕𝑣𝑖
vectorial, producto punto. ⃗ ∙ 𝑣)
(∇ ⃗ ∙ 𝑣) = ∑
(∇
𝜕𝑥𝑖
𝑖=1
Divergencia de un campo 3 3
tensorial tipo I, (e.g. tensor de
⃗ ∙𝝅
[∇ ̿] [∇ ̿ ] = ∑ ⃗⃗⃗
⃗ ∙𝝅 𝛿𝑖 (∑ 𝜋𝑗𝑖 )
esfuerzos moleculares: 𝝅 ̿ ),
𝑖=1 𝑗=1
producto punto.
Divergencia de un campo
tensorial tipo II, (e.g. tensor de
transporte convectivo: 3 3
{𝜌𝑣𝑣 }), producto punto. 𝜕(𝜌𝑣𝑗 𝑣𝑖 )
⃗ ∙ 𝜌𝑣 𝑣]
[∇ ⃗ ∙ 𝜌𝑣𝑣] = ∑ ⃗⃗⃗
[∇ 𝛿𝑖 (∑ )
𝜕𝑥𝑗
𝑖=1 𝑗=1

Lapaciano de un campo vectorial 3 3

(e.g., vector velocidad, 𝑣), 𝜕 2 𝑣𝑖
⃗ ∙∇
[∇ ⃗ 𝑣 ] = ∇2 𝑣 ∇2 𝑣 = ∑ ⃗⃗⃗
𝛿𝑖 (∑ )
producto punto equivalente a un 𝜕𝑥𝑗2
𝑖=1 𝑗=1
producto diádico.

3 3 3
𝜕𝑣𝑘
[∇ ⃗ ] = ∑ ∑ ∑ 𝜀𝑖𝑗𝑘 ⃗⃗⃗
⃗ ×v 𝛿𝑖
𝜕𝑥𝑗
𝑖=1 𝑗=1 𝑘=1
Rotacional de un campo Donde,
vectorial (e.g., vector velocidad, ⃗ ×v
[∇ ⃗] 𝜀𝑖𝑗𝑘 = +1, si ijk = 123, 231 ó 312
𝑣), producto cruz.
𝜀𝑖𝑗𝑘 = −1, si ijk = 321, 132 ó 213
𝜀𝑖𝑗𝑘 = 0, si cualquier índice se repite, e.g.,
112.

Producto escalar entre dos

campos tensoriales (e.g., tensor 3 3
𝜕𝑣𝑗
de esfuerzos cortantes o (𝝉̿: ⃗∇𝑣) ⃗ 𝑣) = ∑ ∑ 𝜏𝑖𝑗
(𝝉̿: ∇
viscosos: 𝝉̿), producto doble 𝜕𝑥𝑖
𝑖=1 𝑗=1
punto.

28
Ejemplos adicionales de operaciones entre escalares, vectores y tensores.
Ejemplo II.2. La ley generalizada de Newton de la viscosidad es una ecuación constitutiva (o
ecuación de transporte) de carácter tensorial que expresa la relación entre el esfuerzo cortante
aplicado a un fluido y su tasa de deformación. La ecuación se expresa de la siguiente manera:
ǂ
𝝉̿ = −𝜇 {∇ ⃗ 𝑣 ) } + (2 𝜇 − 𝜅) (∇
⃗ 𝑣 + (∇ ̿
⃗ ∙ 𝑣 )𝜹 (EII.1)
3

Donde, µ y κ son las viscosidades dinámica (Newtoniana) y volumétrica (dilatacional),

respectivamente. Ambas son escalares y tienen dimensiones ML-1T-1. El símbolo ǂ indica la
transpuesta del tensor gradiente del vector velocidad. ̿
𝜹 es el tensor unitario.
Desarrolle componente por componente la ley generalizada de Newton de la viscosidad
efectuando las operaciones indicadas.
Solución. El problema planteado pide desarrollar explícitamente los componentes de la
ecuación EII-1. Para hacer esto, hay que reconocer: (i) esta ecuación expresa una relación entre
tensores de segundo orden. Por tanto, la ecuación tiene nueve componentes. (ii) La primera
expresión del lado derecho de la igualdad es la suma de dos tensores que son el producto
diádico de los vectores nabla y velocidad. (iii) Para el segundo término de la ecuación al lado
derecho de la igualdad, el carácter tensorial está dado por la realización del producto diádico
entre un escalar y el tensor unitario. (iv) Para expresar los componentes de la ecuación, se
deben usar los resultados presentados en la Tablas 6 y 7. (v) Las expresiones en la tabla 6 son
directamente aplicables a un sistema de coordenadas cartesianas. Por tanto, en lugar de usar
la notación numérica allí sugerida, haremos i = 1, 2, 3 = x, y, z.
Pasos de solución:

(1) Expresar 𝝉̿ de manera explícita:

𝜏𝑥𝑥 𝜏𝑥𝑦 𝜏𝑥𝑧

𝜏
𝝉̿ = ( 𝑦𝑥 𝜏𝑦𝑦 𝜏𝑦𝑧 ) (EII.2)
𝜏𝑧𝑥 𝜏𝑧𝑦 𝜏𝑧𝑧

⃗ 𝑣 } (Tabla 7):
(2) Calcular {∇
𝜕𝑣 𝜕𝑣𝑦
⃗ 𝑣 } = ∑𝑧𝑖=𝑥,𝑦 ∑𝑧𝑗=𝑥,𝑦 ⃗⃗⃗
{∇ 𝛿𝑖 ⃗⃗⃗
𝛿𝑗 𝜕𝑥𝑗 = ∑𝑧𝑖=𝑥,𝑦 ⃗⃗⃗ ⃗⃗⃗⃗𝑥 (𝜕𝑣𝑥 ) + ⃗⃗⃗⃗
𝛿𝑖 [𝛿 𝛿 𝑦 ( ⃗⃗⃗𝑧 (𝜕𝑣𝑧)] = {𝛿
)+𝛿 ⃗⃗⃗⃗𝑥 ⃗⃗⃗⃗ 𝜕𝑣
𝛿𝑥 ( 𝜕𝑥𝑥 ) +
𝑖 𝜕𝑥 𝑖 𝜕𝑥 𝑖 𝜕𝑥 𝑖
𝜕𝑣𝑦 𝜕𝑣 𝜕𝑣 𝜕𝑣𝑦 𝜕𝑣 𝜕𝑣 𝜕𝑣𝑦
⃗⃗⃗⃗
𝛿𝑥 ⃗⃗⃗⃗
𝛿𝑦 ( 𝜕𝑥 ) + ⃗⃗⃗⃗
𝛿𝑥 ⃗⃗⃗
𝛿𝑧 ( 𝜕𝑥𝑧) + ⃗⃗⃗⃗
𝛿𝑦 ⃗⃗⃗⃗
𝛿𝑥 ( 𝜕𝑦𝑥 ) + ⃗⃗⃗⃗
𝛿𝑦 ⃗⃗⃗⃗
𝛿𝑦 ( 𝜕𝑦 ) + ⃗⃗⃗⃗
𝛿𝑦 ⃗⃗⃗ ⃗⃗⃗𝑧 ⃗⃗⃗⃗
𝛿𝑧 ( 𝜕𝑦𝑧) + 𝛿 𝛿𝑥 ( 𝜕𝑧𝑥 ) + ⃗⃗⃗
𝛿𝑧 ⃗⃗⃗⃗
𝛿𝑦 ( 𝜕𝑧 ) +
⃗⃗⃗ 𝜕𝑣
𝛿𝑧 ⃗⃗⃗
𝛿𝑧 ( 𝜕𝑧𝑧 ) } (EII.3)

La expresión obtenida se puede representar en forma matricial así:

𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧
𝜕𝑥 𝜕𝑥 𝜕𝑥
𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑦
⃗ 𝑣} =
{∇ (EII.4)
𝜕𝑦 𝜕𝑦 𝜕𝑧
𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧
( 𝜕𝑧 𝜕𝑧 𝜕𝑧 )

29
 ⃗ 𝑣 }† usando la ecuación EII.4:
(3) Calcular {∇
𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕𝑣𝑥
𝜕𝑥 𝜕𝑦 𝜕𝑧
ǂ 𝜕𝑣𝑦 𝜕𝑣𝑦 𝜕𝑣𝑦
⃗ 𝑣} =
{∇ (EII.5)
𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝑣𝑧 𝜕𝑣𝑧 𝜕𝑣𝑧
( 𝜕𝑥 𝜕𝑦 𝜕𝑧 )

(4) Usando las ecuaciones EII.4 y EII.5, el primer término al lado derecho de la igualdad de la
ecuación EII.1 es:
𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧 𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕𝑣𝑥
𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑦 𝜕𝑣𝑦 𝜕𝑣𝑦 𝜕𝑣𝑦
−𝜇 𝜕𝑦 𝜕𝑦 𝜕𝑧
+ 𝜕𝑥 𝜕𝑦 𝜕𝑧
(EII.6a)
𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧 𝜕𝑣𝑧 𝜕𝑣𝑧 𝜕𝑣𝑧
{( 𝜕𝑧 𝜕𝑧 𝜕𝑧 ) ( 𝜕𝑥 𝜕𝑦 𝜕𝑧 )}

La suma de las dos tensores se hace término por término; tal como se hace con dos matrices,
así las cosas:
𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑥 𝜕𝑣𝑧 𝜕𝑣𝑥
2
𝜕𝑥 𝜕𝑥
+ 𝜕𝑦 𝜕𝑥
+ 𝜕𝑧
𝜕𝑣𝑥 𝜕𝑣 𝜕𝑣𝑦 𝜕𝑣𝑦 𝜕𝑣𝑦
−𝜇 + 𝑦 2 + (EII.6b)
𝜕𝑦 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑧
𝜕𝑣𝑥 𝜕𝑣𝑧 𝜕𝑣𝑦 𝜕𝑣 𝜕𝑣𝑧
{( 𝜕𝑧 + 𝜕𝑥 𝜕𝑧
+ 𝜕𝑦𝑧 2 𝜕𝑧 )}
(5) Para el segundo término a la derecha de la igualdad de la ecuación EII.1, se tiene que el
carácter tensorial de la expresión está dado por la multiplicación de dos escalares por ̿ 𝜹.
Primero, se puede hallar la divergencia del vector velocidad así:

⃗ ∙ 𝑣 ) = ∑𝑧𝑖=𝑥,𝑦 𝜕𝑣𝑖 = 𝜕𝑣𝑥 + 𝜕𝑣𝑦 + 𝜕𝑣𝑧 (EII.7)

(∇ 𝜕𝑥 𝜕𝑥 𝑖 𝜕𝑦 𝜕𝑧

Mientras tanto, el vector unitario es, por definición:

1 0 0
𝜹̿ = (0 1 0) (EII.8)
0 0 1
Por tanto, el término completo es:
1 0 0
2 ̿ = (2 𝜇 − 𝜅) (𝜕𝑣𝑥 + 𝜕𝑣𝑦 + 𝜕𝑣𝑧) (0
⃗ ∙ 𝑣 )𝜹
(3 𝜇 − 𝜅) (∇ 1 0) =
3 𝜕𝑥 𝜕𝑦 𝜕𝑧
0 0 1
2 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧
(3 𝜇 − 𝜅) ( 𝜕𝑥 + 𝜕𝑦 + 𝜕𝑧 ) 0 0
2 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧
0 ( 𝜇 − 𝜅) (
3 𝜕𝑥
+
𝜕𝑦
+
𝜕𝑧
) 0 (
2 𝜕𝑣 𝜕𝑣𝑦 𝜕𝑣𝑧
{( 0 0 (3 𝜇 − 𝜅) ( 𝜕𝑥𝑥 + 𝜕𝑦
+ )
𝜕𝑧 )}
EII.9)
(6) El uso de la ecuación EII.2 y la suma de las ecuaciones EII.6b y EII.9 resulta entonces en una
expresión explícita de la ecuación EII.1, representada en forma material así:

30
 𝜏𝑥𝑥 𝜏𝑥𝑦 𝜏𝑥𝑧
(𝜏𝑦𝑥 𝜏𝑦𝑦 𝜏𝑦𝑧 ) =
𝜏𝑧𝑥 𝜏𝑧𝑦 𝜏𝑧𝑧
𝜕𝑣𝑥 2 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧 𝜕𝑣𝑦 𝜕𝑣𝑥 𝜕𝑣𝑧 𝜕𝑣𝑥
−2𝜇 + ( 𝜇 − 𝜅) ( + + ) −𝜇 ( + ) −𝜇 ( + )
𝜕𝑥 3 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑥 𝜕𝑧
𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑦 2 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧 𝜕𝑣𝑦 𝜕𝑣𝑦
−𝜇 ( + ) −2𝜇 + ( 𝜇 − 𝜅) ( + + ) −𝜇 ( + )
𝜕𝑦 𝜕𝑥 𝜕𝑦 3 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑧 𝜕𝑧
𝜕𝑣𝑥 𝜕𝑣𝑧 𝜕𝑣𝑦 𝜕𝑣𝑧 𝜕𝑣𝑧 2 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧
( −𝜇 ( + ) −𝜇 ( + ) −2𝜇 + ( 𝜇 − 𝜅) ( + + ))
𝜕𝑧 𝜕𝑥 𝜕𝑧 𝜕𝑦 𝜕𝑧 3 𝜕𝑥 𝜕𝑦 𝜕𝑧
(EII.10)
La ecuación EII.10 se presenta en forma de tabla en el apéndice B de BSL, Tabla B1.
Ejemplo II.3. Forma alternativa de hacer un producto diádico entre dos vectores. Efectúe el
producto {∇⃗ 𝑣 } usando: (i) una técnica análoga a la que se hace para multiplicar un vector
columna con un vector fila; (ii) marcando las posiciones en la matriz resultante con los
productos diádicos de los vectores unitarios ⃗⃗⃗⃗
𝜹𝒊 .
Solución
(i) Para el método pedido en (i), se parte de la definición de producto externo que resulta en
el producto diádico de los dos vectores (ver:
https://en.wikipedia.org/wiki/Row_and_column_vectors):
𝑢𝑥 𝑢𝑥 𝑣𝑥 𝑢𝑥 𝑣𝑦 𝑢𝑥 𝑣𝑧
⃗ ⨂𝑣 = (𝑢𝑦 ) (𝑣𝑥
𝑢 𝑣𝑦 𝑣𝑧 ) = (𝑢𝑦 𝑣𝑥 𝑢𝑦 𝑣𝑦 𝑢𝑦 𝑣𝑧 ) (EII.11)
𝑢𝑧 𝑢𝑧 𝑣𝑥 𝑢𝑧 𝑣𝑦 𝑢𝑧 𝑣𝑧

⃗ 𝑣 }:
Aplicando la definición presentada en la ecuación (E2.3-1) al producto {∇
𝜕∎ 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧
𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑥
𝜕∎ 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑦
⃗ ⨂𝑣 =
∇ 𝜕𝑦
(𝑣𝑥 𝑣𝑦 𝑣𝑧 ) = (EII.12).
𝜕𝑦 𝜕𝑦 𝜕𝑧
𝜕∎ 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧
( 𝜕𝑧 ) ( 𝜕𝑧 𝜕𝑧 𝜕𝑧 )

Nótese que en la ecuación EII.12, el producto entre las componentes del vector gradiente y las
de la velocidad no corresponden a una multiplicación convencional sino a la acción del primero
sobre el segundo porque el vector gradiente es un operador. Por tanto, el símbolo ∎
simplemente se reemplaza por la componente del vector velocidad que corresponda.
(ii) En este caso, la posición de cada componente del tensor gradiente de la velocidad se marca
en la representación matricial así:

⃗⃗⃗⃗⃗
𝜹𝒙 ⃗⃗⃗⃗⃗
𝜹𝒙 ⃗⃗⃗⃗⃗
𝜹𝒚 ⃗⃗⃗⃗⃗
𝜹𝒙 ⃗⃗⃗⃗⃗
𝜹𝒛 ⃗⃗⃗⃗⃗
𝜹𝒙
(⃗⃗⃗⃗⃗
𝜹𝒙 ⃗⃗⃗⃗⃗
𝜹𝒚 ⃗⃗⃗⃗⃗
𝜹𝒚 ⃗⃗⃗⃗⃗
𝜹𝒚 ⃗⃗⃗⃗⃗
𝜹𝒚 ⃗⃗⃗⃗⃗
𝜹𝒛 ) (EII.13).
⃗⃗⃗⃗⃗
𝜹𝒙 ⃗⃗⃗⃗⃗
𝜹𝒛 ⃗⃗⃗⃗⃗
𝜹𝒙 ⃗⃗⃗⃗⃗
𝜹𝒛 ⃗⃗⃗⃗⃗
𝜹𝒛 ⃗⃗⃗⃗⃗
𝜹𝒛
El primer índice del producto diádico de los vectores unitario marca el denominador de la
derivada parcial del gradiente mientras que el segundo índice marca la componente del vector
velocidad, por tanto:

31
 ⃗⃗⃗⃗⃗
𝜹𝒙 ⃗⃗⃗⃗⃗
𝜹𝒙 ⃗⃗⃗⃗⃗
𝜹𝒙 ⃗⃗⃗⃗⃗
𝜹𝒚 ⃗⃗⃗⃗⃗
𝜹𝒙 ⃗⃗⃗⃗⃗
𝜹𝒛
𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧
𝜕𝑥 𝜕𝑥 𝜕𝑥
⃗⃗⃗⃗⃗
𝜹𝒚 ⃗⃗⃗⃗⃗
𝜹𝒙 ⃗⃗⃗⃗⃗
𝜹𝒚 ⃗⃗⃗⃗⃗
𝜹𝒚 ⃗⃗⃗⃗⃗
𝜹𝒚 ⃗⃗⃗⃗⃗
𝜹𝒛
⃗ 𝑣} =
{∇ 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧 (EII.14).
𝜕𝑦 𝜕𝑦 𝜕𝑦
⃗⃗⃗⃗⃗
𝜹𝒛 ⃗⃗⃗⃗⃗
𝜹𝒙 ⃗⃗⃗⃗⃗
𝜹𝒛 ⃗⃗⃗⃗⃗
𝜹𝒚 ⃗⃗⃗⃗⃗
𝜹𝒛 ⃗⃗⃗⃗⃗
𝜹𝒛
𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧
( 𝜕𝑧 𝜕𝑧 𝜕𝑧 )
Por supuesto, los productos diádicos de los vectores unitarios se pueden eliminar de la
representación final del tensor. En general, con la práctica, el lector podría escribir este
resultado rápida y mecánicamente.
Ejemplo II.4. La ecuación de conservación de la energía para un fluido que fluye en forma de
láminas o capas paralelas (i.e., flujo laminar) y cuyas densidad (𝜌[=] 𝑀⁄ 3 )y conductividad
𝐿
térmica (𝑘 [=] 𝑀𝐿⁄ 3 ) no cambian; i.e., el fluido es incompresible y su conductividad térmica
𝑇 Θ
es isentrópica, puede escribirse de la siguiente manera cuando no hay acumulación en el
sistema y si se ignoran los efectos de la viscosidad:

(𝑣 ∙ ⃗∇𝑇) = 𝛼∇2 𝑇 (EII.15)

2 2
Aquí, 𝛼 = 𝑘⁄𝜌𝐶 [=] 𝐿 ⁄Θ es la difusividad térmica. Donde, 𝐶𝑝 [=] 𝐿 𝑀⁄ 2 , es la capacidad
𝑝 𝑇 Θ
calorífica (ver: J.M. Smith, H.C. Van Ness, M.M. Abbott & M.T. Swihart, Introduction to Chemical
Engineering Thermodynamics, 8th edition, McGraw-Hill). Debe tenerse en cuenta que la ecuación
EII.15 también requiere que Cp sea constante.
Solución
La ecuación EII.15 es escalar. La operación a la derecha de la igualdad es el producto de la
difusividad térmica con el Laplaciano de la temperatura. El Laplaciano de un campo escalar f
arbitrario se define como:
𝜕2 𝑓 𝜕2 𝑓 𝜕2 𝑓
∇2 𝑓 = 𝜕𝑥 2
+ 𝜕𝑦2 + 𝜕𝑧2 (EII.16)

Así las cosas, para la parte derecha de la ecuación EII.15 se tiene:

𝜕2 𝑇 𝜕2 𝑇 𝜕2 𝑇
𝛼∇2 𝑇 = 𝛼 ( + + ) (EII.17)
𝜕𝑥 2 𝜕𝑦 2 𝜕𝑧 2

Por otro lado, la parte izquierda de la ecuación EII.15 puede desarrollarse así:
𝜕𝑇 𝜕𝑇 𝜕𝑇 𝜕𝑇 𝜕𝑇 𝜕𝑇
(𝑣 ∙ ⃗∇𝑇) = [𝑣𝑥 ⃗⃗⃗⃗
𝛿𝑥 + 𝑣𝑦 ⃗⃗⃗⃗
𝛿𝑦 + 𝑣𝑧 ⃗⃗⃗⃗
𝛿𝑦 ] ∙ [𝜕𝑥 ⃗⃗⃗⃗
𝛿𝑥 + 𝜕𝑦 ⃗⃗⃗⃗
𝛿𝑦 + 𝜕𝑧 ⃗⃗⃗⃗
𝛿𝑦 ] = 𝑣𝑥 𝜕𝑥 + 𝑣𝑦 𝜕𝑦 + 𝑣𝑧 𝜕𝑧 (EII.18).

Por tanto, la forma explícita de la ecuación EII.15 es:

𝜕𝑇 𝜕𝑇 𝜕𝑇 𝜕2 𝑇 𝜕2 𝑇 𝜕2 𝑇
𝑣𝑥 𝜕𝑥 + 𝑣𝑦 𝜕𝑦 + 𝑣𝑧 𝜕𝑧 = 𝛼 (𝜕𝑥 2 + 𝜕𝑦2 + 𝜕𝑧2 ) (EII.19).

Ejemplo II.5. Cuando un fluido se desplaza, la fricción intermolecular de sus constituyentes

causa la disipación de parte de su energía en calor. A este fenómeno se le denomina disipación

32
viscosa. Una manera de estimar la magnitud de la disipación viscosa corresponde a calcular el
producto doble punto entre el tensor de esfuerzos viscosos y el gradiente del vector velocidad:
⃗ 𝑣 ) (ecuación EII.20). Desarrolle explícitamente este producto.
(𝝉̿: ∇
Solución
De acuerdo con la Tabla 6, la operación indicada en EII.20 se puede calcular así:
𝜕𝑣
⃗ 𝑣 ) = ∑𝑧𝑖=𝑥,𝑦 ∑𝑧𝑗=𝑥,𝑦 𝜏𝑖𝑗 𝑖 (EII.21).
(𝝉̿: ∇ 𝜕𝑥 𝑗

𝜕𝑣𝑖 𝜕𝑣𝑖 𝜕𝑣 𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕𝑣𝑦

Entonces, ⃗ 𝑣 ) = ∑𝑧𝑖=𝑥,𝑦 𝜏𝑖𝑥
(𝝉̿: ∇ + 𝜏𝑖𝑦 + 𝜏𝑖𝑧 𝑖 = 𝜏𝑥𝑥 + 𝜏𝑥𝑦 + 𝜏𝑥𝑧 + 𝜏𝑦𝑥 +
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥
𝜕𝑣𝑦 𝜕𝑣𝑦 𝜕𝑣𝑧 𝜕𝑣𝑧 𝜕𝑣𝑧
𝜏𝑦𝑦 𝜕𝑦
+ 𝜏𝑦𝑧 𝜕𝑧
+ 𝜏𝑧𝑥 𝜕𝑥
+ 𝜏𝑧𝑦 𝜕𝑦 + 𝜏𝑧𝑧 𝜕𝑧 (EII.22).

La ecuación (EII.22) corresponde a la expresión (C) en la Tabla A.7-1 de BSL.

33
Ejercicios
1. Indique qué orden tienen los resultados de las siguientes operaciones y represéntelas de
acuerdo con ello.

a) 𝜂𝜸
̿ ; b) 𝜌𝒈
̿ ; c) 𝑣 ∙ 𝑤 ̿ 𝝉̿; f) 𝑣 × 𝑤
⃗⃗ ; d) 𝑣 ∙ 𝝉̿; e)∙ 𝝃: ⃗⃗ ; g) 𝜌𝑣 𝑣 ; h) 𝝉̿: 𝜌𝑣 𝑣

⃗ ×v
2. Se dice que un campo vectorial es irrotacional si se cumple que [∇ ⃗ ] = 0. Verifique si los
siguientes campos de velocidad son irrotacionales:
𝑥 2
⃗ = 𝐴 (1 − (𝐵) ) ⃗⃗⃗
(a) v 𝛿𝑧 ;

𝛼𝑦
1 1 𝑦 1
⃗ = 𝐴 (𝑒 𝛼 ( −
(b) v ) −𝑒𝛽 ( − )) ⃗⃗⃗⃗
𝛿𝑥 ;
𝛼 𝛼2 𝛼𝛽 𝛼2

∆𝑝𝑏2 2𝜇 𝜇−𝜇 𝑧 𝑧 2
(c) v
⃗ = (
2𝜇𝐿 𝜇+𝜇0
+ (𝜇+𝜇0 ) (𝑏) − (𝑏) ) ⃗⃗⃗⃗
𝛿𝑦 ;
0

(d) 𝑣𝑥 = 𝑥 + 𝑦 + 𝑧; 𝑣𝑦 = 𝑥 − 𝑦 + 𝑧; 𝑣𝑧 = 𝑥 + 𝑦 + 3;

(e) 𝑣𝑥 = 𝑥𝑦; 𝑣𝑦 = 𝑦𝑧; 𝑣𝑧 = 𝑦𝑧 + 𝑧 2

3. Verifique las siguientes identidades:

⃗ ∙ [∇
(a) (∇ ⃗ ×v
⃗ ]) = 0

⃗ × [∇
(b) [∇ ⃗ 𝜑]] = 0

4. Considere el siguiente tensor de esfuerzos cortantes y los siguientes vectores velocidad:

𝑥 0 0
𝝉̿ = (0 𝑦 1)
0 1 𝑧
⃗⃗⃗⃗𝑥 + 𝑏𝑥𝛿
𝑣 = 𝑏𝑦𝛿 ⃗⃗⃗⃗𝑦 ;

⃗⃗⃗⃗𝑥 + (𝑥 − 𝑦 + 𝑧)𝛿
⃗ = (𝑥 + 𝑦 + 𝑧)𝛿
𝑢 ⃗⃗⃗⃗𝑦 + (𝑥 + 𝑦 + 𝑎)𝛿
⃗⃗⃗𝑧

𝑤 ⃗⃗⃗⃗𝑥 + 𝑦(𝑥 + 𝑧)𝛿

⃗⃗ = 𝑥(𝑦 + 𝑧)𝛿 ⃗⃗⃗⃗𝑦 − ((𝑥 + 𝑦)𝑧 + 𝑧 2 )𝛿
⃗⃗⃗𝑧

5. Haga las siguientes operaciones:

(a) ∇2 𝑢
⃗ ; ∇2 𝑣 ; ∇2 𝑤
⃗⃗

⃗𝑢
(b) (𝝉̿: ∇ ⃗ 𝑣 ); (𝝉̿: ∇
⃗ ); (𝝉̿: ∇ ⃗𝑤⃗⃗ )

(c) [𝝉̿ ∙ 𝑢
⃗ ]; [𝝉̿ ∙ 𝑣 ]; [𝝉̿ ∙ 𝑤
⃗⃗ ]
6. Determine las componentes del tensor de esfuerzos viscosos que satisface la ley
generalizada de Newton de la viscosidad para los campos de velocidad 𝑢 ⃗⃗ .
⃗ ,𝑣 𝑦 𝑤

34
7. Halle el campo de velocidad para el cual 𝝉̿ satisface la Ley generalizada de Newton de la
viscosidad.
8. Desarrolle de manera explícita la ecuación de continuidad (balance de masa):
𝜕𝜌
⃗ ∙ 𝜌𝑣 ) = 0
+ (∇
𝜕𝑡
9. Compruebe si los campos de velocidad de los ejercicios (a) y (c) cumplen con la ecuación de
continuidad.
10. Desarrolle explícitamente las siguientes expresiones:
(a) Vector densidad de flujo (flux) másico respecto a ejes fijos para una especie i dentro de una
mezcla multicomponente:

𝑛𝑖 = 𝜌𝑤𝑖 𝑣 + ⃗𝑗𝑖
⃗⃗⃗
𝜌𝑖
Donde, 𝑤𝑖 = 𝜌
(fracción másica del componente i en la mezcla), ⃗𝑗𝑖 [=]𝑀𝐿−2 Θ−1 es el vector
densidad de flujo másico para la especie i con respecto al vector velocidad 𝑣 ; ⃗𝑗𝑖 = 𝜌𝑖 (𝑣
⃗⃗⃗𝑖 − 𝑣 ).
(b) Tensor combinado de densidad de flujo de cantidad de movimiento:
̿ = 𝜌𝑣 𝑣 + 𝝅
𝝓 ̿
(c) Vector densidad de flujo de energía
1
̂ + 𝑣 2 ) 𝑣 + 𝑞 + [𝝅
𝑒 = 𝜌 (𝑈 ̿ ∙ 𝑣 ].
2

(d) Ecuación de continuidad para el componente i de una mezcla:

𝐷𝑤𝑖
𝜌 ⃗ ∙ ⃗𝑗𝑖 ) + 𝑟𝑖
= −(∇
𝐷𝑡
Donde, 𝑟𝑖 [=]𝑀𝐿−3 Θ−1, velocidad de reacción de la especie i.
(e) Ecuación de cantidad de movimiento para una mezcla:
𝑛
𝐷𝑣
𝜌 ⃗ ∙ 𝜋̿] + ∑ 𝜌𝑖 ⃗⃗⃗
= −[∇ 𝑔𝑖
𝐷𝑡
𝑖=1

𝑔𝑖 es el vector de aceleración gravitacional para una especie i.

Donde, ⃗⃗⃗
(f) Ecuación de energía para una mezcla:
𝑛
𝐷 1
𝜌 (𝑈̂ + 𝑣 2 ) = −(∇
⃗ ∙ 𝑞 ) − (∇
⃗ ∙ [𝜋̿ ∙ 𝑣 ] + ∑( 𝑛⃗𝑖 ∙ ⃗⃗⃗
𝑔𝑖 )
𝐷𝑡 2
𝑖=1

Referencias: Calculus: Concepts & Context, Stewart; Cálculo, Tom Apóstol.

35