“Año del Fortalecimiento de la Soberanía Nacional”

Dinámica de Fluidos Computacional (CFD)

Avance de Trabajo Final 1

Grupo: 6

Alumno: Cieza Mostacero, Kenji Gianfranco – U21208400

Docente: Taipe Ramos, Raphael Cristian

Lima – Perú

2022
 Tabla de Contenidos
1. Introducción:.....................................................................................................................3
1.1. Objetivos....................................................................................................................4
2. Desarrollo del Trabajo:......................................................................................................5
2.1. Procedimiento para la Tobera 1................................................................................5
2.2. Procedimiento para la Tobera 2..............................................................................17
2.3. Procedimiento para la Tobera 3..............................................................................29
3. Resultados:.....................................................................................................................42
3.1. Resultados de la Tobera 1.......................................................................................42
3.2. Resultados de la Tobera 2.......................................................................................45
3.3. Resultados de la Tobera 3.......................................................................................47
4. Análisis de los Resultados:.............................................................................................49
5. Conclusiones:.................................................................................................................51
1. Introducción:

En este trabajo se hará uso de la herramienta Ansys Fluent para realizar el análisis de las
siguientes toberas convergentes - divergentes:

Figura 1: Tobera 1 para analizar en 3D

Figura 2: Tobera 2 para analizar en 2D

 Figura 3: Tobera 3 para analizar en 2D

1.1. Objetivos

 Aplicar la herramienta Ansys Fluent para analizar diferentes tipos de toberas bajo
diferentes condiciones iniciales.
 Obtener algunos gráficos de los diversos parámetros que se pueden simular con
esta herramienta.
 Analizar y discutir los resultados y gráficos obtenidos de la simulación con este
software.

2. Desarrollo del Trabajo:

 2.1. Procedimiento para la Tobera 1
Abrimos Workbench, hacemos doble click en Fluid Flow (Fluent), configuramos la geometría
en 3D, le ponemos nombre al proyecto y abrimos el DesignModeler de geometría:

Cambiamos las unidades a milimetros:

Creamos un nuevo sketch en el plano XY y con la herramienta polyline dibujamos la forma
de la mitad de la tobera 1:

En dimensiones ponemos las medidas correspondientes:

Usamos el comando Revolve y elegimos como geometría al sketch 1 y como eje al eje X y
luego hacemos click en Generate:

Finalmente cerramos la geometria:

Hacemos doble click en Mesh:

Creamos los nombres para el inlet, outlet y walls de la tobera

Cambiamos el tamaño del elemento del mallado a 3 mm, generamos la malla y cerramos:

Click derecho en mesh y actualizamos la malla, luego hacemos doble click en setup:
Elegimos doble precisión y 8 procesadores para los análisis (dependiendo de la maquina):

Cambiamos el modelo de viscosidad:

Creamos el material benzene-vapor (c6h6):

Cambiamos el material de las celdas al benzene-vapor:

Colocamos la presión en la entrada en 100 GPa:

Verificamos la presión a la salida:

Elegimos el método de solución:

Configuramos todos los residuos en 0.00001:

Luego vamos a inicialización, y hacemos una inicialización estándar:

Finalmente vamos a run calculation y calculamos la solución para 1000 iteraciones:

Esperamos que termine de calcular la solución:

Nos vamos a resultados > graficas > contornos y creamos el plano XY para ver las gráficas
de contornos:
Similarmente creamos graficas de vectores:

Por último, obtenemos los graficos XY de velocidad y presión:

 2.2. Procedimiento para la Tobera 2
Cargamos un nuevo proyecto de Ansys Fluent que nombraremos como Tobera 2,
configuramos la geometría en 2D y abrimos el DesignModeler

Cambiamos a milímetros, creamos un sketch en el plano XY y dibujamos la mitad de la

forma de la tobera 2:
Se colocan las dimensiones correctas de la Tobera 2

Creamos una superficie a partir de este sketch

Dibujamos 2 líneas para separar la superficie en 3 caras:

Usamos FaceSplit con las 2 lineas que creamos anteriormente:

Cambiamos el material a tipo fluido:

Entramos al mallado:
Colocamos los nombres al inlet, outlet, walls y al eje de simetria

Se insertan mallados de cara y de tamaño y se genera la malla

Actualizar la malla y abrir el Setup

Poner doble precisión y 8 procesadores (dependiendo de la maquina):

Verificar la calidad de la malla y cambiar los valores para la solución:

Activar la ecuación de la energía:

Cambiamos el modelo de viscosidad

Cambiamos la densidad y viscosidad del aire, por la densidad de un gas ideal y la ley
Sutherland de viscosidad.
Ponemos como dato inicial la presión a la entrada de 20 atmosfera

Cambiar el courant number (tiempo que una partícula permanece en una celda) a 0.5
Configurar los residuos en 0.000001

Inicializar la solución con una inicialización estándar:

Correr los cálculos para 20000 número de pasos

Una vez la solución converge podemos parar la simulación:

Para las simulaciones, hacemos las gráficas de contornos:

Para que se vea la tobera completa aplicamos el eje X de simetría:

 2.3. Procedimiento para la Tobera 3
Cargamos otro Fluid Flow (Fluent) para hacer la Tobera 3 y abrimos DesignModeler

Cambiamos las unidades a milímetros y dibujamos la forma de la tobera 3 en un sketch en

el plano XY:
Ponemos algunas restricciones para que el dibujo sea simétrico:

Colocamos las dimensiones correspondientes para la tobera 3

Generamos una superficie a partir del sketch:

Dibujamos 2 líneas que corten la tobera a 5 mm de la línea de la garganta:

Se usa Face Split para dividir la superficie en 4 caras y se cierra la geometría:

Entramos al mallado
Se crea el inlet, outlet y walls:

En el mallado se inserta un mallado de caras y 2 sizing

Se genera la malla y cerramos

Actualizar mallado y abrir setup:

Configuramos Fluent con doble precisión y 8 procesadores lógicos (depende de maquina):

Chequeamos la calidad de la malla y cambiamos el Solver a Density-Based

Activamos la ecuación de la energía:

Cambiamos la densidad y la viscosidad del aire:

Configuramos la presión en la entrada a 70 barn = 7MPa

Ponemos los valores de referencia:

Configuramos los residuos en 0.0000001:

Inicializamos la solución:
Calculamos la solución para 600 iteraciones:

Obtenemos las gráficas de contornos para 5 simulaciones:

Graficamos la presión y la velocidad en las paredes vs la posición y cerramos el fluent:

Actualizar la solución y abrir resultados:

Generamos un streamline de la velocidad:

Finalmente animamos esta gráfica:

3. Resultados:

3.1. Resultados de la Tobera 1

Figura 4: Presión Estática en la Tobera 1

Figura 5: Presión Dinámica en la Tobera 1

Figura 6: Magnitud de la Velocidad en la Tobera 1

 Figura 7: Densidad del gas que pasa por la Tobera 1

Figura 8: Turbulencia de la Energía Cinética en la Tobera 1

Figura 9: Vectores de la Presión Estática en la Tobera 1

Figura 10: Vectores de la Turbulencia de la Energía Cinética en la Tobera 1

Figura 11: Presión Estática (Pa) vs Posición (mm) en la Tobera 1

Figura 12: Magnitud de la Velocidad (m/s) vs Posición (mm) en la Tobera 1

3.2. Resultados de la Tobera 2

Figura 13: Residuos de la Simulación para la Tobera 2

Figura 14: Presión Estática en la Tobera 2

Figura 15: Temperatura Estática en la Tobera 2

 Figura 16: Número de Mach del fluido en la Tobera 2

Figura 17: Turbulencia Cinética en la Tobera 2

Figura 18: Densidad del fluido a lo largo de la Tobera 2

3.3. Resultados de la Tobera 3

Figura 19: Presión Total en la Tobera 3

Figura 20: Presión Estática en la Tobera 3

Figura 21: Número de Mach en la Tobera 3

 Figura 22: Densidad en la Tobera 3

Figura 23: Temperatura Estática en la Tobera 3

Figura 24: Turbulencia Cinética en la Tobera 3

 Figura 25: Presión Estática (Pa) vs Posición de las paredes (m)

Figura 26: Numero de Mach vs Posición de las paredes (m)

4. Análisis de los Resultados:

Tobera 1

 Se usa como gas el benceno-vapor (c6h6), no se usa la ecuación de la energía y se usa

el modelo de viscosidad el k-épsilon.
 Se tiene que la presión estática es alta en la entrada de la tobera, disminuye hasta
valores negativos altos en la garganta de la tobera y luego, pasada la garganta vuelve a
aumentar hasta valores cercanos a la presión de salida o atmosférica.
 La magnitud de la velocidad en cambio, empieza en valores bajos y al acercarse a la
garganta se acelera hasta su valor máximo para luego volver a decrecer
 paulatinamente, excepto en la zona afectada por la turbulencia en donde la disminución
es drástica.
 La presión dinámica al estar relacionada con la velocidad tiene el mismo
comportamiento que la magnitud de la velocidad.
 La turbulencia aparece con mas notoriedad en la zona baja de la parte final de la tobera
 La densidad y la temperatura no varían debido a que el análisis se realizo considerando
al gas como un fluido incompresible

Tobera 2

 Se usa como fluido al aire tratándolo como gas ideal y con su viscosidad dependiente de
la temperatura (Sutherland). Además, se usa la ecuación de la energía y el modelo de
viscosidad k-épsilon. Para el solver se usó un tipo basado en la densidad y un tiempo
transitorio.
 La presión estática, temperatura estática y la densidad tienen un comportamiento muy
similar, empiezan en un valor muy alto y al llegar a la garganta disminuye formando un
patrón característico desde ese punto hacia adelante. El comportamiento similar de
estos 3 parámetros se debe a la consideración del aire como gas ideal.
 Por el contrario, la velocidad como numero de mach, empieza en un valor pequeño y a
partir de la garganta este valor aumenta formando un patrón característico.
 Casi no hay turbulencia en este tipo de tobera, solo se aprecia turbulencia en un
pequeño tramo cerca de las paredes al final de la tobera.

Tobera 3

 Se usa como fluido al aire tratándolo como gas ideal y con su viscosidad dependiente de
la temperatura (Sutherland). Además, se usa la ecuación de la energía y el modelo de
viscosidad k-omega. Para el solver se usó un tipo basado en la densidad.
 Hay una pérdida de presión total en la tobera sobre todo en las zonas de mayor
turbulencia.
 Debido a la alta presión en la entrada se genera una gran turbulencia en esta zona y
también a lo largo de las paredes de la tobera.
 La presión estática empieza en un valor bajo en la entrada, probablemente por la
turbulencia existente. Como no hay una zona de garganta, se genera una zona de
mayor presión en el centro de la sección donde cambia la tobera de convergente a
divergente.
 Debido a que se considera al aire como gas ideal, la temperatura estática y la densidad
tienen un comportamiento similar al de la presión estática.
 La velocidad en cambio, comienza con un valor alto en la entrada y va disminuyendo
hasta llegar a la zona de cambio de tobera convergente a divergente y luego vuelve a
aumentar.

5. Conclusiones:

 Se concluye que las soluciones obtenidas son diferentes dependiendo de las

condiciones iniciales que le demos al problema.
 Las soluciones obtenidas también dependen del tipo de solver elegido, del uso o no de
la ecuación de la energía y del modelo de viscosidad elegido.
 Las zonas de turbulencia afectan mucho al gráfico de las presiones y velocidades.
 Las zonas de turbulencia son diferentes para las diferentes formas de las toberas y para
distintas condiciones iniciales.
 En general, cuando la tobera tiene una garganta pronunciada, en esa zona la presión
estática tiende a disminuir y la velocidad tiende a aumentar.