Serie3 FMC

Fazer download em pdf ou txt
Fazer download em pdf ou txt
Você está na página 1de 2

UNIVERSIDADE DE LISBOA - FCUL

FÍSICA DOS MEIOS CONTÍNUOS


Problemas - 3a Série

1. a) Mostre que para um fluido incompressı́vel a equação de Navier-Stokes


Du
= −∇p∗ − η∇ × Ω
Dt
se pode escrever
Du
= −∇p∗ + η∇2 u.
Dt
b) Mostre que um tensor de ordem dois se pode escrever
 
1 1 2 1
Tij = δij Tmm + Tij + Tji − δij Tmm + (Tij − Tji ),
3 2 3 2
onde o primeiro tensor é isotrópico, o segundo é simétrico e o terceiro anti-simétrico. Escreva
os tensores das tensões e das taxas de deformação nesta forma e discuta o significado fı́sico de
cada um deles.
2. Mostre que a dissipação viscosa num fluido de Navier-Stokes em escoamento se escreve (ver Sec.
6.5 do Faber):
 2
η ∂ui ∂uj
+ .
2 ∂xj ∂xi
3. Considere um plano com inclinação α com a horizontal onde temos uma camada de fluido de
espessura d. O fluido, de viscosidade η, tem uma superfı́cie livre e está sujeito a um campo
gravitacional g. a) Calcule a velocidade em função da distância da placa e a taxa de escoamento
para uma espessura constante e uniforme d. b) Discuta as aproximações necessárias para que se
possa considerar o escoamento como estacionário.
4. Considere um escoamento viscoso para baixo num tubo de seção transversal circular com r = a
sob gradiente de pressão constante P = −dp/dz. a) Mostre que
P 2
uz = (a − r2 ), ur = uθ = 0.

b) Calcule a taxa de escoamento no tubo. Este escoamento é conhecido com Hagen-Poiseuille. A
instabilidade formada em números de Reynolds mais elevados constitui um problema importante
em mecânica dos fluidos (veja, por exemplo, www.youtube.com/watch?v=y0WRJtXvpSo).
   
∂u ∂ui ∂ui ∂uj
5. a) Usando ti = Tij nj e Tij = −pδij + η ∂xji + ∂xj
, obtemos ti = −pn i + ηni ∂xj + ∂xi
, onde
ti é a componente i da tensão no elemento de superfı́cie com normal ni . Mostre que

t = −pn + η[2(n · ∇)u + n × (∇ × u)].

b) Use o resultado da alı́nea “a” e as identidades matemáticas necessárias para mostrar que a
força resultante exercida num volume de fluido pelo fluido circundante é
Z Z
tdS = (−∇p + η∇2 u)dV,
S V

1
onde S é a superfı́cie do volume de fluido. Deduza que se o volume for pequeno, a força resultante
excluindo a gravidade é −∇p + η∇2 u por unidade de volume, em acordo com a equação de
Navier-Stokes para fluidos incompressı́veis.

6. a) Verifique que no caso de um escoamento u = [u(y), 0, 0], a tensão se reduz a


 
du
t = η , −p, 0
dy

quando n = (0, 1, 0). b) Mostre que os termos η(∂uj /∂xi + ∂ui /∂xj ) do tensor das tensões são
nulos para um escoamento uniformemente circular u = Ω × x, onde Ω é um vetor constate.

7. Considere um fluido inicialmente em repouso entre duas placas paralelas: uma em y = d/2 com
velocidade constante U e outra em y = −d/2 estacionária. a) Calcule o campo de velocidades.
b) Calcule a taxa de escoamento. c) Calcule a tensão de corte nas duas placas.

8. Um fluido escoa lentamente através de uma bolha esférica de raio a constituı́da de ar. Considere
que a tensão tangencial é nula, tθ = 0, em r = a. a) Mostre que a componente normal do stress
em r = a é tr = (3ηU/a) cos θ. b) Mostre que a força de arrasto na bolha é D = 4πηU a na
direção do escoamento livre. c) Considere agora que, no lugar da bolha de ar, temos uma gota
constituı́da de outro fluido de viscosidade η̄. Mostre que a força de arrasto é
 
η + 3η̄/2
D = 4πηU a .
η + η̄

d) Discuta os limites η̄/η −→ 0 e η̄/η −→ ∞.

9. O problema de um escoamento viscoso 2D através de um cilindro circular de raio a envolve


calcular o campo de velocidades u = [u(x, y), v(x, y), 0] que satisfaz

∇p
(u · ∇)u = − + ν∇2 u, ∇ · u,
ρ
com as condições de fronteira

u = 0 em x2 + y 2 = a2 ; u −→ (U, 0, 0) em x2 + y 2 −→ ∞.

Reescreva esse problema na forma adimensional usando as variáveis adimensionais

x0 = x/a, u0 = u/U, p0 = p/(ρU 2 )

em vez de x, u e p. Sem tentar resolver o problema, mostre que o padrão de linhas de corrente
depende de ν, a e U apenas através da combinação Re = U a/ν. Usando o código fornecido
(code2-cilinder.py), verifique que escoamentos com o mesmo número de Reynolds são geometri-
camente similares (varie estas 3 quantidades de forma a manter Re constante).

10. a) Calcule a velocidade u(y) num escoamento viscoso entre duas placas paralelas estacionárias em
y = 0 e y = L onde o fluido, de viscosidade cinemática ν = η/ρ, escoa na direção x devido a uma
força (e.g., gravidade) f = f ĥ. b) Compare a solução analı́tica com a obtida numericamente
usando o código fornecido (code3-Poiseuille.py) para os mesmos parâmetros. Discuta o que
acontece quando se altera a separação entre as placas L e a viscosidade ν de forma que L/ν
permaneça constante.

Você também pode gostar