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 ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
56<br />
Estas incertezas são obtidas através do uso <strong>de</strong> alguns estimadores <strong>de</strong> erros. Como<br />
aparece <strong>em</strong> Marchi (2001), exist<strong>em</strong> vários métodos para estes cálculos, porém, po<strong>de</strong>m ser<br />
divididos <strong>em</strong> dois grupos, o primeiro estima o erro com base na solução numérica calculada<br />
sobre uma única malha, o segundo grupo consi<strong>de</strong>ra a estimativa <strong>de</strong> erro sobre as soluções<br />
numéricas obtidas sobre diferentes malhas.<br />
2.12.3 Estimador <strong>de</strong> Richardson<br />
O estimador <strong>de</strong> Richardson nos dá a incerteza da solução numérica para casos <strong>em</strong> que<br />
a solução analítica não é conhecida, assim, a incerteza ( U<br />
RI<br />
) <strong>de</strong> uma solução numérica ()<br />
é<br />
obtida <strong>em</strong>pregando a Eq.(2.68), baseado numa solução analítica estimada ( <br />
) .<br />
U ( )<br />
<br />
(2.68)<br />
RI<br />
<br />
Po<strong>de</strong>-se obter o <br />
<br />
por meio da extrapolação <strong>de</strong> Richardson dada por,<br />
( 1<br />
2<br />
)<br />
<br />
<br />
1<br />
<br />
(2.69)<br />
( q p L<br />
1)<br />
on<strong>de</strong> <br />
1<br />
e <br />
2<br />
são as soluções numéricas obtidas com as malhas fina e grossa, cujo tamanho<br />
(h)<br />
dos volumes <strong>de</strong> controle é h<br />
1<br />
e h<br />
2, respectivamente, e q é a razão <strong>de</strong> refino da malha,<br />
dada por,<br />
h<br />
2<br />
q (2.70)<br />
h<br />
1<br />
Com a substituição da Eq. (2.69) na Eq. (2.68) o estimador <strong>de</strong> Richardson resulta <strong>em</strong>,<br />
U<br />
RI<br />
( 1<br />
2<br />
)<br />
( ) <br />
(2.71)<br />
p<br />
( q L<br />
1)<br />
Or<strong>de</strong>m Efetiva:<br />
A or<strong>de</strong>m efetiva é a or<strong>de</strong>m do erro conseguido a partir <strong>de</strong> uma solução analítica<br />
conhecida, mais duas soluções numéricas <strong>em</strong> malhas refinadas com uma razão <strong>de</strong> refino (q).