OPERACIONES DE TRANSFERENCIA DE
CANTIDAD DE MOVIMIENTO
__________________________________________________________________________
CÁLCULO DE PÉRDIDAS EN UN TUBO CON EXPANSIÓN SÚBITA EN FLUJO LAMINAR
Universidad Nacional de Colombia sede Medellín
__________________________________________________________________________
Resumen
Se denomina expansión súbita, cuando un fluido atraviesa de un conducto pequeño a otro
más amplio, su velocidad se reduce de forma repentina provocando turbulencias, las cuales,
ocasionan pérdidas de energía. En este trabajo se desea determinar dichas pérdidas, de
modo que, se modela el comportamiento de un fluido empleando dos métodos, el CFD
utilizando la herramienta Ansys Student Academic y el método de las 2K de Hooper, esto con
el fin de compararlos y comprobar si el valor de K es dado por los tamaños de los dos
conductos y la velocidad del flujo.
Datos
Como no se especifica el material del tubo se usó la guía de tubería ferrosas para obtener
diámetros internos y longitudes:
Tabla 1. Valores de diámetros y espesores para tubería
Diámetro 1: ½ 0.622 in -> 0.015798 m
Diámetro 2: 1 ½ 1.610 in -> 0.04094 m
Datos Simbología Valor Unidades
OPERACIONES DE TRANSFERENCIA DE
CANTIDAD DE MOVIMIENTO
__________________________________________________________________________
DN Sch 40 DN ½ y1½ in
Radio 1 R1 0.0078994 Metro (m)
Longitud 1 L1 0.1579 Metro (m)
Radio 2 R2 0.2047 Metro (m)
Longitud 2 L2 1.0224 Metros (m)
Número de Reynols Re 1100 adimensional
Velocidad flujo agua V 0.06997 m/s
Temperatura T 288.1 Kelvin (ºK)
Presión P 1 atm
Densidad del agua a 25ºC 𝜌 998,2 Kg/m^3
Viscosidad del agua a 25ºC μ 1,003x10-3 Kg/m*s
Tabla 2. Datos utilizados para llevar a cabo los análisis
1. Cálculo de K usando software ANSYS
En primer lugar, para simular la expansión súbita por una tubería con un fluido (agua a
288,16 °K) se requiere establecer condiciones de frontera. Una de esas condiciones iniciales
corresponde a la velocidad del fluido. La relación que permite hallar este parámetro es la que
describe la ecuación de Reynolds, puesto que su valor es conocido y está dado por:
Ec [1] al despejar V obtenemos que: Ec [2]
Donde Re = 1100; D1 =0.015798 m; D2 = 0.04094 m; 𝜌 = 998,2 Kg/m^3; 𝜇 = 1,003x10-3
Kg/m*s. Sustituyendo y reorganizando la Ecuación [2], obtenemos que la velocidad es:
V1 = Velocidad respecto al diámetro (D1)
V2 = Velocidad respecto al diámetro (D2)
V3 =Velocidad respecto al diámetro promedio.
((D1+D2)/2)
OPERACIONES DE TRANSFERENCIA DE
CANTIDAD DE MOVIMIENTO
__________________________________________________________________________
Para los análisis en ANSYS CFD se va a utilizar la velocidad del diámetro más pequeño, es de decir,
diámetro 1 (D1).
Figura 1. Perfil de velocidad de la tubería
Figura 2. Perfil de presión de la tubería
Figura 3. Gradiente de velocidad para 1000 iteraciones
OPERACIONES DE TRANSFERENCIA DE
CANTIDAD DE MOVIMIENTO
__________________________________________________________________________
Figura 4. Gradiente de presión para 1000 iteraciones
Figura 5. Presiones a la entrada y salida
A partir de los resultados obtenidos en la simulación, se determina que 𝑃𝑒 = 1.8191254 [𝑃𝑎]y
𝑃𝑠 = 0.36382509 [𝑃𝑎]. De la ecuación de Bernoulli, se sabe que las pérdidas están dadas por
la relación del cambio de presión en el sistema, la densidad del fluido y la gravedad, teniendo
así para calcular las pérdidas o caídas de presión la ecuación:
Ec [3]
Con este resultado y a partir de la siguiente ecuación se despeja K.
Ec [4]
OPERACIONES DE TRANSFERENCIA DE
CANTIDAD DE MOVIMIENTO
__________________________________________________________________________
2. Cálculo usando método de las dos k
El método de las dos K es una estrategia matemática usada para hallar una constante K que
engloba la relación entre f, L y D del accesorio estudiado. Para hallar las K, se deben
emplear tablas, de acuerdo con la forma del accesorio, su diámetro y el tipo de material,
estiman un valor experimental para sus propiedades geométricas generales. Esta constante
está dada por:
Bernoulli entre A1 y A2, donde Z1 = Z2
Balance de continuidad: (1) Fuerza en términos de área y presión:
(2)
Igualando (1) y (2), y reemplazando en términos de las pérdidas, obtenemos que:
Teorema de Carnot
Ec [5]
Pérdidas causadas por la expansión, con la ecuación analítica
dada para una expansión
súbita:
OPERACIONES DE TRANSFERENCIA DE
CANTIDAD DE MOVIMIENTO
__________________________________________________________________________
3. Cálculo del porcentaje de error
Análisis de resultados
Las pérdidas obtenidas por el método de las 2k son más altas que las calculadas por ANSYS
Student Academic, una posible causa es porque la simulación toma el tubo como ideal,
entonces las pérdidas son menores. Al haber discrepancia entre los dos métodos el
porcentaje de error es elevado.
La simulación muestra un cambio de velocidad en la tubería y el efecto de la relación de
áreas (A1/A2) en las pérdidas de energía que ocurren por ampliación brusca del diámetro del
tubo donde se puede notar que las pérdidas aumentan al decrecer dicha relación, debido a
que se incrementa la perturbación del flujo por el ensanchamiento drástico.
Para seleccionar la ecuación que permite estimar el valor del coeficiente de pérdidas K en los
dispositivos se consideró un modelo que tuviera la mayor simplificación posible, permitiendo
obtener el valor del coeficiente de pérdidas con la mayor facilidad y credibilidad, por ello, en
el método de las 2k, el valor de K depende exclusivamente de la geometría al ser una
expansión súbita.
Conclusiones
La revisión literaria demuestra que el método para tasar el valor del coeficiente K requiere la
utilización de tablas y gráficos, sin embargo, para la realización de este trabajo -expansión
OPERACIONES DE TRANSFERENCIA DE
CANTIDAD DE MOVIMIENTO
__________________________________________________________________________
súbita- solo se precisa el valor de la relación de áreas A1/A2 esta expresión es útil y factible
para obtener el coeficiente K obviando la manipulación de tablas y gráficos.
El programa ANSYS permite establecer parámetros para simular los cambios que afectan el
modelo estudiado, en el caso de la presión se observa que en un primer instante la presión
es mayor, sin embargo, a una mayor relación de expansión se obtiene una caída de presión,
y un aumento de perdidas por turbulencia debido a que el área de recirculación es mayor.
Finalmente, a futuro se espera que las ecuaciones [3] y [4], sean programadas en algoritmos
computacionales que impliquen el análisis y modelamiento como ANSYS para simplificar la
tarea en la obtención de las pérdidas de energía y el valor de K, al incluirlas en su Software
permite ahorrar tiempo en el cálculo de dichos valores .
Referencias
Pérdidas por fricción en accesorios.” n.d.
[Link]
N. Connor, "thermal engineering", [Link], 2020. [Online]. Available:
[Link]
local-definicion/
J.J. Villegas-León, A.A. López-Lambraño, A.A. López-Ramos, C. Fuentes & M. López-
Lambraño (2018) Expresiones analíticas del coeficiente de pérdida K para la ampliación
brusca ó gradual del diámetro, Ribagua, 5:2, 128-142.