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.
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