Enunciado

PA11.8 – Uma bomba-turbina é um dispositivo reversível que usa um reservatório para gerar energia durante o dia e depois bombeia a água de volta para o reservatório durante a noite. Vamos redefinir o Problema P6.62 como uma bomba-turbina. Recorde-se que Δz=36,6 m e que a água escoa através de 610 m de tubo de ferro fundido de 152 mm de diâmetro. Para simplificar, suponha que a bomba opera no PMR (92%) com H*Bomba=61 m e que a turbina opera no PMR (89%) com H*Turbina=30,5 m. Despreze as perdas localizadas. Calcule (a) a potência de eixo, em watts, necessária para a bomba; e (b) a potência, em watts, gerada pela turbina. Para outras leituras técnicas, consulte o URL www.usbr.gov/pmts/hydraulics_lab/pubs/EM/EM39.pdf.

Passo 1

A)

Potência requerida para a bomba, Pe,Bomba.

É a potência necessária que deve-se fornecer a bomba para ela funcionar dentro dos parâmetros operacionais estabelecidos. Ela pode ser determinada pela expressão do rendimento da bomba, cuja equação é

η B o m b a = P h P e , B o m b a

Onde,

-ηBomba é o rendimento da bomba;

-Ph é a potência hidráulica, potência fornecida ao fluido;

-Pe,Bomba é a potência de eixo que é requerida pela bomba;

Vamos discutir um pouco sobre o rendimento...

O rendimento da bomba é a quantidade de potência hidráulica que o equipamento está fornecendo para o fluido, PH, dividido pela quantidade de potência de eixo fornecida para a bomba.

O rendimento da bomba pode ser entendido também da seguinte forma: é a quantidade percentual de energia que o fluido recebeu da bomba, da quantidade total que foi fornecida para a bomba em potência de eixo.

Passo 2

Vamos agora revisar Ph

Ph como foi dito, é a potência hidráulica, e esta representa a quantidade total de energia que a bomba transmitiu para o fluido.

Essa variável, lembrando então, é função das propriedades do fluido e da vazão do fluido. A expressão para ela é a seguinte

P h = ρ g Q H B o m b a

Onde,

-ρ é a massa especifica (densidade) do fluido (água), que vale 998kg/m3 para a água;

-g é constante da gravidade, g=9,81 m/s2;

-Q é a vazão do fluido;

-HBomba é a altura de carga da bomba;

Passo 3

Substituindo, a expressão para a eficiência da bomba agora é

η B o m b a = ρ g Q H B o m b a P e , B o m b a

Isolando Pe,Bomba

P e , B o m b a = ρ g Q H B o m b a η B o m b a                     E q 1

Passo 4

Precisamos da vazão, Q

Podemos achar a vazão pela velocidade de escoamento do fluido, que podemos achar pela equação de Darcy-Weisbach da perda de carga

h L T = h p + h a c

h L T = f L T D T V 2 2 g + V 2 2 g ∑ i K i           E q 2

Onde,

hLT é a perda de carga total do sistema, devido ao atrito do fluido com a tubulação e ao atrito com os acessórios localizados;

-hac é a perda de carga localizada;

-hp é a perda de carga por atrito na tubulação;

-f é o fator de atrito;

-LT é o comprimento total da tubulação;

-DT é o diâmetro interno da tubulação;

-V é a velocidade do fluido na tubulação;

-Ki é uma constante que é especifica para cada tipo de acessório localizado na tubulação;

-g é a constante da gravidade, g=9,81 m/2.

Foi dito para desprezarmos as perdas localizadas, então

h L T = h p = f L T D T V 2 2 g       E q 3

Passo 5

Para expressarmos em termos de vazão, podemos usar a seguinte relação entre vazão e velocidade de escoamento

Q = A V

Onde,

-Q é a vazão volumétrica;

-A é a área da seção transversal do tubo;

-V é a velocidade do fluido.

A área de seção transversal é função do diâmetro da tubulação

A = π 4 D T 2

Assim,

Q = π 4 D T 2 V

Isolando V

V = 4 Q π D T 2

Substituindo na Eq3

h p = f L T D T 16 Q 2 2 π 2 g D T 4

h p = f 8 π 2 L T Q 2 g D T 5           E q 4

Passo 6

Fazendo um balanço de energia mecânica (Eq. de Bernoulli), conseguimos determinar quem é hp.

Relembrando da equação de Bernoulli aplicada para processos com máquinas de fluxo é

H B o m b a = ∆ p ρ g + ∆ V 2 2 g + ∆ z + h L T

Onde,

- HBomba é a altura de carga da bomba, que é a energia que a bomba deve fornecer ao fluido;

-Δp é a variação de pressão do sistema;

-ΔV2 é a variação de energia cinética do sistema;

-Δz é a variação de energia potencial do sistema;

-hLT é a perda de carga total do sistema, devido ao atrito do fluido com a tubulação e ao atrito com os acessórios;

-ρ é a densidade do fluido;

-g é a constante da gravidade, g=9,81 m/s2.

Para hLT

h L T = h p + h a c

Onde,

-hp e hac são a perda de carga devido ao atrito do fluido com a tubulação e a perda de carga devido ao atrito com os acessórios, respectivamente;

Passo 7

Aplicando o balanço de energia mecânica (Eq de Bernoulli) entre o reservatório de baixo (1) e o reservatório de cima (2), da figura P6.62, lá em cima no enunciado

  p 1 ρ g + V 1 2 2 g + z 1 + H B o m b a = p 2 ρ g + V 2 2 2 g + z 2 + h p                 E q 5

Onde,

-p1 e p2 são a pressão no reservatório de baixo e no de cima, respectivamente;

-V1 e V2 são a velocidade do fluido no reservatório de baixo e no de cima, respectivamente;

-z1 e z2 são a altura no nível no reservatório de baixo e no de cima, respectivamente;

-hp é a perda de carga por atrito do fluido na tubulação;

-ρ é a densidade da água, ρ=998 kg/m3;

-g é a constante da gravidade, g=9,81 m/s2.

Na Eq5 não foi adicionado o termo hac, pois foi pedido para desprezar as perdas localizadas.

Isolando HBomba

H B o m b a = p 2 - p 1 p g + V 2 2 - V 1 2 2 g + z 2 - z 1 + h p

Passo 8

Vamos agora analisar a figura P6.62, lá em cima para ver quem é quem da equação de Bernoulli.

Analisando o esquema da figura P6.62, vemos que

p 1 = p 2 = p a t m

Onde, patm é a pressão atmosférica que você pode considerar para as duas pressões, pois os dois reservatórios estão abertos para a atmosfera.

Também vemos que as velocidades em 1 e 2 são referentes as velocidades de nível de fluido nos reservatórios, e como estes são de grande dimensão (D muito alto), então as velocidades são muito pequenas. Assim,

V 2 ≈ 0   e   V 1 ≈ 0

Já para as alturas z1 e z2, temos pela figura que

z 2 - z 1 = 36,6   m

Passo 9

A Eq. da altura de carga da bomba, HBomba agora é

H B o m b a = z 2 - z 1 + h p

Isolando hp

h p = H B o m b a - ( z 2 - z 1 )

Calculando

h p = 61 m - 36,6 m = 24,4   m

Passo 10

Voltando na Eq4

h p = f 8 π 2 L T Q 2 g D T 5 = 24,4   m

Isolando Q

Q 2 = 24,4 π 2 8 f g D T 5 L T

Q = 24,4 π 2 8 f g D T 5 L T 1 2           E q 6

Para encontrarmos Q, precisaremos de f, o fator de atrito

Passo 11

O fator de atrito depende de Reynolds (Re) e da rugosidade relativa, ε/D

  f = f R e , ε D

Pelo problema 6.62, o tubo é de ferro fundido e para este material a rugosidade, ε, vale

ε = 0,00026 m

ε D = 0,00026 m 0,152   m = 0,00171

A equação de Reynolds é

R e = 4 ρ Q π μ D T

Onde,

-Re é o número de Reynolds;

-ρ é a densidade do fluido (água), ρ=998 kg/m3 para a água;

-DT é o diâmetro da tubulação, pelo problema 6.62, DT=152mm=0,152m;

-Q é a vazão do fluido;

-μ é a viscosidade dinâmica do fluido (água), μ=0,001 kg/m.s para a água;

Calculando Re:

R e = 4 998 k g m 3 π 0,001 k g m   s ( 0,152   m ) Q

R e = 8359822,8 Q

Passo 12

Você pode ver que Re precisa de Q, mas não temos Q, que é o que queremos achar. Então, Re depende de Q, f depende de Re e Q depende de f, não podemos dessa forma determinar diretamente o valor de nenhum nem outro. O que temos que fazer é usar um método iterativo para encontrar a solução desse problema. Lembrando que um método iterativo você precisa dar um chute inicial para alguma das variáveis e a partir desse chute inicial você começa a calcular com as equações, pegando os valores calculados e atualizando o chute e recalculando novamente até que seja atingido uma convergência. Essa convergência você pode verificar da seguinte maneira: pegue o valor anterior calculado e diminua pelo novo valor calculado, estipule uma tolerância, ai quando essa diferença for menor que a tolerância você chegou em uma convergência.

Passo 13

  • Algoritmo Resolução pelo método iterativo:

I – chute inicial para f (você pode olhar no diagrama de Moody a faixa de valores de f);

II – calcula Q pela equação

Q = 24,4 π 2 8 f g D T 5 L T 1 2

Q = ( 24,4   m ) π 2 8 f 9,81   m / s 2 0,152   m 5 ( 610   m ) 1 2

Q = 0,006267 f 1 2

III – calcula Re pela equação

R e = 8359822,8 Q

IV – determina novo valor de f usando diagrama de Moody ou a Eq. de Churchill

f = 8 8 R e 12 + A + B - 1,5 1 12

A = 2,457 ln ⁡ 7 R e 0,9 + 0,27 ε D - 1 16 ; B = 37530 R e 16

ε D = 0,00171

V – calcula novo valor de Q pela mesma equação;

Q = 0,006267 f 1 2

VI – calcular a diferença percentual (erro) entre o novo valor de Q com o valor anterior pela expressão

e r r o   p e r c e n t u a l = 100 Q A t u a l - Q A n t e r i o r Q A n t e r i o r

VII – verificar se erro é menor que a tolerância (estipular um valor da tolerância, um valor muito pequeno, tipo uns 10-5%);

VIII – se o erro for maior que a tolerância, você volta a etapa III para uma nova iteração; (você vai continuar as iterações até que o erro seja menor que a tolerância)

Se você fizer este método iterativo à mão, vai levar algum trabalho. A melhor opção é usar um programa computacional, como o MATLAB, ou Python, ou até mesmo o EXCEL.

Para resolver eu usei o EXCEL, este é o mais simples. Abaixo o print do método no EXCEL.

Veja pela imagem acima que na iteração de nº 5 já temos um erro percentual de 0,000078%, e também veja que o valor de Q, praticamente não se altera mais a partir dessa iteração. Portanto, houve uma convergência na iteração 2.

Os valores convergidos foram:

f = 0,0198 ; R e = 3,7 x 10 5 ;   Q = 0,0445 m 3 s

Passo 14

Agora que temos a vazão, podemos calcular a potência requerida pela bomba pela Eq1

P e , B o m b a = ρ g Q H B o m b a η B o m b a

Substituindo os valores

P e , B o m b a = 998 k g m 3 9,81 m s 2 0,0445 m 3 s ( 61   m ) 0,92 = 28901   W

Passo 15

B)

Analogamente para a turbina, você vai fazer literalmente os mesmos passo a passo, como feito para a bomba.

O que terá de diferente é o cálculo da potência de eixo gerada pela turbina, Pe,Turbina

Calculando Pe,Turbina pela expressão de rendimento da turbina.

O rendimento da turbina é a quantidade de potência que o equipamento gerou, Pe, dividido pela quantidade de potência disponível que a turbina tem para receber do fluido, PH, cuja expressão é

η T u r b i n a = P e , T u r b i n a P H

Onde,

-ηturbina é o rendimento da turbina;

-Pe é a potência de eixo de saída da turbina;

-PH é a potência hidráulica disponível, é a potência que a turbina recebe do fluido;

O rendimento da turbina pode ser entendido também da seguinte forma: é a quantidade percentual de energia que a turbina recebeu do fluido, ou melhor, a quantidade percentual de energia que a turbina removeu do fluido, que escoa por dentro da turbina, para gerar a potência de eixo.

Passo 16

Revisando PH

PH como foi dito, é a potência hidráulica disponível, e esta representa a quantidade total de energia disponível que o fluido possui para fornecer à turbina.

Essa variável, lembrando então, é função das propriedades do fluido e da vazão do fluido. A expressão para ela é a seguinte

P H = ρ g Q H T u r b i n a

Onde,

-ρ é a massa especifica (densidade) do fluido;

-g é constante da gravidade, g=9,81 m/s2;

-Q é a vazão do fluido;

-HTurbina é a altura de queda da turbina, a energia que a turbina remove do fluido;

Substituindo, a expressão para a eficiência agora é

η T u r b i n a = P e , T u r b i n a ρ g Q H T u r b i n a

Isolando Pe

P e , T u r b i n a = η T u r b i n a ρ g Q H T u r b i n a           E q 7

Passo 17

Agora é o mesmo processo do que foi feito na questão A) para determinar Q

Então, analogamente

Pela Eq. de Bernoulli

H T u r b i n a = z 2 - z 1 - h p

Veja que na equação de Bernoulli para a bomba chegamos na expressão

H B o m b a = z 2 - z 1 + h p

Para a turbina o termo de hp está com sinal negativo, isso porque a turbina remove energia do fluido para gerar potência de eixo, então essa perda de carga (hp) , isto é, perda de energia do fluido é uma energia não é aproveitada pela turbina. E no caso da bomba, o sinal fica positivo, pois é uma energia que a bomba deve fornecer ao fluido para suprir a perda de carga, a perda de energia que o fluido sofreu.

Continuando, então. Isolando hp

h p = z 2 - z 1 - H T u r b i n a

Calculando

h p = 36,6 m - 30,5 m = 6,1   m

Pela Eq. de Darcy-Weisbach

h p = f 8 π 2 L T Q 2 g D T 5 = 6,1

Vamos isolar Q e chegar em

Q = 6 , 1 π 2 8 f g D T 5 L T 1 2

O fator de atrito é por Reynolds e rugosidade

f = f R e , ε D

ε = 0,00026 m

ε D = 0,00026 m 0,152   m = 0,00171

R e = 4 ρ Q π μ D T

Passo 18

Agora vamos para as iterações, também de forma análoga da parte (A)

Algoritmo Resolução pelo método iterativo:

I – chute inicial para f (você pode olhar no diagrama de Moody a faixa de valores de f);

II – calcula Q pela equação

Q = ( 6,1   m ) π 2 8 f 9,81   m / s 2 0,152   m 5 ( 610   m ) 1 2

Q = 0,0031336 f 1 2

III – calcula Re pela equação

R e = 4 998 k g m 3 π 0,001 k g m   s ( 0,152   m ) Q

R e = 8359822,8 Q

IV – determina novo valor de f usando diagrama de Moody ou a Eq. de Churchill

f = 8 8 R e 12 + A + B - 1,5 1 12

A = 2,457 ln ⁡ 7 R e 0,9 + 0,27 ε D - 1 16 ; B = 37530 R e 16

ε D = 0,00171

V – calcula novo valor de Q pela mesma equação;

Q = 0,0031336 f 1 2

VI – calcular a diferença percentual (erro) entre o novo valor de Q com o valor anterior pela expressão

e r r o   p e r c e n t u a l = 100 Q A t u a l - Q A n t e r i o r Q A n t e r i o r

VII – verificar se erro é menor que a tolerância (estipular um valor da tolerância, um valor muito pequeno, tipo uns 10-5%);

VIII – se o erro for maior que a tolerância, você volta a etapa III para uma nova iteração; (você vai continuar as iterações até que o erro seja menor que a tolerância)

Pela imagem acima

f = 0,0206 ; R e = 1,8 x 10 5 ;   Q = 0,0218 m 3 s

Passo 19

Agora que temos a vazão, podemos calcular a potência requerida pela bomba pela Eq1

P e , T u r b i n a = η T u r b i n a ρ g Q H T u r b i n a

Substituindo os valores

P e , T u r b i n a = 0,89 998 k g m 3 9,81 m s 2 0,0218 m 3 s ( 30,5   m ) = 5800   W

Resposta

A)

P e , B o m b a = 28901   W

B) P e , T u r b i n a = 5800   W

Ver Outros Exercícios desse livro

Exercícios de Livros Relacionados

A bomba para água modelo A-12251 da Worthington Corp., operando no ponto de máximo rendimento, produz 16,2 m de altura de carga a 3.550 rpm, consome 1,1 hp de potência de eixo a 3.200 rpm, e fornece 3
Deseja-se usar uma bomba de 45 hp que produza uma altura de carga de 61 m quando estiver operando no PMR com gasolina a 20 °C a 1.200 rpm. Utilizando as correlações nas Figuras P11.49 e P11.50, determ
Considere que se deseja fornecer 330 l/s de gás propano (peso molecular = 44,06) a 1 atm e 20 °C com uma elevação de pressão de 200 mm de H20 em um único estágio. Determine o tamanho e a rotação aprop
O ventilador real discutido no Problema 11.61 tinha um diâmetro de 6,1 m [Referência 20, p. 339]. Qual seria o diâmetro apropriado para a família de bombas da Figura 11.14 para fornecer 500 m3/s de ar
Um ventilador de ventilação de uma mina, operando a 295 rpm, fornece 500 m3/s de ar ao nível do mar com uma elevação de pressão de 1.100 Pa. Esse ventilador é axial, centrífugo ou misto? Calcule seu d