Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

TEMA 5 Flujo Potencial

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 49

TEMA5:FLUJOPOTENCIAL

ECUACINDE BERNOUILLI
El estudio del flujo sin friccin a travs de un tubo de corriente infinitesimal, como muestra la
figura 4.1a, proporciona una relacin muy utilizada entre la presin, la velocidad y la altura, que
se denomina ecuacin de Bernouilli.

Fig. 4.1 Ecuacin de Bernouilli para flujos sin friccin a lo largo de una lnea de corriente: (a) fuerzas y
flujos, (b) resultante de las fuerzas de presin despus de retar una presin p uniforme.

El volumen de control analizado coincide con un tubo de corriente infinitesimal de rea variable
A(s) y longitud Js, donde s representa la direccin de la lnea de corriente. Las propiedades
(p, p, I) pueden variar con s y con el tiempo, pero consideraremos que son constantes en la
seccin transversal, la cual es suficientemente pequea. El tubo de corriente est inclinado un
ngulo arbitrario 0, de forma que la variacin de altura entre las secciones es Jz = Js scn0. La
figura muestra una friccin inevitable en las paredes del tubo, que despreciamos, adems en el
lmite JA u el tubo de corriente coincide con la lnea de corriente.
La ecuacin de Bernouilli se aplica a lo largo de una lnea de corriente en un flujo no viscoso.
Aplicamos la conservacin de la masa en el volumen de control
J
Jt
__ pJV
vC
_ +m
suI
-m
cnt
= u =
op
ot
JV +Jm
m = pAI y JV = A Js Jm = J(pAI) = -
op
ot
AJs
Escribimos ahora la ecuacin de conservacin de la cantidad de movimiento en la direccin de
la corriente
JF
s
=
J
Jt
__ IpJ
vC
V_ +(m I)
suI
-(m I)
cnt.
=
o
ot
(pI)A Js +J(m I)
Si despreciamos los esfuerzos tangenciales en las paredes (flujo sin friccin), los trminos de
fuerza se deben slo a la presin y la gravedad. La fuerza de gravedad en la direccin de la
corriente es la componente del peso del fluido contenido en el volumen de control.
JF
s,gu
= -Jw scn0 = -gpA Js scn0 = -gpAJz
La fuerza resultante de presin es JF
s,pcs
= -AJp
Sustituimos estos dos trminos en la ecuacin de conservacin de la cantidad de movimiento
JF
s
= -gpAJz -AJp =
o
ot
(pI)A Js +J(m I)
-gpAJz -AJp =
o
ot
(p)IA Js +
o
ot
(I)pA Js +m J(I) +IJ(m )
De la ecuacin de conservacin de masa utilizamos Jm = -
p
t
AJs y m = pAI .
Sustituimos en la ecuacin de conservacin de la cantidad de movimiento
-gpAJz -AJp =
o
ot
(p)IA Js +
o
ot
(I)pA Js +pAIJ(I) -I
op
ot
AJs
-gpAJz -AJp =
o
ot
(I)pA Js +pAIJ(I)
oI
ot
Js +
Jp
p
+IJI +gJz = u
Hemos llegado as a la ecuacin de Bernouilli para flujo no estacionario sin friccin a lo largo de
una lnea de corriente. Esta ecuacin se puede integrar entre dos puntos 1 y 2 a lo largo de una
lnea de corriente
_
oI
ot
Js
2
1
+_
Jp
p
2
1
+
1
2
(I
2
2
-I
1
2
) +g(z
2
-z
1
) = u
Si consideramos el caso estacionario [
v
t
= u e incompresible, densidad constante, la
ecuacin de Bernouilli queda de la forma
p
2
-p
1
p
+
1
2
(I
2
2
-I
1
2
) +g(z
2
-z
1
) = u
p
1
p
+
1
2
I
1
2
+gz
1
=
p
2
p
+
1
2
I
2
2
+gz
2
= ctc.
Esta es pues la ecuacin de Bernouilli para un flujo estacionario incompresible y sin friccin a lo
largo de una lnea de corriente. Esta ecuacin est relacionada con la ecuacin de la
conservacin de la energa en rgimen estacionario. Recordamos que la ecuacin de
conservacin de la energa para un flujo en un tubo de corriente con una entrada y una salida
era de la forma
b
`
1
+
o
2
I
1
2
+gz
1
= [b
`
2
+
o
2
I
2
2
+gz
2
-q +w
s
+w


p
1
p
+
o
2
I
1
2
+gz
1
= _
p
2
p
+
o
2
I
2
2
+gz
2
] +(u
2
-u
1
) -q +w
s
+w



Esta relacin es mucho ms general que la ecuacin de Bernouilli, ya que permite tener en
cuenta la friccin, la transferencia de calor, el trabajo mecnico y el trabajo viscoso, pero
implica una simplificacin importante; asumir constantes las propiedades del flujo en la entrada
y salida del circuito.
La ecuacin de Bernouilli en aplicable en las siguientes condiciones:
- flujo estacionario;
- incompresible, aceptable si el nmero de Mach es inferior a 0.3;
- sin friccin;
- flujo a lo largo de una lnea de corriente;
- sin intercambio de trabajo ni transferencia de calor entre los dos puntos considerados;
Es importante destacar que lneas de corriente distintas pueden tener diferentes constantes de
Bernouilli w
0
=
p
p
+
v
2
2
+gz, dependiendo de las condiciones del flujo.
En la obtencin de la ecuacin de Bernouilli no se consideran tampoco transferencia de calor o
trabajo. La razn bsica de estas restricciones es que en fluidos reales los intercambios de
calor y trabajo estn ligados a efectos de friccin, lo que invalida la hiptesis de flujo sin
friccin. Estos efectos termodinmicos son tenidos en cuenta en la ecuacin de la energa en
un flujo estacionario.
A continuacin se muestran algunas limitaciones prcticas del uso de la ecuacin de Bernouilli.

Fig. 4.2 Ilustracin de las zonas de validez o no validez de la ecuacin de Bernouilli: (a) modelo en un
tnel aerodinmico, (b) hlice, (c) chimenea

Una interpretacin visual muy til de la ecuacin de Bernouilli se obtiene representando dos
lneas de flujo. La lnea de nivel de energa (LNE), tambin conocida como lnea de cargas o
alturas totales que muestra la altura de la constante de Bernouilli b
0
= z +
p
gp
+
v
2
2g
. En un flujo
sin friccin y sin aplicacin de calor o trabajo, la lnea es una lnea de nivel constante. La lnea
de altura motriz (LAM), tambin conocida como lnea de cargas o alturas piezomtricas, indica
el nivel correspondiente a la altura geomtrica ms la de presin z +
p
gp
, esto es la LNE menos
la altura de velocidad
v
2
2g
. La LAM es la altura a la que subira el lquido en un tubo piezomtrico
incorporado al flujo. En el flujo en un canal abierto, la LAM es la superficie libre del agua.

Fig. 4.3 Lneas de nivel de energa y lnea de altura motriz para flujo sin friccin en un conducto
La figura anterior muestra las lneas LNE y LAM para un flujo sin friccin en un conducto. Los
tubos piezomtricos de las secciones 1 y 2 miden la carga de la presin esttica z +
p
gp
, y por
tanto la LAM. Los tubos de pitot de presin de remanso miden la altura total z +
p
gp
+
v
2
gp
, que
corresponde a la LNE. En este caso particular la LEN es constante y la LAN asciende lo que
indica que la velocidad est disminuyendo. En condiciones ms generales la LEN disminuir
lentamente debido a las prdidas por friccin, y descender bruscamente por prdidas
localizadas (una vlvula u obstruccin) o debido a la extraccin de trabajo (en una turbina). La
LEN slo puede ascender si se comunica trabajo (como en una bomba o una hlice). La LAM
sigue el comportamiento de la LEN respecto a prdidas y trabajo motor y asciende o desciende
al disminuir o aumentar la velocidad, respectivamente.
INTRODUCCINA LATEORA DEPERFILES
Una de las aplicaciones bsicas de la mecnica de fluidos se refiere al estudio de las denominadas
superficiessustentadoras.Tpicamente,unasuperficiesustentadoraestconstituidoporuncuerpofino
currentiforme (perfil) que se enfrenta al fluido con un pequeo ngulo de ataque o de incidencia
1
.
Ejemplos de superficies sustentadoras son el ala de un avin, el timn y las palas de la hlice de un
buque,laquilladeunveleroounavelaenceida.
La caracterstica principal de este tipo de cuerpos es que su enfrentamiento a una corriente de fluido
provoca la aparicin de una fuerza, que puede ser descompuesta en dos componentes. Una en la
direccinperpendicularalflujo(direccindeavancedelperfil)quesedenominasustentacin(lift)yotra
enladireccindelflujo(resistenciaodrag).
EnlaFigura 1semuestraelesquemadefuncionamientodeunperfildecuerdaC,endondedebidoa
un flujo de velocidad U, con un ngulo de ataque , se desarrolla una fuerza F que se puede
descomponerenlasustentacinLyenlaresistenciaD(lalongitudlsedenominacuerdadelperfil).

Figura1.Esquemadefuncionamientodeunperfil.
Esposibleexplicardeunamaneraintuitivalaaparicindelafuerzadesustentacinenunperfil,conuna
ciertainclinacin,enfrentadoaunacorrienteuniforme,apartirdelestudiodelaasimetradelflujoque
aparece. Dada su forma y teniendo en cuenta la continuidad (conservacin de masa) del fluido, las
partculas de la parte superior del perfil (cara de succin) deben ir ms rpido que las partculas de la
parte inferior donde su recorrido es menor (cara de presin). Esta aceleracin se traduce en que hay
menorpresinenlapartesuperiorqueenlainferioryestadiferenciadepresioneseslaqueprovocala
sustentacin.Lajustificacindeestefenmenopuedehacerseconmayorsentidofsicobasndoseenla
ecuacindeBernoulli.
En efecto, aunque la ecuacin de Bernoulli es slo aplicable a los fluidos ideales, resulta muy til para
analizar algunos problemas en los que el fenmeno de turbulencia es poco significativo o est muy
localizado.SinosfijamosenlaFigura 2,podemosobservarelflujotpicoalrededordeunperfilconun
ngulodeataquepequeo,obtenidomedianteeltrazadodelneasdecorrienteenunexperimento.Es

1
La mayora de los perfiles utilizados en ingeniera naval son simtricos,aunque su inclinacin(ngulo
deataque)lohaceaparecernosimtricofrenteaunflujouniforme.
evidentequeenestoscasos,elflujopotencialrepresentaconbastanteaproximacinlarealidad
2
.Algo
similarocurreporejemploenelflujoentoberasojets.

Figura2.Ensayoexperimentaldelflujoalrededordeunperfil
En casos como este, es posible calcular la distribucin de presin sobre un perfil, una vez conocida la
velocidadyviceversa
3
,atravsdelaError!Noseencuentraelorigendelareferencia..
Como ltima conclusin, baste confirmar que la Error! No se encuentra el origen de la referencia.,
justificalaaparicindelfenmenodesustentacinenlosperfiles,puesvinculalaaceleracindelfluido
conunadepresin(amayorvelocidad,menorpresinyviceversa).
FLUJO ALREDEDOR DEUN PERFIL
Elflujohabitualdeunfluidoquecirculaalrededordeunperfil,puededividirseendiversasregiones,de
manera similar a como se present en el anlisis del flujo alrededor de una carena (ver Error! No se
encuentraelorigendelareferencia.).
Engeneral,existirunazonaexteriorenlaquelosefectosdelaviscosidadsondespreciablesyunazona
cercanaalperfilyqueseextenderaguasabajodelmismoenlaestosefectossernimportantes,taly
comosemuestraenlaFigura 3.Enelcasodelosperfiles,lazonaocupadaporlacapalmite(lazona
de transicin entre el flujo cerca del perfil y el flujo exterior), aunque dependiente del nmero de
Reynolds,esenlaprcticamuydelgada.Eshabitualquesuvalorenelbordedesalidadelflujoseadel
ordendel35%delacuerdadelperfil.Sinembargo,lasfuerzasqueactansobreelperfildependende

2
Noocurreassiseaumentasuficientementeelngulodeataque.
3
LaconstantequeapareceenlaecuacindeBernoullipuedecalcularsesiavanzamosaguasarribapor
unalneadecorriente,hastaencontrarelflujonoperturbado.Comoesaconstanteeslamismapara
cadalneadecorriente,ycomoenesazonalapresineslahidrosttica,podemosescribirlaigualdad:
2
1
2
u p gz + =
Dondegeslaaceleracindelagravedadyzlaprofundidaddelalneadecorrienteaguasarriba.
manera muy importante de cmo se desarrolle esta capa a lo largo del perfil. Por otra parte, este
desarrollodependedeladistribucindepresinsobreelperfil,yenparticularsobrelosgradientesde
presin a lo largo de la cuerda. Se suele hablar de gradientes de presin favorables, cuando la presin
disminuye a lo largo de la cuerda, facilitando la aceleracin del flujo. A la inversa se denominan
gradientesdepresindesfavorablesalosqueseoponenaladireccindeavancedelflujo.
Esimportanteremarcarqueladistribucindepresindependedemaneramuyimportantedecmose
desarrollaelflujoexterioralacapalmite.

Sipasamosadescribirelflujohabitualenlacapalmite,podemosdecirque,cercadelbordedeentrada
delflujoenelperfil,lacapalmitesemuestraordenada,siguiendounflujolaminar,perorpidamente
se alcanza una zona de transicin, a partir de la cual el movimiento del flujo se vuelve errtico, dando
lugar a un flujo turbulento. En una zona importante del perfil, este flujo turbulento contina, en
promedio,siguiendolasuperficiedelperfil.
Los gradientes de presin favorables ayudan a retrasar el punto de transicin del flujo laminar al
turbulento,mientrasquelosgradientesdepresindesfavorables,facilitanqueestatransicinaparezca
rpidamente.

Figura3.Regionesdeflujolaminaryturbulentoalrededordeunperfil.
Amedidaqueavanzamosendireccindelfluido,elflujoturbulentocontinadesarrollndoseysellega
aunmomentoenquepequeosgradientesdesfavorablesprovocanlaaparicindeunflujoensentido
opuesto,yconllaaparicindelfenmenodeseparacin.Amedidaqueseavanzaenlacuerda,elflujo
se hace ms sensible a separarse, de manera que pequeos gradientes desfavorables son capaces de
provocarestefenmeno.Porelcontrario,cercadelbordedeentrada
4
delperfilsonnecesariosgrandes
gradientesparaprovocarlaaparicindeseparacin.Porotraparte,esimportantehacernotarque,ante
gradienteslosuficientementeintensos,esposiblelaaparicindeseparacininclusoenlazonadeflujo
laminar. Pero, afortunadamente, en esos casos es posible que ante gradientes favorables, el flujo se
reacople, apareciendo una burbuja de flujo separado entre el punto de separacin y el de
reacoplamiento.
El aumento de empuje que se produce al incrementar el ngulo de ataque o la curvatura del perfil,
provoca el aumento de los gradientes favorables en la zona cercana al borde de ataque del perfil y en
cierto momento la aparicin de una burbuja de flujo desprendido. Si continuamos el incremento, la
longitud de la burbuja de separacin aumenta, y llega un momento en que se funde con la zona de
Di i d l Fl j Estela
(Flujo desprendido)
l j l i
Flujoturbulento
separacindelbordedesalida.Duranteesteprocesoelempujellegaatomarunvalormximo,apartir
delcualdisminuye.Enparalelo,laresistenciaaumentademaneradrstica(sedicequeelperfilentraen
prdida).
CURVA CARACTERSTICA DE UNPERFIL
Los valores de las resistencia y sustentacin deunperfildependen delnmero de Reynolds(calculado
con la cuerda del perfil) y del ngulo de ataque del caso. La variacin de estos valores, suele
representarseencurvascomolaquesemuestraenlaFigura 4yenlaFigura 6.Encadaunadeellas
podemosvertrescurvascaractersticasdeunmismoperfilquesecorrespondenaotrostantosnmeros
deReynolds.LascurvasrepresentanlavariacindelcoeficientedesustentacinC
L
ydelcoeficientede
resistenciaC
D
.Estoscoeficientessedefinenpor:
2 2
,
1 1
l l
2 2
L D
L D
C C
U U
= = ,

Figura4.Curvacaractersticadelcoeficientede
sustentacindeunperfil(Rn1>Rn2>Rn3).
Como puede verse, las curvas que aparecen en la Figura 4 tienen un comportamiento prcticamente
linealparapequeosngulosdeataque(porloquesepuedesuponerque,enesazonalasustentacin
depende nicamente del ngulo de ataque y no es sensible a variaciones del nmero de Reynolds),
siendolapendientedelacurvaunacaractersticapropiadelperfil.

-1.5
-1
-0.5
0
0.5
1
1.5
-30 -20 -10 0 10 20 30
C
L
ngulo de Ataque o (grados)
Rn 1
Rn 2
Rn 3 a

Figura5.Evolucindelcoeficientedesustentacin
deunaseccinNACA2412(resultadosdeJavafoil).
Sin embargo cuando el ngulo de ataque o incidencia aumenta (los valores habituales estn alrededor
delos15),elperfilseoponecadavezmsalavancedelfluidoysepropicialaaparicindefenmenos
turbulentosylaseparacindelflujo.

Figura6.Curvacaractersticadelcoeficientede
resistenciadeunperfil(Rn1>Rn2>Rn3).
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
-30 -20 -10 0 10 20 30
C
D
ngulo de Ataque o
Rn 1
Rn 2
Rn 3
A partir de ese punto el aumento del coeficiente de sustentacin con el ngulo de ataque es cada vez
menor,hastallegaraunpunto(ngulodeprdida)enelqueelperfilcomienzaaperdersustentacinde
maneramuyrpida.
La eleccin/diseo de un perfil ptimo ha de tener en cuenta el desarrollo de una distribucin de
presin adecuada que permita un adecuado crecimiento de la capa lmite a lo largo de la cuerda. El
objetivo final ser la optimizacin de la extensin del flujo laminar sobre la superficie de succin,
mediante el desarrollo de gradientes de presin favorables, a la vez que retrasar la aparicin de la
separacin,limitandolosvaloresdelosgradientesdepresinnegativos[
1
].
Por otra parte, la eleccin de un perfil u otro para una aplicacin concreta, debe tener en cuenta los
siguientespuntos:
1. Para ser efectivo, cualquier perfil debera trabajar con ngulos de ataque menores que el
ngulo de prdida. El rango de ngulos de ataque para los cuales el comportamiento del
perfilessatisfactorioesmayorparaperfilesanchosconradiosdecurvaturamayoresenel
bordedeataque.
2. La resistencia del perfil es mnima cuando opera cerca de su ngulo de ataque ideal. El
ngulo de ataque ideal se incrementa a medida que aumenta la curvatura del perfil,
aumentandoparalelamenteelempujegeneradoparaelngulodeataqueideal.
3. El nmero de Reynolds de operacindebe ser adecuadoal perfil en cuestin. Engeneral,
losperfilesslosonadecuadosparatrabajarenunrangodenmerosdeReynolds,quese
correspondeconlosdatosdisponiblesenlascurvascaractersticas.
4. Enoperacinnormal,lasustentacindelperfileslafuerzatil,mientrasquelaresistencia
eselprecioquehayquepagarparaconseguireseempuje.Porello,unposiblecriteriode
eleccinseraaquelperfil,cuyarelacinC
L
/C
D
fueramayorenelrangodeoperacin.
5. El valor de la sustentacin en el punto de diseo debe ser suficiente para una correcta
operacin.

PERFILES TPICOS
Entrelosperfilesmssencillosquesepuedenutilizarestnlosqueseconstruyenmedianteun
segmentodecrculo,elipseoparbolaquedefinelacaradesuccin,mientrasquelacaradepresines
rectaycoincideconlacuerda.Estetipodeperfilespuedenutilizarseenlaspalasdelashlicesdelos
buquesenlasseccionesmsalejadasdelejedebidoatenerunadistribucindepresinenlacarade
succinmsuniformeydisminuyeelpeligrodecavitacin.
Sinembargo,losperfilesmsconocidosfuerondesarrolladosapartirdelosaostreintaporlaNational
Advisory Committee for Aeronautics (NACA), ms tarde conocida como NASA, que se embarc en una
serie de ensayos en perfiles definidos de forma sistemtica y racional. Muchos de estos perfiles NACA
siguenusndoseenlaactualidad[2].
Una de las familias ms usadas son los perfiles NACA de la serie 6 de perfiles de flujo laminar,
desarrolladadurantelaSegundaGuerraMundial.UnejemplopuedeserelperfilNACA65218.Elprimer
dgito indica la serie; el segundo, el punto de mnima presin definido como en dcimos de la cuerda
medidadesdeelbordedeentrada
4
;eltercero,eselcoeficientedesustentacindediseodivididopor
diezylosdosltimosdgitosdanelmximoespesorentantoporcientodelacuerda
5
.
Durantelosaossetenta,laNASAdesarrollunaseriedeperfilesdebajavelocidadconun
comportamientosuperioralosperfilesNACAanteriores.Estasnuevasseriesfuerondiseadascon
ordenadorbasndoseentcnicasnumricas.

4
Elborde de entrada (leading edge) es el extremo delperfil pordondeincide el flujo,mientras que se
denominabordedesalida(trailingedge)alextremoopuestopordondeescapaelflujo.
5
Si el coeficiente de sustentacin de diseo o la posicin del mximo espesor no son valores enteros,
estascantidadessemuestranentreparntesis.

FLUJOSIRROTACIONALES YNO VISCOSOS
Cuando el flujo es no viscoso e irrotacional, la ecuacin de la cantidad de movimiento se
reduce a la ecuacin de Euler:
p
JF
Jt
= pg -Vp
La aceleracin es de la forma:
JF
Jt
=
oF
ot
+(F V)F =
oF
ot
+V_
1
2
I
2
] + F

Donde = rotF es la vorticidad del fluido.
Reorganizamos la ecuacin, dividimos por la densidad o multiplicamos por un vector
desplazamiento arbitrario dr .

_
oF
ot
+V_
1
2
I
2
] + F +
1
p
Vp -g_ Jr = u

El sumando ( F) Jr = u se anula en las condiciones siguientes:

1. La velocidad es cero, que es el caso trivial, no hay flujo (hidrosttica).
2. es cero, caso irrotacional.
3. Jr es perpendicular a F: este flujo es bastante peculiar y resulta poco comn.
4. Jr es paralelo a F: integramos a lo largo de una lnea de corriente.

Si integramos a lo largo de una lnea de corriente, en un flujo compresible no viscoso, y
tomamos por conveniencia g = -gk , la ecuacin anterior queda de la forma
oF
ot
Jr +J _
1
2
I
2
] +
Jp
p
+gJz = u

( ) ( )
} }
= + + +
c
c
2
1
1 2
2
1
2
2
2
1
0
2
1
z z g V V
dp
ds
t
V


Siendo ds el elemento diferencial de longitud a lo largo de la lnea de corriente. Esta ecuacin
es la ecuacin de Bernoulli para flujo no estacionario y no viscoso a lo largo de una lnea de
corriente. Para flujo incompresible estacionario se reduce a
te cons gz V
p
tan
2
1
2
= + +


La constante puede variar de una lnea de corriente a otra, a menos que el flujo sea tambin
irrotacional, 0 = , , en cuyo caso la constante ser la misma en todo el campo fluido.
Mostramos una corriente libre que se aproxima a dos cuerpos prximos entre s,, creando un
flujo interno entre ellos y un flujo externo por encima y por debajo de ellos. En la parte frontal
de los cuerpos hay una regin de gradiente favorable (la presin disminuye a lo largo de la
superficie) y la capa lmite, que estar adherida, ser delgada.


Fig. 4.4 Acoplamiento entre las regiones viscosas y no viscosas de un flujo. La teora potencial no es
aplicable a la zona de capa lmite

En los flujos internos, las capas lmite crecen desde las paredes, y al encontrarse desaparece
el ncleo no viscoso. La teora no viscosa debe funcionar bien en los flujos externos,
especialmente cerca de la parte frontal del cuerpo, hasta que el gradiente de presiones a lo
largo de la superficie se vuelve adverso (la presin aumenta) y la capa lmite se desprende. La
corriente desprendida deflecta y modifica las lneas de corriente del flujo exterior no viscoso,
que interacciona fuertemente con el flujo viscoso cerca de la pared.

POTENCIALDE VELOCIDADES
Si el flujo es irrotacional, se puede demostrar que la velocidad se puede igualar al gradiente de
una funcin escalar.
Si V F = u entonces F = V, donde ) , , , ( t z y x | | = recibe el nombre de potencial de
velocidades.
z
w
y
v
x
u
c
c
=
c
c
=
c
c
=
| | |
; ;

En este caso, la ecuacin de continuidad V F = u, se convierte en la ecuacin de Laplace para
.
v
2
=
o
2

ox
2
+
o
2

oy
2
+
o
2

oz
2
= u

Las lneas o superficies | constante se denominan lneas o superficies equipotenciales del
flujo.

Si el flujo es irrotacional,la ecuacin de Bernoulli no estacionaria se simplifica utilizando el
potencial de velocidades.

oF
ot
Jr =
o
ot
(V) Jr = J _
o
ot
]

A lo largo de cualquier direccin arbitraria la ecuacin de Bernouilli para movimiento irroracional
no estacionario, incompresible queda de la forma.
o
ot
+
p
p
+
1
2
I
2
+gz = ctc. JonJc I = |V|
Las condiciones de contorno para la ecuacin son:
Contornos exteriores,
q
x
,
q

,
q
z
conocidos
Componente de velocidad normal al contorno de la superficie del cuerpo, nula
q
n
=
u JonJc n cs pcrpcnJiculor ol cucrpo
Integrar la ecuacin de Laplace, con valores conocidos para la derivada de en el contorno, es
mucho ms sencillo que utilizar directamente las ecuaciones completas de Navier-Stokes. Hay
muchas tcnicas para encontrar las funciones potenciales que satisfacen la ecuacin de
Laplace, incluyendo la superposicin de funciones elementales, los mtodos numricos
mediante diferencias finitas, elementos finitos y elementos de contorno. Habiendo determinado
(x, y, z, t) mediante alguno de estos procedimientos, se determina F por derivacin.

REPASO DEL CONCEPTO DE FUNCIN DECORRIENTE
Si el flujo est descrito slo por dos coordenadas, existe tambin la funcin de corriente .
Para el flujo plano incompresible en coordenadas cartesianas x, y, la forma apropiada es
u =
o
oy
: = -
o
ox

La condicin de irrotacionalidad nos lleva de nuevo a la ecuacin de Laplace, en este caso
para la funcin de corriente.
C. IrrotocionoliJoJ: V F =
o:
ox
-
ou
oy
=
o
Jx
_-
o
ox
] -
o
oy
_
o
oy
] = u
o
2

ox
2
+
o
2

oy
2
= u
Las condiciones de contorno son de nuevo velocidad conocida en la corriente libre y flujo nulo a
travs de cualquier superficie slida.
Corriente exterior:

x
conocidas
Superficies slidas:
cucpo
= ctc
La ltima condicin se puede interpretar como que cualquier lnea constante puede
interpretarse como la pared de un cuerpo.
Si un flujo es irrotacional y puede describirse mediante slo dos coordenadas, existen tanto la
funcin de corriente como el potencial de velocidades | , y las lneas de corriente y
equipotenciales son ortogonales entre s excepto en los puntos de remanso, donde u y v son
cero. Por ejemplo, para el flujo incompresible en el plano xy , tendamos
y z
v
x y
u
c
c
=
c
c
=
c
c
=
c
c
=
| |
;
Las relaciones anteriores implican la ortogonalidad entre las lneas de corriente y las
equipotenciales.
Una vez determinadas, cualquier conjunto de lneas puede considerarse como las lneas
equipotenciales, y las otras sern lneas de corriente. Ambos conjuntos de lneas son
soluciones de la ecuacin de Laplace y pueden intercambiarse sus papeles.

Fig. 4.5 Las lneas de corriente y equipotenciales son ortogonales y pueden invertirse sus papeles si los
resultados son tiles: (a) flujo no viscoso tpico, (b) lo mismo que (a) pero con los papeles invertidos

Muchas de estas soluciones conviene escribirlas utilizando coordenadas polares (r, 0). Las
expresiones para las componentes de la velocidad en funcin de las derivadas de y
adoptan ahora la forma:
:

=
o
or
=
1
r
o
o0
:
0
=
1
r
o
o0
= -
o
or

Y la ecuacin de Laplace se describe como
1
r
o
or
_r
o
or
] +
1
r
2
o
2

o0
2
= u
GENERACIN DE ROTACIONALIDAD
Un flujo inicialmente irrotacional puede llegar a ser rotacional si:
1- Hay fuerzas viscosas apreciables inducidas por chorros, estelas o paredes slidas. En
este caso la ecuacin de Bernouilli no es vlida en estas regiones viscosas.
2- Hay gradientes de entropa originados por ondas de choque curvadas.
3- Hay gradientes de densidad originados por estratificacin ms que por gradiente de
presin.
4- Hay efectos no inerciales importantes tales como la rotacin de la tierra.
En los casos 2 y 4 la ecuacin de Bernouilli sigue siendo vlida a lo largo de una lnea de
corriente si la friccin es despreciable.
El caso 1, se da cerca de superficies slidas, donde la condicin de no deslizamiento crea una
capa lmite a travs de la cual la velocidad cae a cero (la velocidad ser cero en la superficie de
contacto con el cuerpo), y en chorros y estelas, donde corrientes con distintas velocidades
aparecen separadas por una capa delgada de cortadura muy intensa.
Los flujos internos, como en tubos y conductos, son principalmente viscosos, pues las capas
lmite de las paredes crecen hasta extenderse a todo el conducto. La ecuacin de Bernouilli no
es vlida en estos flujo, a menos que se modifique incluyendo los efectos de las prdidas
viscosas.

Fig. 4. 6 Modelos tpicos de flujo que muestran las regiones viscosas acopladas a regiones no viscosas:
(a) flujo subsnico alrededor de un cuerpo: flujo no viscoso e irrotacional fuera de la capa lmite y la
estela, (b) flujo supersnico alrededor de un cuerpo: flujo rotacional no viscoso fuera de la capa lmite y la
estela

Los flujos externos, como el flujo alrededor de cuerpos sumergidos en una corriente, son
parcialmente viscosos y parcialmente no viscosos; las dos regiones se acoplan en el borde de
la capa de cortadura o capa lmite. Se muestran dos ejemplos en la figura 4.6. La figura 4.6a
muestra un flujo subsnico alrededor de un cuerpo. La corriente incidente uniforme es
irrotacional, pues el rotor de una constante es cero, pero los esfuerzos viscosos originan una
capa de cortadura en torno y aguas abajo del cuerpo. En trminos generales, la capa lmite es
laminar, cerca del borde de ataque del cuerpo y turbulenta en la parte posterior. Generalmente
aparece una regin desprendida cerca del borde de salida, seguida de una estela turbulenta no
estacionaria, que se extiende lejos aguas abajo.
La figura 4.6b muestra un flujo supersnico alrededor de un cuerpo con borde de ataque
redondeado. Generalmente se forma una onda de choque curvada en la parte delantera, y el
flujo aguas abajo es rotacional debido a los gradientes de entropa. En este caso la zona de
separacin es pequea o inexistente, y la estela es generalmente ms delgada.

FLUJOS POTENCIALESPLANOS,SOLUCIONES ELEMENTALES
Ahora estudiaremos algunos flujos incompresibles no viscosos que poseen funcin de corriente
y potencial de velocidades. Se pueden obtener muchas soluciones haciendo uso de la
superposicin de funciones. As que empezaremos buscando las soluciones para flujos
elementales.
Cada uno de los tres flujos elementales es irrotacional e incompresible, por tanto, satisface las
dos ecuaciones de flujo potencial, v
2
= u y v
2
= u . Como ambas son ecuaciones
diferenciales en derivadas parciales lineales, cualquier suma de tales soluciones bsicas, es
una solucin. Algunas de estas composiciones de soluciones resultan muy interesantes y
prcticas.
CORRIENTEUNIFORME UNIDIMENSIONAL
Tenemos una corriente uniforme en la direccin del eje x, F = u|, esta corriente posee una
funcin de corriente y un potencial de velocidades, por lo que puede ser descrita como
u = u =
o
ox
=
o
oy
: = u =
o
oy
= -
o
ox

Se puede integrar cada expresin y despreciar las constantes de integracin, lo cual no afecta
a las velocidades del flujo (ya que el flujo se obtiene derivando o , al derivar se pierden las
constantes).
As, para una corriente uniforme |u tenemos = uy = ux.
Las lneas de corriente son lneas rectas horizontales (y=cte.) y las lneas equipotenciales son
verticales (x=cte.), ortogonales a las lneas de corriente como era de esperar.

Fig. 4.7 Corriente uniforme unidimensional. Las lneas slidas son lneas de corriente y las discontinuas
son lneas equipotenciales.
FUENTE O SUMIDEROBIDIMENSIONAL EN ELORIGEN

Fig. 4.8 Fuente. Las lneas slidas son lneas de corriente y las discontinuas equipotenciales

Supongamos que el eje z (que se adentra en el papel) fuese una tubera delgada a travs de la
cual fluye un caudal total uniformemente largo de su longitud b. Fijando la atencin en el
plano xy, se vera un flujo radial hacia fuera o fuente. Es conveniente tratar el problema en
coordenadas polares. Considerando que no existe velocidad circunferencial, para cualquier
radio r, la velocidad es
:

=

2nrb
=
m
r
=
1
r
o
o0
=
o
or
:
0
= u = -
o
or
=
1
r
o
o0

Integrando, descartando las constantes de integracin, se obtienen las funciones para este flujo
radial simple:
Fuente o sumidero: = m0 = m lnr
Donde m = (2nb) es una constante, positiva para una fuente y negativa para un sumidero.
Las lneas de corriente son rectas radiales (0 constante) y las lneas equipotenciales son
circunferencias (r constante).
TORBELLINO IRROTACIONAL EN EL ORIGEN

Fig. 4.9 Torbellino. Las lneas slidas son lneas de corriente y las discontinuas son equipotenciales

Un torbellino bidimensional en el origen es un movimiento estacionario puramente circulatorio,
:
0
= (r) solamente y :

= u. Se puede demostrar que slo una funcin :


0
(r) es irrotacional,
es decir, V F = u y :
0
= Kr , donde K es una constante. Para este flujo, llamado en
ocasiones vrtice libre, la funcin de corriente y el potencial de velocidades se obtienen de
:

= u =
1
r
o
o0
=
o
or
:
0
=
K
r
= -
o
or
=
1
r
o
o0

Integrando, se determinan las funciones apropiadas
= -Klnr = K0
Donde K es una constante llamada intensidad del torbellino. Las lneas de corriente son
crculos (r constante) y las lneas equipotenciales son rectas radiales (0 constante). El flujo
originado al vaciarse un depsito por un orificio inferior central responde aproximadamente al
patrn del torbellino potencial.
El flujo inducido por un torbellino bidimensional es irrotacional en toda partes, excepto en el
origen, donde la vorticidad V F es infinita. Esto significa que una cierta integral de lnea
denominada circulacin no se anula cuando se integra a lo largo de un circuito que encierra el
ncleo del torbellino.
La circulacin se define como la integral a lo largo de una curva cerrada C, en el sentido
contrario al de las agujas del reloj, de la componente de la velocidad tangente a la curva por la
longitud del arco Js del elemento de curva:
= _ Icoso Js
C
= _ F Js
C
= _ (uJx +:Jy +wJz)
C

De la definicin de , F Jx = V Jx = J para un flujo irrotacional, as que generalmente
en un flujo irrotacional sera igual al valor final menos el inicial de . Puesto que la integral es
cerrada, empezamos y acabamos en el mismo punto, por lo que se obtendra = u, pero no
para un torbellino al ser = K0, ya que cada vez que da una vuelta completa 0 se incrementa
en 2n, de manera que en la curva que encierra el torbellino se cumple
= 2nK
En general es igual a la suma algebraica de las intensidades de todos los torbellinos que hay
en la regin interior de la curva cerrada. Una regin de circulacin finita en una corriente est
sometida a una fuerza de sustentacin proporcional a u y .
Si no hay torbellinos en el campo fluido, la circulacin ser cero alrededor de cualquier curva
cerrada que encierre un nmero arbitrario de fuentes y sumideros.

SUPERPOSICIN DEFLUJOS PLANOSELEMENTALES
Ahora podemos construir una variedad de flujos potenciales interesantes, sin ms que sumar
los potenciales de velocidad y las funciones de corriente de corrientes uniformes, fuentes,
sumideros y torbellinos. La superposicin es vlida porque las ecuaciones bsicas v
2
= u y
q
t
+
p
p
+
1
2
I
2
+gz = ctc., son lineales.
SUMIDEROMS TORBELLINO EN EL ORIGEN
Cuando se superponen un sumidero y un torbellino, ambos centrados en el origen, la
composicin de funcin de corriente y potencial de velocidad resultante es
= m0 -Klnr = mlnr +K0

Fig. 4.10 Superposicin de un sumidero ms un torbellino. Simula un tornado


La representacin est formada por dos familias de espirales logartmicas ortogonales. Este
patrn es una simulacin bastante realista de un tornado (donde el sumidero asciende segn el
eje z en la atmsfera) o bien en el torbellino que se forma en el vaciado rpido de una baera.
CORRIENTEUNIFORME MS FUENTE EN EL ORIGEN: CUERPO SEMIINFINITODE
RANKINE
Si enfrentamos una corriente uniforme en la direccin x frente a una fuente aislada, aparece
una forma de cuerpo semiinfinito. Si la fuente se posiciona en el origen, la superposicin de
funciones de corriente resultante es, expresada en coordenadas polares
= urscn0 +m0


Fig. 4.11 La superposicin de una corriente uniforme ms una fuente forma un cuerpo semiinfinito de
Rankine.
En la familia de curvas, puede apreciarse una forma curva, ms o menos elptica, como si
fuese el contorno de un cuerpo semiinfinito, que recibe el nombre de cuerpo semiinfinito de
Rankine. La forma del cuerpo se corresponde con las lneas de corriente particulares = nm.
La anchura media del cuerpo es, a una cierta distancia aguas abajo igual a nmu.
La superficie superior puede trazarse a partir de la expresin
r =
m(n -0)
uscn0

La nariz del cuerpo, es un punto de estancamiento, I = u, y est situado en (x, y) = (-o, u),
donde o = mu.
FILAINFINITA DETORBELLINOS
Consideramos una fila infinita de torbellinos de la misma intensidad K equiespaciados una
distancia o.

Fig. 4.12 Fila de torbellinos

El torbellino i tiene por funcin de corriente

= -Klnr

, de modo que la funcin de corriente


de la fila infinita es
= -Klnr

=1

Puede demostrarse que esta suma infinita de logaritmos es equivalente a la funcin
= -
1
2
Kln _
1
2
_cosb
2ny
o
-cos
2nx
o
]_
A continuacin se muestran las lneas de corriente

Fig. 4.13 Lneas de corriente para una fila de torbellinos

Se observa la configuracin llamada ojo de gato, con celdas de recirculacin que rodean a los
torbellinos individuales. Por encima de los ojos de gato, el flujo es hacia la izquierda y por
debajo es hacia la derecha. Adems estos flujos hacia izquierda y derecha se hacen uniformes
para |y| o
u =
o
oy
_
||u
=
nK
o

El signo ms corresponde al flujo por debajo de la fila y el menos por encima. Podemos
representar las corrientes uniformes

Fig. 4.14 Flujo de la fila de torbellinos visto desde lejos

Este efecto est inducido por torbellinos, en este ejemplo no hay corriente uniforme hacia la
fila.

CAPADE TORBELLINOS
Si observamos la fila infinita de torbellinos desde lejos, lo que se ve son las corrientes
uniformes hacia la derecha, por arriba, y hacia la izquierda, por abajo y los torbellinos parecen
estar tan prximos los unos a los otros que se ven como una capa de torbellinos contnua. La
intensidad de la capa se define como
y =
2nK
o

En el caso ms general y puede variar con x. La circulacin alrededor de cualquier curva
cerrada que encierre una longitud Jx de la capa, ser
J = u

Jx -u
s
Jx = (u

-u
s
)Jx =
2nK
o
Jx = yJx
Donde los subndices i y s significan inferior y superior, respectivamente. Por tanto, la
intensidad de la capa es la circulacin de la capa por unida de longitud y = JJx. Cuando una
capa de torbellinos est inmersa en una corriente uniforme, y es proporcional a la sustentacin,
por unidad de longitud de cualquier superficie que rodee la capa.
No hay velocidad perpendicular a la capa en la superficie de la misma, por tanto, una capa de
torbellinos puede simular un cuerpo delgado, como una placa o un perfil delgado. Esta es la
base de la teora de perfiles delgados que se utiliza por ejemplo en la modelizacin del flujo
alrededor de las velas.
DOBLETE
Cuando nos situamos lejos del par fuente-sumidero.

Fig. 4.15 Flujo potencial debido a una fuente ms un sumidero de igual intensidad. Las lneas slidas son
lneas de corriente y las discontinuas son equipotenciales

El patrn de flujo comienza a parecerse a una familia de crculos tangentes en el origen.

Fig. 4.16 Un doblete, o par fuente sumidero, constituye el caso lmite de la figura anterior vista desde
lejos. Las lneas de corriente son crculos tangentes al eje x en el origen.

En este lmite, en el que la distancia o se hace muy pequea, la configuracin se llama doblete.
La funcin de corriente para esta configuracin es
= -
zy
x
2
+y
2

El parmetro z se denomina intensidad del doblete.
Las lneas de corriente son crculos con centro en el eje y y tangentes en el origen y se rigen
por la ecuacin
x
2
+_y +
z
2
]
2
= _
z
2
]
2

El potencial de velocidad es =
xx
x
2
+
2
y las lneas equipotenciales cumplen
y
2
+_x -
z
2
]
2
= _
z
2
]
2

Las lneas equipotenciales son crculos tangentes en el origen con sus centros en el eje x. Se
obtienen girando 90 en el sentido de las agujas del reloj la figura anterior, de manera que son
perpendiculares a las lneas de corriente.

FLUJO ALREDEDOR DEUN CILINDRO CONCIRCULACIN
La funcin de corriente para el flujo alrededor de un cilindro circular con circulacin centrado en
el origen es la de una corriente uniforme ms un doblete y un torbellino situados en el origen.
= u

scn0 _r -
o
2
r
_ -Kln
r
o

Dibujamos las lneas de corriente para cuatro valores distintos de la intensidad adimensional
del torbellino K(u

o).

Fig. 4.17 Flujo alrededor de un cilindro circular con circulacin para valores de K(u

o). (a) 0.0, (b) 1.0,


(c) 2.0 y (d) 3.0

En todos los casos, la lnea = u corresponde al crculo de radio r = o, esto es, al cuerpo de
forma cilndrica. Cuando la circulacin = 2nK aumenta, crece la velocidad en la parte inferior
del cilindro y decrece en la parte superior. Las componentes de la velocidad estn dadas por
:

=
1
r
o
o0
= u

cos 0 _1 -
o
2
r
2
_
:
0
= -
o
or
= -u

sen0 _1 +
o
2
r
2
_ +
K
r

La velocidad en la superficie r = o es tangencial
:

(r = o) = u :
0
(r = o) = -2u

sen0 +
K
o

Para valores pequeos de K hay dos puntos de remanso sobre la superficie del cilindro,
situados en los ngulos 0
s
en los que se cumple :
0
= u. Se ha de cumplir por tanto scn0
s
=
K
20

u
.
El caso (a) de la figura anterior corresponde a K = u, 0
s
= u y 18u, esto es un flujo no viscoso
doblemente simtrico alrededor de un cilindro circular con circulacin. El caso (b) corresponde
a
K
(0

u)
= 1, 0
s
= Su y 1Su, y el caso (c) corresponde al caso lmite
K
(0

u)
= 2, en que los dos
puntos de remanso coinciden en el punto ms alto del cilindro, 0
s
= 9u.
Para K > 2u

o la ecuacin scn0
s
=
K
20

u
no es vlida, y slo hay un punto de remanso fuera
del cilindro, como en el caso (d) de la figura.
Para los casos (b) a (d) de flujos alrededor de cilindros hay una fuerza vertical hacia abajo, o
sustentacin negativa, denominada efecto Magnus-Robins,que es proporcional a la velocidad
de la corriente uniforme y a la intensidad del torbellino. Del esquema de las lneas de corriente
se deduce que la velocidad en la parte superior del cilindro es menor que en la parte inferior, y
segn la ecuacin de Bernouilli, la presin es ms alta en la parte superior que en la inferior, lo
que explica que exista esta fuerza hacia abajo. Por supuesto, estas fuerzas no son viscosas, ya
que la teora utilizada es no viscosa.
La velocidad en la superficie es
:

(r = o) = u :
0
(r = o) = -2u

sen0 +
K
o

As para la lnea de corriente que entra se ha de cumplir la ecuacin de Bernouilli
p

+
1
2
pu

= p
s
+
1
2
p _-2u

sen0 +
K
o
]
2

De manera que podemos encontrar la presin en la superficie
p
s
= p

+
1
2
pu

(1 -4scn
2
0 +4[scn0 -[
2
)
[ =
K
0

u
y p

es la presin de la corriente incidente. Si b es la anchura del cilindro


perpendicular al papel, la resistencia D es la integral sobre la superficie de la componente
horizontal de las fuerzas de presin:
= -_ (p
s
-p

)cos0
2n
0
bo J0
Donde (p
s
-p

) =
1
2
pu

(1 -4scn
2
0 +4[scn0 -[
2
) y si resolvemos la integral obtenemos
(cilinJro con circulocion) = u. Este es un caso particular de la paradoja de DAlembert: De
acuerdo con la teora no viscosa, cualquier cuerpo de forma arbitraria inmerso en una corriente
uniforme no tiene resistencia.
DAlembert public este resultado en 1752, indicando l mismo que no concordaba con lo que
ocurra en los flujos de fluidos reales. Esta desafortunada paradoja dio pie a una reaccin
exagerada y todos rechazaron las teoras no viscosas, hasta que Prandtl, en 1904, mostr cul
era el efecto, tan importante en el flujo, de la delgada capa lmite viscosa en la parte posterior
del cuerpo.
La sustentacin I perpendicular a la corriente incidente, tomada positiva hacia arriba, est
dada por la integral de las fuerzas verticales de presin:
I = -_ (p
s
-p

) scn0
2n
0
bo J0
Al igual que antes, sustituimos (p
s
-p

) y como la integral entre 0 y 2 de cualquier potencia


impar de scn0 es cero, la integral se reduce a
I = -
1
2
pu

2
4K
ou

bo _ scn
2
0 J0
2n
0
= -pu

(2nK)b
I
b
= -pu


Esta ecuacin fue generalizada por W.M. Kutta en 1902 e independientemente por N.
Joukowski en 1906, de la siguiente forma:
De acuerdo con la teora no viscosa, la sustentacin por unidad de envergadura de un cilindro
de forma arbitraria inmerso en una corriente uniforme es igual a pu

, donde es la circulacin
total alrededor del cuerpo. La direccin de la sustentacin se obtiene girando 90 la direccin
de la corriente incidente, en el sentido opuesto a la circulacin.
El problema principal del anlisis de perfiles consiste en determinar la circulacin como
funcin de la forma y orientacin de los mismos.
Los flujos analizados en la figura anterior son matemticos, un doblete y un torbellino situado
en el mismo punto, ms una corriente uniforme. Pero se podra conseguir un modelo para su
representacin fsica haciendo girar un cilindro en una corriente. La condicin de no
deslizamiento en un fluido viscoso obliga al fluido en contacto a moverse tangencialmente con
la velocidad :
0
= o, consiguiendo una circulacin neta .
APLICACIN AL ESTUDIO DE PERFILES
La anterior teora puede aplicarse al estudio de perfiles, para determinar la circulacin en
funcin de la forma del perfil y del ngulo de ataque de la corriente incidente o.
CONDICIN DE KUTTA
Aunque la forma del perfil y el ngulo de ataque sean conocidos, la solucin proporcionada por
la teora potencial no es nica, se puede encontrar una familia infinita de soluciones, cada una
de ellas correspondiente a un valor de . En la siguiente figura mostramos tres soluciones para
la corriente alrededor de un perfil que son matemticamente aceptables.


4.18 La condicin de Kutta simula de forma apropiada el flujo alrededor de un perfil; (a) circulacin menor
de la necesaria, el punto de remanso posterior est en la superficie superior, (b) circulacin excesiva, el
punto de remanso posterior est en la superficie inferior, (c) circulacin correcta, la condicin de Kutta
implica que la corriente abandona el borde de salida suavemente.

En el caso (c), que cumple la condicin de Kutta, los flujos superior e inferior se encuentran y
abandonan el borde de salida suavemente. Si el borde de salida es ligeramente redondeado,
habr all un punto de remanso. Si el borde de salida es afilado, como en la mayora de los
perfiles, las velocidades del fluido en las superficies superior e inferior deben ser iguales al
abandonar el perfil. Este razonamiento fsico proporciona el valor apropiado de y se atribuye
generalmente a W. M. Kutta, de ah el nombre de condicin de Kutta. El valor correcto de la
circulacin
Kuttu
depende de la velocidad incidente, del ngulo de ataque y de la forma del
perfil.
PLACAPLANA
Analizamos el perfil ms sencillo, una placa plana, para resolver el problema utilizamos la capa
de torbellinos. Suponemos que la placa se puede modelar como una capa de torbellinos de
intensidad variable y(x). La corriente libre u

forma un ngulo de ataque o con la cuerda de la


placa.

Fig. 4. 19 Solucin para la placa plana con ngulo de ataque utilizando la capa de torbellinos.

De manera experimental se sabe que la placa se sustentar, sabemos que la sustentacin es
hacia arriba, por lo que la circulacin ha de ser en sentido horario. (Recordamos: La direccin
de la sustentacin se obtiene girando 90 la direccin de la corriente incidente, en el sentido
opuesto a la circulacin).
La capa de torbellinos debe originar un flujo hacia la derecha ou = +
1
2
y en la superficie
superior, e igual y en sentido contrario en la superficie inferior.
Se puede calcular la sustentacin total sumando la intensidad de la capa de torbellinos sobre
todo el perfil. Para un perfil de anchura b:
I = pu

b = _ (x)ux
C
0

Una forma alternativa de determinar la sustentacin es a partir del coeficiente adimensional de
presin C
p
en las superficies superior e inferior:
C
p
i,s
=
p
s,
-p

1
2
pu

2
= 1 -
u
,s
2
u

2

La ltima igualdad se obtiene utilizando la ecuacin de Bernouilli. El cuadrado de la velocidad
en la superficie se obtiene combinando la corriente uniforme y las componentes de la velocidad
debidas a la capa de torbellinos.
u
,s
2
= (u

coso ou)
2
+(u

scno)
2
= u

2
2u

ou coso +ou
2
= u

2
_1
2ou
u

]
La sustentacin es la integral de la diferencia de presiones extendida a toda la longitud del
perfil supuesto de anchura b:
I = _ (p
I
-p
u
)b Jx
C
0

C
L
=
I
1
2
pu

2
bC
= _ (C
p
i
-C
p
s
)
Jx
C
1
0
= 2_
y
u

J [
x
C

1
0

La intensidad de la capa y(x) se determina de la condicin de que la velocidad normal :(x) es
cero en y = u, ya que representa una placa slida o superficie de corriente. La velocidad
normal en el punto x inducida por toda la capa es:
:
cupu
= _
-yJx
2n(x
0
-x)
C
0

La corriente uniforme induce una velocidad normal constante en cada punto de la capa dada
por
:
cocntc
= u

scno
Para que la velocidad normal en la superficie sea nula, se ha de cumplir que la suma de las
velocidades normales inducidas por la capa y por la corriente, sea igual a cero:
:
cupu
+:
cocntc
= u _
yJx
2n(x
0
-x)
C
0
= 2nu

scno
La integral debe resolverse para y(x) con la condicin de Kutta y(C) = u (La condicin de Kutta
recordamos que impone que en el borde de salida la velocidad en la superficie superior ha de
ser igual a la velocidad en la parte inferior).
El resultado final que se obtiene, despus de muchos clculos, es:
y(x) = 2u

scno _
C
x
-1]
12

De manera que el coeficiente de presin en la superficie, queda de la forma:
C
p
s,i
= 2scno _
C
x
-1]
12


Dibujamos los coeficientes de presin en funcin de la posicin x en la placa.

Fig. 4. 20 Coeficiente terico de presin sobre las superficies superior e inferior

Observamos que en la superficie superior la presin aumenta continuamente con x, hay un
gradiente adverso de presin.
Representamos la velocidad en la superficie superior u
s
= u

+ou = u

+
1
2
y para varios
ngulos de ataque:

Fig. 4.21 Velocidad en la superficie superior, donde los puntos D indican
el punto de desprendimiento de la capa lmite laminar

Se observa que para ngulos de ataque por encima de o = S, la contribucin ou de la capa es
alrededor de 20% de u

, por lo que se viola la ley de pequeas perturbaciones. En la figura


tambin se muestran los puntos de desprendimiento calculados por el mtodo de Thwaites, y
se observa que a mayores ngulos de ataque, antes se produce el desprendimiento. La
prediccin, aproximadamente correcta, es que la placa plana sufre un desprendimiento masivo
en la superficie superior que provoca la entrada en prdida para o > 6.
El coeficiente de sustentacin del perfil
C
L
= 2_
y
u

J [
x
C

1
0
= 4 scno _ _
C
x
-1]
12
J [
x
C

1
0
= 2nscno = 2no
Tambin es interesante el coeficiente de momento alrededor del borde de ataque (BA) del
perfil, considerado positivo en el sentido de las agujas del reloj:
C
M
BA
=
H
BA
1
2
pu

2
bC
2
= _ (C
p
s
-C
p
i
)
x
C
J [
x
C

1
0
=
n
2
scno =
1
4
C
L

Por lo tanto el centro de presiones (CP), o posicin de la sustentacin resultante, est situado
en el punto un cuarto de la cuerda, independientemente de cul sea el ngulo de ataque.
[
x
C

CP
=
1
4

ANLISIS NUMRICO
Cuando el flujo potencial presenta geometras complicadas o condiciones de corriente
inusuales, el mtodo clsico de superposicin de flujos resulta poco atractivo. En este caso el
enfoque ms apropiado es utilizando el anlisis numrico, del que existen, al menos, tres
mtodos distintos:
1- El mtodo de los elementos finitos (FEM, Finite Element Method)
El mtodo de elementos finitos se puede aplicar a todos los tipos de ecuaciones
diferenciales en derivadas parciales, tanto lineales como no lineales, de la fsica y la
ingeniera. El dominio computacional se divide en pequeas regiones, o celdas,
tpicamente con forma de tringulo o cuadriltero. Las celdas se definen utilizando un
nmero finito de nodos donde queremos calcular las variables de campo, como la
temperatura, la velocidad, presin, funcin de corriente, etc. La solucin en cada celda
se aproxima por una combinacin algebraica de los valores nodales locales. A
continuacin se integran estas funciones aproximadas sobre la celda y se minimiza el
error, para lo que suelen utilizarse funciones de peso. Se obtiene as un conjunto de N
ecuaciones algebraicas para los N valores nodales incgnita. Las ecuaciones nodales
se deben resolver de forma simultnea, invirtiendo una matriz o mediante iteracin.

2- El mtodo de diferencias finitas (FDM, Finite Difference Method)
La idea de este mtodo es aproximar las derivadas parciales que aparecen en la
ecuacin fsica por diferencias entre los valores de la solucin en una serie de nodos
separados entre s, una cierta distancia finita, si bien los nodos no tienen que estar
equiespaciados. La ecuacin original en derivadas parciales se sustituye as por una
serie de ecuaciones algebraicas para los valores nodales. Para el flujo potencial (no
viscoso), estas ecuaciones algebraicas son lineales, pero en general, para el flujo
viscoso son no lineales. Finalmente para obtener los valores nodales se debe iterar o
invertir una matriz.
Como ejemplo podemos escribir la ecuacin de Laplace para la funcin de corriente
haciendo uso del mtodo de diferencias finitas.
o
2

ox
2
+
o
2

oy
2
= u
Primero dividimos el capo fluido utilizando nodos equiespaciados.

Fig. 4.22 Definicin esquemtica de una malla rectangular de diferencias finitas para un
problema bidimensional.

Las derivadas se escriben utilizando una aproximacin algebraica:
o
ox
=
(x +x, y) -(x, y)
x
=
1
x
(
+1,]
-
,]
)
o
oy
=
(x, y +y) -(x, y)
y
=
1
y
(
,]+1
-
,]
)
o
2

ox
2
=
1
x
_
(x +x, y) -(x, y)
x
-
(x, y) -(x -x, y)
x
_
=
1
(x)
2
(
+1,]
-
,]
+
-1,]
)

o
2

oy
2
=
1
y
_
(x, y +y) -(x, y)
y
-
(x, y) -(x, y -y)
y
_
=
1
(y)
2
(
,]+1
-
,]
+
,]-1
)
De manera que la ecuacin de Laplace se transforma en una ecuacin algebraica
lineal.
1
(x)
2
(
+1,]
-
,]
+
-1,]
) +
1
(y)
2
(
,]+1
-
,]
+
,]-1
) = u

Si la malla es cuadrada x = y
1
4
(
+1,]
+
-1,]
+
,]+1
+
,]-1
) =
,]

Este modelo en diferencias finitas de la ecuacin de Laplace indica que cada valor
nodal de la funcin de corriente
,]
es una combinacin lineal de las funciones de
corriente de sus cuatro vecinos ms prximos.
El error numrico en relacin a la solucin exacta de la ecuacin de Laplace es
proporcional al cuadrado del tamao de la celda computacional.
La resolucin del sistema de ecuaciones o la iteracin se puede programar usando
cualquier lenguaje de programacin.

3- A. Mtodos integrales de singularidades distribuidas
B. El mtodo de los elementos de contorno
Esta es una tcnica relativamente nueva para la resolucin de ecuaciones diferenciales
en derivadas parciales es el mtodo de los elementos de contorno (BEM, Boundary
Element Method). En este mtodo no existen elementos interiores, todos los nodos se
sitan en la frontera del dominio.

Fig. 4.23 Elementos de contorno de intensidad constante para flujo potencial plano

Cada elemento es una pequea regin del contorno que rodea al nodo
correspondiente, cuya intensidad puede ser constante o variable.
Cada elemento i tiene una intensidad distinta

y r representa la distancia desde


dicho elemento hasta cualquier otro punto del campo fluido. Sumando los efectos de
todos los elementos e imponiendo las condiciones de contorno apropiadas se obtiene
la solucin final del problema de flujo potencial.
ALAS DEALARGAMIENTO FINITO

4.24 Perfil tpico grueso y con curvatura
Los resultados previos de la teora de perfiles son vlidos para alas bidimensionales o de
alargamiento infinito. Pero todas las alas reales tienen extremos y son, por tanto, de
envergadura finita o relacin de aspecto RA finito, definido como
RA =
b
2
A
p
=
b
C

Siendo b la envergadura o distancia entre extremos o puntas del ala y A
p
es el rea de la
forma en planta del ala vista desde arriba. Los coeficientes de sustentacin y resistencia de un
ala de alargamiento finito dependen fuertemente de la relacin de aspecto y muy poco del rea
de la forma en planta.

4.25 Teora de la lnea sustentadora para alas finitas: (a) sistema de torbellinos real en la estela de un ala;
(b) simulacin del sistema de torbellinos ligados al ala; (c) velocidad vertical inducida en el ala debida a
un elemento infinitesimal de torbellinos desprendidos

Los torbellinos no pueden terminar en el fluido; o bien se extienden hasta los contornos o bien
forman un circuito cerrado. La figura 4.25a muestra cmo los torbellinos que proporcionan la
circulacin alrededor del ala se curvan aguas abajo en los extremos de un ala de alargamiento
finito, alinendose con la corriente para unirse lejos aguas abajo formando un torbellino de
arranque.

Fig. 4.26 Etapas del desarrollo de la sustentacin: (a) sin sustentacin (b) el borde de salida afilado
induce la separacin y se forma un torbellino de arranque: sustentacin pequea (c) el flujo arrastra el
torbellino de arranque, el flujo en el borde de salida es suave (d) el torbellino de arranque arrastrado
aguas abajo
Los torbellinos de mayor intensidad se desprenden de los extremos, pero algunos se
desprenden del interior del ala, como muestra esquemticamente la figura 4.25b. La circulacin
efectiva (y) de los torbellinos desprendidos es cero en los bordes y, generalmente, tiene un
mximo en el plano central, o raz del ala. En 1918 Prandtl model este flujo de forma
satisfactoria reemplazando el ala por una lnea sustentadora y una capa semiinfinita de
torbellinos de intensidad y(x) =
dI
d
, como en la figura 4.25c. Cada torbellino elemental y(p)Jp
induce una velocidad vertical Jw(y) dada por
Jw(y) =
y(p)Jp
4n(y -p)

En el punto y de la lnea sustentadora. El factor 4n es causa de que los torbellinos se extienden
desde 0 hasta .
La velocidad total descendente w(y) inducida por el sistema completo de torbellinos
desprendidos es
w(y) =
1
4n
_
y(p)Jp
(y -p)
1
2
b
-[
1
2
b

Cuando esta velocidad vertical se suma vectorialmente con la corriente incidente u

, el ngulo
de ataque efectivo en esta seccin del ala es
o
c]
= o -o

= tg
-1
w
u

=
w
u


Tomando la aproximacin w u

.
La teora de Prandtl de la lnea sustentadora para alas de alargamiento finito:
(y) = nC(y)u

_o(y) -
1
4nu

_
y(p)Jp
(y -p)
1
2
b
-[
1
2
b
_
Que es una ecuacin integrodiferencial para (y) con las condiciones [
1
2
b = [-
1
2
b = u.
Una vez que se ha resuelto esta ecuacin, la sustentacin del ala y la resistencia inducida
estn dadas por
I = pu

_ (y)Jy
1
2
b
-[
1
2
b

= pu

_ (y)o

(y)Jy
1
2
b
-[
1
2
b

Tenemos aqu un caso en el que la resitencia no es nula en un flujo no viscoso, debido a que la
velocidad vertical hace que la sustentacin se incline hacia atrs un ngulo o

de modo que
proporciona una componente de resistencia paralela a la direccin de la corriente incidente,
J

= JIscno

= JI o

.
La solucin ms sencilla se da en el caso de un ala de forma en planta elptica y sin torsin,
dada por la funcin C(y) = C
0
_1 -[
2
b

2
_
12
.
El rea y alargamiento del ala son
A
p
= _ CJy
1
2
b
-[
1
2
b
=
1
4
nbC
0
RA =
4b
nC
0

La distribucin de circulacin es pues
(y) =
0
_1 -_
2y
b
]
2
_
12

La sustituimos en (y) = nC(y)u

_o(y) -
1
4n0

]
y(q)dq
(-q)
1
2
b
-[
1
2
b
_ e integramos, de manera que se
obtiene
0
=
nC
0
0

u
1+2RA
donde o es constante a lo largo del ala sin torsin.
La sustentacin para el ala elptica es:
I =
1
4
n
2
bC
0
pu

2
o(1 +2RA)
C
L
=
2no
1 +2RA

Si generalizamos esto a un ala finita gruesa con curvatura y de forma en plana
aproximadamente elptica, tenemos:
C
L
=
2nscn(o +[)
1 +2RA

[ = tg
-1
2b
C

EJEMPLO: ANLISIS NUMRICODEVELAS.
El problema de las velas es tridimensional, para resolverlo se utiliza el denominado mtodo de
los paneles. En este mtodo se discretizan las geometras, en nuestro caso las velas, y se
distribuyen elementos singulares en la superficie.
El mtodo ms econmico, computacionalmente hablando, para resolver el problema de una
superficie fina sustentable es utilizar como elemento singular los anillos de vrtice.
Primero se discretiza la geometra por medio de elementos o paneles cuadrilteros con sus
nodos colocados sobre la superficie y luego se distribuyen sobre ella los elementos singulares.
Como la superficie de la vela es fina, consideraremos un espesor nulo, por lo que usaremos los
mismos elementos para la cara de succin y la de presin, y esto se puede lograr gracias al
empleo de singularidades tipo vrtice o dobletes.


Fig. 4.27
Cada anillo de vrtice est formado por cuatro segmentos de vrtice, por lo que la velocidad
inducida en un punto P por cada anillo se puede calcular sumando las velocidades inducidas
en este punto, por cada uno de los segmentos que forman el anillo.
El cuerpo se discretiza en N paneles y la estela en N
w
paneles. Se considera que en cada
punto P slo influyen los paneles que se encuentran a ms de 2.5-5 veces la diagonal de un
panel.
Con el mtodo de los paneles obtenemos el campo de velocidades como si todo el fluido fuera
potencial.
ANEXO:MANUAL DE USUARIO DEJAVAFOIL
Javafoil es un programa de uso libre para clculo de perfiles, basado en:
- Resolucin numrica del flujo potencial mediante un mtodo de paneles de alto orden
- Resolucin numrica de las ecuaciones de capa lmite (mtodo integral)

La herramienta est disponible en: http://www.mh-aerotools.de/airfoils/javafoil.htm
Javafoil se presenta en un entorno de una nica ventana, con diversas pginas o tabs, que se
describen a continuacin.
OPTIONSCARD
This card offers some information about your Java system and it contains a combo box to select
a different country setting. The country setting also affects the decimal separator.Initially, the
language will be selected automagically, based on your system settings (or according to your
command line parameters). The default language is English, but if you prefer your native
language, contact me by eMail to receive a file with the character strings to translate.

You can save and restore the current state of JavaFoil to revert later to a previous project
(seesecuritysettings).
Also you can specify some properties of the fluid where you want the airfoil to operate. The
KINEMATIC VISCOSITYis needed for the calculation of the localReynolds numberand theSPEED OF
SOUNDis needed for theMach number. Currently, these parameters are currently only required
for the Aircraft card.
The aspect ratio is used for an approximate correction of the results on the Polar and Aircraft
cards for a finite wing. First the 3D lift coefficient CLis determined by adapting the 2D Cl. Mach
number and aspect ratio are taken into account. Then the 3D drag coefficient CDis calculated
by adding the induced drag coefficient for a wing with elliptical lift distribution to the Cdof the
airfoil.
Finally, the scripting facility can be used to automate simple command sequences.

GEOMETRY CARD

The Geometry card is used to prepare the coordinates of your current airfoil. This is the
geometry which is used by all tools in JavaFoil.


It shows a list of coordinates and a plot of the resulting airfoil shape. You can enter or paste
arbitrary coordinates into this field and press the "Update View" button to copy the coordinates
into the working airfoil. Remember that the coordinates must be ordered trailing edge > upper >
nose > lower > trailing edge.
Additionally it contains options to create standard airfoils from several families, namely:
- 4-digit series, like NACA 2415, based on
o maximum thickness and its x/c location
o maximum camber and its x/c location
- 5 digit series like NACA 23015, based on
o maximum thickness and its x/c location
o design lift coefficient
o x/c location of maximum camber
- 16-series like NACA 16-412, based on
o maximum thickness and its x/c location
o design lift coefficient
- 6-series like NACA 64-412, based on
o maximum thickness and its x/c location
o design lift coefficient
o camber line specification (a=...)
o A-Modification
Note that the 6-series sections are only approximated, for correct data use
tables from NACA reports.
- TsAGI "B" series airfoils, based on
o maximum thickness
- NPL "EC", "EH" series airfoils, based on
o maximum thickness and its x/c location
o maximum camber and its x/c location
- symmetrical Circular Arc airfoils, based on
o maximum thickness
- symmetrical Double Wedge airfoils, based on
o maximum thickness and its x/c location
- Cambered Plate airfoils, based on
o maximum thickness
o maximum camber and its x/c location
- Van de Vooren conformal mapping airfoil, based on
o maximum thickness
o trailing edge angle
- Newman airfoil, based on
o maximum thickness
- Helmbold-Keine airfoil, based on
o maximum thickness and its x/c location
o trailing edge angel
o radius at center of airfoil
o leading edge radius
Note that only a subset of possible parameters leads to reasonable shapes.
Note about all x-y graphs:
1. You can copy most of the graphs to the clipboard or save them to a file in AutoCad
DXF, Adobe Illustrator encapsulated Postscript or SVG vector graphics format by
pressing the right mouse button while the mouse pointer is located over the graph
(context menu).
Paste the test into a text editor, save as filename.dxf, .eps/.ai or .svg and import into
your favorite graphics, CAD or word processing program.
2. If you have allowed the applet to write files, or if you run JavaFoil standalone, you can
also export the graphics directly to a file, without using the clipboard.
3. Additionally, you can import numerical data into a graph to enhance it. For example you
can import test results into the the polar plot for comparison. Import works only, if the
graph already contains some data, so that the axis system is already set up. The data is
lost when the graph is cleared for a new calculation, though.
MODIFYCARD

Here you find tools to modify the geometry of the current airfoil.


You can create a new distribution of points, change the camber and the thickness or deflect a
plain flap. A scaling and rotating option is also available to transform the airfoil. The result of any
modification is shown at the bottom of the card and is also transferred back to the Geometry
card. Each transformation is executed when you press the button to the left of the
corresponding text entry field.
Before a modification is performed JavaFoil saves a copy of the current geometry on top of a
stack and you can either use the "Undo" button to back up to the previous configuration or
select the desired configuration from the combo box. It is recommended to enter a name which
reflects your modification before pushing the "Modify" button (e.g. "NACA 0010 / F+10" for a 10
flap deflection). To reduce memory overhead, the number of undo steps is currently limited to
10. Each additional modification will drop one saved geometry from the bottom of the stack.
Changing the Number of Points
You can increase or decrease the number of points of the selected airfoil element. If
you specify a negative number of points, JavaFoil uses a special tight spline curve
which produces an almost linear interpolation between the given coordinate points.

Deflection of a simple Flap
Uses the given a flap chord and deflection angle (positive = trailing edge down) to rotate
the coordinate points covered by the flap. No points are added to the convex side of the
airfoil, so any initially coarse distribution might lead to a poor representation of the
flapped airfoil. Also some points on the concave side are deleted if they would move
inside of the flapped airfoil.

Modification of Thickness and Camber
Creates a new shape based on the camber distribution and the thickness distribution of
the current airfoil. Both distributions are derived by analysis of the y-coordinates for a
given x station - the thickness distribution is not applied at right angles to the camber
line. This will introduce a small error where the camber line has a larger inclination. The
extracted distributions are then scaled with the ratio of prescribed to current maximum
value. The x-axis is defined to be the straight line from the trailing edge to the
coordinate point having the largest distance to the trailing edge.

Scaling
Scaling is always performed with respect to the origin of the coordinate system (0/0). X
and Y coordinates are multiplied with the given factor.

Trailing Edge Gap
The upper and lower surface of the airfoil are rotated around the leading edge point to
open or close the trailing edge.

Rotation and Translation x/y
Allows arbitrary modifications to the position of the airfoil.

Duplicate
Creates a copy of the selected airfoil element. You have to move and/or rotate it
afterwards so that it does not overlap its parent.

Delete
Deletes the selected airfoil element.

Flip Y
Mirrors the selected airfoil element along the x-axis (upside down transformation).

Copy (Text)
Copies a text buffer to the clipboard if your system security settings allow this. The
buffer contains the coordinates of the camber distribution (x
c
/y
c
) and the coordinates of
the thickness distribution (x
t
/y
t
) for the current airfoil.
Notes:
1. Changes to the camber are performed by scaling an existing camber line. Thus it is not
possible to camber a symmetrical airfoil, as this would additionally require the
specification of a camber line shape.
2. Changes act on the selected airfoil element(s) (in case of multi-element airfoils).
DESIGNCARD
This card can be used to design an airfoil based on a prescribed target pressure (coefficient)
distribution. Such a method is called "inverse design" - the geometry of the airfoil is a result of
the given pressure distribution.

Pressing the "Setup" button initializes the design procedure by copying the current airfoil. This
airfoil is then analyzed at the given angle of attack and produces an initial target pressure
distribution.
Now you can modify this pressure distribution either by using a smooth distortion of the target
Cp or in a single point mode. Simply grab a point on the target distribution and drag it up or
down to modify the curve.
After you are satisfied with the target Cp-distribution, you can run several design cycles and
check the result. The current solution is overlaid on the previous results to allow for
convergence check.
You can "Redraw" the screen to clear the intermediate results.
Instead of viewing the Cp-distribution in the usual way Cp=f(x/c) you can "unfold" the distribution
by plotting Cp=f(s). Here s is the arc length measured along the airfoil surface (the order is
upper surface, nose, lower surface). This representation makes it easier to modify the leading
edge region. When this viewing mode is active, a slider can be used to enlarge the leading
edge or trailing edge regions.
A relaxation factor helps to stabilize the procedure and to smooth the geometry changes -
usually a value between 10% and 25% is sufficient.
Typically, the number of design steps should be between 10 and 50 steps. You can repeat the
design until JavaFoil has reached the target.
The analysis takes into account the effects of ground proximity as well as multiple
elements. Note that the design procedure only acts on the first element of multiple elements.
You have to change the order on the Geometry card if you want to design another element.
If you want to define the target pressure distribution by numbers, you can use the "Details..."
button to open a window where you can enter or paste x/c, y/c, Cp triples. From these
data,JavaFoil only reads the Cp values, changing x/c or y/c will have no effect.
Notes:
1. When "symmetric Cp modification" is checked, the Cp curves for upper and lower
surfaces are modified in the same manner. This results in a global modification of the
airfoil thickness.
2. When "anti-symmetric Cp modification" is checked, the Cp curves for upper and lower
surfaces are modified with reversed signs. This results in a global modification of the
airfoil camber.
3. When none of both options is checked, you can modify the Cp values for single points,
which is useful for smoothing wavy airfoils.
4. The method is not foolproof and may diverge if the changes are too large. Problems
may occur if modifications close to the leading and trailing edges are performed. This is
because the stagnation point region is very sensitive to small changes in overall
circulation.
5. The tool tries to translate any arbitrary Cp-distribution into a shape. If the target
distribution is "ill posed", it may happen that no shape exists which creates such a
distribution!
6. Markers on the airfoil shape indicate the results of the boundary layer analysis: green:
transition, red: separation.
VELOCITYCARD
This card can be used to calculate the velocity distribution on the surface of the airfoil for
several angles of attack. As usual, all angles are counted from the x-axis of the airfoil.
The graph shows the velocity on both sides of the airfoil and can be used to smooth airfoils. If
the velocity distributions show wiggles and zig-zag waves, any subsequent boundary layer
analysis on the Polars card will probably create unrealistic results. To smooth an airfoil, go back
to the Geometry card and change single y-coordinate values or use the Design card to modify
the velocity distribution directly. Then re-analyze. Yes, this is slow, but possible and probably
better than an automatic global smoother which would smooth over the whole airfoil.

You can display either the distribution of the local velocity v/V or the resulting local pressure
coefficient Cp. If a Mach number other than zero is specified on the Options card, the critical
velocity ratio V* respectively the critical pressure coefficient Cp* is plotted also. The distributions
are corrected for compressibility effects by the Karman-Tsien rule, but one must be aware of the
fact that such corrections are valid only for Mach numbers below approximately 0.7.
Notes:
1. The Mach number from the Options card will be used for a compressibility correction on
all cards. The corrections work reasonably well for Mach numbers below 0.75. If
supersonic flow occurs somewhere on the airfoil surface, the code will not fail, but you
would need a different type of analysis code for transonic and supersonic flow.
2. The table in the upper right of the card lists force and moment coefficients which are
determined for the Reynolds number taken from the boundary layer card. However, the
velocity distribution shown is for the inviscid flow. i.e. will not be affected by Reynolds
number.
3. The table also shows the critical Mach number for each angle of attack. If the onset flow
exceeds this Mach number, supersonic flow will occur somewhere on the airfoil. This
Mach number is also known as the "Critical Mach number". As usual you can copy it
using the context menu.
FLOWFIELD CARD
If you want to get an impression of how the flow around the airfoil looks like, this card is for you.
The panel analysis method works with the surface of an airfoil only, but when the surface
velocity has been determined, potential flow theory can be used to calculate the flow velocity
and direction anywhere in space.

You can specify a regular x-y grid and an angle of attack. After solving for the surface velocity
distribution, an evaluation is performed for each point on the grid.
The buttons in the control bar at the bottom of the card perform the following actions:
Analyze It!
performs the analysis of the flow field at all grid points plus any additional post-
processing like streamlines or iso-Cp lines.

Print...
Sends a copy of the picture to the printer, if your system security settings allow this.

Save...
Creates a text buffer which contains the coordinates and the velocity vector for each
grid point (x, y, v
x
, v
y
and v). This text buffer is saved to a file if your system security
settings allow this. The format is suitable for the commercial plotting
programTecplot (by Amtec Corp.).

Copy (Text)
Like the Save... command, but the text buffer is copied to the clipboard if your system
security settings allow this.

You can select from the following display options:
- tufts ... shows a black "tuft" at each grid point, which is aligned with the local flow
direction.
- colored cells ... colors a rectangle around each grid point according to the local
pressure. This makes for a nice picture of the pressure field around the airfoil and
shows how far the pressure is changed due to the airfoil.
- iso Cp ... plots lines of constant pressure, like altitude lines on a map.
- Cp-Vectors ... plots pressure vectors normal to the airfoil surface.
- stream lines ... black stream lines are drawn starting at the left border. The result
resembles the injection of smoke into a wind tunnel.
- timed stream lines ... the stream lines are dashed to show the distance traveled during
equal time intervals. Models something like a pulsed smoke generator. You will notice
that particles arriving side by side at the leading edge will NOT meet again at the trailing
edge.
- higher accuracy ... this option controls the calculation of the stream lines only. When
selected, a 4th order Runge-Kutta scheme and a smaller time step are used, otherwise
simple forward differencing is employed. The results are more accurate in regions
where the streamlines are highly curved - calculation time is increased by a factor of 5,
though.
Notes:
1. As each evaluation can be time consuming, it is recommended to start with the coarse
default grid and refine depending on the power of your computer and the amount of
time you wish to spend.
2. The rectangular grid is simply placed on top of the airfoil. It is not adapted to the airfoil
shape. While interior points are excluded in the calculation, a coarse grid will create a
rough outline of the airfoil only.
3. The graph shows, how the airfoil affects the whole flow field, not only the immediate
neighborhood of the section.
BOUNDARYLAYERCARD
This card is a source of information for the experts. It shows all important boundary layer
parameters like thickness and shape functions. Additional parameters are available in the
listings. A fixed transition location can be defined on the Polars card (see below).

POLARSCARD
When you have created or imported a sufficient smooth airfoil shape, you can calculate lift and
drag on this card.

After specification of the desired Reynolds number and angle of attack range as well as
selecting a surface roughness you can start the analysis. For each Reynolds number / angle of
attack condition, JavaFoil will first calculate the velocity distribution and then perform a
boundary layer analysis. The resulting lift, drag and moment coefficients as well as location of
transition and separation will be presented in graphs and tables. A transition strip can be
simulated by specifying a transition location x/c for both sides of the airfoil (the default setting of
100% corresponds to natural transition).

REFERENCIAS

1. Sailing Yatch Design. Practice/Theory. Edited by A. Claughton, J. Wellicome & A.


Shenoi. Addison Wesley Longman, Lmtd. 1998.
2. I.H Abbot y A.E. Doenhoff. Theory of Wing Sections. McGraw-Hill Book Company,
Inc. 1949.
3. R. A. Royce. A Rational Prismatic Hull Approach for Planing Hull Analysis. SNAME
Meeting. Cleveland (Ohio). January 1994.
4. F. M. White. Mecnica de Fluidos. Mc. Graw Hill (2008)

También podría gustarte