17.04.2013 Views

Modelagem Física e Computacional de um Escoamento Fluvial

Modelagem Física e Computacional de um Escoamento Fluvial

Modelagem Física e Computacional de um Escoamento Fluvial

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

4.4 INTEGRAÇÃO NA DIREÇÃO z DAS EQUAÇÕES DA CONSERVAÇÃO DA<br />

QUANTIDADE DE MOVIMENTO NAS DIREÇÕES x E y<br />

Partindo da equação (3.44) para a direção x (i = 1), inicialmente sem incluir a<br />

força <strong>de</strong> Coriolis, obtém-se:<br />

∂u<br />

∂t<br />

+ ∂u2<br />

∂x<br />

+ ∂uv<br />

∂y<br />

+ ∂uw<br />

∂z<br />

= Fx<br />

ρ<br />

1 ∂p<br />

−<br />

ρ ∂x<br />

<br />

1 ∂τxx ∂τxy<br />

+ +<br />

ρ ∂x ∂y<br />

<br />

∂τxz<br />

+<br />

∂z<br />

47<br />

. (4.15)<br />

A integração em z dos quatro termos do primeiro membro, aplicando a regra <strong>de</strong> Leibnitz,<br />

fornece:<br />

∂<br />

∂t<br />

a+h<br />

a<br />

∂a<br />

u dz + u| a ∂t − u| a+h<br />

+ ∂<br />

∂x<br />

a+h<br />

a<br />

+ ∂<br />

∂y<br />

u 2 dz + u 2 a<br />

a+h<br />

a<br />

∂(a + h)<br />

∂t +<br />

∂a<br />

∂x − u2 ∂(a + h)<br />

a+h ∂x +<br />

uv dz + uv| a<br />

∂a<br />

∂y − uv| ∂(a + h)<br />

a+h ∂v +<br />

− uw| a + uw| a+h .<br />

Os dois últimos termos <strong>de</strong> cada linha anulam-se através da aplicação das condições <strong>de</strong> con-<br />

torno cinemáticas na superfície livre (equação 4.8) e no fundo (equação 4.9), resultando 5 ,<br />

através <strong>de</strong> (4.11):<br />

O termo Fx<br />

ρ<br />

h ∂U<br />

∂t<br />

+ hU ∂U<br />

∂x<br />

+ hV ∂U<br />

∂y<br />

. (4.16)<br />

<strong>de</strong> (4.15) representa as forças por unida<strong>de</strong> <strong>de</strong> massa na direção x.<br />

Em alguns casos <strong>de</strong> escoamentos <strong>de</strong> gran<strong>de</strong>s dimensões po<strong>de</strong> tornar-se relevante a força<br />

<strong>de</strong> Coriolis:<br />

F Coriolis<br />

x = 2 ρ h Ω V sen(φ) , (4.17)<br />

em que Ω é a velocida<strong>de</strong> <strong>de</strong> rotação angular da Terra e φ a latitu<strong>de</strong> local.<br />

5 A rigor, nesta <strong>de</strong>dução restariam os termos com correlações <strong>de</strong> <strong>de</strong>svios <strong>de</strong> velocida<strong>de</strong>s médias na<br />

vertical ũ 2 dz e ũ˜v dz, que são assimilados pelo mo<strong>de</strong>lo <strong>de</strong> turbulência.

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

Saved successfully!

Ooh no, something went wrong!