Fluent Manejo Del Program

Área de Mecánica de Fluidos. Universidad de Oviedo.
UNIVERSIDAD DE OVIEDO
ÁREA DE MECÁNICA DE FLUIDOS
FLUENT
Manejo del programa y ejemplos para las
prácticas de simulación numérica de
Mecánica de Fluidos
1
DINÁMICA DE FLUIDOS COMPUTACIONAL (CFD)
OBJETIVOS DE ESTAS PRÁCTICAS
Con la realización de las prácticas de simulación numérica de la asignatura Mecánica

de Fluidos, se pretende que los alumnos alcancen los siguientes objetivos.
OBJETIVOS ESPECÍFICOS
- Familiarización con las técnicas de simulación numérica aplicadas a la Mecánica de Fluidos

(técnicas CFD, “Computational Fluid Dynamics”), comprendiendo la importancia de su apli-
cación en la resolución de problemas de esta disciplina.
- Aprendizaje del manejo del programa comercial FLUENT, siendo el alumno capaz de resol-
ver problemas de flujos sencillos con geometrías bidimensionales.
- Exposición de forma adecuada de los resultados de una simulación numérica, interpretándo-

los correctamente mediante la aplicación de criterios científicos y técnicos.
OBJETIVOS ESPECÍFICOS
- Mejora de la confección y presentación de informes científico-técnicos.
- Realización de trabajos en equipo.
2
INTRODUCCIÓN
La dinámica de fluidos computacional o CFD es el arte, más que la técnica, que intenta utilizar los or-
denadores para la simulación del movimiento de los fluidos y, en ocasiones, de otros fenómenos asociados:
transferencia de calor, reacciones químicas, arrastre de sólidos, etc.
Los paquetes de CFD existentes en el mercado son lo suficientemente potentes y fáciles de utilizar como
para que resulte rentable su utilización a nivel industrial. Sus beneficios vienen principalmente de la reducción del
número de ensayos experimentales necesarios y del tiempo de desarrollo.
También pueden proporcionar bastante información complementaria del comportamiento detallado, que
resulta muy difícil conocer experimentalmente. Un valor añadido es poder poner en el producto la etiqueta de
"Diseñado con ayuda del ordenador", y la facilidad para generar dibujos espectaculares, que estimulan la "com-
pra" del producto.
Con la proliferación de programas comerciales, un número creciente de técnicos ha entrado en contacto
con estos métodos. Sin embargo, frecuentemente no se conocen bien las características que tiene el CFD, y por ello,
los resultados a los que se llegan pueden no ser correctos, ni útiles. Por ello, se ha hecho muy importante para el
manejo de CFD, la formación en dinámica de fluidos y el conocimiento de la filosofía, capacidades y limitaciones
del sistema.
Utilidad
El software de CFD busca el cálculo detallado del movimiento de los fluidos por medio de la utilización
del ordenador para la resolución de las ecuaciones matemáticas que expresan las leyes por las que se rigen.
En los resultados de estas técnicas, junto con el movimiento y la presión, pueden obtenerse las variacio-
nes de las propiedades, las fuerzas que ejercen sobre los sólidos adyacentes, los intercambios de energía, etc.
Algunos de los campos de aplicación son:
- Aerodinámica de vehículos terrestres, aviones, entrada y salida de la atmósfera de vehículos espa-
ciales.
- Diseño de motores de combustión, calderas, turbomáquinas (bombas, ventiladores, compresores,
turbinas, etc.)
- Refrigeración de equipos eléctricos y electrónicos.
- Equipos para procesos físicos y químicos: reactores, sedimentadores, mezcladores, intercambiado-
res, eyectores, etc.
- Desarrollo de sistemas de ventilación, calefacción y aire acondicionado.
- Predicción meteorológica.
- Influencia del viento sobre edificios, puentes, etc.
- Dispersión de contaminantes en la atmósfera, ríos y mares.
- Biomedicina: respiración y circulación de la sangre.
- Hidrología y oceanografía: corrientes en ríos, estuarios y océanos.
- Hidrodinámica de buques.
En general, son casos en los que resulta necesario entrar en los detalles de la dinámica de fluidos para
obtener los resultados que se buscan.
Las ecuaciones que definen en cualquier punto del espacio la velocidad y presión de un fluido fueron des-
cubiertas hace más de siglo y medio por el ingeniero francés Claude Navier y el matemático irlandés George Sto-
kes. Estas ecuaciones se derivan directamente de las leyes del movimiento de Newton, y son ecuaciones diferencia-
les en derivadas parciales. Las ecuaciones son las mismas para cualquier situación de flujo de fluidos. La particula-
rización a los casos concretos viene definida por las condiciones de contorno y los valores iniciales. Por ejemplo,
los cuerpos sólidos sumergidos en una corriente, o actuando como contornos exteriores van a definir por exclusión
3
la zona del espacio en la que hay que resolver las ecuaciones, y van a introducir condiciones de contorno en la va-
riable velocidad (velocidad nula en las paredes).
Estas ecuaciones son lo suficientemente complicadas como para que su solución analítica sólo sea posible
en casos muy elementales. La utilización del ordenador para su resolución numérica es lo que ha dado origen a la
dinámica de fluidos computacional o CFD.
Incluso hoy día, la complejidad del cálculo y las limitaciones de los superordenadores más potentes, hacen
que sea absurdo intentar utilizar las técnicas de CFD en los casos en los que otras técnicas han logrado simplifica-
ciones adecuadas, como el cálculo de pérdidas de carga en tuberías y canales, golpe de ariete, diseño de sistemas
oleohidráulicos y neumáticos, etc. También resultan impracticables para procesos muy extensos, como puede ser la
simulación global de una planta química, una depuradora, o incluso un motor de explosión o una turbina de gas en
su conjunto (sí se pueden estudiar por partes).
Algunas ventajas e inconvenientes

Hasta el final de los 60, los ordenadores no alcanzaron velocidades de cálculo suficientes para resolver
casos sencillos, como el flujo laminar alrededor de un obstáculo. Antes de eso, la experimentación constituía el
medio básico de cálculo y desarrollo.
En la actualidad, los ensayos experimentales siguen siendo necesarios para la comprobación de diseños no
excesivamente complejos. Los continuos avances en los ordenadores y algoritmos, permiten una reducción impor-
tante en el número de ensayos necesarios. El diseño típico de un modelo de ala de avión, se hace ahora con 3 ó 4
ensayos en túnel aerodinámico, en vez de los 10 ó 15 que eran necesarios anteriormente.
Las ventajas que proporciona el análisis por CFD se pueden resumir en:
- Reducción sustancial de tiempos y costes en los nuevos diseños.
- Posibilidad de analizar sistemas o condiciones muy difíciles de simular experimentalmente: veloci-
dades hipersónicas, temperaturas muy altas o bajas, movimientos relativos, etc.
- Capacidad de estudiar sistemas bajo condiciones peligrosas o más allá de sus condiciones límite de
funcionamiento, por ejemplo accidentes.
- Nivel de detalle prácticamente ilimitado. Los métodos experimentales son tanto más caros cuanto
mayor es el número de puntos de medida, mientras que los códigos de CFD pueden generar un gran
volumen de resultados sin coste añadido, y es muy fácil hacer estudios paramétricos.
Las técnicas de CFD no son baratas. En primer lugar, son necesarias máquinas de gran capacidad de
cálculo (los investigadores de CFD son usuarios habituales de los ordenadores más potentes que existen), y un
software con precio todavía no accesible al gran público.
En segundo lugar, se necesita personal cualificado que sea capaz de hacer funcionar los programas y
analizar adecuadamente los resultados.
Los desarrollos en el campo del CFD se están acercando cada vez más a los de otras herramientas de
CAE como las de análisis de esfuerzos en sólidos y estructuras. El motivo del retraso es la gran complejidad de
las ecuaciones y, sobre todo, la dificultad de modelizar adecuadamente la turbulencia.
El mayor inconveniente del CFD consiste en que no siempre es posible llegar a obtener resultados sufi-
cientemente precisos, y la facilidad de cometer graves errores de bulto. Esto proviene de:
- La necesidad de simplificar el fenómeno a estudiar para que el hardware y software sea capaz de
tratarlo. El resultado será tanto más preciso cuanto más adecuadas hayan sido las hipótesis y simpli-
ficaciones realizadas.
- La limitación de los modelos existentes para la turbulencia, flujo bifásico, combustión, etc.
Explicación resumida del método

Considérese, por ejemplo, el flujo del aire alrededor de un coche. En teoría, con las ecuaciones de Na-
vier-Stokes, se puede calcular la velocidad y la presión del aire en cualquier punto (en un caso de flujo compre-
sible o en el que haya transferencia de calor, también son variables a considerar la temperatura y la densidad).
Esto permite calcular la resistencia aerodinámica, la adherencia al suelo de los alerones o faldones delanteros, la
adecuada colocación de las tomas de aire, etc.
4
Junto con las ecuaciones, hay que utilizar las condiciones iniciales y de contorno referentes a las variables
y las referentes a la superficie sólida. En este caso, las condiciones referentes a las variables, vienen definidas por la
velocidad del coche y la presión, igual a la atmosférica en los puntos suficientemente alejados. Las condiciones de
las superficies sólidas vienen definidas por la forma, expresada matemáticamente en las coordenadas adecuadas.
Para resolver las ecuaciones, el programa transformará las ecuaciones diferenciales en ecuaciones alge-
braicas, y las resolverá únicamente en un número finito de puntos del espacio. Por eso, lo primero que hay que
hacer es representar el vehículo o, mas bien, su entorno mediante una malla de cálculo; cuanto mayor sea el número
de puntos de esta malla, mayor será la precisión y realismo de la simulación, y más difícil de generar y resolver. En
casos con geometría complicada, esta fase puede ocupar días e incluso semanas.
RESOLUCIÓN NUMÉRICA DE LAS ECUACIONES DE FLUJO
Ecuaciones que describen el movimiento de un fluido

Las ecuaciones que describen un fluido en movimiento se pueden deducir de la ley de conservación de
la masa, y de la ley de conservación de la cantidad de movimiento. Aplicando estas leyes de conservación a un
elemento de fluido, se obtienen las ecuaciones de Navier-Stokes. Simplificadas para un fluido incompresible
son:
r
Continuidad: ∇⋅v = 0
r
dv r r
Cantidad de movimiento: ρ = −∇p + ρg + µ∇ 2 v
dt
Donde ρ es la densidad, y µ es la viscosidad del fluido. Resolviendo estas ecuaciones, se halla la velo-
cidad del fluido, y la presión en cualquier punto del dominio estudiado.
Como se puede ver, son ecuaciones diferenciales en derivadas parciales que sólo en problemas con
geometría y condiciones de contorno sencillas pueden ser resueltas de forma analítica.
Incluso en estos casos, la solución que se obtiene de manera analítica es la del flujo laminar. Por des-
gracia, habitualmente hay que contar con otro factor: la turbulencia, que hace que el flujo real se aparte bastante
de la solución laminar.
Para poder obtener una solución en un caso turbulento complejo, estas ecuaciones han de ser resueltas
por métodos numéricos, con ayuda del ordenador y utilizando modelos de turbulencia, que complican más, si
cabe, el proceso.
Modelos de turbulencia
En la práctica, casi todos los flujos que interesan a científicos e ingenieros son turbulentos.
Para entender en qué consiste la turbulencia pueden resultar ilustrativos unos ejemplos corrientes. Al
abrir un grifo de cocina sólo un poco y el agua que fluye desde el grifo lo hace suave y limpiamente. Este flujo
se llama laminar. Si se abre algo más, el chorro deja de ser transparente y tiene una apariencia... ¡turbulenta!. Lo
mismo se puede ver en el humo de un cigarrillo encendido en el aire en calma. Al salir del cigarrillo el flujo es
laminar. Algo más arriba se transforma en ondulado y difuso.
La turbulencia se podría definir como un fenómeno de inestabilidad intrínseca del flujo. Provoca que el
fluido pase a comportarse de forma aparentemente caótica. De una manera descriptiva podría hablarse de la
formación de torbellinos más o menos aleatorios alrededor de la dirección media del movimiento. Surge cuando
la velocidad del fluido supera un umbral específico, por debajo del cual las fuerzas viscosas amortiguan el com-
portamiento caótico.
El medio que se suele utilizar para determinar si un flujo va a ser turbulento o no es el número de Rey-
nolds. Este número indica la relación, o la importancia relativa, entre las fuerzas inerciales y las fuerzas viscosas
de la corriente. Fuerzas inerciales grandes en relación con las viscosas tienden a favorecer la turbulencia, mien-
tras que, una viscosidad alta, la evita. Dicho de otro modo, la turbulencia aparece cuando el número de Reynolds
supera cierto valor.
5
La turbulencia no es un fenómeno siempre indeseable. En los cilindros de un motor de combustión in-

terna, por ejemplo, la turbulencia mejora el mezclado del combustible y el comburente y produce una combus-
tión más limpia y eficiente. La capa límite laminar es más resistente al desprendimiento que la turbulenta; por
eso en los aviones de aeromodelismo se busca conseguir que el flujo sobre las alas sea turbulento de forma que
puedan volar en circunstancias más extremas sin entrar en pérdida (en aviones de mayor tamaño no existe este
problema porque el flujo es siempre turbulento). Los hoyos de las pelotas de golf hacen pasar la capa límite a
turbulenta debido a que al desprenderse más tarde que la laminar, reducen el tamaño de la estela, disminuyendo
sustancialmente el arrastre. Con esto se consigue lanzar la bola a más del doble de distancia con la misma fuerza.
Los torbellinos, también llamados vórtices, que caracterizan el flujo turbulento son de muy diversos ta-
maños. Estos vórtices se forman y deshacen sin solución de continuidad. Los torbellinos grandes se rompen en
otros menores, éstos en torbellinos más pequeños, y así sucesivamente. Cuando los torbellinos se hacen lo sufi-
cientemente pequeños se disipan en forma de calor debido a la viscosidad.
Las ecuaciones de Navier-Stokes son una representación matemática adecuada de los flujos de fluidos,
incluso turbulentos. Sin embargo, la resolución de estas ecuaciones de forma que se calcule el flujo turbulento
requiere una discretización temporal y de la malla de cálculo tan detallada que en la práctica resulta inviable. En
un artículo publicado por P. Moin y J. Kim (“Simulación de la turbulencia mediante superordenadores”, Investi-
gación y Ciencia, marzo, 1997) se pone el siguiente ejemplo:
La resolución de las ecuaciones de Navier-Stokes para el flujo alrededor de un avión requiere una ma-
lla de cálculo poco espaciada con el fin de discernir los torbellinos más pequeños. Por otra parte, la malla ha
de ser lo bastante grande como para abarcar el avión entero y parte de su entorno. La disparidad de escalas de
longitud en un flujo turbulento -relación entre los tamaños mayor y menor de los torbellinos- puede calcularse
elevando a ¾ el número de Reynolds del flujo. Podemos emplear esta relación para estimar el número de puntos
de malla exigidos en una simulación razonablemente precisa: debido a la existencia de tres dimensiones, el
número es proporcional al cubo de esta relación de escalas de longitud. Por consiguiente, el número de puntos
requerido para una simulación numérica será proporcional al número de Reynolds elevado a 9/4. Con otras
palabras, duplicar el número de Reynolds conduce a multiplicar por cinco, aproximadamente, el número de
puntos en la malla requeridos para simular el flujo.
Consideremos un avión de transporte con un fuselaje de 50 metros y con alas de unos 5 metros de
cuerda (distancia entre los bordes de ataque y de salida). Si la aeronave se encuentra en vuelo de crucero con
velocidad de 250 metros por segundo y altitud de 10.000 metros, se necesitan unos 10.000 billones (1016) de
puntos para simular, con una aproximación razonable, la turbulencia cerca de la superficie.
¿Qué exigencias de cálculo impone tal número de puntos? Una estimación grosera, basada en los ac-
tuales algoritmos y software, sugiere que incluso un superordenador capaz de efectuar un billón (1012) de ope-
raciones en coma flotante por segundo ¡tardaría varios miles de años en calcular el flujo para un segundo de
tiempo de vuelo! Y tal ordenador no existe todavía.
Por fortuna, no se necesita simular el flujo de esta manera para obtener información útil. Incluso, aun-
que fuera factible, se generarían más datos de los que se pueden manejar. Habitualmente lo que interesa son los
efectos de la turbulencia sobre los valores medios de las variables: la velocidad media y la presión media en el
caso del flujo en un conducto; en el caso de un avión, las fuerzas medias de resistencia y sustentación; para el
caso de un motor, los efectos de la turbulencia sobre las relaciones de mezcla entre combustible y comburente;
etc.
Para conseguir esto, las ecuaciones de Navier-Stokes se promedian sobre las escalas de las fluctuaciones
de turbulencia. En la práctica, esto significa que no suele calcularse el movimiento de todos y cada uno de los
pequeños torbellinos. Antes bien, se calculan los torbellinos grandes y se utilizan modelos de turbulencia para
estimar los efectos de los pequeños torbellinos sobre los grandes. Estos métodos dan lugar a un campo de flujo
promediado y simulado que es más uniforme que el flujo real, y, por tanto, reduce drásticamente el número de
puntos de malla y la discretización temporal necesaria para simular el campo.
Los modelos básicos de turbulencia que se utilizan actualmente varían en complejidad desde simples
coeficientes de viscosidad ajustados hasta sistemas completos de ecuaciones adicionales. Algunos de ellos son:
- Modelo de longitud de mezcla.
- Modelo K- Epsilon (con multitud de variantes).
- Modelos de Esfuerzos cortantes de Reynolds.
Todos ellos incluyen coeficientes que hay que ajustar para cada tipo de flujo y que, habitualmente, se
obtienen por vía experimental. Por otra parte, cada modelo tiene sus preferencias: van mejor con flujos abiertos o
6
confinados, con separación o sin ella… En definitiva, las simulaciones de flujos turbulentos promediados poseen
sólo la exactitud de los modelos que incorporan.
Como resumen de estos comentarios sobre la turbulencia, se puede decir que la validez de una simula-
ción numérica depende en gran parte de la modelización de la turbulencia, y que, respecto a ésta, todavía no se
dispone de modelos de resolución perfectos.
DISCRETIZACIÓN DE LAS ECUACIONES DEL FLUJO

Para que un ordenador pueda resolver las ecuaciones que describen el movimiento del fluido, éstas han
de ser transformadas en expresiones algebraicas que sólo contengan números, combinados mediante operaciones
sencillas, tales como sumar, restar y multiplicar.
La transformación de las ecuaciones diferenciales en su análogo numérico, es lo que se llama proceso de
discretización numérica. Hay varias técnicas de discretización, dependiendo de los principios en que se basen. Las
más usadas son: diferencias finitas, volúmenes finitos, y elementos finitos. Aquí se describirá brevemente el método
de los volúmenes finitos, que es el utilizado por el programa FLUENT .
Método de los volúmenes finitos

Este método consiste, fundamentalmente, en convertir las ecuaciones diferenciales en su análogo numé-
rico, mediante una transformación física de las ecuaciones. Por ejemplo, la ecuación de la cantidad de movimien-
to puede ser considerada como una serie de flujos dentro de un volumen de fluido, junto con un término de fuen-
te que es el gradiente de presiones.
Para mostrar la discretización, se crean volúmenes finitos (particiones del espacio) en la dirección x, y se
fija la atención en el volumen n-ésimo por simplicidad. En la figura se puede ver un volumen finito, o celda, típico
donde el centro del volumen, punto P, es el punto de referencia en el que se quiere hallar el análogo numérico de la
ecuación diferencial. Los puntos centrales de las celdas vecinas se denominan W y E (celda oeste y celda este). Las
caras que delimitan la celda a un lado y a otro se denominan w y e (cara oeste y cara este).
w e
W P E
∆y
∆x
Supóngase que se quiere discretizar la derivada segunda siguiente:
∂2U
∂x 2
Hay que hacerlo a partir de las derivadas primeras en las caras que limitan la celda:
 ∂U ∂U 
 − 
∂ U
2
 ∂x e ∂x w 
 2  =
 ∂x  P xe − xw
Estas derivadas se calculan, a su vez, a partir de los valores de la variable U en los centro de las celdas ve-
cinas:
 ∂U  U − UP
  = E
 ∂x  e x E − x P
 ∂U  U − UW
  = P
 ∂x  w x P − x W
Estas expresiones son utilizadas para implementar las derivadas de las ecuaciones en una celda cual-
quiera. De esta manera, se puede hallar el análogo numérico de las ecuaciones diferenciales para cada volumen
7
finito, y después resolver el sistema de ecuaciones algebraicas que resulta. De la misma manera, pueden introdu-
cirse las condiciones de contorno, añadiendo términos fuente en las ecuaciones de las celdas correspondientes.
Resolución de las ecuaciones discretizadas

Una vez discretizadas las ecuaciones diferenciales, lo que se obtiene es un sistema de ecuaciones alge-
braicas, que se puede escribir en notación matricial de la siguiente manera:
A es la matriz del sistema, b es la columna de términos independientes, y x es el vector incógnita. Para
resolver este sistema de ecuaciones, existen dos tipos de métodos: los directos y los iterativos.
Los métodos directos consisten en hallar la inversa de la matriz A. Normalmente estos métodos parten
A∗x = b
de la descomposición de la matriz A en el producto de dos matrices: L*U, una triangular superior y la otra es
triangular inferior, a partir de lo cual todo es muy sencillo. Cuando la matriz es un poco grande estos métodos
no son rentables, computacionalmente hablando.
Los métodos iterativos tratan de hallar la solución partiendo de una solución aproximada y calculando, a
partir de ella, una solución más aproximada, acercándose en cada iteración a la solución real. De este tipo son los
métodos de Jacobi, Gauss-Seidel, TDMA (algoritmo de la matriz tridiagonal), SIP (procedimiento fuertemente
implícito) de Stone, CGM (método del gradiente conjugado), etc. Cada uno de ellos con sus adeptos y variantes,
lo que normalmente quiere decir que van bien en unos casos y no tan bien en otros.
Dependencia del tiempo y no linealidad

Los métodos que se han comentado resuelven sistemas de ecuaciones lineales, con coeficientes constan-
tes, y no se pueden usar directamente para resolver las ecuaciones del flujo, pues no son lineales (hay en ellas
productos de variables). Además, la velocidad y la presión del fluido pueden depender del tiempo, por lo que hay
que resolver todo el sistema para cada instante de tiempo. Por tanto, habrá que discretizar el tiempo, como ya se
hizo con el resto de variables, y resolver el sistema para cada paso temporal.
Para que sea posible emplear los métodos de resolución de sistemas que se han visto, se linealizan las
ecuaciones. Los términos no lineales, como hemos visto, son los productos de las velocidades por sus derivadas.
Por ello, se discretiza la derivada, y se usa el valor de la velocidad existente en la celda de que se trate:
∂u u − u i-1, j
u se convierte en u i +1, j
∂x 2δx
Así se puede resolver iterativamente el sistema lineal de ecuaciones. Una vez resuelto, el valor de la ve-
locidad en la celda i,j habrá cambiado, se vuelve a linealizar el sistema con los nuevos valores de las velocida-
des, y se repite el proceso hasta que se llegue a una solución.
Posteriormente, habrá que incrementar el paso temporal, y resolver el sistema para el nuevo instante de
tiempo. En el caso de que sea un problema estacionario, este último proceso de cambiar el paso temporal, se
omite.
Condiciones de contorno
La solución de un problema depende de las condiciones iniciales y de las condiciones de contorno a las
que está sometido el fluido. Las condiciones de contorno más habituales que suelen aparecer en los problemas
resueltos con CFD son las siguientes:
- Presión fija en un determinado contorno. Se puede fijar presión estática o presión total.
- Velocidad fija en el contorno. Se puede fijar dirección y módulo, o sólo dirección, combinada con la
condición de presión.
Imponiendo estas condiciones, se puede simular el comportamiento de los contornos reales de un fluido,
tales como paredes, entradas en tuberías, salidas al exterior, etc.
Las paredes sólidas son condiciones de contorno especialmente críticas. Esto es debido al comporta-
miento del fluido en las proximidades de la misma. En las superficies sólidas se genera una capa límite de fluido
en la que la velocidad varía rápidamente a medida que nos acercamos a la pared. La simulación directa de la
8
capa límite es extremadamente difícil y requiere unos recursos desproporcionados, por lo que se suelen utilizar
modelos específicos.
Metodología
En los paquetes de software existentes de CFD, el usuario debe especificar las condiciones del problema
que se va a resolver, así como proveer al ordenador de ciertos parámetros de resolución para que el programa
consiga encontrar una correcta solución del problema. Los pasos en los que esta información es introducida en el
ordenador son los siguientes:
1. Especificación de la geometría del problema. Es preciso introducir al programa los contornos del
dominio del problema. Este proceso es relativamente sencillo cuando el modelo es bidimensional,
porque se trata de definir unas líneas, que limitan la zona del problema. Estas líneas pueden ser
muy sencillas en el caso de rectas o circunferencias. Si las curvas son de otro tipo, se introducirán
mediante puntos. Cuando el modelo es tridimensional, la especificación de la geometría requiere
introducir las superficies de contorno del problema, que puede tener mucha complicación si hubiera
que definirlas por puntos.
2. Creación del mallado, o volúmenes en los que van a ser calculadas todas las variables. El mallado
puede ser estructurado o no estructurado. El mallado estructurado está formado por rectángulos o
paralelepípedos, que pueden deformarse para adaptarse a las superficies curvas de la geometría. El
mallado estructurado está formado por triángulos o tetraedros, construidos a partir de la división de
los contornos de la geometría. Este último tipo de mallado es mucho más sencillo de crear que el
anterior.
3. Imposición de las condiciones de contorno, que fijan los valores de ciertas variables en los límites
del dominio. Estas condiciones se le imponen a las variables de la presión estática, presión total,
velocidad, o gradiente de las variables. También hay que imponer condiciones para las variables de
turbulencia que haya, según el modelo de turbulencia que se esté utilizando.
4. Introducción de las condiciones iniciales. En el caso de un problema no estacionario, hace falta de-
finir unos valores de las variables para todos los puntos del dominio, desde los cuales empezará el
programa a calcular las soluciones de los sucesivos pasos temporales. En el caso de un problema
estacionario, es preciso introducirle al programa unos valores iniciales de las variables, de los cua-
les partirá el proceso iterativo.
5. Especificación de las propiedades del fluido, tales como la viscosidad y la densidad. Estas magni-
tudes pueden variar con la temperatura, la presión, etc. En ese caso habría que introducirle al pro-
grama el modelo matemático de variación, con los parámetros correspondientes. También pueden
considerarse varios tipos de problema, según se tenga en cuenta o no la transferencia de calor. En
ese caso, habría que introducirle al programa las propiedades térmicas del fluido de que se trate.
6. Control de los parámetros que afectan a la resolución numérica del problema. Como ya se explicó
antes, se puede actuar sobre el proceso iterativo variando los factores de relajación, y fijando los
criterios de finalización de las iteraciones.
9
AL EMPEZAR CON EL SOFTWARE
Para realizar una simulación numérica de un caso de dinámica de fluidos con un paquete comercial de
CFD, básicamente hay que seguir los siguientes pasos:
- Definir la geometría de la zona que se quiere calcular. Esta zona corresponde con la ocupada por el
fluido. Normalmente, los sólidos que rodean o están contenidos por el fluido actúan como contor-
nos.
- Generar el mallado de la zona(s) de fluido. Donde la geometría requiera más detalle o se esperen
mayores gradientes en las variables, el mallado deberá ser más fino, para que el flujo se calcule en
un mayor número de puntos.
- Seleccionar los modelos que hay que utilizar para el cálculo: flujo bi o tridimensional, compresible
o incompresible, modelo de turbulencia, tipo de fluido, etc.
- Definir las condiciones de contorno en los bordes que rodean a la zona de cálculo.
- Seleccionar los parámetros de resolución para el proceso iterativo.
- Realizar el cálculo realmente dicho.
- Visualizar y analizar la solución obtenida.
El paquete Fluent consta de dos programas que permiten realizar las operaciones anteriores: Gambit y
Fluent 5. Con el primero se define la geometría y se genera el mallado. También hay que especificar en él qué
tipo de condición de contorno tiene cada borde, aunque no se especifican sus parámetros. Con el segundo se
realizan todas las demás operaciones.
Tal como se van a utilizar, el Fluent 5 funciona como un programa normal de Windows NT; el progra-
ma no está instalado en cada ordenador sino que se lee de un servidor. Lo mismo sucede con el Gambit, pero
además, para ejecutarle debe estar instalado en cada ordenador el programa Exceed, pues la salida gráfica del
Gambit está basada en ventanas de Xwindows.
Aparte de esto, en el entorno de NT tienen que estar definidas unas ciertas variables en Panel de con-
trol, Sistema, Entorno, Variables del usuario:
En el caso de que el programa no consiga arrancar es aconsejable comprobar que se tiene acceso al ser-
vidor en Entorno de Red.
FICHEROS QUE SE GENERAN

La salida del mallado desde el programa Gambit al Fluent 5 se realiza a través de un fichero de tipo
.msh, exportando el mallado con el menú:
File
Export
Mesh...
Sin embargo, este fichero no puede volver a ser leído por Gambit. Para poder realizar cambios, guardar
una sesión que no se ha finalizado o como seguridad (no vaya a ser que haya un apagón), conviene guardar el
trabajo con el formato propio de Gambit (que no puede ser utilizado directamente por el programa de resolu-
ción). Esto se realiza en
File
Save as...
Al hacer esto, GAMBIT genera una serie de ficheros que tienen en común el nombre, variando la exten-
sión:
- .dbs:
Fichero básico de datos de Gambit. Es el único fichero realmente necesario. Si se pierde o se dete-
riora este fichero, hay que volver a empezar el trabajo.
10
- .jou:
Listado de todas las instrucciones ejecutadas con el Gambit para el trabajo en concreto.
- .trn:
Listado también de instrucciones, pero más relacionado con operaciones de lectura-escritura.
- .lok:
Cuando se está utilizando Gambit se genera este fichero para bloquear su uso desde otros progra-
mas. Al salir de Gambit se borra automáticamente. En el caso de que se salga del programa de for-
ma no estándar (apagando el ordenador o cerrando Exceed), no se borra y bloquea la utilización
posterior del caso. Para desbloquearlo basta con borrar este fichero. Un problema especialmente
desconcertante surge cuando lo anterior sucede sin haber dado nombre al caso. Gambit utiliza por
defecto un caso llamado default_id y la existencia previa del fichero default_id.lok impide abrir el
programa. La solución, como se ha comentado, es borrar este fichero.
Además de estos ficheros, el programa crea unos subdirectorios llamados Gambit.(número), que se
pueden borrar sin especial remordimiento.
Por defecto, estos ficheros se guardan en el directorio definido como HOME, aunque no es necesario
que estén ahí.
El programa Fluent 5 lee el fichero .msh
File
Read
Case
y genera dos tipos de ficheros de salida: .cas y .dat.
File
Write
Case // una de las tres
Data
Case & Data
En el primero está contenido el mallado, las condiciones de contorno y todas las definiciones que se rea-
lizan antes de ejecutar los cálculos. En el segundo están contenidos exclusivamente los datos numéricos del
resultado. Si se quiere guardar el trabajo para continuar en otro momento o por seguridad, es suficiente con
guardar el caso (.cas). Una vez terminada la ejecución se pueden guardar también los datos (.dat). Si se quiere
examinar un caso terminado con anterioridad es preciso leer los dos ficheros. El fichero de datos puede ocupar
mucho espacio; en caso de saturación, es suficiente con guardar el fichero de caso, pues, a partir de él, se pueden
regenerar los datos sin más que poner el programa a calcular.
Respecto al lugar de almacenamiento, para trabajar es aconsejable usar el disco duro del ordenador que
se esté utilizando (en principio en el subdirectorio \temp), pues así se minimizan los posibles errores de red y se
evita el bloqueo por sobrepasar la capacidad del disquete o cuenta. Los ficheros .dbs, .msh, .cas y, si cabe, .dat
habría que guardarlos en la cuenta personal al finalizar, e incluso hacer una copia de seguridad en disquete.
Estos apuntes están estructurados en dos grandes bloques (GAMBIT y FLUENT) constando cada uno
de ellos de:
a) Explicación genérica de las opciones de cada programa.
b) Resolución de un ejemplo concreto.
Finalmente se añade un caso de aplicación : estudio del flujo en una tobera convergente-divergente.
11
GAMBIT
De cómo utilizar el GAMBIT para generar figuras geométricas, crear el ma-

llado interior y definir las condiciones de contorno. Un ejemplo al detalle.
INTRODUCCIÓN
GAMBIT es el generador de geometría, mallado y condiciones de contorno usado por Fluent. Tiene ca-
pacidad para generar mallados tanto estructurados como no estructurados.
En el caso bidimensional, la geometría debe estar contenida en un plano y estará formada por una o va-
rias superficies cerradas (caras). En el caso tridimensional constará de uno o varios volúmenes.
El mallado consiste en la división de las caras o volúmenes en celdas. Se pueden generar mallados tanto
de triángulos como de cuadriláteros (tetraedros, hexaedros, prismas en 3D), aunque aquí sólo se van a comentar
los mallados bidimensionales triangulares. El programa también tiene capacidad para leer mallados o geometrías
generados por otro software de CAE.
Se maneja íntegramente con el ratón a partir de los menús. La ayuda es bastante buena.
El esquema de un proceso típico incluiría los siguientes pasos:
- Generar la geometría
- Crear puntos
- Unir los puntos con líneas
- Crear la superficie
- Generar el mallado
- Poner nodos en las líneas
- Ejecutar la acción de mallado
- Poner las condiciones de contorno
- Guardar el fichero
Téngase en cuenta que el programa tiene muchas más capacidades de las que se van a explicar, por lo
que no se detallan todas las opciones.
GENERAR LA GEOMETRÍA
Al arrancar el programa (pinchando el icono correspondiente), aparece la pantalla principal del
GAMBIT, con la que se trabajará siempre.
12
En la parte superior se tiene un menu principal con las siguientes opciones:
La opción File permite guardar y recuperar un fichero, así como importar geometrías realizadas con
otros programas y exportar el fichero de mallado, una vez terminado todo el proceso, con la opción Export...
Mesh..., para poder después recuperarlo con el programa de resolución.
La opción Solver permite escoger el programa con el que se va a resolver el caso: FLUENT 5.
Las operaciones realizadas, así como los errores aparecen en la ventana inferior izquierda.
La parte superior derecha muestra las tres operaciones básicas: geometría , mallado y condicio-
nes de contorno :
13
Cada una de estas operaciones se describirá a continuación. En la geometría aparece la entrada de pun-
tos y la creación de líneas, superficies y volúmenes:
Los iconos que llevan una punta de flecha en la parte inferior, poseen más información. Para ver todas
las opciones se pulsa el botón derecho del ratón. Para escoger una opción se pulsa el botón izquierdo.
Cada vez que se pasa el ratón por encima de un icono, en la parte inferior de la pantalla principal se ob-
tiene una descripción de la operación que realiza.
Para introducir los puntos por medio de sus coordenadas geométricas se pulsa el primer icono del apar-
tado Vertex apareciendo (dibujo superior izquierdo).
Se pueden introducir en coordenadas cartesianas, cilíndricas y esféricas cambiando la opción Type. Se
introducen los valores en Global y se le puede asignar un nombre al punto o dejar que se asigne uno automáti-
camente con la denominación vertex.número. Si los datos son correctos se pulsa Apply, si no lo son se pueden
corregir directamente y después se pulsa Apply. Si se ha pulsado Apply y los datos son incorrectos se va al icono
Anular y se borra la última operación. También se puede ir al icono Borrar, se selecciona el punto y se pulsa
Apply.
14
A medida que se van introduciendo los puntos, estos se van mostrando en pantalla. Se puede cambiar la
escala y la orientación usando de los iconos que aparecen en la parte inferior derecha o el ratón:
El icono Encajar en pantalla permita ajustar el dibujo al tamaño de la ventana. El icono Cambiar orien-
tación de los ejes permite colocar el dibujo con respecto a una orientación determinada. Es muy útil para volver-
lo a colocar después de haber mareado la vista con el ratón. El icono visualizar atributos permite ver en el dibujo
los puntos, líneas, mallado, etc. así como los nombres. En general conviene verlo todo salvo que las cosas se
compliquen.
Arrastrando con el botón izquierdo del ratón en la pantalla, se gira la vista con respecto a un eje conte-
nido en el plano del dibujo y perpendicular a la dirección de arrastre. Arrastrando con el botón derecho en hori-
zontal, se gira el dibujo respecto a un eje perpendicular a la pantalla. Arrastrando en vertical se hace zoom. Apre-
tando la tecla Control y seleccionando un área con el botón izquierdo, se aumenta la zona seleccionada.
Una vez dibujados todos los puntos se procede a la creación de líneas. Seleccionando el segundo botón
del apartado Geometry:
En el dibujo anterior izquierdo se observan los tipos de líneas que se pueden crear. Las instrucciones
son bastante sencillas. Los puntos se pueden seleccionar en la pantalla por medio de la combinación Shift +
Botón izquierdo ratón (el punto cambia de color) o en una lista que aparece en la flecha vertical de la opción
Vértices. Cuando un elemento está seleccionado, cambia su color en el dibujo. El resto es similar a lo comentado
para los puntos.
Shift + botón derecho anula la selección. Si se aprieta el botón izquierdo mientras se mantiene el dere-
cho, cambia la funcionalidad de los botones respecto a la tecla Shift, lo cual puede ser un poco desconcertante.
Se reconoce el cambio porque el cursor pasa de ser una flecha a un circulo con un punto en medio. Repitiendo la
operación se deshace el entuerto.
El siguiente paso sería crear la/s superficie/s. Para ello se selecciona el tercer icono de la sección Geo-
metry (dibujo anterior centro)
El icono superior izquierdo del apartado Face contiene los métodos de generación de caras. El más ele-
mental, Wireframe, genera una cara a partir de un conjunto de líneas que definan una superficie cerrada.
Las líneas que forman cada superficie se seleccionan de la misma manera que los puntos para las líneas.
Con estas superficies se pueden realizar operaciones booleanas usando el tercer icono de la sección Fa-
ce (dibujo superior derecho). Permite la unión (Unite), la resta (Subtract) y la intersección (Intersect) con la
opción de mantener también la/s superficie/s originales.
El quinto icono posee la opción de renombrar las superficies, aunque conviene tener en cuenta que se
pierden los nombres cuando se realiza cualquiera de la operaciones anteriores.
15
GENERAR EL MALLADO
Una vez generada la geometría bidimensional, el siguiente paso será realizar el mallado, primero con
todas las líneas y después con las superficies.
Se selecciona el icono de mallado Mesh (el segundo del apartado Operation), a continuación, el icono
de línea (segundo del apartado Mesh), y después el icono de mallar (primero de Edge). La línea se puede selec-
cionar de la lista que aparece pulsando la flecha de la opción Edges o seleccionándola directamente como se ha
comentado anteriormente. La pantalla que aparece es la del dibujo inferior izquierdo.
En la sección Spacing se puede seleccionar el colocar los nodos cada cierta distancia (Interval size) o
bien poner un número fijo de nodos (Interval count). Se mallará con la opción Interval Count. La sección Gra-
ding permite distribuir los nodos de manera no uniforme a lo largo de la línea.
Para hacerlo efectivo se pulsa Apply. En otros iconos de la sección Edge, se pueden anular las últimas
instrucciones, se pueden borrar mallados en líneas y se dispone de información acerca de lo realizado.
Para mallar las superficies, se selecciona el icono de mallado de caras (tercero de Mesh) y después el
icono de mallar (primero de Face). La superficie se puede seleccionar de la lista que aparece pulsando la flecha
de la opción Faces o seleccionándola con Shift + Botón izquierdo en una de las líneas laterales de la cara.
En el apartado Scheme se puede seleccionar el tipo de celdas: cuadrangulares o triangulares. En los
ejercicios propuestos se utilizará únicamente mallados triangulares. La información es la del dibujo inferior
derecho. La cara se mallara de acuerdo con el mallado que poseen todos sus lados. Para ello se desactivará la
opción Spacing – Apply. Es posible mallar directamente las caras sin especificar los nodos de las líneas laterales
– en este caso, el apartado Spacing permite conseguir un mallado más o menos denso – pero habitualmente, con
el mallado anterior de las líneas se obtienen mallados de superficies más adecuados.
ESPECIFICAR TIPOS DE CONDICIONES DE CONTORNO

Si los pasos anteriores han sido correctos, la última tarea con este programa será el indicar los tipos de
condiciones contorno y el medio en que se desenvuelve.
Se especifican primero los tipos de contorno (dibujo inferior izquierdo). Se selecciona el tercer icono de
la sección Operation y, a continuación, el primero de la sección Zones. Si la geometría es bidimensional, las
condiciones de contorno hay que ponerlas en las líneas, por lo que en el apartado Entity hay que seleccionar
Edges (por defecto está Faces, para mallados tridimensionales).
Con el ratón o en la lista emergente (flecha hacia arriba) se escoge la línea – o líneas, en el caso de que
queramos que varias líneas tengan la misma condición de contorno – y en el apartado Type se selecciona el tipo
de condición: wall, velocity inlet, pressure outlet, etc. El apartado Name permite indicar un nombre para el con-
16
torno en lugar del asignado automáticamente. Estas condiciones de contorno pueden ser modificadas posterior-
mente con el programa FLUENT 5.1.
Si en el apartado Action está marcado Add, al pulsar Apply se asigna la condición. También se puede
modificar (Modify) o borrar (Delete) contornos. Por defecto está activado Add.
Se dispone de una ventana con el nombre y el tipo de todos los contornos asignados.
Ya sólo queda indicar el medio, fluido o sólido (dibujo superior derecho). El segundo icono de la sec-
ción Zones abre la ventana que permite elegir los tipos de continuo; habitualmente las zonas malladas en las que
se resuelve el problema deberán ser zonas de fluido. Únicamente se utilizan zonas sólidas cuando se quiere re-
solver un problema asociado de transmisión de calor... En los casos bidimensionales las entidades serán caras
(Faces).
Una vez realizadas estas operaciones, solamente queda escribir en disco el fichero con el mallado gene-
rado, para abrirlo posteriormente con el programa de resolución. Esto se realiza en el menú File:
File
Export
Mesh...
Por defecto la terminación de los ficheros de mallado es msh.
A lo largo del proceso y al acabar, conviene guardar el trabajo con el formato propio de Gambit (que no
puede ser utilizado directamente por el programa de resolución), pues los ficheros msh no pueden ser vueltos a
leer por Gambit. Esto se realiza en
File
Save as...
ALGUNAS PRECISIONES SOBRE GEOMETRÍA Y MALLADO

Quizá convenga hacer unas consideraciones respecto a la filosofía de la generación del mallado con el
Gambit. Lo que realmente utiliza el programa de resolución es el mallado. La geometría únicamente tiene una
función de apoyo: las superficies que no se mallen no se utilizarán (los puntos y líneas no se utilizan en ningún
caso). Incluso los nodos que se colocan en las líneas para realizar el mallado de las superficies tampoco tienen
17
ninguna utilidad ulterior. Las líneas únicamente van a servir para asignar condiciones de contorno pero, para
esto, deben formar parte del borde de las superficies. Una línea aislada no tiene ninguna función.
Por otra parte ha de tenerse presente la idea de que una superficie está definida por una serie de líneas
que deben ser continuas entre si (donde acaba una empieza la siguiente) y formar un conjunto cerrado. Una línea
que corte parcial o totalmente la superficie, realmente no formará parte de ella; no será utilizada en el mallado, ni
en la resolución.
Incorrecto Correcto
Si se quiere dividir una superficie, hay que hacer dos. Esto no quiere decir que el flujo no pase de una a
otra, basta con definir la línea(o líneas) que separa las dos superficies como condición de contorno interior.
También se puede hacer con varias líneas seguidas y definir unos tramos como pared y otros como interior, blo-
queando el paso sólo en parte.
Incorrecto Correcto
Un caso especial son los objetos dentro de una superficie, el procedimiento es crear la superficie global,
crear la superficie objeto y sustraer el objeto de la superficie global, generando así una superficie con un agujero.
En los bordes del objeto (que si es sólido no hace falta mallarlo, ni siquiera definirlo más que para restar) se
pondrán condiciones de contorno. Si no se sustrae el objeto de la superficie global, ésta tendrá una parte super-
puesta con el objeto y el mallado no será coherente o bien no tendrá definido el objeto. El objeto debe ser una
superficie; no es suficiente con que sea una línea.
Incorrecto Correcto
Un punto puede formar parte de varias líneas, así como una línea puede formar parte de varias superfi-
cies. Dos superficies adyacentes pueden estar formadas, en su borde común, por la misma línea o por líneas
distintas (en este caso las dos líneas estarán superpuestas en el dibujo). La única diferencia es que al asignar
condiciones de contorno habrá que asignar la condición en una línea o en dos.
Estas ideas son semejantes en casos tridimensionales sin más que cambiar superficies por volúmenes y
líneas por superficies.
18
ALGUNOS CONSEJOS E INDICACIONES

Hay que procurar mantener los casos lo más sencillos posible e ir probando una sola cosa cada vez. Los
primeros ejemplos deberían tener pocos nodos (y celdas, por consiguiente), de forma que no se tenga que utilizar
mucho tiempo de trabajo del ordenador antes de descubrir que uno se ha equivocado. Un corolario de esta regla
dice que si un caso complicado funciona a la primera, el error es más grave de lo que parece.
Por razones de claridad compensa poner nombres de pocos caracteres a puntos y líneas.
En caso de duda conviene separar en varias las condiciones de contorno. Es más fácil imponer las mis-
mas condiciones a varios tramos en el Fluent que tener que redefinir la geometría. No tiene mucha importancia el
tipo de condición de contorno que se imponga dentro del Gambit; se puede modificar posteriormente.
Un caso de flujo interno es aquel en el que el flujo está delimitado por paredes sólidas, por ejemplo un
tramo de tubería, una tobera, una habitación, etc. Habitualmente las condiciones de contorno serán las paredes,
una o varias entradas y una o varias salidas. Las que dan más problemas son las salidas: no todos los tipos son
compatibles o funcionan correctamente según las condiciones del flujo. También puede haber una línea de sime-
tría, en los casos que tengan un eje principal de simetría; esto reduce a la mitad los cálculos (sólo hay que gene-
rar y calcular media geometría). Únicamente hay que tener cuidado de que, tanto la geometría como el flujo,
sean realmente simétricos.
Los casos de flujo externo son los que tienen zonas sin delimitar, por ejemplo un vehículo, un perfil de
ala, un objeto sólido sumergido en una corriente, etc. En estos hay que definir los límites del mallado encerrando
el objeto con líneas de manera un poco arbitraria. Las condiciones de contorno impuestas en estas zonas son
habitualmente entradas y salidas, debiendo estar lo suficientemente lejos del objeto origen del estudio como para
que las condiciones sean realistas, es decir, para que la presencia del objeto no influya en las velocidades de
entrada y las presiones o en las direcciones de salida. En caso contrario se estaría forzando a la solución y el
resultado no sería realista. Un efecto típico son las estelas: la presencia de un cuerpo se deja sentir aguas abajo a
distancias mucho mayores que aguas arriba. Lejos del cuerpo el mallado puede ser mucho más basto que en las
cercanías debido a que las variaciones de las variables son mucho menores.
En la generación de la geometría, el número de celdas que va a tener la simulación viene definido por
los nodos que se ponen en las líneas de contorno: muchos nodos implican muchas celdas pequeñas y al contrario.
Hasta que se tenga una cierta experiencia con los programas, es mejor poner pocos nodos. Compensa rehacer la
geometría para conseguir mejor definición o simular más adecuadamente el flujo que empezar con mallados que
tarden mucho en calcularse.
Si dos paredes opuestas están muy próximas, por ejemplo el espacio entre la parte inferior de un coche y
la carretera, aunque se esté haciendo un mallado basto, seguramente habrá que poner una mayor densidad de
nodos que en otras partes. Una distancia entre nodos del mismo orden de magnitud que la separación entre pare-
des, daría solamente una o dos filas de celdas entre ellas. Por otra parte, el espaciado de los nodos debe ser simi-
lar en las dos paredes, para evitar que las celdas queden demasiado deformadas.
Para simular adecuadamente la capa límite, sobre todo en los sitios donde existe desprendimiento, el
mallado debe ser de un orden de magnitud similar al espesor típico de ésta. Es conveniente no preocuparse mu-
cho por conseguirlo en las primeras pruebas, pues el número de celdas necesario es elevado.
Un objeto sólido en el interior de un fluido se define practicando un agujero con la forma del sólido en
la zona fluida. El objeto sólido no tiene que estar mallado. Se supone que, al ser sólido, el fluido no lo atraviesa.
No hacer falta, por tanto, calcular el flujo dentro de él. Para que el flujo rodee al objeto (inexistente) se utilizan
condiciones de contorno tipo pared en todas las líneas que dan la forma del objeto.
En cualquiera de las operaciones descritas, si uno se equivoca, se puede usar la tecla de deshacer .
Sin embargo esta tecla deja de estar activa cada vez que se guarda algún dato intermedio (Save). Esto último es
recomendable si se quieren evitar problemas desagradables.
19
UN EJEMPLO COMPLETO
Se simulará un coche introducido en un túnel de viento. Se utilizarán formas reducidas de los coman-
dos. Los textos detrás de la doble barra (//) son comentarios.
Para los nombres de los puntos se utilizará una sola letra: A – F contorno exterior; K – R coche.
// Abrir GAMBIT
File
New // crear una nueva sesión
ID: ruta de acceso al fichero coche // identificación de la sesión. Nombre para los ficheros
Title: coche // aquí se guardarán los ficheros cada vez que se selecciones file - save
Accept
Solver
FLUENT 5 // programa de resolución que se va a utilizar
// Definición de los puntos
En Operation selecciónese el icono Geometry (los nombres aparecen en la ventana Description)
En Geometry selecciónese el icono puntos (Vertex)
En Vertex Selecciónese el icono Create Vertex (From Coordinates dentro de las opciones desplega-
bles con el botón derecho)
// Introdúzcase el punto A (0, 0). Salvo que se haya definido otro sistema de geometría, las coordenadas
Globales y Locales coinciden. El tipo será Cartesianas (por defecto). La coordenada z se mantendrá 0 en todo el
caso, por ser bidimensional. Conviene recordar que, aunque se haya puesto nombre al punto, o posteriormente a
líneas o superficies, si éste forma parte de un conjunto que se va a unir a otro conjunto, el programa lo renombra
automáticamente. Si los puntos se van de escala úsese la opción autoescalar comentada anteriormente para
ver la figura completa.
// de la misma forma se introduce el resto de los puntos
Pto X Y Pto X Y
B 1.750 4.250 L 3.000 0.500
C 6.000 6.000 M 3.100 0.900
D 10.25 4.250 N 3.900 1.000
E 12.00 0.000 O 4.100 1.500
F 6.200 0.000 P 5.700 1.500
G 2.800 0.000 Q 6.000 0.600
K 3.100 0.300 R 5.800 0.300
// Selecciónese la tecla Close del menú Create Real Vertex.
// El programa pone nombres automáticamente si uno se ha olvidado de hacerlo.
20
// Definición de las líneas

En Geometry selecciónese el icono líneas (Edges)
En Edge selecciónese el icono Create edge. Dentro de las opciones de éste (botón derecho), hay que
escoger NURBS para las líneas ABC y CDE (también se podría hacer con ARC, siempre y cuando lo puntos
sean sólo tres), y Straight para el resto.
// Creación de la línea ABC.
Selecciónense los puntos A, B y C de la lista o con Shift + Botón izquierdo y pulse Apply.
// de la misma forma se hace con el resto de las líneas
Lin Ptos
CDE C,D,E
EF E,F
FG F,G
GA G,
KL K,L
LM L,M
MN M,N
NO N,O
OP O,P
PQ P,Q
QR Q,R
RK R,K
// Al acabar selecciónese Close para cerrar la ventana.
21
En el dibujo mostrado se ha empleado la opción ARC en lugar de NURBS.
Recuérdese que se puede hacer ZOOM (SHIFT + botón derecho del ratón y se selecciona un área de-
terminada). Si el dibujo se gira se puede volver a orientar con el botón de ejes.
Se puede activar o desactivar la escritura de los nombres de los puntos usando la tecla visualizar atribu-
tos , parte inferior derecha.
Si al hacer Zoom, por error se gira el dibujo se puede restablecer su posición usando la tecla cambiar
orientación ejes , parte inferior derecha.
22
// Generación de superficies
Se definirán dos superficies, una será el perímetro externo y la otra el coche. Posteriormente se realizará
la operación booleana de la substración.
// Creación de la superficie externa:

En Geometry seleccionar Face.
En Face seleccionar Create Face con la opción (botón derecho) Wireframe (generación de una super-
ficie mediante las líneas que forman su perímetro).
Se seleccionan, tal como se ha indicado, todas las líneas que configuran esa superficie Dibujo inferior
izquierdo).
// Poner nombre en Label.

// Creación de la superficie del coche:
Se repite el mismo proceso con las líneas del coche.
// Substracción (dibujo superior derecho)

Dentro de la sección Face, seleccionar el icono Boolean Operations y dentro de las opciones de esta
(botón derecho) la de Substract (es la que activa el submenú dibujado anteriormente, Substract Real Faces).
Face: TUNEL
Subtract:
Face: COCHE
En las dos Face se desactiva Retain. De esta forma creará una nueva superficie, la del tunel con el agu-
jero del coche, borrando las originales (esta opción está por defecto).
//Cerrar menú con Close
// Generación de los nodos
De Operation selecciónese el icono mallado (Mesh)

De Mesh selecciónese el icono líneas (Edges)
De Edges selecciónese el icono de mallado de líneas (Mesh Edges)
23
// Aunque se utilizará la opción Interval count (especificar el número de nodos), los números elegidos
suponen divisiones cada 0.1 m en el coche y cada 0.5 m en los bordes exteriores salvo debajo del coche que
serán cada 0.1 m. Como se puede ver es bastante basto.
// Mallado de la línea ABC

En el gráfico se observan los valores que hay que introducir. Aunque está activado el apartado Gra-
ding, con los valores elegidos no tiene ningún efecto. Se desactiva Apply y se generará el mallado automática-
mente.
// Se repite el proceso para el resto de las curvas

Lin Num Divisiones
ABC 20
CDE 20
EF 14
FG 36
GA 8
KL 4
LM 6
MN 10
NO 6
OP 18
PQ 10
QR 6
RK 28
// Cerrar con Close el menú
La líneas malladas se
ponen de color azul y se ven los
puntos en que se ha dividido la
línea.
Se puede observar que
las líneas que definen el coche y
las más próximas tiene mayor
densidad de puntos. Su finalidad
es que el mallado que se genere
se más denso en esa parte para
que pueda recoger más informa-
ción.
24
// Generación del mallado
De Mesh selecciónese el icono superficies (Face)

De Face selecciónese el icono de mallado de caras (Mesh Faces)
Se selecciona la superficie (TUNEL).

Se realizará mallado triangular: Elements Tri // Triángulos
Type Pave // No hay otra
Spacing: se desactiva Apply (se han mallado las líneas).
Options: se activan Mesh y Remove old mesh (no hay un mallado anterior, pero puede usarse esta op-
ción para corregir). Se desactiva Remove lower mesh (si no, se borrarían los nodos de las líneas; no pasaría
nada si se ha hecho todo lo demás correctamente).
Al seleccionar Apply en la parte

superior de la pantalla aparece una barra con
el % que se ha mallado y luego aparece el
dibujo de la malla, que será similar al mos-
trado en la figura.
25
// Especificación de los contornos (dibujo inferior izquierdo)

Se definen las condiciones de contorno: líneas que serán entrada, salida de flujo y paredes.
De Operation selecciónese Zones. De Zones selecciónese Specify Boundary Types.
En el apartado Entity hay que elegir Edges.
Después hay que definir: Línea ABC como Type: velocity-inlet
// Puede ocurrir que, por cuestiones de resolución de pantalla, no aparezca visible esta operación. Se selecciona
entonces otra opción cualquiera y al volver a cambiar la condición de contorno ya aparecen visibles todas las
opciones.
Línea CDE como pressure-outlet
Resto de líneas como wall. En estas se pueden definir varias líneas como la misma codición de contor-
no, pero al menos hay que definir dos: la carretera por una parte y el coche por otra, ya que habrá que asignarles
valores diferentes.
Tal como se ha comentado, y se puede ver en la figura, se han perdido los nombres de algunos
elementos por que se han hecho con ellos operaciones booleanas. Para saber cuál el nombre asignado por el
programa a cada uno de ellos se selecciona cada elemento y se activará su nombre.
Conviene dar nombres que sean fácilmente identificables.
En el programa FLUENT 5.1 se hablará acerca de estas condiciones.

// Especificación del medio (dibujo superior derecho)
De Zones selecciónese el icono Specify Continuum Types.
Como Entity hay que elegir Faces (por defecto está volumen).
Como sólo hay una superficie, sólo hay que especificar un medio, que en este caso es fluido. En la figu-
ra derecha se pueden ver los pasos que hay que seguir.
/ Para terminar nos despedimos
File
Save // Guarda los ficheros con el formato de Gambit.
File
Export
Mesh // Guarda el fichero con el formato que acepta el programa de resolución.
File
Exit // Esta es la forma correcta de salir del Gambit.
En el directorio de trabajo aparecen guardados todos los ficheros coche.*
26
FLUENT 5
De cómo utilizar el Fluent para calcular una posible solución de flujo. De

cómo leer un mallado, definir el fluido y características del flujo, condiciones de con-
torno, parámetros de la solución y de cómo verlo todo.
INTRODUCCIÓN
El programa llamado Fluent 5 es el que realiza el resto del proceso: definir los parámetros de ejecución,
realizar las iteraciones de cálculo y representar la solución obtenida. Parte de un mallado no estructurado, creado
en principio con el GAMBIT, aunque admite mallados generados por otros programas.
Se maneja íntegramente con el ratón a partir de los menús de la parte superior y, lo mismo que en el
GAMBIT, en la pantalla principal van apareciendo mensajes acerca de las instrucciones ejecutadas.
La ayuda en línea es muy completa, suple totalmente los manuales. Hay un botón ad hoc en cada una de
las ventanas que se abran y es muy recomendable consultarla siempre que se tenga alguna duda en un paso o
selección.
Las unidades que el programa utiliza son las del Sistema Internacional, metros para las coordenadas,
metros por segundo para las velocidades, pascales para las presiones (relativas), grados kelvin para las tempera-
turas. Se pueden cambiar para algunas cosas, aunque no es aconsejable.
El proceso de operaciones que hay que realizar es el siguiente:
- Leer mallado
- Definición de modelos y fluido
- Especificar las condiciones de contorno
- Fijar los parámetros de resolución para el proceso iterativo
- Ver los resultados
- Análisis de la solución
Al abrir el programa se elige entre la ejecución en 2D o 3D (las opciones 2ddp y 3ddp son para realizar
los cálculos con doble precisión), después se abre la ventana principal:
27
LEER EL MALLADO
Para leer el mallado se usa el menú File:
File
Read
Case…
Se lee como Case, aunque un caso propiamente dicho incluye el mallado y las demás condiciones y pa-
rámetros que se imponen para la ejecución. Los ficheros Data son los que guardan la información de la solución.
Estos son los formatos para guardar y volver a leer los trabajos realizados.
Al leer el mallado Coche.msh aparece lo siguiente:
> Reading "E:\Gambit\coche.msh"...

1790 nodes.
20 2D velocity-inlet faces, zone 2.
20 2D pressure-outlet faces, zone 3.
8 2D wall faces, zone 4.
4998 2D interior faces, zone 16.
3394 triangular cells, zone 1.
Done.
Building...
grid,
materials,
interface,
zones,
default-interior
wall.14
wall.13
wall.12
wall.11
wall.10
28
wall.9
wall.8
wall.7
wall.5
wall.4
wall.3
salida
entrada
fluido
Done.
Además de leer el mallado adjudica parámetros por defecto a las condiciones de contorno y al fluido.
En el menú Grid se puede ver la información relativa al mallado leído
Grid
Info
Size
Grid Size
Level Cells Faces Nodes Partitions
0 3394 5184 1790 1
1 cell zone, 14 face zones.
También en Grid se puede cambiar la escala del dibujo con Scale. Hay que tener en cuenta que el
Gambit trabaja con unidades del Sistema Internacional.
La opción de Smooth/Swap… puede servir para intentar mejorar las características de Skew de un ma-
llado no muy bueno generado con un programa diferente de Gambit (si se ha hecho con Gambit es mejor esfor-
zarse por hacerlo bien desde el principio).
MODELOS
Después de leer el mallado, lo primero que hay que hacer es definir los modelos que se van a utilizar.
Define
Models
dentro de los modelos están las siguientes opciones:
29
Solver
El sistema de solución (Solver) se deja el que tiene por defecto
(Segregated). El sistema Coupled es más apropiado para el caso de flujos
compresibles, sobre todo si aparecen velocidades supersónicas. La opción
Formulation se deja Implicit. La opción Velocity Formulation se deja tam-
bién la que está por defecto (Absolute), por el momento no se van a utilizar
zonas con movimiento relativo.
Permite elegir tres tipos de espacio (en 2D): bidimensional sim-
ple, axisimétrico y axisimétrico con una componente de giro. Los dos
últimos casos son bidimensionales en el fondo, aunque las áreas y algunos
otros pequeños detalles se calculan de manera diferente. Para el coche se
dejará en 2D. También permite definir si el caso es estacionario o no. En el
caso estacionario se hace un solo cálculo iterativo, en el no estacionario, se
parte de unas condiciones iniciales para todas las variables y se va calcu-
lando su evolución a lo largo del tiempo; el programa realiza una iteración para cada instante de tiempo fijado.
En el caso del coche, el proceso de aceleración sería no estacionario aunque, por ahora, se va a resolver el caso
de velocidad constante, estacionario: steady.
Viscous
Se puede seleccionar entre resolver el caso de forma laminar o turbulenta y para este último, entre va-
rios modelos. Dentro de cada modelo se pueden ajustar varias constantes y otros parámetros. No se va a profun-
dizar aquí en este tema; baste saber que cada modelo puede ser más apropiado en unos casos que en otros: si el
flujo es externo o interno, si hay estelas o capas de cortadura, desprendimiento de la capa límite, etc. Con las
constantes cambian las características del modelo y en el caso de tener datos experimentales o soluciones de
casos similares se puede ajustar un poquito. Por el momento se elegirá el modelo k-epsilon (2 eqn), Standard
con las constantes por defecto.
Las opciones de Near Wall Treatment se refieren a diferentes modelos para simular los efectos de la
capa límite. Téngase en cuenta que a menos que el mallado sea lo suficientemente fino en los contornos de las
superficies sólidas (varias capas de celdas en el espesor de la capa límite) , no se calcularían bien de forma dire-
cta las velocidades y pérdidas de la capa límite. Estos modelos permiten no tener que definir el mallado hasta
esos extremos, aunque si se quiere tener un buen resultado, éste debe ser de un orden de magnitud similar al
espesor de la propia capa límite.
30
Energy
Hay que activar a opción Energy cuando haya transferencia de ca-
lor dentro del fluido y/o del fluido con las paredes. Por supuesto, si el flujo
es compresible tiene que estar activada. Un flujo de un gas en el que las
velocidades esperadas en algún punto sean mayores de un tercio de la del
sonido hay que tratarlo como compresible; si son inferiores, se puede con-
siderar incompresible.
En realidad, más que utilizar un modelo, lo que hace el programa es añadir la ecuación de la energía a
las de continuidad y cantidad de movimiento que ya resuelve siempre.
Para el coche no hace falta activar la resolución de esta ecuación.
El resto se refiere a modelos para flujos especiales, en los que se quiere tener en cuenta la radiación,
existen varias especies, para una fase sólida transportada por el fluido o para la dispersión de contaminantes.
MATERIALES
Ahora hay que definir el/los fluidos/sólidos que se vayan a utilizar:
Define
Materials
Por defecto el programa sólo tiene seleccionado un único sólido (aluminio) y un único fluido (aire). Si
se fuese a trabajar con otro fluido (agua, por ejemplo) se tendría que seleccionar de una base de datos que trae el
programa y copiarlo a la lista de fluidos. Una vez realizado se vería que los fluidos disponibles serían aire y
agua.
Si se quiere añadir un material se pulsa Database, se selecciona el material deseado y se pulsa Copy, se
selecciona otro, Copy… y Close.
Los materiales sólidos sólo son de aplicación cuando hay transferencia de calor con las paredes y se
quiere que el programa calcule esta transferencia en función de las propiedades de la pared.
Respecto a los fluidos, si no se ha elegido el modelo de transferencia de calor, únicamente aparecen las
propiedades de densidad y viscosidad. La densidad es constante y para la viscosidad se puede elegir entre varias
leyes. Si se especifica un valor o se cambia algo en propiedades hay que pulsar Change/Create para que surta
efecto.
31
Para casos con gases en los que haya que resolver la transferencia de calor, lo habitual será utilizar, no
una densidad constante, si no la obtenida a partir de la ecuación de los gases perfectos. Esto se selecciona dentro
de las opciones (flecha hacia abajo) de Density. Al activar la ecuación de la energía, aparecen también, dentro de
las propiedades, la capacidad calorífica y la conductividad térmica. También aparece la posibilidad de utilizar
leyes por tramos o polinomios del cambio de las propiedades en función de la temperatura, en vez de mantener-
las constantes.
CONDICIONES DE OPERACION
En el apartado de Operating conditions aparece la opción
de definir la presión de referencia (operating pressure) este valor es
respecto al cual se calculan los valores de presión en los resultados;
normalmente se pondrá aquí la presión atmosférica, de forma que se
obtengan valores de presión relativa.
Si es importante la presencia de la gravedad en un fluido se
activa la opción Gravity. Esto sucede, por ejemplo, en casos bidi-
mensionales cuando el plano de trabajo es vertical y hay cambios de
densidad por transferencia de calor (o por coexistir varias especies o
fases...).
En el caso que se está trabajando no hay transferencia de calor, el fluido es aire y el flujo incompresible
(hasta unos 300 km/hora), por lo que se pueden dejar las opciones mínimas que vienen por defecto en este apar-
tado y el anterior.
CONDICIONES DE CONTORNO
Las condiciones de contorno, junto con la geometría, son las que definen la solución de las ecuaciones
diferenciales que representan el flujo, por lo que son uno de los factores más críticos para que el resultado sea
realista.
En Fluent se definen las condiciones de contorno en zonas. Existen tres grupos de zonas y cada uno
puede tener unos tipos de condiciones de contorno diferentes. Si no se asigna un nombre manualmente, las zonas
se nombran con el tipo de condición de contorno, que tienen asignado en cada momento, y un número. El tipo de
condición se puede cambiar, pero sólo dentro de las que son válidas para su grupo.
1. Contornos exteriores: entradas y salidas, paredes, ejes y planos de simetría (pressure inlet, velocity
inlet, pressure outlet, outflow, wall, axis, symmetry).
2. Zonas de celdas internas: Internal cell zones (fluid & solid).
3. Caras internas: Internal face boundaries (wall & interior)
La relación anterior no es exhaustiva, pero las zonas citadas son las más importantes. En todo este apar-
tado es muy conveniente consultar la ayuda en línea. También es la forma de ir aprendiendo las características de
otros tipos de condiciones de contorno.
En las zonas de celdas internas se especifica el fluido (o sólido) que está ocupando ese grupo de celdas.
En un caso sencillo sólo existirá una zona de celdas interior. Se pueden obtener varias separando la geometría en
varias zonas cerradas (por ejemplo, en la geometría generada para el coche existían originariamente dos superfi-
cies: la total y el coche; al hacer la substracción se dejó solamente la resultante y las líneas que definen el coche
son ahora contornos exteriores, pero se podría haber mantenido la zona del coche...).
Aparte de especificar el tipo de fluido, también se pueden definir en ellas fuentes de calor, de masa, de
cantidad de movimiento… El moving reference frame, que permite que unas zonas estén en movimiento relati-
vo con respecto a otras, es otro concepto distinto.
Cada zona de celdas interna tiene una zona de caras internas asociada. Estas tienen que ir definidas co-
mo interior y son de uso interno del programa, se pueden olvidar. Existe también otro tipo de caras internas que
surge cuando se ha definido una línea dentro de la zona de celdas (el contorno del coche si no se hubiera borrado
la zona interior). Estas caras se utilizan habitualmente para separar zonas, pero sin bloquear el paso del fluido -
se les asigna entonces condición de contorno interior.
32
Los contornos exteriores son el meollo de la cuestión. Aunque se les puede asignar cualquiera de las
condiciones válidas para este grupo, en beneficio de la explicación se van a separar en entradas-salidas y pare-
des.
Las paredes imponen en primer lugar la condición de que no las puede atravesar el flujo, con lo que la
componente normal de velocidad es nula. En segundo lugar la condición de no deslizamiento, que impone que la
componente tangencial de velocidad sea también nula. Esta segunda condición viene originada por el rozamiento
del fluido con una superficie sólida y es la que genera la capa límite, que se trata con el modelo correspondiente.
Las condiciones de planos y ejes de simetría permiten reducir la geometría de cálculo cuando la geometría tiene
una simetría especular o axial. Afectan al flujo de manera que se comporte como si hubiera otra zona exactamen-
te simétrica o axisimétrica. En el fondo, estas condiciones se reducen a paredes en las que no existe rozamiento y
no se aplica, por tanto, la segunda condición.
Las entradas y salidas son las que van a provocar normalmente el flujo y habitualmente son imprescin-
dibles, al menos una de cada (un ejemplo de lo contrario puede ser el estudio de la circulación del aire en una
habitación cerrada con una fuente de calor: el flujo se produce por convección natural). Básicamente hay dos
tipos de entrada: velocity inlet y pressure inlet y dos tipos de salida: pressure outlet y outflow. En la velocity
inlet se impone la velocidad de entrada. En la pressure inlet se impone la presión total y la dirección de la velo-
cidad, la magnitud de la velocidad de entrada queda entonces a expensas de la presión estática que vaya calcu-
lando el programa. En la pressure outlet lo que se fija es la presión estática de salida. En la condición de out-
flow no se fija ninguna variable, el programa fuerza que las derivadas de las variables normales a este contorno
sean nulas. Se van a estudiar ahora con más detalle algunas de estas condiciones.
Velocity inlet
La condición de contorno
velocity inlet se utiliza para imponer
una velocidad de entrada constante.
También hay que fijar otras variables
escalares: temperatura, parámetros de
turbulencia… pero no la presión.
Esta se ajustará al valor necesario
para asegurar que se da la velocidad
exigida. Suele utilizarse cuando se
conoce el caudal de entrada, aunque
también es muy adecuado para flujos
externos cuando se modelizan en
movimiento relativo, por ejemplo: un
vehículo moviéndose a una cierta
velocidad, para que el caso sea esta-
cionario se simula con el vehículo
fijo y el fluido avanzando hacia él.
Esta condición no es apro-
piada para flujos compresibles por-
que puede llevar a valores absurdos de las condiciones de estancamiento.
Excepcionalmente se puede utilizar para definir la velocidad en una salida, porque es el dato que se co-
noce o se quiere imponer. Si se hace esto, hay que tener especial cuidado en que se cumpla la continuidad global
En este caso no es de aplicación la definición de las variables escalares (se toman los valores de la parte interior).
En el panel de especificación de esta condición aparecen las siguientes opciones:
- Zone name: se puede cambiar el nombre.

- Velocity Specification Method: permite elegir entre tres formas de definir la velocidad de entrada:
- Magnitud y dirección
- Componentes según los ejes cordenados
- Perpendicular a la superficie de entrada
- Reference Frame: es indiferente mientras no se utilicen zonas con movimiento relativo.
33
- Temperature: en caso de tener activada la transferencia de calor, junto con la velocidad hay que fi-
jar también la temperatura.
- Turbulence Specification Method: para flujo no laminar. Permite definir la turbulencia de entrada
de varias formas. Habitualmente son necesarios datos experimentales de caso similares (en caso de
duda, déjense los valores por defecto).
Pressure inlet
Se utiliza para fijar la presión total de
entrada, junto con las otras variables escalares
(temperatura, variables de turbulencia). Es muy
útil cuando se conoce la presión pero no el caudal
o la velocidad, lo que suele ser frecuente en flujos
internos (por ejemplo, una válvula entre un depó-
sito a presión y la atmósfera; se conoce la presión
total de entrada: la del depósito, y la presión está-
tica de salida: atmosférica, pero el caudal va a
depender de las pérdidas que provoque). Es apro-
piada tanto para flujo compresible como incom-
presible. Aunque no hay que definir la magnitud
de la velocidad, si hay que especificar la direc-
ción.
- Gauge Total Pressure: presión total (relativa) de entrada.

- Supersonic / Initial Gauge Pressure: presión estática; sólo se usa para iniciar variables y en el ca-
so de que la entrada llegue a ser supersónica durante las iteraciones.
Outflow
Este tipo de condición de contorno se utiliza en la sali-
da cuando la distribución de velocidad y la presión no se cono-
cen (el caudal sí se debe conocer, está fijado por la entrada). No
hace falta definir nada. Representa una salida donde el flujo está
completamente desarrollado, cuando no hay cambio de sección
y los perfiles de velocidad han llegado a ser constantes. Esto se
basa en que la condición que en realidad se está imponiendo es
que sean nulos los gradientes normales de todas las variables, excepto la presión; los gradientes tangenciales no
tienen por qué ser nulos. Los valores de las variables a los que se llega son extrapolación de los valores en el
interior.
No es apropiada cundo el flujo es supersónico ni cuando se utilizan pressure inlets. En la práctica, un
caso con esta condición no llega a converger si a través de esta salida existe recirculación, y funcionará mal si
está cerca de una estela o similar.
La opción de Flow Rate fraction sirve para repartir el caudal entre varias salidas con outflow, aunque
lo normal es que la salida sea única.
Pressure outlet
Fija la presión estática en la salida. Es
obligatorio utilizar esta condición cuando se ha
fijado una pressure inlet. No se puede mezclar con
la condición de outflow. Se especifican también
una serie de valores (temperatura y variables de
turbulencia) por si existe flujo inverso (backflow)
en la solución final o durante las iteraciones. En
realidad no es nada aconsejable, para hacer una
buena simulación, que a través de una condición de
34
contorno haya flujo inverso, habría que extender más allá el dominio de cálculo.
Condiciones correctas de entrada – salida

La correcta selección de condiciones de contorno de entrada y salida depende del régimen de flujo:
- Flujo incompresible (obviando la temperatura) y compresible subsónico.
- Entrada: presión y temperatura total (pressure inlet). Salida: presión estática (pressure ou-
tlet).
- Entrada: velocidad y temperatura total (velocity inlet). Salida: presión estática (pressure ou-
tlet).
- Entrada: velocidad y temperatura total (velocity inlet). Salida: outflow.
- Se pueden mezclar entradas de velocidad y presión total (velocity inlet & presure inlet). Sali-
da: presión estática (pressure outlet).
- Se pueden usar las condiciones de velocidad y de presión total como salidas.
- Flujo compresible transónico: sónico en el interior pero con entradas y salidas subsónicas.
- Entrada: presión y temperatura total (pressure inlet). Salida: presión estática (pressure ou-
tlet).
- Flujo supersónico en entradas o salidas: se pueden utilizar como entradas pressure inlet y mass
flow inlet; como salidas pressure outlet y pressure far field. Para estos casos es necesario utilizar
el solver de tipo coupled. En estos apuntes no se van a tratar estos casos.
En flujo incompresible, subsónico o transónico se pueden usar sin problemas varias entradas y/o salidas
del mismo tipo, con las restricciones anteriores.
Recuérdese que siempre que se especifica la presión se habla de presión relativa respecto a la presión de
operación (operating pressure) definida en Operating Conditions.
Wall
Las paredes definen los lími-
tes físicos de la geometría. En cual-
quier caso, se impone la condición de
que el flujo no atraviese la pared.
Normalmente también se impone la
condición de velocidad tangencial
nula, aunque se pueden simular pare-
des que se mueven con respecto a la
zona de cálculo definiendo la veloci-
dad (Moving wall) o bien el esfuerzo
cortante (Specified shear).
En el caso del coche, se puede hacer
que la pared que simula la carretera se
mueva hacia atrás a la velocidad relati-
va del coche, para simular correctamen-
te el paso del aire por debajo del vehí-
culo.
La opción de Wall Rough-
ness permite corregir las leyes de la
pared para simular correctamente la
capa límite sobre superficies rugosas.
Como en muchos otros casos, para aplicar estos valores correctamente suele ser necesario contar con datos expe-
rimentales. De forma simplificada, se puede mantener la constante de rugosidad (Roughness Constant) en 0.5,
que va bien para muchos casos, y poner en la altura de rugosidad (Roughness Height) un valor estimado.
Cuando se activa la transferencia de calor, se ponen en las paredes las condiciones térmicas de contor-
no. Hay varias opciones: fijar el flujo de calor, fijar la temperatura en la pared, o utilizar algún modelo (convec-
35
ción, radiación) que calcule el flujo de calor a partir de una cierta temperatura exterior. Si se utiliza Heat Flux,
un valor de cero es simplemente la condición adiabática. Valores positivos introducen calor. La opción de Tem-
perature fija el valor de la temperatura en la pared; se introduce o extrae el calor necesario para que permanezca
constante.
Las opciones de Wall Thickness y Material se utilizan únicamente cuando la pared separa dos zonas
de cálculo, simulando la transferencia de calor por medio de una ley unidimensional (ver la ayuda para más
información).
Symmetry & Axis

Se utilizan en casos que tienen una simetría especular o
axial y permiten definir solamente la mitad de la geometría.
Esta condición actúa en el flujo de manera que se comporte
exactamente igual que si se tuviera también la otra zona. En
realidad funciona mejor, porque, si se define todo el dominio,
no se va a conseguir una simetría perfecta, empezando porque el
mallado no la tendrá. No hay que definir ningún parámetro en
estos tipos de condición.
Entonces, ¿qué CONDICIONES hay que PONER con el coche?

El caso del coche requiere un planteamiento de movimiento relativo. Un objeto se está moviendo con
respecto al aire y al suelo (que están en reposo). En principio, la situación es no estacionaria ya que un punto del
fluido cambia del reposo al movimiento durante el paso del coche, para volver al reposo después. Este plantea-
miento desde un sistema de referencia absoluto es bastante complejo. Si se consideran, en cambio, unos ejes de
referencia que se muevan con la misma velocidad que el coche, el objeto está en reposo, el aire y el suelo se
acercan a él a la misma velocidad a la que el objeto se movía en el
sistema de referencia absoluto, pero con sentido contrario. En esta
nueva situación, las velocidades del fluido en todos los puntos perma-
necerán constantes en el tiempo, con lo que se ha conseguido una
situación estacionaria. Al ser constante la velocidad del sistema de
referencia, no aparecen fuerzas inerciales, por lo que la distribución de
presiones y las fuerzas que provocan sobre las superficies sólidas son
las mismas. El campo de velocidades absoluto se puede obtener, si se
quiere, restando a cada punto la velocidad del sistema de referencia.
Este planteamiento es mucho más sencillo para el cálculo.
En cuanto a las condiciones de contorno, el coche está parado,
el aire tiene que entrar a la velocidad a la que avanza el coche en el
sistema absoluto, pero con sentido contrario y la carretera avanza hacia
el coche a la misma velocidad. Las condiciones de salida podrían ser de
outflow, pero, tal como se ha definido la geometría, en la parte superior
por encima del coche el flujo tiende a salir del contorno y detrás del
coche a volver a entrar, lo que invalida esta condición. Es mejor
poner una condición de presión estática constante (atmosférica) en
la parte posterior del dominio de cálculo.
Al abrir la ventana de las condiciones de contorno, si se
han seguido los pasos indicados hasta el momento, tendrá una
forma similar a la de la figura. En caso de duda sobre a qué zona
corresponden se puede comprobar con Display Grid. En esta
ventana de Boundary Conditions hay una lista de zonas y otra de
tipos. Cuando se selecciona una determinada zona, en la lista de
tipos aparecen los tipos de condiciones de contorno compatibles
con esa zona. Se muestra sombreado el tipo actualmente elegido
para esa zona.
Pulsando dos veces en el tipo sombreado o pulsando Set,
se abre la ventana característica de esa condición. Si se desea
cambiar el tipo basta con pinchar en el tipo deseado; pide confir-
mación y abre la ventana característica.
36
Si se selecciona la zona tunel y se pulsa set, aparece el panel donde se puede cambiar el nombre, etc. El
fluido por defecto es aire, así que no hace falta cambiarlo –tampoco se podría a menos que se hubiera seleccio-
nado otro fluido en la base de datos de materiales. Para el caso del coche no hay que tocar nada en esta condi-
ción.
La condición de default-interior, solo permite cambiar el nombre:
En la zona de entrada del flujo: entrada lo más práctico es elegir el Velocity Specification Method
como Magnitude and Direction. Un valor positivo de la mag-
nitud (40 m/s : 144 km/hora) define la dirección según el eje
positivo de las x, de frente hacia el coche, tal como se busca.
Los valores de los parámetros de turbulencia definen la turbu-
lencia del aire lejos del coche; se dejan por defecto a falta de
otros datos.
La zona de salida, que se había definido como salida
en GAMBIT, no hay que cambiarla. Se pulsa set y aparece su panel.
Por defecto tiene un valor de pre-
sión estática relativa cero (atmosférica). Los
valores de los parámetros de turbulencia
para el caso de Backflow se mantienen igua-
les a los de entrada.
Las paredes son wall.3, wall.4 y
wall.5 que corresponden a la carretera y
wall.6 ...wal13, que corresponden al coche.
En la superficie correspondiente al coche no
hay que definir nada, porque según el
planteamiento
teamiento relativo
relativo
que se
que
ha se
hecho,
ha hecho,
el coche
el está parado y
su rugosidad para este caso se puede despreciar. En la
carretera,
Se puede poner una altura de rugosidad de medio
centímetro (por ejemplo) y marcando Moving Wall, una
velocidad de traslación en la dirección x de 40 m/s (la
misma que la del aire). Es indiferente que sea Relative to
adjacent Cell Zone o Absolute, porque el Motion Type
de la zona fluida es Stationary.
Lo mismo para wall4 y wall.5.
Con esto estaría definido el caso. Hay que ajustar
ahora los parámetros propios de la resolución.
PARÁMETROS DE RESOLUCIÓN
El procedimiento general sería el siguiente:
1. Escoger el esquema de discretización y el método de acoplamiento entre las ecuaciones de presión
y velocidad.
2. Fijar los valores de relajación.
3. Seleccionar y fijar los parámetros del esquema de resolución Multigrid.
4. Valores iniciales.
5. Fijar los valores de convergencia y la monitorización de los residuos.
6. Definir un número de iteraciones y comenzar el cálculo.
Los tres primeros pasos se refieren a opciones del solver. Los valores se pueden cambiar en el menú
Solve
37
Controls
Los valores por defecto son apropiados para la mayoría de los problemas. Por ahora es conveniente no
tocarlos salvo cuando se tienen problemas de convergencia, estabilidad o precisión. En general no es que sean
los valores con los que se obtienen mejores resultados, sino los que más fácilmente obtienen una solución cohe-
rente.
INICIALIZACIÓN
La inicialización de la solución consiste en
asignar valores iniciales a todas las variables en
todos los puntos de cálculo. Estos serán los valores
de los que partirán las iteraciones. Habitualmente
los valores por defecto son buenos, aunque cuanto
más aproximados sean a la solución, menos tardará
en converger. En ocasiones la asignación de valores
adecuados puede ser crítica para evitar la divergen-
cia en las primeras iteraciones. La inicialización se
realiza en el menú
Solve
Initialize
Initialize
Los valores iniciales para cada variable están en la ventana Initial Values. Se pueden fijar valores ma-
nualmente o bien elegir una zona en Compute From para que los calcule automáticamente. En este caso se han
tomado como valores iniciales en todas las celdas el valor de la velocidad de entrada.
Pulsando Init se realiza la inicialización: se asignan los valores a la memoria de cálculo. Pulsando Ap-
ply se pueden guardar los valores pero sin asignarlos a la memoria (también se guardan al pulsar Init) y con
Reset se recuerdan.
La opción Reference Frame solo cambia algo cuando se tienen zonas con distinto movimiento relativo.
Si ya se ha realizado una serie de iteraciones pero no se ha llegado a la convergencia de la solución, se
puede seguir iterando a partir de lo conseguido pero ¡sin inicializar de nuevo!. Si no se ha cambiado nada de los
parámetros de cálculo, inicializar llevaría a repetir inútilmente las iteraciones anteriores. Si se ha cambiado algo,
la última solución puede ser un buen punto de partida, que reduzca mucho el número de iteraciones necesarias, o
quizá no...
CRITERIO DE CONVERGENCIA Y MONITORIZACIÓN

Después de realizar cada iteración el programa calcu-
la los residuos. Estos son una medida de los desequilibrios que
permanecen en las ecuaciones de conservación ( ver Defini-
tion of Residuals en la ayuda, para más información). Si los
cálculos se realizaran con precisión infinita, los residuos ten-
derían a cero cuando la solución fuese convergiendo. En la
práctica se alcanzan antes los errores de redondeo y el número
de iteraciones necesario para reducirlos aumenta desaforada-
mente cuanto más pequeños son. A menudo no es necesario
llegar tan lejos para tener una solución de una precisión cohe-
rente con la simplificaciones hechas en el caso.
El criterio de convergencia por defecto del Fluent es
que los residuos (Scale Residuals) bajen de 0.001, salvo la
ecuación de la energía que es 0.000001. Los Scale Residuals
son los residuos normales divididos por un factor representati-
vo del flujo en el dominio. Este criterio es adecuado para la
38
mayoría de los problemas (más información en Judging Convergence). Si la prueba de convergencia está activa-
da, el cálculo se termina cuando se alcanza el criterio en todas las ecuaciones, aunque no se haya llegado al nú-
mero de iteraciones señalado. Si se detiene la ejecución por llegar al número de iteraciones fijado, pero no se ha
conseguido la convergencia, hay que continuar iterando.
Una forma práctica de saber si la solución marcha adecuadamente es representar gráficamente los valo-
res de los residuos y algunas otras variables en función del número de iteraciones. Cuando la solución tiende a
converger, los residuos van disminuyendo progresivamente y los valores de las variables se hacen constantes.
El valor de los residuos se va guardando en el fichero de datos. Se pueden monitorizar en
Solve
Monitors
Residual
La opción de Plot es la que activa la representación gráfica (es muy entretenido mientras se espera a que
converja, incluso da una buena idea de si lo está haciendo bien). Las ventanas inferiores indican las ecuaciones
en las que se realiza la monitorización y el valor del criterio de convergencia.
Aunque no tienen tanto interés como los residuos, se pueden monitorizar otras variables mientras se rea-
liza la resolución. No tiene tanto interés, porque los valores válidos son los de la solución definitiva, aunque, en
ocasiones, también son un indicativo de que la resolución
está siendo satisfactoria. Por ejemplo, para ver la evolu-
ción de la fuerza de arrastre sobre el coche mientras se
realizan las iteraciones:
Solve
Monitors
Force
Representa un coeficiente de fuerza (arrastre, sustentación o mo-

mento). Para que los valores tengan sentido hay que ajustar los parámetros
del coeficiente en:
Report
Reference Values
Para el arrastre hay que indicar la velocidad de referencia, la densi-

dad (que ya está) y el área frontal, que, como la profundidad está en 1 m,
viene definida por la altura del coche (la profundidad sólo interviene en el
cálculo de áreas; el cálculo es bidimensional, equivalente a profundidad
infinita).
La evolución de otras variables se puede comprobar con
Solve
Monitors
Surface
Antes de iterar conviene guardar el fichero con todas las condicio-
nes de contorno puestas. Se haría con save case y el nombre que tiene por
defecto: coche.cas.
39
A ITERAR
Se fija el número de iteraciones en
solve
Iterate
En Number of Iterations se indica el número de iteracio-

nes. ¿Cuántas hay que poner? Los casos más sencillos llegarán a la
convergencia en unos cuantos cientos de iteraciones, los más compli-
cados necesitarán varios miles... La única indicación sencilla que se
puede dar es que si, después de muchas iteraciones, todavía está lejos
de converger, se tienen problemas.
El número de Reporting Interval indica cada cuántas itera-
ciones se monitorizan los residuos y demás variables. Los casos que necesitan muchas iteraciones se agilizan un
poco monitorizando menos.
UDF Profile Update Interval señala cada cuántas iteraciones se actualizan las funciones definidas por el
usuario. Por el momento no tiene utilidad.
Púlsese Iterate.
Se abren las ventanas...
...
..............................................................................
! 135 solution is converged
Evolución de los residuos: Evolución del coeficiente de arrastre:
Residuals
continuity
x-velocity 1e+01
y-velocity 20.0000
k
epsilon 18.0000
1e+00
16.0000
1e-01
14.0000
1e-02 12.0000
Cd 10.0000
1e-03
8.0000
1e-04 6.0000
4.0000
1e-05
2.0000
1e-06 0.0000
0 20 40 60 80 100 120 140 0 20 40 60 80 100 120 140
Iterations Iterations
Después de un pico inicial, los residuos disminuyen progresivamente y el arrastre permanece constante
en un valor algo inferior a 2, que no es una aberración, al menos a primera vista.
Antes de ver los resultados es muy recomendable guardarlos . Se hace con save date y sirve el nombre
que tiene por defecto: coche.da.t
VER LOS RESULTADOS

Tres son los menús que permiten analizar los resultados: Display, Plot y Report. En el primero se pue-
den ver, entre otros, los mapas de velocidades de forma vectorial y mapas de contornos (presión, temperatura,
líneas de corriente, etc.). Plot se utiliza principalmente para construir gráficos XY que representen, por ejemplo,
40
la distribución de presión sobre una superficie. En Report se consiguen las fuerzas e integraciones variadas de
los diferentes parámetros.
Display
Velocity Vectors
7.11e+01
6.40e+01
5.70e+01
4.99e+01
4.29e+01
3.58e+01
2.88e+01
2.17e+01
1.47e+01
7.61e+00
5.56e-01
Display Vectores de velocidad

Contours
depende de la variable y opción elegida:
1.92e+03
1.41e+03
9.12e+02
4.08e+02
-9.47e+01
-5.98e+02
-1.10e+03
-1.60e+03
-2.11e+03
-2.61e+03
-3.11e+03
Presión estática
1.68e+02
2.74e+02
1.52e+02
2.43e+02
1.35e+02
2.13e+02
1.18e+02
1.82e+02 1.01e+02
1.52e+02 8.46e+01
1.22e+02 6.78e+01
9.12e+01 5.11e+01
6.08e+01 3.44e+01
3.04e+01 1.76e+01
0.00e+00 9.00e-01
(Velocity – Stream Function, sin Filled Turbulence - Turbulent Kinetic Energy
41
3.35e+02
2.07e+03
2.71e+02
1.69e+03
2.06e+02
1.32e+03
1.41e+02
9.37e+02
7.70e+01
5.59e+02
1.25e+01
1.81e+02
-5.20e+01
-1.97e+02
-1.17e+02
-5.75e+02
-1.81e+02
-9.53e+02
-2.46e+02
-1.33e+03
-3.10e+02
-1.71e+03
(Pressure – Total Pressure: indicación pérdidas energía Vorticidad (Velocity – Vorticity Magni-
tude),
Obsérvese el sentido de giro contrario por encima y debajo del coche y en parte superior e inferior de la
estela.
Los gráficos XY se suelen utilizar par ver con mayor exactitud la evolución de una variable a lo largo
de una superficie o contorno:
Plot
XY Plot
wall-3
2.00e+03
1.00e+03
0.00e+00
Static -1.00e+03
Pressure
(pascal)
-2.00e+03
-3.00e+03
-4.00e+03
3 3.5 4 4.5 5 5.5 6
Position (m)
Distribución de presión sobre la superficie del coche:
Para hallar el valor numérico de las fuerzas sobre la superficie del coche se utiliza el menú Report:
Report
Forces
42
Fuerzas en la dirección horizontal y coeficiente de arrastre:

Force vector: (1 0 0)
pressure viscous total pressure viscous total
zone name force force force coefficient coefficient coefficient
------------------------- -------------- ------------------- ----------------- ---------------- ---------------------- ---------------
wall.6 250.83672 0.40921041 251.24593 0.21329653 0.00034796804 0.21364449
wall.7 677.29822 0.19733027 677.49555 0.57593386 0.00016779785 0.57610166
wall.8 54.401051 1.305827 55.706878 0.046259397 0.0011103971 0.04736979
wall.9 318.45258 0.6462431 319.09882 0.27079301 0.00054952644 0.27134253
wall.11 622.005 -0.45319721 621.55181 0.52891582 -0.00038537178 0.52853045
wall.12 213.8665 -0.018207315 213.84829 0.18185927 -1.5482411e-05 0.18184379
wall.13 0 13.590525 13.590525 0 0.011556569 0.011556569
------------------------- --------------- ------------------ ---------------- ------------------ -------------------- ------------------
net 2136.8601 15.677731 2152.5378 1.8170579 0.013331404 1.8303893
Fuerzas en la dirección vertical y coeficiente de sustentación:

Force vector: (0 1 0)
pressure viscous total pressure viscous total
zone name force force force coefficient coefficient coefficient
------------------------- --------------- ------------------ ---------------- ------------------ -------------------- ------------------
wall.6 125.41835 -0.81842226 124.59993 0.10664826 -0.00069593729 0.10595232
wall.7 -169.32439 0.78932172 -168.53506 -0.14398332 0.00067119194 -0.14331213
wall.8 -435.20877 0.16322851 -435.04554 -0.37007549 0.00013879975 -0.36993669
wall.9 -127.38091 1.6156093 -125.7653 -0.1083171 0.0013738174 -0.10694329
wall.11 207.33536 1.359594 208.69495 0.17630558 0.0011561173 0.17746169
wall.12 -142.57756 -0.027311049 -142.60487 -0.12123942 -2.3223681e-05 -0.12126265
wall.13 -1874.7408 0 -1874.7408 -1.5941674 0 -1.5941674
------------------------- --------------- ------------------ ---------------- ------------------ ------------------- ------------------
net -2416.4788 3.0820202 -2413.3967 -2.0548289 0.0026207655 -2.0522081
Aunque pueda no ser fácil, es muy importante indicar correctamente las variables que se representan y
sus unidades. En esta exposición no se hace hincapié en ello, asumiendo que las unidades son las del SI que
utiliza el programa, pero en un informe habría que indicarlo explícitamente en todas las gráficas y figuras.
Como se ve, una vez resueltas las ecuaciones, es muy fácil sacar una gran cantidad de información. La
dificultad está en saber interpretarla y, sobre todo, en asegurarse de que la solución es correcta.
ANÁLISIS DE LA SOLUCIÓN
Si se tiene información experimental o
analítica suficientemente completa del caso que
se está resolviendo, el análisis de la solución se
limita a comparar (y habitualmente ... intentar
justificar las diferencias). Pero esto no suele ser
lo normal. Cuando se estudia un caso del que se
posee poca información hay dos vías: buscar CD 1.12 2.01 1.6
casos similares para comparar al menos cualitativamente y hacer un análisis intrínseco de los detalles del flujo.
Esto último consiste en comprobar que las estelas, capas límites, zonas de presión y turbulencia, gradientes de
velocidad, etc. están en los sitios adecuados; se desprenden, ensanchan y estrechan en los momentos oportunos;
giran en los sentidos adecuados; mantienen unos órdenes de magnitud razonables, etc.
43
El coche. Como no se tienen datos experimentales, se va a comparar con otros casos similares y analizar
algunos los detalles del flujo obtenido.
El coeficiente de arrastre de un coche suele estar en torno a 0.5 - 0.35. La simulación arroja un valor de
1.82.
Lo primero que se puede comentar es que la mayor parte de ese coeficiente, 1.803, es debida al arrastre
de presión (desprendimiento de la capa límite, estela...) y sólo 0.017 a la fricción viscosa; Para las formas y ta-
maños de superficies que se tienen, esto es cualitativamente correcto. Lo segundo es que se ha simulado un caso
bidimensional y, en éstos, los coeficientes de arrastre son mucho mayores. En la bibliografía se pueden obtener
valores típicos para algunas formas, bastante independientes del número de Reynolds debido a las esquinas.
Aunque estos datos son para flujo libre, sin carretera debajo, se puede considerar que el coeficiente obtenido es
bastante razonable.
Una pregunta que surge en este momento es ¿para qué sirve una simulación bidimensional si se aparta
tanto de la realidad tridimensional?. En primer lugar no siempre se aparta tanto, a veces son bastante buenos
modelos. En segundo, aunque los valores no son correctos, las mejoras que se puedan conseguir en una geome-
tría bidimensional, casi seguro que serán mejoras tridimensionales. En tercero, solo se está empezando con el
CFD.
Pasando ahora al análisis intrínseco del flujo, la capa límite, dada la geometría y los pocos nodos sobre
la superficie del coche, no hay ningún sitio donde se haya desarrollado de manera adecuada.
Sí se puede observar que la estela se desprende justo en las esquinas, lo que es correcto dada la brus-
quedad del cambio de pendiente.
7.11e+01
7.11e+01
6.40e+01 6.40e+01
5.70e+01 5.70e+01
4.99e+01 4.99e+01
4.29e+01 4.29e+01
3.58e+01 3.58e+01
2.88e+01 2.88e+01
2.17e+01 2.17e+01
1.47e+01 1.47e+01
7.61e+00 7.61e+00
5.56e-01 5.56e-01
También quedan bien definidas en la estela las dos zonas de recirculación, en sentido contrario en la
parte superior e inferior. Esto se puede ratificar con la vorticidad (dibujo inferior izquierdo).
3.35e+02 7.11e+01
2.71e+02 6.40e+01
2.06e+02 5.70e+01
1.41e+02 4.99e+01
7.70e+01 4.29e+01
1.25e+01 3.58e+01
-5.20e+01 2.88e+01
-1.17e+02 2.17e+01
-1.81e+02 1.47e+01
-2.46e+02 7.61e+00
-3.10e+02 5.56e-01
Las zonas de estancamiento en el morro y en la esquina inferior del parabrisas son claramente visibles
(dibujo superior derecho), aunque esta última está dentro del desprendimiento de la esquina superior del morro.
Así mismo, en la parte superior del parabrisas se produce un desprendimiento por el cambio brusco de forma.
La distribución de presión estática muestra el aumento de presión en la zona de estancamiento del mo-
rro y la disminución en las esquinas con desprendimiento (dibujo inferior izquierdo).
44
En la estela (dibujo superior derecho), la presión es prácticamente atmosférica, lo que es coherente con que el arras-
tre por presión en los cuerpos no se produce por depresión en la parte posterior sino por falta de recuperación de la
presión.
1.92e+03 1.92e+03
1.41e+03 1.41e+03
9.12e+02 9.12e+02
4.08e+02 4.08e+02
-9.47e+01 -9.47e+01
-5.98e+02 -5.98e+02
-1.10e+03 -1.10e+03
-1.60e+03 -1.60e+03
-2.11e+03 -2.11e+03
-2.61e+03 -2.61e+03
-3.11e+03 -3.11e+03
En la parte inferior del coche, la presión es ligeramente inferior a la atmosférica, a causa de la separa-
ción en la esquina inferior del morro y a las mayores velocidades del flujo en esa zona.
1.92e+03
Esta depresión es la responsable de la fuerza
vertical negativa (de adherencia a la carretera) que se ha 1.41e+03
obtenido. Este efecto se ve amplificado cuando se utili- 9.12e+02
zan "spoilers" en la parte inferior del morro. 4.08e+02
-9.47e+01
Para conseguir mayor adherencia también se -5.98e+02
emplean alerones. Los realmente eficaces están un poco
-1.10e+03
separados de la carrocería en la parte superior y tienen
-1.60e+03
forma aerodinámica (perfil de ala invertido). Los alero-
nes en la parte posterior del techo aumentan la adheren- -2.11e+03
cia pero a costa de un gran arrastre. El mismo alerón en -2.61e+03
la parte posterior del maletero tiene una influencia prác- -3.11e+03
ticamente nula por estar situado en la estela. Caso aparte son los alerones en las cabinas de los camiones.
Algunos cambios que se sugieren para mejorar la simulación son:
- Ampliar la zona de cálculo alrededor del coche hasta que su propia presencia no afecte las condi-
ciones en los bordes.
- Mallado más detallado en la superficie del coche.
- Redondear las esquinas de la geometría. Por supuesto, con unas formas mucho más suaves se puede
conseguir una reducción sustancial del coeficiente de arrastre. Un detalle que hay que tener en
cuenta es que, si la parte posterior del techo es muy suave y no se produce el desprendimiento en
una esquina sino en mitad de una superficie curva, el punto de desprendimiento (y el coeficiente de
arrastre, por tanto) va a estar muy influenciado (pero mucho...) por el modelo de turbulencia y los
parámetros elegidos.
45
TOBERA CONVERGENTE DIVERGENTE
Ejemplo de cálculo de flujo compresible transónico utilizando el Fluent.
INTRODUCCIÓN
Este ejemplo consiste en el cálculo de una tobera convergente divergente del tipo de las que se pueden
utilizar en los motores de cohetes. Este tipo de toberas están diseñadas para que partiendo de un flujo subsónico
en la entrada, 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, llegar a sónico en la
garganta y permanecer subsónico después o llegar a sónico y 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 sali-
da. El Fluent con el solver de tipo segregated no es capaz de resolver problemas en los que la entrada o salida
sea supersónica (para estos problemas hay que utilizar el coupled), pero el flujo si puede serlo en el interior. En
este caso se va a simular solamente la tobera con entrada y salida subsónicas.
La tobera tiene un área en la garganta de 300 cm2 y en la salida de 600 cm2. El fluido es aire. La presión
a la salida (absoluta) 100 kPa. La temperatura total a la entrada 3000 K (se supone que viene de una cámara de
combustión). La presión total a la entrada se va a ajustar para buscar los distintos casos.
Con estos datos, la teoría unidimensional de flujo supersónico sin viscosidad, etc., indica que a partir de
una presión total (abs) de 106.7 kPa en la entrada, en la garganta se alcanza flujo sónico. En estas condiciones el
flujo másico es de 2.36 kg/s, la temperatura a la salida de 2952 K y el nº de Mach a la salida de 0.306. Para una
presión de 1065 kPa (abs) el flujo sería supersónico en toda la zona divergente. Con valores intermedios se for-
ma una onda de choque y el flujo sale subsónico.
Se va a simular una tobera de sección circular de forma bidimensional, para ello sólo se generará la
geometría de una sección circunferencial y se utilizará el modelo axisimétrico.
GEOMETRÍA
Con el GAMBIT se traza el contorno de la sección. La curvatura de la parte exterior se define de forma
que la entrada y salida sean lo más axiales posible. Los puntos son:
PUNTO X Y
A 0 0
B 0 0.15
C 0.15 0.14
D 0.3 0.09772 Garganta; radio para sección 300 cm2
E 0.5 0.12
F 0.9 0.1382 Salida; radio para sección 600 cm2
G 0.9 0
Y las líneas definidas con ellos:
NOMBRE PUNTOS TIPO DE LÍNEA
AB A,B Straight
BCDEF B,C,D,E,F URBS
FG F,G Straight
GA G,A Stright
Con estas líneas se define una superficie.
46
MALLADO
Primero se ponen nodos en las líneas:
NOMBRE Nº NODOS
AB 8
BCDEF 40
FG 8
GA 40
Después se genera el mallado de la superficie con celdas triangulares, basándose en los nodos de las lí-
neas. El mallado generado tiene un número de celdas de 538. El número de celdas es bajo, pero esto es sólo un
ejemplo.
Líneas:
AB pressure-inlet entrada
FG pressure-outlet salida
GA axis eje de simetría
BCDEF wall pared
La superficie se define como fluid.
MODELO
Pasando ahora al Fluent, se lee el mallado y se definen los modelos y el fluido que se van a utilizar.
Define - Models:
- Solver:
- Axisymmetric. El dominio axisimétrico permite realizar el cálculo de la tobera de sección
circular de forma bidimensional, pero teniendo en cuenta la corrección necesaria de las
áreas.
- Steady
- Energy: Activado.
- Viscous: Turbulent: Standard k-epsilon. Con los parámetros por defecto.
- Viscous heating activado. Las velocidades van a ser elevadas, aunque, dada la longitud de
la tobera, su efecto va a ser pequeño.
Define - Materials:
Fluido: aire (por defecto). Para la densidad hay que escoger la opción Ideal Gas Law
Define - Operating Conditions:
Dado que la presión en la salida es de 100000 Pa, se fija la Operating Pressure en ese valor, aunque no
es obligatorio. Téngase en cuenta que esta presión es con la que el programa define las presiones relativas (gau-
ge) que utiliza en las condiciones de contorno y en las representaciones; Cuando se realizan cálculos manuales
de flujo supersónico habitualmente se trabaja con presiones absolutas pues son las que hay que poner en las
fórmulas termodinámicas.
47
Para flujo transónico, las condiciones de entrada y salida deben ser presión total y presión estática res-
pectivamente. La pared se va a definir sin rugosidad (aunque sigue existiendo la condición de no deslizamiento)
y sin transferencia de calor, lo que la convierte en adiabática. El eje de simetría es el eje del dominio axisimétri-
co.
Define - Boundary Conditions:
- pared: Wall Rougness: Roughness Height: 0 m. Thermal Conditions: Heat Flux: 0 w/m2.
- salida: Gauge Pressure: 0 Pa (100000 Pa abs). Backflow Total Temperature: 3000 k (con la espe-
ranza de que no la tenga que utilizar mucho, no debería haber flujo inverso en la salida). Backflow
Turb: por defecto.
- entrada: Gauge Total Pressure: se comenta luego. Total Temperature: 3000 K. Direction Specifica-
tion Method: Direction Vector: X = 1, Y = 0; la dirección de entrada según el eje x, dirección axial.
Turbulence: por defecto. Supersonic/Initial Gauge Pressure: igual que Gauge Total Pressure.
La presión total a la entrada (fijada una presión estática a la salida) fija la condición de funcionamiento
de la tobera. Se ha ejecutado el programa con tres presiones distintas: 5000 Pa (105000 abs) que mantiene el
flujo subsónico en toda la tobera; 20000 Pa (120000 abs) con la que se llega a condiciones sónicas en la tobera y
subsónicas en el resto y 80000 Pa (180000 abs) donde el flujo es supersónico en parte de la zona divergente,
produciéndose una onda de choque antes de la salida.
PARÁMETROS DE LA SOLUCIÓN
En este caso no se puede empezar con valores nulos en las variables porque se produce una divergencia
rápida en las primeras iteraciones y el cálculo rompe (si no diverge demasiado rápidamente no se colgará el
ordenador, aunque todo podría ser). En el menú Solve - Initialize - Initialize se puede utilizar Compute From:
pressure inlet-2, que pone la presión y temperatura a los valores de la entrada. También es bueno añadir una
cierta velocidad en la dirección x (200 m/s…) para facilitarle el trabajo y evitar que en las primeras iteraciones
exista flujo inverso en la salida.
En la sección Solve - Controls se dejan los valores por defecto.
En Solve - Monitors - Residuals se puede activar el dibujo de los residuos. Los valores de los criterios
de convergencia se dejan por defecto.
El caso subsónico converge en pocas iteraciones. Los casos transónicos necesitan unos cientos de ellas.
ANÁLISIS DE LOS RESULTADOS

Los dibujos de velocidad en forma de vectores no dan mucha información en este caso, porque el flujo
es bastante axial. En un difusor con mayor gradiente son útiles para comprobar si se produce desprendimiento.
Vectores de velocidad, 105 kPa
5.91e+02
5.52e+02
5.13e+02
4.74e+02
4.34e+02
3.95e+02
3.56e+02
3.16e+02
2.77e+02
2.38e+02
1.99e+02
48
Los mapas del número de Mach (Velocity – Mach number) son bastante explíctos. En el de 105 kPa se
puede ver el aumento hasta la garganta y la disminución posterior. La difusión es mayor cerca de la pared que en
el eje.
Mach 105 kPa

5.51e-01
5.14e-01
4.77e-01
4.40e-01
4.03e-01
3.66e-01
3.30e-01
2.93e-01
2.56e-01
2.19e-01
1.82e-01
En el de 120 kPa se puede ver cómo en la garganta se alcanza flujo sónico en la totalidad de la sección;
en la pared los valores son un poco mayores que la unidad. En el difusor la velocidad vuelve a ser subsónica.
Mach 120 kPa
1.15e+00
1.06e+00
9.69e-01
8.78e-01
7.87e-01
6.96e-01
6.05e-01
5.14e-01
4.23e-01
3.32e-01
2.41e-01
Con la presión de 180 kPa en la entrada, el flujo es totalmente sónico en la garganta y pasa a supersóni-
co en la zona divergente. Hacia la mitad del difusor se forma una onda de choque y el número de Mach baja a
subsónico, disminuyendo posteriormente un poco hasta la salida.
Mach 180 kPa

2.04e+00
1.86e+00
1.68e+00
1.50e+00
1.32e+00
1.14e+00
9.60e-01
7.81e-01
6.01e-01
4.22e-01
2.42e-01
En las gráficas de presión estática se puede ver el mismo comportamiento: disminución de la presión
desde la entrada hasta la garganta, debido a la aceleración del flujo y, para 105 y 120 kPa, aumento hasta la sali-
49
da. En la gráfica de 180 kPa se observa cómo la presión sigue disminuyendo a partir de la garganta en la zona
supersónica, aumentando bruscamente al atravesar la onda de choque.
Presión estática 105 kPa

2.70e+03
8.05e+02
-1.09e+03
-2.98e+03
-4.87e+03
-6.76e+03
-8.65e+03
-1.05e+04
-1.24e+04
-1.43e+04
-1.62e+04

1.54e+04
8.62e+03
1.79e+03
-5.04e+03
-1.19e+04
-1.87e+04
-2.55e+04
-3.24e+04
-3.92e+04
-4.60e+04
-5.28e+04

7.31e+04
5.76e+04
4.20e+04
2.65e+04
1.10e+04
-4.54e+03
-2.01e+04
-3.56e+04
-5.11e+04
-6.66e+04
-8.22e+04
El mapa de temperatura total en el último caso muestra unos valores bastante constantes hasta la onda
de choque, disminuyendo sustancialmente después. La temperatura total da una indicación de la pérdida de ener-
gía que tiene lugar en la onda de choque.
50
Temperatura total 180 kPa

3.00e+03
2.98e+03
2.96e+03
2.94e+03
2.93e+03
2.91e+03
2.89e+03
2.87e+03
2.85e+03
2.83e+03
2.81e+03
Se pueden apreciar mejor las características antes comentadas en los gráficos XY. En la figuras siguien-
tes se representa el número de Mach y la presión estática a lo largo del eje de simetría. En la primera se observa
cómo los valores de 120 y 180 kPa son iguales en la zona convergente. Se podría observar el mismo comporta-
miento en las presiones si en lugar de regular la tobera con la presión total a la entrada se hubiera hecho redu-
ciendo la presión estática a la salida. Una vez que se alcanza el flujo sónico en la garganta, la zona convergente
quedaría bloqueada.
105 kPa 105 kPa
120 kPa 120 kPa
180 kPa 180 kPa
2.2 8e+04
2 6e+04
1.8
4e+04
1.6
2e+04
1.4
1.2 0
Mach Presión
1 estática -2e+04
0.8
-4e+04
0.6
-6e+04
0.4
0.2 -8e+04
0 -1e+05
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
X X
El flujo másico y la temperatura estática media a la salida para los tres casos son:
Caso Caudal másico Temperatura
105 kPa 1.74 kg/s 2963 K
120 kPa 2.6 kg/s 2883 K
180 kPa 3.92 kg/s 2705 K
Si se compara la presión total del segundo caso con el valor que se ha calculado por medio de la teoría
unidimensional no viscosa (120 frente a 106.7 kPa), puede parecer que la diferencia es sustancial, pero no lo es
tanto si se considera que en el caso simulado, aparte de ser bidimensional y tener en cuenta la viscosidad y la
turbulencia, el nº de Mach que se alcanza en la pared –donde el flujo se ve acelerado debido a la curvatura- es
superior a la unidad.
La presión estática en la entrada no coincide con la total ya que la sección de entrada no es infinita y
existe una cierta velocidad –bastante elevada, por cierto-.
La onda de choque real es bastante más brusca que la simulada; el frente de onda tiene un espesor real-
mente pequeño. Se podría conseguir más precisión refinando el mallado en esa zona.
Como ya se ha comentado anteriormente, para el cálculo de flujos compresibles es más adecuado el sol-
ver de tipo coupled en vez de segregated y, dada la brusquedad de las variaciones existentes en las variables,
también es muy conveniente utilizar una discretización de segundo orden para las variables, pero eso es otra
historia.
51
TRATAMIENTO DE GRÁFICOS CON GAMBIT Y FLUENT

GAMBIT
Los gráficos generados con Gambit se pueden enviar a un fichero mediante el procedimiento siguiente:
File
Print Graphics
Se escoge el formato de fichero gráfico. Lo
normal será elegir los formatos TIFF o BMP. En File
Name se escribe el nombre que se le quiera dar al fiche-
ro.
FLUENT
Desde el Fluent se pueden imprimir y copiar
los gráficos a otros documentos. Si se pincha en la parte
superior izquierda de la ventana, aparece un menú con
las opciones habituales, más otras tres para imprimir y
copiar los dibujos.
Print... envía el dibujo a la impresora
Copy to clipboard... lo pone en el portapapeles
Page Setup... define las opciones de impresión
El procedimiento para incluir un gráfico en un
documento de Word es el siguiente: definición de las
opciones (con Page Setup); copiar al portapapeles, y
pegarlo desde el Word.
Son interesantes las opciones de guardar en el
portapapeles como Metafile con la que se consigue una
mejor definición, sobre todo de texto, y la de Reverse
Foreground/Background que hace que el fondo de las figuras sea blanco y
no negro (para ahorrar tóner). Si se va a imprimir en blanco y negro,
conviene utilizar Gray Scale. De esta forma, además de utilizar sólo tonos
grises, gradúa los tonos de más claro a más oscuro. Los gráficos con
vectores de velocidad se pueden imprimir con la opción Monochrome,
pues la magnitud de la velocidad viene ya definida por la longitud de los
vectores.
Otra opción disponible en el Fluent es

guardar los gráficos en ficheros para insertarlos
posteriormente en un documento Word. Para ello,
en el menú de Fluent se realiza lo siguiente:
File
Hardcopy
con lo que aparece la ficha siguiente en la que se
pueden selecionar las mismas opciones que en la
ficha Page Setup, pudiendo almacenar el gráfico
en una gran variedad de formatos.
52
REDACCIÓN DE INFORMES
Una vez realizadas todas las sesiones de prácticas de simulación numérica, debe presentarse un informe
escrito sobre cada uno de los casos resueltos en dichas sesiones. En este apartado se dan algunas indicaciones
para redactar un informe correcto.
Un buen informe debe contener toda la información necesaria para la comprensión de un determinado
problema, describiéndolo perfectamente en todas sus fases: planteamiento, solución adoptada, resultados,
conclusiones, etc. En palabras sencillas, un informe bien redactado debe cumplir la condición siguiente: tras
leerlo, alguién sin conocimientos previos sobre el tema podrá comprenderlo y utilizar sus conclusiones, sin
necesidad de recurrir a otra información adicional
Características generales de un buen informe
A continuación se resumen las características generales que debe reunir un buen informe técnico sobre
cualquier campo o tema:
− Corrección técnica: la solución adoptada al problema planteado debe ser correcta y estar
adecuadamente justificada.
− Adecuado tratamiento de la información gráfica, de forma que sea legible y comprensible
para cualquier persona con una mínima formación técnica (y sin necesidad de utilizar una lupa)
− Presentación en un soporte adecuado: como mínimo, debe presentarse en papel, dentro de un
portafolios o carpeta que sujete las hojas impidiendo su pérdida (¡Ojo! Un simple clip o grapa no cum-
ple esta función). Las diferentes hojas que lo componen deben estar numeradas, apareciendo en la pri-
mera de ellas el título del trabajo, sus autores y la fecha de realización. No es obligatorio utilizar el or-
denador para redactarlo, pero si se escribe a mano debe ser con tinta azul o negra: no es admisible pre-
sentar un informe a lápiz.
− Corrección gramatical y ortográfica.
− Brevedad y concisión: debe incluirse solo la información relevante al caso, eliminando aspec-
tos superfluos o referencias que pueden encontrarse en otros documentos. Estas referencias deben do-
cumentarse adecuadamente, incluyendo un apartado de bibliografía al final del informe.
− Originalidad: el informe debe aportar información útil y novedosa sobre el problema resuelto.
No resulta adecuado, por tanto, copiar en el informe apartados completos de los apuntes o de un libro.
Contenidos mínimos de los informes de estas prácticas
Como mínimo, el informe de cada uno de los casos resueltos en estas prácticas deberá constar de los
siguientes apartados:
1. Objetivos
En este apartado se describirá con claridad el problema planteado y los objetivos que se persiguen en su
resolución.
2. Definición de la geometría
En este aparado se incluirá una lista de los puntos que se han utilizado para introducir la geometría del
problema en el GAMBIT. También se incluirá un dibujo de la geometría y el mallado del problema, que puede
realizarse desde el FLUENT con la opción “Display Grid”.
53
3. Parámetros de resolución
En este apartado se incluirán todos los parámetros que se hayan utilizado en la resolución del problema,
justificando su elección: modelo utilizado, propiedades de los materiales, condiciones de contorno, valores de
referencia (para el cálculo de coeficientes), parámetros de iteración, etc.
4. Exposición y análisis de resultados

En este apartado se incluirán los gráficos más importantes que ilustren la resolución del problema, con
vistas generales y ampliaciones de ciertas zonas en caso necesario. Los gráficos servirán de soporte a explicacio-
nes breves sobre los fenómenos fluidodinámicos más relevantes observados.
5. Conclusiones
En este apartado se citarán las conclusiones que se pueden extraer de la realización del trabajo: no se
trata de repetir párrafos de apartados anteriores, sino de repetir lo más importante.
6. Bibliografía
En este apartado se citará el material de apoyo (libros, apuntes, artículos, www, ...) que se han consultado
y utilizado en la elaboración del informe.
54

Fluent Manejo Del Program

Cargado por

Copyright:

Formatos disponibles

Fluent Manejo Del Program

Cargado por

Información del documento

Descripción 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

Fluent Manejo Del Program

Cargado por

Copyright:

Formatos disponibles

Área de Mecánica de Fluidos. Universidad de Oviedo.

ÁREA DE MECÁNICA DE FLUIDOS

DINÁMICA DE FLUIDOS COMPUTACIONAL (CFD)

OBJETIVOS DE ESTAS PRÁCTICAS

Con la realización de las prácticas de simulación numérica de la asignatura Mecánica

- Familiarización con las técnicas de simulación numérica aplicadas a la Mecánica de Fluidos

- Exposición de forma adecuada de los resultados de una simulación numérica, interpretándo-

- Mejora de la confección y presentación de informes científico-técnicos.

- Realización de trabajos en equipo.

Algunas ventajas e inconvenientes

Explicación resumida del método

RESOLUCIÓN NUMÉRICA DE LAS ECUACIONES DE FLUJO

Ecuaciones que describen el movimiento de un fluido

La turbulencia no es un fenómeno siempre indeseable. En los cilindros de un motor de combustión in-

DISCRETIZACIÓN DE LAS ECUACIONES DEL FLUJO

Método de los volúmenes finitos

Supóngase que se quiere discretizar la derivada segunda siguiente:

Resolución de las ecuaciones discretizadas

Dependencia del tiempo y no linealidad

AL EMPEZAR CON EL SOFTWARE

FICHEROS QUE SE GENERAN

De cómo utilizar el GAMBIT para generar figuras geométricas, crear el ma-

En la parte superior se tiene un menu principal con las siguientes opciones:

ESPECIFICAR TIPOS DE CONDICIONES DE CONTORNO

ALGUNAS PRECISIONES SOBRE GEOMETRÍA Y MALLADO

ALGUNOS CONSEJOS E INDICACIONES

// de la misma forma se introduce el resto de los puntos

// Definición de las líneas

En el dibujo mostrado se ha empleado la opción ARC en lugar de NURBS.

// Creación de la superficie externa:

// Poner nombre en Label.

// Substracción (dibujo superior derecho)

// Generación de los nodos

De Operation selecciónese el icono mallado (Mesh)

// Mallado de la línea ABC

// Se repite el proceso para el resto de las curvas

// Generación del mallado

De Mesh selecciónese el icono superficies (Face)

Se selecciona la superficie (TUNEL).

Al seleccionar Apply en la parte

// Especificación de los contornos (dibujo inferior izquierdo)

En el programa FLUENT 5.1 se hablará acerca de estas condiciones.

De cómo utilizar el Fluent para calcular una posible solución de flujo. De

> Reading "E:\Gambit\coche.msh"...

- Zone name: se puede cambiar el nombre.

- Gauge Total Pressure: presión total (relativa) de entrada.

Condiciones correctas de entrada – salida

Symmetry & Axis

Entonces, ¿qué CONDICIONES hay que PONER con el coche?

CRITERIO DE CONVERGENCIA Y MONITORIZACIÓN

Representa un coeficiente de fuerza (arrastre, sustentación o mo-

Para el arrastre hay que indicar la velocidad de referencia, la densi-

En Number of Iterations se indica el número de iteracio-

VER LOS RESULTADOS

Display Vectores de velocidad

(Velocity – Stream Function, sin Filled Turbulence - Turbulent Kinetic Energy

Distribución de presión sobre la superficie del coche:

Fuerzas en la dirección horizontal y coeficiente de arrastre:

Fuerzas en la dirección vertical y coeficiente de sustentación: