18.07.2014 Views

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

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

50<br />

2.8.2 Estimativas <strong>de</strong> Erros a Priori<br />

As estimativas <strong>de</strong> erro a priori são usadas para estimar a or<strong>de</strong>m assintótica ( pL ) do<br />

erro <strong>de</strong> discretização, proporcionando uma análise qualitativa do erro <strong>de</strong> discretização antes<br />

da obtenção da solução numérica (SUERO, 2006).<br />

O procedimento adotado para estimar o erro <strong>de</strong> discretização a priori é estimar o erro<br />

<strong>de</strong> truncamento ( ε τ<br />

) do mo<strong>de</strong>lo matemático por meio da série <strong>de</strong> Taylor (TANNEHILL et al.,<br />

1997) e, admitir que o erro <strong>de</strong> discretização tenha a mesma forma funcional da Eq. (2.29)<br />

quando o tamanho ( ∆ x)<br />

dos volumes <strong>de</strong> controle ten<strong>de</strong> a zero:<br />

pL<br />

( ) C ∆x<br />

E φ (para ∆x → 0)<br />

(2.30)<br />

= 1<br />

on<strong>de</strong> C<br />

1<br />

representa um coeficiente com valor constante, porém <strong>de</strong>sconhecido.<br />

Portanto, somente com a estimativa <strong>de</strong> erro a priori não é possível obter o valor do<br />

erro <strong>de</strong> discretização (φ ),<br />

E mas é possível obter o valor da or<strong>de</strong>m assintótica ( )<br />

avaliar qual o efeito produzido pela redução <strong>de</strong><br />

pL e também,<br />

∆ x sobre o erro ( E ). Então, conhecendo-se<br />

pL e, pelo menos, duas soluções numéricas em malhas diferentes ( x f<br />

= malha _ fina)<br />

( ∆ = malha _ grossa)<br />

tem-se:<br />

x g<br />

∆ e<br />

pL<br />

( ) C ∆<br />

E φ (2.31)<br />

f<br />

= 1<br />

x f<br />

pL<br />

( ) C ∆<br />

E φ (2.32)<br />

g<br />

= 1<br />

x g<br />

realizando a operação <strong>de</strong> divisão entre as duas equações, tem-se:<br />

E<br />

E<br />

( φ<br />

f<br />

)<br />

( φ )<br />

g<br />

pL<br />

C1∆x<br />

f<br />

=<br />

pL<br />

C1∆x<br />

g<br />

=<br />

⎛ ∆x<br />

⎜<br />

⎝ ∆x<br />

f<br />

g<br />

⎞<br />

⎟<br />

⎠<br />

pL<br />

(2.33)<br />

para o caso do refino da malha ser: ( ) 2<br />

∆ x = ∆ e pL = 2 , tem-se:<br />

f<br />

x g<br />

E<br />

E<br />

( φ<br />

f<br />

)<br />

( φ )<br />

g<br />

=<br />

⎛<br />

⎜<br />

⎜<br />

⎜<br />

⎜<br />

⎝<br />

( ∆x<br />

)<br />

2<br />

g<br />

( ∆x<br />

)<br />

g<br />

⎞<br />

⎟<br />

⎟<br />

⎟<br />

⎟<br />

⎠<br />

2<br />

=<br />

⎛ 1 ⎞<br />

⎜ ⎟⎠<br />

⎝ 2<br />

2<br />

= 4<br />

1<br />

(2.34)

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

Saved successfully!

Ooh no, something went wrong!