1.1 ESCOAMENTO EM CANAIS
O escoamento em condutos livres (ou canais) está em contato com a superfície livre e, por isso, está sujeito à pressão atmosférica (Fig. 1-1). Por outro lado, o escoamento em condutos forçados (ou tubulações) está sujeito apenas à pressão hidráulica que atua no conduto. Devido à sua superfície livre, a análise do fluxo de um canal é geralmente mais complexa que a de um conduto forçado, pois a superfície livre é susceptível a variações espaciais e temporais. Quando o escoamento de um canal não varia no espaço, ele é chamado de uniforme, do contrário ele é variado. Se o escoamento não varia no tempo, ele é chamado de permanente, mas se varia, ele é não permanente.
O escoamento dos condutos forçados é delimitado pela seção transversal das paredes do tubo. Já no canal aberto, o escoamento não tem seção transversal fixa, variando de acordo com a vazão a ser escoada. Nos condutos forçados, a rugosidade varia de acordo com o material utilizado para a sua confecção, entre outros fatores, desde um latão liso a tubos corroídos, enquanto que nos canais, o rugosidade varia de vidro acrílico ou lucite® (um tipo de plástico muito suave), ambos usados em experimentos laboratoriais, até a dos materiais dos leitos e margens de canais naturais, além das planícies de inundação. A pressão hidráulica no centro dos condutos forçados define a linha do gradiente hidráulico ou linha piezométrica (Fig. 1-2). A linha de energia corresponde à V 2/(2g) na seção transversal. A perda de carga, ou perda de carga por atrito, da seção 1 à 2, é representadax por hf.
Nos canais, a profundidade do fluxo, medida a partir do leito, define a elevação da superfície da água, equivalente à linha do gradiente hidráulico de um conduto forçado (veja a figura 1-3). A linha de energia total se refere à parcela de energia cinética V 2/(2g) na seção transversal. A perda de carga, ou perda de carga por atrito, da seção 1 à 2, é representada por hf.
Observa-se que no conduto forçado a água sobe até o nível que define a linha do gradiente hidráulico, associada à pressão hidráulica. Por outro lado, no conduto livre, a superfície da água é a linha do gradiente hidráulico que está sujeita à pressão atmosférica. 1.2 TIPOS DE ESCOAMENTO
Existem dois tipos gerais de seções transversais de canais:
Os canais artificiais (ou feitos pelo homem) são geralmente prismáticos, apresentando uma forma e tamanho constantes, pelo menos para um determinado comprimento de canal. Por outro lado, os canais naturais são tipicamente não prismáticos, isto é, a forma e o tamanho da seção transversal variam ao longo do canal. Várias propriedades geométricas e hidráulicas ajudam a descrever um canal aberto (Fig. 1-4), entre elas:
Nos canais prismáticos, a profundidade do fluxo d é conhecida como y, particularmente quando não pode ser confundida com a altura do nível da água obtida pela leitura da régua limnimétrica. Além disso, a inclinação do lado do canal é frequentemente referida como zH:1V, particularmente quando não pode ser confundida com a elevação do leito. Classificação do escoamento em canais O fluxo de um canal pode ser classificado da seguinte forma:
O fluxo é permanente quando as variáveis hidráulicas (vazão, área hidráulica, velocidade média, profundidade do fluxo, etentre outras) não variam ao longo do tempo. Do contrário, o fluxo é instável. O escoamento permanente é relativamente simples de calcular, quando comparado com o instável. O fluxo é uniforme quando o canal é prismático e as variáveis hidráulicas (Q, A, V, d, etc) são constantes no tempo e no espaço. O fluxo está em equilíbrio quando o canal não é prismático e as variáveis hidráulicas são aproximadamente constantes no tempo e no espaço. O cálculo do fluxo uniforme é relativamente direto quando comparado ao de outros estadotipos de escoamento. O fluxo varia gradualmente quando a vazão Q é constante, mas as outras variáveis hidráulicas (A, V, d, etc) variam gradualmente no espaço. Sob o fluxo gradualmente variado, a distribuição de pressão na direção vertical, normal ao fluxo, é muito próxima da hidrostática, isto é, proporcional à profundidade do fluxo. O fluxo varia bruscamente quando a vazão é constante, mas as outras variáveis hidráulicas (A, V, d, etc.) variam bruscamente no espaço, de modo que uma distribuição de pressão hidrostática não possa ser assumida na direção vertical normal ao fluxo. Embora o cálculo do fluxo variado gradualmente seja mais factível, o cálculo do fluxo variado rapidamente é geralmente mais complexo, baseado em fórmulas empíricas, por falta de uma solução teórica. O fluxo é espacialmente variado quando a descarga Q varia apenas no espaço (ou seja, ao longo do canal), normalmente devido à alguma entrada ou saída lateral. Ocorrência de vários tipos de escoamento O escoamento uniforme ocorre em um canal prismático (Fig. 1-5), enquanto que o escoamento em equilíbrio ocorre em um canal não prismático. O fluxo uniforme instável não existe na natureza, porque o fluxo não pode ser uniforme e instável ao mesmo tempo. A palavra "instável" implica em não equilíbrio; portanto, não existe fluxo de equilíbrio instável.
O fluxo gradualmente variado e constante é representado pelos perfis da superfície da água, também conhecidos como cálculos de água de retorno (ou rebaixamento)(consulte o Capítulo 7). O fluxo instável e gradualmente variado é utilizado nos cálculos dos fluxos de inundação, ou roteamento de inundação (consulte o Capítulo 10). O fluxo constante e bruscamente variado é representado pelo fluxo nos vertedouros ou pelo salto hidráulico. O fluxo instável e bruscamente variado é representado pelo salto hidráulico em movimento, surtos, ondas de rolagem, choques cinemáticos e remansos de maré. A Figura 1-6 mostra ondas rolantes em um canal de irrigação íngreme. O fluxo espacialmente variado ocorre em um canal artificial quando a vazão varia ao longo do canal, devido a extrações laterais da água ou transbordamento do canal.
1.3 ESTADO DE ESCOAMENTO
O estado do fluxo de um canal pode ser descrito em função das velocidades e viscosidades características. Velocidade é a razão do comprimento (distância) ao longo do tempo, com unidades Razões de velocidade Existem três velocidades características no fluxo de um canal:
A velocidade média do fluxo uniforme usando a equação de Manning (unidades SI) é:
Em que: n = coeficiente de atrito de Manning; A velocidade média do fluxo uniforme constante usando a equação de Chezy é:
Em que: C = coeficiente de atrito de Chezy. Em geral, quatro forças estão ativas em um volume de controle no fluxo de canal. Essas forças são causadas por atrito, gravidade, gradiente de pressão (profundidade do fluxo) e inércia. Ondas cinemáticas são aquelas em que o equilíbrio do momento é expresso apenas em termos das forças de atrito e forças gravitacionais (Lighthill and Whitham, 1955). A celeridade das ondas cinemáticas, ou celeridade de Seddon, é (Seddon, 1900; Chow, 1959; Ponce, 1989):
Em que: β = expoente de classificação da área da curva-chave A, que é definida da seguinte forma:
Em que: α = razão da profundidade hidráulica D com a profundidade do fluxo y. As ondas dinâmicas são aquelas em que o equilíbrio do momento é expresso apenas em termos do gradiente de pressão e das forças inerciais. A celeridade das ondas dinâmicas é:
Em que: g = aceleração da gravidade; D = profundidade hidráulica, D = A /T. Da Eq. 1-3, a celeridade relativa das ondas cinemáticas é:
Da Eq. 1-5, o valor absoluto da celeridade relativa das ondas dinâmicas é:
Para canais retangulares (onde, D = d ), ou para os canais de grande largura (onde, D ≅ d ), a celeridade relativa das ondas dinâmicas é:
A equação 1-8 é conhecida como equação de celeridade de Lagrange (relativa), em homenagem à Lagrange (1788), que a derivou pela primeira vez. Número de Froude O número de Froude é definido da seguinte forma (Chow, 1959):
O número Froude caracteriza a condição de:
Sob o escoamento subcrítico, as ondas de superfície (perturbações) podem fluir para a montante, porque a sua celeridade à montante -w é maior que a velocidade média do fluxo -u. Sob o escoamento crítico, as ondas superficiais (perturbações) permanecem estacionárias, porque a sua celeridade (absoluta) w é igual à velocidade média do fluxo u. Sob o escoamento supercrítico, as ondas de superfície (perturbações) podem fluir apenas para a jusante, porque a sua celeridade à montante -w é menor que a velocidade média do fluxo u. Número de Vedernikov O número de Vedernikov é definido da seguinte forma (Vedernikov, 1945; 1946; Powell, 1948; Craya, 1952):
O número de Vedernikov caracteriza os seguintes estados de fluxo:
Sob fluxo estável, a celeridade relativa da onda cinemática v é menor que a celeridade relativa da onda dinâmica w e, portanto, as ondas de superfície (perturbações) são capazes de atenuar (dissipar). Sob fluxo neutro estável, a celeridade relativa da onda cinemática v é igual à celeridade relativa da onda dinâmica w e, portanto, as ondas superficiais (perturbações) não atenuam nem amplificam. A amplificação equivale à dissipação negativa. Sob um fluxo instável, a celeridade relativa da onda cinemática v é maior que a celeridade relativa da onda dinâmica w. Portanto, as ondas de superfície (perturbações) estão sujeitas à amplificação. Na prática, a condição V ~ 1 leva ao desenvolvimento de ondas rolantes, que formam um conjunto de ondas que viajam para rio abaixo, tipicamente em canais artificiais de declive acentuado (Cornish, 1907) (Fig. 1-7).
As três velocidades u, v, e w, levam a apenas duas relações de velocidade independentes, que são os números de Froude (Eq. 1-8) e de Vedernikov (Eq. 1-9). A terceira relação é a celeridade da onda cinemática relativa adimensional, que corresponde ao expoente da classificação da área da curva-chave subtraida de 1:
Dessa forma, pode-se verificar que o expoente β na Eq. 1-4 está em função dos números de Froude e de Vedernikov. O valor de β varia de acordo com o tipo de regime de atrito (laminar, transitório ou turbulento; e Manning ou Chezy turbulento) e com a forma da seção transversal. Para o escoamento laminar, β = 3. Para o escoamento turbulento, sob atrito de fricção: 1 ≤ β ≤ 5/3, dependendo da forma da seção transversal. Para o escoamento turbulento, sob atrito Chezy: ≤ β ≤ 3/2, dependendo da forma da seção transversal. Tipos de forma da seção transversal Existem três formas transversais assintóticas em canais:
Fluxo neutro estável Para a estabilidade neutra: V = 1. Portanto, da Eq. 1-11, o número de Froude correspondente ao fluxo neutro estável é:
A Tabela 1-1 mostra os valores de Fns para os valores selecionados de β. Observa-se que, como β varia de 3 (fluxo laminar) para 1 (canal inerentemente estável), o valor de Fns varia de 1/2 a ∞. Em outras palavras, como β tende a 1, Fns tende a ∞. Na prática, como o atrito possui efetivamente um limite inferior, o número de Froude é restrito a um limite superior, que raramente excede F = 25. Portanto, na maioria dos casos, um valor de β = 1,04 já seria estável para fins práticos.
Conforme o mostrado na Tabela 1-1, os valores de β no fluxo de canal e no transbordamento são limitados no intervalo Relações de viscosidade Existem três viscosidades características no fluxo de canais:
A viscosidade cinemática ν do fluido varia em função da temperatura (Apêndice A). O conceito de difusividade hidráulica νh é devido a Hayami (1951). Hayami combinou as equações de fluxo de canais (Capítulo 10) para desenvolver uma única equação de convecção-difusão, ou seja, uma equação que descreve a convecção (processo de primeira ordem) e a difusão (processo de segunda ordem) de uma onda de inundação. A difusividade hidráulica é definida da seguinte forma:
Em que: qo = vazão da largura unitária; Observa-se que a difusão das ondas de inundação é diretamente proporcional à vazão de largura unitária e inversamente proporcional à inclinação do atrito. A equação 1-13 pode ser expressa em termos de velocidade uo e profundidade do fluxo do da seguinte maneira:
Um valor relacionado de difusividade, independente da inclinação, é:
Em geral, para uma forma de seção transversal arbitrária:
Em que: Ro = raio hidráulico Na teoria das ondas cinemáticas, o comprimento de alcance característico é definido da seguinte forma (Lighthill and Whitham, 1955):
Em que: Lo = comprimento do canal no qual o fluxo de equilíbrio decai no mesmo valor que a sua profundidade. Assim, em termos do comprimento de alcance característico, a difusividade hidráulica é:
De maneira semelhante à difusividade hidráulica, a difusividade das ondas é convenientemente definida da seguinte forma:
Em que: L = comprimento de onda da perturbação. Número de Reynolds O número de Reynolds R pode ser expresso da seguinte forma (Chow, 1959):
O número de Reynolds R descreve o regime de fluxo como:
Sob condições de fluxo estável no fluxo de canal aberto, o fluxo laminar ocorre para R ≤ 500 e fluxo turbulento para R > 2000. O fluxo de transição ocorre na faixa intermediária: 500 < R ≤ 2000. Sob fluxo instável, o fluxo laminar-turbulento misto descrito na Tabela 1-1 é semelhante ao fluxo de transição, apresentando uma faixa comparável de números de Reynolds. Na prática, para a maioria dos casos de escoamento em canais, o regime é do tipo turbulento. Por outro lado, a maior parte dos casos de escoamento em transbordamentos (ou seja; em superfícies livres planas) são do tipo laminar ou misto turbulento-laminar. Número de onda adimensional O número de onda adimensional σ é definido da seguinte forma (Ponce and Simons, 1977):
O número de onda σ descreve a escala de comprimento adimensional da onda, como mostrado na Fig. 1-11, em termos de:
A Figura 11 é aplicável ao caso do Coeficiente de atrito de Chezy em canais hidraulicamente amplos. Sob o escoamento cinemático, descrito no lado esquerdo da Fig. 1-11:
Sob o fluxo dinâmico, descrito no lado direito da Fig. 1-11:
Sob o fluxo misto cinemático-dinâmico, representado na seção central da Fig. 1-11:
Difusividade hidráulica dinâmica
A difusividade hidráulica dinâmica, que, diferentemente da difusividade hidráulica de Hayami da
Para números de Vedernikov baixos, V ⇒ 0, a difusividade hidráulica dinâmica se reduz à expressão de difusividade hidráulica cinemática, isto é, Eq. 1-18. Por outro lado, para números altos de Vedernikov, V ⇒ 1, a difusividade hidráulica dinâmica desaparece. Sob essa condição de fluxo, a ausência total de atenuação das ondas é propícia ao desenvolvimento de ondas rolantes (Figs. 1-6 e 1-7).
1.4 REGIMES DE ESCOAMENTO
Os regimes de fluxo de um canal podem ser do tipo:
Os regimes de escoamento são caracterizados pelo número R de Reynolds, Eq. 1-20. No fluxo de um canal, o regime laminar prevalece para R ≤ 500, o regime transitório para 500 < R ≤ 2000 e o regime turbulento para R > 2000. Os regimes de escoamento variam com a rugosidade da superfície do canal. A Figura 1-12 mostra a relação entre o número de Reynolds R e o fator de atrito de Darcy-Weisbach f para o escoamento em canais lisos. A Figura 1-13 mostra a relação entre R e f para o escoamento em canais rugosos. A fórmula de atrito de Darcy-Weisbach, desenvolvida em paralelo com a do escoamento em tubulações, é:
Em que: hf = perda de carga por atrito;
Observando-se as Fig. 1-12 e 1-13 pode-se concluir que:
Fator de atrito de Darcy-Weisbach para o escoamento de canais A fórmula de Darcy-Weisbach, Eq. 1-23, é estritamente aplicável ao escoamento de condutos forçados. Nessas tubulações, o comprimento de atrito característico é o diâmetro do tubo. Por outro lado, em canais, o comprimento de atrito característico é o raio hidráulico R, isto é, a razão entre a área do escoamento e o perímetro molhado:
Como a área de escoamento de um tubo circular (totalmente cheio) é A = π do2/4, e o perímetro molhado é
Em que: V = velocidade média do escoamento no canal. Em canais com escoamento uniforme, a inclinação da energia, é a mesma do atrito e do leito:
Para uma seção transversal de forma arbitrária, o número de Froude é:
Em que: D = profundidade hidráulica: D = A /T. A equação 1-30 pode ser expressa em relação ao número de Froude, da seguinte maneira:
A equação 1-32 indica a proporcionalidade entre a inclinação da energia e o número de Froude, com o fator de proporcionalidade, sendo uma função do fator de atrito de Darcy-Weisbach e do fator de forma D /R. Para um canal hidraulicamente amplo, no qual D ≅ R, a Eq. 1-32 se reduz para:
Assim, para um canal hidraulicamente amplo, o fator de proporcionalidade entre a inclinação da energia e o número de Froude é apenas uma função do fator de atrito de Darcy-Weisbach f. Essencialmente, para o escoamento de canais, deve-se utilizar o fator modificado de Darcy-Weisbach f, igual a 1/8 do fator de atrito convencional f de Darcy-Weisbach. A equação de Darcy-Weisbach modificada para escoamento em canais é:
A Tabela 1-2 mostra os valores aproximados do fator de atrito f de Darcy-Weisbach e o correspondente fator de atrito modificado f para valores selecionados de R na faixa turbulenta.
QUESTÕES
PROBLEMAS
REFERÊNCIAS
Chow, V. T. 1959. Open-channel Hydraulics. Mc-Graw Hill, New York. Cornish, V. 1907. Progressive waves in rivers. The Geographical Journal. Vol. 29, No. 1, January, 23-31. Craya, A. 1952. The criterion for the possibility of roll wave formation. Gravity Waves, Circular 521, 141-151, National Institute of Standards and Technology, Gaithersburg, Md. Dooge, J. C. I., W. B. Strupczewski, and J. J. Napiorkoswki. 1982. Hydrodynamic derivation of storage parameters in the Muskingum model. Journal of Hydrology, 54, 371-387. Fread, D. 1993. "Flow Routing," Chapter 10 in Handbook of Hydrology, D. R. Maidment, editor, McGraw-Hill, New York. Hayami, I. 1951. On the propagation of flood waves. Bulletin, Disaster Prevention Research Institute, No. 1, December. Lagrange, J. L. de. 1788. Mécanique analytique, Paris, part 2, section II, article 2, p 192. Lighthill, M. J., and G. B. Whitham. 1955. On kinematic waves: I. Flood movement in long rivers. Proceedings, Royal Society of London, Series A, 229, 281-316. Ponce, V, M., and D. B. Simons. 1977. Shallow wave propagation in open-channel flow. Journal of the Hydraulics Division, ASCE, Vol. 103, No. HY12, December, 1461-1476. Ponce, V. M. 1989. Engineering Hydrology: Principles and Practices. Prentice-Hall, Englewood Cliffs, New Jersey. Ponce, V. M. 1991a. The kinematic wave controversy. Journal of Hydraulic Engineering, ASCE, Vol. 117, No. 4, April, 511-525. Ponce, V. M. 1991b. New perspective on the Vedernikov number. Water Resources Research, Vol. 27, No. 7, 1777-1779, July. Ponce, V. M., and P. J. Porras. 1995. Effect of cross-sectional shape on free-surface instability. Journal of Hydraulic Engineering, ASCE, Vol. 121, No. 4, April, 376-380. Powell, R. W. 1948. Vedernikov's criterion for ultra-rapid flow. Transactions, American Geophysical Union, Vol. 29, No. 6, 882-886. Seddon, J. A. 1900. River hydraulics. Transactions, ASCE, Vol. XLIII, 179-243, June. Vedernikov, V. V. 1945. Conditions at the front of a translation wave disturbing a steady motion of a real fluid, Dokl. Akad. Nauk SSSR, 48(4), 239-242. Vedernikov, V. V. 1946. Characteristic features of a liquid flow in an open channel, Dokl. Akad. Nauk SSSR, 52(3), 207-210.
|
Documents in Portable Document Format (PDF) require Adobe Acrobat Reader 5.0 or higher to view; download Adobe Acrobat Reader. |