Informe Aerodinamica (Otiniano Ramos Nelson Jesús)
Informe Aerodinamica (Otiniano Ramos Nelson Jesús)
Informe Aerodinamica (Otiniano Ramos Nelson Jesús)
SAN MARCOS
(Universidad del Perú. Decana de América.)
FACULTAD DE CIENCIAS FISICAS
AERODINAMICA
https://drive.google.com/file/d/1sjHbUbSKWtI_n2kVhwPZVgnL2fv8D1NA/view?
usp=sharing
ALUMNO:
Nelson Jesús Otiniano Ramos
DOCENTE:
Alvarado Torres Wilson Emilio
2021-II
INTRODUCCION
1) Objetivo General
- Analizar por medio de la simulación computacional el comportamiento
aerodinámico de un vehículo eléctrico utilitario, evaluando diferentes
configuraciones de diseño de la carrocería que disminuyan su resistencia al
avance.
2) Objetivos Específicos
- Estructurar los conceptos relevantes en simulación de flujo computacional con el fin
de realizar un análisis aerodinámico para vehículos de transporte terrestre.
- Estudiar el diseño de detalle de la carrocería del vehículo utilitario para indicar los
aspectos aerodinámicos implementados y sus dimensiones generales.
- Establecer las condiciones de frontera asociadas al modelo a estudiar con el
propósito de preparar el modelo virtual para su posterior análisis en CFD.
- Analizar los resultados obtenidos del estudio en CFD, para determinar las
alternativas de diseño aerodinámicamente más recomendables a implementar en
la carrocería del vehículo eléctrico utilitario.
MARCO TEORICO
1) Fundamentos del análisis aerodinámico de vehículos
La aerodinámica tiene gran influencia en el comportamiento global del vehículo, no
juega únicamente un papel decorativo. Cuando un automóvil empieza el movimiento,
se genera flujo de aire por la superficie exterior del mismo. Este flujo da origen a
fuerzas de arrastre y sustentación, las cuales dificultan el movimiento y la estabilidad
del vehículo. En otras palabras, cada punto en la superficie del vehículo estará
expuesto a fuerzas de rozamiento, tangentes a la superficie del cuerpo, debidas a la
viscosidad del aire, y a fuerzas de presión, normales a la superficie, a causa de la
velocidad relativa existente entre el cuerpo y el fluido.
La fuerza de sustentación, está definida como la componente normal a la superficie del
cuerpo, de la cual depende que las altas velocidades no afecten la estabilidad del
vehículo, haciendo que este se separe del suelo. En aeronáutica, por ejemplo, se
busca que la fuerza de sustentación tenga dirección ascendente para propiciar el
despegue de la aeronave, al contrario del caso vehicular, donde se busca que la
dirección de la fuerza de sustentación sea descendente, para dar agarre y permitir
mayor velocidad al vehículo sin perder estabilidad.
El avance tecnológico ha permitido el desarrollo de herramientas como la Dinámica de
Fluidos Computacional (CFD), en la que se agrupan las ecuaciones de continuidad,
momentum y modelos de turbulencia, permitiendo el análisis aerodinámico de cuerpos
inmersos en fluidos. También llamadas ecuaciones de Navier-Stokes, son la base de la
Dinámica de Fluidos Computacional (CFD). Esta herramienta da solución a las
ecuaciones gobernantes de los fluidos aplicándolas en dominios discretizados del
cuerpo en análisis, facilitando la solución de problemas complejos en fluidos 2D y 3D.
A continuación, se describen algunos parámetros necesarios para realizar el análisis
aerodinámico de un vehículo:
a) Aerodinámica
Rama de la mecánica de fluidos que estudia el comportamiento de los gases y de
su interacción con los cuerpos circundantes. Está fundamentada en el principio de
Bernoulli, el cual establece que la energía de un flujo uniforme es contante en su
trayectoria, por lo tanto, el incremento de velocidad conlleva a una compensación o
disminución de la Presión.
𝑉 + 𝑃 = 𝐶𝑡𝑒
Donde:
V: Velocidad
P: Presión
b) Coeficiente de arrastre
Corresponde a la medida relativa sobre cuanta resistencia al avance puede
generar un cuerpo. El Cx depende de la forma, de la rugosidad de la superficie y
del número de Reynolds.
c) Fuerza de Arrastre Aerodinámica
Es causada por los efectos viscosos del aire y la diferencia de presión presente en
diversas zonas de la superficie de la carrocería en movimiento. Su dirección es
opuesta y paralela al movimiento del cuerpo.
1
𝐹𝑋 = ∗ 𝐶𝑋 ∗ 𝜌 ∗ 𝑉 2 ∗ 𝐴
2
Donde:
CX: Coeficiente de Arrastre
A: Área del vehículo
1 2
∗𝜌∗𝑉 : Presión Dinámica
2
Entre las principales causas de arrastre en un cuerpo se encuentra la forma del
vehículo, la cubierta de la rueda, los aditamentos aerodinámicos, como alerones y
alas, los accesorios como espejos, manijas en las puertas y sellos de ventanas.
e) Eficiencia Aerodinámica
Es la relación entre la sustentación y la resistencia
𝐶𝑍
𝜂=
𝐶𝑋
En vehículos, la eficiencia aerodinámica busca conseguir la menor resistencia
posible con el mayor agarre.
f) Arrastre y Potencia
La relación entre la velocidad máxima y la potencia al freno BHP, se basa en la
ecuación de cálculo de la fuerza de arrastre.
𝐶𝑋 ∗ 𝐴 ∗ (𝑉𝑚𝑎𝑥 )3
𝐵𝐻𝑃 =
1.225
donde:
g) Fluido
Un fluido es una sustancia que tiene la propiedad de fluir, se encuentra en
constante deformación y cede ante cualquier fuerza
h) Flujo Laminar
Bajos números de Reynolds indican predominación de las fuerzas viscosas sobre
las fuerzas inerciales. Las partículas del flujo se desplazan formando capas
paralelas, sin que se lleguen a mezclar de forma significativa. Suele presentarse
para números de Reynolds menores o iguales a 2000.
i) Fujo Turbulento
De carácter caótico, este flujo es intrínseco de números de Reynolds altos,
indicando predominancia de las fuerzas de inercia sobre las fuerzas viscosas. En
este tipo de flujo, las partículas se desplazan de forma desordenada y en
trayectorias irregulares, formando vórtices. Suele presentarse para números de
Reynolds mayores o iguales a 4400.
j) Flujo en Transición
Representa la transición entre el flujo laminar y el flujo turbulento.
m) Numero de Mach
Permite expresar la velocidad de un objeto respecto a la velocidad del sonido. Se
define como:
𝑉
𝑀=
𝐶
Donde:
V: Velocidad del Fluido
C: Velocidad del Fluido en el medio en que se desplaza
a) Conservación de la Masa
𝜕𝜌 𝜕𝜌𝑈𝑖
+ =0
𝜕𝑡 𝜕𝑥𝑖
Donde:
Ui(t,xi): Velocidad en la dirección del eje ‘i’
Xi: Coordenada del eje ‘i’
𝜕𝑈𝑖 𝑈𝑖 1 𝜕𝑃 𝜕 2 𝑈𝑖
+ 𝑈𝑖 ∗ =− ∗ +𝑣∗ + 𝑔𝑖
𝜕𝑡 𝜕𝑥𝑗 𝜌𝑟 𝜕𝑥𝑖 (𝜕𝑥𝑗)2
Donde:
P(t,xi) = Presión
𝜌(𝑡, 𝑥𝑖) = Densidad Local
𝜌𝑟 = Densidad de Referencia
𝑔𝑖 = Gravedad en la dirección xi
c) Modelos de Turbulencia
Los modelos de turbulencia son procedimientos computacionales que sirven para
cerrar el sistema de ecuaciones gobernantes del flujo de fluidos.
La turbulencia es una de las características más relevantes de los fluidos y se
encuentra formada por infinidad de escalas turbulentas. Algunas de estas
representan fluctuaciones muy pequeñas, que demandan tiempo y altos recursos
computacionales. El objetivo de los modelos de turbulencia es desarrollar las
ecuaciones que predicen los campos de velocidad, presión y temperatura de un
fluido, eliminando las escalas más pequeñas de turbulencia, permitiendo la
solución de las ecuaciones constituyentes con cargas computacionales más
pequeñas.
2) Aditamentos aerodinámicos
Adicional a los cambios realizados para la preparación de la geometría, se
agregaron deflectores y carenados traseros o cola de bote (ángulos en la parte
posterior del vehículo). Estos elementos buscan reducir el coeficiente de
arrastre modificando el perfil aerodinámico de la carrocería del vehículo.
Ambos elementos, tanto el deflector como el ángulo trasero, serán probados a
diferente inclinación 𝜃 y 𝛼 respectivamente. El deflector será probado a ángulos
𝜃 equivalentes a 30°, 25° y 20°; alrededor de los cuales se presenta la mayor
reducción del arrastre aerodinámico en el cuerpo de Ahmed. El ángulo trasero
será probado a ángulos 𝛼 equivalentes a 10°, 25° y 60°; cada ángulo se
encuentra dentro de uno de los tres regímenes de flujo generados por colas de
bote en la parte posterior del vehículo según Hanfeng & Cia.
Cabe resaltar que, para este estudio, los elementos aerodinámicos son
completamente sólidos, y no fueron modelados como laminas. Las cavidades
de estos aditamentos son un aspecto a considerar en estudios posteriores.
MODELAMIENTO, MALLADO Y SIMULACION EN CFD
1) Geometría
La Ilustración 15 muestra el resultado de las modificaciones realizadas a la geometría
con el fin de adecuar el modelo a una simulación de flujo. Por otro lado, en la
Ilustración 15, también se pueden apreciar en detalle las zonas en las cuales se
agregaron los elementos aerodinámicos, dichos aditamentos y algunas dimensiones de
la geometría.
3) Mallado
La malla inicial se realizó de manera automática, se realizaron los refinamientos y
cambios en la configuración para mejorar la calidad de la malla. En primer lugar, en
este estudio se realizó análisis de independencia de malla con siete diferentes
configuraciones. La configuración de malla seleccionada presentó un criterio de
convergencia menor al 5% en los resultados. Las características de la malla
seleccionada para el desarrollo se muestran en la Tabla 3.
En segundo lugar, a la malla seleccionada se le realizó un refinamiento alrededor del
vehículo, reduciendo el tamaño de los elementos a 0,02 m. El tamaño de la zona de
refinamiento de los elementos puede visualizarse en la Ilustración 20.
Por último, con el fin de obtener mayor precisión en el comportamiento del flujo
alrededor de la superficie del vehículo, se definió una zona más fina o “Inflation”. Esta
configuración, crea una capa de elementos circundantes al contorno del vehículo que
se alojan sobre la superficie del mismo. Los elementos de esta capa son aún más
pequeños y permiten desarrollar de manera precisa la capa límite aerodinámica. En la
tabla 4, se encuentran los parámetros del “inflation” recomendados.
4) Condiciones de Frontera
La Ilustración 23 muestra las condiciones de frontera que fueron abordadas en el
modelo.
Estas condiciones fueron asignadas a las superficies mostradas, las cuales son
enumeradas en la Ilustración 24.
a) La superficie 1 corresponde a la entrada de aire del dominio, cuya condición de
frontera “Velocity Inlet” (Velocidad de entrada) permite establecer la velocidad y la
dirección del fluido en la entrada del flujo, así como todas las propiedades
escalares relevantes del flujo [8]. Esta condición es aplicable sólo a fluidos
incompresibles. Para este caso, la velocidad máxima del vehículo considerada es
de 40 km/h (11 m/s).
b) La condición “Pressure outlet” (Presión de salida) asignada a la superficie 2, exige
definir la presión estática (manométrica) a la salida del dominio, permitiendo
conocer la presión en cualquier punto al interior de este. Este valor es utilizado
mientras el flujo sea sub-sónico. El estudio se realizó con una presión de 1 atm.
c) La condición “No slip Wall” (No deslizamiento) fue impuesta a las superficies 3 y 7.
Esta condición de frontera es utilizada para unir regiones fluidas y sólidas.
También, permite establecer si los muros tienen movimiento, su rugosidad o
condiciones de deslizamiento o no deslizamiento en la superficie, entre otras. En
este caso, las superficies 3 y 7 representan sólidos, y su condición de no
deslizamiento se debe a que son superficies que atraparán las partículas
adyacentes de fluido. Esto genera, debido a la viscosidad, la formación de la capa
límite aerodinámica.
a) Tipo de solver: Hay dos tipos de solver disponibles en fluent para realizar una
simulación: Basado en presión (Pressure-based) y basado en densidad (Density-
based). El solver basado en presión, se desarrolló originalmente para flujos a baja
velocidad e incompresibles, mientras que el solver basado en densidad se utilizó
principalmente para flujos a alta velocidad y compresibles. Sin embargo,
actualmente, estos enfoques han sido reformulados y ampliados para operar en un
amplio rango de condiciones de flujo más allá de su objetivo tradicional. Ambos
métodos obtienen el campo de velocidad a partir de las ecuaciones de momento.
Sin embargo, en el método basado en presión, el campo de presión es hallado
resolviendo una ecuación de corrección de presión que se obtiene de la
manipulación de la ecuación de continuidad y momento, mientras que, en el
método basado en densidad, la ecuación de continuidad es utilizada para obtener
el campo de densidad y el campo de presión es determinado a partir de la
ecuación de estado. Ambos enfoques ahora son aplicables a una amplia gama de
flujos (desde incompresible hasta altamente compresible). No obstante, el
solucionador basado en presión se ha usado tradicionalmente para flujos
incompresibles o levemente compresibles a bajas velocidades, por lo que es el
solver seleccionado para este estudio.
b) Velocity formulation: Existen dos opciones de formulación de la velocidad: absoluta
y relativa. La formulación de velocidad absoluta es recomendable en aplicaciones
donde el flujo no es rotacional en la mayoría del dominio, mientras que la
formulación de velocidad relativa aplica en el caso contrario, cuando la mayoría del
flujo se encuentra en rotación. En este caso, no hay necesidad de definir marcos
de referencia relativos y el fluido no tiene movimiento rotacional, por lo que la
formulación de la velocidad es absoluta.
c) Time: Las simulaciones en estado estable son aquellas donde las características
del flujo no cambian en un intervalo de tiempo relativamente lago. Por lo tanto, no
se requiere información en tiempo real para describir las características del fluido.
Por otro lado, las simulaciones en estado transitorio requieren información en
tiempo real para conocer las características del campo de fluido. Muchos flujos, en
particular los impulsados por la flotabilidad, no tienen una solución de estado
estable y pueden presentar un comportamiento crítico. Para este estudio, teniendo
en cuenta que las propiedades del fluido (aire) no cambian en un intervalo de
tiempo considerablemente amplio, se realizará la simulación en estado estable.
d) Modelos de viscosidad: En primer lugar, todas las simulaciones fueron iniciadas en
flujo laminar para facilitar la convergencia de la simulación. Posteriormente, al cabo
de unas cuantas iteraciones, el flujo se torna turbulento, cambiando manualmente
el modelo de turbulencia. Este estudio se desarrolló utilizando las dos opciones de
modelos viscosos de dos ecuaciones disponibles en fluent: k-epsilon y k-omega.
Estos dos modelos han demostrado ser precisos y se adaptan a las limitaciones
computacionales presentes en este proceso. El primer modelo ejecutado fue el K-
ε, el cual ha sido utilizado en una amplia gama de estudios aerodinámicos y es
recomendable para corrientes libres de flujo. Sin embargo, ninguna de las
configuraciones de este modelo (RNG, Realizable o Standar) en conjunto con las
funciones de pared disponibles lograron llegar a la convergencia de la simulación,
en algunos casos, el resultado del estudio tendía a la divergencia.
Los resultados del modelo K-ε, los cuales no se muestran en el presente
documento, mostraron resultados anómalos, como flujo en reversa y presiones
absolutas negativas, que no permitían dar continuidad al estudio.
El segundo modelo en ser ejecutado fue k- ω. El modelo K-ω resuelve mejor las
proximidades de la capa limite y tiene un mejor comportamiento cerca de las
paredes.
La configuración k- ω estándar fue corrida obteniendo resultados no satisfactorios.
Los resultados de las variables monitoreadas tenían un comportamiento similar al
obtenido con el modelo K-ε. Posteriormente fue ejecutado el modelo k- ω SST, el
cuál es la base del modelo de turbulencia SST (Shear stress transport) de cuatro
ecuaciones. El modelo k- ω SST logró estabilizar y converger los residuales con un
criterio de 10−3 en dos mil iteraciones, exceptuando los residuales de continuidad,
los cuales teniendo en cuenta su pendiente negativa, requieren de un mayor
número de iteraciones y tiempo de cómputo para llegar a la convergencia con este
criterio. El modelo k- ω SST es preciso y logra simular un comportamiento
razonable del fluido alrededor del vehículo. Por esta razón, este último es el
modelo de turbulencia seleccionado para el estudio.
Otro tipo de modelos diferentes al viscoso no fueron considerados y no son
relevantes para lograr el objetivo de este estudio.
e) Materiales: El fluido que circula alrededor del vehículo es aire, cuya densidad es
1,225 kg/m3 y cuya viscosidad es 1,7894 e-5 kg/ms en condiciones normales. El
acoplamiento Presión velocidad: El acoplamiento presión-velocidad se consigue a
partir de la ecuación de continuidad discretizada. Esta última, permite obtener una
ecuación para la presión. Fluent permite escoger entre tres tipos de esquemas de
acoplamiento presión velocidad: SIMPLE, SIMPLEC y PISO. Ansys Fluent
recomienda utilizar los algoritmos SIMPLE o SIMPLEC para simulaciones en
estado estable, mientras que para simulaciones en estado transitorio se
recomienda el esquema PISO. El algoritmo SIMPLE emplea una relación entre la
velocidad y la correlación de la presión para hacer que se cumpla la ecuación de la
conservación de la masa y obtener así, el campo de presiones. Sobre este
esquema se han variaciones y reformulaciones para obtener algoritmos como el
SIMPLEC. Este último es recomendable para flujos complicados que involucran
turbulencia porque acelera la convergencia en problemas en los que el
acoplamiento presión-velocidad es el principal obstáculo para llegar a la solución.
Por otro lado, todas las simulaciones fueron realizadas en primer orden durante las
primeras 100 iteraciones. A partir de este valor, la presión, el momentum, la
energía cinética turbulenta y la tasa de disipación de energía fueron cambiadas a
segundo orden con el fin de obtener mayor precisión en la solución.
RESULTADOS
1) Residuales de convergencia
Como se puede observar en la Ilustración 25, los residuales de convergencia
correspondientes a la velocidad en x, y, z y omega alcanzaron valores convergentes en
dos mil iteraciones con el criterio de 10-3. Los residuales de continuidad y k1, tienen
tendencia decreciente y para propósitos académicos el resultado es aceptable. Debido
al tiempo disponible para la ejecución de las simulaciones, no fue posible lograr un
mayor refinamiento en la malla o incrementar el número de iteraciones para llevar al
criterio de convergencia a todas las variables monitoreadas. Se espera que los
resultados obtenidos puedan ser mejorados en la medida en que los recursos
computacionales sean enriquecidos.
Es importante resaltar que el coeficiente de arrastre del vehículo con deflector está
por debajo del coeficiente de arrastre teórico. Por otro lado, el deflector con ángulo
de inclinación de 25º logro reducir el 29,4% del coeficiente de arrastre presentado
para el vehículo sin elemento aerodinámico, porcentaje de reducción mayor al
encontrado por Chilbule & Cia. Es pertinente entonces, estudiar el comportamiento
del coeficiente de arrastre con otros elementos aerodinámicos y determinar si es
posible reducirlo aún más.
a) Contornos de Presión
La Ilustración 39, Ilustración 40 e Ilustración 41 muestran los contornos de presión
alrededor y sobre el vehículo para los tres diferentes casos del elemento
aerodinámico: ángulo trasero.
Puede observarse en la Ilustración 39 que el ángulo trasero de 10° reduce el
gradiente de presión en la parte posterior del vehículo. En esta región trasera, las
presiones oscilan entre -12,2 Pa y 17,5 Pa utilizando el ángulo de 10°, lo cual
genera un contorno de presión más uniforme que en los otros dos casos.
Es importante resaltar que el vehículo con ángulo trasero de 10° fue el único caso,
de los tres casos simulados con ángulo trasero, en el que se obtuvo disminución
considerable en el coeficiente de arrastre respecto al vehículo sin elementos
aerodinámicos. Este elemento logró reducir el 14,7% del coeficiente de arrastre
respecto al vehículo básico. A pesar de que no fue un porcentaje de reducción
mayor al de los casos con deflector, es recomendable asociar los elementos
aerodinámicos que mostraron los mejores resultados de Cx, y determinar si es
posible reducir aún más este coeficiente a partir de la vinculación de diversos
elementos.
CONCLUSIONES
1) Conclusiones
La recopilación de conceptos relacionados con aerodinámica y dinámica de fluidos
computacional permitió reconocer los parámetros que definirían el arrastre del vehículo
para los casos estudiados. Por otra parte, se observaron los aspectos aerodinámicos
que ya se encontraban implementados en la carrocería y se identificaron las posibles
regiones donde sería necesario implementar otros aditamentos. Además, la inspección
de la carrocería permitió reconocer la necesidad de realizar modificaciones a las
llantas, con el fin de conseguir que el modelo se acercara más al comportamiento real
de un vehículo en movimiento. El modelo de turbulencia y las condiciones de frontera
establecidas prepararon el modelo virtual para el posterior análisis en CFD. La
configuración del solver realizada conllevó a resultados coherentes, no se alejan de los
datos usualmente obtenidos en análisis aerodinámicos para este tipo de automotores.
En este trabajo se logró determinar la alternativa aerodinámicamente más
recomendable a implementar en el vehículo eléctrico utilitario. Los resultados arrojaron
contornos de presión y campos vectoriales de velocidad que permitieron identificar que
el vehículo con deflector de 25° y cola de bote o ángulo trasero de 10° es la
configuración que presenta menor presencia de vórtices y mayor nivel de orden en el
comportamiento del flujo de los casos estudiados. Además, se obtuvieron los valores
exactos de los coeficientes de arrastre para cada caso, con el fin de encontrar
cuantitativamente la opción que representara menor consumo de combustible.
2) Aportes
a) Se realizó el análisis aerodinámico del vehículo en construcción del proyecto
“diseño y construcción virtual de un vehículo eléctrico utilitario”. Este análisis
entrego los componentes aerodinámicos que pueden reducir el coeficiente de
arrastre del vehículo en un 32%, convirtiéndolo no solo en una alternativa funcional
y utilitaria de transporte de carga o pasajeros, sino también, en una alternativa
eficiente desde el punto de vista aerodinámico y amigable con el medio ambiente.
b) Se obtuvieron los contornos de presión y campos vectoriales de velocidad para el
vehículo eléctrico utilitario con y sin elementos aerodinámicos. Estas ilustraciones
permitieron analizar el comportamiento del fluido alrededor del vehículo en
movimiento, para posteriormente, lograr reducir las presiones extremas y el caos
de flujo.
c) Se desarrollaron simulaciones en Ansys Fluent del modelo del vehículo eléctricos
utilitario a partir de las cuales se obtuvieron los coeficientes de arrastre de cada
uno de los casos estudiados en este documento para la velocidad máxima de
operación del vehículo.
BIBLIOGRAFIA
1) María A. Chadid Gutiérrez (2018). Análisis Aerodinámico Mediante Dinámica
de Fluidos Computacional de un Vehículo Eléctrico Utilitario. Facultad de
Ingeniería Mecánica. Universidad Tecnológica de Pereira.