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

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

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

Saved successfully!

Ooh no, something went wrong!