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