verificação de funções de interpolação em advecção-difusão 1d ...

72
Para facilitar o entendimento, a partir deste ponto as ações foram divididas em três
etapas, que foram as seguintes.
3.5.1 Obtenção das Soluções Numéricas das Variáveis de Interesse
Conforme já dito em seção anterior, para cada conjunto de FI utilizada, e aí entendase
conjunto a FI para o termo advectivo e outra para o termo difusivo, foi gerado um sistema
de equações algébricas, que foram resolvidas através do método TDMA.
O programa computacional utilizado para obtenção destes valores foi
Peclet_1Dp_2p2, tendo sido o programa gerado em 09-set-2008. Cujo algoritmo de solução
tem, de uma forma geral, a seguinte rotina :
1) leitura dos dados de entrada (variáveis, condições de contorno, número de volumes, valor
de Pe, função de interpolação utilizada, nome do arquivo de saída);
2) obtenção do tamanho dos volumes da malha;
3) cálculo dos coeficientes e termos fonte de cada volume da malha;
4) resolução do sistema de equações algébricas, obtendo a solução numérica, com o uso do
método TDMA; e
5) visualização dos resultados das variáveis e erros de discretização, obtidos através da
comparação da solução numérica com a solução analítica.
E finalmente, para obtenção da solução numérica final, foram assumidos ainda:
Precisão quádrupla (REAL*16);
O número de Peclet considerado para todos os processos de simulação foi 5;
Para a simulação com aproximação ALFA, considerou-se o coeficiente alfa
como 0,05;
O refinamento da malha foi de ordem 3, iniciando a primeira simulação com 5
nós e a última simulação com 23.914.845 nós. Ou seja, para cada FI obteve-se a
solução numérica em 15 malhas distintas;
O processo acima foi feito para cada FI (10 no total) considerada, totalizando
150 simulações. Nesta fase foram obtidos os valores numéricos das variáveis de
interesse para cada simulação feita. Isto significa que cada variável de interesse
teve 15 resultados (das 15 malhas) para cada função de interpolação;