Origem: Wikilivros, livros abertos por um mundo aberto.
Um líquido ideal é um fluido incompressível e invíscido (ou seja, sua viscosidade é nula). Muitos líquidos de interesse têm um comportamento que pode ser razoavelmente aproximado pelo líquido ideal. As equações básicas apresentam uma forma muito simplificada quando comparadas às equações gerais.
Portanto, para um líquido ideal, teremos
ρ
=
c
o
n
s
t
a
n
t
e
μ
=
0
{\displaystyle \rho \;=\;constante\qquad \mu \;=\;0}
Vamos aplicar essa simplificação às equações básicas em suas diversas formas (para sistemas, integral e diferencial) e obter as equações resultantes. É importante que, como estamos tratando de líquidos, adotaremos a convenção de fazer o eixo Z apontar para baixo, e não para cima. Dessa maneira, o eixo X aponta para a direita e o eixo Y aponta para o observador.
A equação de continuidade em forma diferencial
∂
ρ
∂
t
+
(
∇
⋅
(
ρ
v
→
)
)
=
0
{\displaystyle {\frac {\partial \rho }{\partial t}}\;+\;(\nabla \cdot (\rho {\vec {v}}))\;=\;0}
quando aplicada a um líquido ideal, simplifica-se para
∇
⋅
v
→
=
0
{\displaystyle \nabla \cdot {\vec {v}}\;=\;0}
Conforme visto anteriormente , as equações de Navier-Stokes
(
∂
∂
x
(
−
p
−
2
3
μ
(
∂
v
x
∂
x
+
∂
v
y
∂
y
+
∂
v
z
∂
z
)
+
2
μ
∂
v
x
∂
x
)
+
∂
∂
y
(
μ
∂
v
x
∂
y
+
μ
∂
v
y
∂
x
)
+
∂
∂
z
(
μ
∂
v
x
∂
z
+
μ
∂
v
z
∂
x
)
)
=
{\displaystyle \left({\frac {\partial }{\partial x}}\left(-\;p\;-\;{\frac {2}{3}}\mu \left({\frac {\partial v_{x}}{\partial x}}\;+\;{\frac {\partial v_{y}}{\partial y}}\;+\;{\frac {\partial v_{z}}{\partial z}}\right)\;+\;2\mu {\frac {\partial v_{x}}{\partial x}}\right)\;+\;{\frac {\partial }{\partial y}}\left(\mu {\frac {\partial v_{x}}{\partial y}}\;+\;\mu {\frac {\partial v_{y}}{\partial x}}\right)\;+\;{\frac {\partial }{\partial z}}\left(\mu {\frac {\partial v_{x}}{\partial z}}\;+\;\mu {\frac {\partial v_{z}}{\partial x}}\right)\right)\;=\;}
=
ρ
(
v
x
∂
v
x
∂
x
+
v
y
∂
v
x
∂
y
+
v
z
∂
v
x
∂
z
+
∂
v
x
∂
t
)
{\displaystyle \;=\;\rho \left(v_{x}{\frac {\partial v_{x}}{\partial x}}\;+\;v_{y}{\frac {\partial v_{x}}{\partial y}}\;+\;v_{z}{\frac {\partial v_{x}}{\partial z}}\;+\;{\frac {\partial v_{x}}{\partial t}}\right)}
(
∂
∂
y
(
−
p
−
2
3
μ
(
∂
v
x
∂
x
+
∂
v
y
∂
y
+
∂
v
z
∂
z
)
+
2
μ
∂
v
y
∂
y
)
+
∂
∂
x
(
μ
∂
v
x
∂
y
+
μ
∂
v
y
∂
x
)
+
∂
∂
z
(
μ
∂
v
z
∂
y
+
μ
∂
v
y
∂
z
)
)
=
{\displaystyle \left({\frac {\partial }{\partial y}}\left(-\;p\;-\;{\frac {2}{3}}\mu \left({\frac {\partial v_{x}}{\partial x}}\;+\;{\frac {\partial v_{y}}{\partial y}}\;+\;{\frac {\partial v_{z}}{\partial z}}\right)\;+\;2\mu {\frac {\partial v_{y}}{\partial y}}\right)\;+\;{\frac {\partial }{\partial x}}\left(\mu {\frac {\partial v_{x}}{\partial y}}\;+\;\mu {\frac {\partial v_{y}}{\partial x}}\right)\;+\;{\frac {\partial }{\partial z}}\left(\mu {\frac {\partial v_{z}}{\partial y}}\;+\;\mu {\frac {\partial v_{y}}{\partial z}}\right)\right)\;=\;}
=
ρ
(
v
x
∂
v
y
∂
x
+
v
y
∂
v
y
∂
y
+
v
z
∂
v
y
∂
z
+
∂
v
y
∂
t
)
{\displaystyle \;=\;\rho \left(v_{x}{\frac {\partial v_{y}}{\partial x}}\;+\;v_{y}{\frac {\partial v_{y}}{\partial y}}\;+\;v_{z}{\frac {\partial v_{y}}{\partial z}}\;+\;{\frac {\partial v_{y}}{\partial t}}\right)}
(
∂
∂
z
(
−
p
−
2
3
μ
(
∂
v
x
∂
x
+
∂
v
y
∂
y
+
∂
v
z
∂
z
)
+
2
μ
∂
v
z
∂
z
)
+
∂
∂
x
(
μ
∂
v
x
∂
z
+
μ
∂
v
z
∂
x
)
+
∂
∂
y
(
μ
∂
v
z
∂
y
+
μ
∂
v
y
∂
z
)
)
+
ρ
g
=
{\displaystyle \left({\frac {\partial }{\partial z}}\left(-\;p\;-\;{\frac {2}{3}}\mu \left({\frac {\partial v_{x}}{\partial x}}\;+\;{\frac {\partial v_{y}}{\partial y}}\;+\;{\frac {\partial v_{z}}{\partial z}}\right)\;+\;2\mu {\frac {\partial v_{z}}{\partial z}}\right)\;+\;{\frac {\partial }{\partial x}}\left(\mu {\frac {\partial v_{x}}{\partial z}}\;+\;\mu {\frac {\partial v_{z}}{\partial x}}\right)\;+\;{\frac {\partial }{\partial y}}\left(\mu {\frac {\partial v_{z}}{\partial y}}\;+\;\mu {\frac {\partial v_{y}}{\partial z}}\right)\right)\;+\;\rho g\;=\;}
=
ρ
(
v
x
∂
v
z
∂
x
+
v
y
∂
v
z
∂
y
+
v
z
∂
v
z
∂
z
+
∂
v
z
∂
t
)
{\displaystyle \;=\;\rho \left(v_{x}{\frac {\partial v_{z}}{\partial x}}\;+\;v_{y}{\frac {\partial v_{z}}{\partial y}}\;+\;v_{z}{\frac {\partial v_{z}}{\partial z}}\;+\;{\frac {\partial v_{z}}{\partial t}}\right)}
derivadas do princípio de conservação do momento linear, descrevem a dinâmica de um volume diferencial de fluido, juntamente com a equação de continuidade, mas são extremamente difíceis de resolver. No caso de um líquido ideal, as equações acima se reduzem a
−
∂
p
∂
x
=
ρ
(
v
x
∂
v
x
∂
x
+
v
y
∂
v
x
∂
y
+
v
z
∂
v
x
∂
z
+
∂
v
x
∂
t
)
{\displaystyle -\;{\frac {\partial p}{\partial x}}\;=\;\rho \left(v_{x}{\frac {\partial v_{x}}{\partial x}}\;+\;v_{y}{\frac {\partial v_{x}}{\partial y}}\;+\;v_{z}{\frac {\partial v_{x}}{\partial z}}\;+\;{\frac {\partial v_{x}}{\partial t}}\right)}
−
∂
p
∂
y
=
ρ
(
v
x
∂
v
x
∂
x
+
v
y
∂
v
x
∂
y
+
v
z
∂
v
x
∂
z
+
∂
v
x
∂
t
)
{\displaystyle -\;{\frac {\partial p}{\partial y}}\;=\;\rho \left(v_{x}{\frac {\partial v_{x}}{\partial x}}\;+\;v_{y}{\frac {\partial v_{x}}{\partial y}}\;+\;v_{z}{\frac {\partial v_{x}}{\partial z}}\;+\;{\frac {\partial v_{x}}{\partial t}}\right)}
−
∂
p
∂
z
+
ρ
g
=
ρ
(
v
x
∂
v
x
∂
x
+
v
y
∂
v
x
∂
y
+
v
z
∂
v
x
∂
z
+
∂
v
x
∂
t
)
{\displaystyle -\;{\frac {\partial p}{\partial z}}\;+\;\rho g\;=\;\rho \left(v_{x}{\frac {\partial v_{x}}{\partial x}}\;+\;v_{y}{\frac {\partial v_{x}}{\partial y}}\;+\;v_{z}{\frac {\partial v_{x}}{\partial z}}\;+\;{\frac {\partial v_{x}}{\partial t}}\right)}
Que podem ser escritas também na forma mais concisa
ρ
g
→
−
∇
p
=
ρ
D
v
→
D
t
{\displaystyle \rho {\vec {g}}-\nabla p\;=\;\rho {\frac {D{\vec {v}}}{Dt}}}
Essas equações são chamadas de equações de Euler .
As equações de Euler em coordenadas cilíndricas podem ser obtidas a partir da forma concisa, substituindo-se a definição da derivada direcional nesse sistema
ρ
g
→
−
∇
p
=
ρ
D
v
→
D
t
=
ρ
(
∂
v
→
∂
t
+
v
→
⋅
∇
v
→
)
{\displaystyle \rho {\vec {g}}\;-\;\nabla p\;=\;\rho {\frac {D{\vec {v}}}{Dt}}\;=\;\rho \left({\frac {\partial {\vec {v}}}{\partial t}}\;+\;{\vec {v}}\cdot \nabla {\vec {v}}\right)}
ρ
g
r
−
∂
p
∂
r
=
ρ
(
∂
v
r
∂
t
+
v
r
∂
v
r
∂
r
+
v
θ
r
∂
v
r
∂
θ
+
v
z
∂
v
r
∂
z
−
v
θ
2
r
)
{\displaystyle \rho g_{r}\;-\;{\frac {\partial p}{\partial r}}\;=\;\rho \left({\frac {\partial v_{r}}{\partial t}}\;+\;v_{r}\;{\frac {\partial v_{r}}{\partial r}}\;+\;{\frac {v_{\theta }}{r}}\;{\frac {\partial v_{r}}{\partial \theta }}\;+\;v_{z}\;{\frac {\partial v_{r}}{\partial z}}\;-\;{\frac {v_{\theta }^{2}}{r}}\right)}
ρ
g
θ
−
1
r
∂
p
∂
θ
=
ρ
(
∂
v
θ
∂
t
+
v
r
∂
v
θ
∂
r
+
v
θ
r
∂
v
θ
∂
θ
+
v
z
∂
v
θ
∂
z
+
v
r
v
θ
r
)
{\displaystyle \rho g_{\theta }\;-\;{\frac {1}{r}}\;{\frac {\partial p}{\partial \theta }}\;=\;\rho \left({\frac {\partial v_{\theta }}{\partial t}}\;+\;v_{r}\;{\frac {\partial v_{\theta }}{\partial r}}\;+\;{\frac {v_{\theta }}{r}}\;{\frac {\partial v_{\theta }}{\partial \theta }}\;+\;v_{z}\;{\frac {\partial v_{\theta }}{\partial z}}\;+\;{\frac {v_{r}v_{\theta }}{r}}\right)}
ρ
g
z
−
∂
p
∂
z
=
ρ
(
∂
v
z
∂
t
+
v
r
∂
v
z
∂
r
+
v
θ
r
∂
v
z
∂
θ
+
v
z
∂
v
z
∂
z
)
{\displaystyle \rho g_{z}\;-\;{\frac {\partial p}{\partial z}}\;=\;\rho \left({\frac {\partial v_{z}}{\partial t}}\;+\;v_{r}\;{\frac {\partial v_{z}}{\partial r}}\;+\;{\frac {v_{\theta }}{r}}\;{\frac {\partial v_{z}}{\partial \theta }}\;+\;v_{z}\;{\frac {\partial v_{z}}{\partial z}}\right)}
As equações de Euler em coordenadas esféricas podem ser obtidas da mesma forma:
ρ
g
→
−
∇
p
=
ρ
(
∂
v
→
∂
t
+
v
→
⋅
∇
v
→
)
{\displaystyle \rho {\vec {g}}\;-\;\nabla p\;=\;\rho \left({\frac {\partial {\vec {v}}}{\partial t}}\;+\;{\vec {v}}\cdot \nabla {\vec {v}}\right)}
ρ
g
r
−
∂
p
∂
r
=
ρ
(
∂
v
r
∂
t
+
v
r
∂
v
r
∂
r
+
v
θ
r
∂
v
r
∂
θ
+
v
ϕ
r
sin
(
θ
)
∂
v
r
∂
ϕ
−
v
ϕ
2
+
v
θ
2
r
)
{\displaystyle \rho g_{r}\;-\;{\frac {\partial p}{\partial r}}\;=\;\rho \left({\frac {\partial v_{r}}{\partial t}}\;+\;v_{r}\;{\frac {\partial v_{r}}{\partial r}}\;+\;{\frac {v_{\theta }}{r}}\;{\frac {\partial v_{r}}{\partial \theta }}\;+\;{\frac {v_{\phi }}{r\;\sin(\theta )}}\;{\frac {\partial v_{r}}{\partial \phi }}\;-\;{\frac {v_{\phi }^{2}+v_{\theta }^{2}}{r}}\right)}
ρ
g
θ
−
1
r
∂
p
∂
θ
=
ρ
(
∂
v
θ
∂
t
+
v
r
∂
v
θ
∂
r
+
v
θ
r
∂
v
θ
∂
θ
+
v
ϕ
r
sin
(
θ
)
∂
v
θ
∂
ϕ
+
v
r
v
θ
−
v
ϕ
2
cot
(
θ
)
r
)
{\displaystyle \rho g_{\theta }\;-\;{\frac {1}{r}}\;{\frac {\partial p}{\partial \theta }}\;=\;\rho \left({\frac {\partial v_{\theta }}{\partial t}}\;+\;v_{r}{\frac {\partial v_{\theta }}{\partial r}}\;+\;{\frac {v_{\theta }}{r}}\;{\frac {\partial v_{\theta }}{\partial \theta }}\;+\;{\frac {v_{\phi }}{r\;\sin(\theta )}}\;{\frac {\partial v_{\theta }}{\partial \phi }}\;+\;{\frac {v_{r}v_{\theta }\;-\;v_{\phi }^{2}\cot(\theta )}{r}}\right)}
ρ
g
ϕ
−
1
r
sin
(
θ
)
∂
p
∂
ϕ
=
ρ
(
∂
v
ϕ
∂
t
+
v
r
∂
v
ϕ
∂
r
+
v
θ
r
∂
v
ϕ
∂
θ
+
v
ϕ
r
sin
(
θ
)
∂
v
ϕ
∂
ϕ
+
v
r
v
ϕ
+
v
ϕ
v
θ
cot
(
θ
)
r
)
{\displaystyle \rho g_{\phi }\;-\;{\frac {1}{r\;\sin(\theta )}}\;{\frac {\partial p}{\partial \phi }}\;=\;\rho \left({\frac {\partial v_{\phi }}{\partial t}}\;+\;v_{r}{\frac {\partial v_{\phi }}{\partial r}}\;+\;{\frac {v_{\theta }}{r}}\;{\frac {\partial v_{\phi }}{\partial \theta }}\;+\;{\frac {v_{\phi }}{r\;\sin(\theta )}}\;{\frac {\partial v_{\phi }}{\partial \phi }}\;+\;{\frac {v_{r}v_{\phi }\;+\;v_{\phi }v_{\theta }\cot(\theta )}{r}}\right)}
Equações de Euler em coordenadas de linhas de corrente [ editar | editar código-fonte ]
As equações de Euler em coordenadas de linhas de corrente podem ser obtidas a partir da forma concisa tomando-se um elemento de volume infinitesimal δV e aplicando-se o princípio de conservação do momento linear na direção do fluxo, após considerar ρ constante e μ = 0. Seja o fluxo normal ao eixo Y e Θ, o ângulo que a linha de corrente faz com a horizontal nesse ponto. A pressão pl- na face posterior do volume será
p
l
−
=
p
−
δ
p
l
2
=
p
−
∂
p
∂
l
δ
l
2
{\displaystyle p_{l-}\;=\;p\;-\;{\frac {\delta p_{l}}{2}}\;=\;p\;-\;{\frac {{\frac {\partial p}{\partial l}}\;\delta l}{2}}}
onde p é a pressão no centro do volume. A pressão na face anterior do volume, similarmente, será
p
l
+
=
p
+
δ
p
l
2
=
p
+
∂
p
∂
l
δ
l
2
{\displaystyle p_{l+}\;=\;p\;+\;{\frac {\delta p_{l}}{2}}\;=\;p\;+\;{\frac {{\frac {\partial p}{\partial l}}\;\delta l}{2}}}
e a diferença entre as pressões será, evidentemente
Δ
p
l
=
p
l
+
−
p
l
−
=
∂
p
∂
l
δ
l
{\displaystyle \Delta p_{l}\;=\;p_{l+}\;-\;p_{l-}\;=\;{\frac {\partial p}{\partial l}}\;\delta l}
pressão esta cujo sentido é contrário ao fluxo. O peso do volume terá uma componente na direção do fluxo igual a δm · g · sen Θ, também no sentido contrário a ele. Assim, podemos escrever
−
δ
m
g
sin
θ
−
∂
p
∂
l
δ
l
d
A
l
=
δ
m
a
l
{\displaystyle -\;\delta m\;g\;\sin \theta \;-\;{\frac {\partial p}{\partial l}}\;\delta l\;dA_{l}\;=\;\delta m\;a_{l}}
onde Al é a área transversal ao fluxo e al é a aceleração do volume. Desenvolvendo, teremos
−
ρ
d
l
d
n
d
y
g
∂
z
∂
l
−
∂
p
∂
l
δ
l
d
n
d
y
=
ρ
d
l
d
n
d
y
a
l
⇒
ρ
a
l
=
−
∂
p
∂
l
−
ρ
g
∂
z
∂
l
{\displaystyle -\;\rho \;dl\;dn\;dy\;g\;{\frac {\partial z}{\partial l}}\;-\;{\frac {\partial p}{\partial l}}\;\delta l\;dn\;dy\;=\;\rho \;dl\;dn\;dy\;a_{l}\Rightarrow \;\;\;\rho a_{l}\;=\;-\;{\frac {\partial p}{\partial l}}\;-\;\rho g{\frac {\partial z}{\partial l}}}
Mas, sobre uma linha de corrente
a
l
=
D
v
D
t
=
∂
v
∂
t
+
v
∂
v
∂
l
{\displaystyle a_{l}\;=\;{\frac {Dv}{Dt}}\;=\;{\frac {\partial v}{\partial t}}\;+\;v{\frac {\partial v}{\partial l}}}
Assim
ρ
(
∂
v
∂
t
+
v
∂
v
∂
l
)
=
−
∂
p
∂
l
−
ρ
g
∂
z
∂
l
{\displaystyle \rho \left({\frac {\partial v}{\partial t}}\;+\;v{\frac {\partial v}{\partial l}}\right)\;=\;-\;{\frac {\partial p}{\partial l}}\;-\;\rho g{\frac {\partial z}{\partial l}}}
No caso de fluxo em regime permanente
ρ
v
∂
v
∂
l
=
−
∂
p
∂
l
−
ρ
g
∂
z
∂
l
{\displaystyle \rho v{\frac {\partial v}{\partial l}}\;=\;-\;{\frac {\partial p}{\partial l}}\;-\;\rho g{\frac {\partial z}{\partial l}}}
Procedendo de forma similar para a componente normal à linha de corrente
Δ
p
n
=
p
n
+
−
p
n
−
=
∂
p
∂
n
δ
n
{\displaystyle \Delta p_{n}\;=\;p_{n+}\;-\;p_{n-}\;=\;{\frac {\partial p}{\partial n}}\;\delta n}
−
δ
m
g
cos
θ
−
∂
p
∂
n
δ
n
d
A
n
=
δ
m
a
n
{\displaystyle -\;\delta m\;g\;\cos \theta \;-\;{\frac {\partial p}{\partial n}}\;\delta n\;dA_{n}\;=\;\delta m\;a_{n}}
Desenvolvendo, teremos
−
ρ
d
l
d
n
d
y
g
∂
z
∂
n
−
∂
p
∂
n
δ
n
d
l
d
y
=
ρ
d
l
d
n
d
y
a
n
⇒
ρ
a
n
=
−
∂
p
∂
n
−
ρ
g
∂
z
∂
n
{\displaystyle -\;\rho \;dl\;dn\;dy\;g\;{\frac {\partial z}{\partial n}}\;-\;{\frac {\partial p}{\partial n}}\;\delta n\;dl\;dy\;=\;\rho \;dl\;dn\;dy\;a_{n}\Rightarrow \;\;\;\rho a_{n}\;=\;-\;{\frac {\partial p}{\partial n}}\;-\;\rho g{\frac {\partial z}{\partial n}}}
Mas a aceleração na direção normal à linha de corrente é a aceleração centrípeta
a
n
=
−
v
2
r
{\displaystyle a_{n}\;=\;-\;{\frac {v^{2}}{r}}}
onde r é o raio de curvatura da linha de corrente no ponto considerado. O sinal negativo indica que o elemento de volume está sendo acelerado para dentro da curva. Assim
ρ
v
2
r
=
∂
p
∂
n
+
ρ
g
∂
z
∂
n
{\displaystyle \rho {\frac {v^{2}}{r}}\;=\;{\frac {\partial p}{\partial n}}\;+\;\rho g{\frac {\partial z}{\partial n}}}
Em regiões onde as linhas de corrente são linhas retas, r = ∞, o que implica em an = 0. Nessas regiões, não há variação de presssão entre as linhas de corrente, pois
∂
p
∂
n
=
0
{\displaystyle \;{\frac {\partial p}{\partial n}}\;=\;0}
.