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