Problemas Térmicos Conjugados (CHT): Fundamentos, Modelagem e Desafios Numéricos

Introdução

Os Problemas Térmicos Conjugados (em inglês, Conjugate Heat Transfer – CHT) surgem quando há interação térmica entre um fluido em escoamento e um sólido, de forma que o campo de temperatura do fluido afeta o sólido e vice-versa.

Esse tipo de problema aparece com frequência em engenharia: trocadores de calor, resfriamento eletrônico, aerodinâmica térmica, sistemas de combustão, turbomáquinas e microdispositivos. Diferentemente de problemas puramente térmicos ou puramente fluidodinâmicos, aqui existe uma forte acoplagem fluido–estrutura térmica, o que torna a modelagem e a solução computacional mais desafiadoras.

1- Abordagens para tratar Problemas Térmicos Conjugados

Existem duas abordagens principais para resolver numericamente problemas de transferência de calor conjugada:

Domínio não particionado (monolítico)

Nessa abordagem, o domínio do sólido e do fluido é tratado como um único domínio computacional, resolvendo simultaneamente todas as equações governantes em um sistema único e fortemente acoplado.

Apesar de ser matematicamente elegante, essa estratégia apresenta:

• Maior complexidade de implementação;
• Sistemas lineares muito grandes e mal condicionados;
• Dificuldade de reutilização de códigos especializados (CFD e condução térmica).

Domínio particionado (partitioned approach)

Na abordagem particionada, os domínios do sólido e do fluido são tratados separadamente, cada um com seu próprio conjunto de equações e, muitas vezes, com solvers distintos.

Essa é considerada a melhor abordagem na prática, especialmente em aplicações industriais e acadêmicas avançadas, pelos seguintes motivos:

• Permite o uso de solvers especializados e maduros para fluido e sólido;
• Maior flexibilidade para refinamento de malha em cada domínio;
• Facilita paralelização e escalabilidade;
• Melhor manutenção e reutilização de código.

Por isso, a grande maioria dos frameworks modernos de CHT adota o domínio particionado.

2- Troca de informações no domínio particionado e risco de divergência

No domínio particionado, cada domínio é resolvido de forma isolada, porém existe uma interface fluido–sólido, onde ocorre a troca de informações térmicas durante o processo iterativo.

Nessa interface, ocorre a transferência de:

• Fluxo de calor do fluido para o sólido;
• Temperatura da parede do sólido para o fluido.

O grande desafio é que, se os parâmetros numéricos dessa troca não forem escolhidos adequadamente, o acoplamento pode se tornar instável, levando à divergência da solução.

Isso é especialmente crítico em problemas transientes, com:

• Altos gradientes térmicos;
• Fortes efeitos advectivos no fluido;
• Grandes diferenças de propriedades térmicas entre sólido e fluido.

3- Equação governante do sólido – Difusão de calor 2D transiente

No domínio sólido, o fenômeno dominante é a difusão térmica. A equação governante 2D transiente, com termo fonte, é dada por:

$$\rho_s c_{p,s} \frac{\partial T_s}{\partial t} = \nabla \cdot \left( k_s \nabla T_s \right) + \dot{q}_s$$

onde:

• $\rho_s$ é a densidade do sólido
• $c_{p,s}$​ é o calor específico
• $k_s$​ é a condutividade térmica
• $T_s$​ é a temperatura no sólido
• $\dot{q}_s$​ é o termo fonte volumétrico de calor

4- Equações governantes do fluido incompressível 2D

Para o fluido incompressível, o problema é governado por três equações fundamentais:

Equação da Continuidade (fluido incompressível)

$$\nabla \cdot \vec{u} = 0$$

O que cada termo significa

• $\vec{u} = (u, v)$
  Vetor velocidade do fluido
• $u$: componente da velocidade na direção $x$
• $v$: componente da velocidade na direção $y$

• $\nabla \cdot \vec{u}$(divergente da velocidade)
  Mede quanto o fluido “se expande” ou “se contrai” localmente.

Interpretação física

• Para fluido incompressível, a densidade é constante.
• Logo, o que entra em um ponto deve sair.
• Não há criação nem destruição de massa.

Em palavras simples:
o volume de fluido é conservado ponto a ponto.

Equação do Balanço da Quantidade de Movimento (Navier–Stokes)

$$\rho \left( \frac{\partial \vec{u}}{\partial t} + \vec{u} \cdot \nabla \vec{u} \right) = – \nabla p + \mu \nabla^2 \vec{u} + \vec{f}$$

Essa equação nada mais é do que a Segunda Lei de Newton aplicada a um fluido:

massa × aceleração = soma das forças

Lado esquerdo — termos inerciais (aceleração)

ρ — densidade do fluido

• Mede a massa por volume.
• Multiplica a aceleração para gerar força inercial.

$\dfrac{\partial \vec{u}}{\partial t}$— aceleração local

• Representa a variação da velocidade no tempo em um ponto fixo do espaço.
• Importante em escoamentos transientes.

Exemplo:
O fluido começa parado e passa a se mover → aceleração local ≠ 0.

$\vec{u} \cdot \nabla \vec{u}$ — aceleração convectiva

• Representa a variação da velocidade devido ao movimento do fluido no espaço.
• Surge porque o fluido se desloca para regiões onde a velocidade é diferente.

Exemplo clássico:
Fluxo acelerando dentro de um bocal → mesmo em regime permanente, esse termo existe.

Esse termo é não linear e o grande responsável pela dificuldade numérica da CFD.

Lado direito — forças atuando no fluido

$– \nabla p$ — força de pressão

• Gradiente de pressão gera força.
• O sinal negativo indica que o fluido acelera da alta para a baixa pressão.

É o principal motor de muitos escoamentos internos.

$\mu \nabla^2 \vec{u}$ — força viscosa

• $\mu$: viscosidade dinâmica do fluido.
• $\nabla^2 \vec{u}$: operador laplaciano da velocidade.

Representa:

• Difusão de quantidade de movimento;
• Resistência ao movimento;
• Efeitos de atrito interno.

Esse termo é dominante em:

• Regime laminar;
• Baixos números de Reynolds.

$\vec{f}$— forças de corpo

• Forças que atuam em todo o volume do fluido.
• Exemplos:
  • Gravidade
  • Forças eletromagnéticas
  • Força de Coriolis (em sistemas rotativos)

Equação do Balanço de Energia Térmica (Advecção–Difusão)

$$\rho_f c_{p,f} \left( \frac{\partial T_f}{\partial t} + \vec{u} \cdot \nabla T_f \right) = k_f \nabla^2 T_f + \Phi + \dot{q}_f$$

• $\Phi$ representa o termo de dissipação viscosa
• $\dot{q}_f$​ é o termo fonte térmico no fluido

5- Estratégia correta de acoplamento térmico no domínio particionado

Na solução particionada de problemas CHT, a troca de informações entre os domínios deve ser mista:

• Um domínio fornece fluxo de calor para o outro;
• O outro domínio fornece temperatura de volta.

Nunca se deve trocar apenas fluxo ou apenas temperatura.

Se apenas fluxos forem impostos, o problema fica mal condicionado.
Se apenas temperaturas forem impostas, o sistema perde consistência física.

Essa escolha mista (fluxo ↔ temperatura) garante:

• Conservação de energia na interface;
• Estabilidade numérica;
• Convergência do acoplamento iterativo.

Conclusão

Problemas Térmicos Conjugados representam um dos maiores desafios da Computação Científica aplicada à Engenharia, justamente por envolverem múltiplos fenômenos físicos fortemente acoplados.

A abordagem por domínio particionado é hoje a mais eficiente e robusta, desde que a troca de informações na interface fluido–sólido seja feita de forma correta e estável. Entender as equações governantes, os mecanismos de acoplamento e os riscos de divergência é essencial para quem trabalha com simulações térmicas avançadas.