fabiana de fátima giacomini - PG-Mec Programa de Pós-Graduação ...

141
φ
e
= 1
( − φ )
dφ 1
P
= 2 (D.6)
dx ∆x
e
CONDIÇÕES DE CONTORNO COM VOLUME FICTÍCIO
• Face oeste ( w ):
para P = 0 , tem-se φ φ = φ = 0 , a aproximação para este volume é:
cc
=
w a
φ
cc
φW
+ φP
= ⇒ φW
= −φ
P
(D.7)
2
para P = volumes internos, tem-se:
φ
w
φW
=
+ φ
2
P
dφ
dx
w
φP
−φW
=
∆x
(D.8)
Substituindo φ
W
= −φ
, os coeficientes assumem os valores de:
P
φ
w
= 0
d φ P
= 2 φ
(D.9)
dx ∆x
w
• Face leste ( e ):
para
P = N , tem-se φ φ = φ = 1, a aproximação para este volume é:
cc
=
e b
φ
cc
φE
+ φP
= ⇒ φE
= 2 φcc
− φP
(D.10)
2
para P = volumes internos, tem-se:
φP
+ φE
φe
=
2
dφ
dx
e
φE
− φP
=
∆x
(D.11)
Substituindo φ
E
= 2 φ − φ , os coeficientes assumem os valores de:
cc
P
φ
e
= φ cc
dφ
dx
e
( φ −φ
)
cc P
= 2 (D.12)
∆x