UTILIZAÇÃO DA MODELAGEM COMPUTACIONAL PARA AVALIAR ...

esultar em um contorno irregular, que pode causar instabilidade dentro de um
modelo de elementos finitos.
Já nos algoritmos de “fissura”, os elementos de secagem e os nós não são
removidos do domínio. Em vez disso, os elementos e os nós entram num domínio de
reduzida capacidade de escoamento, geralmente representado conceitualmente por
uma fissura. Ou seja, ocorrendo o secamento dos nós de um elemento da malha o
escoamento continua a um nível subterrâneo e a profundidade da água neste
domínio de escoamento é parametrizada e transformada numa profundidade d‟água
positiva.
Desta forma, o nó ou elemento permanece no domínio da solução, removendo o
problema de fronteiras dos modelos de elementos finitos. Isso elimina os problemas
com a continuidade de massa e melhora a estabilidade global da solução. O método
geralmente exige uma região em que há uma transição gradual entre um elemento
com escoamento totalmente molhado e um elemento totalmente seco, com
capacidade de vazão reduzida. Outros parâmetros podem variar dentro deste
domínio de alagamento e secamento, incluindo a especificação da velocidade zero
dentro do elemento, o aumento da rugosidade e da escala de redução da
capacidade de escoamento.
Nielsen e Alpet (2003) investigaram o efeito da variação dos parâmetros da
porosidade de pântano no desempenho do modelo hidrodinâmico bidimensional
RMA2, baseado em elementos finitos, em situações em que ocorre a inundação de
áreas inicialmente secas, tais como planícies de maré. Estes autores utilizaram o
método da porosidade de pântano, que é um tipo de algoritmo de alagamento e
secamento de “fissura”. Em tal estudo foram realizadas comparações entre
exemplos simples e cálculos teóricos.
A esquematização do método de fissura padrão é mostrada na Figura 5 - a. Neste
esquema os elementos e nós secos da malha nunca são removidos do domínio
computacional, ao invés disso, estes entram num domínio de capacidade de
escoamento reduzido (zona subterrânea de baixa porosidade).
Neste caso a redução na capacidade de escoamento ocorre sobre todo o elemento
e em ambas as direções de escoamento horizontal, e a transição na porosidade
37