18.07.2014 Views

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 ...

SHOW MORE
SHOW LESS

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)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!