verificação de funções de interpolação em advecção-difusão 1d ...
verificação de funções de interpolação em advecção-difusão 1d ...
verificação de funções de interpolação em advecção-difusão 1d ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
39<br />
2.8.5 PLDS (Power Law Difference Sch<strong>em</strong>e)<br />
Esqu<strong>em</strong>a <strong>de</strong>senvolvido por Patankar <strong>em</strong> 1979. Sua base conceitual é a busca por<br />
corrigir algumas inconsistências do esqu<strong>em</strong>a híbrido. A principal <strong>de</strong>las é a consi<strong>de</strong>ração <strong>de</strong><br />
que o efeito da difusão cessa tão logo tenha-se<br />
malha, dado por,<br />
2 Peh<br />
e Pe<br />
h<br />
2, on<strong>de</strong> Pe<br />
h<br />
é o Peclet da<br />
. u.<br />
x<br />
Pe h<br />
,<br />
<br />
(2.16)<br />
Neste esqu<strong>em</strong>a, o efeito da difusão cessa assim que Pe exce<strong>de</strong> o valor 10. Se Pe está<br />
entre zero e <strong>de</strong>z, o fluxo é analisado usando uma expressão polinomial.<br />
Consi<strong>de</strong>rando a possibilida<strong>de</strong> <strong>de</strong> escrever a equação <strong>de</strong> advecção-difusão discretizada<br />
<strong>de</strong> uma forma genérica por (PATANKAR, 1980),<br />
a . a . a . b , (2.17)<br />
P<br />
P<br />
E<br />
E<br />
W<br />
W<br />
po<strong>de</strong>-se <strong>de</strong>finir os valores <strong>de</strong> a E e a W para o esqu<strong>em</strong>a PLDS como,<br />
a<br />
a<br />
E<br />
W<br />
p<br />
D .max 0,<br />
(2.18)<br />
e<br />
5<br />
1<br />
0,1. Pee<br />
<br />
max<br />
Fe<br />
,0<br />
5<br />
1<br />
0,1. Pew<br />
<br />
maxFw<br />
,0<br />
D .max 0,<br />
(2.19)<br />
e<br />
sendo que, as relações D, F e Pe são dadas por,<br />
Pe <br />
F<br />
D<br />
. u<br />
<br />
<br />
x<br />
(2.20)<br />
A proximida<strong>de</strong> <strong>de</strong>ste esqu<strong>em</strong>a ao esqu<strong>em</strong>a exponencial é muito gran<strong>de</strong>, <strong>em</strong> termos <strong>de</strong><br />
resultado, o que leva Patankar (1980) a consi<strong>de</strong>rá-lo como o mais propício <strong>em</strong> probl<strong>em</strong>as <strong>de</strong><br />
advecção-difusão, sendo o esqu<strong>em</strong>a híbrido relegado somente a alguns casos mais específicos.<br />
2.8.6 WUDS<br />
Abreviação <strong>de</strong> Weighted Upstream Differencing Sch<strong>em</strong>e, esta função <strong>de</strong> interpolação<br />
é associada a dois coeficientes, e , que <strong>de</strong>pen<strong>de</strong>m do número <strong>de</strong> Peclet (Pe) e traduz<strong>em</strong> o<br />
peso da advecção e difusão no fenômeno. Como esta seção trata da aproximação do termo