Perda De Carga Em Tubulações

EQUAÇÕES EXPLÍCITAS PARA O FATOR DE ATRITO DA EQUAÇÃO DE DARCY-WEISBACH CARLOS ROGÉRIO DE MELLO1 DANIEL FURTADO FERREIRA2 JACINTO DE ASSUNÇÃO CARVALHO3 RESUMO - O objetivo deste trabalho é modelar equações matemáticas que permitam o cálculo do fator de atrito da equação de Darcy-Weisbach, na forma explícita em f, devido ao fato de os modelos mais usados possuírem soluções matemáticas complexas, em razão da forma implícita de f nos mesmos. Para isto, trabalhou-se com situações de fluxo turbulento de transição e liso, estimando-se modelos específicos para ambas situações, limites de velocidade de 0,5 a 2,0 m/s, viscosidades da água às temperaturas de 20oC a 30oC, diâmetros de 15 a 600 mm, rugosidades relativas de 0,000088 a 0,01. As diversas situações foram estudadas, respectivamente, com as equações de Colebrook-White e Von

Kármán, utilizando-se um programa computacional específico de matemática (MAPLE V Versão 3 para Windows) para obter soluções de f, uma vez que essas equações são tidas como as mais precisas. Para a determinação dos modelos, avaliaram-se diversas variáveis, verificando a significância estatística e estimando seus respectivos parâmetros, através de regressão linear múltipla com o Sistema SAS para Windows, utilizando procedimentos de seleção de variáveis (Backward). Foi possível obter modelos cujas estimativas de parâmetros estatísticos foram de alta qualidade, proporcionando soluções próximas aos valores exatos, caracterizando modelos de elevada precisão.

TERMOS PARA INDEXAÇÃO: Perda de carga, equações matemáticas, fator de atrito.

EXPLICIT EQUATIONS FOR THE DARCY-WEISBACH EQUATION FRICTION COEFFICIENT ABSTRACT - The purpose of this work is to fit explicit mathematic equations for calculating the friction coefficient from Darcy-Weisbach equation, in f explicit form, based on fact that the most used models contain complex mathematics solutions due to their implicit form. For these studies, it was used situations of smooth pipe and transition zone turbulent flow, with the objective to estimate specific models for both situations, considering flow velocity between 0.5 and 2.0 m/s, kinematic viscosity of water at the temperatures between 20 and 30oC, diameters between 15 and 600 mm and relative roughness between

0.000088 and 0.01. Different situations were studied using respectively the equations of Colebrook-White and Von Kármán with the software MAPLE V Release 3 for Windows, once these equations are considered most exact. Several variables were studied in order to determine the models, by estimating their respective parameter and verifying the statistical significance using the software SAS for Windows, with the procedure Backward for selection of variables. With these studies, it was possible to obtain models with high quality statistical parameters, providing solutions close to the exact values, which characterizes models of high precision.

INDEX TERMS: Charge loss, mathematic equations, friction coefficient. INTRODUÇÃO O cálculo de perdas de carga em situações que envolvam fluxo de água em tubulações é fonte constante de estudos, uma vez que esse fator refere-se à perda de energia provocada por atritos que ocorrem entre a água e as paredes das tubulações, como conseqüência da intera-

ção entre viscosidade e rugosidade, sendo refletida nos custos variáveis da instalação (Kamand, 1988). Existem várias formulações desenvolvidas por alguns pesquisadores para o cálculo de perda de carga, destacando-se as equações de Flamant, Hazen-Willians e Darcy-Weisbach. Segundo Gomes (1997), a equação desenvolvida por este último é a mais utilizada no meio

1. Eng. Agrícola, Pós-Graduando em Irrigação e Drenagem, Departamento de Engenharia da UNIVERSIDADE FEDERAL DE LAVRAS (UFLA). Enviar correspondência para D.F.F. (DEX-UFLA). 2. Eng. Agrônomo, Dr., Professor Adjunto do Departamento de Ciências Exatas/ UFLA, tel.: (035) 829-1369. 3. Eng. Agrícola, Dr., Professor Adjunto do Departamento de Engenharia/UFLA, tel.: (035) 829-1482.

366 científico, garantindo maior ajuste dos dados à realidade física. Esta equação, também conhecida como Fórmula Universal de Perda de Carga, possui a seguinte expressão: hf = f ×

L D

×

V

2

(1)

2g

f

em que hf é a perda de carga (m), f é o fator de atrito (adimensional), L é o comprimento da tubulação (m), D é o diâmetro (m), V é a velocidade (m/s) e g é a constante gravitacional (m/s2). O fator de atrito representa a principal dificuldade ao cálculo da perda de carga, pois as formulações propostas são do tipo implícitas, com f em ambos os membros da equação, sendo de difícil resolução, mesmo com o uso de calculadoras programáveis e programas de informática específicos de matemática. Bernardo (1989) propõe o uso de resolução gráfica através de um diagrama conhecido por Diagrama de Moody. Tal procedimento é trabalhoso tanto quanto impreciso, não produzindo resultados satisfatórios. De acordo com Azevedo Netto e Alvarez (1991), o fluxo de água em uma tubulação pode ser classificado em turbulento, laminar ou crítico (transitório). Essa característica é determinada através do cálculo de um parâmetro adimensional denominado Número de Reynolds (Re): V×D (2) Re =

Como a expressão (3) depende do valor de f, considera-se, inicialmente, a situação de fluxo turbulento de transição, trabalhando com a equação de Colebrook - White.  E (4) 1 2,51 

Estimando-se f e observando se a relação 3 proporciona um valor entre 14 e 200, intervalo esse que define o regime turbulento de transição, o valor encontrado de f por (4) estará correto (Yanagi Júnior, 1995). Para valores menores que 14, f deverá ser reestimado pela equação de Von Kármán, que constitui no regime de fluxo turbulento liso, no qual a rugosidade relativa não influencia significativamente na perda de carga, e apenas o Número de Reynolds é relevante (Azevedo Netto e Alvarez, 1991), significando que o fator E/D não tem importância física:  2,51  = −2 × LOG 10    Re× f  f

1

E Re× × f ×  D

(5)

Nos casos de valores acima de 200, regime de fluxo turbulento rugoso, utiliza-se a equação de Nikuradse (Azevedo Netto e Alvarez,1991), a qual é explícita em f: 1

ν

f

em que Re é o número de Reynolds (adimensional) , ν é a viscosidade cinemática (m2/s), V e D definidos na equação 1. Conforme o valor encontrado, trabalha-se com formulações diferentes do fator de atrito, exceção feita ao regime transitório, que não possui formulação. Sendo assim, para um Re menor que 2000, tem-se o regime laminar; entre 2000 e 4000, o transitório, e maior que 4000, o turbulento. A maioria dos escoamentos em tubulações de sistemas de irrigação pressurizados ocorre em regime turbulento (Gomes, 1997). O fluxo turbulento pode ser classificado em três regimes diferentes, o que leva a equações diferentes do fator f (Yanagi Junior, 1995). Trabalha-se com os regimes: turbulento liso, misto e rugoso, o que pode ser caracterizado a partir do seguinte cálculo:

= −2 × LOG 10  +   3,7 × D Re× f 

 E  = −2 × LOG10    3,7 × D 

(6)

Cossolosso e Satto (1996) citam algumas equações explícitas para o fator f, destacando-se a equação de Chen Ning Hsing , que proporciona resultados satisfatórios somente dentro do fluxo turbulento de transição, não gerando precisão desejada para a turbulência lisa, 1 f

  11098 ,  1 5,0452  E 5,8506   E = −2 × LOG10 − × LOG10 ×  + 0,8981    37065 × D  28257 D , Re , ( Re)   

(7)

em que os parâmetros f, E, D e Re estão definidos pelas expressões anteriores. Para o fluxo turbulento liso, destaca-se a equação de Blasius, que é matematicamente simples. Entretanto, há o inconveniente de que essa equação é específica para cada tipo de material, necessitando-se estimar os parâmetros a e b da mesma.

(3)

em que E é rugosidade absoluta (m), Re é o número de Reynolds (adimensional), f é o fator de atrito (adimensional) e D é o diâmetro (m).

f =a×

1 Re b

(8)

Re e f estão definidos nas equações anteriores; a e b são constantes a serem determinadas.

Ciênc. e Agrotec., Lavras, v.23, n.2, p.365-371, abr./jun., 1999

367 Segundo Von Bernuth (1990), a equação de Blasius é altamente precisa para tubos de plástico lisos com pequenos diâmetros e em situações em que o Número de Reynolds for menor que 100.000, o que, em projetos de irrigação, correspondem a um diâmetro máximo de 80mm, gerando, portanto, dificuldade técnica em certos dimensionamentos. Devido às dificuldades matemáticas geradas pelas equações implícitas que exprimem o fator f, este trabalho propõe o desenvolvimento de modelos matemáticos explícitos desse fator, tendo-se como base as equações de Von Kármán e Colebrook-White. MATERIAL E MÉTODOS Para a realização deste trabalho, foram utilizados dados de diâmetro, rugosidade relativa, velocidade, viscosidade cinemática e Número de Reynolds, com valores próximos da realidade prática dos projetos de irrigação pressurizada e sistematizados em pequenos intervalos de variação. Com isto, trabalhou-se com diâmetros de 15 a 600 mm, viscosidades referentes às temperaturas de 20oC a 30oC (Von Bernuth, 1990), velocidades de 0,5 a 2,0 m/s (valores econômicos para dimensionamentos, segundo Bernardo,1989), proporcionando número de Reynolds maiores que 4000, sempre caracterizando fluxo turbulento, e rugosidades relativas (E/D) de 8,8x10-5 a 10-2. Foram consideradas as situações de fluxo turbulento misto e liso, uma vez que as equações implícitas em f (equações 4 e 5) são as que governam essas situações. A fim de obter bons ajustes para os modelos, procurou-se trabalhar em cada situação de diâmetro, velocidade e temperatura, com um intervalo de rugosidades relativas que abrangesse os extremos do fluxo turbulento de transição. Assim, foi investigado, em cada situação, os valores de E/D e de f que proporcionassem resultados da equação 3 igual a 200 e igual a 14, gerando também, alguns dados intermediários, cobrindo, portanto, todas as possibilidades de fluxo turbulento de transição. Essas investigações foram realizadas com a equação de Colebrook - White, que é tida por vários pesquisadores como sendo a mais precisa dentre vários modelos, sendo aceita universalmente (Gomes, 1997). Ela foi solucionada com o Sistema MAPLE V, Versão 3 para Windows, gerando soluções com várias casas decimais, tendo-se adotado, neste caso, 5 casas decimais. Do total de valores combinados dentro dos limites estabelecidos, foram gerados 593 valores para f, os quais foram encontrados com alta precisão. A equação de Von Kármán foi solucionada por meio do referido programa computacional, para a situa-

ção de fluxo turbulento liso, para cada diâmetro, situações de velocidades entre 0,5 e 2,0m/s e viscosidades às temperaturas de 20oC a 30oC, gerando, assim, 418 situações. Após essa etapa, estudou-se o modelo mais adequado para cada situação, utilizando-se o Sistema SAS para Windows (System Analysis Statistical, 1985), tendo as equações de Colebrook-White e Von Kármán como base para determinação do modelo. Trabalhou-se também com regressão linear múltipla dos dados obtidos, através do procedimento Proc Reg, utilizando-se o método Backward para seleção de variáveis. Foram testadas a significância do modelo e dos respectivos parâmetros da regressão (Ho: bi=0), ou seja, o teste em que a variável em questão não trouxe contribuição significativa para a variação do fator de atrito, para a situação de fluxo turbulento de transição e para o fluxo turbulento liso. Como parâmetros estatísticos para averiguação da qualidade do modelo, foram considerados: o coeficiente de determinação ajustado (R2), o coeficiente de variação (CV), significância ao nível de 0,01 de probabilidade para as estimativas dos parâmetros e a média dos desvios em relação aos dados originais. Estes desvios foram calculados seguindo metodologia de Kamand (1988): ∆f ( %) =

fE − f f

× 100

(8)

sendo fE estimados pelos modelos e f oriundos das equações de Colebrook-White e Von Kármán. RESULTADOS E DISCUSSÃO As estimativas dos parâmetros para o modelo relativo ao fluxo turbulento de transição estão apresentados na Tabela 1, com suas respectivas variáveis e significância ao nível de 0,01 de probabilidade. Observase que todas as variáveis testadas foram altamente significativas, demonstrando a sua importância na descrição do modelo. Foram obtidos, coeficiente de determinação ajustado (R2) de 99,85%, coeficiente de variação (CV) de 1,298% e média dos desvios entre os valores estimados e os originais de 0,955%, sendo os valores máximo e mínimo observados, respectivamente, de 6,53% e 0,00059%. Como foram consideradas 5 casas decimais, a maioria dos erros ocorreu na quarta ou quinta casa decimal, e como todos os parâmetros estatísticos são de alta qualidade matemática, o modelo descrito abaixo é capaz de prever, com alto grau de confiabilidade, o fator de atrito, desde que todas as variáveis estejam dentro dos intervalos especificados na metodologia. O modelo estimado está apresentado a seguir.

Ciênc. e Agrotec., Lavras, v.23, n.2, p.365-371, abr./jun., 1999

368  E  1,1098 1   f = 0,011096 − 0,000108 × V − 4, 39543 × 10 −10 × Re + 0,178996 × LOG 10   + Re   D    2

   E  1,1098 1   E  46,161642    − 0,157153 × LOG10  + 0, 004001 ×  LOG10   + +    3, 7 × D  Re Re D        E + 0, 032358 ×   × V  D

TABELA 1 - Estimativas dos parâmetros e sua significância segundo as variáveis do modelo para situação de fluxo turbulento de transição. Variável

Estimativas

Prob > |t| para Ho:bi=0

Intercepto

0,011096

0,0003

V

-0,000108

0,0010

Re

-4,39543x10-10

0,0001

-2xLOG10[(E/D)1,1098+1/Re]

-0,089498

0,0001

{LOG10[(E/D)1,1098 + 1/Re]}2

0,004001

0,0001

-LOG10[E/(3,7xD)]

0,157153

0,0001

1/Re

46,161642

0,0001

(E/D)xV

0,032358

0,0001

Analisando-se a Fig. 1, que caracteriza as situações de diâmetros pequeno, médio e grande, para uma velocidade em torno de 1,5m/s, nota-se uma pequena variação entre os valores de f estimados a partir do modelo desenvolvido para fluxo turbulento de transição e a equação de Colebrook-White, o que confirma o esperado, devido ao alto valor de R2 encontrado. Cossolosso e Satto (1996) citam vários modelos explícitos para essa situação, os quais possuem desvios superiores aos obtidos neste trabalho, como as equações de Barr e Swamee e Jain, que obtiveram desvios, respectivamente, de 1% e 6,22%. A equação de Chen Ning Hsing (equação 7) apresentou um desvio médio de 0,30%, inferior ao obtido com o modelo.

Contudo, as diferenças em relação à equação de Colebrook-White ocorreram, em ambas, apenas na quarta e/ou quinta casa decimal, mostrando que, na prática, ambos os modelos proporcionam resultados equivalentes. Para os casos de fluxo turbulento liso, a Tabela 2 apresenta os parâmetros estimados com suas respectivas significâncias, tendo-se obtido um R2 de 99,98%, um CV de 0,473% e média dos desvios de 0,347%, com valores máximo e mínimo, respectivamente, de 3,14% e 0,0016%, mostrando que o modelo foi satisfatório quanto à sua qualidade em predizer os valores de f, em relação aos valores exatos. O modelo estimado está apresentado a seguir:

f = 0 ,059737 − 0 , 000124 × V + 3,9783 × 10 −9 × Re − 1, 03909 × 10 −15 × Re 2 − 0 ,008767 × LOG 10 ( Re )  50 ,646040  +0 , 000287 × LOG 10 ( V × Re) +     Re

Ciênc. e Agrotec., Lavras, v.23, n.2, p.365-371, abr./jun., 1999

369

0,0405

0,027 Colebrook-White Modelo Estimado D=20mm; V=1,4m/s; 20oC

Colebrook-White Modelo Estimado D=300mm; V=1,57m/s;20oC

0,0255

0,039

0,024

Fator de atrito (f)

Fator de atrito (f)

0,0375

0,036

0,0345

0,0225

0,021

0,0195

0,033 0,018

0,0315

0,03 0,003

0,0165

0,004

0,005

0,006

0,0075

0,009

0,015 0,0005

0,01

0,024

0,00075

0,0009

0,001

0,002

0,00267

Rugosidade Relativa (E/D)

Rugosidade relativa (E/D)

Colebrook-White Modelo Estimado D=600mm; V=1,45m/s; 20oC

0,0225

Fator de atrito (f)

0,021

0,0195

0,018

0,0165

0,015

0,0135

0,012 0,00014

0,00025

0,0005

0,00075

0,001

Rugosidade Relativa (E/D)

FIGURA 1 - Comparação da regressão obtida com a equação de Colebrook-White para os diâmetros de 20, 300 e 600 mm, a 20o C e praticamente à mesma velocidade.

Ciênc. e Agrotec., Lavras, v.23, n.2, p.365-371, abr./jun., 1999

370 Na Fig. 2 estão apresentados os valores estimados pelo modelo proposto e os valores originais de Von Kármán. Verificou-se, novamente, que o modelo de regressão obtido foi altamente satisfatório para predizer o valor de f. Com relação ao uso dos modelos, deve-se verificar a importância do parâmetro E/D, considerando fluxo turbulento de transição. Com o valor de f calculado,

determina-se a sua característica através da equação 3. Se ficar caracterizado fluxo turbulento liso, utiliza-se o modelo correspondente; ficando caracterizado fluxo turbulento de transição, o valor encontrado estará correto; caracterizando fluxo turbulento rugoso, trabalhase com a equação de Nikuradse (equação 6), que não apresenta dificuldades para sua solução.

TABELA 2 - Estimativas dos parâmetros e sua significância para as variáveis do modelo para situação de fluxo turbulento liso. Variável

Estimativas

Prob > |t| para Ho:bi=0

Intercepto

0,059737

0,0001

V

-0,000124

0,0002

Re

3,9783x10-9

0,0001

Re2

-1,03909x10-15

0,0001

-LOG10(Re)

0,008767

0,0001

-LOG10(VxRe)

-0,000287

0,0003

1/Re

50,646040

0,0001

0,0265

Equação de Von-Kármán Modelo estimado V=1,5m/s; 20oC

0,025 0,0235 0,022

Fator f

0,0205 0,019 0,0175 0,016 0,0145 0,013 0,0115 0,01 23347

122525

267327

429703

583168

898515

Número de Reynolds

FIGURA 2 - Comparação do modelo de regressão ajustada com a equação de Von Kármán.

Ciênc. e Agrotec., Lavras, v.23, n.2, p.365-371, abr./jun., 1999

371 Os modelos desenvolvidos não apresentam os inconvenientes de serem específicos para um tipo de tubo (material) podem ser usados para número de Reynolds bem acima de 100.000, não mostrando a magnitude dos desvios apresentados por algumas equações citadas; e o fator f é explícito, não havendo grandes complicadores, principalmente em termos de desenvolvimento de programas de computador para dimensionamentos hidráulicos. Deve-se ressaltar que poderá haver ligeiras diferenças em relação aos dados originais para alguns casos, mas em termos de cálculos hidráulicos em projetos de irrigação, tais diferenças não são significativas. CONCLUSÕES Com base nas soluções obtidas para os valores de f, conclui-se que os modelos desenvolvidos são apropriados para estimar o fator de atrito, desde que sejam usados dentro dos limites estabelecidos. Deve-se ressaltar, também, que esses modelos além de auxiliar o desenvolvimento de programas de computador para cálculos hidráulicos, podem ser usados utilizando-se uma simples calculadora científica, agilizando cálculos do fator f e conseqüentemente de perda de carga, sem perder precisão significativa. REFERÊNCIAS BIBLIOGRÁFICAS AZEVEDO NETTO, J. M.; ALVAREZ, G. A. Manual

de hidráulica. 7. ed. São Paulo: Blücher, 1991. 724p. BERNARDO, S. Manual de irrigação. 5. ed. Viçosa: UFV, 1989. 596p. COSSOLOSSO, M. A.; SATTO, J. Equações explícitas de Colebrook-White. In: CONGRESSO BRASILEIRO DE ENGENHARIA AGRÍCOLA, 25, Anais...Bauru, UNESP, 1996. (CD-ROM). GOMES, H. P. Engenharia de irrigação. 2. ed. Campina Grande: UFPB, 1997. 390p. KAMAND, F. Z. Hydraulic friction factors for pipe flow. Journal of Irrigation and Drainage Engineering, New York, v. 114. n.2. p.311-323, may, 1988. SYSTEM ANALYSES STATISTICAL. Language guide for personal computers. 6. ed. Cary: SAS Institute Inc, 1985. 429p. VON BERNUTH, R. D. Simple and accurate friction loss equation for plastic pipe. Journal of Irrigation and Drainage Engineering, New York, v.116, n.2, p.294-298, Mar./Apr. 1990. YANAGI JUNIOR, T. Dimensionamento e predição da distribuição de água em sistemas de irrigação por aspersores autopropelidos. Lavras: UFLA, 1995. 79p. (Dissertação-Mestrado em Irrigação e Drenagem).

Ciênc. e Agrotec., Lavras, v.23, n.2, p.365-371, abr./jun., 1999