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 ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
41<br />
A condição principal para um esqu<strong>em</strong>a TVD é que, a variação total da solução, dada<br />
pela Eq.2.26, <strong>de</strong>cresça com as iterações (BLAZEK, 2001).<br />
TVD U I <br />
U<br />
)<br />
(2.26)<br />
I<br />
(<br />
1 I<br />
Os esqu<strong>em</strong>as TVD são impl<strong>em</strong>entados como uma média dos fluxos advectivos,<br />
combinados com um termo <strong>de</strong> dissipação adicional. Se o termo <strong>de</strong> dissipação <strong>de</strong>pen<strong>de</strong>r do<br />
sinal da velocida<strong>de</strong>, então é chamado TVD simétrico, caso contrário é um TVD upwind, que<br />
t<strong>em</strong> melhor or<strong>de</strong>m <strong>de</strong> acurácia que o anterior (BLAZEK, 2001). É um esqu<strong>em</strong>a <strong>de</strong> 2ª or<strong>de</strong>m<br />
<strong>em</strong> regiões <strong>de</strong> fluxos suaves, e <strong>de</strong> 1ª or<strong>de</strong>m <strong>em</strong> regiões <strong>de</strong> <strong>de</strong>scontinuida<strong>de</strong> (BLAZEK, 2001).<br />
O esqu<strong>em</strong>a TVD utilizado neste trabalho será do tipo Superbee, proposto por Roe <strong>em</strong><br />
1983. Sua representação, como segue, foi tirada <strong>de</strong> Marchi (1993),<br />
1<br />
<br />
e<br />
(1 e<br />
)<br />
(2.27)<br />
2<br />
r e<br />
no cálculo <strong>de</strong><br />
* *<br />
( <br />
P<br />
W<br />
)<br />
(2.28)<br />
* *<br />
( )<br />
E<br />
<br />
P<br />
<br />
MAX ; MIN(2r<br />
;1); MIN(<br />
r ;2)<br />
(2.29)<br />
e<br />
0<br />
e<br />
e<br />
As Eqs. (2.27) e (2.28) são válidas para u>0. O fator r e introduz uma não-linearida<strong>de</strong><br />
<br />
e<br />
e <br />
e<br />
(MARCHI, 1993). Após calculados os valores <strong>de</strong> “α” para a face leste,<br />
e face oeste, estes <strong>de</strong>v<strong>em</strong> ser substituídos nas Eqs.(2.22) e (2.23) para obtenção das<br />
aproximações com o esqu<strong>em</strong>a TVD.<br />
2.8.9 ADS<br />
O esqu<strong>em</strong>a ADS, ou Adaptable Difference Sch<strong>em</strong>e, é apresentado a seguir pelas Eqs.<br />
*<br />
(2.30), (2.31), (2.32) e (2.33) (MARCHI, 1993). Neste esqu<strong>em</strong>a o valor <strong>de</strong> e<br />
não <strong>de</strong>pen<strong>de</strong> do<br />
sentido do escoamento.<br />
<br />
<br />
e<br />
1 *<br />
*<br />
*<br />
P<br />
.( P<br />
W<br />
)<br />
2<br />
(2.30)<br />
<br />
m<br />
e<br />
1 *<br />
*<br />
.( <br />
P<br />
E<br />
)<br />
2<br />
(2.31)