Universidad Tecnológica de Panamá

Centro Regional de Veraguas

Facultad de Ingeniería Mecánica
Departamento de Energía y Ambiente

PARCIAL II

Nombre: __Melany Flores________ Cédula: 9-759-151__

__Rosa Mosquera_______ 2-739-885__

Indicaciones: desarrolle en grupo máximo de 2 estudiantes cada problema. Sea ordenado,

objetivo y presente la resolución a computadora e impresa.
Fecha de entrega 05/12/2022

Problema N° 1
Tema: fluidos compresibles (valor 50 pts)
A continuación, presente 2 casos de estudios respecto a la aplicabilidad de las siguientes
temáticas:
1. Flujo unidimensional de gases con Fricción: Flujo de Fanno.
2. Flujo unidimensional de gases con transferencia de calor: Flujo de Rayleigh.
3. Problemas de aplicación de Flujo de Fanno y flujo de Rayleigh.
Justifiquen cada uno de estos con artículos científicos.

Problema N° 2
Tema: dinámica computarizada de fluidos (CFD) (valor 50 pts)
Presenten 1 caso de estudio su respectiva solución y aplicación considerando la temática Dinámica
de Fluidos Computacional o CFD por sus siglas en inglés Computational Fluid Dynamics la cual
representa el área de conocimiento que trata sobre la simulación numérica de flujos fluidos para
describir su transferencia de energía, masa y movimiento asociados a la dinámica de fluidos.
Justifiquen el caso propuesto con el uso de artículos científicos.

Nota: no es necesario utilizar el formato I+D. Sin embargo, es indispensable que cumpla
con el marco científico en cuanto a documento final.

MECÁNICA DE FLUIDOS II
MCs. GUILLERMO LÓPEZ TENORIO
 Fluidos compresibles y dinámica computarizada de
fluidos
Compressible fluids and computational fluid dynamics
Melany M. Flores G., Rosa Mosquera 1, Guillermo López2*
Facultad de Ingeniería Mecánica, Centro Regional de Veraguas, Universidad Tecnológica de Panamá
1
1
Facultad de Ingeniería Mecánica, Centro Regional de Veraguas, Departamento de Energía y Ambiente y
Ciencia de los Materiales, Universidad Tecnológica de Panamá

*Autor de correspondencia: guillermo.lopez2@utp.ac.pa

RESUMEN. Este documento hace referencia a dos problemas, el primero sobre flujos compresibles y el
segundo sobre dinámica computarizada de fluidos en el primer problema se aplicarán dos casos de estudio
con respecto a temas puntuales como son el flujo unidimensional de gases con fricción (Flujo de Fanno), flujo
unidimensional de gases con transferencia de calor (Flujo de Rayleigh) y problemas de aplicación sobre estos
tipos de flujos. Para el segundo problema se realizó un caso de estudio en donde se aplicó los conocimientos
adquiridos en el último laboratorio sobre la Dinámica de Fluidos Computacional, explicando la transferencia
de energía, masa y movimiento aplicados a los fluidos, basándose en la influencia del número de Reynolds
sobre las características de un flujo. Es imperante mencionar que para la resolución de ambos problemas se
utilizó Excel, el Software ANSYS y artículos científicos relacionados con los temas como justificación para la
resolución de estos.

Palabras clave. Flujo unidimensional, Reynolds, gases, energía, masa

ABSTRACT. This document refers to two problems, the first on compressible flows and the second on
computerized fluid dynamics. In the first problem, two case studies will be applied with respect to specific
topics such as the one-dimensional flow of gases with friction (Fanno Flow)., one-dimensional flow of gases
with heat transfer (Rayleigh Flow) and application problems on these types of flows. For the second problem,
a case study was carried out where the knowledge acquired in the last laboratory on Computational Fluid
Dynamics was applied, explaining the transfer of energy, mass and movement applied to fluids, based on the
influence of the Reynolds number. about the characteristics of a flow. It is imperative to mention that for the
resolution of both problems Excel, the ANSYS Software and scientific articles related to the topics were used
as justification for their resolution.

Keywords. One-dimensional flow, Reynolds, gases, energy, mass

1. Introducción
En el análisis del flujo compresible, es indispensable considerar las variaciones de la densidad y la
temperatura, es por ello por lo que se aplican en el estudio de este la ecuación de energía como la ecuación
de estado que son necesarias para la resolución de este tipo de flujo compresible. Además, se caracteriza
por medio del parámetro adimensional conocido como número de Mach. En función de valores con altas
velocidades aerodinámicas externas que suelen clasificarse como flujo subsónico que va de 0.3<M<0.8 y
flujo sónico donde M=1, existen otras categorías, pero para el presente documento solo se han de utilizar
las anteriormente mencionadas en la aplicación de los dos casos de estudio planteado en el primer
problema[1].
El flujo Fanno es un modelo que se utiliza en los estudio y diseño de tuberías o toberas, por lo general si
no hubiese la presencia de la fricción, las propiedades de este tipo de flujo permanecerían constantes en el
ducto, pues la fricción genera que un incremento de la entropía y por consiguiente el número de Mach es
considerado subsónico, ya que disminuye al aumentar la fricción.[2]
En cambio, el flujo Rayleigh es aquel flujo estacionario unidimensional que circula por el ducto y tiene
propiedades caloríficas constantes en la presencia de transferencia de calor y cuya fricción se desprecia
generalmente. Todas las consideraciones que se realicen sobre este tipo de flujo afectan las ecuación sobre
conservación de masa, cantidad de movimiento y energía, provocando un modelo matemático que
posibilita la determinación de propiedades como la temperatura, velocidad, densidad, número de Mach,
entre otras.
Es imperativo que al trabajar tanto el flujo unidimensional de gases con fricción (flujo Fanno) y Flujo
unidimensional de gases con transferencia de calor (flujo Rayleigh) se realicen ciertas consideraciones
sobre los mismo, pues las ecuaciones se pueden volver muy complejas y las consideraciones permite que
se pueda simplificar un poco el análisis.[1], [2]
Para el problema en CFD se simula un flujo a través de un cilindro circular, tomando en cuenta la
influencia del número de Reynolds sobre las características del flujo. Como ya se conoce que, si se trata de
un cilindro simétrico, el número de Re a medida que aumenta el flujo se irá separando detrás del cilindro,
lo que provoca un desprendimiento de vórtices, cabe destacar que este desprendimiento es inestable y para
ello con el programa ANSYS se logra crear una simulación dependiente del tiempo donde se puede
apreciar estos efectos.
2. Marco teórico

MECÁNICA DE FLUIDOS II
MCs. GUILLERMO LÓPEZ TENORIO
 En los gases se diferencia entre el flujo con volumen constante (incompresible) y el flujo con cambio de
volumen (compresible). Por lo general al hablar de flujos incompresibles de gases, se hace referencia a los
procesos de flujo donde se contemplan de la misma forma que en el caso del flujo de un líquido[3]. En
caso de cambios por incremento de la presión y temperatura del fluido gaseoso, no se pueden dejar de
mencionar las relaciones físicas entre la presión, temperatura y volumen. Este flujo se denomina
compresible. Por lo general se aplica para gases ideales y de ello se destaca que cuando la velocidad de un
flujo de gas es mayor que el número de Mach 0.3, es necesario considerarse compresible[4].

En este análisis se considera dos casos de estudio uno mediante el flujo de Fanno y el otro sobre el flujo
Rayleigh, pero ambos tienen de principal geometría la de una tubería adiabática de diámetro constante. El
objetivo para analizar es la validación del algoritmo de volumen finito basado en presión cuasi
unidimensional, utilizando como referencia el estudio de Alak Bandyopadhyay y Alok Majumdar [5] que
implementando el programa de simulación de sistemas de fluidos generalizados (GFSSP), para predecir el
flujo compresible con fricción, transferencia de calor. Las predicciones numéricas se comparan con las
soluciones clásicas de flujo compresible anteriormente mencionadas, pues el flujo Fanno proporciona una
solución analítica de flujo comprimible en una tubería larga y delgada don el flujo subsónico que entra
puede obstruirse por la fricción. Por lo cual, el flujo Raleigh da la solución analítica de flujo compresible
sin fricción con transferencia de calor donde el flujo subsónico que entra puede ahogarse en el limite con
el calor añadido al volumen de control.

Actualmente uno de los requisitos principales de simulación de flujo compresible es modelar con suma
precisión la fuerza de inercia que a menudo no se le da el interés debido en muchos códigos de análisis de
flujo de red. Es por ello por lo que la NASA Marshall Space Flight Center desarrollo un código de análisis
de flujo de red que se basó en volumen finito de carácter general: Programa de simulación de sistemas de
fluido generalizado (GFSSP) que es utilizado en su mayoría para el diseño de sistemas de propulsión para
vehículos de lanzamiento y análisis de flujo secundario de turbobomba. Estos casos de estudios
presentados a continuación sobre el flujo Fanno y Rayleigh buscan examinar a profundidad las
predicciones numéricas del GFSSP, es decir verificarlas y validar su veracidad[5].

3. Materiales y Métodos
3.1 Problema 1
En este análisis se considera un tipo principal de geometría: una tubería de diámetro constante.
El efecto producido por la fricción y la trasferencia de calor sobre la presión, la temperatura y el número
de Mach, se estudian utilizando las dos geometrías anteriormente mencionadas. Por lo tanto, se divide en
dos casos distinto, que son:
 Casos sobre flujo Fanno y Rayleigh
Se analiza el flujo Fanno como un flujo con fricción en una tubería adiabática de área constante. Para
este caso se considera que el flujo es unidimensional y la presión, temperatura y el número de Mach
se evalúan de manera analítica y numérica, ver figura uno.

Figura 1. Diagrama esquemático de una tubería de área constante[5].

Consideración: la figura uno es utilizado también en el caso de estudio del Flujo Rayleigh, ya que la
geometría es la misma, una tubería de área constante, pero se diferencian en que este flujo se analiza
sin fricción.
Nótese que la tubería mostrada en la figura uno tiene un factor de fricción constante (f). La longitud
de la tubería se escoge a modo que el flujo se obstruya a la salida de la tubería según lo determinado
por medio de la solución analítica.

 Cálculos matemáticos de referencia

El flujo compresible unidimensional generalizado se puede escribir matemáticamente utilizando las
ecuaciones de conservación[5]. Dichas ecuaciones son aplicables para estudiar un conducto de área
constante y están en forma diferencial:
La ecuación uno presenta la conservación de masa:
dρ dA dV
+ + =0 (1 )
ρ A V
donde:
dρ
=¿ Diferencial de densidad
ρ

MECÁNICA DE FLUIDOS II
MCs. GUILLERMO LÓPEZ TENORIO
 dA
=¿ Diferencial de área
A
dV
=¿ Diferencial de volumen
V
En la ecuación dos se muestra la conservación de momento:
dp γ M 2 fdx γ M 2 dV
+ ∙ + =0 (2)
p 2 D V
V V
M= = (2 b)

√
C p
γ∙
ρ

donde:
V 4 R ṁ √ T
M= =
√
∙
C π γ g c p D2
(2 c )

M= Número de Mach
R=Constante del gas
D= Diámetro
ṁ= Caudal
T= Temperatura
γ =Coeficiente de dilatación adiabática

Se considera la utilización de la definición del número de Mach, con la temperatura de

estancamiento, las ecuaciones uno y dos también, para poder expresar como una sola ecuación
diferencial ordinaria de primer orden[4].

( γ −1 2
) γ M ∙ f + (1+ γ M ) d T −γ M ∙ 1 ∙ dA (3)
[ ]
M 1+ ∙M 2
dM 2 2 0 2
=
dx ( 1+ M )
2 D 2 T0 A dx
Con valor límite, M ( x =0 )=M 1 ( 3b )

dT0
en la ecuación tres se puede determinar a partir de la ecuación cuatro sobre la energía y se
dx
puede expresar como:
q ( πD ) dx=ṁ c p d T 0 (4)

Dado el valor límite de la ecuación 3(b) y la ecuación diferencial de primer orden (ecuación 3) en M
se resuelve para encontrar el número de Mach en cualquier ubicación x. Como esta ecuación no es
lineal en M, se procede a resolver por el método de Runge-Kutta de cuarto orden[4], [5].

A partir del número de Mach, la temperatura y la presión estáticas se pueden determinar por medio
de las ecuaciones 5 y 6:
γ −1 2
T 0 ( x ) 1+ ∙( M ( 0) )
T(x) 2
= (5)
T ( 0) T ( 0 ) 1+ γ−1 2
0 ∙( M ( 0) )
2

 Modelo Numérico
√
p (x) A ( 0 ) M ( 0 ) T ( x )
=
p(0) A ( x ) M ( x ) T ( 0 )
(6)

Al basarnos en el artículo sobre Modelado de flujo compresible con fricción y transferencia de calor
utilizando el programa de simulación de sistemas de fluidos generalizados (GFSSP), se observa que
el programa utilizado emplea una formulación de volumen finito sobre las ecuaciones de
conservación de masa, momento y energía en una red que consta de nodos y ramas. Las ecuaciones
se dan en términos de entropía con generación de entropía debido a la disipación viscosa, en donde la
ecuación de conservación de masa y energía se resuelven para entropías en los nodos y para
presiones, mientras que la ecuación de conservación de momento se resuelve en las ramas para
calcular el caudal. La fricción en la tubería fue modelada por el factor de fricción de Darcy.

La caída de presión en una tubería se expresa como la ecuación 7:

∆ p=k f m2 (7)
donde:
8 fL
k f= (7 b)
ρ u π 2 D 5 gc

MECÁNICA DE FLUIDOS II
MCs. GUILLERMO LÓPEZ TENORIO
 Para la simulación numérica utilizando GFSSP, toda la tubería se divide en varios sectores, y cada
uno de estos sectores está representado por una tubería de diámetro constante. Los diámetros de los
sectores que están adyacentes pueden cambiar dependiendo el caso.
Se considera como fluido el nitrógeno y las soluciones numéricas dadas por GFSSP se observan en
los resultados y discusión, además se comparan con la solución analítica[5].

Se hace la consideración de que para la figura uno, en la simulación numérica se le hace una división
de esta en un número finito de tuberías de longitud no uniforme, como se aprecia en la figura dos. Se
coloca un nodo al final de ambos sectores de la tubería, incluyendo dos nodos limite. Se observan un
total de 21 nodos con nodos agrupados cerca de la entrada y salida en la figura dos, lo que es
suficiente para que tubería de área constante para tanto el caso del flujo Fanno como el de Rayleigh
obtenga una solución independiente de la red[5].

Figura 2. Distribución de nodos no uniforme para las tuberías de área constante[5].

 Caso 1: Flujo Fanno
´Para el flujo Fanno con fricción, pero sin cambio de área ni transferencia de calor se considera que
dT0
la temperatura de estancamiento es constante y =0.
dx
Se considero una longitud de 3207 pulgadas que proviene de la solución analítica de este tipo de
¿
flujo donde se tiene la longitud critica de la tubería L la cual es determinada mediante la ecuación
8:

(
(1+ γ) M 2
)
¿ 2
f L 1−M 1+γ
= + ln (8)
D γM
2
2γ 2+(1+γ ) M
2

donde:
M= Número de Mach de entrada
f= Factor de fricción
D= Diámetro de tubería
L*= Longitud crítica
γ =Coeficiente de dilatación adiabática
 Caso 2: Flujo Rayleigh
La solución analítica para los dos casos se obtiene a partir de la solución numérica generalizada de la
ecuación tres y las soluciones también pueden ser verificadas con las tablas de flujo Fanno y
Rayleigh disponibles en los libros de texto estándar.

3.2 Problema 2
A continuación, se presenta un caso de estudio aplicando la dinámica computarizada de fluidos (CFD) en
donde se verá la influencia del número de Reynolds en las características de flujo y transferencia de calor.
Descripción del caso de estudio
Un cilindro circular infinitamente largo de diámetro D=0.01 m se introduce a un flujo cruzado uniforme
con velocidad U, modelado en dos dimensiones.
Temperatura de 20°C y se utilizará la presión atmosférica.

MECÁNICA DE FLUIDOS II
MCs. GUILLERMO LÓPEZ TENORIO
 Como ya se sabe, el número de Reynolds es adimensional y se utiliza para ampliar las propiedades del
flujo, a continuación, la ecuación siguiente define el número de Re:
ρDU
ℜ=
μ
Donde:
Re: número de Reynolds
ρ : densidad
D: diámetro
μ: viscosidad

Re<5 Régimen de flujo no

separado

5 a 15 Re ≤ 40 Un par de vórtices
en la estela.

40 ≤ Re<90 y 90≤ Re<150 Dos

regímenes, es laminar.

150≤ Re<300 Es el rango de

transición a turbulencia

3x105<Re<3.5x106 La capa límite

está pasando a turbulenta, la estela
es un poco más estrecha.

3.5x106≤ Re<∞ Aquí es turbulento

y la estela es más delgada.

4. Resultados y discusión
4.1 Problema 1

 Caso 1: Flujo Fanno

De la ecuación 6 para el caso de flujo Fanno se obtiene que la longitud critica de la tubería es
básicamente la longitud necesitada para que el flujo se ahogue en la salida (considerando M=1 en la
salida). Con el número de Mach se obtiene una entrada de 0.5 y un factor de fricción de 0.002,
además se tiene un diámetro de tubería de 6 pulgadas, por lo cual la longitud critica es de
aproximadamente 3207 pulgas, como se muestra en la figura uno. Dicha longitud se mantiene fija
para los casos de tuberías de área constante.
En la figura dos se muestra un gráfico de la relación p/p * con diferentes tipos de distribución de
nodos para la solución numérica en comparación con la solución analítica. Una distribución de nodos
no uniforme con un total de (20 volúmenes de control) es suficiente para obtener una solución
independiente de la red.

Figura 2. Distribución de presión para flujo Fanno con varias distribuciones de rejilla[5].

MECÁNICA DE FLUIDOS II
MCs. GUILLERMO LÓPEZ TENORIO
 La gráfica proviene del artículo de Alak Bandyopadhyay y Alok Majumdar [5]. Nótese que cada
distribución de nodos presenta un color, para el color negro se representa una cuadricula no uniforme
con 20 cuadriculas, el color rojo es teórico, el color azul representa cuadriculas uniformes con 40
nodos y el verde cuadriculas uniformes con 20 nodos.

El número de Mach es una cantidad derivada de la solución numérica y se calculó a partir de la

ft ∙ lbf 2
ecuación (2c) de la cual se tiene que R=5.19 , g =32.17 ft /s .
lbmR c
La figura tres muestra un gráfico del número de Mach a lo largo de la dirección axial, y nuevamente
concuerda con la solución analítica. La única diferencia por considerar es en la entrada, lo que ocurre
por el GFSSP, donde no se determina la tasa de flujo másico, sino que se especifica las condiciones
límite de presión y se calcula la tasa de flujo.

Figura 3. Gráfica del número de Mach a lo largo de la longitud de la tubería[5].

La figura cuatro muestra la distribución de temperatura correspondiente al flujo Fanno.

Figura 4. Distribución de temperatura para flujo Fanno con varias distribuciones de cuadricula[5].
Nótese que la distribución en la gráfica roja muestra la teórica, en la de color negro se muestra una
cuadricula no uniforme con 20 nodos, la de color verde presente cuadriculas uniformes con 20 nodos
y la de color azul muestra una cuadricula uniforme con 40 nodos.

 Caso 2: Flujo de Rayleigh

El flujo de Rayleigh sin fricción y transferencia de calor uniforme muestra la distribución de las
figuras 2, 3 y 4 de temperatura, presión y número de Mach para una tasa de entrada de calor de
2088 Btu /¿s en la misma geometría de la tubería que se utilizó para el caso del flujo Fanno. El
número de Mach de entrada se eligió como 0.46 y se ha calculado analíticamente que con ese
número de Mach de entrada y una entrada de calor de 2088 Btu /s hace que el flujo se ahogue en la
salida (es decir, número de Mach=1). Nuevamente se aprecia que los resultados numéricos
concuerdan bastante con la solución analítica presentada en el artículo de Alak Bandyopadhyay y
Alok Majumdar [5].

Figura 7. Distribución de presión para el flujo de Rayleigh[5].

MECÁNICA DE FLUIDOS II
MCs. GUILLERMO LÓPEZ TENORIO
 Figura 8. Número de Mach de la distribución de flujo Rayleigh[5].

Figura 9. Distribución de temperatura para flujo Rayleigh (Q= 2088 Btu/s)[5].

Hay que considerar que tanto el gráfico de temperatura y de distribución de presión presenta el
número de Mach igual a uno. Y si se aumenta la tasa de calor hará que el flujo sea supersónico y, el
articulo utilizado solo considera flujo de tipo subsónico [4], [5].

4.2 Problema 2

Punto de
estancamiento

Figura 10. Para Re=1, se puede ver el coeficiente de fricción en su punto más alto, la velocidad del
fluido se desacelera hasta cero.

MECÁNICA DE FLUIDOS II
MCs. GUILLERMO LÓPEZ TENORIO
 Figura 11. Para Re=100, el coeficiente de arrastre empieza a aumentar, lo que provoca un aumento en la
presión en dirección opuesta al movimiento del flujo.

Figura 12. Re=1 x106, Flujo laminar. El coeficiente de arrastre se reduce un poco y tiende a tener una
diferencia de presión mayor entre el punto de estancamiento y la parte posterior del cilindro.

Figura 13. Re=1 x106, Flujo turbulento. El coeficiente de fricción va aumentando a medida que el flujo
acelera alrededor de la superficie del cilindro, para este caso de turbulento, una vez el flujo alcanza su
velocidad máxima, este comienza a desacelerarse y se refleja en el coeficiente de fricción.

MECÁNICA DE FLUIDOS II
MCs. GUILLERMO LÓPEZ TENORIO
 Figura 14. Simulación en 2D de un flujo en un cilindro.

5. Conclusiones
 En el problema uno se presento dos casos de estudio, el primero sobre el flujo Fanno en un flujo con
fricción en una tubería adiabática de área constante y el otro un flujo sin fricción y transferencia de
calor uniforme en una tubería adiabática de área constante que es el flujo Rayleigh, de ello se
concluye que ambos flujos pueden mantener una relación termodinámicamente por lo cual se puede
considerar que ambos pueden llegar a establecer una condición bajo la cual pueden tener un choque
normal en el flujo compresible de un gas en una tubería , además una vez analizado el efecto de la
fricción, transferencia de calor y el área en el flujo subsónico se pudo verificar la precisión del
algoritmo de volumen finito, donde se basó en la presión utilizando el programa de simulación de
sistemas de fluidos generalizados como referencia.
 Las soluciones numéricas de presión, temperatura y el número de Mach se compararon con la
solución analítica en los dos casos y se concluyo que hay un buen acuerdo tanto en la solución
numérica como en la solución de referencia.
 El flujo alrededor de un cilindro circular ha sido objeto de estudio, porque es un fenómeno muy
complicado, a pesar de ello existen software en el que podemos simular estos casos.
 Se pudo notar que el número de Reynolds determina el tipo de flujo y a medida que este aumenta
cambia de laminar a transitorio y por último a turbulento.

6. Referencias
[1] E. Rivera Chávez, “Flujo Compresible Apuntes de Clase Teoría y problemas resueltos.” [Online]. Available:
http://erivera-2001.com/flujo-compresible.html

[2] “Capitulo-3 FLUJO COMPRESIBLE UNIDIMENSIONAL.”

[3] “UNIVERSIDAD CATÓLICA SANTO TORIBIO DE MOGROVEJO FACULTAD DE INGENIERÍA

ESCUELA DE INGENIERÍA MECÁNICA ELÉCTRICA.”

[4] S. Colle, “Thermodynamic bounds for existence of normal shock in compressible fluid flow in pipes,” An
Acad Bras Cienc, vol. 89, no. 2, pp. 1313–1337, 2017, doi: 10.1590/0001-3765201720160292.

[5] A. Bandyopadhyay and A. Majumdar, “Modeling of Compressible Flow with Friction and Heat Transfer
using the Generalized Fluid System Simulation Program (GFSSP).”

[6] ANSYS Fluent Inc., Fluent 15.0.7 User’s Guide, 2014

[7] Fundamentals of Fluid Mechanics, B. R. Munson, et al (Capítulo 9) ‐Boundary layer characteristics.

9.2.6 Effect 

[8] Frank White. Mecánica de Fluidos, McGraw-Hill, Sexta Edición, México, 2008

[9] E.M. Ziemniak, N.K. Mitra y M. Fiebig, "Numerical Study of Two-Dimensional Vortex
Structure in a Channel Behind a Circular Cylinder", Numerical Methods in Laminar and
Turbulent Flow, Vol. 6, parte 1, pp. 553-563, (1989).

MECÁNICA DE FLUIDOS II
MCs. GUILLERMO LÓPEZ TENORIO