INSTITUTO POLITÉCNICO NACIONAL

ESCUELA SUPERIOR EN INGENIERÍA MECÁNICA Y ELÉCTRICA

UNIDAD PROFESIONAL TICOMÁN
INGENIERÍA AERONAÚTICA

FACTIBILIDAD DEL USO DE SOFTWARE CFD PARA

EL ANÁLISIS DE RENDIMIENTO DE UN TURBO-FAN
CON ACUMULACIÓN DE HIELO

TESIS
PARA OBTENER EL GRADO DE INGENIERO AERONÁUTICO CON ESPECIALIDAD EN
DISEÑO Y CONSTRUCCION

PRESENTA
CINDY BARRIOS ROJAS

ASESOR:
DR. OLIVER MARCEL HUERTA CHÁVEZ
CO-ASESOR:
ING. ISAAC GARIBAY SANDOVAL

MEXICO CD. MX
 FACTIBILIDAD DEL USO DE FLUENT CFD PARA EL
ANÁLISIS DE RENDIMIENTO DE UN TURBO-FAN
CON ACUMULACIÓN DE HIELO

INDICE
OBJETIVO....................................................................................................................................................... 3
JUSTIFICACION .............................................................................................................................................. 3
ALCANCE. ...................................................................................................................................................... 3
METODOLOGIA. ............................................................................................................................................ 3
CAPITULO1. INTRODUCCION........................................................................................................................ 3
1.1 ANTECEDENTES ......................................................................................................................... 3
1.2 PÉRDIDA DE POTENCIA ............................................................................................................. 4
1.3 ACUMULACIÓN DE HIELO ......................................................................................................... 5
CAPITULO 2. MODELOS MULTIFASICOS. .................................................................................................. 5
2.1 FLUJO MULTIFASE EN LA GOTA DE AGUA................................................................................. 5
2.1 MODELO TÉRMICO.................................................................................................................... 6
2.1 MODELO DE TURBULENCIA K-EPSILON..................................................................................... 7
CAPITULO 3. DOMINIO COMPUTACIONAL Y GEOMETRIA ............................................................................ 7
3.1 GEOMETRIA............................................................................................................................... 7
CAPITULO 4. ANÁLISIS DE RESULTADOS ....................................................................................................... 8
4.1 CASO 1: FAN ............................................................................... ¡Error! Marcador no definido.
4.2 CASO 2: STATIONARY WALL ....................................................... ¡Error! Marcador no definido.
4.3 CASO 3: MOVING WALL ............................................................. ¡Error! Marcador no definido.
BIBLIOGRAFIA................................................................................................................................................ 8
OBJETIVO.
Realizar el estudio de Factibilidad del uso de herramientas numéricas CFD para analizar la
perdida de rendimiento del fan del motor rotativo Trent 900 de RR por acumulación de hielo en
regiones clave.

JUSTIFICACIÓN
Es bien sabido que bajo condiciones de clima adversos, el rendimiento de los motores de
turbina de gas se deteriora y en casos extremos podría llegar a causar serios incidentes. Debido
a que el hielo se forma y se desprende continuamente en los alabes del fan, el rendimiento del
motor se ve afectado por la ingesta de este. Asimismo la modificación de la geometría en el
labio altera la dinámica del flujo teniendo repercusiones en el ciclo termodinámico.

Estos efectos están bien documentados y los procedimientos de diseño y certificación han sido
desarrollados para su identificación y erradicación.

ALCANCE.
Realizar el estudio de factibilidad del uso de software CFD para analizar el fenómeno de
acumulación de hielo en la parte frontal del motor, el cual suele presentarse bajo ciertas
condiciones climatológicas, lo que ayudara a saber o entender con más detalle las variables
principales que provocan la formación de hielo en el fan y en el labio de la carcasa del motor.

METODOLOGÍA.

CAPITULO 1. INTRODUCCIÓN
1.1 ANTECEDENTES
Desde 1990 han existido casos de motores Turbofan que presentan perdida de potencia en
altitudes superiores a los 22000 pies (6700 m), el límite más alto conocido en la industria por la
existencia de agua líquida súper-enfriada. A dichas altitudes pueden estar presentes partículas
en forma de cristales de hielo o incluso acumulaciones de cristales que han colisionado con
gotas de agua súper-enfriada para formar partículas más densas y esféricas tales como granizo.
El tamaño de estas partículas de hielo puede abarcar un amplio rango de tamaño, desde
micrómetros hasta centímetros.

Los principales efectos aerotermodinámicos son debido a los cambios de fase naturales en el
flujo de trabajo, el calor y la masa se transfieren entre las diferentes fases, provocando
modificaciones en las propiedades térmicas y de transporte.
 No es posible asegurar que el agua que entra en el motor, siempre será convertida en vapor al
momento que el fluido deja el compresor de alta presión o incluso la cámara de combustión.
Por lo tanto el estado del fluido y la distribución de la sección transversal del agua son de
preocupación a lo largo de la trayectoria del flujo en el motor

El rango de alcance de dichos efectos, dependen del diseño del motor, la cantidad del agua
ingerida y la respuesta del sistema de control.

1.2 PÉRDIDA DE POTENCIA

Por otro lado, la pérdida de potencia es definida como la inestabilidad en el motor tal como la
agitación, paro o extinción de llama dando como resultado una condición de reposo. A grandes
alturas se cree que estos cristales de hielo han causado la mayoría de estos eventos.

Existen distintas formas en la que se presenta la pérdida de potencia, dependiendo de las

características particulares de diseño en cada motor, algunas de estas se resumen en la
siguiente tabla [2].

TIPO DE PERDIDA DESCRIPCIÓN EFECTO

DE POTENCIA
Agitación /Paro El hielo se desprende hacia el compresor, Pérdida de empuje y
provocando una agitación en el motor, éste se temperaturas altas de los
para y la velocidad empieza a decaer reduciendo gases de escape
el flujo del aire, mientras la cámara de
combustión permanece encendida

Apague de Flama El hielo se desprende hacia la cámara de Pérdida de empuje y

combustión apagando la flama. descenso de todos los
parámetros

Daño al motor Los alabes del motor sufren daños al ser El daño producido en los
impactados por el hielo alabes, puede provocar
vibración en el motor o
paro del mismo

Tabla 1. Tipos de Pérdida de potencia que se presenta en el motor

 1.3 ACUMULACIÓN DE HIELO
Para poder entender la acumulación de hielo, primero se deben conocer algunos de los
principios que rigen este fenómeno. Lo primero que se debe considerar, es que existen
diferentes tipos de hielo que se pueden formar, un cristal de hielo con un tamaño típicamente
menor a 0.1mm es considerado como sólido, lo que lo difiere de la gota, ya que ésta se sigue
considerando como líquido [4].
La formación de escarcha y la de cristales, es el proceso que representa el crecimiento del hielo
por medio de la acumulación de gotas de agua súper-enfriada en superficies sólidas.
 La escarcha es una capa lisa de hielo que se extiendo a lo largo de la cuerda del
álabe, causada por la congelación de gotas de agua súper-enfriada al momento de
impactarse con la superficie, su apariencia opaca es debido al aire que queda
atrapado en las gotas de agua.
 El acristalamiento es causada por las gotas que fluyen en la superficie del álabe y se
congelan en una zona diferente de donde impactaron. Su apariencia se caracteriza
por ser irregular y transparente, además de la formación de dos picos (cuernos), por
la congelación de las gotas al ir y venir del punto de impacto.
 Y el mixto se caracteriza por la presencia de ambos, escarcha y cristales.

CAPITULO 2. MODELOS MULTIFASICOS.

2.1 FLUJO MULTIFASE EN LA GOTA DE AGUA.
La producción de cristales de hielo se lleva a cabo por la generación de calor latente, es decir,
la energía producida por la molécula de agua, mientras los átomos de la molécula se organizan
para formar una estructura cristalina.
La energía requerida durante este proceso es el calor latente de fusión, que tiene un valor de
333.7 kJ/kg.
El primer paso para realizar el análisis de la formación de hielo es, determinar la ubicación
donde impactan las gotas de agua súper-enfriada, para esto los parámetros de choque son
definidos de la siguiente manera: h- altura frontal de la corriente libre proyectada o área de la
estructura aerodinámica donde impactan las gotas, Y - los valores del límite superior 𝑌𝑈 e
inferior 𝑌𝐿 del impacto de la gota, 𝑆- la distancia de la envolvente, con el punto de
estancamiento del flujo del alabe, tomando como punto de referencia 0, y las distancias
superior 𝑆𝑈 e inferior 𝑆𝐿 de la envolvente.
Existen dos maneras diferentes de implementar el análisis de impacto de las gotas: por el
método numérico multifase Lagrangiano, y por el método Euleriano. El método Lagrangiano
consiste en calcular la posición de la trayectoria de cada gota desde el punto de origen hasta el
punto de impacto. Mientras que el Euleriano consiste en usar volúmenes de control para
determinar la posición de la gota. Ambos tienen sus puntos fuertes y débiles, el Lagrangiano es
un método más rápido, pero el Euleriano tiene mejor precisión y es fácilmente aplicable a los
códigos de la dinámica de fluidos computacional viscosos.
Cabe mencionar que los códigos computacionales para calcular la trayectoria de la gota dentro
de la turbo-maquinaria del compresor aún siguen en desarrollo.

2.1 MODELO TÉRMICO.

El modelo Térmico más convencional para la formación de hielo fue desarrollado por
Messinger, adicional a este modelo, otros dos modelos térmicos han sido recientemente
propuestos, uno basado en el Problema de Stefan y otro basado en los Medios Porosos de
Transferencia de Calor Multifase. Estos dos últimos son igualmente validos como el modelo
Messinger, pero no tan ampliamente utilizados en la comunidad de la formación de hielo
durante el vuelo.

A continuación se presenta el Modelo Messinger más detallado

𝑚̇ℎ𝑖𝑒𝑙𝑜 = 𝑚̇𝑖𝑚𝑝 + 𝑚̇𝑖𝑛𝑡 + 𝑚̇𝑒𝑥𝑡 + 𝑚̇𝑒𝑣𝑝

𝑚̇𝑖𝑚𝑝 = 𝑣∞ ∗ 𝐶𝐴𝐿 ∗ 𝛽

Donde como se había mencionado anteriormente, CAL es la cantidad de agua líquida por
unidad de volumen (g/m3) y 𝛽 es la eficiencia de la acumulación de partículas.

Por lo tanto, el flujo neto de masa es la suma de las gotas impactadas y el flujo reversible
dentro del volumen de control:

𝑚̇𝑛𝑒𝑡−𝑖𝑛 = 𝑚̇𝑖𝑚𝑝 + 𝑚̇𝑖𝑛

𝑚̇ℎ𝑖𝑒𝑙𝑜
𝑛=
𝑚̇𝑛𝑒𝑡−𝑖𝑛

La fracción de congelación (𝑛) es definida por la cantidad de hielo formado por la entrada de
masa neta.
Dado que las gotas se congelan de forma paralela a la trayectoria, este tipo de acumulación de
hielo tiende a ser muy aerodinámico, aunque su efecto es perjudicial debido a la modificación
de la geometría original del álabe.
La fracción de congelación es menor a 1.0, habrá evaporización y liquido reversible, en cuyo
caso 𝑚̇𝑒𝑥𝑡 + 𝑚̇𝑒𝑣𝑝 > 0.0. Este tipo de hielo es caracterizado por el predominio de agua líquida
durante el proceso de congelación generando hielo visualmente traslucido, similar a los cubos
de hielo. Debido a que el hielo glacial crece a lo largo de la zona que tenga mayor transferencia
de calor.

𝑞̇ 𝑓𝑟𝑒𝑒𝑧𝑒 = 𝑛 ∗ 𝑚̇𝑛𝑒𝑡−𝑖𝑛 ∗ 𝐿𝑓
 2.1 MODELO DE TURBULENCIA K-EPSILON
El denominado modelo k-épsilon estándar es el modelo más simple de los modelos de
turbulencia de dos ecuaciones de transporte. Con el objetivo de potenciar algunas deficiencias
de este modelo, a partir de él se han derivado dos modelos de turbulencia, estos son, el modelo
k−épsilon RNG y el modelo k-épsilon realizable [18].En éste trabajo nos enfocaremos al análisis
de éste último.

Modelo k-épsilon realizable.

En un modelo desarrollado recientemente, en el cual incluye una formulación alterna para la

viscosidad turbulenta y una nueva ecuación de transporte para el ratio de disipación, ε. El
beneficio del modelo realizable es que tiene una resolución superior para flujos que incluyen
rotación, capas limite bajo gradientes de presión fuertes y adversos, separación y recirculación.

Además de que presenta avances, cuando las características del flujo incluyen fuertes
curvaturas, vórtices y rotaciones, éste modelo ha demostrado que da mejores resultados y
comportamiento frente a otras versiones de modelos k-épsilon en el caso de flujos separados y
flujos secundarios complejos.

CAPITULO 3. DOMINIO COMPUTACIONAL Y GEOMETRIA

3.1 GEOMETRIA
Para este análisis se realizó el diseño de un motor Turbofan con base al modelo del motor Trent
900 de la compañía Rolls Royce, dicha geometría es un modelo 3D creado en Solidworks.

Uno de los objetivos de este trabajo es analizar la formación de hielo en la parte interna del
labio de la carcasa del motor; es por eso que al realizar este modelo se utilizó un perfil NACA
0015, el cual es un perfil simétrico, lo que ayudara a que se genere dicha acumulación de hielo
debido a que por su geometría tiende a generarse una zona de succión tanto en el extradós
como en intradós, siendo ésta última donde el hielo se depositara.

METODOLOGIA
Modelo  Ecuación de la Energía - on
 Solidification & Melting - on
 K-epsilon realizable
Condiciones de Frontera  Moving Wall
 Stationary wall
Inicialización de Solución  Inicialización standard

Calculo  Time step 1e-3

 Numero de time step
Tabla 2. Metodología para realizar los análisis
 CAPITULO 4. ANÁLISIS DE RESULTADOS

CAPITULO 4. ANÁLISIS DE RESULTADOS

Para la realización de este estudio se propusieron tres casos, los cuales son presentados en este
capítulo.

4.1 CASO 1: FAN CON FRONTERA MOVIL

El primero caso que se realizó es con la misma geometría del Turbofan presentada anteriormente, en
cuanto a las especificaciones de la malla son las que se encuentran en la Tabla 2, la variación que se
presenta es en los parámetros de la simulación, ya que para éste análisis solamente se produjo el paso
de flujo (aire) en el Turbofan para poder visualizar el paso del fluido, a continuación se muestran algunas
Figuras que nos permiten observar dicha interacción.

Figura 2. Contornos de velocidad-vista frontal Figura 1. Contornos de velocidad-vista trasera

Figura 4. Contornos de presión total-vista trasera Figura 3. Contornos de presión total-vista frontal
Con estas imágenes se puede corroborar que con este análisis y sus casos posteriores, regidos bajo el
mismo criterio de malla, conducen a resultados confiables, debido a que la velocidad máxima por el
análisis es de 403 m/s equivalente a las 2610 RPM.

4.2 CASO 2: FRONTERA FIJA

Para el Segundo caso se aplicó la condición de frontera fija trabajando con la misma metodología de la
Tabla 2.

Figura 5. Contornos de temperatura total en el Fan

Figura 6. Contornos de temperatura total en el labio

Figura 7. Contornos de temperatura total

Al aplicar la condición de frontera fija, se puede observar que tanto en el fan como en la carcasa interna
(labio) se genera un descenso de temperatura en el escalón debido a la zona de estancamiento propia
del mismo. El flujo bifásico entrando a bajas temperaturas no es comprimido lo suficiente por el fan
como para aumentar su temperatura y esta zona de obstrucción propicia que dicho flujo genere una
zona de baja presión cerca del escalón, dando lugar un estado de perpetuo frio en esta región.

4.3 CASO 3: FRONTERA MOVIL

En el caso tres se hizo la variación al aplicar la condición de frontera móvil trabajando con 2610 RPM
siguiendo la metodología de la Tabla 2.
Figura 9. Contornos de temperatura total en el labio Figura 8. Contornos de temperatura total en el Fan

Figura 10. Contornos de temperatura total en el cuerpo.

Con esta condición se puede observar que se sigue generando el mismo efecto en las zonas claves
donde se pretendía visualizar dicho fenómeno, teniendo una ligera variación en cuanto a la temperatura
debido a la velocidad angular que se declaró.

BIBLIOGRAFIA.
[1] Mason G. Jean, Strapp J. Walter, Chow P.,”The Ice Particle Threat to Engines in Flight”,
AIAA-2006-206, Reno Nevada, pp 6-8, Enero 2006.

[2] Mason J., “Engine Power Loss in Ice Crystal Conditions”, Boeing Aero-magazine Engine
Performance and Operability Propulsion Systems Division Artículo.

[3] Nikolaidis T. y Pilidis P., “The effect of water ingestion on an axial flow compressor
performance”, Institution of Mechanical Engineers, Part G, Diciembre 2012.

[4] Hutchings R., ‘’Effects of Supercooled Water Ingestion on Engine Performance”, University
of Tennessee, Knoxville, pp 1-4, August 2011.
[5] AGARD-AR-344, “Ice Accretion Simulation,” Advisory Group of Aerospace Research and
Development (AGARD) por sus siglas en ingles, NATO, 1997.

[6] Denti E. y Giuseppe Bindi F., “Flight in Icing Conditions”, por parte de la DGAC de Francia,
pp 7-9, 1997.

[7] Rios Pabón M.A., “Ice Crystal Ingestion by Turbofans”, Drexel University, Marzo 2012

[8] Libbrecht, K.G., “The physics of snow crystals”, Reports on Progress in Physics, pp. 855,
2005.

[9] Fukusako, S., “Thermophysical properties of ice, snow, and sea ice”. International Journal of
Thermophysics, pp. 353-372, 1990.

[10] Özgen, S., E. Tarhan, and M. Canıbek, “3-D In-flight icing simulations and use of parallel
computing in lagrangian droplet trajectory calculations”. 2011.

[11] Tong, X., E. Luke, ”Robust and Accurate Eulerian Multiphase Simulations of Icing Collection
Efficiency Using Singularity Diffusion Model”, Aplicación de Ingeniera de Dinámica de Fluidos
Computacionales, pp. 483-495, 2010.

[12] Crowe, C.Tl, “Multiphase flows with droplets and particles”, imprenta CRC, 2011.

[13] Ruff, G.; Berkowitz, B., “Users manual for the NASA Lewis Ice Accretion Prediction
Code(LEWICE)(Final Report)”. 1990.

[14] Saeed, F., “State-of-the-Art Aircraft Icing and Anti-Icing Simulation”. Revista ARA (American
Romanian Academy), pp. 2000-2002.

[15] Myers, T.G., “Extension to the Messinger model for aircraft icing”, revista AIAA, pp. 211-
218, 2001

[16] Fedorov, A., V. Fomin, y A. Shul’gin, “Mathematical description of motion and adherence of
droplets to the flow boundaries. Journal of Engineering Physics and Thermophysics”, pp. 1123-
1127, 2008.

[17] Mishima, O. and H.E. Stanley, “The relationship between liquid, supercooled and glassy
wáter”. Universidad Politécnica de Cartagena, pp. 329-335, 1998.

[18] Crevillén D., “Simulación numérica del flujo a travès del tubo diffusor de una turbine de
reacciòn”. Nature, pp. 26-29, 2012.

[19] Fluent Inc., Fluent 6.3 user’s guide, 2006.

[20] Tasdelen V., Hertel J., Fiola R.,” A numerical Investigation on ingested water separation in
modern turbofan engines”, Institute of Aircraft Propulsion Systmes , University of Stuttgart,
pp1-3.

[21] Obermuller M., Schmidt and D.Peitsch K.J., “Some aspects on wet compression-physical
effects and modeling strategies used inn engine performance tool”, Deutscher Luft-und
Raumfahrtkongress, DLRK-2012-281210, Berlin 2012.

[22] White A. and Meacock A.J., “An evaluation of the effects of water injection on Compressor
Performance”, Proceedings of GT2010, ASME Turbo Expo, 2010.

[23] Pilch M. y Erdman C.A, “Use of breakup time data and velocity history data to predict the
maximum size of stable fragments for acceleration-induced breakup of a liquid drop”,
International Journal of Multiphase Flow, 1987.