Equações de Navier-Stokes

Teoria Completa

Comportamento de forças em fluidos newtonianos.

Faaala!!!

Hoje no RespondeAí Repórter vamos ver o comportamento do conjunto de equações mais importante da Mecânica dos Fluidos: as EQUAÇÕES DE NAVIER-STOKES.

Onde vivem? Do que se alimentam? Pra que servem?

Essas e todas perguntas serão resolvidas hoje!

Não, elas não se alimentam dos sonhos de formatura dos alunos e não, elas não servem para nada.

Antes de falarmos por trás delas, vamos conhecer a figurinha, se liga:

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

Calma jovem aprendiz vamos olhar o que elas estão nos dizendo.

Para isso, precisamos responder uma pergunta:

Qual grandeza estamos analisando em Navier-Stokes?

Tá vendo aquele termo ρ d u d t no lado direito da primeira equação?

Lembra lá da Segunda Lei de Newton, aquela tal de F = m a , força igual a massa vezes aceleração?

Lembra que aceleração é a primeira derivada da velocidade no tempo? E que massa é igual ao volume vezes a massa específica ρ ? Em outras palavras:

a = d V d t

m = ρ ∇

Substituindo na nossa expressão:

F = m a

F = ρ ∇   d V d t

Ou, olhando para uma componente só,

F x = ρ ∇ d u d t

E se a gente dividir os dois lados por ? Só pra gente pensar em força por unidade de volume:

F ∇ = ρ ∇ d u d t ∇

F ∇ =   ρ d u d t

Ta vendo aquele termo no lado direito? Não é o mesmo do lado direito da primeira equação de Navier-Stokes?

Acabamos de descobrir o que está acontecendo naquelas três equações: estamos analisando força sobre volume.

Mas tem uma pegadinha, essa aceleração d u d t é uma das componentes do campo de aceleração do fluido. E em MecFlu nós não seguimos o movimento de uma só partícula igual em movimento de sólidos, aqui nós vamos observar o fluido ponto a ponto, isso significa que precisamos trabalhar com campos de velocidade e aceleração.

Nosso campo de velocidade vai ser representado da seguinte forma:

V x , y , z , t = u x , y , z , t i + v x , y , z , t j + w x , y , z , t k

Aceleração em Navier-Stokes

A aceleração então é:

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

A pegadinha entra agora: como u ,   v   e w não dependem só do tempo, devemos utilizar a regra da cadeia, derivando a variação de x , y , z no tempo t também:

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

O mesmo vale para as outras componentes de velocidade.

Por definição:

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

O que substituindo nas nossas expressões da aceleração dá o seguinte:

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

Essas relações são muito importante pra quando a gente começar a fazer exercícios.

Substituindo nas 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 ρ 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

A essa altura, você pode estar pensando...

Então, eu não queria ser chata, mas a resposta é sim. Essas são as equações de Navier-Stokes. Se fosse um hambúrguer, ia ser o X-TUDÃO de MecFlu. Maaaaas, a boa notícia é que em cada exercício você vai conseguir simplificar bastante. :D

Vamos resolver juntos os problemas, mas a equação não costuma ser esse problemão aí todo não. Bora deixar isso mais simples! Para isso, vamos entender os termos.

Entendendo os termos de Navier-Stokes

Agora vamos dar uma olhada no lado esquerdo das equações de Navier-Stokes, usando a primeira como exemplo:

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

Podemos separar em três pedaços

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

O primeiro, lembra um pouco uma força conhecida, aquela tal de força peso m g → ...

Então é exatamente isso! É a força que a gravidade aplica nos nossos escoamentos, só que agora ela está no formato para fluidos: força sobre volume, N m ³ se usarmos o sistema internacional de medidas.

Próximo pedaço:

O que é aquele p ? Não é nada mais que a pressão no fluido, e como estamos no eixo x , queremos saber como ela varia nesse mesmo eixo.

Se olharmos as unidades, veremos que representa a mesma coisa: força sobre volume

Unidade de pressão: P a = N m ²

Unidade de comprimento: m

N / m ² m = N / m ³

O sinal negativo apareceu aí devido a convenção de sinais usados na dedução da expressão.

Terceiro e último pedaço: componentes cisalhantes

Quando arrastamos um bloco no chão sentimos uma resistência entre a superfície de contato do bloco e do chão. Essa força é a força de atrito que tanto conhecemos.

Em fluidos acontece a mesma coisa, só que em todo ele.

Imagine que o fluido são pequenos pedaços infinitesimais formados por algumas moléculas que estão sempre encostados no outro, se arrastando, um passando por cima do outro.

Ou seja, sempre há uma resistência própria de cada fluido a esses movimentos. Por exemplo o ar é muito menos resistente que a água não é?

A gente tem um outro nome pra isso: viscosidade ou atrito interno. Representado pela letra μ .

Todo fluido é viscoso, uns mais, outros menos. A gente ta mais acostumado em chamar o mel de algo viscoso, né mesmo? Mas todos eles apresentam essa propriedade inerente de resistência.

Quando um fluido tem viscosidade constante no tempo e para qualquer estado de tensões , chamamos esse fluido de fluido newtoniano.

As equações de Navier-Stokes só são válidas para escoamentos incompressíveis de fluidos newtonianos.

Agora, bora lá praticar! \o/

Expandir

Exercício Resolvido #1

Fox, McDonald, Pritchard. Introdução à Mecânica dos Fluidos, 8ª ed. Gen LTD Editora, 2014, pp. 293, ex. 5.38

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 ) .

MOSTRAR SOLUÇÃO COMPLETA

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