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.
27<br />
el<strong>em</strong>entar, o menor volume <strong>de</strong> fluido ainda tratado como um meio contínuo, ou seja, que não<br />
é afetado pelas ativida<strong>de</strong>s moleculares / atômicas do fluido.<br />
Po<strong>de</strong>-se obter as equações discretizadas <strong>de</strong> duas formas através do método dos<br />
volumes finitos. A primeira é através <strong>de</strong> um balanço da proprieda<strong>de</strong> no volume el<strong>em</strong>entar. Já<br />
a segunda é por meio <strong>de</strong> integração, sobre o volume <strong>de</strong> controle (CV), no espaço e no t<strong>em</strong>po,<br />
das equações na forma conservativa. Utilizando a integração sobre o volume el<strong>em</strong>entar e<br />
consi<strong>de</strong>rando a Eq.(2.1), obtém-se (VERSTEEG e MALALASEKERA, 1995),<br />
<br />
CV<br />
(<br />
)<br />
dV <br />
t<br />
<br />
CV<br />
<br />
<br />
( V ).<br />
dV<br />
<br />
<br />
CV<br />
<br />
<br />
( )<br />
dV <br />
<br />
CV<br />
S<br />
<br />
dV<br />
(2.2)<br />
Aplicando o teor<strong>em</strong>a da divergência <strong>de</strong> Gauss na Eq. (2.2), chega-se a<br />
<br />
t<br />
<br />
CV<br />
<br />
dV n<br />
( V<br />
).<br />
dA n<br />
( )<br />
dA S dV<br />
(2.3)<br />
ACV<br />
<br />
<br />
<br />
ACV<br />
<br />
A Fig.2.1 mostra o significado físico <strong>de</strong> cada um dos termos <strong>de</strong>sta equação.<br />
<br />
<br />
CV<br />
<br />
Taxa <strong>de</strong><br />
acréscimo<br />
<strong>de</strong> <br />
Taxa <strong>de</strong><br />
<strong>de</strong>scréscimo <strong>de</strong> <br />
<strong>de</strong>vido à advecção<br />
através das<br />
fronteiras<br />
Taxa <strong>de</strong><br />
acréscimo <strong>de</strong> <br />
<strong>de</strong>vido à<br />
difusão através<br />
das fronteiras<br />
Taxa <strong>de</strong><br />
criação /<br />
<strong>de</strong>struição<br />
<strong>de</strong> <br />
Figura 2.1 – Significado físico dos termos da Eq. (2.3)<br />
2.3 DISCRETIZAÇÃO<br />
Para discretização da EDP <strong>de</strong>ve-se primeiramente discretizar o espaço, ou domínio,<br />
que é a região <strong>em</strong> que a mo<strong>de</strong>lag<strong>em</strong> física foi feita. Isto é realizado dividindo-o <strong>em</strong> pontos, ou<br />
volumes <strong>de</strong> controle (CV), processo representado pela Fig.2.2.<br />
A divisão da região, ou discretização, gera uma malha <strong>de</strong> pontos/volumes <strong>de</strong><br />
controle. Os pontos, ou nós, são os locais on<strong>de</strong> serão <strong>de</strong>finidas as variáveis <strong>de</strong> interesse, por<br />
isto é importante que seja feita uma análise pormenorizada na geração da malha.<br />
Verifica-se <strong>em</strong> Schnei<strong>de</strong>r (2007) que a geração da malha po<strong>de</strong> ser <strong>de</strong> quatro formas<br />
básicas: uniforme, não-uniforme, estruturada e não-estruturada.