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 partir da equação (3.45) e sem se utilizar os termos <strong>de</strong> filtragem, a primeira<br />
linha <strong>de</strong> (4.22) resulta:<br />
∂<br />
∂x<br />
a+h<br />
a<br />
2 µ xx ∂u ∂<br />
t dz +<br />
∂x ∂y<br />
a+h<br />
a<br />
2 µ xy<br />
t<br />
∂u<br />
∂y<br />
<br />
∂v<br />
+<br />
∂x<br />
49<br />
dz (4.23)<br />
que, no caso do mo<strong>de</strong>lo RMA2, on<strong>de</strong> consi<strong>de</strong>ra-se µt homogêneo e isotrópico, é ajustada<br />
para:<br />
h<br />
ρ µt<br />
∂2U h<br />
+<br />
∂x2 ρ µt<br />
A tensão na superfície livre, gerada pelo vento, é dada por:<br />
∂2U . (4.24)<br />
∂y2 τ S x = ρ ζ V 2 w cos(ψ) , (4.25)<br />
on<strong>de</strong> ζ é <strong>um</strong> coeficiente empírico, Vw é a velocida<strong>de</strong> do vento e ψ é direção do vento medida<br />
no sentido anti-horário a partir do eixo x. A tensão junto ao fundo po<strong>de</strong> ser avaliada pela<br />
aplicação <strong>de</strong> <strong>um</strong>a fórmula empírica, como a <strong>de</strong> Manning, da seguinte forma:<br />
τ F = ρ g Rh<br />
<br />
n2V 2<br />
u+v<br />
R 4/3<br />
<br />
h<br />
on<strong>de</strong> Rh é o raio hidráulico, n o coeficiente <strong>de</strong> rugosida<strong>de</strong> <strong>de</strong> Manning, e Vu+v a velocida<strong>de</strong><br />
resultante da soma <strong>de</strong> U e V , cujo módulo é √ U 2 + V 2 . Consi<strong>de</strong>rando o raio hidráulico<br />
igual à profundida<strong>de</strong> h e <strong>de</strong>nominando β o ângulo entre Vu+v e U, a tensão junto ao fundo<br />
na direção x será dada por:<br />
τ F x = τ F 2 2 n V<br />
cos(β) = ρ g h<br />
u+v<br />
h 4/3<br />
U<br />
Vu+v<br />
,<br />
= ρ g n2 U<br />
h1/3 √<br />
U 2 + V 2 . (4.26)<br />
Substituindo os termos já integrados na equação (4.15) obtém-se a forma final<br />
da equação da quantida<strong>de</strong> <strong>de</strong> movimento na direção x integrada na profundida<strong>de</strong>, como