Notas Hidr Exp 2019 S2

Hidráulica Experimental
Notas de Aula - Versão 1.8 - 2019/s2

Prof. Milton Dall’Aglio Sobrinho
0 A DICIPLINA HIDRÁULICA EXPERIMENTAL .................................................. 01

0.1 O que são e o que fazem os Fenômenos de Transporte .................... 01
0.2 Relação com as disciplinas do curso .................................................. 03
0.3 Algumas definições iniciais................................................................. 04
1 QUANTIFICAÇÃO DOS ESCOMENTOS ......................................................... 08

1.1 Vazão e Fluxos ................................................................................... 08
1.2 Relação Básica ente Velocidade e Vazão .......................................... 09
1.3 Fluxo de Grandezas Extensivas Transportadas ................................. 11
1.4 Vazão em Seções com Velocidade Variável ...................................... 15
1.5 Exemplos Numéricos ......................................................................... 19
1.6 Exercícios Sugeridos ......................................................................... 21
1.7 Relação Geral Entre Velocidade e Fluxos .......................................... 24
1.8 Exercícios Sugeridos ......................................................................... 30
2 DESCRIÇÃO DOS ESCOAMENTOS ............................................................... 32

2.1 Trajetória de uma Partícula Fluida ................................................. 32
2.2 Velocidade e Aceleração de uma Partícula Fluida ....................... 33
2.3 Linha de Corrente - Um Novo Ponto De Vista .................................... 37
2.4 Velocidade e Aceleração em um Ponto (Análise Euleriana) ............... 39
2.5 Linha de Emissão e Linha de Tempo ............................................... 47
2.6 Perfis de Velocidade ........................................................................ 50
2.7 Sistema e Volume de Controle............................................................ 51
2.8 Classificação dos Escoamentos ....................................................... 53
3 CONSERVAÇÃO DE GRANDEZAS (Equação da Continuidade) ..................... 57

3.1 Conservação da Massa .................................................................... 57
3.2 Misturas Homogêneas - Balanço de Grandeza Extensiva N ............ 62
3.3 Equação Integral do Balanço de Massa ........................................... 68
3.4 Discussão Sobre a Taxa de Variação da Grandeza no V.C. ............... 73
3.5 Exercícios ........................................................................................... 77
4 TRANSFORMAÇÕES DE ENERGIA NOS ESCOAMENTOS ........................... 79

4.1. Equação de Bernoulli ......................................................................... 79
4.2. Conservação da Energia nos Escoamentos ...................................... 81
4.3. Energias e Cargas na Equação de Bernoulli ............................... 85
4.4. Aplicação a Medições de Vazão e Velocidade ................................... 86
4.4.1. Medidor Venturi ....................................................................... 87
4.4.2. Tubo de Pitot ........................................................................... 87
4.4.3. Orifícios de pequenas dimensões ........................................... 88
4.4.4. Bocais em condutos forçados .................................................. 90
4.5. Exercícios ........................................................................................... 91
5 TRANSFORMAÇÃO DE REYNOLDS (Relação Sistema x Volume de Controle) ... 95
5.1 Introdução .......................................................................................... 95
5.2 Do Sistema ao Volume de Controle ............................................... 95
5.3 Balanço Global de Grandezas Extensivas ......................................... 98
6 BALANÇO GLOBAL DE ENERGIA ................................................................... 99

6.1 Aplicação a um V.C. em Regime Permanente ........................ 102
6.2 Problemas isotérmicos: bombas, turbinas hidráulicas e tubulações . 103
6.3 Exemplos Ilustrativos ...................................................................... 105
6.4 Efeito do Atrito nos Escoamentos .................................................... 109
6.4.1. Perdas de Carga em Escoamento em Tubos ........................ 111
6.4.2. Sobre o Fator de Atrito ........................................................... 115
6.5 Efeitos das Bombas e Turbinas sobre as Cargas ............................ 121
6.6 Resumo das Transformações de Energia ........................................ 123
6.7 Exercícios Propostos ....................................................................... 128
7 BALANÇO GLOBAL DE QUANTIDADE DE MOVIMENTO ............................ 130

7.1 Aplicações Elementares: Pás Defletoras ......................................... 131
7.2 Aplicações Elementares: Curvas em Tubulações ............................ 133
7.3 Aplicações Elementares: Perdas em Expansão Brusca .................. 136
7.4 Aplicações Elementares: Estruturas em Canais Abertos ................. 137
7.5 Exercícios Propostos ....................................................................... 140
A N E X O : RESPOSTA DOS EXERCÍCIOS PROPOSTOS ................................. 144
Diferenças em relação à versão 1.7 (de 2016-s1)

– Inclusão do capítulo inicial, sobre a disciplina, sua relação com Fenômenos de Transporte e
a grade curricular e definições básicas sobre mecanismos de fluxos.

– Texto sobre comportamento dinâmico de escoamentos no cap 2
– Exemplo numérico no item 3.2; Exercício proposto no item 3.5

. – Pequenas correções na numeração de figuras e equações.
– Acréscimo de texto sobre hidrometria no cap 1
. – Acrescentados texto e exemplos numéricos no cap. 2.

. – Acrescentados os itens 6.4.1 e 6.4.2 num total de 9 páginas.
– Acrescentado o Anexo com as respostas no corpo do texto, 5 páginas.
CAPÍTULO 0
A DICIPLINA HIDRÁULICA EXPERIMENTAL
No currículo da Engenharia Civil da FE-IS/Unesp aparece a disciplina Hidráulica

Experimental, como a primeira da Área de Hidráulica e Saneamento. Uma forma simples de
definir a disciplina é dizer que ela constitui um subconjunto da matéria Fenômenos de
Transporte. As razões para dividir a matéria Fenômenos de Transporte em duas disciplinas
em nosso currículo são, basicamente, de ordem didática, para facilitar o aprendizado.
Adicionalmente, a divisão permite abordar aspectos de Recursos Hídricos importantes para
o Engenheiro Civil. As partes selecionadas da matéria são as que podem ser empregadas
em recursos hídricos para descrever a quantidade e qualidade dos corpos de água. Desta
forma a disciplina poderia ser definida também como uma introdução à Hidráulica Ambiental,
pelo enfoque que prestigia os escoamentos de água e o transporte de substâncias pelos
escoamentos.
É necessário falar um pouco sobre a matéria Fenômenos de Transporte, na forma
como é adotada nos currículos de engenharia, para poder situar a disciplina Hidráulica
Experimental no contexto mais geral da área de Recursos Hídricos e Qualidade Ambiental,
assim como no contexto da nossa grade curricular.
0.1 O que são e o que fazem os Fenômenos de Transporte
É fundamental para o engenheiro conhecer as leis que regem os fenômenos físicos,

tanto para construir e otimizar o rendimento das máquinas e equipamentos, quanto para
avaliar (quantificar) como suas obras interagem com o meio.
A interação das obras com o meio pode ser radical, quando a atividade do engenheiro
envolve a própria modificação do meio para atingir seus objetivos. Além disso, essas
intervenções podem ser de grande escala, como uma usina hidroelétrica ou uma rodovia.
Em muitas situações a interação com o meio é mais sutil, como no caso das trocas de calor
entre o ar e uma habitação, que podem tornar uma moradia agradável ou insuportável,
dependendo do conhecimento e talento empregados no projeto.
Uma interação importante entre as obras e o meio ambiente ocorre quando os fluidos
se movimentam, transportando energia e massa. Estes fenômenos constituem o campo de
estudo da matéria Fenômenos de Transporte. Em que medida eles são comuns, e qual a
sua influência sobre nós são perguntas que ocorrem naturalmente a quem se inicia no
estudo de FT. Nossa vivência proporciona algumas respostas.
A própria terra, e a vida que ela sustenta, depende da existência da atmosfera e da
hidrosfera transportando e dispersando elementos, tanto os nocivos quanto os necessários à
vida. A circulação dos ventos na atmosfera é um mecanismo que transfere, em forma de
Notas de Hidráulica Experimental – versão 1.8 – 2019/s2 2
calor, a energia recebida do sol. Neste movimento o ar leva grandes quantidades de vapor
d'água (massa) e induz correntes marítimas (quantidade de movimento), contribuindo para a
manutenção das condições climáticas.
A circulação da água é descrita pelo ciclo hidrológico. A energia solar evapora a água,
dando origem às nuvens, que entram na circulação geral da atmosfera, para novamente
precipitarem-se, alimentando rios, vegetação e lençóis subterrâneos. Neste processo, a
água carrega substâncias dissolvidas, provoca erosão e transporta solo em suspensão, e é
utilizada pelo homem como matéria prima e meio de dispersão dos rejeitos de toda espécie
decorrentes da atividade humana.
Passando da escala global ao funcionamento de máquinas e sistemas de engenharia,
notamos que sempre há o transporte de alguma grandeza (massa, calor ou quantidade de
movimento), em que os fluidos exercem papel relevante.
Máquinas, e mesmo organismos vivos, dependem do ar para dissipar o calor
excedente (transporte de calor). Um ventilador aumenta a velocidade desta transferência,
tornando mais suportável o calor de Ilha Solteira. Uma roupa que seca no varal ilustra o
transporte de vapor d´água (massa) pelo ar, fenômeno que também ocorre num lago ou num
campo com vegetação. Aqui também a velocidade do ar tem importância na intensidade do
processo.
Um importante fenômeno que envolve a transferência de massa é a oxigenação dos
rios, responsável pelo processo de autodepuração. Na oxigenação ocorre a incorporação do
ar atmosférico pela água. Este processo é maior nos rios de montanha, com águas agitadas
e ocorre mais lentamente nos rios calmos, com escoamento lento.
Os fluidos em movimento provocam forças nos corpos imersos. Percebemos isso
facilmente quando andamos de bicicleta ou caminhamos em dias de vento forte, ou mesmo
quando mexemos a colher no café para misturar o açúcar. Essas forças resultam da troca de
quantidade de movimento entre os fluidos e corpos sólidos. Esses mesmos fenômenos
ocorrem no interior do fluido, sempre que houver diferenças de velocidade.
É comum ocorrer transporte de mais de uma grandeza. Entre um vento frio e a água
de um reservatório há troca de calor, de massa e de quantidade de movimento, esta última
evidenciada pela formação de ondas. Um bolo quente esfriando troca calor e umidade com o
ar ambiente. O ar aquecido sobe, iniciando um movimento do ar que aumenta a velocidade
da evaporação e do resfriamento do bolo. Um carro quente numa estrada transfere
movimento e calor para o ar circundante.
Finalmente, devemos citar os fenômenos em que há transporte de uma grandeza sem
depender de movimento do fluido na direção do transporte. Exemplo típico é o calor
transmitido no interior de corpos sólidos, de grande importância para a engenharia, por
exemplo, no desempenho térmico das edificações.
0.2 Relação com as disciplinas do curso
As disciplinas de Hidráulica Experimental e Fenômenos de Transporte são as

primeiras da Área de Hidráulica e Saneamento. Como o movimento dos fluidos é muito
importante nos processos de transporte, grande parte da Hidráulica Experimental é
dedicada à medição dos escoamentos e quantificação das grandezas que a água transporta
ao escoar. Essa matéria compõe a Parte 1 da apostila de Fenômenos de Transporte. Com a
base fornecida em Hidráulica Experimental é mais fácil enfrentara segunda disciplina da
grade, Fenômenos de Transporte, com estudo mais detalhado dos escoamentos,
abrangendo as trocas de energia e as forças decorrentes e causadoras do movimento, e a
transferência de massa e calor no ambiente e nos sistemas de engenharia.
As atividades rotineiras do engenheiro hidráulico dependem muito da quantificação

dos escoamentos vista em Hidráulica Experimental. O dimensionamento de tubulações,
canais e sistemas de bombeamento, entre outras, são aplicações estudadas em Hidráulica I
e Hidráulica II, a partir dos conhecimentos teóricos e da experiência adquirida em Hidráulica
Experimental e em Fenômenos de Transporte.
Outra classe de problemas de engenharia envolve a água nos sistemas naturais. Por
exemplo, qual seria o nível máximo das águas de um rio? Essa pergunta é importante para
definir a altura de uma ponte, ou para saber quanta água deve suportar um bueiro
rodoviário. Os projetos das estradas contam com inumeráveis obras hidráulicas deste tipo.
Problemas similares são enfrentados quando se deseja saber qual o volume de água que
deve ser armazenado numa represa para garantir o abastecimento à população. Estes
exemplos, entre tantos outros de quantificação da água disponível ou esperada em um
sistema natural, são tratados em Hidrologia. Na quantificação das diversas fases do ciclo
hidrológico e de sua influência sobre a natureza através de fenômenos como a erosão e
carreamento de solo, a hidrologia utiliza intensivamente conceitos de FT. Os mecanismos
envolvidos no equacionamento da evaporação, por exemplo, e na formação de precipitações
são inteiramente baseados em difusão e advecção de massa, energia e quantidade de
movimento.
A área de recursos hídricos preocupa-se também com a qualidade da água, fator que
a cada dia torna-se mais crítico. Além de levar a água até a população, é necessário garantir
a qualidade que deve ser atingida no tratamento das águas, com custos diferentes para
cada manancial. Depois disso vem a preocupação com a disposição dos resíduos líquidos, e
da sua influência sobre a qualidade disponível para demais aproveitamentos a jusante. O
tratamento de água e de esgotos utiliza intensivamente técnicas desenvolvidas com base
em Fenômenos de Transporte. O estudo da recuperação da qualidade dos rios é
inteiramente baseado em transferência de massa.
0.3. Algumas definições iniciais

Grandeza
É qualquer coisa que pode ser medida fisicamente. Por exemplo, temperatura,
velocidade, massa, energia.
Grandeza Intensiva: o valor da medida não depende da quantidade de massa

considerada. Exemplos: temperatura, velocidade, pressão, massa específica,
concentração.
Grandeza Extensiva: o valor medido depende da quantidade de massa

considerada. Exemplos: quantidade de calor, quantidade de movimento, energia
cinética, volume, massa.
Mecanismos básicos de transporte

Existem dois tipos básicos de transporte de grandezas que dependem do meio para
ocorrer, sendo que um deles depende do movimento do fluido e o outro não.
Advecção - Mecanismo de transporte de massa, calor ou quantidade de movimento,

em que a grandeza é transportada por meio do movimento de fluidos
(Mecanismos Macroscópicos).
Difusão - Mecanismo de transporte de massa, calor ou quantidade de movimento que
ocorre a nível molecular (movimento Browniano) e, portanto, não
depende da movimentação de fluido na direção do transporte
(Mecanismos Microscópicos). Em alguns contextos é usado o termo
“Difusão Molecular” para realçar o mecanismo de transporte.
Radiação - Mecanismo de transporte de calor por ondas eletromagnéticas. É o único
que não depende de meio físico.
Condução - Transporte de calor por difusão é chamado também de condução,
principalmente em meios sólidos.
Convecção - Transporte de calor por advecção é chamado também de convecção.
A Figura 0.1 ilustra o transporte de massa pelo processo de advecção. Imagine que
um poluente seja lançado num escoamento em uma região em que a velocidade é uniforme.
Após o lançamento, no tempo inicial t0, a massa de poluente lançada é transportada com a
mesma velocidade pelo escoamento, de forma que a mancha não se espalha.
V Mancha de poluente no Rio
t0 t1 t2 t3 t4 t5
x
Figura 0.1: Advecção de um lançamento de massa num escoamento com velocidade uniforme.
A posição do centro de massa da mancha de poluentes lançada no instante t0 no

ponto x0 é determinada pela equação cinemática:
𝑥𝐶𝐺 = 𝑥0 + 𝑉 (𝑡 − 𝑡0 )
O exemplo de advecção ilustra ainda a importância do movimento do fluido no

processo de transferência. Por isso, parte do assunto tratado na disciplina Fenômenos de
Transporte integra o corpo de conhecimentos da matéria “Mecânica dos Fluidos”.
A Figura 0.2 ilustra o transporte de massa pelo processo de difusão. Imagine que um
poluente seja lançado num corpo d’água com velocidade nula. Após o lançamento, no tempo
inicial t0, a massa de poluente permanece no mesmo local, mas o tamanho da mancha
aumenta. Esse espalhamento da mancha é resultado da difusão, pois não há movimento
macroscópico da água. O transporte ocorre devido ao movimento aleatório das moléculas,
em escala microscópica, e resulta num transporte da grandeza lançada por difusão. Após
um grande tempo decorrido a tendência é a massa distribuir-se uniformemente por todo o
corpo d’água.
Mancha de poluente em água parada: V = 0
t0 t1 t2 t3
Figura 0.2: Espalhamento de um poluente lançado em água parada, por difusão molecular.
A advecção é muito mais eficiente que a difusão, que depende da movimentação

molecular. Entretanto, a difusão é importante porque existem problemas puramente difusivos
na engenharia, e também porque os processos convectivos podem ser entendidos como a

soma de uma difusão com um transporte mecânico.
Mecanismos combinados de transporte

É muito comum que ocorra transporte pela combinação de vários mecanismos. A
Figura 0.3 ilustra o caso de um lançamento de poluente em que ocorre a advecção
combinada com a difusão.
V Mancha de poluente no Rio
t0 t1 t2 t3
x
Figura 0.3: Espalhamento de um poluente pela combinação de advecção e difusão.
Na Figura 0.3 percebe-se que com o passar do tempo a difusão provoca o

espalhamento da mancha, enquanto que a advecção transporta o centro de massa da
mancha com a velocidade do escoamento.
A Figura 0.4 apresenta uma ilustração da transferência de calor num dissipador com
aletas. As aletas são as superfícies que se projetam da base do dissipador, e são
empregadas para aumentar a área de contato do metal com o ar, facilitando a transferência
do calor. O calor é gerado no chip e passa por condução (difusão) para a base do dissipador.
No interior da base o fluxo de calor é redistribuído, passando para as aletas, por difusão. Da
superfície das aletas o calor é transferido para o ar, por convecção (advecção). O ar
aquecido inicia um movimento de advecção ascendente, transportando para longe o calor
transferido pelas aletas.
Advecção do Ar Aquecido
Convecção
Advecção Advecção
do Ar do Ar
Aquecido Aquecido
Condução
Condução
Fluxo Gerado no Chip
Figura 0.4: Resfriamento de Chip com dissipador: difusão de calor no metal do dissipador e convecção
das aletas para o ar que inicia movimento ascendente advecção de massa e calor.
Difusão Turbulenta
Na maioria dos escoamentos de água e ar envolvidos na transferência de grandezas
é difícil ocorrer a advecção em conjunto apenas com a difusão molecular, porque sendo o ar
e água fluidos com baixa viscosidade, os escoamentos normalmente são turbulentos. As
características da turbulência serão estudadas posteriormente, mas o efeito visível da
turbulência é provocar o espalhamento das manchas de poluente, da mesma forma que a
difusão molecular, mas com uma velocidade muito maior, numa escala macroscópica.
Devido ao efeito provocado ser semelhante ao da Difusão Molecular, essa ação de
espalhamento das manchas pelo escoamento turbulento é frequentemente chamada de
Difusão Turbulenta.
CAPÍTULO 1: QUANTIFICAÇÃO DOS ESCOAMENTOS
1.1 Vazão ou Fluxo de Volume
É muito importante conhecer o volume de fluido que um escoamento transporta. Como os

escoamentos são contínuos é conveniente expressar o volume transportado por unidade de
tempo, ou seja, pelo Fluxo de Volume, FVol, também conhecido como Vazão:
Volume Transportado
FVol  VAZÃO  ( m3 / s )
Tempo decorrido 1.1
A vazão de água transportada por um rio é fundamental em muitos problemas

práticos. Por exemplo, para sabermos se é possível utilizar a água para abastecimento de
uma cidade, ou se o rio comporta o lançamento de esgotos com um determinado nível de
tratamento.
Para medir uma vazão podemos imaginar o experimento representado pela Figura
1.1, conhecido como “método volumétrico direto”. Conhecemos o volume inicial de água no
reservatório e, no instante t = 0, colocamos o recipiente sob o jato de água, parando o
cronômetro ao final de um tempo t qualquer, quando lemos o volume final. A diferença de
volumes fornece o volume escoado durante o intervalo de tempo considerado.
Figura 1.1: Medição de volume transportado pelo escoamento num intervalo de tempo.
Aplicando a definição da equação 1.1 com o volume Vol e com o intervalo de tempo
decorrido t, obtemos o valor da vazão média no período de tempo da medição:
Vol
Q  ( Valor Médio no intervalo t) 1.2
t
Para que a definição seja válida no caso de escoamento variável no tempo, interessa
o valor instantâneo.
Vol
Q  lim
t  0 t
 Q
dVol
dt
(Valor Instantâneo) 1.3
A dimensão do fluxo de volume é [ M3 / T ], e as unidades mais comuns são m 3/s,

m3 /h, L / h, m3 /dia.
Uma vazão só tem sentido quando associada a uma determinada seção. No caso da
Figura 1.1, trata-se da seção de saída do tubo, com área S. Um sinônimo de fluxo é Taxa de
Passagem. Então podemos dizer também que a vazão é a taxa de passagem de volume
através de uma dada superfície
Vazão é um Fluxo de Volume, ou seja, a quantidade de volume
por unidade de tempo que atravessa uma determinada área .
_____ Fluxo de Massa

Em muitas ocasiões é importante conhecer a taxa de transferência de massa através
de uma seção de escoamento. Isso é particularmente verdadeiro no caso de escoamentos
compressíveis.
Dada uma seção qualquer de um escoamento, a quantidade de massa que atravessa
a seção por unidade de tempo é o Fluxo de Massa.
F 
m

 Vol  FM   Q 1.4
M
t t
A dimensão do Fluxo de Massa é [ M / T ], e as unidades são:( Kg/h ), ( ton/h ),
( Kg/s ), ( utm/s ) etc.
1.2 Relação Básica entre Velocidade e Vazão
Nossa experiência cotidiana, por exemplo, com torneiras e mangueiras de jardim,

indica que a vazão é função da velocidade do escoamento. A velocidade do fluido é um dos
fatores principais para definir a capacidade de transporte de grandezas dos escoamentos. A
outra é a área da seção transversal, conforme veremos neste item.
Imagine o escoamento num duto retangular de seção transversal A, transportando
água, com velocidade V uniforme e constante no tempo, conforme esquema da Figura 1.2.
x
 Vol
t=0 t = t A
x
Figura 1.2: Escoamento uniforme num duto retangular – volume que atravessa a seção.
O perfil uniforme significa que qualquer partícula tem a mesma velocidade. Além
disso, o movimento é unidirecional, ou seja, ocorre apenas na direção x. Podemos marcar
uma partícula qualquer com corante, e determinar sua velocidade por meio do deslocamento
registrado num intervalo de tempo t dado:

x
V ; valor médio da velocidade no intervalo t. 1.5
t
Com a velocidade conhecida, é fácil determinar quais partículas serão capazes de
atravessar a sessão “A” num intervalo t. Basta ver que o deslocamento possível nesse
tempo é x = V t . Concluímos que um volume igual ao hachurado irá atravessar a seção
de área A no intervalo t. Então:
Vol Ax
Q   Q V A 1.6
t t
A equação 1.6, embora simplificada, é importantíssima. É empregada na grande
maioria dos cálculos de tubulação, com V igual à velocidade média no tubo. As
simplificações adotadas foram:
– o módulo da velocidade é o mesmo em toda a seção A,
– a direção da velocidade é a mesma em toda a seção A,
– a direção da velocidade é perpendicular à seção A.
A primeira hipótese é equivalente a afirmar que V é a velocidade média na seção. Já
a segunda hipótese é praticamente impossível de ser satisfeita num escoamento real devido
a presença dos contornos sólidos, como as paredes do tubo ou o fundo dos canais.
Posteriormente adaptaremos a equação 1.5 para uso num caso geral.
Exemplo 1.1 : Uma tubulação com 50mm de diâmetro interno abastece um caminhão
tanque de 15.000 l de capacidade com gasolina (  = 860 kg/m3). Sabendo que a velocidade
média no tubo é de 2,0 m/s, pede-se: a) Qual a massa de gasolina transportada; b) Vazão
que sai do tubo; c) Fluxo de massa que entra no tanque; d) Qual o tempo de enchimento
completo do tanque?
SOLUÇÃO:
a) Pode-se usar o valor médio porque a massa é uniformemente distribuída. A partir da
definição de massa específica e sabendo que 1m3 equivale a 1000 litros, obtemos:
kg 1m 3
m   Vol  860 ( 3 ) 15000 (l )( )  12900kg
m 1000l
b) Sabendo que a vazão é a velocidade multiplicada pela área do escoamento, temos:
 d2  0,0502
A   0,00196 m 2
4 4
m m3
Q  VA  2 ( )  0,00196m 2  0,00393  4 litros / s
s s
c) Aplicando a definição do fluxo de massa,
kg m3 kg
FM   Q  869 ( 3
)  0,00393 ( )  3.377
m s s
d) O tempo de enchimento vem da aplicação da definição de vazão:
Vol Vol 15m 3

Q  t    t  3.820s
t Q 0,00393m 3 / s
Exemplo 1.2 : Um fluido com massa específica constante escoa pela redução de diâmetro
de 100mm para 75mm representada na figura. Sabendo que a velocidade no tubo maior é
1m/s, calcule a velocidade no tubo de menor diâmetro.
Volume Constante
V1 Vol1 Vol2 V2
 x1  x2
SOLUÇÃO :
Como o volume de fluido no interior da redução (tracejado na figura) é constante,
deduzimos que o volume trazido pelo tubo de 100mm em cada intervalo de tempo deve ser
igual ao volume que sai pelo tubo menor no mesmo intervalo (Vol 1 = Vol2 ). Mas, pela
definição de vazão é possível calcular os volumes, já que o intervalo de tempo considerado
é o mesmo:
Vol1 Vol 2
Q   ; mas Vol1  x1 A1 e Vol 2  x 2 A2
t t
x1 A1 x2 A2 A1
  V1 A1  V2 A2  V2  V1
t t A2
substituindo as áreas,
A1  d12 / 4 d12 0,12 m
    1,78  V2  1,78
A2  d 22 / 4 d 22 0,0752 s
1.3 Fluxo de Grandezas Extensivas Transportadas

Ao considerarmos um fluido escoando através de uma seção qualquer, podemos
quantificar não só os fluxos de volume e massa do fluido, mas também a quantidade das
grandezas extensivas que o fluido carrega em seu meio.
O Fluxo de uma Grandeza Extensiva, FN, é a quantidade da grandeza N transferida

por unidade de tempo através de uma determinada área.
N
FN 
t 1.7
Um fluxo só tem significado quando associado à área em que ocorre. Essa definição é
válida tanto para fluxos advectivos quanto para difusivos.
As unidades do Fluxo de uma grandeza são as unidades da grandeza divididas pelo

tempo. Por exemplo, um fluxo de massa será expresso em kg/s e um fluxo de calor em
Joules/segundo, ou Watts.
Concentração de Grandezas Extensivas

A concentração de uma grandeza extensiva N qualquer, CN, é a quantidade da
grandeza por unidade de volume do meio considerado. O meio pode ser sólido, líquido ou
gasoso. A concentração média no volume considerado é dada por:
𝛥𝑁
𝐶𝑁 = (valor médio em ∆Vol ) 1.8
𝛥𝑉𝑜𝑙
As unidades da concentração são as unidades da grandeza por metro cúbico. Por

exemplo, para N = quantidade de calor, temos CN dado em J/m3.
O uso mais comum da concentração é para o caso da grandeza N ser a massa de
uma substância dissolvida. Nesse caso tem-se CN dado em kg/m3.
Quando o volume considerado é cada vez menor, a concentração média tende a
um valor pontual (ver item sobre meio contínuo). A concentração num ponto é dada por:
𝛥𝑁
𝐶𝑁 =
𝛥𝑉𝑜𝑙 0 𝛥𝑉𝑜𝑙
𝑁
𝐶𝑁 = ( em um ponto ) 1.9
𝑉𝑜𝑙
Da definição podemos concluir que a massa específica de um meio é uma

concentração da grandeza extensiva N = massa do meio. Por seu uso frequente a massa
específica possui uma notação especial:
= ( ) 𝐶𝑁 (com N=massa do meio)
𝑉𝑜𝑙
Quantidade Específica
A concentração pode ser expressa em termos de quantidade da grandeza por
unidade de massa do meio considerado. Chamada de concentração em massa, ou razão de
mistura média, ou ainda, como adotamos neste curso, a Quantidade Específica de uma
grandeza extensiva N qualquer, ,é dada por:
𝛥𝑁
= (valor médio em ∆m ) 1.10
𝛥
O valor em um ponto, como já vimos, é dado pelo limite para um volume tendendo
a zero:
𝑁
= ( em um ponto ) 1.11
A quantidade específica da grandeza N, para o caso da massa de uma substância A

distribuída no meio (N = mA), dá origem à propriedade chamada de “razão de massa”,
apresentada em (kg/kg). A diferença é que a quantidade específica pode ser definida para
qualquer grandeza extensiva e não somente para a massa da substância dissolvida.
Exemplo 1.3
Com base no conceito de grandeza extensiva, deduzir a equação das concentrações
de Calor, Energia Cinética, e Quantidade de Movimento em um escoamento.
Solução:
a) Concentração de Calor: A quantidade de calor depende da temperatura T e do calor
específico c. Aplicando as definições, temos:
a) Concentração de Calor: A quantidade de calor depende da temperatura T e do
calor específico c. Aplicando as definições, temos:
(𝑇 𝑇 ) 𝑉𝑜𝑙 (𝑇 𝑇 )
𝐶 = =  𝐶 = ( − 0)
𝑉𝑜𝑙 𝑉𝑜𝑙
(𝑇 𝑇 )
 =   = ( − 0)
onde T0 é uma temperatura de referência.
b) Concentração de Energia Cinética:

𝐸𝑐 𝑉2 𝑉2
𝐶𝐸𝑛.𝐶𝑖𝑛é𝑡𝑖 𝑎 = = 𝑉𝑜𝑙 =
𝑉𝑜𝑙 𝑉𝑜𝑙 2 2
𝐸𝑐 𝑉2 𝑉2
𝐸𝑛.𝐶𝑖𝑛é𝑡𝑖 𝑎 = = =
2 2
c) Concentração de Quantidade de Movimento: A quantidade de movimento é uma
⃗ =
grandeza vetorial dada por 𝑃 ⃗ . Assim, sua concentração em volume e massa num
𝑉
escoamento fica:
𝐶𝑄𝑢𝑎𝑛𝑡.𝑀𝑜𝑣 = ⃗
𝑉  𝑄𝑢𝑎𝑛𝑡.𝑀𝑜𝑣 = ⃗V
Fluxo da Grandeza Específica

O Fluxo de uma grandeza extensiva N qualquer pode ser dado em relação à
concentração da grandeza, ou em relação à quantidade específica.
Fluxo de N em função da sua concentração CN :
N N Vol
FN    CN Q  FN  CN VA 1.12
t Vol  t
Fluxo em função da quantidade específica :
FN 
N m
m  t
  FM  FN  VA 1.13
Ilustração: Uma dedução alternativa das equações do fluxo ocorre ao considerar a analogia
entre o escoamento e um trem em movimento, conforme a Figura 1.3. Os vagões equivalem
ao fluido em escoamento e os passageiros nos vagões são análogos às grandezas

extensivas conduzidas pelo escoamento.
Analogia:
Seção S
V
Escoamento = Trem
Fluido = Vagões
Pessoas = Grandeza N
Num.Pessoas Num.Pessoas Num.Vagões

FluxoPessoas  
U .Tempo Vagão U .Tempo
      
definição Quantidade Fluxode
deFluxo Específica Vagões
Figura 1.3: Analogia com trem em movimento para definição do fluxo de grandezas extensivas.
Pensando num vagão como 1 m3 (unidade de volume), ou como 1 kg (un. de massa)

de fluido, obtemos as equações genéricas dos fluxos pela extensão do raciocínio utilizado
para calcular o fluxo de pessoas:
Quant. da Grandeza Quant. de Massa

FN  x
Quant. de Massa Un. de Tempo
  
 FM
Quant. da Grandeza Quant. de Vol.

FN  x
Quant. de Vol. Un. de Tempo
    
C Q
_____ Exemplos de Fluxos de Grandezas Extensivas

Aprendemos nesse item que o fluxo de qualquer grandeza cuja quantidade total no
fluido depende da massa de fluido considerada pode ser descrito em função do fluxo de
massa. Por exemplo, para algumas grandezas extensivas consideradas:

FCALOR  c T (  V A ) FQ.MOV  V ( V A )
V2
FE .CINÉTICA  ( V A ) FENERGIA  e (  V A )
2
Em todos os exemplos vimos que sempre a quantidade específica da grandeza é
multiplicada por uma parte comum que é o Fluxo de Massa Esse termo representa, como já
vimos, a quantidade da grandeza transportada (por Advecção) por unidade de massa do
fluido transportador.
Exemplo 1.4 : Um rio possui vazão de 10m3/s de água com concentração de sólidos totais
de 250mg/L. Calcular: a) o fluxo de massa de sólidos totais e b) a massa de sólidos
transportada pelo rio em um dia.
Solução:
a) uma vez conhecida a concentração da grandeza extensiva (sólidos totais), o fluxo é dado
pela equação 1.12.
mg g L m3 g m3 g
FST  ( 250  0,001 1000 3 )  10  250 3 10  2.500
L mg m s m s s
b) a massa transportada num dia vem da definição do fluxo médio (eq. 1.7):
kg s
mST  FST  t  2,5  86.400  216.000 kg / dia
s dia
1.4 Vazão em Seções com Velocidade Variável
A hipótese de perfil uniforme utilizada na relação básica praticamente nunca ocorre na

prática, devido à influência da viscosidade dos fluidos e ao fenômeno da adesão do fluido
aos contornos sólidos do escoamento. Nos fluidos reais sempre vai existir uma região
próxima aos contornos sólidos com variação pronunciada da velocidade, chamada de
camada limite.
Imagine um trecho de rio retilíneo esquematizada na Figura 1.4, com a seção
transversal com diferentes profundidades. A utilização de traçadores permite concluir que
existe um perfil variável de velocidades, que pode ser aproximado por 3 velocidades
diferentes.
V3
V1 V2 V3
V2
A1 A3
V1 A2
Corte da Seção Transversal
Seção
Figura 1.4: Escoamento com velocidade variável
A vazão pode ser considerada como a soma da contribuição de 3 seções distintas e
independentes, com velocidades e áreas diferentes.
Q  Q1  Q2  Q3  V1 A1 V2 A2  V3 A3
generalizando para um número qualquer de áreas, temos:

n
Q  V
i 1
i Ai 1.14
_____ Discussão sobre Modelo de Medição

O sinal de aproximadamente igual na equação 1.14 surge ao fazermos apenas 3
medições de velocidade para aproximar um perfil real de velocidades que varia
continuamente, conforme ilustrado na Figura 1.5.
V3 A3
V2 A2
perfil
perfil V1
real
aproximado A1
Seção
Figura 1.5: Escoamento Real e Modelo de Medição, com Velocidades Constantes
É claro que o perfil aproximado não representa com perfeição o perfil real de
velocidades, que varia continuamente. Entretanto, ao adotar um perfil composto de apenas 3
velocidades constantes, estamos adotando um modelo de medição que pode ser
suficientemente exato para nossos propósitos.
Sabemos que o perfil real de velocidades não é como descrito pelo modelo
simplificado de medição. Podemos reduzir o erro de modelo fazendo mais medições de
velocidade ao longo da seção transversal, mas o custo das medições adicionais necessárias
pode não ser viável.
O erro de modelo numa medição pode ser aceitável ou não, dependendo de nosso
objetivo. No caso de uma medição de vazão em rios utilizando flutuadores, pode ser
aceitável um modelo bem simplificado, se nosso objetivo for uma estimativa para fins de
anteprojeto.
Pode-se perceber a partir da Figura 1.6, que a divisão da seção em áreas menores e
um maior número de medições de velocidade diminui o erro de modelo.
Seção Seção
A6
A3
A5
A2 A4
A3
perfil perfil
real
A1 real
A2
A1
aproximado com 3 velocidades aproximado com 6 velocidades
Figura 1.6: Aumento do número de medições de velocidade diminui o erro de modelo.
O efeito do aumento do número de n de subáreas consideradas na equação 1.14

pode ser visualizado num gráfico como o da Figura 1.7.
Vazão (m3/s) Q1
Q4
QReal
1 2 3 4 5 6 7 ... j Número de subáreas

Figura 1.7: Aumento do número de medições de velocidade diminui o erro de modelo.
A conclusão que se impõe é que o valor correto surge no limite de uma série de
medições com um número crescente de subdivisões de áreas.
n
QRe al  lim Qn  lim  Vi Ai
n n 
i 1
1.15
O sinal de igualdade na equação 1.15 indica que no limite, para número muito grande
de áreas, deixa de existir o erro de modelo. Matematicamente essa operação é denotada
pelo sinal da integral:
Q   VdA
A 1.16
Pontos Importantes na equação 1.16
1. Velocidade na direção perpendicular a Área

2. “dA” é o elemento diferencial de área: é a maior área em que V pode ser
considerado constante (não é a derivada da função área)
3. O limite da integração A não é operacional, apenas indica que os limites reais
devem cobrir toda a área desejada.
4. Na vazão calculada resta apenas o erro de medição.
_____ Hidrometria – Medição de Vazão em Rios
A determinação de vazões em rios é uma aplicação prática importante das equações

1.4 e 1.16. Uma diferença importante, entre os esquemas apresentados e a utilização na
prática, é a consideração da variação de velocidades com a profundidade, e não só com a
distância da margem do rio.
Um dos métodos de divisão da seção para a medição de vazão de rios é o chamado
método dos dois pontos, em que as velocidades são medidas a 20% e a 80% da
profundidade da subárea considerada. Um exemplo de divisão segundo a técnica de dois
pontos é mostrado na Figura 1.8.
i=1 2 wj n
0,6p
Profundidade pj
Figura 1.8: Divisão da seção segundo a técnica dos dois pontos em cada vertical.
Observa-se na Figura 1.8 que nas seções mais rasas a velocidade foi medida em
apenas um ponto, correspondendo a 60% da profundidade. Adota-se este critério quando a
profundidade for entre 0,15 e 0,60m.
A ANA (Agência Nacional de Águas) considera justificável o método dos dois pontos
quando o fator tempo é importante na medição, por exemplo, em cheias e medições com
grande variação de nível da água.
A ANA considera preferível o método detalhado que, como indica o nome, adota uma
subdivisão mais fina em cada vertical. A Tabela 1.1 indica as recomendações da ANA para
adoção do método detalhado, com o número de medições em cada vertical definido em
função da profundidade.
Tabela 1.1: Tabela de pontos de medição de velocidade no método detalhado. Fonte: ANA
N° de Posição na vertical Velocidade média na Profundidade

pontos (em relação a p ) vertical (vM) (m)
1 0,6 V0,6 0,15 – 0,6
2 0,2 e 0,8p (V0,2+V0,8)/2 0,6 – 1,2
3 0,2; 0,6 e 0,8 (V0,2+ 2V0,6+V0,8)/4 1,2 – 2,0
4 0,2; 0,4; 0,6 e 0,8 (V0,2+2V0,4+2V0,6+V0,8)/6 2,0 – 4,0
6 S*; 0,2; 0,4; 0,6 e F* [VS+2(V0,2+V0,4+V0,6+V0,8)+VF]/10 > 4,0
VS = v na superfície; VF = v no fundo; S = 0,10m; F = fundo, determinado pelo lastro.
Quanto ao número de seções a ANA recomenda 0,05% a 0,025% da largura como

distância entre duas verticais de medição (20 a 40 divisões), e indica 0,30m como a mínima
distância entre seções.
Os medidores tipo ADCP (Acoustic Doppler Current Profiler), ou perfiladores acústicos
de velocidade são cada vez mais utilizados na prática. Conforme seu nome indica, os
ADCPs determinam a velocidade por efeito “Doppler” em feixes de ultrassom, em um grande
número de células em cada vertical. Ao ser conduzido através de uma seção transversal o
equipamento mede o deslocamento e calcula a velocidade média em cada célula. Com as
velocidades da água e o deslocamento do barco o software integra as contribuições de cada
célula para indicar diretamente o valor da vazão ao final da travessia de uma margem a
outra do rio.
A Figura 1.9 apresenta um esquema da divisão da seção utilizada pelos medidores
tipo ADCP, que avaliam a integral da equação 1.16 de forma automática.
Figura 1.9: Divisão da seção em células para integração da vazão pelos medidores tipo ADCPl. Fonte:
ANA (2012), disponível em http://arquivos.ana.gov.br/infohidrologicas
1.5. Exemplos Numéricos

Exemplo 1.5: Deseja-se saber a vazão de um córrego com a seção transversal dada na
Figura 1.10. Tendo em vista os objetivos da medição, julgou-se suficiente a divisão da seção
em duas subáreas, nas quais foram medidas as velocidades seguintes: V 1 = 0,3m/s e V2 =
0,7m/s. Determinar a vazão.
Seção Real
Modelo Adotado
A1 0,35m
A2 0,8m
2,5m 0,9m
Figura 1.10: Seção transversal real e modelo adotado para a medição de velocidade.
Adotando-se o modelo de medição exposto na figura, tem-se:

Q  Q1 + Q2 = V1A1 + V2A2
Q  0,3 (2,5 G 0,35/2) + 0,7 (0,9 G 0,8) = 0,13125 + 0,504 = 0,63525
Resposta: a vazão do rio é aproximadamente 0,6m3/s, ou 600L/s.

Exemplo 1.6: Um rio com seção retangular de 10m de largura com 1,5m de profundidade
possui um perfil de velocidades dado na Figura 1.11. Determinar a vazão.
2,0m/s A z
y 10m
1,5m y
x
Vista Superior Seção Transversal
Figura 1.11: Seção transversal e perfil de velocidades.
Solução:
1. Função da velocidade:
Observa-se que a velocidade varia linearmente com y. Ajustando-se uma reta aos pontos
dados (y = 0, V = 0 e y = 10, V = 2) obtém-se:
V = 0,2y
2. Determinação do elemento diferencial de área:
Analisando a função velocidade percebe-se que V não depende de z. Por isso podemos
adotar o elemento diferencial de área dado no esquema:
z
dy
1,5 y
dA = z dy
Elemento diferencial de área adotado:
3. Solução, usando a equação 1.16
10
y 10 y2
Q   VdA  
A y 0
0,2 y  1,5 dy  0,3
2
 15
0
4. Resposta: A vazão do rio é de 15m3/s.
Exemplo 1.7: Um rio possui seção transversal que pode ser considerada triangular de
acordo com a Figura 1.12. O perfil de velocidades, dado na figura, é o mesmo do exemplo
1.5. Pede-se calcular a vazão.
2,0m/s A
z
10m
y
2m
y
x
Figura 1.12: Seção transversal e perfil de velocidades.
Solução:
1. Função da velocidade:
A velocidade varia linearmente com y como no exemplo anterior:
V = 0,2y
2. Determinação do elemento diferencial de área:
Analisando a função velocidade percebe-se que V não depende de z. Por isso podemos
adotar o elemento diferencial de área dado no esquema, com dimensão finita na vertical:
dy
h dA = h dy
z
y
O elemento de área adotado pode ser escrito apenas em função de y, pois h = 0,2y:
dA = 0,2ydy.
3. Solução, usando a equação 1.16

3 10
y 10 y
Q   A
VdA  
y 0
0,2 y  0,2 y dy  0,04
3
 13,33
0
4. Resposta: A vazão do rio é de 13,3m3/s.
1.6. Exercícios Sugeridos
1.6.1. Um escoamento de água quente a 45°C ( ρ = 995kg/m 3) ocorre com velocidade de

2m/s num tubo com área 0,01m2. A água possui uma concentração de 200mg/L de sólidos
totais. Pede-se: a) calcule a vazão de água; b) calcule o fluxo de massa de água; c) calcule
o fluxo de sólidos totais transportados pela água; d) calcule o fluxo de energia térmica
(quantidade de calor) em relação à temperatura de referência de 0°C transportado pelo
escoamento. Dado: calor específico da água c = 4180J/kg°C.
1.6.2. Um rio recebe a água de um afluente pouco antes de um trecho retilíneo, com a seção
dada na figura. No trecho indicado foram lançados flutuadores e medidas as concentrações
de matéria orgânica em 3 pontos de amostragem ao longo da seção transversal, conforme a
figura com a seção de medição. A escala da seção é dada pelo quadriculado com 0,5m de
lado. Os valores medidos foram: V1 = 0,7m/s, V2 = 1,5m/s, V3 = 2,0m/s; C1 = 200mg/L, C2 =
195mg/L e C3 = 25mg/L.
1 2 3
y
Rio
Afluente Seção de Z
x Seção de Medição
medição
Planta de Situação Y
Pede-se:
a) Adote e justifique um modelo de medição para a seção e para os perfis de
velocidade e de concentração;
b) com o modelo adotado em (a) calcule a vazão de água no canal;
c) com os modelos adotados em (a) calcule o fluxo de massa de matéria orgânica
transportado pelo canal.
1.6.3. O reservatório de acumulação de uma pequena hidrelétrica recebe contribuição de 3

rios, com as vazões e conteúdo de sólidos suspensos médios dados na tabela.
Rio 1 Rio 2 Rio 3 Saída

Vazão (m3/s) 10 20 40 QSAI
CSS (mg/L) 100 2000 1000 100
a) Qual a vazão média de saída?

b) Qual a taxa média de acúmulo de massa de Sólidos Suspensos no reservatório?
c) Sabendo que a massa específica do material sólido depositado no reservatório é de
1600 kg/m3, calcular a perda anual de volume útil do reservatório devido ao acúmulo
de sólidos.
1.6.4. Um meio para determinar a vazão de rios consiste na injeção de substâncias

traçadoras, como sais ou corantes. Numa determinação de vazão em um córrego foram
lançados 2l/s de água com uma concentração de corante fluorescente igual a 5g/l. Numa
seção a jusante, após a completa mistura do traçador, retirou-se uma amostra da água,
obtendo-se uma concentração de corante de 0,2g/l. Qual a vazão do córrego?
1.6.5. A sua equipe executou medidas de velocidade e determinou trajetórias de partículas

num trecho de rio onde se pretende lançar um efluente industrial, obtendo as trajetórias
apresentadas. Observou ainda que as trajetórias de partículas ao longo do tempo
praticamente não variam, podendo-se considerar o escoamento permanente. A velocidade
na seção 1 é praticamente uniforme em toda a seção e igual a 0,4m/s. Com base nessas
informações, responda as seguintes questões:
2
c
a
b
a'
1 3
2,5m 3m
1m
60m 25m 40m
Seção 1 Seção 2 Seção 3
a) Qual a vazão do rio e as velocidades nas seções 2 e 3 ?

b) Se no ponto "a" forem lançados 150l/s de efluente com uma concentração volumétrica de
500mg/l de uma substância poluente inexistente no trecho a montante do rio, qual será a
concentração resultante nos pontos "b", e "c"?
c) Existe a alternativa de lançamento do outro lado do rio, no ponto “a”. Do ponto de vista
dos habitantes da ilha, qual dos pontos é preferível? Justifique sua resposta.
Dica: considere a divisão do escoamento pelas linhas de corrente e mistura completa em
cada seção a jusante.
1.6.6. Dados os perfis de velocidade e de concentração de Cloretos na água do rio da figura,

com seção aproximadamente triangular, pede-se calcular: a) o fluxo de volume em m3/s; b) o
fluxo de massa de cloretos, em gramas por segundo. Dados V Max = 2,0m/s, CMax = 200mg/L
e CMin = 200mg/L.
C Max
A z
VMax C Min
10m
y
2,0m y
x
1.6.7. Um canal retangular de 1 m de profundidade e 3 m de largura, transporta água

salgada (= 1000 kg/m3) com concentração igual a 100 mg/kg, com um perfil de velocidade
dado pela equação abaixo, com V em (m/s) e a cota y em metros, com origem no fundo do
canal:
V  15 
, 1  ( y  1) 2 
.
Pede-se calcular: a) velocidade média; b) fluxo de volume (vazão); c) fluxo de massa de
água no canal; d) fluxo de sal conduzido pelo canal.
1.6.8. Para obter permissão legal para operar, uma indústria comprometeu-se a lançar no
máximo 3 litros por segundo de efluentes com uma concentração máxima de cianetos igual
a 3 miligramas por por litro. Uma associação de defesa ambiental desconfia do cumprimento
da lei pela indústria, mas uma comissão de vistoria formada para investigar o problema não
foi bem recebida pela empresa. Em consequência, você foi consultado para reunir dados
para amparar uma ação legal contra a indústria. Sua equipe fez medições da seção e
velocidade do rio a montante da indústria suspeita, e das concentrações de cianeto acima e
a jusante do ponto de lançamento, obtendo os seguintes dados:
MONTANTE: V = 0,6m/s, A = 6,3m2, CCN = 0,0000mg/l
JUSANTE : CCN = 0,0081mg/l
Determinar se há base legal para processar a indústria.
1.7. Relação Geral Entre Velocidade e Fluxos
No item 1.4 deduzimos o caso de perfil de velocidades variável, perpendicular à seção

considerada. No caso mais geral a seção pode ter forma e inclinação qualquer. Este caso
será trabalhado, em primeiro lugar, transformando a superfície curva em uma superfície
aproximada por várias superfícies planas. Assim, o problema geral de superfície curva se
reduz a uma sucessão de problemas de superfícies planas com inclinação qualquer.
A seguir mostraremos como calcular o fluxo através de uma superfície plana de
inclinação qualquer.
_____ Passo inicial - Definição Vetorial da Área

Uma superfície de inclinação qualquer no espaço pode ser definida pelo seu vetor
área, conforme mostra a Figura 1.9. Um elemento de área dA é definido por um vetor com
módulo dA e direção do versor n, normal à superfície considerada. O sentido do vetor área é
positivo quando se dirige para fora em relação a uma superfície fechada. Se não existir uma
superfície fechada para referência o vetor só possui direção definida.
dA dA
dA dA
dA
Figura 1.9: Definição vetorial da área.

  
dA  dA  n sendo n o versor normal 1.17
_____ Segundo passo – identificar o volume que atravessa a seção

Para isso iniciaremos com uma situação mais simples dada por velocidade constante
na seção, (perfil uniforme) e área com inclinação constante. Imaginemos então um

escoamento com perfil uniforme de velocidades, representado pelas linhas de corrente da
Figura 1.10. A dimensão na direção z é dz.
dx dA

y V dh
ds
x
Seção dA = dsdz
Figura 1.10: Fluxo de volume através de uma seção inclinada em relação à velocidade.
O módulo do vetor área na Figura 1.10 é dado por:

dA  ds dz
O volume é dado pela porção hachurada, que corresponde ao volume de um prisma
cuja seção é um losango de base dx e altura dh.
dVol  dx dh dz
Temos, pela geometria da seção, que
dh  ds cos 
dVol  dx ds dz cos  dVol  dx dA cos 
Lembrando a definição de fluxo e que dx  V dt , vem:
dVol dx
dQ   dA cos  V dA cos
dt dt 1.18
A equação 1.18 indica um produto escalar entre os vetores da velocidade e da área,
de modo que:
 
dQ  V  dA 1.19
_____ Terceiro passo – identificar a Integral de Área
Já vimos no item 1.4 que um perfil qualquer de velocidades pode ser aproximado por
segmentos elementares nos quais a velocidade é constante. A mesma idéia é válida para
dividirmos também uma seção de forma qualquer em vários planos retilíneos.
Assim, qualquer área e qualquer perfil podem ser aproximados, no limite, por uma
sucessão de áreas planas e perfis constantes, sendo que cada uma contribui com uma
vazão elementar, conforme a equação 1.19. Veja a Figura 1.11.
Seção aproximada
V1 dA 1
V
V dA 2
2
V3 dA 3
V4 dA4
Perfil real Seção real Perfil aproximado
Situação real Modelo aproximado
Figura 1.11: Fluxo de volume através de uma seção qualquer.
Podemos repetir o raciocínio utilizado no item 1.4. A vazão total é aproximada por
uma soma que engloba as contribuições de toda a área:
Q  dQ1  dQ 2  dQ 3  dQ 4
A aproximação é exata no limite, quando o número de áreas dA  
n  
Q  V . dA i
i 1
  n  
Q  lim  V .dAi  Q   V . dA 1.20
A
n i 1
O símbolo "A" na integral significa que o somatório das contribuições deve envolver
toda a área A, e não que ela seja a variável de integração. Dependendo da forma da
equação para expressar o elemento diferencial de área dA, que depende da função da
velocidade, poderemos ter que efetuar uma integração simples ou dupla.
Uma vez que estabelecemos o fluxo de volume, fica fácil escrever diretamente a
massa desse volume para encontrarmos a equação do fluxo de massa:
 
FMASSA   V . dA 1.21
A
E, para uma grandeza extensiva N qualquer, vale a expressão geral:

 
FN  A
  V . dA 1.22
A equação 1.21 é a forma mais geral para o fluxo de volume, e a 1.22 sua equivalente
para fluxo de grandeza extensiva qualquer transportada pelo fluido.
_____ Sobre o significado do sinal na equação vetorial
Nas equações 1.20, 1.21 e 1.22, válidas em 2 e 3 dimensões, o sinal indica

diretamente se o fluxo é de entrada ou de saída. Devemos lembrar que o sentido do vetor
área é de dentro para fora, quando são definidas superfícies fechadas. Veja os esquemas a
seguir com a superfície fechada de um Volume de Controle.
A1 1 A2
2
V1
Vol. V2
Área de
Entrada
Controle
Área de Saída
Áreas de Entrada Superfície de Controle

Sinal negativo Áreas de Saída
Sinal positivo
dA V
L.C.
V
V
Áreas Laterais dA
dA
Valor nulo
   
QENTRA    V . dA QSAI    V . dA
AE AS
Figura 1.12: definição de áreas de entrada e saída por meio do ângulo entre os vetores.
Toda vez que o ângulo entre os dois vetores for  > 90° o produto escalar será
negativo. Isso só ocorre nos fluxos de entrada.
Por outro lado, um ângulo  <90° indica uma situação entre os vetores que só ocorre
em áreas de saída.
Exemplo 1.8:
Uma trincheira de drenagem intercepta um aquífero numa seção retangular com 2m de
altura, dada em m2 por A = 50 i  25 j . A velocidade de percolação da água na seção
considerada é dada em m/dia por V = 3,0 i  5 j . Calcule o fluxo de volume (vazão) de
água a ser retirada da trincheira, para que a água não se acumule.
Análise:
Trata-se de um caso de velocidade constante ao longo da área, pois a velocidade não é
função de x ou de y. Assim, a equação 1.20 pode ser aplicada diretamente a toda a área:
   
dQ  V  dA  Q  V A
O esquema a seguir permite visualizar a geometria do problema:
y 
V
12,5
0
z 0 x
12,5
y Axz 5,0
A=Lz
25 2 
A yz Ayz i
x 3,0
q

Área A q Axz j
n
-25
Solução:
     
Q  V  A   3 i  5,0 j   50 i  25 j    150  125   275 m 3 / d
Comentários: A solução é teórica porque, na prática, a abertura da vala e o bombeamento
irão alterar as condições de contorno, mudando as cargas e a direção da velocidade nas
proximidades da abertura. Entretanto, o procedimento serve para ilustrar o cálculo, assim
como permite introduzir a discussão sobre o valor negativo do fluxo. Afinal, o que significa
este sinal?
Exemplo 1.9:
A figura mostra o traço de uma seção plana com y
1m de espessura na direção z, perpendicular ao
papel, submetido a um campo bidimensional de
velocidades dado por V = 200 x i + 50 y j (m/s). 1,5m
Determinar o fluxo de volume que atravessa a
1,0m A1
seção A1 indicada na figura.
x
0,5m 1,5m
Solução:
Inicialmente é necessário definir a área A1.em termos vetoriais. Observe o esquema:
Temos:
 
dA  z ds n
  
e também d A  dAyz i  dAxz j ,
dAyz  z ds senq  z dy
dAxz  z ds cosq  z dx
Portanto,
  
dA  z dy i  z dx j
A contribuição da densidade de fluxo na área dA é dada por:
     
dQ  V  dA  200 x i  50 y j   z dy i  z dx j 
Efetuando o produto escalar e lembrando que z = 1
dQ  200 x dy  50 y dx
O fluxo total é a somatória de todas as contribuições ao longo da área A, dada pela
integral de dQ:
x 1, 5 y 1, 5
Q   dQ  
A x 0, 5 
y 1, 0
200 x dy  50 y dx
A integral dupla não pode ser avaliada porque os limites não estão separados. Mas, ao
longo do limite de integração temos que dy = dx/2. Isto transforma a integral dupla em
simples:
x 1, 5
Q
x 0,5
100 x  50 y  dx
Ainda não pode ser avaliada porque sobre a área y é função de x. Para resolver,
temos que notar que: y  0,75  0,5 x . Assim, o fluxo fica:
1, 5
x2
Q   100x  37,5  25x  dx  75
1, 5
 37,5 x  37,5
0,5 2 0,5
Resposta: A vazão que atravessa a seção A1 é de Q = 37,5 m3/s.

1.8. Exercícios Sugeridos
1.8.1. Na seção definida pela figura foram observados os valores de velocidade dados por: V
= (1 + 0,5xy) i + 0,2xy j + 0,2y k , sendo V(m/s) e x e y em metros. Calcular a vazão através
da seção considerada.
y (m) Figura Ex. 2
(0,2) (1,2)
x (m)
(4,0)
1.8.2. A Figura mostra um trecho de um canal regular com seção parabólica. Sabendo que a
velocidade é dada por V = 0,5 z2 i, e que a cota do fundo é dada por Z = 0,5y2, calcular a
vazão transportada.
1.8.3. Considere uma seção de escoamento paralela ao plano XZ, com 2 lados horizontais e
2 verticais. A seção é quadrada com 2m de lado. Um dos lados horizontais está situado em z
= 0m e y = 5m, entre x = 0m e x = 2m. Outro lado horizontal está situado em z = 2m (y=5m).
Esta seção está num escoamento dado por 𝑉 ⃗ = 𝑦𝑧 + 𝑧2 + 𝑥 ⃗ , sendo V em (m/s)
e ordenadas x, y e z em metros. Calcule a vazão através da superfície.
1.8.4. Um fluido escoa através das seções hachuradas do V.C. mostrado na figura. Pede-se:
a) Sendo a velocidade em m/s dada por 𝑉 ⃗ = + + ⃗ , calcule a vazão total que
entra ou sai do V.C.
b) Com as velocidades do item anterior e sendo a concentração de uma substância
dissolvida na água dada por C (mg/L) = 20 mg/L calcule o fluxo de massa total da
substância que entra ou sai do V.C.
c) Com as velocidades do item anterior e sendo a concentração de uma substância
dissolvida na água dada por C (mg/L) = 20x + 20y, calcule o fluxo de massa total da
d) Sendo a velocidade dada por 𝑉 ⃗ = 𝑦 + 𝑧 + 𝑥 ⃗ . Calcule a vazão total que entra
ou sai do V.C..
1.8.5. Um fluido escoa através da seção hachurada do V.C. mostrado na figura. Pede-se:
a) Sendo a velocidade em m/s dada por 𝑉 ⃗ = + + ⃗ , calcule a vazão que entra
ou sai do V.C.
b) Com as velocidades do item anterior e sendo a concentração de uma substância
dissolvida na água dada por C (mg/L) = 20 mg/L calcule o fluxo de massa da
c) Com as velocidades do item anterior e sendo a concentração de uma substância
dissolvida na água dada por C (mg/L) = 20x + 20z, calcule o fluxo de massa da
d) Sendo a velocidade dada por 𝑉 ⃗ = 𝑦 + 𝑧 + 𝑥 ⃗ , calcule a vazão que entra ou sai
do V.C. pela superfície hachurada.
Figura 1
CAPÍTULO 2
DESCRIÇÃO ELEMENTAR DOS ESCOAMENTOS
A geometria de um escoamento qualquer fica completamente descrita pelo seu

campo de velocidades pontuais. Algumas vezes é mais vantajoso descrever um escoa-
mento por meio de outras características cinemáticas, que serão definidas neste capítulo.
Será demonstrado que existem dois métodos fundamentais para descrever um
escoamento: o método Lagrangeano e o Euleriano. Entender as diferenças entre as duas
formas de abordagem e as descrições e características cinemáticas derivadas de cada
método é um dos objetivos deste capítulo.
O capítulo também pretende que o leitor se familiarize com as técnicas analíticas e
experimentais existentes para descrição dos escoamentos.
2.1 Trajetória De Uma Partícula Fluida
O conceito de trajetória é bastante intuitivo. Imagine que você pode marcar uma
determinada partícula do escoamento, e anotar sua posição ao longo do tempo. O resultado
é uma linha definida como trajetória, ou seja, o lugar geométrico ocupado por uma partícula
ao longo do tempo, mostrada na Figura 2.1.
V(t=0)
t=1 t=2 t=3

t=0 t=4
t=5
y
Figura 2.1: Trajetória de uma partícula
As trajetórias podem ser obtidas na prática por método fotográfico, lançando algumas
partículas no escoamento e fazendo exposições sucessivas do mesmo negativo. As
partículas sólidas lançadas no escoamento assumem a função de um traçador, ou seja, de
uma substância que se move com a mesma velocidade do fluido em seu entorno.
A trajetória pertence a uma partícula, que é acompanhada no decorrer do tempo ao
se deslocar pelo escoamento. Por isso se diz que a trajetória é um Conceito Lagrangeano
de descrição do escoamento.
2.2 Velocidade e Aceleração de uma Partícula Fluida
_____ Velocidade de uma Partícula (Lagrange)

O movimento de uma partícula de fluido em um escoamento, pela abordagem
Lagrangeana, é analisado de forma idêntica a um ponto material com uma trajetória
curvilínea. Imagine um ponto movendo-se entre P1 e P2 num plano x – y (Figura 2.2).
y P1
s1 s
P2
s2
j
x
i
Figura 2.2: Vetor posição e vetor deslocamento entre dois pontos

O vetor posição é s em relação à origem. Em P1:
  
s1  x1 i  y1 j 2.1
 
No ponto P2 o vetor posição após o deslocamento pode ser escrito como s  s ,
 
onde s é o vetor deslocamento. Se o deslocamento
 s ocorrer num intervalo de tempo t,

a velocidade média durante o deslocamento é Vmédia  s / t . Sua direção é a mesma do

deslocamento s sobre a corda P1P2.
A velocidade instantânea é calculada tomando-se intervalos de tempo cada vez mais
curtos. Com isso, o comprimento da corda tende a zero, e a direção tende para a tangente à
curva da trajetória em P1.
P4 P4
P V P
3 3
 s (  t3) V
P2 P
 s (  t2) V 2
P1 P1
 s (  t1)
V
Figura 2.3: Velocidade instantânea como limite das velocidades médias
Na Figura 3.3 observamos que, a velocidade instantânea é sempre tangente à

trajetória, e é dada pelo limite:
  
s ds ds 
V  lim   e 2.2
t  0 t dt dt s
O módulo do vetor é a velocidade escalar da partícula, ds/dt, e a direção é dada pelo

versor tangente à trajetória.
A velocidade pode também ser calculada em coordenadas cartesianas, a partir das
projeções sobre os eixos x e y. Enquanto o deslocamento se dá entre P 1 e P2 as

componentes do deslocamento s movem-se entre x1 e x2 e entre y1 e y2.
x dx 
Vx médio   Vx  
t dt   ds  
 V   Vx i  Vy j 2.3
y dy  dt
Vy médio   Vy  
t dt 
A velocidade ao longo da trajetória é a soma vetorial das componentes em x e y, e é

sempre tangente à trajetória.
O método Lagrangeano é a base empregada nas técnicas de medição de campos de
escoamento por imagem de partículas marcadas por traçadores, chamadas de PIV (Particle
Image Velocimetry). Na medição com PIV, técnicas de computação gráfica são usadas para
superpor duas imagens separadas por um intervalo de tempo t. Na imagem superposta,
técnicas de correlação estatística são empregadas para identificar cada partícula na posição
original com sua imagem deslocada, o que permite determinar os deslocamentos x e y e,
portanto, a velocidade da partícula, com as equações 2.3. Com o crescimento da
capacidade dos computadores, a PIV tem sido muito utilizada atualmente.
Também se usa a notação u = Vx ; v = Vy e w = Vz para os componentes do vetor
velocidade. Com esta notação a velocidade num escoamento tridimensional seria expressa
por:
⃗ = 𝑢 + 𝑣 + 𝑤⃗
𝑉
Em termos físicos, a descrição Lagrangeana da velocidade de uma dada partícula
⃗ (𝑥 𝑦 𝑧 𝑡) que descreve a velocidade de
precisa acompanhar a partícula. Assim, na função 𝑉
cada partícula, as ordenadas x, y e z não são fixas. Por exemplo, a velocidade de uma
partícula fluida, originalmente na posição x0, y0, z0 no instante t0, é expressa pelas funções
abaixo:
𝑢 = 𝑉𝑥 = [𝑥(𝑡) 𝑦(𝑡) 𝑧(𝑡) 𝑡]
𝑣 = 𝑉𝑦 = [𝑥(𝑡) 𝑦(𝑡) 𝑧(𝑡) 𝑡]
𝑤 = 𝑉𝑧 = [𝑥(𝑡) 𝑦(𝑡) 𝑧(𝑡) 𝑡]
_____ Aceleração de uma Partícula (Lagrange)
A abordagem Lagrangeana permite calcular facilmente a aceleração de uma partícula

de fluido em um determinado ponto de sua trajetória ao longo do escoamento.
Quando o vetor velocidade tem sua direção constantemente mudada ao longo de
uma trajetória curva, existe uma aceleração mesmo que o módulo da velocidade seja
constante. Pensando no movimento das componentes x e y, vemos que são movimentos
retilíneos acelerados, cuja soma vetorial compõe o movimento ao longo da curva. Deste
modo a aceleração pode ser calculada a partir das componentes em x e y.
d d
ax  Vx ; a y  Vy 2.4
dt dt

 dV   d  d 
a   a xi  a y j  Vx i  V y j 2.5
dt dt dt

 dV d  
a  Vx i  V y j 2.6
dt dt
_____ Coordenadas Intrínsecas
Usando o sistema de coordenadas intrínsecas a aceleração da partícula terá as

componentes tangencial e normal à trajetória, mostradas na Figura 2.4:
en
es
= + 𝑛 𝑛
s
Figura 2.4: Sistema de coordenadas intrínsecas
_____ Aceleração Tangencial

V s V
a s  lim V s  s
t 0 t s
 V s V 1 V2
V
a s  lim s  lim
t 0 t t 0  s t
V
s
 as 
2  s2
2.7
_____ Aceleração Normal
A Figura 2.5 mostra o deslocamento de uma partícula entre P e P', com velocidade
constante em módulo, numa trajetória curva com raio r:
P V
, q
V
V
r V
r
q s
,
P
, q r
V
Figura 2.5: Aceleração normal numa trajetória circular

  
A variação de velocidade entre P e P' é dada por V  V  V ,. Por semelhança de
triângulos, temos que:
V s

V r
s é o comprimento da corda PP', que é aproximadamente igual ao arco de circunferência
entre P e P'. Este comprimento é percorrido pelo ponto em um intervalo t, ou seja,
PP'  V t :
V V t V V2
  
V r t r
A relação torna-se exata no limite, quando t  0 e q  0. Nessas condições temos
a aceleração normal instantânea.
V V2
an  lim 
t  0 t r
A direção e o sentido são os mesmos de V, ou seja, segundo o raio da curva, no
sentido da circunferência para o centro:
 V2 
an   e 2.8
r n
O versor normal aponta sempre para fora da curva, o que explica o sinal negativo. A
aceleração é sempre dirigida em direção ao centro, portanto em sentido contrário ao versor.
2.3 Linha de Corrente - Um Novo Ponto De Vista
Até aqui consideramos uma partícula de fluido, acompanhando seu deslocamento ao

longo do tempo. Entretanto, em muitas análises de escoamentos interessa descrever o
movimento a partir da observação de um ponto fixo no espaço, em uma seção de
interesse.
A análise do escoamento a partir de um ponto fixo no espaço é

denominada Análise Euleriana, ou método de Euler (1707-1783).
Em oposição ao método de Euler, no item 2.1 consideramos uma partícula fluida,

acompanhando-a no espaço, em seu movimento ao longo do tempo, definindo o conceito de
trajetória.
O ponto de vista que considera uma partícula é chamado de

Análise Lagrangeana, ou método de Lagrange (1736-1813).
Se considerarmos um ponto do escoamento e tentarmos descrever as velocidades de

todas as partículas do escoamento que passam pelo ponto especificado (abordagem
Euleriana), teremos uma descrição diferente das velocidades. Não podemos nos valer da
trajetória, pois não estamos mais acompanhando as partículas e cada partícula que
sucessivamente passa pelo ponto de interesse pode ter uma trajetória diferente. Para
desenvolver o equacionamento Euleriano, portanto, é necessário usar o conceito de Linhas
de Corrente
_____ Linhas de Corrente
Uma representação dos escoamentos pode ser obtida quando se traçam linhas
contínuas que são, em cada ponto, tangentes ao vetor velocidade. Essas linhas são
chamadas de Linhas de Corrente. Podem ser obtidas por meio de uma fotografia do
escoamento, onde se lançou um grande número de partículas visíveis. Com um tempo de
exposição apropriado, cada partícula deixará no negativo um segmento correspondente ao
caminho percorrido durante o tempo de exposição, conforme demonstra o esquema da
Figura 2.6.
trajetória de uma
partícula marcada posição no início
do intervalo
posição no final
do intervalo
Linha de Corrente:
V une pontos na direção
tangente a V
V é tangente ao traço
Figura 2.6: Esquema mostrando a técnica de traçado de linhas de corrente
É importante observar que as linhas de corrente descrevem simultaneamente a

direção instantânea de muitas partículas. A Linha de Corrente pertence ao escoamento, ou
seja, descreve as direções do campo de velocidades num dado instante, ao passo que as
trajetórias pertencem a uma determinada partícula ao longo do tempo.
Para não esquecer

Linhas de corrente (Euler): exposição única da foto, muitas partículas
Trajetória (Lagrange): uma partícula, múltiplas exposições do mesmo
negativo
A Figura 2.7 mostra um exemplo real de visualização num escoamento bidimensional

ao redor de um perfil de asa. Nesse caso a água foi marcada com partículas de pó de
alumínio, deixando os traços brancos que se pode ver na fotografia. Com esse apoio pode-
se traçar facilmente as linhas de corrente do escoamento. Mais exemplos dessa técnica
podem ser vistos no Rui Vieira, cap. 1 vol.2, Cinemática.
Figura 2.7: Exemplo de visualização de escoamento para traçar linhas de corrente. Fonte Rui Vieira.
Como a Linha de Corrente é tangente aos vetores de velocidade de todas as

partículas que definem seu traçado, não há escoamento de fluido através de uma Linha de
Corrente. Uma partícula que se move ao longo de uma Linha de Corrente, em qualquer
tempo, tem um deslocamento s com componentes x, y e z, cuja direção é a mesma do
vetor velocidade V. Assim, temos:
𝑥 = 𝑉𝑥 𝑡 ; 𝑦 = 𝑉𝑦 𝑡 ; 𝑧 = 𝑉𝑧 𝑡
𝑥 𝑦 𝑧
𝑡= = =
𝑉 𝑉 𝑉
No limite para dt  0, obtém-se a equação da Linha de Corrente:
𝑥 𝑦 𝑧
= = 2.9
𝑢 𝑣 𝑤
2.4 Velocidade e Aceleração em um Ponto (Análise Euleriana)

Vamos considerar uma Linha de Corrente e as partículas que passam por um ponto P
de um escoamento com velocidade variável no tempo e no espaço. Imagine, por exemplo as
linhas de corrente no interior de um tubo curvo que drena um reservatório de água. Com a
diminuição do nível na caixa, a velocidade diminui em cada ponto, e, ao passar pela curva,
cada partícula sofre uma aceleração que muda a direção de sua velocidade. Essa situação é
esquematizada na Figura 2.8.
VP ( t ) VP (t+ t)
P P
Q VQ ( t ) Q VQ (t+ t)
L.C. L.C.
(a) - tempo t (b) - tempo t +  t
Figura 2.8: Variações da velocidade num escoamento não permanente numa curva de tubulação.
_____ Aceleração Local:
Ao medirmos a velocidade VP num ponto P em dois intervalos de tempo, podemos

verificar uma variação, conforme o esquema vetorial da Figura 2.9.
Variação Local ( no ponto P)
VP (t+ t)
V
- VP ( t ) aLocal  lim
 t 0 t
V
Figura 2.9: Variação local da velocidade: mesmo ponto, dois instantes de tempo.
A taxa de variação da velocidade com o tempo em um dado ponto do escoamento é

chamada de Aceleração Local. Corresponde a uma aceleração, das partículas que passam
pelo ponto P, que ocorre no decorrer do tempo.
Portanto, a aceleração local é dada por
( ) ()
= lim
t 0
= 2.10
sendo V(P, t) e V(P, t + t) a velocidade de duas partículas que passam pelo ponto P nos
dois instantes de tempo considerados. O limite corresponde à derivada da velocidade em
relação ao tempo. Foi usado o símbolo de derivada parcial em 2.10 porque a velocidade
depende também da localização no espaço.
_____ Aceleração Convectiva:
Podemos também observar dois pontos diferentes no mesmo instante de tempo,

como os pontos P e Q da Figura 2.8 e verificarmos que suas velocidades VP e Vq são
diferentes, seja em módulo ou direção ou em ambos, conforme o esquema vetorial da Figura
2.10.
Variação Convectiva (entre ponto P e Q)
VP ( t )
V
V aConvectiva  lim
VQ ( t )  t 0 t
Figura 2.10: Variação convectiva da velocidade: dois pontos no mesmo instante de tempo.
O esquema vetorial mostra que as partículas do escoamento estão sofrendo uma

variação se deslocarem s no espaço entre os pontos P e Q. A taxa de variação no tempo
sentida pela partícula ao se deslocar no espaço é chamada de Aceleração Convectiva.
A variação da velocidade no espaço é dada por:

V = V ( )− V ( )
A partícula sofre esta variação de velocidade no tempo que levou para percorrer a
distância s entre os pontos P e Q com velocidade V, ou seja, t = s/V:
Portanto, a aceleração convectiva é dada por
() () () ()
= lim
t 0
= lim
t 0
2.11
No limite quando t tende a zero o deslocamento também fica infinitesimal, ou seja, s0, e
a equação 2.11 fica:
() ()
= V lim
s0
( ) = V 2.12
Em que dispensamos a identificação do ponto porque os pontos P e Q coincidem no limite.

É interessante observar que a velocidade com que a partícula se desloca entre os
dois pontos influencia na aceleração que a partícula sofre. Mesmo que exista uma grande
diferença de velocidade entre os dois pontos, se a partícula demorar muito tempo no
deslocamento a aceleração convectiva será pequena.
_____ Aceleração Total: Derivada Substantiva
Uma partícula de fluido no escoamento sente simultaneamente as duas acelerações.

Assim, a aceleração da partícula, medida com variáveis com variáveis Eulerianas, é dada
por:
= 𝑜 𝑎𝑙 + 𝐶𝑜𝑛𝑣 𝑡𝑖𝑣𝑎
= + V 2.13
Para tornar clara a distinção entre o uso de variáveis Lagrangeanas e Eulerianas,

usa-se definir a aceleração com o operador que chamamos de derivada substantiva, definido
a seguir:
= + V 2.14
Em que a notação DV/Dt indica que a derivada é uma operação a ser efetuada com
as velocidades de uma determinada partícula da substância em escoamento, ou seja,
variáveis Lagrangeanas. O segundo membro da equação 2.14 é, conforme deduzimos, a
mesma quantidade (aceleração sentida pela partícula) definida com as variáveis Eulerianas
(velocidades medidas em pontos definidos do espaço).
Aceleração de uma partícula
Euler
V V
a   V
t s
Lagrange Local Convectiva
O resultado da equação 2.15 pode ser deduzido também a partir das regras do
cálculo de funções de várias variáveis, pois V = f (s, t), sendo s a coordenada intrínseca que
define a trajetória.
Assim, segundo o cálculo:
V V V V
V  t  s no limite, dV  dt  ds 2.15
t s t s
A aceleração fica então:
 V dt  V ds V V
a    V 2.16
 t dt  s dt t s
Uma apresentação alternativa de 2.13 ou 2.16 é a seguinte:
V  V2
a   ( ) 2.17
t s 2
_____ Aceleração de uma Partícula em Coordenadas Cartesianas
A extensão da equação 2.16 para o caso da velocidade descrita pelas ordenadas do

sistema cartesiano é direta. Uma partícula de fluido num escoamento tridimensional tem sua
velocidade descrita em termos das componentes cartesianas conforme segue:
⃗ = 𝑢 + 𝑣 + 𝑤⃗
𝑉 2.18
Sendo a velocidade uma função de x, y, z e do tempo, a derivada total é dada por:

⃗
𝑉 ⃗
𝑉 ⃗
𝑉 𝑥 ⃗
𝑉 𝑦 ⃗
𝑉 𝑧
= + + + 2.19
𝑡 𝑡 𝑥 𝑡 𝑦 𝑡 𝑧 𝑡
Sendo dx/dt = u; dy/dt = v e dz/dt = w, as componentes da velocidade da partícula.

Com isso a eq. 2.19 fica:
⃗
𝑉 ⃗
𝑉 ⃗
𝑉 ⃗
𝑉 ⃗
𝑉
⃗⃗ = = + + + 2.20
𝑡 𝑡 𝑥 𝑦 𝑧
A equação 2.20 corresponde no sistema cartesiano à equação 2.16 no sistema

intrínseco de coordenadas, para o caso de escoamento tridimensional. A aceleração pode
ser descrita de forma mais compacta, usando a notação do operador gradiente para
descrever as derivadas espaciais,
⃗
𝑉 ⃗
𝑉 ⃗
𝑉
⃗ 𝑉
⃗ = + + ⃗
𝑥 𝑦 𝑧
Vemos que o termo convectivo pode ser escrito como
⃗
𝑉 ⃗
𝑉 ⃗
𝑉
⃗ = (𝑢 + 𝑣 + 𝑤 ⃗ )
⃗ .⃗ 𝑉
𝑉 + (𝑢 + 𝑣 + 𝑤 ⃗ ) + (𝑢 + 𝑣 + 𝑤 ⃗ ) ⃗
𝑥 𝑦 𝑧
Os versores perpendiculares se anulam no produto escalar, restando:
⃗
𝑉 ⃗
𝑉 ⃗
𝑉
⃗ .⃗ 𝑉
𝑉 ⃗ = + +
𝑥 𝑦 𝑧
Assim, em notação vetorial a equação 2.20 fica:
⃗
𝑉 ⃗
𝑉
⃗ = = ⃗ .⃗ 𝑉
+𝑉 ⃗ 2.21
𝑡 𝑡
As três equações escalares correspondentes à equação 2.20, ou à sua forma
compacta, equação 2.21, são dadas por:
𝑢 𝑢 𝑢 𝑢
𝑥 = + (𝑢 + 𝑣 + 𝑤 )
𝑡 𝑥 𝑦 𝑧
𝑢 𝑢 𝑢 𝑢
𝑥 = + (𝑢 + 𝑣 + 𝑤 )
𝑡 𝑥 𝑦 𝑧
𝑢 𝑢 𝑢 𝑢
𝑥 = + (𝑢 + 𝑣 + 𝑤 )
𝑡 𝑥 𝑦 𝑧
2.22
EXEMPLO 2.1:
O escoamento permanente de água através de um bocal convergente, conforme
mostrado na Figura 2.11, pode ser descrito pela abordagem unidimensional com uma
velocidade média que varia em função de x, u(x).
Figura 2.11: Bocal convergente com variação linear da velocidade.
Supondo que a velocidade varie linearmente entre V 0 e 3V0 ao longo do bocal com
comprimento L, pede-se: a) calcule a aceleração como função de x; b) sendo V 0 = 3 m/s e
L = 1m, calcule a aceleração na entrada e na saída do bocal.
Solução:
O problema é unidimensional na direção x: ⃗ = 𝑉𝑥 = 𝑢
𝑉
Com as condições de contorno do problema, a velocidade é dada por:
𝑥
( ) = 𝑉0 ( + )
Portanto,
𝑢 𝑉0
=
𝑥
Item a) Usando a equação 2.18:
𝑢 𝑢 𝑉02 𝑥
= +𝑢 = ( + )
𝑡 𝑥
Item b) Substituindo os valores numéricos:
Na entrada, x = 0,
2
= ( + ) = 𝑠2
Na saída, x = 1m,
2
= ( + ) = 𝑠2
A aceleração equivale a 1,8g na entrada do bocal, atingindo cerca de 50 vezes a

aceleração da gravidade na saída do bocal. Este exemplo mostra que esforços dinâmicos
elevados podem se desenvolver no interior dos escoamentos, mesmo em regime
permanente.
EXEMPLO 2.2:
Encontre a aceleração de uma partícula no campo Euleriano de velocidade dado por:
⃗ (𝑥 𝑦 𝑧 𝑡) =
𝑉 𝑡 + 𝑥𝑧 + 𝑡𝑦 2 ⃗
Solução: Será aplicada a equação 2.22.
2
Inicialmente, identificamos os componentes: u = 3t ; v = xz e w = ty .
Calculamos agora as derivadas parciais necessárias
⃗
𝑉 ⃗
𝑉 ⃗
𝑉 ⃗
𝑉
= + 𝑦2 ⃗ ; =𝑧 ; = 𝑡𝑦 ⃗ ; = 𝑥
𝑡 𝑥 𝑦 𝑧
Aplicando as derivadas conforme a equação 2.22:
= + 𝑦 2 ⃗ + ( 𝑡)𝑧 + (𝑥𝑧) 𝑡𝑦 ⃗ + (𝑡𝑦 2 )𝑥

Colocando os termos juntos:
= + ( 𝑡𝑧 + 𝑡𝑦 2 𝑥 ) + ( 𝑦 2 + 𝑥𝑧𝑦𝑡) ⃗
EXEMPLO 2.3:
Um escoamento permanente bidimensional ocorre no plano xz, sendo x > 0 , z > 0 e A
uma constante, com a velocidade dada por:
⃗ (𝑥 𝑧) = − 𝑥 + 𝑧 ⃗
𝑉
Determine a equação das Linhas de Corrente e a aceleração das partículas.
Solução: Será aplicada a equação 2.9.

Sendo dz e dx o espaçamento entre dois pontos da L.C, temos:
𝑥 𝑧 𝑧 𝑤
= = = = −
𝑢 𝑤 𝑥 𝑢 −
𝑧 𝑥 𝑧 𝑥
=− ∫ = −∫
𝑧 𝑥 𝑧 𝑥
(𝑧) = − (𝑥) + 𝐶
Fazendo a substituição: ln C = C’:
(𝑧) = (𝐶 𝑥) 𝑧𝑥 = 𝐶
Observa-se que as L.C. formam uma família de hipérboles, cujo aspecto é mostrado
na figura a seguir:
5
C= 1
4
C= 2
3 C= 3
C= 4
2
0
0 1 2 3 4 5 6
Figura 2.12: Linhas de Corrente de escoamento nas proximidades de um canto.
2 2
Fica para o leitor demonstrar as acelerações: ax = A x ; az = A z.
_____ Taxa de variação de outras grandezas
O conceito de variação local e convectiva surge sempre que precisarmos avaliar taxas
de variação no tempo de uma grandeza qualquer usando informações Eulerianas.
A Figura 2.13 ilustra o caso com a temperatura sentida pelos ocupantes de um carro
durante uma viagem de Ilha Solteira para São Carlos. A informação Euleriana disponível é a
variação local das temperaturas medidas nas duas cidades. Qual é a taxa de variação média
no tempo, sentida pelos ocupantes do carro, que viaja com janelas abertas?
Temperatura em ISA (°C) Temperatura em S.Carlos (°C)

30 30
25 Variáveis 25
20 Eulerianas 20
15 15
t (h) t (h)
10 10
17 18 19 20 21 22 23 17 18 19 20 21 22 23
Ilha Solteira Variável Lagrangeana

t = 17h
São Carlos
t = 23h
s = 0 km
T = 30 °C s = 420 km
Taxa de variação da T = 10 °C
temperatura no carro = ?
Figura 2.13: Ilustração de cálculo de taxa de variação Lagrangeana com variáveis Eulerianas.
A taxa de variação média (Lagrangeana) é dada por:

− − 𝐶
= = = − 𝐶
𝑡 −
A taxa de variação local em ISA é dada por:
− − 𝐶
= =
𝑡 −
A velocidade média do carro é
𝑉= =
A taxa de variação convectiva da temperatura é dada por:
− 𝐶 − 𝐶
𝑉 = ( ) ( )=
𝑠
Portanto, também no campo de temperaturas do exemplo, a
Variação Lagrangeana (-20°C/6h) = Local (-10°C/6h) mais a Convectiva (-10°C/6h).
2.5 Linha de Emissão e Linha de Tempo
_____ Linha de Emissão

Se injetarmos continuamente um corante num determinado ponto do escoamento,
obteremos ao fim de algum tempo uma figura chamada de Linha de Emissão.
A Figura 2.14 apresenta duas fotografias obtidas no túnel hidrodinâmico do DEM-
FEIS/UNESP, contendo vórtices que se formam a jusante de obstáculos nos escoamentos.
A seta mostra o sentido do escoamento.
ponto de emissão do corante
(a)
vórtice (b)
Figura 2.14: Linhas de emissão de ponto a jusante de cilindro. Fotos: cortesia do Prof. Edson Del Rio.
Todas as partículas marcadas com o corante passaram em instantes anteriores pelo

ponto de injeção do corante. Portanto, o conjunto de vórtices constitui a linha de emissão do
ponto de injeção do corante, identificado na foto (a) pelo ponto branco.
Uma linha de emissão é o lugar geométrico ocupado pelas partículas que passaram
por um dado ponto do escoamento em instantes anteriores. Cada ponto do escoamento
pode ter uma linha de emissão diferente, e como a linha de emissão é uma representação
instantânea, pode variar ao longo do tempo. A forma mais comum de linhas de emissão que
observamos no dia a dia é proporcionada pelas chaminés de fábricas.Os diferentes padrões
de escoamento identificados pelos vórtices das fotos (a) e (b) ocorrem devido à variação da
velocidade do fluido.
O mesmo tipo de estrutura (vórtices a jusante de um cilindro) pode ser visualizado por
meio das linhas de corrente, conforme o exemplo da Figura 2.15. No caso da figura foi
utilizada a técnica do pó de alumínio, obtendo-se a foto com tempo de exposição
suficientemente longo para que as trajetórias apareçam como traços brancos.
Figura 2.15: Vórtices observados pelas linhas de corrente. Obtida de VIEIRA, R.C.C.S.
Você pode observar também com facilidade as linhas de emissão de um bocal de

mangueira de jardim. Movimente a mão de forma ritmada em um percurso fixo, e divirta-se
com os desenhos que o jato forma no ar. Esses desenhos nada mais são que as linhas de
emissão do bocal da mangueira em movimento. Observe como a trajetória de cada gota de
água em particular é completamente diferente da linha de emissão. A Figura 2.16 mostra um
esquema das linhas de emissão que podem ser obtidas com esse simples experimento.
Posições em tempos anteriores
Trajetória de uma partícula

v
t t t t t
1 2 3 4 5
Oscilação do bocal
Trajetória
Trajetória
Linha de Emissão
Trajetória
Figura 2.16: Exemplo esquemático de Linha de Emissão de um bocal oscilante.
A linha de emissão foi desenhada na Figura 2.16 unindo a posição de diferentes

partículas num mesmo instante. As trajetórias extremas e a central são apresentadas em
tracejado e em pontilhado os espaços percorridos pelas partículas a cada ¼ de ciclo do
bocal. Se a frequência do movimento de oscilação variar as linhas de emissão resultantes
descreverão curvas mais abruptas (para aumento da frequência) ou mais suaves (para
diminuição da frequência de oscilação).
_____ Linha de Tempo

A Linha de Tempo é individualizada marcando-se num determinado instante as
partículas alinhadas segundo algum critério de interesse. Diversas linhas de tempo podem
ser obtidas fotografando-se o escoamento ao longo do tempo. A Figura 2.17 mostra
exemplos de duas linhas de tempo, marcadas com (b) e (c), obtidas com injeção de bolhas
de hidrogênio na linha (a).
(a) (b) (c)
Figura 2.17: Exemplo de linhas de tempo construídas com injeção de bolhas de hidrogênio
A Figura 2.17 mostra como linhas de tempo, marcadas pelo método de bolhas de
hidrogênio, podem ser usadas para determinar a diferença de velocidades num escoamento.
As bolhas são geradas por eletrólise da água que ocorre no contato com um fio submetido a
uma corrente polarizada. As bolhas de hidrogênio são carreadas pelo escoamento, servindo
como traçador. Observa-se na linha (a) o lugar onde inicialmente as partículas foram
marcadas. A corrente elétrica foi fornecida durante um intervalo de tempo conhecido,
gerando muitas linhas de tempo que foram carreadas pelo escoamento. A foto só permite
visualizar com clareza a primeira linha de tempo (c), marcada no início do pulso de corrente,
e a última (b), que recebeu as bolhas em (a) no final do pulso de corrente.
2.6 Perfis de Velocidade

A representação gráfica em escala das velocidades ao longo de uma linha
perpendicular à direção da velocidade dá origem a um perfil. Um perfil pode ser obtido
experimentalmente, ou calculado por meio de equações do escoamento.
A técnica utilizada para determinar experimentalmente um perfil depende da escala e
do tipo de escoamento. Uma experiência simples consiste em medir o perfil de velocidades
do ar próximo à superfície da terra. Nesse caso as velocidades podem ser medidas com
anemômetros de conchas. Um resultado possível desse experimento é representado na
Figura 2.18 (a).
Figura 2.18: Exemplos de perfis verticais de velocidade.
Existem escoamentos em que a velocidade varia segundo duas direções, como, por
exemplo, em um rio. As velocidades são menores perto das margens, aumentando em
direção ao centro do rio. Além disso, considerando um determinado ponto do rio, a
velocidade varia na direção vertical, conforme a Figura 2.18 (b). É menor junto ao fundo,
atinge um máximo em algum ponto intermediário e depois decresce ligeiramente até a
superfície (veja figura). As velocidades em rios e canais são medidas com um equipamento
semelhante ao anemômetro de conchas, chamado molinete fluviométrico.
Escoamentos em escala menor são medidos com outras técnicas. Como exemplo,
podemos medir as velocidades do ar no interior de dutos utilizando a anemometria de fio
quente.
2.7 Sistema e Volume de Controle
Um Sistema é uma quantidade definida de massa. Por isso é normalmente mais

fácil quantificar as grandezas de um Sistema usando variáveis Lagrangeanas. Num
escoamento, um sistema pode mudar de forma, acompanhando o escoamento, mas sempre
contém a mesma massa.
O Volume de Controle (VC) é uma região definida do espaço. A massa é livre para
entrar e sair do VC, que então pode possuir massa variável com o tempo.
As grandezas num Volume de Controle são usualmente mais fáceis de quantificar
com variáveis Eulerianas, medidas com instrumentos instalados em pontos fixos do volume
de controle, por exemplo, na seção de entrada e na de saída. A Figura 2.19 traz de forma
esquemática um escoamento qualquer representado por suas linhas de corrente, e um
volume de controle.
Volume de Controle
L.C.
Superfície de Controle
Figura 2.19: Representação esquemática de volume de controle num escoamento.
O Volume de Controle é separado do resto do escoamento pela Superfície de

Controle. A superfície de controle possui áreas de entrada, de saída e áreas laterais, onde
não ocorre entrada nem saída de fluxo.
O Volume de Controle é conceitualmente diferente do sistema. Imagine um
escoamento, visualizado esquematicamente na Figura 2.20 por meio de suas linhas de
corrente. Existe um Volume de Controle qualquer, que contém inicialmente em seu interior
uma massa que constitui o Sistema sob análise (Fig. 2.20(a)). Podemos pensar que o
volume foi demarcado por uma linha de tempo, acompanhando a forma do sistema no
instante inicial.
Decorrido um intervalo de tempo t, o Sistema, carregando consigo sua massa e suas
grandezas extensivas N, terá se deslocado devido ao escoamento, enquanto que o V.C.
permanece fixo no espaço, conforme ilustrado na Figura 2.20(b). Lembre-se que a Linha de
Tempo pode deformar-se ao acompanhar o escoamento, mas sempre demarcará a mesma

quantidade de fluido inicial (com a mesma massa e mesmas quantidades das grandezas
extensivas N - ver propriedades da linha de tempo).
(a) tempo t (b) tempo t + t

Sistema deslocou-se
1 2 3
V.C. = Sistema V.C. permanece fixo

Figura 2.20: Volume de Controle e Sistema em dois instantes consecutivos num escoamento.
Pense no caso da Figura 2.20 e no problema de como descrever o conteúdo de uma

certa grandeza N no sistema, dispondo apenas de variáveis medidas na entrada, na saída e
em pontos fixos no interior do Volume de Controle. No tempo inicial é fácil, pois o Sistema
coincide com o Volume de Controle. Mas, para descrever o sistema após o intervalo t é
necessário pegar o conteúdo do V.C., subtrair a quantidade que entrou (1) e somar a
quantidade que saiu do VC (3). Este artifício matemático é chamado de Transformação de
Reynolds, ou de Relação Sistema-Volume de Controle, e será visto com mais detalhe no
Capítulo 5.
A Conservação de Massa com medições Eulerianas:
Pense no caso da massa do fluido ou de alguma substância carregada pelo fluido em
escoamento. Podemos dizer, em relação ao sistema, usando variáveis Lagrangeanas, que:
|𝑆 = = 𝑆𝑡 𝑡 − 𝑆𝑡 = 2.23
A equação 2.23 expressa o princípio de conservação da massa, sendo a própria

definição de sistema. Agora vamos substituir a massa do sistema pelas massas que ocupam
o Volume de Controle e as regiões do espaço 1 e 3 nos tempos inicial (t) e final (t+t).
|𝑆 = = 𝑆𝑡 𝑡 − 𝑆𝑡 =( 𝑉𝐶 + − )𝑡 𝑡 −( 𝑉𝐶 )𝑡 2.24
Com isso passamos a determinar a massa do sistema pelas variáveis Eulerianas.

Rearranjando,
|𝑆 = =( 𝑉𝐶 𝑡 𝑡 − 𝑉𝐶 𝑡 ) + 𝑡 𝑡 − 𝑡 𝑡 2.25
|𝑆 = = |𝑉𝐶 + 𝑡 𝑡 − 𝑡 𝑡 2.26
Observando as linhas de tempo, percebemos que a massa da região 3 é a que saiu

do VC e a massa da região 1 é a que entrou no VC durante o intervalo de tempo t.
Portanto, as massas 1 e 3 podem ser relacionadas com os fluxos médios de entrada e saída
do VC.
|𝑉𝐶 + 𝑆 𝑡− 𝐸 𝑡= 2.27
Dividindo por t obtemos o valor médio da taxa de variação da massa no VC:
| + 𝑆 − 𝐸 = 2.28
𝑡 𝑉𝐶
A equação 2.28 expressa a conservação da massa com variáveis Eulerianas,

observadas no interior do VC e nas Superfícies de Entrada e Saída.
2.8 Classificação dos Escoamentos
_____ Escoamentos Unidimensionais, Bidimensionais e Tridimensionais.
Em algumas situações, como no caso de tubos com escoamento em altas

velocidades, o perfil de velocidades é praticamente constante, podendo ser desprezadas as
variações na seção para fins práticos. Nessas condições o escoamento é determinado
apenas pela velocidade média na seção do escoamento, originando os chamados
escoamentos unidimensionais, ou 1D. Num escoamento 1-D basta saber a ordenada da
seção para determinarmos a velocidade no ponto desejado.
Observando os diversos exemplos de escoamento que nos rodeiam, podemos
verificar que em alguns deles, devido a condições de simetria, basta apenas um perfil de
velocidades para descrever o escoamento. Esse é o caso do escoamento a baixas
velocidades em dutos circulares, e do escoamento do ar sobre a superfície plana da Terra.
Não importa ao longo de qual linha foram obtidas as velocidades, os perfis resultantes serão
idênticos. Esses escoamentos são chamados bidimensionais ou 2-D. Para determinar a
velocidade num escoamento 2-D precisamos conhecer duas coordenadas do ponto
desejado.
No caso de um rio o escoamento não fica totalmente determinado apenas com um
perfil vertical, pois os perfis variam conforme a distância da margem. Esses casos definem
os escoamentos tridimensionais, ou 3-D. Para determinar a velocidade num escoamento 3-D
é necessário saber as coordenadas x, y e z do ponto desejado.
_____ Outras Classificações dos Escoamentos
Aprendemos que quando o critério de análise dos escoamentos é o número de

variáveis necessárias à descrição do campo de velocidades obtemos sua classificação em
uni, bi e tridimensionais. Outros critérios podem ser utilizados.
Massa específica
Utilizando a massa específica, podemos dividir os escoamentos em compressíveis,
quando a massa específica varia de uma seção para outra ou incompressíveis, quando
permanece constante. Observe que a classificação pertence ao escoamento, e não ao
fluido: um mesmo fluido pode participar de escoamentos compressíveis e incompressíveis,

dependendo dos gradientes de pressão observados.
Quando a massa específica varia em uma mesma seção do escoamento temos os
chamados escoamentos estratificados, em oposição aos não estratificados, em que a massa
específica é constante na seção.
Tempo
Quando o critério de classificação é o tempo, devemos escolher uma dada seção do
escoamento e observar o que ocorre. Se as grandezas não variarem temos um escoamento
permanente, em oposição aos não-permanentes, também chamados de transientes.
Comparação entre seções
Num dado instante de tempo podemos considerar o comportamento de uma grandeza
em duas seções. O mais comum é utilizarmos a velocidade para esta análise, dando origem
a duas situações: escoamento uniforme, quando não há variação de módulo, direção e
sentido da velocidade, e não uniforme ou variado quando ocorre o oposto.
Comportamento dinâmico
Divide os escoamentos entre Laminares e Turbulentos. Pela sua importância na
análise dos escoamentos a divisão entre laminares e turbulentos é abordada com mais
detalhe a seguir.
_____ Escoamentos Laminares e Turbulentos
O comportamento dinâmico dos escoamentos é caracterizado pelas perturbações

introduzidas pelas forças de inércia e pelas forças viscosas que se opõem às perturbações.
Os escoamentos laminares ocorrem nas situações em que as forças viscosas são
maiores que qualquer força de inércia que aparece tentando tirar as partículas do
alinhamento. Nesses escoamentos os fluidos escoam em camadas perfeitamente definidas,
como lâminas superpostas, e com velocidades também perfeitamente definidas.
A Figura 2.21 mostra um esquema de duas linhas de emissão em escoamento
laminar. As linhas nunca se cruzam, pois o fluido escoa como se fossem lâminas deslizando
umas sobre as outras. Devido a essa característica ordenada, no escoamento laminar só há
difusão de quantidade de movimento a nível molecular entre as camadas adjacentes de
fluido.
Agulhas injetando
corante
Linhas de Emissão
paralelas
a) Escoamento Laminar em baixas velocidades
Agulhas injetando
corante
Linhas de Emissão
não se misturam
b) Escoamento Laminar aproximando-se da transição p/ turbulento
Figura 2.21: Esquema de linhas de emissão em escoamento laminar.
Quando aumenta a velocidade as forças de inércia começam a aumentar em relação

às forças de inércia, e o escoamento passa por uma fase de transição, começando a ficar
cada vez mais desordenado, até chegar ao escoamento turbulento. A fase de transição é
ilustrada pelas linhas de emissão da Figura 2.21-b.
Quando as forças de inércia são grandes em relação às forças provocadas pela
viscosidade os escoamentos apresentam um comportamento desordenado, com mistura de
quantidades macroscópicas de fluido na direção transversal à velocidade média. Nesses
escoamentos a velocidade instantânea apresenta flutuações aleatórias em torno de um valor
médio, e são chamados de escoamentos turbulentos.
A estrutura dos vórtices em um escoamento turbulento pode ser visualizada pela linha
de emissão da Figura 2.22. Foi injetado traçador líquido num escoamento turbulento em um
canalete de laboratório, por meio de uma agulha, observando-se a mistura na direção
transversal ao escoamento provocada pela turbulência.
Ponto de Injeção
Velocidade Média
Figura 2.22: Detalhe de uma linha de emissão em escoamento turbulento de água em canal.
A principal característica dos escoamentos turbulentos é a maior capacidade de

misturar as grandezas, devido ao movimento macroscópico provocado pelos vórtices. A
movimentação dos vórtices é equacionada como uma flutuação de velocidade aleatória, e de
média nula, que se superpõe à velocidade média.
𝑉 = 𝑉+𝑣
_____ Outras Grandezas que descrevem os escoamentos
Outras grandezas são necessárias para caracterizar completamente um

escoamento, e sua enumeração depende do tipo de problema.
Considerando o escoamento de um fluido sem misturas e isotérmico, as outras
grandezas necessárias são a pressão, a massa específica e a cota geométrica, que
definem, conjuntamente com a velocidade, a energia mecânica total.
No caso de um escoamento destinado a resfriar uma determinada peça ou
equipamento, a temperatura passa a ser importante também.
Em escoamentos que envolvem misturas de substâncias, é necessário acrescentar
a concentração para descrever completamente o problema.
CAPÍTULO 3: CONSERVAÇÃO DE GRANDEZAS

Equação da Continuidade
Aplicaremos aqui o conceito familiar de conservação de massa, discutindo como sua

aplicação em fenômenos de transporte é possível na análise Euleriana, ou seja, usando
volumes de controle e não um sistema.
Um caso prático interessante é a operação de reservatórios de geração de energia. É
necessário manter um contínuo controle do volume represado, para enfrentar as épocas
sem chuvas, controlar as cheias e obter um rendimento ótimo de turbinas.
3.1 Conservação da Massa

A região escolhida para a análise é chamada de Volume de Controle. Um Volume de
Controle é uma porção definida do espaço onde se dá o escoamento, com quantidade de
massa que pode variar, conforme a Figura 3.1:
Fe
m As
Ae
Fs
Figura 3.1: Fluxos de massa num volume de controle.
Podemos dizer, para um determinado intervalo de tempo t:

m INICIAL + m ENTRA – m SAI = m FINAL
m FINAL – m INICIAL = m = m ENTRA – m SAI 3.1
A equação 1.1 está ligada a um intervalo de tempo t. Como os escoamentos são
contínuos, é mais conveniente escrever as taxas médias no intervalo de tempo, dividindo a
equação 3.1 por t:
m mE mS Os termos do segundo membro são Fluxos

  de Massa médios no intervalo t.
t t t
m
 F Me  F Ms   Q e   Q s 3.2
t
A equação 3.2 é exata quando as vazões não variam no tempo t ou quando são
usados os valores médios no intervalo de tempo. No caso de fluxos variáveis, é necessário
usar um valor instantâneo, obtido pelo limite da variação da massa quando t tende a zero.
dm
 F Me  F Ms 3.3
dt
A equação 3.3 diz que “A taxa instantânea de variação da massa é igual ao
saldo dos fluxos de entrada e saída”
_____ balanço de volumes
No exemplo da usina hidrelétrica e em muitos casos da prática a massa específica
não varia. Esse é o caso dos chamados escoamentos incompressíveis. Para  constante, o
balanço de massas fica equivalente a um balanço de volumes:
m Vol
   Qe   Qs 3.4
t t
Vol
 Qe  Qs 3.5
t
A equação 3.5 diz que “A taxa média de variação do volume é igual à diferença de
vazões de entrada e saída”
_____ fluxos variáveis no tempo
Quando a vazão varia no tempo, pode-se usar o mesmo raciocínio, mas escrevendo a
fórmula com valores instantâneos.
dVol
 Qe  Qs 3.6
dt
Nesse caso, cada pequeno intervalo diferencial de tempo dt traz uma variação
diferencial no volume dVol.
A variação total num dado tempo finito é o somatório das variações diferenciais ao
longo do intervalo de tempo considerado. Analiticamente, isto se consegue pela integração
da equação diferencial da equação 3.6:
V = ( − )
t t
Vol TOTAL  
t0
dVol  
t0
(Qe  Qs ) dt
3.7
_____ aproximação numérica

A integral da equação 3.7 muitas vezes não tem solução analítica. Para casos em que
a solução analítica não existe ou é inconveniente, pode-se chegar ao resultado por
aproximação numérica. Para isso divide-se o tempo total numa sucessão de intervalos de
tempo finitos t, no qual as vazões são consideradas constantes. Obtemos então:
n n
  Voli   (Q
t0  n  t
VolTOTAL  Vol t0 E  QS )  t 3.8
i 1 i 1
É claro que a equação 3.8 possui um erro, pois as vazões variam continuamente. A
aproximação torna-se mais fina à medida que decresce o intervalo de tempo considerado,
ou seja, aumenta o limite n do somatório de parcelas finitas. A equação só é totalmente
exata no limite para t  0. O limite da série infinita de somas é, como sabemos,
equivalente à integral:
 
Vol TOTAL  lim  Voli  
t
 t 0  
 t0
dVol 3.9
 i 1 
Assim, como em vários casos da prática, a quantidade n de somas a ser efetivamente
realizada depende dos objetivos do cálculo. Usualmente, quanto mais rápida a variação dos
fluxos no tempo, menor deve ser o intervalo de tempo t adotado.
_____ balanço de massas

Se for conveniente lidar com a massa, a equação 3.7 fica:
t t
mTOTAL  t0
eQe dt  
t0
sQs dt 3.10
O segundo membro foi dividido em duas parcelas porque num caso geral a massa
específica pode variar. Já o caso da variação de volume só tem sentido em escoamento
incompressível, e dividir ou não os termos de entrada e saída dos fluxos é apenas uma
questão de conveniência e/ou clareza.
EXEMPLO 3.1:
Um reservatório prismático com área da base Ab = 10m2 possui um volume inicial de 20 m3 e
recebe, durante 1 hora, uma vazão média de 10m 3/h, fornecendo uma vazão média de
30m3/h. Com esses valores médios o reservatório fica completamente vazio após 1 hora.
Determine a dinâmica da variação de nível no reservatório por meio de um gráfico do nivelo
em função do tempo, para os casos seguintes.
___ Caso 1: vazões constantes
Temos apenas um intervalo de tempo para aplicar a equação 3.5:
Vol = (10 – 30) t = -20 t, e a variação do nível é linear, conforme os gráficos.
Q(m3/h) Nível (m)
QS
30 3,0
20 2,0
QE
10 1,0
t (h) t (h)
0 0,5 1,0 0 0,5 1,0

___ Caso 2: vazão constante de entrada e saída de 60 m 3/h durante a última meia hora:
Neste caso temos dois intervalos de 0,5h, em que as vazões são constantes, nos
quais podemos usar a equação 3.5:
Vol1 = (10 – 0) t = 10(m3/h) G 0,5 h = + 5 m3  Vol1 = Vol0 + Vol1 = 25m3
Vol2 = (10 – 60) t = -50(m3/h) G 0,5 h = - 25 m3  Vol2 = Vol1 + Vol2 = 0m3
Usando a equação 3.5 em 2 intervalos de tempo consecutivos chegamos ao mesmo
resultado final do caso 1 (reservatório vazio). Isto é lógico, uma vez que as vazões médias
não variaram.
Observamos que neste caso o nível aumentou durante o primeiro intervalo, como
demonstra o gráfico a seguir.
Q(m3/h)
Nível (m)
QS
60 3,0
2,0
QE 1,0 t (h)
10 t (h)
0 0,5 1,0
0 0,5 1,0
A equação 3.7 fornece a variação total nos dois intervalos de tempo, levando ao
mesmo resultado final:
VolTOTAL = Vol1 + Vol2 = + 5 – 25 = - 20 m3
VolFINAL = VolINICIAL + VolTOTAL = 20 – 20 = 0m3
Comentário: Os casos 1 e 2 ilustram que podemos chegar ao mesmo resultado final, mas
com dinâmicas diferentes. Somente com a divisão em dois intervalos de tempo foi possível
captar a variação real de nível no caso 2.
Assim, para conhecer em detalhe a evolução dos níveis de água no reservatório, com
vazões que variam continuamente, é necessário utilizar intervalos de tempo cada vez
menores. No limite, chega-se à equação 3.9. O exemplo 3.2 mostra um caso prático da
situação de vazões variáveis.
EXEMPLO 3.2:
Um reservatório prismático com área da base Ab = 10m2 recebe uma vazão constante de
10m3/h. A vazão de saída em m3/h é dada por Qs = 10H, sendo H a cota do nível da água
em metros. Considerando que inicialmente a água está na cota H i = 5,0m, Calcular o nível
após decorrido 1 hora.
Qe = 10 m3 /h
dVol = A B dH
H
Qs = 10H m3 /h
Análise: como neste problema a vazão de saída varia continuamente, a solução exata
precisa partir do balanço instantâneo. A variação total de volume será encontrada pela
integração do balanço instantâneo.
Solução:
𝑉𝑜𝑙 𝑉𝑜𝑙
= −  = −
𝑡 𝑡
A equação diferencial resultante não pode ser integrada porque há 3 variáveis (Vol, H e t).
Usando a relação entre volume e altura do nível d’água podemos reduzir a 2 variáveis:
𝑉𝑜𝑙 =
= −  = 𝑡
𝑡 −
Com as variáveis separadas podemos integrar
∫ = ∫ 𝑡
− 0
Fazendo a mudança de variáveis: 𝑢= −  𝑢 = −

𝑢 𝑢 𝑢
−∫ 0 𝑢
= ∫0 𝑡 − 𝑢 |𝑢 0 = 𝑡 |0  0
= −
𝑢
=  𝑢= −
0
− = −  =
Resposta: ao final de 1 hora o nível será de 2,471m e o volume será 24,71m3.
EXEMPLO 3.3:
Resolva o problema 3.2 utilizando método numérico aproximado (eq. 3.7) utilizando t de 6
minutos (0,1h).
Análise: o método numérico considera que os fluxos se mantêm constantes durante o
intervalo. Evidentemente esta hipótese contém um erro, já que o fluxo de saída varia
continuamente. Esta é a razão pela qual a solução numérica contém erros, sendo apenas
uma aproximação da solução.
Para o primeiro intervalo de tempo (i = 1):

Qs = 10  5 = 50m3/h
𝑉𝑜𝑙 = ( − ) = −  Vol1 = 46m3  H1 = 4,6m
Segundo intervalo de tempo (i = 2)
Qs,2 = 10  4,6 = 46m3/h
𝑉𝑜𝑙2 = ( − ) = −  Vol2 = 42,4m3  H2 = 4,24m
Resumindo para os demais intervalos:
Qs,3 = 10  4,24 = 42,4m3/h
𝑉𝑜𝑙2 = ( − ) = −  Vol3 = 39,16m3  H3 = 3,916m
𝑉𝑜𝑙 = ( − ) = −  Vol4 = 36,2m3  H4 = 3,620m
𝑉𝑜𝑙 = ( − ) = −  Vol5 = 33,58m3  H5 = 3,358m
𝑉𝑜𝑙 = ( − ) = −  Vol6 = 31,222m3  H6 = 3,1222m
𝑉𝑜𝑙 = ( − ) = −  Vol7 = 29,0998m3  H7 = 2,91m
𝑉𝑜𝑙 = ( − ) = −  Vol8 = 27,19m3  H8 = 2,719m
𝑉𝑜𝑙 = ( − ) = −  Vol9 = 25,471m3  H9 = 2,5471m
𝑉𝑜𝑙 = ( − ) = −  Vol10 = 23,9239m3  H10 = 2,392m
Resposta: a solução numérica iterativa com intervalos de tempo de 0,1h indica um nível
d’água no reservatório de 2,392m ao fim de 1 minuto.
Obs: A solução numérica apresentou erro de -3,2%, que tende a diminuir com intervalos de
tempo menores.
3.2 Misturas Homogêneas - Balanço de Grandeza Extensiva N
Normalmente a água em escoamento não é pura, mas misturada com várias

substâncias sobre as quais é possível ter informação por meio do balanço de massas. Vimos
no item anterior que a taxa de variação da massa é igual à diferença entre os fluxos de
entrada e saída. Esse raciocínio é válido não só para a massa do fluido em escoamento,
mas também para as outras massas e grandezas dependentes da massa que o fluido
transporta em seu meio.
= −
𝑁
= 𝑁 − 𝑁 3.11
𝑡
Mas os Fluxos responsáveis pelas quantidades Ne e Ns da grandeza extensiva dependem

do fluxo de massa e das concentrações, pois:
=  = 
Então, para valores médios no intervalo de tempo,
𝑁
=  −  3.12
𝑡
Lembrando que quando os fluxos variam no tempo, o balanço deve ser escrito para um
instante de tempo:
𝑁
=  −  3.13
𝑡
EXEMPLO 3.4:
Um meio para determinar a vazão de rios consiste na injeção de substâncias traçadoras,
como sais ou corantes. Numa determinação de vazão em um córrego foram lançados 2L/s
de água com uma concentração de corante fluorescente igual a 5mg/L. Numa seção a
jusante, após a completa mistura do traçador, retirou-se uma amostra da água, obtendo-se
uma concentração de corante de 0,2mg/L. Qual a vazão do córrego?
FN,1 Q2
FN,2
Q1
VC FN,3
Q3
Análise: trata-se de um exemplo de balanços em regime permanente em que ocorre a

conservação da massa (volume) da água do rio e também da massa do traçador no VC
dado pelo trecho de rio. Será usada a eq. 3.13 para N = Volume, e também para N = massa
do traçador.
Solução:
= −
Regime Permanente: = = −
- Grandeza N = Volume
= −  + 2 =
- Grandeza N = massa do traçador:
= −  𝐶 + 𝐶2 2 =𝐶
O traçador não é encontrado naturalmente no rio: C1 = 0. Com os valores numéricos dados

os fluxos de massa do traçador ficam:
+ = ( + )
𝑠 𝑠
Que resolvida fornece Q2 = 48 L/s.

Resposta: a vazão do córrego é de 48 litros por segundo (0,048m 3/s), equivalente a
162m3/h.
EXEMPLO 3.5:
Considere um recipiente com 100 litros de água à temperatura de 20C, recebendo 1 L/s de
água a 80C e com uma vazão de saída de 1L/s. Durante os instantes iniciais a água sai
com temperatura de 20C, mas depois disso a água sai com a temperatura média do
reservatório. Isto quer dizer que a água no reservatório é bem misturada. Analise o
transiente da temperatura da água no reservatório, supondo que a variação da massa
específica da água com a temperatura é desprezível. Calcule a variação da quantidade de
calor armazenada em 10 segundos e a temperatura média da água na caixa ao final deste
período. Adote ρ = 995 kg/m3 e calor específico c = 4,18kJ/kg°C.
Fe
AS
Ae FS
Solução:
Parte a) _____ Análise do transiente
A grandeza extensiva considerada é a quantidade de calor:
𝑁
= ( − 0)   = = ( − 0)
Como a temperatura da água varia continuamente, o fluxo de calor que deixa o reservatório
também varia, de forma que o balanço de energia só é válido instantaneamente, e deve ser
expresso na forma diferencial:
𝑁 ( 𝑇)
𝑡
= −  𝑡
= −
em que a temperatura de referência adotada foi T0 = 0.
O balanço de massa da água no reservatório mostra que = , pois o regime é

permanente ( não há variação da massa no reservatório) e com massa específica constante
(não há variação do volume no reservatório).
Assim, o balanço de energia fica:
𝑇
= ( − )
𝑡
Fazendo =( − ) a equação do balanço pode ser escrita como:
− 𝑡
=
𝑄
𝑉𝑜𝑙
 = − 𝑡
𝑉𝑜𝑙
Lembrando que o tempo de detenção hidráulico, definido como = , é o tempo
𝑄
médio de permanência de cada partícula no reservatório.
A equação acima pode ser integrada facilmente, entre um instante de tempo inicial ti
em que a temperatura do reservatório é Ti e um tempo t qualquer, obtendo-se:
− 𝑡 𝑡
= =
𝑖 − 𝑖
Que é a solução para o transiente de temperatura no reservatório. Quando o tempo tende ao
instante inicial, T tende a Ti e q  1; por outro lado, quando o tempo é muito grande, q0 e
T  Te.
Observe que a solução não depende da massa específica nem do calor específico do
fluido no reservatório, apenas depende da relação entre o volume do reservatório e a vazão
de alimentação, dado pelo tempo de detenção hidráulico. Desta forma a solução vale para
qualquer problema semelhante, podendo ser expressa pelo gráfico adimensionalizado da
figura seguinte.
1,0
0,8
0,6
q / qi
0,4
0,2
0,0
0 1 2 3 4
t/ 
Resposta adimensional da temperatura no reservatório
Parte b) ____ cálculo da temperatura após 10 segundos

Basta substituir os valores numéricos na equação do transiente
− 𝑡 𝑡
= =
𝑖 − 𝑖
𝑄 00
sendo o tempo de detenção dado por: = = = 𝑠 . Assim,
𝑉𝑜𝑙
0⁄ 00
− ( 𝑠) = ( − )
( 𝑠) = − = 𝐶
Resposta: a temperatura da água no reservatório após 10 segundos será de 25,7°C.
EXEMPLO 3.6:
Calcule numericamente um valor aproximado para a temperatura da água do reservatório do
exemplo anterior ao final de 10 segundos.
Solução:
Para utilizar uma aproximação numérica da equação diferencial do balanço de energia,
vamos considerar que os fluxos permanecem constantes ao longo do intervalo de tempo
finito adotado.
Em função dos fluxos de massa e volume, a variação da quantidade de calor armazenada é

calculada por:
= = = ( − ) 𝑡
= ( − )𝐶 𝑠=
𝐶 𝑠
Essa quantidade de calor, acrescentada à massa da caixa, permite calcular o acréscimo de
temperatura:
= = 𝐶
𝐶
Resposta: A temperatura da caixa ao final do intervalo de 10 segundos será de 26,0C,
aproximadamente.
Comentários: Note que o erro da solução numérica aproximada foi de 0,3°C em relação à
solução analítica. Este erro tende a diminuir com a adoção de menores intervalos de tempo.
Se o cálculo for repetido para mais um intervalo de 10 segundos o acréscimo de
temperatura será menor, porque maior quantidade de energia deixa o reservatório, visto que
a água sairá mais quente.
A solução numérica partiu da premissa que a água no reservatório permanece a 20°C
durante todo o intervalo de tempo considerado. Evidentemente, esta simplificação implica
em que a resposta possui certo erro. Entretanto, o erro tende a diminuir com o intervalo de
tempo considerado no cálculo.
EXEMPLO 3.7:
Considere novamente o reservatório do exemplo 3.5, levando em conta a variação da massa
específica com a temperatura. Calcule numericamente a temperatura da água ao final de 10
segundos.
Fe
Dados:  (80C) = 971,8kg/m3
AS  (20C) = 998,2kg/m3
Ae FS c = 4,18 kJ/kgC
Análise: este problema envolve o balanço transiente de massa da água, visto que, com o
aumento da temperatura causado pela entrada de água quente, o reservatório de volume
constante irá conter cada vez menos massa; além disso, a outra grandeza envolvida é a
energia, também com balanço transiente, visto que a saída de água cada vez mais quente
conduzirá cada vez mais energia para fora da caixa e o fluxo de entrada de energia térmica
é constante.
Em vista dessas dificuldades, torna-se mais simples resolver o problema por aproximação
numérica.
Solução:
A grandeza extensiva considerada é a quantidade de calor: N  mc ( T  T0 ). A variação da
quantidade de calor no período é dada por:
= 𝑖𝑛𝑎𝑙 − 𝑛𝑖 𝑖𝑎𝑙 = (𝑡 + 𝑡) − (𝑡)
mas também podemos expressar a variação por meio dos Fluxos:
= 𝐸𝑛𝑡 𝑜𝑢 − 𝑆𝑎𝑖𝑢 = 𝑡− 𝑡
Expressando agora os fluxos de energia em função dos fluxos de massa e volume,
𝑁 = 𝑉 ; = ( − 0)  𝑁 = ( − 0 )𝑉
Considerando a temperatura de referência T 0 como nula e substituindo os valores
numéricos, temos as quantidades de energia que entraram e saíram:
= ( ) ( ) ( 𝐶) ( ) ( ) (𝑠) =
𝐶
= ( ) ( ) ( 𝐶) ( ) ( ) (𝑠) =
𝐶
Conhecendo a quantidade de calor aduzida, o cálculo da variação da temperatura depende

da quantidade de massa no reservatório. A massa final é obtida pelo balanço de massas:
m  me  ms  m  ( Fe  Fs )t
m  ( 971,8 x 0,001  998,2 x 0,001 )60  1,584kg
m f  mi  m  98,416kg
=  = = 𝐶
𝐶
Resposta: Portanto, a temperatura ao final de 10 segundos será 25,9 C.
Para continuar o cálculo da aproximação numérica ao longo do próximo intervalo de

tempo precisamos da massa específica da água à temperatura de 25,9°C.
Continuam valendo as observações dos exemplos anteriores. Ao assumir que as
temperaturas permanecem constantes ao longo de todo intervalo de tempo a solução
numérica carrega um erro de aproximação. Como antes, este erro tende a diminuir com o
intervalo de tempo considerado.
3.3. Equação Integral do Balanço de Massa
Nosso próximo desafio será descrever matematicamente o balanço, de forma que

seja válida para uma situação geral, admitindo-se velocidades variáveis e com inclinação
qualquer em relação às seções de entrada e saída da região de interesse, que é chamada
de Volume de Controle.
_____ Volume de Controle
O Volume de Controle (VC) é uma região definida do espaço. A massa é livre para
entrar e sair do VC, que então pode possuir massa variável com o tempo. As grandezas num
VC são usualmente mais fáceis de quantificar com variáveis Eulerianas. A Figura 3.2 traz de
forma esquemática um escoamento qualquer representado por suas linhas de corrente, e
um volume de controle.
O Volume de Controle é conceitualmente diferente do sistema. Um Sistema é uma
quantidade definida de massa. Por isso é normalmente mais fácil quantificar as grandezas
de um sistema usando variáveis Lagrangeanas. Num escoamento, um sistema pode mudar
de forma, acompanhando o escoamento, mas sempre contém a mesma massa.
Volume de Controle
L.C.
Superfície de Controle
Figura 3.2: Representação esquemática de volume de controle num escoamento.
Um sistema não é muito útil para efetuar balanços de massa em escoamentos. Todas
as equações de balanços usam como base um volume de controle. Isso ocorre porque não
há interesse, por exemplo, em saber “qual água” se encontra num reservatório, e sim
“quanta água?” é a pergunta importante. Salvo raras exceções, toda água é igual do ponto
de vista da solução de problemas práticos. Por exemplo, no caso de hidrelétricas, é
necessário operar o reservatório com segurança e produção ótima de energia elétrica.
Portanto interessa equacionar o problema a partir do reservatório, uma região fixa do espaço
(VC), e não a partir de cada massa de água que escoa no rio (Sistema).
Outra razão, talvez menos óbvia, é que as variáveis que utilizamos para descrever os
escoamentos são Eulerianas. Isso significa que foram medidas em pontos definidos do
espaço, e não em partículas definidas de matéria.
____ Balanço
Para o balanço geral de massa vamos imaginar um VC com fluxos que entram e
saem por várias seções de entrada e saída. Inicialmente vamos simplificar o efeito de todas
as entradas em uma só e de todas as saídas também em uma só.
Vamos pensar agora no efeito dos fluxos de entrada e saída do nosso VC. Se isso
parece difícil com um VC abstrato, podemos pensar no conhecido problema do reservatório,
conforme a Figura 3.3, com o mesmo resultado:
Área Lateral
FM,e
V.C. M
t
Área de Saída
Área de Entrada
FM,s
Figura 3.3: Reservatório atuando como volume de controle num escoamento.
Sabemos que o balanço de massa pode ser expresso por:

M
FM ,E  FM ,S  3.14
t
Sabemos também que os fluxos podem ser expressos pela integral vista no item 3.1.
Mas, para substituirmos sem erro na equação do balanço, é necessário lembrar do sinal
algébrico incluído na integral. O sinal é consequência do produto escalar da velocidade e
área, devido à convenção de sentido para o vetor normal. Veja o esquema da Figura 3.4:
Áreas de Entrada Superfície de Controle

Sinal negativo
L.C.
dA V
V
Áreas Laterais V
dA Áreas de Saída
Valor nulo dA
Sinal positivo
Figura 3.4: Sentido do vetor área em relação a densidades de fluxo de entrada e saída.
Outro resultado interessante da integral é que ao longo da superfície de controle

existem apenas áreas de entrada, de saída ou áreas laterais. Uma área lateral é aquela
onde não há fluxo de entrada ou saída, como se a superfície fosse impermeável. Nas áreas
laterais o vetor velocidade sempre fica perpendicular ao vetor área, fazendo com que o
produto escalar seja nulo.
Com as propriedades da integral em mente, podemos escrever então que:
 
FM , E    V . dA 3.15-a
AE
 
FM , S    V . dA 3.15-b
AS
onde os limites AE e AS nas integrais referem-se às áreas de entrada e saída,

respectivamente.
Utilizando então a notação geral de fluxos de entrada e saída das equações 3.15 na
expressão do balanço, temos:
    M
 
AE

 V . dA 
  AS


 V . dA 
t 3.16
F M ,E F M ,S
Podemos agora levar as integrais para o segundo membro, e também acrescentar um
termo nulo, sem alterar o balanço:
M      
   V . dA    V . dA    V . dA  0 3.17
T AE
  AS
  AL

 FM , E FM ,S NULO
Uma propriedade interessante da integração nos permite somar todos os limites de

integração numa mesma integral, de forma que podemos escrever:
M  
   V . dA  0 3.18
t   AS 
AE AL
SC
A equação 3.18 é válida apenas se os Fluxos permanecerem constantes durante o

intervalo de tempo t considerado. Se a velocidade variar alteram-se os fluxos e o balanço
que escrevemos deixa de ser válido. Para que a ideia fique exata é preciso pensar na
variação de massa que ocorre em um tempo t muito pequeno. Somente nesse caso o fluxo
instantâneo é igual ao médio, mesmo nos regimes transientes.
Escrevendo então o balanço para o caso do limite de t  0 temos:
dM  
dt
  SC  V  dA  0 3.19
A equação 3.19 utiliza valores instantâneos das velocidades e da taxa de variação da

massa. Portanto, continua válida mesmo durante os transientes.
_____ Equacionando a taxa de variação da massa no VC ( = termo dM / dt )
Para encerrar o balanço integral de massa num V.C. falta explicitar como, de uma
forma geral, é calculada a massa contida no VC.
Se imaginarmos uma massa específica uniforme ao longo de todo o volume, fica
muito fácil:
= 𝑉𝑜𝑙 onde Vol é o volume total.
Desta forma a taxa fica:
= ( V ) 3.20
A equação 3.20 descreve a variação de massa dos chamados “modelos

concentrados”, em que apenas um valor de qualquer das grandezas consideradas
descreve toda a massa considerada. O equacionamento concentrado de uma dada
grandeza pode ser adotado para um Volume de Controle ou para um Sistema.
Num caso geral a massa específica pode variar ao longo do volume de
controle. Basta, por exemplo, que varie a temperatura ou a salinidade do fluido para que
cada ponto do fluido tenha um  diferente. Como avaliar a massa total? Nesses casos é
necessário usar os “modelos distribuídos”. O modelo distribuído pode ser necessário para
descrever tanto um V.C. como um Sistema.
Como exemplo, a Figura 3.5 mostra o esquema de um estuário, região onde os rios
deságuam no mar. Estuários são estudados em hidráulica ambiental, pela sua importância
como áreas de reprodução de muitas espécies.
água doce
Rio
0 maré alta
maré baixa
zona de mistura
 = f (x,y,z,t) água salgada
1
Figura 3.5: Corte esquemático de estuário, com variação da massa específica da água.
Os estuários são regiões em que a salinidade da água depende da proporção da

mistura entre as águas do mar e do rio. A mistura é influenciada pela vazão do rio, pela
topografia, pelos ventos e principalmente pelas marés. A massa específica da água varia no
espaço e no tempo, tornando necessário o uso de “modelos distribuídos”.
Uma aproximação razoável para a massa M nos modelos distribuídos pode ser obtida
se dividirmos o volume total em vários pequenos volumes Vol . Esses volumes são tão
pequenos que podemos considerar a massa específica constante em seu interior. A massa
M fica então:
n
M   Vol (Quanto menor o Vol, melhor a aproximação) 3.21
i 1
A medida que aumenta o limite n do somatório, diminui o tamanho dos volumes Vol.
Quanto menor o volume, mais preciso fica o cálculo da massa, já que consideramos 
constante em cada volume e ele pode variar continuamente. Esse processo, no limite, leva-
nos à massa como resultado de uma integração. A integral corresponde ao somatório das
massa de infinitos volumes diferenciais dVol:
n
M  lim (  i Voli )    dVol 3.22
VC
n  i 1
O limite de integração denotado genericamente como “VC” indica que o somatório

deve incluir todo o volume de controle. A especificação dos limites dependerá, em cada
caso, do elemento diferencial de integração “dVol” adotado.
_____ A Equação Final
Colocando todas as nossas considerações em conjunto, chega-se à equação geral do

balanço de massa, em sua forma integral.
BALANÇO GLOBAL DE MASSA
  
t  VC
 dVol   SC
 V . dA  0
3.23
Surpreendente e Elegante!
_____ Equação da Continuidade
Muitos problemas importantes ocorrem em condições de escoamento incompressível.

Nesse caso a massa específica pode sair da integral. O balanço de massa assume a forma
de balanço de volumes, conforme a equação 3.24.
  

t  VC
dVol    SC
V . dA  0
  
t  VC
dVol   SC
V . dA  0 3.24
O balanço integral de volumes da equação 3.24 é conhecido como “Equação da

Continuidade”. Quando a continuidade é aplicável, ou seja, escoamento com massa
específica constante, tanto o balanço de massa como o de volumes dão o mesmo resultado.
_____ Balanço de Substâncias Transportadas
A massa das substâncias transportadas pelo fluido em escoamento pode ser

expressa na forma de balanço integral para um V.C. com o mesmo desenvolvimento usado
para deduzir a equação 3.23. O balanço integral é dado por:
𝑡
∫𝑉𝐶  𝑉𝑜𝑙 + ∫𝑆𝐶  ⃗
𝑉 ⃗⃗⃗⃗⃗ = 3.25
Em que temos: N = mA , a massa da substância A dissolvida na água;

 = dmA / dm, a concentração em massa, kg de A por kg de fluido;
CA =   , a concentração de A em kg por m3 de fluido.
Observe que a equação 3.23 está contida na equação 3.25, para o caso em que a
grandeza é a massa do próprio fluido em escoamento. Tem-se neste caso N = m e  = 1.
3.4 Discussão Sobre a Taxa de Variação da Grandeza no V.C.
Existe a possibilidade de expressar o balanço de massas das equações 3.23 ou 3.25

de outra forma, sem mudar o significado físico, conforme apresentado na equação 3.26. A
diferença reside no termo que contém a derivada no tempo.
𝑁 
| = ∫𝑉𝐶 𝑉𝑜𝑙 + ∫𝑆𝐶  ⃗
𝑉 ⃗⃗⃗⃗⃗ 3.26
𝑡 𝑆 𝑡
A equação 3.25 aponta para inicialmente calcular a quantidade da grandeza N no

V.C. e depois calcular sua taxa de variação no tempo pela derivada.
Já a equação 3.26 indica que primeiro deve-se calcular a taxa de variação local da
grandeza N em cada elemento diferencial de volume, para depois integrar (somar) as
contribuições de todo o V.C.
Fisicamente a ordem de cálculo diferente das duas equações não tem impacto no
resultado final. O exemplo ilustrativo a seguir demonstra essa conclusão.
_____ Exemplo Ilustrativo: caixa de areia
Considere uma caixa de areia, conforme esquema da Figura 3.6. As caixas de areia
são utilizadas na entrada das estações de tratamento de esgoto para remover os sólidos
sedimentáveis, que podem prejudicar as etapas posteriores do tratamento.
comporta
VISTA SUPERIOR
Areia sedimentada
CORTE LONGITUDINAL
Figura 3.6: Esquema de uma caixa de areia com duas câmaras de sedimentação.
A caixa de areia tem duas câmaras que funcionam de forma alternada. Enquanto uma
câmara é usada, a outra deve ser limpa. O problema em questão é determinar o tempo
médio de utilização de cada câmara. Para isso é necessário saber a taxa de acúmulo de
areia na câmara de sedimentação.
O problema pode ser analisado com o balanço da grandeza extensiva N dada pela
massa de sólidos (areia) no esgoto. São definidos:
N = mS , a massa de sólidos sedimentáveis na água;
 = dmS / dm, a concentração de sólidos em kg por kg de fluido;
CS =   , a concentração de sólidos em kg por m3 de fluido.
Com a notação acima definida o balanço de sólidos na caixa de areia é dado pela
equação geral 3.26 da seguinte forma:
⃗ ⃗⃗⃗⃗⃗ =
∫ 𝐶𝑆 𝑉𝑜𝑙 + ∫ 𝐶𝑆 𝑉
𝑡 𝑉𝐶 𝑆𝐶
A resposta pode ser obtida experimentalmente pela avaliação do saldo de fluxos de

massa de sólidos, dado pela integral de superfície (segundo termo do segundo membro), ou
pela avaliação direta da taxa de acúmulo, dada pelo primeiro termo do segundo membro.
_____ Solução pelo saldo de fluxos de entrada e saída

Para o volume de controle da Figura 3.6, adotando-se o modelo concentrado para os
fluxos, pode-se escrever:
− + =
𝑡
−𝐶 +𝐶 =
𝑡
Sendo o regime permanente, a vazão que entra é igual à que sai.

= (𝐶 − 𝐶 )
𝑡
Portanto, é possível conhecer a taxa de variação da massa de areia medindo-se as
concentrações de areia na entrada e saída e a vazão. A medição deve ser mantida por um
período de tempo significativo, em intervalos suficientemente pequenos. No caso em
questão, deve ser previsto uma campanha de medição de 24h, no mínimo, com amostras
obtidas a cada hora, durante um dia típico. Com os dados experimentais tem-se a variação
diária:
𝑡
𝑆 = ∫ 𝑆 = ∫ ( 𝐶𝑠 − 𝐶 ) 𝑡
0
Aproximando numericamente a integral por meio de um somatório:
2
𝑆 ∑(𝐶𝑠 − 𝐶 ) 𝑡
𝑖
Com esse procedimento calculamos a variação diária de massa e, conhecendo a
quantidade máxima de areia, o período de tempo máximo para encher a caixa de areia.
Essa abordagem fornece a resposta com rapidez, mas necessita de um levantamento
experimental de dados custoso.
_____ Solução pela determinação direta do volume de areia acumulado
Para aplicar a técnica de medição direta devem ser realizados dois levantamentos
batimétricos, separados por um período de tempo significativo, que pode ser de dois a três
dias no caso em questão.
O levantamento pode adotar, por exemplo, a técnica de medição de profundidades
em pontos organizados numa malha, conforme o esquema da Figura 3.7.
j
n
Zi,j
A i,j
1
i
0
0 1 2 ... m
Zi,j
H
hi,j
Figura 3.7: Esquema do levantamento batimétrico numa malha de m x n pontos.

Observa-se na figura que é medida a profundidade Zi,j em cada ponto de uma malha
de m x n pontos. A cada profundidade medida corresponde uma altura de areia depositada
hi,j, que pode ser calculada conhecendo-se a profundidade total da caixa H, sendo que
hi,j = H – Zi,j.
Cada ponto representa uma área Ai,j delimitada pela linha de fronteira traçada na
metade da distância até cada um dos pontos vizinhos. Na malha retangular as áreas de
influência de todos os pontos internos são iguais. Nos pontos do contorno sólido a área é
metade e nos cantos, um quarto das áreas interiores.
Método 1: Levantar os volumes sedimentados totais no início e final

Com o levantamento batimétrico é possível aproximar o volume da areia no início e no
final de um período de tempo representativo da seguinte forma:
𝑛 𝑛
𝑆 = ∑∑ 𝑉𝑜𝑙𝑖 = ∑∑ 𝑖 𝑖
𝑖 𝑖
Em que o índice k representa o tempo em que foi realizado o levantamento e a

massa específica da areia sedimentada.
Com as massas de areia no início (índice k) e final do intervalo de tempo (índice k+1),
podemos calcular a taxa de variação da massa acumulada por:
𝑆 𝑆 𝑆 − 𝑆
=
𝑡 𝑡 𝑡
Esta abordagem de cálculo seguiu a sequência estabelecida na equação 3.24:
primeiro são calculadas as massas totais pela integração no V.C. e depois é calculada a
derivada.
Método 2: Levantar a variação de volume em cada ponto e totalizar

Outra forma de abordar o problema consiste em calcular inicialmente a variação de
volume de areia acumulado em cada área Ai,j e posteriormente somar a variação em todas
as áreas. Esta sequência de cálculo é descrita pela equação a seguir:
𝑛 𝑛
𝑆 = 𝑆 − 𝑆 = ∑∑ 𝑉𝑜𝑙𝑖 = ∑∑ ( 𝑖 − 𝑖 ) 𝑖
𝑖 𝑖
A abordagem do segundo método adota a sequência descrita pela equação 3.25:

primeiro é calculada a taxa de variação da grandeza no tempo (derivada), e depois as taxas
são somadas em todo o V.C. (integral) para fornecer a taxa de variação total no volume de
controle.
3.5 Exercícios
1 - Na operação de uma hidrelétrica ocorreram durante um dia as seguintes vazões médias:

Qrio = 1.500 m3/s, Qturbinas = 600 m3/s. Qual a variação de volume armazenado no
período de 24 horas?
2 - Uma hidrelétrica tem um lago alimentado por dois afluentes, e 6 turbinas. As vazões
afluentes observadas no dia foram de Q1 = 600 m3/s e Q2 = 120 m3/s. Cada turbina operou
com vazão de 100 m3/s, mas 2 turbinas ficaram ligadas apenas 6 horas. Qual a variação do
volume do lago?
3 - Um meio para determinar a vazão de rios consiste na injeção de substâncias traçadoras,

como sais ou corantes. Numa determinação de vazão em um córrego foram lançados 1L/s
de água com uma concentração de sal igual a 25mg/L. Numa seção a jusante, após a
completa mistura do traçador, retirou-se uma amostra da água, obtendo-se uma
concentração de 0,013mg/L. Qual a vazão do córrego, sabendo que numa seção a montante
do ponto de injeção foi medida uma concentração de 0,003mg/L?
4 - Em um período de cheias uma hidrelétrica operou com as seguintes vazões afluentes:

Rio A: Q1 = 1.300 m3/s; Rio B: Q2 = 700 m3/s. As vazões efluentes são fluxos turbinados
para produção de energia e fluxos vertidos, que escoam pelos vertedores. A usina possui 8
máquinas e 6 vertedores de superfície. A configuração de operação foi a seguinte: Vazão
turbinada - 4 máquinas com 100 m3/s durante 24 horas ; 2 máquinas com 180 m3/s durante
2 horas; Vazão vertida - 6 vertedores operando continuamente com 200 m 3/s cada um.
Calcule a variação do volume armazenado em 24 horas.
5 - A figura mostra um esquema de chaminé de equilíbrio. As chaminés de equilíbrio são

utilizadas na prática para atenuar as variações de pressão que podem ocorrer durante
transientes em tubulações de alimentação de bombas e turbinas. No instante considerado, a
velocidade no tubo de alimentação da turbina é 0,5m/s a montante da chaminé e de 2,5m/s
a jusante da mesma. Calcule a vazão fornecida e a velocidade instantânea de abaixamento
do nível d'água na chaminé.
Reservatório Chaminé
D = 3m
D = 1m Turbina
6 - Um tanque cilíndrico possui uma área de base igual a 1m 2. A água escoa por um orifício
de 50 mm de diâmetro, segundo a figura. A vazão que escoa em m3/s é dada por Q = 2,7Ah
a) Calcular o nível da água após 1 minuto

b) Quanto tempo leva para escoar 2 m3 do tanque?
5,0m c) Repetir o item (a) usando integração numérica,
com t = 20s.
d) Repetir o item (a) usando solução numérica, com
t = 5s.
e) Quanto tempo leva para escoar 1 m3 do tanque?
7 - Uma represa forma um reservatório de 5x107m3 de capacidade. O lago recebe a

contribuição de dois rios, com as seguintes vazões e concentrações médias de sedimentos:
Rio A: Q = 12 m3/s; Csedimentos = 10g/L
3
Rio B: Q = 3 m /s; Csedimentos = 18g/L
Sabendo que na saída a concentração de sedimentos é 2g/l, e que a massa específica dos
sedimentos é sed = 2,65 g/cm3, determinar qual o tempo de vida estimado para o
reservatório.
8 - O duto da figura tem seção transversal quadrada com 0,1m de lado e descarrega água
por quatro fendas de 0,01m por 1m localizadas em nas faces laterais de uma derivação.
Sabendo que o regime é permanente, que o duto da derivação é fechado em sua
extremidade inferior e com base nas velocidades nas faces dadas na figura, pede-se:
a) vazões nas faces 1 e 3; b) módulo e sentido da velocidade média na seção b.
Seção a b
0.1m
Va = 8 m/s Vb
Velocidades:
0.1m V1 = V2 = 4 - 2Z
V3 = V4 = 2 - Z
Z V3
1.0m
V2 V1
Y
0.01m X V4
CAPÍTULO 4
TRANSFORMAÇÕES DE ENERGIA NOS ESCOAMENTOS
4.1. Equação de Bernoulli
Uma das equações mais importantes da hidrodinâmica é a equação de Bernoulli, que

explica como variam a pressão e os termos de energia potencial e cinética de um fluido em
escoamento. A equação de Bernoulli é deduzida com aplicação do teorema que relaciona a
variação da energia mecânica ao trabalho realizado sobre um sistema.
Imagine um fluido ideal incompressível escoando em regime permanente entre as
seções 1 e 2 de um tubo de corrente, conforme a Figura 4.1. Vamos considerar o Volume de
Controle como o tubo de corrente entre as duas seções, e o Sistema como sendo a masssa
de fluido que no instante inicial da análise encontra-se dentro do volume de controle,
identificado pela região tracejada do VC.
V2
(a) tempo t
F2
V1 III
2
F1
I Z2
Z1 1
x 1
V2
(b) tempo t +  t F2
III II
V1
F1 Z2
x 2
Z1
Figura 4.1: Escoamento de um sistema num tubo em dois instantes de tempo
No instante inicial da análise o Sistema coincide com o VC. Com o passar do tempo a
água deixa o VC, como é mostrado na Figura 4.1 (b), em que a porção tracejada (Sistema)
não mais coincide com o VC.
Aplicaremos o princípio da conservação da energia mecânica ao Sistema:
= = − 4.1
Em que a energia do sistema, em qualquer tempo, é a soma da energia potencial e da

cinética.
O sistema é descrito pelas variáveis Eulerianas por meio da região do espaço
ocupada em cada instante. Assim, temos:
()= ( )+ () 4.2
( + )= ( + )+ ( + ) 4.3
O trabalho sobre o sistema é realizado apenas por forças de pressão, porque as

forças tangenciais são nulas, devido à hipótese de fluido sem viscosidade (ideal). Além
disso, será realizado trabalho apenas nas seções 1 e 2, porque o restante da fronteira do
sistema não se move na direção das forças de pressão.
Na seção 1 o resto do fluido (meio) exerce sobre o sistema uma força p 1A1. Se o
intervalo de tempo considerado for muito pequeno, as pressões podem ser consideradas
constantes, assim com as áreas. Portanto, a força é constante e, ao longo do intervalo de
tempo considerado, esta força desloca o sistema de x1. Como consequência, o sistema
avança através da seção 2, deslocando-se por uma distância x O restante do fluido opõe-
se ao deslocamento do sistema, com a força de pressão p 2A2 . Da mesma forma pode-se
considerar a força constante ao longo do deslocamento.
Das considerações acima, durante o intervalo de tempo considerado, as forças que o
meio exerce sobre o sistema realizam um trabalho dado por:
𝑊 = 𝑝 𝑥 − 𝑝2 2 𝑥2 4.4
Aplicando os resultados 4.2 a 4.4 em 4.1:

− 2 2 2 = ( + )+ ( + )− ( )− () 4.5
Pensando na variação de energia do sistema, vemos que a região III compõe a parte
do sistema cuja energia não variou entre t e t+t. Toda essa parte contém, em cada ponto, a
mesma quantidade de fluido e à mesma velocidade no início e no fim do período. Isto quer
dizer que a soma das energias cinética e potencial das unidades de massa dessa região do
sistema não varia. Portanto, a variação ocorre porque a parte I do sistema desapareceu,
dando origem à parte II.
− 2 2 2 = ( + )− () 4.6
Se o escoamento é incompressível verificarmos, pelo princípio de conservação das

massas, que a massa da parte I é a mesma da parte II.
 1
 Ec1  mV12 Ep1  mg z1
2
m1  m2  m 
 Ec  1 mV 2 Ep2  mg z2
 2 2 2
Reunindo todas as variações na equação 4.6, temos:

1
p1 A1x1  p2 A2x2  m( V22  V12 )  mg( z2  z1 ) 4.7
2
Os termos da equação 4.7 representam energia total. Dividindo a equação pela
quantidade de massa envolvida na variação da energia, temos, lembrando que “Ax” é o
volume de fluido que entrou e saiu do VC:
p1  p2 1 2
 (V2  V12 )  g ( z2  z1 ) 4.8
 2
Em que os termos representam energia por unidade de massa (Nm/kg). Rearranjando
os termos obtemos a Equação de Bernoulli:
p1 V12 p2 V22
  g z1    g z2 4.9
 2  2
Daniel Bernoulli (1700-1782)
A equação de Bernoulli representa com pressões e velocidade em seções definidas

do espaço, portanto variáveis Eulerianas, a variação da energia de um Sistema. As
limitações da análise incluem um intervalo de tempo tendendo a zero, ou seja, uma condição
instantânea. Mas, como o escoamento é em Regime Permanente, a relação fica válida para
qualquer tempo.
Como as seções 1 e 2 podem ser quaisquer, é usual também expressar a equação de
Bernoulli na forma:
Hipóteses Utilizadas
2
p V Escoamento Incompressível
  g z  Cte

Regime Permanente
2 Atrito desprezado
Em um tubo de corrente
Equação de Bernoulli (1738)
4.2. Conservação da Energia nos Escoamentos
Para explorar as consequências da equação de Bernoulli é útil analisar o caso dos

fluidos em escoamento em relação aos corpos rígidos, tirando partido das semelhanças e
discutindo as diferenças entre os dois casos.
A Figura 4.2 mostra o comportamento de uma esfera rolando sem atrito sobre uma
superfície ao se deparar com aumento ou diminuição da cota.
V2
z 2
m V1 t
z1 V2 < V 1
t 0
m V1
z1 V2 > V 1
t 0 V2
z 2
t
Figura 4.2: comportamento de sistema sólido numa rampa
Esse exemplo é bastante familiar. Já nos acostumamos, da observação cotidiana dos

corpos isolados (sistemas), com o fato de que uma rampa transforma energia cinética em
potencial. A conservação da energia, na ausência de atrito, permite escrever:
V12 V22
m  m g z1  m  m g z2 4.10
2 2
_____ Conclusão: Em um corpo isolado a velocidade é livre para variar.
A Figura 4.3 mostra três situações possíveis para um escoamento forçado no interior
de um tubo. A ilustração fala por si só. Nos escoamentos forçados a energia cinética deixa
de ser controlada pela cota, como aconteceria com um corpo isolado.
V2 D2 = D1
Escoamento Forçado
V2 = V1
2
V1
V2
1 D2 < D1
V2 > V1
2
V1
1 V2 D2 > D1
V2 < V1
2
V1
Área determina a Velocidade
1
Figura 4.3: Comportamentos possíveis de um escoamento quando aumenta a cota da tubulação.
A Figura 4.3 deixa claro que nos escoamentos forçados não há como compensar a
variação da cota com a energia cinética. A velocidade é controlada apenas pelo diâmetro da
tubulação.
_____ Conclusão :
Nos escoamentos forçados, a Velocidade não é determinada
pelas variações da cota.
_____ Aperfeiçoando a Analogia Trem-Escoamento com Energia Potencial Elástica
Já vimos que não é possível comparar o escoamento com uma massa isolada, sob
pena de comprometer o entendimento físico do problema. Entretanto, a nossa analogia do
escoamento com um trem, constituído de muitos vagões interligados, é bem mais
aproximada da situação que realmente ocorre. Observe a situação da Figura 4.4.
Na Figura 4.4 a velocidade dos vagões não irá variar quando passarem pela
elevação, pois a velocidade de cada vagão é a mesma do trem, que é constante O que
ocorre então, se a energia Potencial aumentou e a Cinética não pode diminuir para
compensar?
Corpo Isolado: V 1= V 3
V 2 < V1 V2
V
V1 3
V2
Trem: V 1 = V 2 = V 3
V1 V
3
Figura 4.4: Comportamento de corpo isolado e de um trem frente a uma elevação de cota.
A resposta a essa pergunta pode ser dada quando pensamos nas barras de
união entre os vagões, introduzindo mais um componente na analogia, capaz de armazenar
energia. Pensando que nos vínculos como molas, eles se tornam capazes de armazenar
energia potencial sob a forma de deformação elástica.
O potencial elástico das molas de união é análogo, no caso real, à pressão existente
nos escoamentos.
_____ Energia Potencial Elástica
Imagine (Fig. 4.5) uma mola com constante elástica K, sendo comprimida por uma
força F. O trabalho realizado pela força ficará armazenado na mola sob forma de energia
potencial.
F
FK x
xi=0
F dW F dxK xdx
xf x x2
xf
W  
xi
dW  
0
K x dx  K
2
Figura 4.5: mola armazena energia
No caso do trem subindo a colina, a energia cinética permanece constante, mas a

força nos vínculos diminui à medida que o vagão sobe a colina.
Quando um sistema massa-mola é desacelerado a força gerada pela variação da
energia cinética fica armazenada como energia potencial elástica na mola. Esse efeito é
ilustrado pela Figura 4.6.
V
(a) li
V=0
(b)
lf
Figura 4.6: Desaceleração de sistema massa-mola comprime a mola.
No caso da Figura 4.6(a), com a massa deslocando-se livremente, a mola está

totalmente distendida, e toda a energia está na forma cinética.
1
Ec  mV 2 E pot  0 4.11
2
Na Figura 4.6(b) a massa foi desacelerada e a velocidade é nula. Supondo atrito nulo,
toda a energia cinética disponível foi usada para comprimir a mola:
K ( li  lf ) 2 1
Ec  0 E pot   mV 2
4.12
2 2
O caso da Figura 4.6(b) corresponde, na analogia trem-escoamento, ao que ocorre

nos pontos de estagnação, em que toda a energia cinética da água é convertida em
pressão.
Com esses exemplos percebemos que o nosso modelo de analogia entre trem e
escoamento precisa incluir a capacidade de armazenar energia sob a forma de deformação
elástica. A variável que corresponde à mola é a pressão.
4.3. Energias e Cargas na Equação de Bernoulli
Para obter a equação 4.11 dividimos a energia mecânica expressa na equação 4.10
pela quantidade de massa envolvida (m = A x) . Uma forma bastante conveniente da
equação surge ao dividirmos a equação 4.10 pelo peso da massa envolvida, ou seja, W =
mg = g A x . A equação 4.10 fica:
p1 V12 p2 V22
  z1    z2 4.13
 2g  2g
Lembrando que o peso específico  = g.
Os termos da equação 4.13 representam energia por unidade de peso do fluido. As
unidades no sistema SI são (Nm/N) ou simplesmente (m). Pelas unidades percebe-se que
os termos da equação 4.13 são Cargas. Esta nomenclatura vale-se da correspondência
entre pressões e colunas de fluido que exercem a mesma pressão. O termo carga surgiu
para definir a altura da coluna de água, conforme o esquema da Figura 4.7.
dA
Se p1 = p2, dizemos que H é a
coluna equivalente
de fluido altura de coluna equivalente, ou CARGA:
H p1   g H   H
tubo p1
H 
1 2

Figura 4.7: Relação entre pressão e altura equivalente de coluna de fluido ou Carga.
Com base nessa analogia, todos os termos de energia por peso da equação foram
chamados de cargas. Temos:
𝑝
Carga de Pressão: 𝑝 =
𝑉2
Carga Cinética: =
2
Carga Potencial : = 𝑧
Lembrando que as cargas de pressão resultam do trabalho das forças externas sobre
o sistema e as demais compõem a energia mecânica do sistema (por unidade de peso).
Com o uso das cargas a equação de Bernoulli pode ser escrita da forma:
𝑝
𝑇 = + 𝑧 + 𝑉 (m) 4.14
Os dois primeiros termos compõem a Carga Piezométrica, porque correspondem à

cota do nível que a água atinge num piezômetro instalado na seção considerada. Veja o
esquema da Figura 4.8.
Fluido Manométrico
Carga Piezométrica:
p1
Fluido de Trabalho 𝑝
 𝑝𝑧 = + 𝑧
1
Z1
Figura 4.8: Piezômetro – cota do nível do fluido coincide com a carga piezométrica.
4.4. Aplicação a Medições de Vazão e Velocidade
4.4.1. Medidor Venturi

Os medidores tipo Venturi utilizam um estreitamento da seção para provocar aumento
da carga cinética. Manômetros medem a diferença na carga de pressão. Um esquema de
medidor de vazão tipo Venturi é apresentado na Figura 4.9.
Figura 4.9: Medidor tipo Venturi.
Tomando um tubo de corrente entre as seções 1 e 2 e aplicando Bernoulli, temos:

p1 V12 p2 V2
 z1    z2  2
 2g  2g
2 2
p1 p2 V2 V
   1
  2g 2g
Pela continuidade, temos 𝑉 = 𝑉2 2  𝑉2 = 𝑉

2
p1  p2 V 2
A 2

 1
 1
 1
 2g A2
2

Resolvendo em função de V1 temos:
𝑝 − 𝑝2
𝑉 = 2
√
( 2 − )
2
Com a velocidade média pode-se calcular a vazão:

= 𝑉
4.4.2. Tubo de Pitot

O medidor tipo Pitot é formado por um tubo com abertura perpendicular às linhas de
corrente, conforme visto na seção 2 do esquema da Figura 4.10. O fluido é desacelerado no
trajeto entre 1 e 2. No ponto 2, sobre a linha de corrente central do escoamento, a
velocidade é nula, pois o líquido está estático no interior do tubo. O ponto 2 é chamado de
Ponto de Estagnação.
1 2 L.C.
V
a 1 2
L
(a) (b)
Figura 4.10: Medidor tipo Pitot (a) e detalhe do ponto de estagnação (b).
As tomadas de pressão perpendiculares às linhas de fluxo, como a do ponto 2, são

chamadas na prática de “tomadas dinâmicas” de pressão, porque são afetadas pela
velocidade do escoamento. Em contraposição, tomadas de pressão paralelas às linhas de
fluxo, como as do ponto 1, são chamadas de “tomadas estáticas” de pressão, porque a
leitura indicada não é afetada pela velocidade do fluido.
Entre os pontos 1 e 2, ao longo da linha de corrente central, pode ser aplicada a
equação de Bernoulli, supondo que a desaceleração da água ocorre sem perda apreciável
de energia.
p1 V12 p2 V22
 z1    z2 
 2g  2g
Temos z1 = z2 e V2 = 0, pois o ponto 2 é ponto de estagnação.

V12 p2 p1
 
2g  
𝑝2 𝑝
𝑉 = √
O manômetro diferencial de tubo em U permite calcular a diferença de pressões.

𝑝 +  + 𝑀 −  −  = 𝑝2
𝑝2 − 𝑝 = ( 𝑀 −  )
sendo “M “ o peso específico do fluido manométrico e “ “ o peso específico do fluido de
trabalho.
4.4.3. Orifícios de pequenas dimensões

Quando a carga sobre o orifício é grande em relação a seu diâmetro, a velocidade de
saída do fluido é aproximadamente constante e o orifício é chamado de pequenas
dimensões. A Figura 4.11 mostra o esquema de um grande reservatório descarregando
água por um orifício, formando um jato livre. Em um grande reservatório o nível da água
varia muito lentamente.
Veia contraída
V1 Seção
Contraída Distribuição
h de pressões
no orifício
V2
d
Figura 4.11: Orifício de parede delgada em reservatório de grandes dimensões e detalhe da veia
contraída.
A seção contraída é a primeira seção a partir do orifício em que a pressão do jato é

conhecida. Após a seção contraída o jato é livre (p = p atm ) e antes dela a pressão segue
uma distribuição desconhecida, conforme indica o detalhe da Figura 4.11.
Aplicamos a equação de Bernoulli entre um ponto qualquer do reservatório e a seção
contraída, com o referencial de cota no eixo do orifício.
p1 V12 p2 V22
 z1    z2 
 2g  2g
Qualquer que seja a posição do ponto 1 no reservatório, desde que suficientemente

afastado do orifício, temos:
V12
 0 : a carga cinética é desprezível;
2g
p1
 z1  h : carga média sobre o orifício (não depende da posição do ponto 1).

Como na seção 2 o jato é livre, temos :
V22
h  0  0  𝑉2 = √
2g
Entretanto, a velocidade real é menor que a calculada, devido às perdas no processo
de aceleração. Essas perdas são levadas em conta por um fator de correção experimental
da velocidade, chamado de Coeficiente de Velocidade (Cv).
𝑉2 𝑅 𝑎𝑙
𝐶𝑉 =
𝑉2 𝑇 𝑖 𝑎
A área do jato também é menor, devido à curvatura das linhas de corrente ao passar
pelo orifício. O fenômeno da contração do jato é levado em conta por outro coeficiente
experimental, chamado de Coeficiente de Contração (Cc).
𝑆 𝑜 𝐶𝑜𝑛𝑡 𝑎 𝑎
𝐶𝐶 =
𝑂 𝑖 𝑖𝑜
Para orifícios circulares de parede delgada Cv varia entre 0,95 e 0,99 e o Cc = 0,62.
A vazão através do orifício com área Ao fica então:
= 𝑉𝑅 𝑎𝑙 𝑎𝑡𝑜 = 𝐶𝑉 √ 𝐶𝐶 𝑂
É usual combinar os dois coeficientes experimentais num único coeficiente, chamado

de Coeficiente de Descarga (CD).
𝐶 = 𝐶𝐶 𝐶𝑉
E a vazão fica:
= 𝐶 √
O coeficiente de descarga típico de orifícios circulares de parede delgada é 0,61.
4.4.4. Bocais em condutos forçados
Quando um tubo descarrega por meio de um bocal que causa estreitamento da seção
e aumento da velocidade, a vazão descarregada pelo tubo pode ser facilmente determinada
com auxílio de um manômetro. Esta aplicação é ilustrada pela Figura 4.12 que mostra um
bocal na extremidade de um tubo dotado de um piezômetro.
Piezômetro
a
h
Bocal L.C. retilíneas
2 z
1 e paralelas
Jato livre 1
Figura 4.12: Bocal descarregando em jato livre e detalhe da variação de pressão no interior do tubo.
A vazão é determinada aplicando-se a equação de Bernoulli entre os pontos 1 e 2

situados ao longo da linha de corrente que passa pelo eixo do tubo.
A carga de pressão no ponto 1, situado na linha média do tubo pode ser calculada
diretamente com a leitura no piezômetro, tendo em vista que no interior do tubo as linhas de
corrente são retilíneas e paralelas:
𝑝
𝑝 = 𝑝𝑎𝑡 +  =
A equação de Bernoulli fica:

V12 V2 V22 V2
h  z1   0  z2  2   1  h
2g 2g 2g 2g
A equação da continuidade fornece a relação entre as velocidades necessária para
resolver o problema. Pela continuidade, temos:
𝑉 = 𝑉2 2  𝑉2 = 𝑉
2
Resolvendo para V1, a velocidade no tubo:
𝑉 = 2
√
( 2 − )
2
4.5. Exercícios
4.5.1) Um fluido de  = 10.000 N/m3 escoa por um tubo horizontal, com uma redução de
diâmetro de 150mm para 50mm. Sabendo que a pressão na seção 1, antes da redução, é
de 500.000 N/m2, e a velocidade é 2m/s, calcular a pressão na seção 2, após a redução.
4.5.2) Desprezando as perdas no Venturi da figura, calcule a vazão de água transportada.
ar
0.75m
água
D1 = 0,3m
D2 = 0,15m 2
4.5.3) A figura mostra dois Tubos Venturi verticais, transportando água ( = 1000 kg/m3),
instalados lado a lado e interligados por dois circuitos manométricos, com fluido
manométrico de  = 2000 kg/m3 . Sabendo que no ponto 1 a velocidade é 1m/s e a pressão
é 10.000 Pa, pede-se:
a) Calcular a pressão no ponto 2 b) Calcular a pressão no ponto 3
c) Calcular a pressão no ponto 4 d) Calcular a vazão no venturi da esquerda.
3 1 z=1m
2,0m
4 2 z=0,5m
0,25m
Q 50mm
100mm
4.5.4) A figura mostra um reservatório de grandes dimensões, com ar pressurizado na parte

superior, conforme demonstrado pelo manômetro. O reservatório possui um orifício de
0,01m de diâmetro em sua parede lateral, descarregando um jato livre na atmosfera, no qual
foi posicionado um tubo de Pitot, em posição que não obstrui o jato.
Pede-se determinar a vazão que sai pelo orifício e a leitura h, sabendo que o coeficiente de
contração do orifício é 0,66.
Ar
h
0,15m
0,40m
Jato Livre
Hg Água
4.5.5) O tubo de 100mm de diâmetro descarrega água em jato livre através do bocal de
50mm. As perdas de carga entre a saída e o piezômetro são h = 10V2/2g sendo V a
velocidade no tubo. A leitura do piezômetro é h= 2,5m. Calcular a vazão.
h
100mm
2m
50mm
4.5.6) Água de um reservatório de grandes dimensões escoa em regime permanente na

tubulação da figura, que termina num bocal de 75mm descarregando em jato livre. O fluido
no manômetro é mercúrio (d=13,6). Pede-se determinar, desprezando as perdas: a) a vazão;
b) a pressão p indicada pelo manômetro no tubo de 100mm, e c) a carga H.
4.5.7) Água de um reservatório de grandes dimensões escoa em regime permanente na

tubulação da figura, que termina num bocal de 100mm descarregando em jato livre. O fluido
no manômetro é mercúrio (d=13,6). Pede-se determinar, desprezando as perdas: a) a vazão;
b) a pressão p indicada pelo manômetro no tubo de 100mm, e c) a carga H.
4.5.8) Um reservatório de grandes dimensões contém água ( = 9.800 N/m3 ) escoando por
um orifício circular de parede delgada na parede inclinada a 60, sendo que o jato sobe até
um nível 2,0 m acima do orifício. Sabe-se que a área do orifício é 5 cm2, e que seus
coeficientes são CC = 0,60 e CV = 0,9. Calcule a) a carga H no orifício b) a vazão escoada.
H= ?
2m
4.5.9) Na tubulação da figura a água escoa com velocidade de 2,4 m/s no ponto A. Aonde o
nível da água chegará no piezômetro C?
Ex. 4.5.9
Ex. 4.5.10
4.5.10) Se cada manômetro indicar a mesma leitura para uma vazão de 28 l/s, qual deverá
ser o diâmetro da constrição na figura deste problema?
4.5.11) Para uma vazão de ar de 2 m3/s ( = 12,0 N/m3) qual deverá ser a maior seção A2
necessária para que a água se eleve até a abertura do piezômetro? Despreze efeitos de
compressibilidade.
4.5.12) Bombeia-se ar através de um tanque conforme indicado na figura. Desprezando

efeitos de compressibilidade, calcule a velocidade do ar no tubo de 100 mm. A pressão
atmosférica é de 91 kN/m2 e o peso específico do ar é igual a 11 N/m3.
Ex. 4.5.12
Ex. 4.5.13
4.5.13) A água escoa num tubo vertical de 50mm conforme a figura, caindo em jato livre
sobre um disco com 0,30m de diâmetro. No centro do disco está um manômetro diferencial.
O escoamento é axissimétrico, com a água deixando o disco horizontalmente, num jato com
1mm de altura. Calcule a vazão e a deflexão no manômetro.
CAPÍTULO 5: TRANSFORMAÇÃO DE REYNOLDS

Relação Sistema x Volume de Controle
5.1 Introdução
Vamos iniciar com o balanço global de massa (equação 3.22), deduzido no item 3.3.
  
t VC  dVol   SC  V . dA  0 (eq. 3.22)
Alguém com talento para generalizações pode perceber, ao considerar a equação

5.22, que o segundo membro representa uma declaração válida para qualquer Sistema, em
relação à massa, da seguinte forma:
dM
 0 5.1
d t SISTEMA
O primeiro membro surgiu ao aplicar o conceito de conservação de massa a um VC,
ou seja, representa a taxa de variação da massa do V.C. ao longo do tempo. Mas, o que é a
taxa de variação no tempo da massa de um V.C.? Aplicando diretamente a definição do
cálculo obtemos equação 5.2.
𝑀 𝑀 (𝑡 𝑡) 𝑀 (𝑡)
| = 𝑡 0 5.2
𝑡 𝑉𝐶 𝑡
Onde M indica a massa do Volume de Controle e os índices indicam o tempo em que ela é
avaliada.
Pergunta : Podemos demonstrar que a derivada da equação 5.2 é equivalente
às duas integrais do balanço global da equação 3.23?
Se provarmos que a derivada da equação 5.2 é igual às integrais do balanço da

equação 3.22, mostraremos que, na verdade, a equação do balanço relaciona uma
propriedade da massa válida em um sistema à descrição dessa mesma propriedade com
variáveis medidas de forma Euleriana, ou seja, em pontos definidos do espaço. Tentaremos
deixar claro este ponto no próximo item.
5.2 Do Sistema ao Volume de Controle
A relação entre Sistema e Volume de Controle para uma dada lei física é chamada
também de Transformação de Reynolds. É uma ferramenta muito útil, porque todas as leis
da física clássica se aplicam a quantidades de massa definidas. Na Mecânica dos Fluidos,
como já pudemos perceber, é importante que tenhamos leis válidas para uma região
definida do espaço, o Volume de Controle. Logicamente, tudo que for válido para a Massa,
também será válido para qualquer outra grandeza extensiva N. Nessa generalização para
outras grandezas extensivas nosso esforço será amplamente recompensado.
Imagine um escoamento, visualizado esquematicamente na Figura 5.1 por meio de
suas linhas de corrente. Existe um Volume de Controle qualquer, que contém inicialmente
em seu interior uma massa que constitui o Sistema sob análise (Fig. 5.1(a)). Podemos
pensar que o volume foi demarcado por uma linha de tempo, acompanhando a forma do
sistema no instante inicial.
Decorrido um intervalo de tempo t, o Sistema, carregando consigo sua massa e suas
grandezas extensivas N, terá se deslocado devido ao escoamento, enquanto que o V.C.
permanece fixo no espaço, conforme ilustrado na Figura 5.1(b). Lembre-se que a Linha de
Tempo pode deformar-se ao acompanhar o escoamento, mas sempre demarcará a mesma
quantidade de fluido inicial (com a mesma massa e mesmas quantidades das grandezas
extensivas N - ver propriedades da linha de tempo).
(a) tempo t (b) tempo t + t

Sistema deslocou-se
1 2 3
V.C. = Sistema V.C. permanece fixo

Figura 5.1: Volume de Controle e Sistema em dois instantes consecutivos num escoamento.
A situação do sistema em relação ao VC permite definir 3 regiões do espaço,

conforme assinalado na Figura 5.1(b). A massa do sistema no instante t corresponde à
massa contida nas regiões 1 e 2. No instante seguinte a massa do sistema ocupa as regiões
2 e 3. A mesma consideração vale para as quantidades de grandeza N que a massa possui:
[ ] = ( + 2) 5.3
[ ] = ( + 2) 5.3-a
[ ]  = ( 2 + )  5.4
[ ]  = ( 2 + )  5.4-a
Com esse expediente, as equações 5.3 e 5.4 relacionam a massa do sistema por
meio de medições nas regiões do espaço correspondentes, indicadas pelos índices1, 2 e 3.
Observe que as quantidades de massa e da grandeza N devem ser descritas por variáveis
Eulerianas, embora a quantidade total das grandezas seja de um sistema.
A taxa de variação no tempo da grandeza N (ou da massa M) do sistema fica então:
dN N S t t  N S t ( N 2  N 3 ) t  t  ( N 1  N 2 ) t
 lim  lim 5.5
dt S t 0 t t 0 t
Lembrando que N é substituído por M no caso do balanço de massa.
Observe que parte da massa do sistema que antes estava no V.C. foi afastada pelo
fluxo que entra no V.C., levando consigo suas grandezas extensivas. Assim, podemos
associar a grandeza N (ou a massa) da região 1 da Figura 5.1 ao fluxo de entrada da
grandeza N ( ou da massa):
N1  FN ,E t 5.6
Da mesma forma, ao acompanhar o escoamento, o sistema teve parte de sua massa

atravessando a área de saída, identificada pelo volume 3 da Figura 5.1. Portanto, o fluxo de
saída da grandeza N pode ser usado para calcular a quantidade N3.
N 3  FN ,S t 5.7
Esses resultados serão usados para simplificar a equação 5.5:
dN ( N 2 )t t  ( N 2 )t  N 3  N1
 lim 5.8
dt S t 0 t
dN ( N 2 )t t  ( N 2 )t  ( FN ,S  FN ,E )t
 lim 5.9
dt S
t 0 t
É fácil perceber que os fluxos da equação 5.9 não dependem do limite considerado,
porque são fluxos médios no intervalo de tempo. Assim, a derivada da massa no sistema
fica com dois termos, mas apenas um varia no limite para t  0:
dN ( N 2 )t t  ( N 2 )t
 lim  FN ,S  FN , E 5.10
dt S
t 0 t
O próximo passo é verificar que quanto menor o tempo t , menor é o volume que vai
entrar e sair do V.C. devido ao escoamento. Isto implica que no limite o volume 2 tende para
o próprio Volume de Controle, de forma que podemos escrever a Lei Física da Conservação
da Massa da seguinte forma:
𝑁 𝑁
| = | + 𝑁𝑆 − 𝑁𝐸 5.11
𝑡 𝑆 𝑡 𝑉𝐶
Usando as variáveis de interesse nos termos do segundo membro da equação 5.11

obtém-se, para o caso da Massa:
𝑀
| = = ∫𝑉𝐶 ⃗
𝑉𝑜𝑙 + ∫𝑆𝐶 𝑉 ⃗⃗⃗⃗⃗ 5.12
𝑡 𝑆 𝑡
O esquema da Figura 5.2 realça os aspectos interessantes da equação 5.12

Lei Física Mesma Lei descrita com medições

válida para Sistema no Volume de Controle
Variáveis Lagrangeanas Variáveis Eulerianas
Figura 5.2: Transformação de Reynolds relaciona grandezas Lagrangeanas a grandezas Eulerianas
5.3 Balanço Global de Grandezas Extensivas

Vimos que tudo que é válido para a massa, também é válido para as grandezas
extensivas que a massa transporta. Assim, a relação sistema - volume de controle para uma
grandeza extensiva N qualquer pode ser escrita diretamente, a partir das equações 5.11 e
5.12 como:
𝑁
𝑡 𝑆
| =
𝑡
∫𝑉𝐶  𝑉𝑜𝑙 + ∫𝑆𝐶  ⃗
𝑉 ⃗⃗⃗⃗⃗ 5.13
Lembrando que a quantidade específica  da grandeza N é dada por: dN =  dm.

A equação 5.13 é chamada também de “balanço global de grandezas extensivas”.
Observe que não podemos considerar nula a taxa de variação da grandeza no sistema (lado
esquerdo - Lagrange), pois isso é válido somente para a massa, e não para grandezas
extensivas em geral. Isso será visto com detalhe nos itens seguintes. A Figura 5.3 realça os
aspectos de mudança de pontos de vista envolvidos na Transformação de Reynolds.
Lei Física: Mesma Lei descrita com medições
vale para Sistema no Volume de Controle
Variáveis Lagrangeanas Variáveis Eulerianas

Figura 5.3: Transformação de Reynolds para uma grandeza extensiva N qualquer.
CAPÍTULO 6: BALANÇO GLOBAL DE ENERGIA

A energia é uma propriedade das substâncias. Contrariamente ao trabalho ou
calor, o sistema possui uma quantidade determinada de energia, que é grandeza extensiva.
Sabemos da termodinâmica como representar a variação da energia de um sistema durante
um processo qualquer. Neste item vamos desenvolver o mesmo conceito em uma equação
a
válida para um Volume de Controle, aplicando a 1 . Lei a um sistema e usando a relação
Sistema - V.C.
A grandeza em questão é a energia E, de forma que seu valor específico por unidade
de massa (joules por kilograma, no SI) é :
dN V2
   e   gz  u
dm 2 interna
potencial
cinética 6.1
O primeiro princípio da termodinâmica é a lei da conservação da energia, deduzida a

partir de observação, e aplica-se a todos fenômenos físicos. Para um processo qualquer, a
primeira lei diz que a variação da energia é igual à diferença entre o calor fornecido ao
sistema e o trabalho realizado sobre o sistema durante qualquer variação de estado.
Q-
Q+ Sistema
 Q  W   E W+
W-
Processo
Convenção de Sinal
a.
Portanto, unindo a 1 lei da termodinâmica, que descreve a variação da energia num
sistema, ao resultado previsto pela relação Sistema - V.C., ficamos com a seguinte equação:
dE Q W   

t  t 
t VC e  dVol  SC e  V  dA
6.2
dt S
que representa o balanço de energia num V.C.

Vamos agora colocar a equação da energia em termos das grandezas físicas que são
usadas nos cálculos de escoamentos: pressões, velocidades, cotas, potências de bombas e
perdas por atrito viscoso.
O primeiro passo nesse sentido é a divisão do trabalho realizado pelo sistema sobre o
ambiente. Um sistema realiza trabalho sobre o ambiente quando ele se contrai ou expande.
Nesse caso são as forças de contato ao longo da sua superfície as responsáveis pelo
trabalho. Forças de contato originam-se das pressões e das tensões de cisalhamento ao
longo da superfície do sistema.
Além disso, trabalho pode ser extraído ou retirado de um sistema fluido sem que suas
fronteiras se movam, por meio do torque existente em um eixo com pás. É o caso das
bombas e turbinas.
trabalho de eixo (bomba ou turbina)
W =W + W + W
p t e
trabalho de forças tangenciais (cisalhamento)
trabalho de forças normais (pressão) 6.3
______ trabalho das forças de pressão

Imagine um sistema que está se expandindo. Por exemplo, os gases no interior de um
cilindro de motor empurrando o pistão, como na Figura 6.1:
t t +  t
pistão de
área A
Fp Fp
p
x
Figura 6.1: Trabalho realizado pelas forças de pressão
O trabalho realizado pelas forças de pressão na fronteira do sistema durante a

expansão é dado por
 
W p  F  x  W p  p A  x 6.4
A taxa de realização de trabalho é dada por W p / t:

W p  x  
 p A  pV  A 6.5
t t
O resultado da eq. 6.5 é válido quando a pressão e a velocidade são constantes ao
longo da área e durante o intervalo t.
No caso geral devemos usar valores instantâneos da velocidade, e forças atuando em
uma área elementar ⃗⃗⃗⃗⃗ . Essas forças são dadas por ⃗⃗⃗⃗ = 𝑝⃗⃗⃗⃗⃗ Assim o trabalho das forças
ao longo de toda a área móvel é:
W p  
t
 
A
p dA V (trabalho das forças de pressão) 6.6
_____ trabalho das forças tangenciais
Podem existir também tensões de cisalhamento atuando na superfície do sistema.

Mas o trabalho das forças tangenciais será anulado sempre que as áreas forem
perpendiculares ao vetor velocidade. Devemos nos lembrar disso quando escolhermos o
V.C. dos problemas, para facilitar a solução. Veja a Figura 6.2:

p
 
FT
V
WT  FT   x
 
Fp
WT  FT  V  t  0
Figura 6.2: Fronteira móvel em que o trabalho das forças tangenciais é nulo
Portanto, o trabalho das forças de contato fica reduzido apenas ao trabalho das forças
de pressão nas entradas e saídas do V.C., desde que as áreas de entrada e saída sejam
perpendiculares às velocidades. Assim, temos que:
W W p We   We
t

t

t
  AE  AS
p dAV 
t 6.7
As áreas de entrada e saída podem ser substituídas por toda a superfície do V.C.,
pois ao longo das áreas laterais a velocidade é, por definição, perpendicular à área e o
produto escalar ⃗⃗⃗⃗⃗ 𝑉
⃗ é nulo.
W   We
t
  SC
p dA V 
t 6.8
Usando o resultado da equação 6.8 para substituir o termo do trabalho na equação

6.2 vemos que:
Q W p  We   
t

t

t

t VC
e  dVol   SC
e  V  dA 6.9
Q W e  p  
t

t

t VC e  dVol   SC  
 e  V  dA 6.10
Balanço Integral da energia

Os termos da equação do balanço integral da energia, eq. 6.10, representam trabalho
por unidade de tempo. Assim, todos os termos equivalem a taxas de variação da energia no
tempo. No caso do segundo termo do segundo membro essa variação da energia no tempo
decorre da diferença de fluxos de entrada e saída de trabalho das forças de pressão (p/) e
energia interna (e) nas fronteiras do sistema. Suas dimensões são dadas por:
ML
L
 W  FL T2 ML 2
     6.11
 t 
3
T T T
As unidades no SI na equação 6.10 são de joules por segundo, ou Watt.
6.1 Aplicação a um V.C. em Regime Permanente
As situações mais comuns na vida prática são as que envolvem um VC com apenas
uma entrada e uma saída, operando em regime permanente, conforme a Figura 6.3. Vamos
considerar também um trabalho de eixo, já que é muito comum a presença de bombas ou
turbinas nas aplicações hidráulicas.
We
A entrada
A lateral
V1
2
seção 2 : u
1 2
z
seção 1 : u1 2
z
1 V.C. V2
A saída
Figura 6.3: VC para aplicação do balanço integral da energia
Como o regime é permanente, a integral de volume da equação 6.10 desaparece.

Usando os valores médios das grandezas do escoamento nas áreas 1 e 2, de entrada e
saída, podemos efetuar a integral de superfície. Com esses passos a equação da energia
fica:
Q  We
   p2
 e2   2 V2 A2   p1
 e1  1 V1 A1
t t 2 1 6.12
Podemos rearranjar os termos da equação 6.12 e escrever a energia específica “e”

em termos de suas componentes cinética, potencial e interna:
Q p1 V 12 W e p2 V 22
(   gz 1  u 1 )  1V 1 A 1  (   gz 2  u 2 )  2 V 2 A 2
t 1 2 t 2 2
6.13
Se aplicarmos agora o princípio da conservação da massa para escoamentos

permanentes, veremos que VA = VA = FM . Podemos dividir então os dois membros
da equação 6.13 pelo fluxo de massa.
Q p1 V 12 W e p2 V 22
    gz 1  u 1       gz 2  u 2  6.14
 t  VA 1 2  t  VA 2 2
O fluxo de massa multiplicado por um intervalo de tempo t fornece a massa escoada

durante o tempo:
m
 t V A   t  m 6.15
t
Com 6.15 em 6.14 vemos que os termos iniciais dos dois membros da equação 6.14
representam, respectivamente, quantidade de calor e trabalho por unidade de massa
acrescentados ou retirados do V.C.
Q p1 V 12 W e p2 V 22
    gz 1  u 1       gz 2  u 2  6.16
m 1 2 m 2 2
Todos os termos da equação do balanço global da energia representam trabalho, ou

energia, por unidade de massa do escoamento. Suas dimensões são [L 2T-2] e as unidades
no sistema SI são m2/s2.
Partindo da forma básica surgem simplificações para as aplicações mais comuns na
prática, visto que nem todos os termos têm importância igual em todos os problemas.
6.2 Problemas isotérmicos: bombas, turbinas hidráulicas e tubulações
O escoamento dos fluidos reais sempre leva a gradientes de velocidade e ao

aparecimento de tensões de cisalhamento. Sempre existe a geração de calor pelo trabalho
das forças viscosas dissipativas, mesmo que não existam fontes de calor no volume de
controle. Nessas aplicações hidráulicas é conveniente agrupar o termo do calor aos da
energia interna e o balanço de energia, a partir da equação 6.16, fica:
V12  We p2 V22 Q
 gz2  u2  u1  
p1
  gz1    6.17
1 2  m 2 2 m
A ação das forças de atrito viscoso numa partícula que escoa entre as seções 1 e 2
provoca seu aquecimento. Portanto u2 - u1 >0. Entretanto, parte desse calor gerado é
perdido para o ambiente, o que torna o termo do calor (Q/m) negativo pela nossa
convenção. Assim, verifica-se que o último termo da equação 6.17 é sempre positivo.
Uma vez que a transformação de trabalho em calor é um processo irreversível, a
quantidade de energia representada pelo último termo não pode mais ser recuperada para
trabalho útil. Dizemos então que o termo final da equação representa uma perda de energia
do sistema. Trata-se, em outras palavras, da irreversibilidade do processo, ou seja, a
diferença entre o trabalho realizado ao longo de um caminho reversível e o real.
Usando a notação de perdas para o termo das irreversibilidades do escoamento, a
equação da energia fica:
p1 V12  We p2 V22
 gz1     gz2   Perdas1 2 6.18
1 2  m 2 2
Os termos da equação 6.18 representam energia por unidade de massa, ou joules por
kilograma no SI.
Dividindo-se a equação 6.18 pela aceleração da gravidade o resultado não se altera e
a equação fica expressa em termos chamados de Carga, muito utilizados na prática quando
se trata de escoamentos incompressíveis. Em termos de carga, o balanço de energia fica:
𝑉2 2 𝑉22
+ 𝑧 + − 𝑀 = + 𝑧2 + + 2 6.19
2 2
A carga representa energia por unidade de peso do fluido em escoamento, e é

expressa em metros no sistema SI, dimensões [L]. A carga de eixo (He) é chamada também
de altura manométrica da bomba ou turbina. Lembre que a altura manométrica de uma
bomba hidráulica é negativa, pois a bomba adiciona energia ao fluido em escoamento.
Inversamente, a altura manométrica de uma turbina é positiva. O termo das perdas
dissipativas por unidade de peso é chamado de perda de carga, e é sempre positivo.
O termo dissipativo (Perdas) não é calculado com a termodinâmica e sim com
equações independentes, baseadas na dinâmica dos escoamentos. É necessário conhecer
as tensões de cisalhamento no interior do escoamento, o que só é possível em escoamentos
laminares. Essas demonstrações são encontradas em textos mais avançados de mecânica
dos fluidos. No escoamento turbulento, entretanto, é indispensável usar resultados
experimentais para avaliar as perdas por dissipação.
Várias fórmulas empíricas foram ajustadas pelos pesquisadores em hidráulica. O item
6.4 apresenta algumas fórmulas para o cálculo da perda de carga.
6.3 Exemplos Ilustrativos
EXEMPLO 6.3.1: ______ troca de calor em escoamento incompressível

Um aquecedor solar com área exposta ao sol de 4m 2 é utilizado para aquecer água. Em
determinado instante a placa capta a radiação solar com densidade de fluxo de 500W/m 2. A
água a ser aquecida passa através do aquecedor com velocidade de 0,5m/s no interior de
tubos de 10 mm de diâmetro interno. Sabe-se que a água entra no aquecedor com
temperatura de 20°C, e que a eficiência do aquecedor é de 45%, porque apenas parte do
calor captado é transferido para a água. Supondo que a massa específica da água
permaneça constante (ρ = 1000kg/m3), pede-se escrever a equação da energia para o
problema e calcular a temperatura de saída da água. Dado: c = 4185 J/(kg°C).
Análise:
O escoamento é considerado incompressível, visto que a diminuição da massa
específica com a temperatura é desprezada. O volume de controle a ser usado (conjunto de
captador solar e serpentina de tubos que formam o aquecedor) possui apenas uma entrada
e uma saída. Não há bombas ou turbinas no interior do V.C. de modo que o trabalho de eixo
é nulo. Além disso, o regime de escoamento é permanente. Uma hipótese simplificadora
importante.
Solução:
Partindo da equação básica 5.10
− = ∫ V + ∫ ( + ) ⃗V ⃗⃗⃗⃗⃗
Descartando os termos de trabalho de eixo e de variação no tempo
= ∫ ( + ) ⃗V ⃗⃗⃗⃗⃗
Considerando que o VC só tem uma entrada (seção 1) e uma saída (seção 2) de

água, a integral em SC é dividida
= ∫ ( + ) ⃗V ⃗⃗⃗⃗⃗ + ∫2 ( + ) ⃗V ⃗⃗⃗⃗⃗
Considerando que as grandezas são igualmente distribuídas pelo interior do tubo nas
seções de entrada e saída, os termos da energia específica “e” e do trabalho nas fronteiras
(p/ρ) pode sair da integral
= ( + )∫ ⃗ ⃗⃗⃗⃗⃗ + (
V 2
+ 2 ) ∫2
⃗ ⃗⃗⃗⃗⃗
V
Considerando que foi dada a velocidade média nas seções de entrada e saída,
avaliamos as integrais:
2
= ( + )(− V ) + ( + 2 )(+ V2 2)
Sendo o escoamento permanente, V = V2 2 e, sendo incompressível,

V = V2 2 , ou seja, vazão de entrada e saída são iguais. Além disso, como a área do
tubo é constante, a velocidade não se altera. Usando essas informações e desdobrando a

energia específica:
2 2
2 2
= *( + 2+ + 2) − ( + + + )+ V
2 2
Podemos desprezar a diferença de pressão, assim como a diferença de cota entre a

entrada e saída. Assim, a equação do balanço global de energia, com as simplificações do
problema em questão fica:
= (𝑢2 − 𝑢 ) 𝑉 = ( 2 − ) 𝑉
𝑡
Os valores numéricos conhecidos são:
𝑄
𝑡
=  2
 2
= 𝑊
00 2
𝑉 =   2
=
Substituindo os valores numéricos na equação do balanço de energia:

𝑊=   ( 2 − ) 𝐶
𝐶 𝑠
Resolvendo obtém-se: 2 = 𝐶
Resposta: nas condições dadas a água deixará o aquecedor com temperatura de 22,1°C.
Comentário: A temperatura de saída ficou muito baixa. Observando-se a equação deduzida
para o problema percebe-se que, para melhorar o resultado, é possível (a) aumentar a área
do coletor solar ou (b) diminuir a vazão de água através dos tubos.
EXEMPLO 6.3.2: ______ escoamento compressível com troca de calor
Uma turbina a vapor usa 4.600 kg/h de vapor e entrega 700 kW de potência a um gerador
elétrico, conforme a figura. Os dados de entrada e saída do vapor são dados a seguir. Pede-
se calcular a perda de calor através da carcaça da turbina e nos mancais. Calcular a efi-
ciência da turbina. Entrada: h1 = 2790kJ/kg; V1 = 60m/s. Saída: h2 = 2093kJ/kg; V2 = 270m/s.
Análise: Com o volume de controle desenhado na figura, aplicaremos o balanço de energia

em regime permanente.
Nos casos de escoamentos compressíveis com variação de temperatura entre a
entrada e saída do VC, é útil escrever a equação explicitando a entalpia h (h = u + p/.
Nesses casos, após as simplificações para regime permanente em um VC com apenas uma
entrada e uma saída, a equação 5.10 é transformada na eq. 5.14 fica:
Q p1 V 12 W e p2 V 22
    gz 1  u 1       gz 2  u 2 
 t  VA 1 2  t  VA 2 2
Pelo fato do escoamento ser compressível, os termos p/ρ e u são somados para formar a
entalpia do gás, de forma que a equação do balanço fica:
𝑉2 𝑊 𝑉22
+ ( + 𝑧 + ) = + ( + 𝑧2 + 2)
𝑡 𝑉 𝑡 𝑉
Nas condições do problema a diferença de energia potencial pode ser desprezada. Os

demais termos foram fornecidos.
Solução:
a) Perda de calor:
Substituindo os valores numéricos na equação acima, temos:
= 𝑉 = = 𝑠
𝑠
𝑊 𝑉 𝑉
= + [( + )− ( + )]
𝑡 𝑡
= + [( − ) +( − ) ]
𝑡 𝑠 𝑠
= − 𝑊
𝑡
O sinal negativo indica que o calor deixa o volume de controle.
b) Eficiência
O cálculo da eficiência usa o fato de que quando há 100% de eficiência a perda de calor é
nula. Assim, o trabalho reversível é dado pela energia total retirada do vapor, conforme a
equação:
𝑊 𝑣 𝑉 𝑉
− = [( + )− ( + )]
𝑡
𝑊 𝑣
= 𝑊
𝑡
O trabalho reversível pode ser calculado diretamente, somando o trabalho líquido (700kW)
às perdas (147kW).
A eficiência é dada pelo trabalho líquido em relação ao trabalho máximo (reversível)
𝑊
= = =
𝑊 𝑣
Portanto, 82,6% de eficiência da turbina.
EXEMPLO 6.3.3: ______ escoamento incompressível isotérmico

A figura mostra uma bomba que retira água de um grande reservatório por meio de uma
tubulação de 150mm de diâmetro, e descarrega em jato livre num ponto 30m acima do eixo
da bomba.
A pressão relativa na seção 1 (sucção da bomba é de 100 kPa) e a velocidade é 6m/s.

Sabendo que a velocidade no bocal de saída é 10m/s e que as perdas de carga na
tubulação são de 5 m.c.a., calcular a potência fornecida pela bomba. Sabendo que o
rendimento da bomba é de 65%, calcular a potência consumida pela bomba.
Análise:
O escoamento é considerado incompressível e permanente. Adotamos como volume de
controle a água na tubulação entre as seções 1 e 2. O jato descarrega com pressão
atmosférica (jato livre), portanto com pressão relativa nula. O balanço de energia para este
caso é dado pela equação 5.19 em termos de trabalho por unidade de massa:
p1 V12  We p2 V22
  gz1     gz2  Perdas1 2
1 2  m 2 2
𝑃 𝑎 2
Sabendo que = , temos que 𝑃 𝑟 𝑠 2 = 2 ( 2
)
Solução:
Adotando massa específica da água ρ = 1000 kg/m3 e substituindo valores na equação
acima,
2 2
𝑊
+ + − = + + +
𝑊
=− ( )
O termo é negativo porque se trata de uma bomba que realiza trabalho sobre a água.
A potência fornecida pela bomba é calculada multiplicando-se o termo do trabalho de eixo
por unidade de massa pelo fluxo de massa do escoamento.
𝑊
𝑃𝑜𝑡 = 𝑉
𝑃𝑜𝑡 = 𝑥 = 𝑊
𝑃𝑜𝑡 0
A potência consumida : 𝑃𝑜𝑡 𝑜 𝑎 = = = 𝑊
 0
Obs: a potência consumida pela bomba é maior do que a transferida ao líquido devido
principalmente a perdas por atrito nos mancais e por recirculação do líquido no interior da
bomba, das zonas de alta para as de baixa pressão.
6.4 Efeito do Atrito nos Escoamentos
Podemos fazer a analogia entre o escoamento no interior de um tubo e um trem

sendo empurrado nos trilhos por uma locomotiva, conforme o esquema da Figura 6.4.
1 2 3 4 5
V
F
Figura 6.4: Diagrama de forças aplicadas a cada ligação entre vagões

empurrados por uma locomotiva.
Na analogia entre trem e escoamento, cada vagão pode ser imaginado como uma
certa quantidade de massa do fluido em escoamento. As forças transmitidas para cada
vagão pelos vínculos correspondem no escoamento às forças resultantes das pressões que
atuam em cada face dos elementos de massa.
Supondo massas e atrito iguais em todos os vagões, podemos calcular as forças nos
vínculos entre os vagões, que são as barras 1,2,3,4 e 5 vistas na Fig.6.4. Nada impede que
as forças sejam colocadas num gráfico da posição dos vagões.
Na Figura 6.4 a energia para o movimento é fornecida pela locomotiva, e consumida
pelo trabalho dissipativo do atrito ao longo do trem. Com isso as forças nos vínculos entre
vagões diminuem.
No caso de fluido escoando num tubo, as forças nos vínculos são substituídas por
distribuição de tensões ao longo da área de contato entre os elementos de massa, conforme
o esquema da Figura 6.5.
F V
1 2 3 4 5 6
l
(a)
p Diagrama de Forças
F
5-4
F
3-4
P 4
3-4 P
5-4 F
at
(b) p (c)
Figura 6.5: Fluido escoando em um tubo - (a) analogia entre as massas de fluido e os vagões do trem;
(b) isolando o elemento de fluido; (c) diagrama de forças resultantes no elemento 4.
Pensando em cada uma das porções de massa de fluido, representadas pelos blocos
numerados da Figura 6.5, vemos que elas recebem energia para o escoamento ( fornecida
por uma bomba ) e, à medida que passa o tempo, essa energia é dissipada pelo trabalho
das forças de atrito, fazendo com que a pressão diminua. Podemos então equacionar as
diferenças de forças entre os vínculos em função do atrito:
= + 𝑎𝑡
Usando AL = Área lateral e AT = Área de seção transversal vemos que:
𝑎𝑡 = 𝑝
𝑝
= 𝑝 𝑇] 𝑝 = 𝑝 −
𝑇
= 𝑝 𝑇
Observe que o bloco 5 de fluido já esteve anteriormente na posição do bloco 1. Com o

avanço do tempo ele se desloca e perde pressão. Entretanto, outra massa trazida pelo
escoamento ocupa a posição do bloco 1 quando ele se desloca para a posição do bloco 2 e
assim sucessivamente.
Assim, embora as unidades de massa ( blocos ) sejam transportadas no espaço ao
longo do tubo, o perfil de pressões visualizado em um determinado trecho do tubo
permanece constante no tempo. Este raciocínio leva ao perfil de pressões no tubo conforme
a Figura 6.6:
P Perfil de Pressões
observadas no tubo.
x
1 2 3 4 5
Figura 6.6: Diagrama de pressões resultantes no escoamento de fluido num tubo.

devido à viscosidade do fluido (forças de atrito).
Veja que a analogia não é perfeita devido a essa maior liberdade de raciocínio que
temos ao analisar os escoamentos, já que uma unidade de massa pode sempre ser
substituída por outra, num determinado local de análise, à medida que passa o tempo
No trem, as forças permanecem constantes ao longo do tempo em cada vagão
(Sistema – Análise Lagrangeana). Com o fluido, as pressões permanecem constantes em
uma dada região do tubo (Volume de Controle – Análise Euleriana).
6.4.1 Perdas de Carga em Escoamento em Tubos
A energia por unidade de peso do fluido em escoamento, ou carga total, foi definida
pela equação 4,14, assim como seus componentes, a carga piezométrica e a carga cinética.
Os termos da equação da energia podem ser representados graficamente num escoamento
em tubos pelas Linhas de Carga, conforme o esquema da Figura 6.7.
Figura 6.7: Linhas de carga num trecho de tubo: H = carga total, H = perda de carga.
A equação da energia, na forma de cargas da Equação 6.19, pode ser aplicada entre
as seções 1 e 2 do escoamento da Figura 6.7.
p1 V12 P V2
 z1   H M  2  z 2  2  H 12
 2g  2g
Como não há trabalho de eixo a equação fica:

𝑉2 2 𝑉22
+ 𝑧 + = + 𝑧2 + + 2
Desta forma, no caso do escoamento em tubo da Figura 6.7, a perda de carga é

escrita como:
2 2
2 2
2 =( + 2 + ) − ( + + ) 6.20
2 2
As perdas de energia por unidade de peso, ou perdas de carga, numa tubulação com
escoamento forçado, podem ser expressas pela chamada equação universal:
L V2
H 12  f 6.21
D 2g
em que a perda é dada em metros, f é um adimensional chamado fator de atrito, L é o
comprimento em metros entre as seções 1 e 2 e D é o diâmetro da tubulação em metros.
O fator de atrito f é determinado por ajuste de dados experimentais. Sabe-se que f
depende do número de Reynolds do escoamento e da rugosidade relativa /D. A utilização
das fórmulas originais de regressão para o fator de atrito é dificultada porque são equações
implícitas, ou seja, com a incógnita nos dois membros.
Várias aproximações explícitas para o fator de atrito f foram propostas ao longo da
segunda metade do século 20. Uma aproximação bem simples e com faixa de validade
bastante ampla foi apresentada por Swamee e Jain (1976), conforme a equação 6.22.
1 , 325
f  6.22
2
   5 , 74 
 ln   
  3 ,7 D Re 0 , 9 
  
sendo  a rugosidade absoluta da parede do tubo, D o diâmetro interno do tubo e Re o

adimensional número de Reynolds.
O número de Reynolds é dado por:
= = 6.23
Sendo: ρ (kg/m3) a massa específica;  (kg/m.s) a viscosidade e  = /ρ (m2/s) a

viscosidade cinemática.
A fórmula de Swamee e Jain possui um erro de no máximo 1% em relação ao
diagrama de Colebrook- Moody. Os limites de utilização são: 10-6<  / D < 10-2 e
5000 < Re < 108.
Para escoamentos laminares o fator de atrito é :
64
f  ; válida para Re < 2000 6.24
Re
Para facilitar a solução de problemas eventuais que requeiram a determinação do

fator de atrito, pode se utilizada uma solução gráfica, com os valores numéricos do fator de
atrito fornecidos pelo diagrama de Colebrook – Moody, apresentado na Figura 6.8.
Figura 6.8: Diagrama de Colebrook – Moody, apresentando o comportamento do fator de atrito f.
O diagrama de Moody mostra que o fator adimensional de atrito possui um

comportamento relativamente complexo, em função do tipo de escoamento, da rugosidade
do tubo e do número de Reynolds. A explicação deste comportamento é apresentada com
mais detalhes no item 6.4.2, com auxílio da análise dimensional do problema.
Exemplo 6.4.1: Um reservatório de nível constante alimenta uma tubulação de ferro

galvanizado (e = 0,15mm) de 250 mm de diâmetro e 500m de comprimento, que conduz
água ( = 1,1 x 10 –6 m2/s) até um reservatório com o nível de água situado 20m abaixo. Uma
vazão de 0,15m3/s é regulada por um registro no final da tubulação. Pede-se calcular: a) a
perda de carga na tubulação; b) a pressão antes do registro de regulagem.
Análise:
Será necessário calcular a perda de carga na tubulação e, para isso, vamos supor
escoamento turbulento, aplicando a equação de Swamee e Jain (6.22). A carga antes do
registro será calculada pela equação da energia aplicada entre um ponto qualquer da água
parada no reservatório e o ponto imediatamente a montante do registro de saída da
tubulação.
Inicialmente, calculamos o número de Reynolds do escoamento:

VD VD 4Q 4  0,15
Re      694494
  D   1,1106  0,25
Como Re está dentro dos limites de utilização, a eq. 6.22 pode ser usada:
1,325
f  2
 0,0181
  0,15 5,74 
ln  
0,9 
  3,7 * 200 694494 
A perda de carga é dada pela eq. 6.21. Com os valores numéricos temos:
500 8  0,152
HTubo  0,0181   17,28m
0,25 9,8   2  0,255
Resposta item a: Cada Newton de fluido perdeu 17,28 Joules de energia devido ao
atrito viscoso, ao percorrer os 500m de tubo. Essa energia foi transformada em calor.
Resta calcular a pressão a montante do registro, aplicando a equação da energia.
Chamamos o ponto inicial do reservatório de “1”, tomado, por simplicidade, na superfície do
reservatório, e de “2” o ponto imediatamente a montante do registro.
𝑉2 2 𝑉22
+ 𝑧 + − 𝑀 = + 𝑧2 + + 2 (eq. 6.19)
2 2
Como não há trabalho de eixo a equação fica:

𝑉2 2 𝑉22
+ 𝑧 + = + 𝑧2 + + 2
Substituindo os valores numéricos, lembrando que V = QA = 3,06m/s, vem:

p 3,062 p
0  20  0  2  0   17,28  2  20  17,28  0,48  2,24m
 2  9,8 
Tendo a carga equivalente a 2,24m, o cálculo da pressão é imediato:
𝑝2 = = . = . 𝑃
Resposta item b: A pressão imediatamente a montante do registro é 21,5kPa.
Comentários adicionais: A perda de carga no tubo permite calcular a potência total

dissipada pelo atrito viscoso ao longo do comprimento do tubo. Como a perda de carga é a
energia perdida por cada Newton que escoa, basta multiplicá-la pelo fluxo de água em
Newtons por segundo, como segue:
𝑃𝑜𝑡 𝑖 𝑖𝑝𝑎 𝑎 𝑇𝑢 𝑜 = 2 = = . 𝑊
Imediatamente após o registro está a saída do tubo com pressão nula, e a carga
cinética é a mesma antes e depois do registro, porque o escoamento é permanente. Com
isso deduzimos que a perda de carga no registro é igual a 2,24m. Para que a vazão seja
regulada o registro dissipa uma potência dada por:
𝑃𝑜𝑡 𝑖 𝑖𝑝𝑎 𝑎 𝑅 𝑖 𝑡 𝑜 = 𝑅 𝑖 𝑡 𝑜 = = . 𝑊
6.4.2 Sobre o Fator de Atrito
Todas as grandezas da equação 6.21 podem ser medidas experimentalmente, para

determinar o coeficiente de atrito, por meio de ensaios de simples execução. Um esquema
experimental típico é mostrado na Figura 6.9.
1 2
Venturi
x
Figura 6.9: Arranjo experimental para determinação de perdas de carga.
A seção experimental situa-se entre os pontos 1 e 2 da Figura 6.9, sendo alimentada

por água em circuito fechado por meio da bomba centrífuga. A vazão é regulada pelo
registro e medida pelo tubo Venturi. A aplicação da equação da energia entre a seção 1 e 2 é
apresentada na equação 6.20.
2 2
2 2
2 =( + 2+ ) − ( + + )
2 2
Nas condições do ensaio a carga cinética e a carga potencial se mantêm constantes

entre as seções 1 e 2 porque o tubo é horizontal e com diâmetro constante.
Dessa forma, a dissipação de energia é indicada pela diferença de pressão medida
pelo manômetro diferencial.
2
2 = −
Os resultados obtidos nos ensaios mostraram as seguintes características da perda
de carga:
1. A perda de carga não depende da pressão do tubo;
2. A perda de carga é diretamente proporcional ao comprimento do tubo;
3. A perda de carga depende da rugosidade do tubo no escoamento
turbulento, mas não no escoamento laminar;
4. A perda de carga diminui com o aumento do diâmetro e aumenta com a
velocidade do escoamento;
5. A perda de carga depende da viscosidade do fluido.
A partir das informações experimentais e analisando o problema físico podemos

afirmar que o fator de atrito irá depender de 7 variáveis dimensionais, ou seja:
= (𝑉 𝜇 𝜀 𝜀′ ) 6.25
Na relação 6.25 temos:  é uma medida do tamanho da rugosidade e tem dimensão L

(comprimento), ’ é uma medida do arranjo da rugosidade e também tem dimensão L, m é
um fator de forma, que depende do formato das rugosidades e é adimensional. A Figura 6.10
mostra um exemplo de como rugosidades de mesmo tamanho podem ser dispostas em
diferentes arranjos e como a forma da rugosidade pode influenciar no escoamento.
 
’
 
’
Arranjo Fator de forma
Figura 6.10: Ilustração de diferenças de arranjo e de forma de rugosidades de mesmo tamanho.
Além das 3 variáveis descritoras da rugosidade o fator de atrito depende da

velocidade V do escoamento, do diâmetro D do tubo, e do fluido em escoamento,
representado pela massa específica  e pela viscosidade dinâmica . Como f é
adimensional, ele deve ser função das 7 variáveis agrupadas de forma adimensional.
Pensando na situação de tubo liso,  , ’ e m são nulos, de forma que o fator de atrito
depende apenas das 4 primeiras grandezas agrupadas num adimensional. Da análise das
dimensões das grandezas, ou pela aplicação do Teorema  da análise dimensional,
podemos concluir que este adimensional é o número de Reynolds.
Para adimensionalizar os descritores da rugosidade basta dividir pelo diâmetro, ou
seja, usar tamanhos relativos. Assim, a descrição do fator de atrito em termos adimensionais
deve ser uma função do tipo:
𝑉 𝜀 𝜀′
= ( ) 6.26
𝜇
Embora a análise dimensional seja útil para definir os parâmetros adimensionais que
irão descrever o fenômeno em qualquer escala, ela nada fala sobre a forma da função f, que
dever ser determinada experimentalmente.
Blasius, em 1913, foi o primeiro a correlacionar os dados experimentais em
escoamento turbulento, mas apenas para tubos lisos. A fórmula empírica de Blasius
correlacionou os resultados dos ensaios com uma dispersão de  5%, sendo dada por:
0
= 4 6.27
𝑅
O efeito do parâmetro /D, chamado de rugosidade relativa, foi esclarecido por
Nikuradse, em 1933. Nikuradse utilizou tubos de vidro revestidos internamente com grãos de
areia de tamanho uniforme. Os ensaios de Nikuradse demonstraram que o fator de atrito de
tubos rugosos comporta a divisão em 3 regiões, que podem ser observadas na Figura 6.8 e,
de forma mais esquemática, na Figura 6.11. Na primeira região, chamada de escoamento
turbulento hidraulicamente liso, o tubo possui o mesmo fator de atrito de um tubo liso; nesta
região f não depende da rugosidade relativa do tubo. Com o crescimento do número de
Reynolds, o fator de atrito começa a depender tanto do número de Reynolds quanto da
rugosidade; esta segunda região é chamada de escoamento turbulento de transição. Na
terceira região, chamada de escoamento turbulento hidraulicamente rugoso, o fator de atrito
depende apenas da rugosidade relativa.
Com os dados de Nikuradse para tubos lisos e para a região de escoamento
hidraulicamente liso dos tubos rugosos, a equação do fator de atrito é dada por:
1
 0,86 ln (Re f )  0,8
f
6.28
A equação também vale para tubos de pequenas rugosidades relativas, quando o

número de Reynolds é baixo. O escoamento na região em que o fator de atrito de um tubo
rugoso fica igual ao de um tubo liso é chamado de Escoamento Hidraulicamente Liso.
Para a região de turbulência completa, ou de escoamento hidraulicamente rugoso, a
equação de regressão dos dados experimentais é:
1  
 1,14  0,86 ln  
f  D 6.29
Observe que na região de turbulência completa, ou de Escoamento

Hidraulicamente Rugoso, o fator de atrito não depende de Re, e as curvas de regressão
ficam horizontais no diagrama de Moody.
Na zona de transição entre os escoamentos hidraulicamente lisos e rugosos,
entretanto, os resultados de Nikuradse não são válidos para tubos comerciais, que possuem
uma rugosidade composta de diversos tamanhos e formas, dispostos em padrões que
dependem da técnica de fabricação do tubo. O comportamento de tubos rugosos comerciais
foi estudado por Colebrook que apresentou a equação para a zona de transição entre
escoamento hidraulicamente liso e escoamento hidraulicamente rugoso, em 1939.
1   2,51 
  0,86 ln   
 3,7 D Re f  6.30
f  
Finalizando a descrição do comportamento do fator de atrito, para escoamentos
laminares, temos uma equação que pode ser deduzida teoricamente, dada por:
64
f  6.31
Re
A explicação para o comportamento complexo do fator de atrito surge ao observar o
tamanho relativo das rugosidades  em relação à espessura da camada limite .
A Figura 6.11 reproduz as características principais do diagrama de Moody, realçando
as diversas subdivisões de comportamento do fator de atrito.
0,08
0,05
Rugosidade Relativa ( /D)

Zona Crítica
Fator de Atrito f
0,04
0,01
Turbulência Completa
0,03
0,02 0,001
0,0002
0,0001
0,01
103 104 105 106 107 108
Número de Reynolds - Rey
Figura 6.11: Diagrama de Moody realçando as 4 regiões de comportamento diferente do fator de atrito.
A Figura 6.12 ilustra o tamanho das rugosidades em relação ao tamanho da camada

limite laminar, ou da subcamada limite laminar, no caso dos escoamentos turbulentos. Com
o aumento da velocidade, que reflete no aumento do número de Reynolds, a espessura da
camada limite diminui, fazendo com que as rugosidades comecem a interferir no
escoamento. Na região de turbulência completa toda a rugosidade fica exposta.
Laminar Hidraulicamente Liso

Subcamada Limite Amortecedora
Camada Limite Laminar Subcamada Limite Laminar

 '
 
 >>  ’ > 
Transição Hidraulicamente Rugoso
Subcamada Limite Amortecedora

'  '
’   ’ << 
Figura 6.12: Tamanho relativo das rugosidades e camada limite nas diversas regiões do
Diagrama de Moody determinam o comportamento do fator de atrito.
A Figura 6.12 permite compreender também porque os tubos comerciais tem

comportamento diferente dos tubos de rugosidade uniforme na região de transição. Nos
tubos com rugosidade artificial de Nikuradse as asperezas vão emergindo uniformemente da
subcamada laminar, ao passo que nos tubos comerciais as asperezas tem tamanhos
irregulares. Por esta razão os dados de Nikuradse não podem ser empregados para prever o
fator de atrito de tubos comerciais na região de transição. Apenas com os resultados de
Colebrook (eq. 6.30) a descrição do fator de atrito ficou completa.
A Figura 6.13 apresenta um resumo das equações utilizadas para descrever os
resultados experimentais obtidos nos ensaios de Nikuradse e de Colebrook.
Resumo das Fórmulas do Fator de Atrito f
Escoamento Laminar ( Re < 2.300 )

64
f 
Re
Escoamento Turbulento ( Re > 4.000 )
Tubos Lisos e Hidraulicamente Lisos
1
 0,86 ln (Re f )  0,8
f ’ < 1,9 
Tubos Rugosos
Turbulência Completa (hidraulicamente rugoso)
1  
 1,14  0,86 ln  
f  D ’ < 0,008 
Região de Transição (Colebrook-Moody)
1   2,51 
  0,86 ln   
 3,7 D Re f 
f   0,008 < ’ < 1,9 
Figura 6.13: Resumo das equações do fator de atrito e faixa de validade.
Existe uma razoável dificuldade prática de utilização das diversas equações que
descrevem o comportamento do fator de atrito, visto que várias são implícitas, ou seja, com
o fator de atrito aparecendo nos dois membros da equação. Devido à dificuldade de
utilização, o esforço para desenvolver fórmulas explícitas envolveu diversos pesquisadores.
A fórmula de Samee e Jain apresentada na equação 6.22 é um exemplo desse
desenvolvimento. Entretanto, embora seja bem simples e aplicável na maioria das situações
práticas, ela apresenta limitações, ou seja, não reproduz o comportamento completo em
todas as faixas do número de Reynolds.
Dentre as fórmulas válidas para toda a faixa de número de Reynolds, incluindo o
escoamento laminar, podemos citar as de Churchill (1977)1 , Chue (1984)2 , e a de Pereira e
Almeida (1986)3 .
1
CHURCHILL, S.W, 1977. Friction fator equation spans all fluid regimes. Chem. Engng., 84-7, pp.
91-92.
2
CHUE, S.H., 1984. A pipe skin friction Law of universal applicability. Proc. Inst. Civil. Engrs., Part
2, 77 – mar 1977, pp. 43-48.
3
PEREIRA, A.J. e ALMEIDA, A.B., 1986. Formulação explícita e universal da resistência em tubos.
XII Congresso Latino-Americano de Hidráulica. São Paulo, S.P.
A fórmula de Swamee (1990), apresentada a seguir, tem obtido boa aceitação:
6.32
6.5 Efeitos das Bombas e Turbinas sobre as Cargas
Aplicando a analogia trem-escoamento podemos representar um escoamento,

mantido por uma bomba hidráulica num tubo por meio do esquema da Figura 6.14.
ANALOGIA DE BOMBA HIDRÁULICA E TUBULAÇÃO
Sucção
Recalque Saída
V
Bomba Mola comprimida

Mola distendida
Figura 6.14: Representação de uma bomba hidráulica com a analogia entre

escoamento e trem de massas unidas por molas
A bomba pode ser encarada como um mecanismo que coloca as massas no início do
tubo e as comprime em direção à saída. A energia fornecida pela bomba é usada para
acelerar as massas e para comprimir as molas (aumentar a pressão). Com o passar do
tempo, enquanto as massas deslocam-se no tubo, a pressão cai devido ao trabalho das
forças de atrito, conforme deduzido no item 6.1 e 6.2, e as molas vão se distendendo
progressivamente. Na saída, à pressão atmosférica, é como se as molas estivessem
totalmente distendidas, e só há energia cinética.
_____ Variação da Energia com uma Bomba: Linhas de Carga
Num caso geral a energia adicionada pela bomba pode provocar variações na carga
de pressão e na carga cinética. Essas transformações de energia são visualizadas por meio
dos diagramas de carga, que contêm as Linhas de Carga Piezométrica (LP) e Linha de
Carga Total (LC).
A Linha de Carga Piezométrica indica em cada ponto a altura que a água do tubo
alcançaria caso fosse instalado um piezômetro com tomada estática naquele ponto do tubo.
𝑃
𝑃 = + 𝑧
A Linha de Carga Total, como o nome indica, mostra a altura que seria atingida pela
água num piezômetro conectado a um Tubo de Pitot (tomada dinâmica). A distância entre as
duas linhas representa a carga cinética HC.

𝑃 𝑉2
𝑇 = + 𝑧 + = 𝑃 + 𝐶
A Figura 6.15 mostra o caso em que a energia introduzida provoca apenas aumento
da pressão.
LC
Carga (m)
HM LP
Hc x
VS VR
Sucção Bomba Recalque VS = VR
Figura 6.15: Bomba adiciona energia que se manifesta no aumento de pressão, já que a energia cinética
permanece constante.
As bombas usadas em tubos para aumento de pressão são chamadas de “Boosters”.

Normalmente as bombas são instaladas na saída de reservatórios e com um diâmetro
maior na tubulação de sucção porque a energia disponível na sucção não é grande. Esse
caso é ilustrado pela Figura 6.16. Com uma carga cinética menor a pressão pode ser maior,
para uma dada energia total. Este efeito é conveniente para evitar a cavitação, um fenômeno
que ocorre quando a pressão fica próxima à pressão de vapor.
L.C.
Carga (m)
Hc
HM
L.P.
VS < VR
EC,S < EC,R
x
VS VR
Sucção Bomba Recalque

Figura 6.16: Com diâmetros são diferentes, a energia fornecida pela bomba transforma-se em aumento
de pressão e de energia cinética.
A bomba da Figura 6.16 mostra um caso mais comum do que o booster da Figura
6.15. A energia fornecida pela bomba é dividida entre aumento de carga cinética e aumento
de pressão.
O efeito de redistribuição entre formas de energia não depende da bomba existir,
acontece também nas tubulações sempre que varia o diâmetro. Veja o esquema da Figura
6.17.
HTotal
Carga (m)
Hc H
LC
Hc LP
P/
P/
x
V1 V2 V1
Figura 6.17: Quando a velocidade aumenta, a pressão diminui pois a energia disponível é limitada.
O efeito de diminuição de pressão num escoamento quando aumenta a velocidade é

chamado de “efeito Venturi”. A distribuição da energia total disponível no fluido entre as
formas de energia cinética e pressão é descrita pela equação de Bernoulli, já deduzida.
6.6 Resumo das Transformações de Energia
Entre duas seções de um escoamento permanente:
𝐸 𝑛𝑖 𝑖𝑎𝑙 𝐸 𝑜 𝑛 𝑖 𝑎 𝐸 𝑖𝑛𝑎𝑙 𝐸𝑃 𝑖 𝑎
+ = +
𝑈𝑛 𝑃 𝑠𝑜 𝑈𝑛 𝑃 𝑠𝑜 𝑈𝑛 𝑃 𝑠𝑜 𝑈𝑛 𝑃 𝑠𝑜
A relação anterior expressa em palavras a conservação de energia. Numericamente
os termos de energia por unidade de peso terão a dimensão de espaço [L], com unidades
em (m) metros no SI. A energia por unidade de peso é chamada genericamente de Carga
Hidráulica.
_____ Caso de Bombas
A carga de eixo fornecida é chamada de Altura Manométrica da Bomba (HB). A

energia dissipada pelo escoamento entre as duas seções resulta do atrito, e é chamada de
Perda de Carga (H).
Assim, para um trecho de tubulação com início na seção 1 e final na seção 2, e uma
bomba em um ponto qualquer entre a seção final e inicial a equação 6.19 fornece:
2 2
p V p V
1
z  1
H  z   H
2 2
 
1 B 2
2g 2g
 -He = - (-HB)
_____ Caso de Turbinas

Turbinas são máquinas hidráulicas que retiram energia do escoamento. A carga de
eixo retirada por uma turbina é a energia fornecida pela água por unidade de peso do fluido
escoado. Esta carga também é chamada de Altura Manométrica da Turbina (HT), com a
diferença que é um termo negativo no primeiro membro da equação 6.19.
p1 V12 p2 V22
 z1   HT   z2   H
 2g  2g
 -He = - (+HT)
O termo genérico carga de eixo (He) ou ainda altura manométrica (HM) é usado para
designar uma máquina hidráulica, que pode ser uma bomba ou turbina.
Assim, a equação da energia, na presença de uma máquina hidráulica, seja ela uma
bomba ou uma turbina, fica com os termos cujo significado é realçado a seguir.
Carga Inicial Hi Carga Final HF
p1   V12 p V 2
z1  H M  2  z2  2  H
 2g  2g
_____ Relação entre Cargas e Potências
As bombas e turbinas são especificadas pela potência (W) que consomem ou

produzem. A potência é um fluxo de trabalho, ou seja, trabalho por unidade de tempo. A
carga é trabalho por peso de fluido escoado. Assim temos, pelas dimensões:
𝑟 𝑏 𝑙 𝑜 𝑟 𝑏 𝑙 𝑜 𝑃 𝑠𝑜 𝐸𝑠 𝑜 𝑜
𝑃𝑜𝑡ê𝑛 = =
𝑝𝑜 𝑃 𝑠𝑜 𝐸𝑠 𝑜 𝑜 𝑝𝑜
O peso escoado por unidade de tempo (N/s) é dado pelo peso específico do fluido
“” (N/m3) multiplicado pelo fluxo de volume ou vazão “Q” (m 3/s). Assim, da relação entre as
dimensões é fácil verificar que:
= 6.33
Note que a equação 6.33 define a potência retirada do fluido, ou acrescentada ao

fluido. Para obtermos as potências da máquina hidráulica é necessário considerar o seu
rendimento, conforme detalhado a seguir.
____ Bombas
A potência fornecida ao eixo de uma bomba será maior que a acrescentada ao fluido:
𝑃𝑜𝑡𝐸𝑖𝑥𝑜 =
𝑜 𝑎

em que B é o rendimento da bomba (0 < B <1) .
A potência de eixo da bomba é também chamada de potência bruta da bomba, e
corresponde à potência líquida fornecida pelo motor elétrico. A potência consumida pelo
motor elétrico que aciona a bomba é maior que a potência fornecida ao eixo da bomba:
𝑃𝑜𝑡𝐸𝑖𝑥𝑜 𝑜 𝑎
𝑃𝑜𝑡𝑀𝑜𝑡𝑜 = =
𝐸𝑙é𝑡 𝑖 𝑜
𝑀  𝑀
em que M é o rendimento do motor elétrico (0 < M <1)
____ Turbinas
A potência retirada do eixo de uma turbina será menor do que a retirada do fluido.
𝑃𝑜𝑡𝐸𝑖𝑥𝑜 𝑇𝑢 𝑖𝑛𝑎 = 𝑇 𝑇
em que T é o rendimento da turbina (0 < T <1) . A potência de eixo da turbina é chamada

também de potência líquida da turbina.
A potência retirada do gerador elétrico acoplado à turbina é menor que a potência
fornecida pelo eixo da turbina:
𝑃𝑜𝑡𝐺 𝑎 𝑜 𝐸𝑙é𝑡 𝑖 𝑜 = 𝐺 𝑃𝑜𝑡𝐸𝑖𝑥𝑜 𝑇𝑢 𝑖𝑛𝑎 = 𝐺  𝑇 𝑇
em que G é o rendimento do gerador elétrico (0 < T <1)

Exemplo 6.6.1:
Sabendo que as perdas de carga entre as seções 1 e 2 na tubulação da figura são de
5m e que o diâmetro da tubulação na seção 1 é 1,0m e na seção 2 é 0,50m, pede-se: a)
Identifique se a máquina hidráulica é uma bomba ou uma turbina; b) Calcule a potência
retirada ou acrescentada ao fluido.
Máquina
Seção 1: p = 100kPa; V = 0,5m/s; z = 3m
Seção 2: p = 30kPa; z = 2m
 = 9800 N/m3.
1 2
Análise: deve ser aplicado o balanço de energia em forma de cargas, conforme a equação
6.19, para descobrir se a carga da máquina é positiva (bomba) ou negativa (turbina). Antes é
necessário calcular a velocidade na seção 2 pela continuidade.
Solução:
Admitindo regime permanente e escoamento incompressível:
𝑉 2
𝑉2 = = 𝑉 2  𝑉2 = 𝑠
2 2
A equação do balanço de energia fica:

00000 0 2 0000 22
+ + − 𝑀 = + + +
00 2 00 2
− 𝑀 =  𝑀 =  a máquina é uma Turbina
Potência extraída da água (eq. 6.15):

2
𝑃𝑜𝑡 = 𝑀 = = 𝑊
Resposta: a máquina é uma Turbina, que retira uma potência de 11,7kW da água.
Exemplo 6.6.2:
A figura mostra uma tubulação alimentada por um reservatório de grandes dimensões,
que descarrega num jato livre na cota 200,00m, com velocidade 5m/s. As perdas de carga
no percurso total da tubulação são de 25m. A carga total na entrada da tubulação (seção 1)
é de 220m. Informe se a máquina hidráulica é uma bomba ou turbina.
1 Bomba ou Turbina?
2
Análise: trata-se de um problema semelhante ao anterior, exceto que não são fornecidas as
cargas individuais na seção de entrada da tubulação.
Solução:
2
− 𝑀 = + + +
𝑀 = −
Resposta:
A máquina hidráulica é uma bomba, com altura manométrica de 6,3m.
6.7 Exercícios Propostos
6.7.1. Água escoa de um grande reservatório e descarrega em jato livre. Calcule H(m) e a
pressão relativa p (kN/m2) indicada pelo manômetro antes do bocal. Perdas desprezíveis.
6.7.2. A água se escoa na tubulação da figura. Calcule o diâmetro necessário, d, para que
as leituras manométricas sejam as mesmas.
Ex. 6.7.2
Ex. 6.7.3
6.7.3. Calcule a potência da bomba que recalca 120 l/s de água.

6.7.4. Supondo que o bocal divergente permaneça cheio, calcule a potência da bomba.
Ex. 6.7.5
Ex. 6.7.4
6.7.5. Calcule a altura h necessária para produzir uma vazão de 85 l/s e uma potência de 15
kW na turbina. Despreze as perdas.
6.7.6) Considere um recipiente fechado dotado de um tubo de entrada de água de 10mm de

diâmetro e um tubo de saída de 5mm de diâmetro. O recipiente possui um aquecedor que
gera 2.000 W de potência em regime permanente. A massa específica da água é igual a
1000 kg/m3 na entrada e saída e o calor específico é igual a 4180 J/kg. Sabendo que a água
que entra chega com uma temperatura de 25°C, e o perfil de velocidades é dado por:
r2
V (r )  1 
2.5x105 ; V em (m/s) e r em metros, pede-se:
a) V média, Fluxo de Volume e Fluxo de massa que entram no recipiente.
b) Fluxo de calor trazido ao recipiente pela água que entra.
c) Velocidade média de saída
d) Fluxo de calor que a água transporta para fora do recipiente.
e) Temperatura da água na saída.
6.7.7) O escoamento da figura ocorre em regime permanente no sentido de 1 para 2. Sabe-

se que A1 = 1m2; V1= 1m/s; z1 = 10m; p1 = 9800 Pa; V2= 0,5m/s; z2 = -10m; p2 = 450800 Pa.
As perdas de carga totais entre as seções 1 e 2 são de 5m. Pede-se: a) informar se a
máquina hidráulica “M” é uma bomba ou turbina, justificando; b) a altura manométrica da
máquina “M”; c) a potência líquida da máquina; d) a energia consumida (ou fornecida) pela
máquina, em Joules, durante um período de operação de 8 horas, admitindo um rendimento
de 85% na máquina.
6.7.8) Calcular a vazão e a potência fornecida pela bomba. As perdas de carga na tubulação
são dadas e a perda no bocal é desprezível. Dados: ATUBO = 0,1m2; ABOCAL = 0,05m2; Perdas
no tubo HTUBO = 5V2/2g; densidade do mercúrio dHG = 13,6.
Bocal
Bomba Z = 20m
2m
1m
Z = 0m
Hg
CAPÍTULO 7: BALANÇO GLOBAL DE QUANTIDADE DE

MOVIMENTO
A quantidade de movimento de um dado volume de fluido é uma grandeza extensiva.

Quando um escoamento atravessa um volume de controle, a água carrega consigo
quantidades de movimento que entram e saem. O efeito imediato mais importante disso é
que aparecem forças exercidas pelo escoamento sobre o volume de controle.
A grandeza extensiva = 𝑉 ⃗ é um vetor. A quantidade específica é  = 𝑉
⃗.
A lei física que trata da variação da quantidade de movimento num sistema é a
segunda lei de Newton.
| = =
𝑡 𝑆
Aplicando-se a relação sistema-volume de controle para a quantidade de movimento
tem-se:
⃗  = ∫ ⃗V V + ∫ ⃗V ⃗V ⃗⃗⃗⃗⃗ 7.1
O primeiro termo do segundo membro (variáveis Eulerianas) representa a taxa de

variação da quantidade de movimento no VC. Pode ocorrer porque as velocidades variam no
tempo, no interior do VC,m ou porque a quantidade de massa varia no tempo, ou pela
combinação de ambos.
O segundo termo representa o saldo de fluxos de quantidade de movimento (fluxo
que sai menos fluxo que entra).
A equação é resolvida para cada componente da força. Em problemas
bidimensionais, as componentes ficam:
= ∫ V V + ∫ V ⃗V ⃗⃗⃗⃗⃗ 7.2
= ∫ V V + ∫ V ⃗ ⃗⃗⃗⃗⃗
V 7.3
Observe nas equações 7.2 e 7.3 que o termo 𝑉⃗ ⃗⃗⃗⃗⃗ não é decomposto porque é um
escalar.
A avaliação correta do primeiro membro depende, em cada caso, do correto
isolamento do volume de controle, e da consideração de todas as forças, tanto de contato
como de campo ou de ação à distância (normalmente a força peso).
Nos itens seguintes serão desenvolvidos exemplos de aplicação a casos simples de
interesse prático.
7.1 Aplicações Elementares: Pás Defletoras
No cálculo de pás defletoras admite-se a hipótese simplificadora básica de que a pá

não muda a velocidade do jato. Isto implica em ignorar as forças tangenciais devido à
viscosidade no contato entra o fluido e a pá defletora.
A Figura 7.1 mostra um jato livre (p = patm) de água defletido no plano vertical pela pá
fixa, presa por um suporte (não representado na figura). Analise o problema para determinar
as componentes da força transferida pelo jato à pá.
Figura 7.1: Pá defletora fixa.
Inicialmente é necessário definir com clareza o volume de controle adotado, para não
cometer erros no isolamento. Foi escolhido para este problema um V.C. que inclui apenas o
jato de água. Seu isolamento aparece no lado direito da figura.
A pressão atmosférica não atua no contato jato-pá. Mas a área de ação da pressão
atmosférica sobre a superfície superior do jato tem a mesma projeção que a superfície
inferior da pá, de forma que as forças devido à pressão atmosférica vão se anular. Por isso é
mais conveniente usar pressões efetivas. Com isso, todas as forças de pressão se anulam,
com exceção da força F exercida pela pá sobre o jato de fluido.
A força F é resultado da integração de um diagrama de pressões efetivas que se
desenvolve ao longo da pá, como ilustrado na Figura 7.1. Esta distribuição de pressões é
desconhecida, mas o balanço integral não depende desta informação para calcular a força
resultante F. Esta é uma grande vantagem da técnica de balanços globais. Aplicando o
balanço, para regime permanente, temos:
Direção x:
= − = ∫ V ⃗V ⃗⃗⃗⃗⃗
Direção y:
= − = ∫ V ⃗ ⃗⃗⃗⃗⃗
V
Em que W representa o peso da água sobre a pá.
Desenvolvendo a componente x:
− = ∫ V ⃗V ⃗⃗⃗⃗⃗ + ∫ V ⃗V ⃗⃗⃗⃗⃗
2
Como as seções 1 e 2 possuem velocidades constantes, as componentes podem ser
retiradas da integral, assim como a massa específica.
− = (+V ) ∫ ⃗V ⃗⃗⃗⃗⃗ + (+V2 ) ∫2 ⃗V ⃗⃗⃗⃗⃗
− = (+V ) (−V ) + (+V2 ) (+V2 2)

Observe que há duas considerações de sinal a fazer em cada termo. A primeira
referente ao sentido da componente de velocidade, e a segunda referente à integral em área
de entrada (negativa) ou de saída (positiva).
Pela equação da continuidade, V = V2 2 = .
Pela hipótese básica de desprezar a ação do cisalhamento na pá, V1 = V2 .
Portanto, A1 = A2 .
= V( − )
Desenvolvimento semelhante no eixo vertical leva a:
− = V
EXEMPLO 7.1: A figura mostra um bocal com 0,05m2 de área, descarregando um jato livre
com Vj = 10m/s numa pá defletora horizontal. A pá divide o jato pelo meio, e o ângulo
= 30°. Calcule a força necessária para manter a pá no lugar.
Solução: Adotando um VC que inclui a pá e o jato, conforme o esquema a direita da figura,

percebemos que a força devido à pressão atmosférica se anula. A força peso não age no
problema, restrito ao plano horizontal.
=− = ∫ V ⃗V ⃗⃗⃗⃗⃗ = ∫ V ⃗V ⃗⃗⃗⃗⃗ + ∫ V ⃗V ⃗⃗⃗⃗⃗ + ∫ V ⃗V ⃗⃗⃗⃗⃗

2
− = (+ V ) ∫ ⃗V ⃗⃗⃗⃗⃗ + (− V ) ∫ ⃗V ⃗⃗⃗⃗⃗ + (− V ) ∫ ⃗V ⃗⃗⃗⃗⃗

2
(+V ) (+V )
− = (+ V )(−V ) + (− V ) + (− V )
− = (− )− ( )( )
=
O sinal positivo indica que o sentido adotado(-i) é correto.

A repetição do balanço na direção y mostrará que F y = 0, devido à simetria entre os
jatos de saída nas áreas 2 e 3.
7.2 Aplicações Elementares: Curvas em Tubulações
O cálculo de forças em curvas de tubulação com escoamento forçado é similar ao de

pás defletoras, com a adição das forças devido à pressão nas áreas de entrada e saída. A
Figura 7.2 apresenta uma curva com redução em tubulação com paredes finas.
Figura 7.2: Curva com redução.
Adota-se como volume de controle a água no interior da tubulação. O volume de

controle isolado como corpo livre aparece no esquema à direita da figura. As forças de
contato são as forças de pressão nas faces 1 e 2, e as forças decorrentes das distribuições
de pressão (pw) e de tensão de cisalhamento (w) ao longo da parede lateral da curva. A
resultante dessas distribuições desconhecidas é a força resultante que a curva aplica sobre
a água que compõe o VC. Além das forças de contato, já descritas, atua a força de campo
(peso) na direção y.
Para calcular as forças de pressão vamos considerar a pressão efetiva nas seções 1 e 2.
Com isso podemos ignorar a pressão atmosférica atuando sobre a curva.
Aplicando o balanço, para regime permanente, temos:
Direção x:
= − 2 + = ∫ V ⃗ ⃗⃗⃗⃗⃗
V
Direção y:
= − − 2 = ∫ V ⃗V ⃗⃗⃗⃗⃗
Em que W representa o peso da água contida no interior da curva.
Desenvolvendo a componente x:
− 2 2 + = ∫ V ⃗V ⃗⃗⃗⃗⃗ + ∫ V ⃗V ⃗⃗⃗⃗⃗
2
− = (+V ) ∫ ⃗V ⃗⃗⃗⃗⃗ + (+V2 ) ∫2 ⃗V ⃗⃗⃗⃗⃗
− 2 2 + = (+V ) (−V ) + (+V2 q) (+V2 2)
Lembrando que, pela continuidade,

𝑉 = 𝑉2 2 =
A equação do balanço, para o caso especial da curva da figura, fica:
− 2 2 + = (V2 q−V )
O tratamento da direção y é semelhante, levando a:
− − 2 2 = V2 q
EXEMPLO 7.2: A figura mostra uma curva vertical de 120 com redução de 300mm para
200mm. A cota média da seção 2 situa-se 1,5m acima da cota média da seção 1. A pressão
relativa na seção 1 é de 50 kPa. A vazão transportada pela tubulação é de 200 L/s. O
volume interno da curva é de 85 L. Calcule a força exercida pela água sobre a curva.
Solução:
a) Primeiramente, calcular V1 com a vazão dada e V2 pela continuidade

𝑉 = = = 𝑠
𝑉2 = = = 𝑠
2
b) Com as velocidades, calcular p2 com a equação da energia

𝑝 𝑉2 𝑝2 𝑉22
+ 𝑧 + = + 𝑧2 + + 2
2 2
𝑝2
+ + = + + +
𝑝2 = 𝑃
c) Adotar e isolar um VC para calcular o somatório de forças e aplicar o balanço de
quantidade de movimento.
O VC escolhido será a água no interior da curva. As distribuições de tensões tangenciais
e normais ao longo das paredes internas são desconhecidas, mas sua integral sobre a área
resulta na força da curva sobre a água, denotada por F na figura.
Direção x:
+ 2 2 − = ∫ V ⃗V ⃗⃗⃗⃗⃗ + ∫ V ⃗V ⃗⃗⃗⃗⃗
2
+ 2 2 − = (+V ) (−V ) + (−V2 ) (+V2 2)
= + 2 2 + (V + V2 )
𝑥 =
O sinal positivo indica que o sentido adotado (-i) é correto.
Direção z:
− 2 2 − + = ∫ V ⃗V ⃗⃗⃗⃗⃗ + ∫ V ⃗V ⃗⃗⃗⃗⃗
2
− 2 2 − + = (+ ) (−V ) + (+V2 ) (+V2 2)
= 2 2 + V + (+V2 )
= + + =
O sinal positivo indica que o sentido adotado (+k) é correto.
7.3 Aplicações Elementares: Perdas em Expansão Brusca
Quando a água em escoamento forçado encontra um aumento de diâmetro, passa

pelo fenômeno da expansão brusca. O comportamento da água é mostrado na Figura 7.3: o
núcleo do escoamento expande-se gradualmente e surge uma zona morta, com vórtices e
escoamento reverso na região próxima da expansão. Como resultado, há uma perda de
energia que pode ser quantificada com a aplicação do balanço global da quantidade de
movimento.
Adota-se como volume de controle a água entre a seção imediatamente após
expansão (1) e a seção em que a velocidade volta a ser uniforme na seção (2). O VC isolado
é mostrado na figura (b).
Figura 7.3: Expansão brusca em escoamento forçado.
Como primeira hipótese simplificadora, a força de atrito resultante da distribuição de

tensões tangenciais na parede lateral do VC é desprezada.
Como as acelerações resultantes dos vórtices na zona de estagnação são pequenas,
pode-se considerar, como segunda hipótese simplificadora, que a pressão é constante ao
longo da seção 1.
Com as simplificações adotadas, o balanço de quantidade de movimento sobre o VC
fica:
2 − 2 2 = (+V ) (−V ) + (+V2 ) (+V2 2) *
Aplicando a equação da energia entre as seções 1 e 2:

𝑝 𝑉2 𝑝2 𝑉22
+ 𝑧 + = + 𝑧2 + + 2
𝑝 𝑝2 𝑉22 𝑉2
− = − + 2
𝑉22 𝑉2
𝑝 − 𝑝2 = − + 2 **
2 2
Igualando o termo p1-p2 das equações marcadas com asterisco, pode-se obter a perda de
carga:
𝑉22 𝑉2
− 2 = V22 − V 2
= − + 2
2
𝑉22 𝑉2
2 = + ( − )
2
2
2
𝑉22 𝑉2 𝑉 ( ) + 𝑉2 ( − )
2 2
2 = + ( − )=
2
2
2
𝑉 [( ) + − ] 𝑉2 2
2 2
2 = = ( − )
2
2
𝑉2
2 = ( − )
2
O resultado mostra que as perdas aumentam com o quadrado da velocidade e diminuem

quando A1 se aproxima de A2.
7.4 Aplicações Elementares: Estruturas em Canais Abertos
O balanço global de quantidade de movimento pode ser usado para calcular os

esforços provocados pela água em estruturas como comportas e vertedores em
escoamentos livres. A Figura 7.4 ilustra um trecho de canal retangular, com largura L na
direção perpendicular à figura, com a água represada por uma comporta plana com
inclinação qualquer.
V.C. F
W
F1
Fat F2
z
x
FB
Figura 6.4: Estrutura em escoamento livre.
______ Isolando o V.C.

O volume de controle adotado é a água situada entre as seções 1 e 2. O diagrama à
direita mostra o V.C. isolado, com as forças atuantes. O somatório de forças na direção x
fica:
𝑥 = − 2 − 𝑎𝑡 − 𝑥
A força Fx é a componente horizontal da força F que a comporta exerce sobre a água.

A força F é resultante de uma distribuição desconhecida de pressões, ilustrada no diagrama.
Observe que a distribuição de pressões mostra a pressão atmosférica atuando nos dois
extremos da comporta em contato com a água. Sabemos também que F é perpendicular à
face da comporta, porque consideramos desprezíveis as forças tangenciais, originadas da
viscosidade do fluido. Da mesma forma, consideramos desprezível a força Fat causada pelas
tensões tangenciais ao longo das paredes laterais e do fundo do canal.
Nas seções 1 e 2 o diagrama de pressões é triangular (distribuição isostática), porque
as linhas de corrente são retilíneas e paralelas.
______ Força de Pressão

A força F resultante da distribuição isostática de pressões sobre uma área plana retangular é
calculada a seguir, com base na Figura 7.5..
Figura 7.5: Carregamento de pressão isostática sobre superfície plana.
= ∫ = ∫ 𝑝 = ∫0
2
=
______ Aplicando o balanço

Aplicando o balanço de quantidade de movimento (eq. 7.2) na direção x ao V.C.
adotado,
𝑦2 𝑦22
− − = ∫ V ⃗V ⃗⃗⃗⃗⃗ + ∫ V ⃗V ⃗⃗⃗⃗⃗
2 2 2
𝑦2 𝑦22
− − = (+V )(− V ) + (+V2 )(+ V2 2)
2 2
Lembrando que, pela continuidade,

2
𝑉 𝑦 = 𝑉2 𝑦2 = 𝑞 ( )
sendo q = Q/L a vazão por metro de largura do canal, a equação do balanço pode ser
resolvida para fornecer Fx.
Com a componente horizontal e conhecendo a inclinação da comporta calcula-se a
componente vertical, sabendo que a resultante é perpendicular à superfície da comporta.
EXEMPLO 7.3: A figura mostra um trecho de canal retangular com 2m de largura, com uma
comporta vertical controlando a vazão em regime permanente. Calcule a vazão e a força
exercida pela água sobre a comporta, sabendo que o nível de montante é 5 m e o de jusante
2 m. Despreze a ação das tensões de cisalhamento.
1
V1
y1 2
A1 V2
y2 A2
Análise:
A força de atrito no fundo será desprezada, bem como a força vertical da água sobre a
comporta, resultante das tensões tangenciais. A equação da energia entre as seções 1 e 2,
bem como a continuidade, serão utilizadas para determinar a vazão. Com a vazão
conhecida pode ser aplicado o balanço global da quantidade de movimento para o volume
de controle constituído pela água entre as seções 1 e 2.
Solução:
O volume de controle isolado é mostrado no esquema a seguir:
V.C.
Fx
z F1
x F2
O balanço global aplicado ao V.C. isolado acima fornece, conforme a discussão do item 7.4:
𝑦2 𝑦22
− − = (+V )(− V ) + (+V2 )(+ V2 2)
2 2
O balanço de energia entre as seções 1 e 2 fornece:

𝑝 𝑉2 𝑝2 𝑉22
+ 𝑧 + = + 𝑧2 + + 2
Desprezando as perdas e tomando as pressões e cotas na seção média

𝑉2 𝑉22
+ + = + +
A continuidade diz que:

𝑉 𝑦 = 𝑉2 𝑦2  𝑉2 = 𝑉
Portanto,
𝑉2
( 2
− ) =  𝑉 = 𝑠
2
𝑉2 = 𝑠
= 𝑠
Substituindo no balanço de quantidade de movimento:
2 2
( − ) − = ( − )
= 𝑥
O sinal positivo indica que o sentido adotado para a força exercida sobre a água é o correto
(-i).
Resposta: A força exercida pela água sobre a comporta é de 37,8kN, no sentido positivo do
eixo x (+i).
7.5 Exercícios Propostos
7.5.1. Calcule a força resultante da água sobre esta placa com orifício. O tubo tem diâmetro
de 0,3 m e o diâmetro do orifício é de 0,2 m, sendo a vena contracta de 0,16 m.
Ex. 7.5.1 Ex. 7.5.2

7.5.2. Calcule a força exercida pela água sobre esta placa com orifício. Suponha que a água
entre o orifício e a vena contracta pese 18N.
7.5.3. O projétil enche parcialmente o terminal do tubo de 0,3 m. Calcule a força necessária
para manter o projétil na posição indicada, quando a velocidade média no tubo for 6 m/s.
Ex. 6.5.3
7.5.4. Calcule a força horizontal exercida pela água sobre o cone quando a velocidade
média no tubo de 200 mm for igual a 3 m/s.
Ex. 6.5.4
7.5.5. Uma vazão de 1.820 l/min de água ocorre num tubo de 50 mm de diâmetro, o qual se
alarga para 100 mm. Se a pressão no tubo pequeno for igual a 138 kN/m2, calcule o módulo
e o sentido da força horizontal sobre o alargamento.
7.5.6. Um bocal de 100 mm é rebitado (com 6 rebites) numa dobra de uma tubulação
horizontal de 300 mm e descarrega água na atmosfera. Determine a carga de tensão em
cada rebite quando a pressão no tubo for igual a 600 kN/m2. Despreze as forças verticais.
7.5.7. Um tubo cônico divergente se encontra na horizontal, possui comprimento de 0,3 m,
sendo o diâmetro do estrangulamento igual a 75 mm; o diâmetro na saída e igual a 100 mm,
por onde ocorre uma descarga de 28,3 l/s para a atmosfera. Calcule o módulo e o sentido de
cada componente da força que a água exerce sobre o tubo.
7.5.8. Neste escoamento de água no plano horizontal, determine a força sobre as cavilhas
que prendem o conjunto curva-bocal à tubulação. O bocal descarrega em jato livre e as
perdas são desprezíveis.
Ex. 7.5.8 Ex. 7.5.9
7.5.9. Água está fluindo em regime permanente no jato livre com vazão de 50l/s. O jato é
desviado pela placa circular com 0,3m de diâmetro. Um manômetro situado no centro da
placa acusa uma pressão de 980kPa. Determine a velocidade do jato, a força exercida pela
água sobre a placa e a força a ser exercida sobre a placa para mantê-la no lugar.
7.5.10. Quando a bomba indicada na figura pára de funcionar (e não oferece nenhuma
resistência ao escoamento), a força exercida sobre a mola é igual a 672 N. Quando a bomba
está funcionando, a força cresce para 2,24 kN. Qual a potência que a bomba está
fornecendo para a água? Suponha que a superfície da água no tanque não se modifique..
Ex. 7.5.11
7.5.11. Um tanque é montado sobre um carro, conforme a figura. Considere a água que sai
do tanque através de um bocal de 0,05m 2, com velocidade de 10m/s. O nível da água no
tanque é mantido constante pela adição de água por meio de uma tubulação vertical.
Determinar a tensão no arame que mantém o carro estacionário.
7.5.12. A figura mostra um dispositivo que recebe um jato de água com 10m/s de velocidade
e 0,1m2 de área. O jato recebido é desviado e descarrega por duas fendas laterais de 0,1m2
de área cada uma, com a velocidade formando um ângulo de 60° com o plano horizontal. A
massa do dispositivo é de 5kg e 10 litros de a água ficam retidos em seu interior, em regime
permanente. Pede-se: a) calcular a velocidade média de saída em cada fenda lateral; b) a
força indicada no dinamômetro.
y
x
Ex. 7.5.13
60
Ex. 7.5.12
7.5.13. Calcule a força sobre a curva com bocal da figura, situada no plano horizontal e
descarregando com um jato livre. A área do tubo é 1m 2 , a velocidade da água no tubo é
3m/s e na saída do bocal é 9m/s. A pressão na entrada da curva pode ser calculada
sabendo que as perdas de carga entre a curva e o jato livre podem ser desprezadas.
7.5.14. Um bocal de 0,01m2 de área lança um jato no plano horizontal sobre a pá dupla da
figura, com velocidade 5m/s. Sabendo que na pá o jato se divide em dois jatos iguais,
calcule as resultantes Rx e Ry sobre a pá.
Ex. 7.5.14 o 2
60
o
45
5m/s 6m/s
y 1
y
3
o x
60
x Ex. 7.5.15
7.5.15. A figura mostra um jato livre de água de área 0,1 m2 e velocidade 6m/s, atingindo
um anteparo no plano horizontal. Ao bater no anteparo o módulo da velocidade não varia e o
jato é dividido em duas partes. Sabe-se que sai pela seção 3 duas vezes mais água do que
sai pela seção 2. Pede-se:
a) áreas A2 e A3
b) as reações sobre o anteparo
A N E X O - RESPOSTA DOS EXERCÍCIOS PROPOSTOS
CAPÍTULO 1
ITEM 1.6
1.6.1
a) Q = 0,02m3/s; b) FM = 19,9kg/s; c) FN = 4g/s (N = sólidos totais); d) FN = 3743kW.
1.6.2
O mapa mostra que as maiores concentrações de matéria orgânica devem-se, provavelmente, à
contribuição do afluente. O problema possui solução aberta. Uma das possíveis soluções é discutida a
seguir.
a) Adotei divisão da seção em subáreas com velocidades constantes.
 Seção 1 triangular, com menores velocidades devido à pequena profundidade, com base 2m e
profundidade máxima 0,5m. : A1 = 2x0,5/2 = 0,5m2;
 Seção 2 trapezoidal, entre y = 2,0 e y =3,5m, com área de transição entre as velocidades
baixas da margem e o canal principal do rio com as maiores velocidades: A2 = (0,5 + 1,0)/2 x
1,5 = 1,125m2;
 Seção 3 retangular entre y = 3,5m e y = 5,5m, com 1,5m da profundidade, considerando o

núcleo de maiores velocidades do escoamento: A3 = 2 x 1,5 = 3,0 m2.
b) Q = V1A1 + V2A2 + V3A3 = 8,0m3/s;

c) FM = 550g/s.
1.6.3
a) 70m3/s; b) 74 kg/s; c) 1.458.540m3.
1.6.4
Q = 48 L/s (0,048m3/s).
1.6.5 – trajetórias em regime permanente coincidem com as Linhas de Corrente. A seção 1 possui 6
tubos de corrente com igual vazão, visto que a velocidade é uniforme.
1. Q1 = 60m3/s, V2 = 0,533m/s, V3 = 0,5m/s;
2. Cc = 1,875g/m3, Cb = 0;
3. Cc = 0, Cb = 3,75g/m3; A alternativa de lançamento no ponto (a) é melhor do que no ponto (a’),

porque a concentração resultante é menor.
1.6.6
Devido à simetria é possível resolver para metade da seção e dobrar o resultado.
a) Q = 13,33m3/s, b) FN = 1.667g/s (N = massa de cloretos).
1.6.7
a) VM = 1m/s; b) Q = 3m3/s; c) FM = 3.000kg/s ; d) FN = 300g/s (N = massa de sal).
1.6.8
N = massa de cianetos.
Fluxo medido no rio FN = 30,6mg/s;
Fluxo máximo permitido pelo acordo FN = 9mg/s;

A empresa não está cumprindo o acordo.
ITEM 1.8
1.8.1: Q = 0,8m3/s
1.8.2: Q = 3,20m3/s
1.8.3: Q = 1,667m3/s
1.8.4: Considerei Área 1 na direção i; Área 2 na direção j;
a) Q1 = 2m3/s, saída; Q2 = 3m3/s; saída, Total Q = 5m3/s, saída
b) FM = 100g/s, saída
c) FM1= 20g/s, FM2= 30g/s, Total: 50g/s , saída
d) Q1 = 1 m3/s, Q2 = 1,5m3/s
1.8.5:
a) 7,5m3/s; saída
b) 150g/s; saída
c) 150g/s
d) 7,5 m3/s
CAPÍTULO 3 ITEM 3.5

Ex 1 : Vol = 77.760.000 m3.
Ex. 2: Vol = 23.328.008 m3.
Ex.3: Q = 2499 L/s (2,50m3/s).
Ex. 4: Vol = 31.968.000 m3.
Ex. 5: Q = 1,571m3/s; V = 0,222 m/s
Ex. 6: a) h = 3,638m; b) t = 96,4s; e) t = 42,1s
Item c)
t (s) h (m)
0 5,000
20 4,470
40 3,996
60 3,572
t (s) h (m)
0 5,000
5 4,867
10 4,738
15 4,613
20 4,491
25 4,372
30 4,256
35 4,143
40 4,033
45 3,926
50 3,822
55 3,721
Item d) 60 3,622
Item e) numericamente, interpolar linearmente entre 40 e 45segundos:

( − )
𝑡= + ( − )= 𝑠
( − )
Obs: analiticamente, obtém-se t = 42,1s
Ex. 7: t = 29,2 anos
Ex. 8: a) Q1 = Q2 = 0,030m3/s;
b) Q3 = Q4 = 0,015m3/s; Vb = 1m/s, sentido (-i) , portanto vazão de entrada.
CAPÍTULO 4
Ex. 4.5.1 P = 336.735Pa = 337kPa.
Ex. 4.5.2: Q = 0,070m3/s
Ex. 4.5.3: a) P2 = 7.400 Pa ; b) P3 = 29.600Pa ; c) P4 = 4.950 Pa ; d) Q = 0,0156 m3/s.
Ex. 4.5.4: VOrificio = 5,67m/s ; h = 1,64m ;
A vazão depende da área contraída o orifício: Q = 5,67 ASeção Contraída = 0,0073 m3/s.
Ex. 4.5.5: Q = 0,0148m3/s.
Ex. 4.5.6: a) Q = 0,0430m3/s ; b) p = 84.264 Pa; c) H = 11,82m.
Ex. 4.5.7: a) Q = 0,0964m3/s ; b) p = 44.491 Pa; c) H = 7,69m.
Ex. 4.5.8: a) H = 3,29m ; b) Q = 0,00215m3/s.
Ex. 4.5.9: A carga de pressão do ponto B é 1,51m. O nível do piezômetro C ficará 0,31 acima da cota
do ponto A.
Ex. 4.5.10: D = 0,0543m
Ex. 4.5.11: A  0,0219m2.
Ex. 4.5.12: V = 59,7m/s.
Ex. 4.5.13: Q = 0,0082m3/s ; deflexão = 0,397m.
CAPÍTULO 6
6.7.1: P = 32.279 Pa; H = 4,82m
6.7.2: d = 0,236m
6.7.3: Bomba com HM = 16,73m; Pot = 19,7kW.
6.7.4: Pot = 30.093W
6.7.5: h = 23,98m
6.7.6:
a) VM = 0,5m/s; Q = 3,96x10-6 m3/s ; FM = 3,93x10-3 kg/s
b) 410,7W
c) 2,0 m/s
d) 2410,7W
e) 146,75°C ! (a água irá evaporar, mudando as condições do problema)
6.7.7:
a) a) Bomba ; b) HM = 29,96m; c) Pot = 293.625 J/s; d) C = 2.763,5 kWh = 9,95x109 J.
6.7.8: OBS: Nas versões anteriores não foi incluída na proposta do problema a distância entre o
centro do tubo e o nível do mercúrio no lado esquerdo do manômetro (2,0m). Com esse dado, as
respostas ficam:
Q = 0,871m3/s ; b) Pot = 302,9 kW
CAPÍTULO 7
Ex. 7.5.1
Força exercida pela água sobre a placa: Fa  p = (17.723 + 0,0393 Patm ) N (sentido +x); esta resposta
depende da pressão atmosférica local.
Força resultante sobre a placa : F = 17.723 N (sentido -x); esta é a força a ser aplicada pelos vínculos
que unem o tubo à placa (parafusos, rebites, solda, etc). Esta força não depende da pressão
atmosférica local.
Ex. 7.5.2
Força exercida pela água sobre a placa: Fa  p = (322,6 + 0,01767 Patm ) N (sentido -z); esta resposta
Força resultante sobre a placa : F = 322,6 N (sentido +z); esta é a força a ser aplicada pelos vínculos
que unem o tubo à placa (parafusos, rebites, solda, etc). Esta força não depende da pressão
atmosférica local.
7.5.3
Força necessária para manter o projétil F = 6.570N (sentido -x)
7.5.4
Força exercida pela água sobre o cone: Fa  c = (1.477 + 0,031415 Patm ) N (sentido +x); esta resposta
Força resultante sobre o cone : F = 1.477 N (sentido –x); esta é a força a ser aplicada pela barra que
mantém o cone no lugar. Esta força não depende da pressão atmosférica local.
7.5.5
Força aplicada pelo vínculo sobre a expansão F = 1339N (sentido +x); a expansão comprime o tubo
menor.
7.5.6
Tensão em cada rebite T = 5655N (tração)
7.5.7
Componentes da Força exercida pela água sobre o tubo:
Componente X:
Fx a  t = (–17,34 + 0,003436 Patm ) N (sentido –x); esta resposta depende da pressão atmosférica
local.
O sentido coincide com o adotado (–x) quando Fx a  t > 0. Isto ocorre para
Patm > 5049,5Pa, o que abrange a maioria das situações comuns.
Componente Y:
Fy a  t = 18 N (sentido –y)
7.5.8
Força do conjunto curva + bocal sobre as cavilhas FVC  Cavilhas = 8.747 N (sentido –x); O sentido indica
que as cavilhas são tracionadas.
7.5.9
Velocidade do jato : V = 44,27m/s;
Força exercida pela água sobre a placa: Fa  p = (2.213,5 + 0,0707 Patm ) N (sentido +x); esta resposta
Força resultante sobre a placa : F = 2.213,5 N (sentido –x); esta é a força a ser aplicada pelos vínculos
sobre a placa para mantê-la no lugar. Esta força não depende da pressão atmosférica local.
Ex. 7.5.10
Potência da bomba Pot = 8.828W
Ex. 7.5.11
Tensão no arame T = 2.500N (tração no arame)
Ex. 7.5.12
Velocidade média de saída em cada fenda lateral V = 5m/s.
Força indicada no dinamômetro F = 14.477N (força de compressão na mola)
Ex. 7.5.13
Componente X da força sobre a curva Fx = 34.500N (sentido +x).
Componente Y da força sobre a curva: Fy = 23.383N (sentido –y).
Ex. 7.5.14
Componente x da resultante sobre a pá Rx = 375N (sentido –x)
Componente y da resultante sobre a pá Ry = 0.
Obs: a resultante sobre a pá é a força necessária para manter a pá no lugar.
Ex 7.5.15
a) Áreas: A2 = 0,03333... m2; A3 = 0,06666... m2.
b) Reações sobre o anteparo
Componente x: Rx = 4.448N (sentido –x)
Componente y: Ry = 1.551N (sentido –y)
Obs: a reação sobre o anteparo é a força necessária para manter o anteparo no lugar.

Notas Hidr Exp 2019 S2

Enviado por

Dados do documentoclique para ver informações do documento

Direitos autorais:

Formatos disponíveis

Notas Hidr Exp 2019 S2

Enviado por

Dados do documento

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

Notas Hidr Exp 2019 S2

Enviado por

Direitos autorais:

Formatos disponíveis

Hidráulica Experimental

Notas de Aula - Versão 1.8 - 2019/s2

0 A DICIPLINA HIDRÁULICA EXPERIMENTAL .................................................. 01

1 QUANTIFICAÇÃO DOS ESCOMENTOS ......................................................... 08

2 DESCRIÇÃO DOS ESCOAMENTOS ............................................................... 32

3 CONSERVAÇÃO DE GRANDEZAS (Equação da Continuidade) ..................... 57

4 TRANSFORMAÇÕES DE ENERGIA NOS ESCOAMENTOS ........................... 79

6 BALANÇO GLOBAL DE ENERGIA ................................................................... 99

7 BALANÇO GLOBAL DE QUANTIDADE DE MOVIMENTO ............................ 130

A N E X O : RESPOSTA DOS EXERCÍCIOS PROPOSTOS ................................. 144

Diferenças em relação à versão 1.7 (de 2016-s1)

Diferenças em relação à versão 1.6 (de 2015-s2)

Diferenças em relação à versão 1.5 (de 2015-s1)

Diferenças em relação à versão 1.4 (de 2014-s2)

A DICIPLINA HIDRÁULICA EXPERIMENTAL

No currículo da Engenharia Civil da FE-IS/Unesp aparece a disciplina Hidráulica

0.1 O que são e o que fazem os Fenômenos de Transporte

É fundamental para o engenheiro conhecer as leis que regem os fenômenos físicos,

0.2 Relação com as disciplinas do curso

As disciplinas de Hidráulica Experimental e Fenômenos de Transporte são as

As atividades rotineiras do engenheiro hidráulico dependem muito da quantificação

0.3. Algumas definições iniciais

Grandeza Intensiva: o valor da medida não depende da quantidade de massa

Grandeza Extensiva: o valor medido depende da quantidade de massa

Mecanismos básicos de transporte

Advecção - Mecanismo de transporte de massa, calor ou quantidade de movimento,

V Mancha de poluente no Rio

A posição do centro de massa da mancha de poluentes lançada no instante t0 no

O exemplo de advecção ilustra ainda a importância do movimento do fluido no

Mancha de poluente em água parada: V = 0

A advecção é muito mais eficiente que a difusão, que depende da movimentação

na engenharia, e também porque os processos convectivos podem ser entendidos como a

Mecanismos combinados de transporte

V Mancha de poluente no Rio

Na Figura 0.3 percebe-se que com o passar do tempo a difusão provoca o

CAPÍTULO 1: QUANTIFICAÇÃO DOS ESCOAMENTOS

1.1 Vazão ou Fluxo de Volume

É muito importante conhecer o volume de fluido que um escoamento transporta. Como os

A vazão de água transportada por um rio é fundamental em muitos problemas

A dimensão do fluxo de volume é [ M3 / T ], e as unidades mais comuns são m 3/s,

_____ Fluxo de Massa

1.2 Relação Básica entre Velocidade e Vazão

Nossa experiência cotidiana, por exemplo, com torneiras e mangueiras de jardim,

registrado num intervalo de tempo t dado:

Vol Vol 15m 3

1.3 Fluxo de Grandezas Extensivas Transportadas

O Fluxo de uma Grandeza Extensiva, FN, é a quantidade da grandeza N transferida

As unidades do Fluxo de uma grandeza são as unidades da grandeza divididas pelo

Concentração de Grandezas Extensivas

As unidades da concentração são as unidades da grandeza por metro cúbico. Por

Da definição podemos concluir que a massa específica de um meio é uma

A quantidade específica da grandeza N, para o caso da massa de uma substância A

onde T0 é uma temperatura de referência.

b) Concentração de Energia Cinética:

Fluxo da Grandeza Específica

Fluxo de N em função da sua concentração CN :

Fluxo em função da quantidade específica :

ao fluido em escoamento e os passageiros nos vagões são análogos às grandezas

Num.Pessoas Num.Pessoas Num.Vagões