APOSTILA_ESCOAMENTOS_VISCOSOS_2010

ESCOAMENTOS VISCOSOS

PUCRS - DEM – Prof. Alé (2010) 1

APOSTILA DE ESCOAMENTOS VISCOSOS

Prof. Jorge Villar Alé 2º Semestre

2010


Escoamentos Viscosos 2



1.1 ESCOAMENTO VISCOSO E INCOMPRESSÍVEL.............................................................................................5

1.1.1 Conceito de Escoamento Plenamente Desenvolvido...................................................................................5

1.2 DISTRIBUIÇÃO DA TENSÃO DE CISALHAMENTO EM TUBOS.... ......................................................7

1.3 ESCOAMENTO LAMINAR EM TUBULAÇÕES ................... ......................................................................9

1.4 ESCOAMENTO TURBULENTO EM TUBULAÇÕES...............................................................................................12 1.4.1 Tensão de cisalhamento............................................................................................................................13 1.4.2 Distribuição da Velocidade no Escoamento Turbulento..........................................................................14

1.5 EQUAÇÃO DE ENERGIA COM VELOCIDADE MÉDIA .........................................................................................16 1.6 PERDA DE PRESSÃO NO ESCOAMENTO EM TUBULAÇÕES................................................................................17 1.7 PERDA DE CARGA PRINCIPAL .........................................................................................................................18

1.7.1 Perda de Carga Principal - Escoamento Laminar ...................................................................................18 1.7.2 Perda de Carga Principal - Escoamento Turbulento ...............................................................................19 1.7.3 DIAGRAMA DE MOODY.........................................................................................................................20

1.8 MÉTODOS PARA DETERMINAR AS PERDAS DE CARGA SECUNDÁRIAS.................................................................23 1.8.1 Método do comprimento equivalente........................................................................................................23 1.8.2 Método do coeficiente de perda de carga .................................................................................................24

1.9 PERDA DE CARGA EM ELEMENTOS SECUNDÁRIOS..........................................................................................25 1.9.1 Saídas e Entradas Abruptas.....................................................................................................................25 1.9.2 Expansão e Contração Abruptas ..............................................................................................................26 1.9.3 Expansão e Contração Gradual ...............................................................................................................27

1.10 PROBLEMAS TÍPICOS DE ESCOAMENTOS EM TUBOS.......................................................................................28 1.10.1 Determinação da Vazão.......................................................................................................................28 1.10.2 Determinação do Diâmetro da Tubulação...........................................................................................28

1.11 RESUMO DA TENSÃO DE CISALHAMENTO NAS PAREDES................................................................................29 1.12 CONCEITO DE DIÂMETRO HIDRÁULICO .........................................................................................................30 2.1 ESCOAMENTOS TURBULENTOS - TENSOES DE REYNOLDS...........................................................32 2.2 REPRESENTACOES SEMI-EMPIRICAS DAS TENSOES DE REYNOLDS............................................36 2.3 CONCEITO DE MISTURA DE PRANDTL.................................................................................................36 3.1 PERFIL DE VELOCIDADES EM ESCOAMENTOS TURBULENTOS – TUBOS....................................40 3.2 PERFIL DE VELOCIDADES EM ESCOAMENTOS TURBULENTOS – TUBOS LISOS.................................42 3.3 PERFIL DE VELOCIDADES EM ESCOAMENTOS TURBULENTOS – TUBOS RUGOSOS.........................47 3.4 PERFIL DE VELOCIDADES EM ESCOAMENTOS TURBULENTOS – LEI EXPONENCIAL......................48 4.1 ESCOAMENTO VISCOSO EXTERNO: CONCEITOS DE CAMADA LIMITE .............................................................51 4.2 ESCOAMENTO EM TORNO DE CORPOS............................................................................................................51

4.2.1 Efeito do Número de Reynolds no Escoamento Externo...........................................................................51 4.3 ESCOAMENTO SOBRE PLACA PLANA ..............................................................................................................52

4.3.1 Forças Viscosas Predominantes – Reynolds muito baixo - Re≈≈≈≈0,1 ..........................................................52 4.3.2 Forças Viscosas Moderadas – Reynolds baixo - Re≈≈≈≈10 ..........................................................................52 4.3.3 Forças de Viscosas Confinadas – Reynolds Alto - Re≈≈≈≈107 ......................................................................53

4.4 CARACTERÍSTICAS DA CAMADA LIMITE .........................................................................................................54 4.5 ESPESSURA DA CAMADA LIMITE ....................................................................................................................55 4.6 ESPESSURA DE DESLOCAMENTO.....................................................................................................................55 4.7 ESPESSURA DA QUANTIDADE DE MOVIMENTO ...............................................................................................56 4.8 COEFICIENTE DE ARRASTO EM PLACAS PLANAS ............................................................................................57 4.9 COEFICIENTE DE ARRASTO E FORÇA DE ARRASTO PELA TENSÃO DE CISALHAMENTO ...................................58 4.10 EQUAÇÕES DE BLASIUS – PLACA PLANA – CAMADA LIMITE LAMINAR .......................................................59 ESPESSURA DA CAMADA LIMITE ...................................................................................................................................59 ESPESSURA DE DESLOCAMENTO DA CAMADA LIMITE....................................................................................................59 4.11 COEFICIENTE DE ARRASTO LOCAL OU COEF. DE TENSÃO DE CISALHAMENTO.................................................59 4.12 COEFICIENTE DE ARRASTO MÉDIO..................................................................................................................59 4.13 TRANSIÇÃO DE ESCOAMENTO LAMINAR PARA TURBULENTO.........................................................................60 4.14 CAMADA LIMITE TURBULENTA EM PLACA PLANA .........................................................................................61

4.14.1 Coeficiente de Arrasto Local - Coeficiente de Tensão de Cisalhamento .............................................61 4.14.2 Coeficiente de Arrasto Médio ..............................................................................................................61

4.15 ESPESSURA DA CAMADA LIMITE - ESCOAMENTO TURBULENTO.....................................................................63 4.16 RESUMO DAS EQUAÇÕES DA CAMADA LIMITE EM PLACA PLANA ..................................................................64



5.1 ESCOAMENTOS EXTERNOS: CAMADA LIMITE EM PLACA PLANA ...............................................67 5.2 RESULTADOS PARA ESCOMANETO LAMINAR:..................................................................................70 5.3 RESUMO DAS EQUACOES DE CAMADA LIMITE EM PLACA PLANA - LAMINAR E

TURBULENTO...........................................................................................................................................................77 5.4 RELACOES BASICAS......................................................................................................................................78 5.5 COEFICIENTE DE ARRASTO EM PLACA PLANA – REGIME LAMINAR E TURBULENTO ......................................79 6.1 ESCOAMENTOS EXTERNOS CAMADA LIMITE EM PLACA PLANA .........................................................81 7.1 FORÇAS AERODINÂMICOS DE SUSTENTAÇÃO E ARRASTO.......................................................................................88

5.5.1 Coeficiente de Arrasto ..............................................................................................................................89 7.2 ESCOAMENTO SOBRE CILINDROS - EFEITO DA VISCOSIDADE...................................................................................91 7.3 ESCOAMENTO NÃO VISCOSO NUM CILINDRO...........................................................................................................92 7.4 ESCOAMENTO VISCOSO NUM CILINDRO : EFEITO DO GRADIENTE ADVERSO DE PRESSÃO......................................94 7.5 SUSTENTAÇÃO AERODINÂMICA .............................................................................................................................99 7.6 RELAÇÃO ENTRE COEFICIENTE DE PRESSÃO E SUSTENTAÇÃO..............................................................................101 7,7 CURVA DE SUSTENTAÇÃO VERSUS ÂNGULO DE ATAQUE......................................................................................102

5.5.2 Influência da Velocidade Induzida na Força de Arrasto ........................................................................105 5.5.3 Velocidade mínima de vôo .....................................................................................................................107



1.1 ESCOAMENTO VISCOSO E INCOMPRESSÍVEL

1.1.1 Conceito de Escoamento Plenamente Desenvolvido

Consideramos no estudo o escoamento viscoso interno num tubo com fluido incompressível (Fig.1). Se o tubo

estivesse imerso num reservatório (ou na saída de um reservatório) a velocidade U0 na entrada poderia ser considerada como uniforme. A medida que o fluido entra no tubo os efeitos viscosos provocam aderência do fluido às paredes do tubo. Esta é conhecida como condição de não deslizamento. Assim, o fluido em contato com as paredes sempre terá velocidade nula ao longo de todo o comprimento da tubulação.

Figura 1 Região de entrada em um tubo

A medida que o fluido escoa para dentro do tubo (na direção x) se desenvolve uma camada limite, devido ao

efeito das forças de cisalhamento das paredes, que retardam o escoamento. A medida que avança para o interior do tubo tal efeito aumenta. Os efeitos viscosos são importantes dentro da camada limite. Na região do núcleo não atingida pela camada limite os efeitos viscosos são desprezíveis.

Considerando que o escoamento é incompressível a velocidade na linha central do tubo aumenta com a distância a partir da entrada satisfazendo a equação da continuidade. O perfil de velocidades u(r,x) muda conforme aumenta a camada limite. Contudo como a seção do tubo é constante a velocidade média deve ser a mesma em qualquer seção:

∫= AduA

utotal

rr1

Como na região de entrada a velocidade é uniforme também é verdadeiro que u=U0: Numa determinada

posição x a camada limite atinge a linha central da tubulação e o perfil de velocidade não muda com a posição x que encontramos no tubo.



Comprimento de entrada L Distância da entrada até o local onde a camada limite atinge a linha central (de simetria) do tubo (x=L). A partir deste ponto o perfil de velocidade é plenamente desenvolvido significando que seu formato não varia mais na direção de x. • Para x > L o perfil de velocidade não varia mais com x, nesse caso denomina-se perfil de velocidades plenamente

desenvolvido. Posição no tubo Perfil de velocidades Na entrada do tubo x=0. u=Uo = constante Na região de desenvolvimento x ≤ L u=u(r,x) Na região plenamente desenvolvida x > L u=u (r)

• O formato do perfil plenamente desenvolvido depende se o regime de escoamento é laminar ou turbulento. Para Escoamento Laminar o comprimento de entrada é função do número de Reynolds:

Re06,006,0 =≈µ

ρ DV

D

L

Onde ρ é a massa especifica do fluido (kg/m3), V é a velocidade média do escoamento (m/s), D é o diâmetro interno da tubulação (m) e µ é a viscosidade dinâmica do fluido (Pa.s). Considerando que o escoamento é laminar até

2300Re< podemos estimar o comprimento de entrada neste caso.

DL

DxL

DL

140

230006,0

Re06,0

≈≈≈

O escoamento laminar plenamente desenvolvido ocorrerá para L > 100 D

Para escoamento turbulento a mistura entre camadas de fluido aumenta rapidamente a camada limite (mais rápido que a laminar). A experiência mostra que a velocidade torna-se plenamente desenvolvida para

DL )40...25(≈

Dependendo das características do escoamento turbulento podem ser encontrados casos em que o escoamento atinge um perfil de velocidades plenamente desenvolvido para valores de L ≅80D. Para estimar-se o comprimento L num escoamento turbulento pode ser dotada a expressão:

( ) 6/1Re4,4≈D

L



1.2 DISTRIBUIÇÃO DA TENSÃO DE CISALHAMENTO EM TUBOS

No escoamento permanente plenamente desenvolvido num tubo horizontal, seja laminar ou turbulento, a queda de pressão é equilibrada pelas forças de cisalhamento nas paredes do tubo.

Figura 2 Volume de controle para análise da tensão de cisalhamento

Aplicando a equação da quantidade de movimento na direção x

∫∫ +∂∂=+

scvc

bxsx AdVudVut

FFrr

ρρ

Hipóteses: (1) Tubo horizontal FBx=0 (2) Escoamento permanente. (3) Escoamento incompressível. (4) Escoamento plenamente desenvolvido. Desta forma 0=sxF . Para o elemento de fluido da Fig. 8.2 o balanço de forças é dado por:

02

0222

2

22

=+∂∂−

=+

∂∂+−

∂∂−=

rdxrdx

x

p

rdxrdx

x

ppr

dx

x

ppF

rx

rxsx

πτπ

πτππ

obtendo-se finalmente:

x

prrx ∂

∂=2

τ Válido para escoamento Laminar ou Turbulento

desta forma a tensão de cisalhamento no fluido varia linearmente na direção transversal ao tubo, de zero na linha de centro até um máximo na parede. Denominando tensão de cisalhamento na parede como τw, e sabendo que a variação da pressão ao longo do tubo é constante

x

pRRrrxw ∂

∂−=−== 2

ττ



A expressão fica negativa (-) já que se considera a tensão de cisalhamento na parede com a mesma magnitude da tensão do fluido porém agindo em sentido contrário.

como ( )

cteL

P

L

pp

x

p =∆−=−

=∂∂ 12

Figura 3 Perda de presão numa tubulação

substituída na equação anterior obtém-se a equação que relaciona a tensão de cisalhamento na parede com a queda de pressão em tubos válida para escoamento laminar ou turbulento.

L

pRw

∆=2

τ ou L

pDw

∆=4

τ

A distribuição da tensão de cisalhamento é mostrada na figura abaixo. É representada como uma função linear do tipo

crrx =τ onde a constante c=∆P/2L.

Figura 4 Perfil de velocidade e de tensão de cisalhamento em tubulações

Desta forma podemos relacionar a queda de pressão com a tensão de cisalhamento na parede

wD

Lp τ4=∆

Uma pequena tensão de cisalhamento na parede pode produzir uma grande diferença de pressão quando a tubulação for muito longa (L/D >> 1). Obs: As equações da tensão de cisalhamento obtidas aqui são válidas para escoamento laminar e turbulento já que a dedução foi realizada independente destes regimes de escoamento.



1.3 ESCOAMENTO LAMINAR EM TUBULAÇÕES Perfil de Velocidades No escoamento laminar unidimensional a tensão de cisalhamento é dada por:

dr

durx µτ =

Explicitando desta equação a velocidade:

drdu rxτµ1=

substituindo o termo da tensão de cisalhamento: rL

Prx

∆=2

1τ

rdrL

Pdu

∆=2

11

µ

integrando

∫∫

∆=R

r

rdrL

Pdu

2

11

µ

R

r

r

L

Pu

∆=

22

11 2

µ

{ }22

4

1rR

L

Pu −∆=

µ

ou também:

−∆=22

14 R

r

L

PRu

µ

Esta equação representa o perfil de velocidades para escoamento laminar em tubos.



Vazão Volumétrica A vazão volumétrica ou simplesmente vazão no elemento de fluido da Fig. 8.2 é dada por:

rdrudQ π2=

∫∫ =R

rdrudQ0

2π

substituindo a velocidade u=u(r) pelo termo deduzido anteriormente: { }22

4

1rR

L

Pu −∆=

µ

{ } ∫−∆=R

rdrrRL

PQ

0

22 24

πµ

{ }∫ −∆=R

rdrrRL

PQ

0

2224

πµ

−∆=

422

4

44 RR

L

PQ π

µ

42

4

4R

L

PQ π

µ∆=

4

8R

L

PQ

µπ∆=

ou em função do diâmetro:

L

PDQ

µπ128

4∆= (Equação de Hagen - Poussiulle)

Velocidade Média

2

4

D

Q

A

QV

π==

Substituindo a expressão de Hagen-Poussiulle:

L

PD

DV

µπ

π 128

4 4

2

∆=

L

PDV

µ32

2∆= Ou também L

PRV

µ8

2∆=



Velocidade Máxima Sabemos que o perfil de velocidades num escoamento laminar é dada por:

−∆=2

2 14 R

rR

L

Pu

µ

A velocidade máxima ocorre na linha central do tubo, isto é para r=0.

L

PRu

µ4

2

max

∆=

Relação entre Velociade Máxima e Velocidade Média:

2

4

82

2

max =∆

∆

=

L

PR

L

PR

V

u

µ

µ

Vu 2max = (para escoamento Laminar)

Perfil de Velocidades em Função da Velocidade Máxima

−∆=2

2 14 R

rR

L

Pu

µ

−=2

max 1R

ruu

O perfil de velocidades num escoamento laminar é parabólico

Figura 5 Perfil de velocidade para escoamento laminar numa tubulação



1.4 Escoamento Turbulento em Tubulações A natureza do escoamento nos tubos pode ser laminar ou turbulento. Tais regimes são dependentes do valor do número de Reynolds.

µρ DV=Re

Onde ρ é a massa específica do fluido (kg/m3), V é a velocidade média do escoamento (m/s), D é o diâmetro interno da tubulação (m) e µ é a viscosidade dinâmica do fluido (Pa.s). Fluido Laminar: O fluido escoa em camadas (lâminas) não existe mistura macroscópica das camadas adjacentes. Escoamento Turbulento: Manifestam-se pequenas flutuações da velocidade de alta freqüência superpostas ao movimento predominante. Medindo a componente da velocidade x num local fixo da tubulação podemos observar na Fig.8.6 o comportamento da velocidade para o caso laminar e turbulento. No escoamento turbulento a velocidade instantânea (u) é tão uniforme que sua velocidade é a mesma

uu = Observa-se que no caso do escoamento turbulento existe uma componente aleatória de flutuação da velocidade instantânea (u´). Desta forma a velocidade instantânea é dada pela soma algébrica velocidade média mais a componente de flutuação:

uuu +=

Figura 6 Variação da velocidade num escoamento laminar e turbulento unidimensional

No caso do escoamento real tridimensional a natureza do escoamento é mais complicada já que a velocidade manifesta três componentes de flutuação, sendo a velocidade instantânea dada como:

´

´

´

www

vvv

uuu

+=+=+=



1.4.1 Tensão de cisalhamento

No escoamento laminar unidimensional a tensão de cisalhamento é dada por:

dy

duyx µτ =

Conhecido o perfil de velocidades, podemos através da sua derivada (du/dy), determinar as tensões de cisalhamento no escoamento. Para escoamento turbulento não se tem uma relação direta como no caso do escoamento laminar, mesmo com velocidade média unidimensional. As flutuações aleatórias da velocidade tridimensional u´, v´, w´ transportam quantidade de movimento aumentando a tensão de cisalhamento efetiva. Desta forma não existe uma relação universal entre o campo de tensões e da velocidade no caso do escoamento turbulento. No caso do escoamento turbulento para determinar as tensões de cisalhamento utilizam-se teorias semi-empíricas e de dados experimentais. Neste caso a tensão de cisalhamento se expressa como sendo formada por uma componente laminar e outra turbulenta.

turbulentoarla τττ += min

onde dy

udlam µτ = ´´vuturb ρτ −=



1.4.2 Distribuição da Velocidade no Escoamento Turbulento

(a) Lei Exponencial Empírica Num escoamento turbulento o perfil de velocidades não pode ser deduzido da maneira como foi realizado para o escoamento laminar, devido a que não podemos utilizar a lei de Newton para relacionar a tensão de cisalhamento com o gradiente de velocidades.

Figura 7 Perfil de Velocidades num escoamento turbulento

Num escoamento turbulento adotam-se perfis de velocidades obtidos de relações empíricas. Por exemplo a lei exponencial empírica considera um perfil do tipo:

n

R

ruru

/1

max 1)(

−=

Tal equação não pode ser aplicada próxima à parede (R=0) já que o gradiente de velocidade é infinito. Contudo pode ser utiliza para y/R < 0,004 sendo y= R – r. O termo n depende do número de Reynolds como mostra a Fig. 8.8. O valor para n=7 é geralmente utilizado com precisão razoável em muitas situações reais. Também podemos utilizar a expressão 96.1log(Re)85.1 −=n

Figura 8 Expoente n do perfil da lei exponencial de velocidade turbulento



A Fig.8.9 mostra um perfil turbulento utilizando a expressão exponencial com n=6 e n=10. Para comparação também mostra-se o perfil laminar de velocidade. Observa-se que os perfis turbulentos são muito mais “achatados” que os laminares. O achatamento aumenta com o número de Reynolds isto é´, com o aumento de n.

Figura 9 Perfil de velocidades num tubo

A razão entre a velocidade média ( Vouu ) e a velocidade máxima (Umax) para um perfil exponencial de velocidade é dada por:

( )( )121

2 2

max ++=

nn

n

U

u

(b) Distribuição da Velocidade Considerando Fator de Atrito O fator de atrito ( f ) pode ser determinado para escoamentos em regime laminar e turbulento. O expoente n pode ser determinado no caso de escoamento turbulento como:

fkn

8=

onde k=0,41 é denominada constante de von karman. No caso de escoamento turbulento podemos também utilizar a seguinte expressão para determinar o perfil de velocidades em função do fator de atrito ( f )

( ){ }RyffVu /log15,243,11 10++=

onde y=r - R.

A velocidade máxima é dada como:

{ }fVu 43,11max +=



1.5 Equação de Energia com Velocidade Média

Considerando escoamento em regime permanente incompressível uma análise de energia entre duas seções, que incluem dissipação e/ou ganhos adicionais de energia, pode ser representada como:

2

222

1

211

22z

g

u

g

phHHz

g

u

g

pLTRA ++=−−+++

ρρ

onde HA representa a energia adicionada, HR, representa a energia retirada do sistema e hLT representa a dissipação de energia. Num problema em particular nem todos os termos de energia são utilizados. Nos escoamentos viscosos o perfil de velocidade numa dada seção não pode ser uniforme. É conveniente portanto utilizar a velocidade média, para tal é necessário definir o coeficiente de fluxo de energia cinética (α). Aplicando a equação de energia numa tubulação entre os pontos 1 e 2, onde não existem dispositivos mecânicos (HA =0 e HR =0):

2

22

22

1

21

11

22z

g

u

g

phz

g

u

g

pLT ++=−++ α

ρα

ρ

o coeficiente de energia cinética é definido como

2

3

Vm

dAVA

&

∫=

ρα

� No caso de escoamento laminar: α=2.

� No caso de escoamento turbulento:

( )( )

++

=nn

n

U

U

233

2 23

maxα

Por ex. para os números de Reynolds considerados

Re n αααα

4,0x103 6 1,08

3,2x106 10 1,03

Se observa que α≅1. Desta forma para a maioria dos casos de engenharia nos cálculos de perda de carga considera-se α=1.

Observação: No material do Fox, McDonald a energia e perda de carga é apresentada como energia por unidade de massa (J/kg). No nosso caso é dada como energia por unidade de peso (J/N), ou metro de coluna de fluido (m).



1.6 Perda de Pressão no Escoamento em Tubulações A variação de pressão num duto resulta da variação da elevação, da velocidade e do atrito. Escoamento sem atrito A variação de pressão pode ser determinada aplicando a Eq. de Bernoulli.

),( VZfP →∆ já que hLT=0. Escoamento real com atrito a variação de pressão pode ser determinada aplicando a Eq. da Energia

),,( LThVZfP →∆ � O atrito origina uma diminuição da pressão. � Causa uma perda de pressão comparada com o caso de escoamento sem atrito.

Figura 10 Perda de carga em sistema de bombeamento

Perda de carga Total A perda de carga em tubulações é dada por duas parcelas.

accLLT hhh +=

Perda de Pressão ou de Carga Principal: (hL) � Devido ao atrito no escoamento plenamente desenvolvido entre pontos da tubulação com área constante. Perda de Carga Secundária - (hac) � Devido ao escoamento através de acessórios como válvulas, joelhos, registros e em porções do sistema de área

variável tais como saídas de reservatórios, bocais convergentes e divergentes. � A perda de carga na entrada ou saída de uma tubulação é considerada como perda de carga secundária.



1.7 Perda de Carga Principal � Transformação da energia cinética para energia térmica por efeitos viscosos. Consideremos um escoamento plenamente desenvolvido numa tubulação de comprimento L. Analisando uma tubulação com área constante A1=A2 e desta forma pela Eq. da continuidade u1=u2

( ) ( )2121 zz

g

pphL −+

−=

ρ

No caso de uma tubulação horizontal (z1=z2).

( )g

P

g

pphL ρρ

∆=−

= 21

1.7.1 Perda de Carga Principal - Escoamento Laminar

Utilizando a expressão da velocidade média

∆=L

PDV

µ32

2

e desta forma: 2

32

D

LVP

µ=∆

substituindo esta última expressão na equação da perda de carga

=∆=gD

LV

g

PhL ρ

µρ

1322

Podemos expressar esta equação em função do Número de Reynolds

µρ DV=Re explicitando a viscosidade dinâmica:

Re

DVρµ =

=

=gD

LV

gD

LVDVhL

132

Re

132

Re

2

2 ρρ

expressando em função da energia cinética

g

V

D

LhL 2Re

64 2

= Escoamento laminar



1.7.2 Perda de Carga Principal - Escoamento Turbulento

� No caso de escoamento turbulento não existem expressões que permitam avaliar analiticamente a queda de

pressão. � Se utiliza análise dimensional e correlações de dados experimentais. Analisando o caso de escoamento turbulento plenamente desenvolvido a queda de pressão é função das seguintes variáveis.

),,,,,( µρεφ VLDP =∆ D diâmetro da tubulação L, comprimento da tubulação, V, Velocidade média, ε, rugosidade absoluta, ρ massa específica, µ, viscosidade dinâmica. Aplicando-se análise dimensional se obtém uma expressão da forma:

=∆

DD

L

VDV

P ερ

µφρ

,,2

como o termo é dado por g

PhL ρ

∆= podemos explicitar a variação de pressão (∆P) e substituir a mesma na equação

do análise dimensional.

==DD

L

V

gh

V

gh LL εφρρ

,Re,22

experimentos mostram que a perda de carga é diretamente proporcional a L/D. Para que a perda de carga seja obtida adimensionalizada em relação à energia cinética se introduz o termo 1/2 na equação ficando como:

=DD

L

g

V

hL εφ Re,

2

2

A função φ é conhecida como fator de atrito ou coeficiente de atrito.

=D

fεφ Re,

desta forma se obtém a equação da perda de carga que representa a energia dissipada por unidade de peso do fluido escoando.

g

V

D

LfhL 2

2

= Equação de Darcy-Weisbach.

O fator de atrito determina-se experimentalmente. Utiliza-se o Diagrama de Moody.



1.7.3 DIAGRAMA DE MOODY

Para determinar o fator de atrito se utiliza o Diagrama de Moody. Para tal deve-se ter o valor do número de

Reynolds e a rugosidade relativa ε/D. A rugosidade absoluta ε depende do tipo de material da tubulação e do seu acabamento. Representa o valor médio das alturas da rugosidade da parede interna da tubulação. A Tabela dada mostra os valores da rugosidade absoluta para os materiais típicos de tubulações industriais utilizadas para o escoamento de fluidos.

Figura 11 Representação da rugosidade absoluta em tubulações

Tabela 1 Rugosidade absoluta (mm) de tubulações industriais Material Rugosidade absoluta (mm)

Aço, revestimento asfalto quente 0,3 a 0,9 Aço, revestimento esmalte centrifugado 0,011 a 0,06 Aço enferrujado ligeiramente 0,15 a 0,3 Aço enferrujado 0,4 a 0,6 Aço muito enferrujado 0,9 a 2,4 Ferro galvanizado novo, com costura 0,15 a 0,2 Ferro galvanizado novo, sem costura 0,06 a 0,15 Ferro fundido revestido com asfalto 0,12 a 0,20 Ferro fundido com crostas 1,5 a 3,0 PVC e Cobre 0,015 Cimento-amianto novo 0,05 a 0,10

Fonte: - Equipamentos Industriais e de Processo - (Macintyre)



O diagrama de Moody apresenta uma zona laminar (Re < 2000), uma zona crítica (Re de 2000 e 4000) uma zona

de transição e uma zona inteiramente rugosa. Nestas zonas o fator de atrito f apresenta diferentes dependências em relação ao número de Reynolds (Re) e em relação a rugosidade relativa ε/D as quais são resumidas a seguir: 1. Na zona laminar fator de atrito f é independente da rugosidade ε/D e inversamente proporcional ao número de Re 2. Na zona crítica o fator de atrito apresenta aumentos bruscos. 3. Na zona de transição para um determinado Re o fator de atrito f diminui conforme a rugosidade relativa ε/D

diminui. 4. Na zona de transição, para uma determinada rugosidade relativa ε/D o fator de atrito f diminui ao aumentar o Re

até alcançar a região inteiramente rugosa. 5. Dentro da zona inteiramente rugosa, para uma determinada rugosidade relativa ε/D, o fator de atrito f, se mantém

praticamente como um valor constante independente do Re. 6. Na zona de transição, conforme diminui a rugosidade relativa ε/D o valor do Re no qual inicia a região plenamente

turbulenta começa a aumentar

Figura 12 Representação do Diagrama de Moody



I - Escoamento Laminar O fator de atrito para escoamento laminar pode ser obtido igualando a equação

g

V

D

LfhL 2

2

= com a equação da perda de carga laminar g

V

D

LhL 2Re

64 2

= se obtém:

Re

64=f válido para Re < 2500

� No escoamento laminar o fator de atrito ( f ) é função somente do número de Reynolds. � Independe da rugosidade da tubulação. II - Escoamento com Tubos Hidraulicamente Lisos Nesta região pode utilizar-se a Eq. de Blasius ou a Eq. de Drew Koo e McAdams

( ) 4/1Re

316,0=f Eq. de Blasius 4000 < Re < 105

32,0Re5,00056,0 −+=f Eq. de Drew Koo e McAdams 105 < Re < 3x106

III - Escoamento Turbulento com Tubos Hidraulicamente Semi-Rugosos Permite determinar o fator de atrito para escoamento turbulento:

+−=

f

D

f Re

51,2

7,3

/log0,2

1 ε Equação de Colebrook 5,0x103 < Re < 1x108

Como tal equação é do tipo transcendente deve ser utilizado um procedimento iterativo para determinar f. Uma alternativa é utilizar uma equação explícita:

2

9,0Re

74,5

7,3

/log25,0

−

+= Df

ε Equação Explícita 5,0x103 < Re < 1x108

Utilizando a Eq. acima se encontram valores de f com margem de erro de +-1% comparados com os obtidos com a Eq. de Colebrook, para: ε/D de 1,0x10-4 (0,0001) até 1x10-6 (0,000001) IV - Escoamento Turbulento com Tubos Hidraulicamente Rugosos O fator de atrito depende unicamente da rugosidade relativa e pode ser determinado pela equação:

−=7,3

/log2

1 D

f

ε Equação de Von Karman



1.8 Métodos para Determinar as Perdas de Carga Secundárias

1.8.1 Método do comprimento equivalente

Os acessórios são todos aqueles elementos que existem numa tubulação através dos quais o fluido escoa, tais como curvas, bocais, registros e válvulas. Cada um destes elementos produz uma dissipação de energia que é avaliada pela perda de carga (hac) definida como:

g

V

D

Lfh eq

ac 2

2

= (m)

O comprimento equivalente em metros de canalização retilínea (Leq) é tabelado segundo o tipo de acessório, o material utilizado e o diâmetro da tubulação. Se substituirmos um certo acessório por uma tubulação retilínea com o comprimento igual ao comprimento equivalente (com igual material e diâmetro) ambos originariam a mesma perda de carga. A tabela abaixo mostra o comprimento equivalente adimensional (Leq/D) de diversos acessórios.

Figura 13 Representação do comprimento equivalente em acessórios

Tabela 2 Perda de carga localizada Tipo de Acessório Comprimento Equivalente

(Leq/D) Válvula de globo aberta 340 Válvula de gaveta aberta 8 3/4 aberta 35 1/2 aberta 160 1/4 aberta 900 Válvula tipo borboleta aberta 45 Válvula de esfera aberta 3 Válvula de retenção tipo globo 600 Válvula de retenção tipo em ângulo 55 Válvula de pé com crivo: de disco móvel 75 Cotovelo padronizado 900 30 Cotovelo padronizado 450 16 Te padronizada fluxo direto 20 Te padronizada fluxo ramal 60

Válvula globo

Válvulas tipo borboleta

Te com flanges

Figura 14 acessórios utilizados em instalações industriais



1.8.2 Método do coeficiente de perda de carga

Uma outra forma de representar a perda de carga nos acessórios (hac) é definindo a mesma na forma:

g

VKhac 2

2

= (m)

Onde K é o coeficiente de perda de carga e V a velocidade média. O coeficiente de perda de carga será maior quanto mais abruto seja o elemento originando zonas de recirculação de fluxo e altos níveis de turbulência, aumentando desta forma a energia dissipada. A tabela mostra o coeficiente de perda e carga de diversos elementos.

Tabela 3 Coeficiente de perda de carga de acessórios

Tipo de Acessório K Tipo de Acessório K

Ampliação Gradual 0,20* Junção 0,40 Bocais 2,75 Medidor venturi 2,5 Comporta aberta 1,00 Redução gradual 0,15 Controlador de vazão 2,50 Registro de ângulo aberto 5,0 Cotovelo 900 0,9 Registro de gaveta aberto 0,20 Cotovelo 450 0,4 Registro de globo aberto 10,0 Crivo 0,75 Saída de canalização 1,00 Curva 90 0,4 Tê passagem direta 0,6 Curva 45 0,20 Tê saída de lado 1,30 Curva 22,5 0,10 Tê saída bilateral 1,80 Entrada normal em canalização 0,50 Válvula de pé 1,75 Entrada de borda 1,0 Válvula de retenção 2,50 Existência de pequena derivação 0,03

Velocidade 1,0 * com base na velocidade maior (seção menor) ** Relativa à velocidade de canalização

Igualando as equações de perda de carga por acessórios se obtém:

D

LfK eq=

mostrando a relação entre o coeficiente de perda de carga (K) e o comprimento equivalente (Leq).

Curva de 900

Joelho de 900 Registro de gaveta

Válvula de pé com crivo Figura 15 Exemplo de diversos acessórios utilizados em instalações industriais



1.9 Perda de Carga em Elementos Secundários

1.9.1 Saídas e Entradas Abruptas

Quando o fluido escoa de um tubo para um reservatório sua velocidade cai bruscamente até próximo de zero. A perda de carga para este caso é igual à energia cinética dissipada. K=1.

(a) saída de tubos K=1 (b) entrada de tubos: K depende do tipo de entrada

Figura 16 Representação de escoamento na saída e na entrada de tubos Entrada Abruta de um Reservatório para um Tubo No escoamento dado entre um reservatório e uma tubulação, a velocidade passa de um valor muito baixo para um valor elevado. O coeficiente de perda de carga depende do tipo de união entre o tubo e o reservatório. Três casos típicos apresentam diferentes perda de carga:

( a ) Entrada com tubo para dentro K=1,0 ( b ) Entrada com cantos vivos K=0,5 ( c ) Entrada com cantos arredondados K conforme os dados da tabela abaixo:

r/D 0,02 0,06 ≥0,15 K 0,28 0,15 0,04

(a) tubo para dentro K=1

(b) cantos vivos K=0,5

(c) cantos arredondados

Figura 17 Entrada com (a) tubo para dentro (b) cantos vivos e (c) cantos arredondados



1.9.2 Expansão e Contração Abruptas

Expansão abrupta Numa expansão abrupta o fluido escoa de um tubo de seção menor para um outro de seção maior. A velocidade cai abruptamente formando-se uma região de turbulência e recirculação de fluxo a qual provoca uma perda de carga proporcional à relação das seções dos tubos. A perda de carga localizada é determinada pela expressão: Onde V é a velocidade média do tubo menor.

(a) Contração abrupta (b) Expansão abrupta

Figura 18 Contração abrupta e expansão abrupta Contração Abrupta Neste tipo de elemento, a perda de carga é originada pela contração da linhas de corrente formando uma veia contracta e regiões de recirculação de fluxo.

0

0.1

0.2

0.3

0.4

0.5

0.6

0 0.2 0.4 0.6 0.8 1

A2/A1

K

(a) Contração abrupta

0

0.2

0.4

0.6

0.8

1

0 0.2 0.4 0.6 0.8 1

A1/A2

K

(b) Expansão abrupta

Figura 19 Coeficiente de perda de carga para contração e expansão abrupta

Para determinar a perda de carga com estas relaçoes se utiliza a velocidade correspondente a seção de menor diâmetro. O mesmos é valido para avaliar a perda de carga em peças com expansão o contração gradual como visto no proximo item.



1.9.3 Expansão e Contração Gradual

A expansão gradual é obtida com uma peça de transição unindo um tubo de menor diâmetro com outro de maior diâmetro permite uma menor dissipação de energia do que uma transição abrupta direta entre dois tubos de diferente diâmetro. O coeficiente de perda de carga (K) depende da relação de diâmetros (D2/D1) e do ângulo do cone. Obtém-se uma perda de carga mínima adotando-se um ângulo do cone de 70 .

(a) Contração gradual

(b) Expansão gradual

Figura 20 Contração gradual e expansão gradual

Figura 21 Perda de carga em expansão gradual

Contração Gradual Da mesma forma que numa contração brusca a perda de carga depende da relação de diâmetros e do ângulo da contração.

Tabela 4 Coeficiente de perda de carga (K) de contração gradual de tubos

Angulo da contração - θ A2/A1 10o 15 o a 40 o 50 o a 60 o 90 o 120 o 150 o 180 o

0,10 0,05 0,05 0,08 0,19 0,29 0,37 0,43 0,25 0,05 0,04 0,07 0,17 0,27 0,35 0,41 0,50 0,05 0,05 0,06 0,12 0,18 0,24 0,26

Obs. Válido para tubos redondos e retangulares. Fonte: Fox



1.10 Problemas Típicos de Escoamentos em Tubos A variação de pressão entre dois pontos de uma tubulação depende basicamente das variáveis envolvidas na Eq. da Energia.

1.10.1 Determinação da Vazão

Q = φ (L,hL,D) 1. Escrever a Eq. da energia introduzindo as grandezas conhecidas 2. Expressar a perda de carga em função da velocidade e do fator de atrito hL =φ (V,f) 3. Explicite a velocidade em função do fator atrito V= φ(f) 4. Expresse o número de Reynolds em função da velocidade Re =φ (V) 5. Calcule a rugosidade relativa ε/D. 6. Selecione um valor inicial do fator de atrito f=fo tomando como referência o valor da rugosidade relativa ε/D e

admitindo um Re na faixa turbulenta. 7. Calcule a velocidade em função do fator de atrito assumido Vcal=φ(f) 8. Calcule o Re com a nova velocidade Re =φ (Vcal) 9. Com Recal e ε/D obtenha um novo valor do fator de atrito f= fcal. 10. Se fcal ≠ f Adote f= fcal e repita o procedimento a partir do passo 7 até convergir o valor da fator de atrito. A solução do problema é encontrada quando o fator de atrito converge, determinado a vazão com a velocidade final calculada.

1.10.2 Determinação do Diâmetro da Tubulação

D = φ (L,Q, hL) 1. Explicite da Eq. da energia a perda de carga. 2. Expresse a vazão em função da velocidade e do diâmetro na Eq, da perda de carga. 3. Explicitar o diâmetro da Eq. da perda de carga ficando uma expressão na forma: D=(C1f)0,2 4. Expresse o número de Reynolds como função do diâmetro Re= C2/D. 5. Adote um valor inicial do fator de atrito f=f0 (por exemplo f0 =0,02) 6. Calcule o diâmetro pela expressão obtida: D=(C1f)0,2 7. Calcule o número de Reynolds pela expressão: Re= C2/D.

1. Calcule a rugosidade relativa ε/D. 2. Com Re e ε/D determine um novo valor do fator de atrito fcal. 3. Se fcal ≠ f adote f= fcal e repita o procedimento a partir do passo 7 até convergir o valor da fator de atrito. A solução do problema é encontrada quando o fator de atrito converge, determinado o diâmetro com o fator de atrito final.



1.11 Resumo da Tensão de Cisalhamento nas Paredes A tensão na parede no escoamento laminar e turbulento é dada por:

4

D

L

Pw

∆=τ

tal valor representa a tensão de cisalhamento máxima τw =τmax A Eq. de Darcy-Weisbach também é válida para escoamento laminar e turbulento

g

V

D

LfhL 2

2

=

A tensão de cisalhamento em função do fator de atrito (f) para regime laminar ou turbulento é obtida igualando-se as duas expressões anteriores obtendo-se

24

2Vfw ρτ = válida para escoamento laminar ou turbulento

A tensão de cisalhamento para qualquer posição r do duto é dada como:

R

rmaxττ = válida para escoamento laminar ou turbulento

r=0 é no centro da tubulação e r=R na parede da tubulação Perfil de velocidades e tensão de cisalhamento para escoamento Laminar e Turbulento

Figura 22 Escoamento laminar e turbulento: perfil de velocidades e tensão de cisalhamento



1.12 Conceito de Diâmetro Hidráulico Os equacionamentos de perda de carga estudados neste capítulo também podem ser aplicados a tubulações com seções não circulares utilizando a definição de diâmetro hidráulico (Dh) :

P

ADh

4=

Onde A é a área da seção transversal do tubo P é o perímetro molhado, que é o comprimento da parede em contato com o fluido. A equação acima para um duto circular A=πD2/4 e P=πD e desta forma Dh=D.

Figura 23 Diversas geometrias de tubulações

A Fig. 23 mostra diversas geometrias de seções transversais de tubos que podem ser utilizados nas aplicações industriais. Devido as limitações de espaço nas instalações de ar condicionado se utilizam freqüentemente dutos retangulares. Em trocadores de calor podem ser utilizados tubos achatados, hexagonais, ovais e outros com cilíndricos concêntricos para escoamento anular. Em canais de regadios, rios, córregos, canais de represamento e calhas o fluido não preenche totalmente a seção transversal do duto, isso deve ser considerado para determinar corretamente o perímetro molhado.




2. ESCOAMENTOS TURBULENTOS - TENSOES DE REYNOLDS



2.1 ESCOAMENTOS TURBULENTOS - TENSOES DE REYNOLDS Objetivo:

• Deduzir as tensões aparentes ou Tensões de Reynolds para escoamento turbulento utilizando as equações da conservação da massa e Eq. de Navier Stokes considerando fluido incompressível

• Apresentar as equações de Navier-Stokes para escoamento turbulento.

Definimos as equações que da conservação da massa e Eq. de Navier Stokes

• Equação valida para escoamento laminar e turbulento • Fluido incompressível (massa especifica constante) e viscosidade constante • Sem iteração térmica. • Condições de não-deslizamento e condições conhecidas na entrada e saída

Equação da conservação da massa

0=∂∂+

∂∂+

∂∂

z

w

y

v

x

u

(1)

Equação de Navier Stokes

VgpDt

VD rrr

2∇++−∇= µρρ (2)

Utilizando as grandezas escalares: Equação de Navier Stokes Componentes escalares

wgz

p

z

ww

y

wv

x

wu

t

w

vgy

p

z

vw

y

vv

x

vu

t

v

ugx

p

z

uw

y

uv

x

uu

t

u

z

y

x

2

2

2

∇++∂∂−=

∂∂+

∂∂+

∂∂+

∂∂

∇++∂∂−=

∂∂+

∂∂+

∂∂+

∂∂

∇++∂∂−=

∂∂+

∂∂+

∂∂+

∂∂

µρρ

µρρ

µρρ

(3)

No escoamento turbulento velocidade instantânea e definida como sendo a soma da media temporal mais a componente de flutuação

'

'

'

'

ppp

www

vvv

uuu

+=

+=

+=

+=

(4)



Media Temporal A media temporal u da função ),,,( tzyxu e definida como:

∫=

T

oudt

Tu

1 ∫=

T

ovdt

Tv

1 ∫=

T

owdt

Tw

1 ∫=

T

opdt

Tp

1

(5)

T

Período de calculo da media, o qual deve ser maior que o período das flutuações. Para escoamentos turbulentos em gases e água T=5seg. e um período apropriado.

Media da Flutuação

( ) 01' =−= ∫

T

odtuu

Tu

(6)

Media do Quadrado da Flutuação

01 2'2' ≠= ∫

T

odtu

Tu

(7)

Media do Produto das Flutuações 0'' ≠vu 0'' ≠wu 0'' ≠wv (8)

Utilizando as velocidades instantâneas e introduzidas na conservação da massa:

( ) ( ) ( )0

'''

=∂+∂+

∂+∂+

∂+∂

z

ww

y

vv

x

uu

0'''

=

∂∂+

∂∂+

∂∂+

∂∂+

∂∂+

∂∂

z

w

y

v

x

u

z

w

y

v

x

u

Tomando valores médios se obtém: Equação da Conservação da Massa para escoamento Turbulento. 0=

∂∂+

∂∂+

∂∂

z

w

y

v

x

u

(9)

Equação de Navier Stokes para Escoamento Turbulento Substituindo as definições de velocidades instantâneas se obtém:

( ) ( ) ( )

( ) ( ) ( )

( ) ( ) ( )

∂∂+

∂∂+

∂∂−∇++

∂∂−=

∂∂+

∂∂+

∂∂+

∂∂

∂∂+

∂∂+

∂∂−∇++

∂∂−=

∂∂+

∂∂+

∂∂+

∂∂

∂∂+

∂∂+

∂∂−∇++

∂∂−=

∂∂+

∂∂+

∂∂+

∂∂

z

w

y

wv

x

wuwg

z

p

z

ww

y

wv

x

wu

t

w

z

wv

y

v

x

vuvg

y

p

z

vw

y

vv

x

vu

t

v

z

wu

y

vu

x

uug

x

p

z

uw

y

uv

x

uu

t

u

z

y

x

2'''''2

''2'''2

''''2'2

ρµρρ

ρµρρ

ρµρρ

(10)

Observa-se comparando a Eq.10 com a Eq,1 que no caso do escoamento turbulento surge forcas adicionais devido à turbulência denominada forcas aparentes. As tensões associadas a estas forcas são chamadas de tensões aparentes ou tensões de Reynolds



Forcas Aparentes Associadas ao Escoamento Turbulento.

( ) ( ) ( )

( ) ( ) ( )

( ) ( ) ( )

∂∂+

∂∂+

∂∂−=

∂∂+

∂∂+

∂∂−=

∂∂+

∂∂+

∂∂−=

z

w

y

vw

x

uwdf

z

wv

y

v

x

uvdf

z

wu

y

vu

x

udf

Tz

Ty

Tx

2'''''

''2'''

''''2'

ρ

ρ

ρ

(11)

Estas forcas aparentes estão relacionadas as tensões do escoamento turbulento

∂∂+

∂∂

+∂

∂−=

∂∂

+∂

∂+

∂∂

−=

∂∂+

∂∂

+∂

∂−=

zyxdf

zyxdf

zyxdf

TzzTzyTzxTz

TyzTyyTyxTy

TxzTxyTxxTx

σττ

τστ

ττσ

(12)

Estas tensões são denominadas tensões de Reynolds as quais são conseqüência das flutuações da velocidade. ( ) ( ) ( )

( ) ( ) ( )( ) ( ) ( )

−=

2'''''

''2'''

''''2'

wvwuw

wvvuv

wuvuu

TzzTzyTzx

TyzTyyTyx

TxzTxyTxx

ρστττστττσ

(13)

Desta forma as tensões no escoamento turbulento podem ser consideradas como sendo a soma das parcelas das tensões laminares mais as tensões turbulentas ou tensões de Reynolds: Tensões de para Escoamento turbulento

Txxxxxx σσσ +=

Tyyyyyy σσσ +=

Tzzzzzz σσσ +=

Txyxyxy τττ +=

Tyzyzyz τττ +=

Tyzxzxzx τττ +=

(14)



Para escoamento com fluido incompressível a s tensões laminar e turbulenta são definidas pelas equações a seguir:

Txxxxxx σσσ +=

x

upxx ∂

∂+−= µσ 2 ( )2'uTxx ρσ −=

Tyyyyyy σσσ +=

y

vpyy ∂

∂+−= µσ 2 ( )2'vTyy ρσ −=

Tzzzzzz σσσ +=

z

wpzz ∂

∂+−= µσ 2 ( )2'wTzz ρσ −=

Txyxyxy τττ +=

∂∂+

∂∂=

x

v

y

uxy µτ ( )''vuTxy ρτ −=

Tyzyzyz τττ +=

∂∂+

∂∂=

z

v

y

wyz µτ ( )''wvTyz ρτ −=

Tyzxzxzx τττ +=

∂∂+

∂∂=

x

w

z

uzx µτ ( )''uwTzx ρτ −=

A Eq. 10 também pode ser representada como:

( ) ( ) ( )( ) ( ) ( )

( ) ( ) ( )

−∂∂

∂∂+

−

∂∂

∂∂+

−∂∂

∂∂++

∂∂−=

∂∂+

∂∂+

∂∂+

∂∂

−∂∂

∂∂+

−

∂∂

∂∂+

−∂∂

∂∂++

∂∂−=

∂∂+

∂∂+

∂∂+

∂∂

−∂∂

∂∂+

−

∂∂

∂∂+

−∂∂

∂∂++

∂∂−=

∂∂+

∂∂+

∂∂+

∂∂

2'''''

''2'''

''''2'

wz

w

zuw

y

w

yuw

x

w

xg

z

p

z

ww

y

wv

x

wu

t

w

wvz

v

zv

y

v

yuv

x

v

xg

y

p

z

vw

y

vv

x

vu

t

v

wuz

u

zvu

y

u

yu

x

u

xg

x

p

z

uw

y

uv

x

uu

t

u

z

y

x

ρµρµρµρρ

ρµρµρµρρ

ρµρµρµρρ

A equação acima pode ser escrita em forma mais compacta como: Equação de Navier-Stokes para Escoamento turbulento

zyxg

z

ww

y

wv

x

wu

t

w

zyxg

z

vw

y

vv

x

vu

t

v

zyxg

z

uw

y

uv

x

uu

t

u

zzzyzxz

yzyyyx

y

xzxyxxx

∂∂

+∂

∂+

∂∂

+=

∂∂+

∂∂+

∂∂+

∂∂

∂∂

+∂

∂+

∂∂

+=

∂∂+

∂∂+

∂∂+

∂∂

∂∂

+∂

∂+

∂∂

+=

∂∂+

∂∂+

∂∂+

∂∂

σττρρ

τστρρ

ττσρρ

(15)



2.2 REPRESENTACOES SEMI-EMPIRICAS DAS TENSOES DE REYNOLDS

• Ainda não existe um modelo de turbulência geral e completo que descreva como varia a tensão de cisalhamento num campo de escoamento incompressível viscoso e turbulento.

• Existe uma grande dificuldade em determinar as tensões de Reynolds ou que representa não conhecer a

viscosidade turbulenta efetiva. Duas soluções semi empíricas podem ser descritas:

(a) Conceito de Comprimento de Mistura de Prandtl em 1925 (b) Conceito de Viscosidade Turbulenta Efetiva de Boussinesq.

Considerando um dos termos para um escoamento numa direção predominante:

( )''vudy

udturblam ρµτττ −=+=

demonstra-se que o produto ( )''vu e sempre negativo (-):

( ) ( )'''' vudy

udvu

dy

ud ρµρµτ +=−−=

Desta forma num escoamento turbulento, a tensão total (lam + turb) e sempre maior que no escoamento laminar. Parcela lamτ : Dominante na sub-camada viscosa (muito fina: 0,1% do Raio da tubulação)

Parcela Turbτ : Dominante na camada externa ou camada turbulenta. (100 a 1000 maior que lamτ nesta camada)

• Transições dos efeitos laminar e turbulento ocorrem na camada de superposição ou amortecedora.

2.3 CONCEITO DE MISTURA DE PRANDTL

• As partículas de fluido viajam de camada para camada. • Neste transporte percorrem um caminho com comprimento de mistura l

Pode ser mostrado que as flutuações das velocidades são relacionadas por:

dy

udlu 1

' = e dy

udlv 2

' =

Onde l1 e l2 são os comprimentos de mistura para o transporte da quantidade de movimento. Definimos também que:



212 lllm =

O comprimento de mistura e definido como

kylm =

onde y e a distancia normal a parede e k e a constante de Von karman (k=0,4). Estudos posteriores mostram que ml

não apresenta um valor constante. Desta forma a tensão turbulenta e expressa como:

( )dy

ud

dy

udlvu mTurb

2'' ρρτ =−=

Tensão de Reynolds utilizando hipótese de Prandtl.

dy

ud

dy

udlmTurb

2ρτ =

Boussinesq define a viscosidade turbulenta efetiva, representando a tensão turbulenta como:

dy

udTurb ητ =

Tensão de Reynolds utilizando hipótese de Boussinesq

dy

udTurb ητ =

Onde η representa a viscosidade turbulenta efetiva que podemos relacionar a expressão de Prandtl por: Viscosidade turbulenta efetiva.

dy

udlm

2ρη =

Desta forma se obtém uma representação da tensão turbulenta como:

Tensão de Cisalhamento para escoamento turbulento

( )dy

udηµτ +=

dy

ud

dy

udturblam ηµτττ +=+=



A viscosidade efetiva turbulenta (η ) e relacionada nestas equações com a difusividade turbilhonar conhecida como denominada viscosidade cinemática aparente. Viscosidade cinemática aparente.

ρηε =m

Tensão de Reynolds utilizando a viscosidade cinemática aparente.

dy

udmTurb ρετ =

Desta forma a tensão de cisalhamento total num escoamento turbulento pode ser dada por: Tensão de Cisalhamento para escoamento turbulento

( )dy

udmρεµτ +=




3. PERFIL DE VELOCIDADES TUBULAÇOES

EM ESCOAMENTOS TURBULENTOS



3.1 PERFIL DE VELOCIDADES EM ESCOAMENTOS TURBULENTOS – TUBOS Num escoamento turbulento em dutos o perfil de velocidade cresce desde a parede até um máximo no centro da tubulação. Este escoamento pode ser divido em três regiões principais:

• Uma subcamada laminar ou viscosa muito próxima da parede • Uma camada intermediaria ou de superposição • Uma camada turbulenta externa (na região central da tubulação).

A natureza do escoamento e portando do perfil de velocidade e totalmente diferente nestas três camadas:

• Na subcamada a viscosidade do fluido e um parâmetro significativo e a massa especifica não. • Na camada externa a massa especifica e um parâmetro significativo e a viscosidade não.

Escoamento turbulento num tubo (a) tensão de cisalhamento e (b) velocidade média. Sabemos que num escoamento turbulento a tensão de cisalhamento e composta por uma parcela de tensão laminar e uma turbulenta.

( )''vudy

udturblam ρµτττ −=+=

A equação pode ser representada como:

( )dy

udTµµτ +=

Onde µ representa a viscosidade absoluta do fluido e µT representa a viscosidade aparente ou efetiva.

dy

udlmT

2ρµ =



dy

ud

dy

udl

dy

udm

+= 2ρµτ

Visto desta forma podemos colocar que: Sub-camada laminar ou viscosa (região da parede)

Tturblam µµττ >>>>

Camada de amortecedora ou de superposição: Tturblam µµττ ≅≅

Camada turbulenta externa:

µµττ >>>> Tlamturb



3.2 PERFIL DE VELOCIDADES EM ESCOAMENTOS TURBULENTOS – TUBOS LISOS Objetivo:

Determinar o perfil de velocidade num escoamento turbulento numa tubulação considerada lisa e um fluido com propriedades constantes.

Utiliza-se como equação básica a expressão da tensão turbulenta:

Para equacionar o perfil de velocidade e utilizado o conceito de velocidade de atrito. Velocidade de Atrito

ρτWu =*

onde Wτ e a tensão de cisalhamento na parede e ρ a massa especifica do fluido. Alem disto são introduzidas duas grandezas adimensionais:

• Velocidade adimensional +u

• Distancia a partir da parede adimensional +y Velocidade de Adimensional

*u

uu =+

Distancia da parede adimensional

v

yuy

*

=+

rRy −= Representa a distancia normal a partir da parede

ν Viscosidade cinemática do fluido Em termos das camadas: Sub-camada viscosa (região da parede) 5≤+y

Camada de superposição: 305 ≤< +y

dy

ud

dy

udl

dy

udm

+= 2ρµτ



Camada externa: 30>+y

SOLUCAO

1. dy

ud

dy

udl

dy

udm2ρµτ +=

2. sabemos que kylm =

3. ( )2

2

+=

dy

udky

dy

ud ρµτ

4. Dividimos a Equação pela massa especifica

5. ( )2

2

+=

dy

udky

dy

ud

ρµ

ρτ

6. O termo νρµ = representa a viscosidade cinemática:

7. ( )2

2

+=

dy

udky

dy

udv

ρτ

8. ( )dy

ud

dy

udkyv

+= 2

ρτ

9. ( )dy

ud

dy

udkyvu

+= 22*

Equação Deduzir as Equações Especificas.

( )dy

ud

dy

udkyvu

+= 22*



Caso No1 – Escoamento na Sub-camada Viscosa (Região da Parede): Como: turblam ττ >> o termo que representa a tensão viscosa da equação torna-se nulo.

( )dy

udvu =2*

dyv

uud

2*

=

sendo que : *uduud += e que *u

dydy

ν+

= *

2**

u

dy

v

uudu

ν++ = Se obtém:

++ = dydu

Integrando

cyu += ++

Nas condições de contorno na parede: 0 0 por tanto 0 0 ==== ++ yeuuypara e c=0. Lei da Parede Sub-camada Laminar ou Viscosa

++ = yu para 5≤+y



Caso No2 – Camada turbulenta (Região de Externa): Como: lamTrub ττ >> o termo que representa a tensão laminar torna-se nulo.

Podemos utilizar as relações adimensionais:

( )2

22*

=

dy

udkyu

dy

udkyu =*

O termo: *2*

*

*

*

*u

dy

duy

dy

udu

u

vy

u

dy

udu

u

vy

dy

udy +

++

+

++

+

++

=

=

=

νν

** udy

dukyu +

++=

+

++ =

ky

dydu

Integrando:

cyk

u += ++ ln1

Onde c e uma constante que depende da rugosidade da tubulação. Para paredes consideradas lisas, na literatura se encontra c=5 ou também c=5,5. Lei da Logarítmica Camada Externa plenamente turbulenta.

5,5ln5,2 += ++ yu para 30>+y

Pode ser mostrado que integrando a equação anterior se obtém a velocidade media do perfil de velocidades. Lei da Logarítmica Velocidade Media

34,1ln5,2*

+≅ +yu

Vmedia para 30>+y

Caso No3 – Região de Superposição: Neste caso adota-se um perfil de ajuste logarítmico do tipo. Lei da Logarítmica Camada de Superposição.

05,3ln0,5 −= ++ yu para 305 ≤< +y



RESUMO DAS EQUACOES DO PERFIL DE VELOCIDADES Sub-camada viscosa (região da parede) ++ = yu para 5≤+y

Camada de superposição: 05,3ln0,5 −= ++ yu para 305 ≤< +y

Camada externa:

5,5ln5,2 += ++ yu para 30>+y

Perfil de velocidade turbulenta num tubo liso

Escoamento turbulento num tubo (a) tensão de cisalhamento e (b) velocidade média.

y+=u*y/ν

u+



3.3 PERFIL DE VELOCIDADES EM ESCOAMENTOS TURBULENTOS – TUBOS RUGOSOS

• Uma superfície e considerada hidraulicamente lisa quando as saliências da superfície ( ε) ou rugosidade for muito menor que a espessura da sub-camada viscosa (δV)).

• Define-se o parâmetro Parâmetro de Rugosidade

*uνεε =+

• Estudos em tubos em escoamento turbulento utilizando rugosidade areia para aumentar artificialmente a rugosidade permitem concluir que a as superfícies podem ser classificadas em função do parâmetro:

Hidraulicamente Lisa:

• Sem efeito da rugosidade sobre o atrito

50 ≤≤ +ε

Transitórias

• Efeito moderado do numero de Reynolds

705 ≤< +ε

Completamente Rugosa

• A subcamada viscosa e totalmente destruída e o atrito dependem do numero de Reynolds.

70>+ε

Resultados mostram que para escoamento em tubos rugosos, a lei logarítmica da velocidade para escoamento plenamente turbulento e dado por:

Lei da Logarítmica da Velocidade

• Tubos Rugosos • Camada Externa plenamente

turbulenta.

5,8ln5,2 +=+

εy

u para 70>+ε

Integrando esta equação se obtém a velocidade media do perfil de velocidades na tubulação:

Velocidade media

• Tubos Rugosos • Camada Externa plenamente

turbulenta.

5,8ln5,2*

+=εD

u

Vmedia para 70>+ε



3.4 PERFIL DE VELOCIDADES EM ESCOAMENTOS TURBULENTOS – LEI EXPONENCIAL

Uma alternativa para descrever a distribuição de velocidade num escoamento turbulento numa tubulação e dada pela lei exponencial:

Lei Exponencial

n

R

ruu

/1

max 1

−=

onde o expoente n e uma função do numero de Reynolds e da rugosidade do material e varia de 5 a 10.

Para tubos lisos:

Re 4x103 105 106 > 2x106

n 6 7 9 10

Podemos também utilizar uma expressão aproximada. 96.1Relog85.1 −=n

O expoente n esta relacionado com o fator de atrito pela equação empírica: Expoente n

f

n1=



Lembrando que tanto para escoamentos laminar e turbulento o atrito esta relacionada com:

Equações validas para Escoamento Laminar e turbulento

2

8

media

W

Vf

ρτ

= 4

D

L

PW

∆=τ g

V

D

LfhL 2

2

=

• O perfil de velocidades da lei exponencial não poder ser utilizado para determinar a tensão de cisalhamento na parede já que:

• ∞=

parededy

ud

• Para determinar Wτ deve-se relacionar o fator de atrito com a tensão de cisalhamento com as equações

apresentadas acima.

Pode-se obter a velocidade media em função da velocidade máxima pela integração da velocidade:

A

dAu

A

QV

R

media

∫== 0 ( ) max

2

20

)12(1

22)(

unn

n

R

rdrru

V

R

media ++==

∫

π

π

Relação entre a velocidade media e velocidade máxima.

( ) )12(1

2 2

max ++=

nn

n

u

Vmedia




4. Escoamento Viscoso Externo: Conceitos de Camada Limite



4.1 Escoamento Viscoso Externo: Conceitos de Camada Limite

Quando um corpo se move através de um fluido existe um interação entre este e o fluido. Tal interação pode

ser descrita por forças que atuam na interface fluido-corpo. Estas forças se devem aos efeitos viscosos e aos efeitos de pressão. Em Engenharia, para avaliar os efeitos globais é mais interessante representar estas forças em função da denominada força de arrasto que atua na direção do escoamento e a força de sustentação que atua na direção normal ao escoamento denominada sustentação. O arrasto e sustentação podem ser obtidos pela integração das tensões de cisalhamento e as forças normais ao corpo. No Cap.11 são abordadas as forças de sustentação e arrasto para escoamentos externos viscosos sobre superfícies curvas tais como cilindros e aerofólios. No presente capítulo é abordado o escoamento externo sobre placas planas.

4.2 Escoamento em Torno de Corpos A característica do escoamento em torno de um corpo depende de vários parâmetros como: forma do corpo,

velocidade, orientação e propriedades do fluido que escoa sobre o corpo. Os parâmetros mais importantes para descrever o escoamento sobre um corpo são o número de Reynolds e número de Mach.

4.2.1 Efeito do Número de Reynolds no Escoamento Externo

O número de Reynolds (Re= ρ VD/µ) representa a relação entre os efeitos de inércia e os efeitos viscosos. • Sem os efeitos viscosos (µ=0) , o número de Reynolds é infinito. • Por outro lado na ausência de todos os efeitos de inércia (ρ=0) o número de Reynolds é nulo.

Qualquer escoamento real apresenta um número de Reynolds entre esses limites. A natureza do escoamento varia muito se Re >>1 ou se Re <<1. A maioria dos escoamentos que nos são familiares estão associados a objetos de tamanho moderado com comprimento característico da ordem de 0,01m a 10m. As velocidades ao longe destes escoamentos (água e ar) apresentam ordem de grandeza de 0,01m/s até 100m/s. Desta forma o Re destes escoamento está entre 10 < Re < 109. Escoamentos com Re > 100 são controlados por efeitos de inércia. Escoamentos com Re < 1 são controlados por efeitos viscosos. A maioria dos escoamentos são controlados por efeitos de inércia.

Antes de Prandtl a Mecânica dos Fluidos evoluiu com resultados teóricos e experimentais que diferiam. Prandtl introduziu o conceito de camada limite fornecendo o elo entre teoria e prática. Prandtl mostrou que muitos escoamentos viscosos podem ser analisados considerando duas regiões: uma próxima das fronteiras sólidas e outra cobrindo o restante. Apenas na região muito delgada adjacente a fronteira sólida (camada limite) o efeito da viscosidade é importante. Na região fora da camada limite o efeito da viscosidade é desprezível e o fluido pode ser tratado como não-viscoso. Em muitas situações reais, a camada limite desenvolve-se sobre uma superfície sólida plana. Exemplo disso é escoamento sobre cascos de navios e de submarinos, asas de aviões e movimentos atmosféricos sobre terreno plano. Estes casos podem ser ilustrados pelo caso mais simples analisando uma placa plana. Tal caso será estudado a seguir.

Figura 1 Camada limite sobre uma placa plana (espessura exagerada)



4.3 Escoamento sobre Placa Plana

4.3.1 Forças Viscosas Predominantes – Reynolds muito baixo - Re≈0,1 As Fig. .2 a Fig.4 mostram três tipos de escoamentos sobre uma placa plana que tem comprimento L. Para o caso em que Re≈≈≈≈0,1 (Fig..2) os efeitos viscosos são predominantes afetando o escoamento uniforme. Devemos nos afastar consideravelmente da placa plana para alcançar uma região do escoamento que tem sua velocidade alterada em menos de 1%. A região afetada pelos efeitos viscosos é bastante ampla quando o número de Reynolds do escoamento é baixo.

Figura 2 Escoamento sobre placa plana com efeitos viscosos predominantes

4.3.2 Forças Viscosas Moderadas – Reynolds baixo - Re≈10

Com o aumento do Re no escoamento (por ex. aumento de Uoo), neste caso Re≈≈≈≈10, a região onde os efeitos viscosos são importantes se torna menor em todas as direções, exceto a jusante da placa (Fig. 3). Se observa que as linhas de corrente são deslocadas da posição original do escoamento uniforme, mas o deslocamento não é grande como na situação referente ao Re≈≈≈≈0,1.

Figura 3. Placa Plana - Efeitos viscosos moderados



4.3.3 Forças de Viscosas Confinadas – Reynolds Alto - Re≈107

Para escoamento com Re muito alto (Re≈≈≈≈107) predominam os efeitos das forças de inércia. Os efeitos das forças viscosas são praticamente desprezíveis em todos os pontos, exceto naqueles muito próximos da placa plana e na região de esteira localizada a jusante da placa (Fig. 4). Como a viscosidade do fluido não é nula o fluido adere à superfície sólida (condição de não escorregamento). Desta forma a velocidade varia desde zero na superfície da placa até um valor Uoo, na fronteira de uma região muito fina denominada camada limite. Essa região conhecida como camada limite (δ) é sempre muito menor que o comprimento da placa. A espessura desta camada aumenta na direção do escoamento e é nula no borda de ataque da placa. O escoamento na camada limite pode ser laminar ou turbulento. Se define a espessura da camada limite δ como a distância da superfície ao ponto em que a velocidade situa-se dentro de 1% da velocidade de corrente livre.

Figura 4 Placa Plana - Efeitos de inércia importantes

As linhas de corrente fora da camada limite são aproximadamente paralelas àplaca plana. O leve

deslocamento das linhas de corrente externas (fora da camada limite) se deve ao aumento da espessura da camada limite na direção do escoamento e é nula no bordo de ataque da placa. A existência da placa plana tem pouco efeito nas linhas de corrente externas tanto na frente, acima ou abaixo da placa. Por outro lado, a região de esteira é provocada por efeitos viscosos. Camada Limite – Prandtl

O físico alemão Prandtl (1875-1953) realizou um dos grandes avanços na Mecânica dos Fluidos, em 1903, concebendo a idéia da camada limite na qual define – Uma região muito fina dentro da camada limite e adjacente à superfície do corpo onde os efeitos viscosos são muito importantes, onde a componente axial da velocidade varia rapidamente com a distância y. Uma região fora da camada limite denominada região de escoamento potencial onde o fluido se comporta como se fosse um fluido não viscosos, ou investido onde as forças de cisalhamento são desprezíveis. Certamente a viscosidade dinâmica é a mesma em todo o campo de escoamento. Desta forma a importância relativa de seus efeitos (devido aos gradientes de velocidade) é diferente fora e dentro da camada limite. Os gradientes de velocidades normais ao escoamento são relativamente pequenos fora da camada limite e o fluido se comporta como se fosse não viscoso.



4.4 Características da Camada Limite

O tamanho da camada limite e a estrutura do escoamento nela confinado variam muito. Parte desta variação é provocada pelo formato do objeto onde se desenvolve a camada limite. A seguir se analisa o efeito da camada limite para o caso de um fluido viscoso e incompressível sobre uma placa plana de comprimento infinito (x varia de 0 a infinito). Se o Re é muito alto somente o fluido confinado na camada limite sentirá a presença da placa. Exceto na região fora da camada limite a velocidade será essencialmente igual a velocidade de corrente livre V=Ui. Para uma placa finita, o comprimento L pode ser utilizado como comprimento característico. No caso da placa plana de comprimento infinito definimos o Rex= Ux/ν. Se a placa é longa o Re é alto, apresentando uma camada limite exceto na região muito pequena próxima da borda da placa. A presença da placa é sentida em regiões muito finas da camada limite e da esteira.

Figura 10.5. Efeito rotacional de partículas de fluido dentro da camada limite

Consideremos o escoamento de uma partícula de fluido no campo de escoamento. Quando a partícula entra na camada limite começa a distorcer devido ao gradiente de velocidade do escoamento – a parte superior da partícula apresenta uma velocidade maior do que na parte inferior. O elemento de fluido não tem rotação fora da camada limite mas começa a rotar quando atravessa a superfície fictícia da camada limite e entre na região onde os efeitos viscosos são importantes. O escoamento é irrotacional fora da camada limite O escoamento é rotacional dentro da camada limite. A partir de uma certa distância x do bordo de ataque, o escoamento na camada limite torna-se turbulento e as partículas de fluido tornam-se extremadamente distorcidas devido a natureza irregular da turbulência. Uma das características da turbulência é o movimento de misturas produzido no escoamento. Esta mistura é devido a movimentos irregulares de porções de fluido que apresentam comprimentos que variam da escala molecular até a espessura da camada limite. Quando o escoamento é laminar a mistura ocorre somente em escala molecular. A transição do escoamento de laminar para turbulência ocorre quando o Re atinge um valor critico (Rec). Placa Plana: • Rec varia de 2x105 até 3x106 . É função da rugosidade da superfície e da intensidade da turbulência. � Considera-se o valor crítico igual a Rec=5x105. ( 500.000) � Considera-se que a camada limite é turbulenta quando Rex > 3x106 ( 3.000.000) � Laminar Re < 5x105 Turbulento Re > 3,0 x10 6



4.5 Espessura da Camada Limite Na camada limite a velocidade muda de zero na superfície da placa até o valor da velocidade de corrente livre na fronteira da camada limite. Desta forma o perfil de velocidades u=u(x,y) que satisfaz as condições de contorno: V=0 em y =0 e V≈U00 em y =δ. Matematicamente como fisicamente o perfil de velocidade não apresenta nenhuma singularidade. Isto é, u tende a Uoo quando mais nos afastamos da placa (não é necessário que u seja precisamente igual a U00 em y=δ). Se define a espessura da camada limite δ como a distância da superfície ao ponto em que a velocidade situa-se dentro de 1% da velocidade de corrente livre.

∞== U,uy 990 onde δ

4.6 Espessura de Deslocamento A Fig. 10.6 mostra dois perfis de velocidade para escoamento sobre uma placa plana: um (Fig. 6a) considerando perfil uniforme de velocidade (sem atrito) e outro (Fig.6b) com viscosidade no qual a velocidade na parede é nula.

Figura 6. Camada limite e conceito de espessura de deslocamento

Devido à diferença de velocidade U – u dentro da camada limite, a vazão através da seção b – b é menor do que aquela na seção a – a . Se deslocarmos a placa plana na seção a – a de uma quantidade δ* , as vazões pelas seções serão idênticas. Esta distância é denominada espessura de deslocamento.

Figura 7 Perfil de velocidade para definir espessura de deslocamento



A definição é verdadeira se

( )∫∞

−=0

* bdyuUbUδ

onde b é a largura da placa. Desta forma:

∫∞

−=0

* 1 dyU

uδ

A espessura de deslocamento representa o aumento da espessura do corpo necessário para que a vazão do escoamento uniforme fictício seja igual a do escoamento viscoso real. Também representa o deslocamento das linhas de corrente provocado pelos efeitos viscosos. Tal idéia permite simular a presença da camada limite no escoamento pela adição de uma espessura de deslocamento da parede real e tratar o escoamento sobre o corpos mais espessos como se fossem não viscoso.

4.7 Espessura da Quantidade de Movimento A diferença de velocidades existente na camada limite U – u, provoca uma redução do fluxo da quantidade de movimento na seção b – b mostrado na Fig.7 . O fluxo é menor do que aquele na seção a – a da mesma figura. Esta diferença de fluxo de quantidade de movimento na camada limite, também conhecida como déficit do fluxo da quantidade de movimento no escoamento real é dada por:

( ) ( )∫ ∫∞

−=−0

dyuUubdAuUu ρρ

por definição estas integrais representam o déficit do fluxo da quantidade de movimento numa camada limite de velocidade uniforme U e espessura θ. Assim,

( )∫∞

−=0

2 dyuUubbU ρθρ

∫∞

−=0

1 dyU

u

U

uθ

as três definições de espessura de camada limite δ , δ*e θ são utilizadas nas análises de camada limite. A hipótese da camada limite ser fina é essencial para o desenvolvimento do modelo de escoamento. Esta hipótese, na análise do escoamento sobre uma placa plana garante que δ seja muito menor que x (δ <<x) e também que (δ* << x) e (θ << x) onde x é a distância em relação ao bordo de ataque da placa. Como ordem de grandeza se utilizam:

=δ

δδ

8

13

1min

*

turbulento

arla

=δ

δθ

10

17

1min

turbulento

arla



4.8 Coeficiente de Arrasto em Placas Planas O coeficiente de arrasto ou de resistência de um corpo é dado por:

DfDpD CCC +=

CDf representa o coeficiente de tensão de cisalhamento.

AU

FC Df

Df2

2

1∞

=ρ

onde A representa a área superficial ou área molhada. Por exemplo numa placa paralela ao escoamento A=bxL onde b é a largura da placa. O termo CDp representa o coeficiente de arrasto por pressão.

AU

FC Dp

Dp2

2

1∞

=ρ

Neste caso A pode representar projeção num plano normal da área do corpo. Por exemplo num cilindro A=DxL O coeficiente de arrasto total é assim definido:

AU

FC D

D2

2

1∞

=ρ

onde FD= FDp + FDf

No caso de uma placa perpendicular ao fluxo a tensão de cisalhamento não contribui para a força de resistência. O coeficiente de arrasto deve-se unicamente ao arrasto por pressão. Desta forma CD= CDp.

Figura .8 Placa plana perpendicular ao fluxo

No caso de uma placa plana paralela ao escoamento o arrasto deve-se unicamente ao atrito superficial. Desta forma CD= CDf.

CD=CDp



4.9 Coeficiente de Arrasto e Força de Arrasto pela Tensão de Cisalhamento

Considerando que o perfil de velocidade u(x,y) da camada limite seja conhecido. A tensão de cisalhamento τw na parede que atua ao longo da superfície em qualquer posição x é determinado a partir da definição:

0

),(

=∂∂=

y

w y

yxuµτ

Desta forma conhecendo a distribuição de velocidades na camada limite, pode-se determinar a força de cisalhamento, devido ao escoamento que está atuando sobre a superfície sólida. Como a equação anterior não é muito prática para aplicações de Engenharia define-se a tensão de cisalhamento ou força de arrasto local como função do coeficiente de arrasto local Cf. Também denominado coeficiente de tensão de cisalhamento (Cx no texto de Ozisik).

2

2∞=

UC fw

ρτ

onde ρ é a massa específica do fluido e U00 a velocidade de corrente livre. Desta forma conhecendo o coeficiente da tensão de cisalhamento Cf podemos determinar a força de arrasto exercida pelo fluido que está escoando sobre a placa plana. Igualando as equações anteriores se obtém:

0

2

),(2

=∞ ∂∂=

y

f y

yxu

UC

ν

o coeficiente local de arrasto poderá ser determinado se o perfil de velocidade u(x,y) na camada limite for conhecido. O valor médio do coeficiente da tensão de cisalhamento CDf de x=0 até x=L é definido como:

∫=

=L

x

fDf dxCL

C0

1

determinado o CDf podemos calcular a força de arrasto FD atuando sobre a placa de x=0 até x=L numa largura da placa b (lembrando que a área superficial é A=bxL).

2

2∞= U

bLCF DD

ρ

Obs. Para placa plana paralela à direção do escoamento. CD=CDf.



4.10 Equações de BLASIUS – Placa Plana – Camada Limite Laminar

Para o caso de placa plana existem diferentes soluções para determinar a espessura da camada limite, espessura de deslocamento da camada limite, coeficiente de arrasto local e coeficiente de arrasto médio. Em 1908 Blasius, discípulo de Prandtl, obtém a solução exata da camada limite numa placa plana (gradiente de pressão nulo) considerando: � Escoamento em regime laminar. � Escoamento permanente � Escoamento bidimensional � Escoamento incompressível Soluções aproximadas foram também determinadas para tal problema considerando o perfil de velocidades como um polinômio de segundo grau, de terceiro grau e de quarto grau. A seguir são apresentadas as equações denominadas exatas, determinadas por Blasius, válidas para escoamento laminar Re < 5,0x105 até 1,0x106

Figura 9. Esquema de placa plana

Espessura da camada Limite

x

x

Re

5=δ

Espessura de deslocamento da camada limite

∞

=U

vx73,1*δ ou também

xRe

73,1* =δ

4.11 Coeficiente de arrasto local ou Coef. de tensão de cisalhamento

2/1Re

664,0

x

fC =

4.12 Coeficiente de arrasto médio

Podemos determinar o coeficiente de arrasto médio CDf integrando Cf de x=0 até x=L ∫=

=L

x

fDf dxCL

C0

1

2/12/1 Re

328,1

Re

664,02

LL

DfD CC === onde ReL= VL/ν

v

xUx

00Re =

x

L



4.13 Transição de Escoamento Laminar para Turbulento

O valor de Re de transição é uma função muito complexa de vários parâmetros como rugosidade da superfície, curvatura da superfície e intensidade das perturbações existentes no escoamento. No caso do ar a transição de escoamento laminar para turbulento, na camada limite sobre uma placa plana, ocorre para Rec na faixa de 2x105 a 3x106. Para efeitos práticos utiliza o valor fixo Rec=5x105 que na verdade corresponde ao limite inferior da região de transição. O processo de transição envolve instabilidade do campo de escoamento. Pequenas perturbações impostas sobre a camada limite, como vibrações na placa, rugosidade da superfície, pulsações no escoamento principal aumentam ou diminuem a instabilidade dependendo do lugar onde a perturbação for introduzida: � Se a perturbação ocorre em Rex < Rec são amortecidas fazendo com que a camada limite retorne ao regime

laminar � Se a perturbação ocorre em Rex > Rec irão crescer transformando o escoamento em regime turbulento. � A mudança do escoamento laminar para turbulento também provoca uma mudança na forma do perfil de

velocidades.

Figura 10 Perfis de velocidades em placa plana - regime laminar, transição e turbulento (ar).

Observa-se na Fig.10 que o perfil turbulento de velocidades é mais plano apresentando um alto gradiente de velocidade na parede. Trata-se do escoamento de ar com uma velocidade de corrente livre de 27m/s. • Numa placa plana a camada limite será sempre turbulenta para Re > 4,0 x106



4.14 Camada Limite Turbulenta em Placa Plana

4.14.1 Coeficiente de Arrasto Local - Coeficiente de Tensão de Cisalhamento

A partir de dados experimentais Schilichting apresentou a seguinte correlação para o coeficiente de arrasto

local para placa plana lisa.

7x

52,0 10Re5x10 para do váli Re0592,0 <<= −xfC

para número de Reyndols altos, recomenda-se a correlação de Schultz-Grunow:

( ) 9x

7-2,584 10Re10 para válido Relog370,0 <<= xfC

4.14.2 Coeficiente de Arrasto Médio

Para uma camada limite que é inicialmente laminar e passa por uma transição em algum ponto sobre a placa plana, o coeficiente de arrasto turbulento deve ser ajustado para levar em conta o escoamento laminar no comprimento inicial.

Figura 11. Placa plana com região laminar, transição e turbulenta.

Consideremos um escoamento na camada limite sobre uma placa plana que seja:

Laminar na região entre 0 ≤ x ≤ xc e Turbulenta na região xc < x ≤ L . Coeficiente de Arrasto Local Os coeficientes locais de atrito em cada uma das duas regiões são:

(laminar) xx0 em Re664,0 c5,0 ≤≤= −

xfc - Eq. de Blasius

o)(turbulent xem Re0592,0 c2,0 Lxc xf ≤<= − - Eq. de Schilichting

Coeficiente de Arrasto Médio O coeficiente de arrasto médio CDf (igual a CD em placa plana) na região inteira

+

=

+=

∫∫

∫∫

−−

−−

dxxv

Udxx

v

U

Lc

dxcdxcL

c

L

xc

xc

D

L

xc xTurb

xc

xLamD

2,0

2,0

00

0

5,0

5,0

00

0

059,0664,01

1

efetuando a integração se obtém:



L

ccLDC

Re

Re328,1Re074,0Re074,0

5,08,02,0 −−= − válido para Rec ≤ ReL < 107 ,

Definição do Número de Reynolds Total e Crítico e Local

v

LUL

00Re = é o número de Reynolds para o comprimento total (L) da placa plana.

v

xU cc

00Re = é o número de Reynolds crítico da transição do escoamento laminar para turbulento

v

xUx

00Re = é o número de Reynolds crítico em qualquer posição da placa

Forma Geral do Coeficiente de Atrito Médio O coeficiente médio de arrasto CD sobre a região onde o escoamento é parcialmente laminar e parcialmente turbulento, depende do valor do número de Reynolds crítico Rec . Por isto a Eq. anterior é especificada de maneira mais compacta.

7Lc

2,0 10Re Re para válido Re

Re074,0 <<−= −

LLD

BC

Válida quando existe a camada limite turbulenta com camada laminar anterior. O termo B é dada como:

( )arLaDoDTurbulent CCcB minRe −=

o qual depende do número de Reynolds crítico ( Rec ) e das características do arrasto plenamente laminar e turbulento. Para diversos número de Rec o valor de B é dado a seguir.

=

=

=

=

=

6c

5c

5c

5c

101Re para 3340

105Re para 1740

103Re para 1050

102Re para 700

x

x

x

x

B

• No caso em que B=0 corresponde ao escoamento turbulento começando desde x=0 e desta forma se utiliza e

equação para regime turbulento denominada Eq. de Karman-Prandtl:

72,0 10ReRe para Re074,0 <<= −LcLD C

Para altos número de Reynolds 107 < Re < 109 se utiliza a seguinte expressão:

( ) LL

DCRe

1710

Relog

455,058,2

−= para 107 < Re < 109



4.15 Espessura da Camada Limite - Escoamento Turbulento Utilizando uma relação empírica para a tensão de cisalhamento na parede na forma:

2

5/1

0

0296,0 ∞∞=

=

∂∂= U

xU

v

y

u

y

w ρµτ

são obtidas expressões que permitem avaliar a espessura da camada limite turbulenta para placa plana. 1. Para a camada limite plenamente turbulenta, começando da borda de ataque da placa (x=0).

5/1Re381,0)( −= xx

xδ

2. No caso em que a espessura da camada limite é laminar até o ponto em que Rec=5x105, e então se torna

plenamente turbulenta.

7515/1 10Re5x10 para válidaRe10256Re381,0)( <<−= −−

xxxx

xδ

A Fig. 12 apresenta graficamente um resumo os tipos de coeficiente de arrasto médios para placa plana lisa.

Figura 12 Coeficiente de arrasto para placa plana lisa.



4.16 Resumo das Equações da Camada Limite em Placa Plana

Como ser observa existe uma grande quantidade de equações que podem ser utilizadas para avaliar o coeficiente de arrasto médio em placas planas. O uso de cada equação dependerá do regime de escoamento. A placa poderá apresentar escoamento plenamente laminar, escoamento plenamente turbulento ou se na placa plana existe uma região com escoamento laminar e posteriormente uma região com escoamento turbulento. A seguir, para simplificar, podemos utilizar as seguintes relações em exercícios específicos. I - Camada Limite Laminar Perfis de Velocidade

Linear δy

U

u = Parabólico 2

2

−=δδyy

U

u Senoidal

=δ

π y

U

u

2sen

Equação de Blasius 3x105 < Rex < 5x105

2/1Re

664,0

x

fC = 2/1Re

328,1

L

DC =

x

w

U

Re

332,0 2ρτ =

xx

x

Re

5)( =δ

xx

x

Re

73,1)(*

=δ )(346,0)(* xx δδ = )(

7

1)( xx δθ =

II - Camada Limite Turbulenta (escoamento turbulento desde a borda de ataque)

Perfis de Velocidade Exponencial 7/1

=δy

U

u

Equação de Kárman – Prandtl (Rec=5x105)

72,0 10ReRe para Re074,0 <<= −LcLD C

75/1 10ReRe para Re0594,0 <<= −

Lcxf C

Equação de H. Schlichting

( )9

L7

58,210 Re 10 para

Relog

455,0 <<=L

DC

4/120233,0

=δ

νρτU

Uw

5/1Re381,0)( −= xx

xδ

8

)()(* x

xδδ = )(

72

7)( xx δθ =

III – Camada Limite Turbulenta com Camada Laminar Anterior

7L

52,0 10Re 5x10 para Re

1700Re074,0 <<−= −

LLDC

( )9

L7

58,210Re 10 para

Re

1700

Relog

455,0 <<−=LL

DC

7515/1 10Re5x10 para Re10256Re381,0)( <<−= −−

xxxx

xδ



Exemplo - Comparação das Variáveis para Camada Limite Laminar e Turbulenta.

Água (ρ=1000kg/m3 e ν=1x106 m2/s) escoa com velocidade de U=1,0 m/s sobre uma placa plana de L=1m. Avalie a espessura da camada limite δ(x), a espessura de deslocamento δ*(x) e espessura de quantidade de movimento θ(x) e a tensão de cisalhamento na parede τw(x) para x=L. (a) Considere que é mantido escoamento laminar em toda a placa. (b) Considere que a camada limite é provocada, de modo que se torna turbulenta a partir na borda de ataque. 1. Determinamos o número de Reynolds para x=L.

νUL

L =Re 66

10101

11Re == −x

xL

(a) Considerando Equações para Regime Laminar

xx

x

Re

5)( =δ

xx

x

Re

73,1)(*

=δ )(

7

1)( xx δθ =

mmx

510

516

==δ mmxx 73,15346,0)(* ==δ mmx 71,07

5)( ==θ

000664,010

664,0

Re

664,06

===x

fC ( )

2

2

3

2

332,02

11000000664,0

2 m

N

m

kgx

UC fw === ∞ρτ

(b) Considerando Equações para Regime Turbulento

5/1Re381,0)( −= xx

xδ

8

)()(* x

xδδ = )(

72

7)( xx δθ =

( ) mmxx 2410381,01)(5/16 == −δ mmx 3

8

24)(* ==δ mmx 34,224

72

7)( ==θ

( ) 00375,0100,0594 Re0594,05/165/1 === −−

xfC ( )

2

22

875,12

1100000375,0

2 m

Nx

UC fw === ∞ρτ

Espessura da camada limite 8,45

24

)min(

)( ==arLa

Turbulento

δδ

Espessura de deslocamento da camada limite 73,173,1

3

)min(

)(*

*

==arLa

Turbulento

δδ

Espessura da quantidade de movimento 86,371,0

74,2

)min(

)( ==arLa

Turbulento

θθ

Tensão de cisalhamento na parede. 63,5332,0

87,1

)min(

)(==

arLa

Turbulento

w

w

ττ

Obs: Existe um crescimento maior das variáveis na camada limite turbulenta devido a uma tensão de cisalhamento na parede mais alta.




5. ESCOAMENTOS EXTERNOS - CAMADA LIMITE EM PLACA PLANA



5.1 ESCOAMENTOS EXTERNOS: CAMADA LIMITE EM PLACA PLANA

• REGIME LAMINAR < Re 5x105 • TRANSICAO ReC 5x105 (Pode variar entre 2x105 ate 3x106 segundo tipo de rugosidade) • TURBULENTA > Re 3x106

EQUACAO INTEGRAL DA QUANTIDADE DE MOVIMENTO DE VON KARMAN

Estudamos o escoamento numa placa plana lisa submetida a uma velocidade de corrente livre U∞ paralela a placa. A placa não bloqueia o escoamento sendo a única resistência ao escoamento e dado ao cisalhamento. Devido a condição de não escorregamento provoca uma desaceleração brusca das partículas do fluido e estas retardam as partículas vizinhas surgindo uma espessura de camada cisalhante desacelerada denominada camada limite de espessura )(xy δ= . Para determinar a forca de arrasto sobre a placa deve-se realizar a integração das tesões viscosas ao longo da parede. Largura da placa: b Comprimento da placa: L Espessura da camada limite: δ Velocidade de corrente livre: ∞U

Volume de controle: (1) Região de entrada do fluido no VC : )0,0(),( =yx ate ),0(),( hyx =

• Nesta região existe a uma velocidade de corrente livre iUV ˆ1 ∞=r

(2) Região de linha de corrente externa no VC : ),0(),( hyx = ate ),(),( hLyx =

• Nesta região não existe fluido atravessando as fronteiras: 02 =Vr

(3) Região de saída do fluido no VC : ),0(),( Lyx = ate )0,(),( Lyx =

• Nesta região existe a uma velocidade de corrente livre iyuV ˆ)(3 =r

(4) Região sobre a placa plana no VC : )0,0(),( =yx ate )0,(),( Lyx =

• Nesta região não existe fluido atravessando as fronteiras: 04 =Vr



1.1 Aplicando a Eq. da Conservação da Massa:

0=+∀∂∂

∫∫SCVC

AdVdt

rrρρ

Considerando regime permanente:

( ) ( ) ( ) ( )∫∫∫∫∫ +++=43

3

2

2

1

1

AAAASC

AdVAdVAdVAdVAdVrrrrrrrrrr

ρρρρρ

( ) ( )∫∫∫ +=3

3

1

1

AASC

AdVAdVAdVrrrrrr

ρρρ

0)(0

=+− ∫∞

δρρ bdyyubhU

Relação de velocidades.

∫=∞

δ

0)( dyyuhU

( 1)

1. 2 Aplicando a Eq. da quantidade de movimento. – Determinação da Forca de Arrasto.

∫∫∑ +∀∂∂=

SCVC

AdVVdVt

Frrrrr

ρρ

Aplicando na direção x em regime permanente:

∫∑ =SC

x AdVuFrr

ρ

As únicas forcas agindo são as forcas de superfície por cisalhamento. A pressão e constante e desta forma a forca de pressão resultante e nula.

Ax FF −=∑

( ) ( ) ( ) ( )∫∫∫∫ +++=−43

3

2

2

1

1

AAAA

A AdVuAdVuAdVuAdVuFrrrrrrrr

ρρρρ

( ) ( )∫∫ +=−3

3

1

1

AA

A AdVuAdVuFrrrr

ρρ

bdyuuAuuFA 2

0

2111 ρρδ

∫+−=−

Forca de Arrasto local FA(x) dyyubbhUF oA ∫−=

δ

ρρ0

22 )(

( 2)

Obs. Trata-se da forca de arrasto para uma posição x da placa plana FA(x) já que a espessura da camada limite depende de x. )(xδδ =



1. 3 Espessura da Quantidade de Movimento. Arranjando a Eq. (2) na forma:

( ) dyyubhUbUFA ∫−= ∞∞

δ

ρρ0

2 )(

Substituindo (I) em (II)

dyyubdyyubUFA ∫∫ −= ∞

δδρρ

0

2

0)()(

dyuUubF

dyuuUbF

dyyudyyuUbF

A

A

A

)(

)(

)()(

0

2

0

0

2

0

−=

−=

−=

∫

∫

∫∫

∞

∞

∞

δ

δ

δδ

ρ

ρ

ρ

∫

∫

∞∞

∞∞

∞

−=

−=

δ

δ

ρ

ρ

0

0

)1(

)(

dyU

uuUbF

dyuUU

UubF

A

A

−= ∫∞∞

∞

δρ

0

2 )1( dyU

u

U

ubUFA

O termo entre parêntesis e denominado Espessura da Quantidade de Movimento: Espessura da Quantidade de Movimento.

∫∞∞

−=δ

θ0

)1( dyU

u

U

u

(3)

Desta forma a forca de arrasto e dada por: Forca de Arrasto Local.

θρ 2∞= bUFA

(4)



1. 4 Tensão de Cisalhamento na Parede Sendo a espessura da parede considerada como uma medida do arrasto total da placa. Von Karman notou que o arrasto também e equivalente a integral da tensão de cisalhamento.

∫=x

WA dxxbxF0

)()( τ

)(xbdx

dFW

A τ=

dx

dbU

dx

dFA θρ 2∞=

Tensão de cisalhamento na parede de uma Placa Plana. Valida para escoamento Laminar e Turbulento.

dx

dUx

θρτ 2)( ∞=

(5)

5.2 RESULTADOS PARA ESCOMANETO LAMINAR: 2.1 ) Espessura do Momento da Quantidade de Movimento Para o escoamento laminar Von karman considerou que o perfil de velocidades tivesse um formato aproximadamente parabólico ajustado pela expressão:

−= ∞ 2

22),(

δδyy

Uyxu Valido para ),(0 yxy δ≤≤

(6)

Utilizando esta expressão na definição da espessura da quantidade de movimento:

∫∞∞

−=δ

θ0

)1( dyU

u

U

u

Se obtém:

∫ ≅

+−

−=

δδ

δδδδθ

0 2

2

2

2

15

221

2dy

yyyy

Espessura do Momento da Quantidade de Movimento . Solução de Von Karman. Regime Laminar

)(15

2)( xx δθ =

(7)



2.2 Espessura da Camada Limite Laminar Da mesma forma podemos determinar a tensão de cisalhamento na parede:

δµµτ ∞

=

≅∂∂=

U

y

ux

y

W

2)(

0

Igualando as expressões da tensão de cisalhamento na parede:

dx

dUxW

θρτ 2)( ∞= e δ

µτ ∞=U

x2

)(

com:

dx

d

dx

d δθ15

2=

se obtém:

δµδρ ∞

∞ =U

dx

dU

2

15

22

dxU

d∞

= νδδ 15

Integrando de 0 a x, considerando 0=δ em x=0.

xU ∞

= νδ 152

1 2

Espessura da Camada Limite. Solução de Von Karman. Regime Laminar

xx

x

Re

5,5)( =δ

(8)

Esta solução da espessura da camada limite e somente 10% maior que a solução exata da espessura da camada limite numa placa plana em regime laminar.



2.3 Espessura de deslocamento da Camada Limite Laminar.

Utilizando a figura observa-se que a linha de corrente externa desvia-se uma distancia )(* xδ para satisfazer a conservação da massa entre a entrada e saída.

0=+∀∂∂

∫∫SCVC

AdVdt

rrρρ

Considerando fluido incompressível em regime permanente.

∫=∞

δρρ

0)( bdyyubhU

onde: )()( * xhx δδ += . Cancelando b e ρ e substituindo esta expressão na anterior:

( )

( )

∫

∫

∫∫

∫

∫

∫

−=

−=

−=

−=

=−

=−

∞

∞∞

∞∞

∞∞

∞∞

∞

δ

δ

δδ

δ

δ

δ

δ

δ

δ

δδ

δδ

δδ

0

*

0

*

00

*

0

*

0

*

0

*

1 dyU

u

dyuUU

udydyUU

udyUU

udyUU

udyU

Espessura de Deslocamento da Camada Limite.

∫

−=

∞

δδ

0

* 1 dyU

u

(9)



2.4 Espessura de deslocamento da CL - Laminar. Utilizando o perfil de velocidades para camada limite laminar:

−=

∞2

22),(

yyU

yxu δδ

∫

−=δ

δ0

* 1 dyU

u

∫

+−=

δ

δδδ

0 2

2* 2

1 dyyy

3

332

2

21

*

02

32*

0 2

2*

δδ

δδδδδ

δ

δδδ

δ

δ

=

+−=

+−=

+−= ∫

yyy

dyyy

Espessura de Deslocamento da Camada Limite. Solução Aproximada. Regime Laminar

3

)()(* x

xδδ =

(10)

Utilizando a aproximação de Von Karman para a espessura da CL

xx

x

Re

5,5)( =δ

Espessura de Deslocamento da Camada Limite. Solução Aproximada. Regime Laminar

xx

x

Re

83,1)(*

=δ

(11)



2.5 Coeficiente de Arrasto Local CL – Laminar A tensão de cisalhamento na parede esta relacionada com o coeficiente de arrasto superficial local por:

fW CUx 2

2

1)( ∞= ρτ

Igualando com a expressão da tensão de cisalhamento na parede:

δµτ ∞=

UxW

2)(

δµρ ∞

∞ =U

CU Df

2

2

1 2

Considerando nesta equação a solução da espessura da camada limite:

x

xRe

5)( =δ

x

f

UCU

Re

52

2

1 2 ∞∞ =

µρ

xx

xfC

Re

73,0

Re

Re

5,5

4 ==

Coeficiente de Arrasto Superficial Local. Solução de Von Karman Regime Laminar

x

fCRe

73,0=

(12)



2.6 Coeficiente de Arrasto Total CL Laminar

∫=L

WA dxxbLF0

)()( τ

fW CUx 2

2

1)( ∞= ρτ

x

fCRe

73,0=

∫ ∞=L

fA dxCUbLF0

2

2

1)( ρ

∫ ∞=L

x

A dxUbLF0

2

Re

73,0

2

1)( ρ

∫−

∞∞=

L

A dxxU

vUbLF

0

2/12

2

173,0)( ρ

LA

LUbL

U

vUbLF

Re

73,02

2

173,0)(

22/12 ∞

∞∞ ==

ρρ

bLCULF DA2

2

1)( ∞= ρ

L

DCRe

73,02=

)(2 LCC fD =

Coeficiente de Arrasto Superficial Total Solução de Von Karman. Regime Laminar

L

DCRe

46,1=

(13)

Também:

∫=L

fD dxxCL

C0

)(1



3. DEMOSTRACAO:

∫ ≅

+−

−=

δδ

δδδδθ

0 2

2

2

2

15

221

2dy

yyyy

Resolvendo primeiro o termo:

+−

−

2

2

2

2 21

2

δδδδyyyy

)(15

2)(

15

2

53

5

54

4

3

5

2

2

54

4

3

5

2

2

452

452

22422

12

4

5

3

4

2

32

4

5

3

4

2

32

0 4

4

3

3

2

2

4

4

3

3

2

2

4

4

3

3

2

2

3

3

2

2

2

2

2

2

xx

yyyy

dyyyyy

yyyy

yyyyyyyyyy

o

δθ

δδδδδθ

δδ

δδ

δδ

δδθ

δδδδθ

δδδδθ

δδδδ

δδδδδδδδδδ

δ

δ

=

=−+−=

−+−=

−+−=

−+−=

−+−

−+−+−=

+−

−

∫

\



5.3 RESUMO DAS EQUACOES DE CAMADA LIMITE EM PLACA PLANA - LAMINAR E TURBULENTO ESCOAMENTO LAMINAR Rec < 5x105 Equação Solução Von Karman Solução Exata (Blasius) Espessura da Camada Limite

xx

x

Re

5,5)( =δ

xx

x

Re

5)( =δ

Espessura de deslocamento da Camada Limite

3

)()(* x

xδδ =

xx

x

Re

83,1)(*

=δ

89,2

)()(* x

xδδ =

xx

x

Re

73,1)(*

=δ

Espessura da Quantidade de Movimento )(

15

2)( xx δθ = )(

7

1)( xx δθ =

Coeficiente de Local de Arrasto

x

fCRe

73,0= x

fCRe

664,0=

Coeficiente de Arrasto médio. fD CC 2=

L

DCRe

46,1=

fD CC 2=

L

DCRe

328,1=

ESCOAMENTO TURBULENTO Rec > 5x105

Equação Turbulento

5x105 < ReL < 107

Turbulento com Laminar Anterior 5x105 < ReL < 107

(Transição)

Placa Plana Rugosa

Espessura da Camada Limite ( ) 5/1Re

381,0)(

xx

x =δ

( ) xxx

x

Re

10256

Re

381,0)(5/1

−=δ

Espessura de deslocamento da Camada Limite

8

)()(* x

xδδ =

8

)()(* x

xδδ =

Espessura da Quantidade de Movimento

)(72

7)( xx δθ = )(

72

7)( xx δθ =

Coeficiente de Local de Arrasto ( ) 5/1Re

0594,0

x

fC = Antes de xc x

fCRe

664,0=

Após xc ( ) 5/1Re

0594,0

x

fC =

5,2

log58,187,2−

+=εx

C f

Coeficiente de Arrasto médio. ( ) 5/1Re

074,0

L

DC = ( ) LL

DCRe

1700

Re

074,05/1

−= 5,2

log62,189,1−

+=εL

CD

Turbulento 107 < ReL < 109

Turb. Com Laminar anterior 107 < ReL < 109

Coeficiente de Arrasto médio. ( ) 58,2Relog

455,0

L

DC = ( ) LL

DCRe

1700

Relog

455,058,2

−=



5.4 RELACOES BASICAS Tensão de cisalhamento na parede. (Escoamento Laminar)

0

)(=∂

∂=y

W y

ux µτ

( 1 )

Tensão de cisalhamento na parede.

)(2

1)( 2 xCUx fW ∞= ρτ

( 2 )

Coeficiente de Arrasto Local.

nx

f

AxC

Re)( =

( 3 )

Forca de Arrasto local.

∫=x

WA dxxbxF0

)()( τ

( 4 )

Forca de Arrasto Local.

)()( 2 xbUxFA θρ ∞=

(5)

Coeficiente de arrasto médio ou total.

∫=L

fD dxxCL

C0

)(1

( 6 )

Forca de Arrasto da Placa.

DA ACUF 2

2

1∞= ρ

( 7 )

Espessura da Quantidade de Movimento.

∫∞∞

−=δ

θ0

)1()( dyU

u

U

ux

(8)


∫

−=

∞

δδ

0

* 1)( dyU

ux

(9)

Espessura da Camada Limite.

nx

A

x

x

Re

)( =δ

(10)



5.5 Coeficiente de Arrasto em Placa Plana – Regime Laminar e Turbulento

Fonte: White Mecânica do Fluidos 4ª Edição. 2002




6. ESCOAMENTOS EXTERNOS - CAMADA LIMITE EM PLACA PLANA

ESCOAMENTOS TURBULENTOS



6.1 ESCOAMENTOS EXTERNOS CAMADA LIMITE EM PLACA PLANA ESCOAMENTO TURBULENTO No caso de escoamento turbulento sobre placa plana não existe uma teoria exata e sem varias aproximações computacionais utilizando vários modelos empíricos. No presente material será adotada uma solução simplificada utilizando equacionamento integral com apoio de equação empírica. Determinação da Espessura da Quantidade de Movimento: Para o escoamento turbulento considera-se como valida o perfil de velocidades exponencial: Perfil de velocidades exponencial: n=7 Re < 107 n=8 107 < Re < 108 n=9 108 < Re < 109

7/1

),(

=∞ δ

y

U

yxu valido para ),(0 yxy δ≤≤

( 1 )

Espessura da Quantidade de Movimento.

∫∞∞

−=δ

θ0

)1( dyU

u

U

u

( 2 )

δδδδδ

δδθ

δδθ

δδθ

δδθ

δ

δ

δ

72

7

9

7

8

7

9

7

8

7

9

7

8

7

)1

7/2

7/9

7/1

7/8

07/2

7/9

7/1

7/8

0

7/27/1

0

7/17/1

=−=−=

−=

−

=

−

=

∫

∫

yy

dyyy

dyyy

Espessura da Quantidade de Movimento. Escoamento Turbulento.

)(72

7)( xx δθ =

( 3 )



Determinação da Espessura da Camada Limite No material de Fox e Macdonald se utiliza uma expressão empírica com base a resultado de escoamento em tubulações representando a tensão de cisalhamento na parede dada por:

4/1

20233,0)(

=

∞∞ δ

ρτU

vUxW

A tensão de cisalhamento na parede e relacionada com a espessura de deslocamento quantidade de movimento.

dx

dUx

θρτ 2)( ∞=

Foi determinado o termo que para escoamento turbulento:

dx

x

dx

xd )(

72

7)( δθ =

Igualando as expressões da tensão de cisalhamento:

dx

dU

U

vU

δρδ

ρ72

70233,0 2

4/1

2∞

∞∞ =

dx

d

U

v δδ 72

70233,0

4/1

=

∞

dx

d

U

v δδ 72

710233,0

4/1

4/1

=

∞

δδ ddxU

v 4/1

4/1

0233,07

72 =

∞

dxU

vd

4/1

4/1 0233,07

72

=

∞

δδ

Integrando:

∫∫

=

∞

xdx

U

vd

0

4/1

0

4/1 0233,07

72δδδ

ctexU

v +

=

∞

4/1

4/5 0233,07

72

5

4δ



Considerando que para 0=δ 0=x se obtém que cte=0.

5/1

5/1

5/4

5/1

5/44/1

Re

382,0)(

382,0)(

382,0)(

0233,07

72

4

5)(

xx

x

xU

v

x

x

xU

vx

xU

vx

=

=

=

=

∞

∞

∞

δ

δ

δ

δ

Espessura da Camada Limite Escoamento Turbulento.

5/1Re

382,0)(

xx

x =δ

( )



Determinação do Coeficiente Local da Atrito: A tensão de cisalhamento na parede e dada pelas expressões:

4/1

20233,0)(

=

∞∞ δ

ρτU

vUxW

fW CUx 2

2

1)( ∞= ρτ

Igualando as e explicitando o coeficiente local do arrasto na placa:

4/1

22 0233,02

1

=

∞∞∞ δ

ρρU

vUCU f

4/1

0466,0

=

∞δU

vC f

δ∞

=

U

vC f

4

0466,0

Substituindo a expressão da espessura da camada limite:

5/1Re

382,0)(

xx

x =δ

=

∞ 382,0

Re

0466,0

5/14

xf

xU

vC

382,0

1

Re

Re

0466,0

5/14

x

xfC=

382,0

Re

0466,0

5/44 −

=

xfC

Extraindo a rais quarta:

( )5/1

4/1

5/1

Re0593,0382,0

Re0466,0 −

−

== xx

fC

Espessura da Camada Limite Escoamento Turbulento.

5/1Re0593,0 −= xfC

( )



Determinação do Coeficiente de Arrasto Médio:

dxCL

CL

x fDf ∫ ==

0

1

( )dxL

CL

x xDf ∫ =

−=0

5/1Re0593,01

( )dxxU

v

LC

L

xDf ∫ =

−−

∞

=

0

5/1

5/11

0593,0

5/4

5/1

4

510593,0 L

U

v

LCDf

−

∞

=

LLU

v

LCDf

5/1

5/1

4

510593,0 −

−

∞

=

5/1

4

50593,0

−

∞

=

LU

vCDf

( ) 5/1Re074,0 −= LDfC

Coeficiente de Arrasto Médio. Escoamento Turbulento.

( ) 5/1Re074,0 −= LDfC

( )



(C) Espessura de deslocamento da CL.


∫

−=

∞

δδ

0

* 1 dyU

u

( )

Utilizando a equação da distribuição da velocidade para camada limite turbulenta:

7/1),(

=∞ δ

y

U

yxu

Pode ser mostrado que: Espessura de Deslocamento da Camada Limite.

8

)()(* x

xδδ =

( )




7. Escoamento Viscoso Externo:

Forças Aerodinâmicas



Escoamento Viscoso Externo: Forças Aerodinâmicas

7.1 Forças Aerodinâmicos de Sustentação e Arrasto Num escoamento externo quando o corpo se movimento através do fluido se manifesta uma interação fluido-corpo resultando em forças que podem ser descritas em função da tensão de cisalhamento na parede (τw) provocada pelos efeitos viscosos e uma tensão normal provocada pela distribuição de pressão (p).

Figura 1 Forças aerodinâmicas sobre um corpo

• A componente da força resultante que atua na direção normal ao escoamento é denominada força de sustentação

(Lift, L ou FL). • A componente da força resultante que atua na direção do escoamento é denominada força de arrasto. (Drag, D ou

FD) . Consideremos um elemento diferencial localizado na superfície do corpo em estudo. As componente x e y da força que atua no pequeno elemento de área dA são:

θτθ coscos dApdAdF wx +=

θτθ cossen dApdAdF wy +−=

O arrasto e a sustentação podem ser determinados pela integração das tensões de cisalhamento e das tensões normais ao corpo. A força de sustentação é dada por:

θτθ cossen dApdAdFF wyL +−== ∫

A força de arrasto é dada por:

∫∫∫ +== dAdApdFF WxD θτθ sencos

Para determinar esta força é necessário determinar o formato do corpo e as distribuições da tensão de cisalhamento na parede e da distribuição de pressão ao longo da superfície do corpo.



5.5.1 Coeficiente de Arrasto

Na forma adimensional esta força é definida pelo coeficiente de arrasto como:

AU

FC D

D2

2

1∞

=ρ

onde

O coeficiente de arrasto ou de resistência de um corpo é dado por:

DfDpD CCC +=

onde CDf representa o coeficiente de tensão de cisalhamento.

AU

FC Df

Df2

2

1∞

=ρ

A representa a área superficial ou área molhada. Por exemplo, numa placa paralela ao escoamento A=bL onde b é a largura da placa e L o comprimento da placa. O termo CDp representa o coeficiente de arrasto por pressão.

AU

FC Dp

Dp2

2

1∞

=ρ

Neste caso A pode representar projeção num plano normal da área do corpo. Por exemplo num cilindro A=DL , onde D é o diâmetro do cilindro No caso de uma placa perpendicular ao fluxo a tensão de cisalhamento não contribui para a força de resistência. O coeficiente de arrasto deve-se unicamente ao arrasto por pressão. Desta forma CD= CDp.

Figura 2 Placa plana perpendicular ao fluxo

Como foi visto no Cap.10, no caso de uma placa plana paralela ao escoamento, o arrasto se deve unicamente ao atrito superficial. Desta forma CD= CDf.

CD=CDp



Na Tab.1 se são dados os valores do coeficiente de arrasto para diferentes corpos rombudos entre eles,

esferas, semi-esferas, cilindros, placas planas, aerofólios; também é dado o coeficiente de arrasto de corpos típicos como asas de avião e automóveis. Cabe salientar que estes são valores de referência. Um estudo mais apurado deverá ser realizado para projetos de sistemas específicos.

Tabela 1 Coeficiente de Arrasto para diferentes tipos de corpos Corpos rombudos CD Esfera rugosa 0.40 Esfera lisa 0.10 Semi-esfera oca oposta à corrente 1.42 Semi-esfera oca com face para a corrente 0.38 Semi-cilindro oco oposto a corrente 1.20 Semi-cilindro oco com face para a corrente 2.30 Placa plana 90° 1.17 Placa plana comprida a 90° 1.98 Roda girando oca h/D=0.28 0.58 Corpos afinados CD Placa Plana Laminar 0.001 Placa Plana Turbulenta 0.005 Aerofólio valor mínimo 0.006 Aerofólio próximo do estol 0.025 Asa em escoamento subsônico mínimo 0.05 Automóveis CD Avião de transporte subsônico 0.016 Avião supersônico M=2.5 0.025 Barcos 0.4-1.2 Helicópteros 0.3 -0.4 Carro de esporte 0.4 -0.5 Carro Econômico 0.5 Camioneta e caminhão 0.6-0.7 Trator e Trailers 0.7-0.9 Pessoas CD Homem em pé 1.0 - 1.3 Esquiador 1.2 - 1.3 Skier 1.0 - 1.1 Outros Fios e cabos 1.0 - 1.3 Prédio Empire State 1.3 - 1.5 Torre de Eiffel 1.8 - 2.0



7.2 Escoamento sobre cilindros - Efeito da viscosidade Número de Reynolds Muito Baixo Para Reynolds baixo (Re < 0.1) o escoamento apresenta uma grande região onde os efeitos viscosos são importante. As linhas de corrente são praticamente simétricas com comportamento muito similar na parte anterior e posterior do cilindro. Este tipo de escoamento pode ser estudado utilizando a teoria de escoamentos potenciais.

Figura 3 Escoamento com baixo Re

Número de Reynolds Moderado Para escoamento em regime moderado (Re≅50) a região onde os efeitos viscosos são importantes se torna menor a montante do cilindro. A jusante a região viscosa aumenta. O escoamento perde sua simetria. Forma-se uma bolha de separação atrás do cilindro existindo um escoamento em sentido contrário ao fluxo principal.

Figura 4 Escoamento para Re moderado

Número de Reynolds Alto No caso de escoamento com número de Reynolds alto (Re > 105) a área afetada pelas forças viscosas é concentrada na parte de atrás do cilindro. Na parte frontal do cilindro se desenvolve uma camada muito fina de fluido onde os efeitos viscosos são importante. Na parte frontal, após a separação, o escoamento torna-se turbulento originando-se uma região com emissão de vórtices.

Figura 5 Escoamento para Re alto Nestas regiões o fluido apresenta gradientes consideráveis de velocidade. Como a tensão de cisalhamento é proporcional a estes gradientes, os efeitos viscosos são significativos. Fora da camada limite e da região de vórtices o fluido se comporta como se fosse um fluido não viscoso. Cabe salientar que a viscosidade dinâmica permanece a mesma em todo o campo do escoamento já que o fluido é o mesmo.



7.3 Escoamento não viscoso num cilindro

Do estudo do escoamento da camada limite numa placa plana sabemos que a fronteira da camada limite tenderá ao valor da velocidade de corrente livre (Voo) admitida a jusante da placa. Neste caso aplicando a Eq. de Bernoulli podemos constatar que não existe variação da pressão ao longo da placa. No caso do escoamento sobre um cilindro isto é bem diferente. Consideremos um escoamento não viscoso sobre um cilindro. Neste tipo de escoamento as linhas de corrente formadas em torno do corpo são simétricas e a linha de corrente que atinge o ponto de estagnação contorna o cilindro aderida ao mesmo. Devido à curvatura do cilindro a velocidade do fluido que contorna o cilindro (U) é diferente da velocidade de corrente livre e dependente da posição angular. Neste caso aplicando a Eq. de Bernoulli pode ser constatado que existe uma variação da pressão dependente da variação da velocidade que contorna o cilindro.

Figura. 6 Esquema de escoamento não viscoso

Consideremos que a montante do cilindro a corrente livre não perturbada apresenta uma velocidade Voo e uma pressão Poo. Podemos aplicar a Eq. de Bernoulli que contorna o cilindro considerando um ponto a montante do cilindro e outro sobre a superfície da mesma com pressão p e velocidade U=U(θ).

g

U

g

p

g

V

g

p

22

22

+=+ ∞∞

ρρ

Para analisar a distribuição de pressão utilizamos na forma adimensional definindo o coeficiente de pressão (Cp):

2

2

1∞

∞−=

V

ppcp

ρ

Explicitando o termo (p - poo) da Eq. de Bernoulli e substituída na Eq. de Cp se obtém:

2

1

−=

∞V

Ucp

A equação obtida mostra a dependência da distribuição de pressão em função da velocidade do fluido que contorna o cilindro. Para escoamento não viscoso a solução teórica (potencial) da distribuição de pressão é dada como:

θ2sen41−=pc

Da mesma forma a velocidade ao longo da superfície é dada por:



θθ sen2)( ∞= VU

a qual pode ser representada na forma adimensional

θθθ sen2)(

)(* ==∞V

UU

Utilizando estas expressões podemos graficar a distribuição de Cp e do perfil de velocidades em torno do

cilindro. A pressão é simétrica em relação ao semi-plano vertical atingindo seu máximo nos pontos de estagnação A e F. Observa-se que a velocidade nos pontos de estagnação (θ=0 e θ=1800) é nula (U*(θ) =0), alcançando seu máximo em θ=900 sendo sua magnitude o dobro da velocidade de corrente livre (U=2Voo).

Figura 7 Distribuição do coeficiente de pressão e da velocidade tangencial

Considerando o escoamento não viscoso o arrasto por atrito será nulo (CDf =0). Devido à simetria da distribuição de pressão em torno ao cilindro o arrasto por pressão é nulo (CDp=0). Dados experimentais mostram que sempre existirá um arrasto no cilindro mesmo tratando-se de fluidos com viscosidade muito pequena. Isto nos leva ao denominado Paradoxo de d´Alambert o qual especifica que o arrasto num corpo é sempre nulo para escoamento não viscoso, porém o arrasto num corpo imerso num fluido viscoso não é nulo.

Cilindros : Escoamento não viscoso CDp = CDf=0



7.4 Escoamento viscoso num cilindro : Efeito do Gradiente Adverso de Pressão

Numa placa plana paralela ao escoamento a camada limite se desenvolve num campo de escoamento onde a pressão permanece constante. Isto significa que o gradiente de pressão é nulo. No caso de geometrias mais complexas, ou placa plana com inclinação, o campo de pressão deixa de ser uniforme. No caso de um cilindro na camada limite se desenvolve um gradiente de pressão devido à variação da velocidade da corrente livre que contorna a fronteira da camada limite.

Consideremos uma partícula de fluido, que escoa dentro da camada limite, que viaja do ponto A para o ponto F. Tal partícula está submetida à mesma distribuição de pressão das partículas de fluido próximas, porém fora da camada limite. Contudo, devido aos efeitos viscosos, a partícula localizada dentro da camada limite sofre perdas de energia. Sendo assim a partícula não tem energia suficiente para vencer o gradiente adverso de pressão quando escoa de C para F. Considera-se que a partícula de fluido quando chega em C não tem quantidade de movimento suficiente para vencer o gradiente de pressão adverso. Se define gradiente de pressão adverso quando a pressão aumenta no sentido do escoamento ou ∂p/∂x > 0 Se define gradiente de pressão favorável quando a pressão diminui no sentido do escoamento ou ∂p/∂x < 0

Figura 8 Escoamento com gradiente adverso de pressão sobre um cilindro

Observando o perfil de velocidades dentro da camada limite (Fig. ) vemos que no ponto D, onde ocorre a separação do escoamento, o gradiente de velocidade e a tensão de cisalhamento na parede são nulos. Após este ponto se origina um escoamento reverso dentro da camada limite.

0 e 0 separação de ponto No 0

==

∂∂

=w

yy

u τ

Atualmente as soluções computacionais conseguem identificar nos escoamento viscosos a separação da camada limite e a emissão de vórtices tal como representado na Fig. 11.9.



Figura 9 Solução numérica (CFD) do escoamento num cilindro com emissão de vórtices

Devido aos efeitos da separação da camada limite a pressão média na metade traseira do cilindro é muito

menor que na metade dianteira. Isto origina um arrasto (CD) devido principalmente à parcela de arrasto por pressão (CDp) já que o arrasto por efeitos viscosos (CDf) pode ser muito pequeno. O arrasto por pressão é denominado também arrasto por forma devido a sua dependência da forma do objeto. Cilindros : Escoamento viscoso CDp >> CDf Pelo efeito da separação da camada limite podemos compreender o paradoxo de d`Alambert. No escoamentos sobre um corpo submerso, mesmo para fluido com pequena viscosidade, se manifestará uma força de arrasto, a qual é, geralmente independente da magnitude da viscosidade do fluido. Dependência do Regime de Escoamento A localização do ponto de separação, a largura da esteira de vórtices originados na parte traseira do corpo e a distribuição de pressão na superfície do corpo dependem da natureza do escoamento, seja ele laminar ou turbulento. A energia cinética e a quantidade de movimento associadas ao escoamento na camada limite turbulenta são maiores do que as associadas ao escoamento na camada limite laminar. Isto se deve basicamente ao seguinte: (1) O perfil de velocidade na camada limite é mais uniforme no caso do escoamento turbulento que no caso do

escoamento laminar. (2) A energia associada com os movimentos turbulentos aleatórios é maior que a laminar. Desta forma o descolamento da camada limite turbulenta desenvolvida em torno de um cilindro descola numa posição posterior daquela da camada limite laminar tal como se observa na figura.

(a) Laminar (b)Turbulento

Figura 10 Separação do escoamento laminar e turbulento.



Perfil de velocidades na camada limite no cilindro analisado

Figura 11 Distribuição de pressão em cilindro escoamento não viscoso e viscoso

Efeito do Regime de Escoamento no Arrasto de Esferas e Cilindros Como mostra a Fig. 12 existe uma dependência do coeficiente de arrasto nos cilindros e esferas lisas, muito semelhantes em função do número de Reynolds. Para escoamento com baixo número de Reynolds o arrasto é função de 1/Re. Para escoamentos moderados (103 a 105) o coeficiente de arrasto tem comportamento constante. Quando o Re atinge o valor crítico a camada limite se torna turbulenta e existe um queda abrupta do arrasto. Para determinar o coeficiente de arrasto (CD) numa esfera lisa podemos também utilizar as equações sugeridas por Chow:

Re ≤ 1 1 < Re ≤ 400 400 <Re ≤ 3x105 3x105 < Re ≤ 2x106 Re > 2x106

Re

24=DC ( ) 646,0Re

24=DC

5,0=DC ( ) 4275,0Re

000366,0=DC

18,0=DC



Figura 12 Coeficiente de arrasto em função do número de Reynols para cilindros e esferas

A estrutura típica do escoamento segundo o número de Reynolds é mostrada na Fig. 13. Para baixo número de Reynolds (Re≅0,1) se observa o escoamento típico (A) sem separação. A medida que Reynolds aumenta (Re≅10) se origina uma região de separação na parte traseira do corpo (B). A formação de vórtices oscilantes (C) se origina (Re≅100), conhecidos como vórtices de Von Karman. Para maiores Re se produz a configuração do escoamento laminar (D) no qual o arrasto é quase constante. Posteriormente quando se alcança o Re crítico (≅3x105) o escoamento torna-se turbulento (E) no qual o ponto de separação desloca-se para a parte traseira do perfil originando-se queda brusca do arrasto.

Figura 11.13 Tipos de escoamentos associados aos pontos indicados no gráfico



Efeito da Rugosidade Superficial no Arrasto de Esferas e Cilindros

Geralmente o arrasto aumenta com o aumento da rugosidade superficial nos corpos delgados como os perfis aerodinâmicos. Isto se deve a que o escoamento se torna turbulento. Nesta condições a maior contribuição para o arrasto total se deve ao arrasto por atrito (CDf) que é muito maior no escoamento turbulento que no escoamento laminar. Por outro lado, como se observa na figura abaixo, nos corpos rombudos, como um cilindro circular ou esferas, o aumento da rugosidade superficial pode causar uma diminuição do arrasto total. Para uma esfera lisa quando o Re atinge o valor crítico (Re≅3x105), a camada limite se torna turbulenta. Nesta condição a esteira atrás da esfera fica mais estreita. Isto origina uma diminuição significativa do arrasto por pressão (CDp) e um leve aumento do arrasto por atrito (CDf). A combinação desta duas parcelas de arrasto (CDp + CDf) fornece um arrasto total menor que nas condições de escoamento laminar. O aumento da rugosidade superficial pode conseguir que a camada limite se torne turbulenta para um Re mais baixo e com isto conseguir um arrasto total menor. Esta é, por exemplo, a técnica utilizada nas bolas de golfe que apresentam uma rugosidade artificial exagerada para conseguir um escoamento turbulento com menor Re (≅ 4x104) e diminuir assim o arrasto. Desta forma com uma tacada a bola pode alcançar maiores distâncias percorridas comparadas com o caso de uma esfera lisa.

Figura 14 Efeito da rugosidade no coeficiente de arrasto em esferas lisas

Figura 15 Diferença do escoamento de uma esfera lisa e uma bola de golfe.



7.5 Sustentação Aerodinâmica A sustentação é a componente da força aerodinâmica perpendicular ao movimento do fluido. Tal força é a responsável pelo vôo dos aviões e princípio de acionamento de muitos tipos de turbomáquinas. Nos aviões, por exemplo, as asas apresentam um formato aerodinâmico (Fig) cuja seção é denominado aerofólio ou perfil aerodinâmico. Estes são projetados para produzir sustentação com a menor força de resistência possível.

Figura 16 Detalhe de seção transversal de uma asa definindo um aerofólio

Um aerofólio apresenta uma borda de ataque e uma borda de fuga. Denomina-se corda ( c ) a linha que une a borda de ataque com borda de fuga. A linha curva que é sempre simétrica às superfícies superior e inferior denomina-se linha de camber ou linha média. Um perfil aerodinâmico é simétrico quando a linha da corda e a linha de camber são retas coincidentes. O formato de um aerofólio apresenta uma curvatura que atinge seu máximo indicada pela espessura máxima.

Figura 17 Nomenclatura básica de um aerofólio



Figura 8 Detalhe das forças de sustentação e arrasto num aerofólio

Um perfil aerodinâmico quando submetido a uma corrente de fluido com velocidade V∞ apresenta uma força resultante ( R ou FR) que é formada por duas componentes. Uma componente denominada força de sustentação (L ou FL) que atua perpendicular à velocidade e uma força de arrasto (D ou FD) que atua paralela à velocidade. O ângulo de ataque ( α ) é o ângulo formado entre a linha da corda e a velocidade de corrente livre. A força de sustentação é apresentada na forma adimensional como:

p

L

AV

LC

2

2

1∞

=ρ

Onde CL é o coeficiente de sustentação L a força de sustentação V∞ a velocidade de corrente livre e Ap a área projetada máxima da asa. Ap=cb onde c é a corda do aerofólio e b a envergadura da asa. Da mesma forma define-se o coeficiente de arrasto como:

p

D

AV

DC

2

2

1∞

=ρ

Onde CD é o coeficiente de arrasto e D a força de arrasto. Num perfil aerodinâmico o arrasto total origina-se pelo arrasto devido à pressão CDf, o arrasto devido ao atrito (superficial) CDf e o arrasto induzido CDi por efeitos de envergadura finita. Geralmente nos aerofólios o arrasto superficial é o mais importante. Isto pode se inverter para relações t/c maiores que 25% onde t é a espessura máxima do perfil e c a corda do mesmo. Aerofólio: Geralmente CDf >> CDp



7.6 Relação entre Coeficiente de Pressão e Sustentação A sustentação depende de vários parâmetros entre eles o formato do aerofólio, o número de Reynolds e o ângulo de ataque do perfil. Num corpo pode ser determinada quando se conhece a distribuição de pressão em torno do corpo. Na forma adimensional a distribuição de pressão é dada por:

2

2

1∞

∞−=

V

ppcp

ρ

cp é denominado coeficiente de pressão que é a diferença entre a pressão estática local e a pressão estática de corrente livre adimensionalizada pela pressão dinâmica da corrente livre. Na figura abaixo mostra-se a curva típica da distribuição de pressão em torno de um aerofólio. A parte inferior do aerofólio apresenta uma pressão maior que na parte superior. Geralmente isto se apresenta trabalhando com o eixo de cp negativo tal como mostrado. O ponto de estagnação ocorre próximo da borda de ataque. Neste local a velocidade V=0. Para escoamento incompressível Cp=0 neste ponto. Quando a corda é unitária a sustentação é relacionada com o coeficiente de pressão:

( )∫

−=1

0 c

xdCCC pspiL

Onde Cpi é o coeficiente de pressão da superfície inferior e Cps representa coeficiente de pressão da superfície superior.

Figura 19 Distribuição do coeficiente de pressão num aerofólio



7,7 Curva de Sustentação versus Ângulo de Ataque. Um aerofólio com um determinado ângulo de ataque originará uma distribuição de pressão tal como mostrada na figura acima. Para graficar o comportamento da sustentação versus o ângulo de ataque de um perfil aerodinâmico devemos previamente avaliar a distribuição de pressão para cada angulo desejado e posteriormente graficar o resultado. Uma curva típica deste resultado pode ser observada na figura abaixo. Como se aprecia existe uma região em que a sustentação aumenta linearmente com o ângulo de ataque até alcançar a sustentação máxima (CLmax). Nesta região o escoamento apresenta-se suave sem separação da camada limite. Após este máximo o gradiente adverso de pressão provoca a separação do escoamento na superfície superior do aerofólio originando-se um esteira turbulenta. Nestas condições o aerofólio entra em estol o que significa que perde sustentação e ocorre aumento do arrasto. O ângulo em que se origina este fenômeno denomina-se ângulo de estol.

Figura 20 Curva típica de sustentação aerodinâmica versus ângulo de ataque

Um aerofólio simétrico apresentará uma curva de CL versus α que passa pela origem. Isto é para α=00 a sustentação CL=0. No caso de perfis assimétricos para α=00 o aerofólio apresenta um sustentação, contudo existirá um ângulo tal que terá sustentação nula tal como mostrado na figura abaixo.

Figura 21 Curva de sustentação para aerofólios simétricos e assimétricos

Um aerofólio é uma seção de asa que, para efeitos de análise de escoamento, considera-se como bidimensional. Trata-se portanto de uma asa de envergadura infinita. Quando se estudam perfis com envergadura finita devem ser considerados os efeitos tridimensionais provocados pelas pontas das asas, as quais reduzem a sustentação e aumentam o arrasto. Num aerofólio de envergadura finita são originados vórtices de fuga devido a que a pressão média na superfície inferior é maior que a pressão média na superfície superior. Esta diferença de pressão se manifesta



perto das pontas na qual o fluido tende a escoar da parte superior para a parte inferior. Como a asa está em movimento para jusante do aerofólio formam-se estes vórtices de fuga tal como mostrados na figura abaixo.

Figura 22 Circulação e feito de vórtices de fuga num perfil de envergadura finita

Os efeitos de envergadura são correlacionados utilizando a definição da razão de aspecto

pA

b

ojetadaÁrea

asadaocomprimentdoQuadradoarAspectodeRazão

2

Pr

)( ==

onde b é a envergadura e Ap a área projetada. Se o comprimento da corda é constante tal como numa asa retangular, esta relação fica simplificada como ar = b/c. As asas compridas são mais eficientes que as asas curtas devido às perdas das pontas são menos significativas. O efeitos de pontas também origina um arrasto induzido o qual deve ser determinado e adicionado ao arrasto por atrito e por pressão do aerofólio.

Figura 23 Definição de envergadura e área planiforme de uma asa

A relação sustentação/arrasto (L/D) é um parâmetro importante que mede a qualidade aerodinâmica de um perfil. Quanto maior esta relação maior será a eficiência do perfil. Seções modernas de baixo arrasto atingem L/D em torno de 400. Um planador de alto desempenho com ar=40 pode ter um L/D=40. Um avião típico (ar≅12) pode ter L/D≅20.



Figura 24 Definição do ângulo de ataque geométrico efetivo e induzido

Figura 25 Efeito da envergadura finita na sustentação aerodinâmica

Nos perfis com envergadura finita as velocidades dirigidas para baixo reduzem o ângulo de ataque efetivo em proporção ao coeficiente de sustentação.

iefec ααα +=

onde αefec é o ângulo efetivo numa asa com envergadura finita, αi é o ângulo de ataque induzido por efeito da velocidade para baixo originada pelos vórtices de fuga. Isto origina uma redução da inclinação da curva da sustentação como observado na figura. Da teoria de fluido incompressível o ângulo induzido é determinado como:

ar

CLi π

α =

A inclinação da curva de sustentação para um aerofólio com envergadura infinita é definida como coeficiente de inclinação:

αd

dCa L

o =

Desta forma a sustentação pode ser avaliada para uma asa de envergadura infinita em função de ao curva utilizando a relação:



eefctoL aC α=

( )ioL aC αα −=

Figura 26 Determinação da sustentação para aerofólios de envergadura finita

7.71. Influência da Velocidade Induzida na Força de Arrasto

Numa asa de envergadura finita os vórtices de fuga (Fig.11.22) originam velocidades para baixo que provocam um aumento do coeficiente de arrasto CD , o qual pode ser avaliado como:

DiDD CCC += ∞

onde CD∞ é o coeficiente de arrasto da seção considerada um perfil de envergadura infinita e CDi é o arrasto induzido que pode ser avaliado pela expressão:

ar

CCC L

iLDi πα

2

==

A Fig. 27 mostra a as curvas típicas de sustentação e arrasto para um perfil aerodinâmico em função do ângulo de ataque. Observa-se na curva de sustentação o comportamento linear de CL até o alcançar ângulo de estol (α≅150). Após este ângulo o aerofólio entra em estol, observando-se um queda brusca de CL e um aumento acentuadado do coeficiente de arrasto.



(a ) Sustentação (b) Arrastro

Figura 27 Curvas típicas de sustentação e arrasto para um aerofólio

Efeito da compressibilidade Para corpos perfilados para escoamentos com número de M<0,5 os efeitos de compressibilidade no coeficiente de arrasto não são significativos. Já para escoamentos com M alto o coeficiente de arrasto é fortemente dependente do número de Mach, como se observa no exemplo da figura.

Figura 28 Efeito da compressibilidade no escoamento de aerofólios



7.7.2 Velocidade mínima de vôo

Nas condições de estado de vôo constante (condições de cruzeiro) a sustentação (FL) deve ser igual ao peso da aeronave (W).

AVLW 2

2

1∞== ρ

Figura 29 Equilíbrio do peso e da sustentação num avião em velocidade de cruzeiro

A velocidade mínima (Vmin) de vôo é obtida quando CL=CLmax.

AC

WV

L maxmin

2

ρ=

Desta forma a velocidade mínima de aterrissagem pode ser reduzida pelo aumento de CLmax. ou pelo aumento da área da asa. Os flapes são partes móveis da borda de fuga de uma asa que podem ser prolongados num aterrissagem e decolagem com a finalidade de aumentar a área efetiva da asa.



ANEXO

EQUACOES BASICAS DE MECANICA DOS FLUIDOS



8. EQUACOES BASICAS DE MECANICA DOS FLUIDOS

Equação da conservação da massa 0)( =m

dt

d

Equação da Quantidade de Movimento (2ª Lei de Newton)

FVmdt

d rr=)(

Equação do Momento da Quantidade de Movimento

( ){ } FrrVmdt

d rrrr×=×

Equação da Conservação da Energia dt

dW

dt

dQE

dt

d −=)(

Onde: m Massa do fluido

Vr

Vetor de velocidade da partícula de fluido

rr

Vetor posição da partícula de fluido

Fr

Vetor das forcas agindo sobre a partícula de fluido

E Energia total Q Calor W Trabalho



FORMAS INTEGRAIS DAS EQUACOES DO MOVIMENTO As equações integrais podem ser descritas a partir de uma equação geral reconhecendo os efeitos externos e termos característicos.

∫ ∫+∀∂∂=

vc sc

ext AdVdt

Err

ξρξρ

Conservação da massa: 0=extE 1=ξ Quantidade de Movimento:

BSext FFErr

+= Vr

=ξ Momento da Quantidade de Movimento: eixoBSext TFrFrE

rrrrr+×+×= Vr

rr ×=ξ

Equação da Energia dt

dW

dt

dQEext −= e=ξ

Onde e representa a energia total por unidade de massa e E/m

++= int2

2

1ugzVe

sendo uint a energia especifica interna (energia por unidade de massa). As forcas que agem em fluidos são basicamente as forcas de superfície e as forcas de campo. As forcas de superfície são formadas pelas forcas por efeito o de tensões normais ou de pressão e das tensões tangenciais ou de cisalhamento.

∫∫ +=+=AASSpS dApdAFFF ττ

rrr

A forca de campo dada por:

∫∫∫ ∀=∀==vcvcvcB dgdBdmBF ρρ rrrr

As forcas de campo e de superfície podem ser representadas pelas suas componentes:

idFidFidFFd SpzSpySpxSpˆˆˆ ++=

r

idFidFidFFd zSySxSSˆˆˆ

ττττ ++=r



FORMA VETORIAL DO CAMPO DE VELOCIDADES O vetor de posição ou de deslocamento de uma partícula de fluido e dado por:

krjrirr zyxˆˆˆ ++=r

A velocidade e uma função vetorial da posição e do tempo com três componentes u,v e w sendo cada componente um campo escalar

ktzyxwjtzyxvitzyxutrV ˆ),,,(ˆ),,,(ˆ),,,(),( ++=r

Outras grandezas podem ser determinadas manipulando matematicamente o campo de velocidades, denominadas propriedades cinemáticas: Propriedades Cinemáticas: Vetor de Deslocamento

dtV∫=rr

r

Aceleração

dt

Vdr

r =a

Vazão em Volume

∫= AdVrr

Q

Vetor rotação – Velocidade Angular

Vrr

×∇=2

1ω

2ª Lei de Newton aplicada a Fluidos.

amFrr

=

Apresenta-se para fluidos em movimentos definindo a aceleração substancial da partícula de fluido.

=

Dt

VDmF

rr

Onde

( )VVt

V

Dt

VD rrrr

∇+∂∂=



Para estudar o movimento dos fluidos devemos conhecer algumas regras básicas assim como operadores específicos. Regra da Cadeia. Seja uma variável de f, Dependente de coordenadas espaciais e do tempo de f(x,y,z,t), Para obter uma derivada temporal escalar da mesma variável pode-se aplicar a regra da cadeia.

dt

dz

z

f

dt

dy

y

f

dt

dx

x

f

t

f

dt

df

∂∂+

∂∂+

∂∂+

∂∂=

Gradiente ou Operador Nabla As variáveis de cinemática dos fluidos podem ser manipuladas escritas de modo mais compacto quando se utiliza o operador denominado Gradiente o Operador Nabla definido como. Gradiente k

zj

yi

xˆˆˆ

∂∂+

∂∂+

∂∂=∇

O produto deste operador com um vetor velocidade resulta no divergente do vetor . Por exemplo, o divergente do vetor velocidade e dado por: Divergente da Velocidade z

w

y

v

x

uV

∂∂+

∂∂+

∂∂=∇

r

Conservação da massa escoamento compressível e incompressível em regime não-permanente: Eq. da conservação da massa.

( ) 0=∇+∂∂

Vt

rρρ

z

w

y

v

x

u

tV

t ∂∂+

∂∂+

∂∂+

∂∂=∇+

∂∂ ρρρρρρ r

No caso em que o escoamento é em regime permanente, com fluido incompressível. Escoamento Incompressível Regime permanente.

0=∇Vr



ACELERAÇÃO DE UMA PARTÍCULA DE FLUIDO A aceleração de uma partícula de fluido e dada por:

kajaia zyxˆˆˆa ++=

r

a qual pode ser determinara em função do vetor velocidade : dt

Vdr

r=a

kdt

dwj

dt

dvi

dt

du

dt

Vd ˆˆˆ ++=r

Utilizando a regra da cadeia para cada componente u,v,w:

dt

dz

z

u

dt

dy

y

u

dt

dx

x

u

t

u

dt

du

∂∂+

∂∂+

∂∂+

∂∂=

Como se trata de uma partícula especifica. dt

dww

dt

dyv

dt

dxu ===

z

uw

y

uv

x

uu

t

u

dt

du

∂∂+

∂∂+

∂∂+

∂∂=

De modo compacto podemos representar esta equação como:

)(uVt

u

dt

du ∇+∂∂=

r

Aplicando o mesmo procedimento para o componente u, v e e w encontramos as seguintes expressões:

z

uw

y

uv

x

uu

t

u

dt

du

∂∂+

∂∂+

∂∂+

∂∂=

z

vw

y

vv

x

vu

t

v

dt

dv

∂∂+

∂∂+

∂∂+

∂∂=

z

ww

y

wv

x

wu

t

w

dt

dw

∂∂+

∂∂+

∂∂+

∂∂=

)(uVt

u

dt

du ∇+∂∂=

r

)(vVt

v

dt

dv ∇+∂∂=

r

)(wVt

w

dt

dw ∇+∂∂=

r



A aceleração total de uma partícula e denominada também aceleração substancial ou material Aceleração total de uma partícula

VVt

V

Dt

VD rrrr

∇+∂∂=

Aceleração total de uma partícula

z

Vw

y

Vv

x

Vu

t

V

Dt

VD

∂∂+

∂∂+

∂∂+

∂∂=

rrrrr

Derivada substancial

zw

yv

xu

tDt

D

∂∂+

∂∂+

∂∂+

∂∂=

Aceleração Local

∂∂

t

Vr

• Trata-se de uma aceleração que ocorre no tempo.

• Ocorre em escoamentos transientes e em regime permanente.

• E nula para escoamento em regime permanente.

Aceleração Convectiva

∂∂+

∂∂+

∂∂

z

Vw

y

Vv

x

Vu

rrr

• Aceleração que se manifesta em escoamentos com mudanças de geometria. • Escoamentos em regime permanente podem ter grandes acelerações

convectivas devido a mudanças de geométrica.



ROTACIONAL O rotacional e o produto do operador Nabla por uma função vetorial. O rotacional da velocidade e dado por:

ky

u

x

vj

x

w

z

ui

z

v

y

w

wvu

zyx

kji

V ˆˆˆ///

ˆˆˆ

∂∂−

∂∂+

∂∂−

∂∂+

∂∂−

∂∂=∂∂∂∂∂∂=×∇

r

Desta forma, o vetor da velocidade angular (vetor rotação) local como:

Vrr

×∇=2

1ω

kji zyxˆˆˆ ωωωω ++=

r

ky

u

x

vj

x

w

z

ui

z

v

y

w ˆ2

1ˆ2

1ˆ2

1

∂∂−

∂∂+

∂∂−

∂∂+

∂∂−

∂∂=ωv

Vorticidade Defini-se a vorticidade como duas vezes o valor da rotação

Vrrr

×∇== ωζ 2 A vorticidade e o rotacional estão associados com escoamentos viscosos os quais apresentam tensões de cisalhamento. Escoamento Irrotacional 0=ωr Um escoamento não viscoso não apresenta tensões de cisalhamento, portanto e denominado irrotacional. Desta forma 0=ωr . Significa que suas componentes também devem ser nulas.

0

0

0

=

∂∂−

∂∂

=

∂∂−

∂∂

=

∂∂−

∂∂

y

u

x

v

x

w

z

u

z

v

y

w

Escoamento Irrotacional 0=ωr

0=×∇ Vr

Escoamento Rotacional 0≠ωr

0≠×∇ Vr





ANEXO

GRAFICOS



APOSTILA_ESCOAMENTOS_VISCOSOS_2010

Documents

Transcript of APOSTILA_ESCOAMENTOS_VISCOSOS_2010