Hybrid numerical scheme for simulation of transient flow in a pluvial water reservoir

Latest revision as of 16:22, 23 February 2021

Resumo

Este trabalho apresenta a simulação hidráulica do escoamento em um reservatório, utilizando o método dos volumes finitos para resolver o sistema de equações que modelam o fluxo bidimensional em águas rasas, negligenciando as tensões tangenciais. Para resolver as equações, utilizou-se um esquema híbrido de volumes finitos. A solução do sistema de equações linearizadas foi obtida por implementação computacional do método iterativo de Gauss-Seidel. Para ilustrar a aplicação do esquema numérico em engenharia hidráulica, apresenta-se o estudo do escoamento em um reservatório subterrâneo com pilares internos, cujo fluxo é gerado pela abertura de dois portões de saída. O código desenvolvido permite a análise temporal do volume e do fluxo no reservatório, a interpretação gráfica do fenômeno físico investigado bem como o cálculo da precisão do modelo. Os resultados obtidos mostram a interpretação física esperada, boa concordância com os dados da literatura, boa precisão, estabilidade e baixo custo computacional.

Palavras chave: Equações de águas rasas, volumes finitos, interpolação híbrida, estabilidade numérica, simulação de leito seco, reservatório de contenção

Abstract

This work presents the numeric simulation of flow in a reservoir, using the finite volume method for solver the system of equations that model the two-dimensional flow in shallow water, neglecting the tangential tensions. The solution of the system of linearized equations was obtained by computational implementation of the Gauss-Seidel iterative method. To illustrate the application of the numerical scheme in hydraulic engineering, it was considered the study of flow in an underground reservoir with internal pillars, whose flow is generated by the opening of two outlet gates. The employed code allowed the temporal analysis of the volume and the flow in the reservoir, the graphical interpretation of the investigated physical phenomenon as well as, the calculation of the precision of the model. The results obtained show the expected physical interpretation, good agreement with literature data, good precision,stability and low computational cost.

Keywords: Shallow water equations, finite volumes, hybrid interpolation, numerical stability, dry bed simulation, containment reservoir

1. Introdução

A ocorrência ou agravamento de inundações nas grandes cidades, bem como a quebra de barragens, tornaram-se problemas cada vez mais recorrentes, causando mortes e fortes impactos ambientais. Nesse cenário, fica evidente a necessidade de conhecimento prévio do comportamento dos fluxos de superfície livre e simulações de abertura de comportas de barragens para manejo de recursos hídricos, avaliação de riscos e usos para obtenção de melhor disponibilidade visando a preservação ambiental. Para tanto, é necessário utilizar metodologias que possam descrever as variáveis hidrodinâmicas envolvidas, como campo de velocidade do fluxo, turbulência e mudanças temporais na altura da superfície da água [1,2]. Em [3] sugere-se o uso de modelos matemáticos para entender os padrões de fluxo em corpos de água. Modelos matemáticos adequados são ferramentas úteis no gerenciamento de recursos hídricos, pois reduzem o tempo de análise, custos e riscos envolvidos em levar dados experimentais em fluxos reais auxiliando na tomada de decisões e na definição de ações corretivas diante de vazamentos iminentes. Matematicamente, os fluxos de superfície livre podem ser descritos de várias maneiras. Uma abordagem mais simples e muito aplicada é dada pelas Equações de Águas Rasas, que representam uma simplificação das Equações de Navier-Stokes e é aplicável aos fluxos nos quais as tensões e acelerações verticais podem ser negligenciadas sem prejuízo do resultado. Essas equações são derivadas por meio da integração das equações governantes (Navier-Stokes e Continuidade) na direção vertical. Sua derivação está fora do escopo deste trabalho, mas pode ser vista em detalhes em [4] . Esse sistema de equações, embora definido em um domínio bidimensional, fornece a altura da coluna de água e as componentes de velocidade em um fluxo instável. No entanto, a obtenção de uma solução analítica para o sistema de equações de águas rasas é, na maioria dos casos, uma tarefa difícil, pois trata-se de um conjunto de equações diferenciais parciais hiperbólicas não lineares com termos fontes. O que requer um tratamento numérico robusto para evitar instabilidades numéricas. Ao longo das últimas três décadas, muitos avanços foram realizados na solução numérica das equações de águas rasas para simulação de escoamentos unidimensionais aplicados a hidráulica fluvial em rios, canais e tubulações [5,6,7,8,9], e modelos bidimensionais aplicados à análise de escoamentos causados pela abertura instantânea de comportas e ruptura de barragens [10,11,12,13,14,15]. Em [10] utilizaram esquemas explícitos de diferenças finitas, na discretização das equações de águas rasas bidimensionais, adicionando condições de estabilidade e viscosidade artificial para suavizar as oscilações de alta frequência. Para validar, os esquemas discretos foram aplicados em um problema típico de engenharia hidráulica: abertura parcial de uma barragem e passagem de uma onda de inundação através de uma contração do canal, adotando a razão entre o nível de água a jusante $(h_{r})$ e à montante $(h_{l})$ como sendo, $h_{r}/h_{l}=0.5$ . Os resultados dos diferentes esquemas foram comparados e a ocorrência de oscilações numéricas foi observada, evidenciando a necessidade de se utilizar esquemas implícitos. O estudo de caso apontado em [10] foi analisado por [15] usando um esquema de elementos finitos de alta ordem e uma melhoria na relação $h_{r}/h_{l}$ para $0,05$ foi obtida, buscando simular a condição de leito seco em uma determinada região do domínio. Em [12] usaram dados experimentais em um evento de quebra de barragem para verificar a precisão, estabilidade e confiabilidade de uma solução numérica das equações de águas rasas obtidas por um esquema de volumes finitos bidimensional. O modelo numérico permitiu uma precisão de segunda ordem, tanto no espaço como no tempo, acarretando porém um alto custo computacional, apresentando algumas discrepâncias entre os resultados reais dos medidores e os resultados numéricos nos pontos pesquisados. Os autores destacam a necessidade de simular a propagação de uma onda de inundação para um fundo seco e áspero e otimizar o algoritmo para reduzir o tempo de computação. Um esquema de volumes finitos não estruturado para malhas triangulares para simular fluxos bidimensionais instáveis em águas rasas com frentes molhadas/secas sobre topografia complexa é apresentado em [11]. Utiliza-se em [11] para o esquema temporal o método explícito de Runge-Kutta e para o espacial a malha não estruturada. Comparações satisfatórias foram obtidas com dados mensurados e numéricos. Todavia, os autores sugerem a incorporação de um procedimento de escalonamento num esquema implícito temporal, e uma implementação paralela do algoritmo, a fim de reduzir o tempo computacional. A eficiência de um esquema implícito de volumes finitos é explorada em [13] para simulação de fluidos em águas rasas bidimensionais, com frentes secas/molhadas. Adotou-se em [13] um tratamento complexo para a discretização espacial, fazendo uso de uma malha não estruturada flexível para obter um melhor ajuste à geometrias irregulares. No entanto, é feito ressalvas quanto ao tempo computacional para obtenção da solução por esquemas implícitos, pois embora um método implícito normalmente exija menos etapas da solução, cada uma exige mais tempo computacional do que as soluções obtidas por um esquema explícito. Os autores apontam a necessidade de se utilizar um solucionador de matrizes que com um número reduzido de iterações alcance a convergência. Neste contexto, o objetivo principal deste trabalho é formular um esquema híbrido e totalmente implícito de volumes finitos para discretização das equações de águas rasas, a fim de obter uma solução numérica estável com um algoritmo iterativo otimizado, com convergência rápida e baixo custo computacional, para aplicação na simulação de fluxos bidimensionais transitórios em topografias complexas (com declives e sujeitos a resistência ao fluxo). Condições iniciais descontínuas representando leitos úmidos/secos foram adotadas. Para tanto, o trabalho está organizado da seguinte forma: A Seção 2 apresenta as equações hiperbólicas, que regem a natureza das águas rasas, considerando as perdas de atrito e declives da topografia do domínio computacional e detalhando a configuração numérica destas equações; na Seção 3, o modelo é validado e uma comparação é feita com os resultados obtidos com o trabalho apresentado em [10]. A robustez da solução numérica é verificada com uma extrapolação na relação $h_{r}/h_{l}$ para $0.00002$ , na Seção 4; o esquema é aplicado na análise de um problema hidráulico real na Seção 5 e; as conclusões são apresentadas na Seção 6.

2. Modelo matemático e solução numérica

Na análise do escoamento bidimensional transiente em um canal, cuja topografia do fundo é descrita por uma superfície $z_{0}(x,y)$ , busca-se determinar o nível do fluido de superfície livre $h(x,y,t)$ e o campo de velocidade do fluxo $U(u,v)$ para cada instante de tempo $t$ , conforme Figura 1.

Figura 1. Notação adotada na análise de fluxo em águas rasas

2.1 Sistema de Equações Governantes

O sistema de equações diferenciais que descreve o escoamento hidrostático de um fluido incompressível em águas rasas com superfície superior livre é obtido pela integração ao longo da profundidade das equações tridimensionais que modelam a quantidade de movimento e conservação de massa num escoamento, assumindo distribuição de velocidade uniforme na direção vertical. Como em Mahmood and Yevjevich [3], adota-se as seguintes simplificações: (i) a distribuição de pressão na vertical é puramente hidrostática; (ii) não há contribuições laterais ao escoamento no reservatório analisado; (iii) a massa específica e a viscosidade do fluido no reservatório não sofrem variações significativas ao longo do tempo e do espaço; (iv) as tensões tangenciais são desprezíveis. Deste modo, as equações governantes do modelo hidrodinâmico, também ditas equações de águas rasas ou equações bidimensionais de Saint Venant, podem ser escritas na forma conservativa (com os fluxos dentro do sinal da derivada):

(1)

(2)

(3)

onde $x$ e $y$ representam as componentes cartesianas nas direções longitudinal e transversal, respectivamente; $h$ é profundidade da água no reservatório; $u$ e $v$ são as componentes cartesianas da velocidade; $g$ é a aceleração da gravidade, $t$ é o tempo, e $S_{y}$ e $S_{x}$ são os declives do fundo do reservatório definidos como:

(4)

enquanto os declives de resistência ao escoamento (perda de energia), nas direções $x$ e $y$ , respectivamente, são dadas por:

(5)

onde $c_{M}(m/s^{\frac {1}{3}})$ é o coeficiente de rugosidade de Manning, que quantifica o atrito do fluxo com o fundo do canal. As equações de águas rasas (1) - (3) não possuem solução analítica fechada, portanto, uma abordagem numérica é empregada para resolvê-las de forma discreta e linearizada, onde o domínio é dividido em um número finito de células resultando em um sistema linear algébrico que é resolvido por um algoritmo computacional apropriado codificado em FORTRAN 90. Neste estudo, diferentes condições de contorno e topografia de fundo foram impostas para verificar a robustez e aplicabilidade do código. Na etapa inicial, o nível da água é sempre conhecido e todas as velocidades são iguais a zero para todo o domínio. Essas condições estão detalhadas nas Subseções 3.1 e 3.2.

2.2 Discretização do domínio e linearização das equações governantes

A solução numérica do modelo hidrodinâmico descrito pelas equações (1)-(3) é obtida por meio do método dos volumes finitos (MVF). Os fundamentos matemáticos e físicos do MVF são bem discutidos por [16,17,18,19,20]. De acordo com [20], o MVF consiste na integração espacial e temporal das equações diferenciais na forma conservativa (com todos os fluxos envolvidos pelo operador derivada) em um volume genérico de controle. O sistema de equações (1)-(3) é integrado espacialmente (sobre um volume de controle retangular) e no tempo (de $t$ a $t+\Delta t$ ) e é empregado para a integração temporal a seguinte aproximação numérica:

(6)

Adota-se aqui uma notação mais compacta:

(7)

E assim, a Eq. (6) pode ser reescrita como:

(8)

Neste artigo, adota-se a formulação totalmente implícita, $(\theta =1)$ , ou seja, todas as derivadas são avaliadas no nível de tempo $t+\Delta \,t$ . Com uma implementação mais complexa, a formulação totalmente implícita é limitada apenas pela precisão e não sofre instabilidades numéricas devido à mudança de coeficientes de sinal pela variação do intervalo de tempo. E assim, permite o uso de etapas de tempo maiores com menor tempo de computação. Esse é um método incondicionalmente estável no tempo. Obviamente, esta estabilidade pode ser alterada devido às não linearidades do acoplamento e ao comportamento dos termos de fontes, como a perda de energia devido à declividade, o que pode gerar oscilações numéricas. Portanto, faz-se necessário um tratamento cuidadoso na obtenção de uma solução numérica estável e também precisa. A derivada temporal é avaliada com uma aproximação de Euler regressiva de $1^{a}$ ordem:

(9)

onde $\phi$ representa uma variável. Na integração espacial é utilizado o Teorema Fundamental do Cálculo para reescrever as integrais das equações (1)-(3) como um sistema de equações constituído pelos valores das variáveis avaliadas nos pontos de integração $w,e,s,n$ , que por sua vez, são avaliados por funções de interpolação. Para o cálculo destes valores, adota-se uma combinação dos esquemas de diferenciação central de $2^{a}$ ordem (CDS - Central Differencing Scheme) e um esquema upwind (UDS - Upwind Differencing Scheme) de $1^{a}$ ordem. Na interpolação central, os valores das variáveis nas faces do volume de controle são dados pelo valor médio dos centróides vizinhos, isto é,

(10)

No esquema upwind, os valores das variáveis nas faces do volume de controle são dados considerando primeiramente a direção do fluxo que transporta alguma propriedade do centróide vizinho para a face onde o fluxo é avaliado, por exemplo,

(11)

O sistema (1) - (3) apresenta produtos de variáveis envolvendo $u$ , $v$ e $h$ . Esta característica é a natureza não-linear das equações de Saint-Venant e uma das razões pelas quais não existe uma solução analítica para problemas complexos de águas rasas, e assim, o sistema deve ser linearizado para permitir que o método numérico gere um sistema linear resolvido por qualquer solucionador linear. A linearização é feita considerando no produto de variáveis apenas uma delas como a variável ativa, e as outras são constantes ou variáveis secundárias calculadas a partir da condição inicial da iteração de looping do tempo anterior e/ou do último passo de tempo, por exemplo,

(12)

onde neste exemplo, o índice $0$ representa os valores disponíveis do cálculo anterior e esses valores são usados para calcular o coeficiente da matriz para a variável ativa $h$ . Naturalmente, outras equações estão disponíveis para resolver $u$ e $v$ como a variável ativa (ou variável principal), e estas são empregadas para calcular cada coeficiente, assim, quando uma nova solução é alcançada, um novo cálculo é iniciado e novos coeficientes são definidos, e uma nova solução é executada até a convergência. Na combinação CDS/UDS, o esquema UDS é escolhido para a variável principal e a aproximação CDS para a variável secundária conhecida, por exemplo na integral da Eq. (1) o UDS é usado para a variável principal $h$ , e para $u$ e $v$ é aplicada a interpolação CDS usando os valores nodais vizinhos conhecidos $u^{0},\,v^{0}$ do nó $P$ . Todos os valores nas faces são funções dos valores nodais $P,W,E,S$ , usando coeficientes $\alpha ^{,}s$ que representam o sinal das componentes $u$ e $v$ nas faces do volume de controle, isto é, a direção do vetor de velocidade:

(13)

(14)

(15)

(16)

onde,

(17)

(18)

Portanto, a partir da formulação totalmente implícita e da interpolação híbrida CDS/UDS, a equação algébrica discreta para a Eq. (1) é dada por:

(19)

onde,

(20)

(21)

(22)

(23)

(24)

(25)

Para integrar a Eq. (2) aplica-se o Teorema Fundamental do Cálculo considerando a formulação totalmente implícita. Deste modo, o lado esquerdo desta equação é reescrita em função das variáveis calculadas em $P$ e nas faces $w,e,s,n$ , conforme a expressão $($ 26 $)$ :

\left[(uh)_{P}-(uh)_{P}^{0}\right]{\frac {\Delta x\Delta y}{\Delta t}}+\left(\Delta y\left[\left(u^{2}h+{\frac {1}{2}}gh^{2}\right)_{e}-\left(u^{2}h+{\frac {1}{2}}gh^{2}\right)_{w}\right]+\Delta x\left[(uvh)_{n}-(uvh)_{s}\right]\right)_{t+\Delta t}.

(26)

Na expressão $($ 26 $)$ o termo $u^{2}$ é escrito como $(u^{0}\,u)$ onde $u^{0}$ é conhecido e $u$ é a variável ativa a ser calculada. A interpolação UDS é usada para avaliar $u$ , enquanto a interpolação CDS para avaliar $u^{0}$ nas faces. Por exemplo,

(27)

com $\alpha _{e}$ definido na equação (17). O termo quadrático da expressão (26) da variável secundária $h$ é calculado usando a interpolação CDS com os valores disponíveis da Eq. (19) e assim,

(28)

A mesma abordagem é usada nos termos advectivos da expressão (26). O produto $(uvh)$ é dividido em $(hv^{0})u$ com UDS para variável ativa $u$ e CDS para as variáveis secundárias, $h$ atualizada e $v^{0}$ . Para a integração do lado direito da Eq. (2), que representa a resistência à declividade, o valor de $h$ é calculada no centróide $P$ , pois fisicamente a fonte representa a média integral no volume discreto centrado em $P$ , isto é, o fonte não flui, é um termo volumétrico. Portanto, $h=h_{P}$ e a magnitude da velocidade é constante. Dessa forma, a integral do lado direito da Eq. (2) é reescrita da seguinte forma:

\left[gh_{P}(-1)(z_{e}-z_{w})\Delta y+(-1)\left({\frac {gc_{M}^{2}}{h_{P}^{1/3}}}{\sqrt {(u_{P}^{0})^{2}+(v_{P}^{0})^{2}}}\right){\frac {u_{e}+u_{w}}{2}}\Delta x\Delta y\right]_{t+\Delta t},

(29)

onde os valores da variável principal $u$ nas faces são obtidos por UDS. Assim, a forma discreta linearizada da Eq. (2) é dada por:

(30)

onde,

(31)

(32)

(33)

(34)

(35)

(36)

com

(37)

A mesma abordagem é usada para o cálculo de $v$ , integrando a Eq. (3). Assim, são obtidos três sistemas lineares a serem resolvidos:

(38)

onde $A^{\phi }$ representa a matriz de coeficientes da variável $\phi$ ( $h,u$ ou $v$ ), enquanto $B^{\phi }$ é o vetor que contém as condições de contorno e termos fontes.

3. O algoritmo

O código de simulação computacional foi programado em Linguagem Fortran, compilado na versão Fortran Power Station 4.0 licenciada pela Microsoft Corporation. A escolha por esse programa deve-se à sua velocidade de execução e alta capacidade de armazenamento de dados. Neste trabalho, o domínio físico retangular $D=[x_{a},x_{b}]\times [y_{a},y_{b}]$ é dividido em $N_{x}$ e $N_{y}$ volumes de controle com comprimentos dados por $\Delta x,\Delta y$ nas direções $x$ e $y$ , respectivamente. Da mesma forma, o espaço temporal $[t_{0},t_{f}]$ é dividido em $N_{t}$ passos de tempo uniformes iguais a $\Delta t$ . A malha cartesiana ordenada e estruturada é varrida de acordo com o par $(i,j)$ da esquerda para a direita aumentando o comprimento $(i-1)\Delta x$ na direção $x$ , e de baixo para cima aumentando o comprimento $(j-1)\Delta y$ na direção $y$ , começando no centróide $(1,1)$ que representa o ponto $(x_{a}+\Delta x/2,y_{a}+\Delta y/2)$ . Desta forma, faz-se a seguinte correspondência:,

(39)

Após a discretização do domínio, são lidas a topografia da parte inferior, representada pela superfície $z_{0}(x,y)$ , e as condições iniciais $h_{0},u_{0}$ e $v_{0}$ . O sistema linear é resolvido por um método de Gauss-Seidel de acordo com a direção varrida anteriormente descrita. Ele itera até que um dos dois critérios seja atingido: o número máximo residual $tol$ ou máximo de iterações $itmax$ . Se for atingido resíduo $R<tol$ em um número $k$ de iterações e $2k<itmax$ , então continua-se o processo iterativo até $2k$ iterações. Ou seja, ao atingir a tolerância do resíduo máximo, faz-se mais $k$ iterações, isso torno o resíduo bem inferior ao residual máximo permitido. Para este trabalho, o número máximo de iterações para cada intervalo de tempo é fixado em $itmax=1000$ iterações, enquanto a tolerância para o resíduo máximo global permitido é igual a $tol=10^{-11}$ , sendo o resíduo global $R$ definido como:

(40)

onde $N_{x}N_{y}$ é o número total de volume de controles, e os resíduos $r_{h},r_{u},r_{v}$ são a soma dos resíduos individuais de cada volume de controle,calculados incrementalmente de acordo com a Eq. (41):

(41)

onde $d_{\phi }$ representa a diferença absoluta entre os dois lados da equação linearizada para a variável $\phi$ . O algoritmo completo é apresentado abaixo. O código também permite usar a última solução para inicializar uma nova execução com novas condições de contorno, se necessário.

Declare as matrizes principais $h,u$ e $v$ e as variáveis secundárias $h_{0},v_{0}$ e $u0.$

Leia os dados do modelo: tamanho da malha espacial e temporal, coeficiente de atrito.

Calcule o comprimento dos passos espaciais e temporal

Faça a leitura da topografia $(z_{0}).$

Leia as condições iniciais $h_{0},u_{0}$ e $v_{0}.$

Calcule em cada passa de tempo as variáveis principais $h,u,v$

Se $1\leq k\leq N_{t}$ então

Calcular as variáveis no contorno

Se $2\leq i\leq N_{x}-1$ e $2\leq j\leq N_{y}-1$ então

Interpolar $h,u,v,$ e $z_{0}$ nos pontos de integração nas faces do volume de controle

Calcular os coeficientes do sistema linear de $h$

Calcular o resíduo $r_{h}$

Se $h(i,j)>z(i,j)$ então

Calcular $h(i,j)$

Calcular a coluna de água $H(i,j)=h(i,j)-z(i,j)$

Atualizar os valores de $h$ nas faces dos volumes de controle

Calcular os coeficientes do sistema linear de $u$

Calcular o resíduo $r_{u}$

Calcular $u(i,j)$

Atualizar os valores de $u$ nas faces dos volumes de controle

Calcular os os coeficientes do sistema linear de $v$

Calcular o resíduo $r_{v}$

Calcular $v(i,j)$

$\,\,\,$ caso contrário

Faça $h(i,j)=z(i,j)$ e $r_{h},u,v$ iguais a zero

$\,\,\,$ (o nível da água é calculado até $h=z_{0}$ .)

$\,\,\,$ (fim da malha varrida)

Teste de Convergência

Atualização de $u,v$ e $h$ no tempo $(u_{0}=u,v_{0}=v,h_{0}=h)$

$\,\,\,$ (fim do loop temporal)

Faça o pós-processamento (tabelas e gráficos)

4. Validação da solução numérica

Com o objetivo de validar a abordagem numérica apresentada na seção anterior, esta seção apresenta uma comparação com os resultados obtidos para um modelo de abertura instantânea de comportas apresentado em [10]. Este problema consiste em um domínio quadrado com $200\,m$ de comprimento com uma barragem interna com $10\,m$ de largura, disposta na faixa de $90-100\,m$ , que divide o domínio como é mostrado na Figura 2. Para simulação da abertura de comportas, considera-se uma abertura não simétrica e instantânea com $75\,m$ de largura entre $30$ e $105\,m$ . Essa abertura causa uma liberação sutil da coluna de água a jusante, resultando em um fluxo com gradientes espaciais e temporais acentuados.

Figura 2. Domínio do modelo

Todos os parâmetros numéricos e físicos são mostrados na Tabela (1). De acordo com [10], neste problema a resistência ao escoamento é desprezada e sua parte inferior é plana.

Tabela 1. Parâmetros para simulação de abertura da comporta

Tamanho da malha	Tamanho dos elementos	$z_{0}$	Passo de tempo	Coeficiente de Manning	$g$
$200\,m\,\times \,200\,m$	$5\,m\,\times \,5\,m$	$0$	$0.01\,s$	$0.0$	$9.81\,m/s^{2}$

4.1 Condições iniciais

Para o tempo inicial, o fluido é quiescente, ou seja, $u=v=0$ para todo o domínio, enquanto a coluna d'água inicial é dada por:

(42)

com o nível da água à montante $h_{l}=10\,m,\,$ e a jusante $h_{r}=5\,m$ . Consequentemente, a razão entre os níveis de água é $h_{r}/h_{l}=0.5\,m.$

4.2 Condições de Contorno

Por uma questão de clareza, esta seção é apresentada aqui para mostrar de forma prática como as condições de contorno são tratadas neste trabalho usando para isso, o caso de validação como exemplo. As abordagens apresentadas aqui também são aplicadas a outros casos de maneira similar, respeitando, é claro, a configuração geométrica de cada um. Como apontado em [21], existem pelo menos três maneiras de lidar com as condições de contorno: realizar um balanço para cada volume de controle na fronteira, chamado simplesmente de equilíbrio de contorno; considerar um meio volume na fronteira com seu centróide pertencente ao próprio limite; e, por último, considerar volumes fictícios em torno dos limites, isso significa que o valor de limite prescrito deve ser igual à interpolação (qualquer que seja) entre os centróides do volume real de controle e o novo volume de controle fictício. Esta última abordagem é preferida, pois a primeira implica um tratamento matemático distinto para cada fronteira, o que torna o código mais complexo de acordo com a geometria. O caminho de $2^{a}$ ordem, o meio volume, não é conservador porque o valor é imposto, em vez disso, como resultado de um balanço de fluxo. A abordagem de volume fictício permite que toda a malha seja constituída por volumes inteiros, o que é fácil de codificar e tem a mesma natureza conservadora para todos os volumes. Embora este artifício aumente o número de incógnitas do problema, acarretando custo computacional, é um procedimento útil com codificação fácil. Na Figura 3 é apresentado um esboço deste conceito.

Figura 3. Volumes fictícios ao redor da malha real

Utiliza-se a Condição de Neumann nula para o nível da água nas paredes externas e nos contornos da represa. Desse modo,

(43)

(44)

onde $w$ e $e$ são os pontos dos contornos esquerdo e direito, respectivamente; e $s$ e $n$ são os pontos dos contornos inferior e superior, respectivamente, como descrito na Figura 4.

Figura 4. Esquema da condição de contorno para o nível

do fluido

Todas as derivadas nas faces fronteiras ( $w,e,n,$ e $s$ ) podem ser interpoladas como funções dos centróides vizinhos destacados na Figura 4: um externo fictício ( $W,E,S$ ou $N$ ) e um interno $(P)$ pertencente ao domínio. Consequentemente,

(45)

(46)

(47)

(48)

o que resulta em,

(49)

Para as variáveis de fluxo ( $u$ e $v$ ), as paredes são consideradas impermeáveis e todas as componentes de velocidade normal aos contornos são nulas. Esta condição de Dirichlet é imposta considerando o esquema de diferenciação central, ou seja,

(50)

(51)

(52)

(53)

e portanto,

(54)

(55)

Por outro lado, as componentes da velocidade tangencial são calculados considerando derivada nula para cada uma, ou seja, uma Condição de Neumann nula para cada componente tangencial, e então

(56)

(57)

(58)

(59)

e portanto,

(60)

(61)

4.3 Resultados

O algoritmo proposto e a formulação numérica devem ser validados para determinar sua precisão e robustez. Isso é feito por meio da aplicação do modelo e do código numérico a um problema representativo de águas rasas: o problema da abertura instantânea de comportas, amplamente utilizado para testar códigos numéricos ou modelos matemáticos que tratam de vazões superficiais sujeitas a mudanças abruptas, mudanças no tempo e no espaço, e a abertura instantânea de comportas tem esses detalhes físicos e é amplamente estudada na literatura [10,14,22,23,24,25].

4.3.1 Nível da água

A Figura 5 mostra a superfície livre de água após $7.1$ segundos após a abertura da comporta. Os resultados são comparados com aqueles obtidos em [10]. Observa-se na Figura 5 uma frente de onda após a abertura da eclusa e uma onda de depressão no lado esquerdo à barreira. Os níveis de superfície computados estão de acordo com os apresentados por [10] para o mesmo instante de tempo.

	]
(a) Resultado obtido pelo esquema híbrido proposto	(b) Resultado obtido em [10]
Figura 5. Comparação do comportamento da superfície livre $3D$ após $7.1\,s$ após a abertura da comporta

A Tabela 2 mostra os valores de nível para algumas posições do domínio para fins de comparação quantitativa.

Tabela 2. Valores de nível de água usando a relação

yx	$0$	$25$	$50$	$75$	$100$	$125$	$150$	$175$	$200$
$0$	$9.99932$	$9.96352$	$9.43066$	$8.76246$	$5.06748$	$7.08372$	$5.31831$	$5.00090$	$5.0$
$25$	$9.99497$	$9.88268$	$9.18459$	$8.70866$	$5.16149$	$6.34225$	$6.43242$	$5.01876$	$5.0$
$50$	$9.98848$	$9.77139$	$8.67804$	$8.08645$	$7.28272$	$6.84739$	$7.01455$	$5.08294$	$5.0$
$75$	$9.98760$	$9.75073$	$8.51860$	$7.88839$	$7.52336$	$6.87814$	$7.10481$	$5.09955$	$5.0$
$100$	$9.99177$	$9.83175$	$8.98309$	$8.49358$	$5.82968$	$6.54152$	$6.75214$	$5.04260$	$5.0$
$125$	$9.99851$	$9.95297$	$9.51005$	$9.17919$	$5.68769$	$5.94785$	$5.60476$	$5.00278$	$5.0$
$150$	$9.99996$	$9.99796$	$9.93089$	$9.72443$	$5.48812$	$5.30022$	$5.00670$	$5.00001$	$5.0$
$175$	$9.99999$	$9.99997$	$9.99829$	$9.98574$	$5.00626$	$5.00135$	$5.00001$	$5.00000$	$5.0$
$200$	$10.0000$	$9.99999$	$9.99999$	$9.99985$	$3.75001$	$5.00000$	$5.00000$	$5.00000$	$5.0$

Duas linhas de amostra foram extraídas em $7.1$ segundos antes da abertura da comporta $(x=75\,m)$ e após a abertura $(x=105\,m)$ e os dados foram comparados a [10], onde são confrontados outros dois esquemas numéricos: MacCormack e Gabutti (descritos em [26] e [27], respectivamente), e ambos os métodos são de precisão de segunda ordem e requerem um cálculo mais complexo, uma vez que uma etapa do preditor deve ser avaliada em $t+\Delta t/2.$ A Figura 6 apresenta a comparação dos perfis transversais do nível da água após $7.1$ segundos da abertura da barreira interna, e os dados estão em boa concordância com [10]. Algumas diferenças podem ser observadas antes e depois da abertura da barragem, principalmente com a onda de depressão antes do portão, mas os níveis estão na mesma ordem e com o mesmo comportamento físico. Além disso, observa-se algumas elevações locais na frente de onda.

$Comparação do perfil do nível de água em duas linhas de amostra em x=115\,m e x=75\,m$

Figura 6. Comparação do perfil do nível de água em duas linhas de amostra em

e

A Figura 7 mostra os perfis longitudinais em $y=70\,m$ e $y=150\,m$ a $7.1$ segundos. A linha $h(x,150)$ está quebrada porque intercepta a parede da barragem, enquanto a linha $h(x,75)$ passa pela abertura da comporta. Pode ser visto na linha $h(x,150)$ que a onda frontal após a abertura da eclusa mostra uma elevação de pico que é devida à diferença de velocidade entre a frente de onda e a onda de depressão em movimento a jusante. A frente de onda é desacelerada pelo corpo de água em repouso na frente dela. Embora esse comportamento não seja mostrado nos resultados de [10], a posição da frente de onda e o início da onda de depressão são coincidentes. O perfil $h(x,150)$ tem uma melhor correspondência, uma vez que nesta posição o comportamento do fluido é mais suave e menos dinâmico devido à presença da parede. É evidente que o local mais sensível a ser resolvido é a região de abertura da comporta devido ao abrupto colapso da coluna de água.

$Comparação das linhas de amostra longitudinais do nível de água em y=70\,m e y=150\,m$

Figura 7. Comparação das linhas de amostra longitudinais do nível de água em

e

As mudanças no nível da água podem ser monitoradas por um período de tempo, em um ponto fixo do domínio. Na Figura 8, pode-se ver a amostragem do nível de água em $P(115,70)$ , localizado na frente da abertura da comporta que permite capturar o avanço da frente de onda, durante $7.1$ segundos.

Figura 8. Comparação do tempo de amostragem do nível da água em

Pode ser visto em ambos os resultados que a frente de onda leva $1$ segundo após a abertura para alcançar o ponto de teste $P$ , e depois de $3$ segundos, o nível da água está próximo de $7\,m$ em ambos os casos. Novamente, há um comportamento mais dinâmico em comparação com os resultados de [10], que é mais suave com oscilações brandas devido a utilização de termos dissipativos artificiais, enquanto os resultados obtidos pelo esquema híbrido proposto apresentam uma onda de pico que tende a ser atenuada com o tempo, de acordo com o fenômeno físico esperado.

4.3.2 Campo de Velocidade

A Figura 9 compara os resultados de [10] e deste trabalho para o campo de velocidade no instante de tempo igual a $7.1$ segundos. Os vetores do campo de velocidade estão em muito boa concordância apresentando as mesmas direções, magnitudes (representadas pelos tamanhos das flechas) com um espalhamento pela mesma área ao redor da abertura da comporta.


(a) Campo de velocidade usando o esquema híbrido	(b) Campo de velocidade apresentado em [10]
Figura 9. Comparação do campo de velocidade em torno da abertura da comporta em $7.1$ segundos

4.4 Considerações

Apesar das oscilações locais ou da suavidade, os resultados estão em boa concordância geral com as frentes de onda e ondas de depressão de acordo com os dados da literatura, o que reflete a habilidade da abordagem numérica proposta em lidar e capturar o comportamento físico do fluxos de águas rasas que validam o modelo e a abordagem numérica. É importante enfatizar que, em [10] artifícios numéricos foram usados como uma viscosidade artificial para lidar com as oscilações da superfície livre, o que torna o resultado de [10] mais difusivo e quase sem flutuações locais no nível da superfície. Por outro lado, neste trabalho, nenhuma abordagem similar foi empregada para lidar com uma possível instabilidade numérica devido à aproximação de tempo de alta ordem. O esquema de tempo implícito proposto de $1^{a}$ ordem neste trabalho mostrou-se capaz de capturar o comportamento físico e os gradientes de tempo rígidos sem etapas de tempo proibitivas. E a formulação UDS em associação com a abordagem CDS foi capaz de lidar com mudanças abruptas nos campos de variáveis sem instabilidades ou oscilações não físicas. Todas as oscilações de superfície livre têm um significado físico quando as diferenças de velocidade entre as porções distintas da onda são contabilizadas porque essas diferenças compactam a onda e aumentam sua amplitude. Como apontado em [10], as oscilações numéricas puras não são amortecidas com o tempo, em vez disso, elas tendem a ser amplificadas com gradientes rígidos, o que faz com que o sistema divirja. Porém fazendo uso do esquema numérico híbrido aqui apresentado, a solução se mantêm suave com resíduos muito baixo, inferiores a tolerância máxima de $10^{-11}$ . Mesmo para grandes intervalos de tempo, não há crescimento no valor destes resíduos, atingindo valores inferiores a tolerância máxima em poucas iterações.

5. Verificação numérica

Uma validação numérica assegura que o modelo é capaz de fornecer uma resposta física correta sob certas condições iniciais e de contorno. No entanto, mais testes são feitos para descobrir se, sob outro cenário, o modelo computacional é capaz de prever o comportamento físico do sistema. Esta etapa, juntamente com a validação anterior, permite definir se o modelo computacional é útil como uma ferramenta preditiva para qualquer fluxo de água superficial sob condições diversas.

5.1 Teste da condição de leito seco

Para avaliar a robustez do código, o modelo da Seção 3 é considerado, fazendo a razão $h_{r}/h_{l}$ tender a zero, isto é, $h_{r}$ é considerado com uma fina película líquida, como uma “wet condition”. Esta situação é geralmente uma situação física difícil de resolver devido às singularidades que podem surgir nos termos fontes de atrito, e o possível mau condicionamento da matriz de coeficientes devido a valores nulos que surgirão ao longo dos cálculos quando $h$ tende a zero. Nesta situação, as equações de águas rasas tendem a enfraquecer seu significado físico. Em vista disso, $h_{r}$ não pode ser exatamente zero, e o código deve ser capaz de lidar com essa situação. Na Figura 10, pode-se observar o nível da água para os primeiros $7.1$ segundos com $h_{r}/h_{l}=0.00002$ .

Figura 10. Visão tridimensional da superfície livre em

segundos após a abertura da comporta com

Os valores dos níveis são reproduzidos na Tabela 3, com todas as unidades em metros.

Tabela 3. Valores de nível de água usando

yx	$0$	$25$	$50$	$75$	$100$	$125$	$150$	$175$	$200$
$0$	$9.9984$	$9.9274$	$9.1889$	$8.4453$	$0.0303$	$0.8589$	$1.2135$	$0.0002$	$0.0002$
$25$	$9.9904$	$9.8118$	$8.9461$	$8.3330$	$0.4597$	$1.1416$	$1.1143$	$0.0583$	$0.0002$
$50$	$9.9789$	$9.6532$	$8.3391$	$7.3257$	$5.4764$	$3.3378$	$1.9742$	$0.4390$	$0.0002$
$75$	$9.9773$	$9.6221$	$8.1432$	$6.9206$	$5.6710$	$3.7372$	$2.1506$	$0.6072$	$0.0002$
$100$	$9.9848$	$9.7398$	$8.7161$	$8.0375$	$3.4672$	$1.9603$	$1.4454$	$0.1380$	$0.0002$
$125$	$9.9969$	$9.9163$	$9.3453$	$8.9770$	$0.1084$	$0.5429$	$0.6798$	$0.0058$	$0.0002$
$150$	$9.9999$	$9.9950$	$9.8806$	$9.6130$	$0.0002$	$0.0235$	$0.0028$	$0.0002$	$0.0002$
$175$	$9.9999$	$9.9999$	$9.9960$	$9.9728$	$0.0002$	$0.0002$	$0.0002$	$0.0002$	$0.0002$
$200$	$10.000$	$9.9999$	$9.9999$	$9.9996$	$0.00015$	$0.0002$	$0.0002$	$0.0002$	$0.0002$

5.2 Considerações

Em Fennema and Chaudhry [10], Fernández-Pato et al. [13] e Sheu and Fang [15] e são apontadas as dificuldades em lidar com modelo que representam leitos quase secos, devido as possíveis singularidades. O esquema numérico híbrido CDS/UDS aqui proposto mostrou-se eficaz e com baixo custo computacional para simular escoamentos transientes em águas rasas para leitos quase secos, o que evidencia a acurácia e robustez do esquema.

6. Escoamento em um reservatório de contenção subterrâneo: uma aplicação do modelo completo

Os reservatórios subterrâneos são projetados para zonas urbanas populosas, onde não há lugar para armazenar água da chuva. Eles são construídos com uma rede de pilares de sustentação e canais de drenagem na parte inferior, como ilustrado na Figura 11.

Figura 11. Reservatório de águas pluviais subterrâneas

Para testar a solução numérica obtida para as equações completas das águas rasas, analisou-se o fluxo em um reservatório com obstáculos como colunas estruturais, e drenos controlados por saídas estreitas como fendas na barragem. A topografia de fundo foi considerada com inclinações não nulas e com resistência ao escoamento no fundo do canal, com coeficiente de atrito de Manning $N=0.013m/s^{1/3}$ , que quantifica a resistência ao escoamento no fundo do canal para superfícies de cimento, apresentado em [28]. O fundo tem a presença de canais que conduzem a água pelas aberturas de saída até uma área aberta sem paredes. Pode-se dizer que o problema considerado é uma composição do problema de abertura instantânea de comportas, para leitos molhado/seco, com topografia mais complexa e configuração de domínio. A Figura 12 representa a configuração do domínio físico. O domínio consiste em um reservatório retangular de $12m\times 24m$ de dimensão com nove pilares de sustentação.

Figura 12. Domínio físico do reservatório subterrâneo

Como ilustrado na Figura 13, a parte inferior do reservatório possui dois canais de drenagem com seção trapezoidal e um declive sutil em direção às fendas da barragem.


(a) representação tridimensional do fundo	(b) Perfil frontal do fundo do reservatório
Figura 13. Topografia do fundo do reservatório

As configurações da topografia do fundo são descritas por $z_{0}$ na Eq.(62).

z_{0}(x,y,t)=\left\{{\begin{array}{lll}0.2,\,\,\,\,\,\,\,\,\,\,{\textrm {se}}\,\,\,\,0\leq \,x\,\leq 0.5\,\,\,{\textrm {e}}\,\,\,\,0\leq \,y\,\leq 8,\\\\1.2-2x,\,\,\,\,\,\,\,\,\,\,{\textrm {se}}\,\,\,\,0.5<\,x\,<0.6\,\,\,{\textrm {e}}\,\,\,\,0\leq \,y\,\leq 8,\\\\0.0,\,\,\,\,\,\,\,\,\,\,{\textrm {se}}\,\,\,\,0.6\leq \,x\,\leq 0.9\,\,\,{\textrm {e}}\,\,\,\,0\leq \,y\,\leq 8,\\\\-1.8+2x,\,\,\,\,\,\,\,\,\,{\textrm {se}}\,\,\,\,\,0.9<\,x\,<1.0\,\,\,{\textrm {e}}\,\,\,\,0\leq \,y\,\leq 8,\\\\0.2,\,\,\,\,\,\,\,\,\,\,{\textrm {se}}\,\,\,\,1.0\leq \,x\,\leq 11.0\,\,\,{\textrm {e}}\,\,\,\,0\leq \,y\,\leq 8,\\\\22.0-2x,\,\,\,\,\,\,\,\,\,\,{\textrm {se}}\,\,\,\,11.0<\,x\,<11.1\,\,\,{\textrm {e}}\,\,\,\,0\leq \,y\,\leq 8,\\\\0.0,\,\,\,\,\,\,\,\,\,\,{\textrm {se}}\,\,\,\,11.1\leq \,x\,\leq 11.4\,\,\,{\textrm {e}}\,\,\,\,0\leq \,y\,\leq 8,\\\\-22.8+2x\,\,\,\,\,\,\,\,\,\,\,{\textrm {se}}\,\,\,\,11.4<\,x\,<11.5\,\,\,{\textrm {e}}\,\,\,\,0\leq \,y\,\leq 8,\\\\0.2,\,\,\,\,\,\,\,\,\,\,{\textrm {se}}\,\,\,\,11.5\leq \,x\,\leq 12.0\,\,\,{\textrm {e}}\,\,\,\,0\leq \,y\,\leq 8,\\\\0.0,\,\,\,\,\,\,\,\,\,\,{\textrm {se}}\,\,\,\,0\leq \,x\,\leq 12.0\,\,\,{\textrm {e}}\,\,\,\,8.15\leq \,y\,\leq 24.\\\end{array}}\right.

(62)

No reservatório, o nível da água está inicialmente em repouso e tem $2.5\,m$ de profundidade, enquanto na zona de inundação, há apenas uma camada fina de água, ou seja, “wet condition”, a altura da água é tão pequena quanto possível, como foi testado na Seção 5.1. Neste caso, a razão fixada é $h_{r}/h_{l}=0.0002$ , como mostrado na Figura 14. Todas as condições de contorno para o nível da água são as mesmas do problema apresentado na Seção 4.

Figura 14. Condições iniciais para o nível da água

Com relação aos componentes de velocidade, é necessário dizer que: como as paredes do reservatório, os pilares e paredes laterais do domínio externo são considerados impermeáveis; o fluxo da água que é armazenada no reservatório ocorre apenas nas saídas da parede norte do reservatório e na extremidade norte do domínio externo, como indicado na Figura 12. Mais precisamente, a componente de velocidade $u$ nas paredes internas leste-oeste e nos contornos internos é nula, enquanto a componente de velocidade $v$ é considerada nula no limite sul e nos limites internos norte-sul. Por outro lado, nos contornos norte-sul para a componente $u$ e nos contornos leste-oeste para o componente $v$ respectivamente, é adotada a condição de Neumann nula, como na Seção 4. Para o domínio externo, uma vez que existem apenas paredes laterais, o fluxo é livre para todos os pontos da extremidade norte. Essa condição de contorno é descrita pelas expressões:

(63)

(64)

sendo $C=1.0(2g)^{\frac {1}{2}}$ um coeficiente arbitrário, $h(m)$ o nível da água e $g(m/s^{2})$ a aceleração da gravidade. Para a simulação computacional do modelo proposto utilizamos uma malha estruturada, de espaçamento $\Delta \,x=\Delta y=0.05\,m,$ no espaço bidimensional, com intervalo de tempo de $\Delta \,t=0.001\,s$ . Como há uma topografia especificada na parte inferior do reservatório, é necessário estabelecer uma condição compatível com a física do modelo: se o nível da água atingir o fundo do reservatório, o fluxo será considerado nulo neste ponto. Mais precisamente se, $h(x,y,t)>z_{0}(x,y,t)$ os cálculos para $h,u,v$ são feitos para o próximo passo de tempo, caso contrário faz-se: $h(x,y,t)=z_{0}(x,y,t)$ e as componentes de velocidade são anuladas neste ponto. O volume de fluido armazenado no interior do reservatório foi calculado em cada ponto de tempo. E foi estabelecido como critério de parada no código, o volume mínimo $V=0.01\,m^{3}.$

6.1 Resultados

Devido às aberturas na parede norte do reservatório, há um declínio no nível da água dentro do reservatório e ondas de elevação no domínio externo, especialmente nas proximidades das aberturas. A Figura 15 é uma representação em perspectiva da superfície da água em todo o domínio para $t=21\,s$ e $t=70\,s$ após a abertura das saídas.

$t=21\,s$	$t=70\,s$
(a) $t=21\,s$	(b) $t=70\,s$
Figura 15. Visão tridimensional do nível da água após o início do despejo

A representação gráfica do campo vetorial é usada para a análise da intensidade e direção do fluxo, em cada ponto da malha. A Figura 16 mostra as direções e a intensidade do fluxo nos primeiros $4\,s$ do fluxo. Observa-se que a direção do fluxo ocorre de dentro para fora do reservatório, a partir das saídas, com gradientes mais acentuados próximos às aberturas.

$t=1\,s$	$t=2\,s$
(a) $t=1\,s$	(b) $t=2\,s$
$t=3\,s$	$t=4\,s$
(c) $t=3\,s$	(d) $t=4\,s$
Figura 16. Campo de velocidade nos primeiros $4$ segundos

Para $t=15\,s$ observa-se que ainda não há fluxo na região estreita à frente da parede norte do reservatório, nas seções onde não há aberturas. Por sua vez, nas regiões à frente das aberturas, o fluxo é mais acentuado e o campo de velocidade assume um perfil hiperbólico no domínio externo, como na Figura 17a. Por outro lado, como mostra a Figura 17b, a magnitude do vetor velocidade torna-se menos acentuada e mais uniforme, após $70\,s$ de simulação do fluxo.

$t=15\,s$	$t=70\,s$
(a) $t=15\,s$	(b) $t=70\,s$
Figura 17. Campo de velocidade após o início do despejo

O perfil transversal do nível da água pode ser analisado em um determinado instante de tempo. Na Figura 18, o perfil $h(0,75,y)$ é representado para $t=10\,s$ . Uma inclinação significativa no nível da água é observada na faixa de $8\,m<y<8.15\,m$ , o que era esperado devido à condição inicial descontínua imposta.

$Perfil do nível de água em uma linha de amostra em x = 0.75\,m para t = 10\,s$

Figura 18. Perfil do nível de água em uma linha de amostra em

para

O código implementado foi utilizado para análises como: variação de volume no reservatório em função do tempo, o comportamento do fluxo ao redor dos pilares e na vizinhança das paredes. Em todas essas análises, os resultados obtidos são satisfatórios e consistentes com a física do problema.

7. Conclusões

Foi proposto um método de volume finito implícito de $1^{a}$ ordem com uma discretização híbrida CDS/UDS para resolver as equações completas de águas rasas. O algoritmo mostrou-se capaz de lidar com fortes problemas instáveis como a abertura instantânea de comportas, vazão em topografias não uniformes com condição dry/wet, com e sem condições de saídas. Os resultados mostram boa concordância com a literatura com dinâmica de frente de onda mais complexa. Essas diferenças locais são explicadas devido à ausência de artifício numérico, como a viscosidade artificial, que amortece o comportamento da superfície livre e suaviza as ondas. A formulação numérica proposta tem a vantagem de ser mais fácil de codificar em comparação com aquelas disponíveis na literatura e, além disso, o apelo físico para a conservação de fluxo que é característica do método de volumes finitos, torna-se uma opção com a mesma precisão física, porém simples e rápida para ser implementada. A convergência do método iterativo de Gauss-Seidel, mesmo impondo máximo residual muito pequeno, é obtida com poucas iterações, o que garante baixo custo computacional. Como ferramenta preditiva, o modelo computacional tem mostrado bons resultados desde a configuração simples até a mais complexa, com precisão e robustez. A solução numérica obtida pode ser usada para simular fluxos sobre topografias reais encontradas na natureza.

Declaração de disponibilidade de dados

Todos os dados, modelos ou códigos gerados ou utilizados durante o estudo estão disponíveis pelo autor correspondente por solicitação. Lista de itens:

o algoritmo utilizado para obtenção da solução numérica híbrida proposta;
o arquivo de dados da Tabela 2;
o arquivo de dados da Tabela 3;
os dados usados para gerar as figuras.

Agradecimentos

Este estudo foi financiado em parte pela Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - Brasil (CAPES) - Finance Code 001. Os autores agradecem também ao Departamento de Matemática da Universidade Tecnológica Federal do Paraná (UTFPR), Campus Campo Mourão e ao Programa de Pós-Graduação em Métodos Numéricos (PPGMNE) da Universidade Federal do Paraná (UFPR).

Referências

[1] Imbonden D M. The lakes handbook, immnology and limnectic ecology. Blackwell Scienc Ltd., 2004.

[2] Reynolds C. The ecology of freshwater phytoplankton. Cambridge University Press, 1984.

[3] Mahmood K., Yevjevich V. Unsteady flow in open channels. Water Resources Publications, 1st Edition, 1975.

[4] Vreugdenhil C.B. Numerical Methods for Shallow-Water flow. Springer Science+Business Media Dordrecht, 1994.

[5] García-Navarro P., Alcrudo F., Priestley A. An implicit method for water flow modelling in channels and pipes. J. Hydraul. Eng., 32:721–742, 1994.

[6] Burguete J., García-Navarro P. Implicit schemes with large time step for non-linear equations: Application to river flow hydraulic. Int. J. Numer. Methods , 46:607–636, 2004.

[7] Cantero-Chinchilla F.N., Castro-Orgaz O., Dey S., Ayuso J.L. Nonhydrostatic dam break flows. I: Physical equations and numerical schemes. J. Hydraul. Eng., 142(12):1–19, 2016.

[8] Soler J., Bladé E., Sánchez-Juny M. Ensayo comparativo entre modelos unidimensionales y bidimensionales en la modelización de la rotura de una balsa de materiales sueltos erosionables. Revista International de Métodos Numéricos para Cálculo y Diseño en Ingeniería, 28(2):103–105, 2012.

[9] Brufau P., García-Navarro P. Esquemas de alta resolución para resolver las ecuaciones de shallow water. Revista International de Métodos Numéricos para Cálculo y Diseño en Ingeniería, 16(4):493–512, 2000.

[10] Fennema R.J., Chaudhry M.H. Explicit Methods for 2-D transient free surface flows. J. Hydraul. Eng., 116(8):1013–1034, 1990.

[11] Nikolos I.K., Delis A.I. An unstructured node-centered finite volume scheme for shallow water flows with wet/dry fronts over complex topography. Comput. Methods. Appl. Mech. Engrg., 198:3723–3750, 2009.

[12] Valiani A., Caleffi V., Zanni A. Case study: malpasset dam-Break simulation using a two-dimensional finite volume method. J. Hydraul. Eng., 128(5):460–472, 2002.

[13] Fernández-Pato J., Morales-Hernández M., García-Navarro P. Implicit finite volume simulation of 2D shallow water flows in flexible meshes. Comput. Methods. Appl. Mech. Engrg., 328:1–25, 2018.

[14] Melo A.R., Gramani L.M., Kaviski E. High-order CE/SE scheme for dam-break flow simulation. Revista International de Métodos Numéricos para Cálculo y Diseño en Ingeniería, 34(1):1–12, 2018.

[15] Sheu T.W.H., Fang C.C. High resolution finite-element analysis of shallow water equations in two dimensions. Comput. Methods. Appl. Mech. Engrg., 190:2581–2601, 2001.

[16] Hirsch J.J. Water waves: The mathematical theory with applications. Interscience Publishers, Edition, 2007.

[17] Blazek J. Computational fluid dynamics: principles and applications. Elsevier Science Ltd, 1st Edition, 2001.

[18] Versteeg H.K., Malalasekera W. An Introduction to computational fluid dynamics. The finite volume method. Pearson Education Limited, 2nd Edition, 2007.

[19] Maliska C.R. Transferência de calor e mecânica dos fluidos computacional. LTC, 2nd Edition, 2013.

[20] Patankar, S. (1980) Numerical heat transfer and fluid flow. Hemisphere Publishing Corporation, Edition

[21] Lemos, E. M. D. (2011) Implementação de um método de volumes finitos de ordem superior com tratamento multibloco aplicado à simulação de escoamentos em fluidos viscoelásticos. Phd. thesis, UFRJ, Rio de Janeiro, Brasil

[22] Rezende, R. V. P. and Almeida, R. A. and Souza, A. A. U. and Souza, S. M. A. G. U. (2015) A two-fluid model with a tensor closure model approach for free surface flow simulations, Volume 122. Chemical Engineering Science 596–613

[23] Hirt, C. W. and Nichols, B. D. (1981) Volume of fluid (vof) method for the dynamics of free boundaries, Volume 39. Journal of ComputationalPhysics 201–225

[24] Jahanbakhsh, E. and Panahi, R., Seif, M. S.(2007) Numerical simulation fo three-dimensional interfacial flows, Volume 17. Int. J. Numer. Methods Heat Fluid Flow 384–404

[25] Villota, A., Codina, R. (2018) Approximation of the shallow water equations with higher order finite elements and variational multiscale methods, Volume 34 (1). Revista International de Métodos Numéricos para Cálculo y Diseño en Ing. 1–25

[26] MacCormack, R. W. (1969) The effect of viscosity in hypervelocity impact cratering, Volume. American Institutue Aeronautics Astronautics 69–354

[27] Gabutti, B. (1983) On two upwind finite-difference schemes for hyperbolic equations in non-conservative form, Volume 11. Comput. Fluids 3 207–230

[28] Porto, R. M. (2006) Hidráulica Básica. EESC-USP, 4 th Edition

Latest revision as of 16:22, 23 February 2021

Resumo

Abstract

1. Introdução

2. Modelo matemático e solução numérica

2.1 Sistema de Equações Governantes

2.2 Discretização do domínio e linearização das equações governantes

3. O algoritmo

4. Validação da solução numérica

4.1 Condições iniciais

4.2 Condições de Contorno

4.3 Resultados

4.3.1 Nível da água

4.3.2 Campo de Velocidade

4.4 Considerações

5. Verificação numérica

5.1 Teste da condição de leito seco

5.2 Considerações

6. Escoamento em um reservatório de contenção subterrâneo: uma aplicação do modelo completo

6.1 Resultados

7. Conclusões

Declaração de disponibilidade de dados

Agradecimentos

Referências

Document information

Document Score

Share this document

Keywords

claim authorship

@@ Line 24: / Line 24: @@
 |}
 -->
+==Resumo==
-==Abstract==
 Este trabalho apresenta a simulação hidráulica do escoamento em um reservatório, utilizando o método dos volumes finitos para resolver o sistema de equações que modelam o fluxo bidimensional em águas rasas, negligenciando as tensões tangenciais. Para resolver as equações, utilizou-se um esquema híbrido de   volumes finitos. A solução do sistema de equações linearizadas foi obtida por  implementação computacional do método iterativo de Gauss-Seidel. Para ilustrar a aplicação do esquema numérico em engenharia hidráulica, apresenta-se o estudo do escoamento em um reservatório subterrâneo com pilares internos, cujo fluxo é gerado pela abertura de dois portões de saída. O código desenvolvido permite a análise temporal do volume e do fluxo no reservatório,  a interpretação gráfica do fenômeno físico investigado bem como o cálculo da precisão do modelo. Os resultados obtidos mostram a interpretação física esperada, boa concordância com os dados da literatura, boa precisão, estabilidade e baixo custo computacional.
-'''keywords'''
+'''Palavras chave''': Equações de águas rasas,  volumes finitos,  interpolação híbrida,  estabilidade numérica,  simulação de leito seco,  reservatório de contenção
-Equações de águas rasas,  Volumes Finitos,  Interpolação híbrida,  Estabilidade numérica,  Simulação de leito seco,  Reservatório de contenção.
+==Abstract==
-[label1]Eloy Kaviski
+This work presents the numeric simulation of flow in a reservoir, using the finite volume method for  solver the system of equations that model the two-dimensional flow in shallow water, neglecting the tangential tensions.  The solution of the system of linearized equations was obtained by computational implementation of the Gauss-Seidel iterative method. To illustrate the application of the numerical scheme in hydraulic engineering, it was considered the study of flow in an underground reservoir with internal pillars, whose flow is generated by the opening of two outlet gates. The employed code allowed the temporal analysis of the volume and the flow in the reservoir, the graphical interpretation of the investigated physical phenomenon as well as, the calculation of the precision of the model. The results obtained show the expected physical interpretation, good agreement with literature data, good precision,stability and low computational cost.
-==HYBRID NUMERICAL SCHEME FOR SIMULATION OF TRANSIENT FLOW IN A PLUVIAL WATER RESERVOIR==
+'''Keywords''': Shallow water equations, finite volumes, hybrid interpolation, numerical stability, dry bed simulation, containment reservoir
-'''Abstract'''
+==1. Introdução==
-This work presents the numeric simulation of flow in a reservoir, using the finite volume method for  solver the system of equations that model the two-dimensional flow in shallow water, neglecting the tangential tensions.  The solution of the system of linearized equations was obtained by computational implementation of the Gauss-Seidel iterative method. To illustrate the application of the numerical scheme in hydraulic engineering, it was considered the study of flow in an underground reservoir with internal pillars, whose flow is generated by the opening of two outlet gates. The employed code allowed the temporal analysis of the volume and the flow in the reservoir, the graphical interpretation of the investigated physical phenomenon as well as, the calculation of the precision of the model. The results obtained show the expected physical interpretation, good agreement with literature data, good precision,stability and low computational cost.
+A ocorrência ou agravamento de inundações nas grandes cidades, bem como a quebra de barragens, tornaram-se problemas cada vez mais recorrentes, causando mortes e fortes impactos ambientais. Nesse cenário, fica evidente a necessidade de conhecimento prévio do comportamento dos fluxos de superfície livre e simulações de abertura de comportas de barragens para manejo de recursos hídricos, avaliação de riscos e usos para obtenção de melhor disponibilidade visando a preservação ambiental. Para tanto, é necessário utilizar metodologias que possam descrever as variáveis hidrodinâmicas envolvidas, como campo de velocidade do fluxo, turbulência e mudanças temporais na altura da superfície da água <span id='citeF-1'></span><span id='citeF-2'></span>[[#cite-1|[1,2]]].  Em <span id='citeF-3'></span>[[#cite-3|[3]]] sugere-se o uso de modelos matemáticos para entender os padrões de fluxo em corpos de água. Modelos matemáticos adequados são ferramentas úteis no gerenciamento de recursos hídricos, pois reduzem o tempo de análise, custos e riscos envolvidos em levar dados experimentais em fluxos reais auxiliando na tomada de decisões e na definição de ações corretivas diante de vazamentos iminentes.  Matematicamente, os fluxos de superfície livre podem ser descritos de várias maneiras. Uma abordagem mais simples e muito aplicada é dada pelas Equações de Águas Rasas, que  representam uma simplificação das Equações de Navier-Stokes e é aplicável aos fluxos nos quais as tensões e acelerações verticais podem ser negligenciadas sem prejuízo do resultado. Essas equações são derivadas por meio da integração das equações governantes (Navier-Stokes e Continuidade) na direção vertical. Sua derivação está fora do escopo deste trabalho, mas pode ser vista em detalhes em <span id='citeF-4'></span>[[#cite-4|[4]]] .  Esse sistema de equações, embora definido em um domínio bidimensional, fornece a altura da coluna de água e as componentes de velocidade em um fluxo instável. No entanto, a obtenção de uma solução analítica para o sistema de equações de águas rasas é, na maioria dos casos, uma tarefa difícil, pois trata-se de um conjunto de equações diferenciais parciais hiperbólicas não lineares com termos fontes. O que requer um tratamento numérico robusto para evitar instabilidades numéricas.   Ao longo das últimas três décadas, muitos avanços foram realizados na solução numérica das equações de águas rasas para simulação de escoamentos unidimensionais aplicados a hidráulica fluvial em rios, canais e tubulações
+<span id='citeF-5'></span><span id='citeF-6'></span><span id='citeF-7'></span><span id='citeF-8'><span id='citeF-9'></span>[[#cite-5|[5,6,7,8,9]]], e modelos bidimensionais aplicados à análise de escoamentos causados pela abertura instantânea de comportas e ruptura de barragens <span id='citeF-10'></span><span id='citeF-11'></span><span id='citeF-12'></span><span id='citeF-13'><span id='citeF-14'><span id='citeF-15'></span>[[#cite-10|[10,11,12,13,14,15]]].  Em <span id='citeF-10'></span>[[#cite-10|[10]]] utilizaram esquemas explícitos de diferenças finitas, na discretização das equações de águas rasas bidimensionais, adicionando condições de estabilidade e viscosidade artificial para suavizar as oscilações de alta frequência. Para validar, os esquemas discretos foram aplicados em um problema típico de engenharia hidráulica: abertura parcial de uma barragem e passagem de uma onda de inundação através de uma contração do canal, adotando a razão entre o nível de água a jusante <math> (h_r) </math> e à montante <math> (h_l ) </math> como sendo, <math> h_r / h_l = 0.5 </math>. Os resultados dos diferentes esquemas foram comparados e a ocorrência de oscilações numéricas foi observada, evidenciando a necessidade de se utilizar esquemas implícitos.  O estudo de caso apontado em <span id='citeF-10'></span>[[#cite-10|[10]]] foi analisado por <span id='citeF-15'></span>[[#cite-15|[15]]] usando um esquema de elementos finitos de alta ordem e uma melhoria na relação <math> h_r / h_l </math> para <math> 0,05 </math> foi obtida, buscando simular a condição de leito seco em uma determinada região do domínio. Em <span id='citeF-12'></span>[[#cite-12|[12]]] usaram dados experimentais em um evento de quebra de barragem para verificar a precisão, estabilidade e confiabilidade de uma solução numérica das equações de águas rasas obtidas por um esquema de volumes finitos bidimensional. O modelo numérico permitiu uma precisão de segunda ordem, tanto no espaço como no tempo, acarretando porém um alto custo computacional, apresentando algumas discrepâncias entre os resultados reais dos medidores e os resultados numéricos nos pontos pesquisados. Os autores destacam a necessidade de simular a propagação de uma onda de inundação para um fundo seco e áspero e otimizar o algoritmo para reduzir o tempo de computação.  Um esquema de volumes finitos não estruturado para malhas triangulares para simular fluxos bidimensionais instáveis em águas rasas com frentes molhadas/secas sobre topografia complexa é apresentado em <span id='citeF-11'></span>[[#cite-11|[11]]]. Utiliza-se em <span id='citeF-11'></span>[[#cite-11|[11]]] para o esquema temporal o método explícito de Runge-Kutta e para o espacial a malha não estruturada. Comparações satisfatórias foram obtidas com dados mensurados e numéricos. Todavia, os autores sugerem a incorporação de um procedimento de escalonamento num esquema implícito temporal, e uma implementação paralela do algoritmo, a fim de reduzir o tempo computacional.   A eficiência de um esquema implícito de volumes finitos é explorada em <span id='citeF-13'></span>[[#cite-13|[13]]] para simulação de fluidos em águas rasas bidimensionais, com frentes secas/molhadas. Adotou-se em <span id='citeF-13'></span>[[#cite-13|[13]]] um tratamento complexo para a discretização espacial, fazendo uso de uma malha não estruturada flexível para obter um melhor ajuste à geometrias irregulares. No entanto, é feito ressalvas quanto ao tempo computacional para obtenção da solução por esquemas implícitos, pois embora um método implícito normalmente exija menos etapas da solução, cada uma exige mais tempo computacional do que as soluções obtidas por um esquema explícito. Os autores apontam a necessidade de se utilizar um solucionador de matrizes que com um número reduzido de iterações alcance a convergência.   Neste contexto, o objetivo principal deste trabalho é formular um esquema híbrido e totalmente implícito de volumes finitos para discretização das equações de águas rasas, a fim de obter uma solução numérica estável com um algoritmo iterativo otimizado, com convergência rápida e baixo custo computacional, para aplicação na simulação de fluxos bidimensionais transitórios em topografias complexas (com declives e sujeitos a resistência ao fluxo). Condições iniciais descontínuas representando leitos úmidos/secos foram adotadas. Para tanto, o trabalho está organizado da seguinte forma: A Seção [[#2 Modelo Matemático e Solução Numérica|2]] apresenta as equações hiperbólicas, que regem a natureza das águas rasas, considerando as perdas de atrito e declives da topografia do domínio computacional e detalhando a configuração numérica destas equações; na Seção [[#3 O algoritmo|3]], o modelo é validado e uma comparação é feita com os resultados obtidos com o trabalho apresentado em <span id='citeF-10'></span>[[#cite-10|[10]]]. A robustez da solução numérica é verificada com uma extrapolação na relação <math> h_r / h_l </math> para <math>0.00002</math>, na Seção [[#4 |4]]; o esquema é aplicado na análise de um problema hidráulico real  na Seção [[#5 Verificação Numérica|5]] e; as conclusões são apresentadas na Seção [[#6 Escoamento em um reservatório de contenção subterrâneo: uma aplicação do modelo completo|6]].
-''Keywords'': Shallow water equations, Finite volumes, Hybrid interpolation, Numerical stability, Dry bed simulation, Containment reservoir.
+==2. Modelo matemático e solução numérica==
-==1 Introduco==
+Na análise do escoamento bidimensional transiente em um canal, cuja topografia do fundo é descrita por uma superfície <math>z_0(x,y)</math>, busca-se determinar o nível do fluido de superfície livre <math>h(x,y,t)</math> e o campo de velocidade do fluxo <math>U(u,v)</math> para cada instante de tempo <math>t</math>, conforme  [[#img-1|Figura 1]].
-A ocorrência ou agravamento de inundações nas grandes cidades, bem como a quebra de barragens, tornaram-se problemas cada vez mais recorrentes, causando mortes e fortes impactos ambientais. Nesse cenário, fica evidente a necessidade de conhecimento prévio do comportamento dos fluxos de superfície livre e simulações de abertura de comportas de barragens para manejo de recursos hídricos, avaliação de riscos e usos para obtenção de melhor disponibilidade visando a preservação ambiental. Para tanto, é necessário utilizar metodologias que possam descrever as variáveis hidrodinâmicas envolvidas, como campo de velocidade do fluxo, turbulência e mudanças temporais na altura da superfície da água <span id='citeF-1'></span><span id='citeF-2'></span>[[#cite-1|[1,2]]].  Em <span id='citeF-3'></span>[[#cite-3|[3]]] sugere-se o uso de modelos matemáticos para entender os padrões de fluxo em corpos de água. Modelos matemáticos adequados são ferramentas úteis no gerenciamento de recursos hídricos, pois reduzem o tempo de análise, custos e riscos envolvidos em levar dados experimentais em fluxos reais auxiliando na tomada de decisões e na definição de ações corretivas diante de vazamentos iminentes.  Matematicamente, os fluxos de superfície livre podem ser descritos de várias maneiras. Uma abordagem mais simples e muito aplicada é dada pelas Equações de Águas Rasas, que  representam uma simplificação das Equações de Navier-Stokes e é aplicável aos fluxos nos quais as tensões e acelerações verticais podem ser negligenciadas sem prejuízo do resultado. Essas equações são derivadas por meio da integração das equações governantes (Navier-Stokes e Continuidade) na direção vertical. Sua derivação está fora do escopo deste trabalho, mas pode ser vista em detalhes em <span id='citeF-4'></span>[[#cite-4|[4]]] .  Esse sistema de equações, embora definido em um domínio bidimensional, fornece a altura da coluna de água e as componentes de velocidade em um fluxo instável. No entanto, a obtenção de uma solução analítica para o sistema de equações de águas rasas é, na maioria dos casos, uma tarefa difícil, pois trata-se de um conjunto de equações diferenciais parciais hiperbólicas não lineares com termos fontes. O que requer um tratamento numérico robusto para evitar instabilidades numéricas.   Ao longo das últimas três décadas, muitos avanços foram realizados na solução numérica das equações de águas rasas para simulação de escoamentos unidimensionais aplicados a hidráulica fluvial em rios, canais e tubulações  [GarciaNavarro1994,Burguete2004,Cant2016,Soler2012,Brufau2000], e modelos bidimensionais aplicados à análise de escoamentos causados pela abertura instantânea de comportas e ruptura de barragens <span id='citeF-8'></span><span id='citeF-9'></span><span id='citeF-10'></span><span id='citeF-11'></span>[[#cite-8|[8,9,10,11]]].  Em <span id='citeF-8'></span>[[#cite-8|[8]]] utilizaram esquemas explícitos de diferenças finitas, na discretização das equações de águas rasas bidimensionais, adicionando condições de estabilidade e viscosidade artificial para suavizar as oscilações de alta frequência. Para validar, os esquemas discretos foram aplicados em um problema típico de engenharia hidráulica: abertura parcial de uma barragem e passagem de uma onda de inundação através de uma contração do canal, adotando a razão entre o nível de água a jusante <math> (h_r) </math> e à montante <math> (h_l ) </math> como sendo, <math> h_r / h_l = 0.5 </math>. Os resultados dos diferentes esquemas foram comparados e a ocorrência de oscilações numéricas foi observada, evidenciando a necessidade de se utilizar esquemas implícitos.  O estudo de caso apontado em <span id='citeF-8'></span>[[#cite-8|[8]]] foi analisado por <span id='citeF-12'></span>[[#cite-12|[12]]] usando um esquema de elementos finitos de alta ordem e uma melhoria na relação <math> h_r / h_l </math> para <math> 0,05 </math> foi obtida, buscando simular a condição de leito seco em uma determinada região do domínio. Em <span id='citeF-10'></span>[[#cite-10|[10]]] usaram dados experimentais em um evento de quebra de barragem para verificar a precisão, estabilidade e confiabilidade de uma solução numérica das equações de águas rasas obtidas por um esquema de volumes finitos bidimensional. O modelo numérico permitiu uma precisão de segunda ordem, tanto no espaço como no tempo, acarretando porém um alto custo computacional, apresentando algumas discrepâncias entre os resultados reais dos medidores e os resultados numéricos nos pontos pesquisados. Os autores destacam a necessidade de simular a propagação de uma onda de inundação para um fundo seco e áspero e otimizar o algoritmo para reduzir o tempo de computação.  Um esquema de volumes finitos não estruturado para malhas triangulares para simular fluxos bidimensionais instáveis em águas rasas com frentes molhadas/secas sobre topografia complexa é apresentado em <span id='citeF-9'></span>[[#cite-9|[9]]]. Utiliza-se em <span id='citeF-9'></span>[[#cite-9|[9]]] para o esquema temporal o método explícito de Runge-Kutta e para o espacial a malha não estruturada. Comparações satisfatórias foram obtidas com dados mensurados e numéricos. Todavia, os autores sugerem a incorporação de um procedimento de escalonamento num esquema implícito temporal, e uma implementação paralela do algoritmo, a fim de reduzir o tempo computacional.   A eficiência de um esquema implícito de volumes finitos é explorada em <span id='citeF-11'></span>[[#cite-11|[11]]] para simulação de fluidos em águas rasas bidimensionais, com frentes secas/molhadas. Adotou-se em <span id='citeF-11'></span>[[#cite-11|[11]]] um tratamento complexo para a discretização espacial, fazendo uso de uma malha não estruturada flexível para obter um melhor ajuste à geometrias irregulares. No entanto, é feito ressalvas quanto ao tempo computacional para obtenção da solução por esquemas implícitos, pois embora um método implícito normalmente exija menos etapas da solução, cada uma exige mais tempo computacional do que as soluções obtidas por um esquema explícito. Os autores apontam a necessidade de se utilizar um solucionador de matrizes que com um número reduzido de iterações alcance a convergência.   Neste contexto, o objetivo principal deste trabalho é formular um esquema híbrido e totalmente implícito de volumes finitos para discretização das equações de águas rasas, a fim de obter uma solução numérica estável com um algoritmo iterativo otimizado, com convergência rápida e baixo custo computacional, para aplicação na simulação de fluxos bidimensionais transitórios em topografias complexas (com declives e sujeitos a resistência ao fluxo). Condições iniciais descontínuas representando leitos úmidos/secos foram adotadas. Para tanto, o trabalho está organizado da seguinte forma: A seção [[#2 Modelo Matemático e Solução Numérica|2]] apresenta as equações hiperbólicas, que regem a natureza das águas rasas, considerando as perdas de atrito e declives da topografia do domínio computacional e detalhando a configuração numérica destas equações; na seção [[#3 O algoritmo|3]], o modelo é validado e uma comparação é feita com os resultados obtidos com o trabalho apresentado em <span id='citeF-8'></span>[[#cite-8|[8]]]. A robustez da solução numérica é verificada com uma extrapolação na relação <math> h_r / h_l </math> para <math>0.00002</math>, na seção [[#4 |4]]; o esquema é aplicado na análise de um problema hidráulico real  na seção [[#5 Verificação Numérica|5]] e; as conclusões são apresentadas na seção [[#6 Escoamento em um reservatório de contenção subterrâneo: uma aplicação do modelo completo|6]].
+<div id='img-1'></div>
+{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 40%;"
-==2 Modelo Matemático e Solução Numérica==
-Na análise do escoamento bidimensional transiente em um canal, cuja topografia do fundo é descrita por uma superfície <math>z_0(x,y)</math>, busca-se determinar o nível do fluido de superfície livre <math>h(x,y,t)</math> e o campo de velocidade do fluxo <math>U(u,v)</math> para cada instante de tempo <math>t</math>, conforme Fig. [[#img-1|1]].     <div id='img-1'></div>
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
 |-
-|[[Image:Draft_Coelho_632154385-fig1.png|300px|Notação adotada na análise de fluxo em águas rasas]]
+|style="padding:10px;text-align:center;"|[[Image:Draft_Coelho_632154385-fig1.png|300px|Notação adotada na análise de fluxo em águas rasas]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figure 1:''' Notação adotada na análise de fluxo em águas rasas
+| colspan="1" style="padding:10px;"| '''Figura 1'''. Notação adotada na análise de fluxo em águas rasas
 |}
 ===2.1 Sistema de Equações Governantes===
-O sistema de equações diferenciais que descreve o escoamento hidrostático de um fluido incompressível em águas rasas com superfície superior livre é obtido pela integração ao longo da profundidade das equações tridimensionais que modelam a quantidade de movimento e conservação de massa num escoamento, assumindo distribuição de velocidade uniforme na direção vertical. Como em <span id='citeF-3'></span>[[#cite-3|[3]]], adota-se as seguintes simplificações: (i) a distribuição de pressão na vertical é puramente hidrostática; (ii) não há contribuições laterais ao escoamento no reservatório analisado; (iii) a massa específica e a viscosidade do fluido no reservatório não sofrem variações significativas ao longo do tempo e do espaço; (iv) as tensões tangenciais são desprezíveis. Deste modo, as equações governantes do modelo hidrodinâmico, também ditas equações de águas rasas ou equações bidimensionais de Saint Venant, podem ser escritas na forma conservativa (com os fluxos dentro do sinal da derivada):
+O sistema de equações diferenciais que descreve o escoamento hidrostático de um fluido incompressível em águas rasas com superfície superior livre é obtido pela integração ao longo da profundidade das equações tridimensionais que modelam a quantidade de movimento e conservação de massa num escoamento, assumindo distribuição de velocidade uniforme na direção vertical. Como em Mahmood and Yevjevich <span id='citeF-3'></span>[[#cite-3|[3]]], adota-se as seguintes simplificações: (i) a distribuição de pressão na vertical é puramente hidrostática; (ii) não há contribuições laterais ao escoamento no reservatório analisado; (iii) a massa específica e a viscosidade do fluido no reservatório não sofrem variações significativas ao longo do tempo e do espaço; (iv) as tensões tangenciais são desprezíveis. Deste modo, as equações governantes do modelo hidrodinâmico, também ditas equações de águas rasas ou equações bidimensionais de Saint Venant, podem ser escritas na forma conservativa (com os fluxos dentro do sinal da derivada):
 <span id="eq-1"></span>
@@ Line 78: / Line 74: @@
 {| style="text-align: left; margin:auto;width: 100%;"
 |-
-| style="text-align: center;" | <math> \frac{\partial }{\partial t}(uh) +\frac{\partial }{\partial x}(u^2h+\frac{1}{2}gh^2)+\frac{\partial }{\partial y}( uvh)= gh(S_{0}^x - S_{f}^x),  </math>
+| style="text-align: center;" | <math> \frac{\partial }{\partial t}(uh) +\frac{\partial }{\partial x}\left(u^2h+\frac{1}{2}gh^2\right)+\frac{\partial }{\partial y}(uvh)= gh(S_{0}^x - S_{f}^x),  </math>
 |}
 | style="width: 5px;text-align: right;white-space: nowrap;" | (2)
@@ Line 124: / Line 120: @@
 ===2.2 Discretização do domínio e linearização das equações governantes===
-A solução numérica do modelo hidrodinâmico descrito pelas equações ([[#eq-1|1]]) - ([[#eq-3|3]]) é obtida por meio do método dos volumes finitos (MVF). Os fundamentos matemáticos e físicos do MVF são bem discutidos por <span id='citeF-13'></span>[[#cite-13|[13]]], <span id='citeF-14'></span>[[#cite-14|[14]]], <span id='citeF-15'></span>[[#cite-15|[15]]] e <span id='citeF-16'></span>[[#cite-16|[16]]].  De acordo com <span id='citeF-17'></span>[[#cite-17|[17]]], o MVF consiste na integração espacial e temporal das equações diferenciais na forma conservativa (com todos os fluxos envolvidos pelo operador derivada) em um volume genérico de controle.  O sistema de equações ([[#eq-1|1]]) - ([[#eq-3|3]]) é integrado espacialmente (sobre um volume de controle retangular) e no tempo (de <math>t</math> a <math>t + \Delta t</math>) e é empregado para a integração temporal a seguinte aproximação numérica:
+A solução numérica do modelo hidrodinâmico descrito pelas equações ([[#eq-1|1]])-([[#eq-3|3]]) é obtida por meio do método dos volumes finitos (MVF). Os fundamentos matemáticos e físicos do MVF são bem discutidos por <span id='citeF-16'></span>[[#cite-16|[16,17,18,19,20]]].  De acordo com <span id='citeF-20'></span>[[#cite-20|[20]]], o MVF consiste na integração espacial e temporal das equações diferenciais na forma conservativa (com todos os fluxos envolvidos pelo operador derivada) em um volume genérico de controle.  O sistema de equações ([[#eq-1|1]])-([[#eq-3|3]]) é integrado espacialmente (sobre um volume de controle retangular) e no tempo (de <math>t</math> a <math>t + \Delta t</math>) e é empregado para a integração temporal a seguinte aproximação numérica:
 <span id="eq-6"></span>
@@ Line 176: / Line 172: @@
 |}
-onde <math>\phi </math> representa  uma variável.  Na integração espacial é utilizado o Teorema Fundamental do Cálculo para reescrever as integrais das equações ([[#eq-1|1]]) - ([[#eq-3|3]]) como um sistema de equações constituído pelos valores das variáveis avaliadas nos pontos de integração <math>w, e, s, n </math>, que por sua vez, são avaliados por funções de interpolação. Para o cálculo destes valores, adota-se uma combinação dos esquemas de diferenciação central de <math>2^a</math> ordem ''(CDS - Central Differencing Scheme)'' e um esquema ''upwind (UDS - Upwind Differencing Scheme)'' de <math>1^a</math> ordem.  Na interpolação central, os valores das variáveis nas faces do volume de controle são dados pelo valor médio dos centróides vizinhos, isto é,
+onde <math>\phi </math> representa  uma variável.  Na integração espacial é utilizado o Teorema Fundamental do Cálculo para reescrever as integrais das equações ([[#eq-1|1]])-([[#eq-3|3]]) como um sistema de equações constituído pelos valores das variáveis avaliadas nos pontos de integração <math>w, e, s, n </math>, que por sua vez, são avaliados por funções de interpolação. Para o cálculo destes valores, adota-se uma combinação dos esquemas de diferenciação central de <math>2^a</math> ordem ''(CDS - Central Differencing Scheme)'' e um esquema ''upwind (UDS - Upwind Differencing Scheme)'' de <math>1^a</math> ordem.  Na interpolação central, os valores das variáveis nas faces do volume de controle são dados pelo valor médio dos centróides vizinhos, isto é,
 <span id="eq-10"></span>
@@ Line 374: / Line 370: @@
 {| style="text-align: left; margin:auto;width: 100%;"
 |-
-| style="text-align: center;" | <math> \left[(uh)_P - (uh)_P^0\right]\frac{\Delta x \Delta y}{\Delta t}+ \left(\Delta y\left[(u^2h+\frac{1}{2}gh^2)_e- (u^2h+\frac{1}{2}gh^2)_w\right]+\Delta x\left[(uvh)_n - (uvh)_s\right]\right)_{t+\Delta t}. </math>
+| style="text-align: center;" | <math> \left[(uh)_P - (uh)_P^0\right]\frac{\Delta x \Delta y}{\Delta t}+ \left(\Delta y\left[\left(u^2h+\frac{1}{2}gh^2\right)_e- \left(u^2h+\frac{1}{2}gh^2\right)_w\right]+\Delta x\left[(uvh)_n - (uvh)_s\right]\right)_{t+\Delta t}. </math>
 |}
 | style="width: 5px;text-align: right;white-space: nowrap;" | (26)
@@ Line 527: / Line 523: @@
 onde <math>A^\phi </math> representa a matriz de coeficientes da variável <math>\phi </math> (<math>h,u</math> ou <math>v</math>), enquanto <math>B^\phi </math> é o vetor que contém as condições de contorno e termos fontes.
-==3 O algoritmo==
+==3. O algoritmo==
 O código de simulação computacional foi programado em Linguagem Fortran, compilado  na versão Fortran Power Station 4.0 licenciada pela Microsoft Corporation. A escolha por esse programa deve-se à sua velocidade de execução e alta capacidade de armazenamento de dados.    Neste trabalho, o domínio físico retangular <math> D = [x_a, x_b] \times [y_a, y_b] </math> é dividido em <math> N_x </math> e <math> N_y </math> volumes de controle com comprimentos dados por <math> \Delta x, \Delta y </math> nas direções <math> x </math> e <math> y </math>, respectivamente. Da mesma forma, o espaço temporal <math> [t_0, t_f] </math> é dividido em <math> N_t </math> passos de tempo uniformes iguais a <math> \Delta t </math>.    A malha cartesiana ordenada e estruturada é varrida de acordo com o par <math>(i, j)</math> da esquerda para a direita aumentando o comprimento <math>(i - 1)\Delta x</math> na direção <math>x</math>, e de baixo para cima aumentando o comprimento <math>(j - 1) \Delta y</math> na direção <math>y</math>, começando no centróide <math>(1, 1)</math> que representa o ponto <math>(x_a + \Delta x/2, y_a + \Delta y/2)</math>. Desta forma, faz-se a seguinte correspondência:,
@@ Line 544: / Line 540: @@
 <span id="eq-40"></span>
-{| class="formulaSCP" style="width: 100%; text-align: left;"
+ {| class="formulaSCP" style="width: 100%; text-align: left;"
 |-
 |
@@ Line 568: / Line 564: @@
 onde <math>d_\phi </math> representa a diferença absoluta entre os dois lados da equação linearizada para a variável <math>\phi </math>.      O algoritmo completo é apresentado abaixo. O código também permite usar a última solução para inicializar uma nova execução com novas condições de contorno, se necessário.
-{| style="margin: 1em auto;border: 1px solid darkgray;"
+{| style="margin: 1em auto;border: 1px solid darkgray;font-size:90%;"
 |-
 |
-:   Declare as matrizes principais <math>h, u</math> e <math>v</math> e as variáveis secundárias <math>h_0, v_0</math> e <math>u0.</math>   Leia os dados do modelo: tamanho da malha espacial e temporal, coeficiente de atrito.   Calcule o comprimento dos passos espaciais e temporal   Faça a leitura da topografia <math>(z_0).</math>   Leia as condições iniciais <math>h_0,u_0</math> e <math>v_0.</math>   Calcule em cada passa de tempo as variáveis principais <math>h,u,v</math> <math> 1 \leq k \leq N_t</math>    Calcular as variáveis no contorno <math> 2\leq i \leq N_{x} -1 </math> e <math> 2\leq j \leq N_{y} -1 </math>     Interpolar <math>h,u,v,</math> e <math>z_0</math> nos pontos de integração nas faces do volume de controle    Calcular os coeficientes do sistema linear de <math>h</math>    Calcular o resíduo <math>r_h</math> <math>h(i,j)> z(i,j)</math>     Calcular <math>h(i,j)</math>    Calcular a coluna de água <math>H(i,j)=h(i,j)-z(i,j)</math>    Atualizar os valores de <math>h</math> nas faces dos volumes de controle    Calcular os coeficientes do sistema linear de <math>u</math>    Calcular o resíduo <math>r_u</math>    Calcular <math>u(i,j)</math>    Atualizar os valores de <math>u</math> nas faces dos volumes de controle    Calcular os os coeficientes do sistema linear de <math>v</math>    Calcular o resíduo <math>r_v</math>    Calcular <math>v(i,j)</math>        Faça <math>h(i,j)=z(i,j)</math> e <math>r_h,u,v</math> iguais a zero <math>\,\,\,</math> (o nível da água é calculado até <math>h=z_0</math>.) <math>\,\,\,</math> (fim da malha varrida)    Teste de Convergência     Atualização de <math>u,v</math> e <math>h</math> no tempo <math>(u_0=u,v_0=v,h_0=h)</math> <math>\,\,\,</math> (fim do ''loop'' temporal)        Faça o pós-processamento (tabelas e gráficos)
+:
+Declare as matrizes principais <math>h, u</math> e <math>v</math> e as variáveis secundárias <math>h_0, v_0</math> e <math>u0.</math>
+Leia os dados do modelo: tamanho da malha espacial e temporal, coeficiente de atrito.
+Calcule o comprimento dos passos espaciais e temporal
+Faça a leitura da topografia <math>(z_0).</math>
+Leia as condições iniciais <math>h_0,u_0</math> e <math>v_0.</math>
+Calcule em cada passa de tempo as variáveis principais <math>h,u,v</math>
+Se <math> 1 \leq k \leq N_t</math>  então
+Calcular as variáveis no contorno
+Se <math> 2\leq i \leq N_{x} -1 </math> e <math> 2\leq j \leq N_{y} -1 </math> então
+Interpolar <math>h,u,v,</math> e <math>z_0</math> nos pontos de integração nas faces do volume de controle
+Calcular os coeficientes do sistema linear de <math>h</math>
+Calcular o resíduo <math>r_h</math>
+Se <math>h(i,j)> z(i,j)</math> então
+Calcular <math>h(i,j)</math>
+Calcular a coluna de água <math>H(i,j)=h(i,j)-z(i,j)</math>
+Atualizar os valores de <math>h</math> nas faces dos volumes de controle
+Calcular os coeficientes do sistema linear de <math>u</math>
+Calcular o resíduo <math>r_u</math>
+Calcular <math>u(i,j)</math>
+Atualizar os valores de <math>u</math> nas faces dos volumes de controle
+Calcular os os coeficientes do sistema linear de <math>v</math>
+Calcular o resíduo <math>r_v</math>
+Calcular <math>v(i,j)</math>
+<math>\,\,\,</math> caso contrário
+Faça <math>h(i,j)=z(i,j)</math> e <math>r_h,u,v</math> iguais a zero
+<math>\,\,\,</math> (o nível da água é calculado até <math>h=z_0</math>.)
+<math>\,\,\,</math> (fim da malha varrida)
+Teste de Convergência
+Atualização de <math>u,v</math> e <math>h</math> no tempo <math>(u_0=u,v_0=v,h_0=h)</math>
+<math>\,\,\,</math> (fim do ''loop'' temporal)
+Faça o pós-processamento (tabelas e gráficos)
 |}
-==4 ==
+==4. Validação da solução numérica ==
-Validação da solução numérica    Com o objetivo de validar a abordagem numérica apresentada na seção anterior, esta seção apresenta uma comparação com os resultados obtidos para um modelo de abertura instantânea de comportas apresentado em <span id='citeF-8'></span>[[#cite-8|[8]]]. Este problema consiste em um domínio quadrado com <math>200\,m </math> de comprimento com uma barragem interna com <math>10\,m </math> de largura, disposta na faixa de <math> 90-100\,m</math>, que divide o domínio como é mostrado na Fig. ([[#img-2|2]]) .  Para simulação da abertura de comportas, considera-se uma abertura não simétrica e instantânea com <math>75\,m</math> de largura entre <math>30</math> e <math>105\,m</math>. Essa abertura causa uma liberação sutil da coluna de água a jusante, resultando em um fluxo com gradientes espaciais e temporais acentuados.     <div id='img-2'></div>
+Com o objetivo de validar a abordagem numérica apresentada na seção anterior, esta seção apresenta uma comparação com os resultados obtidos para um modelo de abertura instantânea de comportas apresentado em <span id='citeF-10'></span>[[#cite-10|[10]]]. Este problema consiste em um domínio quadrado com <math>200\,m </math> de comprimento com uma barragem interna com <math>10\,m </math> de largura, disposta na faixa de <math> 90-100\,m</math>, que divide o domínio como é mostrado na  [[#img-2|Figura 2]].  Para simulação da abertura de comportas, considera-se uma abertura não simétrica e instantânea com <math>75\,m</math> de largura entre <math>30</math> e <math>105\,m</math>. Essa abertura causa uma liberação sutil da coluna de água a jusante, resultando em um fluxo com gradientes espaciais e temporais acentuados.
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+<div id='img-2'></div>
+{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 50%;"
 |-
-|[[Image:Draft_Coelho_632154385-fig2.png|396px|Domínio do modelo]]
+|style="padding:10px;"|[[Image:Draft_Coelho_632154385-fig2.png|496px|Domínio do modelo]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figure 2:''' Domínio do modelo
+| colspan="1" style="padding-bottom:10px;" | '''Figura 2'''. Domínio do modelo
 |}
-Todos os parâmetros numéricos e físicos são mostrados na Tabela ([[#table-1|1]]). De acordo com <span id='citeF-8'></span>[[#cite-8|[8]]], neste problema a resistência ao escoamento é desprezada e sua parte inferior é plana.
-{|  class="floating_tableSCP wikitable" style="text-align: center; margin: 1em auto;min-width:50%;"
-|+ style="font-size: 75%;" |<span id='table-1'></span>Table. 1 Parâmetros para simulação de abertura da comporta
+Todos os parâmetros numéricos e físicos são mostrados na Tabela ([[#table-1|1]]). De acordo com <span id='citeF-10'></span>[[#cite-10|[10]]], neste problema a resistência ao escoamento é desprezada e sua parte inferior é plana.
-|- style="border-top: 2px solid;"
-|     Tamanho da malha
+<div class="center" style="font-size: 75%;">'''Tabela 1'''. Parâmetros para simulação de abertura da comporta</div>
-|  Tamanho dos elementos
-| <math>z_0</math>
+<div id='table-1'></div>
-|  Passo de tempo
+{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"
-|  Coeficiente de Manning
+|-style="text-align:center"
-| <math>g</math>
+!     Tamanho da malha  !!  Tamanho dos elementos !!  <math>z_0</math> !!  Passo de tempo !! Coeficiente de Manning  !!  <math>g</math>
-|- style="border-top: 2px solid;border-bottom: 2px solid;"
+|- style="border-bottom: 1px solid;text-align:center;"
 | <math>200\,m\,\times \,200\,m</math>
 | <math>5\,m\,\times \,5\,m</math>
@@ Line 602: / Line 660: @@
 | <math>0.0</math>
 | <math>9.81\,m/s^2</math>
 |}
-===4.1 Condições Iniciais===
+===4.1 Condições iniciais===
 Para o tempo inicial, o fluido é quiescente, ou seja, <math> u = v = 0</math> para todo o domínio, enquanto a coluna d'água inicial é dada por:
@@ Line 624: / Line 681: @@
 ===4.2 Condições de Contorno===
-Por uma questão de clareza, esta seção é apresentada aqui para mostrar de forma prática como as condições de contorno são tratadas neste trabalho usando para isso, o caso de validação como exemplo. As abordagens apresentadas aqui também são aplicadas a outros casos de maneira similar, respeitando, é claro, a configuração geométrica de cada um.    Como abordagem teórica, existem pelo menos três maneiras de lidar com as condições de contorno: realizar um balanço para cada volume de controle na fronteira, chamado simplesmente de equilíbrio de contorno; considerar um meio volume na fronteira com seu centróide pertencente ao próprio limite; e, por último, considerar volumes fictícios em torno dos limites, isso significa que o valor de limite prescrito deve ser igual à interpolação (qualquer que seja) entre os centróides do volume real de controle e o novo volume de controle fictício. Esta última abordagem é preferida, pois a primeira implica um tratamento matemático distinto para cada fronteira, o que torna o código mais complexo de acordo com a geometria. O caminho de <math>2^a</math> ordem, o  meio volume, não é conservador porque o valor é imposto, em vez disso, como resultado de um balanço de fluxo.     A abordagem de volume fictício permite que toda a malha seja constituída por volumes inteiros, o que é fácil de codificar e tem a mesma natureza conservadora para todos os volumes. Embora este artifício aumente o número de incógnitas do problema, acarretando custo computacional, é um procedimento útil com codificação fácil. Na Fig. ([[#img-3|3]]) é apresentado um esboço deste conceito.     <div id='img-3'></div>
+Por uma questão de clareza, esta seção é apresentada aqui para mostrar de forma prática como as condições de contorno são tratadas neste trabalho usando para isso, o caso de validação como exemplo. As abordagens apresentadas aqui também são aplicadas a outros casos de maneira similar, respeitando, é claro, a configuração geométrica de cada um.    Como apontado em <span id='citeF-21'></span>[[#cite-21|[21]]], existem pelo menos três maneiras de lidar com as condições de contorno: realizar um balanço para cada volume de controle na fronteira, chamado simplesmente de equilíbrio de contorno; considerar um meio volume na fronteira com seu centróide pertencente ao próprio limite; e, por último, considerar volumes fictícios em torno dos limites, isso significa que o valor de limite prescrito deve ser igual à interpolação (qualquer que seja) entre os centróides do volume real de controle e o novo volume de controle fictício. Esta última abordagem é preferida, pois a primeira implica um tratamento matemático distinto para cada fronteira, o que torna o código mais complexo de acordo com a geometria. O caminho de <math>2^a</math> ordem, o  meio volume, não é conservador porque o valor é imposto, em vez disso, como resultado de um balanço de fluxo.     A abordagem de volume fictício permite que toda a malha seja constituída por volumes inteiros, o que é fácil de codificar e tem a mesma natureza conservadora para todos os volumes. Embora este artifício aumente o número de incógnitas do problema, acarretando custo computacional, é um procedimento útil com codificação fácil. Na  [[#img-3|Figura 3]] é apresentado um esboço deste conceito.
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+<div id='img-3'></div>
+{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 40%;"
 |-
-|[[Image:Draft_Coelho_632154385-fig3.png|330px|Volumes fictícios ao redor da malha real]]
+|style="padding:10px;"|[[Image:Draft_Coelho_632154385-fig3.png|430px|Volumes fictícios ao redor da malha real]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figure 3:''' Volumes fictícios ao redor da malha real
+| colspan="1" style="padding-bottom:10px;"| '''Figura 3'''. Volumes fictícios ao redor da malha real
 |}
 Utiliza-se a Condição de Neumann nula para o nível da água nas paredes externas e nos contornos da represa. Desse modo,
@@ Line 655: / Line 716: @@
 |}
-onde <math>w</math> e <math>e</math> são os pontos dos contornos esquerdo e direito, respectivamente; e <math>s</math> e <math>n</math>  são os pontos dos contornos inferior e superior, respectivamente, como descrito na  Fig. ([[#img-4|4]]).     <div id='img-4'></div>
+onde <math>w</math> e <math>e</math> são os pontos dos contornos esquerdo e direito, respectivamente; e <math>s</math> e <math>n</math>  são os pontos dos contornos inferior e superior, respectivamente, como descrito na   [[#img-4|Figura 4]].
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+<div id='img-4'></div>
+{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 50%;"
 |-
-|[[Image:Draft_Coelho_632154385-fig4.png|360px|Esquema da condição de contorno para o nível h do fluido]]
+|style="padding:10px;"|[[Image:Draft_Coelho_632154385-fig4.png|460px|Esquema da condição de contorno para o nível h do fluido]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figure 4:''' Esquema da condição de contorno para o nível <math>h</math> do fluido
+| colspan="1" style="padding-bottom:10px;"| '''Figura 4'''. Esquema da condição de contorno para o nível <math>h</math> do fluido
 |}
-Todas as derivadas nas faces fronteiras (<math>w,e,n,</math> e <math>s</math>) podem ser interpoladas   como funções dos centróides vizinhos destacados na Fig. 4: um externo fictício (<math>W, E, S</math> ou <math>N</math>) e um interno <math>(P)</math> pertencente ao domínio.  Consequentemente,
+Todas as derivadas nas faces fronteiras (<math>w,e,n,</math> e <math>s</math>) podem ser interpoladas   como funções dos centróides vizinhos destacados na Figura 4: um externo fictício (<math>W, E, S</math> ou <math>N</math>) e um interno <math>(P)</math> pertencente ao domínio.  Consequentemente,
 <span id="eq-45"></span>
@@ Line 863: / Line 928: @@
 ===4.3 Resultados===
-O algoritmo proposto e a formulação numérica devem ser validados para determinar sua precisão e robustez. Isso é feito por meio da aplicação do modelo e do código numérico a um problema representativo de águas rasas: o problema da abertura instantânea de comportas, amplamente utilizado para testar códigos numéricos ou modelos matemáticos que tratam de vazões superficiais sujeitas a mudanças abruptas, mudanças no tempo e no espaço, e a abertura instantânea de comportas tem esses detalhes físicos e é amplamente estudada na literatura [Fennema1990,Adriano2018,Rezende:2015, Hirt:1981, Jahanbakhsh2007,Villota2017].
+O algoritmo proposto e a formulação numérica devem ser validados para determinar sua precisão e robustez. Isso é feito por meio da aplicação do modelo e do código numérico a um problema representativo de águas rasas: o problema da abertura instantânea de comportas, amplamente utilizado para testar códigos numéricos ou modelos matemáticos que tratam de vazões superficiais sujeitas a mudanças abruptas, mudanças no tempo e no espaço, e a abertura instantânea de comportas tem esses detalhes físicos e é amplamente estudada na literatura <span id='citeF-10'></span><span id='citeF-14'></span><span id='citeF-22'></span><span id='citeF-23'><span id='citeF-24'><span id='citeF-25'></span>[[#cite-22|[10,14,22,23,24,25]]].
 ====4.3.1 Nível da água====
-A Fig. ([[#img-5|5]]) mostra a superfície livre de água após <math>7.1</math> segundos após a abertura da comporta. Os resultados são comparados com aqueles obtidos em <span id='citeF-8'></span>[[#cite-8|[8]]]. Observa-se na Fig. ([[#img-5|5]])  uma frente de onda após a abertura da eclusa  e uma onda de depressão no lado esquerdo à barreira. Os níveis de superfície computados estão de acordo com os apresentados por <span id='citeF-8'></span>[[#cite-8|[8]]] para o mesmo instante de tempo. <div id='img-5a'></div>
+A Figura [[#img-5|5]] mostra a superfície livre de água após <math>7.1</math> segundos após a abertura da comporta. Os resultados são comparados com aqueles obtidos em <span id='citeF-10'></span>[[#cite-10|[10]]]. Observa-se na  [[#img-5|Figura 5]]  uma frente de onda após a abertura da eclusa  e uma onda de depressão no lado esquerdo à barreira. Os níveis de superfície computados estão de acordo com os apresentados por <span id='citeF-10'></span>[[#cite-10|[10]]] para o mesmo instante de tempo.
+<div id='img-5a'></div>
 <div id='img-5b'></div>
 <div id='img-5'></div>
 {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
 |-
-|[[Image:Draft_Coelho_632154385-fig5a.png|534px|Resultado obtido pelo esquema híbrido proposto]]
+|style="padding:10px;"|[[Image:Draft_Coelho_632154385-fig5a.png|534px|Resultado obtido pelo esquema híbrido proposto]]
-|[[Image:Draft_Coelho_632154385-fig5b.png|496px|Resultado obtido em <span id='citeF-8'></span>[[#cite-8|[8]]]]]
+|style="padding:10px;"|[[Image:Draft_Coelho_632154385-fig5b.png|496px|Resultado obtido em <span id='citeF-8'></span>[[#cite-8|[8]]]]]
 |- style="text-align: center; font-size: 75%;"
 | (a) Resultado obtido pelo esquema híbrido proposto
-| (b) Resultado obtido em <span id='citeF-8'></span>[[#cite-8|[8]]]
+| (b) Resultado obtido em <span id='citeF-10'></span>[[#cite-10|[10]]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="2" | '''Figure 5:''' Comparação do comportamento da superfície livre <math>3D</math> após <math>7.1 \, s </math> após a abertura da comporta
+| colspan="2" style="padding-bottom:10px;padding-top:10px;"| '''Figura 5'''. Comparação do comportamento da superfície livre <math>3D</math> após <math>7.1 \, s </math> após a abertura da comporta
 |}
-A Tabela ([[#4.3.1 Nível da água|4.3.1]]) mostra os valores de nível para algumas posições do domínio para fins de comparação quantitativa.
-{|  class="floating_tableSCP wikitable" style="text-align: center; margin: 1em auto;min-width:50%;"
-|+ style="font-size: 75%;" |<span id='table-2'></span>Table. 2 Valores de nível de água  usando a relação <math>h_r/h_l=0.5</math>
+A  [[#table-2|Tabela 2]] mostra os valores de nível para algumas posições do domínio para fins de comparação quantitativa.
-|- style="border-top: 2px solid;"
-| style="border-right: 2px solid;" |         yx
+<div class="center" style="font-size: 75%;">'''Tabela 2'''. Valores de nível de água  usando a relação <math>h_r/h_l=0.5</math></div>
-| style="border-left: 2px solid;" | <math>0</math>
-| <math>25</math>
+<div id='table-1'></div>
-| <math>50</math>
+{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"
-| <math>75</math>
+|-style="text-align:center"
-| <math>100</math>
+!   yx !!  <math>0</math> !!  <math>25</math> !!  <math>50</math> !!  <math>75</math> !!  <math>100</math> !!  <math>125</math> !!  <math>150</math> !!  <math>175</math> !!  <math>200</math>
-| <math>125</math>
+|- style="text-align:center;"
-| <math>150</math>
+|  <math>0</math>
-| <math>175</math>
+|  <math>9.99932</math>
-| <math>200</math>
-|- style="border-top: 2px solid;"
-| style="border-right: 2px solid;" | <math>0</math>
-| style="border-left: 2px solid;" | <math>9.99932</math>
 | <math>9.96352</math>
 | <math>9.43066</math>
@@ Line 905: / Line 968: @@
 | <math>5.00090</math>
 | <math>5.0</math>
-|-
+|-style="text-align:center;"
-| style="border-right: 2px solid;" | <math>25</math>
+|  <math>25</math>
-| style="border-left: 2px solid;" | <math>9.99497</math>
+| <math>9.99497</math>
 | <math>9.88268</math>
 | <math>9.18459</math>
@@ Line 916: / Line 979: @@
 | <math>5.01876</math>
 | <math>5.0</math>
-|-
+|-style="text-align:center;"
-| style="border-right: 2px solid;" | <math>50</math>
+|  <math>50</math>
-| style="border-left: 2px solid;" | <math>9.98848</math>
+|  <math>9.98848</math>
 | <math>9.77139</math>
 | <math>8.67804</math>
@@ Line 927: / Line 990: @@
 | <math>5.08294</math>
 | <math>5.0</math>
-|-
+|-style="text-align:center;"
-| style="border-right: 2px solid;" | <math>75</math>
+|  <math>75</math>
-| style="border-left: 2px solid;" | <math>9.98760</math>
+|  <math>9.98760</math>
 | <math>9.75073</math>
 | <math>8.51860</math>
@@ Line 938: / Line 1,001: @@
 | <math>5.09955</math>
 | <math>5.0</math>
-|-
+|-style="text-align:center;"
-| style="border-right: 2px solid;" | <math>100</math>
+|  <math>100</math>
-| style="border-left: 2px solid;" | <math>9.99177</math>
+|  <math>9.99177</math>
 | <math>9.83175</math>
 | <math>8.98309</math>
@@ Line 949: / Line 1,012: @@
 | <math>5.04260</math>
 | <math>5.0</math>
-|-
+|-style="text-align:center;"
-| style="border-right: 2px solid;" | <math>125</math>
+|  <math>125</math>
-| style="border-left: 2px solid;" | <math>9.99851</math>
+|  <math>9.99851</math>
 | <math>9.95297</math>
 | <math>9.51005</math>
@@ Line 960: / Line 1,023: @@
 | <math>5.00278</math>
 | <math>5.0</math>
-|-
+|-style="text-align:center;"
-| style="border-right: 2px solid;" | <math>150</math>
+|  <math>150</math>
-| style="border-left: 2px solid;" | <math>9.99996</math>
+|  <math>9.99996</math>
 | <math>9.99796</math>
 | <math>9.93089</math>
@@ Line 971: / Line 1,034: @@
 | <math>5.00001</math>
 | <math>5.0</math>
-|-
+|-style="text-align:center;"
-| style="border-right: 2px solid;" | <math>175</math>
+|  <math>175</math>
-| style="border-left: 2px solid;" | <math>9.99999</math>
+|  <math>9.99999</math>
 | <math>9.99997</math>
 | <math>9.99829</math>
@@ Line 982: / Line 1,045: @@
 | <math>5.00000</math>
 | <math>5.0</math>
-|- style="border-bottom: 2px solid;"
+|- style="text-align:center;"
-| style="border-right: 2px solid;" | <math>200</math>
+| <math>200</math>
-| style="border-left: 2px solid;" | <math>10.0000</math>
+|  <math>10.0000</math>
 | <math>9.99999</math>
 | <math>9.99999</math>
@@ Line 993: / Line 1,056: @@
 | <math>5.00000</math>
 | <math>5.0</math>
 |}
-Duas linhas de amostra foram extraídas em <math>7.1</math> segundos antes da abertura da comporta <math>(x = 75\, m) </math> e após a abertura <math> (x = 105\, m) </math> e os dados foram comparados a <span id='citeF-8'></span>[[#cite-8|[8]]], onde são confrontados outros dois esquemas numéricos: MacCormack e Gabutti (descritos em <span id='citeF-22'></span>[[#cite-22|[22]]] e <span id='citeF-23'></span>[[#cite-23|[23]]], respectivamente), e ambos os métodos são de precisão de segunda ordem e requerem um cálculo mais complexo, uma vez que uma etapa do preditor deve ser avaliada em <math> t + \Delta t/2.</math>  A Fig. ([[#img-6|6]]) apresenta a comparação dos perfis transversais do nível da água após <math>7.1</math> segundos da abertura da barreira interna, e os dados estão em boa concordância com <span id='citeF-8'></span>[[#cite-8|[8]]]. Algumas diferenças podem ser observadas antes e depois da abertura da barragem, principalmente com a onda de depressão antes do portão, mas os níveis estão na mesma ordem e com o mesmo comportamento físico. Além disso, observa-se algumas elevações locais na frente de onda. <div id='img-6'></div>
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+Duas linhas de amostra foram extraídas em <math>7.1</math> segundos antes da abertura da comporta <math>(x = 75\, m) </math> e após a abertura <math> (x = 105\, m) </math> e os dados foram comparados a <span id='citeF-10'></span>[[#cite-10|[10]]], onde são confrontados outros dois esquemas numéricos: MacCormack e Gabutti (descritos em <span id='citeF-26'></span>[[#cite-26|[26]]] e <span id='citeF-27'></span>[[#cite-27|[27]]], respectivamente), e ambos os métodos são de precisão de segunda ordem e requerem um cálculo mais complexo, uma vez que uma etapa do preditor deve ser avaliada em <math> t + \Delta t/2.</math>  A  [[#img-6|Figura 6]] apresenta a comparação dos perfis transversais do nível da água após <math>7.1</math> segundos da abertura da barreira interna, e os dados estão em boa concordância com <span id='citeF-10'></span>[[#cite-10|[10]]]. Algumas diferenças podem ser observadas antes e depois da abertura da barragem, principalmente com a onda de depressão antes do portão, mas os níveis estão na mesma ordem e com o mesmo comportamento físico. Além disso, observa-se algumas elevações locais na frente de onda.
+<div id='img-6'></div>
+{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 45%;"
 |-
-|[[Image:Draft_Coelho_632154385-fig6.png|276px|Comparação do perfil do nível de água em duas linhas de amostra em x=115\,m e x=75\,m]]
+|style="padding:10px;"|[[Image:Draft_Coelho_632154385-fig6.png|476px|Comparação do perfil do nível de água em duas linhas de amostra em x=115\,m e x=75\,m]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figure 6:''' Comparação do perfil do nível de água em duas linhas de amostra em <math>x=115\,m</math> e <math>x=75\,m</math>
+| colspan="1" style="padding-bottom:10px;"| '''Figura 6'''. Comparação do perfil do nível de água em duas linhas de amostra em <math>x=115\,m</math> e <math>x=75\,m</math>
 |}
-A Fig. ([[#img-7|7]]) mostra os perfis longitudinais em <math> y = 70 \, m </math> e <math> y = 150 \, m </math> a <math> 7.1 </math> segundos. A linha <math> h (x, 150) </math> está quebrada porque intercepta a parede da barragem, enquanto a linha <math> h (x, 75) </math> passa pela abertura da comporta. Pode ser visto na linha <math> h (x, 150) </math> que a onda frontal após a abertura da eclusa mostra uma elevação de pico que é devida à diferença de velocidade entre a frente de onda e a onda de depressão em movimento a jusante. A frente de onda é desacelerada pelo corpo de água em repouso na frente dela. Embora esse comportamento não seja mostrado nos resultados de <span id='citeF-8'></span>[[#cite-8|[8]]], a posição da frente de onda e o início da onda de depressão são coincidentes. O perfil <math>h (x, 150)</math> tem uma melhor correspondência, uma vez que nesta posição o comportamento do fluido é mais suave e menos dinâmico devido à presença da parede. É evidente que o local mais sensível a ser resolvido é a região de abertura da comporta devido ao abrupto colapso da coluna de água.  <div id='img-7'></div>
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+A  [[#img-7|Figura 7]] mostra os perfis longitudinais em <math> y = 70 \, m </math> e <math> y = 150 \, m </math> a <math> 7.1 </math> segundos. A linha <math> h (x, 150) </math> está quebrada porque intercepta a parede da barragem, enquanto a linha <math> h (x, 75) </math> passa pela abertura da comporta. Pode ser visto na linha <math> h (x, 150) </math> que a onda frontal após a abertura da eclusa mostra uma elevação de pico que é devida à diferença de velocidade entre a frente de onda e a onda de depressão em movimento a jusante. A frente de onda é desacelerada pelo corpo de água em repouso na frente dela. Embora esse comportamento não seja mostrado nos resultados de <span id='citeF-10'></span>[[#cite-10|[10]]], a posição da frente de onda e o início da onda de depressão são coincidentes. O perfil <math>h (x, 150)</math> tem uma melhor correspondência, uma vez que nesta posição o comportamento do fluido é mais suave e menos dinâmico devido à presença da parede. É evidente que o local mais sensível a ser resolvido é a região de abertura da comporta devido ao abrupto colapso da coluna de água.
+<div id='img-7'></div>
+{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 45%;"
 |-
-|[[Image:Draft_Coelho_632154385-fig7.png|282px|Comparação das linhas de amostra longitudinais do nível de água em y=70\,m e y=150\,m]]
+|style="padding:10px;"|[[Image:Draft_Coelho_632154385-fig7.png|482px|Comparação das linhas de amostra longitudinais do nível de água em y=70\,m e y=150\,m]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figure 7:''' Comparação das linhas de amostra longitudinais do nível de água em <math>y=70\,m</math> e <math>y=150\,m</math>
+| colspan="1" style="padding-bottom:10px;"| '''Figura 7'''. Comparação das linhas de amostra longitudinais do nível de água em <math>y=70\,m</math> e <math>y=150\,m</math>
 |}
-As mudanças no nível da água podem ser monitoradas por um período de tempo, em um ponto fixo do domínio. Na Fig. ([[#img-8|8]]), pode-se ver a amostragem do nível de água em <math>P(115,70)</math>, localizado na frente da abertura da comporta que permite capturar o avanço da frente de onda, durante <math>7.1</math> segundos. <div id='img-8'></div>
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+As mudanças no nível da água podem ser monitoradas por um período de tempo, em um ponto fixo do domínio. Na  [[#img-8|Figura 8]], pode-se ver a amostragem do nível de água em <math>P(115,70)</math>, localizado na frente da abertura da comporta que permite capturar o avanço da frente de onda, durante <math>7.1</math> segundos.
+<div id='img-8'></div>
+{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 50%;"
 |-
-|[[Image:Draft_Coelho_632154385-fig8.png|276px|Comparação do tempo de amostragem do nível da água em P(115,70)]]
+|style="padding:10px;"|[[Image:Draft_Coelho_632154385-fig8.png|476px|Comparação do tempo de amostragem do nível da água em P(115,70)]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figure 8:''' Comparação do tempo de amostragem do nível da água em <math>P(115,70)</math>
+| colspan="1" style="padding:10px;"| '''Figura 8'''. Comparação do tempo de amostragem do nível da água em <math>P(115,70)</math>
 |}
-Pode ser visto em ambos os resultados que a frente de onda leva <math>1</math> segundo após a abertura para alcançar o ponto de teste <math>P</math>, e depois de <math>3</math> segundos, o nível da água está próximo de <math>7\,m </math> em ambos os casos. Novamente, há um comportamento mais dinâmico em comparação com os resultados de <span id='citeF-8'></span>[[#cite-8|[8]]], que é mais suave com oscilações brandas devido a utilização de termos dissipativos artificiais, enquanto os resultados obtidos pelo esquema híbrido proposto apresentam uma onda de pico que tende a ser atenuada com o tempo, de acordo com o fenômeno físico esperado.
+Pode ser visto em ambos os resultados que a frente de onda leva <math>1</math> segundo após a abertura para alcançar o ponto de teste <math>P</math>, e depois de <math>3</math> segundos, o nível da água está próximo de <math>7\,m </math> em ambos os casos. Novamente, há um comportamento mais dinâmico em comparação com os resultados de <span id='citeF-10'></span>[[#cite-10|[10]]], que é mais suave com oscilações brandas devido a utilização de termos dissipativos artificiais, enquanto os resultados obtidos pelo esquema híbrido proposto apresentam uma onda de pico que tende a ser atenuada com o tempo, de acordo com o fenômeno físico esperado.
 ====4.3.2 Campo de Velocidade====
-A Fig. ([[#img-9|9]]) compara os resultados de <span id='citeF-8'></span>[[#cite-8|[8]]] e deste trabalho para o campo de velocidade no instante de tempo igual a <math> 7.1 </math> segundos. Os vetores do campo de velocidade estão em muito boa concordância apresentando as mesmas direções, magnitudes (representadas pelos tamanhos das flechas) com um espalhamento pela mesma área ao redor da abertura da comporta.  <div id='img-9a'></div>
+A  [[#img-9|Figura 9]] compara os resultados de <span id='citeF-10'></span>[[#cite-10|[10]]] e deste trabalho para o campo de velocidade no instante de tempo igual a <math> 7.1 </math> segundos. Os vetores do campo de velocidade estão em muito boa concordância apresentando as mesmas direções, magnitudes (representadas pelos tamanhos das flechas) com um espalhamento pela mesma área ao redor da abertura da comporta.
+<div id='img-9a'></div>
 <div id='img-9b'></div>
 <div id='img-9'></div>
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;"
 |-
-|[[Image:Draft_Coelho_632154385-fig9a.png|471px|Campo de velocidade usando o esquema híbrido]]
+|style="padding:10px;"|[[File:Coelho_et_al_2019b_8647_571px-Draft_Coelho_632154385-fig9a.png|250px]]
-|[[Image:Draft_Coelho_632154385-fig9b.png|471px|Campo de velocidade apresentado em  <span id='citeF-8'></span>[[#cite-8|[8]]]]]
+|style="padding:10px;"|[[File:Coelho_et_al_2019b_3874_571px-Draft_Coelho_632154385-fig9b.png|250px]]
 |- style="text-align: center; font-size: 75%;"
 | (a) Campo de velocidade usando o esquema híbrido
-| (b) Campo de velocidade apresentado em  <span id='citeF-8'></span>[[#cite-8|[8]]]
+| (b) Campo de velocidade apresentado em  <span id='citeF-8'></span>[[#cite-10|[10]]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="2" | '''Figure 9:''' Comparação do campo de velocidade em torno da abertura da comporta em <math>7.1</math> segundos
+| colspan="2" style="padding-bottom:10px;padding-top:10px;"| '''Figura 9'''. Comparação do campo de velocidade em torno da abertura da comporta em <math>7.1</math> segundos
 |}
 ===4.4 Considerações===
-Apesar das oscilações locais ou da suavidade, os resultados estão em boa concordância geral com as frentes de onda e ondas de depressão de acordo com os dados da literatura, o que reflete a habilidade da abordagem numérica proposta em lidar e capturar o comportamento físico do fluxos de águas rasas que validam o modelo e a abordagem numérica.  É importante enfatizar que, em <span id='citeF-8'></span>[[#cite-8|[8]]] artifícios numéricos foram usados como uma viscosidade artificial para lidar com as oscilações da superfície livre, o que torna o resultado de <span id='citeF-8'></span>[[#cite-8|[8]]] mais difusivo e quase sem flutuações locais no nível da superfície. Por outro lado, neste trabalho, nenhuma abordagem similar foi empregada para lidar com uma possível instabilidade numérica devido à aproximação de tempo de alta ordem.    O esquema de tempo implícito proposto de <math>1^a</math> ordem neste trabalho mostrou-se capaz de capturar o comportamento físico e os gradientes de tempo rígidos sem etapas de tempo proibitivas. E a formulação ''UDS'' em associação com a abordagem ''CDS'' foi capaz de lidar com mudanças abruptas nos campos de variáveis sem instabilidades ou oscilações não físicas. Todas as oscilações de superfície livre têm um significado físico quando as diferenças de velocidade entre as porções distintas da onda são contabilizadas porque essas diferenças compactam a onda e aumentam sua amplitude. Como apontado em <span id='citeF-8'></span>[[#cite-8|[8]]], as oscilações numéricas puras não são amortecidas com o tempo, em vez disso, elas tendem a ser amplificadas com gradientes rígidos, o que faz com que o sistema divirja. Porém fazendo uso do esquema numérico híbrido aqui apresentado, a solução se mantêm suave com resíduos muito baixo, inferiores a tolerância máxima de <math>10^{-11}</math>. Mesmo para grandes intervalos de tempo, não há crescimento no valor destes resíduos, atingindo valores inferiores a tolerância máxima em poucas iterações.
+Apesar das oscilações locais ou da suavidade, os resultados estão em boa concordância geral com as frentes de onda e ondas de depressão de acordo com os dados da literatura, o que reflete a habilidade da abordagem numérica proposta em lidar e capturar o comportamento físico do fluxos de águas rasas que validam o modelo e a abordagem numérica.  É importante enfatizar que, em <span id='citeF-10'></span>[[#cite-10|[10]]] artifícios numéricos foram usados como uma viscosidade artificial para lidar com as oscilações da superfície livre, o que torna o resultado de <span id='citeF-10'></span>[[#cite-10|[10]]] mais difusivo e quase sem flutuações locais no nível da superfície. Por outro lado, neste trabalho, nenhuma abordagem similar foi empregada para lidar com uma possível instabilidade numérica devido à aproximação de tempo de alta ordem.    O esquema de tempo implícito proposto de <math>1^a</math> ordem neste trabalho mostrou-se capaz de capturar o comportamento físico e os gradientes de tempo rígidos sem etapas de tempo proibitivas. E a formulação ''UDS'' em associação com a abordagem ''CDS'' foi capaz de lidar com mudanças abruptas nos campos de variáveis sem instabilidades ou oscilações não físicas. Todas as oscilações de superfície livre têm um significado físico quando as diferenças de velocidade entre as porções distintas da onda são contabilizadas porque essas diferenças compactam a onda e aumentam sua amplitude. Como apontado em <span id='citeF-10'></span>[[#cite-10|[10]]], as oscilações numéricas puras não são amortecidas com o tempo, em vez disso, elas tendem a ser amplificadas com gradientes rígidos, o que faz com que o sistema divirja. Porém fazendo uso do esquema numérico híbrido aqui apresentado, a solução se mantêm suave com resíduos muito baixo, inferiores a tolerância máxima de <math>10^{-11}</math>. Mesmo para grandes intervalos de tempo, não há crescimento no valor destes resíduos, atingindo valores inferiores a tolerância máxima em poucas iterações.
-==5 Verificação Numérica==
+==5. Verificação numérica==
 Uma validação numérica assegura que o modelo é capaz de fornecer uma resposta física correta sob certas condições iniciais e de contorno. No entanto, mais testes são feitos para descobrir se, sob outro cenário, o modelo computacional é capaz de prever o comportamento físico do sistema. Esta etapa, juntamente com a validação anterior, permite definir se o modelo computacional é útil como uma ferramenta preditiva para qualquer fluxo de água superficial sob condições diversas.
@@ Line 1,044: / Line 1,122: @@
 ===5.1 Teste da condição de leito seco===
-Para avaliar a robustez do código, o modelo da seção [[#3 O algoritmo|3]] é considerado, fazendo a razão <math> h_r / h_l </math> tender a zero, isto é, <math> h_r </math> é considerado com uma fina película líquida, como uma ''“wet condition”''. Esta situação é geralmente uma situação física difícil de resolver devido às singularidades que podem surgir nos termos fontes de atrito, e o possível mau condicionamento da matriz de coeficientes devido a valores nulos que surgirão ao longo dos cálculos quando <math>h</math> tende a zero. Nesta situação, as equações de águas rasas tendem a enfraquecer seu significado físico. Em vista disso, <math> h_r </math> não pode ser exatamente zero, e o código deve ser capaz de lidar com essa situação. Na Fig.([[#img-10|10]]), pode-se observar o nível da água para os primeiros <math> 7.1 </math> segundos com <math> h_r / h_l = 0.00002 </math>.  <div id='img-10'></div>
+Para avaliar a robustez do código, o modelo da Seção [[#3 O algoritmo|3]] é considerado, fazendo a razão <math> h_r / h_l </math> tender a zero, isto é, <math> h_r </math> é considerado com uma fina película líquida, como uma ''“wet condition”''. Esta situação é geralmente uma situação física difícil de resolver devido às singularidades que podem surgir nos termos fontes de atrito, e o possível mau condicionamento da matriz de coeficientes devido a valores nulos que surgirão ao longo dos cálculos quando <math>h</math> tende a zero. Nesta situação, as equações de águas rasas tendem a enfraquecer seu significado físico. Em vista disso, <math> h_r </math> não pode ser exatamente zero, e o código deve ser capaz de lidar com essa situação. Na  [[#img-10|Figura 10]], pode-se observar o nível da água para os primeiros <math> 7.1 </math> segundos com <math> h_r / h_l = 0.00002 </math>.
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+<div id='img-10'></div>
+{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 50%;"
 |-
-|[[Image:Draft_Coelho_632154385-fig10.png|300px|Visão tridimensional da superfície livre em 7.1 segundos após a abertura da comporta com h<sub>r</sub>/hₗ=0.00002]]
+|style="padding-bottom:10px;"|[[Image:Draft_Coelho_632154385-fig10.png|450px|Visão tridimensional da superfície livre em 7.1 segundos após a abertura da comporta com h<sub>r</sub>/hₗ=0.00002]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figure 10:''' Visão tridimensional da superfície livre em <math>7.1</math> segundos após a abertura da comporta com <math>h_r/h_l=0.00002</math>
+| colspan="1" style="padding-bottom:10px;"| '''Figura 10'''. Visão tridimensional da superfície livre em <math>7.1</math> segundos após a abertura da comporta com <math>h_r/h_l=0.00002</math>
 |}
-Os valores dos níveis são reproduzidos na Tabela 3, com todas as unidades em metros.
-{|  class="floating_tableSCP wikitable" style="text-align: center; margin: 1em auto;min-width:50%;"
-|+ style="font-size: 75%;" |<span id='table-3'></span>Table. 3 Valores de nível de água usando <math>hr/hl=0.00002</math>
+Os valores dos níveis são reproduzidos na [[#tab-3|Tabela 3]], com todas as unidades em metros.
-|- style="border-top: 2px solid;"
-| style="border-right: 2px solid;" |     yx
+<div class="center" style="font-size: 75%;">'''Tabela 3'''. Valores de nível de água usando <math>hr/hl=0.00002</math></div>
-| style="border-left: 2px solid;" | <math>0</math>
-| <math>25</math>
+<div id='tab-3'></div>
-| <math>50</math>
+{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;"
-| <math>75</math>
+|-style="text-align:center"
-| <math>100</math>
+!  yx !!   <math>0</math>  !!  <math>25</math> !!  <math>50</math> !!  <math>75</math> !!  <math>100</math> !!  <math>125</math>   !!  <math>150</math> !!  <math>175</math> !!<math>200</math>
-| <math>125</math>
+|- style="text-align:center;"
-| <math>150</math>
+| <math>0</math>
-| <math>175</math>
+| <math>9.9984</math>
-| <math>200</math>
-|- style="border-top: 2px solid;"
-| style="border-right: 2px solid;" | <math>0</math>
-| style="border-left: 2px solid;" | <math>9.9984</math>
 | <math>9.9274</math>
 | <math>9.1889</math>
@@ Line 1,076: / Line 1,152: @@
 | <math>0.0002</math>
 | <math>0.0002</math>
-|-
+|-style="text-align:center;"
-| style="border-right: 2px solid;" | <math>25</math>
+|<math>25</math>
-| style="border-left: 2px solid;" | <math>9.9904</math>
+|  <math>9.9904</math>
 | <math>9.8118</math>
 | <math>8.9461</math>
@@ Line 1,087: / Line 1,163: @@
 | <math>0.0583</math>
 | <math>0.0002</math>
-|-
+|-style="text-align:center;"
-| style="border-right: 2px solid;" | <math>50</math>
+|  <math>50</math>
-| style="border-left: 2px solid;" | <math>9.9789</math>
+|  <math>9.9789</math>
 | <math>9.6532</math>
 | <math>8.3391</math>
@@ Line 1,098: / Line 1,174: @@
 | <math>0.4390</math>
 | <math>0.0002</math>
-|-
+|-style="text-align:center;"
-| style="border-right: 2px solid;" | <math>75</math>
+|  <math>75</math>
-| style="border-left: 2px solid;" | <math>9.9773</math>
+|  <math>9.9773</math>
 | <math>9.6221</math>
 | <math>8.1432</math>
@@ Line 1,109: / Line 1,185: @@
 | <math>0.6072</math>
 | <math>0.0002</math>
-|-
+|-style="text-align:center;"
-| style="border-right: 2px solid;" | <math>100</math>
+| <math>100</math>
-| style="border-left: 2px solid;" | <math>9.9848</math>
+|  <math>9.9848</math>
 | <math>9.7398</math>
 | <math>8.7161</math>
@@ Line 1,120: / Line 1,196: @@
 | <math>0.1380</math>
 | <math>0.0002</math>
-|-
+|-style="text-align:center;"
-| style="border-right: 2px solid;" | <math>125</math>
+|  <math>125</math>
-| style="border-left: 2px solid;" | <math>9.9969</math>
+|  <math>9.9969</math>
 | <math>9.9163</math>
 | <math>9.3453</math>
@@ Line 1,131: / Line 1,207: @@
 | <math>0.0058</math>
 | <math>0.0002</math>
-|-
+|-style="text-align:center;"
-| style="border-right: 2px solid;" | <math>150</math>
+| <math>150</math>
-| style="border-left: 2px solid;" | <math>9.9999</math>
+|  <math>9.9999</math>
 | <math>9.9950</math>
 | <math>9.8806</math>
@@ Line 1,142: / Line 1,218: @@
 | <math>0.0002</math>
 | <math>0.0002</math>
-|-
+|-style="text-align:center;"
-| style="border-right: 2px solid;" | <math>175</math>
+|  <math>175</math>
-| style="border-left: 2px solid;" | <math>9.9999</math>
+|  <math>9.9999</math>
 | <math>9.9999</math>
 | <math>9.9960</math>
@@ Line 1,153: / Line 1,229: @@
 | <math>0.0002</math>
 | <math>0.0002</math>
-|- style="border-bottom: 2px solid;"
+|- style="text-align:center;"
-| style="border-right: 2px solid;" | <math>200</math>
+|  <math>200</math>
-| style="border-left: 2px solid;" | <math>10.000</math>
+|  <math>10.000</math>
 | <math>9.9999</math>
 | <math>9.9999</math>
@@ Line 1,164: / Line 1,240: @@
 | <math>0.0002</math>
 | <math>0.0002</math>
 |}
 ===5.2 Considerações===
-Em  <span id='citeF-8'></span>[[#cite-8|[8]]], <span id='citeF-12'></span>[[#cite-12|[12]]] e <span id='citeF-11'></span>[[#cite-11|[11]]] são apontadas as dificuldades em lidar com  modelo que representam leitos quase secos, devido as possíveis singularidades.   O esquema numérico híbrido ''CDS/UDS'' aqui proposto mostrou-se eficaz e com baixo custo computacional para simular escoamentos transientes em  águas rasas para leitos quase secos, o que evidencia a acurácia e robustez do esquema.
+Em  Fennema and Chaudhry <span id='citeF-10'></span>[[#cite-10|[10]]], Fernández-Pato  et al. <span id='citeF-13'></span>[[#cite-13|[13]]] e Sheu and Fang <span id='citeF-15'></span>[[#cite-15|[15]]] e  são apontadas as dificuldades em lidar com  modelo que representam leitos quase secos, devido as possíveis singularidades.   O esquema numérico híbrido ''CDS/UDS'' aqui proposto mostrou-se eficaz e com baixo custo computacional para simular escoamentos transientes em  águas rasas para leitos quase secos, o que evidencia a acurácia e robustez do esquema.
-==6 Escoamento em um reservatório de contenção subterrâneo: uma aplicação do modelo completo==
+==6. Escoamento em um reservatório de contenção subterrâneo: uma aplicação do modelo completo==
-Os reservatórios subterrâneos são projetados para zonas urbanas populosas, onde não há lugar para armazenar água da chuva. Eles são construídos com uma rede de pilares de sustentação e canais de drenagem na parte inferior, como ilustrado na Fig. ([[#img-11|11]])   <div id='img-11'></div>
+Os reservatórios subterrâneos são projetados para zonas urbanas populosas, onde não há lugar para armazenar água da chuva. Eles são construídos com uma rede de pilares de sustentação e canais de drenagem na parte inferior, como ilustrado na  [[#img-11|Figura 11]].
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+<div id='img-11'></div>
+{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 60%;"
 |-
-|[[Image:Draft_Coelho_632154385-fig11.png|360px|Reservatório de águas pluviais subterrâneas]]
+|style="padding:10px;"|[[Image:Draft_Coelho_632154385-fig11.png|460px|Reservatório de águas pluviais subterrâneas]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figure 11:''' Reservatório de águas pluviais subterrâneas
+| colspan="1" style="padding-bottom:10px;"| '''Figura 11'''. Reservatório de águas pluviais subterrâneas
 |}
-Para testar a solução numérica obtida para as equações completas das águas rasas, analisou-se o fluxo em um reservatório com obstáculos como colunas estruturais, e drenos controlados por saídas estreitas como fendas na barragem.  A topografia de fundo foi considerada com inclinações não nulas e com resistência ao escoamento no fundo do canal, com coeficiente de atrito de Manning <math>N = 0.013 m/s^{1/3}</math>, que quantifica a resistência ao escoamento no fundo do canal para superfícies de cimento, apresentado em <span id='citeF-24'></span>[[#cite-24|[24]]].   O fundo tem a presença de canais que conduzem a água pelas aberturas de saída até uma área aberta sem paredes. Pode-se dizer que o problema considerado é uma composição do problema de abertura instantânea de comportas, para leitos molhado/seco, com topografia mais complexa e configuração de domínio.  A Fig. ([[#img-12|12]]) representa a configuração do domínio físico. O domínio consiste em um reservatório retangular de <math> 12 m \times 24  m </math> de dimensão com nove pilares de sustentação. <div id='img-12'></div>
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+Para testar a solução numérica obtida para as equações completas das águas rasas, analisou-se o fluxo em um reservatório com obstáculos como colunas estruturais, e drenos controlados por saídas estreitas como fendas na barragem.  A topografia de fundo foi considerada com inclinações não nulas e com resistência ao escoamento no fundo do canal, com coeficiente de atrito de Manning <math>N = 0.013 m/s^{1/3}</math>, que quantifica a resistência ao escoamento no fundo do canal para superfícies de cimento, apresentado em <span id='citeF-28'></span>[[#cite-28|[28]]].   O fundo tem a presença de canais que conduzem a água pelas aberturas de saída até uma área aberta sem paredes. Pode-se dizer que o problema considerado é uma composição do problema de abertura instantânea de comportas, para leitos molhado/seco, com topografia mais complexa e configuração de domínio.  A  [[#img-12|Figura 12]] representa a configuração do domínio físico. O domínio consiste em um reservatório retangular de <math> 12 m \times 24  m </math> de dimensão com nove pilares de sustentação.
+<div id='img-12'></div>
+{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 40%;"
 |-
-|[[Image:Draft_Coelho_632154385-fig12.png|600px|Domínio físico do reservatório subterrâneo]]
+<!--|style="padding:10px;"|[[Image:Draft_Coelho_632154385-fig12.png|600px|Domínio físico do reservatório subterrâneo]]-->
+|style="padding:10px;"|[[File:Coelho_et_al_2019b_3585_Image 2.png|200px|Domínio físico do reservatório subterrâneo]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figure 12:''' Domínio físico do reservatório subterrâneo
+| colspan="1" style="padding-bottom:10px;"| '''Figura 12'''. Domínio físico do reservatório subterrâneo
 |}
-Como ilustrado na Fig. ([[#img-13|13]]), a parte inferior do reservatório possui dois canais de drenagem com seção trapezoidal e um declive sutil em direção às fendas da barragem.  <div id='img-13a'></div>
+Como ilustrado na  [[#img-13|Figura 13]], a parte inferior do reservatório possui dois canais de drenagem com seção trapezoidal e um declive sutil em direção às fendas da barragem.
+<div id='img-13a'></div>
 <div id='img-13b'></div>
 <div id='img-13'></div>
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 85%;"
 |-
-|[[Image:Draft_Coelho_632154385-fig13a.png|471px|representação tridimensional do fundo]]
+|style="padding¡:10px;"|[[Image:Draft_Coelho_632154385-fig13a.png|571px|representação tridimensional do fundo]]
-|[[Image:Draft_Coelho_632154385-fig13b.png|471px|Perfil frontal do fundo do reservatório]]
+|style="padding:10px;"|[[Image:Draft_Coelho_632154385-fig13b.png|571px|Perfil frontal do fundo do reservatório]]
 |- style="text-align: center; font-size: 75%;"
 | (a) representação tridimensional do fundo
 | (b) Perfil frontal do fundo do reservatório
 |- style="text-align: center; font-size: 75%;"
-| colspan="2" | '''Figure 13:''' Topografia do fundo do reservatório
+| colspan="2" style="padding-bottom:10px;padding-top:10px;"| '''Figura 13'''. Topografia do fundo do reservatório
 |}
 As configurações da topografia do fundo são  descritas por <math>z_0</math> na Eq.([[#eq-62|62]]).
@@ Line 1,213: / Line 1,301: @@
 |}
-No reservatório, o nível da água está inicialmente em repouso e tem <math> 2.5\,m</math> de profundidade, enquanto na zona de inundação, há apenas uma camada fina de água, ou seja, “wet condition”, a altura da água é tão pequena quanto possível, como foi testado na Seção 5.1. Neste caso, a razão fixada é <math> h_r / h_l = 0.0002 </math>, como mostrado na Fig. ([[#img-14|14]]). Todas as condições de contorno para o nível da água são as mesmas do problema apresentado na Seção [[#4 |4]].  <div id='img-14'></div>
+No reservatório, o nível da água está inicialmente em repouso e tem <math> 2.5\,m</math> de profundidade, enquanto na zona de inundação, há apenas uma camada fina de água, ou seja, “wet condition”, a altura da água é tão pequena quanto possível, como foi testado na Seção 5.1. Neste caso, a razão fixada é <math> h_r / h_l = 0.0002 </math>, como mostrado na  [[#img-14|Figura 14]]. Todas as condições de contorno para o nível da água são as mesmas do problema apresentado na Seção [[#4 |4]].
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+<div id='img-14'></div>
+{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 60%;"
 |-
-|[[Image:Draft_Coelho_632154385-fig14.png|300px|Condições iniciais para o nível da água]]
+|style="padding:10px;"|[[Image:Draft_Coelho_632154385-fig14.png|450px|Condições iniciais para o nível da água]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figure 14:''' Condições iniciais para o nível da água
+| colspan="1" style="padding-bottom:10px;" | '''Figura 14'''. Condições iniciais para o nível da água
 |}
-Com relação aos componentes de velocidade, é necessário dizer que: como as paredes do reservatório, os pilares e paredes laterais do domínio externo são considerados impermeáveis; o fluxo da água que é armazenada no reservatório ocorre apenas nas saídas da parede norte do reservatório e na extremidade norte do domínio externo, como indicado na Fig. ([[#img-12|12]]). Mais precisamente, a componente de velocidade <math>u</math> nas paredes internas leste-oeste e nos contornos internos é nula, enquanto a componente de velocidade <math>v</math> é considerada nula no limite sul e nos limites internos norte-sul. Por outro lado, nos contornos norte-sul para a componente <math>u</math> e nos contornos leste-oeste para o componente <math>v</math> respectivamente, é adotada a condição de Neumann nula, como na Seção [[#4 |4]].  Para o domínio externo, uma vez que existem apenas paredes laterais, o fluxo é livre para todos os pontos da extremidade norte. Essa condição de contorno é descrita pelas expressões:
+Com relação aos componentes de velocidade, é necessário dizer que: como as paredes do reservatório, os pilares e paredes laterais do domínio externo são considerados impermeáveis; o fluxo da água que é armazenada no reservatório ocorre apenas nas saídas da parede norte do reservatório e na extremidade norte do domínio externo, como indicado na  [[#img-12|Figura 12]]. Mais precisamente, a componente de velocidade <math>u</math> nas paredes internas leste-oeste e nos contornos internos é nula, enquanto a componente de velocidade <math>v</math> é considerada nula no limite sul e nos limites internos norte-sul. Por outro lado, nos contornos norte-sul para a componente <math>u</math> e nos contornos leste-oeste para o componente <math>v</math> respectivamente, é adotada a condição de Neumann nula, como na Seção [[#4 |4]].  Para o domínio externo, uma vez que existem apenas paredes laterais, o fluxo é livre para todos os pontos da extremidade norte. Essa condição de contorno é descrita pelas expressões:
 {| class="formulaSCP" style="width: 100%; text-align: left;"
@@ Line 1,246: / Line 1,338: @@
 ===6.1 Resultados===
-Devido às aberturas na parede norte do reservatório, há um declínio no nível da água dentro do reservatório e ondas de elevação no domínio externo, especialmente nas proximidades das aberturas.  A Fig. ([[#img-15|15]])  é uma representação em perspectiva da superfície da água em todo o domínio para <math> t = 21 \, s </math> e <math> t = 70 \, s </math> após a abertura das saídas. <div id='img-15a'></div>
+Devido às aberturas na parede norte do reservatório, há um declínio no nível da água dentro do reservatório e ondas de elevação no domínio externo, especialmente nas proximidades das aberturas.  A  [[#img-15|Figura 15]]  é uma representação em perspectiva da superfície da água em todo o domínio para <math> t = 21 \, s </math> e <math> t = 70 \, s </math> após a abertura das saídas.
+<div id='img-15a'></div>
 <div id='img-15b'></div>
 <div id='img-15'></div>
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 80%;"
 |-
-|[[Image:Draft_Coelho_632154385-fig15a.png|582px|t=21\,s]]
+|style="padding:10px;"|[[Image:Draft_Coelho_632154385-fig15a.png|620px|t=21\,s]]
-|[[Image:Draft_Coelho_632154385-fig15b.png|582px|t=70\,s]]
+|style="padding:10px;"|[[Image:Draft_Coelho_632154385-fig15b.png|620px|t=70\,s]]
 |- style="text-align: center; font-size: 75%;"
 | (a) <math>t=21\,s</math>
 | (b) <math>t=70\,s</math>
 |- style="text-align: center; font-size: 75%;"
-| colspan="2" | '''Figure 15:''' Visão tridimensional do nível da água após o início do despejo
+| colspan="2" style="padding:10px;"| '''Figura 15'''. Visão tridimensional do nível da água após o início do despejo
 |}
-A representação gráfica do campo vetorial é usada para a análise da intensidade e direção do fluxo, em cada ponto da malha. A Fig. ([[#img-16|16]]) mostra as direções e a intensidade do fluxo nos primeiros <math> 4 \, s </math> do fluxo. Observa-se que a direção do fluxo ocorre de dentro para fora do reservatório, a partir das saídas,  com gradientes mais acentuados próximos às aberturas.  <div id='img-16a'></div>
+A representação gráfica do campo vetorial é usada para a análise da intensidade e direção do fluxo, em cada ponto da malha. A  [[#img-16|Figura 16]] mostra as direções e a intensidade do fluxo nos primeiros <math> 4 \, s </math> do fluxo. Observa-se que a direção do fluxo ocorre de dentro para fora do reservatório, a partir das saídas,  com gradientes mais acentuados próximos às aberturas.
+<div id='img-16a'></div>
 <div id='img-16b'></div>
 <div id='img-16c'></div>
@@ Line 1,266: / Line 1,364: @@
 {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
 |-
-|[[Image:Draft_Coelho_632154385-fig16a.png|474px|t=1\,s]]
+|style="padding-top:10px;"|[[Image:Draft_Coelho_632154385-fig16a.png|474px|t=1\,s]]
-|[[Image:Draft_Coelho_632154385-fig16b.png|474px|t=2\,s]]
+|style="padding-top:10px;"|[[Image:Draft_Coelho_632154385-fig16b.png|474px|t=2\,s]]
 |- style="text-align: center; font-size: 75%;"
 | (a) <math>t=1\,s</math>
@@ Line 1,278: / Line 1,376: @@
 | (d) <math>t=4\,s</math>
 |- style="text-align: center; font-size: 75%;"
-| colspan="2" | '''Figure 16:''' Campo de velocidade nos primeiros <math>4</math> segundos
+| colspan="2" style="padding-bottom:10px;"| '''Figura 16'''. Campo de velocidade nos primeiros <math>4</math> segundos
 |}
-Para <math> t = 15 \, s </math> observa-se que ainda não há fluxo na região estreita à frente da parede norte do reservatório, nas seções onde não há aberturas. Por sua vez, nas regiões à frente das aberturas, o fluxo é mais acentuado e o campo de velocidade assume um perfil hiperbólico no domínio externo, como na Fig. ([[#img-17|17]])a. Por outro lado, como mostra a Fig. ([[#img-17|17]])b, a magnitude do vetor velocidade torna-se menos acentuada e mais uniforme, após <math> 70 \, s </math> de simulação do fluxo. <div id='img-17a'></div>
+Para <math> t = 15 \, s </math> observa-se que ainda não há fluxo na região estreita à frente da parede norte do reservatório, nas seções onde não há aberturas. Por sua vez, nas regiões à frente das aberturas, o fluxo é mais acentuado e o campo de velocidade assume um perfil hiperbólico no domínio externo, como na  [[#img-17|Figura 17]]a. Por outro lado, como mostra a  [[#img-17|Figura 17]]b, a magnitude do vetor velocidade torna-se menos acentuada e mais uniforme, após <math> 70 \, s </math> de simulação do fluxo.
+<div id='img-17a'></div>
 <div id='img-17b'></div>
 <div id='img-17'></div>
 {| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
 |-
-|[[Image:Draft_Coelho_632154385-fig17a.png|471px|t=15\,s]]
+|style="padding:10px;"|[[Image:Draft_Coelho_632154385-fig17a.png|471px|t=15\,s]]
-|[[Image:Draft_Coelho_632154385-fig17b.png|471px|t=70\,s]]
+|style="padding:10px;"|[[Image:Draft_Coelho_632154385-fig17b.png|471px|t=70\,s]]
 |- style="text-align: center; font-size: 75%;"
 | (a) <math>t=15\,s</math>
 | (b) <math>t=70\,s</math>
 |- style="text-align: center; font-size: 75%;"
-| colspan="2" | '''Figure 17:''' Campo de velocidade após o início do despejo
+| colspan="2" style="padding-bottom:10px;"| '''Figura 17'''. Campo de velocidade após o início do despejo
 |}
-O perfil transversal do nível da água pode ser analisado em um determinado instante de tempo. Na Fig. ([[#img-18|18]]), o perfil <math> h (0,75, y) </math> é representado para <math> t = 10 \, s </math>. Uma inclinação significativa no nível da água é observada na faixa de <math> 8 \, m <y <8.15 \, m </math>, o que era esperado devido à condição inicial descontínua imposta.  <div id='img-18'></div>
-{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 100%;max-width: 100%;"
+O perfil transversal do nível da água pode ser analisado em um determinado instante de tempo. Na  [[#img-18|Figura 18]], o perfil <math> h (0,75, y) </math> é representado para <math> t = 10 \, s </math>. Uma inclinação significativa no nível da água é observada na faixa de <math> 8 \, m <y <8.15 \, m </math>, o que era esperado devido à condição inicial descontínua imposta.
+<div id='img-18'></div>
+{| class="floating_imageSCP" style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: 40%;"
 |-
-|[[Image:Draft_Coelho_632154385-fig18.png|240px|Perfil do nível de água em uma linha de amostra em x = 0.75\,m para t = 10\,s ]]
+|style="padding:10px;"|[[Image:Draft_Coelho_632154385-fig18.png|440px|Perfil do nível de água em uma linha de amostra em x = 0.75\,m para t = 10\,s ]]
 |- style="text-align: center; font-size: 75%;"
-| colspan="1" | '''Figure 18:''' Perfil do nível de água em uma linha de amostra em <math>x = 0.75\,m</math> para <math>t = 10\,s</math>
+| colspan="1" style="padding-bottom:10px;"| '''Figura 18'''. Perfil do nível de água em uma linha de amostra em <math>x = 0.75\,m</math> para <math>t = 10\,s</math>
 |}
 O código implementado foi utilizado para análises como: variação de volume no reservatório em função do tempo, o comportamento do fluxo ao redor dos pilares e na vizinhança das paredes. Em todas essas análises, os resultados obtidos são satisfatórios e consistentes com a física do problema.
-==7 Conclusões==
+==7. Conclusões==
 Foi proposto um método de volume finito implícito de <math>1^a</math> ordem com uma discretização híbrida CDS/UDS para resolver as equações completas de águas rasas. O algoritmo mostrou-se capaz de lidar com fortes problemas instáveis como a abertura instantânea de comportas, vazão em topografias não uniformes com condição ''dry/wet'', com e sem condições de saídas.  Os resultados mostram boa concordância com a literatura com dinâmica de frente de onda mais complexa. Essas diferenças locais são explicadas devido à ausência de artifício numérico, como a viscosidade artificial, que amortece o comportamento da superfície livre e suaviza as ondas.  A formulação numérica proposta tem a vantagem de ser mais fácil de codificar em comparação com aquelas disponíveis na literatura e, além disso, o apelo físico para a conservação de fluxo que é característica do método de volumes finitos, torna-se uma opção com a mesma precisão física, porém  simples e rápida para ser implementada. A convergência do método iterativo de Gauss-Seidel, mesmo impondo máximo residual muito pequeno, é obtida com poucas iterações, o que garante baixo custo computacional.  Como ferramenta preditiva, o modelo computacional tem mostrado bons resultados desde a configuração simples até a mais complexa, com precisão e robustez.  A solução numérica obtida pode ser usada para simular fluxos sobre topografias reais encontradas na natureza.
-==8 Declaração de disponibilidade de dados==
+==Declaração de disponibilidade de dados==
 Todos os dados, modelos ou códigos gerados ou utilizados durante o estudo estão disponíveis pelo autor correspondente por solicitação.  Lista de itens:
@@ Line 1,319: / Line 1,427: @@
 Este estudo foi financiado em parte pela Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - Brasil (CAPES) - Finance Code 001.   Os autores agradecem também ao Departamento de Matemática da Universidade Tecnológica Federal do Paraná (UTFPR), Campus Campo Mourão e ao Programa de Pós-Graduação em Métodos Numéricos (PPGMNE) da Universidade Federal do Paraná (UFPR).
-==9 Referências==
+==Referências==
-===BIBLIOGRAPHY===
+<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]]]''' Imbonden, D. M. (2004) "The Lakes Handbook, Limmnology and Limnectic Ecology". Blackwell Scienc Ltd
+[1] Imbonden D M.  The lakes handbook, immnology and limnectic ecology. Blackwell Scienc Ltd., 2004.
 <div id="cite-2"></div>
-'''[[#citeF-2|[2]]]''' Reynolds, C. (1984) "The ecology of freshwater phytoplankton". Cambridge University Press
+[2] Reynolds C.  The ecology of freshwater phytoplankton. Cambridge University Press, 1984.
 <div id="cite-3"></div>
-'''[[#citeF-3|[3]]]''' Mahmood, K. and Yevjevich, V. (1975) "Unsteady flow in open channels". Water Resources Publications, 1 th Edition
+[3] Mahmood K., Yevjevich V. Unsteady flow in open channels. Water Resources Publications, 1st Edition, 1975.
 <div id="cite-4"></div>
-'''[[#citeF-4|[4]]]''' Vreugdenhil, C. B. (1994) "Numerical Methods for Shallow-Water flow". Springer Science+Business Media Dordrecht
+[4] Vreugdenhil C.B. Numerical Methods for Shallow-Water flow. Springer Science+Business Media Dordrecht, 1994.
 <div id="cite-5"></div>
-'''[[#citeF-5|[5]]]''' García-Navarro, P. and Alcrudo, F. and Priestley, A. (1994) "An implicit method for water flow modelling in channels and pipes", Volume 32. J. Hydraul. Eng. 721&#8211;742
+[5] García-Navarro P., Alcrudo F., Priestley A.  An implicit method for water flow modelling in channels and pipes. J. Hydraul. Eng., 32:721&#8211;742, 1994.
 <div id="cite-6"></div>
-'''[[#citeF-6|[6]]]''' Burguete, J. and García-Navarro, P. (2004) "Implicit schemes with large time step for non-linear equations: Application to river flow hydraulic", Volume 46. Internat. J. Numer. Methods Fluids 607&#8211;636
+[6] Burguete J., García-Navarro P.  Implicit schemes with large time step for non-linear equations: Application to river flow hydraulic. Int. J. Numer. Methods , 46:607&#8211;636, 2004.
 <div id="cite-7"></div>
-'''[[#citeF-7|[7]]]''' Cantero-Chinchilla, F. N. and Castro-Orgaz, O. and Dey, S. and Ayuso, J. L. (2016) "Nonhydrostatic Dam Break Flows. I: Physical Equations and Numerical Schemes", Volume 142. J. Hydraul. Eng. 12 1&#8211;19
+[7] Cantero-Chinchilla F.N., Castro-Orgaz O., Dey S., Ayuso J.L.  Nonhydrostatic dam break flows. I: Physical equations and numerical schemes. J. Hydraul. Eng., 142(12):1&#8211;19, 2016.
 <div id="cite-8"></div>
-'''[[#citeF-8|[8]]]''' Fennema, R. J. and Chaudhry, M.H. (1990) "Explicit Methods for 2-D transient free surface flows", Volume 116. J. Hydraul. Eng. 8 1013&#8211;1034
+[8] Soler J., Bladé E., Sánchez-Juny M. Ensayo comparativo entre modelos unidimensionales y bidimensionales en la modelización de la rotura de una balsa de materiales sueltos erosionables. Revista International  de Métodos Numéricos para Cálculo y Diseño en Ingeniería, 28(2):103&#8211;105, 2012.
 <div id="cite-9"></div>
-'''[[#citeF-9|[9]]]''' Nikolos, I.K. and Delis, A.I. (2009) "An unstructured node-centered finite volume scheme for shallow water flows with wet/dry fronts over complex topography", Volume 198. Comput. Methods. Appl. Mech. Engrg 3723&#8211;3750
+[9] Brufau P., García-Navarro P.  Esquemas de alta resolución para resolver las ecuaciones de shallow water. Revista International  de Métodos Numéricos para Cálculo y Diseño en Ingeniería, 16(4):493&#8211;512, 2000.
 <div id="cite-10"></div>
-'''[[#citeF-10|[10]]]''' Valiani, A. and Caleffi, V. and Zanni, A. (2002) "Case Study: Malpasset Dam-Break Simulation using a Two-Dimensional Finite Volume Method", Volume 128. J. Hydraul. Eng. 5 460&#8211;472
+[10] Fennema R.J.,  Chaudhry M.H.  Explicit Methods for 2-D transient free surface flows. J. Hydraul. Eng., 116(8):1013&#8211;1034, 1990.
 <div id="cite-11"></div>
-'''[[#citeF-11|[11]]]''' Fernández-Pato, J. and Morales-Hernández, M. and García-Navarro, P. (2018) "Implicit finite volume simulation of 2D shallow water flows in flexible meshes", Volume 328. Comput. Methods. Appl. Mech. Engrg 1&#8211;25
+[11] Nikolos I.K.,  Delis A.I. An unstructured node-centered finite volume scheme for shallow water flows with wet/dry fronts over complex topography. Comput. Methods. Appl. Mech. Engrg., 198:3723&#8211;3750, 2009.
 <div id="cite-12"></div>
-'''[[#citeF-12|[12]]]''' Sheu, T. W. H. and Fang, C. C. (2001) "High resolution finite-element analysis of shallow water equations in two dimensions", Volume 190. Comput. Methods. Appl. Mech. Engrg 2581&#8211;2601
+[12] Valiani A., Caleffi V., Zanni A. Case study: malpasset dam-Break simulation using a two-dimensional finite volume method. J. Hydraul. Eng., 128(5):460&#8211;472, 2002.
 <div id="cite-13"></div>
-'''[[#citeF-13|[13]]]''' Hirsch, J. J. (2007) "Water Waves: The Mathematical Theory with Applications", Volume. Interscience Publishers, Edition
+[13] Fernández-Pato J.,  Morales-Hernández M., García-Navarro P.  Implicit finite volume simulation of 2D shallow water flows in flexible meshes. Comput. Methods. Appl. Mech. Engrg., 328:1&#8211;25, 2018.
 <div id="cite-14"></div>
-'''[[#citeF-14|[14]]]''' Blazek, J. (2001) "Computational fluid dynamics: principles and applications.". Elsevier Science Ltd, 1 th Edition
+[14] Melo A.R.,  Gramani L.M.,  Kaviski E.  High-order CE/SE scheme for dam-break flow simulation. Revista International  de Métodos Numéricos para Cálculo y Diseño en Ingeniería, 34(1):1&#8211;12, 2018.
 <div id="cite-15"></div>
-'''[[#citeF-15|[15]]]''' Versteeg, H.K. and Malalasekera, W. (2007) "An Introduction to Computational Fluid Dynamics. The finite volume method". Pearson Education Limited, 2 th Edition
+[15] Sheu T.W.H.,  Fang C.C. High resolution finite-element analysis of shallow water equations in two dimensions. Comput. Methods. Appl. Mech. Engrg., 190:2581&#8211;2601, 2001.
 <div id="cite-16"></div>
-'''[[#citeF-16|[16]]]''' Maliska, C. R. (2013) "Transferência de Calor e Mecânica dos Fluidos Computacional". LTC, 2 th Edition
+[16] Hirsch J.J. Water waves: The mathematical theory with applications. Interscience Publishers, Edition, 2007.
 <div id="cite-17"></div>
-'''[[#citeF-17|[17]]]''' Patankar, S. (1980) "Numerical heat transfer and fluid flow". Hemisphere Publishing Corporation, Edition
+[17] Blazek J. Computational fluid dynamics: principles and applications. Elsevier Science Ltd, 1st Edition, 2001.
 <div id="cite-18"></div>
-'''[[#citeF-18|[18]]]''' Melo, A. R. (2001) "Titulo kkkkkkkk", Volume 190. Comput. Methods. Appl. Mech. Engrg 2581&#8211;2601
+[18] Versteeg H.K.,  Malalasekera W. An Introduction to computational fluid dynamics. The finite volume method. Pearson Education Limited, 2nd Edition,  2007.
 <div id="cite-19"></div>
-'''[[#citeF-19|[19]]]''' Rezende, R. V. P. and Almeida, R. A. and Souza, A. A. U. and Souza, S. M. A. G. U. (2015) "A two-fluid model with a tensor closure model approach for free surface flow simulations.", Volume 122. Elsevier. Chemical Engineering Science 596&#8211;613
+[19] Maliska C.R.  Transferência de calor e mecânica dos fluidos computacional. LTC, 2nd Edition, 2013.
 <div id="cite-20"></div>
-'''[[#citeF-20|[20]]]''' Hirt, C. W. and Nichols, B. D. (1981) "Volume of fluid (VOF) method for the dynamics of free boundaries", Volume 39. ScienceDirect. Journal of Computational Physics 201&#8211;225
+[20] Patankar, S. (1980) Numerical heat transfer and fluid flow. Hemisphere Publishing Corporation, Edition
 <div id="cite-21"></div>
-'''[[#citeF-21|[21]]]''' Jahanbakhsh, E. and Panahi, R. and Seif, M. S. (2007) "Numerical simulation of three-dimensional interfacial flows.", Volume 17. Int. J. Numer. Methods Heat Fluid Flow 384&#8211;404
+[21] Lemos, E. M. D. (2011) Implementação de um método de volumes finitos de ordem superior com tratamento multibloco aplicado à simulação de escoamentos em fluidos viscoelásticos. Phd. thesis, UFRJ, Rio de Janeiro, Brasil
+<div id="cite-22"></div>
+[22] Rezende, R. V. P. and Almeida, R. A. and Souza, A. A. U. and Souza, S. M. A. G. U. (2015) A two-fluid model with a tensor closure model approach for free surface flow simulations, Volume 122. Chemical Engineering Science 596&#8211;613
+<div id="cite-23"></div>
+[23] Hirt, C. W. and Nichols, B. D. (1981) Volume of fluid (vof) method for the dynamics of free boundaries, Volume 39. Journal of ComputationalPhysics 201&#8211;225
+<div id="cite-24"></div>
+[24] Jahanbakhsh, E. and Panahi, R., Seif, M. S.(2007) Numerical simulation fo three-dimensional interfacial flows, Volume 17. Int. J. Numer. Methods Heat Fluid Flow 384&#8211;404
+<div id="cite-25"></div>
+[25] Villota, A., Codina, R. (2018) Approximation of the shallow water equations with higher order finite elements and variational multiscale methods, Volume 34 (1). Revista International  de Métodos Numéricos para Cálculo y Diseño en Ing. 1&#8211;25
-<div id="cite-22"></div>
+<div id="cite-26"></div>
-'''[[#citeF-22|[22]]]''' MacCormack, R. W. (1969) "The effect of viscosity in hypervelocity impact cratering", Volume. American Institutue Aeronautics Astronautics 69&#8211;354
+[26] MacCormack, R. W. (1969) The effect of viscosity in hypervelocity impact cratering, Volume. American Institutue Aeronautics Astronautics 69&#8211;354
 <div id="cite-23"></div>
-'''[[#citeF-23|[23]]]''' Gabutti, B. (1983) "On two upwind finite-difference schemes for hyperbolic equations in non-conservative form", Volume 11. Comput. Fluids 3 207&#8211;230
+[27] Gabutti, B. (1983) On two upwind finite-difference schemes for hyperbolic equations in non-conservative form, Volume 11. Comput. Fluids 3 207&#8211;230
-<div id="cite-24"></div>
+<div id="cite-28"></div>
-'''[[#citeF-24|[24]]]''' Porto, R. M. (2006) "Hidráulica Básica". EESC-USP, 4 th Edition
+[28] Porto, R. M. (2006) Hidráulica Básica. EESC-USP, 4 th Edition
+</div>