Métodos

numéricos para
transferencia
de calor y masa
Transferencia de calor II

Esta foto de Autor desconocido está bajo licencia CC

 Necesidades de predecir
A nuestro alrededor ocurren una gran cantidad de procesos, naturales o forzados, o una
mezcla de ambos, estos procesos afectan nuestras vidas en muchos sentidos, por tanto es
importante conocer y tener la capacidad de predecir estos fenómenos.

Es indispensable observar que los procesos que gobiernan la producción de potencia,

como ejemplo, son los mismos que gobiernan los procesos de calentamiento o
acondicionamiento de aire.

Las industrias química y metalúrgica utilizan componentes como hornos,

intercambiadores de calor, condensadores, reactores, etc. En el diseño de máquinas
eléctricas y circuitos electrónicos, la transferencia de calor es un factor importante.
La predicción nos permite conocer el comportamiento de uno o varios
elementos antes de ponerlos a funcionar en las condiciones reales no
controladas.

Para poder predecir un comportamiento necesitamos conocer el

fenómeno que gobierna, podemos de esta forma observar el
comportamiento de varias propuestas y decidir con seguridad y
eficiencia la mejor solución.
La predicción de una situación física dada consiste
en obtener los valores de las variables relevantes
que gobiernan los procesos de interés, por ejemplo
una descripción completa en una cámara de
combustión debe predecirse los valores de la
velocidad, la presión, la temperatura, la
concentración de especies químicas relevantes, etc.;
también pueden obtenerse los esfuerzos, el flujo de
calor, y flujo de masa.

La predicción debe establecer como estas cantidades

cambian en respuesta a los cambios de geometría,
flujos, propiedades de los fluidos, etc.
Métodos de predicción
Podemos mencionar dos métodos principales, uno, investigación experimental, la
otra, cálculo teórico. Veamos una breve comparación entre las dos:

Investigación experimental

Cálculo teórico
 Investigación
experimental

La información más confiable acerca de un

proceso físico se obtiene a través de la medición
directa. Una investigación experimental, involucra
equipos de escala completa que pueden utilizarse
para predecir como copias idénticas de los
equipos se desempeñarían bajo las mismas
condiciones

Estas pruebas de escala completa es en muchos

casos excesivamente caras y frecuentemente
imposibles. La alternativa es realizar
experimentos en modelos a escala.
 Investigación
experimental

La información resultante, debe ser extrapolada a

la escala completa, por regla general, hacer esto
no está disponible muy a menudo. Además, los
modelos a escala no siempre simulan todas las
características medidas en la escala completa;
frecuentemente características importantes,
como la combustión o la ebullición, son omitidas
en las pruebas a escala.

Esto reduce la utilidad de los resultados de la s

pruebas. Finalmente, hay que considerar que
pueden existir situaciones de difícil medición, y
los instrumentos de medición no están libres de
error.
 Cálculo
teórico
Para un proceso físico de interés, un modelo matemático
principalmente consiste en un juego de ecuaciones diferenciales. Si
los métodos de las matemáticas clásicas son usados para resolver
estas ecuaciones, entonces, solo una pequeña porción de todos los
fenómenos pueden ser resueltos.

Además, estas soluciones suelen contener series infinitas, funciones

especiales, ecuaciones trascendentales, etc. Donde la solución
numérica puede presentarse como una gran tarea.

Afortunadamente, el desarrollo de métodos numéricos y la

disponibilidad de computadoras cada vez más potentes nos permiten
resolver casi cualquier modelo matemático de cualquier problema
práctico.
 Ventajas del cálculo teórico
a) Bajo costo. La ventaja más importante de la predicción computacional, es el
bajo costo. En muchas aplicaciones, el costo del uso de máquina es, por
mucho, menor al costo de la investigación experimental correspondiente.
Además, los costos de los insumos para la experimentación tienden a
aumentar, mientras los costos de computadoras y su uso van en disminución.

b) Rapidez. La investigación computacional, puede ser desarrollada con rapidez.

Un diseñador, puede estudiar la implicación de cientos de diferentes
configuraciones para acercarse al diseño óptimo. Por el otro lado, es fácil
imaginar el largo tiempo que es necesario para experimentar con la misma
cantidad de prototipos que la investigación computacional.
 Ventajas del cálculo teórico
c) Información completa. Una solución en computadora, proporciona información
detallada y completa. Se puede proporcionar los valores de todas las variables
relevantes en la región de interés. Contrario a lo que sucedería en la experimentación,
que además presenta el inconveniente de que los medidores causan cambios en los
flujos y por tanto en las propias mediciones y más si los puntos de medición son
varios.

d) Capacidad para simular condiciones reales. En el cálculo teórico, condiciones

naturales pueden ser fácilmente simuladas. No es necesario recurrir a modelos en
escala. Para un programa de computadora, no es problema tener dimensiones grandes
o pequeñas, o temperaturas muy altas o muy pequeñas, o sustancias tóxicas o
inflamables, o procesos lentos o rápidos.
 Ventajas del cálculo teórico
e) Capacidad de simular condiciones ideales. Los métodos de predicción,
son algunas veces utilizados para estudiar fenómenos básicos, tanto como
aplicaciones de ingeniería compleja. En el estudio de fenómenos básicos,
se presta atención a pocos parámetros esenciales, y se eliminan todas las
características irrelevantes. Además muchas idealizaciones son deseables,
como por ejemplo, la bidimensionalidad, densidad constante, superficie
adiabática, etc. Son situaciones que experimentalmente es muy difícil de
obtener.
 Desventajas del cálculo
teórico
a) Un análisis en computadora depende de un modelo matemático
y del método numérico, contrario a la investigación
experimental la cual observa la realidad.

b) Todo modelo matemático tiene que ser validado, o sea los

resultados de una situación dada tiene que estar de acuerdo con
situaciones ya conocidas y resueltas. Una perfecta y satisfactoria
técnica numérica puede producir resultados sin valor si un
inadecuado modelo matemático es utilizado.
 a) Problemas para los cuales una adecuada
descripción matemática puede obtenerse
(calor en conducción, flujo laminar, flujos
turbulentos simples, capa límite, etc.)

Los problemas
prácticos los
b) Problemas para los cuales una adecuada
podemos descripción matemática no puede
establecerse (flujos turbulentos complejos,
dividir en dos ciertos flujos no Newtonianos, algunas
reacciones en combustiones violentas,
algunos flujos de dos fases, etc.)
 Desventajas
del grupo a) Si el objetivo de la predicción es bastante limitado (como
encontrar la caída de presión en un aparato muy complicado), la
solución por computadora puede no ser más barato que la
experimentación. Para problemas difíciles que involucren
geometrías complejas, nolinearidades fuertes, variaciones en las
propiedades del fluido sensibles, etc., una solución numérica
puede ser trabajosa de obtener y hasta excesivamente cara.

Esta foto de Autor desconocido está bajo licencia

CC BY-ND Finalmente, cuando el problema matemático admite más de una
solución, a veces es difícil determinar cual de la solución en
computadora corresponde a la realidad.

Las desventajas mencionadas tienden a disminuir a medida que

la investigación progresa.
 Desventajas
del grupo b) Las desventajas del grupo b), son todas las del
grupo a), adicionando la no certeza de que los
resultados obtenidos corresponden a la realidad,
en estos casos, algunos experimentos son
altamente deseables.

Avances en la investigación de métodos

matemáticas causan que casos del grupo b) se
muevan hacia el grupo a). Estas investigaciones
consisten en proponer un modelo, trabajar con
el en algunos casos y comparar los resultados
Esta foto de Autor desconocido está bajo licencia CC BY-ND
con datos experimentales.
La solución numérica de procesos de transferencia de
calor, flujo de fluidos, etc. Se puede comenzar cuando
Descripción
las leyes gobernantes de estos procesos han sido
expresados en forma matemática, generalmente en matemática de
forma de ecuaciones diferenciales.
En principio nuestro propósito es familiarizarnos con la
un fenómeno
forma y significado de estas ecuaciones. Observaremos
que las ecuaciones relevantes tienen una forma
común, como primer paso para formular un
físico.
procedimiento de solución general.
Individualmente las ecuaciones diferenciales expresan
cierto principio de conservación. Cada ecuación utiliza
cierta cantidad física como variable dependiente y con
ella se realiza un balance incluyendo los factores que
influencian la variable.
La variable dependiente en estas ecuaciones
diferenciales son generalmente propiedades
específicas.
Los términos de una ecuación diferencial del
tipo que nos interesa, denotan la influencia
sobre un volumen de control básico.
Suponga que J denota un flujo que influencia
a una variable dependiente típica φ.
Consideremos un volumen de control de
dimensiones dx, dy, dz.
El flujo Jx (que es la componente de J en la
dirección x) entra por un área de dimensiones
dyּdz, mientras el flujo saliente en la cara
opuesta es Jx+(∂Jx/∂ x)dx.
Siendo el flujo neto (∂Jx/∂x)dxּdyּdz sobre las
superficies de estas caras.
Considerando el flujo en todas las caras
tenemos un flujo neto por unidad de volumen
= (∂Jx/∂x)+(∂Jy/∂y) +(∂Jz/∂z) = div J = flujo de J
Otro ejemplo expresado en términos de
volumen unitario es la rapidez de
cambio ∂(ρφ)/∂t, ρφ significa la
cantidad de una propiedad extensiva
contenida en un volumen unitario.

Una ecuación diferencial es la

compilación de algunos términos,
donde cada uno implica un balance o
una conservación.
 Para un flujo estable a baja velocidad con
disipación debido a la viscosidad
despreciable, la ecuación de la energía puede
ponerse como

div(ρuh)=div(k gradT)+Sh
Ecuación de la
donde h es la entalpía específica, k es la
energía conductividad térmica, T es la temperatura, y
Sh es la rapidez con que se “genera” calor por
unidad de volumen. El término div(k grad T)
representa la influencia de la conducción de
calor dentro del fluido, de acuerdo a la ley de
conducción de Fourier.
Para gases ideales y para sólidos y líquidos podemos
poner
c grad T = grad h,

donde c es el calor específico a presión constante, por

tanto obtenemos
div(ρuh)=div(k/c gradh)+Sh

si c es constante, la relación entre T y h se simplifica

h=cT
Esta foto de Autor desconocido está bajo licencia CC BY

div(ρuT)=div(k/c gradT)+Sh/c

La situación de calor de conducción estable se obtiene

cuando u=0,
div(k gradT)+Sh=0
 La ecuación diferencial gobernante para la
conservación del momentum en una
dirección dada para un fluido Newtoniano
puede ser escrita en forma similar; sin
embargo la complicación surge al considerar
la ley de viscosidad de Stoke.
Ecuación de
∂(ρu)/∂t+div(ρuu)=div(μ gradu)- ∂P/∂x+Bx
momentum
μ :viscosidad
P : presión
Bx : fuerzas externas por unidad de volumen
 +∂(ρφ)/∂t+div(ρuφ)=div(Γ gradφ)+S

+Γ:coeficiente de difusión
Ecuación +S: término fuente
diferencial +Los cuatro términos de esta ecuación son :
+-         El término de estado inestable
general +-         El término de convección
+-         El término de difusión
+Y el término fuente
Para el planteamiento de las ecuaciones, haremos ciertas Estas consideraciones se
consideraciones para simplificar el análisis: hacen necesarias, en algunos
casos, porque sin ellas no se
1. Modelo bidimensional. obtiene mayor claridad pero
2. Flujo laminar. si mayor trabajo, como en el
caso de considerar modelo
3. Flujo incompresible. bidimensional; en otros
casos, considerar un fluido
4. Fluido Newtoniano. como newtoniano o que
5. Hipótesis de Boussinesq. cumpla la hipótesis de
Boussinesq es una buena
6. Trabajo de expansión – compresión despreciable. aproximación a la realidad
para los casos que nos van a
7. Fluido no participante a la radiación. interesar.
 Consideremos un volumen de control
Ecuación de infinitesimal con forma de paralelepípedo
rectangular de dimensiones dx, dy y dz, el
continuidad cual está siendo atravesado por un fluido
que se mueve con una velocidad con
componentes u, v y w en las direcciones x,
y y z respectivamente.
Por una de las caras del volumen de
control, ubicada en la dirección x y paralela
al plano y – z, el fluido entra a éste con una
velocidad con componente en x igual a ux.
En un intervalo de tiempo dt, ingresa al
volumen de control una cantidad de fluido
igual a
ρx dy dz ux dt
donde ρx significa la densidad del fluido en
la posición x.
Por una de las caras del volumen de
control, ubicada en la dirección x y paralela
al plano y – z, el fluido entra a éste con una
velocidad con componente en x igual a ux.
En un intervalo de tiempo dt, ingresa al
volumen de control una cantidad de fluido
igual a
ρx dy dz ux dt
donde ρx significa la densidad del fluido en
la posición x.
Por la cara opuesta, ubicada en la posición
x + dx, en el mismo intervalo de tiempo,
sale una cantidad de fluido igual a

-ρx+dx dy dz ux+dx dt
La acumulación de fluido, dentro del
volumen de control, durante este
intervalo de tiempo, debido a este flujo
es
(ux ρx- ux+dx ρx+dx )dy dz dt
y la rapidez de acumulación de fluido,
dentro del volumen de control, debido a
este flujo es
(ux ρx- ux+dx ρx+dx )dy dz

esta misma relación puede ponerse

como
[∂(ρu)/∂x]dx dy dz
Considerando los flujos en las tres direcciones (velocidades u, v y w), y
tomando en cuenta que dentro del volumen de control puede haber una
fuente o un sumidero de fluido, la ecuación general puede escribirse
como

Rapidez de acumulación de fluido =

= [∂(ρu)/∂x+∂(ρv)/∂y+∂(ρw)/∂z]dx dy dz+dmg/dt

donde el término dmg/dt es la rapidez de generación de masa dentro del volumen de control a través
de una fuente o un sumidero.
Al considerar el caso bidimensional, y suponiendo que no hay generación de masa, esta ecuación se
reduce a
Rapidez de acumulación de fluido = [∂(ρu)/∂x+∂(ρv)/∂y]dx dy