Manual Pract

Área de Mecánica de Fluidos CURSO 2006-07
PRÁCTICAS DE
FLUIDODINÁMICA COMPUTACIONAL
PRÁCTICA Nº 1
SIMULACIÓN DEL FLUJO EN UNA TUBERÍA
1. DESCRIPCIÓN DEL PROBLEMA

2. GENERACIÓN DE LA GEOMETRÍA EN GAMBIT
3. GENERACIÓN DEL MALLADO EN GAMBIT
4. CONDICIONES DE CONTORNO Y EXPORTACIÓN DESDE GAMBIT
5. INICIALIZACIÓN DEL CASO EN FLUENT: CASO LAMINAR.
6. RESOLUCIÓN
7. ANÁLISIS DE RESULTADOS
8. REFINADO DEL MODELO
9. RESULTADOS FINALES
10. CASO TURBULENTO
PRÁCTICA No. 1 - FLUENT
SIMULACIÓN DEL FLUJO EN UNA TUBERÍA Área de Mecánica de Fluidos
1.- DESCRIPCIÓN DEL PROBLEMA
Considérese un fluido moviéndose a lo largo de una tubería de sección circular constante. El

diámetro de la tubería es D y la longitud de la misma L. La velocidad del fluido a la entrada es V.
Considérese así mismo que la velocidad es constante a lo largo de la sección de entrada. El
fluido descarga contra la presión atmosférica a la salida del tubo a una presión de 1 atm. Tómese un
fluido de densidad, ρ=1 kg/m3, y de viscosidad, μ=2x10-3 kg/(ms).
El número de Reynolds basado en el diámetro de la tubería se calcula como:
ρVD
Re =
μ
V
El ejercicio a resolver consiste en la realización de una simulación del comportamiento de este flujo
mediante el programa FLUENT. En concreto, han de mostrarse los siguientes resultados:
• Mostrar gráficamente la distribución de velocidad en la línea media.
• Obtener el coeficiente de fricción en las paredes del conducto.
• Mostrar gráficamente el perfil de velocidad a la salida.
• Mostrar la longitud de entrada necesaria hasta que se obtiene perfil desarrollado.
NOTA: Han de validarse correctamente los resultados. Los valores de las propiedades del flujo no
corresponden a ningún fluido típico. Han sido escogidos para que el flujo sea laminar. El parámetro
clave en este caso es el número de Reynolds.
Análisis preliminar.
Se trata de ver cómo la capa límite viscosa va creciendo desde la entrada a lo largo de la
tubería. De hecho, ha de crecer hasta que rellene la tubería por completo (para ello, la tubería ha de
tener la longitud suficiente). Cuando esto ocurre, el flujo está completamente desarrollado (fully
developed) y ya no hay variación del perfil de velocidad en la dirección del tubo de ahí en adelante
(ver gráfico de la portada).
Existe una solución analítica exacta de este tipo de flujo, que habréis visto en las clases
teóricas. Ahora vamos a tratar de comparar dicha solución con los resultados numéricos que nos de
el programa. Así que no es mala idea revisar los apuntes de teoría para validar los perfiles de
velocidad y los coeficientes de fricción que obtengamos en la simulación.
A continuación se explican detalladamente los pasos para crear la geometría y hacer la malla
en GAMBIT (el programa de generación de casos de FLUENT) para a continuación leer el caso en
el solver FLUENT y ejecutar la solución del ejercicio.
¡Buena suerte!
1
2.- GENERACIÓN DE LA GEOMETRÍA EN GAMBIT
Estrategia a seguir para crear la geometría.
Para crear el rectángulo que simulará la sección longitudinal de la tubería, vamos a crear en
primer lugar los 4 vértices (vertex) en las 4 esquinas. A continuación se unirán los vértices
correspondientes mediante líneas para crear los lados (edges) del rectángulo. Por último, crearemos
la superficie (face) que corresponde al área encerrada por los 4 lados. En el punto 3 de esta guía, se
mallará la superficie del rectángulo.
En casos tridimensionales, antes de mallar habría que generar volúmenes a partir de las
superficies. Por tanto, nótese la estructura jerárquica que presenta el programa: vertex -> edges ->
faces -> volumes.
Crear un directorio de trabajo.
Ha de crear una carpeta “Tubería” dentro del subdirectorio que sea conveniente.
Utilizaremos esa carpeta de trabajo para almacenar los archivos y resultados que se vayan
generando a lo largo de esta sesión.
Inicie GAMBIT.
Inicio > Ejecutar. Escriba cmd y presione al enter. Se abre una ventada de MS2 en la que
debe escribir: gambit –id tuberia y luego ejecute pulsando intro. Si esto no funciona,
escriba en la ventana el path completo:
c:\fluent.inc\ntbin\ntx86\gambit –id tuberia
Todas estas órdenes sirven para abrir el programa. Además, al indicarle como identificador
(-id) el nombre “tubería”, le estamos indicando al GAMBIT que utilice ese nombre como prefijo del
archivo que vamos a construir a lo largo de esta sesión.
Interfaz del GAMBIT.
La interfaz de GAMBIT se compone de las siguientes partes:
• Barra principal. Observe que el nombre tuberia aparece tras el ID en la barra.
• Barra de herarramientas. A lo largo de este caso vamos a emplear muchas veces este panel.
Fíjese que cada uno de los botones superiores, al ser seleccionados, dan paso a diferentes
“sub-botones”. Por ejemplo, los “sub-controles” de Geometría son los que aparecen en la
captura de pantalla:
2
• Controles globales. Controlan la apariencia del modelo, las vistas o el zoom:
En este menú aparecen controles tales como Fit to Screen (ajustar a pantalla) o Undo
(deshacer) que son de gran utilidad mientras se crea la geometría y la malla.
• Pantalla gráfica. Es la ventana en la cual vamos a ver los progresos de nuestro modelo.
3
• Panel de descripción de GAMBIT. Este panel contiene la descripción de los botones y

objetos que está señalando el ratón en todo momento. Mueve el ratón sobre distintos botones
para observar el texto descriptivo correspondiente.
• Ventana de trascripción de GAMBIT. Es la ventana en la que van quedando reflejados los

sucesivos comandos que son ejecutados por el usuario. Si se hace clic sobre la pequeña
flecha de la derecha, arriba se puede minimizar y maximizar esta ventana.
Selección del SOLVER.
Especifique que la malla que va a crear es para ser utilizada con FLUENT 6.0
Main Menu > Solver > FLUENT 5/6
Verifique que ha quedado correctamente registrado a través de la ventana de trascripción, en
la cual puede leer: Command> solver select “FLUENT 5/6”
NOTA: Los tipos de condiciones de contorno que pueden elegirse en el siguiente paso, dependen
del solver que se haya seleccionado en esta opción.
Vamos a asumir que el flujo es AXISIMÉTRICO (esto es, tiene simetría circunferencial con
todas las secciones longitudinales presentando las mismas características). Por tanto, el dominio del
problema es:
D
0≤r ≤ ,0≤ x≤ L
2
donde r y x son las coordenadas radiales y axiales respectivamente.
Origen de coordenadas.
Se colocará el origen del sistema de coordenadas en la esquina inferior izquierda del

rectángulo. Las coordenadas de las esquinas se muestran en la figura inferior:
4
(L,D/2)
(0,D/2)
(0,0) (L,0)
En primer lugar se crearán los cuatro vértices en las cuatro esquinas, juntándose los vértices
adyacentes para crear las líneas del rectángulo. De esta forma se creará el área subtendida por el
rectángulo.
Creando vértices.
Encuentre los botones que se indican a continuación, señalando cada uno de ellos con el
ratón y observando la descripción que se hace de ellos en la ventana de descripción.
Operation Toolpad > Geometry Command Button > Vertex Command Button > Create
Vertex
NOTA: Como puede comprobar, el botón de “Create Vertex” ya está seleccionado por defecto.
Se crea el vértice que representa la esquina inferior izquierda del rectángulo. Para ello, junto
a x: introduzca el valor 0; junto a y: introduzca el valor 0 y junto a z: introduzca el valor 0. Clic
Apply para aplicar. Con esto, se crea el vértice (0,0,0).
En la ventana de transcripción, GAMBIT indica que ha creado el vértice “vertex.1”. Los

vértices (si no se les dan nombre en el campo label) son creados de forma correlativa en el orden en
que son creados.
Repita el proceso para crear los 3 vértices restantes:
Vertex.2: (0,D/2,0)
Vertex.3: (L,D/2,0)
Vertex.4: (L,0,0)
5
Como puede ver, en problemas bidimensionales, la coordenada z se deja siempre con el

valor 0 por defecto.
Operation Toolpad > Global Control > Fit to Window Button
Este comando ajusta los 4 vértices del rectángulo que se acaba de crear a la ventana gráfica
para que podamos verlos en pantalla (se muestra sólo un trozo de la pantalla en la captura
siguiente):
Creando lados.
Ahora han de conectarse entre sí las parejas apropiadas de puntos que constituyan inicio y
fin de un lado del rectángulo.
IMPORTANTE: Para seleccionar cualquier objeto en GAMBIT, presione la tecla
MAYÚSCULAS (Shift) y a la vez haga clic en el objeto que quiera seleccionar. Cuando un objeto
queda seleccionado, aparece en rojo en la pantalla gráfica. Luego, suelte el botón de
MAYÚSCULAS. Se pueden chequear los vértices que han sido seleccionados si se hace clic sobre
la flechita junto al campo Vertices:
Haciendo esto, se abre una ventana desplegable donde aparece una relación de todos los
objetos que han sido seleccionados. Además, los vértices se pueden poner y quitar de la lista
utilizando los campos Available (disponibles) y Picked (seleccionados) y las flechas de izquierda y
derecha para llevarlos de una lista a otra y viceversa.
6
Una vez que se haya seleccionado la pareja correcta de vértices, es el momento de hacer clic
en Apply para generar el lado. Repita el proceso tres veces más hasta que disponga de los cuatro
lados del rectángulo.
Crear la cara.
Operation Toolpad > Geometry Command Button > Face Command Button >
Form Face
Para crear la cara que queda encerrada por los 4 lados es preciso seleccionar los 4 lados y
que éstos formen un lazo cerrado. Si no es así, el programa es incapaz de generar la superficie. Para
seleccionar los lados puede presionar la tecla MAYÚSCULAS y a la vez ir haciendo clic sobre
cada línea y al final soltar la tecla una vez queden seleccionadas las 4 líneas en rojo. También es
posible seleccionar utilizando la ventana desplegable que aparece junto a la flechita de la derecha:
Con esto, se activa la ventana de lista de lados. Haciendo clic en All-> se seleccionan todos
los lados de golpe. Haga clic en Close.
7
Finalmente haga clic en Apply para crear la cara.
3.- GENERACIÓN DEL MALLADO EN GAMBIT
A continuación se detallan los pasos para crear una malla en la cara rectangular con 250
divisiones en la dirección axial y 5 divisiones en la dirección radial. Se mallarán en primer lugar los
cuatro lados y a continuación se mallará la cara. El espaciado en la malla se especifica a través del
mallado de los lados.
Mallado de los lados.
Operation Toolpad > Mesh Command Button > Edge Command Button > Mesh
Edges
Haga MAYÚSCULAS + clic como antes o utilice la ventana desplegable asociada a la lista
de lados para seleccionar los dos lados verticales. Para acceder más fácilmente a los lados, se puede
hacer un zoom ventana de la zona de interés. Para ello, manteniendo pulsada la tecla de CTRL, se
hace clic y se arrastra el ratón sobre una determinada zona (la zona que quiero ver ampliada) y
continuación se suelta el botón de CTRL. Para retornar a la vista completa del modelo, basta con
volver a clicar en: Control Toolpad > Fit to Window Button.
Una vez que el eje vertical ha sido seleccionado, seleccione la opción Interval Count del
botón desplegable que por defecto pone Interval Size dentro de la ventana de mallado de lados.
Introduzca el valor de 5 como número de intervalos de malla.
Haga clic en Apply. Los nodos de la malla aparecen ahora en el lado, indicando que se han
tomado 5 divisiones intermedias.
8
Repita el proceso con los lados horizontales, pero empleando un intervalo de 250 nodos.
Ahora que los lados están mallados, podemos proceder a crear una malla 2D para la cara.
Mallar la cara.
Operation Toolpad > Mesh Command Button > Face Command Button > Mesh
Faces
Haga MAYÚSCULAS + clic en la cara para seleccionar y luego haga clic en Apply.
4.- CONDICIONES DE CONTORNO Y EXPORTACIÓN DESDE GAMBIT
Crear los tipos de condiciones de contorno.
El esquema adjunto nos muestra qué tipo de condición de contorno hemos de fijar en el
modelo. El lado izquierdo es la entrada de la tubería, el lado derecho es la salida, el superior será la
pared y el inferior un eje axisimétrico.
Wall
Inlet Outlet
Axis
Operation Toolpad > Zones Command Button > Specify Boundary Types Command
Button
Con esto se activará el panel de especificación de condiciones de contorno. Para empezar,

vamos a fijar que el lado izquierdo es la entrada del fluido. Para ello, en Entity: seleccionamos la
opción Edges para que GAMBIT sepa que queremos seleccionar un lado (por defecto trae activo
Face).
Ahora seleccionamos el lado izquierdo como siempre, haciendo MAYÚSCULAS + clic

sobre el lado. El lado seleccionado aparece en el recuadro amarillo que está bajo el botón de Entity.
Además, podemos darle nombre dentro del campo Name:, por ejemplo inlet.
Para fijar el Type:, seleccionamos VELOCITY_INLET.
9
Hacemos clic en Apply. Con eso, la condición de contorno queda carga dentro del recuadro
Name/Type, tal como se ve en la captura siguiente:
A continuación se repite el proceso con el resto de lados según la siguiente tabla:
Lado (edge) Nombre (name) Tipo (type)

Izquierdo Inlet VELOCITY_INLET
Derecho Outlet PRESSURE_OUTLET
Arriba Wall WALL
Abajo Centerline AXIS
Si se ha hecho todo correctamente, la lista final de condiciones de contorno debería figurar

de la siguiente manera en el programa:
10
Salvar y exportar.
Main Menu > File > Save As… Especifique el directorio donde quiere guardar el fichero de
GAMBIT. Note que la extensión es *.dbs
Main Menu > File > Export > Mesh… Con esta instrucción se exporta el fichero para que
pueda ser interpretado por Fluent a continuación. Escriba tuberia.msh como nombre de archivo
(File Name). Seleccione Export 2D Mesh, ya que efectivamente el caso es bidimensional. Haga
clic en Accept.
Compruebe que un archivo tuberia.msh ha sido creado en su directorio de trabajo.
Cierre el programa GAMBIT y pase a abrir el programa FLUENT.
5.- INICIALIZACIÓN DEL CASO EN FLUENT. CASO LAMINAR
Ejecutar el programa FLUENT.
Start > Programs > Fluent Inc > FLUENT 6.0 o bien ejecute el icono que debe aparecer en
el escritorio del PC. Seleccione la opción 2ddp dentro de la lista de opciones y luego haga clic en
Run.
La opción 2ddp se utiliza para seleccionar el solver bidimensional de doble precisión. En el
solver de doble-precisión, cada número decimal está representado por 64 bits, en vez de los 32 bits
que se emplean en el solver de precisión sencilla. La opción de doble precisión no aumenta sólo la
precisión sino también el rango máximo y mínimo que pueden alcanzar las magnitudes que se van a
representar. Obviamente, esta opción requiere de más memoria de cálculo.
Importar la malla.
Main Menu > File > Read > Case…

Utilice el explorer que nos aparece para acceder a su directorio de trabajo y seleccionar el
fichero tuberia.msh. Obviamente, ese es el fichero que fue creado con el programa anterior
GAMBIT. Al cargarlo, FLUENT nos da información sobre las características del mallado.
Compruebe el número de nodos, caras (de los distintos tipos) y celdas. Hay 1250 celdas
rectangulares (quadrilateral cells). Obviamente es lo esperable puesto que se usaron 5 divisiones
radiales y 250 en la dirección axial, así que el número total de celdas no podía ser otro que
5x250=1250.
Además, es interesante ver las zonas que se han cargado. Podemos comprobar cómo
aparecen las cuatro condiciones de contorno que definimos con anterioridad: inlet, outlet, wall y
centerline.
11
Validar y dibujar la malla.
En primer lugar, hemos de chequear el mallado para asegurarnos que no ha habido errores ni
en la generación, ni en la exportación ni en la importación.
Main Menu > Grid > Check
En caso de que hubiesen errores, deberían ser detectados por FLUENT en esta operación.
Compruebe el resultado del chequeo. Si todo esta correcto, proceda a comprobar el tamaño del
mallado:
Main Menu > Grid > Info > Size
Los siguientes resultados deberían ser proporcionados por el programa:
A continuación dibuje la malla: Main Menu > Display > Grid …

Asegúrese que los 5 objetos que aparecen bajo el epígrafe Surfaces están seleccionados. A
continuación haga clic en Display. Una ventana gráfica se abre y la malla aparece dibujada en ella.
Ya puede hacer clic en el botón Close de la ventana Grid Display para ganar espacio en el
escritorio. Como puede comprobar, la ventana gráfica no se cierra y permanece visible.
Las más importantes operaciones que se pueden hacer dentro de la ventana gráfica se
resumen a continuación:
• Traslación. El mallado se puede mover en cualquier dirección (centrar el dibujo completo o
solo una zona determinada) si se mantiene pulsado el botón izquierdo del ratón y a la vez
se mueve el ratón en la dirección que se quiera.
• Zoom In. Mantenga pulsado el botón central del ratón y arrastre una ventana desde la
esquina superior izquierda a la esquina inferior derecha de la zona que quiere observar.
Luego suelte y se ejecutará el zoom de esa zona, ajustándolo al tamaño máximo de la
ventana gráfica.
12
• Zoom Out. Mantenga pulsado el botón central del ratón y arrastre una ventana en
cualquier lado desde la esquina inferior derecha a la esquina superior izquierda (inverso al
anterior). Luego suelte.
Practique estas operaciones para obtener una vista detallada de la entrada tal y como se
muestra en la figura siguiente:
Obviamente, puede mostrar sólo aquellas partes del mallado que le interese. Para ello, basta
con seleccionar sólo las superficies que quiere mostrar:
Por comodidad, el botón que está junto a Surfaces selecciona todos los objetos, mientras
que el botón deselecciona todas a la vez.
Cierre la ventana gráfica para continuar.
Definición del SOLVER.
Main Menu > Define > Models > Solver
13
Elija la opción de Axisymmetric bajo el cuadro Space. Se emplearán las otras opciones que
por defecto indica el programa, como el tipo segregated, implicit formulation, steady flow y
absolute velocity formulation. Haga clic en Ok.
Main Menu > Define > Models > Viscous

La opción de flujo laminar está establecida por defecto. Compruebe que efectivamente sus
datos corresponden con una situación de flujo laminar. En ese caso, no necesita cambiar nada en
este menú. Haga clic en Cancel.
Main Menu > Define > Models > Energy

En caso de flujo incompresible, la ecuación de la energía queda desacoplada de las
ecuaciones de continuidad y del momento. Es preciso resolver la ecuación de la energía únicamente
si estamos interesados en conocer la distribución de temperatura. En este caso, sin embargo, no
tratamos ningún aspecto de transmisión de calor. Por tanto, se deja la ecuación de la energía
desactivada y salimos haciendo clic en Cancel.
Definición de las propiedades del fluido.
Main Menu > Define > Materials

Cambie la propiedad Density a valor 1.0 y la de Viscosity a 2x10-3. Estos son los valores que
ya fijamos al principio del ejemplo. Los tomaremos como valores constantes:
Haga clic en Change/Create.
Definición de Condiciones de Operación.
Main Menu > Define > Operating Conditions …

Para todo tipo de flujos, FLUENT utiliza internamente presión manométrica (relativa a la
atmosférica). De todos modos, este modo de proceder exige que en algún momento se fije el valor
14
absoluto de la presión a la cual quedan referidos todos los cálculos a partir de ese momento. Se
emplerará el valor por defecto de 1 atm (101325 Pa) como Operating Pressure. Haga clic en
Cancel para mantener el valor por defecto.
Definición de Condiciones de Contorno.
A continuación se fijarán los valores de velocidad a la entrada y presión a la salida.
Main Menu > Define > Boundary Conditions …

Obsérvese que en este panel aparecen las 4 cuatro condiciones de contorno que definimos en
GAMBIT como zonas en la izquierda de la ventana de boundary conditions. La zona centerline
debe estar seleccionada por defecto. Compruebe que así es, así como que el tipo seleccionado es
axis. Después haga clic en Set… Fíjese cómo no hay nada que especificar para este tipo de
condición. A continuación, haga clic en Ok.
Muévase a lo largo de la lista y seleccione ahora la condición inlet. FLUENT indica que el
tipo de esta condición de contorno es velocity-inlet. Recuerde que el tipo de condición fue ya fijada
en GAMBIT. De todas formas, es posible cambiar aquí el tipo de condición que se fijó en GAMBIT
si es necesario. Basta con seleccionar otro tipo de condición en esta ventana y validar la selección.
Haga clic en Set…
Introduzca el valor de velocidad V que tenga en los datos de partida en el campo Velocity
Magnitude. Haga clic en Ok. Esta operación fija el valor de velocidad del flujo que entra al
dominio desde la izquierda.
La presión (absoluta) que se tiene a la salida es 1 atm. Puesto que la presión de operación se
ha dejado fijada en 1 atmósfera, la presión manométrica a la salida (outlet gauge pressire) ha de ser
cero. Seleccione outlet bajo el campo Zone. El Type de esta condición de contorno es pressure-
outlet. Haga clic en Cancel para mantener el valor por defecto.
15
Finalmente, haga clic en wall y asegúrese que su tipo es wall. Vaya haciendo clic por todas
las pestañas y compruebe que sólo parámetros sobre la ecuación del momento son accesibles. Esto
no será así en posteriores modificaciones, así que recuerde la localización de estas opciones para
después. Haga clic en Ok.
Haga clic en Close para cerrar el menú de Boundary Conditions.
6.- RESOLUCIÓN
Para la resolución se empleará un esquema de discretización de segundo orden.
Main Menu > Solve > Controls > Solution…

Cambie en la opción de Momentum a Second Order Upwind.
Haga clic en Ok.
Inicialización de la solución.
Hay que inicializar el campo de velocidades y presiones para llevar a cabo el método
iterativo de resolución. En este caso, vamos a fijarle a todo el dominio los mismos valores que se
tienen en la condición de contorno de entrada.
Main Menu > Solve > Initialize > Initialize…

En el menu de Solution Initialization que aparece, elegimos inlet dentro del campo Compute
From. De esta manera, la componente axial de la velocidad se fijará al mismo valor de velocidad V
que se tiene de dato de entrada. Asimismo, la componente radial se fijará en 0 y se supondrá una
presión manométrica de valor también 0. Estos valores son tomados de la propia condición de
velocidad a la entrada.
Haga clic en Init. Con esto, se completa la inicialización del modelo.
16
Fijar criterios de convergencia.
Mientras resuelve las ecuaciones de forma iterativa, FLUENT va mostrando en pantalla un

valor de residuo asociado a cada ecuación de gobierno del flujo que está resolviendo. El residuo es
una medida de lo bien que la solución actual cumple las ecuaciones de gobierno discretizadas. En
nuestro caso, vamos a iterar la solución hasta que los valores del residuo se sitúen por debajo del
valor umbral de 10-6.
Main Menu > Solve > Monitors > Residual…

Cambie el valor de residuo bajo el campo Convergence Criterion al valor 1e-06 tanto para
la ecuación de continuidad (continuity) como para las componentes de la ecuación del momento (x-
momentum, y-momentum).
Además, bajo el cuadro Options, seleccione la opción Plot. Con esto, conseguirá que
FLUENT muestre los residuos en pantalla mientras va resolviendo el modelo. Mantenga también
activa la opción Print (irá escribiendo los valores en pantalla).
17
Haga clic en Ok. Con esto se termina todo el proceso de inicialización y definición del
modelo a resolver. Ahora, salve el trabajo:
Main Menu > File > Write > Case&Data …

Escriba el nombre tuberia.cas como Case File. Haga clic en Ok. Compruebe que
efectivamente FLUENT crea dicho archivo en su directorio de trabajo. Si por cualquier razón, cierra
el programa FLUENT, podrá recuperar todo su trabajo simplemente leyendo el archivo
tuberia.cas.
Proceso Iterativo hasta obtener convergencia.
Comience a calcular, fijando un total de 300 iteraciones:
Main Menu > Solve > Iterate…

En la ventana Iterate Window, cambie el número de iteraciones a 300. Haga clic en Iterate.
Los residuos en cada iteración se muestran en pantalla gráficamente y en la ventana de comandos de
FLUENT:
En este caso, con los valores que se hayan fijado de V,D y L, los residuos caen por debajo
de 1e-06 en 97 iteraciones:
Salve la solución a un archivo tipo *.dat. (archivo de datos). Los ficheros tipo *.cas
contienen la información del modelo (geometría, malla, condiciones de contorno…) mientras que
los de tipo *.dat contienen los campos de velocidad, presión, etc…
18
Main Menu > File > Write > Data…

Introduzca tuberia.dat en el campo Data File y haga clic en Ok. Compruebe
nuevamente que dicho archivo queda almacenado en su directorio de trabajo. Nuevamente, podrá
recuperar la solución del modelo siempre que quiera.
7.- ANÁLISIS DE RESULTADOS
En este apartado se mostrará cómo generar gráficos y mapas de vectores y de contornos con
el objeto de poder presentar los resultados que nos proporciona FLUENT.
Velocidad en la línea central (centerline)
Se va a hacer un gráfico que muestre la evolución de la componente axial de la velocidad a

lo largo del eje de la tubería (centerline).
Main Menu > Plot > XYPlot…

Asegúrese de que la opción Position on X Axis está activada dentro del campo Options.
Además, X se fija con el valor 1 e Y con el valor 0 para indicar el vector director de la línea que se
quiere evaluar (centerline); todo esto en el recuadro Plot Direction. De esta manera, FLUENT
dibuja la coordenada x en la abcisa del gráfico.
Bajo Y Axis Function, seleccione Velocity… y después, dentro de la subcategoría, elija la

componente axial, Axial Velocity.
Es importante fijarse en que X Axis Function e Y Axis Function describen el eje x y el eje y
del gráfico respectivamente, y no las direcciones x e y de la tubería.
Finalmente, seleccione centerline bajo la categoría Surfaces puesto que efectivamente lo

que queremos es que nos de el perfil de velocidad a lo largo de dicha línea. Con esto se completa la
preparación de los parámetros de dibujo.
Haga clic en Plot. Se abre una nueva ventana gráfica mostrando el gráfico de tipo xy qus
acabamos de solicitar:
19
Podemos observar que en el gráfico se alcanza un valor constante de velocidad a partir de

una determinada distancia desde la entrada. Esa es la zona de flujo completamente desarrollado.
Para apreciar mejor ese valor, vamos a cambiar el rango de los ejes del gráfico. En el menú
de Solution XY Plot, hacemos clic en Axes… en Options deseleccionamos el AutoRange. Ahora,
bajo Range están accesibles los máximos y mínimos tanto de x como de y. Introduzca un valor
mínimo de x y un máximo de x por la zona donde se comienza a estabilizar el perfil. En el caso de
la figura, elegimos entre 10 y 20. Elija usted en función de los resultados que haya obtenido
siguiendo el mismo criterio.
Además, active las grid lines que ayudarán a estimar por dónde el flujo se vuelve
completamente desarrollado. Para ello, haga clic en Major Rules y Minor Rules en Options. Luego
haga clic en Apply.
Ahora repita lo mismo pero acotando el rango del eje Y. En nuestro ejemplo fijaremos el
rango de y entre 3.6 y 4 para el Minimum y Maximum en Range. Seleccione también Major Rules
y Minor Rules en Options. Luego haga clic en Apply. Finalmente, haga clic en Close.
20
Si observamos la gráfica con los nuevos límites, podemos determinar de forma más precisa
que la zona completamente desarrollada comienza sobre los 17 metros para una velocidad de 3.84
m/s.
Guardar el gráfico.
En la ventana de Solution XY Plot, active la casilla de Write to File en Options. Entonces el

botón Plot se convierte en Write. Haga clic en él. A continuación escriba el nombre del fichero,
vel.xy como el XY File Name y haga clic en Ok. Compruebe que se ha creado ese archivo en su
directorio de trabajo.
A continuación guarde una imagen del mismo gráfico: Mantenga las ventanas de Solution
XY Plot y la propia ventana gráfica donde está el gráfico y vaya a:
Main Menu > File > Hardcopy…
En Format elija una de estas dos opciones:
• TIFF. Proporciona una imagen de gran resolución. Desgraciadamente, suele ocupar mucho
espacio en disco, así que emplee esta opción si dispone de suficiente memoria en disco.
• JPEG. Es una imagen de menor tamaño y compatible con cualquier programa gráfico y
visor. Sin embargo, la calidad de la imagen es menor.
Una vez seleccionado el tipo de imagen, haga clic en Save… Nombre la imagen con la
extensión apropiada: vel.tif o vel.jpg. Como siempre, compruebe que el archivo ha sido
creado en su directorio de trabajo.
Coeficiente de fricción en la pared.
La ayuda de FLUENT nos proporciona la definición del coeficiente de fricción:
21
El cálculo de ese coeficiente de fricción depende de unos valores de referencia que es

necesario fijar en el modelo. Para ello:
Main Menu > Report > Referente Values…

Seleccionamos inlet dentro de la casilla de Compute From para indicarle al programa que
tome como valores de referencia aquellos que se tienen a la entrada. Compruebe que toma los
valores correctos (si no es así, introdúzcalos vía teclado) y haga clic en Ok.
Vuelva a la ventana de Solution XY Plot y deseleccione la opción Write to File, puesto que
ahora pretendemos realizar un dibujo nuevo (solo guardaremos cuando el gráfico sea definitivo). El
resto de Options y la Direction Plot se dejan como estaban pues seguimos queriendo hacer un
gráfico longitudinal.
Bajo la opción de Y Axis Function, haga clic en Wall Fluyes y luego elija la subcategoría
Skin Friction Coefficient. En Surfaces elija wall y deseleccione centerline pinchando sobre ella.
Así mismo, resetee los rangos de los ejes. Retome Axes… y fije de nuevo AutoRange tanto para X
como para Y. Haga clic en Apply, luego en Close y finalmente haga Plot:
22
En este gráfico se observa que la zona completamente desarrollada se inicia alrededor de los
9 metros y que el coeficiente en esa zona es aproximadamente 0.38. Si se compara este valor con el
coeficiente que se obtiene de forma analítica se vería una gran discrepancia. Además, vemos que la
zona desarrollada se predice de forma muy diferente con esta variable que con el perfil de velocidad
anterior. Sobre esto volveremos más adelante.
Salve los datos de este plot haciendo Write to File en options y guardelo haciendo clic en
Write bajo el nombre cf.xy.
Perfil de velocidad.
A continuación se mostrará la velocidad a la salida del tubo en función del radio del tubo.
Para hacer esto, es necesario fijar como eje y del gráfico el eje y del tubo (la dirección radial).
Para poner en el eje y del gráfico la coordenada radial de la tubería, hay que deseleccionar la
opción Position on X Axis en options y elegir en su lugar Position on Y Axis. Además, para fijar
que la dirección radial será la que se quiera dibujar hay que cambiar en Plot Direction la X a 0 y la
Y a 1. Para dibujar la componente axial en el eje x del gráfico seleccione en X Axis Function la
categoría Velocity y luego la subcategoría Axial Velocity. Como se quiere pintar esto en la salida,
elija outlet dentro de Surfaces y deseleccione todo el resto. Recuerde cambiar los eje x e y a
Autorange. Finalmente, desactive Write to File (si ha quedado activado previamente) en Options
para ver el gráfico y luego haga clic en Plot
23
Como se puede comprobar, presenta un perfil parabólico, tal y como era de esperar.
Ahora, proceda a guardar este gráfico como en anteriores ocasiones. Para ello, simplemente,
pulse en Write to File en Options y luego haga clic en Write… A continuación guardelo con el
nombre profile.xy para XY File y haga clic en Ok.
Además, para ver cómo va evolucionando el perfil de velocidades en la zona de desarrollo,

vamos a añadir los perfiles en las posiciones x/D=3 y x/D=6 en el dibujo anterior. Para ello hay que
crear un par de líneas en esas posiciones utilizando la herramienta line/rake.
Main Menu > Surface > Line/Rake

Crearemos una línea recta que vaya de las posiciones (x0,y0)=(3D,0) a (x1,y1)=(3D,D/2).
Para ello, seleccione Line Tool en Options. Introduzca los valores apropiados para x0, y0, x1, y1.
Nombre line1 en el campo New Surface Name. Haga clic en Create.
Para poder ver la línea que acaba de crearse, ejecute:

Main Menu > Display > Grid …
Observe que line1 aparece en la lista de superficies. Seleccione todas las superficies
excepto el default-interior. Haga clic en Display. De esta forma se dibujan todas las superficies
pero no así la malla del modelo. Haga zoom en la región próxima a la entrada para ver la línea
recién creada en x/D=3. Como ve, line1 es la línea blanca vertical en la figura adjunta:
24
De forma análoga se crea la línea para x/D=6, esta vez definiendo las posiciones
(x0,y0)=(6D,0) a (x1,y1)=(6D,D/2). Para ello, seleccione Line Tool en Options. Introduzca los
valores apropiados para x0, y0, x1, y1. Nombre line2 en el campo New Surface Name. Haga clic
en Create.
Finalmente, ya es posible dibujar los perfiles de velocidad en las posiciones x/D=3 y x/D=6
junto con el perfil a la salida. En la ventana Solution XY Plot utilice los mismos parámetros que en
el caso anterior. En Surfaces, añada line1 y line2 a la sección outlet. Seleccione Node Values en
Options y haga clic en Plot. Los símbolos que pinta el FLUENT para cada curva son distintos. De
todas formas, es posible configurar estos símbolos accediendo al botón Curves. Se deja al alumno la
iniciativa para que configure e investigue estos controles por su cuenta. Finalmente, se obtendría un
gráfico como éste:
Se puede ir observando el paulatino desarrollo de la capa límite desde la entrada hasta la

zona completamente desarrollada. Sin embargo, se ha elegido un mallado que no es lo
suficientemente fino cercano a la pared, lo cual imposibilita una buena descripción de la generación
de la capa límite (apenas hay diferencias entre las curvas cerca de la pared).
25
Vectores de velocidad
Es posible dibujar los vectores de velocidad en todo el dominio o bien en una determinada
superficie. Vamos a dibujar los vectores justo en la zona de entrada para ver cómo se desarrolla el
flujo en esa zona.
Main Menu > Display > Vectors …> Display

Haremos un zoom en la región cercana a la entrada. La longitud y el color de las flechas
representan la magnitud de la velocidad. Cada vector es más visible si se hacen las flechas un poco
más cortas de la siguiente forma: cambiar Scale a 0.4 en el menú Vectors y en Style elija filled-
arrow. Luego haga clic en Display. Además, es posible reflejar el dibujo a partir del eje del modelo
para obtener una vista completa de la sección longitudinal. Para ello:
Main Menu > Display > Views

Bajo la opción de Mirror Planes solamente la superficie axis está accesible, puesto que es la
única condición de simetría existente en este caso. Selecciónela y haga clic en Apply. Cierre la
ventana Views. Habrá obtenido algo similar a esto:
Los vectores de velocidad proporcionan una bonita imagen de cómo el flujo va

desarrollándose en el interior de la tubería. Al ir creciendo la capa límite, el flujo que va pegado a la
pared se ve frenado cada vez más, como consecuencia del efecto de la fricción viscosa.
8.- REFINADO DEL MODELO
Es muy importante asegurarse de que los resultados que se obtienen son completamente
independientes del tipo de malla que se ha utilizado para resolverlos. Para saber si los resultados
que hemos obtenido son realmente independientes de la malla, se ha de rehacer el caso empleando
un mallado diferente y comparar los resultados. Por tanto, a continuación se van a rehacer los
cálculos previos empleando una malla de 15x250 y otra de 25x250. A continuación se detallan los
pasos que han de hacerse para generar esas nuevas mallas rápidamente.
26
Modificar las mallas con GAMBIT.
El caso anterior con una malla de 5x250 celdas está almacenado en el directorio de trabajo
con el nombre de tuberia.cas. Vamos a generar dos copias de dicho archivo con el explorador
de windows y cambiarlas de nombre: tuberia2.cas y tuberia3.cas.A continuación,
vamos a cargar con el GAMBIT el archivo tuberia2.cas y procederemos a hacerle las
modificaciones pertinentes para dotarle de una malla de 15x250. Después repetiremos el proceso
con el caso tuberia3.cas.
Se procede a explicar el caso de tuberia2.cas. Hay que borrar el mallado de la cara,

modificar la distribución de nodos de las líneas verticales y finalmente remallar la cara. Para
remallar la cara ejecutamos:
Operation Toolpad > Mesh Command Button > Face Command Button > Delete Face
Meshes
En la ventana Delete Face Meshes que aparece, hemos de deseleccionar la opción Remove
Unused Lower Mesh. Con esta opción, le indicamos a FLUENT que elimine la malla de la
superficie, pero que no borre la malla de los lados que la componen. Puesto que vamos a remallar
sólo dos lados, no es necesario rehacer la malla de los lados horizontales. Seleccione la única cara
que tiene el modelo como siempre con MAYÚSCULAS + clic sobre ella y luego haga clic en
Apply.
Para cambiar el número de nodos de 5 a 15, seleccionamos:

Operation Toolpad > Mesh Command Button > Edge Command Button > Mesh Edges
Seleccione los dos tramos verticales como siempre pulsando la tecla de MAYÚSCULAS y
a la vez haciendo clic sobre ellos. Tras soltar la tecla de MAYÚSCULAS, seleccione el Interval
Count en el botón de Interval Size dentro de la opción Spacing. Introduzca el valor numérico de
15.
Asegúrese de que la opción Remove Old Mesh está activada. De esta forma, garantizaremos
que la antigua distribución de nodos queda borrada antes de aplicar la nueva. Haga clic en Apply.
27
Recuerde que puede hacer zoom siempre que quiera manteniendo pulsada la tecla CTRL y
arrastrando una ventana sobre la zona que quiera aumentar de tamaño, y después soltando la tecla
de CTRL. Haga esto y asegúrese de que realmente el GAMBIT ha colocado 15 nodos en las líneas.
A continuación rehaga la malla:

Operation Toolpad > Mesh Command Button > Face Command Button > Mesh Faces
Haga MAYÚSCULAS + clic sobre la superficie en la ventana gráfica y luego haga clic en
Apply.
Finalmente, salve y exporte la malla de tuberia2.dbs y tuberia2.msh.

Main Menu > File > Save
Main Menu > File > Export > Mesh…
Teclee tuberia2.dbs en el campo FileName :. Seleccione Export 2D Mesh y haga clic
en Accept.
RECUERDE: Repita este último epígrafe de modificación de la malla para generar un caso
de 25x250 celdas y guardelo con el nombre tuberia3.dbs
Análisis de los casos con malla más fina
Repita los pasos 5 y 6 de este guión para los mallados refinados de 15x250 y 25x250 celdas.
(Es una tarea repetitiva pero constituye una buena práctica en CFD).
Una vez que tenga los casos tuberia2.cas y tuberia3.cas bien convergidos y
solucionados pasaremos a comparar los resultados de los distintos tipos de malla entre si.
Obviamente, a mayor número de celdas, más iteraciones son necesarias para conseguir la
convergencia. En el caso de este ejemplo, tuberia2.cas convergió tras 153 iteraciones
mientras que tuberia3.cas lo hizo tras 267.
Cargue ahora el caso resuelto tuberia2.cas y obtenga los ficheros de datos

vel2.xy; cf2.xy y profile2.xy siguiendo el mismo proceso que se mostró en el
apartado 7 del presente guión.
A continuación vamos a comparar los 3 casos sobre un mismo gráfico. Para ello cargamos
ahora en FLUENT el caso resuelto tuberia3.cas. Nuevamente vamos a dibujar el perfil
longitudinal de velocidad, el coeficiente de fricción y el perfil parabólico a la salida, pero se van a
cargar también los ficheros anteriores para verlos conjuntamente.
Una vez cargado el caso tuberia3.cas dibuje el perfil de velocidad a lo largo del eje
centerline como ya sabe. A continuación cargue los ficheros vel.xy y vel2.xy del directorio
28
de trabajo. Para ello, una vez que ha dibujado el perfil del caso tuberia3.cas, haga clic en
Load File… dentro de la ventana Solution XY Plot. Seleccione vel.xy dentro de su directorio y
luego haga clic en Ok. Repita la operación con el archivo vel2.xy y finalmente haga clic en Plot.
Ajuste la escala del gráfico para centrar la zona donde comienza a haber flujo
completamente desarrollado. También puede jugar con la apariencia de las curvas, pulsando en
Curves… para conseguir que se dibujen como líneas sólidas. En definitiva, ha debido obtener algo
parecido a lo siguiente:
La línea roja representa el caso de la malla 5x250, mientras que los casos verde y blanco
corresponden a las mallas refinadas 15x250 y 25x250. El valor de la velocidad en la zona
completamente desarrollada para la malla más fina es prácticamente 3.995 m/s. La línea verde (com
menos malla) también da unos excelentes resultados: 3.975 m/s. Sin embargo, la malla original
demuestra un valor más alejado y un menor grado de exactitud: 3.84 m/s. Nótese que el valor
analítico exacto de este problema es 4 m/s.
NOTA: El valor teórico analítico, corresponde siempre con 2 veces la velocidad media en la
tubería calculada ésta como caudal partido de la sección de la tubería. Ese valor coincide con
el valor de velocidad constante a la entrada del tubo, así que la velocidad máxima (Vmax) se
obtiene fácilmente como dos veces la velocidad V de partida.
Guarde esta imagen para mostrarla en el informe posterior que habrá de entregar. Para ello,
utilice la opción Hardcopy… que ya se detalló en el apartado 7.
Repetimos la comparación entre los 3 mallados, esta vez empleando el coeficiente de

fricción. Comience liberando los datos de velocidad que cargó previamente. Para ello, dentro de la
ventana Solution XY Plot, pulse el botón Free Data. Después, dibuje el coeficiente de fricción del
caso tuberia3.cas y cargue los ficheros cf.xy y cf2.xy. Finalmente, obtendrá un gráfico
como el siguiente.
El mallado más fino (línea blanca) proporciona en este caso un valor del coeficiente de
fricción de 0.399 en la zona completamente desarrollada, muy cercano al valor teórico de 0.4 que se
obtendría para este caso. El mallado más grueso, proporciona un valor menos exacto: 0.384.
29
NOTA: El valor del coeficiente de fricción, definido como en el apartado anterior, se puede
calcular directamente conocido el número de Reynolds. Así, de forma teórica, se sabe que
16
Cf = (Nótese que f = 4Cf donde f es el coeficiente de fricción típico que se emplea
Re
64
normalmente. Así, f = )
Re
Guarde esta imagen para mostrarla en el informe posterior que habrá de entregar. Para ello, utilice
la opción Hardcopy… que ya se detalló en el apartado 7.
NOTA: Finalmente, para realizar la comparativa de los perfiles de velocidad a la salida,

vamos a incluir la comparativa con el perfil teórico. Como es sabido el perfil teórico se obtiene
analíticamente como:
⎛ r2 ⎞
v = Vmax ⎜ 1 − 2 ⎟
⎝ R ⎠
Vamos a crear una función que responda a esta ecuación. Para ello, empleamos:
Main Menu > Define > Custom Field Functions …
A continuación se introduce la fórmula anterior (tecla a tecla), teniendo en cuenta que Vmax
corresponde a dos veces la V dato, que en R se introduce el radio del conducto y para introducir la
variable r hay que especificarle la coordenada y que emplea FLUENT. Para ello elegimos la
categoría Grid… en Field Functions y luego la subcategoría Y-Coordinate. Para introducirla en la
fórmula (campo Definition), hacemos clic en Select. Una vez introducida la fórmula, le damos
nombre, v y a continuación hacemos clic en Define para validar.
30
A continuación hacemos un XY Plot que muestre la variable v. Le decimos al programa que

la dibuje en outlet. Para ello, se ha de elegir dentro de categoría la opción Custom Field
Function… y en subcategoría, la única disponible (es la única que hemos creado): v. Finalmente, la
guardamos en disco haciendo clic en Write to File y le damos como nombre profileth.xy
Finalmente, hacemos la comparativa entre los casos para el perfil de velocidad a la salida,
añadiendo el perfil teórico que se acaba de generar. Proceda como en los gráficos anteriores para
llevar a cabo la comparativa, dibujando el perfil del caso cargado y luego cargue los ficheros
profile.xy, profile2.xy de los otros casos. y profileth.xy del teórico. Dibujando
y escalando, debería llegar a obtener lo siguiente:
Como se puede comprobar, la discretización con 25 nodos proporciona prácticamente el

mismo perfil que la curva teórica. Sin embargo, el mallado más basto no da tan buenos resultados.
Guarde esta imagen para mostrarla en el informe posterior que habrá de entregar. Para ello,
utilice la opción Hardcopy… que ya se detalló en el apartado 7.
9.- RESULTADOS FINALES
Cada alumno ha de presentar al final de la realización de esta práctica un informe con los
resultados que ha obtenido a partir de sus datos personalizados. En dicho informe se han de incluir
los siguientes puntos:
• Datos personalizados
• Gráfica comparando los perfiles longitudinales de las tres mallas estudiadas.
• Gráfica comparando el coeficiente de fricción de las tres mallas estudiadas.
• Gráfica comparando el perfil parabólico de la velocidad de las tres mallas estudiadas con el
perfil teórico.
IMPORTANTE: En cada gráfica hay que razonar los resultados y compararlos con lo que el
análisis teórico prevé.
Además de estos resultados, cuyo proceso de obtención ya fue explicado anteriormente, hay
que estudiar los siguientes puntos:
31
• Para el caso tuberia3.cas, que es el más exacto, hay que ejecutar de nuevo el
modelo para dos nuevos valores del número de Reynolds. En particular, un caso tendrá
un número de Reynolds que sea dos veces el valor original y el otro caso para un número
de Reynolds que sea la mitad del valor original. Para conseguir esto, basta con modificar
el valor de la viscosidad, μ, y ajustarlo para conseguir ese nuevo valor del número de
Reynolds. Nómbrelos, por ejemplo, como tuberia3-200%Re.cas y tuberia3-
50%Re.cas,ejecútelos y salvelos en su directorio de trabajo.
• Calcular el valor de la longitud de la zona de entrada hasta que se tiene flujo

completamente desarrollado utilizando la expresión analítica y compararlo con los
resultados que nos da FLUENT en los tres casos:
Le
≈ 0.06 Re ( caso laminar )
D
• Mostrar las gráficas del perfil de velocidad para los tres valores del número de Reynolds
y del coeficiente de fricción. Explique brevemente el comportamiento que se observa y
razónelo.
• Dibuje también los vectores de velocidad en la zona de entrada para los tres casos de
número de Reynolds. Siga las instrucciones que ya se detallaron en el apartado anterior.
• Finalmente, ¿Qué ocurre con el perfil de velocidad parabólico al cambiar el número de

Reynolds? ¿Se ve una tendencia similar a lo que ocurre con el perfil longitudinal de
velocidad o no?
32
10.- RESOLUCIÓN EN CASO TURBULENTO
En esta última parte se vuelve a analizar el flujo en la tubería, manteniendo los mismos
valores de velocidad V, diámetro D y longitud L del modelo, pero modificando el valor de la
viscosidad del fluido. En este caso, lo hacemos mucho menos viscoso (100 veces menos que antes),
fijando un valor de μ=2x10-5 kg/(ms). Con esto nos aseguramos que el flujo supere el valor crítico
entre 2000 y 4000 y que por tanto estemos en un caso completamente turbulento.
Un flujo turbulento contiene multitud de fluctuaciones de una escala muy pequeña en el

tiempo. Por esta razón, resolver estas fluctuaciones con un código numérico es muy costoso
(cuando no una quimera). Para hacer frente a este problema, las variables del flujo, como la
velocidad o la presión, se promedian en el tiempo con el objeto de eliminar esas fluctuaciones.
Desgraciadamente, las ecuaciones de gobierno promediadas en el tiempo no presentan un cierre
completo, esto es, introducen términos adicionales (medias de fluctuaciones) que necesitan ser
modelizados mediante el empleo de algún tipo de modelo de turbulencia. Además, tampoco existe
un modelo de turbulencia universal que tenga validez en cualquier tipo de flujo o problema.
Normalmente, en función del tipo de flujo que se tenga (flujo confinado, flujo en canal, flujo
abierto, etc…) se utilizan unos modelos u otros, que además han de ser calibrados en función de los
parámetros del fluido y del flujo en cuestión. En el ejercicio que se propone a continuación, usted
adquirirá cierto conocimiento y soltura en el uso del modelo k-épsilon. Pero en el mundo real, hay
que actuar con gran precaución: hay que estar muy seguro que el modelo que se emplea y los
parámetros que se introducen para él son correctos. Es decir, no se puede utilizar FLUENT como si
fuese una caja negra, en la cual se mete un modelo y él te da unos resultados. Hay que estar seguro
que los resultados son apropiados. El modelo k-ε resuelve dos ecuaciones adicionales para hacerse
cargo de la turbulencia: una para el transporte de la energía cinética turbulenta, k, y otra para
disipación turbulenta, ε. Estas dos ecuaciones han de resolverse de forma acoplada junto con la
ecuación de continuidad, la de momento y la de energía. Por tanto, los cálculos en caso de flujo
turbulento resultan siempre mucho más complejos y con un consumo de recursos computacionales
mucho mayores que en caso de flujo laminar.
Refinado de la malla para modelo turbulento.
Comenzamos cargando cualquiera de los casos que hemos generado con anterioridad para
caso laminar en GAMBIT. A continuación crearemos un mallado “especial”, gradual desde la pared
hacia el centro de la tubería, que llamaremos tuberiaturb.dbs.
En primer lugar hay que deshacer la malla de los lados verticales. Para ello, basta con repetir
los pasos que se detallaron en el apartado 8. Ahora bien, para remallarlos, vamos a cambiar la
distribución lineal de veces anteriores por una distribución gradual. Así, para resolver el importante
gradiente de velocidad en la capa límite de un flujo turbulento, utilizaremos un mallado mucho más
fino cerca de la pared, apilando los nodos de malla contra dicha pared.
Para cada línea vertical (inlet y outlet), especificaremos la longitud de la pared al primer
nodo de la malla con un valor de 0.001 (1 mm). Además, vamos a fijar un total de 30 nodos para
cubrir el radio de la tubería. En GAMBIT, cada línea tiene asociada un vector director, que se
muestra como una flecha roja. Hemos de fijar que la flecha apunte desde la pared (basta para ello,
seleccionar el lado y hacer clic si el sentido por defecto es el contrario). Así, la distancia de la pared
al primer nodo se convierte en la “primera distancia” (First Length) , mientras que la distancia del
último nodo al eje es la “última distancia” (Last Length). Por tanto, en Type, hemos de elegir First
Length y fijar el valor de 0.001. Recordemos:
33
Operation Toolpad > Mesh Command Button > Edge Command Button
Seleccionamos los lados verticales haciendo MAYÚSCULAS + clic. Fíjese que aparecen
las flechas rojas apuntando en una dirección (normalmente en la dirección en que se creo el lado).
Si no las ve, emplee la utilidad zoom para acercarse (CTRL + arrastrar ventana). Asegúrese de que
realmente apuntan hacia el eje…si no es así, cambie el sentido haciendo clic sobre ellas mientras
mantiene pulsada la tecla de mayúsculas. También puede cambiar el sentido si una vez que la ha
seleccionado, hace clic en Reverse. Elija la opción First Length en Type y especifique 30 nodos en
Interval Count. Luego haga clic en Apply.
Una vez fijados los nodos de cada uno de los lados, repita los pasos habituales para generar
la malla de la superficie. Si ha hecho todo correctamente, el mallado que obtenga ha de tener el
siguiente aspecto:
Finalmente, salve y exporte la malla como ya sabe. Puede emplear los nombres
tuberiaturb.dbs y tuberiaturb.msh
Al igual que en el caso laminar, vamos a generar un segundo mallado más fino. Para ello,
repita el mismo proceso de este apartado, pero emplee un mallado con 60 nodos en la dirección
radial, con una First Length de 0.5 mm. El resto de parámetros manténgalo igual. Guarde el
34
segundo mallado con otro nombre, por ejemplo: tuberiaturb2.dbs y

tuberiaturb2.msh.
Una vez esté todo definido, y haya corroborado que todos estos archivos se han creado en su
directorio de trabajo, cierre el GAMBIT y pase a FLUENT.
Activación del modelo de turbulencia en FLUENT.
Como siempre, inicie FLUENT y cargue la primera de las mallas que acaba de crear:
tuberiaturb.msh. A continuación prepare el caso exactamente igual que se mostró en el
apartado 5, comprobando que la malla no presenta errores, comprobando el tamaño, etc…. Las
únicas diferencias que debe introducir ahora son:
• Ya no estamos en flujo laminar, así que hemos de definir un modelo de turbulencia.
Emplearemos un modelo k-ε (2 eqn). Al seleccionar, la ventana se expande, dando acceso a
nuevos parámetros. En la opción de Near-wall treatment, seleccione Enhanced Wall
Treatment con el objeto de obtener un mejor resultado.
• Recuerde que hemos fijado una nueva viscosidad para el fluido. Por tanto, cuando la defina
dentro de Materials…, no olvide escribir ahora 2x10-5.
• En la condición de contorno de velocidad, aparecen nuevas variables que hemos también de
fijar. En principio, igual que otras veces, en Velocity Magnitude fije el valor V que tenga
como dato. Para dar un valor a las variables turbulentas, seleccione la opción Intensity and
Hydraulic Diameter dentro de la categoría Turbulence Specification Method. Introduzca
unos valores típicos de 1 para la Turbulence Intensity (intensidad de turbulencia del 1%) e
introduzca el valor D del diámetro de la tubería como Hydraulic Diameter. Después haga
clic en Ok para validar.
35
• En la condición de contorno de presión a la salida, seguimos manteniendo el valor 0 como

valor de gauge pressure. Aparece un nuevo campo que pide valores sobre Backflow, es
decir, flujo de retorno inverso que pudiera producirse a la salida. En principio, esto no tiene
por qué ocurrir en este caso, así que no es preciso modificar los valores que incluye por
defecto.
• Para obtener una buena resolución, no olvide cambiar la discretización a segundo orden de
todas las ecuaciones que resuelve el modelo. Para ello: Main Menu > Solve > Controls >
Solution… Una vez en ese cuadro, cambie la discretización de Momentum, Turbulence
Kinetic Energy y Turbulence Dissipation Rate (haga scroll en la barra de desplazamiento
para ver esta última) a Second Order Upwind. Luego haga clic en Ok.
Este orden de discretización se refiere a la exactitud con que se resuelven los términos
convectivos en las ecuaciones. En FLUENT, la discretización de segundo orden siempre es
apropiada para las ecuaciones turbulentas. Normalmente, este esquema de segundo orden
proporciona mejores resultados que la de primer orden, si bien una de primer orden siempre es más
robusta en la convergencia. Si el segundo orden no converge, pruebe a comenzar a resolver con un
primer orden y luego haga el cambio a segundo orden más adelante.
• Finalmente, no olvide cambiar el valor umbral de los residuos a 1e-06 también a las
variables turbulentas.
Guarde el trabajo realizado (Main Menu > File > Write > Case…) e inicialice luego como
siempre, tomando los valores a la entrada como datos de partida (Compute From… inlet).
Ejecute la solución, fijando un total de 500 iteraciones para que converja. Si en 500 iteraciones
no hay convergido, ejecute más iteraciones hasta que el modelo converja a la solución. Guarde
finalmente también los datos (Main Menu > File > Write > Data…).
Una vez que ya tenga el caso tuberiaturb.cas corrido, convergido y guardado

(compruebe siempre que el archivo ha quedado realmente almacenado en su directorio de
trabajo), proceda haciendo exactamente lo mismo para obtener tuberiaturb2.cas con la
otra malla. A continuación, pasamos a analizar los resultados.
36
Resultados preliminares.
Cualquier flujo turbulento está claramente afectado por la existencia de contornos sólidos
(paredes). El modelo de turbulencia k-ε tiene validez en aquellas zonas que se encuentran alejadas
de dichas paredes, mientras que requiere de un tratamiento especial en las zonas más cercanas a las
paredes. El modelo que se emplea para simular el flujo en la zona de pared es especialmente
sensible a la calidad (o no) de la malla junto a la pared, esto es, a la densidad de malla existente
describiendo la capa límite. Existe un parámetro adimensional que describe perfectamente la
bondad de la densidad de malla empleada, conocido como y+. Aunque no vamos a entrar en detalle
sobre este parámetro, sí vamos a dejar constancia aquí de una regla a tener en cuenta:
• O se escoge una resolución de malla cerca de la pared bastante grosera, de forma que y+>30
y se fija una ley de pared logarítmica (opción Standard Wall Functions en el Near-Wall
Treatment del panel Viscous Model).
• O bien se escoge una resolución de malla muy fina, de forma que y+<5 (cercano a 1,
normalmente) y entonces se hace un tratamiento especial cercano a la pared (opción
Enhanced Wall Treatment en el Near-Wall Treatment del panel Viscous Model).
Cargamos en FLUENT el primer caso: tuberiaturb.cas.

Para poder calcular el valor de y+ de forma correcta es preciso fijarle a FLUENT los valores
de referencia idóneos:
Main Menu > Report > Reference Values …
Seleccionamos como siempre inlet en la casilla de Compute From. Estos valores de
referencia son los que FLUENT empleará para adimensionalizar la distancia del centro de la
primera celda adyacente a la pared (recordemos que es 1 mm) para conocer y+. Haga clic en Ok.
A continuación, dibujamos la distribución de y+ en las celdas adyacentes a la pared para ver

si cumplimos los criterios anteriores:
Main Menu > Plot > XY Plot
Asegúrese de que Position on X Axis está activa en Options y que en X tenemos el valor 1 y
en Y el valor 0. Selecciona la categoría Turbulence dentro de Y Axis Function y luego
seleccione la subcategoría Wall Y Plus. Lógicamente, seleccione la surface wall y luego haga
clic en Plot.
37
Como puede ver en el gráfico, el valor de y+ se encuentra entre 2.6 y 2.8 (olvidando la zona
anómala de la entrada) por lo que efectivamente es menor que 5 y por tanto, la elección de
Enhanced Wall Treatment ha sido acertado:
MUY IMPORTANTE: El valor del y+ sólo es accesible una vez se ha resuelto el modelo.
Por tanto, es un valor que sólo se puede comprobar a posteriori (y+ depende de los valores de
velocidad en las celdas adyacentes a la pared, cálculo que se obtiene tras haber iterado la solución
hasta la convergencia).
Cargue ahora el caso tuberiaturb2.cas y observe qué obtenemos al hacer la misma

gráfica. ¿Qué valores de y+ se obtienen?
• Obtenga el perfil de velocidad a lo largo de la línea media de la tubería (centerline).

Compare los resultados de las dos mallas generadas para el caso turbulento. Deberá
presentar un gráfico similar al que se presenta a continuación (para obtenerlo, proceda de
forma análoga a como ya lo hizo en la parte laminar de esta práctica).
38
Como vemos, en este caso de ejemplo la zona de flujo completamente desarrollado, empieza
sobre los 15 metros. Compruebe qué valor obtiene en los resultados y compárelo con el que
predice la siguiente correlación empírica en caso de flujo turbulento:
Le
≈ 4, 4 Re1 6
D
Además, la velocidad en la zona desarrollada llega a ser 2.38m/s, mucho menos que los casi
4 m/s del caso laminar… ¿por qué es esto así?
• Obtenga el coeficiente de fricción. Compare los dos casos y a la vez compruebe qué valor
nos daría la siguiente correlación empírica (también puede compararlo mediante el diagrama
de Moody):
f = 0.316 Re−1 4 4000 < Re < 105 (ec. de Blausius)
NOTA: Tenga en cuenta que el Cf que vemos en FLUENT es igual a f = 4Cf
• Finalmente, muestre el perfil de velocidad a la salida. Compare ambas mallas y compárelo

cualitativamente con el que se obtuvo en el caso de flujo laminar. Observando el gradiente,
¿cuál es mayor? Y por tanto, ¿qué se puede decir acerca de la mezcla de la capa límite en los
casos laminares y turbulentos?
39
Por último, vamos a analizar el efecto de la variación del número de Reynolds en el caso
turbulento
• Para el caso tuberiaturb.cas, que presenta tan buenos resultados como el caso
más mallado, hay que ejecutar de nuevo el modelo para dos nuevos valores del número
de Reynolds. En particular, un caso tendrá un número de Reynolds que sea dos veces el
valor original y el otro caso para un número de Reynolds que sea la mitad del valor
original. Para conseguir esto, basta con modificar el valor de la viscosidad, μ, y ajustarlo
para conseguir ese nuevo valor del número de Reynolds. Nómbrelos, por ejemplo, como
tuberiaturb-200%Re.cas y tuberiaturb-50%Re.cas,ejecútelos y
salvelos en su directorio de trabajo.
• Calcular el valor de la longitud de la zona de entrada hasta que se tiene flujo
completamente desarrollado utilizando la expresión analítica y compararlo con los
resultados que nos da FLUENT en los tres régimenes de Reynolds.
• Mostrar las gráficas del perfil de velocidad para los tres valores del número de Reynolds
y del coeficiente de fricción. Explique brevemente el comportamiento que se observa y
razónelo.
• Dibuje también los vectores de velocidad en la zona de entrada para los tres casos de
número de Reynolds. Siga las instrucciones que ya se detallaron anteriormente.
• Finalmente, ¿Qué ocurre con el perfil de velocidad parabólico al cambiar el número de
Reynolds? ¿Se ve una tendencia similar a lo que ocurre con el perfil longitudinal de
velocidad o no?
40
PRÁCTICAS DE
PRÁCTICA Nº 2 - FLUENT
SIMULACIÓN DEL FLUJO ALREDEDOR

DE UN CILINDRO

5. INICIALIZACIÓN DEL CASO EN FLUENT: CASO LAMINAR.
6. RESOLUCIÓN
SIMULACIÓN DEL FLUJO ALREDEROR DE UN CILINDRO Área de Mecánica de Fluidos
El objetivo de esta práctica es ilustrar tanto la preparación como la resolución del flujo no
estacionario alrededor de un cilindro. Se pretende analizar el fenómeno de desprendimiento de
vórtices que aparece bajo determinadas condiciones de flujo.
El flujo alrededor de un cilindro es uno de los clásicos problemas de la Mecánica de Fluidos.

En principio, la geometría circular sugiere la existencia de un patrón de flujo estacionario y
claramente simétrico, lo cual se cumple en caso de valores del número de Reynolds muy bajo:
cualquier perturbación que se introduzca aguas arriba es amortiguada por las fuerzas viscosas. Sin
embargo, para números de Reynolds más grandes, las perturbaciones no pueden ser neutralizadas,
lo cual conlleva a importantes fenómenos periódicos aguas abajo del cilindro, conocidos como
“desprendimiento de vórtices” (vortex shedding).
En este ejemplo se va a mostrar cómo hacer las siguientes tareas en FLUENT:

¾ Crear la geometría y la malla en el programa de pretratamiento GAMBIT.
¾ Preparar un caso no estacionario para resolverlo.
¾ Resolver una simulación dependiente del tiempo.
¾ Determinar el coeficiente de sustentación en función del tiempo.
¾ Analizar el desprendimiento de vórtices y generar animaciones que lo describan.
Problema tipo.
Consideramos un cilindro de diámetro unitario. El dominio computacional tiene una

extensión de 11.5 veces el diámetro aguas arriba y 20 diámetros aguas abajo del cilindro. La
extensión transversal del dominio se completa con 12.5 diámetros hacia cada lado. Se fijará un
número de Reynolds, basado en el diámetro del cilindro, de 150.
Fig. 1. Flujo alrededor de un cilindro.
¡Buena suerte!
1
Para crear una malla refinada en la zona de influencia del cilindro, se va a definir un
rectángulo interior de 3 diámetros de ancho por 3 diámetros de alto. A continuación, se creará un
dominio rectangular exterior según las medidas fijadas en el enunciado para completar el modelo.
Este rectángulo exterior estará fragmentado en varios trozos para facilitar la malla progresiva que se
planea para este caso. Se van a crear en primer lugar los vértices (vertex) tanto del cilindro como de
las zonas rectangulares. A continuación se unirán los vértices correspondientes mediante líneas para
crear los lados (edges) de las superficies. Por último, crearemos las superficies (face) que
corresponderán a las distintas zonas encerradas. En el punto 3 de esta guía, se mostrará cómo mallar
adecuadamente todas las superficies.
faces -> volumes.
Ha de crear una carpeta “Cilindro” dentro del subdirectorio que sea conveniente.
Inicie GAMBIT.
debe escribir: gambit –id cilindro y luego ejecute pulsando intro. Si esto no funciona,
c:\fluent.inc\ntbin\ntx86\gambit –id cilindro
(-id) el nombre “cilindro”, le estamos indicando al GAMBIT que utilice ese nombre como prefijo
del archivo que vamos a construir a lo largo de esta sesión.
• Barra principal. Observe que el nombre cilindro aparece tras el ID en la barra.
2
3


Vamos a asumir que el flujo es BIDIMENSIONAL y que el centro del cilindro es el origen
de coordenadas (0,0). Por tanto, el dominio del problema es:
−11.5 ≤ x ≤ 20, − 12.5 ≤ y ≤ 12.5
donde x e y son las coordenadas horizontal y vertical respectivamente.
Se colocará el origen del sistema de coordenadas en el centro del cilindro. Se van a crear dos
rectángulos, uno dentro del otro, antes de definir el cilindro Las cotas de los vértices de esos
rectángulos se muestran a continuación:
4
(-11.5,12.5) (-1.5,12.5) (1.5,12.5) (20,12.5)
(-1.5,1.5) (1.5,1.5) (20,1.5)

(-11.5,1.5)
(-11.5,-1.5) (20,-1.5)
(-1.5,-1.5) (1.5,-1.5)
(-1.5,-12.5) (1.5,-12.5)
(-11.5,-12.5) (20,-12.5)
Así mismo, se van a crear 4 vértices que correspondan al cilindro (bastaría con 3, pero se va
a partir en 4 trozos para facilitar el mallado posterior). Los puntos caen sobre las diagonales del
rectángulo pequeño:
(-1.5,1.5) (1.5,1.5)
(-√2/4, √2/4) (√2/4, √2/4)
(-√2/4, -√2/4) (0,0) (√2/4, -√2/4)
(-1.5,-1.5) (1.5,-1.5)
Conviene crear también el punto (0,0), origen de coordenadas.
Creando vértices.
Vertex
5

que son creados.
Repita el proceso para crear los 20 vértices restantes (incluyendo el origen de coordenadas).

Este comando ajusta los 4 vértices más exteriores del rectángulo grande de la ventana
gráfica para que podamos ver todos los puntos en pantalla (se muestra la pantalla en la captura
siguiente):
6
Creando lados.
fin de un lado de cada rectángulo.
7
en Apply para generar el lado. Repita el proceso hasta que disponga de todos los lados que se
muestran en la figura:
Ahora vamos a crear los lados del cilindro. Para ello, haciendo clic con el botón de la
derecha sobre el icono de crear lado, descubrirá un panel desplegable. Seleccione la opción de Arc.
Especifique el centro del arco como el origen (0,0) y vaya creando los 4 arcos seleccionando
los vértices correspondientes dos a dos. Debe tener activo el método de Center y End-Points.
Especifique en cada casilla los vértices apropiados y asegúrese de tener en la opción Arc la de arco
de menos de 180º. Finalmente haga clic en Apply.
8
Finalmente cierre el dominio creando los lados que faltan. Para ello, vuelva a la opción
Straight haciendo clic con el botón derecho del ratón en el icono de Create edge. Al final, ha de
obtener en la zona del cilindro lo que se ve en la siguiente captura:
NOTA: Como puede observar, el vértice situado en el origen de coordenadas sigue de color blanco,
puesto que no pertenece a ninguna entidad de jerarquía superior. El resto de vértices sí pertenecen
ahora a algún tipo de línea y por esa razón el programa los colorea de amarillo.
Crear las caras.
Form Face
Para crear las caras que quedan encerradas por 4 lados es preciso seleccionar los 4 lados y
9
Con esto, se activa la ventana de lista de lados. Seleccionando los lados que deben formar
una cara y haciendo clic en ---> se seleccionan todos los lados. Haga clic en Close.
Finalmente haga clic en Apply para crear la cara. Repita el proceso hasta que obtenga todas
las caras necesarias: 12.
NOTA: Observe que las caras adjuntas comparten el mismo lado.
10
A continuación se detallan los pasos para crear las mallas de todas las caras. Iremos
mallando desde el cilindro hacia los extremos del modelo siguiendo un mallado regular. Se
mallarán en primer lugar los lados en diagonal que se apoyan sobre el cilindro, así como los arcos
que configuran el propio cilindro, utilizándose un patrón progresivo con el objeto de optimizar la
malla. A continuación se mallarán las caras externas, utilizando celdas estructuradas. El espaciado
en la malla se especificará a través del mallado de los lados.
Edges
de lados para seleccionar los cuatro lados diagonales del modelo. Para acceder más fácilmente a los
lados, se puede hacer un zoom ventana de la zona de interés. Para ello, manteniendo pulsada la tecla
de CTRL, se hace clic y se arrastra el ratón sobre una determinada zona (la zona que quiero ver
ampliada) y continuación se suelta el botón de CTRL. Para retornar a la vista completa del modelo,
basta con volver a clicar en: Control Toolpad > Fit to Window Button.
Una vez que ha seleccionado uno de ellos, seleccione la opción Interval Count del botón
desplegable que por defecto pone Interval Size dentro de la ventana de mallado de lados. Introduzca
el valor de 40 como número de intervalos de malla. Para que los nodos se vayan separando
gradualmente de la superficie del cilindro, elija la opción First Length en Type. Introduzca el valor
0.02. Si los nodos siguen el efecto contrario, haga clic en Reverse para conseguir el efecto de que se
vayan separando cada vez más entre si según se alejan del cilindro.
11
Haga clic en Apply. Los nodos de la malla aparecen ahora en los lados diagonales, con el
siguiente aspecto:
Para el resto de lados que configuran las 4 superficies que rodean al cilindro, introduzca un
mallado equiespaciado (el que el programa fija por defecto) con 35 nodos en cada lado. Debería
obtener algo así:
12
A continuación se coloca la malla en el resto de lados. Salvo para los cuatro lados
horizontales aguas abajo del cilindro, se seguirá siempre un mallado con 35 nodos con un First
Length de valor 0.7, de forma que los nodos se vayan “abriendo” según nos alejamos del cilindro.
Para esos 4 lados finales, utilizaremos 50 nodos con First Length de 0.95. Si ha estado todo
correcto, deberíamos obtener un mallado de líneas como se muestra en la figura:
Mallar las caras.
Faces
Haga MAYÚSCULAS + clic en la cara que quiera seleccionar y luego haga clic en Apply.
Como las superficies son todas regulares de cuatro lados, las iremos mallando como tipo Quad –
Map (esto es, como celdas rectangulares), opción seleccionado por defecto.
Finalmente, se obtiene un dominio con 16,975 celdas que debería tener el aspecto que se
13
modelo. Los tres lados de la izquierda serán la entrada de velocidad, los tres lados de la derecha
serán la salida, los superiores y los inferiores se definirán como simétricos y el cilindro será una
pared. El resto de lados se usaron de apoyo para hacer la malla, pero no tienen entidad como
condiciones de contorno.
Symmetry
Inlet
Wall Outlet
Symmetry
14
Button

Face).

Lados (edge) Nombre (name) Tipo (type)

Izquierda Inlet VELOCITY_INLET
Derecha Outlet PRESSURE_OUTLET
Arriba Top SYMMETRY
Abajo Bottom SYMMETRY
Cilindro Cylinder WALL
15

Salvar y exportar.
pueda ser interpretado por Fluent a continuación. Escriba cilindro.msh como nombre de
archivo (File Name). Seleccione Export 2D Mesh, ya que efectivamente el caso es bidimensional.
Haga clic en Accept.
Compruebe que un archivo cilindro.msh ha sido creado en su directorio de trabajo.
5.- INICIALIZACIÓN DEL CASO EN FLUENT. CASO LAMINAR
Run.
Importar la malla.
16
fichero cilindro.msh. Obviamente, ese es el fichero que fue creado con el programa anterior
rectangulares (quadrilateral cells).
aparecen las cinco condiciones de contorno que definimos con anterioridad: inlet, outlet, top,
bottom y cylinder.
mallado:

17
ventana gráfica.
Practique estas operaciones para obtener una vista detallada del cilindro tal y como se
18

Elija la opción de Unsteady bajo el cuadro Time. Se emplearán las otras opciones que por
defecto indica el programa, como el tipo segregated, implicit formulation, space-2D y absolute
velocity formulation. Haga clic en Ok.

La opción de flujo laminar está establecida por defecto. Compruebe que efectivamente sus
datos corresponden con una situación de flujo laminar. En ese caso, no necesita cambiar nada en
este menú. Haga clic en Cancel.

si estamos interesados en conocer la distribución de temperatura. En este caso, sin embargo, no
tratamos ningún aspecto de transmisión de calor. Por tanto, se deja la ecuación de la energía
desactivada y salimos haciendo clic en Cancel.

Cambie la propiedad Density a valor 150 y la de Viscosity a 1. Estos son los valores que nos
garantizan un número de Reynolds de 150. Los tomaremos como valores constantes:
19
Haga clic en Change/Create.

emplerará el valor por defecto de 1 atm (101325 Pa) como Operating Pressure. Haga clic en
Cancel para mantener el valor por defecto.

Obsérvese que en este panel aparecen las 5 condiciones de contorno que definimos en
GAMBIT como zonas en la izquierda de la ventana de boundary conditions. La zona bottom debe
estar seleccionada por defecto. Compruebe que así es, así como que el tipo seleccionado es
symmetry. Después haga clic en Set… Fíjese cómo no hay nada que especificar para este tipo de
Haga clic en Set…
20
Introduzca el valor de velocidad 1 m/s en el campo Velocity Magnitude. Haga clic en Ok.
Esta operación fija el valor de velocidad del flujo que entra al dominio desde la izquierda.
Finalmente, haga clic en wall y asegúrese que su tipo es wall. Vaya haciendo clic por todas
las pestañas y compruebe que sólo parámetros sobre la ecuación del momento son accesibles. Esto
no será así en posteriores modificaciones, así que recuerde la localización de estas opciones para
después. Haga clic en Ok.
6.- RESOLUCIÓN

Cambie en la opción de Momentum a Second Order Upwind.
Además, se seleccionará como Pressure-Velocity Coupling el método PISO. PISO permite
el uso de pasos temporales altos sin que se vea afectada la estabilidad de la solución. Por tanto, se
recomienda su empleo cuando se resuelven problemas transitorios.
Deseleccione la opción Skewness-Neighbor Coupling
Haga clic en Ok.
21

From. De esta manera, la componente axial de la velocidad se fijará al mismo valor de velocidad V
que se tiene de dato de entrada. Asimismo, la componente radial se fijará en 0 y se supondrá una
presión manométrica de valor también 0. Estos valores son tomados de la propia condición de
velocidad a la entrada.


Mantenga el valor de residuo bajo el campo Convergence Criterion al valor 1e-03 tanto
para la ecuación de continuidad (continuity) como para las componentes de la ecuación del
momento (x-momentum, y-momentum).
En el campo Plotting especifique que los muestre en la Window 1.
22
Haga clic en Ok.
Monitorización del coeficiente de sustentación del cilindro.
Mientras va iterando, vamos viendo la evolución del coeficiente de sustentación sobre el

cilindro. Se han de introducir los valores de referencia apropiados para que la adimensionalización
de fuerzas y momentos sobre el cilindro sea correcta.
Seleccionamos inlet dentro de la casilla de Compute From para indicarle al programa que
tome como valores de referencia aquellos que se tienen a la entrada. Compruebe que toma los
valores correctos (si no es así, introdúzcalos vía teclado) y haga clic en Ok.
23
A continuación active la monitorización del coeficiente de sustentación (lift) sobre el

cilindro:
Main Menu > Solve > Monitors > Force …
Seleccione la opción Lift dentro del cuadro desplegable de Coefficient. A continuación,
seleccione cylinder dentro del cuadro de Wall Zones. Haga clic en Plot y en Print dentro de
Options. Seleccione como Plot Window la número 2. Haga clic en Apply y cierre el panel.
Visualización del proceso de desprendimiento de vórtices.
FLUENT permite la visualización de contornos o gráficas de las variables mientras está

iterando. De esta forma, podemos almacenar una animación del flujo durante el transitorio de
simulación. Para ello, hay que definir una animación de la siguiente forma:
Main Menu > Solve > Animate > Define…

Indicamos que queremos realizar una secuencia: Animation Sequences se fija a 1. A
continuación, se decimos que la actualiza cada 10 pasos temporales: Every: 10 y When: Time Step.
A continuación se hace clic en el botón Define.
En el nuevo panel elegimos como Window la 0 y en Display Type seleccionamos la opción

Contours. A continuación pulsamos en Properties…
Seleccionamos en categoría la opción Velocity… y en la subcategoría la opción Velocity
Magnitude. En la opción Levels introducimos el valor 50. Y dejamos sin seleccionar la opción
Filled. Hacemos clic en Display y ajustamos el zoom para ver con detalle la zona posterior del
cilindro. Como paso final, hacemos Ok en el panel de Animation Sequence y terminamos.
24
Con esto se termina todo el proceso de inicialización y definición del modelo a resolver.
Ahora, salve el trabajo:

Escriba el nombre cilindro.cas como Case File. Haga clic en Ok. Compruebe que
cilindro.cas.

En la ventana Iterate, fije el valor de 0.2 s para el tamaño de paso temporal (Time Step
Size). Fije también en 600 el número total de pasos temporales que van a realizarse (Number of
Time Steps). Finalmente, introduzca un máximo de 25 iteraciones por paso temporal en caso de no
llegar antes a la convergencia y que actualice el valor de los residuos por pantalla cada 5
iteraciones.
El Número de Strouhal.
Cuando en el flujo alrededor de un cilindro se produce un desprendimiento de vórtices, el

número de Strouhal ronda el valor 0.2. Si se pretende capturar de forma adecuada tal
desprendimiento en la simulación, es necesario fijar al menos 20 o 25 pasos temporales por cada
ciclo en el desprendimiento.
fD
El número de Strouhal se define como: St = donde f es la frecuencia del
U
desprendimiento, D es el diámetro del cilindro y U es la velocidad del flujo. En este ejemplo, como
D = 1 y U = 1, entonces la frecuencia del desprendimiento es de 0.2 Hz. Luego el tiempo
característico del desprendimiento de los vórtices es:
1 1
T= = = 5 s.
f 0.2
25
Si queremos tener del orden de 25 pasos temporales por cada ciclo, entonces es obvio que:
T 5
Δt = = = 0.2 s.
N 25
y de ahí la elección de paso temporal que se ha realizado.
NOTA: En función del número de Reynolds, la morfología del desprendimiento y el valor del
número de Strouhal varían. En las figuras adjuntas (“Viscous Fluid Flow”, F. White, Ed. Mc-Graw
Hill), se observa el efecto del número de Reynolds sobre el desprendimiento de los vórtices de Von
Kárman (nombre habitual que reciben estas estructuras en honor al autor del primer artículo
explicando su origen y comportamiento).
26
Una vez que todos los parámetros están fijados, haga clic en Iterate. Espere un rato a que el
programa ejecute los 600 pasos temporales. Al final, habrá obtenido un resultado como el siguiente:

Introduzca cilindro.dat en el campo Data File y haga clic en Ok. Compruebe
En este apartado se mostrará cómo generar gráficos y animaciones con el desprendimiento

de vórtices con el objeto de poder presentar los resultados que nos proporciona FLUENT.
Cálculo de la frecuencia de desprendimiento
Para determinar con exactitud la frecuencia de desprendimiento de vórtices, vamos a fijarnos

en la evolución del coeficiente de sustentación. Para ello, vamos a:
Main Menu > Solve > Monitors > Surface…
En el panel de Force Monitors hacemos clic en Axes…en Options deseleccionamos el
AutoRange. Ahora, bajo Range están accesibles los máximos y mínimos tanto de x como de y.
Introduzca un valor mínimo de x y un máximo de x por la zona donde se comienza a estabilizar el
coeficiente de sustentación. En este caso, elegimos entre 100 y 120. En lo que se refiere a la y,
dejamos activo el AutoRange. Además, active las grid lines que ayudarán a estimar por dónde el
flujo se vuelve completamente desarrollado. Para ello, haga clic en Major Rules y Minor Rules en
Options. Luego haga clic en Apply.
De nuevo en el panel de Force Monitors hacemos clic en Curves… En Marker Style, en
Symbol escoja la opción @. Haga clic en Apply. Salga haciendo Close y dibuje el coeficiente en la
ventana 2 haciendo clic en Plot.
27
En esta gráfica se aprecia perfectamente el periodo de cada oscilación. Aproximadamente

hemos obtenido un valor de 6.25 segundos. Por tanto, el número de Strouhal asociado a este
D 1
desprendimiento será: St = f = f = = 0.16 . Este valor cuadra bastante bien con el gráfico
U T
logarítmico que relaciona el número de Strouhal con el número de Reynolds.
Animación del desprendimiento
Vamos a crear un archivo *.mpeg con la animación que hemos ido realizando mientras
convergía el caso. Para ello:
Main Menu > Solve > Animate > Playback…
Seleccione la opción MPEG dentro del campo Write/Record Format. Haga clic en Write y
cierre el panel. Con esta operación, se creará una animación en formato MPEG en el directorio de
trabajo que podrá ser visualizada usando, por ejemplo, el Windows Media Placer u otro reproductor
cualquiera.
28
Contornos de vorticidad
Se va a hacer un gráfico que muestre la vorticidad tras el cilindro. La vorticidad es una

variable que representa perfectamente la presencia vórtices y torbellinos.
Main Menu > Display > Contours…
Seleccione en la categoría Velocity… la subcategoría Vorticity Magnitude. Haga clic en las
opciones Filled, Node Values y Golbal Range. Deseleccione el resto de opciones. Introduzca los
valores mínimo y máximo a 0.0001 y 2 respectivamente y deje en el campo Levels el valor 50.
Haga clic en display.
La figura muestra claramente el desprendimiento de vórtices:
EJERCICIO PROPUESTO: Efecto del Número de Reynolds.
Repita el ejercicio para números de Reynolds de 100 y 200 y compare los resultados.
29
UNIVERSIDAD DE OVIEDO MEDIDA Y MODELIZACIÓN
Área de Mecánica de Fluidos DE FLUJOS
PRÁCTICA 3:
SIMULACIÓN DE UN INCENDIO EN UN
TÚNEL DE FUNICULAR
• Introducción:
Se trata de analizar la propagación del humo de un incendio en el interior de un
túnel con una pendiente fuerte.
El incendio de un vehículo supone una generación importante de humo a una
temperatura elevada. Cuando este incendio se produce en el interior de un túnel, el
humo queda confinado, tendiendo a acumularse en la parte superior y avanzando en el
túnel con velocidades de propagación que dependen de la intensidad del incendio. En
pocos minutos, una parte más o menos amplia del túnel puede quedar en condiciones de
no supervivencia.
Una fuerte inclinación del túnel provoca un efecto de chimenea, por el cual el
humo se desplaza mayoritariamente por el tramo ascendente de túnel y además se
genera una corriente ascendente de aire fresco que aviva el incendio.
En este ejercicio se busca realizar una simulación bidimensional no estacionaria
de este fenómeno. Como resultados gráficos hay que obtener:
¾ Mapas de contornos de temperatura y velocidad en distintos instantes de
tiempo, tanto en las proximidades del incendio como en toda la longitud del
túnel.
¾ Gráficas XY de la distribución de temperatura y velocidad en una línea a una
altura media de la sección del túnel en toda su longitud, superponiendo en la
misma gráfica las curvas correspondientes a distintos instantes temporales, de
forma que se aprecie su evolución.
• Características de la simulación:
¾ El flujo debe ser modelizado como no estacionario, bidimensional, turbulento y con
la ecuación de la energía y las fuerzas gravitatorias activadas.
¾ La geometría se generará como horizontal. La pendiente se simulará imponiendo la
inclinación al vector gravedad. Utilícense las dimensiones del esquema.
¾ La generación del mallado se hará poniendo nodos aproximadamente cada 0.5
metros en la zona alejada del incendio y cada 0.2 m en el vagón y en la zona del
túnel cercana. Se debe procurar que el número de celdas no sea mucho mayor de
5000.
¾ Condiciones de contorno:
-1-
9 Entrada (abertura inferior del túnel): presión total de 0 Pa (pressure inlet, hay
que definirla como pressure outlet en Gambit y cambiarla en el Fluent)
9 Salida (abertura superior del túnel): presión estática de 0 Pa (pressure outlet)
9 Techo y suelo del túnel: pared con flujo de calor nulo (wall)
9 Inferior y laterales del vagón: paredes con flujo de calor nulo
9 Superior del vagón: entrada de velocidad (velocity inlet). Esta condición es la
que simula el incendio. En esta modelización simplificada no se simula el
humo, solo se introduce aire muy caliente.
• Datos:
¾ Dimensiones del túnel y vagón:
4a5m
150 a 250 m
vector gravedad
α 1.5 a 2 m
0.5 m
10 m
¾ Magnitud del incendio: velocidad de entrada de 2 m/s (“normal to boundary”) y

temperatura de 900 ºK (equivalente a un incendio de 15 MW).
¾ Pendiente del túnel: 40 º.
¾ Características de la simulación temporal:
9 Time step: 1 s.
9 Nº of time steps: 10 pasos cada vez que se calcule. Se deben guardar los
datos cada 10 segundos. En el caso de que los cambios sean pequeños, se
puede aumentar el tiempo a 20 segundos o más.
-2-
SIMULACIÓN DEL FLUJO DE AIRE EN UNA TUBERÍA.

APLICACIÓN: CALIBRACIÓN DE UNA PLACA ORIFICIO.
OBJETO:
El objeto es la obtención del caudal de aire que circula por una tubería de diámetro D mediante una placa
orificio. Cuando una placa orificio se emplea para medir el caudal, es necesario calibrarla, es decir,
obtener su coeficiente de descarga. En placas orificio el coeficiente de descarga varía con la geometría de
la placa y con el número de Reynolds del flujo. Se tratará de obtener la variación del coeficiente de
descarga para tres valores diferentes del número de Reynolds, una vez fijada la geometría.
DOMINIO:
3D 10D
D/2
W d/2
COTAS FIJAS: D = 20 cm; W = 0.5 cm
TRABAJO A REALIZAR:
Como caso de estudio se considerará un modelo estacionario, bidimensional y axisimétrico, según el
esquema de la figura. Los datos específicos para cada alumno son el coeficiente de contracción β definido
como el cociente entre el diámetro de la placa orificio d y el diámetro de la tubería D. En primer lugar se
debe generar la geometría indicada en la figura, mediante el programa GAMBIT. A continuación se
ejecutará el caso en el FLUENT para 3 situaciones distintas: caudal fijado por la velocidad de entrada
correspondiente a tres números de Reynolds diferentes, en concreto, Re = 104, Re = 105 y Re = 106. Por
último se elaborará un informe, en el que se indiquen:
1. Datos de partida.
2. Dominio con mallado.
3. Detalles del flujo en régimen permanente (mapas de presión estática, vectores velocidad y
líneas de corriente) para: i) zona de entrada, ii) zona inmediatamente anterior a la placa, iii) zona
inmediatamente posterior a la placa.
4. Valores del coeficiente de descarga de la placa orificio para los tres números de Reynolds.
Comparación con las predicciones de las fórmulas de uso común.
PASOS A SEGUIR PARA LA SIMULACIÓN:

A) GAMBIT: se seguirá el procedimiento habitual de definir vértices, aristas y caras (o faces). Sólo se
supondrá una cara. El mallado de todas las aristas se efectuará con el parámetro INTERVAL SIZE = 0.25
cm (se pueden hacer todas a la vez). El mallado de la cara se efectuará sobre celdas cuadrangulares,
opción MAP. Como condiciones de contorno se impondrán:
1. Entrada: Velocity Inlet.
2. Salida: Pressure Outlet.
3. Tubería: Wall.
4. Placa Orificio: Wall.
5. Eje de la tubería: Axis.
B) FLUENT (2d): Exportación de la geometría del GAMBIT (*.msh) al FLUENT. A continuación se han
de fijar los parámetros característicos de la simulación, teniendo en cuenta que se trata de una simulación
estacionaria de un flujo axisimétrico. Los modelos y opciones que deben seleccionarse se detallan a
continuación:
1
1.-Grid Î Scale Î Pasar las dimensiones de la geometría a metros.

2.-Define Î Models Î Solver Î Axisymmetric.
3.-Define Î Models Î Viscous Î Reynolds Stress Model. Dejar los parámetros por defecto.
4.-Define Î Boundary Conditions:
- Velocity Inlet: introducir el valor de velocidad correspondiente.
- Pressure Outlet: introducir el valor de presión (=0).
5.-Solve Î Monitors Î Residual: activar Plot.
6.-Solve Î Monitors Î Surface: para controlar la velocidad de salida del agua durante el cálculo:
- Surface monitors = 1.
- Activar Plot.
- Pulsar: Define Î Report type: Integral; Report of: Velocity; Surface: salida.
7.-Solve Î Initialize Î Initialize. Dejar los parámetros por defecto y pulsar Init.
8.-Solve Î Iterate:
- Number of Iterations = 1000.
- Reporting Interval: 10.
- UDF Profile Update Interval: 1.
FORMULACIÓN TEÓRICA:
La medida del caudal en conducciones cerradas consiste en la determinación de la cantidad de volumen
que circula por la conducción por unidad de tiempo (en el SI, la unidad es el m3/s). Los instrumentos que
llevan a cabo la medida de un caudal se denominan, habitualmente, caudalímetros o medidores de caudal.
Un tipo habitual de medidores de caudal son los medidores de presión diferencial. Esta clase de
medidores presenta una modificación (habitualmente una reducción) de la sección de paso del fluido,
dando lugar a que el fluido aumente su velocidad, lo que origina un aumento de su energía cinética y, por
consiguiente, su presión tiende a disminuir en una proporción equivalente, creando una diferencia de
presión estática entre las secciones aguas arriba y aguas abajo del medidor.
Principio de funcionamiento: en la figura se muestra una placa de orificio colocada transversalmente en

un conducto horizontal; la placa posee un orificio centrado con el conducto de sección transversal A0;
aguas abajo de la restricción las líneas de corriente convergen y forman un área de flujo mínima AC,
denominada "contracción de la vena". Se instalan dos tomas de presión estática, una a cada lado de la
placa.
Se aplican la ecuación de continuidad y la energía (para flujo incompresible estacionario) entre los puntos
1 y el plano de la contracción de la vena:
v1 A1 = v C A C
1 1
p1 + ρ g z1 + ρ v 2 = pC + ρ g zC + ρ v C2
2 1 2
Si se combinan ambas ecuaciones y se despeja la velocidad en la contracción v C :
2
2 (p1 − pC )
vC =
(
ρ ⎛⎜1 − A C / A1 )
2⎞
⎟
⎝ ⎠
El caudal teórico se obtiene multiplicando esa velocidad por la sección:
AC 2 (p1 − pC )
QT = vC AC =
1 − ( A C / A1 )
2 ρ
El caudal real es inferior al teórico debido a que en la realidad hay pérdidas de carga. Además, lo que se
puede medir fácilmente es la presión en 2, y el área de la contracción de la vena no es conocida, por lo
que, en la expresión anterior, conviene reemplazar A C por A 0 . Estos efectos se tienen en cuenta mediante
la introducción de un coeficiente de descarga, CD. Por tanto el caudal real será:
CD A 0 2 (p1 − p2 )
Q=
1 − ( A 0 / A1 )
2 ρ
d A0
Para sección transversal circular, conviene introducir la relación de diámetros: β= = , y utilizar la
D A1
anterior expresión como:
2 (p1 − p2 )
Q = CD A 0
(
ρ 1 − β4 )
vD
El coeficiente de descarga CD depende de β y del número de Reynolds: Re = . Dicho coeficiente de
ν
descarga puede obtenerse como el cociente entre el caudal real Q y el caudal teórico QT. Para medir la
diferencia de presión creada por la placa de orificio, es necesario conectar dos tomas, una en la zona
aguas arriba y otra en la zona aguas abajo de la placa. En la figura se muestran dos disposiciones posibles
de las tomas: en las bridas y en la tubería. Estas últimas se sitúan a una distancia D aguas arriba y a una
distancia D/2 aguas abajo de la placa (siendo D el diámetro interior de la tubería).
Una expresión ajustada empíricamente para calcular el coeficiente de descarga de una placa orificio
como la considerada en este trabajo, es la siguiente:
CD = 0.5959 + 0.0312β2.1 − 0.1840β8 + 0.09β4 − 0.01685β3 + 91.71β2.5 Re−0.75

4
(1 − β )
3
PRÁCTICAS DE
SIMULACIÓN DEL FLUJO IDEAL Y

COMPRESIBLE EN UNA TOBERA

5. INICIALIZACIÓN DEL CASO EN FLUENT: FLUJO IDEAL.
6. RESOLUCIÓN
8. EJERCICIOS PROPUESTOS
SIMULACIÓN DEL FLUJO IDEAL Y COMPRESIBLE EN UNA TOBERA Área de Mecánica de Fluidos
El objetivo de esta práctica es ilustrar la puesta en marcha y resolución de un flujo

axisimétrico y compresible a través de una tobera.
El flujo a través de una tobera convergente-divergente es uno de los problemas-tipo que se

suele utilizar para modelizar flujo compresible mediante herramientas numéricas de CFD. La
predicción de la onda de choque en el flujo muestra uno de los efectos más característicos del flujo
compresible. Una descripción adecuada de este fenómeno constituye un importante reto dentro de
esta disciplina de la mecánica de fluidos. Además, con vistas a resolver los grandes gradientes de
presión que se establecen, es habitual el uso de esquemas numéricos especiales así como de mallas
muy finas, con refinados graduales en zonas muy localizadas.
En general, esta práctica nos ayudará a aprender:
• Generar la geometría de una tobera convergente-divergente.

• Simular un flujo compresible.
• Calcular resultados para flujo ideal (inviscid).
• Comparar resultados con diversos modelos.
D1 D2
Figura. Tobera convergente-divergente.
La figura muestra la sección longitudinal de una tobera convergente-divergente, simétrica

respecto del eje horizontal. La longitud de la dicha tobera (L) es 0.6 m. El radio de entrada de la
tobera (D1/2) es 0.1 m y el de la salida (D2/2) es de 0.12 m. El cociente entre el área de la garganta y
el de la entrada es 0.5625. Inicialmente, la diferencia de presión en la tobera entre la entrada y la
salida serán 0.12 MPa.
En primer lugar, se resolverá el problema utilizando un modelo de turbulencia. A

continuación haremos la resolución para flujo ideal, que utilizaremos para comparar con las
fórmulas teóricas.
¡Buena suerte!
1
Dada la simetría axial de la tobera, se va a simular un flujo axisimétrico, de modo que

únicamente será necesario recrear la mitad de la geometría de la figura anterior. Además, para
facilitar el mallado, se va a dividir el dominio en dos partes: uno antes de la garganta de la otbera y
otro después. Se crearán los 4 vértices exteriores y un spline para reproducir la geometría
convergente-divergente de la tobera. Se unirán los vértices correspondientes mediante líneas para
crear los lados (edges) y después crearemos las dos superficies (face) que corresponden a la sección
longitudinal de la tobera. En el punto 3 de esta guía, se mallarn las superficies implementadas.
faces -> volumes.
Ha de crear una carpeta “Tobera” dentro del subdirectorio que sea conveniente. Utilizaremos
esa carpeta de trabajo para almacenar los archivos y resultados que se vayan generando a lo largo de
esta sesión.
Inicie GAMBIT.
debe escribir: gambit –id tobera y luego ejecute pulsando intro. Si esto no funciona,
c:\fluent.inc\ntbin\ntx86\gambit –id tobera
(-id) el nombre “tobera”, le estamos indicando al GAMBIT que utilice ese nombre como prefijo del
• Barra principal. Observe que el nombre tobera aparece tras el ID en la barra.
2
3


Vamos a asumir que el flujo es AXISIMÉTRICO (esto es, tiene simetría circunferencial con
todas las secciones longitudinales presentando las mismas características). Por tanto, el dominio del
problema es:
0 ≤ r ≤ 0.12 , 0 ≤ x ≤ 0.6
donde r y x son las coordenadas radiales y axiales respectivamente.
Se colocará el origen del sistema de coordenadas en la esquina inferior izquierda de la

tobera. Las coordenadas de las esquinas se muestran en la figura inferior:
4
11 12 13 B
01 02 03 04 09 10
A 05 06 07 08
D C
O
En primer lugar se crearán los vértices en las cuatro esquinas, juntándose los vértices
adyacentes para crear las líneas rectas. Luego se unirán los 13 puntos del spline para crear la
sección variable de la tobera. Se romperá en dos ese lado alabeado y se crearán dos superficies
independientes (antes y después de la garganta).
Creando vértices.
Vertex

que son creados.
Repita el proceso para crear todos los vértices restantes:
5
Vértice
Coord-X Coor-Y
(vertex)
A 0 0.1
01 0.025 0.0995
02 0.05 0.09825
03 0.1 0.0935
04 0.15 0.0875
05 0.2 0.0815
06 0.25 0.0768
07 0.3 0.075
08 0.35 0.0782
09 0.4 0.0866
10 0.45 0.0976
11 0.5 0.1086
12 0.55 0.1169
13 0.575 0.1192
B 0.6 0.12
C 0.6 0
D 0.3 0

Este comando ajusta los 4 vértices más externos que se acaban de crear a la ventana gráfica
siguiente):
6
Creando lados.
fin de un lado del rectángulo.
en Apply para generar el lado. Repita el proceso hasta que disponga de los cinco lados rectos
(Straight) del dominio.
Para crear el spline, haga clic con el botón derecho del ratón sobre el icono Operation
Toolpad > Geometry Command Button > Create Edge . Se despliegan varias opciones, y
ha de seleccionar la opción NURBS. Al hacer clic, nos aparece una nueva ventana pidiéndonos los
lados que van a conformar el spline. Introduzca los 15 vértices por orden estricto desde el punto A
al punto B, pasando por los 13 puntos intermedios.
7
Finalmente, para romper el spline recién creado en dos trozos, debemos utilizar el comando
split:
Operation Toolpad > Geometry Command Button > Edge Command Button >
Split Edge. Seleccionamos el lado del spline en la casilla Edge. En el campo Split With pinchamos
y elegimos la opción Vertex. A continuación señalamos el punto 07 como y hacemos clic en Apply.
De esta forma:
Crear las caras.
Form Face
Para crear cada cara que queda encerrada por 4 lados es preciso seleccionar los 4 lados y que
éstos formen un lazo cerrado. Si no es así, el programa es incapaz de generar la superficie. Para
Con esto, se activa la ventana de lista de lados. Haciendo clic en --> se seleccionan aquellos
lados que se quiera y se seleccionan de golpe. Haga clic en Close.
Finalmente haga clic en Apply para crear las caras.
8
A continuación se detallan los pasos para crear una malla en la cara previa a la garganta con
100 divisiones en la dirección axial y 60 divisiones en la dirección radial; y otra malla en la cara
posterior de la garganta con 200 divisiones en la dirección axial y 60 también en la radial. Se
mallarán en primer lugar los lados y a continuación se mallarán las caras. El espaciado en la malla
se especifica a través del mallado de los lados.
Edges
de lados para seleccionar los dos lados verticales. Para acceder más fácilmente a los lados, se puede
hacer un zoom ventana de la zona de interés. Para ello, manteniendo pulsada la tecla de CTRL, se
hace clic y se arrastra el ratón sobre una determinada zona (la zona que quiero ver ampliada) y
continuación se suelta el botón de CTRL. Para retornar a la vista completa del modelo, basta con
volver a clicar en: Control Toolpad > Fit to Window Button.
Una vez que el eje vertical ha sido seleccionado, seleccione la opción Interval Count del
botón desplegable que por defecto pone Interval Size dentro de la ventana de mallado de lados.
Introduzca el valor de 30 como número de intervalos de malla.
Haga clic en Apply. Los nodos de la malla aparecen ahora en el lado, indicando que se han
tomado 60 divisiones intermedias.
9
A continuación se incluye una tabla donde se dan todos los parámetros necesarios para
mallar adecuadamente todos los lados del modelo. Observará que al hacer MAY+clic sobre un lado,
éste se coloreará en rojo con una flecha que indica el sentido del vector director del lado. Si se
utiliza un patrón de nodos progresivo puede ser necesario que la flecha tenga un sentido
determinado. Para cambiarlo, basta con hacer clic en la tecla Reverse.
En la tabla adjunta se muestra el número de nodos (Interval Count) que debe fijar para los
lados, el tipo de progresión de malla (Type), el parámetro que controla la progresión (Ratio) y el
sentido de la flecha:
Lado Interval Tipo (type) Ratio Sentido

(edge) Count
OA 30 Successive Ratio 1 Indiferente
A-07 50 Successive Ratio 1.02 De 07 a A
07-D 30 Successive Ratio 1 Indiferente
DO 50 Successive Ratio 1.02 De D a O
07-B 100 Successive Ratio 1.001 De 07 a B
BC 30 Successive Ratio 1 Indiferente
CD 100 Successive Ratio 1.001 De D a C
Como se podrá ver, el objetivo de la malla progresiva es apilar nodos en la zona próxima a
la garganta, donde se espera encontrar mayores gradientes.
Mallar la cara.
Faces
Haga MAYÚSCULAS + clic en la cara para seleccionar y luego haga clic en Apply. Se
emplearán Elements de tipo Quad y el campo Type se seleccionará la opción Map. Esta opción me
devuelve malla estructurada. Haga clic en Apply para mallar.
modelo. El lado izquierdo es la entrada de la tobera, el lado derecho es la salida, el superior será la
pared y el inferior un eje axisimétrico.
10
B
A Wall
Pressure
Pressure
Outlet
Inlet Axis C
O
Button

Face).

Para fijar el Type:, seleccionamos PRESSURE_INLET.
11

Izquierdo Inlet PRESSURE_INLET
Derecho Outlet PRESSURE_OUTLET
Arriba (2 líneas) Nozzle WALL
Abajo (2 líneas) Centerline AXIS

Salvar y exportar.
pueda ser interpretado por Fluent a continuación. Escriba tobera.msh como nombre de archivo
clic en Accept.
Compruebe que un archivo tobera.msh ha sido creado en su directorio de trabajo.
5.- INICIALIZACIÓN DEL CASO EN FLUENT. FLUJO IDEAL.
Run.
12
Importar la malla.

fichero tobera.msh. Obviamente, ese es el fichero que fue creado con el programa anterior
rectangulares (quadrilateral cells). Obviamente es lo esperable puesto que se usaron 30 divisiones
radiales y 150 en total en la dirección axial, así que el número total de celdas no podía ser otro que
30x150=4500.
aparecen las cuatro condiciones de contorno que definimos con anterioridad: inlet, outlet, nozzle y
centerline.
mallado:
13

ventana gráfica.
Practique estas operaciones para obtener una vista detallada de la garganta tal y como se
14

Elija la opción de Axisymmetric bajo el cuadro Space. Cambie en tipo de solver de
segregated a Coupled. Se emplearán las otras opciones que por defecto indica el programa, como el
tipo implicit formulation, steady flow y absolute velocity formulation. Haga clic en Ok.
Main Menu > Define > Models > Inviscid

La opción de flujo laminar está establecida por defecto. Cambie la opción a Inviscid (flujo
ideal). Haga clic en Ok

Cambie la propiedad Density al tipo ideal-gas. Haga clic en Change/Create. De forma
automática el programa activa la ecuación de la energía (comprobarlo en Main Menu > Define >
Models > Energy).
15

En el caso de flujo compresible, se recomienda que la presión de operación se fije en 0 para
minimizar errores debidos a fluctuaciones de presión:
A continuación se fijarán los valores de presión a la entrada y presión a la salida.

Obsérvese que en este panel aparecen las 4 cuatro condiciones de contorno que definimos en
GAMBIT como zonas en la izquierda de la ventana de boundary conditions. La zona centerline
debe estar seleccionada por defecto. Compruebe que así es, así como que el tipo seleccionado es
axis. Después haga clic en Set… Fíjese cómo no hay nada que especificar para este tipo de
tipo de esta condición de contorno es pressure-inlet. Recuerde que el tipo de condición fue ya fijada
Haga clic en Set…
16
Introduzca el valor de 220000 Pa en el campo de Gauge Total Pressure así como el valor de
210000 Pa para la Supersonic/Inicial Gauge Pressure. Introduzca el valor de 300K para la
temperatura de entrada. Haga clic en Ok. Esta operación fija las condiciones del flujo que entra al
dominio desde la izquierda.
La diferencia entre esas dos presiones (absoluta y relativa), las emplea para determinar el
número de Mach y la velocidad del flujo a la entrada según la expresión (válida para aire):
= (1 + 0.2Ma 2 )
P0 3.5
La presión (absoluta) que se tiene a la salida es presión atmosférica. Puesto que la presión de
operación se ha dejado fijada en cero, la presión manométrica a la salida (outlet gauge pressire) ha
de ser 100000 Pa aprox. Seleccione outlet bajo el campo Zone. El Type de esta condición de
contorno es pressure-outlet. Introduzca el valor de presión en Gauge Pressure y mantenga también
300K para la Backflow Total Temperature. Haga clic en Ok para validar los cambios.
Finalmente, haga clic en nozzle y asegúrese que su tipo es wall. Vaya haciendo clic por
todas las pestañas y compruebe que…¡no se puede cambiar nada!. Esto es así porque ha elegido el
modelo de tipo inviscid. Haga clic en Ok.
6.- RESOLUCIÓN

Mantenga en la opción de Flow el tipo Second Order Upwind.
17
Haga clic en Ok.

From. De esta manera, la componente axial de la velocidad se fijará en función de las condiciones
de presión de entrada. El resto de valores son tomados de la propia condición de velocidad a la
entrada.

18

Mantenga el valor de residuo bajo el campo Convergence Criterion en el valor 1e-03 tanto
momento (x-momentum, y-momentum) y de la energía (energy).
Haga clic en Ok. Con esto se termina todo el proceso de inicialización y definición del
modelo a resolver. Ahora, salve el trabajo:

Escriba el nombre tobera.cas como Case File. Haga clic en Ok. Compruebe que
tobera.cas.

En la ventana Iterate Window, cambie el número de iteraciones a 2000. Fije a 10 el número
de iteraciones que han de ejecutarse para que se refresquen en pantalla los residuos (Reporting
Interval). Haga clic en Iterate. Los residuos cada 10 iteraciones se muestran en pantalla
gráficamente y en la ventana de comandos de FLUENT:
Al cabo de las 2000 iteraciones, el caso no ha convergido. Podemos ver cómo la solución es
inestable puesto que el proceso iterativo ha sido incapaz de llegar a una solución estable:
19
Para intentar obtener la convergencia, vamos a fijar un esquema de primer orden en la

resolución de las ecuaciones del flujo. Para ello, volvemos a Solve > Controls > Solution y
elegimos la opción Fisrt Order Upwind. Hacemos clic en Ok. A continuación reanudamos el
proceso iterativo con 500 iteraciones más…
A las 250 iteraciones, el programa avisa de que se ha llegado al criterio de convergencia.
Ahora, volvemos a cambiar el esquema de discretización a segundo orden y volvemos a iterar.
Fijemos otras 500 iteraciones más…
Al cabo de 320 iteraciones, la solución converge finalmente con un esquema de segundo
orden:
En este proceso hemos aprendido que no siempre se obtiene la convergencia desde una
discretización de orden superior de forma directa. A veces hay que partir de esquemas menos
exactos para ir dando “pequeñas mejoras” al modelo numérico de forma que se garantice la
convergencia final.
Como se puede observar, un modelo acoplado de las ecuaciones bajo condición de flujo
compresible, con altos números de Mach presenta una estabilidad complicada por lo que es
necesario echar mano de los parámetros de resolución (relajación, número de Courant, orden de
discretización) para mejorar el proceso de convergencia.
20
Salve la solución a un archivo tipo *.dat (archivo de datos). Los ficheros tipo *.cas

Introduzca tobera220kPa.dat en el campo Data File y haga clic en Ok. Compruebe
Toberas convergentes-divergentes: un breve apunte.
El cálculo del flujo en una tobera convergente-divergente es de aplicación en el campo de la

aeronaútica, puesto que es frecuente su utilización en los motores de los cohetes. Este tipo de
toberas está diseñado para que partiendo de un flujo subsónico a la entrada, éste se acelere hasta
sónico en la garganta, pasando a supersónico en la parte divergente. Dependiendo de la relación de
presiones entre entrada y salida, el flujo puede ser totalmente subsónico, o puede llegar a ser sónico
en la garganta y permanecer subsónico después o pasar a supersónico. Si pasa a supersónico, puede
formarse una onda de choque en la zona divergente y salir subsónico o bien continuar supersónico
hasta la salida.
La idea que hay detrás es que en flujo subsónico, para acelerar la corriente, es necesario ir
reduciendo progresivamente la sección (utilizando un difusor). Sin embargo, en flujo supersónico es
exactamente al revés: la aceleración se obtiene en secciones divergentes (tobera). Esto es debido a
la ecuación (1D e inviscid) que rige este proceso, que en función de que el número de Mach sea
mayor o menor que 1 hace que se produzca un cambio de signo:
dV dA 1 V
= ; Ma =
V A Ma − 1
2
c
En función de la diferencia de presiones entre entrada y salida, se puede acelerar lo

suficiente el flujo como para alcanzar condiciones sónicas en la garganta (Ma=1). Si no se alcanzan,
la tobera tendrá régimen subsónico en toda su extensión. Si se alcanzan, el flujo puede presentar
onda de choque o no. El siguiente gráfico resume perfectamente todas las posibilidades, mostrando
el gráfico de presión sobre el eje de la tobera:
21
Es muy interesante observar cómo una vez que se alcanzan el régimen sónico en la garganta,
el flujo másico máximo que puede dar la tobera se bloquea y ya no aumenta aunque se siga
incrementando la diferencia de presiones entre entrada y salida (situaciones D a I).
En las situaciones D y E se forma una onda de choque en el interior del tramo divergente,
que no es más que una irreversibilidad habitual en flujo supersónicos. Estas ondas de choque son
muy delgadas (espesor del orden de la micra) y se comportan como discontinuidades en el campo
fluido. Empleando teoría de flujo isentrópico (sin fricción en las paredes y sin transmisión de calor
–inviscid-), se pueden determinar las propiedades de la onda de choque según las siguientes
expresiones:
22
Expresiones como gas perfecto para el aire

(γ=1.4):
p2 1 + 1.4 Ma12
=
p1 1 + 1.4 Ma22
V1 2.4 Ma12 ρ
= = 2
V2 0.4Ma2 + 2 ρ1
2
0.4Ma12 + 2
Ma22 =
2.8Ma12 − 0.4
Distribución del número de Mach en la tobera.
Vamos a representar la distribución del número de Mach en la tobera para identificar la

situación en la que nos encontramos (flujo subsónico en toda la tobera –casos A,B-, flujo
supersónico en la parte divergente con onda de choque –casos D,E-, o sin onda de choque –F-).
Para ello:
Main Menu > Display > Contours … Elegimos Contours of Velocity – Mach Number. En
Levels ponemos el valor 20 y activamos la opción Filled. Hacemos clic en Display.
Vemos que para ese salto de 120 kPa entre entrada y salida, el FLUENT predice la
generación de una onda de choque próxima a la sección de salida de la tobera (caso E del gráfico).
Perfil longitudinal de presiones (centerline)
Se va a hacer un gráfico que muestre la evolución de la presión (adimensionalizada por el

valor de la presión total a la entrada, P0 ) a lo largo del eje de la tobera (centerline).
Primero se crea la función de presión adimensionalizada. Para definir esa expresión en
FLUENT, utilizamos una función definida por el usuario:
23
Main Menu > Define > Custom Field Functions…Vamos seleccionando cada una de las
variables que necesitamos dentro de las casillas Field Functions, buscándolas y luego haciendo clic
en Select para introducirlas. En este caso se necesita: Pressure… Static Pressure. Para terminar le
damos el nombre p-adim y hacemos clic en Define.
NOTA: Se ha dividido por el valor de 220 kPa a la entrada.
Para dibujar el gráfico acudimos a:

quiere evaluar (centerline); todo esto en el recuadro Plot Direction. De esta manera, FLUENT
dibuja la coordenada x en la abcisa del gráfico.
Bajo Y Axis Function, seleccione Custom Field Functions… y después, dentro de la
subcategoría, elija la función recién creada, p-adim.
Es importante fijarse en que X Axis Function e Y Axis Function describen el eje x y el eje y
del gráfico respectivamente, y no las direcciones x e y de la tobera.
Finalmente, seleccione centerline bajo la categoría Surfaces puesto que efectivamente lo
que queremos es que nos de el perfil de presión a lo largo de dicha línea. Con esto se completa la
Haga clic en Plot. Se abre una nueva ventana gráfica mostrando el gráfico de tipo xy que
24
Podemos observar cómo en la zona divergente (x>0.3 m) aparece el salto de presión

asociado a la onda de choque. Para apreciar mejor el gráfico, vamos a cambiar el aspecto de la
línea. En el menú de Solution XY Plot, hacemos clic en Curves… en Line Style elegimos como
Pattern el correspondiente a una línea continua. Así mismo, en Marker Style eliminamos el
marcador en la opción Symbol. Salga haciendo clic en Apply.
Además, active las grid lines dentro del panel Axes… Para ello, haga clic en Major Rules en
Options tanto para X como para Y. Luego haga clic en Apply.

p-220kPa.xy como el XY File Name y haga clic en Ok. Compruebe que se ha creado ese
archivo en su directorio de trabajo.
25
extensión apropiada: p-220kPa.tif o p-220kPa.jpg. Como siempre, compruebe que el
archivo ha sido creado en su directorio de trabajo.
Distribución de la densidad del fluido a lo largo de la tobera.
Vamos a representar la distribución de la densidad en la tobera. Para ello:

Main Menu > Display > Contours … Elegimos Density – Density. En Levels ponemos el
valor 20 y activamos la opción Filled. Hacemos clic en Display.
El flujo está a lata presión a la entrada, por lo que la densidad es alta. A lo largo de la tobera
se expansiona, acelerándose y como consecuencia la densidad se reduce. En el bloqueo, la densidad
recupera el valor de descarga, a presión atmosférica.
Por tanto, a la vista de estos mapas, podemos caracterizar los resultados que nos da
FLUENT para las condiciones del aire antes y después del bloqueo sónico:
Ma1 = 2.45 ; Ma2 = 0.57

ρ1 = 0.35 ; ρ 2 = 1.2
p1 = 15kPa ; p2 = 100kPa
Compruebe si estos valores cuadran con las fórmulas teóricas para flujo isentrópico 1-D.
26
Vectores de velocidad
superficie. Vamos a dibujar los vectores justo en la zona de salida para ver cómo se desarrolla el
flujo tras el bloqueo sónico.

Haremos un zoom en la región cercana a la salida. La longitud y el color de las flechas
representan la magnitud de la velocidad. Cada vector es más visible si se hacen las flechas un poco
más cortas de la siguiente forma: cambiar Scale a 0.4 en el menú Vectors y en Style elija filled-
arrow. Luego haga clic en Display. Además, es posible reflejar el dibujo a partir del eje del modelo
para obtener una vista completa de la sección longitudinal. Para ello:
Main Menu > Display > Views

Bajo la opción de Mirror Planes solamente la superficie axis está accesible, puesto que es la
única condición de simetría existente en este caso. Selecciónela y haga clic en Apply. Cierre la
ventana Views. Habrá obtenido algo similar a esto:
Obsérvese que se pasa de velocidades cercanas a los 600 m/s antes de la onda de choque a
valores cercanos a los 80 m/s en la sección de salida (1/8).
Flujo másico
Finalmente, vamos a calcular el flujo másico que circula por la tobera. Para ello:
Main Menu > Report > Surface Integrals… Elegimos la opción Mass Flow Rate en la
casilla Report Type del panel que nos aparece y seleccionamos la sección de entrada: inlet, dentro
de Surfaces. Hacemos clic en compute.
En este caso, obtenemos un valor de: 9.07 kg/s.
27
8.- EJERCICIOS PROPUESTOS
Vamos a analizar otras situaciones de presión a la entrada. Para ello, cambiaremos en Main
Menu > Define > Boundary Conditions… las condiciones de presión en la condición inlet.
Se plantean ejecutar los siguientes casos:
Gauge Total Supersonic/Inital

Caso
Pressure (Pa) Gauge Pressure (Pa)
Tobera102kPa 102000 Pa 101000 Pa
Tobera103.5kPa 103500 Pa 102000 Pa
Tobera103.95kPa 103950 Pa 102000 Pa
Para cada caso, debe obtener:

• Distribución del número de Mach en la tobera (Contours of…).
• Gráfica del perfil longitudinal de presión adimensionalizada en el eje de la tobera (en
formato *.xy).
• Caudal másico en la salida de la tobera.
Gráficas conjuntas.
Una vez ejecutado el último caso, tobera500kPa.cas, dibuje el perfil de presión

adimensional a lo largo del eje centerline como ya sabe. A continuación cargue los ficheros
p-102kPa.xy , p-103.5kPa.xy, … que ha ido creando anteriormente. Para ello, una vez que
ha dibujado el perfil del caso tobera500kPa.cas, haga clic en Load File… dentro de la
ventana Solution XY Plot. Seleccione p-102kPa.xy dentro de su directorio y luego haga clic en
Ok. Repita la operación con el resto de archivos *.xy y finalmente haga clic en Plot.
28
Compare el gráfico que acaba de obtener con el ilustrado anteriormente por los libros de
texto (reproducido de “Mecánica de Fluidos”, F.M.White, Ed. McGraw-Hill, cap.9. “Flujo
compresible”).
Identifique cada uno de los casos que hemos corrido con los de la tipología A-I que se han
visto en teoría.
Observe que el punto p* p0 en la garganta de la tobera, ronda el valor 0.56, exactamente el
mismo que la relación de áreas entre la garganta y la superficie de entrada.
Para terminar, vamos a representar el gráfico del caudal másico adimensionalizado por su
máximo como función de la presión adimensional a la entrada de la tobera.
Por tanto, el caudal másico lo adimensionalizamos por su máximo, que viene dado por la
expresión:
0.6847 p0 A*
mmax =
&
( RT0 )
0.5
donde p0 es la presión total a la entrada de cada caso, A* es el área de la garganta, R es la

constante de gas ideal para el aire (287 m2/s2K) y T0 es la temperatura de estancamiento a la entrada
(300K en todos los casos). Haciendo el cociente del flujo másico obtenido por el FLUENT,
adimensionalizado en cada caso por esta expresión, y representándolo respecto de la presión a la
salida (100000 Pa) dividida por la p0 de cada caso, se obtiene:
1.1
0.9
0.8
m& m& max
0.7
0.6
0.5
0.2 0.4 0.6 0.8 1
pb p0
29
PRÁCTICAS DE
SIMULACIÓN DE LA MEDIDA DEL CAUDAL EN

FLUJO EN CANALES MEDIANTE VERTEDERO

5. INICIALIZACIÓN DEL CASO EN FLUENT: VERTEDERO ANCHO
6. RESOLUCIÓN
9. VERTEDERO ESTRECHO
PRÁCTICA No. 6
SIMULACIÓN DE LA MEDIDA DE CAUDAL SOBRE UN VERTEDERO Área de Mecánica de Fluidos
Se pretende estudiar numéricamente el flujo de agua en un canal al pasar por un vertedero.

Se considerarán las dos posibilidades: vertedero ancho y vertedero estrecho. La altura del vertedero
será h y la longitud del mismo L. Así mismo, se considerará un canal infinitamente ancho (muy
ancho en comparación con su altura), justificándose plenamente la elección de un modelo
bidimensional para su resolución.
Tras pasar por el vertedero, la lámina de agua descargará contra una sección de salida que se
supondrá a presión atmosférica. Tómese como densidad del agua, ρ=1000 kg/m3, y como
viscosidad, μ=0.001 kg/(ms).
Al ser flujo con superficie libre, se empleará un modelo adicional (VOF) que es capaz de
predecir el comportamiento de la superficie libre del agua y de tratar un problema bifásico; dos
fases coexistirán en el modelo: agua y aire. Puesto que el número de Reynolds será muy alto
(siempre del orden de 106), la fricción tiene poco peso en el análisis y gana total relevancia el
número de Froude:
V
Fr =
gy
donde V es la velocidad del flujo, g es la gravedad e y es la profundidad del agua en el canal.
El ejercicio a resolver consiste en la realización de una simulación del comportamiento de este flujo
mediante el programa FLUENT. En concreto, han de mostrarse los siguientes resultados:
• Mostrar gráficamente el régimen estacionario del flujo al paso por el vertedero.
• Obtener el perfil de velocidad a lo largo del vertedero.
• Determinar la altura de la lámina de agua sobre el vertedero. Compararla con la que predice
la formulación teórica.
• Diagrama longitudinal del número de Froude sobre el vertedero.
NOTA: Han de validarse correctamente los resultados. El flujo es turbulento, pero el número de
Reynolds no tiene relevancia frente al Froude.
Análisis preliminar.
Se trata de ver cómo se comporta el flujo de agua en un canal abierto al paso por un
vertedero. Precisamente, estos vertederos se emplean para poder estimar el caudal de agua que pasa
por el canal, en función de la altura de la superficie libre del agua justo antes de pasar por la
1
PRÁCTICA No. 6
garganta del vertedero. Utilizando una formulación unidimensional basada en la ecuación de

Bernouilli, es relativamente fácil encontrar la relación entre dichas variables. A posteriori, se
introducen una serie de correcciones empíricas para tener en cuenta efectos como la contracción de
vena, la fricción o el efecto de paredes laterales en caso de vertederos no muy anchos.
La práctica se centrará en obtener una solución numérica del flujo que habrá de ser validada
con la formulación empírica que seguro habrá sido vista en las clases teóricas. Por tanto, no es mala
idea revisar los apuntes de teoría para validar los resultados que obtengamos en la simulación.
A continuación se explican detalladamente los pasos para crear la geometría y hacer la malla
en GAMBIT (el programa de generación de casos de FLUENT) para a continuación leer el caso en
el solver FLUENT y ejecutar la solución del ejercicio.
¡Buena suerte!
Para crear el dominio que simulará la sección longitudinal del canal, vamos a crear en
primer lugar los vértices (vertex) imprescindibles de nuestra sección. A continuación se unirán los
vértices correspondientes mediante líneas para crear los lados (edges) del dominio. Por último,
crearemos la superficie (face) que corresponde al área encerrada por todos los lados. En el punto 3
de esta guía, se mallará la superficie del canal.
faces -> volumes.
Ha de crear una carpeta “Canal” dentro del subdirectorio que sea conveniente. Utilizaremos
esta sesión.
Inicie GAMBIT.
debe escribir: gambit –id canal y luego ejecute pulsando intro. Si esto no funciona, escriba
en la ventana el path completo:
c:\fluent.inc\ntbin\ntx86\gambit –id canal
(-id) el nombre “canal”, le estamos indicando al GAMBIT que utilice ese nombre como prefijo del
• Barra principal. Observe que el nombre canal aparece tras el ID en la barra.
2
PRÁCTICA No. 6
3
PRÁCTICA No. 6


Vamos a asumir que el flujo es BIDIMENSIONAL (esto es, se supone un canal lo

suficientemente ancho -ancho semiinfinito- para poder adoptar esta restricción). Por tanto, el
dominio del problema es:
4
PRÁCTICA No. 6
0 ≤ x ≤ LT , 0 ≤ y ≤ H T
donde x e y son las coordenadas longitudinal y transversal del dominio. Se adoptará una
longitud máxima del dominio, LT = 25 m, así como una altura máxima, HT = 10 m, lo
suficientemente grandes como para que las condiciones de contorno que se fijen no puedan prefijar
la solución.
Se colocará el origen del sistema de coordenadas en la esquina inferior izquierda del

rectángulo. Las coordenadas de los puntos que es necesario crear se incluyen en el dibujo:
(0,10) (25,10)
(12,h) (12+L,h)
(0,h)
(0,0) (12,0) (12+L,0) (25,0)
En primer lugar se crearán los nueve vértices de la región, juntándose los vértices
adyacentes para crear las líneas del dominio. De esta forma se creará el área subtendida por la
sección representativa del canal.
Creando vértices.
Vertex
5
PRÁCTICA No. 6

que son creados.
Repita el proceso para crear los 8 vértices restantes:
Vertex.2: (12,0,0)
Vertex.3: (12,h,0)
Vertex.4: (12+L,h,0)
Vertex.5: (12+L,0,0)
Vertex.6: (25,0,0)
Vertex.7: (25,10,0)
Vertex.8: (0,10,0)
Vertex.9: (0,h,0)

OJO: NO OLVIDE CREAR EL ÚLTIMO VÉRTICE, pues será necesario para definir la
zona de entrada al dominio.
Este comando ajusta los 4 vértices del rectángulo que se acaba de crear a la ventana gráfica
siguiente):
6
PRÁCTICA No. 6
Creando lados.
fin de un lado del dominio.
en Apply para generar el lado. Repita el proceso ocho veces más hasta que disponga todos los lados
del dominio.
7
PRÁCTICA No. 6
NOTA: Fíjese que el lado vertical de la izquierda ha de tener dos lados (edge.8 y edge.9).
Como la creación de los lados se ha empezado desde el vértice (0,0), yendo en sentido contrario a
las agujas del reloj, los dos últimos son los lados 8 y 9. Si usted ha seguido otro orden, tendrá
nombres distintos para esos lados.
Crear la cara.
Form Face
Para crear la cara que queda encerrada por los lados es preciso seleccionar los 9 lados y que
Con esto, se activa la ventana de lista de lados. Haciendo clic en All-> se seleccionan todos
los lados de golpe. Haga clic en Close.
Finalmente haga clic en Apply para crear la cara.
8
PRÁCTICA No. 6
A continuación se detallan los pasos para crear una malla en la cara rectangular con una
distancia constante de 0.1 metros tanto en la dirección longitudinal como en la transversal para
todas las celdas. Puesto que la malla es completamente regular a lo largo de todo el dominio y la
geometría del dominio es claramente cuadrangular, es posible mallar directamente la superficie sin
tener que mallar cada lado individualemente. Además, cualquier longitud del dominio es un
múltiplo entero del ancho y alto de las celdas (0.1 metros), por lo que el programa va a ser capaz de
crear un número entero de celdas.
Mallar la cara.
Faces
Haga MAYÚSCULAS + clic en la cara para seleccionar. Seleccione ahora la opción
Interval Size. Introduzca el valor 0.1 como tamaño constante de celda. Luego haga clic en Apply.
modelo. El lado izquierdo contiene la entrada del canal, el lado inferior de la derecha es la salida, el
lado vertical que finaliza el vertedero se considerará como pared. El resto de lados también los
fijaremos como pared (por conveniencia).
Wall
Inlet Wall Wall

Outlet
Button
9
PRÁCTICA No. 6

Face).
Ahora seleccionamos el “edge_9” como siempre, haciendo MAYÚSCULAS + clic sobre el

lado. El lado seleccionado aparece en el recuadro amarillo que está bajo el botón de Entity.
Hacemos clic en Apply. Con eso, la condición de contorno queda cargada dentro del
recuadro Name/Type, tal como se ve en la captura siguiente:

Edge_9 Inlet VELOCITY_INLET
Edge_5 Outlet PRESSURE_OUTLET
Edge_1, 2, 3 Solera WALL
Edge_4 Vertical WALL
Edge_6, 7, 8 Pared WALL
NOTA: Obsérvese que es posible asignar a varios lados una misma condición de contorno
(agrupando lados por una misma condición). De esta forma, los lados 1, 2, 3, por ejemplo,
mantienen una misma condición de contorno.
10
PRÁCTICA No. 6
Si se ha hecho todo correctamente, la lista final de condiciones de contorno debería figurar de la

siguiente manera en el programa:
Salvar y exportar.
pueda ser interpretado por Fluent a continuación. Escriba canal.msh como nombre de archivo
clic en Accept.
Compruebe que un archivo canal.msh ha sido creado en su directorio de trabajo.
5.- INICIALIZACIÓN DEL CASO EN FLUENT. VERTEDERO ANCHO
Run.
Importar la malla.

fichero canal.msh. Obviamente, ese es el fichero que fue creado con el programa anterior
11
PRÁCTICA No. 6
Compruebe el número de nodos, caras (de los distintos tipos) y celdas. Dependiendo de la
altura h que tenga como dato, tendrá mayor o menor número de celdas (más o menos dominio que
habrá sido mallado en GAMBIT). En el caso de este ejemplo, se obtienen 23400 celdas.
aparecen las cinco condiciones de contorno que definimos con anterioridad: inlet, outlet, solera,
vertical y pared.
mallado:

12
PRÁCTICA No. 6

ventana gráfica.
Practique estas operaciones para obtener una vista detallada del vertedero tal y como se
que el botón deselecciona todas a la vez. Cierre la ventana gráfica para continuar.
13
PRÁCTICA No. 6

Elija la opción de 2D bajo el cuadro Space por defecto. También se emplearán las otras
opciones que por defecto indica el programa, como el tipo segregated, implicit formulation, steady
flow y absolute velocity formulation. Haga clic en Cancel.
Definición de las propiedades de los fluidos y las fases.

Necesitamos cargar el material “air” (aire) y el material “water-liquid” (agua). El aire
siempre está cargado por defecto. Para añadir el agua, hacemos clic en FLUENT Database… En la
ventana desplegable que obtenemos, buscamos el material water-liquid, lo seleccionamos y damos a
Copy.
14
PRÁCTICA No. 6
En la ventana de Materials comprobamos que se ha cargado y que tenemos tanto el aire

como el agua accesibles, pulsando la persiana desplegable en Fluent Fluid Materials. Si está todo
correcto, salga haciendo clic en Close.
A continuación, vamos a definir las dos fases a las cuales les asociaremos esos dos
materiales. En primer lugar es preciso activar el modelo multiespecie VOF que nos permitirá
calcular la interfaz agua-aire en nuestro modelo. Para ello:
Main Menu > Define > Models > Multiphase
Seleccionamos la segunda opción, Volume of Fluid. Al hacer clic, se expande la ventana
dejando accesibles nuevos parámetros. En VOF Scheme, seleccionamos la primera opción,
Implicit. Al hacer esta elección, garantizamos una simulación estacionaria. Finalmente, activamos
la opción Implicit Body Force dentro del aparatado Body Force Formulation. Por supuesto,
dejamos por defecto el número de fases en 2.
Main Menu > Define > Phases

A continuación vamos a definir los fluidos que corresponden a cada fase. Al abrir esta
ventana, nos aparecen dos fases: phase-1 y phase-2. Seleccionamos la primera y damos a Set… La
nueva ventana deplegable nos permite elegir el fluido de esa fase y darle nombre. Seleccionamos
water-liquid y le damos como nombre agua. Hacemos clic en Ok y a continuación repetimos el
proceso para el aire.
15
PRÁCTICA No. 6

emplerará el valor por defecto de 1 atm (101325 Pa) como Operating Pressure. Además, esa
presión de referencia es válida para el aire, así que habrá que asociarla a un punto del dominio
donde siempre vaya a haber aire; por ejemplo, la posición (25,10). Introduzca en X(m) el valor 25 y
en Y(m) el valor 10.
Puesto que es un flujo con superficie libre, el número de Froude cobra gran relevancia. Es
necesario por tanto definir la gravedad. Hacemos clic en Gravity y fijamos el valor de -9.81 en la
casilla de Y dentro de Gravitational Acceleration .
Definición del modelo de turbulencia.
Main Menu > Define > Models > Viscous…

El caso que estamos planteando transcurre bajo condiciones de flujo turbulento.
Emplearemos un modelo k-ε (2 eqn). Al seleccionarlo, la ventana se expande, dando acceso a
nuevos parámetros. En la opción de Near-Wall Treatment, deje seleccionada la opción Stándar
Wall Functions. Haga clic en Ok.
16
PRÁCTICA No. 6

Obsérvese que en este panel aparecen las cinco condiciones de contorno que definimos en
GAMBIT como zonas en la izquierda de la ventana de boundary conditions.
en GAMBIT. Además, es posible elegir la Phase, donde nos aparecen 3 opciones: mixture, agua y
aire. En mixture, es donde damos el valor numérico de la condición de contorno. Seleccionamos
mixture y a continuación hacemos clic en Set… Nos encontramos con una ventana donde fijamos el
Q
valor de la velocidad del flujo a la entrada. Ese valor lo determinamos como: V = . En nuestro
bh
caso, como el flujo es bidimensional, usamos un caudal por unidad de ancho (q), con lo que
directamente, la velocidad será el caudal dividido por la altura de la boca de entrada (datos
particulares de cada alumno). En el ejemplo que estamos mostrando, hemos obtenido un valor de
1.5 m/s.
Para dar un valor a las variables turbulentas, dejamos la opción K and Epsilon por defecto
dentro de la categoría Turbulence Specification Method. Mantenga los valores típicos de 1 para la
Turbulence Kinetic Energy (energía cinética turbulenta) y para la Turbulence Dissipation Rate
(tasa de disipación turbulenta). Haga clic en Ok para validar.
Para terminar es preciso indicarle al modelo que por esa entrada SÓLO ENTRARÁ AGUA,
así que en la ventana anterior, en la opción de Phase, elija aire y luego haga clic en Set…Hay que
decirle que no entre nada de aire, así que en Volume Fraction (porcentaje en volumen de entrada
de aire) dejamos 0 por defecto (luego el 100% de lo que entre será aire). Damos Ok y terminamos.
17
PRÁCTICA No. 6
Para la condición de contorno de presión a la salida, hay que indicarle que es presión
atmosférica para el aire. Puesto que la presión de operación se ha dejado fijada en 1 atmósfera, la
presión manométrica a la salida (outlet gauge pressire) ha de ser cero. Seleccione outlet bajo el
campo Zone. El Type de esta condición de contorno es pressure-outlet. Con la opción de mixture
dentro del campo Phase, haga clic en Set… para fijar el valor. Por defecto ya está en 0, así que lo
dejamos y damos a Cancel para salir. Es necesario decirle que por ahí no puede entrar agua (sí salir,
pero no entrar), así que seleccionamos dentro de Phase la opción aire y damos a Set. Acto seguido,
le especificamos en la nueva ventana que Backflow Volume Fraction para el aire es 1 (luego será
cero para el agua, la cual no podrá retornar al dominio).
Finalmente, el resto de condiciones tipo Wall se dejan por defecto.

6.- RESOLUCIÓN

Cambie las opciones de Momentum, Volume Fraction, Turbulence Kinetic Energy y
Turbulence Dissipation Rate a Second Order Upwind. Además, fije en la discretización de
Pressure la opción Body Force Weighted.
Haga clic en Ok.
18
PRÁCTICA No. 6
iterativo de resolución. En este caso, vamos a fijar una zona ya inundada con agua para acelerar la
resolución del problema. Dicha zona comprende desde la entrada hasta la cara anterior del
vertedero, cubriendo la altura h del mismo.

En el menu de Solution Initialization que aparece, dejamos los valores que nos muestra por
defecto. De esta manera, las componentes de la velocidad se fijarán al valor 0. Asimismo, la presión
relativa se fijará en 0 y se supondrán unos valores unitarios de turbulencia.
MUY IMPORTANTE: Moviendo la barra de desplazamiento hacia abajo, queda accesible la
opción aire Volume Fraction. En ella especifique el valor 1 (en vez del valor 0 que pone por
defecto) para llenar de aire todo el dominio en la inicialización.
Ahora hay que “inundar” de agua la zona previa al vertedero. Para ello, definimos esa zona
en FLUENT y a continuación “parcheamos” la inicialización anterior. Para definir esa zona,
hacemos:
Main Menu > Adapt > Region…
Le especificamos a FLUENT la extensión del dominio que vamos a querer “inundar”. En
XMax(m) introducimos el valor 12 y en YMax(m) se introduce el valor h de dato que tenga cada
alumno.
MUY IMPORTANTE: A continuación se hace clic en Mark para “marcar” (=definir) esa
zona. Si por error pulsa en Adapt, FLUENT refinará la malla de esa zona, incrementando el coste de
cálculo del proceso iterativo sin que ello sea necesario. Antes de refinar, le preguntará Ok to Adapt
grid? … a lo cual debe responder no.
Utilice únicamente la opción Mark.
Finalice haciendo clic en Close.
19
PRÁCTICA No. 6
Main Menu > Solve > Initialize > Patch

Para parchear, seleccionamos la region que acabamos de “marcar” y le decimos que la llene
de agua. Para ello, en Registers to Patch, elegimos el hexahedron-r0 (es el nombre automático que
le ha dado a la región que marcamos en la ventana anterior) y en Phase elegimos la opción aire.
Seleccionamos la Variable Volume Fraction y fijamos el valor 0 para el aire. Finalmente,
ejecutamos haciendo clic en Patch.
Para comprobar que efectivamente hemos parcheado bien el modelo, vamos a hacer una
representación de las fases de agua y aire en el modelo. En principio, gracias a las operaciones que
se han realizado, el modelo debería estar repleto de agua (fracción de volumen de agua=1 (100%))
en la zona antes del vertedero y el resto lleno de aire (fracción de volumen de agua=0 (0%)). Vamos
a comprobarlo:
Main Menu > Display > Contours… En la categoría de Contours of elija la opción Phases
(normalmente es la última de todas). En la subcategoría, sólo admite una opción: Volume Fraction.
Después en el campo Phase elija agua (En principio da igual elegir agua o aire, pues obviamente
siempre la fase de agua es igual a uno menos la del aire y la del aire es uno menos la del agua -
complementarias). Haga clic en Filled y luego en Display. Si todo está correcto, deberá obtener una
mapa como éste:
20
PRÁCTICA No. 6
Para salir, cerramos la ventana gráfica y pulsamos Close en el menú Display.

valor de residuo asociado a cada ecuación de gobierno del flujo. El residuo es una medida de lo bien
que la solución actual cumple las ecuaciones de gobierno discretizadas. En nuestro caso, vamos a
iterar la solución hasta que los valores del residuo se sitúen por debajo del valor umbral de 10-6.

Cambie el valor de residuo bajo el campo Convergence Criterion al valor 1e-06 para la
ecuación de continuidad (continuity). El resto puede dejarlas en el valor 1e-03 por defecto.
21
PRÁCTICA No. 6
Monitorización
Mientras vaya iterando, el algoritmo VOF va a ir calculando la interfaz agua-aire (esto es, la
superficie libre) en el avance del agua hacia la salida. Una vez que quede perfectamente
desarrollado el flujo sobre el vertedero, el agua caerá hacia la salida. Por tanto, el caudal másico que
se pueda computar a la salida ha de quedar igualado con el caudal másico de agua que se introduce
en el modelo (por ecuación de continuidad).
Debido a que la convergencia de este tipo de modelos no es excesivamente buena, en este
caso vamos a fijarnos en la igualdad entre el flujo de entrada y de salida para determinar que la
solución ha llegado al régimen estacionario. Para ello, es preciso saber cuál es el caudal másico de
entrada e ir monitorizando el de salida mientras FLUENT itera para ver si se va obteniendo la
solución asintótica del caso.
Para saber el caudal másico de entrada hacemos:
Main Menu > Report > Fluxes… En Options elegimos la opción Mass Flow Rate y en
Boundaries elegimos la opción inlet. Finalmente, pulsamos en Compute.
En este caso, hemos obtenido un valor de 2994.6 kg/s. Conocido este dato, cerramos la
ventana haciendo clic en Close.
Para preparar la monitorización del caudal másico a la salida, definimos una función que nos
lo calcule a la salida. Para ello:
Main Menu > Define > Custom Field Functions… Creamos una function “flujo” que sea el
producto de la componente y de la velocidad (la vamos a evaluar a la salida, donde la componente y
es la que me da idea del flujo que pasa por esa superficie) por la densidad del agua y por la fracción
de volumen del agua. Cada una de estas variables las buscamos en los menús de categoría y
subcategoría y las introducimos en la ecuación usando el botón Select. Al final, le damos nombre a
la función y damos a Define. Por último, salimos haciendo clic en Close.
22
PRÁCTICA No. 6
Para monitorizar mientras itera, vamos a:

Main Menu > Solve > Monitors > Surface… Activamos una ventana de monitorización
fijando el valor 1 en Surface Monitors. Hacemos clic en Plot y Print para ir viéndolo on-line cada
Iteration y acto seguido pinchamos en Define… para seleccionar dónde queremos monitorizar la
variable que acabamos de crear.
Elegimos como Report Type la opción Integral y hacemos un Report of de una Custom
Field Function. Elegimos la única que está activa, la que acabamos de llamar flujo. Finalmente,
elegimos la Surface outlet para que nos la compute ahí. Para terminar seleccionamos Ok.
Haga clic nuevamente en Ok. Con esto se termina todo el proceso de inicialización y
definición del modelo a resolver. Ahora, salve el trabajo:

Escriba el nombre canal.cas como Case File. Haga clic en Ok. Compruebe que
canal.cas.
23
PRÁCTICA No. 6

En la ventana Iterate Window, cambie el número de iteraciones a 500 (Number of
Iterations: 500). Fije que el programa las actualice en pantalla de 10 en 10: Reporting Interval: 10.
Haga clic en Iterate. Los residuos en cada iteración se muestran en pantalla gráficamente y en la
ventana de comandos de FLUENT:
Al cabo de 500 iteraciones, la solución aún no ha convergido… El flujo a la salida no se ha

igualado al de entrada. Si hacemos ahora un contorno de la fase de agua, veríamos algo como esto:
La lámina de agua se ha desarrollado sobre el vertedero y comienza a fluir hacia la salida,

pero el proceso iterativo aún no ha terminado. Por tanto, mandaremos al programa que ejecuta más
iteraciones. Por ejemplo otras 500….Cuando el gráfico de caudal a la salida alcance una forma
asintótica igualándose al de entrada, dibujaremos la fase de agua para ver el flujo sobre el vertedero
completamente desarrollado.
Tras 2000 iteraciones, hemos llegado a una buena solución:
24
PRÁCTICA No. 6
En la iteración 2000, a parte del valor de los residuos, observamos que la variable que
estamos monitorizando alcanza -2990.3 (muy próximo a los 2994.6). La diferencia es un error del
0.14%, perfectamente asumible, por lo que convenimos en decir que la solución está
suficientemente convergida. Cada alumno ha de hacer un razonamiento similar en su caso
particular.

Introduzca canal.dat en el campo Data File y haga clic en Ok. Compruebe nuevamente
que dicho archivo queda almacenado en su directorio de trabajo. Nuevamente, podrá recuperar la
solución del modelo siempre que quiera.
Mostrar gráficamente el régimen estacionario
Se trata de representar en un contorno con la fase de agua, cómo queda la superficie libre del
agua al paso por el vertedero cuando se ha obtenido flujo desarrollado sobre él. Para ello,
simplemente:
Main Menu > Display > Contors …

Al igual que antes, en la categoría de Contours of elija la opción Phases (normalmente es la
última de todas). En la subcategoría, sólo admite una opción: Volume Fraction. Después en el
campo Phase elija agua (En principio da igual elegir agua o aire, pues obviamente siempre la fase
de agua es igual a uno menos la del aire y la del aire es uno menos la del agua -complementarias).
Haga clic en Filled y luego en Display.
25
PRÁCTICA No. 6
Utilice las utilidades del zoom para obtener una imagen precisa del vertedero:
En este tipo de gráficos, no podemos apreciar coordenadas ni alturas alcanzadas, así que
vamos a determinar la superficie libre de la lámina para luego representarla en un diagrama XY:
Main Menu > Surface > Iso-Surface…
Vamos a crear la línea que separa las fases de agua y de aire. Para ello, elegimos en Surface
of Constant, la opción Phases y en la subcategoría, Volume Fraction. En Phase elegimos agua (o
aire, da igual) y especificamos el valor 0.5, que por convenio elegimos como separación en la
interfase. La damos el nombre superficie-libre y pulsamos en Create.
Para comprobar que ha sido creado, vamos a Main Menu > Display > Grid y dibujamos el
dominio, eliminando el default-interior (la malla) y añadiendo la nueva superficie-libre. Después
hacemos clic en Display. Si todo ha estado bien ejecutado, podremos ver la interfaz que se acaba de
crear. A continuación cerramos todas las ventanas gráficas.
Para dar coordenadas a esa superficie, vamos a crear un gráfico XY…
quiere evaluar (a lo largo de la dirección longitudinal); todo esto en el recuadro Plot Direction. De
esta manera, FLUENT dibuja la coordenada x en la abcisa del gráfico.
26
PRÁCTICA No. 6
Bajo Y Axis Function, seleccione Grid… y después, dentro de la subcategoría, elija la

coordenada transversal, Y-Coordinate. Es importante fijarse en que X Axis Function e Y Axis
Function describen el eje x y el eje y del gráfico respectivamente.
Finalmente, seleccione superficie-libre bajo la categoría Surfaces puesto que efectivamente
lo que queremos es que nos de el perfil en altura de dicha curva. Con esto se completa la
Haga clic en Plot. Se abre una nueva ventana gráfica mostrando el gráfico de tipo xy que
Para apreciar mejor la zona del vertedero, vamos a cambiar el rango de los ejes del gráfico.
En el menú de Solution XY Plot, hacemos clic en Axes… en Options deseleccionamos el
AutoRange. Ahora, bajo Range están accesibles los máximos y mínimos tanto de x como de y.
Introduzca un valor mínimo de x y un máximo de x por la zona donde se comienza a estabilizar el
perfil. En el caso de la figura, elegimos entre 10 y 20. Elija usted en función de los resultados que
haya obtenido siguiendo el mismo criterio.
Además, active las grid lines que ayudarán a estimar por dónde el flujo se vuelve
completamente desarrollado. Para ello, haga clic en Major Rules y Minor Rules en Options. Luego
haga clic en Apply.
27
PRÁCTICA No. 6
Ahora repita lo mismo pero acotando el rango del eje Y. En nuestro ejemplo fijaremos el
rango de y entre 2 y 4 para el Minimum y Maximum en Range. Seleccione también Major Rules y
Minor Rules en Options. Luego haga clic en Apply.
Por último, vamos a representar la superficie libre como una línea continua, en vez de una
linea de puntos. Para ello, haga clic en Curves y seleccione el tipo de línea fina en Line Style -
Pattern para un tamaño de linea 2 (Weight) y elimine los puntos, eligiendo el símbolo vacío como
Symbol en Marker Style. Lugo haga clic en Apply.
Si observamos la gráfica con los nuevos límites, podemos determinar de forma más precisa
que la zona prácticamente horizontal en el vertedero tiene una cota de 2.82 metros (yc=0.82 m, ver
figura de portada). Además, aguas arriba, el nivel llega a los 3.58 m (H=1.55m).

superficie.xy como el XY File Name y haga clic en Ok. Compruebe que se ha creado ese
archivo en su directorio de trabajo.
28
PRÁCTICA No. 6
extensión apropiada: superficie.tif o superficie.jpg. Como siempre, compruebe que
el archivo ha sido creado en su directorio de trabajo.
Perfil de velocidad a lo largo del vertedero.
A continuación se mostrará la velocidad a lo largo del vertedero. Se mostrará el perfil justo

al principio del vertedero, a un 25% de su longitud total, al 50%, al 75% y al final en la zona de
descarga. Para ello, hay que crear esas líneas donde se quiere mostrar el perfil de velocidad,
utilizando la herramienta line/rake.
Main Menu > Surface > Line/Rake

Crearemos una línea recta que vaya de las posiciones (x0,y0)=(12,h) a (x1,y1)=(12,10). Para
ello, seleccione Line Tool en Options. Introduzca los valores apropiados para x0, y0, x1, y1.
Nombre linea0% en el campo New Surface Name. Haga clic en Create.
Para poder ver la línea que acaba de crearse, ejecute: Main Menu > Display > Grid …Una
vez esté seguro de se ha creado, cierre la ventana gráfica y el cuadro de diálogo Display.
De forma análoga se crearán el resto de las líneas. Para el 25%, defina las posiciones
(x0,y0)=(12+0.25*L,h) a (x1,y1)=(12+0.25*L,10). Para ello, seleccione Line Tool en Options.
Introduzca los valores apropiados para x0, y0, x1, y1. Nombre linea25% en el campo New
Surface Name. Haga clic en Create. Repita con el resto:
Linea50%: (x0,y0)=(12+0.5*L,h) (x1,y1)=(12+0.5*L,10)

Linea75%: (x0,y0)=(12+0.75*L,h) (x1,y1)=(12+0.75*L,10)
29
PRÁCTICA No. 6
Linea100%: (x0,y0)=(12+L,h) (x1,y1)=(12+L,10)

En definitiva, ha de haber obtenido:
superficie. Vamos a dibujar los vectores sobre las superficies que acabamos de crear para ver cómo
se desarrolla el flujo en esa zona.

Haremos un zoom en la zona del vertedero. La longitud y el color de las flechas representan
la magnitud de la velocidad.
NOTA: Puesto que el campo de presiones y de velocidad son comunes para ambas fases, para poder
dibujar los vectores de velocidad únicamente del agua, vamos a definir una función que sólo tenga
validez en la zona de agua. Para ello, vamos primero a: Main Menu > Define > Custom Field
Functions…y creamos la función velocidad-agua que al ser producto de la fracción de agua (entre
0 y 1 en la interfaz y 1 en la masa de agua) por la velocidad, anula el campo de velocidades en todas
las zonas del dominio llenas de aire. Una vez creada, damos a Define y luego salimos con Close.
Volviendo al Display de los vectores, elegimos en categoría Color By la opción Custom

Field Functions y en subcategoría, la recién creada velocidad-agua (fíjese que también tengo la
función flujo que se definió para analizar la convergencia de la solución).
30
PRÁCTICA No. 6
Cada vector es más visible si se hacen las flechas un poco más cortas de la siguiente forma:
cambiar Scale a 0.4 en el menú Vectors y en Style elija filled-arrow. Deseleccione la opción
Global Range, pero sí mantenga activa el Auto-Range. Seleccione en Surfaces todas las líneas
recién creadas (Linea0%, Linea25%, Linea50%, Linea75% y Linea100%) y luego
haga clic en Display. Si los vectores son muy pequeños, rehaga el dibujo aumentando la opción
Scale, por ejemplo a valor 10.
Los vectores de velocidad proporcionan una bonita imagen de cómo el flujo va acelerándose
al ir llegando a la descarga del vertedero.
Se puede también obtener un gráfico XY con los perfiles de la velocidad en esos planos. En
primer lugar vamos a “cortar” las líneas que creamos antes al 0, 25, 50, 75 y 100% del vertedero, de
forma que terminen donde la superficie libre de la lámina de agua. Para conseguir esto, hacemos lo
siguiente: Main Menu > Surface > Iso-clip
Con este comando, le indicamos al programa que se quede con la parte de la línea que está
comprendida (clip) entre el valor 1 de fracción de agua y el valor 0.5 (la superficie libre). Al tramo
de la linea0% que cumple esto le llamamos: perfil0%.Hacemos clic en Clip y generamos el
trozo de linea. Luego, repetimos esto para todas las líneas…
31
PRÁCTICA No. 6
Volvemos a la ventana de generación de gráficos: Main Menu > Plot > XY Plot. Y
especificamos los parámetros para que se dibuje el perfil de velocidad sobre los trocitos de línea que
acabamos de definir:
Ajustando las líneas para que sean líneas sólidas en el gráfico (editándolas en Curve) y
centrando el gráfico, obtenemos finalmente:
Comparación de la altura de la lámina en el vertedero con la formulación

teórica.
El análisis del flujo sobre el vertedero ancho se analiza de forma muy apropiada si se tiene
en cuenta que sobre él se establece un tramo que presenta flujo crítico (Fr ~ 1), unidimensional, tal
y como se representa en la figura de la portada. Aplicando Bernouilli desde aguas arriba hasta la
cresta del vertedero, se obtiene (“Mecánica de Fluidos”, cap.10, White, Ed. McGraw-Hill):
32
PRÁCTICA No. 6
V12 V2
+ Y + H ≈ c + Y + yc
2g 2g
2 H V12
Teniendo en cuenta que V = gyc , se obtiene despejando: yc ≈
c
2
+ . Y finalmente,
3 2g
suponiendo que la energía cinética aguas arriba (flujo muy lento) es despreciable, se obtiene una
pequeña relación entre la altura crítica sobre el vertedero y la altura de la lámina agua justo antes
2H
del vertedero (ver figura de portada): yc ≈ .
3
Con todo esto, es posible integrar el caudal por unidad de ancho a la salida del vertedero y
relacionarlo con estas alturas, con el fin de poder estimar el caudal a partir de la medida de dichas
alturas. Así:
32
1 ⎛2⎞ ⎛ V12 ⎞
q = gyc3 ≈ ⎜ ⎟ 2 g ⎜ H + ⎟
3⎝3⎠ ⎜ 2 g ⎟⎠
⎝
En nuestro caso, obtuvimos una altura H=1.55 m. Operando todo esto en la última ecuación:
q=3.28 m2/s, que multiplicando por la densidad del agua, nos da: 3280 kg/s… Como podemos ver
se está cometiendo un error del 8.7% entre la estimación teórica y lo que nos predice el programa.
La discrepancia tiene mucho que ver con la forma en que la superficie libre fue definida,
fijando el criterio de fase=0.5 como límite de separación. Sería mucho más recomendable introducir
un mallado más fino en esa zona para que el algoritmo recalculase la interfaz de forma más precisa.
Diagrama longitudinal del número de Froude sobre el vertedero
Los perfiles de velocidad que hemos obtenido en el último XY Plot muestran cómo la
velocidad cae de forma artificial en la interfaz debido a la definición de la función velocidad-agua.
Para evitar esto, vamos a definir una nueva superficie libre, esta vez para un valor de Volume
Fraction de agua de 0.9, yendo como antes a: Main Menu > Surfaces > Iso-Surface. Vamos a
nombrarla como superficie-libre0.9
Si hacemos un Display-Grid de las dos superficies libres que tenemos vemos que:
33
PRÁCTICA No. 6
Obviamente, la generada para el 90% queda por dentro de la del 50%. Ahora vamos a
quedarnos únicamente con la zona de la superficie del 90% que está sobre el vertedero. Para ello,
hacemos otra vez una Iso-clip, haciendo clip con la X-Coordinate en Grid entre las distancias 12 y
12+L. Al trocito lo llamamos lámina, y sobre él vamos a representar el número de Froude.
Ahora vamos a crear la función que representa al número de Froude. Para ello, vamos a
Main Menu > Define > Custom Field Functions e introducimos el Froude como:
V
Fr =
g ( y − h)
Es una definición un poco particular para tener en cuenta que en el origen de alturas en el
FLUENT no está a la cota de la cresta del vertedero, sino en el suelo (de ahí que a la Y-Coordinate
se le tenga que restar la altura h del vertedero). Además, como FLUENT ejecutará la función en
todo el dominio, para evitar que salgan valores imaginarios en la raíz (FLUENT da error), hemos
introducido el valor absoluto. Nótese que debido a la resta de h, en realidad esta función sólo es
correcta sobre el vertedero, exactamente para el trocito lámina que acaba de ser definido.
En definitiva, yendo a Main Menu > Plot > XY Plot, dibujamos sobre la superficie lámina,
la función froude de las Custom Field Functions:
34
PRÁCTICA No. 6
Efectivamente, en la zona de vertedero, el flujo se acelera, pasando de flujo subcrítico a flujo

supercrítico (ronda un valor de Fr = 1.4). Para el cálculo del caudal teórico, se supuso que el flujo
se movía sobre el vertedero bajo condiciones críticas (Fr=1). Nuevamente, este gráfico queda muy
influenciado por la superficie libre que se utilice (porcentaje de fase donde se considera que está la
interfaz).
• Datos personalizados.
• Número de celdas obtenidas e imagen del mallado alcanzado.
• Representación del estado asintótico alcanzado mediante el gráfico Caudal-Saliente VS Nº
de iteraciones.
• Representación gráfica del estado estacionario.
• Vectores velocidad y perfil de velocidad en las posiciones 0, 25, 50, 75 y 100% a lo largo
del vertedero.
• Comparación del caudal simulado con FLUENT con el caudal estimado teóricamente a
partir de las alturas alcanzadas por la lámina de agua en la simulación.
• Representación del número de Froude a lo largo del vertedero del canal.
IMPORTANTE:
Todos estos resultados, cuyo proceso de obtención ya fue explicado anteriormente, hay que
presentarlos debidamente razonados.
35
PRÁCTICA No. 6
9.- VERTEDERO ESTRECHO
En esta última parte se vuelve a analizar el flujo sobre un vertedero en un canal, pero ahora
vamos a introducir un vertedero estrecho.
Todos los datos se mantienen respecto del caso anterior, excepto la longitud L, que ahora se
va a fijar al valor 0.2 m en todos los casos.
Vertedero estrecho en GAMBIT.
Comenzamos de cero en GAMBIT, generando un nuevo caso que se llame canal2.dbs.

A continuación seguimos paso por paso lo descrito en los apartados 2 y 3, si bien teniendo en
cuenta que ahora L=0.2 m. Si todo está correcto, deberíamos llegar a obtener un dominio como éste:
Finalmente, salve y exporte la malla como ya sabe. Puede emplear los nombres
canal2.dbs y canal2.msh
Una vez esté todo definido, y haya corroborado que todos estos archivos se han creado en su
directorio de trabajo, cierre el GAMBIT y pase a FLUENT.
Ejecución del vertedero estrecho en FLUENT.
Puesto que únicamente se ha variado la geometría del vertedero, tanto el modelo como las
condiciones de contorno que se definieron para el caso del vertedero ancho son exactamente iguales
que las que se han de introducir en este nuevo caso. Por tanto, revisite los puntos 5 y 6 del presente
guión paso a paso, repitiendo el mismo procedimiento que antes.
Una vez completada esa tarea, guarde el trabajo realizado (Main Menu > File > Write >
Case&Data…) bajo el nombre canal2.cas y canal2.dat). Ejecute la solución, fijando un
total de 500 iteraciones para que converja. Si en 500 iteraciones no ha convergido, ejecute más
iteraciones hasta que el modelo converja a la solución. Guarde finalmente también los datos (Main
Menu > File > Write > Data…).
Una vez que ya tenga el caso canal2.cas corrido, convergido y guardado (compruebe
siempre que el archivo ha quedado realmente almacenado en su directorio de trabajo), pasemos al
análisis de estos nuevos resultados.
36
PRÁCTICA No. 6
Justificación del modelo de turbulencia.
Cualquier flujo turbulento está claramente afectado por la existencia de contornos sólidos
(paredes). El modelo de turbulencia k-ε tiene validez en aquellas zonas que se encuentran alejadas
de dichas paredes, mientras que requiere de un tratamiento especial en las zonas más cercanas a las
paredes. El modelo que se emplea para simular el flujo en la zona de pared es especialmente
sensible a la calidad (o no) de la malla junto a la pared, esto es, a la densidad de malla existente
describiendo la capa límite. Existe un parámetro adimensional que describe perfectamente la
bondad de la densidad de malla empleada, conocido como y+. Aunque no vamos a entrar en detalle
sobre este parámetro, sí vamos a dejar constancia aquí de una regla a tener en cuenta:
• O se escoge una resolución de malla cerca de la pared bastante grosera, de forma que y+>30
y se fija una ley de pared logarítmica (opción Standard Wall Functions en el Near-Wall
Treatment del panel Viscous Model).
• O bien se escoge una resolución de malla muy fina, de forma que y+<5 (cercano a 1,
normalmente) y entonces se hace un tratamiento especial cercano a la pared (opción
Enhanced Wall Treatment en el Near-Wall Treatment del panel Viscous Model).
En estos ejemplos con flujo bifásico siempre se ha adoptado una discretización muy grosera,
por lo que hemos adoptado desde un principio el uso de Standard Wall Functions. Vamos a
comprobar en este ejemplo con vertedero estrecho si en los contornos sólidos de la solera
tenemos efectivamente un y+>30.
Para poder calcular el valor de y+ de forma correcta es preciso fijarle a FLUENT los valores
de referencia idóneos:
Seleccionamos inlet en la casilla de Compute From. Estos valores de referencia son los que
FLUENT empleará para adimensionalizar la distancia del centro de la primera celda adyacente a la
pared (recordemos que es 0.1 m) para conocer y+. Haga clic en Ok.
A continuación, dibujamos la distribución de y+ en las celdas adyacentes a la pared para ver
si cumplimos los criterios anteriores:
Main Menu > Plot > XY Plot
Asegúrese de que Position on X Axis está activa en Options y que en X tenemos el valor 1 y
en Y el valor 0. Selecciona la categoría Turbulence dentro de Y Axis Function y luego seleccione
la subcategoría Wall Y Plus. Lógicamente, seleccione la surface solera y luego haga clic en Plot.
37
PRÁCTICA No. 6
Como puede ver en el gráfico, el valor de y+ se encuentra por encima de 1000 a lo largo de
la solera previa al obstáculo, así que la malla es muy grosera. Por tanto, la elección de Standard
Wall Fuctions queda justificada.
MUY IMPORTANTE: El valor del y+ sólo es accesible una vez se ha resuelto el modelo.
Por tanto, es un valor que sólo se puede comprobar a posteriori (y+ depende de los valores de
velocidad en las celdas adyacentes a la pared, cálculo que se obtiene tras haber iterado la solución
hasta la convergencia).
Vamos a ver ahora los resultados que hemos obtenido:
Resultados preliminares.
En el caso del ejemplo, se puede establecer que el vertedero estrecho alcanza la

convergencia después de unas 1000 iteraciones. En la figura adjunta se muestra el gráfico con los
residuos y la evolución del caudal saliente:
El caudal saliente alcanza un valor asintótico de unos -2957 kg/s. Comparado con el de
entrada (2994.6 kg/s), en este caso se introduce un mayor error que antes: 1.2%. Hacemos a
continuación un contour de la fase de agua para observar cómo queda la lámina de agua en el
vertedero:
38
PRÁCTICA No. 6
En la figura se aprecia cómo el flujo no tiene el suficiente impulso al pasar sobre el

vertedero como para separarse de la pared. Como consecuencia, el agua rebosa el vertedero y
resbala por la pared vertical hacia la salida. Al crearse una nueva capa límite en la pared vertical del
rebosadero, el cálculo no es tan preciso como en el caso del vertedero ancho.
Para evitar este fenómeno, es muy habitual “airear” el vertedero en la realidad. Para ello, por
los laterales, se suele inyectar una corriente de aire, evitando que la lámina de agua se pegue a la
pared (ver figura de la portada: ventilation). En este caso, vamos a simular la inyección de aire
cambiando la condición de contorno que hemos denominado como vertical. Entrando en Main
Menu > Define > Boundary Conditions, la seleccionamos y cambiamos su tipo wall por uno
velocity-inlet bajo la condición mixture. Introducimos un valor típico de velocidad de aire, por
ejemplo, 0.5 m/s, y validamos dando Ok. A continuación, le decimos que por ahí sólo puede entrar
aire, así que eligiendo en Phase la opción aire, tras dar a Set, imponemos un valor de Volume
Fraction de 1. Finalizamos dando Ok.
A continuación, utilizando la solución que tenemos como punto de partida, iteramos con las
nuevas condiciones para ver qué ocurre. Comenzamos damos 100 iteraciones, y una vez que se
hayan ejecutado, hacemos un Display-Contours de la fase de agua. Si todo ha sido fijado
correctamente, obtendríamos algo parecido a esto:
39
PRÁCTICA No. 6
Se puede apreciar cómo la lámina de agua ha sido obligada a separarse de la pared. A

continuación iteramos otras 900 iteraciones y vamos controlando el gráfico de caudales para ver si
tras el transitorio inicial, la solución camina hacia el estado asintótico.
La solución indica que el caudal saliente es 2994.2 kg/s…lo cual supone únicamente un
error de 0.07% respecto al caudal de entrada. Ahora, con aireación, la lámina de agua sobre el
vertedero queda:
Para completar la gráfica, vamos a determinar la cota H que alcanza el agua sobre el
vertedero aguas arriba (utilizando un gráfico XY Plot) y compararemos esto con la formulación
teórica para obtener el caudal en función de las cotas de agua (“Mecánica de Fluidos”, cap.10,
White, Ed. McGraw-Hill).
Para crear el gráfico XY Plot de la superficie libre, creamos en primer lugar con el comando
Iso-Surface dicha superficie de Volume Fraction de agua 0.9 y a continuación dibujamos un
gráfico de Y-Coordinate en XY Plot. Así:
40
PRÁCTICA No. 6
La altura de la lámina alcanza aproximadamente 3.13 m. Por tanto, H=1.13 m (ver figura de
la portada).
Por otro lado, en el caso de un vertedero estrecho, siguiendo razonamientos similares al caso
de vertedero ancho, se puede llegar a obtener una expresión del caudal por unidad de ancho en
función de la altura del agua, aguas arriba del vertedero:
⎛2⎞
q ≈ 0.81⎜ ⎟ 2 g ⋅ H 3 2
⎝3⎠
Operando esta expresión para una altura de 1.13 m, obtenemos: 2.87 m2/s… Esto, es 2870
kg/s… En este caso, la aproximación es notable: se ha cometido únicamente un error del 4%.
los siguientes puntos en el caso de vertedero estrecho:
• Datos personalizados.
• Número de celdas obtenidas e imagen del mallado alcanzado.
• Representación del contorno de la fase agua cuando no se ha aireado.
• Representación gráfica del estado estacionario una vez se ha aireado: diagrama de fase de
agua y contorno XY de la superficie libre.
• Comparación del caudal simulado con FLUENT con el caudal estimado teóricamente a
partir de las alturas alcanzadas por la lámina de agua en la simulación.
IMPORTANTE:
Todos estos resultados, cuyo proceso de obtención ya fue explicado anteriormente, hay que
presentarlos debidamente razonados.
41
PRÁCTICAS DE
SIMULACIÓN DE UN INCENDIO EN
UN TUNEL DE CARRETERA

5. INICIALIZACIÓN DEL CASO EN FLUENT
6. RESOLUCIÓN
8. EFECTO DE LA PENDIENTE Y
DEL GRADIENTE DE PRESIONES
SIMULACIÓN DE UN INCENDIO EN UN TUNEL DE CARRETERA Área de Mecánica de Fluidos
El objetivo de esta práctica es ilustrar la resolución de un flujo no estacionario multi-especie

en un dominio tridimensional. Se pretende analizar el fenómeno de dispersión y transporte del
humo generado por un incendio que tiene lugar en el interior de un túnel de carretera.
Se comenzará mostrando las estrategias más habituales empleadas para la generación de

dominios tridimensionales en CFD. A continuación se darán las instrucciones necesarias para poner
en marcha el modelo multiespecie. Será necesario emplear una simulación transitoria mientras se va
alcanzando el estado asintótico para la dispersión del humo generado en el incendio. Se analizará la
influencia de algunos parámetros en la evolución y composición de la distribución final de
temperaturas y CO2 en el túnel, tales como la diferencia de presión entre las bocas del túnel o el
efecto de la pendiente.

¾ Crear una geometría tridimensional y mallarla en el programa de pretratamiento GAMBIT.
¾ Preparar un caso multiespecie, no estacionario para su resolución.
¾ Determinar las distribuciones finales de temperatura y CO2 en el interior del túnel.
¾ Analizar los resultados bajo criterios reales de seguridad en túneles y generar animaciones
que describan el estado transitorio en el túnel.
Problema tipo.
Consideramos un túnel típico de carretera, con una longitud total de 50 metros. Se creará
una zona (centrada transversalmente en el túnel pero más cerca de una de las bocas) como
localización de origen del incendio. Debido a la simetría transversal, se simulará únicamente la
mitad del túnel, lo cual permitirá ahorrar un número importante de celdas. La potencia calorífica
que se liberar es de 5 MW, valor característico para el caso de deflagración de un turismo
convencional.
Fig. 1. Geometría del túnel.
¡Buena suerte!
1
Se comenzará creando la sección transversal del túnel. A continuación, se extruirá

longitudinalmente dicha sección para obtener el dominio tridimensional. En ese punto, habrá que
introducir la sección plana que simulará la localización del incendio. Acto seguido, será necesario
incluir un par de postizos a la salida de ambas bocas del túnel con el objetivo de no imponer
condiciones de contorno ficticias en dichas bocas. Para facilitar la creación de la geometría
tridimensional, en lugar de seguir la filosofía jerárquica de GAMBIT, procuraremos utilizar todas
las utilidades que incorpora el programa para generar volúmenes de forma rápida y sencilla.
Ha de crear una carpeta “Tunel” dentro del subdirectorio que sea conveniente. Utilizaremos
esta sesión.
Inicie GAMBIT.
debe escribir: gambit –id tunel y luego ejecute pulsando intro. Si esto no funciona, escriba
en la ventana el path completo:
c:\fluent.inc\ntbin\ntx86\gambit –id tunel
(-id) el nombre “tunel”, le estamos indicando al GAMBIT que utilice ese nombre como prefijo del
2
3


El flujo es TRIDIMENSIONAL. Vamos a colocar el origen de coordenadas (0,0) justo en

uno de los vértices de la superficie del incendio. Veremos que el dominio del problema quedará
comprendido entre:
0 ≤ x ≤ 5, 0 ≤ y ≤ 6, − 20 ≤ z ≤ 50
donde x e y definen las coordenadas transversales y z define el eje del túnel.
Se colocará el origen del sistema de coordenadas en una de las esquinas de la sección del
túnel, justo como un vértice de la superficie del incendio. Se va a crear la sección del túnel a partir
de un arco de circunferencia que representa la bóveda. El centro de dicho arco se encuentra sobre el
plano de simetría del modelo:
4
C A (0,0,0)
B (5,0,0)
C (0,6,0)
O (0,1.2,0)
O
A B
NOTA: El punto O sólo se crea por ser el centro del arco de circunferencia. Por tanto, la
línea AC será una línea única (no es necesario que el punto O divida esa línea en dos trozos).
Creando vértices.
Vertex
5

que son creados.
Repita el proceso para crear los 3 vértices restantes (puntos B,C y O).
Como puede ver, como aún estamos definiendo una sección bidimensional, la coordenada z
se deja siempre con el valor 0 por defecto.
Este comando ajusta los 3 vértices más exteriores de la sección en la ventana gráfica para
que podamos ver todos los puntos en pantalla (se muestra la pantalla en la captura siguiente):
Creando lados.
Ahora han de conectarse entre sí los tres vértices de la sección del túnel.
6
en Apply para generar el lado. Repita el proceso hasta que disponga de los lados rectos.
Ahora vamos a crear el arco de circunferencia. Para ello, haciendo clic con el botón de la
derecha sobre el icono de crear lado, descubrirá un panel desplegable. Seleccione la opción de Arc.
Especifique el centro del arco en el puntoO y seleccione como End Points los puntos B y C.
Debe tener activo el método de Center y End-Points. Especifique en cada casilla los vértices
apropiados y asegúrese de tener en la opción Arc la de arco de menos de 180º. Finalmente haga clic
en Apply.
7
Si todo se ha hecho correctamente, el resultado que obtendrá será:
NOTA: Como puede observar, el vértice situado en el punto O sigue de color blanco, puesto que no
pertenece a ninguna entidad de jerarquía superior. El resto de vértices sí pertenecen ahora a algún
tipo de línea y por esa razón el programa los colorea de amarillo.
Crear la cara.
Form Face
Para crear las caras que represente la sección del túnel es preciso seleccionar los 3 lados y
8
Con esto, se activa la ventana de lista de lados. Haciendo clic sobre el botón AllÆ , se
seleccionan todos los lados. Haga clic en Close.
Extrusión de la cara.
En primer lugar vamos a extruir la cara recién creado según la dirección del eje z. La vamos
a extruir exactamente la longitud que ocupará la zona del incendio a lo largo: 1 metro.
Operation Toolpad > Geometry Command Button > Volume Command Button >
Stitch Faces . Haciendo clic con el botón derecho del ratón, desplegamos varias opciones y
seleccionamos de entre todas ellas la de extruir: Sweep Face:
En el menú desplegable que nos aparece seleccionamos la cara que tenemos para extruir y
elegimos en el campo Path la opción vector. Por defecto nos fija el vector unitario en la dirección z.
Si quisiésemos otra dirección u otro valor de longitud a extruir, se haría clic en Define y se
especificaría ahí lo necesario. Para hacer la extrusión hacemos clic en Apply.
NOTA: Como podrá observar, las líneas son ahora verdes, lo cual significa que pertenecen a una
entidad jerárquica superior: un volumen. Para ver una perspectiva del volumen, haga clic sobre la
ventana gráfica con el botón izquierdo del ratón y muévase a izquierda o derecha con el ratón,
manteniendo pulsado dicho botón izquierdo.
9
Al hacer modelos tridimensionales, las utilidades que implementa el ratón son muy útiles
para centrar o hacer zoom. Practique las siguientes opciones:
• Manteniendo pulsado el BOTÓN IZQUIERDO del ratón y arrastrando el puntero por la

pantalla gráfica podrá hacer giros tridimensionales de las vistas.
• Manteniendo pulsado el BOTÓN CENTRAL del ratón y arrastrando el puntero por la
pantalla gráfica podrá hacer centrar la geometría en la ventana gráfica.
• Manteniendo pulsado el BOTÓN DERECHO del ratón y arrastrando el puntero por la
pantalla gráfica podrá hacer zoom dinámico de la vista.
• Manteniendo pulsada la tecla de CTRL, se hace clic y se arrastra el ratón sobre una
determinada zona (la zona que quiero ver ampliada) y continuación se suelta el botón de
CTRL, conseguirá hacer zoom ventana.
Generación de la zona de incendio.
Para marcar la superficie desde la cual se van a liberar los humos del incendio, es preciso
crear la cara del incendio y “restársela” al volumen recién creado.
Form Face . Haciendo clic con el botón derecho del ratón, desplegamos varias opciones y
seleccionamos de entre todas ellas la de extruir: Sweep Edges: . Elegimos el eje del túnel y lo
extruimos en la dirección de x positiva con una magnitud de 1 metro. (Ver capturas siguientes):
10
Hacemos clic en Apply y así conseguimos crear una cara que está “encima” del suelo del
volumen. Simplemente, nos basta con restar esa cara al volumen. Para ello, utilizamos el comando
Split:
Split Volume . Elegimos el volumen que tenemos creado y bajo la opción Split With elegimos
Face(Real) y ejecutamos haciendo clic en Apply
Al hacer esta operación, la cara inferior del volumen nos aparece partida en dos trozos.
Extrusión del resto de túnel.
Se procede a continuación a extruir cada una de las dos caras transversales hacia cada lado.
La posterior la extruiremos 10 metros en el sentido de z negativo, mientras que la cara anterior la
extruiremos 39 metros en el sentido de z positivo. Empleando el comando extrusión (Sweep Faces)
de forma similar a como se mostró antes, se obtiene finalmente:
11
Inclusión de postizos en ambas bocas.
Vamos a crear dos caras adicionales en las dos bocas del túnel recién creado. Para ello
definiremos 3 puntos en cada una de las caras adicionales que luego habrá que unir mediante líneas.
Finalmente, crearemos las dos caras a partir del lazo cerrado que conformarán los lados.
Los nuevos puntos que es necesario crear son:
D (0,11,40) ; E (10,11,40) ; F (10,0,40)
G (0,11,-10) ; H (10,11,-10) ; I (10,0,-10)
Uniendo esos puntos y creando las caras, debería obtener:
12
Para terminar, haga una extrusión de las dos caras recién creadas así como de las dos bocas
del túnel de valor 10 metros en cada uno de los sentidos. Esto es:
A continuación se detallan los pasos para crear las mallas de todos los volúmenes. Iremos
mallando desde la sección transversal en la zona de incendio y lo extruiremos (cooper)
longitudinalmente hasta llegar a los postizos. Se utilizará un patrón progresivo con el objeto de
optimizar la malla en la zona próxima al incendio. A continuación se mallarán los postizos,
empleando también técnicas de extrusión de malla.
Mallado de la zona túnel.
Edges
de lados para seleccionar los cuatro lados de la sección transversal. Para acceder más fácilmente a
los lados, se puede hacer un zoom ventana de la zona de interés. Para ello, manteniendo pulsada la
tecla de CTRL, se hace clic y se arrastra el ratón sobre una determinada zona (la zona que quiero
ver ampliada) y continuación se suelta el botón de CTRL. Para retornar a la vista completa del
modelo, basta con volver a clicar en: Control Toolpad > Fit to Window Button
En la tabla adjunta se muestra el número de nodos (Interval Count) que debe fijar para cada
uno de los lados:
13
Lado Interval Count

A-Punto 6
Punto-B 12
B-C 18
A-C 18
Respecto a los lados longitudinales, vamos a definir tres patrones diferentes, en función de
que pertenezcan al volumen del incendio (central) o a los volúmenes de 39 y 10 metros de largo
respectivamente. De esa manera:
• Lados longitudinales del volumen del incendio: equiespaciados con Interval Count=5
• Lados longitudinales del volumen de 39 m: Interval Count=39 con Sucessive Ratio=1.064
• Lados longitudinales del volumen de 10 m: Interval Count=15 con Sucessive Ratio=1.112
Procedemos a mallar los tres volúmenes. Haciendo sucesivas extrusiones de la malla de la

sección transversal, evitamos tener que mallar todas las caras de los volúmenes. Comenzamos
mallando la sección transversal:
Faces
Haga MAYÚSCULAS + clic en la sección. Se emplearán Elements de tipo Quad y el
campo Type se seleccionará la opción Pave. Esta opción me devuelve malla estructurada. Haga clic
en Apply para mallar
A continuación, vamos a mallar el volumen central:

Operation Toolpad > Mesh Command Button > Volume Command Button >
Mesh Volume
Seleccione el volumen central. Por defecto, el programa selecciona la opción más adecuada:
Cooper. En esta opción hará una extrusión de la malla desde una cara fuente (Sources) hasta su
opuesta. Como le hemos dicho que tenga 5 celdas en longitudinal, ignora el Spacing que fija por
defecto.
14
Repetimos con los dos volúmenes adyacentes de 39 y 10 metros de longitud una vez que se
ha mallado el volumen central. Finalmente, la malla del túnel debería tener el siguiente aspecto:
Mallado de los postizos.
Para mallar los postizos vamos a seguir empleando una estrategia de extrusión. Para ello, es
necesario mallar en primer lugar las dos caras adyacentes a las bocas del túnel. Se seleccionan
ambas caras y se mallan directamente, empleando un Interval Size de 1 y una malla de tipo Quad –
Pave. Esto es:
15
Finalmente se mallan los 4 volúmenes que restan empleando directamente un Cooper,

fijando como valor de Spacing: 2 metros en Interval Size.
modelo. La zona del incendio será la entrada de velocidad para los gases, el techo de los postizos
serán condiciones de salida, el plano de simetría de la geometría se definirá como tal y el resto de
superficies serán paredes.
16
Outlet
Symmetry (north-portal)
Outlet
(south-portal)
Wall
Velocity-inlet
(inlet-fire)
Button
vamos a fijar que la zona del incendio es la entrada de humos. Para ello, en Entity: seleccionamos la
opción Faces para que GAMBIT sepa que queremos seleccionar una superfice (ya por defecto trae
activo Face).
Ahora seleccionamos la superficie como siempre, haciendo MAYÚSCULAS + clic sobre

ella. El área seleccionada aparece en el recuadro amarillo que está bajo el botón de Entity. Además,
podemos darle nombre dentro del campo Name:, por ejemplo inlet-fire.
Para fijar el Type:, seleccionamos VELOCITY_INLET. Hacemos clic en Apply. Con eso,
la condición de contorno queda carga dentro del recuadro Name/Type, tal como se ve en la captura:
17
A continuación se repite el proceso con el resto de superficies que se muestran en el dibujo

anterior. Agrupándose por condiciones de contorno, si se ha hecho todo correctamente, la lista final
de condiciones de contorno debería figurar de la siguiente manera en el programa:
Salvar y exportar.
pueda ser interpretado por Fluent a continuación. Escriba tunel.msh como nombre de archivo
(File Name). Desactive la opción Export 2D Mesh, ya que efectivamente el caso es tridimensional.
Compruebe que un archivo tunel.msh ha sido creado en su directorio de trabajo.
5.- INICIALIZACIÓN DEL CASO EN FLUENT.
Run.
La opción 3ddp se utiliza para seleccionar el solver tridimensional de doble precisión. En el
Importar la malla.

18
fichero tunel.msh. Obviamente, ese es el fichero que fue creado con el programa anterior
hexagonales (hexahedral cells).
NOTA: Debido al algoritmo que emplea GAMBIT para la generación de la malla, puesto que ésta
no fue definida perfectamente ortogonal, puede haber discrepancias en el número de celdas de unos
casos a otros.
Por otro lado, es interesante ver las zonas que se han cargado. Podemos comprobar cómo
aparecen las cuatro condiciones de contorno que definimos con anterioridad: inlet-fire, north-portal
y south-portal y symmetry. Además, por defecto, ha creado una pared wall que engloba todo el
resto de superficies exteriores del modelo.
mallado:

19
Asegúrese que los 5 objetos que aparecen bajo el epígrafe Surfaces están seleccionados
(todos menos el default-interior). A continuación haga clic en Display. Una ventana gráfica se abre
y la malla aparece dibujada en ella. Ya puede hacer clic en el botón Close de la ventana Grid
Display para ganar espacio en el escritorio. Como puede comprobar, la ventana gráfica no se cierra
y permanece visible.
ventana gráfica.
Practique estas operaciones para obtener una vista en el interior del túnel tal y como se
20

21

si estamos interesados en conocer la distribución de temperatura. En este caso, sí nos interesa. Por
tanto, se activa la ecuación de la energía y salimos haciendo clic en Ok

El flujo en el interior del túnel será turbulento. Seleccionamos como modelo de turbulencia
la opción k-épsilon. Ya que la malla es muy gruesa, seleccionamos también la opción de Standard
Wall Functions, con el objeto de que introduzca en el tratamiento cercano a los contornos sólidos la
ley logarítmica de pared. Haga clic en Ok.
Definición del modelo multiespecie.
Main Menu > Define > Models > Species > Transport & Reaction
Active la opción Species Transport. A continuación seleccione en el campo Mixture
Material la combinación carbon-monoxide-air.
Finalizamos haciendo clic en Ok.
22

Vamos a editar este “fluido” multiespecie. Para ello, hacemos clic en Edit…, botón que se
encuentra junto al campo names en Misture Species. Aparece un nuevo cuadro de diálogo que nos
permite añadir y eliminar componentes de la mezcla multiespecie que tenemos por defecto.
Añadimos el aire a la mezcla y eliminamos el O2, el CO, el H2O y el N2, de forma que el panel
original se transforma en el final:
Haga clic en Ok. Nuevamente sobre el panel de Materials, seleccione en Density la opción
ideal-gas. Termine haciendo clic en Change/Create para validar todos los cambios.

empleará el valor por defecto de 1 atm (101325 Pa) como Operating Pressure.
Además, es preciso definir el valor de la gravedad para tener en cuenta el efecto chimenea,
por lo que hay que activar la casilla Gravity. Vamos a empezar suponiendo que el túnel no tiene
pendiente. Por tanto:
23
Además, se ha fijado un valor de 300K como valor de referencia para la temperatura. Así
mismo, la densidad de referencia que use el modelo ha de corresponder con la del fluido a dicha
temperatura. Vamos a tomar el valor de 1.225 kg/m3 como válido. Haga clic en Ok para terminar.

GAMBIT como zonas en la izquierda de la ventana de boundary conditions. La condición
symmetry debería estar seleccionada por defecto. Si no es así, claque en ella. Compruebe que el tipo
seleccionado es symmetry. Después haga clic en Set… Fíjese cómo no hay nada que especificar
para este tipo de condición. A continuación, haga clic en Ok.
Muévase a lo largo de la lista y seleccione ahora la condición inlet-fire. FLUENT indica que
el tipo de esta condición de contorno es velocity-inlet. Recuerde que el tipo de condición fue ya
fijada en GAMBIT. De todas formas, es posible cambiar aquí el tipo de condición que se fijó en
GAMBIT si es necesario. Basta con seleccionar otro tipo de condición en esta ventana y validar la
selección.
Haga clic en Set…
Velocidad de los humos
Para determinar la velocidad a la que entran los humos por esa superficie, hemos de realizar
unos pequeños cálculos preliminares.
Se fijan unas condiciones estándar para el incendio que queremos modelizar. En este caso
(deflagración de un vehículo), dichas condiciones son:
• Potencia calorífica del incendio: 5 MW.
• Temperatura de los gases calientes: 400ºC.
• Concentración de CO2 de los gases de combustión: 0.05
24
Este modelo multiespecie no resuelve ecuaciones de combustión. Lo que hace es estudiar la

evolución de los gases de combustión que son emitidos desde la superficie del incendio. Con los
datos anteriores es posible determinar el caudal másicos de los gases calientes. Para ello,
imponemos que:
W&comb = m& g C p ΔT
donde W&comb es la potencia liberada en la combustión (los 5 MW), m& g es el caudal másico de
gases que queremos determinar, Cp es la capacidad calorífica del aire (1005 kJ/KgK) y ΔT es la
diferencia de temperatura entre los gases calientes y la temperatura local. (673 K – 300 K). Por
tanto:
5 ⋅106
m& g = ≈ 13.3 kg / s
1005 ⋅ ( 673 − 300 )
P
La densidad de los gases calientes la estimamos mediante la ley de gas ideal: = RT .
ρ
Luego:
101300
ρg = ≈ 0.52 kg / m3
287 ⋅ 673
Por tanto, el caudal volumétrico de gas que entra en el dominio será:
m& 13.3
Q& g = g = = 25.6 m3 / s
ρ g 0.52
Ahora bien, como el modelo es simétrico, todo es caudal entrará tanto por la superficie de
entrada del modelo como por su “imagen” simétrica. Por tanto, para definir la velocidad de entrada,
hay que dividir ese caudal de gases entre dos veces la sección de entrada y así:
Qg 25.6 m3 / s
Vg = = = 12.8 m / s
Sg 2 x1 m 2
Introduzca el valor de velocidad 12.8 m/s en el campo Velocity Magnitude.

Los valores de las variables turbulentas se dejan por defecto.
Finalmente, no olvide cambiar el valor de la temperatura de los gases a la entrada al valor
673 K. Asimismo, ha de fijar el valor 0.05 para el CO2 dentro del apartado Species Mass Fractions
Una vez haya introducido todos estos valores, haga clic en Ok.
25
6.- RESOLUCIÓN

En el campo Discretization, cambie la opción de Pressure a Body Force Weighted, puesto
que ofrece mejor convergencia que el resto de posibilidades. Cambie el resto de parámetros (6) de
Discretization a Second Order Upwind. Además, utilice como Pressure-Velocity Coupling el
método SIMPLEC.
Haga clic en Ok.

26
iterativo de resolución. En este caso, vamos a fijarle a todo el dominio valores nulos para todas las
variables.



En el campo Plotting especifique que los muestre en la Window 1 (por defecto trae la 0).
Haga clic en Ok.
27
Monitorización del flujo de humos en las salidas del dominio.
Mientras va iterando, vamos viendo la evolución del flujo de humo en las salidad. Es de
esperar que se observe un comportamiento asintótico que prediga que tras un cierto tiempo de
estabilidad, la cantidad de humo que sale por ambos portales ha de ser igual al que entra desde l
posición del incendio.
Main Menu > Solve > Monitors > Surface… Hacemos clic tanto en Plot como en Print y
definimos al programa que queremos ver los resultados cada Time Step (Every). A continuación
hacemos clic en el botón Define.
Seleccionamos en Report Type: Mass-flow Rate y le decimos al programa que nos los
muestre en todas las superficies externas: inlet-fire, north-portal y south-portal en la ventana 2
(Plot Window)
Visualización del proceso de difusión de humos.


continuación, se decimos que la actualiza cada paso temporal: Every: 1 y When: Time Step. A
continuación se hace clic en el botón Define.

28
Seleccionamos en categoría la opción Species… y en la subcategoría la opción Mass

fraction of CO2. En la opción Levels introducimos el valor 64. Y dejamos sin seleccionar la opción
Filled. Hacemos clic en Draw grid para ver el dominio y buscamos una vista como la anterior.
Hacemos clic en Display y ajustamos el zoom para ver con detalle el túnel. Ya que el túnel es largo,
podemos estirar la ventana a lo ancho para mejorar la vista Como paso final, hacemos Ok en el
panel de Animation Sequence y terminamos.
NOTA: Para dar renderizado al grid hay que ir a Main Menu > Display > Options y activar la
opción Lights On. Dejamos por defecto el renderizado de tipo Gouraud.
Si todo está correcto, deberíamos poder ver:
Como aún no ha comenzado el incendio, no vemos contornos de CO2 por el dominio

(únicamente se intuye en la entrada de humos debido al valor fijado por la condición de contorno.

Escriba el nombre tunel.cas como Case File. Haga clic en Ok. Compruebe que
tunel.cas.
29
Proceso iterativo para la convergencia.
Main Menu > Solve > Iterate. En la ventana desplegable hemos de imponer un paso
temporal del orden de 1 segundo. Se fijarán 50 iteraciones como número máximo por paso temporal
y 60 pasos temporales para completar 1 minuto completo de simulación de estado transitorio. Le
indicaremos también que refresque los valores de los residuos cada 10 iteraciones.
Una vez que todos los parámetros están fijados, haga clic en Iterate. Espere un rato a que el
programa ejecute los 60 pasos temporales. Al final, habrá obtenido un resultado como el siguiente:
En la evolución del flujo por todas las superficies podemos apreciar el primer instante en
que el humo comienza a salir por la boca más cercana así como el instante final en que también
alcanza la otra boca. Al final, el resultado asintótico tiende a cero, pues obviamente estamos
evaluando en la gráfica el flujo neto de todos los caudales que entran y salen del dominio.
30
7.- ANÁLISIS DE RESULTADOS.
Animación del movimiento del humo.
cualquiera.
En la animación vemos que el penacho de humo sobre la superficie de gases se establece de

forma instantánea. Esto se debe a que hemos fijado un valor de velocidad (12.8 m/s) muy grande
desde el principio. En las situaciones reales se suele fijar una rampa de forma que pasemos de
velocidad nula a los 12.8 m/s de forma progresiva en un determinado período de tiempo (esto es, el
coche no deflagra de forma instantánea…se va calentando y aumentando de potencia).
La animación es muy interesante porque nos muestra cómo se va estratificando el humo

desde la bóveda hacia abajo y cómo se escapa por las dos bocas del túnel.
Distribución de temperaturas.
Se muestra a continuación la distribución de temperaturas en la bóveda del túnel y en la

pared de simetría. Es interesante conocer la temperatura máxima que puede llegar a tenerse en la
cúpula para saber si se verán comprometidas las resistencias térmicas de los materiales de
construcción. Así mismo, una distribución de temperaturas permite saber los umbrales de
resistencia que personas atrapadas en el incendio podrían tolerar.
Main Menu > Display > Contours. Elegimos la categoría Temperatures (subcategoría
Static Temeprature) y seleccionamos que nos los muestre en primer lugar en la superficie de
symmetry. A continuación, desactivamos ésta y se lo pedimos sobre las paredes, wall. Hacemos clic
en la opción Filled y desactivamos que dibuje el grid (Draw grid). Si previamente hemos
desactivado el renderizado Gouraud del panel Options, obtendríamos como resultado:
31
Vamos a determinar las zonas peligrosas para las personas que pudiesen quedar atrapadas en
el interior. Vamos a fijar como máxima temperatura admisible el valor de 60ºC. Por encima de esos
valores, se considerará que no se puede sobrevivir. Vamos a definir la iso-superficie de 60ºC dentro
del modelo. Para ello:
Main Menu > Surface > Iso-surface… Elegimos una superficie de temperatura constante
dentro del plano de simetría: Surface of Constant Temperature (Static Temperature) From
Surface symmetry. También se podría coger la envolvente tridimensional de toda la zona fluid pero
es suficientemente representativo observarlo en el plano central. Le fijamos el valor de 333K (60ºC)
y la llamamos límite. Finalmente, hacemos clic en Apply.
A continuación vamos a representar una gráfica XY con esa curva: Main Menu > Plot >
XYPlot… Siguiendo la dirección z del plano de simetría (Plot Direction), representamos la
coordenada Y de esa curva de isotemperatura a 60ºC. (Y Axis Function: Grid… Y-Coordinate…
en Surfaces “limite”)
32
Variando los límites de los ejes y cambiando el aspecto de la curva, obtenemos un gráfico
como el siguiente:
Como vemos, ronda los 3.5 metros de alto, así que en la mayor parte del túnel (salvo la zona
próxima al fuego), la gente podría transitar buscando la salida sin riesgo elevado de quemaduras
(otro cantar es por asfixia).
Distribución de velocidades.
Vamos a mostrar por último patrones de flujo en la sección transversal del túnel. Para ello,
crearemos una sección transversal justo en el centro del túnel (posición z=13 m). Acudimos
nuevamente a Main Menu > Surface > Iso-surface… Elegimos la zona de fluido Fluid y
especificamos que se cree una superficie de valor constante de coordenada z: Surface of Grid… Z-
Coordinate From Zones Fluid. Le fijamos el valor de 2 m y la llamamos z01. Finalmente, hacemos
clic en Apply.
Repetimos el proceso para los valores de z, 6m, 10m y 20m, renombrando cada uno de ellos
como z02, z03 y z04 respectivamente. Ahora vamos a representar la componente z de la velocidad
sobre esos 4 planos. Como ya sabemos:
33
Main Menu > Display > Contours… Se elige la opción de Contours of… Velocity tomando
como subcategoría Z-velocity. Dejamos desactivada la opción Filled y ponemos un total de 15
Levels. Seleccionamos las 4 superficies (z01,…,z04) que acabamos de crear. Desactivamos la
opción Global Range y dejamos activa la de AutoRange. Dibujamos el grid para apreciar mejor la
perspectiva y hacemos clic en Apply.
Con este tipo de representaciones gráficas podemos entender mejor los patrones de flujo y el
comportamiento del humo al pasar por la bóveda del túnel.
Si repetimos el mismo gráfico de antes pero seleccionando además la superficie de simetría,

los contornos cambian completamente: observamos el flujo de gases con valor negativo de
componente z de velocidad que se mueve hacia la boca más cercana al fuego.
34
8.- EFECTO DE LA PENDIENTE Y DEL GRADIENTE DE PRESIONES
A continuación vamos a repetir el ejemplo de antes pero variando primero la pendiente del
túnel y luego introduciendo una diferencia de presiones entre las bocas.
Efecto de la pendiente.
Para introducir este efecto, vamos a modificar las componentes de la aceleración (es más
fácil esto que rehacer la geometría con inclinación). Vamos a darle al túnel una pendiente del 7.5%
en el sentido hacia la boca más cercana al túnel (north portal). Para ello:
North Portal
South Portal
gz = 0.735 m/s2
7.5%
4.3º
gy = -9.782 m/s2
g = 9.81 m/s2
Vamos al panel de Main Menu > Define > Operating Conditions y en Gravity introducimos
los valores que aparecen en el esquema:
A continuación, inicializamos desde cero el caso y volvemos a ejecutar el modelo durante

60s. Mostramos los flujos a la salida mientras itera así como la animación de la fracción de CO2 en
el plano de simetría. Para evitar problemas, renombramos la nueva animación como sequence-2.
En este caso la solución asintótica se alcanza aproximadamente a los dos minutos. Se

aprecia perfectamente el efecto del tiro del túnel debido a la pendiente. Incluso el penacho de humo
35
es claramente desviado hacia la boca norte. En un principio el humo inunda todo el túnel pero con
la estratificación establecida, el humo tiende a ir todo por la boca norte.
Puede comprobar cómo las distribuciones de velocidad son completamente distintas al caso
anterior. La zona de afección del incendio (en cuanto a temperaturas por encima de 60ºC y zonas de
gran concentración de dióxido de carbono) también cambia sustancialmente.
Efecto de la diferencia de presiones entre bocas.
Imaginemos que a los dos minutos de la situación anterior cambian bruscamente las
condiciones atmosféricas externas. Basta con se produzca una diferencia de 1 milibar entre los
extremos del túnel, para que se genere un importante depresión. Imaginemos que esa depresión de
100 Pa aparece en la boca sur. Introduzca esa variación en el modelo en las condiciones de contorno
de South Portal:
Main Menu > Define > Boudary Conditions… Seleccione south portal y haga clic en
Set… Fije como Gauge Pressure el valor -100 Pa (una depresión respecto al otro extremo del túnel
de 1 mbar).
36
Itere el caso de la solución anterior durante un minuto más a partir del cambio en las
condiciones de contorno y observe los cambios que se producen en el comportamiento del humo en
el interior del túnel. ¿Qué está ocurriendo en el túnel?
37
PRÁCTICAS DE
SIMULACIÓN DE LA DISPERSIÓN DE
CONTAMINANTES EN UN RÍO

6. RESOLUCIÓN
SIMULACIÓN DE LA DISPERSIÓN DE CONTAMINANTES EN UN RIO Área de Mecánica de Fluidos
El objetivo de esta práctica es ilustrar la resolución de un flujo no estacionario multi-especie

en un dominio tridimensional. Se pretende analizar el fenómeno de dispersión y transporte de un
vertido contaminante en el curso de un río.
Se comenzará mostrando las estrategias más habituales empleadas para la generación de

dominios tridimensionales en CFD. A continuación se darán las instrucciones necesarias para poner
en marcha el modelo multiespecie. Será necesario emplear una simulación transitoria mientras se va
alcanzando el estado asintótico para la dispersión del contaminante en el río. Se analizará la
influencia de algunos parámetros en la evolución de la dispersión del vertido, tales como la forma
del curso fluvial o el uso de sumideros para controlar el vertido.

¾ Crear una geometría tridimensional y mallarla en el programa de pretratamiento GAMBIT.
¾ Preparar un caso multiespecie, no estacionario para su resolución.
¾ Determinar las distribuciones finales de contaminante a lo largo del río.
¾ Generar animaciones que describan el estado transitorio del vertido y el efecto de medidas
de control en la difusión del contaminante.
Problema tipo.
Consideraremos un trozo característico de un río, con una longitud total de unos 425 metros.
Se definirá el curso del río a través de una serie de puntos así como la sección transversal del
torrente fluvial. La sección se extruirá a lo largo del curso del río con el objeto de obtener un
volumen que representará la masa de agua. Se incluirán un par de superficies que servirán como
puntos de desagüe. Por uno de esos desagües se introducirá en el modelo el vertido contaminante
mientras que el otro se empleará como punto de recogida de agua.
El caudal del río se estima en unos 100 m3/s. El vertido contaminante será de 0.1 m3/s.
QVERTIDO: 10 m3/s
QRÍO: 100 m3/s

Fig. 1. Geometría del río
1
Se definirá en primer lugar el curso del río mediante una serie de puntos. Al inicio del cauce
se creará la sección transversal del río. A continuación, se extruirá longitudinalmente dicha sección
para obtener el dominio tridimensinal. En ese punto, habrá que introducir las secciones planas que
simularán los desagües tanto de vertido como de captación. Para facilitar la creación de la geometría
tridimensional, en lugar de seguir la filosofía jerárquica de GAMBIT, procuraremos utilizar todas
las utilidades que incorpora el programa para generar volúmenes de forma rápida y sencilla.
Ha de crear una carpeta “Rio” dentro del subdirectorio que sea conveniente. Utilizaremos
esta sesión.
Inicie GAMBIT.
debe escribir: gambit –id rio y luego ejecute pulsando intro. Si esto no funciona, escriba en
la ventana el path completo:
c:\fluent.inc\ntbin\ntx86\gambit –id tunel
(-id) el nombre “rio”, le estamos indicando al GAMBIT que utilice ese nombre como prefijo del
2
3


El flujo es TRIDIMENSIONAL. Vamos a colocar el origen de coordenadas (0,0,0) justo en

uno de los vértices de la sección transversal del cauce fluvial.
Se colocará el origen del sistema de coordenadas en una de las esquinas de la sección

transversal del río, justo como un vértice de la superficie libre del río. Se va a crear la sección del
río a partir de una serie de puntos que unidos con un spline representarán la solera del cauce.
Creando vértices.
4
Vertex
Apply para aplicar. Con esto, se crea el vértice O (0,0,0).

que son creados.
Repita el proceso para crear los vértices que definirán el curso del río:
A (25,0,0) B (50,10,0) C (75,25,0) D (100, 0,0) E (75,-20,0)

F (50,-40,0) G (55,-60,0) H (75,-65,0) I (100,-65,0) J (125,-65,0)
Como puede ver, como estamos definiendo una línea bidimensional, la coordenada z se deja
siempre con el valor 0 por defecto.
Este comando ajusta los vértices más exteriores del cauce en la ventana gráfica para que
podamos ver todos los puntos en pantalla (se muestra la pantalla en la captura siguiente):
5
Creando el curso fluvial.
Ahora han de conectarse entre sí los vértices. Para ello, emplearemos líneas rectas (Straight)
para unir los puntos OA y IJ. El resto de puntos se unirán utilizando la opción NURBS.
6
Es muy importante que los vértices del NURBS se elijan en el orden correcto para que el
spline siga el camino que nosotros le indicamos.
Una vez que se haya seleccionado todos los vértices, es el momento de hacer clic en Apply
para generar el spline. Repita el proceso hasta que disponga de los lados rectos.
Recuerde que la opción NURBS aparece al dar sobre edge con el botón derecho del ratón.
NOTA: Como puede observar, sólo los vértices inicial y final del spline se han vuelto amarillos.
Ambos pertenecen al NURBS como punto inicial y final del mismo. El resto se han empleado para
que el programe calcule la interpolación más cercana que pasase por ellos… pero estrictamente no
pertenecen al spline recién creado y por eso siguen de color blanco.
7
Crear la cara transversal.
Para crear la sección perpendicular del río comenzamos introduciendo una serie de puntos:
P (0,0,-1) Q (0,6,-3) R (0,10,-1) S (0,10,0)
Creamos cuatro líneas: la superficie libre del río, uniendo con una línea recta el origen de
coordenadas y el punto S. La solera del cauce, uniendo con un spline los puntos PQR. Finalmente,
hacemos dos líneas verticales uniendo OP y RS. Para poder seleccionar los puntos recién creados
así como ver de forma óptima la geometría conviene rotar las vistas y hacer zoom dinámicos.
Al hacer modelos tridimensionales, las utilidades que implementa el ratón son muy útiles
para centrar o hacer zoom. Practique las siguientes opciones:
• Manteniendo pulsado el BOTÓN IZQUIERDO del ratón y arrastrando el puntero por la

pantalla gráfica podrá hacer giros tridimensionales de las vistas.
• Manteniendo pulsado el BOTÓN CENTRAL del ratón y arrastrando el puntero por la
pantalla gráfica podrá hacer centrar la geometría en la ventana gráfica.
• Manteniendo pulsado el BOTÓN DERECHO del ratón y arrastrando el puntero por la
pantalla gráfica podrá hacer zoom dinámico de la vista.
• Manteniendo pulsada la tecla de CTRL, se hace clic y se arrastra el ratón sobre una
determinada zona (la zona que quiero ver ampliada) y continuación se suelta el botón de
CTRL, conseguirá hacer zoom ventana.
Las líneas que ha debido obtener serían
Ahora hacemos la cara:

Form Face
8
Seleccionamos los 4 últimos lados que acabamos de crear y que forman un lazo cerrado.
Para seleccionarlos basta con presionar la tecla MAYÚSCULAS y a la vez ir haciendo clic sobre
Con esto, se activa la ventana de lista de lados. Haciendo clic sobre el los últimos y a
continuación al botónÆ , se seleccionan los lados que queremos. Haga clic en Close.
Como resultando, debería haber obtenido:
Extrusión de la cara.
Para obtener el volumen que representará la masa de agua del cauce fluvial vamos a extruir
la cara recién creada según la dirección del NURBS definido a partir de la línea media del río. La
vamos a extruir siguiendo la curva alabeada que representa el cauce del río, definiéndola
perpendicular al propio trayecto.
9
Stitch Faces . Haciendo clic con el botón derecho del ratón, desplegamos varias opciones y
seleccionamos de entre todas ellas la de extruir: Sweep Face:
En el menú desplegable que nos aparece seleccionamos la cara que tenemos para extruir y
elegimos en el campo Path la opción Edge. En la opción Type se elige el caso Perpendicular para
que se vaya extruyendo perpendicularmente a la línea media. Para hacer la extrusión hacemos clic
en Apply.
NOTA: Como podrá observar, las líneas son ahora verdes, lo cual significa que pertenecen a una
entidad jerárquica superior: un volumen. Para ver una perspectiva del volumen, haga clic sobre la
ventana gráfica con el botón izquierdo del ratón y muévase a izquierda o derecha con el ratón,
manteniendo pulsado dicho botón izquierdo.
10
Ya hemos extruido el primer volumen. Ahora repetimos la operación dos veces más con el
objeto de generar los volúmenes de río asociados a los otros dos trozos del trayecto estudiado.
Lógicamente, utilizamos la cara del volumen que se encuentra al inicio del spline para extruir.
Al final, deberá haber obtenido:
Generación de las zonas de desagüe y recogida.
Para marcar la superficie desde la cual se van a liberar los vertidos, es preciso crear ambas
caras y “restárselas” al volumen recién creado.
Form Face . Haciendo clic con el botón derecho del ratón, desplegamos varias opciones y
seleccionamos de entre todas ellas la de extruir: Sweep Edges: . Elegimos el lado vertical de la
margen izquierda del río, común a los dos primeros volúmes y lo extruimos en la dirección de x
negativa con una magnitud de 2 metros. (Ver capturas siguientes):
11
Hacemos clic en Apply y así conseguimos crear una cara que está “encima” del suelo del
volumen. Simplemente, nos basta con restar esa cara al volumen. Para ello, utilizamos el comando
Split:
Split Volume . Elegimos el volumen que tenemos creado y bajo la opción Split With elegimos
Face(Real) y ejecutamos haciendo clic en Apply
Al hacer esta operación, la cara lateral del volumen nos aparece partida en dos trozos.
Repetimos la operación con los otros tres lados verticales que separan los volúmenes. De
esta forma generamos 4 bocas en total para todo el dominio. OJO: Las líneas del volumen 3 debe
extruirlas en la dirección de x positiva. Actuando correctamente, el resultado final sería:
12
Para terminar, vamos a romper el volumen central en dos trozos. Tenemos que hacer esto
para poder mallar de forma óptima el río
MUY IMPORTANTE: Hay que pensar siempre en la malla, no en la geometría. La geometría sirve
como punto de apoyo de la malla…pero realmente el mallado es sobre lo que se discretizarán las
ecuaciones en FLUENT.
Entonces, en primer lugar creamos un plano, centrado en el diedro XZ, para romper el
volumen (Split). Creamos la cara con unas longitudes determinadas: 50 de ancho y 300 de alto bajo
la dirección ZX Centered.
Para hacer el split, simplemente: Operation Toolpad > Geometry Command Button >
Volume Command Button > Split Volume . Elegimos el volumen central y bajo la opción
Split With elegimos Face(Real), seleccionamos el plano recién creado y ejecutamos haciendo clic
en Apply
Para “limpiar” el modelo, podemos borrar esa cara de apoyo: Operation Toolpad >
Geometry Command Button > Face Command Button > Delete Face .
13
A continuación se detallan los pasos para crear las mallas de todos los volúmenes. Iremos
mallando desde la sección transversal en la zona de entrada al río y lo extruiremos (cooper)
longitudinalmente hasta llegar al final. Se utilizará un patrón progresivo con el objeto de optimizar
la malla en la zona próxima a las bocas de desagüe.
Mallado del primer volumen.
Edges
los lados, se puede hacer un zoom ventana de la zona de interés. Para ello, manteniendo pulsada la
tecla de CTRL, se hace clic y se arrastra el ratón sobre una determinada zona (la zona que quiero
ver ampliada) y continuación se suelta el botón de CTRL. Para retornar a la vista completa del
modelo, basta con volver a clicar en: Control Toolpad > Fit to Window Button
En la tabla adjunta se muestra el número de nodos (Interval Count) que debe fijar para cada
uno de los lados, así como el tipo de malla progresiva a emplear:
Lado Interval Count Type Ratio

Verticales 7 Sucessive Ratio 0.9
Horizontales 20 Bell-shaped 0.7
Los nodos de los lados verticales deben tender a acumularse hacia la solera del volumen.
Además, debe mallar las 5 secciones que se ven en el modelo para que al hacer la extrusión de
malla se mantengan los patrones progresivos, así como las líneas verticales de las bocas.
Al final, debe obtener lo siguiente:
Respecto a los lados longitudinales, vamos a definir varios patrones diferentes, en función
de que pertenezcan al volumen central partido o a los dos volúmenes que tienen las bocas de
vertido. Las distribuciones a adoptar se detallan a continuación:
14
• Para el volumen central partido en dos trozos:
Interval Count: 40
Sucessive Ratio: 1.09
Interval Count: 40
Interval Count: 60
Interval Count: 60
• Para los lados longitudinales de las bocas: Interval Count=10, equiespaciados.

• Para los lados longitudinales restantes de los volúmenes de entrada y salida: Interval
Count=20 con Sucessive Ratio=0.9
Como puede ver, la malla siempre se hace más fina hacia las bocas de entrada y salida de
vertidos.
15
Procedemos a mallar todos los volúmenes. Como se han mallado todas las líneas, podemos
hacer directamente un mallado hexagonal (sin extruir) del volumen. Comenzamos mallando el
primer volumen:
Operation Toolpad > Mesh Command Button > Volume Command Button >
Mesh Volume
Seleccione el primer volumen. Por defecto, el programa selecciona la opción más adecuada.
En este caso, lo más cómodo es emplear la opción Hex – Map. El Spacing que marca por defecto
será ignorado, pues ya ha sido especificada por nosotros la distribución de malla en todos los lados
del dominio.
Hacemos clic en Apply. Después, repetimos con los otros volúmenes adyacentes. Al final, la
malla del río debería tener el siguiente aspecto:
Para apreciar la malla con más claridad, es posible utilizar un renderizado de las superficies.
De esta forma, en vez de renderizar con un modelo de alambre, hacemos un Shaded.
16
Operation Toolpad > Global Control > Render Model
modelo. La zona del incendio será la entrada de velocidad para los gases, el techo de los postizos
serán condiciones de salida, el plano de simetría de la geometría se definirá como tal y el resto de
superficies serán paredes.
Inlet
Inlet
Wall
Symmetry
Wall
Outlet
17
Button
vamos a fijar que la zona del incendio es la entrada de humos. Para ello, en Entity: seleccionamos la
opción Faces para que GAMBIT sepa que queremos seleccionar una superfice (ya por defecto trae
activo Face).
Ahora seleccionamos la superficie como siempre, haciendo MAYÚSCULAS + clic sobre

ella. El área seleccionada aparece en el recuadro amarillo que está bajo el botón de Entity. Además,
podemos darle nombre dentro del campo Name:, por ejemplo inlet-rio.
Para fijar el Type:, seleccionamos MASS-FLOW-INLET. Hacemos clic en Apply. Con eso,

anterior. Agrupándose por condiciones de contorno, si se ha hecho todo correctamente, la lista final
de condiciones de contorno debería figurar de la siguiente manera en el programa:
De las cuatro bocas, la boca-1L (primera a la izquierda) se fija como mass-flow-inlet. El

resto se dejan como wall por el momento.
18
Salvar y exportar.
pueda ser interpretado por Fluent a continuación. Escriba tunel.msh como nombre de archivo
(File Name). Desactive la opción Export 2D Mesh, ya que efectivamente el caso es tridimensional.
Compruebe que un archivo rio.msh ha sido creado en su directorio de trabajo.
Run.
Importar la malla.

fichero rio.msh. Obviamente, ese es el fichero que fue creado con el programa anterior
hexagonales (hexahedral cells).
NOTA: Debido al algoritmo que emplea GAMBIT para la generación de la malla, puesto que ésta
no fue definida perfectamente ortogonal, puede haber discrepancias en el número de celdas de unos
casos a otros.
aparecen las cuatro condiciones de contorno que definimos con anterioridad: inlet-rio, outlet-rio
sup-libre, boca-1l, boca-1r, boca-2l y boca-2r. Además, por defecto, ha creado una pared wall que
engloba todo el resto de superficies exteriores del modelo.
19
mallado:

Asegúrese que los 9 objetos que aparecen bajo el epígrafe Surfaces están seleccionados
(todos menos el default-interior). A continuación haga clic en Display. Una ventana gráfica se abre
y la malla aparece dibujada en ella. Ya puede hacer clic en el botón Close de la ventana Grid
Display para ganar espacio en el escritorio. Como puede comprobar, la ventana gráfica no se cierra
y permanece visible.
20
ventana gráfica.
Practique estas operaciones para obtener una vista de la zona de vertido 1L, tal y como se
21


El flujo en el interior del túnel será turbulento. Seleccionamos como modelo de turbulencia
la opción k-épsilon. Ya que la malla es muy gruesa, seleccionamos también la opción de Standard
Wall Functions, con el objeto de que introduzca en el tratamiento cercano a los contornos sólidos la
ley logarítmica de pared. Haga clic en Ok.
22
Definición del modelo multiespecie.
Main Menu > Define > Models > Species > Transport & Reaction
Active la opción Species Transport. Finalizamos haciendo clic en Ok.

Vamos a editar este “fluido” multiespecie. Para ello, primero hemos de introducir como
material fluido el agua así como otra “agua” que hará las veces de contaminante.
En Material Type seleccionamos la opción fluid y a continuación hacemos clic en Fluent
Database. Nos sale un nuevo cuadro de diálogo. En él, elegimos de nuevo la opción fluid dentro de
material type. Buscamos el agua líquida (water-liquid [h2o(l)]), y hacemos clic en Copy. A
continuación hacemos clic en Copy de nuevo. El programa nos advierte que el agua ya está incluida
y que si queremos definir un nuevo material hay que cambiarle el nombre. En new-name lo
renombramos como contaminante y luego damos a Ok. Cerramos el DataBase. De vuelta en el
panel Materials en Material Type seleccionamos ahora Mixture y hacemos clic en Edit…, botón
que se encuentra junto al campo names en Mixture Species. Aparece un nuevo cuadro de diálogo
que nos permite añadir y eliminar componentes de la mezcla multiespecie que tenemos por defecto
(mixture template). Añadimos las dos materiales líquidos que acabamos de cargar (water-liquid y
contamante) y eliminamos el O2, el H2O y el N2, de forma que el panel original se transforma en el
final:
23
Haga clic en Ok. Nuevamente sobre el panel de Materials, seleccione en Density la opción
volume-weighted-mixing-law. Termine haciendo clic en Change/Create para validar todos los
cambios.
NOTA: Fíjese que la fórmula del agua es h2o, mientras que la del contaminante ha quedado como
h2o-new.

si estamos interesados en conocer la distribución de temperatura. En este caso, no nos interesa. Por
tanto, se mantiene desactivada la ecuación de la energía y salimos haciendo clic en Ok

empleará el valor por defecto de 1 atm (101325 Pa) como Operating Pressure.
En este modelo no es preciso definir el valor de la gravedad para tener en cuenta su efecto,
puesto que no estamos modelizando la superficie libre del río (la imponemos como condición de
contorno de simetría). Por tanto:
A continuación se fijarán los valores de flujo a la entrada del río y para la boca de
contaminante.

GAMBIT como zonas en la izquierda de la ventana de boundary conditions. Haga clic en la
condición de sup-libre. Compruebe que el tipo seleccionado es symmetry. Después haga clic en
Set… Fíjese cómo no hay nada que especificar para este tipo de condición. A continuación, haga
clic en Ok.
24
Muévase a lo largo de la lista y seleccione ahora la condición inlet-rio. FLUENT indica que
el tipo de esta condición de contorno es mass-flow-inlet. Recuerde que el tipo de condición fue ya
selección.
Haga clic en Set…
Introduzca el valor de flujo másico 100000 kg/s en el campo Mass Flow-Rate [kg/s]. Este
valor sale de multiplicar los 100 m3/s del caudal del río por la densidad del agua (aprox. 103 kg/m3).
En Direction Specification Method elija la opción Normal to boundary.
Asimismo, mantenga el valor 0 para el h2o-new dentro del apartado Species Mass
Fractions. Una vez haya introducido todos estos valores, haga clic en Ok.
Ahora especifique las condiciones de contorno para la boca-1l. Introduzca en mass flow-rate
el valor 10000 kg/s y fije la dirección del flujo perpendicular a la condición (Normal to boundary).
En el apartado Species Mass Fractions, introduzca el valor 1 para el contaminante (h2o-new).
El resto de condiciones se pueden dejar por defecto.

25
6.- RESOLUCIÓN

En el campo Discretization, cambie la opción de Pressure a Body Force Weighted, puesto
que ofrece mejor convergencia que el resto de posibilidades. Cambie el resto de parámetros (4) de
Discretization a Second Order Upwind. Además, utilice como Pressure-Velocity Coupling el
método SIMPLEC.
Haga clic en Ok.
iterativo de resolución. En este caso, vamos a fijarle a todo el dominio valores nulos para todas las
variables.

26

En el campo Plotting especifique que los muestre en la Window 1 (por defecto trae 0).
Haga clic en Ok.
Monitorización del flujo a la salida del dominio.
Mientras va iterando, vamos viendo la evolución del flujo de agua y contaminante en las
salidas. Es de esperar que se observe un comportamiento asintótico que prediga que tras un cierto
tiempo de estabilidad, la cantidad de contaminante que entra es igual al que sale.
Main Menu > Solve > Monitors > Surface… Hacemos clic tanto en Plot como en Print y
definimos al programa que queremos ver los resultados cada Time Step (Every). A continuación
hacemos clic en el botón Define.
Seleccionamos en Report Type: Flow Rate y le decimos al programa que nos los muestre en
todas las superficies outlet-rio y boca-1l en la Plot Window 2.
27
Visualización del proceso de difusión del contaminante.


continuación, se decimos que la actualiza cada paso temporal: Every: 1 y When: Time Step. A
continuación se hace clic en el botón Define.

Seleccionamos en categoría la opción Species… y en la subcategoría la opción Mass
fraction of h2o(l)-new. Le decimos que los muestre sobre la surface sup-libre. En la opción Levels
introducimos el valor 64 y desactivamos la opción Filled. Hacemos clic en Draw grid para ver el
dominio y buscamos una vista como la anterior. Hacemos clic en Display y ajustamos el zoom para
ver con detalle el túnel. Ya que el túnel es largo, podemos estirar la ventana a lo ancho para mejorar
la vista Como paso final, hacemos Ok en el panel de Animation Sequence y terminamos.
NOTA: Para dar renderizado al grid hay que ir a Main Menu > Display > Options y activar la
opción Lights On. Dejamos por defecto el renderizado de tipo Gouraud.
Si todo está correcto, deberíamos poder ver algo así:
28
Como aún no ha comenzado el vertido, no vemos contornos de contaminante por el dominio

(todo es valor cero).

Escriba el nombre rio.cas como Case File. Haga clic en Ok. Compruebe que
el programa FLUENT, podrá recuperar todo su trabajo simplemente leyendo el archivo rio.cas.
Main Menu > Solve > Iterate. En la ventana desplegable hemos de imponer un paso
temporal del orden de 5 segundos. Se fijarán 50 iteraciones como número máximo por paso
temporal y 24 pasos temporales para completar 1 minuto de simulación de estado transitorio. Le
indicaremos también que refresque los valores de los residuos cada 10 iteraciones.
29
7.- ANÁLISIS DE RESULTADOS.
Al cabo de dos minutos, hemos encontrado la siguiente situación:
El flujo de contaminante a la salida comienza a igualarse con el de entrada por lo que el neto
de ambos tiende asintóticamente a cero. Sin embargo, la distribución de contaminante aún no se ha
establecido por completo en el río. En los meandros, aún tiene que establecerse un patrón
estacionario para las líneas de corriente del contaminante:
Vamos a simular otros 3 minutos, en espera de que la masa de agua en el río se estabilice
por completo. Después de esos tres minutos, veremos cómo queda la distribución de fracción de
contaminante en el agua… Aplicamos entonces otros 36 pasos temporales.
30
Al cabo de 5 minutos, vemos que se ha alcanzado un estado cuasi-estacionario. Hay ciertas

oscilaciones del caudal de contaminante: esto puede indicar que en ciertas zonas (meandros) se
almacenan “bolsas” de vertido que cuando alcanzan cierto valor crítico son arrastradas por la
corriente principal. Por tanto, tenemos una fluctuación a la salida como consecuencia de este efecto.
La distribución final de vertido es:
Animación del transporte de contaminante.
cualquiera.
31
En la animación vemos cómo se va generando una acumulación de contaminante en la parte

interna del meandro, que oscila y que seguramente es responsable de esa oscilación de masa
contaminada en la salida del dominio.
El contaminante se difumina por toda la corriente de agua, por eso a la salida las
concentraciones son muy bajas, pero están completamente extendidas entre el caudal de agua
circulante.
Distribuciones a varias alturas
Hasta ahora hemos analizado el comportamiento en la superficie libre, pero el contaminante

se transporte también por debajo de la superficie. Vamos a crear planos a distintas cotas para ver la
distribución de la fracción de masa del vertido en el río. Para ello, creamos los planos a cota z = -0.5
y cota z = - 2 metros.
Main Menu > Surface > IsoSurface… En Surface of Constant elegimos la opción Grid…
Z-Coordinate y en Iso Values tomamos el valor -0.5. Le decimos que obtenga esa superficie desde
la zona fluid, marcada en el panel. La renombramos como z0.5. Repetimos el procedimiento para la
cota -2 metros.
Ahora representamos la distribución de contaminante entre los valores 0 y 0.5 con un

Display > Contours en la superficie libre y en las otras dos cotas y los comparamos:
32
Superficie libre
Cota – 0.5 m
Cota – 2 m
¿Qué conclusiones se pueden obtener de estas distribuciones?
33
Distribuciones de velocidad.
También es interesante estudiar los patrones de velocidad ya que éstos explican el

comportamiento del arrastre del contaminante. Vamos a observar los patrones en un plano
transversal, situado a x=100 metros del origen.
Para ello, creamos el plano como siempre, con la utilidad IsoSurface… En Surface of
Constant elegimos la opción Grid…X-Coordinate y en Iso Values tomamos el valor 75. Le
decimos que obtenga esa superficie desde la zona fluid, marcada en el panel. La renombramos
como x75. Hacemos un contorno de velocidad absoluta en ese plano.
Vamos a crear una imagen fotorrealista en la que crearemos la IsoSurface de concentración

de contaminante de valor 0.5 y sobre ella dibujaremos vectores velocidad. Esa imagen la
fundiremos con la que acabamos de crear. Para hacer todo esto:
1) Hacemos la Iso-Surface de Species… Mass fraction of h2o-new from zones fluid para
valor 0.5
2) Como tenemos la imagen anterior aún abierta, vamos a Main Menu > Display > Scene y
activamos la opción Overlays y pulsamos en Apply.
3) Hacemos un Display > Vectors sobre la superficie que se acaba de crear ahora. Antes de
hacer Display especificamos que los vectores sean de tipo filled-arrow y en Vector
Options fijamos que tengan una Scale Head de 0.4. En Scale fijamos el valor 4.
4) Finalmente hacemos Display:
34
Parada del vertido.
Podemos reanudar el análisis suponiendo que se para el vertido de pronto. Analizaremos el

tiempo que se tarda en “limpiar” por completo el río de contaminante. Para ello, en condiciones de
contorno, introducimos un valor cero de entrada de caudal de vertido en la superficie boca-1l.
Ejecutamos con la animación activada igual que antes, otros 2 minutos de simulación (24 time
steps).
Al final del proceso hemos obtenido:
Como vemos tras dos minutos prácticamente hemos limpiado el curso del río. La última
acumulación de contaminante está próxima a la salida.
Podemos hacer una interesante animación del proceso de limpieza de agua contaminada.
EJERCICIO PROPUESTO: Otras situaciones.
Se propone al alumno que estudie otras posibilidades. Por ejemplo, se podría estudiar la
situación con el vertido de contaminante saliendo por la boca de la otra margen del río.
Habría que buscar el estado estacionario de esa situación y a continuación plantear un
nuevo cese de vertido.
En ese punto, se propone que desde las bocas cercanas a la salida se aspire caudal con el
objeto de evitar que la gran concentración que se desplaza hacia aguas abajo siga su curso. Para
ello, basta definir unas en esas superficies una condición de contorno de velocity-inlet pero con
signo contrario. Se puede fijar una velocidad en cada una de ellas de forma que el caudal (la sección
de las bocas es de 2 m2) aspirado sea el doble del que había de vertido (20 m3/s).
¿Es efectiva esa aspiración? ¿lo es más desde una margen que desde la otra?
35
PRÁCTICAS DE
SIMULACIÓN DE LA COMBUSTIÓN
GASEOSA EN UN QUEMADOR

6. RESOLUCIÓN Y ANÁLISIS DE RESULTADOS
7. EJERCICIOS PROPUESTOS
SIMULACIÓN DE LA COMBUSTIÓN GASEOSA EN UN QUEMADOR Área de Mecánica de Fluidos
En esta práctica se va a considerar la modelización con FLUENT de una cámara de

combustión cilíndrica. Los objetivos que se plantean son los siguientes:
¾ Introducción a los problemas de mezcla y combustión gaseosa con herramientas CFD.

¾ Estudio de la llama de difusión turbulenta en un quemador axisimétrico de metano/aire.
¾ Selección de diversos modelos para la resolución:
o Elección del modelo para la fracción de mezcla y la reacción química.
o Calor específico constante y variable
¾ Postproceso de los resultados.
Problema tipo.
Vamos a considerar una cámara de combustión cilíndrica como la que se muestra en un

corte en la figura 1. Un chorro de metano penetra en la cámara a una temperatura de 300 K y una
velocidad de 80 m/s a través de una pequeña tobera cilíndrica. Coaxialmente penetra por su exterior
un chorro de aire a 300 K y 0.5 m/s. Dado que el número de Reynolds basado en el chorro de
metano a la entrada de la cámara es Re ∼ 28000, y que las corrientes de ambos gases no están
premezcladas, se formará una llama de difusión turbulenta, cuyo estudio es el objetivo de la
práctica.
La combustión del sistema metano/aire será tratada usando un modelo de química finita,
mediante una reacción global de un paso, de la forma: CH4 + 2O2 → CO2 + 2H2O. Supondremos
que hay conversión total del combustible a CO2 y H2O. Para definir las reacciones debemos
especificar los coeficientes estequiométricos, las entalpías de formación y los parámetros que
controlan la velocidad de reacción.
Una buena aproximación para combustión turbulenta es considerar que la velocidad viene
determinada por la mezcla turbulenta de ambos gases. Estamos suponiendo que la velocidad de
reacción es mucho mayor que la velocidad característica de los torbellinos turbulentos que provocan
la mezcla. Por tanto, en esta aproximación se ignora la cinética química y se supone que la reacción
se produce instantáneamente una vez que combustible y comburente entran en contacto. En
FLUENT disponemos de dicho modelo de reacción limitada por la mezcla turbulenta bajo el
nombre de eddy dissipation model.
Paredes adiabáticas, 300K
Aire
T=300K
U=0.5 m/s
0.225 m
4 mm Metano
T=300K ; U=80 m/s Bujía
Eje axisimétrico
10 mm
1.8 m
Fig. 1. Esquema del problema.

1
Para simular el dominio del quemador se creará una única superficie 2D. Para ello,
empezaremos creando los vértices (vertex) de los extremos. Se introducirán vértices adicionales en
la pared del fondo para garantizar que la malla sea ortogonal. A continuación se unirán los vértices
correspondientes mediante líneas para crear los lados (edges) del rectángulo. Por último, crearemos
la superficie (face) que corresponde al área encerrada por los lados. En el punto 3 de esta guía se
mallará la superficie del rectángulo.
faces -> volumes.
Ha de crear una carpeta “Quemador” dentro del subdirectorio que sea conveniente.
Inicie GAMBIT.
debe escribir: gambit –id quemador y luego ejecute pulsando intro. Si esto no funciona,
c:\fluent.inc\ntbin\ntx86\gambit –id quemador
(-id) el nombre “quemador”, le estamos indicando al GAMBIT que utilice ese nombre como prefijo
del archivo que vamos a construir a lo largo de esta sesión.
2
3


El flujo es BIDIMENSIONAL. Vamos a colocar el origen de coordenadas (0,0) justo en la

esquina inferior izquierda de la tobera. Veremos que el dominio del problema quedará comprendido
entre:
0 ≤ x ≤ 1.8, 0 ≤ y ≤ 0.225
donde x e y definen las coordenadas horizontal y vertical respectivamente.
Se colocará el origen del sistema de coordenadas en la esquina izquierda de la tobera, justo

como el vértice de inicio de la misma. Se señalan en el esquema adjunto los puntos que es preciso
introducir, con sus cotas asociadas.
4
G(1.8,0.225)
F(0,0.225)
E(0,0.009) D(0.01,0.009) H(1.8,0.009)

B(0,0.005) C(0.01,0.005) I(1.8,0.005)
A(0,0) J(1.8,0)
Creando vértices.
Vertex

que son creados.
Repita el proceso para crear los 9 vértices restantes (puntos B,C,D,E,F,G,H,I y J).
5
Como puede ver, como estamos definiendo un dominio bidimensional, la coordenada z se

deja siempre con el valor 0 por defecto.
Este comando ajusta los 4 vértices más exteriores de la sección en la ventana gráfica para
que podamos ver todos los puntos en pantalla (se muestra la pantalla en la captura siguiente):
Los puntos que definen la tobera de entrada para el metano aparecen en la vista
completamente superpuestos por culpa de la vista global. Para verlos en detalle, se puede hacer un
zoom ventana de la zona de interés. Para ello, manteniendo pulsada la tecla de CTRL, se hace clic
y se arrastra el ratón sobre una determinada zona (la zona que quiero ver ampliada) y continuación
se suelta el botón de CTRL. Para retornar a la vista completa del modelo, basta con volver a clicar
en: Control Toolpad > Fit to Window Button. De esa forma:
6
Creando lados.
Ahora han de conectarse entre sí los vértices del dominio.

en Apply para generar el lado. Repita el proceso hasta que disponga de todos los lados que se ven en
la siguiente figura: todos los lados exteriores más los lados CH y DI (en total 10 lados) .
NOTA: Todos los vértices pertenecen ahora a algún tipo de línea y por esa razón el programa los
colorea de amarillo.
Crear las caras.
Form Face
Para crear las caras que representen el quemador es preciso seleccionar todos los lados y que
7
cada línea y al final soltar la tecla una vez queden seleccionadas las líneas en rojo. También es
Con esto, se activa la ventana de lista de lados. Eligiendo en Available los lados precisos y
haciendo clic sobre el botón --> , se seleccionan todos los lados. Haga clic en Close y finalmente
haga clic en Apply para crear la cara.
Repita esta operación para crear las tres caras del modelo.
A continuación se detallan los pasos para crear una malla en la cara rectangular con 60
divisiones en horizontal y 30 en vertical. Se mallarán en primer lugar los lados y a continuación se
mallarán las caras. El espaciado en la malla se especifica a través del mallado de los lados.
Edges
los lados, se puede hacer un zoom ventana de la zona de interés. Para ello, recuerde mantener
pulsada la tecla de CTRL, se hace clic y se arrastra el ratón sobre una determinada zona (la zona
que quiero ver ampliada) y continuación se suelta el botón de CTRL. Para retornar a la vista
completa del modelo, basta con volver a clicar en: Control Toolpad > Fit to Window Button
A continuación se incluye una tabla donde se dan todos los parámetros necesarios para
mallar adecuadamente los 12 lados del modelo. Observará que al hacer MAY+clic sobre un lado,
éste se coloreará en rojo con una flecha que indica el sentido del vector director del lado. Si se
utiliza un patrón de nodos progresivo puede ser necesario que la flecha tenga un sentido
determinado. Para cambiarlo, basta con hacer clic en la tecla Reverse.
8
En la tabla adjunta se muestra el número de nodos (Interval Count) que debe fijar para los
lados, el tipo de progresión de malla (Type), el parámetro que controla la progresión (Ratio) y el
sentido de la flecha:
Lado Interval Tipo (type) Ratio Sentido

(edge) Count
AB 5 Successive Ratio 1 Indiferente
BC 5 Successive Ratio 1 Indiferente
CD 5 Successive Ratio 1 Indiferente
DE 5 Successive Ratio 1 Indiferente
EF 20 Exponent 0.2 De F a E
FG 60 First Length 0.1275 De G a F
GH 20 Exponent 0.2 De G a H
HI 5 Successive Ratio 1 Indiferente
IJ 5 Successive Ratio 1 Indiferente
JA 60 First Length 0.1275 De J a A
CI 55 First Length 0.1275 De I a C
DH 55 First Length 0.1275 De H a D
Como se podrá ver, el objetivo de la malla progresiva es apilar nodos en la zona de entrada
de los gases, donde se espera encontrar mayores gradientes.
Mallado de las caras.
Faces
Haga MAYÚSCULAS + clic en cada cara para ir seleccionándolas y luego haga clic en
Apply. Se emplearán Elements de tipo Quad y el campo Type se seleccionará la opción Map. Esta
opción me devuelve malla estructurada. Haga clic en Apply para mallar.
En detalle, sobre la zona de entrada de gases:
9
modelo. Por la izquierda se tienen las entradas de los flujos. Abajo, se fijará una condición de eje
axisimétrico. A la salida colocaremos una condición de salida de presión para los gases calientes.
La zona de la tobera y la linea horizontal superior serán paredes..
Pared (Wall)
Inlet-aire
Oulet
Tobera
Inlet-gas Eje (Axis)
Button
vamos a fijar que el lado izquierdo es la entrada de aire. Para ello, en Entity: seleccionamos la
Face).
Ahora seleccionamos el lado como siempre, haciendo MAYÚSCULAS + clic sobre él. El
lado seleccionado aparece en el recuadro amarillo que está bajo el botón de Entity. Además,
podemos darle nombre dentro del campo Name:, por ejemplo inlet-aire.
Para fijar el Type:, seleccionamos VELOCITY_INLET. Hacemos clic en Apply. Con eso,
10

anterior. Agrupándose por condiciones de contorno:

AB Inlet-gas VELOCITY_INLET
BC,CD,DE Tobera WALL
EF Inlet-aire VELOCITY-INLET
FG Pared WALL
GH,HI,IJ Outlet PRESSURE-OUTLET
JA Eje AXIS
Si se ha hecho todo correctamente, la lista final de condiciones de contorno debería figurar de la

siguiente manera en el programa:
11
Salvar y exportar.
pueda ser interpretado por Fluent a continuación. Escriba quemador.msh como nombre de
archivo (File Name). Active la opción Export 2D Mesh, ya que efectivamente el caso es
bidimensional. Haga clic en Accept.
Compruebe que un archivo quemador.msh ha sido creado en su directorio de trabajo.
Run.
Importar la malla.

fichero quemador.msh. Obviamente, ese es el fichero que fue creado con el programa anterior
rectangulares (quadrilateral cells). No podía ser de otra manera ya que se usaron 60 divisiones
horizontales por 30 divisiones verticales (menos las 25 celdas de la pared de la tobera).
aparecen las seis condiciones de contorno que definimos con anterioridad: inlet-aire, inlet-gas,
tobera, pared, outlet y axis. Además, por defecto, ha creado una zona fluid que engloba todas las
celdas interiores del modelo.
12
mallado:

ventana gráfica.
13
Practique estas operaciones para obtener una vista de la zona anterior a la placa tal y como
se muestra en la figura siguiente:
14

El problema que vamos a resolver es estacionario, bidimensional y axisimétrico. Elija la
opción Axisymmetric bajo el campo Space. Mantenga la opción de Steady bajo el cuadro Time. Se
emplearán las otras opciones que por defecto indica el programa, como el tipo segregated, implicit
formulation y absolute velocity formulation. Haga clic en Ok.

Para resolver la turbulencia vamos a emplear un modelo k-ε con Standard Wall Functions
para resolver el flujo cercano a los contornos sólidos (la malla que estamos empleando es
extremadamente gruesa, por lo que es seguro que nuestro y+ será muy alto). Salimos haciendo clic
en Ok.
.
15

Como estamos analizando un proceso de combustión, es preciso resolver la ecuación de la
energía para caracterizar la distribución de temperatura. Por tanto, se activa la ecuación de la
energía y salimos haciendo clic en Ok
Definición de las especies químicas. Modelo de reacción.
Main Menu > Define > Models > Species > Transport & Reaction…Para indicar que
queremos calcular las especies químicas presentes accedemos a la opción Species transport. El
panel se extiende para presentar más opciones. En el cuadro Mixtures Material buscamos en la lista
disponible en el desplegable la mezcla metano-aire (methane-air) y la seleccionamos. Vemos que el
número de especies volumétricas se ajusta automáticamente a 5, siendo éste el número de especies
químicas presentes en la mezcla metano-aire: CH4, O2, CO2, H2O y N2.
Seguimos en el panel de Species Model. En el cuadro Reactions activamos la opción
Volumetric para permitir la reacción química en el interior del volumen considerado (el quemador).
El panel se modifica y aparece un nuevo cuadro de opciones titulado Turbulence-Chemistry
Interaction. En él activamos la opción Eddy-Dissipation correspondiente al modelo explicado en la
introducción del guión.
Finalmente en el cuadro de opciones, dejamos la opción Diffusion Energy Source que viene
activada por defecto. Hacemos clic en Ok para hacer efectivos los cambios y cerrar el panel.
Al salir del panel de Species Model nos aparece un panel de información indicando que el
material o los métodos han cambiado y que es necesario confirmar las propiedades antes de
continuar. También podemos ver que en la ventana de comandos aparecen varios mensajes
informativos indicando que es necesario suministrar datos para las propiedades físicas del material
definido. Hacemos clic en Ok para cerrar el panel de información y vamos a proceder a especificar
las propiedades físicas de la mezcla.
16
Una vez establecidas las especies químicas presentes, hemos de especificar sus propiedades
termodinámicas, así como las opciones del modelo de reacción. Esto podemos hacerlo en el panel
Main Menu > Define > Materials… En la opción Mixture Species, si pulsamos el botón Edit…
para comprobar de qué especies se compone la mezcla. En la opción Reaction aparece cargado el
modelo Eddy-Dissipation que elegimos anteriormente. Haciendo clic en Edit… accedemos a los
parámetros que controlan la cinética química del problema:
Las opciones Arrhenius Rate están desactivas ya que en este modelo hemos indicado que la
química se supone mucho más rápida que la mezcla turbulenta, por lo que no es necesario
especificar parámetros relacionados con la cinética química. Para los valores A y B del subpanel
Mixing Rate dejamos los valores por defecto. Compruebe cómo el resto de parámetros es el
adecuado para nuestro problema: una reacción volumétrica con dos reactivos (reactants) y dos
productos (products) con unos coeficientes estequiométricos (Stoich. Coefficient) que
efectivamente siguen la reacción de combustión del metano (CH4 + 2O2 → CO2 + 2H2O). Por tanto,
mantenga los parámetros como están y salga pulsando Ok.
De vuelta en el panel Materials, en la opción Density escogemos incomprensible-ideal-gas.
Moviendo la barra deslizadora del cuadro Properties podemos comprobar que el resto de
propiedades son constantes por defecto. De momento, vamos a dejarlas así, aunque más adelante
haremos cálculos con capacidad calorífica variable.
Una vez establecidas todas las propiedades y opciones de reacción, pulsamos el botón
Change/Create en el panel Materials para activar todos los cambios realizados. Finalice haciendo
clic en Close para cerrar el panel.
17

empleará el valor por defecto de 1 atm (101325 Pa) como Operating Pressure. Haga clic en Cancel
para mantener el valor por defecto.
Main Menu > Define > Boundary Conditions…

GAMBIT como zonas en la izquierda de la ventana de boundary conditions. Haga clic sobre la
condición de contorno eje. Compruebe que el tipo seleccionado para ella es Axis. Después haga clic
en Set… Fíjese cómo no hay nada que especificar para este tipo de condición. A continuación, haga
clic en Ok.
18
Muévase a lo largo de la lista y seleccione ahora la condición inlet-aire. FLUENT indica

que el tipo de esta condición de contorno es velocity-inlet. Recuerde que el tipo de condición fue ya
selección.
Haga clic en Set…Introduzca el valor de velocidad 0.5 m/s en el campo Velocity
Magnitude. Mantenga el valor por defecto de Temperatura a la entrada, 300 K. Para especificar la
turbulencia del aire de entrada seleccione en el desplegable Turbulence Specification Method el
método Intensity and Hydraulic Diameter. Imponga una intensidad de la turbulencia del 10% y un
diámetro hidráulico de 0.4 m. Además, hemos de especificar la composición de la mezcla entrante,
para lo cual introducimos el valor 0.23 (fracción másica del oxígeno en el aire) en el cuadro
correspondiente al O2. Acepte presionando Ok.
Repita los mismos pasos para la velocity-inlet con el metano (condición inlet-gas). Recuerde
introducir un valor de 80 m/s para el chorro de metano así como una temperatura de 300 K.
Seleccionamos el mismo método de especificación de la turbulencia que en la entrada de aire, con
una intensidad turbulenta del 10% pero ahora con un diámetro hidráulico de 0.01 m. Fije la
composición de metano (CH4) al 100% colocando el valor 1 en la casilla ch4 del Species Mass
Fractions.
En la condición de Outlet, resulta apropiada una condición de contorno del tipo pressure-
outlet. Mantenga el valor de 0 gauge pressure por defecto así como la Backflow Total temperature
en 300K. Para las variables turbulentas seleccione la opción habitual y fije un 10% de intensidad y
0.45 m de diámetro hidráulico. En Species Mass Fractions fije el valor 0.23 para el oxigeno y 0
para el resto de especies
Por último, vamos a fijar las condiciones de las paredes. En la pared horizontal de arriba,
bajo la pestaña Thermal, elija la opción Temperature en Thermal Options e imponga una
temperatura constante de 300K. Sin embargo, en la condición tobera, seleccione la opción Heat
Flux y mantenga el valor de 0 W/m2 para simular condiciones adiabáticas. Valide los cambios
haciendo clic en Ok.

19
6.- RESOLUCIÓN Y ANÁLISIS DE RESULTADOS
Para la resolución se emplearán discretizaciones de primer orden que proporcionan una

mejor convergencia para modelos complejos.

En el campo Discretization, mantenga todas las opciones de Momentum, Turbulence,
Energy y Species Fraction en First Order Upwind. Además, mantenga como Pressure-Velocity
Coupling el método SIMPLE.
Haga clic en Ok.


momento (x-momentum, y-momentum), de la turbulencia y de las especies. Deje los residuos de la
energía (energy) en el criterio 1e-06.
En el campo Plotting especifique que los muestre en la Window 0 (por defecto).
Haga clic en Ok.
Con esto se termina todo el proceso de inicialización y definición del proceso a resolver.
20
Main Menu > File > Write > Case&Data… Guardamos el caso y los datos en un archivo
quemador.cas que podremos recuperar cuando queramos. FLUENT genera un archivo *.cas
en el que almacena toda la información relativa al modelo: geometría, modelos y algoritmos de
cálculo, valores de condiciones de contorno, parámetros y criterios de convergencia… Por otro
lado, los valores de las variables en las celdas (los datos), los guarda en otro archivo
complementario de extensión *.dat.
Compruebe que una vez hecha la operación de escritura, se han creado realmente dos
archivos nuevo en su directorio de trabajo bajo el nombre quemador.cas.
6.1.- SOLUCIÓN FRÍA
Comenzaremos el análisis de nuestro problema calculando el campo fluido que se produciría

sin la existencia de reacción química. Para ello, vamos a desactivar la reacción de combustión:
Main Menu > Define > Models > Species > Transport&Reaction… y deshabilite la opción
Volumetric Reactions. Haga clic en Ok para validar.
Ahora hemos de preparar el comienzo del cálculo para lo cual inicializamos la solución.
Main Menu > Solve > Initialize > Initialize… En el subpanel Compute From
seleccionamos la opción All Zones. Haga clic en Init.
Main Menu > Solve > Iterate. En la ventana desplegable fijamos que haga 500 iteraciones
para empezar y que nos vaya actualizando los resultados cada 10 de ellas:
Hacemos clic en Iterate. A las 530 iteraciones, el caso converge:
21
quemador-frio.cas que podremos recuperar cuando queramos.
Resultados de la solución fría
Podemos analizar el campo fluido representando líneas de corriente, que nos darán idea de la
trayectoria de las partículas desde la entrada. Para ello, hacemos: Main Menu > Display >
Contours… Seleccionamos en Contours Of la opción Velocity y debajo Stream Function.
Dejamos sin seleccionar ninguna Surface y desactivamos la opción Filled. Hacemos clic en
Display.
El modelo multiespecie que se utiliza aquí tiene campos de velocidad y presiones únicos, de
forma que vemos líneas de corriente tanto del aire como del metano. Obsérvese cómo el chorro de
metano “arrastra” aire en la zona próxima a la tobera (ingiere el aire). Lógicamente, esa capa de
cortadura junto con el efecto progresivo del mezclado por torbellinos turbulentos permite una
mezcla eficiente de de ambas corrientes (y una buena combustión en el quemador). Se crea
adicionalmente una zona de recirculación muy extensa cercana a la pared superior como
consecuencia de la diferencia de cantidad de movimiento de ambos chorros.
22
Vamos a analizar ahora cada especie por separado. Para ello, vamos a visualizar los
contornos de concentración de combustible (ch4) y de oxidante (o2). Además, para apreciar más
fácilmente la forma de los chorros, vamos a activar la vista simétrica del modelo: Main Menu >
Display > Views… En el campo Mirror planes haga clic sobre la condición eje. A continuación,
haga clic en Apply y podrá observar cómo se duplica la vista de forma simétrica respecto al eje del
quedamos. A continuación representa la fracción de masa de ch4 en el quemador. Main Menu >
Display > Contours… Seleccionamos en Contours Of la opción Species y debajo Mass Fraction of
ch4. Dejamos sin seleccionar ninguna Surface y desactivamos la opción Filled. Eleve el núnero de
Levels a 50. Hacemos clic en Display.
El chorro penetra con gran impulso y poco a poco se va abriendo hasta que se dispersa por
completo dentro de la cámara. La concentración de o2 y n2 es complementaria a esta visualización
(el alumno puede generar la distribución y comprobarlo por sí mismo).
Finalmente, para entender mejor el proceso de dispersión del chorro de metano, se van a
crear una serie de planos transversales, progresivos, sobre los cuales se representarán perfiles de
velocidad longitudinal y concentraciones de metano. Para crear esos planos vamos a:
Main Menu > Surface > Iso-Surface… Seleccionamos en Surface of Constant la opción
Grid y debajo X-Coordinate. Fijamos el valor 0.9 en la casilla Iso-Values y la nombramos como
x0.9 en New Surface Name. Para crear la línea hacemos clic en Apply.
Repetimos el proceso para crear líneas en x=0.45, x=0.225, x=0.1125, x=0.05625 y
x=0.028125.
Una vez finalizada la generación de esos planos, pasamos a dibujar los perfiles. Para ello
crearemos una serie de gráficas mediante el comando XYPlot.
Main Menu > Plot > XY Plot…
Comenzamos activando la opción Position on Y Axis dentro del campo Options. Vamos a
fijar a X con el valor 0 mientras que Y ha de tener el valor 1, pues vamos a representar un gráfico
líneas verticales. Estos cambios en el recuadro Plot Direction harán que FLUENT dibuje la
coordenada y en la ordenada del gráfico.
Bajo X Axis Function seleccione la opción Velocity y después dentro de la subcategoría
elija Velocity Magnitude.
Es importante recordar que X Axis Function e Y Axis Function describen el eje x y el eje y
del gráfico respectivamente y no las direcciones de los ejes de coordenadas del dominio.
Finalmente seleccione outlet y las 6 líneas que acaba de crear bajo la categoría Surfaces
puesto que efectivamente lo que se quiere es que se de el perfil de velocidades a lo largo de dichas
líneas. Con esto se completa la preparación de los parámetros de dibujo.
Por otro lado, para observar mejor los perfiles de velocidad, vamos a utilizar una escala
logarítmica el eje X del gráfico. En el menú Solution XY Plot hacemos clic en Axes… y activando
23
la casilla X, seleccionamos la casilla Log. Podemos activar a su vez las grid lines para que sea más
fácil determinar el límite de la capa límite. Para ello, hacemos clic en Major Rules y Minor Rules
en Options. Luego hacemos clic en Apply. (NOTA: Para activar las grid lines del eje Y debe
seleccionar luego el eje Y, activar las dos opciones de Major Rules y Minor Rules y clicar en
Apply.
Por último, para mejorar la visualización de los perfiles, haga clic en Curves y seleccione
que pinte líneas en el Pattern de cada Curve. Cambie también los estilos de los marcadores
(Marker Style) y hágalos un poco más pequeños (Size=0.2). En definitiva:
¿Qué comentarios podría hacer sobre estos perfiles? ¿a qué es debido ese repunte de la
velocidad en la curva amarilla (x=0.9) cerca de la pared superior?
Repita el gráfico pero mostrando ahora la fracción másica de ch4…Simplemente, cambie en

X Axis Function a la opción Species y después dentro de la subcategoría elija Mass fraction of
ch4. Para apreciar mejor los resultados, edite de nuevo el eje X y deseleccione la opción
AutoRange. Ahora, bajo Range están accesibles los máximos y mínimos de x. Mantenga el valor
máximo de x en 1 y el mínimo en 1e-06.
Explique brevemente estos resultados.

24
6.2.- SOLUCIÓN CON LLAMA A CP CONSTANTE.
Resulta instructivo realizar primero los cálculos con calor específico constante. Para ello, las
opciones que se habían especificado al principio de la práctica son válidas, pudiendo proceder
directamente al cálculo. Esto es, vaya a Main Menu > Define > Models > Species >
Transport&Reaction… y habilite la opción Volumetric Reactions. Así mismo, compruebe que en
CP tenemos seleccionada la opción constant. Introduzca el valor de 1006.4 J/KgK, propio del aire.
Haga clic en Ok para validar. Por otro lado, es necesaria la modificación de los parámetros de
inicialización de modo que partamos de una condición inicial apropiada para el caso con reacción
química. Para ello:
Nueva inicialización de la solución.
Main Menu > Solve > Initialize > Initialize… En el subpanel Compute From
seleccionamos la opción All Zones. Además, vamos a cambiar unos parámetros: Fije una
Temperature inicial de 2000K y una fracción másica de metano, ch4, de 0.2. Haga clic en Init.
La elección de un valor inicial de 2000K para la temperatura proporciona una “chispa
numérica” que permite el comienzo de la reacción química entre las corrientes de metano y aire.
Fijar factores de sobrerelajación.
Al trabajar con problemas que involucren reacción química, resulta conveniente disminuir
los factores de sobrerelajación que posee por defecto FLUENT para las especies químicas. Acceda
a Main Menu > Solve > Controls > Solution… y disminuya los factores asociados a todas las
especies químicas a un valor de 0.9 en el cuadro de nombre Under-Relaxation Factors. El resto de
parámetros los dejamos con sus valores por defecto y presionamos Ok.
para empezar y que nos vaya actualizando los resultados cada 10 de ellas. Hacemos clic en Iterate.
A las 420 iteraciones, el caso converge.
quemador-caliente.cas que podremos recuperar cuando queramos.
Resultados de la solución con llama con CP constante.
Si como en el caso sin llama visualizamos las líneas de corriente del flujo, veremos:
25
Ahora, como resultado de la formación de la llama, la zona de recirculación ha desaparecido

completamente y todas las partículas fluidas de ambas corrientes se dirigen directamente hacia la
sección de salida.
La visualización del campo de isotermas en el quemador nos proporciona información
valiosa sobre la llama. Para dibujarlas, hacemos un Display > Contours de Static Temperature:
Notemos que la temperatura máxima en la cámara es ahora de unos 3000K, valor que resulta
excesivamente alto. El motivo de esta sobreestimación es que estamos despreciando tanto los
efectos de la variación del calor específico de las especies con la temperatura como la radiación
(desactivada en el modelo).
Es extremadamente interesante representar las curvas de nivel de concentración de las
distintas especies, especialmente de aquellas involucradas en la reacción química: CH4, O2, CO2 y
H2O. Podemos conocer el índice de inquemados y la distribución de los gases de combustión a la
salida del quemador.
Compruebe cómo en este modelo sencillo, estequiométrico, esas distribuciones coinciden
con la distribución de la llama.
El calor generado en la combustión se puede calcular como el flujo neto de calor que sale
por la superficie de salida: Main Menu > Report > Fluxes … Elegimos la opción de Total Heat
Transfer Rate a través de la boundary outlet y hacemos clic en Compute: 29.7 kW. Como es un
modelo axisimétrico, FLUENT ya tiene en cuenta que el valor numérico que da el cálculo en esa
sección de salida hay que integrarlo para 2π. Compare ese valor con el estimado para la reacción de
combustión del metano, si su poder calorífico es de 9350 kcal/m3.
26
6.3.- SOLUCIÓN CON LLAMA A CP VARIABLE.
Para el rango de temperaturas que se alcanzan en la llama, proponer un modelo de calor

específico constante no es realista, de modo que hemos de comparar los resultados de calor
específico constante de la sección anterior con otros más realistas que incluyan la variación del
calor específico con la temperatura.
Podemos seleccionar el modelo de calor específico en el panel Main Menu > Define >
Materials… y seleccionamos el modelo mixing-law en la barra CP (j/kg-k). Este modelo realiza un
promedio local, con la fracción másica de cada componente de la mezcla gaseosa como peso.
Presione Change/Create para activas los cambios. Ahora hemos de especificar la dependencia con
la temperatura del calor específico de cada espeficie. Para ello accedemos a la barra Material Type,
eligiendo Fluid. Entonces aparece un panel que nos permite acceder a las propiedades
termodinámicas de cada especie de la mezcla por separado. Seleccionamos co2 en la barra Fluid
Materials. En la casilla CP escogemos piecewise-polynomial, con lo que se abre el panel Piecewise
Polynomial Profile. Presionamos Ok para aceptar los valores por defecto para los coeficientes del
polinomio (FLUENT incorpora una base de datos con las magnitudes termodinámicas de las
sustancias más comunes). Para activar los cambios en las propiedades del CO2, presionamos
Change/Create. Los últimos pasos hemos de hacerlos con el resto de especies químicas presentes
(CH4, N2, O2 y H2O) recordando presionar el botón Change/Create para cada especie química.
Finalmente podemos presionar Close para cerrar el panel.
Iniciando el proceso de convergencia.
Seguimos la misma estrategia de inicialización que en el caso anterior con llama a CP

constante. Por tanto: Main Menu > Solve > Initialize > Initialize… En el subpanel Compute
From seleccionamos la opción All Zones. Además, vamos a cambiar unos parámetros: Fije una
Temperature inicial de 2000K y una fracción másica de metano, ch4, de 0.2. Haga clic en Init.
Se mantendrán los factores de sobrerelajación de 0.9 para las especies. Finalmente, iteramos:
para empezar y que nos vaya actualizando los resultados cada 10 de ellas. Hacemos clic en Iterate.
A las 400 iteraciones, el caso converge.
quemador-cp.cas que podremos recuperar cuando queramos.
Resultados de la solución con llama con CP constante.
Procedemos como en el caso anterior a la visualización de las líneas de corriente, las

isotermas y las curvas de nivel de concentración de especies.
27
La solución fluida no cambia con respecto al caso anterior. ¿Por qué? ¿es esperable este
comportamiento? Vamos a ver qué ocurre ahora con la distribución de temperaturas:
Nótese que el pico de temperatura ha bajado hasta un valor aproximado de 2300K, respecto
a los 3000K en el caso anterior. Además, la distribución de la llama ha cambiado y la zona “fría”
del centro de la llama ha aumentado notablemente respecto al caso anterior.
Para comprobar hasta qué punto varía el calor específico de la mezcla en el interior de la
cámara es interesante obtener una gráfica que muestre sus curvas de nivel. Para ello:
Main Menu > Display > Contours of… y seleccionamos Properties – Specific Heat (CP).
Presionamos Display:
Efectivamente, se alcanzan valores de CP muy grandes en la proximidad de la tobera de

entrada de metano (prácticamente 3000 J/kgK), mientras que ya en la zona de cola de llama, los
valores retoman el valor próximo para el aire a condiciones normales (1000 j/kgK).
Se deja al alumno su iniciativa para que represente cada una de las especies presentes en la
combustión y observe cómo sus distribuciones en el quemador son complementarias y/o
equivalentes.
Nuevamente, el calor generado en la combustión se puede calcular como el flujo neto de
calor que sale por la superficie de salida: Main Menu > Report > Fluxes … Elegimos la opción de
Total Heat Transfer Rate a través de la boundary outlet y hacemos clic en Compute: 30.7 kW. Ha
aumentado respecto al caso anterior, ¿por qué? ¿es más ajustado a lo esperable teóricamente?
Vamos a mostrar perfiles de concentraciones de metano y de CO2 utilizando la utilidad XY

Plot, tal y como se mostró en la parte 6.1 de este guión:
28
Observe las diferencias con el caso sin llama. ¿Hay inquemados a la salida? ¿Más o menos
que en el caso frío? Veamos ahora la distribución de CO2:
Finalmente, podemos hacer una valoración de las fracciones de cada especie en la sección de
salida…Hacemos Main Menu > Report > Surface Integrals… En Report Type elegimos la opción
Mass-Weighted Average y clicamos en la surface outlet. A continuación, elegimos en Field
Variable la opción Species… Mass fraction of ch4. Hacemos clic en Compute. Repetimos la
secuencia para el resto de especies. Y así obtenemos:
% en masa a
Especie
la salida
CH4 0.06 %
O2 5.33 %
CO2 11.47%
H2O 9.39%
N2 69.65%
29
La suma de todas las especies nos da un valor a la salida de: 95.9%. ¿Por qué no se obtiene
un valor del 100%?
7.- EJERCICIOS PROPUESTOS.
Se plantea al alumno que analice cómo se modifican los diversos parámetros (especialmente
a la salida) al variar las velocidades de entrada del combustible y del comburente.
• Manteniendo la entrada del metano a 80 m/s, modifique la velocidad del aire a 2 m/s y
luego a 0.1 m/s. ¿mejora alguna de esas dos configuraciones el caso analizado en el guión
con 0.5 m/s? ¿Es de esperar una relación de velocidades óptimas en función de la relación
de secciones de entrada de ambos compuestos? (razónelo pensando en la relación aire-
metano que es necesario para la combustión estequiométrica. ¿Qué pasaría si hay excedente
o falta de aire para la combustión?).
• Introduzca al modelo el efecto de la velocidad de la reacción química (Arrhenius rate) para
un caso y analice la nueva solución. ¿Hay cambios sustanciales como para hacer que se
tenga en cuenta? ¿o estuvo bien definida desde el principio la hipótesis de “reacción química
instantánea”?
• Modifique los parámetros de turbulencia a la entrada del metano y del aire…Auméntelos al
15% y luego redúzcalos al 5% ¿Mejora la combustión o empeora y en qué sentido? ¿qué
parámetro definiría como el más adecuado para caracterizar la combustión?
30

Manual Pract

Cargado por

Copyright:

Formatos disponibles

Manual Pract

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Manual Pract

Cargado por

Copyright:

Formatos disponibles

Área de Mecánica de Fluidos CURSO 2006-07

SIMULACIÓN DEL FLUJO EN UNA TUBERÍA

1. DESCRIPCIÓN DEL PROBLEMA

1.- DESCRIPCIÓN DEL PROBLEMA

Considérese un fluido moviéndose a lo largo de una tubería de sección circular constante. El

2.- GENERACIÓN DE LA GEOMETRÍA EN GAMBIT

Estrategia a seguir para crear la geometría.

Crear un directorio de trabajo.

Interfaz del GAMBIT.

La interfaz de GAMBIT se compone de las siguientes partes:

• Barra principal. Observe que el nombre tuberia aparece tras el ID en la barra.

• Controles globales. Controlan la apariencia del modelo, las vistas o el zoom:

(deshacer) que son de gran utilidad mientras se crea la geometría y la malla.

• Panel de descripción de GAMBIT. Este panel contiene la descripción de los botones y

• Ventana de trascripción de GAMBIT. Es la ventana en la que van quedando reflejados los

Selección del SOLVER.

donde r y x son las coordenadas radiales y axiales respectivamente.

Se colocará el origen del sistema de coordenadas en la esquina inferior izquierda del

En la ventana de transcripción, GAMBIT indica que ha creado el vértice “vertex.1”. Los

Como puede ver, en problemas bidimensionales, la coordenada z se deja siempre con el

Operation Toolpad > Global Control > Fit to Window Button

Finalmente haga clic en Apply para crear la cara.

3.- GENERACIÓN DEL MALLADO EN GAMBIT

Mallado de los lados.

4.- CONDICIONES DE CONTORNO Y EXPORTACIÓN DESDE GAMBIT

Crear los tipos de condiciones de contorno.

Con esto se activará el panel de especificación de condiciones de contorno. Para empezar,

Ahora seleccionamos el lado izquierdo como siempre, haciendo MAYÚSCULAS + clic

A continuación se repite el proceso con el resto de lados según la siguiente tabla:

Lado (edge) Nombre (name) Tipo (type)

Si se ha hecho todo correctamente, la lista final de condiciones de contorno debería figurar

5.- INICIALIZACIÓN DEL CASO EN FLUENT. CASO LAMINAR

Ejecutar el programa FLUENT.

Main Menu > File > Read > Case…

Validar y dibujar la malla.

A continuación dibuje la malla: Main Menu > Display > Grid …

Definición del SOLVER.

Main Menu > Define > Models > Solver

Main Menu > Define > Models > Viscous

Main Menu > Define > Models > Energy

Definición de las propiedades del fluido.

Main Menu > Define > Materials

Haga clic en Change/Create.

Definición de Condiciones de Operación.

Main Menu > Define > Operating Conditions …

Definición de Condiciones de Contorno.

A continuación se fijarán los valores de velocidad a la entrada y presión a la salida.

Main Menu > Define > Boundary Conditions …

Haga clic en Close para cerrar el menú de Boundary Conditions.

Para la resolución se empleará un esquema de discretización de segundo orden.

Main Menu > Solve > Controls > Solution…

Haga clic en Ok.

Main Menu > Solve > Initialize > Initialize…

Haga clic en Init. Con esto, se completa la inicialización del modelo.

Fijar criterios de convergencia.

Mientras resuelve las ecuaciones de forma iterativa, FLUENT va mostrando en pantalla un

Main Menu > Solve > Monitors > Residual…

Main Menu > File > Write > Case&Data …

Proceso Iterativo hasta obtener convergencia.