Questão 6: Demonstre, em detalhe, que ao admitirmos escoamento irrotacional, podemos, a partir da equação da continuidade para um fluido incompressível, obter a seguinte expressão para o potencial de velocidade:
\nabla^{2} \varphi=0, \text { em que } \vec{u}=\nabla \varphi
Explique qual as vantagens associadas a essa aproximação e ainda descreva como os campos de velocidade e pressão podem ser obtidos.

Passo 1

O enunciado pede que eu demonstre, com detalhes, que ao assumir escoamento irrotacional e incompressível, a equação da continuidade leva à seguinte equação para o potencial de velocidade:
\nabla^{2}\varphi=0,\quad\text{onde}\quad \vec{u}=\nabla\varphi.
Além da demonstração, é pedido que eu explique as vantagens dessa aproximação (potencial) e descreva como, na prática, obtém-se os campos de velocidade e de pressão a partir de \varphi.
Informações dadas / condições do problema:
  • O fluido é incompressível (condição de continuidade).
  • Admite-se escoamento irrotacional.
  • Define-se um potencial escalar \varphi tal que \vec{u}=\nabla\varphi (essa é a relação a ser usada).

Passo 2

Vou começar lembrando e explicando a equação da continuidade para fluido incompressível.
A equação da continuidade (condição de conservação de massa) para um fluido incompressível está escrita por:
\nabla\cdot\vec{u}=0.
Explicação didática:
  • Aqui \vec{u}(\vec{x},t) é o campo vetorial de velocidade do fluido (cada ponto do espaço e tempo tem um vetor velocidade).
  • O operador \nabla\cdot é a divergência; fisicamente, \nabla\cdot\vec{u} mede a taxa de expansão volumétrica local do fluido. Para incompressível, essa taxa é zero: volume local não muda com o tempo.
  • Em coordenadas cartesianas, se \vec{u}=(u_x,u_y,u_z), então \nabla\cdot\vec{u}=\partial_x u_x+\partial_y u_y+\partial_z u_z.

Passo 3

Agora explico o que significa escoamento irrotacional e a consequência topológica necessária.
Definição e interpretação:
  • Escoamento irrotacional significa que a vorticidade é nula:
\nabla\times\vec{u}=0.
  • O operador \nabla\times é o rotacional (curl); que ser zero quer dizer que não há “rotação local” do fluido em torno de pontos (não há vorticidade).
  • Resultado importante da análise vetorial: se \nabla\times\vec{u}=0 numa região simplesmente conexa (sem “buracos”), então existe uma função escalar \varphi(\vec{x},t) tal que
\vec{u}=\nabla\varphi.
Explicação da condição de existência:
  • A exigência de domínio simplesmente conexo evita problemas com campos cuja circulação em torno de laços não nulos impede a existência global de um potencial escalar. Em muitos problemas físicos (domínio externo a um corpo sem topologia complicada) essa condição é satisfeita.

Passo 4

Mostro por que \vec{u}=\nabla\varphi é coerente com a irrotacionalidade.
Propriedade vetorial útil:
  • O rotacional do gradiente de qualquer função é nulo:
\nabla\times(\nabla\varphi)=0,
o que confirma que \vec{u}=\nabla\varphi implica escoamento irrotacional. Isso é só uma verificação da consistência: gradiente tem curl zero.

Passo 5

Agora aplico a equação da continuidade à representação por potencial e obtenho a equação desejada.
Começo com a continuidade:
\nabla\cdot\vec{u}=0.
Substituo \vec{u}=\nabla\varphi:
\nabla\cdot(\nabla\varphi)=0.
Definição do operador laplaciano:
  • O operador \nabla\cdot(\nabla\cdot) é o laplaciano \nabla^{2}. Assim obtemos:
\nabla^{2}\varphi=0.
Conclusão lógica:
  • Portanto, a partir da incompressibilidade e da hipótese de irrotacionalidade, o potencial de velocidade \varphi satisfaz a equação de Laplace (é função harmônica) na região do fluido onde as hipóteses valem.

Passo 6

Comento propriedades matemáticas e físicas importantes da equação \nabla^{2}\varphi=0 e as condições de contorno que aparecem em problemas reais.
Propriedades e tipos de contorno:
  • A equação de Laplace é linear e homogênea; por isso vale o princípio da superposição: soma de soluções é solução.
  • Condições de contorno típicas:
    • Condição de Dirichlet: especifica \varphi na fronteira (por exemplo, potencial em infinito ou em uma superfície).
    • Condição de Neumann: especifica a derivada normal \partial_n \varphi=\vec{n}\cdot\nabla\varphi na fronteira; fisicamente, \partial_n\varphi é a componente normal da velocidade, ou seja, velocidade normal imposta na fronteira.
  • Teoremas de unicidade: com condição de Dirichlet bem colocada, a solução de Laplace é única. No caso de condição de Neumann a solução é única até uma constante (o potencial pode ser somado por uma constante sem alterar \vec{u}).

Passo 7

Vantagens práticas da aproximação potencial (escoamento irrotacional + incompressível).
Vantagens explicadas:
  • Redução do problema: de um campo vetorial desconhecido \vec{u} (três componentes) passa-se para uma única função escalar \varphi. Isso simplifica muito a análise e os cálculos.
  • Linearidade: \nabla^{2}\varphi=0 é linear; permite superposição de soluções elementares (fontes, sumidouros, dipolos, fluxo uniforme) para construir soluções mais complexas.
  • Arsenal teórico: funções harmônicas têm propriedades úteis (média, máximo-mínimo, expansão em funções harmônicas, representações integrais com funções de Green), o que facilita análises analíticas e numéricas.
  • Eficiência numérica: muitos métodos (método de elementos de contorno, BEM) são particularmente eficientes para Laplace porque só discretizam a fronteira.
  • Boa aproximação: para escoamentos externos de fluidos de alta razão de Reynolds (fluido quase inviscido fora das camadas limite), a aproximação potencial costuma ser qualitativamente correta longe das superfícies onde há separação/viscosidade.
Limitações (importante saber onde não usar):
  • Não representa vorticidade: não serve para escoamentos com giro ou regiões onde a viscosidade gera vorticidade (ex.: camadas-limite com atrito, separação).
  • Não captura efeitos viscosos e perda de energia por atrito.
  • Pode falhar perto de superfícies onde ocorre descolamento (separation).

Passo 8

Como obter, passo a passo, os campos de velocidade e pressão a partir de \varphi.
  1. Resolver a equação de Laplace:
  • Resolver \nabla^{2}\varphi=0 na região do fluido com condições de contorno apropriadas (em frente a um corpo: velocidade normal nula na superfície do corpo, condições no infinito etc.). Técnicas: separação de variáveis (casos simétricos), séries de funções harmônicas, funções de Green, métodos numéricos (diferenças finitas, elementos de contorno).
  1. Calcular a velocidade:
\vec{u}=\nabla\varphi.
  • Isso significa calcular as derivadas parciais de \varphi. Em coordenadas cartesianas \vec{u}=(\partial_x\varphi,\partial_y\varphi,\partial_z\varphi).
Exemplo simples (para fixar a ideia):
  • Escoamento uniforme na direção x: \varphi=U x. Verifica-se \nabla^{2}\varphi=0 e \vec{u}=(U,0,0).
  1. Calcular a pressão (assumindo fluxo invíscido e, portanto, uso da equação de Euler/Bernoulli):
  • Equação de Euler (sem viscosidade) é
\rho\bigl(\partial_t\vec{u}+(\vec{u}\cdot\nabla)\vec{u}\bigr)=-\nabla p + \rho\vec{g},
onde \rho é a densidade e \vec{g} a aceleração da gravidade (que pode ser escrita como -\nabla(gz) se z for a vertical).
  • Para escoamento irrotacional temos a identidade vetorial (usar que \nabla\times\vec{u}=0):
(\vec{u}\cdot\nabla)\vec{u}=\nabla\!\left(\tfrac{1}{2}|\vec{u}|^{2}\right).
  • Substituindo \vec{u}=\nabla\varphi e reorganizando, obtém-se a forma integral de Bernoulli (válida globalmente para escoamento irrotacional):
\nabla\!\left(\partial_t\varphi+\tfrac{1}{2}|\nabla\varphi|^{2}+\frac{p}{\rho}+gz\right)=0.
Isto implica que, em cada instante,
\partial_t\varphi+\tfrac{1}{2}|\nabla\varphi|^{2}+\frac{p}{\rho}+gz=f(t),
onde f(t) é uma função apenas do tempo (constante espacialmente; para escoamento permanente f(t) é constante). Variantes usuais:
  • Caso estacionário ( \partial_t\varphi=0):
\tfrac{1}{2}|\nabla\varphi|^{2}+\frac{p}{\rho}+gz=\text{const}.
  • Assim, conhecendo \varphi (logo \nabla\varphi), obtém-se p por Bernoulli até uma constante de referência.
Explicação das variáveis:
  • \rho: densidade do fluido.
  • p: pressão.
  • g: aceleração da gravidade.
  • z: coordenada vertical.
  • |\nabla\varphi|^{2}=\vec{u}\cdot\vec{u}=u_x^2+u_y^2+u_z^2.

Passo 9

Observações práticas e resumo do procedimento para um problema típico de potencial.
Procedimento prático:
  • (i) Formular condições de contorno físicas (por exemplo, impermeabilidade na superfície de um corpo: \vec{n}\cdot\vec{u}=\vec{n}\cdot\nabla\varphi=0; comportamento no infinito; velocidade prescrita em entradas/saídas).
  • (ii) Resolver \nabla^{2}\varphi=0 com essas condições (método analítico se houver simetria; caso contrário, numérico, preferencialmente método de contorno).
  • (iii) Calcular \vec{u}=\nabla\varphi.
  • (iv) Calcular p por Bernoulli (escolhendo a forma estacionária ou não estacionária adequada).
Dicas úteis:
  • Para escoamentos externos, é comum representar a solução como combinação linear de soluções elementares (fluxo uniforme + dipolos + vórtices + fontes) e ajustar coeficientes para satisfazer condições de contorno.
  • Em problemas com condição de velocidade normal conhecida na parede, a formulação em termos de Neumann para \varphi é natural; atenção à unicidade (potencial fixado até constante).

Passo 10

Resumo final pedagógico e uma motivação para o exame.
Resumo didático:
  • Começando de \nabla\cdot\vec{u}=0 (incompressível) e \nabla\times\vec{u}=0 (irrotacional), introduzimos \vec{u}=\nabla\varphi. Aplicando divergência obtemos \nabla^{2}\varphi=0.
  • A solução de \varphi dá diretamente \vec{u} e, via Bernoulli (sob hipótese de fluido invíscido), fornece p.
  • Vantagem prática: redução do problema e acesso a muitas técnicas analíticas e numéricas. Limitação: não captura vorticidade e efeitos viscosos.
Você consegue! Na prova, proceda exatamente como aqui: declare as hipóteses, escreva \vec{u}=\nabla\varphi, aplique divergência e conclua \nabla^{2}\varphi=0; depois diga como extrair \vec{u} e p com Bernoulli — isso é o essencial que os examinadores esperam ver.

Resposta

\nabla^{2}\varphi=0
\vec{u}=\nabla\varphi
Pressão pela equação de Bernoulli para fluxo irrotacional:
\partial_t\varphi+\tfrac{1}{2}|\nabla\varphi|^{2}+\frac{p}{\rho}+gz=f(t),
em particular, no caso estacionário,
\tfrac{1}{2}|\nabla\varphi|^{2}+\frac{p}{\rho}+gz=\text{const}.

Exercícios relacionados