2- Apresente formalmente as equações de balanço com as condições iniciais e de contorno. Explique o significado de cada parcela nos balanços e das condições.

Passo 1

O enunciado pede que eu “apresente formalmente as equações de balanço com as condições iniciais e de contorno” e que eu explique o significado de cada parcela nos balanços e das condições.
Informações dadas (implícitas no enunciado): devo fornecer as equações de balanço gerais (conservação) e explicar cada termo e os tipos de condições iniciais e de contorno que costumam ser usados. Não há figura nem dados numéricos adicionais.

Passo 2

Primeiro vamos escrever a forma integral geral de um balanço para uma grandeza escalar extensiva por unidade de massa \varphi num volume de controle V com fronteira \partial V. A forma integral é pedagógica porque explicita fluxos através da superfície e fontes internas:
\frac{d}{dt}\int_{V} \rho \varphi \, dV +\int_{\partial V} \rho \varphi\, \mathbf{v}\cdot\mathbf{n}\, dS = -\int_{\partial V} \mathbf{J}_\varphi\cdot\mathbf{n}\, dS + \int_{V} s_\varphi \, dV.
Explicação de cada parcela:
  • \displaystyle \frac{d}{dt}\int_{V} \rho \varphi \, dV: taxa de acumulação da quantidade \varphi dentro do volume V. Aqui \rho é densidade (unidade: massa/volume) e \rho\varphi é a densidade volumétrica da grandeza conservada.
  • \displaystyle \int_{\partial V} \rho \varphi\, \mathbf{v}\cdot\mathbf{n}\, dS: fluxo convectivo através da superfície — transporte de \varphi por movimento do fluido com velocidade \mathbf{v}. \mathbf{n} é a normal externa.
  • \displaystyle -\int_{\partial V} \mathbf{J}_\varphi\cdot\mathbf{n}\, dS: fluxo não-convectivo (difusivo ou por condução) através da superfície; o sinal negativo segue a convenção de que um fluxo no sentido externo reduz a quantidade no volume. \mathbf{J}_\varphi é o vetor fluxo difusivo (por exemplo, fluxo de massa por difusão ou fluxo de calor).
  • \displaystyle \int_{V} s_\varphi \, dV: fonte volumétrica interna (geração ou consumo) de \varphi por unidade de volume.

Passo 3

Transformando a forma integral em forma local (diferencial) usando o teorema da divergência e assumindo regularidade, obtemos a forma conservativa local:
\frac{\partial}{\partial t}(\rho \varphi) + \nabla\cdot(\rho \varphi\, \mathbf{v}) = -\nabla\cdot\mathbf{J}_\varphi + s_\varphi.
Interpretação local:
  • \displaystyle \frac{\partial}{\partial t}(\rho \varphi): variação temporal local da densidade volumétrica de \varphi.
  • \displaystyle \nabla\cdot(\rho \varphi\, \mathbf{v}): divergência do fluxo convectivo — representa saída/entrada local de \varphi por convecção.
  • \displaystyle -\nabla\cdot\mathbf{J}_\varphi: contribuição local de difusão/condutividade.
  • \displaystyle s_\varphi: geração local.
Exemplo de lei constitutiva (frequentemente usada):
  • Para difusão linear: \displaystyle \mathbf{J}_\varphi = -\Gamma \nabla \varphi, onde \Gamma é coeficiente difusivo (ex.: coeficiente de difusão ou condutividade térmica); a lei inverte o gradiente: fluxo sai onde \varphi é maior.

Passo 4

Agora as equações fundamentais de fluxo contínuo (as mais comuns em provas): massa, quantidade de movimento (momentum) e energia.
  1. Continuidade (conservação de massa). Forma diferencial:
\frac{\partial \rho}{\partial t} + \nabla\cdot(\rho \mathbf{v}) = 0.
Significado:
  • Termo \partial_t \rho: variação local da densidade.
  • Termo \nabla\cdot(\rho\mathbf{v}): divergência do fluxo de massa.
    Caso incompressível (densidade constante \rho= const): simplifica para
\nabla\cdot\mathbf{v} = 0,
que significa “volume defluido conservado” — não há variação local de volume.

Passo 5

  1. Balanço de quantidade de movimento (equação de Navier–Stokes/Cauchy em forma geral). Forma convectiva/diferencial útil:
\rho\left(\frac{\partial \mathbf{v}}{\partial t} + \mathbf{v}\cdot\nabla\mathbf{v}\right) = \nabla\cdot\boldsymbol{\sigma} + \rho \mathbf{g}.
Aqui:
  • \rho(\partial_t \mathbf{v} + \mathbf{v}\cdot\nabla\mathbf{v}): termo de inércia; \partial_t \mathbf{v} é a aceleração local e \mathbf{v}\cdot\nabla\mathbf{v} é a aceleração convectiva (mudança de velocidade ao mover-se no campo).
  • \nabla\cdot\boldsymbol{\sigma}: forças de superfície por unidade de volume, com \boldsymbol{\sigma} tensor de tensões (stress). Decompõe-se tipicamente em
\boldsymbol{\sigma} = -p\,\mathbf{I} + \boldsymbol{\tau},
onde p é a pressão (tensão normal) e \boldsymbol{\tau} é o tensor de tensões viscosas (tensão de cisalhamento).
  • \rho\mathbf{g}: força volumétrica (por exemplo, gravidade).
Para um fluido newtoniano isotrópico (modelo comum):
\boldsymbol{\tau} = \mu\left(\nabla\mathbf{v} + (\nabla\mathbf{v})^T\right) - \tfrac{2}{3}\mu (\nabla\cdot\mathbf{v})\,\mathbf{I},
com \mu viscosidade dinâmica. Assim fica explícito o termo viscoso (difusivo de momento).

Passo 6

  1. Balanço de energia. Existem formas em termos de energia total ou energia interna. Forma em termos de energia interna específica e (por unidade de massa) é:
\rho\left(\frac{\partial e}{\partial t} + \mathbf{v}\cdot\nabla e\right) = \boldsymbol{\sigma}:\nabla\mathbf{v} - \nabla\cdot\mathbf{q} + \rho r.
Explicação:
  • \rho(\partial_t e + \mathbf{v}\cdot\nabla e): variação material da energia interna por unidade de volume.
  • \boldsymbol{\sigma}:\nabla\mathbf{v}: trabalho das tensões por unidade de volume; inclui contribuição -p\nabla\cdot\mathbf{v} (compressão) e dissipação viscosa \boldsymbol{\tau}:\nabla\mathbf{v}.
  • -\nabla\cdot\mathbf{q}: divergência do fluxo de calor; \mathbf{q} é o vetor fluxo térmico (ex.: lei de Fourier \mathbf{q}=-k\nabla T, com k condutividade térmica).
  • \rho r: fonte volumétrica de energia (ex.: aquecimento interno, reações químicas).
Observação: forma total de energia inclui também energia cinética \tfrac12|\mathbf{v}|^2; às vezes usamos a equação de energia total para fechar modelos compressíveis.

Passo 7

Condições iniciais (CI) — o que precisamos dar em t=0:
  • Para problemas dependentes do tempo, deve-se fornecer o estado inicial em todo o domínio \Omega:
\rho(\mathbf{x},0)=\rho_0(\mathbf{x}),\qquad \mathbf{v}(\mathbf{x},0)=\mathbf{v}_0(\mathbf{x}),\qquad T(\mathbf{x},0)=T_0(\mathbf{x}),
ou para demais campos escalares \varphi(\mathbf{x},0)=\phi_0(\mathbf{x}). Estas funções iniciais têm de ser compatíveis com certas BCs (por exemplo, se a parede tem \mathbf{v}=0 então \mathbf{v}_0 deve satisfazer isso na fronteira).

Passo 8

Condições de contorno (BC) — tipos e significado físico (as que aparecem com mais frequência):
  • Dirichlet (valor imposto): especifica o valor do campo na fronteira. Ex.: velocidade fixa numa parede (no–slip)
\mathbf{v}(\mathbf{x}\in\partial\Omega,t)=\mathbf{v}_b(\mathbf{x},t).
Significado: a própria variável é fixada na fronteira (ex.: temperatura na parede).
  • Neumann (fluxo imposto): especifica o fluxo normal do campo. Ex.: fluxo de calor prescrito
-\mathbf{q}\cdot\mathbf{n} = q''(\mathbf{x},t)\quad\text{em }\partial\Omega.
Significado: especifica quanto da grandeza entra/saí pela fronteira (ex.: calor por unidade de área).
  • Robin (misto / convectivo): combinação linear valor–fluxo, comum em troca convectiva com ambiente
- \mathbf{q}\cdot\mathbf{n} = h\big(T - T_\infty\big)
onde h é coeficiente convectivo e T_\infty temperatura ambiente. É útil quando não sabemos T na outra face, só o ambiente fora.
  • Condição de tração (para momentum): especifica a tensão/força por unidade de área na fronteira:
\boldsymbol{\sigma}\cdot\mathbf{n} = \mathbf{t}_b(\mathbf{x},t).
Ex.: parede sujeita a uma pressão externa ou tração conhecida.
  • Condição impermeabilidade (parede sólida): componente normal da velocidade zero:
\mathbf{v}\cdot\mathbf{n} = 0.
  • Condição no–slip (fluido viscoso): velocidade tangencial na parede igual à velocidade da parede (geralmente zero):
\mathbf{v} = \mathbf{0}\quad\text{na parede}.
  • Periódicas: valor em uma face igual ao da face oposta; usadas em problemas repetitivos.
Observações práticas:
  • A escolha e o número de condições dependem do tipo da PDE: para equações parabolóides (ex.: difusão) precisa-se de CI + BC apropriadas; para equações elípticas (ex.: Poisson), CI não se aplica mas BCs são essenciais; para equações hiperbólicas (ex.: onda, advecção pura) as condições devem seguir a direção das características (mais delicado).
  • Em problemas físicos é comum misturar tipos: valor de temperatura numa parte da fronteira, fluxo numa outra, e condição de isolamento noutra.

Passo 9

Exemplo concreto e simples (1D, escalar advecção-difusão):
PDE:
\frac{\partial \varphi}{\partial t} + u\frac{\partial \varphi}{\partial x} = D\frac{\partial^2 \varphi}{\partial x^2} + s(x,t),
com u velocidade (constante), D difusividade, s fonte. CI:
\varphi(x,0)=\phi_0(x),\quad 0<x<L.
BCs possíveis:
  • Dirichlet: \varphi(0,t)=\phi_L(t),\; \varphi(L,t)=\phi_R(t).
  • Neumann: \partial_x\varphi(0,t)=g_0(t),\; \partial_x\varphi(L,t)=g_L(t).
  • Robin (convectivo): -D\partial_x\varphi(L,t)=h\big(\varphi(L,t)-\phi_\infty\big).
Interpretação rápida: termo \partial_t\varphi acumula, u\partial_x\varphi transporta sem trocar entre vizinhos, D\partial_{xx}\varphi suaviza a solução (difusão), s cria/consome \varphi.

Passo 10

Observações finais e dicas para a prova (estilo amigo):
  • Sempre identifique a grandeza conservada e escolha se usa forma integral (bom para interpretar fluxos) ou forma diferencial (bom para manipular analiticamente).
  • Anote claramente: (i) quais campos precisam de CI; (ii) que tipo de BC cada campo aceita; (iii) verifique compatibilidade entre CI e BC na fronteira no instante inicial.
  • Preste atenção nos sinais: fluxos externos que saem diminuem a quantidade no volume (sinal negativo na forma integral).
  • Se a questão pedir simplificações (fluido incompressível, regime estacionário, sem viscosidade), aplique-as explicitamente nas equações para justificar termos nulos.
  • Em problemas numéricos/teóricos, escrever a forma conservativa \partial_t(\rho\varphi)+\nabla\cdot(\rho\varphi\mathbf{v}+\mathbf{J}_\varphi)=s_\varphi ajuda a garantir conservação.

Resposta

\text{Forma integral geral (para } \varphi\text{):}\quad \frac{d}{dt}\int_{V} \rho \varphi \, dV +\int_{\partial V} \rho \varphi\, \mathbf{v}\cdot\mathbf{n}\, dS = -\int_{\partial V} \mathbf{J}_\varphi\cdot\mathbf{n}\, dS + \int_{V} s_\varphi \, dV.
\text{Forma local (diferencial):}\quad \frac{\partial}{\partial t}(\rho \varphi) + \nabla\cdot(\rho \varphi\, \mathbf{v}) = -\nabla\cdot\mathbf{J}_\varphi + s_\varphi.
\text{Continuidade (massa):}\quad \frac{\partial \rho}{\partial t} + \nabla\cdot(\rho \mathbf{v}) = 0 \quad\text{(incompressível: }\nabla\cdot\mathbf{v}=0\text{).}í
\text{Momentum (Cauchy / Navier–Stokes):}\quad \rho\!\left(\frac{\partial \mathbf{v}}{\partial t} + \mathbf{v}\cdot\nabla\mathbf{v}\right) = \nabla\cdot\boldsymbol{\sigma} + \rho \mathbf{g}, \qquad \boldsymbol{\sigma} = -p\mathbf{I} + \boldsymbol{\tau}.
\text{Energy (interna):}\quad \rho\left(\frac{\partial e}{\partial t} + \mathbf{v}\cdot\nabla e\right) = \boldsymbol{\sigma}:\nabla\mathbf{v} - \nabla\cdot\mathbf{q} + \rho r.
\text{Leis constitutivas típicas:}\quad \mathbf{J}_\varphi=-\Gamma\nabla\varphi,\qquad \mathbf{q}=-k\nabla T,\qquad \boldsymbol{\tau}=\mu\big(\nabla\mathbf{v}+(\nabla\mathbf{v})^T\big)-\tfrac{2}{3}\mu(\nabla\cdot\mathbf{v})\mathbf{I}.í
\text{Condições iniciais (ex.):}\quad \rho(\mathbf{x},0)=\rho_0(\mathbf{x}),\;\mathbf{v}(\mathbf{x},0)=\mathbf{v}_0(\mathbf{x}),\;T(\mathbf{x},0)=T_0(\mathbf{x}).çõ
\text{Condições de contorno (exemplos):}\quad \begin{cases} \text{Dirichlet: } & \varphi|_{\partial\Omega}=\phi_b \;\;(\text{valor imposto}),\\[2pt] \text{Neumann: } & -\mathbf{J}_\varphi\cdot\mathbf{n}=j_b \;\;(\text{fluxo imposto}),\\[2pt] \text{Robin: } & -\mathbf{J}_\varphi\cdot\mathbf{n}=h(\varphi-\phi_\infty)\;\;(\text{misto/convectivo}),\\[2pt] \text{No–slip (momentum): } & \mathbf{v}=\mathbf{0}\;\text{na parede},\\[2pt] \text{Tração: } & \boldsymbol{\sigma}\cdot\mathbf{n}=\mathbf{t}_b. \end{cases}çõçã]

Exercícios relacionados