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.

A maioria das rotinas <strong>de</strong> elementos finitos foram expressamente escritas para<br />

tirar proveito da natureza <strong>de</strong> banda da matriz <strong>de</strong> coeficientes. A largura <strong>de</strong> banda <strong>de</strong><br />

<strong>um</strong>a linha da matriz <strong>de</strong> coeficientes é <strong>de</strong>finida como o número <strong>de</strong> colunas entre o primeiro<br />

coeficiente não-nulo nessa linha e a diagonal da matriz. A máxima largura <strong>de</strong> banda é a<br />

maior largura entre todas as linhas da matriz. Esquemas <strong>de</strong> solução baseados na máxima<br />

largura <strong>de</strong> banda são bastante simples <strong>de</strong> se programar.<br />

As posições dos coeficientes não-nulos na matriz K, e portanto as larguras <strong>de</strong><br />

banda, <strong>de</strong>pen<strong>de</strong>m somente da or<strong>de</strong>nação das incógnitas. Na maioria dos algoritmos <strong>de</strong><br />

solução <strong>de</strong> matrizes <strong>de</strong> banda, essa or<strong>de</strong>m é, <strong>de</strong> certa forma, baseada na n<strong>um</strong>eração dos<br />

nós da malha. Em aplicações <strong>de</strong> elementos finitos, a máxima largura <strong>de</strong> banda B é<br />

tipicamente calculada como:<br />

B = (R + 1) NDF , (6.13)<br />

on<strong>de</strong> R é a maior diferença entre o número dos nós <strong>de</strong> <strong>um</strong> mesmo elemento da malha e<br />

NDF é o número <strong>de</strong> incógnitas (graus <strong>de</strong> liberda<strong>de</strong>) em cada nó.<br />

Rotinas <strong>de</strong> solução baseadas na máxima largura <strong>de</strong> frente da matriz <strong>de</strong> coeficientes<br />

são também usadas em aplicações <strong>de</strong> elementos finitos. A largura <strong>de</strong> frente <strong>de</strong> <strong>um</strong>a linha<br />

da matriz <strong>de</strong> coeficientes é <strong>de</strong>finida como o número <strong>de</strong> colunas ativas naquela linha. Uma<br />

coluna j é dita ativa em <strong>um</strong>a certa linha i se j for maior do que i e houver ao menos <strong>um</strong><br />

coeficiente não-nulo <strong>de</strong>ssa coluna em <strong>um</strong>a linha com índice k maior que i.<br />

Esquemas <strong>de</strong> solução baseados na máxima largura <strong>de</strong> frente são diretamente rela-<br />

cionados ao processo <strong>de</strong> agrupamento dos elementos. Técnicas <strong>de</strong> solução frontal po<strong>de</strong>m<br />

ser baseadas tanto na <strong>de</strong>composição da matriz <strong>de</strong> coeficientes quanto em <strong>um</strong> procedimento<br />

<strong>de</strong> fatoração e agrupamento <strong>de</strong> elemento por elemento. Se a solução frontal for baseada<br />

nas equações completamente agrupadas, a máxima largura <strong>de</strong> frente é também ditada<br />

pelo esquema <strong>de</strong> n<strong>um</strong>eração dos nós. Todavia, se o processo <strong>de</strong> agrupamento <strong>de</strong> elemento<br />

por elemento for utilizado, a máxima largura <strong>de</strong> frente <strong>de</strong>pen<strong>de</strong> apenas da seqüência em<br />

que os elementos são processados e não tem relação com a n<strong>um</strong>eração dos nós.<br />

79

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

Saved successfully!

Ooh no, something went wrong!