Repositorio Institucional de Documentos

Repositorio de la Universidad de Zaragoza – Zaguan http://zaguan.unizar.es Proyecto Fin de Carrera Ingeniería Industrial Simulación del flujo en una etapa de compresor centrífugo Autor Javier Azcona Pérez Director D. Francisco Alcrudo Sánchez Departamento de Ciencia y Tecnología de Materiales y Fluidos Área de Mecánica de Fluidos Centro Politécnico Superior Zaragoza, Diciembre 2011 Simulación del flujo en una etapa de compresor centrífugo - 1 - Contenido 1.- Introducción ......................................................................................................................... - 2 - 2.- Planteamiento y objetivo ..................................................................................................... - 2 - 3.- Geometría del compresor .................................................................................................... - 3 - 4.- Descripción del flujo ideal en el rodete ............................................................................... - 6 - 5.- Ecuaciones .......................................................................................................................... - 11 - 6.- Método de resolución ........................................................................................................ - 13 - 7.- Generación de los álabes ................................................................................................... - 19 - 8.- Mallado .............................................................................................................................. - 24 - 9.- Ejemplo de cálculo ............................................................................................................. - 26 - 10.- Conclusiones .................................................................................................................... - 32 - 11.- Trabajo futuro .................................................................................................................. - 33 - 11.- Bibliografía ....................................................................................................................... - 34 - 12.- Anexos .............................................................................................................................. - 35 - Simulación del flujo en una etapa de compresor centrífugo - 2 - 1.- Introducción Las turbo máquinas radiales llevan siendo usadas más de un siglo, si bien, su primer uso fue la impulsión hidráulica para dar paso posteriormente a los ventiladores, “blowers” y compresores. En los primeros desarrollos de propulsores para aviones, el tipo de compresor usado era el radial pero al poco tiempo fue desplazado por los compresores axiales ya que funcionaban mejor con los rangos de flujo másico que manejan este tipo de propulsores. A causa de esto, el desarrollo del compresor radial se dejó a un lado. Sin embargo, la necesidad de fabricar helicópteros militares impulsados por una pequeña turbina de gas hizo que la investigación y desarrollo del compresor radial volviese a cobrar importancia. Las ventajas que un compresor radial nos da respecto a la configuración axial son las siguientes: • La variación de entalpía que es capaz de suministrar una etapa radial es varias veces superior a la que nos proporciona una etapa axial por lo que se economiza espacio. • Las eficiencias con las que trabajan las máquinas radiales de estas magnitudes son mejores que las de su homólogo axial siempre que el Reynolds del flujo sea lo suficientemente pequeño • Menor coste de fabricación del compresor radial al contar únicamente con una sola etapa. Todo esto hace que estas máquinas sean las idóneas para propulsión de helicópteros, pequeñas turbinas de gas para vehículos terrestres, impulsión de aire para un sistema de ventilación a gran escala o la turbo inyección en motores diesel entre otros. 2.- Planteamiento y objetivo El estudio de este trabajo se centra en el análisis de los flujos que atraviesan un compresor centrífugo de un conjunto de mini turbina de gas. El análisis del flujo se hará sobre un diseño del compresor radial en vez de sobre el axial ya que la configuración radial presenta varias ventajas sobre la axial cuando el tamaño que se requiere que alcance la máquina es pequeño. Se ha simulado únicamente la parte en la que el flujo ya va prácticamente radial al completo debido a que el análisis ha sido bidimensional. También se estima que esta zona es donde se proporciona la práctica totalidad del gradiente de presión ya que es la zona donde el efecto centrífugo se hace más presente y por tanto, es la zona representativa del transcurso del flujo a través del rodete. Simulación del flujo en una etapa de compresor centrífugo - 3 - 3.- Geometría del compresor El compresor toma el aire de la atmosfera ambiente en la que se encuentra a través de la toma de aire. Tras esta, se encuentra con el rodete el cual podemos decomponer en dos partes diferenciadas. La primera parte es el inductor o “inducer”,el cual es casi completamente axial y se encarga de adaptar y guiar el flujo hacia la parte radial donde se realiza el aporte mayoritario de energía. La segunda parte es la zona radial en la cual el fluido sufre la mayor parte de la compresión que se adquiere en el rodete debido a la fuerza centrífuga. Tras su paso por el rodete, el fluido atraviesa el difusor, el cual tiene como cometido tranformar toda la energía cinética posible que sale del rodete en energía en forma de presión. El difusor puede tanto estar hueco como montar álabes fijos/móviles o incluso ambos. Por último tenemos la voluta la cual se encarga de recoger el flujo que sale del difusor y llevarlo a la salida del compresor. Figura 1; Partes de un compresor Simulación del flujo en una etapa de compresor centrífugo - 4 - Este proyecto se centra en la simulación de una aproximación 2D radial del flujo que atraviesa el rodete de un compresor debido a que, de todos los fenomenos que producen la compresión en el rodete de un compresor radial, la fuerza centrífuga es, con mucha diferencia, la que nos da el salto de presión requerido en la salida. El estudio se ha centrado en la zona radial del rodete, la cual es la encargada de aportar el gradiente de presión en esta zona casi en su totalidad. Además, se ha tomado la profundidad del canal tanto en la entrada como en la salida como constante debido a las restricciones del estudio 2D de forma que:   =  (3.1) La Figura 2 representa el corte transversal de un compresor centrífugo en el que se ha indicado la situación y designación de las diferentes variables cinemáticas y geométricas del mismo. En la zona izquierda de la figura se puede observar la vista frontal del corte donde se representan las velocidades tanto del punto de entrada como del de salida. Asi mismo, en la zona derecha de la figura se puede ver el lateral del corte y la representación tanto de los radios de entrada y de salida (R 1 , R 2 ) como de la profundidad del canal de paso del rodete designados por b 1 y b 2 (Como se ha dicho anteriormente, b 1 y b 2 seran iguales aunque en la figura se representen como diferentes). Cabe destacar que el sentido de giro tomado es siempre anti horario. El flujo tanto en la entrada como en la salida queda definido por los triangulos de velocidades donde v será la velocidad absoluta, u la velocidad periférica y w la velocidad relativa. Además, β, ángulo formado entra la tangente al esqueleto del álabe y la tangente al radio interior o Figura 2; Corte transversal de un compresor Simulación del flujo en una etapa de compresor centrífugo - 5 - exterior del rodete (según sea β en la entrada o en la salida respectivamente), es dado por la geometría del álabe y el ángulo α orienta la velocidad absoluta respecto a la velocidad periférica. Se dice que el flujo entra al rodete sin prerotación cuando α 1 vale  , es decir, cuando el flujo entra perpendicular a la velocidad periférica. En este estudio se ha fijado esta condición a la entrada por lo que el triángulo de velocidades resultante sera el siguiente: Figura 3; Triángulo de velocidades en la entrada Sin embargo, del triángulo de velocidades a la salida lo único que se puede fijar es β 2 ya que este ángulo viene dado directamente por la geometría del álabe. Figura 4; Triángulo de velocidades en la salida v 1 u 1 w 1 β 1 α 1 v 2 u 2 w 2 β 2 α 2 Simulación del flujo en una etapa de compresor centrífugo - 6 - 4.- Descripción del flujo ideal en el rodete Para comenzar con la simulación de los diferentes casos, ha sido necesario crear un programa que, mediante la aplicación de la “Teoría 1D” y asumiendo que la compresión sea isentrópica, nos de unos valores aproximados de la solución real para así poder fijar unas condiciones de contorno con cierto criterio. Se fijará T atm (Temperatura atmosférica), P atm (Presión atmosférica), M 1 (Mach en el punto 1), b 1 , b 2 (Profundidad del canal), R 1 , R 2 (Radios interior y exterior del rodete respectivamente), Ω (Velocidad de giro del rodete), R g , γ (Constantes del gas) y β 2 (Angulo β en el punto 2). Supondremos que el flujo entra al rodete sin prerotación por lo que el triángulo de velocidades a la entrada será como se muestra en la Figura 3. Por otra parte, el triángulo de velocidades a la salida tendrá la forma que anteriormente se ha indicado en la Figura 4. El esquema del rodete así como la situación de las variables viene dado por la siguiente figura en la cual se sitúan la entrada y la salida de nuestra simulación con sus respectivos radios, marcadas por los numeros 1 (Entrada) y 2 (Salida), además de las condiciones a las cuales se toma el gas. Figura 5; Situación de las variables de estudio P atm T atm Simulación del flujo en una etapa de compresor centrífugo - 7 - A partir de la imposición de los triángulos de velocidades y los datos de entrada fijados previamente, se procede al cálculo previo de las condiciones de contorno según las siguientes ecuaciones. En primer lugar, calculamos las velocidades periféricas en el punto de entrada y en el de salida mediante la siguiente expresión:   =  (4.1)   =  (4.2) A continuación calculamos la densidad atmosférica a partir de los datos de entrada:   =       (4.3) Tomando las variables anteriores y asumiendo el valor del número de Mach a la entrada (M 1 ) se puede obtener la temperatura (  ), la presión (  ) y la densidad (  ) en la entrada de acuerdo a las siguientes expresiones:   =  (1+     )  (4.4)   =  (1+     ) ! ! " (4.5)   =  (1+     ) " ! " (4.6) Además, serán de utilidad posteriormente las siguientes variables: Velocidad del sonido en la entrada #  =$% &   (4.7) Velocidad absoluta en la entrada '  =  #  (4.8) Sección de entrada (  =2*    (4.9) Gasto másico +  =  (  '  (4.10) Angulo,- en la entrada -  =tan 1 " 2 " (4.11) Coeficientes # 3 y # 1 del aire # 3 = &  (4.12) Simulación del flujo en una etapa de compresor centrífugo - 8 - # 1 = 4 5  (4.13) Entalpía en la entrada ℎ  =# 3   − 1 "8  (4.14) Velocidad relativa en la entrada 9  =$'  +:  (4.15) Mach relativo en la entrada  ; = < " 4 " (4.16) Temperatura y presión de remanso a la entrada en el sistema relativo (rotante)  =; =  (1+    ;  ) (4.17)  =; =  (1+    ;  ) ! ! " (4.18) Aplicando la ecuación de Euler en teoría 1D y la aproximación de compresión ideal e isentrópica obtenemos las siguientes expresiones: ℎ  −ℎ  = 2 88 2 "8  + < "8 < 88  (4.19)  8 > 8! =  " > "! (4.20) De la constancia del gasto másico, la velocidad relativa en la salida queda de la siguiente forma: 9  = < 8 ?@AB 8 = C D  8 E 8 ?@AB 8  > 8 (4.21) Haciendo uso de 4.20 podemos relacionar la entalpía en la salida con la densidad en la misma de la siguiente manera: ℎ  =# 3   = 45 8 >8 = F5 FG F5 FG FG" >"!>8! >8 =  " >"!   (4.22) Sustituyendo 4.21 y 4.22 en 4.19 obtenemos la siguiente expresión con   como incógnita:   " > "!   −ℎ  = 2 88 2 "8  + < "8  −  H C D  8 E 8 ?@AB 8 I  > 88 (4.23) La ecuación anterior relaciona la densidad a la salida con el resto de las variables a la entrada y los parámetros geométricos y cinemáticos del rodete para cada valor del número de Mach a la Simulación del flujo en una etapa de compresor centrífugo - 15 - b. Término fuente: El lado derecho de la ecuación 6.2 da cuenta de las fuerzas másicas debidas a la rotación del sistema. La evaluación de la integral de volumen es realizada con respecto a un punto como se muestra a continuación:  j ‚,ƒ epXf j ‚,ƒ =p €,• ·f €,• =p~{ €,• ,| €,• , €,• •f €,• (6.7) donde se muestra que el término fuente, H, depende tanto de las variables independientes (de posición) como de las variables de flujo U. c. Integración temporal: Una vez discretizadas espacialmente las ecuaciones se obtiene lo que se conoce como sistema semidiscreto que no es más que un conjunto de ecuaciones diferenciales ordinarias que describen la evolución temporal de los valores promedio en las celdas de la malla. Estas ecuaciones deben ser integradas en el tiempo mediante un método apropiado. El código empleado utiliza una discretización temporal implícita de primer orden incondicionalmente estable en el tiempo en el caso lineal. La discretización integral da lugar a un sistema implícito de ecuaciones no lineales que es linealizado de forma aproximada antes de ser resuelto por métodos matriciales. La derivada temporal se aproxima de la siguiente forma: H Wk W I €,• •‘P ≈ k‚,ƒ ’“"k‚,ƒ ’ ∆ (6.8) Con lo que el método completo se expresa como sigue: k‚,ƒ ’“"k‚,ƒ ’ ∆ =•  j‚,ƒ ∑ (st ∗ h) <‡ ·X( <‡ +p €,• ˆ‡‰ – •‘P (6.9) El parámetro θ indica el nivel de tiempo en el que se han evaluado los términos espaciales, es decir, los flujos y el término fuente. Si θ es cero, el método se reduce a un Euler explícito, mientras que si es uno, se trata de un Euler implícito, siendo ambos de primer orden de precisión. Si θ toma como valor ½ , el método se corresponde a la regla del trapecio que teóricamente es de segundo orden. En problemas no lineales muestra una estabilidad muy precaria por lo que, aunque el código de cálculo permite escoger cualquier θ, en la práctica la elección θ = 1 es la más estable y por ende, la que mejor funciona. La ecuación 6.9 representa la expresión del método en su forma no lineal. Para obtener el conjunto de variables a tiempo n+1 en todos los puntos de la malla es necesario resolver un sistema no lineal de tantas ecuaciones como celdas tenga ésta. El código de Simulación del flujo en una etapa de compresor centrífugo - 16 - cálculo hace uso de una linealizacion previa para poder aplicar métodos de resolución de sistemas lineales. En particular, métodos de relajación. k‚,ƒ ’“"k‚,ƒ ’ ∆ •1−—·∆˜H W™ Wk I €,• • –=•  j‚,ƒ ∑ (st ∗ h) <‡ ·X( <‡ +p €,• ˆ‡‰ – • (6.10) d. Condiciones de contorno: Hay basicamente cuatro tipos de condiciones de contorno a imponer en este problema: condición de entrada, condición de pared, condición de periodicidad y condición de salida. d.1 Entrada: Las condiciones a la entrada vienen físicamente determinadas por las propiedades de remanso y el ángulo de la corriente en el sistema relativo. De acuerdo con el procedimiento indicado en la sección 4, el Mach absoluto a la entrada, M 1 , fija directamente el gasto másico, Q m . Esto combinado con la condición de no prerrotación y el giro del rodete proporciona las condiciones de remanso y el ángulo de la corriente en el sistema relativo. Merece la pena destacar que, para un régimen de giro dado, sólo existe un valor del gasto másico, y por tanto del M 1 , que corresponde a un flujo de entrada alineado con la dirección de los álabes. Este sería el punto de funcionamiento nominal. Cualquier otro valor del gasto conduce a flujo no alineado con los álabes a la entrada lo que conlleva desprendimiento, disipación y una pérdida del rendimiento. Desde el punto de vista numérico, se imponen la presión y temperatura de remanso a la entrada junto con el ángulo de incidencia. La variable restante se obtiene por medio de extrapolación desde el interior del dominio de integración. Este procedimiento, al no fijar el Mach a la entrada, no determina el gasto másico, pues ambos son productos de la solución obtenida. En realidad, las condiciones a la entrada terminan siendo fijadas a convergencia por la condición de salida. Lo dicho anteriormente se mantiene siempre que no se alcancen condiciones supersónicas en la entrada o en la salida. d.2 Pared: Corresponden a este tipo las condiciones en los contornos delimitados por los álabes. El código de cálculo permite imponer condiciones de deslizamiento o de no deslizamiento. Las primeras se dan para flujo ideal (viscosidad nula) y las segundas para flujo de un fluido real (con viscosidad). d.3 Periódicas: Se impone periodicidad en las lineas del mallado anteriores y posteriores al alabe que conforman los canales de entrada y salida. Simulación del flujo en una etapa de compresor centrífugo - 17 - d.4 Salida: En el arco de salida es necesario imponer una única condición de contorno si el flujo es subsónico de acuerdo con la teoría de las características. Tres variables deben ser por tanto extrapoladas desde el interior del mallado. Aunque son posibles distintas opciones,la presión estática en salida se ha mostrado como la opción más razonable ya que proporciona resultados estables en la mayoría de los casos y está directamente relacionada con los parametros de funcionamiento del compresor. Una alternativa es imponer el Mach relativo en salida, M 2r , que resulta más estable pero cuya interpretación práctica es más compleja. Así pues, la presión de salida fija a convergencia el número de Mach a la entrada, M 1 , lo que a su vez establece el gasto másico. Figura 7; Condiciones de contorno La Figura 7 nos muestra los tipos de condiciones de contorno que se usan y la aplicación de cada una de ellas a las diferentes paredes de nuestro dominio. Como podemos observar, la condición de pared, marcada en verde, se aplica únicamente a la zona correspondiente a los álabes ya que es la única zona que el flujo no debe atravesar. Las paredes naranjas son aquellas que conectarán con un dominio adyacente. Las paredes naranjas de la zona superior se hallan conectadas con las de la zona inferior de forma que el flujo que atraviese unas sea el Pared Periódicas Periódicas Salida Entrada Simulación del flujo en una etapa de compresor centrífugo - 18 - mismo que atraviesa las otras. En azul se encuentra designada la pared a la cual se le aplican las condiciones de entrada y en rojo tendremos aquella en la que se fijaran las condiciones de contorno en salida. Simulación del flujo en una etapa de compresor centrífugo - 19 - 7.- Generación de los álabes Una parte importante de este trabajo consiste en la construcción de una herramienta destinada a la generación de la malla sobre la que realizar los cálculos. Como paso previo a la generación de la malla es necesario la definición de la geometría de los alabes. Se ha optado por definirlos como arcos de circunferencia que intersectan a los círculos de entrada y salida a unos ángulos prefijados γ 1 y γ 2 . Estos ángulos son los complementarios de los correspondientes a los anteriormente designados como β 1 y β 2 . De esta forma se genera la geometría del álabe para unos determinados ángulos de incidencia y salida que son escogidos por el usuario. La Figura 8, a continuación, muestra la disposición geométrica utilizada para construir la línea del álabe. Figura 8; Designación de las variables del problema general La construcción de la geometría consta de tres circunferencias, dos de las cuales son concéntricas (corresponden en el dibujo a las designadas por los radios R 1 y R 2 ) y otra decalada cuya intersección con las dos anteriores forma el llamado “esqueleto” del álabe (arco verde en la Figura 8). La geometría de las circunferencias que determinan el radio menor y el radio mayor del rodete están completamente resueltas ya que los datos conocidos serán tanto el radio menor, R 1 , de la circunferencia C 1 como el radio mayor, R 2 , de la circunferencia C 2 , ambas centradas en el origen. Fijaremos el punto de corte de la circunferencia decalada, C 3 , con C 1 en el punto P 1 (R 1 ,0). Simulación del flujo en una etapa de compresor centrífugo - 20 - Además de estos datos, también fijaremos los ángulos que formará el álabe respecto a la dirección radial en los puntos de intersección entre circunferencias siendo γ 1 el ángulo en P 1 y γ 2 el ángulo formado en P 2 , punto de intersección entre C 2 y C 3 . Según el esquema planteado en la Figura 8, el desarrollo de la solución es el siguiente: Figura 9; Detalle % , =,š+T (7.1)   = ?@A( " ›) ?@Aœ (7.2) Despejando ω obtenemos: š=,sin  (sin(*−%  −T)  "  8 ) (7.3) Por lo que volviendo a 7.1 y sustituyendo allí 7.3 tendremos: %  =sin  (sin(*−%  −T)  "  8 ),,+T (7.4) De la anterior ecuación se despejará α. Puesto que la ecuación es trascendente ha sido necesario programar una rutina de cálculo iterativo para resolverla. Una vez obtenido α se puede hallar θ, ángulo con el cual situaremos P 2 , mediante la siguiente ecuación: Simulación del flujo en una etapa de compresor centrífugo - 21 - θ=%  +α−sin  (sin(*−%  −T)  "  8 ) (7.5) Tras haber hallado θ, P 2 quedara definido por: {  =  cos(−—) (7.6) P 2 (x 2 ,y 2 ) |  =  sin(−—) (7.7) Con P 1 , P 2 , γ 1 y γ 2 conocidos, determinamos el centro, (x c3 ,y c3 ), y radio, R 3 , de C 3. Para ello, hallamos la intersección de dos rectas que pasa por P 1 y P 2 y que sean perpendiculares a la tangente de C 3 en dichos puntos: Figura 10; Designación de las variables para el cálculo de centro y radio de C 3 |=Ÿ  {+h  (7.8) |=Ÿ  {+h  (7.9) Simulación del flujo en una etapa de compresor centrífugo - 22 - Definimos las pendientes de las rectas como sigue: Ÿ  =  ¡A( " ) (7.10) Ÿ  =  ¡A(P 8 ) (7.11) El término independiente es definido como: h  =|  + m " ¡A( ") (7.12) h  =|  + m 8 ¡A(P 8 ) (7.13) Una vez sustituidas las pendientes y los términos independientes en las ecuaciones 7.8 y 7.9 igualamos ambas para despejar { 4¢ : Ÿ  { 4¢ +h  =,Ÿ  { 4¢ +h  (7.14) { 4¢ = (o8o")H £8 ¤¥¦(U“!8)£" ¤¥¦!"I " ¤¥¦!"" ¤¥¦(U“!8) (7.15) Con { 4¢ hallado, resolvemos la coordenada | 4¢ ,y, a continuación, hallamos el radio correspondiente a la circunferencia de centro ({ 4¢ ,| 4¢ ) que pasa por el punto P 1 y P 2 cuyo arco definirá el esqueleto de nuestro álabe: | 4¢ = m F§ ¡A " +H|  − m " ¡A " I (7.16)  ¢ =V~({  −{ 4¢) +(|  −| 4¢) • (7.17) Simulación del flujo en una etapa de compresor centrífugo - 23 - Figura 11; Esqueleto de los álabes de un rodete El resultado final es una geometría similar a la mostrada en la Figura 11 totalmente determinada. El arco obtenido representa la línea de curvatura media del álabe. Si este es de espesor infinitesimal solo resta generar el mallado entre dos álabes consecutivos. Si se consideran alabes de espesor finito, la herramienta creada es capaz de añadir una distribución parabólica de espesores sobre el esqueleto del alabe según la especificación del usuario. Una vez construido el álabe de espesor finito se procede al mallado del espacio comprendido entre dos alabes consecutivos de forma análoga al caso de espesor infinitesimal. El proceso de generación de la malla se describe en la sección siguiente. Simulación del flujo en una etapa de compresor centrífugo - 24 - 8.- Mallado Para mallar el espacio entre dos álabes, se ha optado por mallas de espaciado progresivo que acumulan nodos en las zonas más críticas y espacian más los mismos en las zonas donde no se requiere tanta resolución. En concreto, esta malla hace un mallado de espaciado progresivo tanto en sentido angular como en sentido radial. Tanto en la entrada/salida como en las paredes adyacentes a los álabes, la resolución de la malla ha de ser muy buena por lo que este espaciado aumenta desde los extremos hasta el punto medio del recorrido y, una vez llega a este punto, el crecimiento se invierte para que la malla vuelva a tener un espaciado pequeño en su otro extremo. La expresión que se ha utilizado para el cálculo del factor de crecimiento entre dos celdas consecutivas, ε, ha sido la de la suma de n términos de una progresión geométrica: ¨=© ; ’  ; (8.1) donde: ª=1+« (8.2) El término ©  corresponde al primer término de la progresión. Este valor será un dato de entrada para el caso del mallado del canal, tanto en sentido angular como en sentido radial. Sin embargo, el término ©  para el mallado radial en las zonas adheridas tanto a la entrada como a la salida del canal, heredará el valor del espaciado de los dos primeros nodos en sentido radial del canal para la zona que añadimos a la entrada, y de los dos últimos nodos en sentido radial del canal para la zona que añadimos a la salida del mismo. Figura 12; γ1 = 20, γ2 = 41, n=171, m=151, álabes curvados hacia delante Simulación del flujo en una etapa de compresor centrífugo - 31 - Figura 23; Resolución del Mach en el canal de la simulación 2 Figura 24; Líneas de corriente en la simulación 2 La Figura 24 nos muestra las líneas de corriente que atraviesan el canal. Como se puede ver, a diferencia del caso anterior, la fuerza de coriolis ayuda a guiar el flujo por lo que este sigue en todo momento el perfil de los álabes de forma ordenada. A la vista de esto, esta configuración de los álabes ayuda que el flujo se mantenga ordenado por lo que es propensa a dar menos problemas de convergencia debido a la formación de vórtices en la salida del rodete. Simulación del flujo en una etapa de compresor centrífugo - 32 - 10.- Conclusiones Este trabajo ha consistido en el estudio y posterior elaboración de un método para la simulación del flujo bidimensional en una etapa de un compresor centrífugo. El estudio se ha centrado en la zona radial del rodete por ser ésta donde se realiza la mayor parte del intercambio energético y ser susceptible del análisis bidimensional propuesto. Esta aproximación tiene la limitación de asumir un canal de paso de ancho constante. Para plantear correctamente el cálculo en la zona de álabes ha sido preciso realizar una estimación de las características de la corriente de aproximación al mismo. Para ello se ha realizado un estudio unidimensional del flujo en el rodete asumiendo flujo ideal, lo que implica entropía constante. Aunque dicho planteamiento es sencillo conceptualmente, las ecuaciones involucradas precisan de solución iterativa. Por tanto se ha acometido la construcción de una aplicación informática que las resuelva. Las hipótesis de flujo unidimensional e isentrópico se alejan considerablemente de la realidad, y sus resultados no son comparables con los de la solución numérica bidimensional obtenida en este trabajo. Sin embargo la aproximación unidimensional isentrópica representa una ayuda considerable a la hora de estimar la operación del rodete previa a la simulación numérica y permite imponer las condiciones de contorno con cierta fiabilidad en lugar de tenerse que recurrir a un tedioso proceso de prueba y error. Una parte considerable de este trabajo se ha dedicado a la generación de mallas destinadas a la simulación numérica. Se ha desarrollado un programa de mallado específico para la geometría característica de rodetes radiales. La aplicación informática permite seleccionar a voluntad todos los parámetros geométricos del rodete, entre los que se deben resaltar los ángulos de entrada y salida así como el espesor del álabe. El mallado, de tipo estructurado, incorpora canales de entrada y salida al rodete a voluntad del usuario así como compresión y expansión de la malla en las zonas de interés (bordes de ataque y salida, zonas de pared). Finalmente se ha adaptado un programa de cálculo de flujo compresible existente en el Area de Mecánica de Fluidos a las características particulares del flujo en un rodete radial. Para ello ha sido necesario estudiar la estructura del código y diseñar y ensamblar nuevas rutinas de cálculo. En particular se ha implementado el cálculo en sistema rotatorio incorporando las fuerzas centrífuga y de Coriolis y se han modificado las condiciones de contorno. Ha sido necesario también tener en cuenta el trabajo realizado por la fuerza centrífuga en la ecuación de la energía total. Se han mostrado ejemplos de cálculo para las dos configuraciones posibles de los álabes de rodete: Curvados hacia delante y hacia detrás. En ambas se han obtenido resultados realistas y compatibles con las estimaciones hechas mediante la aplicación desarrollada para flujo isentrópico unidimensional. Simulación del flujo en una etapa de compresor centrífugo - 33 - 11.- Trabajo futuro Dada la complejidad del flujo en máquinas de compresión radial, el alcance de este proyecto ha tenido que ser necesariamente modesto, ciñéndose a la puesta a punto de la metodología y aplicaciones necesarias para su estudio y simulación. Quedan por realizar algunas tareas importantes para conseguir disponer de una herramienta potente y flexible, pero sobre todo robusta. En este sentido es necesario dedicar una gran cantidad de trabajo a buscar ejemplos de geometrías reales para utilizarlas como casos de estudio y proceder a una verificación o validación sistemática de los resultados de la simulación. No sólo para los puntos de operación nominal sino también fuera de punto y así obtener curvas características. Para alcanzar dicho objetivo sería necesario primero optimizar los parámetros y estrategias de integración numérica. Los pases realizados muestran un gran potencial de optimización que no se ha podido investigar por falta de tiempo. Con el fin de conseguir la mayor fidelidad posible entre la simulación y el flujo real es necesario dedicar un esfuerzo considerable a implementar y ajustar la modelización de la turbulencia. El código del que se ha partido tiene incorporados modelos clásicos de turbulencia de pared que no se han utilizado en el presente trabajo. No obstante sería necesario incorporar modelos más sofisticados y comprobar su desempeño en el sistema rotante. Finalmente, abordar el flujo tridimensional sería el objetivo final más ambicioso, pero ello obliga a desarrollar un código desde cero, tarea que no parece que se pueda encuadrar en el marco de un proyecto fin de carrera sino posiblemente en un trabajo de tesis doctoral. Dentro de las limitaciones de la metodología puesta a punto es posible extender de manera considerable su rango de aplicación. Una prolongación inmediata de este trabajo es el estudio del flujo en el difusor del compresor. La hipótesis de flujo bidimensional empleada se adapta mejor a la zona del difusor que a la del rodete. Otro campo interesante de aplicación inmediata serían las turbinan radiales centrípetas cuya sólo requeriría la inversión del sentido del flujo. Es de esperar que la simulación de una turbina sea más sencilla que la de un compresor debido a que el gradiente de presión es siempre favorable en la primera. Simulación del flujo en una etapa de compresor centrífugo - 34 - 11.- Bibliografía S.L. Dixon, Fluid Mechanics and Themodynamics of Turbomachinery, Third Edition, Pergamon Press 1978 D.G. Wilson, The Design of High Efficiency Turbomachinery and Gas Turbines, Mit Press 1984 Area de Mecánica de Fluidos, Universidad de Zaragoza, Apuntes de Máquinas Hidraulicas y Transporte y Distribucion de Fluidos H.Krain, Review of Centrifugal Compressor’s Application and Development, ASME Journal of Turbomachinery, Vol. 127, pp.25-34 M.Marconcini, F.Rubeccini, A.Arnone,S.Ibaraki, Numerical Investigation of a Transonic Centrifugal Compressor, ASME Journal of Turbomachinery, Vol. 130, pp.1-9 Simulación del flujo en una etapa de compresor centrífugo - 35 - 12.- Anexos Los anexos de este trabajo serán incluidos en formato digital dentro de un DVD que quedará adherido a la contraportada. Los anexos incluyen: • Ejecutable del mallador • Código fuente del mallador en C++ • Ejecutable del preprocesador • Código fuente del preprocesador en C++ • Librerías necesarias