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

124
Considerando então a forma genérica da Eq.(A.4.1), temos para os nós P=1 a P=N:
a
E
1
1
1e
e .
Pe.
x
aW
1
w .
Pe.
x
2
2
a
P
a a 2 ( ). Pe.
x
b 0 (A.4.39)
E
W
e
w
P
A.4.6.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.
A.4.7 TVD / CDS-2
Consideraremos neste caso a FI TVD para o termo advectivo e CDS-2 para o termo
difusivo. Para definição do coeficiente
representadas pelas Eqs. (A.4.40) a (A.4.42)
r e
e
precisamos inicialmente calcular as relações
* *
(
P
W
)
(A.4.40)
* *
( )
E
P
MAX ; MIN(2r
;1); MIN(
r ;2)
(A.4.41)
e
0
e
e
1
e
(1 e
)
(A.4.42)
2
Sendo que:
Se P 0 ou P N 0
(A.4.43)
e
A.4.7.1
Nós P=1 a P=N
Aplicando as aproximações do ADS para os volumes reais (P=1 a P=N), com a
equação (A.3.3), temos:
1 1 1
1
E
P
P
W
Pe. (
e).
P
(
e).
E
(
w).
W
(
w).
P
0
2 2
2
2
x
x
(A.4.44)