Modelagem Física e Computacional de um Escoamento Fluvial
Modelagem Física e Computacional de um Escoamento Fluvial
Modelagem Física e Computacional de um Escoamento Fluvial
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
mantida (a vazão <strong>de</strong> entrada no domínio iguala-se à vazão <strong>de</strong> saída), as vazões calculadas<br />
em seções internas do domínio po<strong>de</strong>m <strong>de</strong>sviar-se dos valores <strong>de</strong> entrada e saída. GEE<br />
e MacARTHUR (1978) apresentam <strong>um</strong>a avaliação dos erros <strong>de</strong> continuida<strong>de</strong> em função<br />
do refinamento da malha <strong>de</strong> elementos finitos. Para <strong>um</strong>a malha <strong>de</strong> 131 elementos em<br />
<strong>um</strong>a região <strong>de</strong> cerca <strong>de</strong> 10x10 km 2 os autores verificaram erros <strong>de</strong> até 21% da vazão em<br />
seções internas do mo<strong>de</strong>lo, valor que cai para 4% com <strong>um</strong>a malha <strong>de</strong> 189 elementos. O<br />
trabalho salienta a relevância da questão dos erros <strong>de</strong> continuida<strong>de</strong> em <strong>um</strong> mo<strong>de</strong>lo <strong>de</strong><br />
qualida<strong>de</strong> da água, preocupação citada também por outros pesquisadores. Em geral, a<br />
questão da continuida<strong>de</strong> é <strong>de</strong> menor importância em <strong>um</strong> estudo hidrodinâmico, on<strong>de</strong> as<br />
pressões e velocida<strong>de</strong>s são as variáveis <strong>de</strong> maior interesse. GEE e MacARTHUR (1978)<br />
apontam, também, que as velocida<strong>de</strong>s são mais afetadas pelos erros <strong>de</strong> continuida<strong>de</strong> do<br />
que as pressões.<br />
KING (1978) arg<strong>um</strong>enta que esse problema aparente do método dos elementos<br />
finitos po<strong>de</strong> ser utilizado como <strong>um</strong>a vantagem, aproveitando-se a verificação da continui-<br />
da<strong>de</strong> como <strong>um</strong>a medida da precisão da solução em várias regiões do mo<strong>de</strong>lo e da possível<br />
necessida<strong>de</strong> <strong>de</strong> refinamento da malha.<br />
6.4 SOLUÇÃO DO SISTEMA DE EQUAÇÕES NÃO-LINEARES DE ELEMENTOS<br />
FINITOS<br />
6.4.1 Classificação dos Métodos para Solução <strong>de</strong> Sistemas Não-Lineares<br />
A solução <strong>de</strong> <strong>um</strong> sistema <strong>de</strong> equações lineares na forma genérica<br />
73<br />
K U = f , (6.9)<br />
on<strong>de</strong> K é a matriz <strong>de</strong> coeficientes globais, U é o vetor coluna das incógnitas nos nós e<br />
f é <strong>um</strong> vetor coluna com funções conhecidas, po<strong>de</strong> ser obtida sem dificulda<strong>de</strong>s através<br />
<strong>de</strong> diversos métodos. A solução torna-se mais difícil para <strong>um</strong> sistema não linear, em que