Grado en Ingeniería Mecánica

Trabajo Fin de Grado

Caracterización de gotas en
atomizadores pressure-swirl
PRESSURE-SWIRL

Autor

Álvaro Lloret Sáez

Tutor/a en Empresa

Alicia Muñoz Agulló

Tutor/a en UPV

Marcos Carreres Talens

25 de mayo de 2024
RESUMEN

RESUMEN

BLA BLA BLA

Palabras Clave:

Página 1 de 57
RESUM

RESUM

BLA BLA BLA

Paraules Clau:

Página 2 de 57
ABSTRACT

ABSTRACT

BLAH BLAH BLA

Key Words:

Página 3 de 57
TABLA DE CONTENIDOS

Tabla de Contenidos

I Memoria 5

II Pliego de Condiciones 54

III Presupuesto 56

Página 4 de 57
Parte I

Memoria
MEMORIA

Índice de la Memoria

1 Objetivo del proyecto 10

1.1 Objeto y justificación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.2 Antecedentes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.3 Estructura del proyecto . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.3.1 Memoria . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.3.2 Pliego de condiciones . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.3.3 Presupuesto . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13

2 Marco Teórico 14
2.1 Turbinas de gas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.1.1 Cámara de combustión . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.2 Inyectores . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.3 Combustión . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.3.1 Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.3.2 Combustión de Spray . . . . . . . . . . . . . . . . . . . . . . . . . . . 17

3 Atomización 19
3.1 Definición de Atomización . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
3.2 Factores que influencian en la atomización . . . . . . . . . . . . . . . . . . . 19
3.2.1 Propiedades del liquido . . . . . . . . . . . . . . . . . . . . . . . . . . 20
3.2.2 Condiciones Ambiente . . . . . . . . . . . . . . . . . . . . . . . . . . 20
3.2.3 Números adimesnionales . . . . . . . . . . . . . . . . . . . . . . . . . 21
3.3 Clasificación atomizadores . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
3.3.1 Atomizador tipo Pressure-Swirl . . . . . . . . . . . . . . . . . . . . . 23
3.4 Procesos Básicos de la atomización . . . . . . . . . . . . . . . . . . . . . . . 25
3.4.1 Formación estática de gotas . . . . . . . . . . . . . . . . . . . . . . . 26
3.4.2 Rotura de gotas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
3.4.2.1 Rotura de gotas en aire en movimiento . . . . . . . . . . . . 27
3.4.3 Desintegración de lámina líquida . . . . . . . . . . . . . . . . . . . . 28
3.5 Características de los sprays . . . . . . . . . . . . . . . . . . . . . . . . . . . 29

4 Herramientas utilizadas. CFD y software Basilisk 33

4.1 CFD . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
4.1.1 ¿Que es el CFD? . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
4.1.2 Funcionamiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
4.1.3 Ecuaciones del fluido . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
4.1.3.1 Ecuación de conservación de la energía . . . . . . . . . . . . 37
4.1.4 Discretización . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
4.1.5 Solucionar de las ecuaciones discretizadas . . . . . . . . . . . . . . . . 39
4.1.6 Modelos de turbulencia . . . . . . . . . . . . . . . . . . . . . . . . . . 40
4.1.7 Direct numericla simulation, DNS . . . . . . . . . . . . . . . . . . . . 41

Página 6 de 57
MEMORIA

4.2 Basilisk . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
4.2.1 Estructura general código Basilisk . . . . . . . . . . . . . . . . . . . . 43
4.2.2 Softwares adicionales . . . . . . . . . . . . . . . . . . . . . . . . . . . 45

5 Desarrollo del trabajo 48

5.1 Geometría del atomizador . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
5.2 Condiciones iniciales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
5.2.1 Dominio y mallado . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
5.2.2 Condiciones de contorno . . . . . . . . . . . . . . . . . . . . . . . . . 52

II Pliego de Condiciones 54

III Presupuesto 56

Página 7 de 57
MEMORIA

Índice de las Figuras de la Memoria

2.1 Esquema TG [?] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14

2.2 ciclo real sencillo Bryton (no regenerativo): (a) en el palno pv; (b) en el plano
Ts ó hs.[?] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.3 Esquema de un atomizador pressure-swirl[?]. . . . . . . . . . . . . . . . . . . 16
2.4 Flujo de atomizador pressure-swirl totalmente desarrollado [?]. . . . . . . . . 17

3.1 Ejemplos de atomizadores de rotación a)guías rectas, b)guías curvas[?]. . . . 22

3.2 Ejemplos de atomizadores según su tipo a)Atom. de presión, b)Atom. de ro-
tación c) Atom. de doble fluido[?]. . . . . . . . . . . . . . . . . . . . . . . . . 23
3.3 Tipos de atomización en atopmizadores tipo pressure-swirl[?]. . . . . . . . . 24
3.4 Esquema de la forma de un atomizador tipo pressure-swirl simplex[?]. . . . . 25
3.5 Fotografías del desarrollo del flujo de un atomizador swirl al ir aumentando
la presión desde 0 [?] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
3.6 Ilustración tipos de gotas[?] . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.7 Fotografía que ilustra el modo wave principalmente [?] . . . . . . . . . . . . 29
3.8 Ejemplo definición ángulo de cono en inyector diésel [?]. . . . . . . . . . . . . 31
3.9 Ejemplo PDF de que muestra la probabilidad de encontrar gotas, en función
de su diámetro, a diferentes distancia del orificio de salida. [?] [?]. . . . . . . 32

4.1 Volumen mínimo, partícula de fluido [?]. . . . . . . . . . . . . . . . . . . . . 36

4.2 Malla 2D o subdivisión de volumen de métodos finitos [?]. . . . . . . . . . . 38
4.3 Ejemplo discretización elemento unidimensional [?]. . . . . . . . . . . . . . . 39
4.4 Ejemplo código #include . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
4.5 Ejemplo código #define . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
4.6 Ejemplo código condiciones de contorno . . . . . . . . . . . . . . . . . . . . 45
4.7 Ejemplo código función de lanzamiento . . . . . . . . . . . . . . . . . . . . . 46
4.8 Ejemplo función para obtener capturas de imagen de la simulación. . . . . . 47

5.1 Corte lateral atomizador pressure-swirl . . . . . . . . . . . . . . . . . . . . . 48

5.2 Atomizador tipo presure swirl comercial, Danfoss OD-H 030H8103 oil nozzle 49
5.3 Esquema datos geométricos del atomizador Danfoss OD-H 030H8103 oil nozz-
le . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
5.4 Ejemplo estructura de árbol 2D, llamada quadtree al dividir en 4 en vez de 8
[?] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52

Página 8 de 57
MEMORIA

Índice de Tablas de la Memoria

5.1 Tabla de datos geométricos del atomizador atomizador Danfoss OD-H 030H8103
oil nozzle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51

Página 9 de 57
Objetivo del proyecto

En la actualidad el estudio del proceso de atomización se ha vuelto muy importante para el

desarrollo y mejora de sistemas de inyección en sistemas de combustión. Lo que se busca con
estos estudios es aumentar el rendimiento de las maquinas térmicas que de las que forman
parte, así como la reducción de gases contaminantes mejorando y optimizando la combustión
mediante una correcta mezcle de aire/combustible.

En adición a su importancia en sistemas de combustión se suma la utilidad de este pro-

ceso en otras múltiples áreas, como pueden ser la medicina o la prevención de incendios.

Debido a la importancia que tiene este proceso han sido cada vez más los métodos de estudio
que se han desarrollado para poder caracterizarlo. No obstante, la complejidad que tiene el
desarrollo de métodos experimentales para estudiarlo, debido a tratarse de flujo bifásico, ha
promovido el desarrollo de métodos de cálculo para su caracterización.

Hoy en día gracias al uso de de los avances en computación se ha desarrollado la diná-

mica de fluidos computacional o CFD (computational fluid dynamics) que permite el estudio
de estos fenómenos de manera menos compleja y con diferentes grados de precisión en fun-
ción del objetivo del mismo.

El uso de técnicas de mayor resolución conocidas como DNS permiten la caracterización

de las gotas que se generan durante la atomización , pudiendo ver como varían las mismas
en función de los propiedades físicos del combustible y el medio al que sobre el que se inyecta.

Gracias a esto se va a desarrollar un modelo de para un atomizador de tipo preasure-swirl

con una características geométricas y de funcionaminto predefinidas, que nos permita la ca-
racterización de las gotas que se generan en función de las propiedades físicas del combustible
usado.
MEMORIA

1.1. Objeto y justificación

La combustión de combustibles líquidos en motores diésel, en turbinas de aerorreacción, o
en encendedores de chispa dependen en gran medida de una correcta atomización del com-
bustible para un óptimo funcionamiento y para que la combustión se lo mas limpia posible.

Como se demuestra en diversas publicaciones, la atomización no depende exclusivamente

de un único factor, siendo muchos los que influyen para romper la vena líquida del fluido,
ni influyendo los mismos factores ni de igual manera en los diversos tipos de atomizadores,
pudiendo afectar no solos las condiciones en las que se inyecta el combustible sino que tam-
bién pueden afectar las propiedades físicas del combustible o del fluido en el que se inyecta.

Sabiendo lo anterior en este Trabajo de Fin de Grado se ha optado por estudiar la variación
en las características de la población de gotas,en un atomizador de tipo preasure-swirl, me-
diante la variación de propiedades físicas del combustible, de manera individual, sin variación
en las condiciones del fluido en el que se inyecta o en las condiciones en las que se inyecta.

1.2. Antecedentes
En los últimos años ha habido y hay diversos proyectos y trabajos académicos que intentan
aclarar como funciona correctamente la atomización y como le influyen diferentes parámetros
físicos.

Hay algunos que sobre el mismo tema, se centran más en la validación experimental del
código de Salvador F.J., Martí-Aldaraví P., Lozano A., Taghavifar H. y,Nemati A. [?], usan-
do en ambos casos DNS.

En otros estudios que usan cálculos de flujo para atomización se estudia el efecto de otros
parámetros como es el caso del trabajo de fin de grado de Lozano A. [?], donde estudia el
efecto de la turbulencia en la inyección en toberas.

Los trabajos y artículos aquí mencionados, siguen un objetivo común, todos buscan au-
mentar la comprensión sobre el proceso de atomización y como este es afectado por diversos
parámetros y a la vez como este afecta a otros procesos, entre otros el de combustión, no
obstante estos son distintos entre ellos ya sea en los parámetros elegidos para el estudio o en
lo que se busca obtener en cada uno de ellos.

Página 11 de 57
MEMORIA

1.3. Estructura del proyecto

El proyecto actual se estructura en tres documentos o partes totalmente diferenciables entre
si. El proyecto esta conformado por una Memoria, un Pliego de Condiciones y un
Presupuesto. No se ha considerado necesario un apartado de planos por la naturaleza del
trabajo.

1.3.1. Memoria
Este documento engloba el desarrollo del proyecto en su totalidad. Se incluye un marco
teórico para facilitar la comprensión del mismo y darle contexto, así como los procesos que
se han seguido durante le proyecto y los resultados y conclusiones obtenidos del mismo. Este
documento se subdivide a su vez en las siguientes 7 partes:

1. Objetivo del proyecto: Se hace una breve introducción al proyecto y se pone en

contexto el mismo explicando la importancia del estudio del proceso de atomización y
los parámetros que afectan al proceso.

2. Marco Teórico: Se da contexto teórico sobre los procesos y aspectos que afectan al
proyecto mediante breves explicaciones teóricas. Algunos de los temas son la inyección
de combustible, la combustión por difusión y turbinas para motores de aviación.

3. Atomiazación: Se indaga más en la teoría conocida sobre el proceso de atomización

y las propiedades físicas que le afectan para la correcta comprensión del trabajo.

4. Herramientas utilizadas: CFD y software Basilisk: Se explica de manera clara

y detallada en que consiste el CFD (Computational Fluid Dynamics) y se explica el
porque de la utilización como la herramienta principal, el software Basilisk.

5. Descripción del trabajo: Se explica como se ha procedido en el trabajo para obtener

los resultados que se buscaban analizar.

6. Análisis de resultados: Se exponen los resultados que se han obtenido de manera

clara y concisa, junto al análisis de los mismos.

7. Conclusiones y futuros desarrollos: Tras el análisis de los resultados se aportan

conclusiones que sirvan de base para futuros estudios y proyectos sobre el mismo tema
o similares.

1.3.2. Pliego de condiciones

En este documento se plantean las limitaciones y exigencias que se han tenido durante la
realización del proyecto, siendo estas de tipo técnico y legal, añadiendo también información
relacionada con el entorno y el material usado.

Página 12 de 57
MEMORIA

1.3.3. Presupuesto
En este último documento se va exponer el coste teórico de la realización del trabajo, tanto
por horas de mano de obra, como por material físico y digital utilizado, llegando al valor
económico hipotético de un proyecto de características similares a las de este.

Página 13 de 57
Marco Teórico

En esté apartado se va a proceder a establecer un contexto teórico sobre el proyecto con

el fin de exponer la importancia que tiene la caracterización de las gotas generadas en un
atomizador de tipo pressure-swirl.En primer lugar se establecerá dentro de las turbinas de
gas que importancia tiene la inyección de combustible y la atomización del combustible. Pos-
teriormente se va a pasar a la explicación de la combustión difusión, que es la que se produce
en general en este tipo de turbinas, y la importancia que tiene la correcta atomización en la
misma. A continuación se procederá a explicar en que consiste la inyección de combustible,
en especial las motores de turbina de gas para aviación.

2.1. Turbinas de gas

Las Turbinas de gas, en su expresión esquemática más simples,como se puede ver en la
figura 2.1 se trata de un motor de combustión interna de ciclo abierto, compuesto por un
compresor, una cámara de combustión y una turbina de gas propiamente dicha.

Figura 2.1: Esquema TG [?]

El funcionamiento de la turbina de gas presentada en el esquema es el siguiente:

1. Entra gas atmosférico por la unidad compresor que es accionada mediante el grupo
turbina de gas
2. Los gases comprimidos del compresor introducen en la cámara de combustión donde se
inyecta y quema el combustible, aumentado la temperatura sin prácticamente variación
de presión.
MEMORIA

3. Tras esto la mezcla de gases atmosféricos junto a los productos de la combustión a alta
temperatura y presión pasan al grupo turbina, donde cede potencia a la misma para
accionar el grupo compresor y lo sobrante al eje, en caso de tener uno.

4. por ultimo los gases de escape se expulsan fuera del grupo turbina [?].

2.1.1. Cámara de combustión

En las turbinas de gas el principal aporte de energía se da en la cámara de combus-
tión,observable en la figura 2.2, donde la combustión se representa entre los puntos 2-3.
Debido a lo comentado, los procesos y los parámetros que afectan a la combustión son de
muy alta importancia para el óptimo funcionamiento de la turbina de gas en su conjunto.
También se puede obtener la importancia de la combustión mediante la observación de como

Figura 2.2: ciclo real sencillo Bryton (no regenerativo): (a) en el palno pv; (b) en el plano Ts
ó hs.[?]

se obtiene el rendimiento térmico en estas máquinas de manera teórica, siendo el rendimiento

térmico ηt y definiendo:

Qa = Calor absorbido en la combustión.

Wt = trabajo que absorbido por la turbina.

Wc = Trabajo realizado por el compresor.

ηt = (Wt − Wc )/Qa (2.1)

La manera en la que se atomiza el combustible no únicamente afecta a la calidad de la
combustión si no que obliga a modificar los parámetros de la cámara de combustión o bien
esta limita de alguna manera como se tiene que dar al a atomización.

Página 15 de 57
MEMORIA

2.2. Inyectores
Los inyectores son parte integral de la cámara de combustión dentro de casi cualquier motor
térmico y por lo tanto componente obligatorio de las turbinas de gas. También tienen otras
aplicaciones a parte de la combustión, como pueden ser aerosoles o algunas aplicaciones mé-
dicas.

En la combustión la función de los mismos es la de introducir el combustible en la cámara

de combustión, además de atomizarlo para facilitar la evaporización combustible y para que
posteriormente se de la mezcla aire-combustible. El más común utilizado para turbinas de
gas es el inyector centrífugo o de swirl [?].Están formados por un orificio de salida, el cual
esta precedido por una cámara de swirl en la cual se introduce el combustible por conduc-
tos o aberturas tangenciales. Gracias a los conductos laterales se hace girar el combustible
creando un vórtice en el interior de la cámara de swirl pegando el combustible a las paredes
de la cámara. Por el orificio de salida el fluido sale en forma de lámina anular, el ángulo de
salida varía en función del ámbito en el que se vaya a utilizar el inyector[?]. El atomizador
que se va a estudiar en este proyecto es del tipo de los mencionados anteriormente, se hará
una descripción más detallada del atomizador estudiado en el capitulo 5: Descripción del
trabajo.

Figura 2.3: Esquema de un atomizador pressure-swirl[?].

2.3. Combustión
2.3.1. Introducción
Se entiende como combustión una reacción química exotérmica de oxidación que se da entre
algunos componentes de los combustibles, con el oxígeno de la atmósfera, aunque el combus-
tible pueda ser sólido, líquido o gaseoso los componentes del mismo que reaccionan con el
comburente tienen que estar en estado gaseoso. En las TG se usan en general combustibles de
origen orgánico, que están compuestos en su mayoría por carbono C, hidrógeno H, oxigeno
O y azufre S.[?].

Página 16 de 57
MEMORIA

Figura 2.4: Flujo de atomizador pressure-swirl totalmente desarrollado [?].

Durante la combustión se dan una serie de procesos, en primer lugar se da una disocia-
ción de los componentes del combustible, posteriormente estos reaccionan con el comburente
generando los que se conoce como radicales, tras esto comienza la reacción de oxidación pro-
piamente dicha y es donde se produce la principal liberación de energía, además de generarse
los productos finales de la reacción de combustión. [?]

Dentro de los tipos de combustión, en función de cuando se haya realizado el proceso de

combustión se pueden distinguir dos tipos de combustión:
Llama premezclada: Esta se dá cuando el combustible y el aire se introducen por el
mismo flujo, es decir como su propio nombre indica la mezcla combustible-aire se da
previa a la reacción de combustión.
Llama de difusión: Esta se da cuando el combustible se introduce de manera inde-
pendiente al aire, o directamente sobre un medio en el cual hay aire estático, estando
por lo tanto combustible y comburente separado, y dándose la reacción química de
combustión y la mezcla de combustible y comburente de manera simultanea.[?]
En el caso estudiado se da un llama difusiva, ya que el atomizador introduce unicamente
el combustible en la cámara de combustión de la turbina de gas, y es en la misma cámara
de combustible donde posteriormente se mezcla el combustión con el aire que viene del
compresor y se realizará la combustión, como se explica en el punto 2.1.

2.3.2. Combustión de Spray

Muchas e importantes aplicaciones, entre las que se encuentran la cámara de combustión de
las turbinas de gas, conllevan la combustión de combustible sólido o líquido, usando gene-

Página 17 de 57
MEMORIA

ralmente aire como comburente. Entendiendo spray como cualquier sistema con un numero
de partículas tal que solo una descripción estadística del mismo se posible [?], por lo que
en el caso de la utilización de que se trate de un líquido, el spray, es el sistema que se
formaría tras la atomización de dicho fluido. A esta combustión sería mas correcto llamarla
combustión de dos fases, ya que una fase condensada, la líquida, se introduce en una fase gas.

Dentro de este tipo de combustión se puede diferencia la siguiente secuencia de proceses:[?]

1. Inyección y atomización del combustible líquido.

2. Mezcla de gotas de combustible con el gas oxidante.

3. Transmisión de calor a las gotas provocando su vaporización.

4. Mezcla vapor de combustible con el gas.

5. Ignición de la mezcla, o por acción de un actuador externo o por autoencendido de la

mezcla.

Como se puede observar la atomización del combustible es parte fundamental para este
tipo de combustión, ya que las características físicas de las gotas que se generen durante la
atomización afectarán al resto del proceso y por lo tanto a la calidad final de la combustión.

Página 18 de 57
Atomización

Durante todo el capitulo anterior se ha destacado la importancia de la atomización para el

proceso de combustión dentro de la cámara de combustión de una turbina de gas. En este
capitulo por lo tanto incide más en la importancia del mismo y se explica como tal de que se
trata el proceso de atomización, los parámetros que afectan al proceso y las características
de los sprays generados por la atomización.

3.1. Definición de Atomización

Se entiende como atomización la transformación de un volumen líquido en un spray y otras
dispersiones físicas de pequeñas partículas en una atmósfera gaseosa[?]. Como se a comen-
tado en el capitulo anterior este fenómeno es de gran interés industrial, no unicamente en el
campo de la combustión como previamente se ha mencionado, si no que tiene gran utilidad
en otros ámbitos, ya sea en la industria química, por la necesidad de que reactivos estén
en condiciones especificas para su correcta reacción, la textil o en la fabricación de piezas,
pinturas y capas protectoras para ciertos materiales necesitan de la correcta atomización de
sus componentes para conseguir los acabados uniformes deseados.

En el proceso de atomización un chorro o una lámina líquidos son desintegrados, la de-

sintegración se da bien por la exposición a un flujo de aire o gas a alta velocidad o bien
por el aporte de energía mecánica mediante un aparato de rotación o de vibración [?]. La
naturaleza aleatoria del proceso junto a las múltiples maneras en las que se puede dar el
proceso de atomización complica en gran medida el estudio de este proceso, de ahí que no
aun no haya ninguna teoría que unifique el proceso de atomización, debido a esto el uso de
CFD ha ido en aumento en el estudio de los factores que afectan a la atomización, gracias a
esto se pueden obtener datos, no solo del resultado final de la atomización si no que también
se pueden obtener datos de como se desarrollan las diferentes fases del proceso.

3.2. Factores que influencian en la atomización

Cualquier atomizador, pese a la complejidad del mismo y de su tipo, depende de las ca-
racterísticas físicas del fluido atomizado, las características físicas del medio en el que se
va a atomizar y de la geometría y dimensiones del atomizador empleado. En los siguientes
MEMORIA

apartados se va a explicar como afectan lo factores previamente nombrados, a excepción de

la geometría, está se explicará en el apartado 3.3.1., ya que es donde se indaga más en las
características del tipo de atomizador que se va a estudiar.

3.2.1. Propiedades del liquido

El flujo y el spray tienen sus características condicionadas en gran parte por la densidad,
la viscosidad y la tensión superficial del líquido atomizado. En la practica, desligar la den-
sidad de cualquiera de las otras dos variable es prácticamente imposible, debido a esto su
relación a de ser tomada de manera cautelosa, además la poca variación de esta cualidad
entre los líquidos disminuye en gran medida su importancia. Lo antes mencionado junto a la
baga bibliografía que relaciona la densidad con el tamaño medio de gota indica que la poco
importancia que tiene[?].

El aumento del área de líquido a lo largo de la atomización esta considerado un indica-

dor de la calidad de la misma . Al salir el flujo de la aguja sale como un cilindro o como una
lámina, pero tras la atomización se desintegrarán hasta convertirse en gotas, estas tienen en
total mayor área que la lámina o el cilindro inicial, se usa un factor de multiplicación como
el indicar del nivel de atomización que se ha explicado antes. La tensión superficial del
líquido es importante para lo antes mencionado, ya que representa las fuerzas que resisten
la formación de nuevas superficies, es decir evita la partición del fluido en partes de menor
volumen. //// La viscosidad del fluido atomizado basas su importancia en que afecta tanto
a la creación de nuevas superficies, como lo hace la tensión superficial, como al flujo en el
atomizador y al patrón generado en el spray.La viscosidad esta estrechamente relzacionada
con el número de Reynolds, por lo que la variación de esta, desemboca en la variación del
tamaño de gota y del tiempo que tarda en aparecer las gotas. Pese a que el efecto de la
viscosidad dentro del flujo del atomizador es importante, este es complejo y complicado de
estudiar ya que el efecto varía en función del tipo de atomizador, en el caso de los tipo
pressure-swirl el incremento de viscosidad estrecha el ángulo de ángulo de atomización ge-
neralmente. En cualquier caso, a menor viscosidad la calidad de la atomización aumentará[?].

Para el análisis de cualquier caso de atomización es importante no solo las características

del fluido anteriormente descritas, ya que otras, como la temperatura a la que se encuentra,
afectan a las anteriormente mencionadas.

3.2.2. Condiciones Ambiente

El aire o gas ambiente en el que se atomiza el fluido varia en presión y temperatura, esta
variación es bastante grande en sistemas de combustión. Para los atomizadores tipo pressure-
swirl la densidad del fluido ambiente es bastante importante ya que la con el aumento de
la misma el ángulo de dispersión disminuye hasta un ángulo mínimo a partir del que no
disminuye más.

La presión ambiente también afecta a la creación de gotas, ya que si se aumenta desde el

Página 20 de 57
MEMORIA

valor de presión atmosférica el tamaño medio de gotas va aumentando hasta llegar un punto
de presión máxima donde el tamaño medio de gota comienza a disminuir lentamente[?].

3.2.3. Números adimesnionales

Los números adimensionales, son parámetros que permiten ralzcionar las características fí-
sicas de los fluidos atomizados con las condiciones de trabajo a las que opera el atomizador.
Estos son bastante imporatnes ya que nos permiten un análisis y control de la atomización
de manera global y sin depender de unidades. Entre estos parámetros destacan el número
de Reynolds, el número de Weber y el coeficiente de densidades combustible/aire.

El número de Reynolds, nos indica el nivel de trubulencia del fluiod, indicando por una
parte si se trata de turbulento o laaminar, además del nivel de turbulencia del fluido, este
este estrechamente relacionado con la viscosidad, tanto cinemática como dinámica, como se
ha comentado en el apartado anterior; pero además también se realciona con la velocidad del
fluido y la longitud carácterisitcas del atomizador. Un aumento de este provoca un aumento
de la turbulencia del fluido y por lo tanto un mejor atomización.
ρDv Dν
Re = =
µ µ
Re → N úmerodeReynolds
ρ → Densidaddelf luido
(3.1)
D → diámetrocaracterísitico
v → velocidaddelf luido
µ → V iscosidaddinámica
ν → viscosidadcinemática
El número de Weber, nos relaciona la fuerzas inerciales que actuan sobre fluido con la
tensión superfical del mismo, teniendo en cuenta esto se puede indicar que el número de
weber relaciona la velocidad a la que se introduce el fluido con la tensión superficial del
mismo, siendo esto muy importante ya que la varaición del mismo con un caudal constante
provaocará que la varía el tamaño de gota míninmo así como la cantidad de gotas. En la
ecuación 3.2 se obsrva la realción que existe entre el número de Weber y distintos parámetros,
observandose que al reducirlo aumentara la tensión superficial y por lo tanto le costará más
generar gotas.
ρv 2 l
We =
σ
(3.2)
l → longitudcaracterística
σ → tensiónsuperf icialdelf luido
El cociente de densidades combustible/aire, como su mismo nombre indica realciona las
densidades de los fluidos que interactuan durante la atomización, esto permite predecir como
va a afectar la densidad del fluido atomizado incluso aunque la densidad del medio sobre el
que se efectua la atomización no sea la misma siempre.

Página 21 de 57
MEMORIA

3.3. Clasificación atomizadores

Para la generación de un spray, todo lo que se necesita es una velocidad relativa bastante
alta, entre el fluido atomizado y el gas ambiente que sobre el que se atomiza. Teniendo esto
en cuenta se puede clasificar el tipo de atomización como aquellos que inyectan combustible
a alta velocidad sobre aire o gas a poca o nula velocidad, y los que utilizan el aire para
atomizar el combustible.

Dentro de los que inyectan el combustible a alta velocidad podemos distinguir:

Figura 3.1: Ejemplos de atomizadores de rotación a)guías rectas, b)guías curvas[?].

Atomizadores de presión: El líquido se introduce por una pequeña apertura me-

diante una alta presión, una vez salido del agujero la energía que se le aporta al flujo
mediante la presión de inyección se transforma en energía cinética. Dentro de estos
encontramos los tipo pressure-swirl, explicados más a detalle en el apartado 3.3.1.
Atomizadores de rotación: En estos atomizadores se le introduce energía cinética
de mediante un aparto que haga rotar a alta velocidad, tras esto el líquido se descarga
por la periferia del atomizador, se ve un ejemplo de est tipo de atomizadores en la
figura 3.1.
Dentro de los que se introduce el fluido a baja velocidad y se les introduce a un flujo de aire
a alta velocidad se pueden distinguir:
Atomizadores asistidos por aire:El atomizador está expuesto a una corriente de
aire a alta velocidad. Dentro se pueden distinguir, la mezcla interna del combustible y
el flujo de aire, de tal manera que se mezcla el combustible y el aire a alta velocidad
antes de de salir del orificio de la descarga, de mezcla externa, se expulsa el combustible
y el aire a alta velocidad los arrastra. Estos son bastante ineficientes desde el punto de
vista energética, no obstante genera gotas de menor tamaño que lo de presión simple.

Página 22 de 57
MEMORIA

Atomizadores Airblast: Estos se distinguen de los anteriores en que los asistidos

por aire usan un pequeña cantidad de aire a muy alta velocidad y los segundos usan
altas cantidades de aire a bajas velocidades, en comparación a los anteriores.

Por último quedan algunos atomizadores que quedan fuera de los anteriormente menciona-
dos, estos son poco utilizados ya que sirven para aplicaciones muy concretas, entro estos se
encuentran los electro-estáticos, los ultrasónicos, los sónnicos (o de silbido) o los vibratorios,
entre otros.

Figura 3.2: Ejemplos de atomizadores según su tipo a)Atom. de presión, b)Atom. de rotación
c) Atom. de doble fluido[?].

3.3.1. Atomizador tipo Pressure-Swirl

Dentro de los atomizadores de presión se pueden encontrar varios tipos entre los que se va a
destacar los atomizadores conocidos como pressure-swirl o simplex. Estos se componen por
una cámara de swirl o de vórtice seguida por una un orificio de salida por donde sale el fluido
en forma de cono. El ángulo que tiene el cono de flujo que sale por el orificio de salida puede
variar en gran medida, de 30o y 180o . La atomización en estos atomizadores es generalmente
buena, pero se obtiene la atomización de gotas de menor tamaño con grandes presiones de
entrega y con ángulos de gran apertura. Se pueden diferenciar diversos tipos básicos según
la forma del fluido que generan a la salida del atomizador.

Solid cone Spray : Las gotas generadas a la salida del atomizador se distribuyen
uniformemente a lo largo de todo el volumen del cono generado.

Hollow cone Spray: En estos las gotas en vez de distribuirse a lo largo de todo el
volumen, se concentran en el el exterior del mismo, como se muestra en la figura 3.3.

Los atomizadores de tipo hollow cone generan una atomización mejor que los de tipo solid
cone, y su distribución radial, con el hueco dentro, es bastante favorable en la numerosas
aplicaciones, sobretodo en los sistemas de combustión.

Página 23 de 57
MEMORIA

Figura 3.3: Tipos de atomización en atopmizadores tipo pressure-swirl[?].

Los tipos de más básicos de hollow cone tienen unos puertos de entrada tangencial a la
cámara de swirl, la entrada de combustible por estos puertos le infiere al combustible una
alta velocidad angular al combustible al entrar a la cámara de swirl, creando ya dentro de
la cámara el cono vacío que los caracteriza. A la salida del orificio el líquido afectado por las
fuerzas axiales y radiales genera una lámina de combustible, cuyo ángulo de dispersión está
definido por los componentes axiales y tangenciales de la velocidad que gobiernas la lámina
del combustible.

El desarrollo de esta lámina y del por lo tanto del flujo de combustible, se desarrolla
mediante un serie de pasos cuando la presión aumenta hasta el punto de funcionamiento
desde cero.

1. Comienza a generar pequeñas gotas.

2. Se genera un haz distorsionado

3. Se genera un cono a la salida del orificio pero se debido a las tensiones superficiales
fuerzan al cono a cerrarse formando una burbuja cerrada.

4. Se abre la burbuja en forma de flor con el borde irregular, desde el borde empieza a
desintegrarse la lámina en gotas relativamente grandes.

5. La curva en forma de flor se endereza, generando la lámina cónica característica de este

tipo de atomizador. Tras esto la se reduce el grosor de la lámina y consecuentemente
comienza a desintegrase en ligamentos y posteriormente en gotas.

En la figura 3.5 se pude observar las distintas fases descritas.

Página 24 de 57
MEMORIA

Figura 3.4: Esquema de la forma de un atomizador tipo pressure-swirl simplex[?].

Figura 3.5: Fotografías del desarrollo del flujo de un atomizador swirl al ir aumentando la
presión desde 0 [?]

3.4. Procesos Básicos de la atomización

El proceso de atomización se puede definir como, convertir un volumen de líquido que se
desintegra en pequeñas gotas.La manera en que se da esto es mediante la oposición, hasta

Página 25 de 57
MEMORIA

superarla, de la acción unificadora de la tensión superficial mediante la acción de agentes

internos y externos. En caso de que no existan las acciones que se opongan a la tensión su-
perficial, la tensión superficial tiende darle forma de esfera al volumen de agua, la viscosidad
del líquido estabiliza la masa de líquido ejerciendo oposición a la variación de geometría.
Entre las acciones que crean las inestabilidades en las gota las fuerzas aerodinámicas sobre
la superficie de líquido son principales para la formación de las gotas.

Las gotas creadas en primera instancia, las más grandes que se generan, tienen mucha inesta-
bilidad por lo que de manera consecuente se volverán ha desintegrar en gotas más pequeñas.
Debido a lo previamente explicado se pude ver que el tamaño final de las gotas que se generan
depende tanto de la primera atomización y como de las atomizaciones que siguen a esta.

3.4.1. Formación estática de gotas

Previamente a la explicación de los procesos mas complejos se explicará la formación gotas
en su caso más simple. Este caso es el cuasi-estático de una gota colgando. Esto se vería
claramente en cualquier aparato que gotee o que tenga una pequeña lamina de agua que de
hacia abajo, en esta la gravedad realiza la acción destabilizante sobre el volumen de agua, en
el momento que esta acción supera a la tensión superficial una pequeña parte del volumen se
desprende y genera una gota. Este es el método más básico de atomización y el mas común
en la naturaleza y aunque no sea útil en aplicaciones industriales, ayuda al entendimiento
de como funciona la desintegración del fluido en gotas.

3.4.2. Rotura de gotas

En el caso más común en la industria, que es el que se da en el caso de estudio en este trabajo,
se da la interacción del fluido con el aire atomizándolo. Cuando la atomización se da por la
interacción del líquido con el aire hay muchos mecanismos involucrados, para conocer todas
las maneras en que las gotas se pueden formar al entrar en contacto con el aire se haría
falta saber a la perfección como actúan las fuerzas aerodinámicas sobre las gotas, esto es
extremadamente complejo ya que con la acción misma de las fuerzas aerodinámicas deforma
las gotas y por lo tanto cambian las acciones que se generan sobre sus superficies. Estas
deformaciones también modifican las fuerzas internas, lo que aumenta aun más lo complejo
de una solución matemática rigurosa.

Para comprender bien como funciona el mecanismo de generación de gotas hay que defi-
nir lo siguiente; la presión interna sobre cualquier punto de la superficie, pI , es una constante
y es sufciente para balancear la presión ejercida por las fuerzas aerodinámicas, pA y las
presión generada por la tensión superficial, pσ , de tal manera que

pI = cte = pA + pσ (3.3)

Viendo esto, siempre que el aumento o la disminución de las presiones, aerodinámicas o de las
tensiones superficiales, sea compensado por la otra y mantengan constante la presión interna,

Página 26 de 57
MEMORIA

no se generan gotas. No obstante, cuando pA es grande, en comparación a pσ , un cambio

en pA no puede ser compensado por pσ y por lo tanto pI no se puede mantener constante.
Esto produce que la gota comience a deformarse hasta que se desintegre en gotas de menor
tamaño. La presiones debida a las tensión superficial se pude definir en gotas esféricas como:
4σ
pσ = (3.4)
D
Es por lo tanto observable que cada vez que se desintegra una gota pσ aumenta para las
nuevas gotas, esté aumento de presión provoca que una vez que pσ sea los suficientemente
grande pueda asumir las variaciones de de pA evitando por lo tanto la aparición de nuevas
gotas.

Estas observaciones ayudan a definir lo que se conoce como tamaño crítico de gota, siendo
este el menor tamaño de gota que puede alcanzar el fluido en unas condiciones. Esto establece
además que el tiempo de formación de gotas va aumentando a medida que va reduciendose
el tamaño de las gotas. Al tiempo de formación gotas también afecta la viscosidad, ya que
como se explicado previamente esta se opone a la deformación de la gota. Además respecto
al tiempo de creación de gotas, aunque pA inicialmente sea suficiente para desintegrar la
gotas, si pA disminuye antes de que se forme la gota, si esta disminuye lo suficiente la rotura
no se dará.

3.4.2.1. Rotura de gotas en aire en movimiento

El estudio de este caso se ha realizado en muchos casos, y es lo mismo que si en vez de
mover una gran flujo de aire sobre un fluido, se mueve a alta velocidad un fluido a través de
aire. Multitud fotografías de alta velocidad [?] indican que las maneras en que se atomiza
un volumen de líquido depende del flujo de aire relativo que es genera a su alrededor. No
obstante se han definido 3 tipos básico[?]:

1. Lenticular: la gota se aplana por los polos de tal manera que adopta la forma de un
elipsoide. Se cree que esta elipsoide se acaba convirtiendo en un toroide estrechándose
hasta desintegrarse.

2. Cigar-shaped: la gota se va alargando de forma cilíndrica hasta que se rompe los

ligamentos generando nuevas gotas.

3. Amorfa: En esta las acciones internas y externas le generan protuberancias que poco
a poco irán separándose del la gota original.

Resumiendo la desintegración de las gotas sometidas a un chorro de aire, o que fluyen a

alta velocidad sobre aire, depende de la presión dinámica, la tensión superficial y la fuerzas
viscosas.

Página 27 de 57
MEMORIA

Figura 3.6: Ilustración tipos de gotas[?]

3.4.3. Desintegración de lámina líquida

La manera en la que se desintegra el fluido en un atomizador, depende en gran parte de la
manara en la que sale el flujo por el orificio de salida, en jets o en lámina, como el atomizador
estudiado en el trabajo se genera un flujo de tipo lamina se va explicar.

El desarrollo de la lámina esta principalmente controlada por la velocidad a la salida del

orificio y de las propiedades físicas del fluido. Tendiendo lo anterior en cuenta se puede en-
tender que para permitir que la lamina se expanda, es necesario una velocidad mínima para
vencer a la tensión superficial que la contrae, esta velocidad en el caso de los pressure-swirl
viene dada de por la presión a la que se introduce el combustible. Aumentar la velocidad
aumenta la longitud de la lámina, que acaba en un borde donde las fuerzas inerciales y las
tensiones superficiales se equilibran.

Se pueden distinguir tres modos o maneras en las que se pueden desintegrar la láminas
llamadas

Rim

Wave

Perforated-sheet

Modo rim, en este modo la fuerzas provocan que en el extremo libre de la lamina se genere
un borde grueso, este borde se desintegra igual que un jet, en este las fuerzas disruptivas y
cohesión interactuan entre si de tal manera que provocan inestabilidades en el flujo que se
manifiestan como oscilaciones, provocando que finalmente se desintegre en gotas.

Modo perfirated-sheet, en este tipo de desintegración aparecen agujeros empujando el

fluido que se encontraba originalmente dentro de estos al borde de los mismos, estos aguje-
ros aumentan de tamaño y cantidad hasta que los bordes de los agujeros continuos colapsan
entre si formando ligamentos que posteriormente producen las gotas al desintegrarse.

Modo wave, se generan una onda de movimiento en la lamina tal y que en los puntos

Página 28 de 57
MEMORIA

de la lamina que se encuentran a longitud de onda o mitad de longitud de onda, partes de

las láminas son despegadas sin llegar al borde de la lámina, estos fragmentos arrancados se
desintegran por la presión aerodinámica o la turbulencia interna de la gota.

Generalmente las láminas que se generan en los atomizadores pueden presentar los tres
modos de desintegración, muchas veces se presentan dos modos de manera simultanea, in-
fluyendo cada uno de los modos al tamaño medio de gota y la distribución de las mismas en
función de la importancia que tenga cada uno de los modos en el proceso de atomización.

Figura 3.7: Fotografía que ilustra el modo wave principalmente [?]

3.5. Características de los sprays

En gran parte de las aplicaciones las funciones que desempeña el atomizadores no son solo
desintegrar el fluido en gotas de pequeño tamaño si no descargarr las gotas generadas sobre
el aire cumpliendo ciertas condiciones de simetría y uniformidad.

Como se ha explicado previamente en el apartado 3.3.1. en los atomizadores tipo pressure-

swirl pueden generar dos tipos de spray, el primero y menos común es el sólido y el segundo
y más usado el hueco. En ambos caso las gotas llevan la misma dirección que lleva el flujo
antes de desintegrarse hasta que por la resistencia del aire, las gotas que están en la parte
más externa del spray van perdiendo su momento tal que se detienen formando una neblina

Página 29 de 57
MEMORIA

de pequeñas gotas, tras esto estas gotas detenidas se mueven según las gravedad.

Los sprays de los atomizadores de tipo pressure-swirl son bastante dependientes de la geo-
metría del atomizador en lo que a la distribución de gotas se refiere, ya que el ángulo que
forma la lámina inicialmente depende en gran parte de la geometría interna de la aguja, sus
condiciones de trabajo y las propiedades del líquido, por lo que al aumentar el ángulo de
salida aumentará, generalmente, la exposición al aire que ha sido descrita previamente.

Algunas de las principales propiedades que presentan los sprays y que pueden ser de in-
terés son las mostradas a continuación.

Dispersión
A la dispersión es posible dar valor cuantitativo para un instante conocido, siempre que se
conozca el volumen de líquido sea conocido. Se puede dar la siguiente definición para el grado
de dispersión, relación de volumen del spray con el volumen de líquido contenido en el[?], se
puede definir pues como cuan disgregado está el líquido del spray.

Una buena dispersión del spray ayuda a que se mezclen el aire y el líquido y por lo tan-
to ayuda a que la evaporación sea más rápida y efectiva.En el caso de los atomizadores
pressure-swirl está controlada por otras características del spray, algunas de estas son, el
ángulo de cono, el tamaño medio de gotas y la distribución de las gotas.

Penetración
Se define penetración como la distancia máxima que alcanza el spray cuando se inyecta en
aire estancado, aunque también se puede definir como distancia máxima que alcanza al in-
yectarse en aire con flujo cruzado, esto varia en función del medio en el que se inyecte el
combustible[?].

La magnitud de esta característica está principalmente controlada por dos fuerzas, una
perteneciente al líquido y otra al aire que lo rodea, siendo la primera la energía cinética
a la salida del atomizador que lleva el líquido y la resistencia aerodinámica que presente
el aire circundante al líquido. Pese a que a la salida del atomizador la velocidad se muy
alta, el aumento en el superficie del spray, en adición a la disminución progresiva de la ener-
gía cinética, provoca que la trayectoria de las gotas en última instancia este gobernada por
la gravedad o por la fuerza que ejerce el movimiento de los gases, si es suficientemente grande.

Como se pude intuir, a mayor es el ángulo de salida en un atomizador menor será la pene-
tración del spray ya que habrá mayor superficie de contacto entre el aire y el líquido y por
lo tanto perderá más energía cinética.

Página 30 de 57
MEMORIA

Figura 3.8: Ejemplo definición ángulo de cono en inyector diésel [?].

Ángulo de cono
Debido a que la lámina que se crea en los atomizadores de tipo pressure-swirl se curva la
definición y la medición del ángulo de cono es bastante complicada, generalmente esta esta
definida por dos lineas recta que definen el contorno del spray desde la salida del atomizador
hasta una distancia especifica del atomizador. Pese a lo complicado de la definición del ángulo
de cono, este es muy importante para conocer el comportamiento exterior del spray, debido
a esto se han desarrollado algunos protocolo estándares para su medición, esto supone que
la aportar cualquier medición hecha del ángulo de cono hay que aportar una descripción
detallada de cual ha sido el método que se ha utilizado.

Patrón
La simetría que tiene el patrón que presenta el spray es una característica es bastante práctica
en diversas aplicaciones, ya que si el spray no es simétrico en aplicaciones relacionadas a la
combustión no permite un correcta mezcla aire y gas.Si la falta de simetría no es muy
grande no se puede detectar fácilmente a simple vista por lo que es necesario el desarrollo de
métodos específicos de control para cada aplicación. El patrón de se caracterizar de manera
radial y circunferencial, para cada una de las maneras de caracterizar el patrón hay métodos
experimentales distintas.

Distribución de Gotas
Una manera de englobar algunas de las caracterísiitcas que se han mencionado en los apar-
tados anteriores es lo que se conoce como distribución de las gotas. Para el estudio de las
misma se busca obtener lo que se conoce como función densidad probabilística, o probability
density fuction,PDF, como se conocen en inglés. En las PDFs se presenta la probabilidad de
que haya una gota en a una distancia específica del orificio de salida, también se puede hacer
para un rango de tamaños de gotas así como presentar varios tamaños en una sola gráfica

Página 31 de 57
MEMORIA

Figura 3.9: Ejemplo PDF de que muestra la probabilidad de encontrar gotas, en función de
su diámetro, a diferentes distancia del orificio de salida. [?] [?].

comparativa. Estas gráficas son bastante interesante y aportan bastante información ya que
pueden mostrar por una parte cual sería la penetración del flujo, así como mostrar donde se
encuentran la mayoría de gotas. La manera de representar estas funciones es bastante varia-
da y varía en gran medida en función de lo que se quiera mostrar. Se muestra un ejemplo de
estas funciones en la figura 3.9. En la realización del trabajo se utilizan estas PDFs para la
caracterización de la atomización en el atomizador estudiado, así como para comparar como
afectan los distintos parámetros estudiados en el mismo.

Página 32 de 57
Herramientas utilizadas. CFD y soft-
ware Basilisk

En la realización de este proyecto se han usado diversas herramientas informáticas, la utilidad

de estas ha sido realizar el pre-proceso de los datos, el cálculo de las soluciones y el post-
proceso. Todas los softwares que se han usado en este proyecto son de uso de libre. Para el
procesado de datos, en pre-proceso y post-proceso, se ha usado Python, que es un lenguaje
libre de código que permite la lectura y escritura de archivos, además de que se permite
hacer cálculos complejos[?]. Para el cálculo del flujo que se necesita para caracterizar las
gotas que se generan se utiliza un software conocido como basilisk que pertenece al grupo de
herramientas informáticas de CFD (Computational Fluid Dynamics) o dinámica de fluidos
computacional. Los cálculos realizados por CFD son los más característicos para realizar
las simulaciones de fluidos, de ahí que se vaya a explorar más en profundo en las próximas
secciones.

4.1. CFD
4.1.1. ¿Que es el CFD?
CFD o dinámica de fluidos computacional es el análisis de sistemas que envuelven el flujo
de fluidos, la transmisión de calor y los fenómenos asociados, entre lo que se destacan las
reacciones químicas, mediante el uso de simulaciones basadas en la computación.

En los últimos años el desarrollo del CFD ha sido para equipararlo en capacidad de cálculo
a otras herramientas de cálculo, tales como las de resistencia de estructuras, ya que debido a
la complejidad de los cálculos del flujo ha provocado que el desarrollo de estas herramientas
se vean rezagadas en comparación a otras.

Actualmente el CFD se emplea en múltiples ámbitos de utilización, algunos de estos ám-

bitos son los siguientes:

Cálculo aerodinámico para aeronaves y otros vehículos.

Hidrodinámica para barcos y otras naves acuáticas.

MEMORIA

Análisis de combustión y proceso de que se dan dentro de motor de combustión interna.

Análisis de flujo interna dentro de turbo máquinas.

Investigar procesos en reacciones químicas y como se desarrollan los productos de las

reacciones.

Estudio de la distribución de contaminantes y partículas.

Teniendo en cuenta que hay una gran variedad de programas de código libre de CFD que
permiten realizar todo tipo de cálculos, el coste de realizar estos cálculos dependen por una
parte del tiempo invertido por parte de un operario para programar la simulación así como el
tiempo invertido en el cálculo, ya que en muchos casos es necesario el uso de computadoras
de gran potencia de cálculo.[?].

No obstante a lo previamente comentado el uso del CFD presenta ventajas respecto al uso
de métodos únicamente experimentales, entre las que se destacan:

Reducción en gran medida de tiempos de diseño y los costes de nuevos diseños.

Posibilidad de estudio de sistemas donde su control en experimentos sea dificultoso o

imposible.

Posibilidad de estudiar un fenómeno en condiciones aleatorios.

Cálculo virtualmente sin límite en detalle del resultado.

En oposición a estas ventajas el uso de CFD conlleva algunas desventajas:

Necesidad de gran inversión en hardware y software.

Necesidad de personal especializado en tanto en el ámbito de estudio, así como en

programación.

4.1.2. Funcionamiento
Los códigos de CFD se estructura en función de algoritmos numéricos. Todos los códigos que
se pueden encontrar tienen tres elementos principales que facilitan la organización[?]:

pre-procesamiento

solver

post-procesamiento

Pre-procesamiento: Se trata de la introducción de los datos de los problemas de flujo

propiamente dicho, se establecen datos desde los que parte el problema y se transforman
para que el solver los procese. Dentro de este se establece la geometría, llamada dominio,
se genera la malla, que se trata de la división del dominio en pequeños dominios que no se

Página 34 de 57
MEMORIA

solapen creando lo que se conoce como malla, se disponen las propiedades físicas y químicas
que se quieren simular, las propiedades de los fluidos que van a ser estudiados y las
condiciones de contorno, aquellas condiciones que limitan el dominio.

Solver: Se refiere al algoritmo numérico que calcula la solución al problema introducido

en el pre-proceso, existen diferentes técnicas de para hacer los cálculos, pero en general
se pueden distinguir tres pasos, integración de las ecuaciones al rededor de todo el volumen
de control, discretización, solución de las ecuaciones que se generan durante la discretización.

Post-proceso: Se trata del procesado de las soluciones generadas por el solver para que
sean accesibles de manera cómoda por el usuario, ya sea para analizarlas o utilizarlos a futu-
ro por otros programas o visualizarlos en softwares de visualización. Dentro de las maneras
que se visualizan se incluyen:
Geometría del dominio y la disposición de la malla, esto puede ser interesante antes de
lanzar una solución, para comprobar que es tal y como queremos.

Gráficos vectoriales de los parámetros de estudio.

Gráficos 2D y 3D de las superficies.

Salida de los gráficos organizados por colores.

4.1.3. Ecuaciones del fluido

Las ecuaciones que gobiernan el comportamiento de los fluidos representan las leyes físicas
de conservación de:
Masa de fluido.

Segunda ley de Newton.

Primera ley de la termodinámica.

El fluido se considera de un continuo y para escalas macroscópicas, igual o mayores a 1 µm,
no se tiene en cuenta la estructura de la materia y los movimientos molecular. Debido a
esto se describe el comportamiento del fluido en términos macroscópicos. Se considerara una
partícula o punto del fluido como el elemento más pequeño cuyas propiedades macroscópicas
no se ven afectadas por moléculas individuales.

La manera en que se define una partícula es considerar un volumen de lados, δx,δy y δz.
Como se puede ver en la figura 4.1 los lados generan unas caras la cuales se nombre N, S,
E, W, T y B, que se nombran así como norte, sur, este, oeste, arriba y abajo en inglés.
En la figura 4.1 se observa también las direcciones positivas y las coordenadas del centro
de la partícula (x,y,z). El estudio sistemático de los cambios en masa, momento y energía
del fluido que atraviesa sus límites, las caras del volumen, proporcionan las ecuaciones de
conservación del fluido.

Página 35 de 57
MEMORIA

Figura 4.1: Volumen mínimo, partícula de fluido [?].

Todas las propiedades de los fluidos son funciones del espacio y el tiempo, por lo que
ρ(x,y,z,t), p(x,y,z,t), T(x,y,z,t) y u(x,y,z,t), para los vectores de densidad, presión, temperatu-
ra y velocidad. A partir de ahora las anteriores magnitudes se escribirán sin las coordenadas
para facilitar la lectura.

Ecuación de conservación de la masa

La ecuación de la continuidad, o de conservación de la masa expresa la imposibilidad de
crear o destruir masa y por lo tanto su conservación a lo largo del tiempo.
δρ → − −
+ ▽(→
ρ ůu) = 0 (4.1)
δt
Si se busca el estudio de un volumen y no de una partícula en concreto se integra la ecuación
anterior obteniendo:
d Z Z
ρdV + ρů(→
−
us − →
−
ue )ů→
−
n dS = 0 (4.2)
dt Vc Sc

para cual:
d R
ρdV : Indica la variación de masa dentro del volumen de control.
dt Vc
→
− → − → −
Sc ρů( us − ue )ů n dS : Indica el flujo neto, lo que entra o sale, a través de la superficie
R

de control.

Ecuación de consevación de la cantidad de movimiento

Con esta ecuación se expresa la segunda ley de Newton. Esta ecuación indica que el sumatorio
de las fuerzas aplicadas sobre un fluido modifica su movimiento, acelerando, decelerando o

Página 36 de 57
MEMORIA

frenando dicho fluido.

∂(ρ→
−
u) →− − → dp →
− −
→
+ ▽(ρ→
u ů−
u ) = − + ▽τ + ρfm (4.3)
∂t dx
Los términos de la izquierda hacen referencia a la variación de movimiento mientras que
el termino de la izquierda hace referencia a las fuerzas que afectan al fluido. Integrando la
ecuación se obtiene:
d Z → Z Z Z
−
→
ρ−
u dV + ρ→
−
u (→
−
u −→
−
uc )ů→
−
n dS = τ ů→
−
n dS + ρfm dV (4.4)
dt Vc Sc Sc Vc

donde:
d R →
ρ− u dV : Indica el cambio de cantidad de movimiento dentro del volumen de
dt Vc
control.

ρ→
−
u (→−u −→ −
uc )ů→
−
R
Sc n dS: Indica el intercambio neto de materia con el exterior por las
superficies de control.

τ ů→
−
R
Sc n dS: Fuerzas superficiales
R −
→
Vc ρfm dV : fuerzas relacionadas a la masa.

4.1.3.1. Ecuación de conservación de la energía

Con esta ecuación se muestra que la cantidad de energía de un sistema aislado no puede
aumentar ni disminuir, puede eso sí transformarse en un tipo distinto de energía. Esta ecua-
ción tiene en cuenta la energía contenida en el sistema y la energía que entra y sale por las
superficies del volumen.
∂ρh → − →
− − → − →
− −
→−
+ ▽ů(ρh→
−
u )) = −p ▽ů→
u + ▽(σů→
−
u ) + ▽(k ▽ (T )) + ρfm →
u +Q (4.5)
∂t
El termino de la izquierda representan la energía interna y su variación, mientras que los de
la derecha muestran las diferentes manera de aportar manera al sistema o que el sistema ceda
energía, mostrando en los dos primero el trabajo ejercido sobre las superficies por unidad de
tiempo, el tercero el calor intercambiado a través de la superficie por conducción, el cuarto
representa el trabajo de las fuerzas másicas y el último el calor intercambiado por rediación
y reacciones químicas.

Estas tres ecuaciones son las mismas para todos los casos de estudio de fluidos. Las con-
diciones y los parámetros que se imponen al sistema cuando se realiza el estudio son las que
simplifican y particularizan las ecuaciones para cada caso. La complejidad de las ecuaciones
imposibilita la solución analítica, exceptuando casos muy básicos, esto provoca que se solu-
cionen mediante el uso de técnicas numéricas, aplicadas mediante computación.

Página 37 de 57
MEMORIA

Figura 4.2: Malla 2D o subdivisión de volumen de métodos finitos [?].

Pese a esto, las limitaciones que aun presentan los superordenadores hace que para aquellos
casos en que se hayan realizado simplificaciones y con estas se hayan obtenido soluciones
correctas, no se use el CFD. Las simulaciones de elemento de gran tamaño, por ejemplo
naves industriales enteras, así como de periodos de tiempo largos, no se pueden simular
en su totalidad, provocando que se simulen unicamente algunas partes, por separado, del
sistema[?].

4.1.4. Discretización
El uso de ordenadores para resolver las ecuaciones que de los fluidos obliga a que estas sean de
tipo algebraico y por lo tanto hay que convertir las ecuaciones diferenciales antes planteadas
(4.1, 4.2 y 4.5) en algebraicas. Esta conversión de ecuaciones diferenciales en algebraicas
se conoce como discretización numérica, las maneras de discretización más utilizadas son,
diferencias finitas, volúmenes finitos y elementos finitos. Para este proyecto se va explicar
brevemente el método de volúmenes finito ya que el software utilizado principalmente usa el
método de volúmenes finitos.

Método de volúmenes finitos

El método de volúmenes finito se vasa en convertir la ecuaciones diferenciales que gobiernan
el flujo mediante la integración de las mismas para un volumen de control, los cuales son
volúmenes finitos. Como se ve para simplificar, el volumen se subdivide en x e y para generar
volúmenes finitos y definidos. Como se pude observar en la figura 4.2 se puede observar una
celda o volumen con el punto central P, que es el que se quiere obtener como análogo a la
solución de la ecuación diferencial. El resto de puntos marcados, E,N,S y W, son los puntos
equivalentes a P en las celdas colindantes y las paredes o límites e,n,s y w son las superficie
que limitan el volumen finito donde se encuentra el punto P.

Página 38 de 57
MEMORIA

Par explicar como se realiza el cálculo de la derivada en P para una magnitud se va ex-
plicar para un caso con una sola dimensión para mayor facilidad:
∂ 2U
(4.6)
∂x2
Pues para obtener los valores del punto P se hace en función de los datos de las superficies
adyacentes:

Figura 4.3: Ejemplo discretización elemento unidimensional [?].

∂2U ∂2U
+
!
∂ 2U ∂x2e ∂x2w
= (4.7)
∂x2 P
xe + x w
Esta por su parte se obtienen por los datos de los puntos adyacentes a la superficie tal que:
∂ 2U UE − UP
= (4.8)
∂xe2 xE − xP
Gracias a estos términos se pueden implementar en cualquier celda para encontrar el análogo
numérico de las derivadas de para cada volumen, de esta manera e implementando unas
condiciones de partida se puede resolver el sistema de ecuaciones que se genera.

4.1.5. Solucionar de las ecuaciones discretizadas

Una vez discretizadas las ecuaciones, se consigue un sistema de ecuaciones algebraicas que
se pueden escribir de manera matricial esquemática como:
Aůx = b (4.9)
Siendo A los coeficientes y b los términos independientes.

La complejidad y el tamaño de estos sistemas de ecuaciones dependen de las dimensio-

nes del problema„ el numero de nodos que tenga la malla y las discretización hecha. Aunque
se pueda usar cualquier método de cálculo para resolver el sistema la potencia de calculo
informático disponible es lo que limita el método usado en gran medida.

Dentro de esos métodos se pueden englobar dos familias o grupos de métodos, siendo estos:

Página 39 de 57
MEMORIA

Métodos Directos.
Métodos indirectos o iterativos.
Los métodos directos son aquellos que basan la resolución de las ecuaciones mediante
métodos como el de la regla de la matriz inversa de Gramer y la eliminación Gaussian. En
estas se puede previamente saber el orden de cuantas operaciones se van a realizar, siendo
en un sistema de N ecuaciones con N incógnitas el orden de operaciones a realizar es N 3 .
Es necesario el almacenaje simulataneo de todos los coeficientes N 2 conjunto de ecuaciones
en la memoria principal [?].

Los métodos indirectos o iterativos como su propio nombre indica en se va en la aplica-

ción de algoritmos relativamente simples hasta alcanzar una convergencia tras un número de
repeticiones. Algunos de los ejemplos más utilizadas son el de Jacobi y el de Gauss-Seidel.
En estos métodos no se puede predecir el número de operaciones que se van a realizar y tam-
poco se puede asegura la convergencia de la solución si no se cumplen ciertas condiciones. Su
principal ventaja es que solo los coeficientes que no sean 0 son necesarios que se almacenen [?].

El método de volúmenes finitos genera un gran número de ceros normalmente. Debido a

lo antes mencionado y que los sistemas generados de problemas de CFD realistas generan un
gran número de ecuaciones, de 100.000 hasta 1 millón de ecuaciones, los métodos iterativos
son en su mayoría mas económicos que los directos. Planteando un sistema de ecuaciones
como este:
2x1 + x2 + x3 = 7
−x1 + 3X2 − x3 = 2 (4.10)
x1 − x2 + 2x3 = 5
Se aísla cada una de las en cada ecuación una de las incógnitas: n sistema de ecuaciones
como este:
x1 = (7 − x2 − x3 )/2
x2 = (2 + x1 + x3 )/3 (4.11)
x3 = (5 − x1 + x2 )/2
Se comienza dando valores supuestos a los diversas incógnitas y cambiando el resto hasta
que se llega a una solución. Tras esto se pasa al siguiente instante temporal.

4.1.6. Modelos de turbulencia

Todos los flujos que se encuentran en la práctica de la ingeniería, tanto los más simples,
como flujo en tuberías o chorros de dimensiones, como los mas complejos tridimensionales,
se vuelven inestables, y por lo tanto turbulentos al alcanzar cierto número de reynolds.

Es observable pues que al alcanzar alto número de Reynolds el flujo se vuelve turbulen-
to, siendo un flujo turbulento aquel que tiene un estado caótico y aleatorio del movimiento,

Página 40 de 57
MEMORIA

cambiando continuamente la velocidad y la presión a lo largo del tiempo.

La turbulencia que aparecen en los flujos, causa que se generen los remolinos de gran va-
riedad de rangos y escalas temporales, que interactúan de manera compleja y dinámica. La
importancia de la turbulencia en el estudio de flujos ha provocado el desarrollo de diferentes
métodos numéricos. Estos métodos se clasifican en tres categorías:

Modelos de turbulencia de Reynolds-avaraged Navier-Stokes (RANS): En

estos modelos se pone el foco en el promedio del flujo y por lo tanto el efecto en el
flujo promedio de la turbulencia. Previamente a la aplicación de los métodos de Navier-
Stokes las ecuaciones son promediadas en el tiempo. Las exigencias computacionales
necesarias para este tipo de modelos no son muy elevados por lo que han sido los más
usados en las últimas décadas.

Large eddy simlation (LES): Se trata de de una forma de cálculo intermedio de

turbulencia, estudiando la evolución de remolinos de gran tamaño. El método incluye
el filtrado de las ecuaciones inestables de Navier-Stokes previamente a la computación,
de tal manera que los remolinos mas pequeños no son calculados mientras que los
de mayor tamaño si que lo son. Estos modelos son interesantes, pese al aumento del
volumen de almacenaje y cálculo en comparación a los modelos de RANS, ya que sirve
como aproximación a problemas de geometría mas compleja.

Direct numerical simulation (DNS): Estas simulaciones, calculan el flujo medio a

la vez que calculan todas las variaciones turbulentas de velocidad. Las ecuaciones de
Navier-Stokes inestables son calculadas mediante una malla espacial suficientemente
fina. Este método solo se suele usar en investigación debido al alto coste en computación
que supone.

4.1.7. Direct numericla simulation, DNS

En este proyecto el método utilizado de cálculo es el conocido como DNS, debido a esto se
explicará más a fondo los como funciona este modelo de turbulencia.

Para este modelo las ecuaciones de continuidad instantánea y las ecuaciones de Navier-
Stokes para un flujo turbulento genera un sistema de 4 ecuaciones con 4 incógnitas, DNS usa
estas ecuaciones como punto de partida para generar una solución con una malla lo suficien-
temente pequeña y con saltos temporales lo suficientemente pequeños como para resolver los
torbellinos mas pequeños y los cambios más rápidos.

Los principales beneficios del DNSs son los siguientes[?]:

Detalles precisos de los parámetros de las turbulencias y su transporte en cualquier

punto del volumen de control, permitiendo crear y validar modelos de turbulencia
nuevos de manera más económica.

Página 41 de 57
MEMORIA

Se pueden generan resultados que no se podrían medir con experimentos convencionales

y además se pueden probar estructuras teóricas de turbulencias.

Se pueden también evaluar y probar técnicas experimentales avanzadas y calibrar al-

gunos aparatos de medida.

Se puede investigar turbulencias fundamentales que no se podrían dar en experimentos

de la realidad incluyendo o excluyendo aspectos concretos de los flujos físicos.

En oposición a los múltiples beneficios que se obtinenen de los DNSs como principal des-
ventaja se encuentra la gran complejidad de la solución directa de las ecuaciones del flujo
debido al gran rango de tamaños y escalas temporales debido a la aparición de remolinos
en los flujos turbulentos. Para la resolución de todos los tamaños de turbulencias que apa-
recen con un numero de Reynols pequeños de 104 requerirá una malla de aproximadamente
109 celdas, siendo el orden generalmente de Re9/4 celdas. Además en el caso de las escalas
temporarels se calcula que el numero de time steps es igual a Re1/2 , siendo para un Re de
104 necesario 100 time steps, no obstante en la practica, se suelen usar un número de steps
mayores para asegurar la obtención de estadísticas mas fiables.

Problemas Numéricos
Debido a que esta sección se realiza una breve introducción al CFD y dentro del mismo al
DNS no se va indagar en los detalles del funcionamiento de los métodos usados en DNS, no
obstante en este apartado se va a hace un resumen de las especificaciones específicas para
DNS.

Discetrización Espacial, actualmente se usan varios métodos de desicretización espa-

cial, algunos de estos son los métodos de elementos espectrales, donde se aprovecha
la flexibilidad geométrica del método de elementos finitos con la convergencia de del
método espectral, y los métodos de diferencias finitas de orden superior, siendo
estos segundo los mas usado a día de hoy para geometrías complejas, aunque estos hay
que prestar bastante atención a que el método sea estable para evitar que se disipen
algunos torbellinos turbulentos.

Discretización Temporal, como se menciona antes, hay un rango muy grande de

escalas temporales. Los saltos temporales grandes que se usan generalmente en sistemas
de CFD no son útiles para DNS, estos se debe a que se necesita una gran resolución
para caracterizar todas las variaciones y la disipación de energía a lo largo del tiempo.
Se ha de ajustar los steps tal que una partícula no recorra mas de una celda por step
temporal para conseguir la resolución deseada.

Página 42 de 57
MEMORIA

4.2. Basilisk
Para la realización de las simulaciones que se realizan en este trabajo se ha seleccionado el
software llamado Basilisk. Basilisk se trata de un softwarelibre, es decir es gratuito. Este
es útil para hacer cálculos de complejos de fluidos , ya sean problemas de flujo externo o
interno, combustión o transferencia de calor, por lo que lo hace una muy buena opción frente
a programas de pago y que tienen prestaciones similares.

El código que usa Basilisk se vasa en el lenguaje C++, los códigos de Basilisk se escri-
ben sobre archivos de texto los que se computaran después con software desde la terminal
del ordenador. El lenguaje que usa Basilisk es más compacto que el C++ ya que tiene libre-
rías propias para trabajar cada caso, además reservar algunas palabras para algunas variables.

A parte de los paquetes incluidos en Basilisk, se usan otros paquetes y programas para
el post-proceso y pre-proceso de los datos, como Python y ParaView.

4.2.1. Estructura general código Basilisk

Para describir como funciona un script o código se coge un ejemplo de script de los usados
para la realización de los cálculos en este trabajo. Antes de eso como se ha dicho antes hay
ciertas palabras para algunas variables, siendo estas:
u v w → componentes de la velocidad

p → presión

t → tiempo

N → tamaño de malla

L0 → tamaño físico
Al principio del código se establecen las librerías que se van a utilizar, para ello se utiliza el
comando #include y después se pone el nombre de la librería que se vaya a utilizar, en la
figura 4.4 se ven distintas librerías que se incluyen.

Tras esto se usa el comando # define para definir las constantes numéricas que se vayan
a utilizar, se ve un ejemplo en la figura 4.5
A continuación se define las condiciones de contorno del volumen de control. En la figura 4.6
se establece las distribuciones de los componentes de la velocidad, de la presión y de la frac-
ción volumétrica de fluido, indicando que entre corchetes a que lado del volumen pertenece.

En la siguiente parte del código se establece la función principal del código, introduciendo
las características de los fluidos y de la malla original.

Página 43 de 57
MEMORIA

include.png

Figura 4.4: Ejemplo código #include

define.png

Figura 4.5: Ejemplo código #define

A partir de aquí ya se establecen lo que se llaman eventos, que son funciones que se realizaran
una acción en concreto en la simualción,como guardar los datos deseados en un archivo .txt
como pueden ser el sacar imagenes y guardarlas. Estas funciones siempre comenzaran con
event y después el nombre de la función.

Página 44 de 57
MEMORIA

contorno.png

Figura 4.6: Ejemplo código condiciones de contorno

4.2.2. Softwares adicionales

Durante la realización del trabajo se han utilizado softwares a parte de basislisk, usando
el software de cfd libre de OpenFoam, y software para el pre-procesado de los datos para
basilsik y para el post-procesado de los resultados obtenidos.

OpenFoam: OpenFoam es un software libre para la realización de cálculos CFD tanto

para estudiantes como para proyectos profesionales, en este proyecto se ha utilizado se
han utilizado para obtener los datos a la salida del atomizador mediante la simulación
del flujo interno del atomizador.

Spyder: Spyder es un software de programación que usa una el lenguaje de programa-

ción conocido como python. Este software tiene múltiples aplicaciones, como el desarro-
llo web, uso educativo y de análisis tratamiento de datos numéricos. En este proyecto
se ha usado para el tratamiento de datos numéricos durante el pre-proceso y el post-
proceso.

Página 45 de 57
MEMORIA

int.png

Figura 4.7: Ejemplo código función de lanzamiento

Paraview: Paraview es un software libre de post-proceso de datos para simulaciones

de cfd. Este se ha usado para sacar imágenes y datos para comprobar de manera visual
el que las simulaciones transcurren correctamente.

Página 46 de 57
MEMORIA

snap.png

Figura 4.8: Ejemplo función para obtener capturas de imagen de la simulación.

Página 47 de 57
Desarrollo del trabajo

En este capitulo, se centra en la realización integra del proyecto, es decir, no se van a centrar
en las soluciones si no que se va a exponer los pasos que se han seguido para la obtención de
los resultados finales. Se describen diferentes aspectos como la geometría del atomizador, el
mallado, el volumen de control, la configuración de las simulaciones hechas con el software
Basilisk. También se hará comentarios sobre el pre-proceso y post-proceso realizados con
python y paraview.

5.1. Geometría del atomizador

Como se ha comentado en apartados anteriores, el atomizador utilizado es el conocido como
atomizador tipo pressure-swirl, que se caracteriza con por tener unos canales tangenciales
por donde se introduce el líquido a lo que se conoce como cámara de swirl, con forma cónica
y a partir de aquí salen por el orficio de salida.

swirlesq.png

Figura 5.1: Corte lateral atomizador pressure-swirl

El atomizador que se ha investigado en el proyecto es el Danfoss OD-H 030H8103 oil nozzle

MEMORIA

comercial, en la figura 5.2 se ven imágenes del modelo del atomizador. Este modelo de ato-
mizador tiene un flujo medio de 1,35 kg/h y el ángulo de dispersión es de 80o . El atomizdor
posee tres puertos tangenciales por los que entra el combustible. Como se observa en la figura
5.2 el atomizdor se componen por varias partes metálicas que forman las cámaras del interior
del atomizador.

modeloatom.png

Figura 5.2: Atomizador tipo presure swirl comercial, Danfoss OD-H 030H8103 oil nozzle

Para obtener los datos geométricos de este modelo de atomizador se ha obtenido mediante
diversas técnicas para medirlos. Los datos geométricos obtenidos mediante estas técnicas se
exponen en la figura 5.3 y la tabla 5.1 .

Para esta simulación es muy importante el D0 ya que será la sección por la que se introduce
el fluido al volumen de control, por lo que se va a dar el valor de está, que es de D0 = 300µm.

5.2. Condiciones iniciales

El estudio se centra en la distribución de las gotas a la salida del atomizador, por lo que es
necesario un buen pre-proceso para poder obtener tanto los datos a la salida del atomizador,
así como para obtener un mallado lo más pequeño posible en el volumen de estudio.

Página 49 de 57
MEMORIA

Figura 5.3: Esquema datos geométricos del atomizador Danfoss OD-H 030H8103 oil nozzle

5.2.1. Dominio y mallado

Basilik puede usar mallas cartesianas adaptativa, por eso se usa los que se llama Adaptive
Mesh Refinement(AMR) lo que permite poder resolver con alto detalle problemas con dis-
tintas escalas de resolución. Estas mallas no son fijas para cada timestep, sino que varían
en cada uno adaptándose en función de parámetros previamente establecidos. En el caso de
estudio se ha utilizado como parámetros la diferencia entre celdas de los componentes de
velocidad y de la fracción volumétrica de fluido, para esto se ha establecido unos valores
umbral que indican la máxima diferencia que puede haber entre celdas, y en cada timestep
la malla se refina en diferentes zonas hasta alcanzar los valores umbral en toda la malla,
provocando que haya zonas con celdas de mayor tamaño que otras. Para evitar que refine en
exceso se le impone un nivel máximo de refinamiento en el que dejará de refinar aunque no
haya alcanzado el valor umbral.

Definir el dominio y la malla es fundamental para el gasto computacional necesario y para la

resolución conseguida. La malla utilizada es estructura cartesiana de árbol u octree la cual
se utiliza ya que permiten realizar la AMR, esta estructura cada vez que aumenta un nivel
de resolución divide la celda en 8 celdas [?].

Una vez definida la estrucutra de la malla lo que sigue es definir el tamaño de celda mínimo,
siendo de un 1µm de lado. Este valor se obtenido mediante la premisa que en DNS se ha de

Página 50 de 57
MEMORIA

Cuadro 5.1: Tabla de datos geométricos del atomizador atomizador Danfoss OD-H 030H8103
oil nozzle

Parámetro Unidades Valor

D0 µm 300
Dt−top µm 180
Dt−bottom µm 165
Ds−top µm 210
Ds−bottom µm 500
Dswirl µm 780
Dtotal µm 3375
ht−top µm 240
hs−top µm 540
hswirl µm 620
htotal µm 2400
Rn µm 308
Rt µm 190
o
ψt 6
o
ψs 61
o
ψslot 30

solucionar las todos los niveles de energía hasta la escala de Kolmogorov, debido a esto se
3/4
ha realizado una aproximación a la misma tal que η ∼ LRf uel , obteniendo que la escala de
Kolmogorov η = 0, 95µm, por lo que como se ha comentado antes se usa dx = 1µm [?].A
parte en más literatura se establece la siguiente relación para la elección de dx, dex/µ = 2,1,
en este caso se a utilizado un tamaño menor debido a que gracias al uso de AMR no se tiene
un aumento excesivo del coste de computación al usar un tamaño menor[?],[?]. Para la malla
adaptativa se han escogido los siguientes valor umbrales como límites, siendo la variación
máxima aceptada para cada componente de la velocidad de 1m/s y la variación máxima de
fracción 1e-6. Estos valores han sido escogidos de simulaciones anteriores donde se observo
que estas permiten obtener una malla lo suficientemente fina en las zonas donde interactua
el fluido con el aire.

Para la elección del dominio se ha tenido en cuenta que este tiene que ser lo suficientemente
pequeño para que el numero de celdas no aumente de manera excesiva el coste computacio-
nal, y por otra parte de ha de ser suficientemente grande como para que el flujo se desarrolle
por completo, que permita que se validen los resultados con estudios previos y para tener
suficientes datos para obtener estadísticas sobre las gotas generadas en la atomización pri-
maria. Se han comprobado diversos tamaños de dominio hasta que se ha seleccionado un
dominio cubico de 8 mm de lado. Este dominio se tiene que aumentar ligeramente, tal que
lx,y,z =8,192, debido a como funciona el refinamiento. En cada nivel de refinamiento se divide
en dos cada dirección de la celda, es decir lx,y,z,ref = lx,y,z /2nref , por ejemplo para un punto

Página 51 de 57
MEMORIA

octree.png

Figura 5.4: Ejemplo estructura de árbol 2D, llamada quadtree al dividir en 4 en vez de 8 [?]

que tiene un nivel de refinamiento de 10 las celdas tendrán un lado de lx,y,z /210 = 8µm,
entiendo esto el aumento de tamaño del dominio se debe a que si usamos un dominio de
8mm no se podría conseguir el tamaño mínimo de celda definido, de 1µm, en cambio si el
tamaño es de 8.192 mm con nref =13 se alcanzaría este tamaño de celda, ya que 213 = 8192.
Siguiendo este razonamiento se ha establecido el tamaño del dominio y el nivel máximo de
refinamiento que salen en la tabla . El número de celdas para t=0 es de 2563 , este se ha
puesto exclusivamente para inicializar la simulación ya que el AMR refinará las zonas que
se necesiten.

5.2.2. Condiciones de contorno

En esta sección se le va a establecer las condiciones de contorno así como parámetros previos
necesarios para la simulación.

Página 52 de 57
MEMORIA

El número de CFL o numero de Courant-Friedrichs-Lewy, es menor que uno para asegu-

rar la estabilidad del cálculo, en consecuencia de tener un valor constante de CFL el dt varia
a lo largo de la simulación. Por otra parte el tiempo de simulación total, se ha establecido
para conseguir que haya datos de estacionarios como para poder procesarlos, estableciendo
este tiempo en 400 µs obtenidos de simulaciones previas.

En lo que respecta a las condiciones de contorno se han usado tres tipos de condiciones
de contorno, inflow, wall y outflow, estas se pueden observar la imagen .

Inflow: Esta condición esta puesta en el orificio de salida del atomizador y introduce
valores de velocidad variable u(t) y de la fracción volumétrica de fluido f(t), los cuales
se obtienen del procesamiento de datos de simulaciones previas del flujo interno hechas
en LES, por otra parte la presión p se intruduce mediante una condición de controno
Neumann, de gradiente cero.

Wall: En la cara del cubo donde se encuentra el orifico de salida del atomizador se
impone esta condición de contorno que impone una condición de Drichlet de cero para
la velocidad y la fracción de volumen ya que ambos parámetros se introducen por el
orificio de salida, en cambio la presión tendrá un condición de contorno de Neumand
de gradiente cero al igual que la de inflow.

Página 53 de 57
 Parte II

Pliego de Condiciones
PLIEGO DE CONDICIONES

Índice del Pliego de Condiciones

Página 55 de 57
 Parte III

Presupuesto
PRESUPUESTO

Índice del Presupuesto

Página 57 de 57