INTRODUCCIÒN

Una discretización (descomponer un elemento complejo en varios elementos más simples

como líneas o planos) parte necesariamente de un modelo continuo.

Obtener una distribución continúa de presiones y hasta de algo más simple que son los
vórtices (basadas en una o varias funciones matemáticas) es una operación sumamente
complicada, por eso utilizamos el método de vórtices discretos, este método es muy
eficiente para poder estudiar perfiles en general y además obtener resultados preliminares
problemas relacionados con interacciones bidimensionales en perfiles multielementos y la
interacción entre el suelo y el perfil.

El perfil alar representa el ADN de la plataforma del ala y otorga las características
fundamentales del desempeño aerodinámico en 2-D, medido en coeficiente de sustentación
y coeficiente de resistencia. El propósito de estudiar los perfiles NACA es realizar una
breve introducción a la nomenclatura general NACA de un perfil alar, sus diferencias
geométricas y su valoración aerodinámica.
MARCO TEÒRICO
Vòrtices
El vòrtice se observa como la rotación de un fluido alrededor de un eje de tal forma que a
mayor distancia del eje, la velocidad de este disminuye de forma proporcional:
ω. r = K
Γ
donde ω = 2π.r

Figura:
En general podemos afirmar que los vòrtices se forman cuando entre dos capas de fluido
existe una diferencia de velocidades, o lo que es igual, cuando existe una discontinuidad de
velocidad entre dos capas fluidas (un caso en particular sería que una de las capas estè en
reposo y otra a cierta velocidad); al siguiente figura muestra una interacción entre una masa
de aire en movimiento A con una masa en reposo B, en el caso de la segunda situación, la
masa en reposo sería aquella que està justo detrás de la circunferencia.

Figura:
Perfiles Naca

La NACA fue una agencia federal de Estados Unidos que se fundó el 3 de marzo de 1915
para emprender, fomentar, e institucionalizar las investigaciones aeronáuticas. Este comité
estaba encargado de realizar pruebas y diseños empíricos de perfiles alares. La agencia se
disolvió el 1 de octubre de 1958 para darle paso a la creación de la NASA (National
Aeronautics and Space Administration).

Estas dos instituciones a nivel mundial formalizaron trabajos y estudios en túneles de

viento durante el período definido entre 1920-1960. Este trabajo resultó en la creación de
familias de perfiles alares y el desarrollo racional de una nomenclatura para la escogencia
del perfil basado en su desempeño aerodinámico. El estudio tiene tanta trascendencia e
importancia, que muchos de los perfiles o standard "Airfoils" son utilizados en aeronaves
hoy en día.

Fundamentos teóricos de la nomenclatura de un perfil alar

Cabe anotar que antes de 1930, los diseñadores aeronáuticos no tenían ningún tipo de
herramienta para seleccionar un perfil alar, por tal motivo, su diseño y escogencia era
personalizado y no obedecía a un análisis dimensionado o consistente en el campo
analítico. (Anderson, 1999).

Esta situación cambio de manera radical cuando la NACA decidió adoptar un estudio
racional y sistemático para el diseño de los perfiles alares aeronáuticos. Este estudio se
basó en pruebas exhaustivas en túneles de viento para el análisis de las características y
propiedades inherentes de cada perfil alar. El estudio requería de una nomenclatura
efectiva para lograr su clasificación de manera ordenada y detallada con el fin de
suministrarle al usuario final una validación completa del perfil alar acorde a su misión.

Concepto geométrico del perfil alar.

La nomenclatura NACA define el concepto geométrico del perfil alar. A continuación se

muestra en la Figura 2, el estándar de lo que se considera un perfil aerodinámico con sus
definiciones geométricas.

La Figura anterior, corresponde al ADN y la definición del desempeño aerodinámico en

análisis 2D. Su definición y por ende su característica más importante es la línea de camber
conocida como "Mean Camber Line" o línea de curvatura media. (Anderson, 1999).

Esta línea de curvatura media es la característica de diseño más importante del perfil alar.
Se define por un grupo de puntos posicionados de manera equidistante entre el extradós y
el intradós con cambio en la longitud de cuerda, medidos perpendicularmente desde la
cuerda.

El punto de inicio y el punto de finalización de la línea de curvatura media están definidos

como el borde de ataque y el borde de salida. La línea recta que conecta el borde de ataque
y el borde de salida es la línea de cuerda del perfil alar y la distancia precisa medida en
esta línea se conoce como la cuerda.

 Borde de Ataque. Es el punto más delantero del perfil.

 Borde de Salida. Es el punto más trasero del perfil.
 Línea de la cuerda: Es la línea recta que pasa por el borde de ataque y por el borde
de salida.
 Cuerda: es la línea recta que une el borde de ataque con el borde de salida. Es una
dimensión característica del perfil.
 Línea de Curvatura Media. Línea equidistante entre el extradós y el intradós. Esta
línea “fija” la curvatura del perfil. Si la línea de curvatura media “queda” sobre la
cuerda (como en la figura 2) se dice que la curvatura es positiva, si queda por
debajo, negativa, y si va por debajo y por arriba, doble curvatura.
 Ordenada máxima de la línea de curvatura media. Es la máxima distancia entre la
línea de curvatura media y la cuerda del perfil. El valor suele darse en porcentaje
de la cuerda.
 Posición de la curvatura máxima: es la distancia media a partir del borde de ataque,
en porcentaje de la cuerda, donde se encuentra la ordenada máxima de la línea de
curvatura media.
 Espesor máximo: es la distancia máxima entre el extradós e intradós, medida
perpendicularmente a la cuerda. Es una característica importante, que se expresa
en porcentaje de la cuerda, el valor varía desde un 3% en los perfiles delgados hasta
un 18% en los más gruesos.
 Posición del espesor máximo: es la distancia paralela a la cuerda, medida desde el
borde de ataque hasta la ordenada donde existe el espesor máximo del perfil
 Radio de Curvatura del Borde de Ataque. Define la forma del borde de ataque y es
el radio de un circulo tangente al extradós e intradós, y con su centro situado en la
línea tangente de curvatura media y que pasa por el borde de ataque
-El camber de un perfil alar es la distancia máxima que existe entre la línea de cuerda y la
línea de curvatura media o "mean camber line" medida perpendicularmente desde la
cuerda. El camber, la forma de la línea de curvatura media y la distribución de espesor del
perfil alar en esencia representan las variables de diseño más importantes porque controlan
las características de sustentación y los momentos del perfil alar.

-El radio de la nariz del borde de ataque tiene un radio de alrededor de un 1% de la longitud
de cuerda en aplicaciones subsónicas. Este radio en efecto es el radio de un círculo centrado
en una línea tangente al camber del borde de ataque conectando puntos en tangente entre
la superficie superior y la superficie inferior. La magnitud del radio del borde de ataque
tiene consideraciones en el tipo de pérdida o "stall" del perfil. Cabe anotar que con
incremento en magnitud del radio la pérdida de sustentación o "stall" es más dócil,
mientras que con la reducción de la magnitud del radio, la pérdida de sustentación es más
dramática y la pendiente es más aguda.

Por último, la distribución del espesor del perfil alar tiene efectos aerodinámicos
importantes relacionados con el coeficiente máximo de sustentación y la capa límite
laminar. Con incremento en espesor, el coeficiente máximo de sustentación aumenta y si
el espesor disminuye el coeficiente máximo de sustentación disminuye.

De igual manera, hay un fenómeno de separación asociado con el incremento en espesor

que incrementa la gradiente de presión adversa en el perfil alar. Por este motivo, es
necesario tener en cuenta hasta qué punto es bueno incrementar el espesor del perfil con
el objetivo de incrementar el coeficiente de sustentación máximo y de igual forma
minimizar la separación de la capa límite.

Clasificaciòn de los perfiles aerodinàmicos

Perfil alar NACA 4 dígitos

NACA decidió generar la nomenclatura para la primera familia de perfiles alares a principios
de 1930. La primera familia de los NACA airfoils se conoce como la familia de perfiles de 4
dígitos.

En base a los valores de m, p y t se pueden calcular las coordenadas para un perfil completo
siguiendo los siguientes pasos:

1.- Elegir valores de x desde 0 hasta la cuerda máxima (c).

2.- Calcular las coordenadas de la línea de combadura media introduciendo los valores de m
y p dentro de las siguientes ecuaciones para cada una de las coordenadas de x.
𝑚
𝑦𝑐 = 𝑝2 (2𝑝𝑥 − 𝑥 2 ) … 𝑑𝑒𝑠𝑑𝑒 𝑥 = 0 ℎ𝑎𝑠𝑡𝑎 𝑥 = 𝑝

𝑚
𝑦𝑐 = (1−𝑝)2 ((1 − 2𝑝) + 2𝑝𝑥 − 𝑥 2 ) … 𝑑𝑒𝑠𝑑𝑒 𝑥 = 𝑝 ℎ𝑎𝑠𝑡𝑎 𝑥 = 𝑐

Donde :

x: Coordenadas a lo largo de la longitud del perfil, desde 0 hasta c, las cuales se posicionan
en la cuerda.

y: Coordenadas por encima y debajo de la línea extendiéndose a lo largo de la longitud del

perfil, estas son tanto para las coordenadas del espesor 𝑦𝑡 como para las coordenadas de la
combadura 𝑦𝑐 .

t: Espesor máximo del perfil en dècimas de la cuerda (por ejemplo, en un 15% de espesor
del perfil debe ser 0.15).
m: Corbadura máxima del perfil en dècimas de la cuerda.

p: Posiciòn de la combadura màxxima a lo largo de la cuerda en dècimas de la cuerda.

3.- Calcular la distribución de espesores por encima (+) y por debajo (-) de la línea media
introduciendo el valor de t dentro de la siguiente ecuación para cada una de las coordenadas
de x.
𝑡
𝑦𝑡 = 0.2 (0.2969√𝑥 − 0.1260𝑥 − 0.3516𝑥 2 + 0.2843𝑥 3 − 0.1015𝑥 4 )

4.- Determinar las coordenadas finales de la superficie superior (𝑥𝑢 , 𝑦𝑢 ) y de la superficie

inferior (𝑥𝐿 , 𝑦𝐿 ) empleando las siguientes relaciones:

𝑥𝑢 = 𝑥 − 𝑥𝑡 𝑠𝑒𝑛𝜃

𝑦𝑢 = 𝑦𝑐 − 𝑦𝑡 𝑐𝑜𝑠𝜃

𝑥𝐿 = 𝑥 − 𝑦𝑡 𝑐𝑜𝑠𝜃

𝑦𝐿 = 𝑦𝑐 − 𝑦𝑡 𝑐𝑜𝑠𝜃

𝑑𝑦
Donde: 𝜃 = arctan ( 𝑑𝑥𝑐)

A continuación se presenta la descripción de dicha nomenclatura utilizando como ejemplo el

perfil alar NACA 2412 de 4 dígitos (Ver Figura 3).
Figura: Perfil Alar Naca de 4 Dígitos Nomenclatura Autor / Figura del Perfil Alar y Su
Geometría. Fuente: Carmona, A. (1999). Aerodinámica y Actuaciones del Avión, p. 79.
Madrid: Thomson-Paraninfo.

De acuerdo a la definición NACA, el primer dígito da el camber máximo (m) en porcentaje

de cuerda. Esto significa que el camber máximo del perfil 2412 es 2% de la longitud de
cuerda. El segundo dígito representa la localización del camber máximo (p) en la décima de
cuerda. En este caso, para el perfil 2412 el máximo camber está localizado en el 40% de la
longitud de la cuerda medido desde el borde de ataque. Los últimos dos dígitos otorgan el
máximo espesor en porcentaje de la cuerda (t). Cabe anotar, que el máximo espesor de un
perfil alar de cuatro dígitos se ubica por defecto al 30% de la cuerda medido desde el borde
de ataque. Por ejemplo, el perfil NACA 2412 tiene un máximo grosor de 12% de longitud de
cuerda @ 30%C.

Es importante mencionar el perfil alar simétrico comúnmente utilizado en superficies de

control como el estabilizador horizontal y el estabilizador vertical (Rudder y Elevador). Este
perfil simétrico se define en la familia de 4 dígitos, pero los primeros dos dígitos están
designados como ceros. El perfil simétrico no tiene camber, solo tiene espesor. Por ejemplo,
un perfil simétrico se define como se indica en la Figura 4.
Figura: Comparación Perfil NACA 2412

La Figura 5 realiza una comparación ilustrativa de la diferencia geométrica entre un perfil

NACA 2412 y un perfil simétrico NACA 0012 con el mismo espesor.

Figura: Comparación Perfil NACA 2412 y Perfil NACA 0012. Fuente: Anderson, J. (1999).
Aircraft Performance and Design. New York: McGraw-Hill.

Perfil alar NACA 5 dígitos

A mediados de 1930 nació la segunda familia de los perfiles alares NACA. Esta familia se
denominó NACA de 5 dígitos. Este tipo de perfil alar surgió de experimentos y estudios
empíricos con el objetivo de incrementar el coeficiente de sustentación máximo. Los
resultados demostraron que si acercaban el punto de camber máximo en dirección al borde
de ataque se obtenía incremento en el coeficiente deEl primer dígito cuando se multiplica por
3/2 otorga el sustentación local. Por tal motivo, se crea el perfil alar NACA de 5 dígitos como
por ejemplo el NACA 23012 (Ver Figura 6).

El primer dígito cuando se multiplica por 3/2 otorga el coeficiente de sustentación de diseño
en décimas. El coeficiente de sustentación de diseño es un índice del incremento en camber
y por ende, si el camber del perfil incrementa, el coeficiente de sustentación de diseño se
incrementa. En este caso, el coeficiente de sustentación de diseño para el perfil NACA 23012
es de 0.3.
Figura: Perfil Alar NACA de 5 Dígitos Nomenclatura Autor / Figura del Perfil Alar y Su
Geometría. Fuente: Carmona, A. (1999). Aerodinámica y Actuaciones del Avión, p. 79.
Madrid: Thomson-Paraninfo.

El segundo y tercer dígito forman un número que multiplicado por ½ otorga la ubicación del
camber máximo en porcentaje de la longitud de línea de cuerda tomada desde el borde de
ataque. En este caso se tiene el camber máximo al 15% de longitud de cuerda. Los últimos
dos dígitos dan el grosor máximo del perfil alar que en este caso es del 12% o 0.12c de la
longitud de cuerda.

Perfil alar NACA serie 1

El perfil alar serie 1 representó el primer intento de diseñar teóricamente un perfil alar con
un tipo de distribución de presión deseado, mantener una presión de gradiente favorable y
una capa límite laminar en el proceso.

A continuación se presenta el ejemplo de nomenclatura con el perfil NACA 16-212. (Ver

Figura 7).
Figura: Perfil Alar NACA de 1 Dígito Nomenclatura Autor / Figura del Perfil Alar y Su
Geometría. Fuente: Carmona, A. (1999). Aerodinámica y Actuaciones del Avión, p. 79.
Madrid: Thomson-Paraninfo.

El primer dígito identifica la serie, el segundo dígito provee la ubicación del punto de mínima
presión en décimas de cuerda medidas desde el borde de ataque. Cabe anotar que la ubicación
del punto de mínima presión es analizado no para el perfil con camber, sino para un perfil
simétrico con la misma distribución de espesor del perfil analizado y al ángulo de ataque
correspondiente a cero sustentación o ángulo de ataque cero para el perfil simétrico. Para el
perfil NACA 16-212 este fenómeno ocurre a 0.6c o 60% de la longitud de cuerda. El tercer
dígito da el coeficiente de sustentación de diseño en décima de cuerda. En este caso el
coeficiente de sustentación de diseño es de 0.2. Finalmente, los últimos dos dígitos otorgan
el máximo grosor en porcentaje de la cuerda, que para el NACA 16-212 sería del 12%.

Perfil Alar NACA Serie 6

A finales de 1930 y principios de 1940 la NACA empieza con los estudios de perfiles alares
con el propósito de fomentar el flujo laminar sobre su superficie minimizando el coeficiente
de resistencia. Este nuevo alcance de la NACA se le llama "Laminar Flow Airfoil" o perfil
alar de flujo laminar. Este estudio involucra a una nueva familia de perfiles denominada serie
6 o NACA 6-series airfoil.
A diferencia de los objetivos del perfil alar de 1 dígito, El NACA serie 6, tiene como objetivo
minimizar la resistencia asegurando el flujo laminar y reduciendo la gradiente de presión
adversa, obtener un Mach number crítico deseado para la reducción de resistencia en régimen
transónico y por último aumentar las características de máxima sustentación.

A continuación se presenta el ejemplo de nomenclatura con el perfil alar NACA 641-212 que
adicionalmente posee un subíndice donde indica gradientes de presión favorables por arriba
y por abajo del coeficiente de sustentación de diseño.

Figura: Perfil Alar NACA de 6 Dígitos Nomenclatura Autor / Figura del Perfil Alar y Su
Geometría. Fuente: Carmona, A. (1999). Aerodinámica y Actuaciones del Avión, p. 79.
Madrid: Thomson-Paraninfo.

El primer dígito identifica la serie, el segundo dígito provee la ubicación del punto de mínima
presión en décimas de cuerda medidas desde el borde de ataque.

Cabe anotar que la ubicación del punto de mínima presión es analizado no para el perfil con
camber, sino para un perfil simétrico con la misma distribución de espesor del perfil
analizado y al ángulo de ataque correspondiente a cero sustentación o ángulo de ataque cero
para el perfil simétrico.

Para el perfil NACA 641-212 este fenómeno ocurre a 0.4c o 40% de la longitud de cuerda.
Adicionalmente, se debe tener un cuenta el tipo de línea de camber o "NACA Mean Line"
cuando se evalúa la nomenclatura. Para este caso en particular, el tipo de línea de camber a
= 0.4 indica que se mantiene una presión uniforme sobre el perfil alar hasta el 40% de la
longitud de la cuerda. Usualmente, el tipo de línea de camber o "a" tiene valores iguales o
mayores al punto de ubicación de mínima presión.

El tercer dígito da el coeficiente de sustentación de diseño en décimas de cuerda. En este caso

el coeficiente de sustentación de diseño es de 0.2. Finalmente, los últimos dos dígitos otorgan
el máximo grosor en porcentaje de la cuerda, que para el NACA 641-212 sería del 12%.

La importancia de los perfiles alares de la familia NACA Serie 6, se deriva de la habilidad

para minimizar la resistencia, especialmente el coeficiente de resistencia con variación del
coeficiente de sustentación.

Usualmente, estos perfiles exhiben en el gráfico de la curva de resistencia un "Drag Bucket"

o una precipitación en el coeficiente de resistencia Cd en el rango de los ángulos de ataque
bajos. Esta precipitación o drag bucket es una característica típica de los perfiles alares
NACA serie 6.

El subíndice de este perfil indica que tiene un rango operativo de Cl= ±0.1 del coeficiente de
sustentación de diseño. En comparación con un perfil de la familia de 4 dígitos, el NACA
serie 6 tiene una reducción de 33% en el coeficiente de resistencia mínima (Abbot &
Doenhoff, 1949). La Figura 9 ilustra el típico gráfico del resultado de las investigaciones y
pruebas en túneles de viento de la NACA en cuanto al perfil de flujo laminar serie 6 NACA
641- 212.

Adicionalmente, se ilustra la variación del coeficiente de resistencia del perfil con variación
del coeficiente de sustentación o en otras palabras, variación de ángulo de ataque. Es visible
el "drag bucket" para números de Reynolds de 3.0x10^6 hasta 9.0x10^6 y coeficiente de
sustentación entre 0.1 y 0.3. Anderson (2001) menciona que la primera aeronave militar en
utilizar este tipo de perfil alar fue el North American P-51 Mustang versión H
Figura . Variación de Coeficiente de Resistencia con Coeficiente de Sustentación Perfil
NACA 64-212. Fuente: Abbot & Doenhoff. (1949). Theory of Wing Sections, p. 571. New
York: Dover Publications.

Teorìa General de Perfiles Delgados

Esta teoría maneja 3 hipótesis básicas:

La viscosidad casi no afecta el flujo de aire sobre el perfil que no ha entrado en

pérdida
La operación de perfil se realiza a bajos ángulos de ataque
Las fuerzas resultantes por la presión son ligeramente influenciadas por el espesor del
perfil
Para que un perfil se considere delgado, este debe tener una curvatura máxima de 2% de la
cuerda y un espesor no mayor a 12% de la cuerda.

Como sabemos la diferencia de velocidades del aire entre el intradós y el extradós del
perfil produce la sustentación, la diferencia de presiones que resulta, puede ser
reemplazada infinitesimalmente por un conjunto de vórtices distribuidos a lo largo de la
línea media del perfil superpuesto sobre el flujo de aire, cuya fuerza rotacional será:

c
Γ=∫0 γ ( s) ds

Donde γ (s ) es la distribución de vórtices sobre la línea media y ds es un elemento

infinitesimal de dicha distribución.

Figura 1. Esquema de la distribución de los vórtices a lo largo de la línea media

El campo de velocidades en la cuerda media entonces resulta de la suma de las velocidades
de los vórtices y de la velocidad del viento. La suma de estas velocidades debe resultar una
velocidad tangente a la línea media en cualquier punto y por lo tanto la velocidad
perpendicular a dicha línea en cualquier punto debe ser cero. Adicionalmente la unión de
los flujos en la parte posterior debe ser suave para que se cumpla el teorema de Kutta.

Existe en la línea media entonces un punto P genérico donde los vórtices de un elemento
diferencial ds influyen generando una velocidad perpendicular a un vector r generado entre
el elemento y el punto.
γdsCos(β)
dV = − 2πr

El componente perpendicular en dicho punto de la velocidad de flujo libre, el cual es

opuesto a la velocidad previamente obtenida será entonces:
𝑈 = 𝑈∞ 𝑠𝑒𝑛(𝛼 − 𝛽)
Donde:
𝑑𝑦
𝛽 = arctan (𝑑𝑥 )

Descripción del Método de Vórtices Discretos en Flujos Bidimensionales

Este método se basa en crear pequeños vórtices ubicados uno a uno en líneas rectas que
sigan a la línea de curvatura media (obviamente mientras más elementos la curva se
parecerá más al conjunto de elementos) que sumados equivalen a la sustentación total del
perfil, al conocer todos los datos (sean de geometría, velocidad del viento y ángulo de
ataque) las únicas incógnitas serían las fuerzas de sustentación de los vórtices.

De manera similar a la del modelo discreto del perfil, la ecuación se crea en base a la
condición de capa límite que dice que el fluido no puede atravesar un perfil, por lo que en
un determinado punto de control se analiza esta condición relacionando la componente
normal al elemento del flujo de aire y la suma sobre dicho elemento de control de la
influencia de cada vórtice creado (obviamente también la componente normal).
A continuación, detallamos el procedimiento para llegar a la solución
-Creación del perfil: Generalmente se utilizan perfiles conocidos normados, una de las
normas son la norma NACA que se emplean para distintos propósitos (crear perfiles
nuevos, etc.) se pueden utilizar estos perfiles como referencia.

-Ubicación de la línea media: Se obtiene de la unión de los puntos equidistantes del

intradós y extradós referenciados a la perpendicular de la cuerda.

-Discretización de la línea media: Los elementos de discretización son líneas rectas de

igual longitud (esto es un procedimiento flexible ya que si se quiere conocer más sobre
algún punto en especial del perfil se pueden utilizar elementos más pequeños en esa zona
(por ejemplo, en las zonas donde se predice el inicio de la separación de capa límite) y más
grandes en las zonas que no se quieren analizar con detenimiento.

-Ubicación de los puntos de vórtices: Los vórtices se ubicarán a h/4 del inicio del
elemento (siendo h la longitud total del elemento y el inicio del elemento el punto más
cercano al borde de ataque).

-Ubicación de los puntos de control: Los puntos de control se ubicarán a 3h/4 del inicio
del elemento.

-Relación de los vórtices: Cada vórtice tendrá influencia sobre cada elemento explicado
matemáticamente de la siguiente manera:
𝑢 Γ 0 1 𝑥𝑖 − 𝑥𝑗
𝑐𝑖𝑗 = [ ] = 2𝜋𝑟𝑖 2 [ ] [𝑦 − 𝑦 ]
𝑤 𝑖𝑗 −1 0 𝑖 𝑗
donde:
𝑟𝑖𝑗2 = (𝑥𝑖 − 𝑥𝑗 )2 + (𝑦𝑖 − 𝑦𝑗 )2
i=1,2,3,…,m son los diferentes puntos de control de cada elemento
j=1,2,3,…, n son los diferentes vòrtices creados
obviamente m=n
Esta ecuación determina la componente normal de la influencia de cada vórtice sobre cada
elemento de control.
Condicòn de capa lìmite: Como ya explicamos, la condición de capa establece el no flujo a
través del perfil. Esa condición de no flujo se describe matemáticamente al igualar
escalarmente la influencia de todos lo vòrtices sobre un punto de control con el flujo de aire
de la siguiente manera:
𝑎11 . Γ1 + 𝑎12 . Γ2 + 𝑎13 . Γ3 + ⋯ + 𝑎1𝑁 . Γ𝑁 + (𝑈∞ cos(𝛼) , 𝑈∞ sen(𝛼)). 𝑛̂1 = 0
Donde: 𝑎𝑖𝑗 es:

𝑐𝑖𝑗
𝑎𝑖𝑗 = ( ) . 𝑛̂𝑖
Γ𝑗

y 𝑛̂𝑖 es el vector normal al elemento.

[(𝑦 −𝑦 )(𝑥 −𝑥 )]
𝑛̂ = [(𝑦𝐴−𝑦𝐵 )(𝑥𝐴−𝑥𝐵 )] Donde A y B son los puntos inicial y final del elemento.
𝐴 𝐵 𝐴 𝐵

Esta ecuación se da para cada elemento juntando y ordenando todas las ecuaciones lineales
obtenemos una ecuación matricial similar a esta:
𝑈1 𝑎1,1 𝑎1,2 𝑎1,3 𝑎1,4 Γ1
𝑈 𝑎2,1 𝑎2,2 𝑎2,3 𝑎2,4 Γ2
[ 2 ] = [𝑎 𝑎3,2 𝑎3,3 𝑎3,4 ] [Γ3 ]
𝑈3 3,1
𝑈4 𝑎 4,1 𝑎4,2 𝑎4,3 𝑎4,4 Γ4

Donde:
𝑈𝑖 = (𝑈∞ cos(𝛼) , 𝑈∞ sen(𝛼)). 𝑛̂𝑖
Una vez establecida la ecuación matricial y aplicando teoría de matrices obtendremos la
fuerza Γ de los vòrtices, con dichos resultados y la ecuación de Kutta se aplica el método:
∆𝐿𝑗 = 𝜌. 𝑈∞ . Γ𝑗

La sustentaciòn total serà la suma de los ∆𝐿′𝑠. En base a esos resultados ya se puede
obtener además los momentos de giro en el perfil y por lo tanto determinar la estabilidad y
predecir el comportamiento del mismo
A continuación, se mostrarà el algoritmo general para n puntos.
ℎ
𝑥(𝑖) = 4 + ℎ ∗ (𝑖 − 1)
x: representa la distancia al vòrtice
3ℎ
𝑦(𝑖) = + ℎ ∗ (𝑗 − 1)
4

y: representa la distancia al punto de referencia

Γ
𝑉(𝑠𝑒𝑛𝛼 − 𝛽) = ∑ 2πh∗(0.5+(j−1))

La anterior ecuación determina la componente normal de la influencia de cada vòrtice sobre

cada elemento de control.
Γ
Si ∑ Vh = 𝑌, entonces se tiene:
𝑌∗1
2𝜋(𝑠𝑒𝑛𝛼 − 𝛽) = 0.5+(𝑗−𝑖)

𝑓 =𝑌∗𝐴
𝑌
Donde 𝐴𝑖𝑗 = 0.5+(𝑗−𝑖) y 𝑓 = 2𝜋(𝑠𝑒𝑛𝛼 − 𝛽𝑃𝑅 )

Finalmente
𝑌 = 𝑓 ∗ 𝐴−1
Una vez establecida la ecuación matricial y aplicando teoría de matrices obtendremos la
fuerza 𝛾 de los vòrtices, con dichos resultados y la ecuación de Kutta aplicada al método:
∆𝐿𝑗 = 𝜌. 𝑈∞ . Γ𝑗

Finalmente sumar todos los ∆𝐿′𝑠 para obtener la sustentaciòn.

DESCRIPCIÒN DEL PROBLEMA
Se tomarà un perfil naca delgado al azar, debido a que encontramos el diseño en la
literatura, tenemos el naca 2310, que es similar al perfil Robin 23010
Es decir que nuestro tipo perfil NACA 4 dígitos finalmente será:

NACA 2 3 1 0
En base a este perfil obtenemos los puntos del extradós e intradós tomando como referencia
el eje x como la cuerda del perfil y el eje y como el primer punto del borde de ataque en el
eje x.
x= [1.000000 0.975528 0.904508 0.793893 0.654508 0.500000 0.345492 0.206107
0.095492 0.024472 0.000000 0.024472 0.095492 0.206107 0.345492 0.500000
0.654508 0.793893 0.904508 0.975528 1.000000];

y= [0.000000 0.004292 0.015977 0.032048 0.048776 0.062419 0.069546 0.066134

0.049080 0.024708 0.000000 -0.018448 -0.027668 -0.030052 -0.029715 -0.025684
-0.019035 -0.011960 -0.005808 -0.001544 0.000000];

Y representamos los puntos de la cuerda media mediante el promedio de las coordenadas

“y” del intradós y el extradós.
Representamos estos puntos en la siguiente tabla:

X y(extradòs) y(intradòs) YLm

0.0000000 0.0000000 0.0000000 0.0000000
0.0244720 0.0247080 -0.0184480 0.0031300
0.0954920 0.0490800 -0.0276680 0.0107060
0.2061070 0.0661340 -0.0300520 0.0180410
0.3454920 0.0695430 -0.0297150 0.0199140
0.5000000 0.0624190 -0.0256840 0.0183675
0.6545080 0.0487760 -0.0190350 0.0148705
0.7938930 0.0320480 -0.0119600 0.0100440
0.9045080 0.0159770 -0.0058080 0.0050845
0.9755280 0.0042920 -0.0015440 0.0013740
1.0000000 0.0000000 0.0000000 0.0000000

Observamos estos puntos en el siguiente gràfico:

 Naca 2310
0.0800000
0.0600000
0.0400000
0.0200000
y

0.0000000
0.0000000
-0.0200000 0.2000000 0.4000000 0.6000000 0.8000000 1.0000000 1.2000000

-0.0400000
x

y(extradòs) y(intradòs) YLm

Nuestro objeto de estudio es la línea media. Las coordenadas de la línea

media ajustamos a una curva mediante el método de mínimos cuadrados para una
curva de grado 3 de la siguiente manera.

𝑓(𝑥) = 𝑎0 + 𝑎1 𝑥 + 𝑎2 𝑥 2 + 𝑎3 𝑥 3
Creamos nuestro sistema de ecuaciones y hallamos los valores de cada sumatoria:

Tenemos la siguiente tabla:

Finalmente reemplazamos los dato en la matriz anterior y obtenemos:

Resolviendo la matriz tenemos los siguientes valores:

𝑎0 = 0.0006
𝑎1 = 0.1153
𝑎2 = −0.1947
𝑎3 = 0.0794
Finalmente, nuestra curva nos quedarà como:
𝑓(𝑥) = 0.0006 + 0.1153𝑥 − 0.1947𝑥 2 + 0.0794𝑥 3
Como resultado tenemos la siguiente tabla:

x f(x)
0.0000000 0.0006
0.0244720 0.003306
0.0954920 0.009904
0.2061070 0.016788
0.3454920 0.020469
0.5000000 0.0195
0.6545080 0.014921
0.7938930 0.009152
0.9045080 0.004356
0.9755280 0.001503
1.0000000 0.0006
Y el siguiente gràfico:
 Distribuciòn de Presiones
0.025

0.02
y = 0.0794x3 - 0.1947x2 + 0.1153x + 0.0006
R² = 1
0.015
f(x)

0.01

0.005

0
0.0000000 0.2000000 0.4000000 0.6000000 0.8000000 1.0000000 1.2000000
x

Algoritmo en programa Matlab