17.04.2013 Views

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

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!