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

122
Considerando então a forma genérica da Eq.(A.4.1), temos para o nó P=N+1, temos:
a
E
0
a
W
w
1
2
aP
1
b 1
(A.4.27)
2
w
P
A.4.5 PLDS
Consideraremos neste caso a FI PLDS para os termo advectivo e difusivo.
Inicialmente definiremos numericamente o termo que tem a função de balancear a
intensidade de cada termo no fenômeno. Desta forma:
( 1
0,1. Pe)
5
(A.4.28)
Pe
A.4.5.1
Nós P=1 a P=N
Aplicando as aproximações do PLDS para os volumes reais (P=1 a P=N), com a
equação (A.3.3), temos:
Pe. P E P W P W
x
x
E P P W
0
(A.4.29)
Considerando então a forma genérica da Eq.(A.4.1), temos para os nós P=1 a P=N:
a E
Pe. x.
1
Pe. x.
1
1
1
2. 2
aP aE
aW
Pe. x.
b 0
(A.4.30)
a W
P
A.4.5.2
Nós P=0 e P=N+1
E finalmente para os nós P=0 e P=N+1, temos os coeficientes dados pelas relações
das Eqs.(A.4.4) e (A.4.5) respectivamente.