Guide pratique sur la simulation électrothermique (partie B) : maillage thermique, conditions aux limites et ajustement aux données expérimentales.

1. Introduction : portée et objectifs de la partie B.

Dans Partie A vous avez effectué une simulation électrothermique, identifié les principaux mécanismes de génération de chaleur (échauffement dû au transport des porteurs de charge, échauffement dû à la recombinaison et pertes parasites), et vous avez appris à déterminer les valeurs de la température du réseau en fonction de la tension appliquée. Cette première partie était axée sur quoi de la chaleur est produite et où il apparaît sur l’appareil. Dans ce cas, Partie Bnous nous concentrons sur l’autre moitié du problème : comment la chaleur est dissipée par l’appareil.Cela est régi par le maillage thermique et, point essentiel, par conditions aux limites thermiquesLes conditions aux limites décrivent la manière dont l’appareil est monté, refroidi ou équipé d’un dissipateur thermique dans la pratique, et elles ont souvent un impact important sur l’augmentation de température qui en résulte. Deux simulations présentant un comportement électrique identique et des termes de génération de chaleur identiques peuvent produire des profils de température très différents, simplement en raison de différentes hypothèses concernant les conditions thermiques aux limites.

Vous comprendrez donc pourquoi la simulation électrothermique utilise généralement un domaine thermique plus vaste que le domaine électrique, comment des conditions aux limites efficaces sont utilisées pour représenter les dissipateurs thermiques macroscopiques sans les mailler explicitement, et pourquoi les propriétés des matériaux dépendant de la température sont précalculées sur une grille de températures discrètes. Ce tutoriel continue d’utiliser un exemple de diode organique, car il est basé sur des données expérimentales et inclut des effets d’auto-échauffement mesurables, mais le processus est entièrement général. Les mêmes principes s’appliquent aux semi-conducteurs inorganiques, aux dispositifs à couches minces, à l’électronique de puissance, aux photodétecteurs et à toute structure dans laquelle la dissipation de puissance électrique influe sur le transport par le biais de la température.

À la fin de la partie B, vous devriez avoir compris ce qui suit :

• pourquoi le réseau thermique est différent du réseau électrique,
• comment les conditions aux limites thermiques régissent l’extraction de chaleur.
• comment la discrétisation de la température est utilisée en interne pour optimiser les performances, et
• pourquoi il est parfois nécessaire de prendre en compte l’échauffement propre afin d’obtenir des résultats conformes aux données expérimentales.

2. Maillage thermique et conditions aux limites.

Commencez dans la fenêtre principale de simulation et accédez au ruban « Thermique ». Il s’agit du même ruban que celui présenté dans la partie A : il contient les commandes pour activer ou désactiver les sources de chaleur et les commandes de configuration principales (paramètres thermiques, maillage thermique et conditions aux limites). Le ruban est illustré dans ??.

2.1 Le treillis thermique et pourquoi il ne s’agit pas d’un treillis électrique.

Cliquez sur « Maillage thermique » pour ouvrir l’éditeur de maillage thermique. Au premier abord, il ressemble aux éditeurs de maillage que vous avez déjà utilisés dans OghmaNano, mais il permet de résoudre un problème physique différent : la diffusion thermique. L’élément essentiel est que la simulation électrothermique utilise deux maillages, car les problèmes électriques et thermiques se posent généralement à des échelles de longueur fondamentalement différentes :

• Le maillage électrique permet de résoudre les équations de transport (dérive-diffusion + Poisson + recombinaison/pièges). Il est généralement appliqué là où l’effet électrique est le plus marqué : dans la couche active et aux interfaces avoisinantes, où les champs et les densités de porteurs varient rapidement. Ces variations peuvent se produire à l’échelle du nanomètre ou du micron.
• Le maillage thermique permet de résoudre l’équation de la chaleur pour les réseaux cristallins. La chaleur ne se limite pas à la couche active : elle se propage dans toute la structure, jusqu’aux contacts, aux substrats, à l’emballage et vers tout élément qui dissipe la chaleur de l’appareil. Ces voies de propagation sont généralement de l’ordre du millimètre à du centimètre.

C’est cette inadéquation qui explique pourquoi l’utilisation d’un seul « maillage partagé » est généralement une approche inadéquate. Le problème électrique nécessite une résolution élevée dans la zone active, tandis que le problème thermique nécessite un domaine suffisamment vaste pour représenter le trajet du flux de chaleur. Par exemple :

• Une diode à couches minces peut avoir une épaisseur de couche active de 100 à 300 nm, mais le substrat peut avoir une épaisseur de 0,5 à 1 mm (en verre), et le bloc de montage ou le support peut mesurer plusieurs centimètres.
• Un contact métallique peut n’avoir qu’une épaisseur de quelques dizaines de nanomètres sur le plan électrique, mais il peut néanmoins être le principal conducteur de chaleur sur le plan latéral, car sa conductivité thermique est élevée par rapport aux couches organiques ou d’oxyde.

Dans cet exemple, vous pouvez constater que la hauteur Y du domaine thermique couvre l’épaisseur totale de l’ensemble du dispositif, et non pas seulement la région électriquement active. C’est typique : le transfert thermique implique toutes les couches qui conduisent la chaleur, et pas seulement celles où circule le courant. Ensuite, nous configurons la manière dont la chaleur peut se dissiper de ce domaine thermique en utilisant des conditions aux limites. Le maillage définit la région dans laquelle la chaleur peut se diffuser ; les conditions aux limites définissent ce qui se passe aux limites de cette région.

2.2 Conditions aux limites : faces isolantes et dissipateur thermique efficace.

Ouvrez l’éditeur des conditions aux limites, comme illustré dans ??. Dans cet exemple, la plupart des faces sont définies avec une condition de Neumann (généralement indiquée par « Neumann (==0) »). Physiquement, une condition aux limites de Neumann spécifie le flux thermique normal à la frontière. Lorsqu’elle est définie à zéro, elle impose ce qui suit :

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

c.-à-d. qu’aucun flux de chaleur ne traverse cette limite. Ces surfaces sont traitées comme des surfaces isolantes sur le plan thermique. Cela ne signifie pas que l’ensemble du dispositif est isolé thermiquement ; cela indique simplement au solveur thermique que ces limites ne font pas partie du trajet prévu pour l’évacuation de la chaleur. L’exception ici est la surface ymax, qui est définie comme étant un dissipateur thermique. Cette limite spécifie une température du dispositif au niveau de ymax (ici, environ 300 K), ainsi qu’une conductivité thermique effective du dissipateur et une longueur du dissipateur (de l’ordre du millimètre dans cet exemple).

D’un point de vue conceptuel, il s’agit d’un modèle efficace pour l’évacuation de la chaleur. Il permet à la chaleur de s’échapper de l’ensemble simulé des dispositifs vers un dissipateur thermique dont la géométrie n’est pas explicitement définie, sans pour autant imposer l’inclusion d’un objet macroscopique dans le maillage thermique. Cet aspect est important, car les dissipateurs thermiques réels sont généralement beaucoup plus grands que les dispositifs à couches minces. Définir explicitement un maillage pour un dissipateur thermique dont les dimensions varient de quelques millimètres à quelques centimètres, avec une résolution inférieure au micron, serait inefficace sur le plan informatique et physiquement inutile. Il convient de souligner que, dans les problèmes de thermique, les conditions aux limites déterminent en grande partie la température atteinte par le dispositif. Un système d’évacuation de chaleur inefficace peut entraîner une augmentation très rapide de la température, tandis qu’un système efficace peut maintenir la température proche de la température ambiante, même en cas de dissipation de puissance importante.

Une analogie utile consiste à considérer l’appareil comme une baignoire qui se remplit d’eau. Les termes liés à la production de chaleur représentent le robinet : ils injectent de l’énergie dans le système. Les conditions aux limites thermiques représentent le trou d’évacuation :

• Si la vanne est complètement ouverte (ce qui permet une excellente dissipation de la chaleur), le niveau d’eau reste bas.
• Si la buse est partiellement obstruée (ce qui limite l’évacuation de la chaleur), le niveau de l’eau monte.
• Si le bouchon est fermé (et qu’il n’y a donc pas d’évacuation de la chaleur), la baignoire finira par déborder.

Dans cette analogie, le niveau de l’eau représente la température du réseau. Par conséquent, l’objectif des conditions limites thermiques n’est pas purement esthétique ou secondaire : elles décrivent l’environnement physique qui détermine en fin de compte si le dispositif fonctionne à une température basse, modérée ou excessive, pouvant entraîner une défaillance.

2.3 Points de discrétisation de la température (pourquoi la température est discrétisée)

L’éditeur de maillage thermique affiche également une plage de températures et un certain nombre de points de température. Il ne s’agit pas du maillage thermique spatial, mais de la grille de température utilisée pour les précalculs. De nombreuses parties du modèle électrique/thermique nécessitent l’évaluation répétée de grandeurs dépendant de la température pendant la résolution couplée, et il serait coûteux en termes de calcul de les recalculer à partir de zéro à chaque température intermédiaire. Par exemple, en fonction des statistiques choisies et du modèle de densité d’états, le solveur peut avoir besoin de relations dépendant de la température pour les densités de porteurs, les fonctions quasi-thermodynamiques et les grandeurs associées. Au lieu d’évaluer ces grandeurs « en temps réel » à chaque mise à jour de la boucle couplée \(T_L\), OghmaNano précalcule des tableaux de référence à un nombre fini de points de température, puis effectue une interpolation entre ces points pendant l’exécution. Dans cet exemple, les tableaux sont générés à 7 points entre 290 K et 350 K.

Cela fait de la plage de températures un choix à effectuer lors de la modélisation : elle doit englober de manière confortable les températures prévues pendant l’auto-échauffement. Si l’appareil chauffe au-delà de la température maximale, vous risquez d’obtenir des artefacts d’extrapolation (ou une saturation, selon les paramètres), ce qui n’est généralement pas souhaitable. En règle générale, choisissez une plage qui couvre les températures de fonctionnement prévues, en laissant une marge.

3. Paramètres thermiques (conductivité thermique et temps de relaxation)

Dans la simulation électrothermique, les conditions aux limites définissent la manière dont la chaleur peut s’échapper de l’appareil, tandis que les paramètres des matériaux définissent la manière dont la chaleur se propage à l’intérieur de la structure. En pratique, l’un des paramètres les plus importants qui contrôlent l’augmentation de la température et la formation de points chauds est la conductivité thermique \(k\) (parfois notée \(\kappa\)).

Pour afficher ou modifier ces paramètres, ouvrez l’onglet « Thermique » et cliquez sur « Paramètres thermiques » (le bouton marqué \(k\) / \(\kappa\)). Cela ouvre l’éditeur de paramètres thermiques, qui affiche les paramètres par couche, comme illustré dans ??.

L’éditeur affiche trois paramètres essentiels pour chaque calque :

• La conductivité thermique \(k\) détermine l’efficacité avec laquelle la chaleur se diffuse à travers cette couche. Les couches à faible \(k\) agissent comme des obstacles thermiques et peuvent entraîner une augmentation importante de l’élévation de température du réseau prévue.
• Temps de relaxation des électrons et temps de relaxation des trous : ces paramètres ne sont nécessaires que lorsque l’on utilise le modèle hydrodynamique/de transport d’énergie, dans lequel les températures des porteurs peuvent différer de la température du réseau cristallin. Dans le modèle standard de la chaleur du réseau cristallin, ils ne sont pas nécessaires, car les températures des porteurs sont liées à \(T_L\).

Dans la plupart des processus de travail, les conductivités thermiques sont sélectionnées parmi les valeurs publiées (ou mesurées lorsque cela est possible), puis affinées uniquement si une justification expérimentale claire est disponible. En association avec les conditions aux limites, \(k\) est l’un des principaux facteurs déterminant la température absolue du dispositif et la structure spatiale des points chauds.

4. Comment fonctionne le solveur électro-thermique couplé (et pourquoi il est plus lent).

Une simulation électrothermique dans OghmaNano est constituée de deux solveurs couplés :

• un solveur électrique entièrement couplé (modèles de dérive-diffusion et de Poisson, incluant la recombinaison et les pièges lorsque ces fonctions sont activées), et
• un solveur thermique qui résout l’équation de diffusion thermique pour le champ de température du réseau \(T_L\) en utilisant des termes de génération de chaleur calculés à partir de la solution électrique.

Pour une tension d’alimentation donnée, le solveur n’exécute pas ces étapes de manière séquentielle. Il effectue plutôt une boucle externe auto-cohérente :

1. résoudre le problème du transport électrique en utilisant le champ de température actuel.
2. calculer les termes liés à la production de chaleur (transport, recombinaison, pertes parasites et, éventuellement, effets optiques),
3. résoudre le problème de diffusion thermique afin de mettre à jour \(T_L\), et
4. Répétez l’opération jusqu’à ce que les résidus électriques et thermiques convergent.

Cette approche itérative est la raison pour laquelle les simulations électrothermiques sont plus lentes que les simulations purement électriques : pour chaque point de fonctionnement, la résolution électrique de Newton et la résolution thermique peuvent être exécutées plusieurs fois avant que le système couplé ne se stabilise. Le résultat est une solution multi-physique véritablement cohérente, plutôt qu’un simple calcul électrique auquel on ajoute la température comme étape de post-traitement.

Le couplage est essentiel, car les problèmes électriques et thermiques se manifestent généralement à des échelles de longueur physiques très différentes. Le transport électrique est principalement régi par la région active, où les champs, les densités de porteurs et les taux de recombinaison varient rapidement. Le transport thermique doit prendre en compte l’ensemble du trajet du flux de chaleur, y compris les contacts, les substrats, l’encapsulation et les dissipateurs thermiques. C’est pourquoi OghmaNano propose des contrôles indépendants pour le maillage thermique et les conditions aux limites thermiques.

Il est également important de reconnaître que le modèle électrothermique implique naturellement plusieurs températures. En général, le modèle fait la distinction entre la température du réseau \(T_l\), qui décrit la température du réseau atomique, la température des électrons \(T_e\) et la température des trous \(T_h\).

Les électrons et les trous sont des populations de porteurs de charge mobiles qui, en théorie, peuvent avoir leurs propres énergies locales, tout comme ils ont leurs propres niveaux quasi-Fermi. Dans le modèle électro-thermique standard utilisé ici, les températures des électrons et des trous sont liées à la température du réseau, de sorte que tous les sous-systèmes partagent un champ de température unique.

Dans les conditions extrêmes, OghmaNano prend également en charge un modèle de transport hydrodynamique dans lequel les températures des électrons et des trous peuvent différer de la température du réseau cristallin. Ce modèle avancé est décrit ailleurs et n’est requis que dans des situations spécifiques. Pour la majorité des simulations de dispositifs, y compris celle présentée dans ce tutoriel, le modèle électrothermique basé sur le réseau cristallin est à la fois approprié et suffisant.

5. Auto-échauffement et comparaison avec les résultats expérimentaux.

Dans cet exemple, le dispositif a été comparé à des données expérimentales à l’aide d’un modèle électrothermique. Aucun ajustement n’est effectué dans ce tutoriel. Au lieu de cela, la comparaison existante est utilisée pour illustrer comment l’auto-échauffement modifie la réponse électrique du dispositif et pourquoi un traitement électrothermique combiné est nécessaire pour des courants de fonctionnement plus élevés.

L’outil de comparaison est accessible via la barre d’outils d’automatisation, comme illustré dans ??. L’ouverture de l’outil d’ajustement aux données expérimentales affiche les données expérimentales de la caractéristique courant-tension de cet appareil, comme illustré dans ??. Cette courbe représente le comportement mesuré de l’appareil sous polarisation.

Effectuer une itération. Après une seule simulation électrothermique, la courbe JV simulée est superposée aux données expérimentales, comme indiqué dans ??. Dans ce cas, la concordance est étroite, ce qui indique que le modèle prend en compte les principaux phénomènes de transport et les aspects thermiques du dispositif.

On comprend mieux le rôle de l’auto-échauffement en répétant la même comparaison, mais en désactivant le modèle thermique :

1. Retournez à la bande thermique.
2. Désactivez le modèle thermique.
3. Effectuez une nouvelle itération de la comparaison.

La courbe I-V obtenue ne correspond plus aux données expérimentales, comme le montre la figure ??. La divergence est la plus marquée aux tensions de polarisation élevées, où les densités de courant sont importantes et l’autochauffement dû au transport des porteurs devient significatif.

Cette comparaison illustre un résultat essentiel de la modélisation électrothermique : une fois que les dispositifs fonctionnent dans des conditions de forte densité de courant, la dissipation de puissance électrique a un impact sur le transport en augmentant la température. Les simulations électriques à température fixe ne peuvent pas prendre en compte cet effet, contrairement au modèle électrothermique couplé.