Simulação numérica do escoamento superficial por meio do método do reticulado de Boltzmann utilizando junção de sub-bacias

Latest revision as of 12:20, 1 October 2021

Resumo

Investigou-se o escoamento superficial em bacia hidrográfica por meio da simulação numérica unidimensional utilizando o método do reticulado de Boltzmann (LBM). Desenvolveu-se um modelo computacional onde a bacia hidrográfica é representada pela junção de nove sub-bacias. Para isto, foram estabelecidas duas funções de distribuição de equilíbrio obtidas por meio da expansão de Chapmann-Enskog e utilizando o reticulado D1Q5, uma adequada para o escoamento na superfície da bacia e outra para o canal principal, obtendo a profundidade da água na superfície da bacia e a área da seção transversal nos canais. Além disso, estabeleceu-se a condição de contorno na passagem de fluxo de uma sub-bacia para outra, levando em consideração a conservação da massa e, para se obter uma simulação mais próxima da realidade, considerou-se uma vazão inicial que ocorre no rio (escoamento de base), em cada trecho de canal. Os resultados numéricos obtidos pelo LBM foram comparados com dados medidos em campo.

Palavras-chaves: Método do reticulado de Boltzmann, escoamento em bacias hidrográficas, escoamento superficial, modelo onda cinemática, junção de sub-bacias

Abstract

The watershed surface runoff was investigated through the one-dimensional numerical simulation using the Lattice Boltzmann Method (LBM). A computational model was developed where the watershed is represented by the junction of nine sub-basins. For this, two equilibrium distribution functions were established through the Chapman-Enskog Expansion on a D1Q5 lattice, one suitable for flow on the basin surface and another for the main channel, obtaining the water depth on the basin surface and the channels cross-sectional area. In addition, the boundary condition was established in the flow passage from one sub-basin to another, taking into account the mass conservation and, in order to obtain a simulation closer to reality, it was considered an initial river flow (baseflow) of each channel stretch. The numerical results obtained by the LBM were compared with data measured in field.

Keywords: Lattice Boltzmann method, runoff in watersheds, overland flow, kinematic wave model, sub-basin junction

1. Introdução

O comportamento do escoamento de um fluido pode ser, satisfatoriamente, avaliado por meio de métodos numéricos. Simulações numéricas, análises teóricas, técnicas experimentais e dados coletados, são ferramentais que se complementam na resolução de problemas em hidrodinâmica [1]. Em particular, o LBM se destaca por capturar as características físicas de problemas, até mesmo detalhes microscópicos, e as apresenta de forma macroscópica [2]. A abordagem do LBM é diferenciada quando comparada aos métodos numéricos tradicionais, pois fornece uma maneira indireta para a solução das equações governantes do escoamento. O LBM não utiliza discretização das equações macroscópicas, o método baseia-se em modelos microscópicos e equações que governam a cinética em um nível mesoscópico. No LBM a dinâmica macroscópica de um fluido é resultado do comportamento coletivo de partículas microscópicas e não se altera com os detalhes subjacentes referentes as interações moleculares do fluido [3]. O método tem sido aplicado nas mais diversas áreas de pesquisa e tornou-se de grande potencial no estudo da dinâmica de fluidos computacional. Escoamentos envolvendo águas rasas [2,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20], o modelo onda cinemática [6,21,22,23,25,24] e várias outras aplicações como em meios porosos [26,30,27,28,29] foram abrangidos pelo LBM e as simulações apresentaram resultados satisfatórios.

As equações de Saint-Venant governam o escoamento em superfície livre. Porém, para sua resolução, elaboradas técnicas numéricas e uma grande quantidade de informações hidráulicas são necessárias. Na busca em reduzir a quantidade de dados necessários e a dificuldade numérica de resolução das equações diferenciais, são empregadas simplificações das equações de Saint-Venant como o modelo onda cinemática, também denominado modelo hidráulico de propagação de ondas de cheias [31].

O modelo onda cinemática é aplicado para descrever escoamentos superficiais em bacias hidrográficas e cursos de água naturais. Foco de interesse neste estudo, as bacias hidrográficas são um sistema natural que exige monitoramento no sentido de prever seu comportamento em ações como precipitações extremas, cheias, estiagens, entre outros. Com a análise destes eventos, é possível gerenciar os recursos hídricos e evitar a degradação do solo, oferecendo assim planejamento adequado para o espaço rural e urbano [1].

Utilizando o LBM com aperador de colisão BGK, a proposta deste estudo é obter a simulação numérica do escoamento superficial na bacia do rio Chopim, representado pelas equações do modelo onda cinemática. A bacia do rio Chopim abrange uma área de ${\textstyle 7.500}$ km ${\textstyle ^{2}}$ e está localizada dentro da bacia do rio Iguaçu, a qual tem grande importância na geração de energia por meio de usinas hidrelétricas e pela grande demanda de recursos hídricos. Considera-se, nesse estudo, um trecho do rio Chopim com comprimento de aproximadamente ${\textstyle 281}$ km com área de drenagem em torno de ${\textstyle 3.205}$ km ${\textstyle ^{2}}$ na superfície da bacia.

Lembramos que um artigo, publicado em 2018, dos mesmos autores deste artigo, fornece um estudo menos amplo que o apresentado aqui. No artigo [24] obteve-se a simulação numérica unidimensional do escoamento superficial em uma bacia hidrográfica natural, por meio do LBM. De acordo com as caractrísticas da bacia hidrográfica estudada, foi possivel obter bons resultados na simulação considerando no modelo computacional a superfície da bacia segmentada em apenas dois planos laterais e um canal principal. No novo artigo que estamos apresentando, analizamos outra bacia hidrográfica, que devido as suas características, foi necessário segmentá-la em 9 sub-bacias. Com esta nova configuração, foi necessário estabelecer condições de contorno na passagem de fluxo de uma sub-bacia para outra, levando em consideração a conservação da massa. No artigo [24] determinou-se duas funções de distribuição de equilíbrio obtidas por meio da expansão de Chapmann-Enskog em escalas de tempo e utilizando o reticulado D1Q3. Neste novo artigo, as funções de equilíbrio foram determinadas para o reticulado D1Q5. Além destas novas contribuições, para se obter uma simulação mais próxima da realidade, considerou-se uma vazão inicial que ocorre no rio (escoamento de base), em cada trecho de canal, antes de ocorrer a precipitação.

As quatro principais contribuições deste artigo são: 1) Aplicação do LBM para obter a simulação numérica do escoamento superficial da bacia do rio Chopim representada pela junção de nove sub-bacias formadas por planos laterais e um segmento de canal. Em cada sub-bacia, a água que escoa dos planos entra no canal sob a forma de contribuição lateral. Na saída de cada sub-bacia, obtém-se um hidrograma, que foi utilizado para ajustar as condições de contorno na junção de dois segmentos de canal. A simulação produziu o hidrograma total o qual fornece a variação da vazão em relação ao tempo correspondente a resposta da bacia hidrográfica após a precipitação efetiva do evento simulado; 2) Determinação de duas funções de distribuição de equilíbrio obtidas por meio da expansão de Chapmann-Enskog em escalas de tempo e utilizando o reticulado D1Q5, uma adequada para o escoamento na superfície da bacia e outra para o canal principal, obtendo a profundidade da água na superfície da bacia e área da seção transversal nos canais; 3) Determinação de condição de contorno na passagem de fluxo de uma sub-bacia para outra, levando em consideração a conservação da massa; 4) Para se obter uma simulação mais próxima da realidade, considerou-se uma vazão inicial que ocorre no rio (escoamento de base), em cada trecho de canal, antes de ocorrer a precipitação. Essa vazão entra no modelo computacional como uma condição inicial. Também é considerada como uma condição de contorno ao ser adicionada na contribuição lateral dos planos, em cada passo de tempo.

Este artigo está estruturado da seguinte forma. Na seção 2 tem-se as equações governantes do escoamento. Em seguida, na seção 3, apresenta-se o método numérico que será utilizado e os quatro elementos principais que o caracterizam: Equação do LBM e o operador de colisão BGK, o reticulado D1Q5, a expensão multi-escala e as funções distribuição de equilíbrio. Tem-se, também, nesta seção informações sobre condições de contorno e estabilidade. O estudo de caso é apresentado na seção 4 considerando as características da bacia do rio Chopim, resultados da simulação e discussões. Em seguida, na seção 6, apresentam-se as conclusões deste trabalho.

2. Equações de escoamento

As equações de Saint-Venant unidimensionais formam um sistema de equações não-lineares composto pelas equações da continuidade (1) e da quantidade de movimento (2). Elas governam o escoamento onde há superfície livre, com a suposição de que a componente vertical e transversal da velocidade do escoamento pode ser desprezada em relação aos componentes longitudinais [32],

(1)

(2)

onde ${\textstyle Q}$ é a vazão, ${\textstyle A}$ é a área da seção transversal, ${\textstyle h}$ é a profundidade do escoamento, ${\textstyle g}$ é a aceleração da gravidade, ${\textstyle t}$ é o tempo, ${\textstyle x}$ é a coordenada espacial, ${\textstyle S_{0}}$ é a declividade do fundo e ${\textstyle S_{f}}$ representa a declividade da linha de energia.

No escoamento em bacias hidrográficas utiliza-se o modelo onda cinemática como simplificação das equações de Saint-Venant. Este modelo, considera a equação da continuidade e a equação da quantidade de movimento, desprezando os termos de pressão e inércia. Dessa forma, a equação da quantidade de movimento resulta [33],

(3)

O modelo onda cinemática baseia-se, principalmente, na equação da continuidade e faz-se uma aproximação da equação da quantidade de movimento por uma fórmula de fluxo uniforme. Definindo a declividade da linha de energia ${\textstyle S_{f}}$ , com uma fórmula de fluxo uniforme, a equação (3) pode ser representada para o escoamento em canal ou em escoamento superficial por uma relação de potência da forma [34],

(4)

em que ${\textstyle \alpha }$ e ${\textstyle \lambda }$ são coeficientes determinados pelas características do escoamento. Este modelo não-linear considera a variabilidade dos parâmetros de acordo com a vazão. Assim, o escoamento em trechos de canais da bacia é dado pelas equações,

(5)

onde ${\textstyle q_{l}}$ é a contribuição lateral, ${\textstyle \beta _{c}}$ e ${\textstyle m_{c}}$ são parâmetros a serem determinados.

As condições iniciais e de contorno nos canais são determinadas por hidrogramas de vazão resultantes do escoamento na superfície e em trechos de canais a montante. A superfície da bacia é aproximada por planos, dessa forma, considerando que a equação da continuidade (1) descreva o escoamento em um canal prismático e retangular, as equações de escoamento para os planos que descrevem a superfície da bacia são as seguintes,

(6)

onde ${\textstyle q}$ é a vazão por unidade de largura, ${\textstyle i_{e}}$ é a precipitação efetiva e ${\textstyle \beta _{s}}$ e ${\textstyle m_{s}}$ são parâmetros a serem determinados.

As condições iniciais e de contorno são representadas por,

(7)

(8)

onde ${\textstyle L}$ é o comprimento de declive.

A aplicação do modelo onda cinemática em escoamento superficial difere da sua aplicação em rios, apenas pelo meio onde ocorre o escoamento. Quanto a rugosidade, utiliza-se a resistência hidráulica de Manning, assim ${\textstyle \beta _{c}=\beta _{s}={\sqrt {S_{0}}}/n_{r}}$ e ${\textstyle m_{c}=m_{s}=5/3}$ , onde ${\textstyle n_{r}}$ é o coeficiente de rugosidade de Manning.

3. Método do reticulado de Boltzmann

De acordo com Chen and Doolen [3], o LBM é um esquema numérico baseado em equações cinéticas formuladas em uma escala mesoscópica que simula a dinâmica do fluido em uma escala macroscópica. A dinâmica deste método é dada pela equação do reticulado de Boltzmann, a qual mostra como ocorrem as distribuições das micropartículas.

O LBM considera uma dinâmica molecular de partículas fictícias em que o espaço, o tempo e as velocidades são discretas. Essas partículas viajam de um ponto ao outro na malha em tempos discretos, se encontrando nestes pontos ao final de cada passo de tempo e trocam quantidade de movimento e energia. No que se refere aos módulos de velocidades, eles assumem valores contínuos. Por outro lado, o LBM pode ser visto como uma forma simplificada da equação cinética de Boltzmann, em que somente são mantidos os detalhes moleculares essenciais para se recuperar o comportamento macroscópico correto [35].

Quatro elementos principais caracterizam o LBM: a equação governante, o operador de colisão, o reticulado e a função distribuição de equilíbrio [2].

3.1 Equação do LBM e o operador BGK

Inicialmente apresenta-se em equação (9), a equação governante do LBM com aperador de colisão BGK, denominada equação do reticulado de Boltzmann. A forma do operador de colisão foi determinada por Bhatnagar et al. [36], com o objetivo de simplificar a equação cinética de Boltzmann. A equação do reticulado de Boltzmann, possui as etapas de propagação e colisão. Na etapa de propagação, as partículas se movem de um nó da malha para um dos seus nós vizinhos, com direção dada pela velocidade. Na etapa de colisão, as partículas que chegam no mesmo nó interagem entre si e mudam suas direções conforme as diretrizes do operador de colisão [37]

(9)

(10)

onde ${\textstyle \Omega _{\alpha }}$ é o operador de colisão, ${\textstyle f_{\alpha }}$ é a função distribuição de partículas e representa a probabilidade de uma partícula seguir determinada direção, ${\textstyle f_{\alpha }^{eq}}$ é a função distribuição de equilíbrio, ${\textstyle {\vec {x}}}$ é o posição do nó na malha, ${\textstyle \Delta t}$ é o incremento no tempo, ${\textstyle {\vec {e_{\alpha }}}}$ são as possíveis direções de movimento na malha, ${\textstyle {(\Delta t)^{2}}{g_{\alpha }}}$ é o termo de força e ${\textstyle \tau }$ é o parâmetro de relaxamento.

3.2 Reticulado D1Q5

No LBM o reticulado tem a função de representar os pontos da malha e determinar as direções de movimento das partículas. As direções finitas e determinadas para o movimento das partículas definem num modelo microscópico para a dinâmica molecular. Na escolha do reticulado a ser utilizado é essencial observar a sua simetria. Ela é necessária para que o reticulado possa representar as equações macroscópicas [2].

Encontra-se em Qian et al. [30], uma família de reticulados denominada ${\textstyle DkQm}$ , onde ${\textstyle k}$ indica o espaço ${\textstyle k}$ -dimensional e ${\textstyle m}$ as direções de movimento para as distribuições de partículas. O modelo de reticulado D1Q5 é unidimensional com quatro direções de velocidades não nulas e uma para velocidade nula. As partículas do fluido se movimentam para dois nós vizinhos a esquerda e dois nós vizinhos a direita. A Figura 1 mostra o reticulado D1Q5 com velocidades ${\textstyle {\vec {e}}_{0}=0}$ , ${\textstyle {\vec {e}}_{1}=-e}$ , ${\textstyle {\vec {e}}_{2}=-2e}$ , ${\textstyle {\vec {e}}_{3}=e}$ e ${\textstyle {\vec {e}}_{4}=2e}$ , sendo ${\textstyle e=\Delta x/\Delta t}$ a velocidade no reticulado, em que ${\textstyle \Delta x}$ é o espaçamento da malha.

Figura 1. Reticulado D1Q5

Fonte: Zhang e Yan [38]

3.3 Expansão multi-escala

Uma expansão multi-escala deve ser aplicada para fazer a ligação entre a escala mesoscópica, onde se insere o LBM, e a escala macroscópica das equações governantes. A expansão mais utilizada é a expansão de Chapmann-Enskog [39], a qual considera escalas de tempo e espaço para que a partir da equação do reticulado de Boltzmann seja derivada e equação governante do escoamento.

Para que não hajam funções de distribuição negativas associadas a equação de onda cinemática, o que contraria as leis da física, deve-se atribuir a expansão multi-escala apenas a coordenada de tempo [22]. Dessa forma, aplicando-se cinco escalas a coordenada temporal, as formas diferenciais para as coordenadas de espaço e tempo são como seguem,

{\begin{aligned}&{\dfrac {\partial }{\partial t}}={\dfrac {\partial }{\partial {t_{0}}}}+\varepsilon {\dfrac {\partial }{\partial {t_{1}}}}+{\varepsilon ^{2}}{\dfrac {\partial }{\partial {t_{2}}}}+{\varepsilon ^{3}}{\dfrac {\partial }{\partial {t_{3}}}}+{\varepsilon ^{4}}{\dfrac {\partial }{\partial {t_{4}}}}+O({\varepsilon ^{5}}),\\&{\dfrac {\partial }{\partial x}}={\dfrac {\partial }{\partial x}}+O(\varepsilon ),\end{aligned}}

(11)

em que ${\textstyle t_{0},\,t_{1},\,t_{2},\,t_{3}}$ e ${\textstyle t_{4}}$ são escalas de tempo e ${\textstyle \varepsilon }$ é o número de Knudsen. O número de Knudsen é uma medida adimensional definido como a razão entre o comprimento do caminho livre médio molecular e uma escala de comprimento fisicamente representativa.

Expandindo ${\textstyle f_{\alpha }}$ em torno da função de distribuição de equilíbrio ${\textstyle f_{\alpha }^{eq}}$ com ${\textstyle \varepsilon }$ pequeno, tem-se,

(12)

em que ${\textstyle f_{\alpha }^{(0)}}$ é função distribuição de equilíbrio,

(13)

e, a função distribuição de não equilíbrio é descrita na forma,

(14)

A expansão de ${\textstyle {f_{\alpha }}({\vec {x}}+{{\vec {e}}_{\alpha }}\Delta t,\,t+\Delta t)}$ em série de Taylor de quinta ordem, considerando ${\textstyle \Delta t=\varepsilon }$ , comparada com a equação (9) de evolução do LBM, resulta em,

(15)

em que,

(16)

Com base nas escalas de tempo e na expansão de ${\textstyle f_{\alpha }}$ , equações (11) e (12), respectivamente, e comparando todas as ordens de ${\textstyle \varepsilon }$ , determinam-se as equações discretas do reticulado de Boltzmann de diferentes ordens de magnitude,

(17)

(18)

(19)

(20)

(21)

em que ${\textstyle p_{1},\,p_{2},\,p_{3}}$ e ${\textstyle p_{4}}$ são polinômios em relação ao parâmetro de relaxamento ${\textstyle \tau }$ , descritos por,

{\begin{aligned}&{p_{1}}(\tau )=\tau ,\\&{p_{2}}(\tau )={\frac {1}{2}}-\tau ,\\&{p_{3}}(\tau )={\tau ^{2}}-\tau +{\frac {1}{6}},\\&{p_{4}}(\tau )=-{\tau ^{3}}+{\frac {3}{2}}{\tau ^{2}}-{\frac {7}{12}}\tau +{\frac {1}{24}}.\end{aligned}}

(22)

Por meio das equações (17), (18), (19), (20) e (21) determinam-se os momentos da função distribuição de equilíbrio, e a partir das equações geradas da expansão dos momentos sobre o reticulado, determinam-se as expressões da função distribuição de equilíbrio nas direções do reticulado.

3.4 Função distribuição de equilíbrio

A função distribuição de equilíbrio juntamente com a equação do Reticulado de Boltzmann (9) desempenham um papel essencial que é recuperar a equação macroscópica do fluido. A seguir, recupera-se a equação da onda cinemática unidimensional (5) e determina-se a sua função distribuição de equilíbrio.

A equação da onda cinemática é um caso especial da equação de Burgers. A equação de Burgers é uma equação de convecção-difusão não-linear que representa um modelo simplificado das equações de Navier-Stokes. É utilizada para testar a eficiência de vários esquemas numéricos e para descrever o escoamento por meio de uma onda de choque que viaja em um fluido viscoso [40].

A equação de Burgers unidimensional pode ser escrita da seguinte forma [22,40],

(23)

onde ${\textstyle v({\vec {x}},t)}$ é uma função que pode representar uma quantidade real como a altura, a vazão ou a área da seção transversal, ${\textstyle \gamma }$ é o coeficiente de difusão, ${\textstyle \beta }$ e ${\textstyle m}$ são parâmetros a serem determinados, e ${\textstyle F}$ representa o efeito externo ou força externa. Assume-se ${\textstyle \gamma }$ nulo para obter a equação da onda cinemática, modelo no qual não é considerada a difusão.

Assumindo que a função distribuição ${\textstyle f_{\alpha }}$ está próxima do equilíbrio localmente, devem ocorrer as seguintes condições de conservação,

(24)

em consequência, as somas sobre a função distribuição de não equilíbrio têm as seguintes restrições,

(25)

Combinando as equações (17) a (21), obedecendo as restrições em (24) e (25), em vista das diferentes escalas de tempo e da equação de Burgers, equação (23), determinam-se os momentos da função distribuição de equilíbrio de zero a quarta ordem, descritos na equação (26). Neste artigo, assume-se que a força externa não varia no espaço e no tempo, portanto são nulas as derivadas nos termos que envolvem ${\textstyle g_{\alpha }}$ do lado direito das equações (20) e (21)

{\begin{aligned}&\sum \nolimits _{\alpha }{f_{\alpha }^{(0)}}=v,\\&\sum \nolimits _{\alpha }{{{\vec {e}}_{\alpha }}f_{\alpha }^{(0)}}=\beta {v^{m}},\\&\sum \nolimits _{\alpha }{{\vec {e}}_{\alpha }^{\,2}f_{\alpha }^{(0)}}={\dfrac {{\beta ^{2}}{m^{2}}}{2m-1}}{v^{2m-1}},\\&\sum \nolimits _{\alpha }{{\vec {e}}_{\alpha }^{\,3}f_{\alpha }^{(0)}}={\dfrac {{\beta ^{3}}{m^{3}}}{3m-2}}{v^{3m-2}},\\&\sum \nolimits _{\alpha }{{\vec {e}}_{\alpha }^{\,4}f_{\alpha }^{(0)}}={\dfrac {{\beta ^{4}}{m^{4}}}{4m-3}}{v^{4m-3}}.\end{aligned}}

(26)

Desenvolvendo as equações em (26) nas direções das velocidades do reticulado D1Q5, obtém-se,

{\begin{bmatrix}1&1&1&1&1\\0&-e&-2e&e&2e\\0&e^{2}&4e^{2}&e^{2}&4e^{2}\\0&-e^{3}&-8e^{3}&e^{3}&8e^{3}\\0&e^{4}&16e^{4}&e^{4}&16e^{4}\end{bmatrix}}{\begin{bmatrix}f_{0}\\f_{1}\\f_{2}\\f_{3}\\f_{4}\end{bmatrix}}={\begin{bmatrix}v\\\beta v^{m}\\{\dfrac {{\beta ^{2}}{m^{2}}}{2m-1}}{v^{2m-1}}\\{\dfrac {{\beta ^{3}}{m^{3}}}{3m-2}}{v^{3m-2}}\\{\dfrac {{\beta ^{4}}{m^{4}}}{4m-3}}{v^{4m-3}}\end{bmatrix}}

(27)

A resolução do sistema formado pelas equações em (27) determina a função distribuição de equilíbrio,

{\begin{aligned}&f_{1}^{eq}={\frac {1}{6}}\left({-4{C_{1}}+4{C_{2}}+{C_{3}}-{C_{4}}}\right),\\&f_{2}^{eq}={\frac {1}{24}}\left({2{C_{1}}-{C_{2}}-2{C_{3}}+{C_{4}}}\right),\\&f_{3}^{eq}={\frac {1}{6}}\left({4{C_{1}}+4{C_{2}}-{C_{3}}-{C_{4}}}\right),\\&f_{4}^{eq}={\frac {1}{24}}\left({-2{C_{1}}-{C_{2}}+2{C_{3}}+{C_{4}}}\right),\\&f_{0}^{eq}=v-\left({f_{1}^{eq}+f_{2}^{eq}+f_{3}^{eq}+f_{4}^{eq}}\right),\end{aligned}}

(28)

onde,

(29)

Tomando o valor médio de ${\textstyle F}$ em cada direção de velocidade do reticulado D1Q5, o termo de força é determinado pela expressão,

(30)

em que ${\textstyle F}$ pode representar a precipitação quando se trata do escoamento na superfície da bacia, ou representar a contribuição lateral quando se trata do escoamento em canais.

3.5 Condições de contorno

A representação adequada das características físicas do problema, por meio das condições iniciais e de contorno, é um fator crucial para a estabilidade e precisão das simulações. A Figura 2 ilustra uma malha unidimensional, com as direções de velocidade do reticulado D1Q5, para discretização do domínio que representa o escoamento na superfície e nos canais da bacia hidrográfica.

Figura 2. Malha computacional

As condições de contorno macroscópicas para ambos o escoamento superficial e o escoamento em canal são conhecidas na fronteira esquerda, porém desconhecidas na fronteira direita. O comportamento macroscópico nas fronteiras é alcançado pela correta aplicação das funções de distribuição de partículas. Com esse objetivo, as funções de distribuição de equilíbrio são utilizadas na fronteira esquerda, e para a fronteira direita, as funções de distribuição são determinadas por extrapolação [22].

No reticulado D1Q5 são aplicadas condições de contorno nos dois primeiros e nos dois últimos nós da malha [21]. Na fronteira esquerda, tem-se,

{\begin{aligned}&{\bar {\bar {f_{0}}}}(x_{0})=f_{0}^{eq}(x_{0}),\\&{\bar {\bar {f_{1}}}}(x_{0})=f_{1}^{eq}(x_{0}),\\&{\bar {\bar {f_{2}}}}(x_{0})=f_{2}^{eq}(x_{0}),\\&{\bar {\bar {f_{3}}}}(x_{0})=f_{3}^{eq}(x_{0}),\\&{\bar {\bar {f_{4}}}}(x_{0})={v^{*}}-({\bar {\bar {f_{0}}}}(x_{0})+{\bar {\bar {f_{1}}}}(x_{0})+{\bar {\bar {f_{2}}}}(x_{0})+{\bar {\bar {f_{3}}}}(x_{0})),\end{aligned}}

{\begin{aligned}&{\bar {\bar {f_{0}}}}(x_{1})={\bar {f_{0}}}(x_{1}),\\&{\bar {\bar {f_{1}}}}(x_{1})={\bar {f_{1}}}(x_{2}),\\&{\bar {\bar {f_{2}}}}(x_{1})={\bar {f_{2}}}(x_{3}),\\&{\bar {\bar {f_{3}}}}(x_{1})={\bar {f_{3}}}(x_{0}),\\&{\bar {\bar {f_{4}}}}(x_{1})=2{\bar {\bar {f_{4}}}}(x_{2})-{\bar {\bar {f_{4}}}}(x_{3}).\end{aligned}}

(31)

Na fronteira direita, aplica-se,

{\begin{aligned}&{\bar {\bar {f_{0}}}}(x_{n-1})={\bar {f_{0}}}(x_{n-1}),\\&{\bar {\bar {f_{3}}}}(x_{n-1})={\bar {f_{3}}}(x_{n-2}),\\&{\bar {\bar {f_{4}}}}(x_{n-1})={\bar {f_{4}}}(x_{n-3}),\\&{\bar {\bar {f_{1}}}}(x_{n-1})={\bar {f_{1}}}(x_{n}),\\&{\bar {\bar {f_{2}}}}(x_{n-1})=2{\bar {\bar {f_{2}}}}(x_{n-2})-{\bar {\bar {f_{2}}}}(x_{n-3}),\end{aligned}}

{\begin{aligned}&{\bar {\bar {f_{0}}}}(x_{n})={\bar {f_{0}}}(x_{n}),\\&{\bar {\bar {f_{3}}}}(x_{n})={\bar {f_{3}}}(x_{n-1}),\\&{\bar {\bar {f_{4}}}}(x_{n})={\bar {f_{4}}}(x_{n-2}),\\&{\bar {\bar {f_{1}}}}(x_{n})=2{\bar {\bar {f_{1}}}}(x_{n-1})-{\bar {\bar {f_{1}}}}(x_{n-2}),\\&{\bar {\bar {f_{2}}}}(x_{n})=2{\bar {\bar {f_{2}}}}(x_{n-1})-{\bar {\bar {f_{2}}}}(x_{n-2}),\end{aligned}}

(32)

em que ${\textstyle {\bar {\bar {f}}}}$ é a função distribuição de velocidade das partículas depois da propagação, ${\textstyle {\bar {f}}}$ é a função distribuição de velocidade das partículas antes da propagação, ${\textstyle v^{*}}$ é a variável macroscópica na fronteira, podendo ser substituída pela profundidade do escoamento ${\textstyle h}$ ou pela área da seção transversal ${\textstyle A}$ .

No modelo computacional, as características físicas em cada trecho de canal são mantidas constantes durante a simulação, isto é, não são alterados os valores do comprimento, largura, declividade e coeficiente de Manning do segmento de canal. O que pode implicar em seções transversais diferentes na junção dos canais, levando a necessidade de implementação de condição de contorno. Para que haja a conservação de massa, a vazão de um segmento de canal a montante para um segmento a jusante deve permanecer a mesma. Considerando a Figura 3, e o cálculo da vazão em uma seção transversal dada pela fórmula de Manning, a vazão no início do canal 2 é determinada por,

(33)

em que ${\textstyle P=b+2h}$ é o perímetro molhado do canal retangular de largura ${\textstyle b}$ .

Figura 3. Junção de canais das sub-bacias

Escrevendo o perímetro molhado em termos da área molhada e da largura do canal, e substituindo na equação (33), tem-se,

(34)

O valor de ${\textstyle A_{2}}$ é calculado pelo método iterativo de Newton-Raphson e considerado no LBM por meio da condição de contorno da equação (31), onde substitui-se ${\textstyle v^{*}}$ por ${\textstyle A_{2}}$ . Utiliza-se a variável área ao invés da vazão, pois a função distribuição de equilíbrio é dependente da área, e por meio da área obtém-se a vazão em uma seção transversal. Quanto ao escoamento na superfície da bacia, ${\textstyle v^{*}}$ é sempre nula para cumprir a condição de contorno ${\textstyle h(0,t)=0}$ em todos os passos de tempo.

3.6 Estabilidade

De acordo com Sterling e Chen (1996) [41], não é possível garantir a estabilidade do LBM e isto ocorre devido a grande quantidade de parâmetros que impedem a sua caracterização completa. No entanto, pode-se seguir algumas condições necessárias, que, se forem satisfeitas simultaneamente, obtém-se a estabilidade da simulação [2]. A primeira expressão de (35) relaciona o parâmetro de relaxamento ${\textstyle \tau }$ e a viscosidade cinemática ${\textstyle \nu }$ , dessa relação resulta ${\textstyle \tau >1/2}$ [41]. A magnitude da velocidade física ${\textstyle {\vec {v}}}$ resultante do fluido deve ser menor do que a velocidade no reticulado, conforme a segunda expressão. Na terceira expressão tem-se a condição para celeridade e na quarta a restrição para escoamentos subcríticos [2]

(35)

Como o LBM é um método intrinsecamente viscoso [42], utiliza-se com frequência, na simulação numérica, uma viscosidade mais elevada do que a viscosidade cinemática do fluido e os resultados estão em concordância com as simulações disponíveis na literatura [2,5,42,43].

4. Estudo de caso: Escoamento em uma bacia hidrográfica

A bacia hidrográfica consiste em uma área bem definida de captação natural da água de precipitação que faz convergir o escoamento até um curso de água. Na bacia hidrográfica, o escoamento é dividido em superficial e em canal. No escoamento superficial a água desloca-se pela superfície da bacia até encontrar um curso de água. Este deslocamento ocorre como resultado da água precipitada que não foi interceptada pela cobertura vegetal e pela parte que não infiltrou no solo. Devido a grande heterogeneidade espacial, a superfície da bacia é representada por planos onde ocorre o escoamento superficial de pouca profundidade. A contribuição lateral nos canais é dada principalmente pela precipitação que ocorre sobre cada plano [33].

Dependendo do tamanho e das características da bacia hidrográfica, pode-se representá-la por uma forma em “V”, segmentada por dois planos e um segmento de canal, como mostra a Figura 4. Neste estudo, utilizam-se várias destas formas para representar as sub-bacias que compõem a superfície da bacia hidrográfica.

Figura 4. Representação de uma sub-bacia

4.1 Escoamento na bacia hidrográfica do rio Chopim

A bacia do rio Iguaçu (Figura 5), com superfície aproximada de ${\textstyle 70.800}$ km ${\textstyle ^{2}}$ , tem grande importância na geração de energia por meio de usinas hidrelétricas e pela grande demanda de recursos hídricos no abastecimento público, industrial e agrícola. A região ainda é caracterizada por um grande número de Unidades de Conservação Ambiental, e abriga importantes corredores de biodiversidade [44].

Figura 5. Localização da bacia hidrográfica do rio Iguaçu

Fonte: Franco (2017) [45]

Dentro da bacia do rio Iguaçu está localizada a Bacia do rio Chopim que abrange uma área de ${\textstyle 7.500}$ km ${\textstyle ^{2}}$ e desenvolve-se basicamente no sentido sudoeste-noroeste, aproximadamente entre os paralelos ${\textstyle 25\,^{\circ }30'}$ e ${\textstyle 26\,^{\circ }40'}$ de latitude sul e os meridianos ${\textstyle 51\,^{\circ }30'}$ e ${\textstyle 53\,^{\circ }20'}$ de longitude oeste. Seus limites estão entre as bacias do rio Uruguai ao sul e com bacias de afluentes do rio Iguaçu nas outras direções [46].

O rio Chopim é o mais importante afluente da bacia do baixo rio Iguaçu. Suas nascentes estão localizadas em altitudes que superam os 1.200 m. A extensão total do curso principal do rio é da ordem de 450 km. Da nascente do rio até em torno do km 209 a inclinação do leito é de aproximadamente 2,9 m/km, e desde este ponto até a foz, a inclinação é de aproximadamente 1,1 m/km [47].

Considera-se nesse estudo, um trecho do rio Chopim localizado entre suas nascentes no município de Palmas, passando pela estação pluviométrica Salto Claudelino até a estação fluviométrica Porto Palmeirinha. Este trecho do rio tem comprimento de aproximadamente ${\textstyle 281}$ km com área de drenagem na superfície da bacia em torno de ${\textstyle 3.205}$ km ${\textstyle ^{2}}$ . As características das estações estão listadas na Tabela A.

A Figura 6 ilustra os principais rios da bacia do rio Iguaçu, suas estações pluviométricas, meteorológicas e fluviométricas, e destaca a localizacão da bacia do rio Chopim.

Figura 6. Localização das estações pluviométricas, meteorológicas e fluviométricas da bacia hidrográfica do rio Iguaçu

Fonte: Adaptado de Uda (2016) [48]

A Figura 7 mostra a bacia do rio Chopim com as divisões de sub-bacias conforme Kaviski e Gonçalves (1996) [49] definidas de acordo com as características físicas da bacia.

Figura 7. Bacia hidrográfica do rio Chopim

Fonte: Adaptado de Kaviski e Gonçalves (1996) [49]

5. Resultados e discussões

Levando-se em conta a delimitação das sub-bacias e os trechos de canais apresentados na Figura 7, fez-se a segmentação da bacia do rio Chopim em nove sub-bacias até a estação Porto Palmeirinha. Cada sub-bacia é composta por um trecho de canal e dois planos laterais de escoamento, com exceção da primeira sub-bacia que é representada por apenas um plano lateral. A numeração das sub-bacias foi mantida a mesma que a utilizada na sequência para os trechos de canais. Essa numeração indica o início do escoamento até a estação fluviométrica Porto Palmeirinha. Os dados das características físicas da segmentação da bacia são apresentados na Tabela B para os segmentos planos, e na Tabela C para os segmentos de canais.

Nos trechos de canais não são conhecidos os valores da declividade e do coeficiente de Manning. Para complementar a descrição das características físicas das sub-bacias hidrográficas, listam-se na Tabela C os parâmetros ${\textstyle {\bar {\alpha }}}$ e ${\textstyle {\bar {\beta }}}$ , que descrevem as curvas de descarga. De acordo com Kaviski e Gonçalves (1996) [49], a relação entre área da seção transversal e vazão é descrita pela equação,

(36)

onde, ${\textstyle {\bar {\alpha }}}$ e ${\textstyle {\bar {\beta }}}$ são estimados pelo método dos mínimos quadrados utilizando valores do nível da água e da descarga, e também da elevação do fundo do canal e a área, medidos em seções transversais.

Ao reescrever a equação (36) como vazão em função da área, tem-se,

(37)

dessa forma, os parâmetros ${\textstyle \beta _{c}}$ e ${\textstyle m_{c}}$ da equação (5) são determinados por,

(38)

Na junção de segmentos de canais (Figura 3), a área da seção transversal é encontrada diretamente pela equação (36).

Por não estarem disponíveis os dados de declividade e coeficiente de Manning dos trechos de canais da bacia do rio Chopim, houve a necessidade de adaptar o modelo para utilizar os parâmetros ${\textstyle {\bar {\alpha }}}$ e ${\textstyle {\bar {\beta }}}$ . A vantagem em utilizar esses parâmetros é a de aplicar o LBM em um canal com seção transversal estimada pelas características físicas naturais do rio.

Os parâmetros da simulação no LBM para os planos são ${\textstyle \Delta x=10}$ m, ${\textstyle \Delta t=10}$ s e velocidade na malha ${\textstyle e=1}$ m/s. Com relação aos canais, ${\textstyle \Delta x=50}$ m, ${\textstyle \Delta t=10}$ s e velocidade na malha ${\textstyle e=5}$ m/s. Em ambos, o parâmetro de relaxamento utilizado é ${\textstyle \tau =12}$ . O número de iterações realizadas foi de ${\textstyle 90.000}$ , o que corresponde a um tempo de simulação de ${\textstyle 250}$ horas ou ${\textstyle 10,42}$ dias.

O período de chuva analisado na bacia é dado no mês 07 de 1987 nos dias 07 e 08. Neste mês não houve chuva além dos dias citados. A Tabela D contém a intensidade média de precipitação em cada sub-bacia. Para se obter uma simulação mais próxima da realidade, considerou-se uma vazão inicial que ocorre no rio (escoamento de base), em cada trecho de canal, antes de ocorrer a precipitação. Essa vazão entra no modelo computacional como uma condição inicial. Também é considerada como uma condição de contorno ao ser adicionada na contribuição lateral dos planos, em cada passo de tempo. A Tabela E, lista as vazões iniciais na seção transversal a montante e a jusante em cada trecho de canal. Utilizou-se interpolação linear para determinar os valores de vazões iniciais a cada ${\textstyle \Delta x=50}$ m no comprimento dos canais no instante inicial da simulação. Nos demais passos de tempo, o escoamento de base é incluído pela equação (39),

(39)

em que ${\textstyle Q_{j}}$ e ${\textstyle Q_{m}}$ são, respectivamente, a vazão a jusante e a vazão a montante em um segmento de canal de comprimento ${\textstyle L}$ . A quantidade ${\textstyle q_{base}}$ é somada a contribuição lateral e uniformemente distribuída no comprimento do canal.

Os resultados obtidos pelo LBM são comparados com dados observados nas duas estações pertencentes ao trecho analisado. Nas Figuras 8 e 9, tem-se os hidrogramas simulado e observado nas estações Salto Claudelino e Porto Palmeirinha referentes a precipitação ocorrida no mês 07/1987.

Figura 8. Hidrograma na sub-bacia 5 - Estação Salto Claudelino

Figura 9. Hidrograma na sub-bacia 9 - Estação Porto Palmeirinha

Os resultados obtidos mostram boa concordância entre a simulação por meio do LBM e os dados observados nas estações, tendo em vista que os dados observados foram coletados a cada 12 horas por leitura feita de forma visual em réguas ou escalas linimétricas fixadas nas margens dos rios para registro das variações no nível da água. Essas leituras são normalmente realizadas por moradores da região. Ainda, assume-se que a precipitação é igualmente distribuída nas últimas 12 horas.

Para quantificar a concordância entre os valores observados e simulados pelo LBM, considerou-se o coeficiente de correlação linear de Pearson ${\textstyle r}$ [50],

r={\dfrac {n\sum \limits _{i=1}^{n}{{X_{i}}{Y_{i}}-\sum \limits _{i=1}^{n}{X_{i}}\sum \limits _{i=1}^{n}{Y_{i}}}}{\sqrt {\left[{n\sum \limits _{i=1}^{n}{X_{i}^{2}-{{\left({\sum \limits _{i=1}^{n}{X_{i}}}\right)}^{2}}}}\right]\left[{n\sum \limits _{i=1}^{n}{Y_{i}^{2}-{{\left({\sum \limits _{i=1}^{n}{Y_{i}}}\right)}^{2}}}}\right]}}},

(40)

onde X e Y são as variáveis em estudo, respectivamente, vazão observada e vazão simulada, ${\textstyle n}$ é o número de pares das observações.

O coeficiente de correlação linear ${\textstyle r}$ é adimensional e tem variação no intervalo ${\textstyle [-1,1]}$ . Se ${\textstyle r=-1}$ tem-se a correlação linear negativa perfeita, enquanto que se ${\textstyle r=1}$ , tem-se a correlação linear positiva perfeita. Para ${\textstyle r=0}$ não há correlação linear entre as variáveis. Esse coeficiente pode ser avaliado quantitativamente como descrito na Tabela 1.

Tabela 1. Classificação do coeficiente linear de Pearson

Intervalo	Classificação
${\textstyle 0<\|r\|<0,3}$	Existe fraca correlação linear
$0,3<\|r\|<0,6$	Existe moderada correlação linear
$0,6<\|r\|<0,9$	Existe forte correlação linear
$0,9<\|r\|<1,0$	Existe correlação linear muito forte

Fonte: Marques e Marques [50]

As Figuras 10 e 11 mostram os gráficos de dispersão das vazões observadas pelas vazões simuladas pelo LBM. A reta tracejada, em ambos os gráficos, representa o ajuste perfeito com coeficiente ${\textstyle r=1}$ .

Figura 10. Dispersão das vazões observadas e simuladas na sub-bacia 5 - Estação Salto Claudelino

Figura 11. Dispersão das vazões observadas e simuladas na sub-bacia 9 - Estação Porto Palmeirinha

O coeficiente de correlação linear de Pearson calculado na sub-bacia 5 foi de ${\textstyle r=0,88}$ e na sub-bacia 9, de ${\textstyle r=0,97}$ . Conforme a Tabela 1, existe forte correlação linear entre as variáveis na sub-bacia 5 e correlação linear muito forte para a sub-bacia 9. Conclui-se que há boa concordância entre as vazões observadas e simuladas. Nota-se na Figura 10 que um ponto está distante da reta, o que faz o coeficiente ${\textstyle r}$ ser um pouco menor.

Para fazer a verificação da significância estatística, foi aplicado o teste de hipóteses com a distribuição t de Student, supondo que os desvios padrões populacionais são desconhecidos e com variâncias equivalentes. A hipótese a ser testada é ${\textstyle H_{0}:\mu _{1}-\mu _{2}=0}$ , isto é, não existe diferença significativa entre as médias populacionais das vazões observadas e simuladas. Os resultados obtidos, para as sub-bacias 5 e 9, estão resumidos na Tabela F.

Os testes realizados mostram que a estatística do teste ${\textstyle t}$ é menor do que o valor do crítico bi-caudal, portanto aceita-se a hipótese ${\textstyle H_{0}}$ , ou seja, não existe diferença significativa, ao nível de significância de ${\textstyle 5\%}$ entre as vazões observadas e vazões simuladas.

6. Conclusão

Esta pesquisa utilizou o LBM com aperador de colisão BGK para obter a simulação numérica do escoamento superficial da bacia do rio Chopim. Computacionalmente a bacia foi representada pela junção de nove sub-bacias dadas por planos laterais e um segmento de canal. A simulação produziu o hidrograma total o qual fornece a variação da vazão em relação ao tempo correspondente a resposta da bacia hidrográfica após a precipitação efetiva do evento simulado. Para que este resultado fosse alcançado, determinou-se duas funções de distribuição de equilíbrio obtidas por meio da expansão de Chapmann-Enskog em escalas de tempo e utilizando o reticulado D1Q5, uma adequada para o escoamento na superfície da bacia e outra para o canal principal, obtendo a profundidade da água na superfície da bacia e área da seção transversal nos canais. Também foi necessário determinar a condição de contorno na passagem de fluxo de um a sub-bacia para outra, levando em consideração a conservação da massa.

Obteve-se uma simulação mais próxima da realidade com a consideração de escoamento de base, em cada trecho de canal, antes que ocorra a precipitação. Essa vazão entra no modelo computacional como uma condição inicial. Também é considerada como uma condição de contorno ao ser adicionada na contribuição lateral dos planos, em cada passo de tempo. Utilizou-se interpolação linear para determinar os valores de vazões iniciais a cada 50 m no comprimento dos canais no instante inicial da simulação. Nos demais passos de tempo, o escoamento de base é incluído como uma condição de contorno. Os resultados obtidos pelo LBM foram comparados com dados observados nas duas estações pertencentes ao trecho analisado e mostram boa concordância, tendo em vista que os dados observados foram coletados a cada 12 horas por leituras feitas de forma visual em réguas fixadas nas margens dos rios para registro das variações no nível da água. Ainda, considera-se que a precipitação é igualmente distribuída nas últimas 12 horas.

Verificou-se que LBM é adequado e preciso na simulação do escoamento em uma bacia hidrográfica representado pelas equações de onda cinemática. Além disso, as condições iniciais e de contorno são de fácil implementação.

Referências

[1] Gribbin J.E. Introduction to hydraulics and hydrology with applications for stormwater management. 3rd Edition, Cengage Learning, USA, 2008.

[2] Zhou J.G. Lattice Boltzmann method for shallow water flows. Springer, New York, 2004.

[3] Chen S., Doolen G.D. Lattice Boltzmann method for fluid flows. Annual Review of Fluid Mechanics, 30(1):329-364, 1998.

[4] Zhou J.G. A lattice Boltzmann model for the shallow water equations. Computer Methods in Applied Mechanics and Engineering, 191:3527-3539, 2002.

[5] Zhou J.G., Liu H., Shafiai S., Peng Y., Burrows R. Lattice Boltzmann method for open-channel flows. Engineering and Computational Mechanics, 163:243-249, 2010.

[6] van Thang P., Chopard B., Lefèvre L., Ondo D.A., Mendes E. Study of the 1D lattice Boltzmann shallow water equation and its coupling to build a canal network. Journal of Computational Physics, 229(19):7373-7400, 2010.

[7] Peng Y., Zhou J.G., Burrows R. Modelling solute transport in shallow water with the lattice Boltzmann method. Computers & Fluids, 50:181-188, 2011.

[8] Peng Y., Zhou J.G., Zhang J.M., Burrows R. Modeling moving boundary in shallow water by lbm. International Journal of Modern Physics, 24(1):1250094, 2013.

[9] Budinski L. Mrt lattice Boltzmann method for 2d flows in curvilinear coordinates. Computers & Fluids, 94:288?301, 2014.

[10] Peng Y., Zhou J.G., Zhang J.M., Liu H. Lattice Boltzmann modeling of shallow water flows over discontinuous beds. International Journal for Numerical Methods in Fluids, 75:608-619, 2014.

[11] La Rocca M., Montessori A., Prestininzi P., Succi S. A multispeed discrete Boltzmann model for transcritical 2d shallow water flows. Journal of Computational Physics, 284:117-132, 2015.

[12] Li S., Huang P., Li J. A modified lattice Boltzmann model for shallow water flows over complex topography. International Journal for Numerical Methods in Fluids, 77(8):441-458, 2015.

[13] Prestininzi P., Montessori A., La Rocca M., Sciortino G. Simulation of arrested salt wedges with a multi-layer shallow water lattice Boltzmann model. Advances in Water Resources, 96:282-289, 2016.

[14] Zhang C.-Z., Cheng Y.-G., Wu J.-Y., Diao W. Lattice Boltzmann simulation of the open channel flow connecting two cascaded hydropower stations. Journal of Hydrodynamics, 28(3):400-410, 2016.

[15] Galina V., Cargnelutti J., Kaviski E., Gramani L.M., Lobeiro A.M. Simulação de onda de maré por meio do método do reticulado de Boltzmann. In: I Simpósio de Métodos Numéricos em Engenharia, Universidade Federal do Paraná, Curitiba, Brasil, 2016.

[16] Zhao Z.-M., Huang P., Li S.-T. Lattice Boltzmann model for shallow water in curvilinear coordinate grid. Journal of Hydrodynamics, 29(2):251-260, 2017.

[17] De Rosis A. A central moments-based lattice Boltzmann scheme for shallow water equations. Computer Methods in Applied Mechanics and Engineering, 319:379-392, 2017.

[18] Cargnelutti J., Galina V., Kaviski E., Gramani L.M., Lobeiro A.M. Two-dimensional numerical simulation of channel flow with submerged obstacles using the lattice Boltzmann method. Revista Internacional de Métodos Numéricos para Cálculo y Diseño en Ingeniería, 34(1), 20, 2018.

[19] Cargnelutti J., Galina V., Kaviski E., Gramani L.M., Lobeiro A.M. Simulation of the two-dimensional flow of the initiation channel of the Itaipu hydroelectric power plant by the lattice Boltzmann method. Revista Internacional de Métodos Numéricos para Cálculo y Diseño en Ingeniería, 34(1), 25, 2018.

[20] Yang G., Kwok C., Sobral Y. The effects of bed form roughness on total suspended load via the lattice Boltzmann method. Applied Mathematical Modelling, 63:591-610, 2018.

[21] Zhang X., Feng J., Zhang D., Liu N. Comparison of lattice Boltzmann method and preissmann implicit difference method in application to overland flow. Nongye Jixie Xuebao/Transactions of the Chinese Society of Agricultural Machinery, 2014.

[22] Zhang X., Feng J., Yang T. Lattice Boltzmann method for overland flow studies and its Eeperimental validation. Journal of Hydraulic Research, 53(5):561-575, 2015.

[23] Liu N. The numerical simulation of one-dimensional overland flow by lattice Boltzmann method. In: 5th International Conference on Advanced Design and Manufacturing Engineeringe, China, 2015.

[24] Galina V., Cargnelutti J., Kaviski E., Gramani L.M., Lobeiro A.M. Application of lattice Boltzmann method for surface runoff in watershed. Revista Internacional de Métodos Numéricos para Cálculo y Diseño en Ingeniería, 34(1), 10, 2018.

[25] Liu H., Wang H., Liu S., Hu C., Ding Y., Zhang J. Lattice Boltzmann method for the Saint-Venant equations. Journal of Hydrology, 524:411-416, 2015.

[26] Benzi R., Succi S., Vergassola M. The lattice Boltzmann equation: Theory and applications. Physics Reports, 222(3):145-197, 1992.

[27] Nestler B., Aksi A., Selzer M. Combined lattice Boltzmann and phase-field simulations for incompressible fluid flow in porous media. Mathematics and Computers in Simulation, 80(7):1458-1468, 2010.

[28] Sobieski W. Numerical investigations of tortuosity in randomly generated pore structures. Mathematics and Computers in Simulation, 166:1-20, 2019.

[29] Andaluz M.E., Galarza V.V., Vera A.R. On hydraulic tortuosity variations due to morphological considerations in 2d porous media by using the lattice Boltzmann method. Mathematics and Computers in Simulation, 169:74-87, 2020.

[30] Qian Y.H., D'Humières D., Lallemand P. Lattice BGK models for Navier-Stokes equation. Europhysics Letters, 17(6):479-484, 1992.

[31] Porto R.M. Hidráulica básica. 4th Edition, EESC-USP, São Carlos, 2006.

[32] Chaudhry M.H. Open-channel flow. Springer, New York, USA, 2008.

[33] Tucci C.E.M. Modelos hidrológicos. Editora da Universidade Federal do Rio Grande do Sul, Porto Alegre, RS - Brasil, 1998.

[34] Miller J.E. Basic concepts of kinematic-wave models. United States Government Printing office, Washington, 1984.

[35] Golbert D.R. Método de lattice Boltzmann em hemodinâmica computacional: interações fluido-estrutura e modelos acoplados 1D-3D. Ph.D. Thesis, Laboratório Nacional de Computação Científica, Petrópolis, RJ - Brasil, 2013.

[36] Bhatnagar P.L., Gross E.P., Krook M. A model for collision processes in gases. I. Small amplitude processes in charged and neutral one-component systems. Physical Review, 94(3):511-525, 1954.

[37] Guo Z., Zheng C., Shi B. Discrete lattice effects on the forcing term in the lattice Boltzmann method. Physical Review E, 65(4):1-6, 2002.

[38] Zhang J., Yan G. Lattice Boltzmann method for one and two-dimensional burgers equation. Physica A: Statistical Mechanics and its Applications, 387(19–20):4771-4786, 2008.

[39] Chapman S., Cowling T.G. The mathematical theory of non-uniform gases. An account of the kinetic theory of viscosity, thermal conduction and diffusion in gases. 3rd Edition, Cambridge University Press, London, 1970.

[40] Zhang J., Yan G. A lattice Boltzmann model for the Burgers-Fisher equation. Chaos: An Interdisciplinary Journal of Nonlinear Science, 20:023129, 2010.

[41] Sterling J.D., Chen S. Stability analysis of lattice Boltzmann methods. Journal of Computational Physics, 123:196-206, 1996.

[42] La Rocca M., Adduce C., Sciortino G. Development of a lattice Boltzmann method for two-layered shallow-water flow. International Journal for Numerical Methods in Fluids, 70:1048-1072, 2012.

[43] Zhou J.G., Liu H. Determination of bed elevation in the enhanced lattice Boltzmann method for the shallow-water equations. Physical Review E, 88(2):1-6, 2013.

[44] Secretaria de Estado do Meio Ambiente e Recursos Hídricos SEMA, Bacias hidrográficas do Paraná, Brasil, 2010.

[45] Franco A.C.L. Calibração do modelo swat com evapotranspiração proveniente de sensoriamento remoto e vazão observada. Master's thesis, Florianópolis - SC, 2017.

[46] Chavasse D.I., Seoane R.S. Assessing and predicting the impact of el Nino Southern Oscillation (ENSO) events on runoff from the Chopim river basin, Brazil. Hydrological Processes, 23(22):3261-3266, 2009.

[47] Burian P.P. Do estudo de impacto embiental à avaliação ambiental estratégica - ambivalências do processo de licenciamento ambiental do setor elétrico. Ph.D. Thesis, Universidade Estadual de Campinas, Campinas, SP - Brasil, Fevereiro 2006.

[48] Uda P.K. Evapotranspiração real da bacia do rio iguaçu por meio do modelo metric. Ph.D. Thesis, Universidade Federal de Santa Catarina, Florianópolis, 2016.

[49] Kaviski E., Gonçalves L.F.A. Análise hidrológica e matemática de operação de reservatórios tech. rep. Centro de Hidráulica e Hidrologia Professor Parigot de Souza - CEHPAR, Curitiba - PR, 1996.

[50] Marques J.M., Marques M.A.M. Estatística básica para os cursos de engenharia. Domínio do Saber, 2005.

Apêndice

Tabela A. A estações hidrométricas no rio Chopim

Nome	Código	Município	Latitude	Longitude	Altitude (m)	Área de contribuição (km ${\textstyle ^{2}}$ )
Porto Palmeirinha	65927000	Coronel Vivida	-26,03	-52,63	501	3.390
Salto Claudelino	65925000	Clevelândia	-26,28	-52,30	797	1.660

Fonte: Franco [45]

Tabela B. Caracterização física dos planos das 9 sub-bacias do rio Chopim - pd e pe são, respectivamente, os planos do lado direito e do lado esquerdo no sentido do escoamento

Sub-Bacia	Plano	Comprimento $L$ (m)	Largura $b$ (m)	Declividade $S_{0}$ (m/m)	Manning $n_{r}$ (m $^{-1/3}$ s)	Área (km $^{2}$ )
1	pd	-	-	-	-	-
1	pe	6.940	26.830	0,0065	0,1200	186,20
2	pd	3.080	47.500	0,0266	0,0840	146,30
2	pe	6.940	47.500	0,0065	0,1200	329,65
3	pd	3.300	30.170	0,0146	0,1044	99,56
3	pe	6.940	30.170	0,0065	0,1200	209,38
4	pd	2.800	26.580	0,0102	0,2832	74,42
4	pe	9.250	26.580	0,0129	0,2832	245,87
5	pd	2.800	41.670	0,0102	0,2832	116,68
5	pe	6.780	41.670	0,0022	0,1020	282,52
6	pd	2.800	13.000	0,0102	0,2832	36,40
6	pe	6.780	13.000	0,0022	0,1020	88,14
7	pd	5.180	65.000	0,0170	0,0972	336,70
7	pe	9.670	65.000	0,0085	0,0972	628,55
8	pd	10.120	22.000	0,0113	0,0984	222,64
8	pe	4.870	22.000	0,0129	0,1020	107,14
9	pd	6.720	8.180	0,0127	0,1032	54,97
9	pe	4.870	8.180	0,0129	0,1020	39,84

Fonte: Kaviski e Gonçalves [49]

Tabela C. Caracterização física dos canais das 9 sub-bacias do rio Chopim

Sub-Bacia	Comprimento do canal $L$ (m)	Parâmetro ${\bar {\alpha }}$	Parâmetro ${\bar {\beta }}$	Área de drenagem (km $^{2}$ )	Área de drenagem acumulada (km $^{2}$ )
1	26.830	0,04113596	1,59495580	186,21	186,21
2	47.500	0,85500771	0,87796837	475,86	662,07
3	30.170	1,01710784	0,86807752	308,94	971,01
4	26.580	1,45162356	0,83154649	320,33	1.291,34
5	41.670	2,55791569	0,76867062	399,15	1.690,49
6	13.000	2,47863913	0,77061510	124,52	1.815,01
7	65.000	2,32144856	0,77466167	965,27	2.780,28
8	22.000	2,14755607	0,77947050	329,82	3.110,10
9	8.180	2,09342670	0,78104723	94,78	3.204,88

Fonte: Kaviski e Gonçalves [49]

Tabela D. Intensidade média de precipitação no rio Chopim, mês 07/1987 - pd e pe são, respectivamente, os planos do lado direito e do lado esquerdo no sentido do escoamento

Sub-Bacia	Plano	07/07/1987 (m/h) $\cdot$ 10 $^{-3}$		08/07/1987 (m/h) $\cdot$ 10 $^{-3}$
Sub-Bacia	Plano	7h	19h	7h	19h
1	pd	-	-	-	-
1	pe	0	0	1,47400098	0,00174585
2	pd	0	0,65031194	2,24260054	0,00242411
2	pe	0	0	1,47400098	0,00174585
3	pd	0	0	1,67223974	0
3	pe	0	0	1,47400098	0,00174585
4	pd	0	0	0,87767508	0
4	pe	0	0	1,15756772	0
5	pd	0	0	0,87767508	0
5	pe	0	0	0,99927618	0
6	pd	0	0	0,87767508	0
6	pe	0	0	0,99927618	0
7	pd	0	0	0,90541231	0
7	pe	0	0	1,15131540	0
8	pd	0	0	1,08201453	0
8	pe	0	0	1,26797648	0,00923398
9	pd	0	0,72930031	1,47131551	0,02406490
9	pe	0	0	1,26797648	0,00923398

Fonte: Kaviski e Gonçalves [49]

Tabela E. Vazões iniciais e coeficientes de interpolação da vazão nos trechos de canais do rio Chopim, mês 07/1987

Canal	Vazão inicial na seção tranversal (m $^{3}$ /s)		Coeficientes da interpolação linear $Q(x)=ax+b$		Vazão inicial como contribuição lateral (m $^{2}$ /s) $\cdot$ 10 $^{-4}$
Canal	Montante	Jusante	$a\,\cdot 10^{-4}$	$b$
1	0	4,91134977	1,83054408	0	1,83054408
2	4,91134977	17,46209717	2,64226261	4,91134977	2,64226261
3	17,46209717	25,61030197	2,70076394	17,46209717	2,70076394
4	25,61030197	34,05889511	3,17855272	25,61030197	3,17855272
5	34,05889511	44,58636093	2,52638969	34,05889510	2,52638969
6	44,58636093	48,07283783	2,68190530	44,58636093	2,68190530
7	48,07283783	75,09886169	4,15784983	48,07283783	4,15784983
8	75,09886169	84,33335876	4,19749867	75,09886168	4,19749867
9	84,33335876	86,98714447	3,24423681	84,33335877	3,24423681

Fonte: Kaviski e Gonçalves [49]

Tabela F. Teste t de Student para as sub-bacias no rio Chopim

Quantidades estatísticas	Sub-bacia 5 Salto Claudelino		Sub-Bacia 9 Porto Palmeirinha
Quantidades estatísticas	Observado	Simulado	Observado	Simulado
Média	80,04	76,58	142,12	139,02
Variância	1.493,98	3.021,13	3.376,02	4.864,39
Número de observações	21		21
Graus de liberdade	40		40
Estatística do teste	0,2361		0,1567
t bi-caudal	2,0211		2,0211
Nível de significância	$5\%$		$5\%$

Latest revision as of 12:20, 1 October 2021

Resumo

Abstract

1. Introdução

2. Equações de escoamento

3. Método do reticulado de Boltzmann

3.1 Equação do LBM e o operador BGK

3.2 Reticulado D1Q5

3.3 Expansão multi-escala

3.4 Função distribuição de equilíbrio

3.5 Condições de contorno

3.6 Estabilidade

4. Estudo de caso: Escoamento em uma bacia hidrográfica

4.1 Escoamento na bacia hidrográfica do rio Chopim

5. Resultados e discussões

6. Conclusão

Referências

Apêndice

Document information

Document Score

Share this document

Keywords

claim authorship

@@ Line 1: / Line 1: @@
-==Simulação Numérica do Escoamento Superficial por meio do Método do Reticulado de Boltzmann utilizando Junção de Sub-bacias==
-'''V. Galina J. Cargnelutti E. Kaviski L. M. Gramani A. M. Lobeiro'''
 ==Resumo==
 Investigou-se o escoamento superficial em bacia hidrográfica por meio da simulação numérica unidimensional utilizando o método do reticulado de Boltzmann (LBM). Desenvolveu-se um modelo computacional onde a bacia hidrográfica é representada pela junção de nove sub-bacias. Para isto, foram estabelecidas duas funções de distribuição de equilíbrio obtidas por meio da expansão de Chapmann-Enskog e utilizando o reticulado D1Q5, uma adequada para o escoamento na superfície da bacia e outra para o canal principal, obtendo a profundidade da água na superfície da bacia e a área da seção transversal nos canais. Além disso, estabeleceu-se a condição de contorno na passagem de fluxo de uma sub-bacia para outra, levando em consideração a conservação da massa e, para se obter uma simulação mais próxima da realidade, considerou-se uma vazão inicial que ocorre no rio (escoamento de base), em cada trecho de canal.  Os resultados numéricos obtidos pelo LBM foram comparados com dados medidos em campo.
-Rev. Int. Mét. Num. Cálc. Dis. Ing.
+<!--Vanderlei Galina
-Vanderlei Galina
 Departamento de Matemática, Universidade Tecnológica Federal do Paraná (UTFPR)
 Rua Cristo Rei, Vila Becker CEP 85902-490 - Toledo PR, Brasil  Telefone Geral +55 (45) 3378-6863
@@ Line 24: / Line 18: @@
 Adilandri Mércio Lobeiro Departamento de Matemática Universidade Tecnológica Federal do Paraná (UTFPR) - Campo Mourão PR, Brasil
-e-mail: adilandri@gmail.com
+e-mail: adilandri@gmail.com-->
+'''Palavras-chaves''': Método do reticulado de Boltzmann, escoamento em bacias hidrográficas, escoamento superficial, modelo onda cinemática, junção de sub-bacias
-''Palavras-chaves:''
+==Abstract==
-Método do reticulado de Boltzmann. Escoamento em bacias hidrográficas. Escoamento superficial. Modelo onda cinemática. Junção de sub-bacias.
-==NUMERICAL SIMULATION OF SURFACE  FLOW THROUGH THE LATTICE  BOLTZMANN METHOD USING SUB-BASIN JUNCTION==
 The watershed surface runoff was investigated through the one-dimensional numerical simulation using the Lattice Boltzmann Method (LBM). A computational model was developed where the watershed is represented by the junction of nine sub-basins. For this, two equilibrium distribution functions were established through the Chapman-Enskog Expansion on a D1Q5 lattice, one suitable for flow on the basin surface and another for the main channel, obtaining the water depth on the basin surface and the channels cross-sectional area. In addition, the boundary condition was established in the flow passage from one sub-basin to another, taking into account the mass conservation and, in order to obtain a simulation closer to reality, it was considered an initial river flow (baseflow) of each channel stretch. The numerical results obtained by the LBM were compared with data measured in field.
-''Keywords:''
+'''Keywords''': Lattice Boltzmann method, runoff in watersheds, overland flow, kinematic wave model, sub-basin junction
-Lattice Boltzmann method. Runoff in watersheds. Overland flow. Kinematic wave model. Sub-basin junction.
+==1. Introdução==
-==1 Introdução==
 O comportamento do escoamento de um fluido pode ser, satisfatoriamente, avaliado por meio de métodos numéricos. Simulações numéricas, análises teóricas, técnicas experimentais e dados coletados, são ferramentais que se complementam na resolução de problemas em hidrodinâmica <span id='citeF-1'></span>[[#cite-1|[1]]]. Em particular, o LBM se destaca por capturar as características físicas de problemas, até mesmo detalhes microscópicos, e as apresenta de forma macroscópica <span id='citeF-2'></span>[[#cite-2|[2]]]. A abordagem do LBM é diferenciada quando comparada aos métodos numéricos tradicionais, pois fornece uma maneira indireta para a solução das equações governantes do escoamento. O LBM não utiliza discretização das equações macroscópicas, o método baseia-se em modelos microscópicos e equações que governam a cinética em um nível mesoscópico. No LBM a dinâmica macroscópica de um fluido é resultado do comportamento coletivo de partículas microscópicas e não se altera com os detalhes subjacentes referentes as interações moleculares do fluido <span id='citeF-3'></span>[[#cite-3|[3]]]. O método tem sido aplicado nas mais diversas áreas de pesquisa e tornou-se de grande potencial no estudo da dinâmica de fluidos computacional. Escoamentos envolvendo águas rasas <span id='citeF-2'></span><span id='citeF-4'></span><span id='citeF-5'></span><span id='citeF-6'></span><span id='citeF-7'></span><span id='citeF-8'></span><span id='citeF-9'></span><span id='citeF-10'></span><span id='citeF-11'></span><span id='citeF-12'></span><span id='citeF-13'></span><span id='citeF-14'></span><span id='citeF-15'></span><span id='citeF-16'></span><span id='citeF-17'></span><span id='citeF-18'></span><span id='citeF-19'></span><span id='citeF-20'></span>[[#cite-2|[2,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20]]], o modelo onda cinemática <span id='citeF-6'></span><span id='citeF-21'></span><span id='citeF-22'></span><span id='citeF-23'></span><span id='citeF-25'></span><span id='citeF-24'></span>[[#cite-6|[6,21,22,23,25,24]]] e várias outras aplicações como em meios porosos <span id='citeF-26'></span><span id='citeF-30'></span><span id='citeF-27'></span><span id='citeF-28'></span><span id='citeF-29'></span>[[#cite-26|[26,30,27,28,29]]] foram abrangidos pelo LBM e as simulações apresentaram resultados satisfatórios.
@@ Line 54: / Line 43: @@
 Este artigo está estruturado da seguinte forma. Na seção [[#2 Equações de escoamento|2]] tem-se as equações governantes do escoamento. Em seguida, na seção [[#3 Método do Reticulado de Boltzmann|3]], apresenta-se o método numérico que será utilizado e os quatro elementos principais que o caracterizam: Equação do LBM e o operador de colisão BGK, o reticulado D1Q5, a expensão multi-escala e as funções distribuição de equilíbrio. Tem-se, também, nesta seção informações sobre condições de contorno e estabilidade. O estudo de caso é apresentado na seção [[#4 Estudo de caso: Escoamento em uma bacia hidrográfica|4]] considerando as características da bacia do rio Chopim, resultados da simulação e discussões. Em seguida, na seção [[#6 Conclusão|6]], apresentam-se as conclusões deste trabalho.
-==2 Equações de escoamento==
+==2. Equações de escoamento==
 As equações de Saint-Venant unidimensionais formam um sistema de equações não-lineares composto pelas equações da continuidade [[#eq-1|(1)]] e da quantidade de movimento [[#eq-2|(2)]]. Elas governam o escoamento onde há superfície livre, com a suposição de que a componente vertical e transversal da velocidade do escoamento pode ser desprezada em relação aos componentes longitudinais <span id='citeF-32'></span>[[#cite-32|[32]]],
@@ Line 66: / Line 55: @@
 |-
 | style="text-align: center;" | <math>\frac{{\partial Q}}{{\partial x}} + \frac{{\partial A}}{{\partial t}} = {0}  </math>
+|}
 | style="width: 5px;text-align: right;white-space: nowrap;" | (1)
+|}
+{| class="formulaSCP" style="width: 100%; text-align: left;"
+|-
+|
+{| style="text-align: left; margin:auto;width: 100%;"
 |-
 | style="text-align: center;" | <math>  \frac{{\partial Q}}{{\partial t}} + \frac{\partial }{{\partial x}}\left({\frac{{{Q^2}}}{A}} \right)+ gA\frac{{\partial h}}{{\partial x}} = gA({S_o} - {S_f}),   </math>
@@ Line 108: / Line 104: @@
 {| style="text-align: left; margin:auto;width: 100%;"
 |-
-| style="text-align: center;" | <math>\frac{{\partial A}}{{\partial t}} + \frac{{\partial Q}}{{\partial x}} = {q_l}</math>
+| style="text-align: center;" | <math>\frac{{\partial A}}{{\partial t}} + \frac{{\partial Q}}{{\partial x}} = {q_l}\quad ; </math> <math> \quad   Q = {\beta _c}{A^{{m_c}}}     , </math>
-|-
-| style="text-align: center;" | <math>   Q = {\beta _c}{A^{{m_c}}}     , </math>
 |}
 | style="width: 5px;text-align: right;white-space: nowrap;" | (5)
@@ Line 125: / Line 119: @@
 {| style="text-align: left; margin:auto;width: 100%;"
 |-
-| style="text-align: center;" | <math>\frac{{\partial h}}{{\partial t}} + \frac{{\partial q}}{{\partial x}} = {i_e}</math>
+| style="text-align: center;" | <math>\frac{{\partial h}}{{\partial t}} + \frac{{\partial q}}{{\partial x}} = {i_e}\quad ; </math> <math> \quad   q = {\beta _s}{h^{{m_s}}}, </math>
-|-
-| style="text-align: center;" | <math>   q = {\beta _s}{h^{{m_s}}}     , </math>
 |}
 | style="width: 5px;text-align: right;white-space: nowrap;" | (6)
@@ Line 142: / Line 134: @@
 |-
 | style="text-align: center;" | <math>h(0,t) = 0, \quad t > 0</math>
+|}
 | style="width: 5px;text-align: right;white-space: nowrap;" | (7)
+|}
+{| class="formulaSCP" style="width: 100%; text-align: left;"
+|-
+|
+{| style="text-align: left; margin:auto;width: 100%;"
 |-
 | style="text-align: center;" | <math>  h(x,0) = 0, \quad 0 \le x \le L </math>
@@ Line 153: / Line 152: @@
 A aplicação do modelo onda cinemática em escoamento superficial difere da sua aplicação em rios, apenas pelo meio onde ocorre o escoamento. Quanto a rugosidade, utiliza-se a resistência hidráulica de Manning, assim <math display="inline">\beta _c = \beta _s = \sqrt{S_0} /n_r</math> e <math display="inline">m_c = m_s = 5/3</math>, onde <math display="inline">n_r</math> é o coeficiente de rugosidade de Manning.
-==3 Método do Reticulado de Boltzmann==
+==3. Método do reticulado de Boltzmann==
-De acordo com Chen and Doolen (1998) <span id='citeF-3'></span>[[#cite-3|[3]]], o LBM é um esquema numérico baseado em equações cinéticas formuladas em uma escala mesoscópica que simula a dinâmica do fluido em uma escala macroscópica. A dinâmica deste método é dada pela equação do reticulado de Boltzmann, a qual mostra como ocorrem as distribuições das micropartículas.
+De acordo com Chen and Doolen <span id='citeF-3'></span>[[#cite-3|[3]]], o LBM é um esquema numérico baseado em equações cinéticas formuladas em uma escala mesoscópica que simula a dinâmica do fluido em uma escala macroscópica. A dinâmica deste método é dada pela equação do reticulado de Boltzmann, a qual mostra como ocorrem as distribuições das micropartículas.
 O LBM considera uma dinâmica molecular de partículas fictícias em que o espaço, o tempo e as velocidades são discretas. Essas partículas viajam de um ponto ao outro na malha em tempos discretos, se encontrando nestes pontos ao final de cada passo de tempo e trocam quantidade de movimento e energia. No que se refere aos módulos de velocidades, eles assumem valores contínuos. Por outro lado, o LBM pode ser visto como uma forma simplificada da equação cinética de Boltzmann, em que somente são mantidos os detalhes moleculares essenciais para se recuperar o comportamento macroscópico correto <span id='citeF-35'></span>[[#cite-35|[35]]].
@@ Line 161: / Line 160: @@
 Quatro elementos principais caracterizam o LBM: a equação governante, o operador de colisão, o reticulado e a função distribuição de equilíbrio <span id='citeF-2'></span>[[#cite-2|[2]]].
-===3.1 Equação do LBM e o Operador BGK===
+===3.1 Equação do LBM e o operador BGK===
-Inicialmente apresenta-se em [[#eq-9|(9)]], a equação governante do LBM com aperador de colisão BGK, denominada equação do reticulado de Boltzmann. A forma do operador de colisão foi determinada por Bhatnagar, Gross e Krook <span id='citeF-36'></span>[[#cite-36|[36]]], com o objetivo de simplificar a equação cinética de Boltzmann.  A equação do reticulado de Boltzmann, possui as etapas de propagação e colisão. Na etapa de propagação, as partículas se movem de um nó da malha para um dos seus nós vizinhos, com direção dada pela velocidade. Na etapa de colisão, as partículas que chegam no mesmo nó interagem entre si e mudam suas direções conforme as diretrizes do operador de colisão <span id='citeF-37'></span>[[#cite-37|[37]]].
+Inicialmente apresenta-se em equação [[#eq-9|(9)]], a equação governante do LBM com aperador de colisão BGK, denominada equação do reticulado de Boltzmann. A forma do operador de colisão foi determinada por Bhatnagar et al. <span id='citeF-36'></span>[[#cite-36|[36]]], com o objetivo de simplificar a equação cinética de Boltzmann.  A equação do reticulado de Boltzmann, possui as etapas de propagação e colisão. Na etapa de propagação, as partículas se movem de um nó da malha para um dos seus nós vizinhos, com direção dada pela velocidade. Na etapa de colisão, as partículas que chegam no mesmo nó interagem entre si e mudam suas direções conforme as diretrizes do operador de colisão <span id='citeF-37'></span>[[#cite-37|[37]]]
 <span id="eq-9"></span>
 {| class="formulaSCP" style="width: 100%; text-align: left;"
@@ Line 193: / Line 191: @@
 No LBM o reticulado tem a função de representar os pontos da malha e determinar as direções de movimento das partículas. As direções finitas e determinadas para o movimento das partículas definem num modelo microscópico para a dinâmica molecular. Na escolha do reticulado a ser utilizado é essencial observar a sua simetria. Ela é necessária para que o reticulado possa representar as equações macroscópicas <span id='citeF-2'></span>[[#cite-2|[2]]].
-Encontra-se em Qian ''et al.'' <span id='citeF-30'></span>[[#cite-30|[30]]], uma família de reticulados denominada <math display="inline">DkQm</math>, onde <math display="inline">k</math> indica o espaço <math display="inline">k</math>-dimensional e <math display="inline">m</math> as direções de movimento para as distribuições de partículas. O modelo de reticulado D1Q5 é unidimensional com quatro direções de velocidades não nulas e uma para velocidade nula. As partículas do fluido se movimentam para dois nós vizinhos a esquerda e dois nós vizinhos a direita. A Figura [[#img-1|1]] mostra o reticulado D1Q5 com velocidades <math display="inline">\vec e_0=0</math>, <math display="inline">\vec e_1=-e</math>, <math display="inline">\vec e_2=-2e</math>, <math display="inline">\vec e_3=e</math> e <math display="inline">\vec e_4=2e</math>, sendo <math display="inline">e=\Delta x / \Delta t</math> a velocidade no reticulado, em que <math display="inline">\Delta x</math> é o espaçamento da malha.
+Encontra-se em Qian ''et al.'' <span id='citeF-30'></span>[[#cite-30|[30]]], uma família de reticulados denominada <math display="inline">DkQm</math>, onde <math display="inline">k</math> indica o espaço <math display="inline">k</math>-dimensional e <math display="inline">m</math> as direções de movimento para as distribuições de partículas. O modelo de reticulado D1Q5 é unidimensional com quatro direções de velocidades não nulas e uma para velocidade nula. As partículas do fluido se movimentam para dois nós vizinhos a esquerda e dois nós vizinhos a direita. A  [[#img-1|Figura 1]] mostra o reticulado D1Q5 com velocidades <math display="inline">\vec e_0=0</math>, <math display="inline">\vec e_1=-e</math>, <math display="inline">\vec e_2=-2e</math>, <math display="inline">\vec e_3=e</math> e <math display="inline">\vec e_4=2e</math>, sendo <math display="inline">e=\Delta x / \Delta t</math> a velocidade no reticulado, em que <math display="inline">\Delta x</math> é o espaçamento da malha.
 <div id='img-1'></div>
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
 |-
-|[[Image:Draft_Galina_237452870-Reticulado_D1Q5.png|360px|Reticulado D1Q5]]
+|style="padding:10px;"| [[Image:Draft_Galina_237452870-Reticulado_D1Q5.png|360px|Reticulado D1Q5]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figura 1:''' Reticulado D1Q5
+| colspan="1" | '''Figura 1'''. Reticulado D1Q5
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Fonte:''' Zhang e Yan (2008) <span id='citeF-38'></span>[[#cite-38|[38]]].
+| colspan="1" style="padding-bottom:10px;"| '''Fonte:''' Zhang e Yan <span id='citeF-38'></span>[[#cite-38|[38]]]
 |}
 ===3.3 Expansão multi-escala===
@@ Line 219: / Line 215: @@
 {| style="text-align: left; margin:auto;width: 100%;"
 |-
-| style="text-align: center;" | <math>\dfrac{\partial }{{\partial t}} = \dfrac{\partial }{{\partial {t_0}}} + \varepsilon \dfrac{\partial }{{\partial {t_1}}} + {\varepsilon ^2}\dfrac{\partial }{{\partial {t_2}}} + {\varepsilon ^3}\dfrac{\partial }{{\partial {t_3}}} + {\varepsilon ^4}\dfrac{\partial }{{\partial {t_4}}} + O({\varepsilon ^5}),</math>
+| style="text-align: center;" | <math>\begin{align}& \dfrac{\partial }{{\partial t}} = \dfrac{\partial }{{\partial {t_0}}} + \varepsilon \dfrac{\partial }{{\partial {t_1}}} + {\varepsilon ^2}\dfrac{\partial }{{\partial {t_2}}} + {\varepsilon ^3}\dfrac{\partial }{{\partial {t_3}}} + {\varepsilon ^4}\dfrac{\partial }{{\partial {t_4}}} + O({\varepsilon ^5}),\\
-|-
+& \dfrac{\partial }{{\partial x}} = \dfrac{\partial }{{\partial x}} + O(\varepsilon ),\end{align}     </math>
-| style="text-align: center;" | <math>   \dfrac{\partial }{{\partial x}} = \dfrac{\partial }{{\partial x}} + O(\varepsilon ),     </math>
 |}
 | style="width: 5px;text-align: right;white-space: nowrap;" | (11)
@@ Line 273: / Line 268: @@
 {| style="text-align: left; margin:auto;width: 100%;"
 |-
-| style="text-align: center;" | <math>\varepsilon {\Delta }{f_\alpha } + \frac{{{\varepsilon ^2}}}{2}{\Delta ^2}{f_\alpha } + \frac{{{\varepsilon ^3}}}{6}{\Delta ^3}{f_\alpha } </math>
+| style="text-align: center;" | <math>\varepsilon {\Delta }{f_\alpha } + \frac{{{\varepsilon ^2}}}{2}{\Delta ^2}{f_\alpha } + \frac{{{\varepsilon ^3}}}{6}{\Delta ^3}{f_\alpha } </math><math>+ \frac{{{\varepsilon ^4}}}{24}{\Delta ^4}{f_\alpha } = - \frac{1}{\tau }\left({{f_\alpha } - f_\alpha ^{eq}} \right)+ {\varepsilon ^2}{g_\alpha }, </math>
-|-
-| style="text-align: center;" | <math>+ \frac{{{\varepsilon ^4}}}{24}{\Delta ^4}{f_\alpha } = - \frac{1}{\tau }\left({{f_\alpha } - f_\alpha ^{eq}} \right)+ {\varepsilon ^2}{g_\alpha }  , </math>
 |}
 | style="width: 5px;text-align: right;white-space: nowrap;" | (15)
@@ Line 305: / Line 298: @@
 |-
 | style="text-align: center;" | <math>O({\varepsilon ^0}) : f_\alpha ^{(0)} = f_\alpha ^{eq}   ,</math>
+|}
 | style="width: 5px;text-align: right;white-space: nowrap;" | (17)
+|}
+{| class="formulaSCP" style="width: 100%; text-align: left;"
+|-
+|
+{| style="text-align: left; margin:auto;width: 100%;"
 |-
 | style="text-align: center;" | <math>  O({\varepsilon ^1}) : \Delta \,f_\alpha ^{(0)} =  - \frac{1}{\tau }f_\alpha ^{(1)}   ,</math>
+|}
 | style="width: 5px;text-align: right;white-space: nowrap;" | (18)
+|}
+{| class="formulaSCP" style="width: 100%; text-align: left;"
+|-
+|
+{| style="text-align: left; margin:auto;width: 100%;"
 |-
 | style="text-align: center;" | <math>  O({\varepsilon ^2}) : \frac{{\partial f_\alpha ^{(0)}}}{{\partial {t_1}}} + {p}_{2}{\Delta ^2}f_\alpha ^{(0)} =  - \frac{1}{\tau }f_\alpha ^{(2)} + {g_\alpha }  ,</math>
+|}
 | style="width: 5px;text-align: right;white-space: nowrap;" | (19)
+|}
+{| class="formulaSCP" style="width: 100%; text-align: left;"
 |-
-| style="text-align: center;" | <math>  O({\varepsilon ^3}) : \frac{{\partial f_\alpha ^{(0)}}}{{\partial {t_2}}} + p_{3} {\Delta ^3}f_\alpha ^{(0)} + 2p_2 \Delta \frac{{\partial f_\alpha ^{(0)}}}{{\partial {t_1}}} =  - \frac{1}{\tau }f_\alpha ^{(3)} </math>
+|
+{| style="text-align: left; margin:auto;width: 100%;"
 |-
-| style="text-align: center;" | <math>-\tau \Delta {g_\alpha }  ,  </math>
+| style="text-align: center;" | <math>  O({\varepsilon ^3}) : \frac{{\partial f_\alpha ^{(0)}}}{{\partial {t_2}}} + p_{3} {\Delta ^3}f_\alpha ^{(0)} + 2p_2 \Delta \frac{{\partial f_\alpha ^{(0)}}}{{\partial {t_1}}} =  - \frac{1}{\tau }f_\alpha ^{(3)} -\tau \Delta {g_\alpha },  </math>
+|}
 | style="width: 5px;text-align: right;white-space: nowrap;" | (20)
+|}
+{| class="formulaSCP" style="width: 100%; text-align: left;"
 |-
-| style="text-align: center;" | <math>  O({\varepsilon ^4}) : \frac{{\partial f_\alpha ^{(0)}}}{{\partial {t_3}}} + p_{4} {\Delta ^4}f_\alpha ^{(0)} + 3p_{3} {\Delta ^2}\frac{{\partial f_\alpha ^{(0)}}}{{\partial {t_1}}} </math>
+|
+{| style="text-align: left; margin:auto;width: 100%;"
 |-
-| style="text-align: center;" | <math>+ 2p_2 \Delta \frac{{\partial f_\alpha ^{(0)}}}{{\partial {t_2}}} + p_{2} \frac{{\partial ^2 f_\alpha ^{(0)}}}{{\partial t_1^{2}}} =  - \frac{1}{\tau }f_\alpha ^{(4)} </math>
+| style="text-align: center;" | <math>  O({\varepsilon ^4}): \frac{{\partial f_\alpha ^{(0)}}}{{\partial {t_3}}} + p_{4} {\Delta ^4}f_\alpha ^{(0)} + 3p_{3} {\Delta ^2}\frac{{\partial f_\alpha ^{(0)}}}{{\partial {t_1}}} + 2p_2 \Delta \frac{{\partial f_\alpha ^{(0)}}}{{\partial {t_2}}} </math>
 |-
-| style="text-align: center;" | <math> + p_2 \tau \Delta ^2 {g_\alpha } - \tau \frac{{\partial g_\alpha }}{{\partial {t_1}}},    </math>
+| style="text-align: center;" | <math> + p_{2} \frac{{\partial ^2 f_\alpha ^{(0)}}}{{\partial t_1^{2}}} =  - \frac{1}{\tau }f_\alpha ^{(4)}  + p_2 \tau \Delta ^2 {g_\alpha } - \tau \frac{{\partial g_\alpha }}{{\partial {t_1}}},    </math>
-| style="width: 5px;text-align: right;white-space: nowrap;" | (21)
 |}
+| style="width: 5px;text-align: right;white-space: nowrap;" | (21)
 |}
@@ Line 335: / Line 352: @@
 {| style="text-align: left; margin:auto;width: 100%;"
 |-
-| style="text-align: center;" | <math>{p_1}(\tau ) = \tau , </math>
+| style="text-align: center;" | <math>\begin{align} & {p_1}(\tau ) = \tau , \\
-|-
+& {p_2}(\tau ) = \frac{1}{2} - \tau , \\
-| style="text-align: center;" | <math>   {p_2}(\tau ) = \frac{1}{2} - \tau , </math>
+&   {p_3}(\tau ) = {\tau ^2} - \tau  + \frac{1}{6}, \\
-|-
+& {p_4}(\tau ) =  - {\tau ^3} + \frac{3}{2}{\tau ^2} - \frac{7}{{12}}\tau  + \frac{1}{{24}}.\end{align}       </math>
-| style="text-align: center;" | <math>   {p_3}(\tau ) = {\tau ^2} - \tau  + \frac{1}{6}, </math>
-|-
-| style="text-align: center;" | <math>   {p_4}(\tau ) =  - {\tau ^3} + \frac{3}{2}{\tau ^2} - \frac{7}{{12}}\tau  + \frac{1}{{24}}.       </math>
 |}
 | style="width: 5px;text-align: right;white-space: nowrap;" | (22)
@@ Line 348: / Line 362: @@
 Por meio das equações [[#eq-17|(17)]], [[#eq-18|(18)]], [[#eq-19|(19)]], [[#eq-20|(20)]] e [[#eq-21|(21)]] determinam-se os momentos da função distribuição de equilíbrio, e a partir das equações geradas da expansão dos momentos sobre o reticulado, determinam-se as expressões da função distribuição de equilíbrio nas direções do reticulado.
-===3.4 Função Distribuição de Equilíbrio===
+===3.4 Função distribuição de equilíbrio===
 A função distribuição de equilíbrio juntamente com a equação do Reticulado de Boltzmann [[#eq-9|(9)]] desempenham um papel essencial que é recuperar a equação macroscópica do fluido. A seguir, recupera-se a equação da onda cinemática unidimensional [[#eq-5|(5)]] e determina-se a sua função distribuição de equilíbrio.
@@ Line 354: / Line 368: @@
 A equação da onda cinemática é um caso especial da equação de Burgers. A equação de Burgers é uma equação de convecção-difusão não-linear que representa um modelo simplificado das equações de Navier-Stokes. É utilizada para testar a eficiência de vários esquemas numéricos e para descrever o escoamento por meio de uma onda de choque que viaja em um fluido viscoso <span id='citeF-40'></span>[[#cite-40|[40]]].
-A equação de Burgers unidimensional pode ser escrita da seguinte forma <span id='citeF-40'></span><span id='citeF-22'></span>[[#cite-40|[40,22]]],
+A equação de Burgers unidimensional pode ser escrita da seguinte forma <span id='citeF-22'></span><span id='citeF-40'></span>[[#cite-22|[22,40]]],
 <span id="eq-23"></span>
@@ Line 362: / Line 376: @@
 {| style="text-align: left; margin:auto;width: 100%;"
 |-
-| style="text-align: center;" | <math>\frac{{\partial v}}{{\partial t}} + \frac{\partial }{{\partial x}}\left({\beta {v^m}} \right)= \gamma \frac{{{\partial ^2}v}}{{\partial x^2}} + F  , </math>
+| style="text-align: center;" | <math>\frac{{\partial v}}{{\partial t}} + \frac{\partial }{{\partial x}}\left({\beta {v^m}} \right)= \gamma \frac{{{\partial ^2}v}}{{\partial x^2}} + F , </math>
 |}
 | style="width: 5px;text-align: right;white-space: nowrap;" | (23)
@@ Line 383: / Line 397: @@
 em consequência, as somas sobre a função distribuição de não equilíbrio têm as seguintes restrições,
 <span id="eq-25"></span>
 {| class="formulaSCP" style="width: 100%; text-align: left;"
@@ Line 395: / Line 408: @@
 |}
-Combinando as equações [[#eq-17|(17)]] a [[#eq-21|(21)]], obedecendo as restrições em [[#eq-24|(24)]] e [[#eq-25|(25)]], em vista das diferentes escalas de tempo e da equação de Burgers, equação [[#eq-23|(23)]], determinam-se os momentos da função distribuição de equilíbrio de zero a quarta ordem, descritos na equação [[#eq-26|(26)]]. Neste artigo, assume-se que a força externa não varia no espaço e no tempo, portanto são nulas as derivadas nos termos que envolvem <math display="inline">g_\alpha </math> do lado direito das equações [[#eq-20|(20)]] e [[#eq-21|(21)]].
+Combinando as equações [[#eq-17|(17)]] a [[#eq-21|(21)]], obedecendo as restrições em [[#eq-24|(24)]] e [[#eq-25|(25)]], em vista das diferentes escalas de tempo e da equação de Burgers, equação [[#eq-23|(23)]], determinam-se os momentos da função distribuição de equilíbrio de zero a quarta ordem, descritos na equação [[#eq-26|(26)]]. Neste artigo, assume-se que a força externa não varia no espaço e no tempo, portanto são nulas as derivadas nos termos que envolvem <math display="inline">g_\alpha </math> do lado direito das equações [[#eq-20|(20)]] e [[#eq-21|(21)]]
 <span id="eq-26"></span>
@@ Line 403: / Line 416: @@
 {| style="text-align: left; margin:auto;width: 100%;"
 |-
-| style="text-align: center;" | <math>\begin{array}{lcl}\sum \nolimits _\alpha  {f_\alpha ^{(0)}}  &=& v, \\   \sum \nolimits _\alpha  {{{\vec e}_\alpha }f_\alpha ^{(0)}}  &=& \beta {v^m }, \\	   \sum \nolimits _\alpha  {\vec e_\alpha ^{\,2} f_\alpha ^{(0)}}  &=& \dfrac{{{\beta ^2}{m^2}}}{{2m - 1}}{v^{2m - 1}}, \\	   \sum \nolimits _\alpha  {\vec e_\alpha ^{\,3} f_\alpha ^{(0)}} &=& \dfrac{{{\beta ^3}{m^3}}}{{3m - 2}}{v^{3m - 2}}, \\	   \sum \nolimits _\alpha  {\vec e_\alpha ^{\,4} f_\alpha ^{(0)}} &=& \dfrac{{{\beta ^4}{m^4}}}{{4m - 3}}{v^{4m - 3}}.  \end{array}   </math>
+| style="text-align: center;" | <math>\begin{align} & \sum \nolimits _\alpha  {f_\alpha ^{(0)}}   = v, \\    & \sum \nolimits _\alpha  {{{\vec e}_\alpha }f_\alpha ^{(0)}}  = \beta {v^m }, \\
+&  \sum \nolimits _\alpha  {\vec e_\alpha ^{\,2} f_\alpha ^{(0)}}   = \dfrac{{{\beta ^2}{m^2}}}{{2m - 1}}{v^{2m - 1}}, \\
+&  \sum \nolimits _\alpha  {\vec e_\alpha ^{\,3} f_\alpha ^{(0)}}  = \dfrac{{{\beta ^3}{m^3}}}{{3m - 2}}{v^{3m - 2}}, \\
+&  \sum \nolimits _\alpha  {\vec e_\alpha ^{\,4} f_\alpha ^{(0)}}  = \dfrac{{{\beta ^4}{m^4}}}{{4m - 3}}{v^{4m - 3}}.  \end{align}   </math>
 |}
 | style="width: 5px;text-align: right;white-space: nowrap;" | (26)
@@ Line 409: / Line 425: @@
 Desenvolvendo as equações em [[#eq-26|(26)]] nas direções das velocidades do reticulado D1Q5, obtém-se,
 <span id="eq-27"></span>
 {| class="formulaSCP" style="width: 100%; text-align: left;"
@@ Line 429: / Line 444: @@
 {| style="text-align: left; margin:auto;width: 100%;"
 |-
-| style="text-align: center;" | <math>f_1^{eq} = \frac{1}{6}\left({ - 4{C_1} + 4{C_2} + {C_3} - {C_4}} \right), </math>
+| style="text-align: center;" | <math>\begin{align} & f_1^{eq} = \frac{1}{6}\left({ - 4{C_1} + 4{C_2} + {C_3} - {C_4}} \right), \\
-|-
+& f_2^{eq} = \frac{1}{{24}}\left({2{C_1} - {C_2} - 2{C_3} + {C_4}} \right), \\
-| style="text-align: center;" | <math>   f_2^{eq} = \frac{1}{{24}}\left({2{C_1} - {C_2} - 2{C_3} + {C_4}} \right), </math>
+&  f_3^{eq} = \frac{1}{6}\left({4{C_1} + 4{C_2} - {C_3} - {C_4}} \right), \\
-|-
+& f_4^{eq} = \frac{1}{{24}}\left({ - 2{C_1} - {C_2} + 2{C_3} + {C_4}} \right), \\
-| style="text-align: center;" | <math>   f_3^{eq} = \frac{1}{6}\left({4{C_1} + 4{C_2} - {C_3} - {C_4}} \right), </math>
+&  f_0^{eq} = v - \left({f_1^{eq} + f_2^{eq} + f_3^{eq} + f_4^{eq}} \right), \end{align}    </math>
-|-
-| style="text-align: center;" | <math>   f_4^{eq} = \frac{1}{{24}}\left({ - 2{C_1} - {C_2} + 2{C_3} + {C_4}} \right), </math>
-|-
-| style="text-align: center;" | <math>   f_0^{eq} = v - \left({f_1^{eq} + f_2^{eq} + f_3^{eq} + f_4^{eq}} \right),     </math>
 |}
 | style="width: 5px;text-align: right;white-space: nowrap;" | (28)
@@ Line 469: / Line 480: @@
 em que <math display="inline">F</math> pode representar a precipitação quando se trata do escoamento na superfície da bacia, ou representar a contribuição lateral quando se trata do escoamento em canais.
-===3.5 Condições de Contorno===
+===3.5 Condições de contorno===
-A representação adequada das características físicas do problema, por meio das condições iniciais e de contorno, é um fator crucial para a estabilidade e precisão das simulações. A Figura [[#img-2|2]] ilustra uma malha unidimensional, com as direções de velocidade do reticulado D1Q5, para discretização do domínio que representa o escoamento na superfície e nos canais da bacia hidrográfica.
+A representação adequada das características físicas do problema, por meio das condições iniciais e de contorno, é um fator crucial para a estabilidade e precisão das simulações. A  [[#img-2|Figura 2]] ilustra uma malha unidimensional, com as direções de velocidade do reticulado D1Q5, para discretização do domínio que representa o escoamento na superfície e nos canais da bacia hidrográfica.
 <div id='img-2'></div>
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
 |-
-|[[Image:Draft_Galina_237452870-Malha_D1Q5_fd.png|420px|Malha computacional]]
+|style="padding:10px;"| [[Image:Draft_Galina_237452870-Malha_D1Q5_fd.png|420px|Malha computacional]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figura 2:''' Malha computacional
+| colspan="1" style="padding:10px;"| '''Figura 2'''. Malha computacional
 |}
@@ Line 490: / Line 501: @@
 {| style="text-align: left; margin:auto;width: 100%;"
 |-
-| style="text-align: center;" | <math>\bar{\bar{f_0}}(x_0) = f_0^{eq}(x_0), </math>
+| style="text-align: center;" | <math>\begin{align} & \bar{\bar{f_0}}(x_0) = f_0^{eq}(x_0), \\
-|-
+& \bar{\bar{f_1}}(x_0) = f_1^{eq}(x_0), \\
-| style="text-align: center;" | <math>   \bar{\bar{f_1}}(x_0) = f_1^{eq}(x_0), </math>
+&  \bar{\bar{f_2}}(x_0) = f_2^{eq}(x_0),\\
-|-
+& \bar{\bar{f_3}}(x_0) = f_3^{eq}(x_0), \\
-| style="text-align: center;" | <math>   \bar{\bar{f_2}}(x_0) = f_2^{eq}(x_0), </math>
+&   \bar{\bar{f_4}}(x_0) = {v^*} -  (\bar{\bar{f_0}}(x_0)+\bar{\bar{f_1}}(x_0)+\bar{\bar{f_2}}(x_0)+\bar{\bar{f_3}}(x_0)), \end{align}    </math>
-|-
-| style="text-align: center;" | <math>   \bar{\bar{f_3}}(x_0) = f_3^{eq}(x_0), </math>
-|-
-| style="text-align: center;" | <math>   \bar{\bar{f_4}}(x_0) = {v^*} </math>
-|-
-| style="text-align: center;" | <math>-  (\bar{\bar{f_0}}(x_0)+\bar{\bar{f_1}}(x_0)+\bar{\bar{f_2}}(x_0)+\bar{\bar{f_3}}(x_0)),     </math>
 |}
 |}
@@ Line 510: / Line 515: @@
 {| style="text-align: left; margin:auto;width: 100%;"
 |-
-| style="text-align: center;" | <math>\bar{\bar{f_0}}(x_1) = \bar{f_0}(x_1), </math>
+| style="text-align: center;" | <math>\begin{align} &\bar{\bar{f_0}}(x_1) = \bar{f_0}(x_1), \\
-|-
+&   \bar{\bar{f_1}}(x_1) = \bar{f_1}(x_2), \\
-| style="text-align: center;" | <math>    \bar{\bar{f_1}}(x_1) = \bar{f_1}(x_2), </math>
+& \bar{\bar{f_2}}(x_1) = \bar{f_2}(x_3), \\
-|-
+&  \bar{\bar{f_3}}(x_1) = \bar{f_3}(x_0), \\
-| style="text-align: center;" | <math>    \bar{\bar{f_2}}(x_1) = \bar{f_2}(x_3), </math>
+&  \bar{\bar{f_4}}(x_1) = 2\bar{\bar{f_4}}(x_2)-\bar{\bar{f_4}}(x_3).        \end{align}    </math>
-|-
-| style="text-align: center;" | <math>    \bar{\bar{f_3}}(x_1) = \bar{f_3}(x_0), </math>
-|-
-| style="text-align: center;" | <math>    \bar{\bar{f_4}}(x_1) = 2\bar{\bar{f_4}}(x_2)-\bar{\bar{f_4}}(x_3).          </math>
 |}
 | style="width: 5px;text-align: right;white-space: nowrap;" | (31)
@@ Line 530: / Line 531: @@
 {| style="text-align: left; margin:auto;width: 100%;"
 |-
-| style="text-align: center;" | <math>\bar{\bar{f_0}}(x_{n-1}) = \bar{f_0}(x_{n-1}), </math>
+| style="text-align: center;" | <math>\begin{align} &\bar{\bar{f_0}}(x_{n-1}) = \bar{f_0}(x_{n-1}), \\
-|-
+&\bar{\bar{f_3}}(x_{n-1}) = \bar{f_3}(x_{n-2}), \\
-| style="text-align: center;" | <math>   \bar{\bar{f_3}}(x_{n-1}) = \bar{f_3}(x_{n-2}), </math>
+& \bar{\bar{f_4}}(x_{n-1}) = \bar{f_4}(x_{n-3}), \\
-|-
+&\bar{\bar{f_1}}(x_{n-1}) = \bar{f_1}(x_n),\\
-| style="text-align: center;" | <math>   \bar{\bar{f_4}}(x_{n-1}) = \bar{f_4}(x_{n-3}), </math>
+& \bar{\bar{f_2}}(x_{n-1}) = 2\bar{\bar{f_2}}(x_{n-2})-\bar{\bar{f_2}}(x_{n-3}),   \end{align}  </math>
-|-
-| style="text-align: center;" | <math>   \bar{\bar{f_1}}(x_{n-1}) = \bar{f_1}(x_n), </math>
-|-
-| style="text-align: center;" | <math>   \bar{\bar{f_2}}(x_{n-1}) = 2\bar{\bar{f_2}}(x_{n-2})-\bar{\bar{f_2}}(x_{n-3}),    </math>
 |}
 |}
@@ Line 547: / Line 544: @@
 {| style="text-align: left; margin:auto;width: 100%;"
 |-
-| style="text-align: center;" | <math>\bar{\bar{f_0}}(x_n) = \bar{f_0}(x_n), </math>
+| style="text-align: center;" | <math>\begin{align} &\bar{\bar{f_0}}(x_n) = \bar{f_0}(x_n), \\
-|-
+&  \bar{\bar{f_3}}(x_n) = \bar{f_3}(x_{n-1}), \\
-| style="text-align: center;" | <math>    \bar{\bar{f_3}}(x_n) = \bar{f_3}(x_{n-1}), </math>
+&  \bar{\bar{f_4}}(x_n) = \bar{f_4}(x_{n-2}), \\
-|-
+& \bar{\bar{f_1}}(x_n) = 2\bar{\bar{f_1}}(x_{n-1})-\bar{\bar{f_1}}(x_{n-2}), \\
-| style="text-align: center;" | <math>    \bar{\bar{f_4}}(x_n) = \bar{f_4}(x_{n-2}), </math>
+& \bar{\bar{f_2}}(x_n) = 2\bar{\bar{f_2}}(x_{n-1})-\bar{\bar{f_2}}(x_{n-2}),  \end{align}    </math>
-|-
-| style="text-align: center;" | <math>     \bar{\bar{f_1}}(x_n) = 2\bar{\bar{f_1}}(x_{n-1})-\bar{\bar{f_1}}(x_{n-2}), </math>
-|-
-| style="text-align: center;" | <math>     \bar{\bar{f_2}}(x_n) = 2\bar{\bar{f_2}}(x_{n-1})-\bar{\bar{f_2}}(x_{n-2}),      </math>
 |}
 | style="width: 5px;text-align: right;white-space: nowrap;" | (32)
@@ Line 562: / Line 555: @@
 em que <math display="inline">\bar{\bar{f}}</math> é a função distribuição de velocidade das partículas depois da propagação, <math display="inline">\bar{f}</math> é a função distribuição de velocidade das partículas antes da propagação, <math display="inline">v^*</math> é a variável macroscópica na fronteira, podendo ser substituída pela profundidade do escoamento <math display="inline">h</math> ou pela área da seção transversal <math display="inline">A</math>.
-No modelo computacional, as características físicas em cada trecho de canal são mantidas constantes durante a simulação, isto é, não são alterados os valores do comprimento, largura, declividade e coeficiente de Manning do segmento de canal. O que pode implicar em seções transversais diferentes na junção dos canais, levando a necessidade de implementação de condição de contorno. Para que haja a conservação de massa, a vazão de um segmento de canal a montante para um segmento a jusante deve permanecer a mesma. Considerando a Figura [[#img-3|3]], e o cálculo da vazão em uma seção transversal dada pela fórmula de Manning, a vazão no início do canal 2 é determinada por,
+No modelo computacional, as características físicas em cada trecho de canal são mantidas constantes durante a simulação, isto é, não são alterados os valores do comprimento, largura, declividade e coeficiente de Manning do segmento de canal. O que pode implicar em seções transversais diferentes na junção dos canais, levando a necessidade de implementação de condição de contorno. Para que haja a conservação de massa, a vazão de um segmento de canal a montante para um segmento a jusante deve permanecer a mesma. Considerando a  [[#img-3|Figura 3]], e o cálculo da vazão em uma seção transversal dada pela fórmula de Manning, a vazão no início do canal 2 é determinada por,
 <span id="eq-33"></span>
@@ Line 578: / Line 571: @@
 <div id='img-3'></div>
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
 |-
-|[[Image:Draft_Galina_237452870-juncao_canais.png|409px|Junção de canais das sub-bacias]]
+|style="padding:10px;"|[[Image:Draft_Galina_237452870-juncao_canais.png|309px|Junção de canais das sub-bacias]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figura 3:''' Junção de canais das sub-bacias
+| colspan="1" style="padding:10px;" | '''Figura 3'''. Junção de canais das sub-bacias
 |}
@@ Line 602: / Line 595: @@
 ===3.6 Estabilidade===
-De acordo com Sterling e Chen (1996) <span id='citeF-41'></span>[[#cite-41|[41]]], não é possível garantir a estabilidade do LBM e isto ocorre devido a grande quantidade de parâmetros que impedem a sua caracterização completa. No entanto, pode-se seguir algumas condições necessárias, que, se forem satisfeitas simultaneamente, obtém-se a estabilidade da simulação <span id='citeF-2'></span>[[#cite-2|[2]]].  A primeira expressão de [[#eq-35|(35)]] relaciona o parâmetro de relaxamento <math display="inline">\tau </math> e a viscosidade cinemática <math display="inline">\nu </math>, dessa relação resulta <math display="inline">\tau > 1/2</math> <span id='citeF-41'></span>[[#cite-41|[41]]]. A magnitude da velocidade física <math display="inline">\vec v</math> resultante do fluido deve ser menor do que a velocidade no reticulado, conforme a segunda expressão. Na terceira expressão tem-se a condição para celeridade e na quarta a restrição para escoamentos subcríticos <span id='citeF-2'></span>[[#cite-2|[2]]].
+De acordo com Sterling e Chen (1996) <span id='citeF-41'></span>[[#cite-41|[41]]], não é possível garantir a estabilidade do LBM e isto ocorre devido a grande quantidade de parâmetros que impedem a sua caracterização completa. No entanto, pode-se seguir algumas condições necessárias, que, se forem satisfeitas simultaneamente, obtém-se a estabilidade da simulação <span id='citeF-2'></span>[[#cite-2|[2]]].  A primeira expressão de [[#eq-35|(35)]] relaciona o parâmetro de relaxamento <math display="inline">\tau </math> e a viscosidade cinemática <math display="inline">\nu </math>, dessa relação resulta <math display="inline">\tau > 1/2</math> <span id='citeF-41'></span>[[#cite-41|[41]]]. A magnitude da velocidade física <math display="inline">\vec v</math> resultante do fluido deve ser menor do que a velocidade no reticulado, conforme a segunda expressão. Na terceira expressão tem-se a condição para celeridade e na quarta a restrição para escoamentos subcríticos <span id='citeF-2'></span>[[#cite-2|[2]]]
 <span id="eq-35"></span>
@@ Line 617: / Line 610: @@
 Como o LBM é um método intrinsecamente viscoso <span id='citeF-42'></span>[[#cite-42|[42]]], utiliza-se com frequência, na simulação numérica, uma viscosidade mais elevada do que a viscosidade cinemática do fluido e os resultados estão em concordância com as simulações disponíveis na literatura <span id='citeF-2'></span><span id='citeF-5'></span><span id='citeF-42'></span><span id='citeF-43'></span>[[#cite-2|[2,5,42,43]]].
-==4 Estudo de caso: Escoamento em uma bacia hidrográfica==
+==4. Estudo de caso: Escoamento em uma bacia hidrográfica==
 A bacia hidrográfica consiste em uma área bem definida de captação natural da água de precipitação que faz convergir o escoamento até um curso de água. Na bacia hidrográfica, o escoamento é dividido em superficial e em canal. No escoamento superficial a água desloca-se pela superfície da bacia até encontrar um curso de água. Este deslocamento ocorre como resultado da água precipitada que não foi interceptada pela cobertura vegetal e pela parte que não infiltrou no solo. Devido a grande heterogeneidade espacial, a superfície da bacia é representada por planos onde ocorre o escoamento superficial de pouca profundidade. A contribuição lateral nos canais é dada principalmente pela precipitação que ocorre sobre cada plano <span id='citeF-33'></span>[[#cite-33|[33]]].
-Dependendo do tamanho e das características da bacia hidrográfica, pode-se representá-la por uma forma em “V”, segmentada por dois planos e um segmento de canal, como mostra a Figura [[#img-4|4]]. Neste estudo, utilizam-se várias destas formas para representar as sub-bacias que compõem a superfície da bacia hidrográfica.
+Dependendo do tamanho e das características da bacia hidrográfica, pode-se representá-la por uma forma em “V”, segmentada por dois planos e um segmento de canal, como mostra a  [[#img-4|Figura 4]]. Neste estudo, utilizam-se várias destas formas para representar as sub-bacias que compõem a superfície da bacia hidrográfica.
 <div id='img-4'></div>
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
 |-
-|[[Image:Draft_Galina_237452870-esc_dois_planos.png|390px|Representação de uma sub-bacia]]
+|style="padding:10px;"| [[Image:Draft_Galina_237452870-esc_dois_planos.png|390px|Representação de uma sub-bacia]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figura 4:''' Representação de uma sub-bacia
+| colspan="1" style="padding:10px;" | '''Figura 4'''. Representação de uma sub-bacia
 |}
 ===4.1 Escoamento na bacia hidrográfica do rio Chopim===
-A bacia do rio Iguaçu (Figura [[#img-5|5]]), com superfície aproximada de <math display="inline">70.800</math> km<math display="inline">^2</math>,  tem grande importância na geração de energia por meio de usinas hidrelétricas e pela grande demanda de recursos hídricos no abastecimento público, industrial e agrícola. A região ainda é caracterizada por um grande número de Unidades de Conservação Ambiental, e abriga importantes corredores de biodiversidade <span id='citeF-44'></span>[[#cite-44|[44]]].
+A bacia do rio Iguaçu ([[#img-5|Figura 5]]), com superfície aproximada de <math display="inline">70.800</math> km<math display="inline">^2</math>,  tem grande importância na geração de energia por meio de usinas hidrelétricas e pela grande demanda de recursos hídricos no abastecimento público, industrial e agrícola. A região ainda é caracterizada por um grande número de Unidades de Conservação Ambiental, e abriga importantes corredores de biodiversidade <span id='citeF-44'></span>[[#cite-44|[44]]].
 <div id='img-5'></div>
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
 |-
-|[[Image:Draft_Galina_237452870-localizacao_relevo_bacia_iguacu.png|420px|Localização da bacia hidrográfica do rio Iguaçu]]
+|style="padding:10px;"| [[Image:Draft_Galina_237452870-localizacao_relevo_bacia_iguacu.png|420px|Localização da bacia hidrográfica do rio Iguaçu]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figura 5:''' Localização da bacia hidrográfica do rio Iguaçu
+| colspan="1" | '''Figura 5'''. Localização da bacia hidrográfica do rio Iguaçu
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Fonte:''' Franco (2017) <span id='citeF-45'></span>[[#cite-45|[45]]].
+| colspan="1" style="padding-bottom:10px;" | '''Fonte:''' Franco (2017) <span id='citeF-45'></span>[[#cite-45|[45]]]
 |}
@@ Line 649: / Line 642: @@
 O rio Chopim é o mais importante afluente da bacia do baixo rio Iguaçu. Suas nascentes estão localizadas em altitudes que superam os 1.200 m. A extensão total do curso principal do rio é da ordem de 450 km. Da nascente do rio até em torno do km 209 a inclinação do leito é de aproximadamente 2,9 m/km, e desde este ponto até a foz, a inclinação é de aproximadamente 1,1 m/km <span id='citeF-47'></span>[[#cite-47|[47]]].
-Considera-se nesse estudo, um trecho do rio Chopim localizado entre suas nascentes no município de Palmas, passando pela estação pluviométrica Salto Claudelino até a estação fluviométrica Porto Palmeirinha. Este trecho do rio tem comprimento de aproximadamente <math display="inline">281</math> km com área de drenagem na superfície da bacia em torno de <math display="inline">3.205</math> km<math display="inline">^2</math>. As características das estações estão listadas na Tabela [[#table-B|B]].
+Considera-se nesse estudo, um trecho do rio Chopim localizado entre suas nascentes no município de Palmas, passando pela estação pluviométrica Salto Claudelino até a estação fluviométrica Porto Palmeirinha. Este trecho do rio tem comprimento de aproximadamente <math display="inline">281</math> km com área de drenagem na superfície da bacia em torno de <math display="inline">3.205</math> km<math display="inline">^2</math>. As características das estações estão listadas na  [[#table-A|Tabela A]].
-A Figura [[#img-6|6]] ilustra os principais rios da bacia do rio Iguaçu, suas estações pluviométricas, meteorológicas e fluviométricas, e destaca a localizacão da bacia do rio Chopim.
+A  [[#img-6|Figura 6]] ilustra os principais rios da bacia do rio Iguaçu, suas estações pluviométricas, meteorológicas e fluviométricas, e destaca a localizacão da bacia do rio Chopim.
 <div id='img-6'></div>
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 60%;"
 |-
-|[[Image:Draft_Galina_237452870-bacia_iguaçu_estacoes2.png|450px|Localização das estações pluviométricas, meteorológicas e fluviométricas da bacia hidrográfica do rio Iguaçu]]
+|style="padding:10px;"| [[Image:Draft_Galina_237452870-bacia_iguaçu_estacoes2.png|500px|Localização das estações pluviométricas, meteorológicas e fluviométricas da bacia hidrográfica do rio Iguaçu]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figura 6:''' Localização das estações pluviométricas, meteorológicas e fluviométricas da bacia hidrográfica do rio Iguaçu
+| colspan="1" style="padding-left:10px;padding-right:10px;" | '''Figura 6'''. Localização das estações pluviométricas, meteorológicas e fluviométricas da bacia hidrográfica do rio Iguaçu
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Fonte:''' Adaptado de Uda (2016) <span id='citeF-48'></span>[[#cite-48|[48]]].
+| colspan="1" style="padding-bottom:10px;" | '''Fonte:''' Adaptado de Uda (2016) <span id='citeF-48'></span>[[#cite-48|[48]]]
 |}
-A Figura [[#img-7|7]] mostra a bacia do rio Chopim com as divisões de sub-bacias conforme Kaviski e Gonçalves (1996) <span id='citeF-49'></span>[[#cite-49|[49]]] definidas de acordo com as características físicas da bacia.
+A  [[#img-7|Figura 7]] mostra a bacia do rio Chopim com as divisões de sub-bacias conforme Kaviski e Gonçalves (1996) <span id='citeF-49'></span>[[#cite-49|[49]]] definidas de acordo com as características físicas da bacia.
 <div id='img-7'></div>
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
 |-
-|[[Image:Draft_Galina_237452870-Bacia_chopim_com_subbacias.png|450px|Bacia hidrográfica do rio Chopim]]
+|style="padding:10px;"| [[Image:Draft_Galina_237452870-Bacia_chopim_com_subbacias.png|450px|Bacia hidrográfica do rio Chopim]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figura 7:''' Bacia hidrográfica do rio Chopim
+| colspan="1" | '''Figura 7'''. Bacia hidrográfica do rio Chopim
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Fonte:''' Adaptado de Kaviski e Gonçalves (1996) <span id='citeF-49'></span>[[#cite-49|[49]]].
+| colspan="1" style="padding-bottom:10px;" | '''Fonte:''' Adaptado de Kaviski e Gonçalves (1996) <span id='citeF-49'></span>[[#cite-49|[49]]]
 |}
+==5. Resultados e discussões==
-==5 Resultados e Discussões==
+Levando-se em conta a delimitação das sub-bacias e os trechos de canais apresentados na  [[#img-7|Figura 7]], fez-se a segmentação da bacia do rio Chopim em nove sub-bacias até a estação Porto Palmeirinha. Cada sub-bacia é composta por um trecho de canal e dois planos laterais de escoamento, com exceção da primeira sub-bacia que é representada por apenas um plano lateral. A numeração das sub-bacias foi mantida a mesma que a utilizada na sequência para os trechos de canais. Essa numeração indica o início do escoamento até a estação fluviométrica Porto Palmeirinha. Os dados das características físicas da segmentação da bacia são apresentados na  [[#table-B|Tabela B]] para os segmentos planos, e na  [[#table-C|Tabela C]] para os segmentos de canais.
-Levando-se em conta a delimitação das sub-bacias e os trechos de canais apresentados na Figura [[#img-7|7]], fez-se a segmentação da bacia do rio Chopim em nove sub-bacias até a estação Porto Palmeirinha. Cada sub-bacia é composta por um trecho de canal e dois planos laterais de escoamento, com exceção da primeira sub-bacia que é representada por apenas um plano lateral. A numeração das sub-bacias foi mantida a mesma que a utilizada na sequência para os trechos de canais. Essa numeração indica o início do escoamento até a estação fluviométrica Porto Palmeirinha. Os dados das características físicas da segmentação da bacia são apresentados na Tabela [[#table-C|C]] para os segmentos planos, e na Tabela [[#table-D|D]] para os segmentos de canais.
+Nos trechos de canais não são conhecidos os valores da declividade e do coeficiente de Manning. Para complementar a descrição das características físicas das sub-bacias hidrográficas, listam-se na  [[#table-C|Tabela C]] os parâmetros <math display="inline">\bar{\alpha }</math> e <math display="inline">\bar{\beta }</math>, que descrevem as curvas de descarga. De acordo com Kaviski e Gonçalves (1996) <span id='citeF-49'></span>[[#cite-49|[49]]], a relação entre área da seção transversal e vazão é descrita pela equação,
-Nos trechos de canais não são conhecidos os valores da declividade e do coeficiente de Manning. Para complementar a descrição das características físicas das sub-bacias hidrográficas, listam-se na Tabela [[#table-D|D]] os parâmetros <math display="inline">\bar{\alpha }</math> e <math display="inline">\bar{\beta }</math>, que descrevem as curvas de descarga. De acordo com Kaviski e Gonçalves (1996) <span id='citeF-49'></span>[[#cite-49|[49]]], a relação entre área da seção transversal e vazão é descrita pela equação,
 <span id="eq-36"></span>
@@ Line 716: / Line 708: @@
 {| style="text-align: left; margin:auto;width: 100%;"
 |-
-| style="text-align: center;" | <math>{\beta _c} = \dfrac{1}{{{\bar{\alpha } ^{1/\bar{\beta } }}}},\,\,\,\,\,{m_c} = \dfrac{1}{\bar{\beta } }.   </math>
+| style="text-align: center;" | <math>{\beta _c} = \dfrac{1}{{{\bar{\alpha } ^{1/\bar{\beta } }}}}\qquad,\qquad {m_c} = \dfrac{1}{\bar{\beta } }.   </math>
 |}
 | style="width: 5px;text-align: right;white-space: nowrap;" | (38)
 |}
-Na junção de segmentos de canais (Figura [[#img-3|3]]), a área da seção transversal é encontrada diretamente pela equação [[#eq-36|(36)]].
+Na junção de segmentos de canais ([[#img-3|Figura 3]]), a área da seção transversal é encontrada diretamente pela equação [[#eq-36|(36)]].
 Por não estarem disponíveis os dados de declividade e coeficiente de Manning dos trechos de canais da bacia do rio Chopim, houve a necessidade de adaptar o modelo para utilizar os parâmetros <math display="inline">\bar{\alpha }</math> e <math display="inline">\bar{\beta }</math>. A vantagem em utilizar esses parâmetros é a de aplicar o LBM em um canal com seção transversal estimada pelas características físicas naturais do rio.
@@ Line 727: / Line 719: @@
 Os parâmetros da simulação no LBM para os planos são <math display="inline">\Delta x = 10</math> m, <math display="inline">\Delta t = 10</math> s e velocidade na malha <math display="inline">e=1</math> m/s. Com relação aos canais, <math display="inline">\Delta x = 50</math> m, <math display="inline">\Delta t = 10</math> s e velocidade na malha <math display="inline">e=5</math> m/s. Em ambos, o parâmetro de relaxamento utilizado é <math display="inline">\tau=12</math>. O número de iterações realizadas foi de <math display="inline">90.000</math>, o que corresponde a um tempo de simulação de <math display="inline">250</math> horas ou <math display="inline">10,42</math> dias.
-O período de chuva analisado na bacia é dado no mês 07 de 1987 nos dias 07 e 08.  Neste mês não houve chuva além dos dias citados. A Tabela [[#table-E|E]] contém a intensidade média de precipitação em cada sub-bacia. Para se obter uma simulação mais próxima da realidade, considerou-se uma vazão inicial que ocorre no rio (escoamento de base), em cada trecho de canal, antes de ocorrer a precipitação. Essa vazão entra no modelo computacional como uma condição inicial. Também é considerada como uma condição de contorno ao ser adicionada na contribuição lateral dos planos, em cada passo de tempo. A Tabela [[#table-F|F]], lista as vazões iniciais na seção transversal a montante e a jusante em cada trecho de canal. Utilizou-se interpolação linear para determinar os valores de vazões iniciais a cada <math display="inline">\Delta x = 50</math> m no comprimento dos canais no instante inicial da simulação. Nos demais passos de tempo, o escoamento de base é incluído pela equação [[#eq-39|(39)]],
+O período de chuva analisado na bacia é dado no mês 07 de 1987 nos dias 07 e 08.  Neste mês não houve chuva além dos dias citados. A  [[#table-D|Tabela D]] contém a intensidade média de precipitação em cada sub-bacia. Para se obter uma simulação mais próxima da realidade, considerou-se uma vazão inicial que ocorre no rio (escoamento de base), em cada trecho de canal, antes de ocorrer a precipitação. Essa vazão entra no modelo computacional como uma condição inicial. Também é considerada como uma condição de contorno ao ser adicionada na contribuição lateral dos planos, em cada passo de tempo. A  [[#table-E|Tabela E]], lista as vazões iniciais na seção transversal a montante e a jusante em cada trecho de canal. Utilizou-se interpolação linear para determinar os valores de vazões iniciais a cada <math display="inline">\Delta x = 50</math> m no comprimento dos canais no instante inicial da simulação. Nos demais passos de tempo, o escoamento de base é incluído pela equação [[#eq-39|(39)]],
 <span id="eq-39"></span>
@@ Line 742: / Line 734: @@
 em que <math display="inline">Q_j</math> e <math display="inline">Q_m</math> são, respectivamente, a vazão a jusante e a vazão a montante em um segmento de canal de comprimento <math display="inline">L</math>. A quantidade <math display="inline">q_{base}</math> é somada a contribuição lateral e uniformemente distribuída no comprimento do canal.
-Os resultados obtidos pelo LBM são comparados com dados observados nas duas estações pertencentes ao trecho analisado. Nas Figuras [[#img-8|8]] e [[#img-9|9]], tem-se os hidrogramas simulado e observado nas estações Salto Claudelino e Porto Palmeirinha referentes a precipitação ocorrida no mês 07/1987.
+Os resultados obtidos pelo LBM são comparados com dados observados nas duas estações pertencentes ao trecho analisado. Nas  [[#img-8|Figuras 8]] e [[#img-9|9]], tem-se os hidrogramas simulado e observado nas estações Salto Claudelino e Porto Palmeirinha referentes a precipitação ocorrida no mês 07/1987.
 <div id='img-8'></div>
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
 |-
-|[[Image:Draft_Galina_237452870-Hidro_sub5.png|512px|Hidrograma na sub-bacia 5 - Estação Salto Claudelino]]
+|style="padding:10px;"| [[Image:Draft_Galina_237452870-Hidro_sub5.png|450px|Hidrograma na sub-bacia 5 - Estação Salto Claudelino]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figura 8:''' Hidrograma na sub-bacia 5 - Estação Salto Claudelino
+| colspan="1" style="padding:10px;" | '''Figura 8'''. Hidrograma na sub-bacia 5 - Estação Salto Claudelino
 |}
 <div id='img-9'></div>
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
 |-
-|[[Image:Draft_Galina_237452870-Hidro_sub9.png|512px|Hidrograma na sub-bacia 9 - Estação Porto Palmeirinha]]
+|style="padding:10px;"| [[Image:Draft_Galina_237452870-Hidro_sub9.png|450px|Hidrograma na sub-bacia 9 - Estação Porto Palmeirinha]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figura 9:''' Hidrograma na sub-bacia 9 - Estação Porto Palmeirinha
+| colspan="1" style="padding:10px;" | '''Figura 9'''. Hidrograma na sub-bacia 9 - Estação Porto Palmeirinha
 |}
@@ Line 776: / Line 768: @@
 onde X e Y são as variáveis em estudo, respectivamente, vazão observada e vazão simulada, <math display="inline">n</math> é o número de pares das observações.
-O coeficiente de correlação linear <math display="inline">r</math> é adimensional e tem variação no intervalo <math display="inline">[-1, 1]</math>. Se <math display="inline">r=-1</math> tem-se a correlação linear negativa perfeita, enquanto que se <math display="inline">r=1</math>, tem-se a correlação linear positiva perfeita. Para <math display="inline">r=0</math> não há correlação linear entre as variáveis. Esse coeficiente pode ser avaliado quantitativamente como descrito na Tabela [[#table-1|1]].
+O coeficiente de correlação linear <math display="inline">r</math> é adimensional e tem variação no intervalo <math display="inline">[-1, 1]</math>. Se <math display="inline">r=-1</math> tem-se a correlação linear negativa perfeita, enquanto que se <math display="inline">r=1</math>, tem-se a correlação linear positiva perfeita. Para <math display="inline">r=0</math> não há correlação linear entre as variáveis. Esse coeficiente pode ser avaliado quantitativamente como descrito na  [[#table-1|Tabela 1]].
+<div class="center" style="font-size: 75%;">'''Tabela 1'''. Classificação do coeficiente linear de Pearson</div>
-{|  class="floating_tableSCP wikitable" style="text-align: left; margin: 1em auto;min-width:50%;"
+<div id='table-1'></div>
-|+ style="font-size: 75%;" |<span id='table-1'></span>Tabla. 1 Classificação do coeficiente linear de Pearson
+{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"
-|- style="border-top: 2px solid;"
+|-style="text-align:center"
-| style="text-align: center;border-right: 2px solid;" |    Intervalo
+!Intervalo   !! style="text-align:left;" | Classificação
-| colspan='1' style="text-align: center;border-left: 2px solid;" | Classificação
+|-
-|- style="border-top: 2px solid;"
+| style="text-align: center;" | <math display="inline">0 < |r| < 0,3</math>
-| style="text-align: center;border-right: 2px solid;" | <math display="inline">0 < |r| < 0,3</math>
+| style="text-align: left;" | Existe fraca correlação linear
-| style="border-left: 2px solid;" | existe fraca correlação linear
 |-
-| style="text-align: center;border-right: 2px solid;" | <math>0,3 < |r| < 0,6</math>
+| style="text-align: center;"|  <math>0,3 < |r| < 0,6</math>
-| style="border-left: 2px solid;" | existe moderada correlação linear
+| style="text-align: left;" | Existe moderada correlação linear
 |-
-| style="text-align: center;border-right: 2px solid;" | <math>0,6 < |r| < 0,9</math>
+| style="text-align: center;" | <math>0,6 < |r| < 0,9</math>
-| style="border-left: 2px solid;" | existe forte correlação linear
+| style="text-align: left;" | Existe forte correlação linear
-|- style="border-bottom: 2px solid;"
+|-
-| style="text-align: center;border-right: 2px solid;" | <math>0,9 < |r| < 1,0</math>
+| style="text-align: center;" | <math>0,9 < |r| < 1,0</math>
-| style="border-left: 2px solid;" | existe correlação linear muito forte
+| style="text-align: left;" | Existe correlação linear muito forte
 |}
-<span style="text-align: center; font-size: 75%;">Fonte: Marques e Marques (2005) <span id='citeF-50'></span>[[#cite-50|[50]]].   </span>
+<div class="center" style="font-size: 75%;text-align: center;">'''Fonte''': Marques e Marques <span id='citeF-50'></span>[[#cite-50|[50]]]   </div>
-As Figuras [[#img-10|10]] e [[#img-11|11]] mostram os gráficos de dispersão das vazões observadas pelas vazões simuladas pelo LBM. A reta tracejada, em ambos os gráficos, representa o ajuste perfeito com coeficiente <math display="inline">r=1</math>.
+As [[#img-10|Figuras 10]] e [[#img-11|11]] mostram os gráficos de dispersão das vazões observadas pelas vazões simuladas pelo LBM. A reta tracejada, em ambos os gráficos, representa o ajuste perfeito com coeficiente <math display="inline">r=1</math>.
 <div id='img-10'></div>
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
 |-
-|[[Image:Draft_Galina_237452870-Dispersao_sub5.png|510px|Dispersão das vazões observadas e simuladas na sub-bacia 5 - Estação Salto Claudelino]]
+|style="padding-top:10px;"| [[Image:Draft_Galina_237452870-Dispersao_sub5.png|450px|Dispersão das vazões observadas e simuladas na sub-bacia 5 - Estação Salto Claudelino]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figura 10:''' Dispersão das vazões observadas e simuladas na sub-bacia 5 - Estação Salto Claudelino
+| colspan="1" style="padding:10px;"| '''Figura 10'''. Dispersão das vazões observadas e simuladas na sub-bacia 5 - Estação Salto Claudelino
 |}
 <div id='img-11'></div>
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
 |-
-|[[Image:Draft_Galina_237452870-Dispersao_sub9.png|510px|Dispersão das vazões observadas e simuladas na sub-bacia 9 - Estação Porto Palmeirinha]]
+|style="padding-top:10px;"| [[Image:Draft_Galina_237452870-Dispersao_sub9.png|450px|Dispersão das vazões observadas e simuladas na sub-bacia 9 - Estação Porto Palmeirinha]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figura 11:''' Dispersão das vazões observadas e simuladas na sub-bacia 9 - Estação Porto Palmeirinha
+| colspan="1" style="padding:10px;"| '''Figura 11'''. Dispersão das vazões observadas e simuladas na sub-bacia 9 - Estação Porto Palmeirinha
 |}
-O coeficiente de correlação linear de Pearson calculado na sub-bacia 5 foi de <math display="inline">r=0,88</math> e na sub-bacia 9, de <math display="inline">r=0,97</math>. Conforme a Tabela [[#5 Resultados e Discussões|5]], existe forte correlação linear entre as variáveis na sub-bacia 5 e correlação linear muito forte para a sub-bacia 9. Conclui-se que há boa concordância entre as vazões observadas e simuladas. Nota-se na Figura [[#img-10|10]] que um ponto está distante da reta, o que faz o coeficiente <math display="inline">r</math> ser um pouco menor.
+O coeficiente de correlação linear de Pearson calculado na sub-bacia 5 foi de <math display="inline">r=0,88</math> e na sub-bacia 9, de <math display="inline">r=0,97</math>. Conforme a  [[#table-1|Tabela 1]], existe forte correlação linear entre as variáveis na sub-bacia 5 e correlação linear muito forte para a sub-bacia 9. Conclui-se que há boa concordância entre as vazões observadas e simuladas. Nota-se na  [[#img-10|Figura 10]] que um ponto está distante da reta, o que faz o coeficiente <math display="inline">r</math> ser um pouco menor.
-Para fazer a verificação da significância estatística, foi aplicado o teste de hipóteses com a distribuição ''t'' de Student, supondo que os desvios padrões populacionais são desconhecidos e com variâncias equivalentes. A hipótese a ser testada é <math display="inline">H_0: \mu _1 - \mu _2 = 0</math>, isto é, não existe diferença significativa entre as médias populacionais das vazões observadas e simuladas. Os resultados obtidos, para as sub-bacias 5 e 9, estão resumidos na Tabela [[#table-A|A]].
+Para fazer a verificação da significância estatística, foi aplicado o teste de hipóteses com a distribuição ''t'' de Student, supondo que os desvios padrões populacionais são desconhecidos e com variâncias equivalentes. A hipótese a ser testada é <math display="inline">H_0: \mu _1 - \mu _2 = 0</math>, isto é, não existe diferença significativa entre as médias populacionais das vazões observadas e simuladas. Os resultados obtidos, para as sub-bacias 5 e 9, estão resumidos na  [[#table-F|Tabela F]].
-Os testes realizados mostram que a estatística do teste ''t'' é menor do que o valor do crítico bi-caudal, portanto aceita-se a hipótese <math display="inline">H_0</math>, ou seja, não existe diferença significativa, ao nível de significância de <math display="inline">5%</math> entre as vazões observadas e vazões simuladas.
+Os testes realizados mostram que a estatística do teste <math display="inline">t</math> é menor do que o valor do crítico bi-caudal, portanto aceita-se a hipótese <math display="inline">H_0</math>, ou seja, não existe diferença significativa, ao nível de significância de <math display="inline">5%</math> entre as vazões observadas e vazões simuladas.
-==6 Conclusão==
+==6. Conclusão==
 Esta pesquisa utilizou o LBM com aperador de colisão BGK para obter a simulação numérica do escoamento superficial da bacia do rio Chopim. Computacionalmente a bacia foi representada pela junção de nove sub-bacias dadas por planos laterais e um segmento de canal. A simulação produziu o hidrograma total o qual fornece a variação da vazão em relação ao tempo correspondente a resposta da bacia hidrográfica após a precipitação efetiva do evento simulado. Para que este resultado fosse alcançado, determinou-se duas funções de distribuição de equilíbrio obtidas por meio da expansão de Chapmann-Enskog em escalas de tempo e utilizando o reticulado D1Q5, uma adequada para o escoamento na superfície da bacia e outra para o canal principal, obtendo a profundidade da água na superfície da bacia e área da seção transversal nos canais. Também foi necessário determinar a condição de contorno na passagem de fluxo de um a sub-bacia para outra, levando em consideração a conservação da massa.
@@ Line 831: / Line 824: @@
 Verificou-se que LBM é adequado e preciso na simulação do escoamento em uma bacia hidrográfica representado pelas equações de onda cinemática. Além disso, as condições iniciais e de contorno são de fácil implementação.
-===BIBLIOGRAFÍA===
+==Referências==
+<div class="auto" style="text-align: left;width: auto; margin-left: auto; margin-right: auto;font-size: 85%;">
 <div id="cite-1"></div>
-'''[[#citeF-1|[1]]]'''  J. E. Gribbin, Introduction to Hydraulics and Hydrology with Applications for Stormwater Management, 3rd Edition, Cengage Learning, USA, 2008.
+[[#citeF-1|[1]]]   Gribbin J.E. Introduction to hydraulics and hydrology with applications for stormwater management. 3rd Edition, Cengage Learning, USA, 2008.
 <div id="cite-2"></div>
-'''[[#citeF-2|[2]]]'''  J. G. Zhou, Lattice Boltzmann Method for Shallow Water Flows, Springer, New York, 2004.
+[[#citeF-2|[2]]]   Zhou J.G. Lattice Boltzmann method for shallow water flows. Springer, New York, 2004.
 <div id="cite-3"></div>
-'''[[#citeF-3|[3]]]'''  S. Chen, G. D. Doolen. Lattice Boltzmann Method for fluid flows, ''Annual Review of Fluid Mechanics'', vol. 30, no. 1, pp. 329-364, 1998.
+[[#citeF-3|[3]]]   Chen S.,  Doolen G.D. Lattice Boltzmann method for fluid flows. Annual Review of Fluid Mechanics, 30(1):329-364, 1998.
 <div id="cite-4"></div>
-'''[[#citeF-4|[4]]]'''  J. G. Zhou, A lattice Boltzmann model for the shallow water equations, Computer Methods in Applied Mechanics and Engineering 191 (2002) 3527-3539.
+[[#citeF-4|[4]]]  Zhou J.G. A lattice Boltzmann model for the shallow water equations. Computer Methods in Applied Mechanics and Engineering, 191:3527-3539, 2002.
 <div id="cite-5"></div>
-'''[[#citeF-5|[5]]]'''  J. G. Zhou, H. Liu, S. Shafiai, Y. Peng, R. Burrows, Lattice Boltzmann Method for Open-channel Flows, Engineering and Computational Mechanics 163 (2010) 243-249.
+[[#citeF-5|[5]]]  Zhou J.G.,   Liu H., Shafiai S.,  Peng Y.,  Burrows R. Lattice Boltzmann method for open-channel flows. Engineering and Computational Mechanics, 163:243-249, 2010.
 <div id="cite-6"></div>
-'''[[#citeF-6|[6]]]'''  P. van Thang, B. Chopard, L. Lefèvre, D. A. Ondo, E. Mendes, Study of the 1D Lattice Boltzmann Shallow Water Equation and Its Coupling to Build a Canal Network, Journal of Computational Physics 229 (19) (2010) 7373-7400.
+[[#citeF-6|[6]]]   van Thang P.,  Chopard B.,  Lefèvre L.,  Ondo D.A.,  Mendes E. Study of the 1D lattice Boltzmann shallow water equation and its coupling to build a canal network. Journal of Computational Physics, 229(19):7373-7400, 2010.
 <div id="cite-7"></div>
-'''[[#citeF-7|[7]]]'''  Y. Peng, J. G. Zhou, R. Burrows, Modelling solute transport in shallow water with the lattice boltzmann method, Computers & Fluids 50 (2011) 181?188.
+[[#citeF-7|[7]]]   Peng Y.,  Zhou J.G.,  Burrows R. Modelling solute transport in shallow water with the lattice Boltzmann method. Computers & Fluids, 50:181-188, 2011.
 <div id="cite-8"></div>
-'''[[#citeF-8|[8]]]'''  Y. Peng, J. G. Zhou, J. M. Zhang, R. Burrows, modeling moving boundary in shallow water by lbm, International Journal of Modern Physics C 24 (1).
+[[#citeF-8|[8]]]  Peng Y.,  Zhou J.G.,  Zhang J.M.,  Burrows R. Modeling moving boundary in shallow water by lbm. International Journal of Modern Physics, 24(1):1250094, 2013.
 <div id="cite-9"></div>
-'''[[#citeF-9|[9]]]'''  L. Budinski, Mrt lattice boltzmann method for 2d flows in curvilinear coordinates, Computers & Fluids 94 (2014) 288?301.
+[[#citeF-9|[9]]]   Budinski L. Mrt lattice Boltzmann method for 2d flows in curvilinear coordinates. Computers & Fluids, 94:288?301, 2014.
 <div id="cite-10"></div>
-'''[[#citeF-10|[10]]]'''  Y. Peng, J. G. Zhou, J. M. Zhang, H. Liu, Lattice boltzmann modeling of shallow water flows over discontinuous beds, International Journal for Numerical Methods in Fluids 75 (2014) 608-619.
+[[#citeF-10|[10]]]  Peng Y.,  Zhou J.G.,  Zhang J.M., Liu H. Lattice Boltzmann modeling of shallow water flows over discontinuous beds. International Journal for Numerical Methods in Fluids, 75:608-619, 2014.
 <div id="cite-11"></div>
-'''[[#citeF-11|[11]]]'''  M. La Rocca, A. Montessori, P. Prestininzi, S. Succi, A multispeed discrete boltzmann model for transcritical 2d shallow water flows, Journal of Computational Physics, 284 (2015) 117-132.
+[[#citeF-11|[11]]]   La Rocca M.,  Montessori A.,  Prestininzi P.,  Succi S. A multispeed discrete Boltzmann model for transcritical 2d shallow water flows. Journal of Computational Physics, 284:117-132, 2015.
 <div id="cite-12"></div>
-'''[[#citeF-12|[12]]]'''  S. Li, P. Huang, Li, A modified lattice boltzmann model for shallow water flows over complex topography, International Journal for Numerical Methods in Fluids 77 (2015) 441-458.
+[[#citeF-12|[12]]]   Li S.,  Huang P., Li J. A modified lattice Boltzmann model for shallow water flows over complex topography. International Journal for Numerical Methods in Fluids, 77(8):441-458, 2015.
 <div id="cite-13"></div>
-'''[[#citeF-13|[13]]]'''  P. Prestininzi, a. Montessori, M. La Rocca, G. Sciortino, Simulation of arrested salt wedges with a multi-layer shallow water lattice boltzmann model, Advances in Water Resources 96 (2016) 282-289.
+[[#citeF-13|[13]]]   Prestininzi P.,  Montessori A.,  La Rocca M.,  Sciortino G. Simulation of arrested salt wedges with a multi-layer shallow water lattice Boltzmann model. Advances in Water Resources, 96:282-289, 2016.
 <div id="cite-14"></div>
-'''[[#citeF-14|[14]]]'''  C.-z. Zhang, C. Y.-G., J.-Y. Wu, W. Diao, Lattice boltzmann simulation of the open channel flow connecting two cascaded hydropower stations, Journal of Hydrodynamics 28 (3) (2016) 400-410.
+[[#citeF-14|[14]]]   Zhang C.-Z., Cheng Y.-G.,  Wu J.-Y.,  Diao W. Lattice Boltzmann simulation of the open channel flow connecting two cascaded hydropower stations. Journal of Hydrodynamics, 28(3):400-410, 2016.
 <div id="cite-15"></div>
-'''[[#citeF-15|[15]]]'''  V. Galina, J. Cargnelutti, E. Kaviski, L. M. Gramani, A. M. Lobeiro, Simulação de onda de maré por meio do método do reticulado de boltzmann, in: I Simpósio de Métodos Numéricos em Engenharia, Universidade Federal do Paraná, Curitiba - Brasil, 2016.
+[[#citeF-15|[15]]]   Galina V.,  Cargnelutti J.,  Kaviski E.,  Gramani L.M.,  Lobeiro A.M. Simulação de onda de maré por meio do método do reticulado de Boltzmann. In: I Simpósio de Métodos Numéricos em Engenharia, Universidade Federal do Paraná, Curitiba, Brasil, 2016.
 <div id="cite-16"></div>
-'''[[#citeF-16|[16]]]'''  Z.-m. Zhao, P. Huang, S.-t. Li, Lattice boltzmann model for shallow water in curvilinear coordinate grid, Journal of Hydrodynamics 29 (2) (2017) 251-260.
+[[#citeF-16|[16]]]   Zhao Z.-M.,  Huang P.,  Li S.-T. Lattice Boltzmann model for shallow water in curvilinear coordinate grid. Journal of Hydrodynamics, 29(2):251-260, 2017.
 <div id="cite-17"></div>
-'''[[#citeF-17|[17]]]'''  A. De Rosis, A central moments-based lattice boltzmann scheme for shallow water equations, Computer Methods in Applied Mechanics and Engineering 319 (2017) 379-392.
+[[#citeF-17|[17]]]   De Rosis A. A central moments-based lattice Boltzmann scheme for shallow water equations. Computer Methods in Applied Mechanics and Engineering, 319:379-392, 2017.
 <div id="cite-18"></div>
-'''[[#citeF-18|[18]]]'''  J. Cargnelutti, V. Galina, E. Kaviski, L. M. Gramani, A. M. Lobeiro, Two-dimensional numerical simulation of channel flow with submerged obstacles using the lattice boltzmann method, Revista Internacional de Métodos Numéricos para Cálculo y Diseño en Ingeniería 34 (1), 2018.
+[[#citeF-18|[18]]]  Cargnelutti J., Galina V., Kaviski E., Gramani L.M.,  Lobeiro A.M. Two-dimensional numerical simulation of channel flow with submerged obstacles using the lattice Boltzmann method. Revista Internacional de Métodos Numéricos para Cálculo y Diseño en Ingeniería, 34(1), 20, 2018.
 <div id="cite-19"></div>
-'''[[#citeF-19|[19]]]'''  J. Cargnelutti, V. Galina, E. Kaviski, L. M. Gramani, A. M. Lobeiro, Simulation of the two-dimensional flow of the initiation channel of the itaipu hydroelectric power plant by the lattice boltzmann method, Revista Internacional de Métodos Numéricos para Cálculo y Diseño en Ingeniería 34 (1), 2018.
+[[#citeF-19|[19]]]  Cargnelutti J., Galina V., Kaviski E., Gramani L.M.,  Lobeiro A.M. Simulation of the two-dimensional flow of the initiation channel of the Itaipu hydroelectric power plant by the lattice Boltzmann method. Revista Internacional de Métodos Numéricos para Cálculo y Diseño en Ingeniería, 34(1), 25, 2018.
 <div id="cite-20"></div>
-'''[[#citeF-20|[20]]]'''  G. Yang, C. Kwok, Y. Sobral, The effects of bed form roughness on total suspended load via the lattice boltzmann method, Applied Mathematical Modelling 63 (2018) 591-610.
+[[#citeF-20|[20]]]   Yang G.,  Kwok C.,  Sobral Y. The effects of bed form roughness on total suspended load via the lattice Boltzmann method. Applied Mathematical Modelling, 63:591-610, 2018.
 <div id="cite-21"></div>
-'''[[#citeF-21|[21]]]'''  X. Zhang, J. Feng, D. Zhang, N. Liu, Comparison of lattice boltzmann method and preissmann implicit difference method in application to overland flow, Nongye Jixie Xuebao/Transactions of the Chinese Society of Agricultural Machinery, 2014
+[[#citeF-21|[21]]]   Zhang X.,  Feng J.,  Zhang D.,  Liu N. Comparison of lattice Boltzmann method and preissmann implicit difference method in application to overland flow. Nongye Jixie Xuebao/Transactions of the Chinese Society of Agricultural Machinery, 2014.
 <div id="cite-22"></div>
-'''[[#citeF-22|[22]]]'''  X. Zhang, J. Feng, T. Yang, Lattice Boltzmann Method for Overland Flow Studies and its Experimental Validation, Journal of Hydraulic Research 53 (5) (2015) 561-575.
+[[#citeF-22|[22]]]   Zhang X.,  Feng J.,  Yang T. Lattice Boltzmann method for overland flow studies and its Eeperimental validation. Journal of Hydraulic Research, 53(5):561-575, 2015.
 <div id="cite-23"></div>
-'''[[#citeF-23|[23]]]'''  N. Liu, The numerical simulation of one-dimensional overland flow by Lattice Boltzmann Method, in: 5th International Conference on Advanced Design and Manufacturing Engineeringe, China, 2015.
+[[#citeF-23|[23]]]   Liu N. The numerical simulation of one-dimensional overland flow by lattice Boltzmann method. In: 5th International Conference on Advanced Design and Manufacturing Engineeringe, China, 2015.
 <div id="cite-24"></div>
-'''[[#citeF-24|[24]]]'''  V. Galina, J. Cargnelutti, E. Kaviski, L. M. Gramani, A. M. Lobeiro, Application of lattice boltzmann method for surface runoff in watershed, Revista Internacional de Métodos Numéricos para Cálculo y Diseño en Ingeniería 34 (1), 2018.
+[[#citeF-24|[24]]]  Galina V.,  Cargnelutti J.,  Kaviski E.,  Gramani L.M.,  Lobeiro A.M. Application of lattice Boltzmann method for surface runoff in watershed. Revista Internacional de Métodos Numéricos para Cálculo y Diseño en Ingeniería, 34(1), 10, 2018.
 <div id="cite-25"></div>
-'''[[#citeF-25|[25]]]'''  H. Liu, H. Wang, S. Liu, C. Hu, Y. Ding, J. Zhang, Lattice Boltzmann Method for the Saint?Venant Equations, Journal of Hydrology 524 (2015) 411-416.
+[[#citeF-25|[25]]]   Liu H.,  Wang H.,  Liu S.,  Hu C.,  Ding Y.,  Zhang J. Lattice Boltzmann method for the Saint-Venant equations. Journal of Hydrology, 524:411-416, 2015.
 <div id="cite-26"></div>
-'''[[#citeF-26|[26]]]'''  R. Benzi, S. Succi, M. Vergassola, The lattice boltzmann equation: Theory and applications, Physics Reports 222 (3) (1992) 145-197.
+[[#citeF-26|[26]]]   Benzi R.,  Succi S.,  Vergassola M. The lattice Boltzmann equation: Theory and applications. Physics Reports, 222(3):145-197, 1992.
 <div id="cite-27"></div>
-'''[[#citeF-27|[27]]]'''  B. Nestler, A. Aksi, M. Selzer, Combined lattice boltzmann and phase-field simulations for incompressible fluid flow in porous media, Mathematics and Computers in Simulation 80 (7) (2010) 1458-1468.
+[[#citeF-27|[27]]]   Nestler B.,  Aksi A.,  Selzer M. Combined lattice Boltzmann and phase-field simulations for incompressible fluid flow in porous media. Mathematics and Computers in Simulation, 80(7):1458-1468, 2010.
 <div id="cite-28"></div>
-'''[[#citeF-28|[28]]]'''  W. Sobieski, Numerical investigations of tortuosity in randomly generated pore structures, Mathematics and Computers in Simulation 166 (2019) 1?20.
+[[#citeF-28|[28]]]   Sobieski W. Numerical investigations of tortuosity in randomly generated pore structures. Mathematics and Computers in Simulation, 166:1-20, 2019.
 <div id="cite-29"></div>
-'''[[#citeF-29|[29]]]'''  M. E. Andaluz, V. V. Galarza, A. R. Vera, On hydraulic tortuosity variations due to morphological considerations in 2d porous media by using the lattice boltzmann method, Mathematics and Computers in Simulation 169 (2020) 74?87.
+[[#citeF-29|[29]]]   Andaluz M.E.,  Galarza V.V.,  Vera A.R. On hydraulic tortuosity variations due to morphological considerations in 2d porous media by using the lattice Boltzmann method. Mathematics and Computers in Simulation, 169:74-87, 2020.
 <div id="cite-30"></div>
-'''[[#citeF-30|[30]]]'''  Y. H. Qian,  D. D'Humières, P. Lallemand, Lattice bgk models for navier-stokes equation, Europhysics Letters 17 (6) (1992) 479-484.
+[[#citeF-30|[30]]]   Qian Y.H., D'Humières D.,  Lallemand P. Lattice BGK models for Navier-Stokes equation. Europhysics Letters, 17(6):479-484, 1992.
 <div id="cite-31"></div>
-'''[[#citeF-31|[31]]]'''  R. M. Porto, Hidráulica Básica, 4th Edition, EESC-USP, São Carlos, 2006.
+[[#citeF-31|[31]]]   Porto R.M. Hidráulica básica. 4th Edition, EESC-USP, São Carlos, 2006.
 <div id="cite-32"></div>
-'''[[#citeF-32|[32]]]'''  M. H. Chaudhry, Open-Channel Flow, Springer, New York, USA, 2008.
+[[#citeF-32|[32]]]   Chaudhry M.H. Open-channel flow. Springer, New York, USA, 2008.
 <div id="cite-33"></div>
-'''[[#citeF-33|[33]]]'''  C. E. M. Tucci, Modelos Hidrológicos, Editora da Universidade Federal do Rio Grande do Sul, Porto Alegre, RS - Brasil, 1998.
+[[#citeF-33|[33]]]   Tucci C.E.M. Modelos hidrológicos. Editora da Universidade Federal do Rio Grande do Sul, Porto Alegre, RS - Brasil, 1998.
 <div id="cite-34"></div>
-'''[[#citeF-34|[34]]]'''  J. E. Miller, Basic Concepts of Kinematic-Wave Models, United States Government Printing office, Washington, 1984.
+[[#citeF-34|[34]]]   Miller J.E. Basic concepts of kinematic-wave models. United States Government Printing office, Washington, 1984.
 <div id="cite-35"></div>
-'''[[#citeF-35|[35]]]'''  D. R. Golbert, Método de Lattice Boltzmann em Hemodinâmica Computacional: interações fluido-estrutura e modelos acoplados 1D-3D, Ph.D. thesis,  Laboratório Nacional de Computação Científica, Petrópolis, RJ - Brasil.
+[[#citeF-35|[35]]]   Golbert D.R. Método de lattice Boltzmann em hemodinâmica computacional: interações fluido-estrutura e modelos acoplados 1D-3D. Ph.D. Thesis,  Laboratório Nacional de Computação Científica, Petrópolis, RJ - Brasil, 2013.
 <div id="cite-36"></div>
-'''[[#citeF-36|[36]]]'''  P. L. Bhatnagar, E. P. Gross, M. Krook, A Model for Collision Processes in Gases. I. Small Amplitude Processes in Charged and Neutral One-Component Systems, Physical Review 94 (3) (1954) 511-525.
+[[#citeF-36|[36]]]   Bhatnagar P.L.,  Gross E.P.,  Krook M. A model for collision processes in gases. I. Small amplitude processes in charged and neutral one-component systems. Physical Review, 94(3):511-525, 1954.
 <div id="cite-37"></div>
-'''[[#citeF-37|[37]]]'''  Z. Guo, C. Zheng, B. Shi, Discrete Lattice Effects on the Forcing Term in the Lattice Boltzmann Method, Physical Review E 65 (2002) 1-6.
+[[#citeF-37|[37]]]   Guo Z.,  Zheng C.,  Shi B. Discrete lattice effects on the forcing term in the lattice Boltzmann method. Physical Review E, 65(4):1-6, 2002.
 <div id="cite-38"></div>
-'''[38]'''  J. Zhang, G. Yan, Lattice boltzmann method for one and two-dimensional burgers equation, Physica A (2008).
+[[#citeF-37|[38]]]   Zhang J.,  Yan G. Lattice Boltzmann method for one and two-dimensional burgers equation. Physica A: Statistical Mechanics and its Applications, 387(19–20):4771-4786, 2008.
 <div id="cite-39"></div>
-'''[[#citeF-39|[39]]]'''  S. Chapman, T. G. Cowling, The Mathematical Theory of Non-uniform Gases. An Account of the Kinetic Theory of Viscosity, Thermal Conduction and Diffusion in Gases, 3rd Edition, Cambridge University Press, London, 1970.
+[[#citeF-39|[39]]]  Chapman S.,  Cowling T.G. The mathematical theory of non-uniform gases. An account of the kinetic theory of viscosity, thermal conduction and diffusion in gases. 3rd Edition, Cambridge University Press, London, 1970.
 <div id="cite-40"></div>
-'''[[#citeF-40|[40]]]'''  J. Zhang, G. Yan, A lattice Boltzmann model for the Burgers-Fisher equation, Chaos: An Interdisciplinary Journal of Nonlinear Science  20 (2010) 023129.
+[[#citeF-40|[40]]]   Zhang J.,  Yan G. A lattice Boltzmann model for the Burgers-Fisher equation. Chaos: An Interdisciplinary Journal of Nonlinear Science,  20:023129, 2010.
 <div id="cite-41"></div>
-'''[[#citeF-41|[41]]]'''  J. D. Sterling, S. Chen, Stability analysis of lattice boltzmann methods, Journal of Computational Physics 123 (1996) 196-206.
+[[#citeF-41|[41]]]   Sterling J.D.,  Chen S. Stability analysis of lattice Boltzmann methods. Journal of Computational Physics, 123:196-206, 1996.
 <div id="cite-42"></div>
-'''[[#citeF-42|[42]]]'''  M. La Rocca, C. Adduce, G. Sciortino, Development of a Lattice Boltzmann Method for Two-layered Shallow-water Flow, International Journal for Numerical Methods in Fluids 70 (2012) 1048-1072.
+[[#citeF-42|[42]]]   La Rocca M.,  Adduce C.,  Sciortino G. Development of a lattice Boltzmann method for two-layered shallow-water flow. International Journal for Numerical Methods in Fluids, 70:1048-1072, 2012.
 <div id="cite-43"></div>
-'''[[#citeF-43|[43]]]'''  J. G. Zhou, H. Liu, Determination of bed elevation in the enhanced lattice Boltzmann method for the shallow-water equations, Physical Review E 88 (2013) 1-6.
+[[#citeF-43|[43]]]   Zhou J.G.,  Liu H. Determination of bed elevation in the enhanced lattice Boltzmann method for the shallow-water equations. Physical Review E, 88(2):1-6, 2013.
 <div id="cite-44"></div>
-'''[[#citeF-44|[44]]]'''  S. de Estado do Meio Ambiente e Recursos Hídricos SEMA, Bacias hidrográficas do paraná, Paraná, Brasil (2010).
+[[#citeF-44|[44]]]  Secretaria de Estado do Meio Ambiente e Recursos Hídricos SEMA, Bacias hidrográficas do Paraná, Brasil, 2010.
 <div id="cite-45"></div>
-'''[45]'''  A. C. L. Franco, Calibração do modelo swat com evapotranspiração proveniente de sensoriamento remoto e vazão observada, Master?s thesis, Florianópolis - SC (2017).
+[45]   Franco A.C.L. Calibração do modelo swat com evapotranspiração proveniente de sensoriamento remoto e vazão observada. Master's thesis, Florianópolis - SC, 2017.
 <div id="cite-46"></div>
-'''[[#citeF-46|[46]]]'''  D. I. Chavasse, R. S. Seoane, Assessing and predicting the impact of el nino southern oscillation (enso) events on runoff from the chopim river basin, brazil, Hydrological Processes (23) (2009) 3261-3266.
+[[#citeF-46|[46]]]   Chavasse D.I.,  Seoane R.S. Assessing and predicting the impact of el Nino Southern Oscillation (ENSO) events on runoff from the Chopim river basin, Brazil. Hydrological Processes, 23(22):3261-3266, 2009.
 <div id="cite-47"></div>
-'''[[#citeF-47|[47]]]'''  P. P. Burian, Do estudo de impacto embiental à avaliação ambiental estratégica - ambivalências do processo de licenciamento ambiental do setor elétrico, Ph.D. thesis, Universidade Estadual de Campinas, Campinas, SP - Brasil (Fevereiro 2006).
+[[#citeF-47|[47]]]   Burian P.P. Do estudo de impacto embiental à avaliação ambiental estratégica - ambivalências do processo de licenciamento ambiental do setor elétrico. Ph.D. Thesis, Universidade Estadual de Campinas, Campinas, SP - Brasil, Fevereiro 2006.
 <div id="cite-48"></div>
-'''[48]'''  P. K. Uda, Evapotranspiração real da bacia do rio iguaçu por meio do modelo metric, Ph.D. thesis (2016).
+[[#citeF-37|[48]]]   Uda P.K. Evapotranspiração real da bacia do rio iguaçu por meio do modelo metric. Ph.D. Thesis, Universidade Federal de Santa Catarina, Florianópolis, 2016.
 <div id="cite-49"></div>
-'''[[#citeF-49|[49]]]'''  E. Kaviski, L. F. A. Gonçalves, Análise hidrológica e matemática de operação de reservatórios, Tech. rep., Centro de Hidráulica e Hidrologia Professor Parigot de Souza - CEHPAR, Curitiba - PR (1996).
+[[#citeF-49|[49]]]   Kaviski E.,  Gonçalves L.F.A. Análise hidrológica e matemática de operação de reservatórios tech. rep. Centro de Hidráulica e Hidrologia Professor Parigot de Souza - CEHPAR, Curitiba - PR, 1996.
 <div id="cite-50"></div>
-'''[[#citeF-50|[50]]]'''  J. M. Marques, M. A. M. Marques, Estatística Básica para os Cursos de Engenharia, Domínio do Saber, 2005.
+[[#citeF-50|[50]]]   Marques J.M.,  Marques M.A.M. Estatística básica para os cursos de engenharia. Domínio do Saber, 2005.
+</div>
 ==Apêndice==
+<div class="center" style="font-size: 75%;">'''Tabela A'''. A estações hidrométricas no rio Chopim</div>
-{|  class="floating_tableSCP wikitable" style="text-align: center; margin: 1em auto;min-width:50%;"
+<div id='table-A'></div>
-|+ style="font-size: 75%;" |<span id='table-A'></span>Tabla. A Teste ''t'' de Student para as sub-bacias no rio Chopim
+{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"
-|- style="border-top: 2px solid;"
+|-style="text-align:center"
-| rowspan='2' colspan='1' style="text-align: left;border-right: 2px solid;border-right: 2px solid;" | Quantidades estatísticas
+!  Nome !! Código !! Município !! Latitude !! Longitude !! Altitude <br> (m) !! Área de contribuição <br> (km<math display="inline">^2</math>)
-| colspan='2' style="border-right: 2px solid;border-left: 2px solid;" |
+|-style="text-align:center"
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
+| style="text-align:left" | Porto <br> Palmeirinha
-|-
+|  65927000
-| Sub-bacia 5
+| Coronel Vivida
-|-
+| -26,03
-| Salto Claudelino
+|  -52,63
+|  501
+|  3.390
+|- style="text-align:center"
+|  style="text-align:left" | Salto <br>  Claudelino
+|  65925000
+| Clevelândia
+|  -26,28
+|  -52,30
+|  797
+|  1.660
 |}
-| colspan='2' style="border-left: 2px solid;" |
+<div class="center" style="font-size: 75%;"> '''Fonte''': Franco  <span id='citeF-45'></span>[[#cite-45|[45]]]</div>
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
-|-
-| Sub-Bacia 9
-|-
-| Porto Palmeirinha
-|}
-|-
-| style="border-left: 2px solid;" | Observado
-| style="border-right: 2px solid;" | Simulado
-| colspan='1' style="text-align: left;border-left: 2px solid;" | Observado
-| colspan='1' style="text-align: left;" | Simulado
-|- style="border-top: 2px solid;"
-| style="text-align: left;border-right: 2px solid;" |    Média
-| style="border-left: 2px solid;" | 80,04
-| style="border-right: 2px solid;" | 76,58
-| style="border-left: 2px solid;" | 142,12
-| 139,02
-|-
-| style="text-align: left;border-right: 2px solid;" | Variância
-| style="border-left: 2px solid;" | 1.493,98
-| style="border-right: 2px solid;" | 3.021,13
-| style="border-left: 2px solid;" | 3.376,02
-| 4.864,39
-|-
-| style="text-align: left;border-right: 2px solid;" | Número de observações
-| colspan='2' style="border-right: 2px solid;border-left: 2px solid;" | 21
-| colspan='2' style="border-left: 2px solid;" | 21
-|-
-| style="text-align: left;border-right: 2px solid;" | Graus de liberdade
-| colspan='2' style="border-right: 2px solid;border-left: 2px solid;" | 40
-| colspan='2' style="border-left: 2px solid;" | 40
-|-
-| style="text-align: left;border-right: 2px solid;" | Estatística do teste
-| colspan='2' style="border-right: 2px solid;border-left: 2px solid;" | 0,2361
-| colspan='2' style="border-left: 2px solid;" | 0,1567
-|-
-| style="text-align: left;border-right: 2px solid;" | ''t'' bi-caudal
-| colspan='2' style="border-right: 2px solid;border-left: 2px solid;" | 2,0211
-| colspan='2' style="border-left: 2px solid;" | 2,0211
-|- style="border-bottom: 2px solid;"
-| style="text-align: left;border-right: 2px solid;" | Nível de significância
-| colspan='2' style="border-right: 2px solid;border-left: 2px solid;" | <math>5%</math>
-| colspan='2' style="border-left: 2px solid;" | <math>5%</math>
-|}
+<div class="center" style="font-size: 75%;">'''Tabela B'''. Caracterização física dos planos das 9 sub-bacias do rio Chopim - pd e pe são, respectivamente, os planos do lado direito e do lado esquerdo no sentido do escoamento</div>
+<div id='table-B'></div>
-{|  class="floating_tableSCP wikitable" style="text-align: center; margin: 1em auto;min-width:50%;"
+{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"
-|+ style="font-size: 75%;" |<span id='table-B'></span>Tabla. B Estações hidrométricas no rio Chopim
+|-style="text-align:center"
-|- style="border-top: 2px solid;"
+!Sub-Bacia !! Plano !! Comprimento <br> <math>L</math> (m) !! Largura <br> <math>b</math> (m) !! Declividade <br> <math>S_0</math> (m/m)!! Manning <br> <math>n_r</math> (m<math>^{-1/3}</math>s) !! Área <br>  (km<math>^2</math>)
-| colspan='1' style="border-right: 2px solid;border-right: 2px solid;" | Nome
+|- style="text-align:center"
-| style="border-left: 2px solid;border-right: 2px solid;" | Código
+| rowspan='2' | 1
-| style="border-left: 2px solid;border-right: 2px solid;" | Município
+|  pd
-| colspan='1' style="border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" |
+|  -
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
+|  -
-|-
+|  -
-| Lati-
+| -
-|-
+| -
-| tude
+|-style="text-align:center"
+|  pe
+|  6.940
+|  26.830
+| 0,0065
+|  0,1200
+|  186,20
+|- style="text-align:center"
+| rowspan='2'  | 2
+|  pd
+|  3.080
+| 47.500
+| 0,0266
+|  0,0840
+| 146,30
+|-style="text-align:center"
+|  pe
+|  6.940
+|  47.500
+|  0,0065
+|  0,1200
+|  329,65
+|- style="text-align:center"
+| rowspan='2'  | 3
+|  pd
+|  3.300
+|  30.170
+|  0,0146
+|  0,1044
+|  99,56
+|-style="text-align:center"
+| pe
+|  6.940
+|  30.170
+|  0,0065
+|  0,1200
+| 209,38
+|- style="text-align:center"
+| rowspan='2' | 4
+| pd
+|  2.800
+|  26.580
+|  0,0102
+|  0,2832
+| 74,42
+|-style="text-align:center"
+| pe
+|  9.250
+|  26.580
+|  0,0129
+|  0,2832
+| 245,87
+|- style="text-align:center"
+| rowspan='2' | 5
+| pd
+|  2.800
+|  41.670
+|  0,0102
+|  0,2832
+| 116,68
+|-style="text-align:center"
+| pe
+|  6.780
+|  41.670
+|  0,0022
+|  0,1020
+| 282,52
+|- style="text-align:center"
+| rowspan='2' | 6
+| pd
+|  2.800
+|  13.000
+|  0,0102
+|  0,2832
+| 36,40
+|-style="text-align:center"
+| pe
+|  6.780
+|  13.000
+|  0,0022
+|  0,1020
+| 88,14
+|- style="text-align:center"
+| rowspan='2' | 7
+| pd
+|  5.180
+|  65.000
+|  0,0170
+|  0,0972
+| 336,70
+|-style="text-align:center"
+| pe
+|  9.670
+|  65.000
+|  0,0085
+|  0,0972
+| 628,55
+|- style="text-align:center"
+| rowspan='2' | 8
+| pd
+|  10.120
+|  22.000
+|  0,0113
+|  0,0984
+| 222,64
+|-style="text-align:center"
+| pe
+|  4.870
+|  22.000
+|  0,0129
+|  0,1020
+| 107,14
+|- style="text-align:center"
+| rowspan='2' | 9
+| pd
+|  6.720
+|  8.180
+|  0,0127
+|  0,1032
+| 54,97
+|- style="text-align:center"
+| pe
+|  4.870
+|  8.180
+|  0,0129
+|  0,1020
+| 39,84
 |}
-| colspan='1' style="border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" |
+<div class="center" style="font-size: 75%;">'''Fonte''': Kaviski e Gonçalves  <span id='citeF-49'></span>[[#cite-49|[49]]]</div>
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
-|-
-| Longi-
-|-
-| tude
-|}
-| style="border-left: 2px solid;border-right: 2px solid;" |
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
-|-
-| Altitude
-|-
-| (m)
-|}
-| style="border-left: 2px solid;" |
+<div class="center" style="font-size: 75%;">'''Tabela C'''. Caracterização física dos canais das 9 sub-bacias do rio Chopim</div>
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
-|-
-| Área de
-|-
-| contribuição
-|-
-| (km<math display="inline">^2</math>)
-|}
-|- style="border-top: 2px solid;"
+<div id='table-C'></div>
-| style="text-align: left;border-right: 2px solid;" |
+{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"
-{|  style="text-align: left; margin: 1em auto;min-width:50%;"
+|-style="text-align:center"
-|-
+!Sub-Bacia !! Comprimento do canal <br> <math>L</math> (m) !! Parâmetro <br> <math>\bar{\alpha }</math>  !! Parâmetro <br> <math>\bar{\beta }</math> !! Área de drenagem <br> (km<math>^2</math>) !! Área de drenagem acumulada <br>  (km<math>^2</math>)
-| Porto
+|- style="text-align:center"
-|-
+|    1
-| Palmeirinha
+|  26.830
+|  0,04113596
+|  1,59495580
+|  186,21
+| 186,21
+|- style="text-align:center"
+|    2
+|  47.500
+|  0,85500771
+|  0,87796837
+|  475,86
+| 662,07
+|- style="text-align:center"
+|    3
+|  30.170
+|  1,01710784
+|  0,86807752
+|  308,94
+| 971,01
+|- style="text-align:center"
+|    4
+|  26.580
+|  1,45162356
+|  0,83154649
+|  320,33
+| 1.291,34
+|- style="text-align:center"
+|    5
+|  41.670
+|  2,55791569
+|  0,76867062
+|  399,15
+| 1.690,49
+|- style="text-align:center"
+|    6
+|  13.000
+|  2,47863913
+|  0,77061510
+|  124,52
+| 1.815,01
+|- style="text-align:center"
+|    7
+|  65.000
+|  2,32144856
+|  0,77466167
+|  965,27
+| 2.780,28
+|- style="text-align:center"
+|    8
+|  22.000
+|  2,14755607
+|  0,77947050
+|  329,82
+| 3.110,10
+|- style="text-align:center"
+|    9
+|  8.180
+|  2,09342670
+|  0,78104723
+|  94,78
+| 3.204,88
 |}
-| style="border-left: 2px solid;border-right: 2px solid;" | 65927000
+<div class="center" style="font-size: 75%;">'''Fonte''': Kaviski e Gonçalves  <span id='citeF-49'></span>[[#cite-49|[49]]]</div>
-| style="border-left: 2px solid;border-right: 2px solid;" |
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
-|-
-| Coronel
-|-
-| Vivida
-|}
-| style="text-align: right;border-left: 2px solid;border-right: 2px solid;" | -26,03
-| style="text-align: right;border-left: 2px solid;border-right: 2px solid;" | -52,63
-| style="border-left: 2px solid;border-right: 2px solid;" | 501
-| style="border-left: 2px solid;" | 3.390
-|- style="border-top: 2px solid;border-bottom: 2px solid;"
-| style="text-align: left;border-right: 2px solid;" |
-{|  style="text-align: left; margin: 1em auto;min-width:50%;"
-|-
-| Salto
-|-
-| Claudelino
-|}
-| style="border-left: 2px solid;border-right: 2px solid;" | 65925000
+<div class="center" style="font-size: 75%;">'''Tabela D'''. Intensidade média de precipitação no rio Chopim, mês 07/1987 - pd e pe são, respectivamente, os planos do lado direito e do lado esquerdo no sentido do escoamento</div>
-| style="border-left: 2px solid;border-right: 2px solid;" | Clevelândia
-| style="text-align: right;border-left: 2px solid;border-right: 2px solid;" | -26,28
-| style="text-align: right;border-left: 2px solid;border-right: 2px solid;" | -52,30
-| style="border-left: 2px solid;border-right: 2px solid;" | 797
-| style="border-left: 2px solid;" | 1.660
-|}
+<div id='table-D'></div>
+{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"
-Fonte: Franco (2017) <span id='citeF-45'></span>[[#cite-45|[45]]].
+|-style="text-align:center"
+!rowspan='2' | Sub-Bacia !! rowspan='2' |Plano !! colspan='2' | 07/07/1987 <br> (m/h) <math>\cdot </math> 10<math>^{-3}</math> !! colspan='2' | 08/07/1987 <br>  (m/h) <math>\cdot </math> 10<math>^{-3}</math>
+|-style="text-align:center"
-{|  class="floating_tableSCP wikitable" style="text-align: right; margin: 1em auto;min-width:50%;"
+|  7h
-|+ style="font-size: 75%;" |<span id='table-C'></span>Tabla. C Caracterização física dos planos das 9 sub-bacias do rio Chopim - pd e pe são, respectivamente, os planos do lado direito e do lado esquerdo no sentido do escoamento.
+|  19h
-|- style="border-top: 2px solid;"
+| 7h
-| style="text-align: center;border-right: 2px solid;" |
+|  19h
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
+|- style="text-align:center"
-|-
+| rowspan='2' | 1
-| Sub-
+|  pd
-|-
+|  -
-| Bacia
+|  -
-|}
+| -
+| -
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | Plano
+|-style="text-align:center"
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" |
+|  pe
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
+|  0
-|-
+|  0
-| Compri-
+|  1,47400098
-|-
-| mento
-|-
-| <math>L</math> (m)
-|}
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" |
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
-|-
-| Largura
-|-
-| <math>b</math> (m)
-|}
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" |
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
-|-
-| Declividade
-|-
-| <math>S_0</math> (m/m)
-|}
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" |
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
-|-
-| Manning
-|-
-| <math>n_r</math> (m<math>^{-1/3}</math>s)
-|}
-| colspan='1' style="text-align: center;border-left: 2px solid;" |
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
-|-
-| Área
-|-
-| (km<math>^2</math>)
-|}
-|- style="border-top: 2px solid;"
-| rowspan='2' style="text-align: center;border-right: 2px solid;" | 1
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pd
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" | -
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" | -
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" | -
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" | -
-| colspan='1' style="text-align: center;border-left: 2px solid;" | -
-|-
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pe
-| style="border-left: 2px solid;border-right: 2px solid;" | 6.940
-| style="border-left: 2px solid;border-right: 2px solid;" | 26.830
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,0065
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,1200
-| style="border-left: 2px solid;" | 186,20
-|- style="border-top: 2px solid;"
-| rowspan='2' style="text-align: center;border-right: 2px solid;" | 2
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pd
-| style="border-left: 2px solid;border-right: 2px solid;" | 3.080
-| style="border-left: 2px solid;border-right: 2px solid;" | 47.500
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,0266
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,0840
-| style="border-left: 2px solid;" | 146,30
-|-
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pe
-| style="border-left: 2px solid;border-right: 2px solid;" | 6.940
-| style="border-left: 2px solid;border-right: 2px solid;" | 47.500
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,0065
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,1200
-| style="border-left: 2px solid;" | 329,65
-|- style="border-top: 2px solid;"
-| rowspan='2' style="text-align: center;border-right: 2px solid;" | 3
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pd
-| style="border-left: 2px solid;border-right: 2px solid;" | 3.300
-| style="border-left: 2px solid;border-right: 2px solid;" | 30.170
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,0146
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,1044
-| style="border-left: 2px solid;" | 99,56
-|-
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pe
-| style="border-left: 2px solid;border-right: 2px solid;" | 6.940
-| style="border-left: 2px solid;border-right: 2px solid;" | 30.170
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,0065
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,1200
-| style="border-left: 2px solid;" | 209,38
-|- style="border-top: 2px solid;"
-| rowspan='2' style="text-align: center;border-right: 2px solid;" | 4
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pd
-| style="border-left: 2px solid;border-right: 2px solid;" | 2.800
-| style="border-left: 2px solid;border-right: 2px solid;" | 26.580
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,0102
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,2832
-| style="border-left: 2px solid;" | 74,42
-|-
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pe
-| style="border-left: 2px solid;border-right: 2px solid;" | 9.250
-| style="border-left: 2px solid;border-right: 2px solid;" | 26.580
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,0129
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,2832
-| style="border-left: 2px solid;" | 245,87
-|- style="border-top: 2px solid;"
-| rowspan='2' style="text-align: center;border-right: 2px solid;" | 5
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pd
-| style="border-left: 2px solid;border-right: 2px solid;" | 2.800
-| style="border-left: 2px solid;border-right: 2px solid;" | 41.670
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,0102
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,2832
-| style="border-left: 2px solid;" | 116,68
-|-
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pe
-| style="border-left: 2px solid;border-right: 2px solid;" | 6.780
-| style="border-left: 2px solid;border-right: 2px solid;" | 41.670
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,0022
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,1020
-| style="border-left: 2px solid;" | 282,52
-|- style="border-top: 2px solid;"
-| rowspan='2' style="text-align: center;border-right: 2px solid;" | 6
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pd
-| style="border-left: 2px solid;border-right: 2px solid;" | 2.800
-| style="border-left: 2px solid;border-right: 2px solid;" | 13.000
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,0102
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,2832
-| style="border-left: 2px solid;" | 36,40
-|-
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pe
-| style="border-left: 2px solid;border-right: 2px solid;" | 6.780
-| style="border-left: 2px solid;border-right: 2px solid;" | 13.000
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,0022
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,1020
-| style="border-left: 2px solid;" | 88,14
-|- style="border-top: 2px solid;"
-| rowspan='2' style="text-align: center;border-right: 2px solid;" | 7
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pd
-| style="border-left: 2px solid;border-right: 2px solid;" | 5.180
-| style="border-left: 2px solid;border-right: 2px solid;" | 65.000
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,0170
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,0972
-| style="border-left: 2px solid;" | 336,70
-|-
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pe
-| style="border-left: 2px solid;border-right: 2px solid;" | 9.670
-| style="border-left: 2px solid;border-right: 2px solid;" | 65.000
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,0085
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,0972
-| style="border-left: 2px solid;" | 628,55
-|- style="border-top: 2px solid;"
-| rowspan='2' style="text-align: center;border-right: 2px solid;" | 8
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pd
-| style="border-left: 2px solid;border-right: 2px solid;" | 10.120
-| style="border-left: 2px solid;border-right: 2px solid;" | 22.000
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,0113
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,0984
-| style="border-left: 2px solid;" | 222,64
-|-
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pe
-| style="border-left: 2px solid;border-right: 2px solid;" | 4.870
-| style="border-left: 2px solid;border-right: 2px solid;" | 22.000
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,0129
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,1020
-| style="border-left: 2px solid;" | 107,14
-|- style="border-top: 2px solid;"
-| rowspan='2' style="text-align: center;border-bottom: 2px solid;border-right: 2px solid;" | 9
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pd
-| style="border-left: 2px solid;border-right: 2px solid;" | 6.720
-| style="border-left: 2px solid;border-right: 2px solid;" | 8.180
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,0127
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,1032
-| style="border-left: 2px solid;" | 54,97
-|- style="border-bottom: 2px solid;"
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pe
-| style="border-left: 2px solid;border-right: 2px solid;" | 4.870
-| style="border-left: 2px solid;border-right: 2px solid;" | 8.180
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,0129
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,1020
-| style="border-left: 2px solid;" | 39,84
-|}
-Fonte: Kaviski e Gonçalves (1996) <span id='citeF-49'></span>[[#cite-49|[49]]].
-{|  class="floating_tableSCP wikitable" style="text-align: right; margin: 1em auto;min-width:50%;"
-|+ style="font-size: 75%;" |<span id='table-D'></span>Tabla. D Caracterização física dos canais das 9 sub-bacias do rio Chopim
-|- style="border-top: 2px solid;"
-| style="text-align: center;border-right: 2px solid;" |
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
-|-
-| Sub-
-|-
-| Bacia
-|}
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" |
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
-|-
-| Comprimento
-|-
-| do canal
-|-
-| <math>L</math> (m)
-|}
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" |
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
-|-
-| Parâmetro
-|-
-| <math>\bar{\alpha }</math>
-|}
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" |
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
-|-
-| Parâmetro
-|-
-| <math>\bar{\beta }</math>
-|}
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" |
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
-|-
-| Área de
-|-
-| drenagem
-|-
-| (km<math>^2</math>)
-|}
-| colspan='1' style="text-align: center;border-left: 2px solid;" |
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
-|-
-| Área de dre-
-|-
-| nagem acumu-
-|-
-| lada (km<math>^2</math>)
-|}
-|- style="border-top: 2px solid;"
-| style="text-align: center;border-right: 2px solid;" |    1
-| style="border-left: 2px solid;border-right: 2px solid;" | 26.830
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,04113596
-| style="border-left: 2px solid;border-right: 2px solid;" | 1,59495580
-| style="border-left: 2px solid;border-right: 2px solid;" | 186,21
-| style="border-left: 2px solid;" | 186,21
-|- style="border-top: 2px solid;"
-| style="text-align: center;border-right: 2px solid;" |    2
-| style="border-left: 2px solid;border-right: 2px solid;" | 47.500
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,85500771
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,87796837
-| style="border-left: 2px solid;border-right: 2px solid;" | 475,86
-| style="border-left: 2px solid;" | 662,07
-|- style="border-top: 2px solid;"
-| style="text-align: center;border-right: 2px solid;" |    3
-| style="border-left: 2px solid;border-right: 2px solid;" | 30.170
-| style="border-left: 2px solid;border-right: 2px solid;" | 1,01710784
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,86807752
-| style="border-left: 2px solid;border-right: 2px solid;" | 308,94
-| style="border-left: 2px solid;" | 971,01
-|- style="border-top: 2px solid;"
-| style="text-align: center;border-right: 2px solid;" |    4
-| style="border-left: 2px solid;border-right: 2px solid;" | 26.580
-| style="border-left: 2px solid;border-right: 2px solid;" | 1,45162356
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,83154649
-| style="border-left: 2px solid;border-right: 2px solid;" | 320,33
-| style="border-left: 2px solid;" | 1.291,34
-|- style="border-top: 2px solid;"
-| style="text-align: center;border-right: 2px solid;" |    5
-| style="border-left: 2px solid;border-right: 2px solid;" | 41.670
-| style="border-left: 2px solid;border-right: 2px solid;" | 2,55791569
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,76867062
-| style="border-left: 2px solid;border-right: 2px solid;" | 399,15
-| style="border-left: 2px solid;" | 1.690,49
-|- style="border-top: 2px solid;"
-| style="text-align: center;border-right: 2px solid;" |    6
-| style="border-left: 2px solid;border-right: 2px solid;" | 13.000
-| style="border-left: 2px solid;border-right: 2px solid;" | 2,47863913
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,77061510
-| style="border-left: 2px solid;border-right: 2px solid;" | 124,52
-| style="border-left: 2px solid;" | 1.815,01
-|- style="border-top: 2px solid;"
-| style="text-align: center;border-right: 2px solid;" |    7
-| style="border-left: 2px solid;border-right: 2px solid;" | 65.000
-| style="border-left: 2px solid;border-right: 2px solid;" | 2,32144856
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,77466167
-| style="border-left: 2px solid;border-right: 2px solid;" | 965,27
-| style="border-left: 2px solid;" | 2.780,28
-|- style="border-top: 2px solid;"
-| style="text-align: center;border-right: 2px solid;" |    8
-| style="border-left: 2px solid;border-right: 2px solid;" | 22.000
-| style="border-left: 2px solid;border-right: 2px solid;" | 2,14755607
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,77947050
-| style="border-left: 2px solid;border-right: 2px solid;" | 329,82
-| style="border-left: 2px solid;" | 3.110,10
-|- style="border-top: 2px solid;border-bottom: 2px solid;"
-| style="text-align: center;border-right: 2px solid;" |    9
-| style="border-left: 2px solid;border-right: 2px solid;" | 8.180
-| style="border-left: 2px solid;border-right: 2px solid;" | 2,09342670
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,78104723
-| style="border-left: 2px solid;border-right: 2px solid;" | 94,78
-| style="border-left: 2px solid;" | 3.204,88
-|}
-Fonte: Kaviski e Gonçalves (1996) <span id='citeF-49'></span>[[#cite-49|[49]]].
-{|  class="floating_tableSCP wikitable" style="text-align: right; margin: 1em auto;min-width:50%;"
-|+ style="font-size: 75%;" |<span id='table-E'></span>Tabla. E Intensidade média de precipitação no rio Chopim, mês 07/1987 - pd e pe são, respectivamente, os planos do lado direito e do lado esquerdo no sentido do escoamento.
-|- style="border-top: 2px solid;"
-| rowspan='2' style="text-align: center;border-right: 2px solid;" |
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
-|-
-| Sub-
-|-
-| bacia
-|}
-| rowspan='2' style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | Plano
-| colspan='2' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" |
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
-|-
-| 07/07/1987
-|-
-| (m/h) <math>\cdot </math> 10<math>^{-3}</math>
-|}
-| colspan='2' style="text-align: center;border-left: 2px solid;" |
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
-|-
-| 08/07/1987
-|-
-| (m/h) <math>\cdot </math> 10<math>^{-3}</math>
-|}
-|-
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" | 7h
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" | 19h
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;" | 7h
-| colspan='1' style="text-align: center;" | 19h
-|- style="border-top: 2px solid;"
-| rowspan='2' style="text-align: center;border-right: 2px solid;" | 1
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pd
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" | -
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" | -
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;" | -
-| colspan='1' style="text-align: center;" | -
-|-
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pe
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
-| colspan='1' style="border-right: 2px solid;border-left: 2px solid;" | 1,47400098
 | 0,00174585
-|- style="border-top: 2px solid;"
+|- style="text-align:center"
-| rowspan='2' style="text-align: center;border-right: 2px solid;" | 2
+| rowspan='2' | 2
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pd
+|  pd
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,65031194
+|  0,65031194
-| colspan='1' style="border-right: 2px solid;border-left: 2px solid;" | 2,24260054
+|  2,24260054
 | 0,00242411
-|-
+|-style="text-align:center"
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pe
+|  pe
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| colspan='1' style="border-right: 2px solid;border-left: 2px solid;" | 1,47400098
+|  1,47400098
 | 0,00174585
-|- style="border-top: 2px solid;"
+|- style="text-align:center"
-| rowspan='2' style="text-align: center;border-right: 2px solid;" | 3
+| rowspan='2' | 3
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pd
+|  pd
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| colspan='1' style="border-right: 2px solid;border-left: 2px solid;" | 1,67223974
+|  1,67223974
 | 0
-|-
+|-style="text-align:center"
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pe
+|  pe
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| colspan='1' style="border-right: 2px solid;border-left: 2px solid;" | 1,47400098
+|  1,47400098
 | 0,00174585
-|- style="border-top: 2px solid;"
+|- style="text-align:center"
-| rowspan='2' style="text-align: center;border-right: 2px solid;" | 4
+| rowspan='2' | 4
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pd
+|  pd
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| colspan='1' style="border-right: 2px solid;border-left: 2px solid;" | 0,87767508
+|  0,87767508
 | 0
-|-
+|-style="text-align:center"
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pe
+|  pe
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| colspan='1' style="border-right: 2px solid;border-left: 2px solid;" | 1,15756772
+|  1,15756772
 | 0
-|- style="border-top: 2px solid;"
+|- style="text-align:center"
-| rowspan='2' style="text-align: center;border-right: 2px solid;" | 5
+| rowspan='2' | 5
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pd
+|  pd
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| colspan='1' style="border-right: 2px solid;border-left: 2px solid;" | 0,87767508
+|  0,87767508
 | 0
-|-
+|-style="text-align:center"
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pe
+|  pe
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| colspan='1' style="border-right: 2px solid;border-left: 2px solid;" | 0,99927618
+|  0,99927618
 | 0
-|- style="border-top: 2px solid;"
+|- style="text-align:center"
-| rowspan='2' style="text-align: center;border-right: 2px solid;" | 6
+| rowspan='2' | 6
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pd
+|  pd
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| colspan='1' style="border-right: 2px solid;border-left: 2px solid;" | 0,87767508
+|  0,87767508
 | 0
-|-
+|-style="text-align:center"
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pe
+|  pe
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| colspan='1' style="border-right: 2px solid;border-left: 2px solid;" | 0,99927618
+|  0,99927618
 | 0
-|- style="border-top: 2px solid;"
+|- style="text-align:center"
-| rowspan='2' style="text-align: center;border-right: 2px solid;" | 7
+| rowspan='2' | 7
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pd
+|  pd
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| colspan='1' style="border-right: 2px solid;border-left: 2px solid;" | 0,90541231
+|  0,90541231
 | 0
-|-
+|-style="text-align:center"
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pe
+|  pe
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| colspan='1' style="border-right: 2px solid;border-left: 2px solid;" | 1,15131540
+|  1,15131540
 | 0
-|- style="border-top: 2px solid;"
+|- style="text-align:center"
-| rowspan='2' style="text-align: center;border-right: 2px solid;" | 8
+| rowspan='2' | 8
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pd
+|  pd
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| colspan='1' style="border-right: 2px solid;border-left: 2px solid;" | 1,08201453
+|  1,08201453
 | 0
-|-
+|-style="text-align:center"
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pe
+|  pe
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| colspan='1' style="border-right: 2px solid;border-left: 2px solid;" | 1,26797648
+|  1,26797648
 | 0,00923398
-|- style="border-top: 2px solid;"
+|- style="text-align:center"
-| rowspan='2' style="text-align: center;border-bottom: 2px solid;border-right: 2px solid;" | 9
+| rowspan='2' | 9
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pd
+|  pd
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| style="border-left: 2px solid;border-right: 2px solid;" | 0,72930031
+|  0,72930031
-| colspan='1' style="border-right: 2px solid;border-left: 2px solid;" | 1,47131551
+|  1,47131551
 | 0,02406490
-|- style="border-bottom: 2px solid;"
+|- style="text-align:center"
-| style="text-align: center;border-left: 2px solid;border-right: 2px solid;" | pe
+|  pe
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
+|  0
-| colspan='1' style="border-right: 2px solid;border-left: 2px solid;" | 1,26797648
+|  1,26797648
 | 0,00923398
 |}
-Fonte: Kaviski e Gonçalves (1996) <span id='citeF-49'></span>[[#cite-49|[49]]].
+<div class="center" style="font-size: 75%;">'''Fonte''': Kaviski e Gonçalves  <span id='citeF-49'></span>[[#cite-49|[49]]]</div>
-{|  class="floating_tableSCP wikitable" style="text-align: right; margin: 1em auto;min-width:50%;"
+<div class="center" style="font-size: 75%;">'''Tabela E'''. Vazões iniciais e coeficientes de interpolação da vazão nos trechos de canais do rio Chopim, mês 07/1987</div>
-|+ style="font-size: 75%;" |<span id='table-F'></span>Tabla. F Vazões iniciais e coeficientes de interpolação da vazão nos trechos de canais do rio Chopim, mês 07/1987
-|- style="border-top: 2px solid;"
-| rowspan='2' style="text-align: center;border-right: 2px solid;" | Canal
-| colspan='2' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" |
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
-|-
-| Vazão inicial na
-|-
-| seção tranversal
-|-
-| (m<math>^3</math>/s)
-|}
-| colspan='2' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" |
+<div id='table-E'></div>
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
+{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"
-|-
+|-style="text-align:center"
-| Coeficientes da
+! rowspan='2' | Canal !!  colspan='2' | Vazão inicial na  seção <br>  tranversal  (m<math>^3</math>/s) !! colspan='2' | Coeficientes da interpolação  <br>  linear <math>Q(x)=ax+b</math> !! rowspan='2' style="vertical-align:top;"| Vazão inicial como contribuição  <br> lateral (m<math>^2</math>/s) <math>\cdot </math> 10<math>^{-4}</math>
-|-
+|-style="text-align:center"
-| interpolação linear
+| Montante
-|-
+|  Jusante
-| <math>Q(x)=ax+b</math>
+|  <math>a\, \cdot 10^{-4}</math>
+|  <math>b</math>
+|- style="text-align:center"
+|  1
+| 0
+| 4,91134977
+| 1,83054408
+| 0
+| 1,83054408
+|- style="text-align:center"
+|    2
+| 4,91134977
+| 17,46209717
+| 2,64226261
+| 4,91134977
+| 2,64226261
+|- style="text-align:center"
+|    3
+| 17,46209717
+| 25,61030197
+| 2,70076394
+| 17,46209717
+| 2,70076394
+|- style="text-align:center"
+|    4
+| 25,61030197
+| 34,05889511
+| 3,17855272
+| 25,61030197
+| 3,17855272
+|- style="text-align:center"
+|    5
+| 34,05889511
+| 44,58636093
+| 2,52638969
+| 34,05889510
+| 2,52638969
+|- style="text-align:center"
+|    6
+| 44,58636093
+| 48,07283783
+| 2,68190530
+| 44,58636093
+| 2,68190530
+|- style="text-align:center"
+|    7
+| 48,07283783
+| 75,09886169
+| 4,15784983
+| 48,07283783
+| 4,15784983
+|- style="text-align:center"
+|    8
+| 75,09886169
+| 84,33335876
+| 4,19749867
+| 75,09886168
+| 4,19749867
+|- style="text-align:center"
+|    9
+| 84,33335876
+| 86,98714447
+| 3,24423681
+| 84,33335877
+| 3,24423681
 |}
-| rowspan='2' colspan='1' style="text-align: center;border-left: 2px solid;" |
+<div class="center" style="font-size: 75%;">'''Fonte''': Kaviski e Gonçalves  <span id='citeF-49'></span>[[#cite-49|[49]]]</div>
-{|  style="text-align: center; margin: 1em auto;min-width:50%;"
-|-
-| Vazão inicial como
-|-
-| contribuição lateral
-|-
-| (m<math>^2</math>/s) <math>\cdot </math> 10<math>^{-4}</math>
-|}
-|-
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" | Montante
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" | Jusante
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" | <math>a\, \cdot 10^{-4}</math>
-| colspan='1' style="text-align: center;border-right: 2px solid;border-left: 2px solid;border-right: 2px solid;" | <math>b</math>
-|- style="border-top: 2px solid;"
-| style="text-align: center;border-right: 2px solid;" |    1
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
-| style="border-left: 2px solid;border-right: 2px solid;" | 4,91134977
-| style="border-left: 2px solid;border-right: 2px solid;" | 1,83054408
-| style="border-left: 2px solid;border-right: 2px solid;" | 0
-| style="border-left: 2px solid;" | 1,83054408
-|- style="border-top: 2px solid;"
-| style="text-align: center;border-right: 2px solid;" |    2
-| style="border-left: 2px solid;border-right: 2px solid;" | 4,91134977
-| style="border-left: 2px solid;border-right: 2px solid;" | 17,46209717
-| style="border-left: 2px solid;border-right: 2px solid;" | 2,64226261
-| style="border-left: 2px solid;border-right: 2px solid;" | 4,91134977
-| style="border-left: 2px solid;" | 2,64226261
-|- style="border-top: 2px solid;"
-| style="text-align: center;border-right: 2px solid;" |    3
-| style="border-left: 2px solid;border-right: 2px solid;" | 17,46209717
-| style="border-left: 2px solid;border-right: 2px solid;" | 25,61030197
-| style="border-left: 2px solid;border-right: 2px solid;" | 2,70076394
-| style="border-left: 2px solid;border-right: 2px solid;" | 17,46209717
-| style="border-left: 2px solid;" | 2,70076394
-|- style="border-top: 2px solid;"
-| style="text-align: center;border-right: 2px solid;" |    4
-| style="border-left: 2px solid;border-right: 2px solid;" | 25,61030197
-| style="border-left: 2px solid;border-right: 2px solid;" | 34,05889511
-| style="border-left: 2px solid;border-right: 2px solid;" | 3,17855272
-| style="border-left: 2px solid;border-right: 2px solid;" | 25,61030197
-| style="border-left: 2px solid;" | 3,17855272
-|- style="border-top: 2px solid;"
-| style="text-align: center;border-right: 2px solid;" |    5
-| style="border-left: 2px solid;border-right: 2px solid;" | 34,05889511
-| style="border-left: 2px solid;border-right: 2px solid;" | 44,58636093
-| style="border-left: 2px solid;border-right: 2px solid;" | 2,52638969
-| style="border-left: 2px solid;border-right: 2px solid;" | 34,05889510
-| style="border-left: 2px solid;" | 2,52638969
-|- style="border-top: 2px solid;"
-| style="text-align: center;border-right: 2px solid;" |    6
-| style="border-left: 2px solid;border-right: 2px solid;" | 44,58636093
-| style="border-left: 2px solid;border-right: 2px solid;" | 48,07283783
-| style="border-left: 2px solid;border-right: 2px solid;" | 2,68190530
-| style="border-left: 2px solid;border-right: 2px solid;" | 44,58636093
-| style="border-left: 2px solid;" | 2,68190530
-|- style="border-top: 2px solid;"
-| style="text-align: center;border-right: 2px solid;" |    7
-| style="border-left: 2px solid;border-right: 2px solid;" | 48,07283783
-| style="border-left: 2px solid;border-right: 2px solid;" | 75,09886169
-| style="border-left: 2px solid;border-right: 2px solid;" | 4,15784983
-| style="border-left: 2px solid;border-right: 2px solid;" | 48,07283783
-| style="border-left: 2px solid;" | 4,15784983
-|- style="border-top: 2px solid;"
-| style="text-align: center;border-right: 2px solid;" |    8
-| style="border-left: 2px solid;border-right: 2px solid;" | 75,09886169
-| style="border-left: 2px solid;border-right: 2px solid;" | 84,33335876
-| style="border-left: 2px solid;border-right: 2px solid;" | 4,19749867
-| style="border-left: 2px solid;border-right: 2px solid;" | 75,09886168
-| style="border-left: 2px solid;" | 4,19749867
-|- style="border-top: 2px solid;border-bottom: 2px solid;"
-| style="text-align: center;border-right: 2px solid;" |    9
-| style="border-left: 2px solid;border-right: 2px solid;" | 84,33335876
-| style="border-left: 2px solid;border-right: 2px solid;" | 86,98714447
-| style="border-left: 2px solid;border-right: 2px solid;" | 3,24423681
-| style="border-left: 2px solid;border-right: 2px solid;" | 84,33335877
-| style="border-left: 2px solid;" | 3,24423681
-|}
+<div class="center" style="font-size: 75%;">'''Tabela F'''. Teste ''t'' de Student para as sub-bacias no rio Chopim</div>
-{|  class="floating_tableSCP wikitable" style="text-align: right; margin: 1em auto;min-width:50%;"
+<div id='table-F'></div>
-|+ style="font-size: 75%;" | Vazões iniciais
+{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"
-|}
+|-style="text-align:center"
+!rowspan='2'  | Quantidades estatísticas !! colspan='2' | Sub-bacia 5<br> Salto Claudelino !! colspan='2' | Sub-Bacia 9<br>  Porto Palmeirinha
+|-style="text-align:center"
+|  Observado
+|  Simulado
+|  Observado
+|  Simulado
+|- style="text-align:center"
+|  style="text-align:left" |  Média
+| 80,04
+| 76,58
+| 142,12
+| 139,02
+|-style="text-align:center"
+| style="text-align:left" | Variância
+|  1.493,98
+|  3.021,13
+|  3.376,02
+| 4.864,39
+|-style="text-align:center"
+|  style="text-align:left" |Número de observações
+| colspan='2'  | 21
+| colspan='2' | 21
+|-style="text-align:center"
+|  style="text-align:left" |Graus de liberdade
+| colspan='2' | 40
+| colspan='2'  | 40
+|-style="text-align:center"
+| style="text-align:left" | Estatística do teste
+| colspan='2' | 0,2361
+| colspan='2' | 0,1567
+|-style="text-align:center"
+| style="text-align:left" |''t'' bi-caudal
+| colspan='2'  | 2,0211
+| colspan='2' | 2,0211
+|- style="text-align:center"
+| style="text-align:left" |Nível de significância
+| colspan='2'  | <math>5%</math>
+| colspan='2' | <math>5%</math>
-Fonte:Kaviski e Gonçalves (1996) <span id='citeF-49'></span>[[#cite-49|[49]]].
+|}