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
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
a:<br />
∂ui<br />
∂xi<br />
41<br />
= 0 . (3.42)<br />
Como essa equação é linear, a aplicação <strong>de</strong> qualquer dos métodos <strong>de</strong> promediação<br />
ou filtragem é trivial, resultando em:<br />
∂ui<br />
∂xi<br />
= 0 . (3.43)<br />
A equação acima expressa que a qualquer instante a taxa líquida <strong>de</strong> fluxos <strong>de</strong><br />
entrada e <strong>de</strong> saída n<strong>um</strong> dado vol<strong>um</strong>e <strong>de</strong> controle fixo no espaço é sempre nula (continuida<strong>de</strong><br />
<strong>de</strong> vol<strong>um</strong>e).<br />
3.5.2.2 Equação da conservação da quantida<strong>de</strong> <strong>de</strong> movimento<br />
Antes <strong>de</strong> escrever a equação <strong>de</strong> Navier-Stokes, correspon<strong>de</strong>nte à segunda Lei<br />
<strong>de</strong> Newton n<strong>um</strong> referencial Euleriano, é necessário explicitar <strong>um</strong>a aproximação básica, a<br />
chamada aproximação <strong>de</strong> Boussinesq. Em corpos <strong>de</strong> água rasos, os gradientes <strong>de</strong> <strong>de</strong>nsida<strong>de</strong><br />
são tipicamente da or<strong>de</strong>m <strong>de</strong> <strong>um</strong> por mil, por isso, a substituição da massa específica<br />
real ρ(xi, t) por <strong>um</strong> valor <strong>de</strong> referência constante ρ0 não apresenta problema, exceto no<br />
termo <strong>de</strong> pressão. No termo <strong>de</strong> pressão, conforme comenta CSANADY (1982), embora os<br />
gradientes <strong>de</strong> <strong>de</strong>nsida<strong>de</strong> sejam pequenos, muitas vezes não são <strong>de</strong>sprezíveis. Isso posto,<br />
po<strong>de</strong>-se escrever a seguinte equação do movimento para o escoamento <strong>de</strong> gran<strong>de</strong> escala<br />
ou resolvível:<br />
∂ui<br />
∂t<br />
∂<br />
+ (ui uj) = −<br />
∂xj<br />
1 ∂p<br />
− g δi3 + 2 ɛijk uj Ωk +<br />
ρ0 ∂xi<br />
1 ∂τ ij<br />
, (3.44)<br />
ρ0 ∂xj<br />
on<strong>de</strong> p é o termo <strong>de</strong> pressão. No termo representando as forças <strong>de</strong> Coriolis, Ωk são<br />
as componentes da velocida<strong>de</strong> angular da Terra n<strong>um</strong> sistema <strong>de</strong> coor<strong>de</strong>nadas locais. A<br />
magnitu<strong>de</strong> <strong>de</strong> Ωk é <strong>de</strong> 2π radianos em 24 horas, ou Ω ∼ = 0, 7292 × 10 −4 s −1 .