Tutorial de Simulação Eletrotérmica (Parte B): Malha Térmica, Condições de Contorno e Ajustes Experimentais

1. Introdução: escopo e foco da Parte B

Na Parte A você executou uma simulação eletrotérmica, identificou os principais mecanismos de geração de calor (aquecimento por transporte de portadores, aquecimento por recombinação e perdas parasitas), e aprendeu como extrair saídas de temperatura da rede em função da polarização aplicada. Essa primeira parte focou em qual calor é gerado e onde ele aparece no dispositivo. Nesta Parte B, focamos na outra metade do problema: como o calor é transportado para longe do dispositivo. Isso é governado pela malha térmica e, criticamente, pelas condições de contorno térmicas. As condições de contorno codificam como o dispositivo é montado, resfriado ou acoplado a um dissipador na prática, e elas frequentemente têm um impacto de primeira ordem no aumento de temperatura resultante. Duas simulações com comportamento elétrico idêntico e termos de geração de calor idênticos podem produzir perfis de temperatura muito diferentes puramente devido a diferentes hipóteses de contorno térmico.

Você verá, portanto, por que a simulação eletrotérmica normalmente usa um domínio térmico maior que o domínio elétrico, como condições de contorno efetivas são usadas para representar dissipadores macroscópicos sem malhá-los explicitamente, e por que quantidades de material dependentes da temperatura são pré-computadas em uma grade discreta de temperatura. Este tutorial continua a usar um exemplo de diodo orgânico porque ele é experimentalmente fundamentado e inclui efeitos mensuráveis de autoaquecimento, mas o fluxo de trabalho é completamente geral. Os mesmos conceitos se aplicam a semicondutores inorgânicos, dispositivos de filme fino, eletrônica de potência, fotodetectores e qualquer estrutura em que a dissipação de potência elétrica realimente o transporte por meio da temperatura.

Ao final da Parte B, você deverá entender:

• por que a malha térmica é distinta da malha elétrica,
• como as condições de contorno térmicas controlam a extração de calor,
• como a discretização da temperatura é usada internamente para eficiência, e
• por que o autoaquecimento às vezes deve ser incluído para corresponder aos dados experimentais.

2. Malha térmica e condições de contorno

Comece na janela principal de simulação e navegue até a faixa Térmica. Esta é a mesma faixa apresentada na Parte A: ela contém os controles de ativação das fontes de aquecimento e os controles principais de configuração (parâmetros térmicos, malha térmica e condições de contorno). A faixa é mostrada em ??.

2.1 A malha térmica e por que ela não é a malha elétrica

Clique em Malha térmica para abrir o editor de malha térmica. À primeira vista ele se parece com os editores de malha que você viu em outros lugares no OghmaNano, mas ele resolve um problema físico diferente: difusão de calor. O ponto principal é que a simulação eletrotérmica usa duas malhas, porque os problemas elétrico e térmico geralmente vivem em escalas de comprimento fundamentalmente diferentes:

• A malha elétrica resolve o transporte (drift–diffusion + Poisson + recombinação/armadilhas). Ela é tipicamente focada onde a ação elétrica é mais forte: a camada ativa e as interfaces próximas, onde campos e densidades de portadores variam rapidamente. Essas variações podem ocorrer em escalas de nanômetros a mícrons.
• A malha térmica resolve a equação de calor da rede. O calor não para na camada ativa: ele se espalha por toda a pilha, para os contatos, substratos, encapsulamento e em direção ao que remove o calor do dispositivo. Esses caminhos são comumente de escalas de milímetros a centímetros.

Essa incompatibilidade é a razão pela qual uma única “malha compartilhada” normalmente é a abstração errada. O problema elétrico precisa de alta resolução na região ativa; o problema térmico precisa de um domínio grande o suficiente para representar o caminho de fluxo de calor. Por exemplo:

• Um diodo de filme fino pode ter uma espessura de camada ativa de 100–300 nm, mas o substrato pode ter 0.5–1 mm de espessura (vidro) e o bloco ou estágio de montagem pode ter centímetros.
• Um contato metálico pode ter eletricamente apenas dezenas de nanômetros de espessura, mas pode ser o espalhador lateral de calor dominante termicamente porque sua condutividade térmica é grande em comparação com camadas orgânicas ou de óxido.

Neste exemplo, você pode ver que a altura Y do domínio térmico abrange a espessura total da pilha do dispositivo em vez de estar restrita à região eletricamente ativa. Isso é típico: o transporte térmico envolve toda camada que conduz calor, não apenas as camadas onde a corrente flui. Em seguida, configuramos como o calor pode sair deste domínio térmico usando condições de contorno. A malha define a região na qual o calor pode se difundir; as condições de contorno definem o que acontece nas bordas dessa região.

2.2 Condições de contorno: faces isolantes e um dissipador térmico efetivo

Abra o editor de condições de contorno mostrado em ??. Neste exemplo, a maioria das faces está definida como Neumann (frequentemente mostrado como “Neumann (==0)”). Fisicamente, uma condição de contorno de Neumann especifica o fluxo de calor normal no contorno. Quando ela é definida como zero, impõe:

\[ -k \nabla T \cdot \hat{n} = 0 \]

isto é, nenhum calor flui através desse contorno. Essas faces são tratadas como termicamente isolantes. Isso não significa que o dispositivo seja termicamente isolado no geral; apenas informa ao solucionador térmico que esses contornos não fazem parte do caminho pretendido de remoção de calor. A exceção aqui é a face y_max, que está definida como Dissipador térmico. Esse contorno especifica uma temperatura do dispositivo em y_max (aqui em torno de 300 K), juntamente com uma condutividade do dissipador térmico efetiva e um comprimento do dissipador térmico (da ordem de um milímetro neste exemplo).

Conceitualmente, isso é um modelo efetivo de remoção de calor. Ele permite que o calor flua para fora da pilha de dispositivo simulada em direção a um dissipador não explicitamente malhado, sem forçar a malha térmica a incluir um objeto macroscópico. Isso importa porque dissipadores reais são tipicamente ordens de magnitude maiores do que dispositivos de filme fino: malhar explicitamente um dissipador em escala de milímetros a centímetros com resolução submícron seria computacionalmente ineficiente e fisicamente desnecessário. Vale enfatizar que, em problemas térmicos, as condições de contorno determinam em grande parte quão quente o dispositivo se torna. Um caminho ruim de extração de calor pode fazer a temperatura subir muito rapidamente, enquanto um caminho eficiente pode manter a temperatura próxima à ambiente mesmo sob dissipação significativa de potência.

Uma analogia útil é pensar no dispositivo como uma banheira sendo enchida com água. Os termos de geração de calor são a torneira: eles despejam energia no sistema. As condições de contorno térmicas são o ralo:

• Se o ralo estiver totalmente aberto (excelente extração de calor), o nível de água permanece baixo.
• Se o ralo estiver parcialmente bloqueado (extração de calor limitada), o nível de água sobe.
• Se o ralo estiver fechado (sem extração de calor), a banheira eventualmente transbordará.

Nessa analogia, o nível de água é a temperatura da rede. O propósito das condições de contorno térmicas, portanto, não é cosmético nem secundário: elas codificam o ambiente físico que, em última análise, controla se o dispositivo opera frio, morno ou catastroficamente quente.

2.3 “Pontos de malha” de temperatura (por que a temperatura é discretizada)

O editor de malha térmica também mostra uma faixa de temperatura e um número de pontos de temperatura. Isso não é a malha térmica espacial: é a grade de temperatura usada para pré-computações. Muitas partes do modelo elétrico/térmico exigem que quantidades dependentes da temperatura sejam avaliadas repetidamente durante a solução acoplada, e é computacionalmente caro recomputá-las do zero em toda temperatura intermediária. Por exemplo, dependendo das estatísticas escolhidas e do modelo de densidade de estados, o solucionador pode precisar de relações dependentes da temperatura para densidades de portadores, funções quase-termodinâmicas e quantidades relacionadas de consulta. Em vez de avaliar essas relações “ao vivo” toda vez que o laço acoplado atualiza \(T_L\), o OghmaNano pré-computa tabelas de base em um número finito de pontos de temperatura e então interpola entre elas durante a execução. Neste exemplo, as tabelas são geradas em 7 pontos entre 290 K e 350 K.

Isso faz da faixa de temperatura uma escolha de modelagem: ela deve abranger confortavelmente as temperaturas esperadas durante o autoaquecimento. Se o dispositivo aquecer além da temperatura superior, você corre o risco de artefatos de extrapolação (ou saturação, dependendo das configurações), o que raramente é o desejado. Como regra prática: escolha uma faixa que cubra as temperaturas operacionais previstas com margem.

3. Parâmetros térmicos (condutividade térmica e tempos de relaxação)

Na simulação eletrotérmica, as condições de contorno definem como o calor pode sair do dispositivo, mas os parâmetros do material definem como o calor se move dentro da pilha. Na prática, uma das quantidades mais importantes que controlam o aumento de temperatura e a formação de hotspots é a condutividade térmica \(k\) (às vezes escrita como \(\kappa\)).

Para visualizar ou editar esses parâmetros, abra a faixa Térmica e clique em Parâmetros térmicos (o botão marcado com \(k\) / \(\kappa\)). Isso abre o editor de parâmetros térmicos camada por camada mostrado em ??.

O editor expõe três parâmetros principais para cada camada:

• Condutividade térmica \(k\): governa quão eficientemente o calor se difunde através dessa camada. Camadas com baixo \(k\) atuam como gargalos térmicos e podem aumentar fortemente a elevação prevista da temperatura da rede.
• Tempo de relaxação de elétrons e tempo de relaxação de lacunas: estes só são necessários ao usar o modelo hidrodinâmico / de transporte de energia, no qual as temperaturas dos portadores podem se afastar da temperatura da rede. No modelo padrão de calor da rede eles não são necessários porque as temperaturas dos portadores ficam travadas em \(T_L\).

Na maioria dos fluxos de trabalho, as condutividades térmicas são escolhidas a partir de valores publicados (ou medidos quando disponíveis) e então refinadas apenas se houver uma justificativa experimental clara. Juntamente com as condições de contorno, \(k\) é um dos principais determinantes da temperatura absoluta do dispositivo e da estrutura espacial de hotspots.

4. Como o solucionador eletrotérmico acoplado funciona (e por que ele é mais lento)

Uma simulação eletrotérmica no OghmaNano consiste em dois solucionadores acoplados:

• um solucionador elétrico totalmente acoplado (drift–diffusion e Poisson, incluindo recombinação e armadilhas quando habilitadas), e
• um solucionador térmico que resolve a equação de difusão de calor para o campo de temperatura da rede \(T_L\) usando termos de geração de calor calculados a partir da solução elétrica.

Em uma dada tensão aplicada, o solucionador não executa essas etapas uma vez em sequência. Em vez disso, ele executa um laço externo autoconsistente:

1. resolver o problema de transporte elétrico usando o campo de temperatura atual,
2. calcular os termos de geração de calor (transporte, recombinação, parasita e, opcionalmente, óptico),
3. resolver o problema de difusão térmica para atualizar \(T_L\), e
4. repetir até que tanto os resíduos elétricos quanto os térmicos estejam convergidos.

Esse acoplamento iterativo é a razão pela qual simulações eletrotérmicas são mais lentas do que simulações puramente elétricas: para cada ponto de polarização, a solução de Newton elétrica e a solução térmica podem ambas ser executadas múltiplas vezes antes que o sistema acoplado se estabilize. O resultado é uma solução multifísica genuinamente autoconsistente, em vez de um cálculo elétrico com temperatura adicionada como etapa de pós-processamento.

O acoplamento é essencial porque os problemas elétrico e térmico normalmente vivem em escalas de comprimento físico muito diferentes. O transporte elétrico é dominado pela região ativa, onde campos, densidades de portadores e taxas de recombinação variam rapidamente. O transporte térmico deve levar em conta todo o caminho de fluxo de calor, incluindo contatos, substratos, encapsulamento e dissipadores térmicos. É por isso que o OghmaNano fornece controles independentes para a malha térmica e as condições de contorno térmicas.

Também é importante reconhecer que a estrutura eletrotérmica naturalmente envolve múltiplas temperaturas. Em geral, o modelo distingue entre a temperatura da rede \(T_l\), que descreve a temperatura da rede atômica, a temperatura dos elétrons \(T_e\), e a temperatura das lacunas \(T_h\).

Elétrons e lacunas são populações de portadores móveis que podem, em princípio, ter suas próprias energias locais, assim como têm seus próprios níveis de quasi-Fermi. No modelo eletrotérmico padrão usado aqui, as temperaturas dos elétrons e das lacunas são travadas à temperatura da rede, de modo que todos os subsistemas compartilham um único campo de temperatura.

Em regimes extremos, o OghmaNano também suporta um modelo de transporte hidrodinâmico no qual as temperaturas dos elétrons e das lacunas podem desviar da temperatura da rede. Esse modelo avançado é descrito em outro lugar e só é necessário em situações especializadas. Para a maioria das simulações de dispositivos, incluindo este tutorial, o modelo eletrotérmico baseado na rede é apropriado e suficiente.

5. Autoaquecimento e comparação com experimento

Este dispositivo de exemplo foi comparado com dados experimentais usando o modelo eletrotérmico. Nenhum ajuste é realizado neste tutorial. Em vez disso, a comparação existente é usada para ilustrar como o autoaquecimento altera a resposta elétrica do dispositivo, e por que um tratamento eletrotérmico acoplado é necessário em correntes operacionais mais altas.

O fluxo de trabalho de comparação é acessado pela barra de ferramentas Automação, mostrada em ??. Abrir a ferramenta ajustar ao experimento exibe os dados experimentais de JV para este dispositivo, como mostrado em ??. Essa curva representa o comportamento medido do dispositivo sob polarização.

Execute uma iteração. Após uma única solução eletrotérmica, a curva JV simulada é sobreposta aos dados experimentais, como mostrado em ??. Neste caso, a concordância é próxima, indicando que o modelo captura a física dominante de transporte e térmica do dispositivo.

O papel do autoaquecimento fica claro ao repetir a mesma comparação com o modelo térmico desabilitado:

1. Retorne à faixa Térmica.
2. Desabilite o modelo térmico.
3. Execute novamente uma iteração da comparação.

A curva JV resultante não corresponde mais aos dados experimentais, como mostrado em ??. O desvio é mais pronunciado em polarizações mais altas, onde as densidades de corrente são grandes e o autoaquecimento por transporte de portadores se torna significativo.

Essa comparação demonstra um resultado central da modelagem eletrotérmica: uma vez que os dispositivos operam em regimes de alta densidade de corrente, a dissipação de potência elétrica realimenta o transporte por meio do aumento de temperatura. Simulações elétricas a temperatura fixa não conseguem capturar esse efeito, enquanto o modelo eletrotérmico acoplado consegue.