Equações de Navier-Stokes

Ver Todos os Exercícios

Teoria

Fala aí! Seja bem vindo ao Responde Aí! Bora aprender a Equação de Navier-Stokes?!

Qual é a Equação de Navier-Stokes?

Sem enrolação, a equação de Navier-Stokes é essa daqui:

Equação de Navier-Stokes
Equação de Navier-Stokes

Sim, ela é é assustadora! Mas fica tranquilo, é só uma questão de tempo pra você pegar intimidade! Como falamos lá em cima, a Equação de Navier-Stokes nos descreve o escoamento de um fluido, então nada mais justo do que essa equação levar em conta DIVERSOS fatores, certo?!

Mas fica tranquilo, quando a gente partir pros exercícios você vai ver que é bem comum usarmos hipóteses simplificadas, então a equação fica bem menor e mais simples!

📢 Clique para ver mais:

Equação de Navier-Stokes: Coordenadas Retangulares

A equação que vimos lá em cima leva em conta os todos os eixos, que quando estamos falando de coordenadas retangulares são e . Ela é uma versão simplificada e mais elegante das Equações de Navier-Stokes!

Mas na prática, na hora dos exercícios a gente precisa abrir essa equação. A versão aberta da equação de Navier-Stokes nas coordenadas retangulares tem essa carinha aqui:

Equação de Navier-Stokes: Coordenadas Retangulares
Equação de Navier-Stokes: Coordenadas Retangulares

Como você pode ver, estão presente todos os termos que vimos na equação lá do início! Só que aqui você pode lidar com cada eixo de uma forma separada!

-Beleza, mas o que eu quero saber é como eu uso a equação de Navier-Stokes na prática! Como ela vai aparecer nos exercícios!

Seu pedido é uma ordem! Se liga nesse videozinho que o Responde Aí preparou pra você! Nele a gente explica melhor qual é o significado de cada termo da Equação de Navier-Stokes, como usa-la em exercícios e muito mais! 👇 🤩

Equação de Navier-Stokes: Coordenadas Cilíndricas

Quando escolhemos o sistema de coordenadas cilíndricas para o nosso problema, precisamos usar a equação de Navier-Stokes nessa coordenada. Então temos:

Como usar a Equação de Navier-Stokes

Pra a gente aprender como usar a equação de Navier-Stokes, vamos fazer um exercício juntos!

Considere um escoamento permanente, bidimensional, incompressível, de um fluido newtoniano no qual o campo de velocidade é . Encontre o escoamento satisfaz a conservação da massa e o campo de pressão se a pressão no ponto é igual a .

Você pode conferir a resolução desse exemplo no vídeo abaixo:

Agora vamos praticar com mais exercícios?!

Como usar a Equação de Navier-Stokes

Exercícios Resolvidos

Exercício Resolvido #1

Considere o campo de velocidade no plano x y dado por:

V   = x y 2 i - 1 3 y 3 j + x y k

Determine ( a )   o número de dimensões do escoamento, ( b ) se ele é um possível escoamento incompressível e ( c ) a aceleração de uma partícula fluida no ponto ( x ,   y ,   z )   =   ( 1 ,   2 ,   3 ) .

Passo 1

E aí, tudo bem? Tudo vai ficar tranquilo quando analisarmos esse perfil de velocidade! 😊

Separando a velocidade V em u , v e w (suas componentes nos eixos x, y e z respectivamente), temos isso aí:

V = u i + v j + w k

u = x y 2 v = - 1 3 y 3 w = x y

Repara uma coisa em u , v e w : nenhuma delas depende de z! As componentes da velocidade dependem apenas de x e y. Isso quer dizer que o escoamento está acontecendo justamente no plano xy. Ou seja, é um escoamento bidimensional :D

Então, encontramos a resposta da letra a. O número de dimensões do escoamento é 2.

Passo 2

Para saber se esse escoamento bidimensional é incompressível, precisamos verificar se a Equação da Continuidade é satisfeita.

Equação da Continuidade:

∂ ρ ∂ t + ∂ ∂ x ρ u + ∂ ∂ y ρ v + ∂ ∂ z ρ w = 0

Considerando escoamento permanente ∂ ρ ∂ t = 0 e incompressível ( ρ = c o n s t a n t e ) :

ρ ∂ u ∂ x + ρ ∂ v ∂ y + ρ ∂ w ∂ z = 0

∂ u ∂ x + ∂ v ∂ y + ∂ w ∂ z = 0

Lembrando que:

u = x y 2 v = - 1 3 y 3 w = x y

Aplicando o campo de velocidade na equação da continuidade:

y 2 - 3 3 y 2 + 0 = 0

y 2 - y 2 = 0

0 = 0

Logo, o campo de velocidade considerando fluido incompressível satisfaz a continuidade.

Ou seja, é um escoamento incompressível!

Passo 3

Para calcular a aceleração em um ponto, precisamos encontrar a expressão da aceleração e aplicá-la nesse ponto.

Como conseguir essa expressão? A derivada da velocidade é a aceleração, então vamos derivar nosso campo de velocidade

a x = ∂ u ∂ t + u ∂ u ∂ x + v ∂ u ∂ y + w ∂ u ∂ z a y = ∂ v ∂ t + u ∂ v ∂ x + v ∂ v ∂ y + w ∂ v ∂ z a z = ∂ w ∂ t + u ∂ w ∂ x + v ∂ w ∂ y + w ∂ w ∂ z

Lembrando que fizemos a hipótese de escoamento permanente no Passo 2, as derivadas em relação ao tempo serão nulas

a x = u ∂ u ∂ x + v ∂ u ∂ y + w ∂ u ∂ z a y = u ∂ v ∂ x + v ∂ v ∂ y + w ∂ v ∂ z a z = u ∂ w ∂ x + v ∂ w ∂ y + w ∂ w ∂ z

Além disso, nenhuma componente da velocidade depende de z, então as derivadas em relação a z também são zero

a x = u ∂ u ∂ x + v ∂ u ∂ y a y = u ∂ v ∂ x + v ∂ v ∂ y a z = u ∂ w ∂ x + v ∂ w ∂ y

Tá ficando menor, vamos que vamos! 😊

Agora vamos lembrar do campo de velocidade e finalmente derivar

u = x y 2 v = - 1 3 y 3 w = x y

a x = x y 2 y 2 + - 1 3 y 3 2 x y a y = x y 2 0 + - 1 3 y 3 - y 2 a z = x y 2 y + - 1 3 y 3 x

Organizando

a x = x y 4 + - 2 3 x y 4 a y = 1 3 y 5 a z = x y 3 + - 1 3 x y 3

O enunciado pede a aceleração no ponto (1, 2, 3), logo:

a x = 1 . 2 4 + - 2 3 2 4 a y = 1 3 2 5 a z = 1 . 2 3 + - 1 3 2 3 ⇒ a x = 16 - 32 3 a y = 32 3 a z = 8 - 8 3

a x = 16 3 a y = 32 3 a z = 16 3

Logo, a aceleração no ponto (1, 2, 3) é

a ( 1 ,   2 ,   3 )   = 16 3 i +   32 3 j + 16 3 k

Resposta

(a)2 dimensões

(b)Sim, é um escoamento incompressível

(c) a ( 1 ,   2 ,   3 )   = 16 3 i +   32 3 j + 16 3 k

Exercício Resolvido #2

O escoamento pelo bocal convergente da Figura pode ser aproximado por uma distribuição de velocidades unidimensional

u   ≈   V 0 1 +   2 x L ,   v   ≈ 0   ,   w   ≈ 0

a Encontre uma expressão geral para a aceleração do fluido no bocal.

( b )   Para o caso específico V 0 = 3   m / s e L = 150   m m , calcule a aceleração, em g’s, na entrada e na saída.

Passo 1

Esse problema exige de nós lembrar que a aceleração de um fluido é a primeira derivada da velocidade.

Como nosso escoamento é unidimensional, nossa expressão da aceleração se reduz a um termo.

Ao invés de analisarmos isso tudo:

d u d t = ∂ u ∂ t + u ∂ u ∂ x + v ∂ u ∂ y + w ∂ u ∂ z d v d t = ∂ v ∂ t + u ∂ v ∂ x + v ∂ v ∂ y + w ∂ v ∂ z   d w d t = ∂ w ∂ t + u ∂ w ∂ x + v ∂ w ∂ y + w ∂ w ∂ z

É necessário olhar apenas a primeira equação:

d u d t = ∂ u ∂ t + u ∂ u ∂ x + v ∂ u ∂ y + w ∂ u ∂ z

Passo 2

Devemos agora substituir a expressão:

u = V o 1 + 2 x L

Na equação:

d u d t = ∂ u ∂ t + u ∂ u ∂ x + v ∂ u ∂ y + w ∂ u ∂ z

Logo, resolvendo cada derivada:

∂ u ∂ t = 0

u ∂ u ∂ x = V o 1 + 2 x L V o 2 L = 2 V o ² L 1 + 2 x L

v ∂ u ∂ y = 0

w ∂ u ∂ z = 0

Somando tudo, temos a resposta da Letra A:

d u d t = 2 V o ² L 1 + 2 x L

Passo 3

Na letra B ele nos pede para substituir V o = 3   m / s e L = 150   m m = 0,150   m , na expressão que encontramos para a aceleração na entrada x = 0 e na saída x = L = 0,150   m :

d u d t = 2 V o ² L 1 + 2 x L

Na entrada ( x = 0 ):

d u d t = 2 × 3 2 0,150 1 + 2 × 0 0,150

d u d t = 2 × 3 2 0,150 1 + 0

d u d t = 18 0,150

d u d t = 120 m s 2 = 12,2 g

Na saída ( x = L )

d u d t = 2 V o ² L 1 + 2 L L

d u d t = 2 V o ² L 1 + 2

d u d t = 2 × 3 2 0,150 × 3

d u d t = 360 m s 2 = 36,7 g

Sabendo que:

1   g = 9,80665   m / s ²

Resposta

Letra a)

d u d t = 2 V o ² L 1 + 2 x L

Letra b)

d u d t ( x = 0 ) = 120 m / s ² = 12,2 g

d u d t ( x = L ) = 360 m / s ² = 36,7 g

Exercício Resolvido #3

Um campo de velocidade bidimensional é dado por

V = x 2 - y 2 + x i - 2 x y + y j

em unidades arbitrárias. Em x ,   y = ( 1 ,   2 ) , calcule ( a ) as acelerações a x e a y , ( b ) o componente de velocidade na direção θ = 40 ° , ( c ) a direção da velocidade máxima e ( d ) a direção da aceleração máxima.

Passo 1

Neste exercício, temos um campo bidimensional de velocidade, logo podemos reduzir o númeo de equações de aceleração para duas apenas:

d u d t = ∂ u ∂ t + u ∂ u ∂ x + v ∂ u ∂ y + w ∂ u ∂ z d v d t = ∂ v ∂ t + u ∂ v ∂ x + v ∂ v ∂ y + w ∂ v ∂ z  

Sabemos que a componente da velocidade no eixo x é igual a u e que no eixo y é igual a v , podemos reescrever o campo nos termos de u e v :

u = x 2 - y 2 + x

v = - ( 2 x y + y )

Passo 2

Na letra a, o problema nos pede para calcular a x e a y , que nada mais são:

a x = d u d t = ∂ u ∂ t + u ∂ u ∂ x + v ∂ u ∂ y + w ∂ u ∂ z

a y = d v d t = ∂ v ∂ t + u ∂ v ∂ x + v ∂ v ∂ y + w ∂ v ∂ z

Resolvendo então as derivadas, primeiro para a x :

∂ u ∂ t = 0 u ∂ u ∂ x = x 2 - y 2 + x 2 x + 1 v ∂ u ∂ y = - 2 x y + y - 2 y w ∂ u ∂ z = 0

Para a y :

∂ v ∂ t = 0 u ∂ v ∂ x = x 2 - y 2 + x ( - 2 y ) v ∂ v ∂ y = - 2 x y + y - 2 x - 1 w ∂ u ∂ z = 0

Em x , y = ( 1,2 ) :

a x = x 2 - y 2 + x 2 x + 1 - 2 x y + y - 2 y

a x = 1 2 - 2 2 + 1 2 × 1 + 1 - 2 × 1 × 2 + 2 - 2 × 2

a x = 1 - 4 + 1 2 + 1 - 4 + 2 - 4

a x = - 2 3 - 6 - 4

a x = - 6 + 24

a x = 18

E, mais rápido dessa vez:

a y = x 2 - y 2 + x ( - 2 y ) - 2 x y + y - 2 x - 1

a y = 1 2 - 2 2 + 1 ( - 2 × 2 ) - 2 × 1 × 2 + 2 - 2 × 1 - 1

a y = - 2 - 4 + - 6 - 3

a y = 26

Passo 3

Na letra b, o exercício nos pede a representação de V em uma direção θ = 40 o .

Isso é nada menos que a projeção de V em um vetor que representa essa direção.

Podemos definir um vetor n = ( cos ⁡ 40 o , sin ⁡ 40 o ) unitário.

Substituindo x , y = 1,2 em V :

u = x 2 - y 2 + x = - 2

v = - 2 x y + y = - 6

V = ( - 2 , - 6 )

Realizando o produto escalar

V ∙ n = ( - 2 , - 6 ) ∙ ( cos ⁡ 40 o , sin ⁡ 40 o )

V ∙ n = - 2 cos ⁡ 40 o - 6 sin ⁡ 40 o

V ∙ n ≅ - 5,39

Passo 4

Na letra c, temos que calcular em que direção temos velocidade máxima.

Já sabemos da letra anterior que V = ( - 2 , - 6 )

Percebe que se seguirmos o mesmo princípio da letra b, devemos achar um ângulo θ que maximize V .

Isso lembra máximo e mínimos de cálculo, e é exatamente o que faremos:

Para θ arbitrário:

V ∙ n = ( - 2 , - 6 ) ∙ ( cos ⁡ θ , sin ⁡ θ )

V ∙ n = - 2 cos ⁡ θ - 6 sin ⁡ θ

Vamos chamar essa expressão de f e derivá-la em θ :

f = V ∙ n = - 2 cos ⁡ θ - 6 sin ⁡ θ

d f d θ = 2 sin ⁡ θ - 6 cos ⁡ θ

Para encontrarmos um ponto ótimo devemos igualar a zero:

d f d θ = 0

2 sin ⁡ θ - 6 cos ⁡ θ = 0

2 sin ⁡ θ = 6 cos ⁡ θ

sin ⁡ θ cos ⁡ θ = 6 2

tan ⁡ θ = 3

θ = tan - 1 ⁡ 3

Temos que lembrar que em um cículo trigonométrico dois ângulos pelo menos tem a mesma tangente, e eles estão separados por 180 o .

No nosso problema um deles nos dará a máxima velocidade e outro a mínima, que são iguais em módulo:

f = - 2 cos ⁡ 71 o - 6 sin ⁡ 71 o = - 6,32

E,

f = - 2 cos ⁡ 251 o - 6 sin ⁡ 251 o = + 6,32

Logo V é máxima em módulo nestes dois ângulos.

Passo 5

Finalmente, repetiremos o processo do Passo 4, mas agora para a aceleração:

Em x , y = 1,2 calculamos a = ( 18,26 )

Sendo g = a ∙ n = 18 cos ⁡ θ + 26 sin ⁡ θ

d g d θ = 0

- 18 sin ⁡ θ + 26 cos ⁡ θ = 0

tan ⁡ θ = 26 18

θ = 55,3 o

Resposta

Letra a)

a x = 18 a y = 26

Letra b)

- 5,39

Letra c)

θ = 251 o

Letra d)

θ = 55,3 o

Exercício Resolvido #4

Um campo de escoamento permanente, incompressível, sem atrito, é dado por

V = 2 x y i - y 2 j

em unidades arbitrárias. Seja a massa específica ρ 0 = c o n s t a n t e e despreze a gravidade. Encontre uma expressão para o gradiente de pressão na direçãoo x .

Passo 1

Neste exercício ele nos dá algumas infomações bem importantes:

-Escoamento permanente: toda derivada temporal será igual a zero

-Escoamento incompressível: a massa específica é constante em todo escoamento

-Sem atrito: viscosidade desprezível

-Desprezaremos a gravidade

Agora só precisamos saber quem é nosso gradiente de pressão ∂ p ∂ x

Passo 2

Aplicando as condições do Passo 1 na primeira expressão de Navier-Stokes:

ρ g x - ∂ p ∂ x + μ ∂ 2 u ∂ x 2 + ∂ 2 u ∂ y 2 + ∂ 2 u ∂ z 2 = ρ d u d t

Ou seja,

g x = 0

μ = 0

O que nos resta é

- ∂ p ∂ x =   ρ d u d t

Sabemos que:

d u d t = ∂ u ∂ t + u ∂ u ∂ x + v ∂ u ∂ y + w ∂ u ∂ z

Sendo ∂ u ∂ t = 0

Juntando tudo temos a seguinte expressão para resolver:

- ∂ p ∂ x =   ρ u ∂ u ∂ x + v ∂ u ∂ y + w ∂ u ∂ z

Passo 3

Do enunciado:

u = 2 x y v = - y 2

Aplicando na nossa expressão

- ∂ p ∂ x =   ρ u ∂ u ∂ x + v ∂ u ∂ y + w ∂ u ∂ z

- ∂ p ∂ x =   ρ 2 x y 2 y + - y 2 2 x + 0

- ∂ p ∂ x =   ρ 4 x y ² - 2 x y ²

∂ p ∂ x = -   2 ρ x y ²

Resposta

∂ p ∂ x = -   2 ρ x y ²

Exercício Resolvido #5

Um campo de escoamento incompressível, tridimensional, proposto tem a seguinte forma vetorial:

V   = K x i + K y j - 2 K z k

( a )   Determine se esse campo é uma solução válida para a continuidade e para Navier-Stokes. ( b ) Se g   =   - g k , encontre o campo de pressão p ( x ,   y ,   z ) .

Passo 1

Neste exercício temos as seguintes informações do enunciado:

-O campo de velocidade V :

u = K x v = K y w = - 2 K z

-Escoamento incompressível:

ρ = c o n s t a n t e

-Gravidade apenas em z :

g z = - g

Passo 2

Letra a e b)

Aplicaremos o campo primeiro a Equação da Continuidade:

∂ ρ ∂ t + ∂ ∂ x ρ u + ∂ ∂ y ρ v + ∂ ∂ z ρ w = 0

∂ u ∂ x + ∂ v ∂ y + ∂ w ∂ z = 0

K + K - 2 K = 0

0 = 0

Portanto, esse campo satisfaz a continuidade

Aplicando a Navier-Stokes:

ρ g x - ∂ p ∂ x + μ ∂ 2 u ∂ x 2 + ∂ 2 u ∂ y 2 + ∂ 2 u ∂ z 2 = ρ ∂ u ∂ t + u ∂ u ∂ x + v ∂ u ∂ y + w ∂ u ∂ z ρ g y - ∂ p ∂ y + μ ∂ 2 v ∂ x 2 + ∂ 2 v ∂ y 2 + ∂ 2 v ∂ z 2 = ρ ∂ v ∂ t + u ∂ v ∂ x + v ∂ v ∂ y + w ∂ v ∂ z ρ g z - ∂ p ∂ z + μ ∂ 2 w ∂ x 2 + ∂ 2 w ∂ y 2 + ∂ 2 w ∂ z 2 = ρ ∂ w ∂ t + u ∂ w ∂ x + v ∂ w ∂ y + w ∂ w ∂ z

- ∂ p ∂ x = ρ u ∂ u ∂ x - ∂ p ∂ y = ρ v ∂ v ∂ y - ρ g - ∂ p ∂ z = ρ w ∂ w ∂ z

- ∂ p ∂ x = ρ K ² x - ∂ p ∂ y = ρ K ² y - ρ g - ∂ p ∂ z = ρ 4 K ² z

∂ p ∂ x = - ρ K ² x ∂ p ∂ y = - ρ K ² y ∂ p ∂ z = - 4 ρ K ² z - ρ g

Agora separando as varíáveis de cada uma das equações e integrando:

∫ ∂ p = ∫ - ρ K 2 x   ∂ x

p = - ρ K 2 x 2 2 + f ( y , z ) + p a

Na segunda equação:

∫ ∂ p = ∫ - ρ K 2 y   ∂ y

p = - ρ K 2 y 2 2 + f ( x , z ) + p a

Na terceira:

∫ ∂ p = ∫ - 4 ρ K ² z - ρ g ∂ z

p = - 2 ρ K ² z ² - ρ g z + f x , y + p a

Podemos dizer por comparação:

p x , y , z = - ρ K 2 x 2 2 - ρ K 2 y 2 2 - 2 ρ K 2 z 2 - ρ g z + p a

Ou,

p x , y , z = - ρ K 2 2 x 2 + y 2 + 4 z ² - ρ g z + p a

Resposta

Letra a)É válido para ambos

Letra b) p x , y , z = - ρ K 2 2 x 2 + y 2 + 4 z ² - ρ g z + p a

Exercício Resolvido #6

Considere um escoamento permanente, bidimensional, incompressível, de um fluido newtoniano no qual o campo de velocidade é conhecido: u =   - 2 x y , v = y 2 - x 2 , w = 0 . ( a ) Esse escoamento satisfaz a conservação da massa? ( b ) Encontre o campo de pressão   p ( x ,   y )   se a pressão no ponto   ( x = 0 ,   y = 0 )   é igual a p a .

Passo 1

Vamos lá,

O enunciado já nos fornece as seguintes informações:

-Escoamento permanente, ou seja, independe do tempo;

-Escoamento incompressível, ou seja, ρ = c o n s t a n t e ;

-Fluido newtoniano, ou seja, μ = c o n s t a n t e ;

-O campo de velocidade:

u = - 2 x y v = y ² - x ² w = 0

-Na letra A devemos aplicar este campo na Equação da Continuidade;

-Na letra B, aplicar nas equações de Navier-Stokes para encontrarmos o campo de pressão.

Passo 2

Letra a)

A Equação da Continuidade:

∂ ρ ∂ t + ∂ ∂ x ρ u + ∂ ∂ y ρ v + ∂ ∂ z ρ w = 0

Aplicando as informações do enunciado:

∂ u ∂ x + ∂ v ∂ y = 0

Substituindo o campo de velocidade:

- 2 y + 2 y = 0

0 = 0

Satisfaz portanto a Continuidade.

Passo 3

Letra b)

Utilizaremos as seguintes equações de Navier-Stokes:

ρ g x - ∂ p ∂ x + μ ∂ 2 u ∂ x 2 + ∂ 2 u ∂ y 2 + ∂ 2 u ∂ z 2 = ρ ∂ u ∂ t + u ∂ u ∂ x + v ∂ u ∂ y + w ∂ u ∂ z ρ g y - ∂ p ∂ y + μ ∂ 2 v ∂ x 2 + ∂ 2 v ∂ y 2 + ∂ 2 v ∂ z 2 = ρ ∂ v ∂ t + u ∂ v ∂ x + v ∂ v ∂ y + w ∂ v ∂ z

Simplificando com as informações do enunciado e que a velocidade não varia em z :

ρ g x - ∂ p ∂ x + μ ∂ 2 u ∂ x 2 + ∂ 2 u ∂ y 2 = ρ u ∂ u ∂ x + v ∂ u ∂ y ρ g y - ∂ p ∂ y + μ ∂ 2 v ∂ x 2 + ∂ 2 v ∂ y 2 = ρ u ∂ v ∂ x + v ∂ v ∂ y

Podemos fazer uma suposição de que só exista componente da gravidade no eixo z , ou seja que nosso referencial esteja com o eixo z paralelo a direção da força de gravidade

- ∂ p ∂ x + μ ∂ 2 u ∂ x 2 + ∂ 2 u ∂ y 2 = ρ u ∂ u ∂ x + v ∂ u ∂ y - ∂ p ∂ y + μ ∂ 2 v ∂ x 2 + ∂ 2 v ∂ y 2 = ρ u ∂ v ∂ x + v ∂ v ∂ y

Substituindo u e v nas expressões

- ∂ p ∂ x + μ 0 + 0 = ρ 4 x y ² + 2 x ³ - 2 x y ² - ∂ p ∂ y + μ - 2 + 2 = ρ 4 x ² y + 2 y ³ - 2 x ² y

- ∂ p ∂ x = ρ 2 x y ² + 2 x ³ - ∂ p ∂ y = ρ 2 x ² y + 2 y ³

Devemos então realizer as integrações parciais:

∫ ∂ p = ∫ - ρ 2 x y ² + 2 x ³ ∂ x ∫ ∂ p = ∫ - ρ 2 x ² y + 2 y ³ ∂ y

Resolvendo a primeira:

∫ ∂ p = ∫ - ρ 2 x y ² + 2 x ³ ∂ x

∫ ∂ p = - 2 ρ ∫ x y ² + x ³ ∂ x

p = - 2 ρ x ² y ² 2 + x 4 4 + f y

O mesmo para a segunda:

∫ ∂ p = ∫ - ρ 2 x 2 y + 2 y 3 ∂ y

∫ ∂ p = - 2 ρ ∫ x ² y + y ³ ∂ y

p = - 2 ρ x ² y ² 2 + y 4 4 + f x

Comparando as duas expresses podemos dizer que:

f x = x 4 4 + C

Ou que,

f y = y 4 4 + C

Substituindo f ( y ) ou f ( x ) em uma das expressões, podemos encontrara pressão p :

p = - 2 ρ x ² y ² 2 + x 4 4 + y 4 4 + C

Passo 4

Aplicando a condição de Contorno fornecida p 0,0 = p a :

p ( x , y ) = - 2 ρ x ² y ² 2 + x 4 4 + y 4 4 + C

p 0,0 = p a = - 2 ρ 0 ² 0 ² 2 + 0 4 4 + 0 4 4 + C

p a = - 2 ρ C

C = - p a 2 ρ

Logo a expressão da pressão é a que segue:

p ( x , y ) = - 2 ρ x ² y ² 2 + x 4 4 + y 4 4 - p a 2 ρ

Arrumando um pouco:

p x , y = p a - ρ x ² y ² + x 4 + y 4 2

Resposta

Letra a)

Sim.

Letra b)

p x , y = p a - ρ x ² y ² + x 4 + y 4 2

Exercício Resolvido #7

Considere a esfera de raio R imersa em uma corrente uniforme U 0 , como mostra a Figura. De acordo com a teoria, a velocidade do fluido ao longo da linha de corrente A B é dada por

V = u i =   U 0 1 +   R ³ x ³ i

Encontre ( a ) a posição da aceleração máxima do fluidod ao longo de A B e ( b ) o tempo necessário para que uma partícula de fluido se desloque de A a B .

Passo 1

O enunciado nos fornece o campo de velocidade

V = u i = U o 1 + R ³ x ³ i

e nos pede a aceleração máxima ao longo de x , antes então devemos encontrar a função que define essa aceleração.

Passo 2

Sabemos que:

a = d V d t = i d u d t + j d v d t + k d w d t

Como é um campo unidimensional,

a = d V d t = i d u d t

Se,

d u d t = ∂ u ∂ t + u ∂ u ∂ x + v ∂ u ∂ y + w ∂ u ∂ z

Temos que:

a = ∂ u ∂ t + u ∂ u ∂ x + v ∂ u ∂ y + w ∂ u ∂ z

Substituindo u :

a = 0 + U o 1 + R ³ x ³ ∂ ∂ x U o 1 + R ³ x ³ + 0 + 0

a = U o ² 1 + R ³ x ³ - 3 R ³ x 4

a = - 3 U o ² R ³ 1 x 4 + R ³ x 7

Devemos então encontrar o ponto que a é máximo, ou seja,

∂ a ∂ x = 0

- 3 U o ² R ³ - 4 1 x 5 - 7 R ³ x 8 = 0

- 4 1 x 5 - 7 R ³ x 8 = 0

4 1 x 5 = - 7 R 3 x 8

x 8 x 5 = - 7 R 3 4

x ³ = - 7 R 3 4

x = R - 7 4 3

x ≅ - 1,205 R

Passo 3

Na letra b) queremos calcular quanto tempo leva para o fluido se movimentar de A até B .

Sabemos que

u = d x d t

U o 1 + R ³ x ³ = d x d t

Se separarmos a varíáveis ficamos com,

U o d t = d x 1 + R ³ x ³

Integrando nos corretos limites de integração:

∫ 0 t U o d t = ∫ - 4 R - R d x 1 + R ³ x ³

Tentar resolver esta integral é muito complicado, pois ela não converge para um valor e tende ao infinito. O objetivo da questão era demonstrar que o fluido leva um tempo infinito para chegar ao ponto de estagnação em B ( x = - R ) , onde a velocidade do fluido é zero.

Logo,

t → ∞

Resposta

Letra a)

x ≅ - 1,205 R

Letra b)

t → ∞

Exercício Resolvido #8

Das equações de Navier-Stokes para escoamento incompressível em coordenadas polares encontre o caso mais geral de movimento puramente circulatório v θ ( r ) , v r =   v z = 0 , para escoamento sem escorregamento entre dois cilindros concêntricos fixos, como na Figura.

Equação da continuidade em coordenadas cilíndricas:

1 r ∂ r v r ∂ r + 1 r ∂ v θ ∂ θ + ∂ v z ∂ z = 0

Equações de Navier-Stokes em coordenadas cilíndricas:

∂ v r ∂ t + v r ∂ v r ∂ r + v θ r ∂ v r ∂ θ + v z ∂ v r ∂ z - v θ 2 r = - 1 ρ ∂ p ∂ r + g r + μ ρ 1 r ∂ ∂ r r ∂ v r ∂ r + 1 r ² ∂ ² v r ∂ θ ² + ∂ ² v r ∂ z ² - v r r 2 - 2 r ² ∂ v θ ∂ θ ∂ v θ ∂ t + v r ∂ v θ ∂ r + v θ r ∂ v θ ∂ θ + v z ∂ v θ ∂ z + v r v θ r = - 1 ρ r ∂ p ∂ θ + g θ + μ ρ 1 r ∂ ∂ r r ∂ v θ ∂ r + 1 r ² ∂ ² v θ ∂ θ ² + ∂ ² v θ ∂ z ² - v θ r 2 - 2 r ² ∂ v r ∂ θ ∂ v z ∂ t + v r ∂ v z ∂ r + v θ r ∂ v z ∂ θ + v z ∂ v z ∂ z = - 1 ρ ∂ p ∂ z + g z + μ ρ 1 r ∂ ∂ r r ∂ v z ∂ r + 1 r ² ∂ ² v z ∂ θ ² + ∂ ² v z ∂ z ²

Passo 1

Primeiro de tudo já sabemos que v θ só está em função de r , e que v r = v z = 0 . Já sabemos daí que a continuidade será satisfeita.

Esses dados simplificarão em muito o problema também. Queremos encontrar o valor da função v θ ( r ) . Vamos supor que os efeitos gravitacionais não existem.

Passo 2

Da equação de θ de Navier-Stokes, já aplicando as condições do Passo 1:

0 = - 1 ρ r ∂ p ∂ θ + μ ρ 1 r ∂ ∂ r r ∂ v θ ∂ r - v θ r 2

Por simetria podemos dizer que p não varia com θ :

0 = μ ρ 1 r ∂ ∂ r r ∂ v θ ∂ r - v θ r 2

1 r ∂ ∂ r r ∂ v θ ∂ r - v θ r 2 = 0

Se resolvermos a derivada pela Regra do Produto:

∂ ² v θ ∂ r ² + 1 r ∂ v θ ∂ r - 1 r 2 v θ = 0

A solução desta equação é:

v θ = C 1 r + C 2 r

Passo 3

Se substituirmos as condições de contorno v θ a = 0 devido ao não escorregamento:

v θ a = C 1 a + C 2 a = 0

C 1 = - C 2 a ²

Para v θ b = 0 também por não-escorregamento:

v θ b = C 1 b + C 2 b = 0

- C 2 a ² b + C 2 b = 0

C 2 - b a 2 + 1 b = 0

Supondo que C 2 ≠ 0

- b a 2 + 1 b = 0

1 b = b a 2

a = b

Se isso for verdade, não haveria escoamento entre os cilindros para satisfazer as condições inciais.

Como o escoamento pode existir, devemos reavliar as hipóteses iniciais. Podemos dizer que o escoamento não existe em estado permanente, ou seja, suas condições variam no tempo o que torna o escoamento muito mais complicado de analisar.

Mas isso não excluiu o fato de existir uma solução geral de escoamentos entre cilindros concêntricos.

Resposta

v θ = C 1 r + C 2 r

Exercícios de Livros Relacionados

O campo de velocidade para um sorvedouro plano do tipo vórtice é dado por V → = ( - q / 2 π r ) ê r+( K / 2 π r ) ê θ ê ê , em que q = 2 m³/s/m e K = 1 m³/s/m. A massa específica do fluido é 1000 kg/m
Um líquido incompressível com viscosidade desprezível e massa específica ρ = 1250 kg/m3 escoa em regime permanente através de um tubo de 5 m de comprimento com seção convergentedivergente, cuja área v
Para modelar a distribuição de velocidade na seção curva de entrada de um túnel de vento, o raio de curvatura das linhas de corrente é expresso como R = LR0/2y. Em primeira aproximação, considere que
Repita o Exemplo 6.1, porém com uma consideração um pouco mais realista de que o escoamento é similar a um perfil de vórtice livre (irrotacional), Vθ = c/r (em que c é uma constante), como mostrado na
Uma dimensão de uma placa retangular é o dobro da outra dimensão. O ar a uma velocidade uniforme flui paralelo à placa e uma camada limite laminar se forma em ambos os lados da placa. Qual a orientaçã