Airfoil mesh generator for OpenFOAM, V0.2. Algunos errores corregidos.

En la última entrada tenía listo el esqueleto del código para generar mallas bidimensionales de perfiles para OpenFOAM a partir de Python, pero incluso en las imágenes que colgué podían apreciarse errores en la malla.

Cosas que se han corregido/añadido:

– Se ha corregido el pequeño desplazamiento en el eje X de la malla que puede verse en la imagen. Punto (1)

Errores en la malla del primer código. Puntos en rojo.

– Control dinámico del reparto de celdas sobre el borde de ataque del perfil y el resto del cuerpo. Esto hace que las mallas sean aptas (a priori) para condiciones supersónicas ya que la distribución puede cambiarse para capturar las ondas de choque creadas en el borde de ataque con gran precisión. Punto (2)
– Inclinación dinámica de la estela. La malla se adecua al ángulo de ataque del perfil. De esta forma no es necesario cambiar las condiciones iniciales de la velocidad en los diccionarios de las variables de OpenFOAM. Esto es opcional.

En las imágenes pueden verse los resultados de estos puntos, salvo el de la estela porque no se aprecia bien en esta pantalla (se admiten regalos y donaciones para una nueva, 24″ full HD)

Se han corregido los errores marcados con (1) en la imagen anterior.

Distribución dinámica de la malla sobre el borde de ataque para capturar efectos. Nota: Se ha acentuado con efectos ilustrativos.

Cosas por hacer:

– Encontrar una relación matemática correcta que haga que las transiciones entre la malla del borde de ataque y la correspondiente al resto del perfil sea suave.
– Implementar un control de y+, especificar el y+ que se busca y que el código se encargue de ello.
– Crear una zona adicional de control en el borde de salida.

Tareas en mente.
– Validar la malla creada con algunos resultados experimentales para los perfiles más comunes, Naca 0012 p.e.
– Refinar los parámetros de densidad de la malla mediante un algoritmo de optimización -tengo uno genético listo, pero creo que el de gradiente va a funcional mejor- valorando el peso de la malla (tamaño o número de celdas) y la aproximación de los resultados al modelo real.

Esta entrada fue publicada en Aerodinamica, cfd, matemáticas, OpenFOAM, python y etiquetada , , , , , . Guarda el enlace permanente.

Responder

Introduce tus datos o haz clic en un icono para iniciar sesión:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión / Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión / Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión / Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión / Cambiar )

Conectando a %s