Método de la Red de Vórtices
(VLM – Vortice Lattice Method)
Bibliografía
▪ The Calculation of Aerodynamic Loading on Surfaces of any Shape –
Falkner, V.M. – Reports and Memoranda Nº1910 - 1943
▪ Aerodynamics for Engineer – Bertin, J.J. y Cummings, R.M. – 2009.
•Introducción
Este es un método similar al método de los paneles pero muy simple de
usar, capaz de proveer muy buenos resultados en el cálculo de la
aerodinámica de alas y la interacción entre componentes.
Es el denominado Vortex Lattice Method (VLM).
Este método esta basado en resolver la ecuación de Laplace, y está
sujeto a las mismas restricciones que el método de los paneles.
Compararemos, ahora, las similitudes y diferencias entre ambos métodos.
Similaridades con paneles:
➢las singularidades son ubicadas en la superficie.
➢La condición de tangencia (el flujo no atraviesa la superficie) debe ser satisfecha en los puntos de
control.
➢Para encontrar la intensidad de las singularidades se resuelve un sistema de ecuaciones
algebraicas lineales.
Diferencias con paneles:
➢ esta orientado a hacia los efectos de sustentación, y las formulaciones clásicas ignoran
espesor.
➢ Las condiciones de borde (CB) son aplicadas a la superficie media, no a la real (no son una
solución exacta de la ecuación de Laplace sobre el cuerpo, pero incluye algunas
aproximaciones adicionales, encontramos el Cp, no la distribución en el intradós y extradós.
➢ Las singularidades no están distribuidas sobre la superficie completa.
➢ Esta orientado hacia combinaciones de las superficies sustentadoras delgadas (el método de
los paneles no tiene limitaciones debido al espesor)
El VLM fue formulado por primera vez en la década del ‘30, y fue llamado en su
principio “Vortex Lattice” en 1943 por Faulkner.
El concepto en el cual se basa es extremadamente simple puesto que es puramente
una aproximación numérica.
Hipótesis empleadas:
Flujo estacionario, NO viscoso, irrotacional e incompresible.
Se representa el ala por una superficie plana sobre la
cual se superpone una grilla de vórtices en herradura.
La velocidad inducida por cada vórtice en herradura en
un punto de control específico es calculada usando la
Ley de Biot-Savart.
Se suma para todos los puntos de control del ala para
producir un conjunto de ecuaciones algebraicas lineales
para la intensidad de los vórtices en herradura que
satisfacen la condición de borde de tangencia (el flujo no
cruza el ala)
La intensidad de los vórtices se relacionan con la
circulación del ala y la diferencia de presiones entre las
superficies inferior y superior.
Las presiones se integran para obtener las fuerzas y
momentos totales.
Se aproxima la distribución continua de vórtices del
cuerpo del ala por un número finito de vórtices en
herradura.
• Cada vórtice en herradura se coloca en paneles
trapezoidales.
• Cada uno de estos se ubica al ¼ de cuerda del panel y
esta alineado con la flecha local.
• Los paneles se colocan en la superficie media del ala y
los vórtices de fuga dejan el ala siguiendo la curvatura.
• Sin embargo, se pueden obtener buenos resultados
utilizando la teoría linealizada, donde los vórtices de
fuga salen en línea recta hacia el infinito.
• Pueden ser paralelos a la velocidad de la corriente libre
o al eje x (en este caso se utilizará este último caso,
pues el cálculo de los coeficientes de influencia es más
simple, con una ventaja no cambian con el ángulo de
ataque).
Condiciones:
• El flujo es tangente a la superficie del ala en cada
punto de control de los 2N paneles.
• Esto genera un conjunto de ecuaciones en las cuales
las incognitas son la intensidad de la circulación de los
vórtices.
Ubicación de los puntos de control
Los puntos de control se ubican a la mitad entre los vórtices de fuga y a ¾ de la cuerda del panel.
Supongamos un filamento vorticoso a ¼ de cuerda de intensidad , la velocidad inducida será:
U=
en el punto c a una distancia r. 2 r
Si el flujo debe ser paralelo a la superficie en el punto de control, entonces la incidencia de la
superficie relativa a la corriente libre será:
U
sen = =
U 2rU
Para placa plana
Cl = 2
Reemplazando
1
l = U 2 c.2 = U
2
c c c 3c
r = c/2 +r = + = = 0,75c
4 4 2 4
Velocidad Inducida
La velocidad inducida por un filamento de intensidad n y longitud dl, esta dado por la Ley de Biot
y Savart.
n sendl
con dV =
4r 2
Como para cada vórtice tenemos tres segmentos, aplicamos lo anterior a cada segmento
2 n
V= n send = (cos1 − cos 2 )
4rp 1 4rp
Siendo 1= 0 y 2 =
Reemplazando en la anterior, considerando la dirección del vector velocidad unitario como,
(válida para cualquier dirección)
Esta se puede utilizar para calcular la velocidad inducida por un vórtice en herradura (sin importar
la orientación).
Supongamos, ahora, un panel n cualquiera.
(x,y,z)
La suma de las anteriores para un punto (x, y ,z) nos da la velocidad inducida allí por un vórtice en
herradura del panel n.
velocidad inducida en el punto de control m por el vórtice en el panel n.
Cm,n es el coeficiente de influencia (depende de la geometría del vórtice en el panel n y la
distancia al punto de control del panel m)
La velocidad inducida total en el punto de control m (del panel m) debido a los 2N vórtices será:
(tendremos 2N ecuaciones)
Condiciones de contorno
Entonces si conocemos la intensidad de n de cada vórtice podemos calcular la velocidad
inducida en cada punto del espacio.
Para conocer n debemos utilizar la condición de borde que indica que la superficie es una
línea de corriente. Es decir el flujo es tangente al ala en cada punto de control.
Si el flujo es tangente al ala, la componente normal de la velocidad inducida en cada punto de
control se cancela con la componente normal de la velocidad de la corriente libre.
Para evaluar las componentes de la velocidad inducida debemos introducir la convención de que
los vórtices de fuga son paralelos al eje x (eje del ala).
La relación de tangencia será:
− um sen cos − vm cos sen + wm cos cos
+ U sen( − ) cos = 0
dz
= arctg
m
dx
Para el caso de curvaturas pequeñas y con
ángulos de ataque pequeños.
Se puede aproximar a:
dz
wm − vmtg + U − = 0
dx m
Luego para encontrar las intensidades se debe
resolver el sistema planteado anteriormente.
Para el caso de un ala plana
Siendo los z = 0 para todos los puntos nos quedarán las siguientes expresiones:
Representa la velocidad en la componente normal al plano en el punto de control m debido al
panel n
Valor total
Siendo el ángulo de diedro
nulo y la pendiente nula
Para pequeños ángulos de
ataque
Ejemplo: ala plana con flecha
Por simetría se puede calcular en una semi-ala
(plano xz).
Se calculan solo los puntos de control de una
semi-ala pero considerando todos los vórtices del
ala.
N N
wm = wm.ns + wm.np
n =1 n =1
Se calculan estas para cada panel y luego se
aplica la condición de tangencia en cada uno.
Luego se obtendrá un sistema de ecuaciones a
resolver siendo las incógnitas las intensidades.
1 + 0 , 5b 1 + 0 , 5b Cl c 1 Cl c Cl c
= Cl c i dy −0,5b ( y − )2 d
0 , 5b
C Dv
S − 0 , 5b
i = −
8
i = −
8
0
+ 2
( y − ) ( y + )
2
d
