Método Paneles

Estudio aerodinmico de un perfil NACA mediante el mtodo
de paneles
Adrin Ortas Delgado
NExp: 112087
Grupo: VA-2
1.
Objetivo
El objetivo de este trabajo es aplicar el mtodo de paneles para estimar las caractersticas aerodinmicas de un perfil NACA, en este caso, del NACA 2712. Asimismo, los resultados obtenidos son
comparados posteriormente con otro mtodo numrico (X-Foil).
2.
Dibujo del perfil

Teniendo en cuenta que la lnea media del perfil est determinada por las parbolas:
8
f
2
>
0 x xf
>
< x2f 2xf x x
yc (x) =
>
>
: f 2 (1 2xf ) + 2xf x x2
xf x c
(1 xf )
y particularizando para el perfil que se estudia, su representacin queda de la siguiente forma:
3.
Mtodo de paneles aplicado a perfiles delgados
En este estudio, se ha considerado utilizar 25 paneles distribuidos de forma equidistante, teniendo

la cuerda un valor unitario, de forma que cada panel mide 0,04, hallndose los puntos de control en la
mitad de cada uno de los paneles. La velocidad de la corriente incidente no perturbada tambin tiene
un valor unitario.
3.1.
Distribucin de sustentacin
A continuacin se representa la distribucin de sustentacin en cada panel para los ngulos de

ataque pedidos. Adems, se incluyen los coeficientes CL , CM 4c , CM BA , xcp , representados y tabulados.
Todos los ngulos se dan en grados.
ngulo de ataque
-5
0
5
CL
-0,23
0,32
0,86
CM 4c
-0,09
-0,09
-0,09
CM 0
-0,03
-0,17
-0,30
xcp
-0,14
0,53
0,35
3.2.
Resto de coeficientes aerodinmicos
En la siguiente tabla se recogen otros datos aerodinmicos que caracterizan el perfil en estudio:
CL0
0,32
CL
0,11
CL =0
-2,92
xCA
0,24
CM ca
-0,09
4.
Comparacin con otros mtodos
En las siguientes tablas se muestran los datos obtenidos mediante otro mtodo de clculo numrico,
en este caso se trata del software X-Foil.
En primer lugar se considera el caso sin viscosidad:
-5
0
5
Valores obtenidos
CL
CM 0 CM 4c
-0,24 -0,03 -0,09
0,37 -0,19 -0,10
0,97 -0,34 -0,10
xcp
-0,12
0,51
0,36
Error cometido ( %)
CL
CM 0 CM 4c
xcp
3,36
-8,77
1,14 -12,74
12,64 6,38
8,42
3,07
10,93 11,91 13,73
0,95
En segundo lugar se considera el caso viscoso con un Re = 9106
-5
0
5
CL
-0,23
0,33
0,89
Valores
CD
0,01
0,01
0,01
obtenidos
CM 0 CM 4c
-0,03 -0,09
-0,17 -0,09
-0,31 -0,09
xcp
-0,13
0,51
0,35
Error cometido ( %)
CL
CM 0 CM 4c
xcp
-0,88 -1,31 0,68
0,59
4,02 -2,06 0,91 -2,27
2,60
1,38 -2,33 -1,39
En esta situacin, se calculan los coeficientes de la polar parablica del perfil, as como el CLmax y
el CLmn :
CLmax
CLmn
2,03
-1,64
ngulo de ataque
20
-20
CD = CD0 + aCL + bCL2

CD0
a
b
0,004 -0,005
0,014
5.
Conclusiones
A la vista de la comparacin de los resultados obtenidos por medio del mtodo de paneles con los
hallados por el programa X-Foil, puede decirse que el estudio ha sido satisfactorio, pues el error en el
caso viscoso es menor del 5 %. Por otro lado, el error aumenta en el caso no viscoso hasta un 14 %,
valor que sigue siendo aceptable en primera aproximacin.
De las curvas de las distribuciones de sustentacin a lo largo de la cuerda del perfil, se aprecia una
asntota vertical cerca del borde de ataque. Esto es debido al uso de la Teora Potencial Linealizada,
que presenta singularidades en los puntos de remanso adems de las que aparecen en las esquinas en
una teora potencial. Para el caso particular del ngulo de ataque ideal, esta singularidad del borde de
ataque desaparece, tal y como se puede comprobar en el grfico.
6.
Anexo
A continuacin se recoge el cdigo MATLAB utilizado para la realizacin de este trabajo, del que
se obtienen todos los datos pedidos. Por comodidad y simplicidad en el cdigo, los vectores b se han
calculado en Excel, aprovechando las ecuaciones para determinar la lnea media. Posteriormente, las
coordenadas de los vectores se han aadido al cdigo MATLAB.
%Construyo m a t r i z B
n=25;
B=zeros ( n ) ;
f o r i =1:n
f o r j =1:n
B( i , j )= ( i
end
end
0.5 j ) log ( abs ( ( i
0.5 j ) / ( i
0.5 ( j
1))))
1;
%Construyo m a t r i z C
C=zeros ( n ) ;
f o r i =1:n
f o r j =1:n
i f i==1
C( j , i )=0;
else
C( j , i )=(( i
end
end
end
1) 0.5 ( j ) ) log ( abs ( ( ( i
1) 0.5 ( j ) ) / ( ( i
1) 0.5 ( j
1 ) ) ) ) +1;
A=B+C
%V e c t o r e s b , i m p o r t a d o s d e s d e e x c e l
%Vector b a l f a = 5
b1 = [ 0 . 8 9 6 7 3 4 6 5 7 , 0 . 8 7 6 2 7 7 5 3 5 , 0 . 8 5 5 8 1 3 4 5 1 , 0 . 8 3 5 3 4 2 8 3 5 , 0 . 8 1 4 8 6 6 1 2 , 0 . 7 9 4 3 8 3 7 3 8 ,
0.773896123 ,0.753403707 ,0.732906925 ,0.712406213 ,0.691902005 ,0.671394738 ,
0.650884849 ,0.630372773 ,0.609858947 ,0.589343809 ,0.568827797 ,0.53691329 ,
0.436622042 ,0.325003276 ,0.213525322 ,0.10225792 , 0.008729979 , 0.119370467 ,
0.229596924];
b1=b1 ;
%Vector b a l f a =0
b2 = [ 0 . 3 4 8 4 2 3 3 1 1 , 0 . 3 2 7 9 6 6 1 8 9 , 0 . 3 0 7 5 0 2 1 0 4 , 0 . 2 8 7 0 3 1 4 8 9 , 0 . 2 6 6 5 5 4 7 7 4 , 0 . 2 4 6 0 7 2 3 9 2 ,
0.225584776 ,0.20509236 ,0.184595579 ,0.164094866 ,0.143590659 ,0.123083392 ,
0.102573502 ,0.082061426 ,0.061547601 ,0.041032463 ,0.02051645 , 0.011398056 ,
0.111689305 , 0.22330807 , 0.334786024 , 0.446053426 , 0.557041325 , 0.667681814 ,
0.777908271];
b2=b2 ;
%Vector b a l f a =5
b3 =[ 0.199888035 , 0.220345158 , 0.240809242 , 0.261279857 , 0.281756572 ,
0.302238954 , 0.32272657 , 0.343218986 , 0.363715768 , 0.38421648 , 0.404720687 ,
0.425227954 , 0.445737844 , 0.46624992 , 0.486763746 , 0.507278883 , 0.527794896 ,
0.559709402 , 0.660000651 , 0.771619417 , 0.883097371 , 0.994364772 , 1.105352671 ,
1.2159931 , 1.326219617];
b3=b3 ;
%Vector b a l f a i d e a l
b4 = [ 0 . 4 1 1 2 5 5 1 6 3 , 0 . 3 9 0 7 9 8 0 4 1 , 0 . 3 7 0 3 3 3 9 5 6 , 0 . 3 4 9 8 6 3 3 4 1 , 0 . 3 2 9 3 8 6 6 2 6 , 0 . 3 0 8 9 0 4 2 4 4 ,
0.288416628 ,0.267924212 ,0.247427431 ,0.226926718 ,0.206422511 ,0.185915244 ,
0.165405354 ,0.144893278 ,0.124379453 ,0.103864315 ,0.083348302 ,0.051433796 ,
0.048857453 , 0.160476218 , 0.271954172 , 0.383221574 , 0.494209473 , 0.604849962 ,
0.715076419];
b4=b4 ;
%D i s t r i b u c i n de v o r t i c i d a d e s
gamma1=A\b1 ;
gamma2=A\b2 ;
gamma3=A\b3 ;
gamma4=A\b4 ;
%C o e f i c i e n t e s de s u s t e n t a c i n
c l 1 =(2sum(gamma1)
gamma1 ( 1 ) ) ( 1 / n )
gamma2 ( 1 ) ) ( 1 / n )
gamma3 ( 1 ) ) ( 1 / n )
gamma4 ( 1 ) ) ( 1 / n )
%C o e f i c i e n t e s de momentos l o c a l e s
f o r i =1:n
i f i==n
c l 1 ( i )=gamma1( n ) ;
c l 2 ( i )=gamma2( n ) ;
c l 3 ( i )=gamma3( n ) ;
c l 4 ( i )=gamma4( n ) ;
else
c l 1 ( i )= gamma1( i ) + gamma1( i +1);
end
end
c l 1=c l 1
c l 2=c l 2
c l 3=c l 3
c l 4=c l 4
f o r i =1:n
cm1 ( i )=
cm2 ( i )=
cm3 ( i )=
cm4 ( i )=
end
cl1 ( i
cl2 ( i
cl3 ( i
cl4 ( i
)((
)((
)((
)((
i /n)
i /n)
i /n)
i /n)
1/(2n ) ) ;
1/(2n ) ) ;
1/(2n ) ) ;
1/(2n ) ) ;
%C o e f i c i e n t e s de momentos g l o b a l e s
cM1=sum( cm1 ) ( 1 / n )
cM2=sum( cm2 ) ( 1 / n )
cM3=sum( cm3 ) ( 1 / n )
cM4=sum( cm4 ) ( 1 / n )
Referencias
[1] Meseguer Ruiz, J. Sanz Andrs, A. (2011). Aerodinmica Bsica. Madrid, Espaa: Ibergarceta
Publicaciones.
[2] Ganda Agera, F. Gonzalo de Grado, J. Margot, X. Meseguer Ruiz, J. (2013). Fundamentos de
los mtodos numricos en aerodinmica. Madrid, Espaa: Ibergarceta Publicaciones.

Método Paneles

Cargado por

Copyright:

Formatos disponibles

Método Paneles

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Método Paneles

Cargado por

Copyright:

Formatos disponibles

Estudio aerodinmico de un perfil NACA mediante el mtodo

Dibujo del perfil

y particularizando para el perfil que se estudia, su representacin queda de la siguiente forma:

Mtodo de paneles aplicado a perfiles delgados

En este estudio, se ha considerado utilizar 25 paneles distribuidos de forma equidistante, teniendo

A continuacin se representa la distribucin de sustentacin en cada panel para los ngulos de

Resto de coeficientes aerodinmicos

Comparacin con otros mtodos

En segundo lugar se considera el caso viscoso con un Re = 9106

CD = CD0 + aCL + bCL2

0.5 j ) log ( abs ( ( i

1) 0.5 ( j ) ) log ( abs ( ( ( i

También podría gustarte