Hidráulica ECIV046 EAMB029 Prof. Marllus Gustavo F. P. das Neves ctec.ufal.br/professor/mgn

HIDRÁULICA – Marllus Gustavo F. P. das Neves

HidráulicaHidráulicaECIV046 EAMB029ECIV046 EAMB029

Prof. Marllus Gustavo F. P. das Neves

www.ctec.ufal.br/professor/mgn

Universidade Federal de Alagoas

Centro de Tecnologia


Introdução à Introdução à hidráulicahidráulica


ApresentaçãoApresentação


Como será a disciplina

?

Ementa: introdução, revisão de alguns conceitos da mecânica dos fluidos, cálculo de condutos forçados, perdas lineares e localizadas, temas diversos a respeito dos condutos forçados, hidráulica dos sistemas de recalques, movimentos uniforme e gradualmente variado


Avaliação

3 Provas 2 questões práticas

AB1 maior nota das 3 provas ponderada com minitestes

AB2 média aritmética das 2 provas restantes ponderada com minitestes

Reavaliação prova que repõe menor AB

Final prova abrangendo toda a disciplina


?

Minitestes de conceitos vídeos no canal do Youtube Marllus Gustavo Neves



?

Datas das provas:

prova 1 (10/08/2017)

prova 2 (14/09/2017)

prova 3 (26/10/2017)

Reavaliação 07/11/2017

Final 14/11/2017


Prova 1: Introdução, Revisão de Mecânica dos Fluidos, Escoamento em condutos forçados até aplicações de perda de carga

Prova 2 : Escoamento em condutos forçados até aplicações de perda de carga (continuação), Máquinas hidráulicas e Análise dos sistemas de recalque

Prova 3 : Características básicas dos escoamentos livres, escoamentos uniforme e gradualmente variado


?


BAPTISTA, Márcio B. & COELHO, Márcia M. Lara P. Fundamentos de engenharia hidráulica.

Bibliografia


Vários autores. Hidráulica aplicadaColeção ABRH 8

www.abrh.org.br


Bibliografia

PORTO, Rodrigo de Melo. Hidráulica Básica


Bibliografia

AZEVEDO NETTO, J. M. Manual de Hidráulica


NEVES, Eurico Trindade. Curso de Hidráulica


A engenharia A engenharia hidráulicahidráulica


Da etimologia para o conceito atual

Escolas tradicionais hidráulica experimental e a hidrodinâmica

Desafios

Vejam a vídeo-aula 1


Hidráulica hydros + aulos

Conjunto de técnicas ligadas ao transporte de líquidos, em geral, e da água, em particular

água

condução

Conceito atual área da engenharia correspondente à aplicação dos conceitos da mecânica dos fluidos na resolução de problemas ligados à captação, armazenamento, controle, transporte e uso da água

E para chegar a este conceito?


Mecânica dos fluidos

Líquidos e gases

Líquidos (água)

Estados: sólido, líquido e gasoso


Revisão de alguns Revisão de alguns conceitosconceitos


2.1. Propriedades 2.1. Propriedades Físicas dos FluidosFísicas dos Fluidos


•Forças de massa ou de corpo de superfície

•Esforços Pressão Tensão

•Massa específica e peso específico•Compressibilidade•Viscosidade

Dinâmica Cinemática Fluidos Newtonianos

•Pressão de vapor



ρgγ

Nosso curso = 9.810 N/m3


Para a transformação Kgf N multiplica-se por 9,81

ρgγ


2.2. Classificação 2.2. Classificação dos escoamentosdos escoamentos


• Pressão reinante forçado Livre canais

•Trajetória das partículas Laminar turbulento

•variação no tempo Permanentes transitórios (não-permanentes)

•Direção, módulo e sentido do vetor velocidade Uniforme e uniforme por seção Variado: gradualmente ou bruscamente

•No de coordenadas do campo de velocidade Unidimensional Bidimensional Tridimensional



forçado

livre


2

max Rr

1uu unidimensional

unidimensional e uniforme em cada seção

bidimensional


Equações Equações fundamentais do fundamentais do

escoamentoescoamento


Lei N Nosso curso

Conservação da massa M 1 Continuidade

2ª lei de NewtonQuantidade de movimento

1ª lei da termodinâmica

E e Bernoulli

P

V

SCVC

dAnVηρηρdtdt

dN

N por unidade de massa

vazão em massa através do elemento de área dA

Elemento de massa contido no VC


Equação da Equação da ContinuidadeContinuidade


Lei N Nosso curso

Conservação da massa M 1 Continuidade

2ª lei de NewtonQuantidade de movimento


E e Bernoulli

P

V

SCVC

dAnVηρηρdtdt

dN Equação geral


Lei N

Conservação da massa

M 1

2ª lei de Newton


E e

P

V

SCVC

dAnVρρdt

0

SCVC

dAnVηρηρdtdt

dN

A massa é constante em VC


Supondo escoamento permanente

Fluxo líquido de vazão em massa na superfície de controle é nulo

SCVC

dAnVρρdt

0

0

0SC

dAnVρ


Supondo escoamento permanente

vazão em massa que entra = vazão em massa que sai

TM

m kg/s

0SC

dAnVρ

SC

dAnVρm

vazão em massa

No caso mais simples:


Para o escoamento incompressível constante; VC indeformável forma e tamanho fixos

Vazão em volume (Q) que entra no VC = Qsai

A

dAnVQ

TL

Q3

m3/s, l/s, ft3/s...

Vazão em volume chamada de Vazão


O caso mais simples

1 22V

1V

Esc. permanente incompressível e uniforme em cada seção

21 A222

A111

SC

dAnVdAnVρdAnVρ0

2n

1n


1 22V

1V

0dAnVdAnV21 A

222A

111

0AnVAnV 222111

11AV- 22AV

QAVAV 2211

2n

1n

0dAnVdAnV21 A

222A

111

uniforme por seção


O caso de uma bifurcação escoamento permanente incompressível e uniforme em cada seção

Q1,V1,A1

Q2,V2,A2

Q3,V3,A3

n1

n2

n3

SCVC

dAnVρρdt

0

0


3

1i Aiii

SC i

dAnVρdAnVρ0

Constante na seção

iiiA

iii AnVdAnVi

integral

3

1iiii AnVρ0

V1

n1

x

y

Seção 1


Seção 2

Seção 3

332211 VAVA VA0

x

y

V2

n2

x

yV3

n3

Q1,V1,A1

Q 2,V 2

,A 2

Q3 ,V

3 ,A3

321 QQQ


Equação da Equação da Quantidade de Quantidade de

movimentomovimento


Lei N Conservação da massa

M 1

2ª lei de Newton


E eP

V

(s)massa(s)s ηρdηdmN

SCVC

dAnVηρηρdtdt

dN


SCVC

dAnVρVρdVtdt

Pd

dtPd

R

SCVC

dAnVρVρdVt

R

Resultante das forças no VC


SCVC

SB dAnVρVρdVt

FF

Forças de massa

Forças de superfície

Equação vetorial decompor nas componentes

SCVC

x dAnVρu ρudR

tNa direção x

wv,u,V


Analogamente nas demais

SCVC

y dAnVρv ρvdR

t

SCVC

z dAnVρw ρwdR

t

nVPrestar atenção no

sinal•verifica-se o sinal do produto escalar;•depois o sinal de cada componente de velocidade


Para o caso mais simples Q constante, uniforme por seção, incompressível

xy

1 22V

1V

2n

1n

SC

x dAnVρuR

Na direção x

SC

y dAnVρvR

Na direção y = 0, pois v = 0


21 A2222

A1111x dAnVVdAnVVρR

xy

1 22V

1V

2n

1n

,0,0VV 11 ,0,0VV 22

wv,u,V


21 A2222


1 22V

1V

2n

1n

21 A2222


22221111x AnVVAnVVρR


1 22V

1V

2n

1n

22221111x AnVVAnVVρR

222111x AVVAVV-ρR

Q Q

12x VVρQR


O caso de uma bifurcação

x

y

Q1

Q2

Q3

n1

n2

n3

Precisamos dos ângulos decompor vetores

Regime permanente, incompressível e uniforme em cada seção


3

1i Aiiii

SC i

dAnVVρdAnVVρR

Constante na seção

iiiiA

iiii AnVVdAnVVi

integral


Direção x

333222

111SC

AVcosVAVcosV

AV-VρdAnVVρ

113322SC

QV-cosVQcosVQρdAnVVρ

O termo da direita

Direção y (faça como exercício) senVQsenVQρdAnVVρ 3322

SC


V1

n1

x

y

Seção 1

Seção 2

Seção 3

x

y

V2

n2

x

yV3

n3


Resumindo

senVQsenVQρR 3322y

113322x QV-cosVQcosVQρR

Os lados esquerdos, Rx e Ry, podem ser decompostos, conforme as forças consideradas


n

1i isinalicosiViQρsR α

De forma geral direção s qualquer

•verifica-se o sinal do produto escalar;•depois o sinal de cada componente de velocidade


n

1i isinalicosiViQρsR α

x

yV3

n3

Direção x:Produto escalar +Sentido da componente +Resultado +Direção y:Produto escalar +Sentido da componente -Resultado -


Equação de BernoulliEquação de Bernoulli


Lei N

Conservação da massa

M 1

2ª lei de Newton


E e

P

V

(s)massa(s)s ηρdηdmN

SCVC

dAnVηρηρdtdt

dN


SCVC

dAnVeρd eρtδt

δWδtδQ

Taxa de trabalho realizado pelo VC

zg2V

ue2

interna específica

cinética específica

potencial específic

a +

Transferência da taxa de energia através da SC por causa de diferença de temperatura


SCVC

dAnVeρd eρtδt

δWδtδQ

Devido a eixos de rotação (bombas, turbina)

Devido à ação do cisalhamento agindo em um contorno em movimento (correia móvel)

δtδW

δtδW

cis

eixo

SC

dAnVp Resultante da força devida à

pressão movendo na SC trabalho de escoamento


SCVC

cis

SC

eixo

dAnVeρd eρt

δtδW

dAnVpδt

δWδtδQ

zg2V

ue2

Colocando na integral e separando os termos de energia interna ...

0

Escoamento permanente

fazendo Weixo = 0 e Wcis = 0

00


SC

2

dAnVρgz2V

upδtδQ

SC

2

SC

dAnVρgz2V

pdAnVuρδtδQ

Termos que representam formas de energia não-utilizáveis

perdas


perdasdAnVρgz2V

p0SC

2

Perdas soma de todos os termos representando formas de energia não-utilizáveisTomando agora um caso simples (2 seções)

perdasdAnVρgz2V

p

dAnVρgz2V

p0

A12222

22

2

A11111

21

1


perdasAnVρgz2V

pAnVρgz2V

p0 2222

22

21111

21

1

perdasAρVgz2V

p)Aρ(-Vgz2V

p0 222

22

2111

21

1

Q Q

perdasρQgz2V

pρQgz2V

p 2

22

21

21

1


perdasρQgz2V

pρQgz2V

p 2

22

21

21

1

Dividindo tudo pela vazão mássica e chamando o termo de perdas de H

m m

122

222

1

211 Hz

2gVp

z2gVp

γγ

Relação entre velocidade, pressão e elevação

Hz2gV

γp 2

carga (energia) total por unidade

de peso


122

222

1

211 Hz

2gVp

z2gVp

γγ

Relação entre velocidade, pressão e elevação

Hz2gV

γp 2

carga (energia) total por unidade

de pesoV é a velocidade ao longo de uma LC ou a velocidade média (idealização de perfil uniforme)

MUNSON, B. R., YOUNG, D. F., OKIISHI, T. H. Fundamentos da mecânica dos fluidos

velocidade média


12

22

22

21

11 ΔH

2g

Vz

γ

p

2g

Vz

γ

p

H1

H2


Significado dos termos

12

22

22

21

11 ΔH

2g

Vz

γ

p

2g

Vz

γ

p

z2gV

γp

H2

Energia ou carga de pressão

Carga de posição (energia potencial)Energia ou carga

cinética


Linha de carga efetiva ou linha piezométrica

(LP)

cotas piezométricas (CP) ou cargas piezométricas

Linha de energia (LE)

zγp

CP


z2gV

γ

pH

2

02gV2

1

ca)(manométri 0γp

γp 41


Considerara não uniformidade do perfil de velocidade

Várias trajetórias

Levar em conta este fato coeficientes de não uniformidade


Coeficientes de não uniformidade

Coeficiente de Coriolis

c1

c2

E

Eα

3c1 ρAU

21

E

A

3c2 dAρV

21

E

1α

AU

dAV

EE

3A

3

c1

c2

fator de correção de energia

1,05 ≥ ≥ 1,15

Em correntes muito irregulares 1,10 ≥ ≥ 2,00


Fazendo-se o mesmo com a QM

1AU

dAV

AρU

dAρV

β2

A

2

2A

2

é o fator de correção da QM ou coeficiente de Boussinesq

Escoamentos:turbulentos em condutos forçados > 1,10laminares em condutos forçados > 1,33turbulentos livres 1,02 ≥ ≥ 1,10


Equação Equação fundamental da fundamental da

hidrostáticahidrostática


0 g ρ p -

A equação abaixo estabelece o campo de pressão em um fluido estático

Força de pressão por unidade de volume em um ponto

Força de massa por unidade de volume em um ponto

Variação de Pressão em um Fluido Estático

Escolhendo um eixo de coordenadas no qual o vetor gravidade esteja alinhado com o eixo z...

z gz = -gkg -


Observando as restrições

fluido estático

a gravidade é a única força de massaeixo z vertical

fluido incompressível

hidrostática

Sendo po no nível de referência zo integrando a equação geral

p – po = -ρg(z-zo) = ρg(zo-z)


Se a superfície do corpo fluido for tomada como referência z - zo = - h

p - po = ρghz

h

Equação da hidrostática

pbar é a leitura barométrica local

pbar pabs= pbar+pm

pm pm é a pressão manométrica

zero absoluto de pressão

ou pressão atmosférica local

Níveis de referência para pressão


pbar pab

s

pm

patm padrão

1 atm

101 kPa

760 mmHg

14,696 psi

2.116 lbf/ft2

22,92 in mercúrio

33,94 ft água


h

Elemento fluido imerso em água com a superfície exposta à atmosfera

p - po = ρghDa equação da hidrostática

pm

patm

pm = γhA pressão exercida pelo fluido é a manométrica


ManometriaManometria


Método de medição de pressões a partir de deslocamentos produzidos numa coluna contendo um ou mais fluidos

piezômetro

Manômetro em U

Manômetro diferencial

Manômetro inclinado,...


Manômetro tipo Bourdon

A medida de pressão é relativa o exterior do tubo está sujeito à pressão atmosférica


A pressão em B é a soma da pressão em A com a pressão da coluna h1

A pressão em B’ é a mesma que em B, pois estão no mesmo nível em um mesmo fluido


Cálculo da pressão em BpB - pA = ρ1gh1

pB = γ1h1 + pA

ou

Por outro ladopB = γ2h2 + pc


Isto resulta em

pA = patm + γ2h2 - γ1h1

Se desprezarmos patm, calcularemos somente pressões manométricas


Surgem então as regras práticas

1) Quaisquer 2 ptos na mesma elevação, num trecho contínuo do mesmo líquido, estão à mesma pressão

2) A pressão aumenta à medida que se caminha líquido, para baixoLembrar da variação de

pressão ao mergulhar numa piscina


Forças hidrostáticas Forças hidrostáticas sobre superfícies sobre superfícies

submersassubmersas


kdAAd

Superfícies planas


Não há tensões de cisalhamento força hidrostática é normal ao elemento de superfície dxdydA

Força no elemento dA kpdAFd

Força resultante na área

kFrdAkp-FA

R

kdA Fd

Ou seja A

RR pdAFF


jyixr ccc

A força resultante tem um ponto de aplicação centro de pressão ou empuxo

ARc ypdAFy A

Rc xpdAFxComo achar?

Para um fluido estático e incompressível:p = p0 + gh h = ysen

yh

A0R ydAρgsenθApF


A última integral é o momento de 1ª ordem da superfície em relação ao eixo dos x

A

0R ydAρgsenθApF

ycg é a coordenada y do centro de gravidade (CG). Logo

AyydA cgA

Chamando hcg = ycgsen ApF cgR


módulo da força resultante em uma superfície plana submersa = produto da área pela pressão unitária que atua em seu centro de gravidade

Como achar o ponto de aplicação (xc,yc)?

Tomando a pressão manométrica (p0=patm) p=gh=gysen

ARc ypdAFy A

2cgc dAyAyy

A última integral é o momento de 2ª ordem da superfície em relação ao eixo dos x Ix


Ou sejaAy

Iy

cg

xc

Do teorema dos eixos paralelos e designando Icg o momento de 2ª ordem em relação ao eixo baricêntrico ou do CG

Ay

Iyy

cg

xcgcgc AyII 2

xcgx cg

Para xc, o resultado é semelhante, usando Ixycg, que é o produto de inércia em relação ao par de eixos xy que passa pelo CG


Ay

Ixx

cg

xycgcgc

Axy xydAI

AyxII cgcgxycgxy

Ay

Ix

cg

xyc

Resumindo superfície plana submersa com a superfície livre à pressão atmosférica

ApF cgR

Ay

Iyy

cg

xcgcgc Ay

Ixx

cg

xycgcgc


Superfícies curvas caso mais geral

FR continua sendo normal à superfície, contudo a direção dos elementos de força varia

Determinar as componentes de FR


npdAFd

x

y

zn

j

ik zyx cosθ,cosθ,cosθn

inpdAiFddF xR

1,0,0cosθ,cosθ,cosθpdAdF zyxxR

xxxR pdApdAcosθdF

xx A

xA

RxRx pdAdFF

Da mesma forma FRy e FRz


No plano zy

xcgAxA

xRx AγhhdAγFx

No plano zx

ycgAyRy AγhF

xAx

z

FRx

0

hcgxxcgAxRx AγhF


γFRz

h

Componente z γdγhdAγF

zAzRz


Escoamento em Escoamento em condutos forçadoscondutos forçados


Introdução à perda Introdução à perda de cargade carga


Forçadolivre


Fatores governantes

Canal gravidade

Conduto forçado gravidade em menor grau, gradiente de pressão principal p1 – p2


12

22

22

21

11 ΔH

2g

Vz

γ

p

2g

Vz

γ

p

H1

H2

Perdas de carga (ou de energia)

1221 ΔHHH



Mecanismos atrito com paredes e turbulência

Grau de agitaçãoNo escoamento laminar é um fenômeno microscópico

No escoamento turbulento é um fenômeno macroscópico

νhh UD

μ

ρUDRe

Material asperezasD

Mecanismos e abrangência


Experimento de ReynoldsLaminar x

turbulento

νhh UD

μ

ρUDRe baixa U tem que ser baixa

para o escoamento ser laminar


u

u'

u'uuA

Matematicamente a turbulência se traduz nas flutuações u’ em torno da média

u

uA


O que são as asperezas?

D

altura média da rugosidade

/D rugosidade relativa



Abrangência distribuída ou localizada

Atrito e turbulência (macroscópica) possuem peso grande

turbulência (geometria) influencia mais


Distribuída


http://www.deg.ufla.br/site/_adm/upload/file/HfLocalizada2007.pdf


Tensões nos regimes Tensões nos regimes laminar e turbulentolaminar e turbulento


Características fundamentais estabelecer fórmulas de perdas de carga

-Perfil de velocidade-Distribuição de tensão de cisalhamento

Escoamento laminar analiticamente

Escoamento turbulento analítica e experimentalmente


Escoamento laminar Escoamento laminar plenamente desenvolvidoplenamente desenvolvido


Perda de carga contínua tensões de cisalhamento

Perfil de velocidade

Tipo de regime de escoamento

laminar turbulento

Hagen-Poiseulle

Perfil parabólico Descoberto de forma analítica


Trecho de comprimento L e queda de pressão p

L128pD

Q4

μ

π

4LpD

p

2gV

DL

fp 2

γ

fator de atrito f = 64/Re

Tubo horizontal

Hp

γ D

L4H p

γ

τ


Escoamento turbulento Escoamento turbulento plenamente desenvolvidoplenamente desenvolvido


Perda de carga contínua tensões de cisalhamento

Perfil de velocidade

Tipo de regime de escoamento

laminar turbulento

Perfil não é mais parabólico

Descoberto com a ajuda de experimentos


2gV

DL

fp 2

γ

f fator de atritoContinua valendo

generalizado

y

y = R – r


Tensão de cisalhamento

dydu

Primeiramente semelhante à lei de viscosidade descoberta por Newton

lamturb

Experimentos: maiores que as calculadas dessa forma nova proposição

dyud

tturb dyud lam


Viscosidade turbulenta

dyud

dyud

tvuturb ''

dyud

tlamturb

Não é propriedade do fluido função das variáveis do escoamento médio

Joseph Boussinesq (1877)

Ainda não era um modelo prático modelar em função de variáveis do escoamento médio


Comprimento de mistura lm

2

2

dyud

mturb l

Também função das variáveis do escoamento médio, da distância da parede, etc. na vizinhança da parede é quase proporcional à distância dessa

Tentativa de solucionar matematicamente o modelo da viscosidade turbulenta

Ludwig Prandtl (início do século XIX)


Perfil de velocidade no Perfil de velocidade no regime turbulentoregime turbulento


CamadasAlém da dificuldade para obter (y), descobriu-se que o escoamento pode ser considerado em camadas

Externa ou turbulenta

Parede

Superposição ou transição ou inercial

amortecedora

viscosa ou laminar ou linear ou de parede << 1% D

(y)u

y


Camada viscosa

*

*

yuuu

Lei da parede

5yu

0 *

Localização

y

p*u

onde

Diminui com o aumento de Re

Tensão de cisalhamento na parede

Dominam os efeitos viscososPerfil quase linear e grande gradiente de velocidadeTensão de cisalhamento laminar

*u5


Camada de amortecimento

305 *yu

Localização

Nenhum perfil de velocidade é exato nessa camada


Camada de Superposição

B*yu

lnk1

uu

*

Lei Logarítmica

30*yu

Localização

Efeitos turbulentos ainda não dominantesVelocidade proporcional ao logaritmo de y

5*yu

ln2,5uu

*

k ≈ 0,4 e B ≈

0,5


Camada turbulenta

A lei logarítmica representa satisfatoriamente bem a região acima da camada de superposição

LUDV Lei universal de distribuição de velocidade

5*yu

ln2,5uu

*


yuu

*

5lny2,5uu

*

Camada turbulenta


*yu

y

Camada turbulenta


Parede

Superposição

amortecedora

viscosa ou laminar ou linear ou de parede << 1% D

(y)u

y

5lny2,5uu

*


Camada turbulenta

Cálculo de B com o requisito de que a velocidade máxima num tubo ocorre no centro (r = 0)

B*yu

lnk1

uu

*

y

y = R – rmaxu

r-RR

ln2,5u

u-u

*

max

Déficit de velocidade


r-RR

ln2,5u

u-u

*

max

Lei do déficit de velocidade

Camada turbulenta em tubos

y

y = R – rmaxu


Harpa de NikuradseHarpa de Nikuradse


Esc. hidraulicamente lisos (HL)Escoamentos de transição

(HT) Esc. hidraulicamente rugosos (HR)Numa tubulação pode ocorrer qualquer

um

LUDV

Teorias de

turbulência

Leis de resistênc

ia específic

as

Harpa de Nikuradse

Equação

Universal



Mecanismos atrito com paredes e turbulência

Grau de agitaçãoNo escoamento laminar é um fenômeno microscópico

No escoamento turbulento é um fenômeno macroscópico

νhh UD

μ

ρUDRe

Material asperezasD

Mecanismos e abrangência


Asperezas e camada viscosa

A relação entre e vai ditar o escoamento

Fator de atrito f função de e função de que Re f é função de Re e

Escoamento laminar 2g

VDL

fp 2

γ

onde f = 64/Re



Parede

Superposição

amortecedora viscosa

(y)u

y

Asperezas e camada viscosa

Escoamentos HLEscoamentos de transição (HT)Escoamentos HR


2gV

DL

fp 2

γ

Dε

,DL

Re,Funçãof

2gV

DL

fH2

Generalizado via

análise dimensional

Equação universal ou de Darcy-Weisbach

Do laminar para o turbulento


Utilizou tubos lisos cuja parede interna esteve revestida com grãos de areia esféricos

Harpa de Nikuradse


Harpa de Nikuradse

Fator de atrito

Rugosidade relativa

UD

ReN° de Reynolds

2gV

DL

fH2


fórmula para laminar: f = 64/Re

I – Re < 2.300: escoamento laminar

Harpa de Nikuradse

Regiões


II – 2.300 < Re < 4.000

região crítica f não caracterizado

Harpa de Nikuradse

Regiões


fórmula para lisos: f = F(Re)

III – curva dos tubos lisos: f = F(Re)

Harpa de Nikuradse

Regiões


IV – turbulento de transição

RegiõesHarpa de Nikuradse


fórmula para rugosos: f = F(Re,)

V – turbulento rugoso

f=F(/D)para um tubocom /Dconstante,f é constante

Harpa de Nikuradse Regiões


Desprendimento da curva de tubos lisos com aumento de Re

O aumento da turbulência provoca diminuição de expõe as asperezas da parede

HT HR


Diagrama de Moody e leis Diagrama de Moody e leis de resistênciade resistência


Esc. hidraulicamente lisos (HL)Escoamentos de transição

(HT) Esc. hidraulicamente rugosos (HR)

LUDV

Teorias de

turbulência

Leis de resistênc

ia específic

as

Harpa de Nikuradse

Equação

Universal

2gV

DL

fH2

Dε

,DL

Re,Funçãof


Fórmulas racionaisFórmulas racionais

2gV

DL

fH2


1939 Ciryl F. Colebrook (1910-1997): combinou dados para o escoamento de transição e turbulento, tanto para tubos lisos e rugosos

fRe

2,513,71D

ε2log

f

110

Indicada para a faixa de transição entre os escoamentos lisos e rugosos, no intervalo

198D/ε

fRe14,14

Equação de Colebrook


1942 Hunter Rouse (1906-1996): confirmou a equação de Colebrook e produziu um gráfico

Diagrama de Moody

fReRe,Funçãof

1944 Lewis F. Moody (1880-1953): Recriou o diagrama de Rouse

https://pt.wikipedia.org/wiki/Diagrama_de_Moody


Laminar Lisos

Diagrama de Moody

Transição

Turbulento Rugoso


1976 Swamee-Jain: fórmula explícita

2

0,9Re

5,743,7D

εlog

0,25f

10-6 ≤ /D ≤ 10-2 e 5.103 ≤ Re ≤108

Exata até 2% do diagrama de Moody


1983 S. E. Haaland fórmula explícita equação com resultados dentro de 2% daqueles obtidos pela equação de Colebrook

Re6,9

3,7Dε

1,8logf

11,11

10

As equações explícitas anteriores podem ser usadas como uma boa primeira estimativa de métodos iterativos no uso das equações implícitas


1993 Swamee fórmula explícita equação geral válida para escoamento laminar, turbulento liso, de transição e turbulento rugoso

O gráfico obtido concorda bem com o diagrama de Moody


Fórmulas empíricasFórmulas empíricas

m

n

DQ

J


fórmula de Blasius Curva limite dos tubos HL faixa 3.000 < Re < 105

0,25Re

0,3164f Ajusta-se bem aos resultados

para tubos lisos, como de PVC

Fórmula para o escoamento laminar a partir de Hagen-Poiseulle, lei de Newton e universal


fórmula de Blasius

Laminar


Laminar

fórmula de Blasius


42 DQ

ρgπ64μ

2gU

ρD32μ

J Laminar

5

22

D

Q0,0827f

2gU

Df

J Turbulento rugoso

Fórmula universal

4,75

1,752

0,25 D

Q0,00078f

D2gU

Re

0,316J Turbulento liso

Fórmula de Blasius

Sob esta inspiração, surgem as fórmulas empíricas


Uma das mais utilizadas é a de Hazen-Williams

4,871,85

1,85

DC

Q10,65J

J(m/m), Q(m3/s), D(m)

C coeficiente de rugosidade = F(natureza, estado das paredes)

Recomendada, preliminarmente para•escoamento turbulento de transição•água a 20 oC não considerar o efeito viscoso•em geral D ≥ 4” (0,1m)•aplicação em redes de distribuição de água, adutoras e sistemas de recalque


Material C Material C

Aço corrugado (chapa ondulada)

60 Aço com juntas lock-bar, tubos novos

130

Aço com juntas lock-bar, em serviço

90 Aço galvanizado 125

Aço rebitado, tubos novos

110 Aço rebitado, em uso 85

Aço soldado, tubos novos

130 Aço soldado, em uso 90

Aço soldado com revestimento especial

130 Cobre 130

Concreto, bom acabamento

130 Concreto, acabamento comum

120


Material C Material C

Ferro fundido novo 130 Ferro fundido 15-20 anos de uso

100

Ferro fundido usado 90 Ferro fundido revestido de cimento

130

Madeiras em aduelas 120 Tubos extrudados PVC 150


Fórmula de Fair-Whipple-Hsiao

Projetos de instalações prediais de água fria recomendada pela ABNT para PVC e aço galvanizado, em instalações hidráulico sanitáriasJ(m/m), D(m) e Q(m3/s)

4,88

1,88

DQ

0,002021J

4,75

1,75

DQ

0,0008695J

Aço galvanizado novo conduzindo água fria

PVC rígido conduzindo água fria


Onde se aplicam as Onde se aplicam as fórmulas de perdas fórmulas de perdas

contínuascontínuas

m

n

DQ

J


Escoamentos completamente desenvolvidos

Trecho 1-2 perfil não uniforme camada limite

Seção 1 perfil uniforme

Seção 2 perfil constante final de le

Trecho 2–3 esc. melhor descrito


Região de entrada e escoamento planamente desenvolvido

Trecho 3-4 esc. complexo como na entradaTrecho 4-5 ainda influência da curva

Trecho 5–6 semelhante ao trecho 2-3


Escoamento em um sistema típico


Perda de carga unitária x linha de energia

ângulo de assentamento da tubulação inclinação da LE

2tg1Jcos

JLcos

Htg

2tg1Jtg

Inclinação da LE > J, a não ser que = 0

Para < 15º diferença desprezível tg = 1,04.J


Perda de carga Perda de carga singularsingular


é função das mudanças de forma, de diâmetro, de direção do escoamento ou de combinações destas

são importantes em condutos curtos

Mudanças alargamentos ou estreitamentos, curvas, bifurcações, equipamentos diversos na canalização (válvulas e outras estruturas).As modificações no escoamento por causa das mudna elementos são as chamadas singularidadesNa prática depende somente da geometria, a não ser nos casos de transições graduais. Para Re > 104, é possível ignorar o efeito da viscosidade


A perda de carga singular é avaliada comparando-se o antes e o depois da singularidade

Sem o efeito da singularidade (regime estabelecido)

Hipótese de escoamento unidimensional válida


Zonas com características fortemente tridimensionais

Aumento das tensões de

cisalhamento


Aceleração e aumento de intensidade de turbulência


Redemoinhos às custas da energia


A energia se transforma em calor


O processo de perda é contínuoMas tratamos de maneira discreta


Coeficientes de perda de Coeficientes de perda de carga singularcarga singular


Em geral, a perda de carga singular é expressa da seguinte maneira K coeficiente adimensional,

determinado experimentalmente para Re > 105 e analiticamente para um pequeno número de casos

U velocidade média de referência. Em geral, nas peças em que há mudanças de diâmetro, é tomada na seção de menor diâmetro (velocidade média maior)


Mudanças de diâmetroMudanças de diâmetro


Mudanças bruscas alargamento brusco, contração brusca, entradas e saídas de canalizaçãoMudanças graduais estreitamentos graduais (convergentes) e alargamentos graduais (difusores ou divergentes);


Experimentos: pAB = p1 em média

VAB ~ V1

AAB ~ A2

Para o alargamento brusco

Ocorre a desaceleração do fluido no trecho curto


Aplicando a equação da QM entre as seções AB e 2, desprezando o atrito entre o fluido e a parede da tubulação 12 VVρQR 122221 VVρQApAp

12221 VVρVpp

Aplicando a equação de Bernoulli, levando-se em conta somente a perda singular

12

222

211 Δh

2gV

γp

2gV

γp

12

21

22

21 γΔh2g

V-Vγpp


Igualando 2g

V-VΔh

221

12

A partir da equação da continuidade

2g

VK

A

A1

2g

VΔh

21

1

22

112

D1/D2 = 0 equivale a uma saída livre em um reservatório


No caso de contração brusca Contração do jato Logo após expansão

Despreza-se a perda de carga entre 1 e 0Reduz-se ao

anterior

h no fluxo acelerado 1-0 << h no fluxo desacelerado 0-2


Entre as seções 0 e 2

2g

V-VΔh

220

02

V0 é a velocidade média do jato na seção contraídaO valor de A0 não é conhecido a priori na maior parte dos casos, é obtido em estudos experimentaisDefinindo Cc como coeficiente de contração

2

0c A

AC

2g

VK

2g

V1

C1

Δh22

22

c02


D2/D1 = 0 ou A2/A1 = 0 equivale a uma entrada de reservatório não reentrante e não ajustada


Entradas de canalizaçãoDepende da forma geométrica e do ângulo de inclinação em relação à parede de entrada

O mais comum é a aresta viva 90º lateral ou fundo dos reservatóriosEntrada normal

No caso de aresta viva K=0,5


Bordos Reentrantes Para Re > 104, K=F(/D, b/D)

Ajuste cônico de bordos K=F(,l/D)


l/D > 0,6 aumento de H (distribuída)

Bordos arredondados h é da mesma ordem do caso de bordos cônicos, com a vantagem de precisar de menor comprimento

K menor


Bordos arredondados

r raio de curvatura da superfície de concordância


saídas de canalização

1)Descarga ao ar livre K=1,0

2) Para dentro de um reservatórioSe não houver recuperação de energia cinética

comDifusores esta será perdida K=1,0

Relembrando...


Estreitamentos graduais Minimizar as perdas na transição ou simplesmente para manter o escoamento mais homogêneo Podem ser cônicas ou curvilíneas

h = F(A2/A1 ou D22/D1

2 e L)

Simplicidade de execução

Melhor homogeneização


Alargamentos graduais (difusores) não só o ângulo de abertura é importante

Formas, comprimento do trecho reto antes do difusor, Re, e relação entre áreas

• > 60º ocorrerá o descolamento da camada limite •Até 6º e L < 4(A2/A1), não ocorrerá


Em geral, empregam-se ângulos fracos e guias correntes internas

Minimizam o comprimento

Para ângulos menores que 40º K também é composto de duas partes


Mudanças de direçãoMudanças de direção


Mudanças de direção

bruscas graduais

Para 0º ≤ ≤ 180º K = C1C2 e C1 e C2 dependem de

K depende de R/D e Re


Equipamentos diversosEquipamentos diversos


Equipamentos diversos1.Válvula de gaveta;2.Válvula de pressão;3.Válvula de retenção (posição

horizontal);4.Válvula de pé;5.Crivo


Válvula de gaveta Válvula em que o

elemento vedante é constituído de um

disco circular (ou retangular) queinterrompe a passagem do

escoamento,movimentando-se verticalmenteh = f(X, geometria

interna)

X abertura do disco


Válvula de pressão Fechar o fluxopor completo e frequentemente sistema fechado mais eficiente, mascom mais perda de carga

Sistema de fechamento disco metálico com anel de material vedante ou não anel sob a ação de uma haste é pressionado sobre o corpo da válvulaTipos•haste a 90º com a entrada e com a saída: tipo globo•0º com a entrada 90º com a saída: angulares ou tipo ângulo•45º com a entrada e com a saída:tipo Y


Empregadas geralmente na saída de condutos em instalações domiciliares para o controle de vazão do sistema


Válvula de retenção Evitar o retorno do fluxo quando a bomba pára o seu movimento a do tipo portinhola é a mais usada para diâmetros médios (50mm<D<300mm)


Válvula de pé Base de tubulações de recalque, quando a bomba não estiver afogada, para que a canalização não se esvazie quando a bomba está parada

Crivo Proteger contra entrada de em estações de recalque, antes da válvula de pé geralmente metálico, composto por um de cesto com furos


Tabela geralTabela geral


Diante de tantas fórmulas e tabelas costumam-se utilizar tabelas mais abrangentes


Comprimento Comprimento equivalente de uma equivalente de uma

singularidadesingularidade


A perda de carga localizada pode ser calculada pelo método dos comprimentos equivalentes ou comprimentos virtuais

Le comprimento de um tubo de diâmetro e rugosidade tal que proporciona a mesma perda de carga da singularidade considerada

Impondo a igualdade

ΔHΔh

2gU

DL

f2gU

k2

e2

DL

fk e


O comprimento obtido pela soma do comprimento do conduto L com os comprimentos equivalentes Le a cada singularidade é chamado comprimento virtual Lv

Valores de Le adaptados da NBR 5626/82 são mostrados a seguir


Aço galvanizado ou ferro fundido (m)


PVC rígido ou cobre (m)


Condutos equivalentesCondutos equivalentes


Um conduto é equivalente a outro(s) conduto(s) quando transporta a mesma vazão sob a mesma perda de carga

Através deste conceito, condutos em série ou em paralelo podem ser transformados, para efeito de cálculo, em um conduto simples

Condutos em série: condutos de características distintas, mas colocados na mesma linha e ligados pelas extremidades Conduzem a mesma vazão e a perda de carga total é a soma das perdas em cada um dos condutos individuais


Normalmente adotam-se De e e, calculando Le


Condutos em paralelo: extremidades de montante e de jusante reunidas num mesmo ponto, mas a vazão é distribuída entre eles entre as duas extremidades

sujeitos à mesma perda de carga, uma vez que as diferenças entre CP de montante e jusante são as mesmas


Equação geral

Então...

As perdas de carga localizadas podem ser representadas pelos comprimentos virtuais Lv

Hidráulica ECIV046 EAMB029 Prof. Marllus Gustavo F. P. das Neves ctec.ufal.br/professor/mgn

Documents

Transcript of Hidráulica ECIV046 EAMB029 Prof. Marllus Gustavo F. P. das Neves ctec.ufal.br/professor/mgn