fabiana de fátima giacomini - PG-Mec Programa de Pós-Graduação ...
fabiana de fátima giacomini - PG-Mec Programa de Pós-Graduação ...
fabiana de fátima giacomini - PG-Mec Programa de Pós-Graduação ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
74<br />
I<br />
principal ( )<br />
Para encontrar a or<strong>de</strong>m do erro <strong>de</strong> truncamento a priori para a <strong>de</strong>rivada da variável<br />
φ com o esquema CDS-2, subtrai-se as Eqs. (3.43) e (3.44). Assim, o primeiro<br />
e<br />
termo do segundo membro representa a aproximação do esquema CDS-2 para o termo<br />
difusivo, dado pela Eq. (2.10). Os outros termos da série representam o valor do erro <strong>de</strong><br />
truncamento, dado pela Eq. (3.48).<br />
( φ −φ<br />
)<br />
I E P 1 III 2 1 V 4 1 VII 6<br />
φ<br />
e<br />
= − φe<br />
∆x<br />
− φe<br />
∆x<br />
− φe<br />
∆x<br />
−K (3.47)<br />
∆x<br />
24 1920 322560<br />
I 1 III 2 1 V 4 1 VII 6<br />
[ φ ] = − φ ∆x<br />
− φ ∆x<br />
− φ ∆x<br />
−K<br />
ε τ e<br />
e<br />
e<br />
e<br />
(3.48)<br />
24 1920 322560<br />
on<strong>de</strong> as or<strong>de</strong>ns verda<strong>de</strong>iras do erro <strong>de</strong> truncamento são dadas por<br />
pV<br />
= 2,4,6, K<br />
, sendo a<br />
or<strong>de</strong>m assintótica pL = 2 .<br />
3.4 SOLUÇÃO NUMÉRICA<br />
A proprieda<strong>de</strong> física consi<strong>de</strong>rada para as sete variáveis <strong>de</strong> interesse estudadas pelas<br />
equações <strong>de</strong> Poisson e advecção-difusão foi a temperatura ( T ) e pela equação <strong>de</strong> Burgers foi a<br />
velocida<strong>de</strong> ( u ). A Tab. 3.4 mostra a simbologia atribuída para as variáveis <strong>de</strong> interesse.<br />
Tabela 3.4: Simbologia atribuída às variáveis <strong>de</strong> interesse<br />
VARIÁVEL DE INTERESSE<br />
SÍMBOLO<br />
1) φ em x = 1 2 (valor nodal) T<br />
nod<br />
u<br />
nod<br />
2) φ em x = 1 2 (média arit.) T<br />
med<br />
u<br />
med<br />
3) φ com a regra do retângulo m ret<br />
T ,<br />
u<br />
m , ret<br />
4) φ com a regra do trapézio m trap<br />
T ,<br />
u<br />
m , trap<br />
5) l<br />
1<br />
do erro <strong>de</strong> discretização<br />
E<br />
m<br />
6) dx<br />
dφ<br />
em x = 0 com DDS I<br />
DDS<br />
7) dx<br />
dφ<br />
em x = 0 com DDS-2 I<br />
DDS−2