Capítulo 3

(Baseado no cap. 9 do livro de Fox e McDonald, 9ª ed.)

Escoamento Externo Viscoso Incompressível

Escoamentos externos são escoamentos sobre corpos imersos em fluido sem fronteiras.
Os escoamentos sobre placa plana semi-infinita e sobre cilindro são exemplos de escoamentos
externos. O objetivo neste capítulo é quantificar o comportamento dos fluidos viscosos e
incompressíveis em escoamentos internos.
Diversos fenômenos que ocorrem no escoamento externo sobre corpo imerso são
ilustrados no esboço do escoamento viscoso com alto número de Reynolds sobre aerofólio, Fig.
3.1. O escoamento na camada limite é inicialmente laminar (CLL). A transição para escoamento
turbulento ocorre a alguma distância do ponto de estagnação, dependendo das condições da
corrente livre, rugosidade da superfície e gradiente de pressão (CLT). Os pontos de transição
estão indicados por “T”. Na região de pressão crescente (com gradiente adverso de pressão) a
separação do escoamento poderá ocorrer. O fluido que estava nas camadas limites da superfície
do corpo forma a esteira viscosa atrás dos pontos de separação.
Existe solução exata para o escoamento na camada limite laminar, mas para a camada
limite turbulenta, somente é possível solução aproximada.
O aerofólio da Fig. 3.1 é submetido a força resultante oriunda das forças de cisalhamento
e de pressão que atuam nas superfícies do corpo. A componente da força resultante
perpendicular ao escoamento uniforme a montante, U  , é chamada de sustentação, enquanto a

componente paralela é denominada arrasto. A presença de separação do escoamento impede a

determinação analítica da sustentação e do arrasto.

Fig.3.1 Detalhes do padrão de escoamento viscoso em torno de aerofólio

Prof. Dr. Admilson T. Franco 57

UTFPR - Curso de Engenharia Mecânica
Cap. 3 - Mecânica dos Fluidos 2
 PARTE A
CAMADAS LIMITES

3.1 Conceito de Camada Limite

O conceito de camada limite foi introduzido primeiro por Ludwig Prandtl, alemão estudioso
de aerodinâmica, em 1904, conceito que marcou o começo da era moderna da mecânica dos
fluidos.
Prandtl mostrou que muitos escoamentos viscosos podem ser analisados dividindo-os em
duas regiões, uma perto das fronteiras sólidas (a camada limite), na qual o efeito da viscosidade
é importante, e outra cobrindo o restante do escoamento (escoamento externo), para a qual o
fluido pode ser tratado como invíscido.
Na camada limite, tanto as forças viscosas quanto às de inércia são importantes. Por isso,
para facilitar a caracterização de escoamentos na camada limite, define-se o número de Reynolds
( Re x ) como a razão entre as forças de inércia e as forças viscosas:

U x U x
Re x = = (3.1)
 
Sendo:
U a velocidade do fluido distante da parede ( U  ); [m/s]
 a massa específica do fluido; [kg/m3]
 a viscosidade dinâmica; [Pa.s]
 a viscosidade cinemática (  /  ); [m2/s]
x o comprimento característico [m], dependente da geometria do escoamento, Fig. 3.2.

Sob condições típicas de escoamento, a camada limite apresenta escoamento laminar

desde o início da placa até a região de transição, a partir da qual, o escoamento passa a ser
turbulento. Entre os fatores que afetam a transição da camada limite, está o gradiente de pressão,
a rugosidade superficial, a transferência de calor, as forças de campo e perturbações da corrente
livre. Em geral, a transição ao longo da placa plana correspondente a número de Reynolds de
500.000.
Em muitas situações reais, a camada limite desenvolve-se sobre superfícies amplas
essencialmente planas. Os exemplos incluem escoamentos sobre cascos de navios e de
submarinos, asas de aviões e movimentos atmosféricos sobre terreno plano.
A Fig. 3.2 ilustra o crescimento da camada limite sobre placa plana. Pode-se verificar o
crescimento mais acelerado da espessura da camada limite na região turbulenta que na laminar.

Prof. Dr. Admilson T. Franco 58

UTFPR - Curso de Engenharia Mecânica
Cap. 3 - Mecânica dos Fluidos 2
 T
y
L
x

Fig. 3.2 Camada limite sobre placa plana (espessura vertical exageradamente ampliada)

Espessuras de Camada Limite

A camada limite é a região adjacente à superfície sólida na qual as forças viscosas são
importantes. Três espessuras são definidas, conforme Fig. 3.3.
i. A espessura de perturbação ou simplesmente espessura, , da camada limite, Fig. 3.3b -
distância da superfície ao ponto em que a velocidade é 99% da velocidade da corrente
livre, u/U=0,99.
ii. A espessura de deslocamento, *, Fig. 3.3b - distância pela qual a fronteira sólida teria
que ser deslocada com escoamento sem atrito para fornecer o mesmo déficit de vazão
em massa que existe na camada limite.
Para escoamento incompressível,  =constante e

U  =   (U − u ) dy
*

 
 u   u 
 =  1 −
*
 dy   1 −  dy [m] (3.2)
0  0 
U U

iii. A espessura da quantidade de movimento,  , Fig. 3.3c - espessura da camada de fluido,

com velocidade U , para qual o fluxo de quantidade de movimento é igual ao déficit do
fluxo de quantidade de movimento através da camada limite. Para  =constante

U 2 =   u (U − u ) dy
0

e
 
u  u  u  u 
 = 1 −  dy   1 −  dy [m] (3.3)
0
U  U  0
U  U 

 * e  são espessuras integrais, cujos integrandos desaparecem na corrente livre, pois uU, e
são mais fáceis de avaliar com precisão a partir de dados experimentais, do que a espessura de
perturbação,  , da camada limite.

Prof. Dr. Admilson T. Franco 59

UTFPR - Curso de Engenharia Mecânica
Cap. 3 - Mecânica dos Fluidos 2
 As hipóteses simplificadoras usualmente feitas em análises de engenharia para o
desenvolvimento da camada limite são:
1. u → U para y = 

2. u y → 0 para y = 

3. v  U dentro da camada limite

Os resultados das análises desenvolvidas nas duas próximas seções mostram que a
camada limite é muito delgada comparada com o seu comprimento de desenvolvimento ao longo
da superfície, /L<<1. Então, é também razoável supor que:
4. A variação de pressão através da camada limite delgada é desprezível. A distribuição
de pressão da corrente livre é impressa sobre a camada limite.

Fig. 3.3 Definições das espessuras da camada limite: a) espessura de deslocamento; b)

espessura de perturbação e c) espessura de quantidade de movimento

Exemplo 3.1 – Escoamento de Camada Limite em Canal

O túnel de vento de laboratório possui seção de teste quadrada, com 305 mm de lado. Os
perfis de velocidade de camada-limite são medidos em duas seções transversais e as espessuras
de deslocamento são avaliadas a partir dos perfis medidos. Na seção 1, onde a velocidade da
corrente livre é U1= 26 m/s, a espessura de deslocamento é 1*=1,5 mm. Na seção 2, localizada
a jusante da seção 1, 2*=2,1 mm. Calcule a variação na pressão estática entre as seções 1 e 2.
Expresse o resultado como fração da pressão dinâmica da corrente livre na seção 1. Considere
a atmosfera na condição padrão.

Prof. Dr. Admilson T. Franco 60

UTFPR - Curso de Engenharia Mecânica
Cap. 3 - Mecânica dos Fluidos 2
3.2 Camada Limite Laminar em Placa Plana: Solução Exata
A solução para a camada limite laminar em placa plana horizontal foi obtida por H. Blasius,
aluno de Prandtl, em 1908. Para escoamento bidimensional, permanente, incompressível, com
gradiente de pressão nulo, as equações que governam o movimento reduzem-se a
u v
+ =0 (3.4)
x y

u u  2u
u +v =v 2 Q.M. em x (3.5)
x y y
com condições de contorno na fronteira, Fig. 3.3:
para y = 0, u = 0, v=0
u
para y→ u =U, =0 (3.6)
y
Blasius argumentou que o perfil de velocidade adimensional, u / U , deveria ser similar
para todos os valores de x , quando traçado contra a distância adimensional em relação à parede;
a espessura de camada limite,  , é escolha natural para tornar adimensional a distância da
parede. Então a solução proposta tem a forma
u y
= g ( ) onde  (3.7)
U 
Fundamentado na solução de Stokes, Blasius concluiu que   vx U e estabeleceu

U
=y (3.8)
vx
Introduzindo a função de corrente,  ,onde

 
u= e v=− (5.4)
y x
A definição de função de corrente adimensional é expressa como:

f ( ) = (3.9)
vU
Faz-se f ( ) a variável dependente e  a variável independente na Eq. (3.5).
As componentes da velocidade são dadas por

   df U df
u= = = vxU =U (3.10)
y  y d vx d

  f 1 vU   df  1 1  1 vU 
v=− = −  vxU + f  = −  vxU −  + f
x  x 2 x   d  2 x  2 x 
1 vU  df 
v=  −
2 x  d
f (3.11)

Prof. Dr. Admilson T. Franco 61
UTFPR - Curso de Engenharia Mecânica
Cap. 3 - Mecânica dos Fluidos 2
 Derivando as componentes da velocidade, também pode ser mostrado que

u U d2 f u d2 f  2u U 2 d 3 f
=−  2 ; = U U vx ; =
x 2 x d y d 2 y 2 vx d 3
Substituindo essas expressões na Eq. (3.5)

d3 f d2 f
2 + f =0 (3.12)
d 3 d 2
com condições de fronteira:
df
para  = 0, f = =0
d
df
para  → , =1 (3.13)
d
As equações diferenciais parciais de segunda ordem que preveem o crescimento da
camada limite em placa plana (Eqs. (3.4) e (3.5)) foram transformadas na equação diferencial
não linear de terceira ordem (Eq. (3.12)) com condições de fronteira dadas pela Eq. (3.13).

Os valores numéricos de f , df d e d f d na Tabela 3.1 foram calculados com

2 2

microcomputador, usando integração numérica de Runge-Kutta de 4ª ordem.

Tabela 3.1 A função f ( ) para a camada limite laminar em placa plana com

ângulo de incidência zero

U u
=y f f '= f ''
vx U

0 0 0 0,3321
0,5 0,0415 0,1659 0,3309
1,0 0,1656 0,3298 0,3230
1,5 0,3701 0,4868 0,3026
2,0 0,6500 0,6298 0,2668
2,5 0,9963 0,7513 0,2174
3,0 1,3968 0,8460 0,1614
3,5 1,8377 0,9130 0,1078
4,0 2,3057 0,9555 0,0064
4,5 2,7901 0,9795 0,0340
5,0 3,2833 0,9915 0,0159
5,5 3,7806 0,9969 0,0066
6,0 4,2796 0,9990 0,0024
6,5 4,7793 0,9997 0,0008
7,0 5,2792 0,9999 0,0002
7,5 5,7792 1,0000 0,0001
8,0 6,2792 1,0000 0,0000

Prof. Dr. Admilson T. Franco 62

UTFPR - Curso de Engenharia Mecânica
Cap. 3 - Mecânica dos Fluidos 2
 Da Tabela 3.1, verifica-se que para  = 5,0, u U = 0,992 . Com a espessura de camada
limite,  , é definida como o valor de y para o qual u U = 0,99 , a Eq. (3.8) resulta em

5, 0 5, 0 x  ( x) 5, 0
 = → = (3.14)
U vx Re x x Re x

A tensão de cisalhamento na parede pode ser expressa como

u  d2 f 
w =  =  
y  y =0
U U vx
d 2  =0

Então
0,332 U 2
 w = 0,332U U x = (3.15)
Re x
e o coeficiente de tensão de cisalhamento na parede
w 0, 664
Cf = = (3.16)
1 Re x
U 2
2
1/ 2
A espessura de camada limite aumenta segundo x , e a tensão de cisalhamento na
1/ 2
parede e o coeficiente de atrito superficial variam de acordo com 1/ x .

3.3 A Equação Integral da Quantidade de Movimento

A solução exata de Blasius fornece expressão para a espessura da camada limite laminar,
 ( x ) , e para a tensão de cisalhamento na parede,  w ( x ) . A solução analítica fechada para o
perfil de velocidade não era possível, então a solução numérica foi necessária.
Como soluções exatas para camadas limites turbulentas não existem, técnicas de
soluções aproximadas são necessárias. O objetivo é desenvolver uma equação que preveja a
como a camada limite cresce em função da distância ao longo do corpo. Deduz-se uma relação
que pode ser aplicada tanto ao escoamento laminar quanto ao turbulento; além disso, a relação
não é restrita aos escoamentos com gradientes de pressão nulos.
Considere o escoamento permanente, incompressível, bidimensional, sobre superfície
sólida. A espessura da camada limite,  , cresce dependente do aumento da distância x . Para a
análise, escolhe-se o volume de controle diferencial, de comprimento dx , largura w e altura

 ( x ) , conforme mostrado na Fig. 3.4. A velocidade da corrente livre é U ( x ) .

Prof. Dr. Admilson T. Franco 63

UTFPR - Curso de Engenharia Mecânica
Cap. 3 - Mecânica dos Fluidos 2
 Fig. 3.4 Volume de controle diferencial em camada limite.

a. Equação da Conservação da Massa

Aplicando a Equação da Conservação da Massa ao VC da Fig. 3.4
= 0 (1)
 (3.12)
0=   d  + SC V dA
t VC
Considerações: (1) Escoamento permanente
(2) Escoamento bidimensional

Portanto: 0=  V
SC
dA → 0 = mab + mbc + mcd

ou mbc = −mcd − mab

Avaliando esses termos para o volume de controle de largura w .

(Note que a componente velocidade, u, e a espessura da camada-limite, δ, o limite superior na integral, dependem de x.)
Prof. Dr. Admilson T. Franco 64
UTFPR - Curso de Engenharia Mecânica
Cap. 3 - Mecânica dos Fluidos 2
   
mbc = −     u dy  dx  w
 x  0  

b. Equação de Balanço de Quantidade de Movimento

Aplicando a componente x da equação da quantidade de movimento ao volume de
controle abcd

= 0 ( 3) = 0 (1)
 (4.18a)
FS x + FBx =  u  d  + SC u V dA
t VC

Suposição: (3) FBx = 0

Então: FS x = + mf ab + mf bc + mf cd

sendo mf a componente x do fluxo da quantidade de movimento.

O fluxo de quantidade de movimento na direção x , através da superfície de controle,

    
 u V dA = −    u  u dy  dx − U    u dy  dx  w

 x  0 
 
x  0 
 
SC

Agora que se tem expressão adequada para o fluxo de quantidade de movimento na

direção x , através da superfície de controle, consideremos as forças superficiais que atuam
sobre o volume de controle na direção x . Reconhecemos que as forças normais com

Prof. Dr. Admilson T. Franco 65

UTFPR - Curso de Engenharia Mecânica
Cap. 3 - Mecânica dos Fluidos 2
componentes não nulas na direção x atuam sobre as três superfícies de controle. Além disso,
uma força cisalhante atua sobre a superfície ad . Como o gradiente de velocidade tende a zero
na borda da camada limite, a força de cisalhamento atuando ao longo de bc é desprezível.

Fig. 3.5 Volume de controle diferencial.

Somando as componentes na direção x de todas as forças atuando sobre o volume de

controle,

Prof. Dr. Admilson T. Franco 66

UTFPR - Curso de Engenharia Mecânica
Cap. 3 - Mecânica dos Fluidos 2
 0 0
 dp 1 dp 1   dp 
FSx = −  dx − dx d −  w dx − d wdx  w = −  dx −  wdx  w
 dx 2 dx 2   dx 
Substituindo as expressões para  u V
SC
dA e FS x na equação de balanço da

quantidade de movimento em x ,

 dp        

 
 −  dx −  dx  w =   u  udy  dx − U    udy  dx  w
 x  0 x  0
w
 dx    
e dividindo a equação por wdx, resulta
 
dp  
−  w =  u  udy − U
x 0
−  udy
dx x 0

O gradiente de pressão, dp dx pode ser determinado pela aplicação da equação de

Bernoulli ao escoamento não viscoso externo a camada limite; resultando em:
dp dx = − UdU dx .

Usando as definições de espessura de deslocamento,  , (Eq. (3.2)), e espessura de

*

quantidade de movimento,  , (Eq. 3.3), obtém-se

w d
= (U 2 ) +  *U
dU
(3.17)
 dx dx
A Eq. (3.17) é a equação integral de balanço da quantidade de movimento, representando
a equação diferencial ordinária para a espessura de camada limite, desde que seja admitida forma
adequada para o perfil de velocidade e que a tensão de cisalhamento na parede possa ser
relacionada com as outras variáveis. Determinada a espessura da camada limite, as espessuras
de quantidade de movimento e deslocamento, e a tensão de cisalhamento na parede podem ser
calculadas.
A Eq. (3.17) fica restrita a escoamento permanente, incompressível e bidimensional sem
a presença de forças de campo paralelas à superfície.
Para fazer uso da equação para estimar a espessura de camada limite como função de
x , deve-se:
1. Propor uma primeira aproximação para a distribuição de velocidade de corrente livre, U ( x ) , o quê pode
ser feita a partir da teoria de escoamento não viscoso. A pressão na camada limite é relacionada com a

velocidade de corrente livre U ( x ) , através da equação de Bernoulli.

2. Admitir forma razoável para o perfil de velocidade dentro da camada limite.
3. Relacionar a tensão de cisalhamento na parede com o campo de velocidade.

Prof. Dr. Admilson T. Franco 67

UTFPR - Curso de Engenharia Mecânica
Cap. 3 - Mecânica dos Fluidos 2
3.4 Uso da Equação Integral da Quantidade de Movimento para Escoamento com Gradiente
de Pressão Nulo
Para o caso especial de placa plana com gradiente de pressão zero, a pressão p e a
velocidade U da corrente livre são ambas constantes, de modo que U(x) = U = constante →
dp dx = 0 .
integral da quantidade de movimento reduz-se então a

d d u  u 
dx 0 U
 w = U 2 = U 2 1 −  dy (3.18)
dx  U 
Note que u U é adimensional e  é função de x apenas. Consequentemente, é

conveniente mudar a variável de integração de y para y 

y
= e dy =  d

e a equação integral da quantidade de movimento para gradiente de pressão zero se torna

d d u  
1
u
dx 0 U
 w = U = U 2
2
1 −  d (3.19)
dx  U 
Para se resolver esta equação para a espessura de camada limite com função de x :

1. Admitir uma distribuição de velocidade na camada limite - qualquer relação funcional da forma

u  y
= f   = f ( )
U  
(a) A distribuição de velocidade admitida precisaria satisfazer certas condições físicas de contorno:
Para y = 0, u=0
Para y =, u =U
u
Para y =  , =0
y
(b) Note que, admitida a distribuição da velocidade, o valor numérico da integral na Eq. (3.19) é;
1
u u 
 U 1 − U  d = 
0
= constante = 

e a equação integral da quantidade de movimento torna-se

d
 w = U 2 
dx
2. Obter expressão para w em termos de  , permitindo resolver para  ( x ) , como ilustrado a seguir.

Prof. Dr. Admilson T. Franco 68

UTFPR - Curso de Engenharia Mecânica
Cap. 3 - Mecânica dos Fluidos 2
3.5.1 Escoamento Laminar
Para escoamento laminar sobre placa plana, suposição razoável para o perfil de

velocidade é o polinômio em y : u = a + by + cy 2
As condições de contorno na camada limite são:
Para y = 0, u=0
Para y =  , u =U
u
Para y =  , =0
y
Avaliando as constantes a , b e c resulta
2
u  y  y
= 2   −   = 2 −  2 (3.20)
U    
A tensão de cisalhamento na parede é dada por

u 
w =  
y  y =0

Substituindo o perfil de velocidade admitido, Eq. (3.20), na expressão para  w , resulta

u  U (u / U ) U d (u / U )
w =   = =
y  y =0  ( y /  ) y /  =0
 d  =0

ou

U d (2 −  2 ) U 2U
w = = d (2 − 2 ) =
 d  =0
  =0 
Agora pode-se de aplicar a equação integral da quantidade de movimento

d u  
1
u
 w = U 2  1 −
dx 0 U  U
 d

(3.19)

Substituindo  w e u U ,obtém-se

15
 d = dx
U
que é a equação diferencial para  . Integrando novamente resulta

2 15
= x+c
2 U
Se for admitido que  = 0 para x = 0 , o que é equivalente a admitir camada limite laminar a partir
do bordo de ataque da placa, então c = 0 , assim

30 x  ( x) 30  5, 48
= ou = = (3.21)
U x Ux Re x

Prof. Dr. Admilson T. Franco 69

UTFPR - Curso de Engenharia Mecânica
Cap. 3 - Mecânica dos Fluidos 2
 É notável constatar que a Eq. (3.21) apresenta erro de somente em 10% (a constante é
maior que 5,0) em comparação com a solução exata (Seção 3.3). A Tabela 3.2 resume resultados
correspondentes calculados com o uso dos perfis de velocidade aproximados e lista os resultados
obtidos da solução exata.
Conhecida a espessura de camada limite, todos os detalhes do escoamento podem ser
determinados. O coeficiente de tensão de cisalhamento na parede, ou de “atrito superficial”, é
definido como:
w
Cf = (3.22)
1
U 2
2
w 2 (U  ) 4 x 1 Re x
e Cf = = = =4
1
U 2
1
U 2 Ux  Re x 5, 48
2 2
Finalmente,
0, 730
Cf = (3.23)
Re x

Como a variação de  w é conhecida, o arrasto viscoso sobre a superfície pode ser

avaliado por integração sobre a área da placa plana.

Exemplo 3.2- Camada Limite Laminar sobre Placa Plana: Solução aproximada usando
Perfil de Velocidade Senoidal
O problema ilustra a aplicação da forma integral da equação de balanço de quantidade de
movimento para escoamento de camada limite laminar sobre placa plana.
Encontre expressões para:
(a) A taxa de crescimento de δ como função de x.
(b) A espessura de deslocamento, δ*, como função de x.
(c) A força de atrito total sobre a placa de comprimento L e largura b.

u  y 
= sin  
U 2

Prof. Dr. Admilson T. Franco 70

UTFPR - Curso de Engenharia Mecânica
Cap. 3 - Mecânica dos Fluidos 2
 Tabela 3.2 Resultados do Cálculo do Escoamento de Camada Limite Laminar sobre Placa
Plana com Ângulo de Incidência Zero Fundamentado em Perfis de Velocidade Aproximados

3.5.2 Escoamento Turbulento

Detalhes do perfil de velocidade turbulento em camadas limites com gradiente nulo de
pressão são muito semelhantes à aqueles para escoamento turbulento em tubos e canais.
Para escoamento turbulento de camada limite, adapta-se a expressão desenvolvida para
escoamento em tubos
0,25
 v 
 w = 0, 0332 V 2   (2.39)
 RV 
Para perfil de potência 1 7 em tubo, a Eq. (2.24) fornece V U = 0,817 . Substituindo para

V = 0,817U e R =  na Eq. (2.39),

1
 v  4
 w = 0, 0233U  
2
(3.25)
 U 
Agora pode-se aplicar a equação integral de balanço de quantidade de movimento

d u  u 
1

dx 0 U  U 
 w = U 2 1 −  d (3.19)

Substituindo para  w e u U e integrando,

d 1 7 7 d
1/ 4
 v 
1
0, 0233   =   (1 − 1 7 ) d =
 U  dx 0 72 dx

Desse modo, obtém-se a equação diferencial para  :

1
1
v 4
 d = 0, 240   dx
4

U 
1/ 4
4 54 v
Integrando resulta em:  = 0, 240   x+c
5 U 
Prof. Dr. Admilson T. Franco 71
UTFPR - Curso de Engenharia Mecânica
Cap. 3 - Mecânica dos Fluidos 2
 Se for admitido que   0 para x = 0 (isso é equivalente a admitir escoamento turbulento
a partir da borda de ataque), então c = 0 e
15
 v  4/ 5
 = 0,382   x
U 
ou
1/5
 ( x)  v  0,382
= 0,382   = (3.26)
x U x  Re1/5
x

Usando a Eq. (3.25), obtém-se o coeficiente de atrito superficial em termos de  :

w
1/ 4
 v 
Cf = = 0, 0466  
1
U 2  U 
2
Substituindo para  ,
w 0, 0594
Cf = = (3.27)
1
U 2 Re1/x 5
2
Experimentos mostram que a Eq. (3.27) prediz muito bem o atrito superficial turbulento

em placa plana para 5  10  Re x  10 . A concordância é notável tendo em vista a natureza

5 7

aproximada da análise .
A tensão de cisalhamento na parede é muito maior na camada limite turbulenta que na
camada limite laminar, razão primária para o desenvolvimento mais rápido das camadas limites
turbulentas.

Exemplo 3.3 – Camada Limite Turbulenta sobre Placa Plana: Solução Aproximada com
Perfil de Velocidade de Potência 1/7.
Água escoa a U = 1 m/s sobre placa plana, com L = 1 m na direção do escoamento. A
camada-limite se torna turbulenta na borda de ataque. Avalie a espessura de perturbação, δ, a
espessura de deslocamento, δ*, e a tensão de cisalhamento de parede, w, para x = L. Compare
com os resultados nessa posição para o escoamento mantido laminar. Considere o perfil de lei
de potência 1/7 para a velocidade na camada-limite turbulenta.

Prof. Dr. Admilson T. Franco 72

UTFPR - Curso de Engenharia Mecânica
Cap. 3 - Mecânica dos Fluidos 2
3.5 Gradientes de Pressão no Escoamento de Camada Limite
A equação integral de balanço de quantidade de movimento para a qual o gradiente de
pressão é zero

d d u  u 
dx 0 U
 w = U 2 = U 2 1 − dy (3.19)
dx  U 
indica que a tensão de cisalhamento na parede é equilibrada pelo decréscimo na quantidade de
movimento do fluido. A espessura da camada limite continua a crescer e o fluido perto da parede
fica sendo continuamente retardado (perdendo quantidade de movimento). “Será o fluido próximo
a parede eventualmente levado ao repouso?”. Ou seja, “Para o caso em que dp dx = 0 é possível

que u y ) y =0 = 0 ?”

Ao considerar as distribuições de tensão de cisalhamento para as placas planas, verifica-

se que para escoamento laminar
 w ( x ) constante
=
U 2 Re x
e para escoamento turbulento
 w ( x ) constante
=
U 2 Re1/x 5

Lembrando que  w = u / y y =0

pode-se então dizer que para qualquer comprimento

finito de placa u / y y =0
nunca será zero.

O ponto na fronteira sólida em que u / y y =0

= 0 é definido como o ponto de separação.

Consequentemente, para dp / dx = 0 , nunca ocorrerá separação do escoamento; a camada de

fluido na vizinhança da superfície sólida não pode ser trazida à velocidade zero.
O gradiente de pressão é dito ser adverso se a pressão aumenta no sentido de
escoamento, se ( dp / dx  0) . Quando dp / dx  0 , a pressão diminui no sentido de escoamento,
e o gradiente de pressão é dito ser favorável.
Considere o escoamento através de canal de seção transversal variável, mostrado na Fig.
3.6.

Prof. Dr. Admilson T. Franco 73

UTFPR - Curso de Engenharia Mecânica
Cap. 3 - Mecânica dos Fluidos 2
 Escoamento
reverso

Fig. 3.6 Escoamento de camada limite com gradiente de pressão (espessura de camada limite
exagerada para maior clareza)

O gradiente de pressão adverso, dp / dx  0 , é condição necessária para a separação,

ou seja, a separação não pode ocorrer a menos que dp / dx  0 .

O perfil turbulento é muito mais obtuso do que o laminar. Para a mesma velocidade de
corrente livre, o fluxo de quantidade de movimento dentro da camada limite turbulenta é maior do
que dentro da camada limite laminar (Fig.3.7b). A separação ocorre quando a quantidade de
movimento de camadas de fluido adjacentes perto da superfície é reduzida a zero pela ação
combinada de forças viscosas e de pressão. Conforme mostrado na Fig. 3.7b, a quantidade de
movimento do fluido próximo da superfície é significativamente maior para o perfil turbulento.
Assim, a camada turbulenta possui maior capacidade de resistir a separação com a condição de
gradiente de pressão adverso.
Os gradientes de pressão adversos causam importantes mudanças nos perfis de
velocidade para ambos os escoamentos de camada limite, turbulento e laminar. Soluções
aproximadas para escoamento com gradiente de pressão diferente de zero podem ser obtidas da
equação integral da quantidade de movimento
w d
= (U 2 ) +  *U
dU
(3.17)
 dx dx
Expandindo o primeiro termo, pode-se escrever
w d dU
=U2 + ( * + 2 ) U
 dx dx
ou

w C f d  dU
= = + ( H + 2) (3.28)
U 2
2 dx U dx

Prof. Dr. Admilson T. Franco 74

UTFPR - Curso de Engenharia Mecânica
Cap. 3 - Mecânica dos Fluidos 2
 *
sendo H = o “fator de forma” do perfil de velocidade. O fator de forma aumenta para gradiente

de pressão adverso. Para escoamento de camada limite turbulento, H aumenta de 1, 3 para

gradiente de pressão zero para aproximadamente 2,5 na separação. Para escoamento laminar

com gradiente de pressão zero, H = 2, 6 ; na separação H = 3,5 .

A distribuição de velocidade de corrente livre, U ( x ) , deve ser conhecida antes que a Eq.

(3.28) possa ser aplicada. Como dp / dx = − UdU / dx , especificar U ( x ) é equivalente a

especificar o gradiente de pressão. Pode-se obter a primeira aproximação de U ( x ) da teoria do

escoamento ideal para o escoamento não viscoso nas mesmas condições.

Fig. 3.7 Perfis adimensionais para escoamento de camada limite sobre placa plana.

Prof. Dr. Admilson T. Franco 75

UTFPR - Curso de Engenharia Mecânica
Cap. 3 - Mecânica dos Fluidos 2