Soluções Exatas para as Equações de Navier-Stokes

As equações de Navier-Stokes são notoriamente não lineares. Então, como alguém pode resolvê-las exatamente?

Simetria. Esse é todo o truque.

Quando a geometria de um escoamento é suficientemente simples (um tubo reto, duas placas planas, um plano infinito), a velocidade só pode apontar em uma direção e variar ao longo de uma ou duas coordenadas. Em muitas dessas configurações simétricas, o termo não linear ou desaparece ou se simplifica tanto que as equações colapsam para uma EDP linear. Em casos estacionários, você frequentemente fica com uma EDO que pode resolver com lápis e papel.

Aqui está a intuição. As equações de Navier-Stokes descrevem todos os movimentos possíveis de fluidos. Mas se você forçar o fluido a uma situação muito ordenada, sem turbilhonamento, sem caos, tudo marchando em uma direção, a maior parte da complexidade da equação se torna irrelevante. A parte difícil de Navier-Stokes é o ciclo de retroalimentação onde o fluido empurra a si mesmo. Nesses escoamentos simétricos, não há nada contra o que empurrar. O termo não linear de auto-advecção simplesmente desaparece.

Essas soluções exatas não são curiosidades. Elas são a base da educação em mecânica dos fluidos, os benchmarks para códigos numéricos, e o ponto de partida para entender quando e como escoamentos reais se desestabilizam.

Uma solução exata para as equações de Navier-Stokes é um campo de velocidade $u$ e pressão $p$ satisfazendo

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0$$

em forma fechada, sem aproximação numérica.

O mecanismo chave é o desaparecimento ou simplificação do termo de advecção não linear $(u \cdot \nabla)u$. Para escoamentos paralelos (escoamentos onde $u = (u(y,t),\, 0,\, 0)$ em coordenadas cartesianas ou $u = (0,\, 0,\, u(r,t))$ em cilíndricas) o campo de velocidade é livre de divergência automaticamente e o termo de advecção desaparece identicamente, porque a velocidade não tem componente na direção do escoamento e nem variação ao longo dela. A equação de momento então se reduz a uma EDP linear, ou para escoamentos estacionários, uma EDO.

Quando o escoamento também é estacionário ($\partial_t u = 0$), o que resta é

$$\nu \Delta u = \nabla p,$$

uma equação de Poisson cujas soluções são as soluções exatas clássicas de escoamento viscoso: escoamento de Poiseuille, escoamento de Couette, e seus parentes.

Para escoamentos paralelos não estacionários ($\partial_t u \neq 0$ mas a advecção ainda desaparece), a equação se torna uma equação de difusão $\partial_t u = \nu\, \partial_{yy} u$, que produz os problemas primeiro e segundo de Stokes.

O escoamento de Poiseuille (também chamado escoamento de Hagen-Poiseuille) é a solução exata mais importante para as equações de Navier-Stokes, e é a que a maioria dos engenheiros encontra primeiro.

Imagine água fluindo constantemente através de um tubo longo e reto. Uma diferença de pressão constante entre as duas extremidades impulsiona o escoamento. As paredes do tubo não se movem, então o fluido tocando a parede está preso a velocidade zero (essa é a condição de não deslizamento). Mais longe das paredes, o fluido acelera. O centro se move mais rápido.

O perfil de velocidade? Uma parábola. Zero na parede. Máximo no centro. Curva suave entre eles. Corte o tubo e olhe a seção transversal: é uma tigela invertida.

O fato quantitativo chave é este: A vazão total escala com a quarta potência do raio do tubo. Quarta potência. Dobre o raio e você não obtém o dobro do escoamento, ou mesmo o quádruplo. Você obtém dezesseis vezes o escoamento. Essa é a lei de Hagen-Poiseuille, e ela explica por que até mesmo um pequeno estreitamento arterial pode sufocar o suprimento de sangue.

Suposições: O escoamento de Poiseuille assume que o fluido é incompressível e newtoniano (viscosidade constante), que o escoamento é estacionário e totalmente desenvolvido (não ainda acelerando desde a entrada), e que o escoamento é laminar. Suave. Ordenado. Na prática, o escoamento em tubo transita para turbulência a um número de Reynolds de aproximadamente 2.300, que é uma observação empírica que ninguém conseguiu derivar da teoria subjacente.

Considere escoamento estacionário, totalmente desenvolvido, axialmente simétrico em um tubo circular de raio $R$, impulsionado por um gradiente de pressão uniforme $dp/dx < 0$ na direção axial $x$. Em coordenadas cilíndricas $(r, \theta, x)$ o campo de velocidade tem a forma $u = (0,\, 0,\, u(r))$.

A condição de incompressibilidade $\nabla \cdot u = 0$ é satisfeita automaticamente. O termo de advecção desaparece porque $u$ não depende de $x$. A equação de momento axial se reduz a

$$0 = -\frac{dp}{dx} + \mu\left(\frac{d^2 u}{dr^2} + \frac{1}{r}\frac{du}{dr}\right),$$

ou equivalentemente

$$\frac{1}{r}\frac{d}{dr}\left(r\frac{du}{dr}\right) = \frac{1}{\mu}\frac{dp}{dx}.$$

Como $dp/dx$ é constante, esta é uma EDO em $r$. Integrando duas vezes com condições de contorno $u(R) = 0$ (não deslizamento) e $du/dr|_{r=0} = 0$ (simetria) dá o perfil de velocidade parabólico:

$$u(r) = \frac{1}{4\mu}\left(-\frac{dp}{dx}\right)(R^2 - r^2).$$

A velocidade máxima na linha central é $u_{\max} = \frac{R^2}{4\mu}\left(-\frac{dp}{dx}\right)$, e a velocidade média é $\bar{u} = u_{\max}/2$.

Integrando sobre a seção transversal produz a lei de Hagen-Poiseuille para a vazão volumétrica:

$$Q = \frac{\pi R^4 \Delta p}{8 \mu L},$$

onde $\Delta p > 0$ é a queda de pressão ao longo do comprimento do tubo $L$.

A dependência de $R^4$ é a característica definidora: pequenas mudanças no raio produzem grandes mudanças na vazão. Em hemodinâmica, é por isso que o fluxo sanguíneo é tão sensível à estenose arterial.

Validade: Esta solução se aplica a escoamento incompressível, newtoniano, estacionário, totalmente desenvolvido, laminar. A transição para turbulência ocorre experimentalmente em $\text{Re} = \bar{u} D / \nu \approx 2{,}300$, onde $D = 2R$ é o diâmetro do tubo. Este limiar é uma observação empírica; não há teorema prevendo-o a partir das equações de Navier-Stokes. Para uma discussão sobre número de Reynolds e a transição laminar-turbulenta, veja Número de Reynolds, Turbulência, e Por Que Pequenas Escalas Importam.

O escoamento de Couette é a solução exata para fluido preso entre duas placas paralelas quando uma placa se move e a outra permanece parada. É uma das soluções exatas mais simples em mecânica dos fluidos.

Imagine um baralho de cartas deitado sobre uma mesa. Arraste a carta superior lateralmente e as cartas abaixo também se deslocam, cada uma um pouco menos que a de cima. É isso. A velocidade varia linearmente de zero na placa inferior até a velocidade da placa superior, e não há nada mais do que isso.

Um perfil de velocidade em linha reta. Placa inferior: estacionária. Placa superior: em movimento. Tudo entre elas apenas interpola linearmente, sem gradiente de pressão necessário, sem configuração complicada, o movimento impulsionado puramente pela fronteira móvel arrastando o fluido.

As coisas ficam mais interessantes quando você também aplica um gradiente de pressão ao longo do canal, porque então você está combinando escoamento impulsionado por cisalhamento e impulsionado por pressão no mesmo espaço. O perfil de velocidade se deforma em uma parábola sobreposta ao perfil linear, às vezes chamado escoamento plano de Poiseuille-Couette, e dependendo da força do gradiente de pressão em relação à velocidade da placa, você pode até obter refluxo perto de uma parede.

Remova completamente a placa móvel, mantenha ambas as paredes estacionárias, e deixe uma diferença de pressão fazer todo o trabalho. Esse é o escoamento plano de Poiseuille, o análogo de placa plana do escoamento em tubo. Parabólico. Mais rápido no meio. Zero em ambas as paredes.

Considere escoamento estacionário entre duas placas paralelas infinitas separadas por um espaço $h$. Seja $y$ a coordenada perpendicular às placas, com a placa inferior em $y = 0$ e a placa superior em $y = h$. O escoamento é paralelo: $u = (u(y),\, 0,\, 0)$.

Escoamento de Couette simples. A placa superior se move a velocidade $U$ na direção $x$; a placa inferior é estacionária; $dp/dx = 0$. A equação de momento se reduz a $d^2u/dy^2 = 0$, produzindo

$$u(y) = U\frac{y}{h}.$$

Esta é a solução exata não trivial mais simples para as equações de Navier-Stokes: um perfil de velocidade linear impulsionado inteiramente pela condição de contorno.

Escoamento plano de Poiseuille. Ambas as placas são estacionárias; um gradiente de pressão constante $dp/dx < 0$ impulsiona o escoamento. A equação de momento se torna

$$\mu \frac{d^2 u}{dy^2} = \frac{dp}{dx},$$

com $u(0) = u(h) = 0$. A solução é

$$u(y) = \frac{1}{2\mu}\left(-\frac{dp}{dx}\right)y(h - y),$$

um perfil parabólico simétrico em torno da linha central $y = h/2$.

Escoamento geral de Couette-Poiseuille. Combinando uma placa superior móvel com um gradiente de pressão dá a superposição

$$u(y) = U\frac{y}{h} + \frac{1}{2\mu}\left(-\frac{dp}{dx}\right)y(h - y).$$

O primeiro termo é a contribuição impulsionada por cisalhamento; o segundo é a contribuição impulsionada por pressão. Dependendo do sinal e magnitude de $dp/dx$ em relação a $\mu U / h^2$, o perfil pode ser monótono, ter um máximo interior, ou até exibir refluxo perto de uma parede.

Os escoamentos de Poiseuille e Couette são estacionários. Nada muda no tempo. Os problemas de Stokes são as soluções exatas não estacionárias mais simples, e eles revelam algo belo: a viscosidade faz o momento difundir através de um fluido, assim como o calor se difunde através de um sólido.

Primeiro problema de Stokes (também chamado problema de Rayleigh): imagine uma vasta piscina parada de fluido repousando acima de uma placa plana. No tempo zero, a placa subitamente começa a deslizar lateralmente a velocidade constante. O fluido logo ao lado da placa é arrastado imediatamente, mas o fluido mais distante leva tempo para perceber. Uma camada limite suave cresce para fora da placa, ficando mais espessa conforme o tempo passa.

A velocidade a qualquer altura acima da placa depende da razão dessa altura para um comprimento característico de difusão $\sqrt{\nu t}$, onde $\nu$ é a viscosidade e $t$ é o tempo decorrido. Fluido mais viscoso? O movimento se espalha para cima mais rápido.

Segundo problema de Stokes: mesma configuração, mas agora a placa oscila para frente e para trás senoidalmente em vez de se mover a velocidade constante. A oscilação só penetra uma distância finita no fluido. Mais acima, o fluido mal percebe. A amplitude do movimento decai exponencialmente com a altura, criando uma fina camada limite oscilatória. Este é o mecanismo por trás das camadas limite oscilatórias: uma placa oscilante coloca o fluido próximo em movimento, mas a perturbação morre exponencialmente com a distância da placa.

Ambos os problemas de Stokes envolvem um fluido semi-infinito ($y > 0$) acima de uma placa plana em $y = 0$. O escoamento é paralelo: $u = (u(y,t),\, 0,\, 0)$. O termo de advecção desaparece, e a equação governante é a equação de difusão unidimensional:

$$\frac{\partial u}{\partial t} = \nu \frac{\partial^2 u}{\partial y^2}.$$

Primeiro problema de Stokes (problema de Rayleigh). Condição inicial: $u(y,0) = 0$ para $y > 0$. Condição de contorno: $u(0,t) = U$ para $t > 0$. Campo distante: $u \to 0$ quando $y \to \infty$.

A variável de similaridade $\eta = y / (2\sqrt{\nu t})$ reduz a EDP a uma EDO. A solução é

$$u(y,t) = U\,\operatorname{erfc}\!\left(\frac{y}{2\sqrt{\nu t}}\right),$$

onde $\operatorname{erfc}$ é a função erro complementar. A espessura da camada limite cresce como $\delta \sim \sqrt{\nu t}$, a marca registrada da difusão.

Segundo problema de Stokes. A placa oscila: $u(0,t) = U\cos(\omega t)$. Buscando uma solução da forma $u(y,t) = \operatorname{Re}[\hat{u}(y)\, e^{i\omega t}]$ dá

$$\hat{u}(y) = U\exp\!\left(-(1+i)\frac{y}{\delta_s}\right), \qquad \delta_s = \sqrt{\frac{2\nu}{\omega}},$$

de modo que

$$u(y,t) = U\exp\!\left(-\frac{y}{\delta_s}\right)\cos\!\left(\omega t - \frac{y}{\delta_s}\right).$$

A amplitude decai exponencialmente com a profundidade de penetração $\delta_s$, e a fase atrasa linearmente: uma onda transversal cuja energia é inteiramente dissipada pela viscosidade. Frequências mais altas penetram menos profundamente.

Os escoamentos de Poiseuille, Couette e Stokes recebem a maior parte da atenção. Eles não são as únicas soluções exatas, no entanto. Nem de longe.

• Vórtice de Taylor-Green: um padrão decrescente de vórtices girando em duas dimensões com estrutura vortical genuína. Testou um código CFD? Você provavelmente o executou contra isso. É o benchmark que todos alcançam primeiro, e tem sido assim por décadas.
• Escoamento de Jeffery-Hamel: escoamento em um canal em forma de cunha que converge ou diverge. Ele captura como o fluido acelera em um espaço estreitante ou desacelera em um expansor.
• Escoamento de ponto de estagnação de Hiemenz: fluido batendo de frente em uma parede plana, desacelerando a zero na superfície e desviando lateralmente. Vento atingindo um prédio. Um jato atingindo uma placa.

Simetria. Cada uma dessas explora uma simetria geométrica específica para tornar as equações tratáveis, e elas importam em contextos especializados, mas para mecânica dos fluidos introdutória cotidiana Poiseuille e Couette ainda fazem todo o trabalho pesado.

Além dos escoamentos paralelos, várias outras famílias de soluções exatas exploram simetrias específicas ou estruturas auto-similares:

• Vórtice de Taylor-Green. Em 2D, $u = (A\cos(ax)\sin(by),\, -B\sin(ax)\cos(by))$ com $aA = bB$ (incompressibilidade) e decaimento temporal exponencial $e^{-\nu(a^2+b^2)t}$. Esta é uma solução exata para as equações completas de Navier-Stokes 2D, incluindo o termo não linear (que acaba sendo um gradiente e é absorvido na pressão). É um caso de validação padrão para códigos DNS.
• Escoamento de Jeffery-Hamel. Escoamento puramente radial estacionário $u_r(r,\theta)$ em uma cunha de meio-ângulo $\alpha$. A função de corrente $\psi = \nu f(\theta)$ satisfaz uma EDO não linear de terceira ordem em $\theta$. Soluções existem tanto para canais convergentes quanto divergentes, com rica estrutura de bifurcação em números de Reynolds mais altos.
• Escoamento de ponto de estagnação de Hiemenz. Uma solução de similaridade 2D para escoamento incidindo em uma parede plana. A função de corrente $\psi = \sqrt{a\nu}\, x\, f(\eta)$, $\eta = y\sqrt{a/\nu}$, reduz Navier-Stokes à EDO de Hiemenz: $f''' + ff'' - f'^2 + 1 = 0$ com $f(0)=f'(0)=0$, $f'(\infty)=1$.

Podemos resolver as equações de Navier-Stokes exatamente em todos esses casos. Então por que ainda há um problema em aberto de um milhão de dólares?

Truques. Cada uma delas depende de um truque. A geometria é escolhida tão cuidadosamente que a parte mais difícil da equação (o termo não linear) ou desaparece inteiramente ou se reduz a algo gerenciável, e o problema se torna solucionável precisamente porque foi drenado de tudo que torna Navier-Stokes difícil. Escoamento em tubo? Unidimensional. Couette? Uma linha. O vórtice de Taylor-Green esconde sua não linearidade dentro da pressão.

O problema do Prêmio do Milênio pergunta sobre escoamentos tridimensionais gerais. Sem truques de simetria. Sem geometria simplificadora. Dados iniciais suaves livres de divergência, interação não linear completa através de todas as escalas. Nessa configuração, ninguém provou que as soluções sempre permanecem suaves, e ninguém provou que elas explodem. Nós genuinamente não sabemos.

Então essas soluções exatas nos dizem algo, mas nem de perto o suficiente. Elas provam que as equações podem produzir soluções suaves explícitas quando você lhes dá forte simetria para se apoiar. A questão de um milhão de dólares é se a suavidade se mantém sempre, para dados iniciais suaves livres de divergência arbitrários nas formulações 3D padrão, ou se em algum lugar na violência completa da turbulência algo dá catastrófica e irreversivelmente errado.

Toda solução exata discutida acima alcança tratabilidade eliminando ou trivializando o termo de advecção não linear $(u \cdot \nabla)u$. Em escoamentos paralelos ele desaparece identicamente. No vórtice de Taylor-Green é um gradiente absorvido na pressão. Em soluções de similaridade ele se reduz via mudança de variáveis a um problema de dimensão inferior.

O problema do Prêmio do Milênio Clay diz respeito ao problema de valor inicial para as equações de Navier-Stokes incompressíveis 3D com dados iniciais suaves e rapidamente decrescentes arbitrários, precisamente o regime onde nenhuma dessas simplificações se aplica.

A questão é se as soluções para

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \quad \nabla \cdot u = 0, \quad u(x,0)=u_0(x)$$

em $\mathbb{R}^3$ permanecem em $C^\infty(\mathbb{R}^3 \times [0,\infty))$ para todo tempo, dados dados iniciais suaves, livres de divergência, rapidamente decrescentes. Soluções exatas não abordam isso porque elas habitam subespaços do espaço de soluções onde a dinâmica não linear completa nunca se engaja.

Em resumo: soluções exatas fornecem soluções suaves explícitas em regimes altamente simétricos. O problema em aberto é se a boa colocação se estende a escoamentos 3D irrestritos. Para a formulação precisa, veja Existência e Suavidade de Navier-Stokes.

Para entender as próprias equações e o que cada termo significa, comece com O Que São as Equações de Navier-Stokes?

Para ver como essas equações são construídas a partir de primeiros princípios, leia Como as Equações de Navier-Stokes São Derivadas.

Para entender quando escoamentos laminares como Poiseuille e Couette se desestabilizam em turbulência, leia Número de Reynolds, Turbulência, e Por Que Pequenas Escalas Importam.

Para entender a questão de um milhão de dólares que soluções exatas não podem responder, leia O Problema do Prêmio do Milênio: Existência e Suavidade.

Próximas páginas recomendadas neste site:

• O Que São as Equações de Navier-Stokes? (o sistema governante e seus termos)
• Como as Equações de Navier-Stokes São Derivadas (do balanço de momento e da lei constitutiva newtoniana à EDP)
• Número de Reynolds, Turbulência, e Por Que Pequenas Escalas Importam (a transição dos regimes laminares descritos aqui para escoamento totalmente turbulento)
• O Problema do Prêmio do Milênio: Existência e Suavidade (a formulação precisa do problema em aberto que soluções exatas não resolvem)