Comment un motif de fentes gravées au laser dans un film Ginestium™ apprend à la chaleur à filer dans une direction et à rester bloquée dans l'autre — et comment lire, sans être thermicien, ce que l'outil calcule pour vous.
Sur une puce EFFIBLUE, on veut souvent que la chaleur d'un point chaud s'évacue vite dans une direction (vers un radiateur, un bord de boîtier) sans se répandre dans l'autre (vers un composant sensible voisin). Un film métallique plein conduit pareil dans toutes les directions — inutile pour ce genre de tri directionnel.
La solution : découper le film Ginestium™ avec des fentes très fines au laser, selon un motif régulier. Bien choisi, ce motif laisse la chaleur filer presque librement dans un sens et la freine fortement dans l'autre. Le simulateur calcule, pour un motif et des réglages donnés, à quel point ce tri fonctionne — avant de graver quoi que ce soit pour de vrai.
L'outil répond à une question concrète : « si je grave ce motif de fentes, avec cette épaisseur de film et sur ce substrat, quelle sera la conduction dans le sens que je veux favoriser (x), et dans celui que je veux bloquer (y) ? » Il donne aussi deux indicateurs de bon sens pratique : le design tiendra-t-il mécaniquement, et les dilatations différentielles film/substrat risquent-elles de le fatiguer au fil des cycles thermiques.
C'est un outil d'exploration avant fabrication : on y règle un motif, on regarde le résultat, on ajuste, et une fois satisfait, on exporte le dessin au format utilisé par la machine de gravure laser.
Imaginez le film vu du dessus, avec des bandes de matière pleine séparées par de fines fentes régulières (remplies d'air, de vide, ou d'un autre isolant). Dans le sens des bandes, la chaleur circule presque sans obstacle : c'est le sens x, celui qu'on veut favoriser. Pour aller dans l'autre sens, elle doit systématiquement traverser une fente — un isolant — ce qui la freine fortement : c'est le sens y, celui qu'on veut bloquer.
Le résultat de ce calcul, ce sont deux nombres : k_xx (conductivité effective dans le sens favorisé) et k_yy (conductivité effective dans le sens bloqué). Leur rapport k_xx / k_yy est l'anisotropie obtenue — plus il est grand, mieux le tri fonctionne.
Le film ne flotte pas dans le vide : il repose sur un substrat EFFIBLUE nettement plus épais. La chaleur peut donc partir dans deux directions différentes, et l'outil les traite avec deux modèles distincts.
Concrètement : k_xx système et k_yy système tiennent compte de cette dilution par le substrat, tandis que la R_th hors-plan répond à « si la chaleur doit traverser l'empilement verticalement, quelle résistance totale rencontre-t-elle ? ».
Dans un film très mince, la chaleur ne circule plus tout à fait comme dans un bloc massif : les porteurs de chaleur (phonons) sont gênés par la proximité des bords du film. Plus le film est fin par rapport à leur « libre parcours moyen » (une distance caractéristique propre à chaque matériau), plus sa conductivité réelle chute sous la valeur du matériau massif. La case « Corriger k selon le confinement » applique cette correction si on la coche — elle est désactivée par défaut, à activer volontairement pour les films très fins.
Un modèle n'est utile que si l'on sait ce qu'il suppose. Voici les hypothèses qui sous-tendent chaque calcul.
Le calcul suppose un motif de fentes qui se répète à l'identique à perte de vue (une seule « cellule » est résolue puis dupliquée mentalement). En bord réel de puce, le motif s'arrête : l'outil ne modélise pas ces effets de bord.
k_xx et k_yy sont des conductivités en régime établi (le système a eu le temps de se stabiliser). Ce simulateur ne dit rien de la vitesse à laquelle la chaleur se propage dans le temps — c'est le rôle du comparateur de diffusion thermique, un outil compagnon. Conséquence directe, détaillée au nota bene (§10) : la capacité calorifique des matériaux n'intervient nulle part dans ce calcul.
Les ponts de jonction mesurent souvent quelques micromètres — trop fins pour la grille de calcul. Leur effet de court-circuit thermique local est ajouté par une formule analytique calibrée (validée à moins de 5 % d'écart contre un solveur haute résolution), pas résolu maille par maille comme le reste du motif.
La correction de confinement suppose un transport de chaleur « diffusif » — valable pour un film avec un intérieur suffisamment épais. En dessous de 5 nm, cette hypothèse devient discutable (le transport devient plutôt bidimensionnel/balistique) ; l'outil ne désactive pas la correction en silence, il signale qu'on sort du domaine de validité.
La « fragilité mécanique » et le « risque CTE » sont des scores qualitatifs simples (fraction de matière évidée, densité de ponts, écart de dilatation thermique × amplitude de cyclage), pensés pour repérer vite une configuration à risque — pas le résultat d'une simulation mécanique par éléments finis ni d'un essai de fatigue réel.
La résistance d'interface (R_th interface) est un nombre unique fourni par l'utilisateur, pas un modèle physique de contact — pratique pour explorer sa sensibilité, mais à recaler sur une mesure réelle avant conclusion définitive. Elle n'agit que sur le chemin vertical : voir le nota bene (§10).
L'outil ne masque jamais un cas limite : il calcule quand même, mais prévient.
k_yy devient peu fiable. Augmenter la résolution de grille ou élargir le pas.Paramètres (JSON) — exporte tous les réglages de l'outil (matériau, géométrie, substrat, ponts…), pas les résultats de calcul. Utile pour reprendre une session plus tard ou la transmettre à un collègue, qui la réimporte à l'identique.
Gerber (.gbr) — exporte le dessin vectoriel réel du motif de gravure, dans le format attendu par la machine de gravure laser. C'est l'aboutissement pratique de la démarche : une fois le motif validé virtuellement, ce fichier passe directement en fabrication.
Cette section est la référence technique du guide : elle décrit les formules réellement exécutées par le moteur serveur, dans l'ordre où elles s'enchaînent, et justifie chaque choix de modèle. Elle n'est pas nécessaire pour se servir de l'outil.
Si la correction d'effet de taille est cochée, la conductivité saisie est d'abord rabattue :
C'est la règle de Matthiessen appliquée au transport de chaleur. Deux mécanismes gênent indépendamment les phonons : les collisions internes au matériau (libre parcours Λ) et les rebonds sur les parois de la structure (parcours limité à L_c). Les taux de collision s'ajoutent, donc les inverses des longueurs s'ajoutent : 1/Λ_eff = 1/Λ + 1/L_c, soit Λ_eff = Λ·L_c/(Λ + L_c). Comme la conductivité est proportionnelle au libre parcours moyen, k_eff/k = Λ_eff/Λ = L_c/(L_c + Λ).
Le minimum dans L_c dit lequel des deux confinements mord : l'épaisseur du film, ou la paroi de matière restée entre deux fentes. Sur ce simulateur, contrairement au balayage paramétrique, les deux sont plausibles — c'est pourquoi le minimum est pris explicitement.
Le même facteur est appliqué à la conductivité hors-plan k_z : le confinement supprime le transport dans les deux directions.
Le rapport k_x/k_y vaut bien A, ce qui est demandé. Mais surtout, le produit k_x·k_y reste égal à k_eff² : la moyenne géométrique est conservée. Tourner le bouton d'anisotropie ne change donc pas subrepticement la quantité totale de chaleur que le matériau sait porter, seulement sa répartition entre les deux directions. Avec A = 1, on retrouve exactement k_eff dans les deux sens.
Une gravure à profondeur partielle laisse de la matière sous la fente. Sur l'épaisseur du film, la fraction p du haut est remplie par le milieu de fente et la fraction (1 − p) du bas est encore du film plein. Vu du transport latéral, ces deux couches conduisent côte à côte — en parallèle — d'où la moyenne arithmétique pondérée par l'épaisseur, et non la moyenne harmonique. À gravure traversante (p = 1), on retrouve k_gap.
Le motif est réduit à une cellule élémentaire qui se répète : largeur = période de ponts (ou pas × rapport d'aspect pour l'hexagonal), hauteur = 2 × pas. Chaque maille reçoit k_fente si elle tombe dans une fente, sinon k_x et k_y. On résout ensuite un problème de conduction pure sur cette cellule, deux fois : une fois dans le sens x, une fois dans le sens y.
Entre deux mailles voisines, la chaleur traverse la demi-maille de l'une puis la demi-maille de l'autre : ce sont deux résistances en série, et la conductivité équivalente de deux moitiés égales est leur moyenne harmonique. Prendre la moyenne arithmétique laisserait la chaleur « fuir » à travers une frontière matière/fente et surestimerait systématiquement k_yy — c'est-à-dire précisément la grandeur que l'outil est censé rendre petite. Autre propriété utile : la moyenne harmonique s'annule dès qu'une des deux conductivités est nulle, ce qui fait d'une maille sous vide poussé une vraie barrière.
(5) est le bilan de flux en volumes finis : à l'équilibre, ce qui entre dans une maille par ses quatre faces en ressort, donc sa température est la moyenne de ses voisines pondérée par les conductances de face. Balayer toutes les mailles en réutilisant les valeurs déjà mises à jour, c'est la méthode de Gauss-Seidel ; elle converge, mais lentement sur une grille fine. Le facteur ω > 1 (sur-relaxation successive) extrapole au-delà de la correction calculée et réduit le nombre d'itérations d'environ un ordre de grandeur. ω doit rester dans ]1 ; 2[ pour que le schéma reste stable ; 1,9 est la valeur retenue.
Conditions aux limites : température imposée à 1 et 0 sur les deux faces perpendiculaires au sens testé, périodicité sur les deux autres — c'est la traduction numérique de l'hypothèse « motif infini ». Le flux total est mesuré toutes les 50 itérations et le calcul s'arrête dès que sa variation relative passe sous 2×10⁻³, avec un plafond de 6 000 itérations.
Comme l'écart de température imposé vaut exactement 1 sur toute la longueur de la cellule, ce flux est la conductivité effective, sans normalisation supplémentaire. Le maillage est dimensionné pour donner plusieurs mailles en travers de la fente (44 à 220 mailles en hauteur, 24 à 220 en largeur) et plafonné à 46 000 mailles au total ; au-delà, la grille est réduite homothétiquement, ce qui déclenche l'avertissement « fente sous-résolue ».
Pour un empilement lamellaire infini de bandes alternées, ce ne sont pas des approximations mais les solutions exactes : le long des bandes les conductances s'ajoutent, en travers ce sont les résistances. Ce sont aussi les bornes de Wiener — aucune microstructure de même fraction évidée ne peut dépasser k_xx^idéal ni descendre sous k_yy^idéal. D'où leur double usage : affichage immédiat (aucun solveur à lancer) et contrôle de cohérence du résultat numérique, qui doit toujours tomber entre les deux.
Le solveur de l'étape 4 tourne sans les ponts. Leur effet est réintroduit ensuite, analytiquement :
Un pont fait quelques micromètres, sur une grille dont les mailles en font plusieurs dizaines. Le résoudre directement demanderait de raffiner la grille partout — un coût d'environ ×100 en nombre de mailles pour un détail localisé. La correction (8) fait mieux avec moins : elle récupère la conductivité de la bande de fente telle que le solveur l'a effectivement calculée, y ajoute le pont comme un court-circuit en parallèle dans cette bande seule, puis réassemble. Le facteur 0,88 traduit le décalage en quinconce : des ponts alternés ne s'alignent pas en un chemin continu et court-circuitent donc un peu moins. L'ensemble est calibré à moins de 5 % d'écart contre un solveur haute résolution.
k_xx n'est pas touché : un pont ajoute de la matière en travers des bandes, il ne change rien à la conduction le long des bandes.
Film et substrat sont deux nappes accolées : latéralement elles conduisent en parallèle, chacune apportant sa conductance k·t. Diviser la somme par l'épaisseur totale reconvertit le tout en conductivité équivalente, directement comparable à celle du film seul. C'est aussi l'explication mécanique de la dilution : avec t_f de quelques centaines de nanomètres contre t_s de l'ordre du millimètre, le terme du substrat pèse des milliers de fois plus, et l'anisotropie du film se noie dans la moyenne. Le rainurage du substrat existe pour rendre le substrat lui-même anisotrope — seule manière de faire remonter K_xx/K_yy.
Verticalement, la chaleur traverse les trois couches l'une après l'autre : les résistances surfaciques s'ajoutent, sans pondération. C'est le seul endroit du modèle où la résistance d'interface intervient — voir le nota bene (§10). Noter que k_z est bien la conductivité hors-plan du matériau, indépendante de celle dans le plan : pour un matériau lamellaire comme le graphène, l'écart atteint un facteur 100, et utiliser la valeur dans le plan rendrait R_z absurdement optimiste.
(12) a un sens physique clair : c'est la déformation différentielle accumulée entre film et substrat sur un cycle thermique. 2 500 ppm, c'est 0,25 % d'écart de dilatation — l'ordre de grandeur où la délamination de films minces est documentée. Le score est donc grossier mais ancré.
(11) n'a pas cette prétention : c'est un classement empirique qui oppose la matière retirée (φ, qui fragilise) à la densité de ponts (b, qui consolide), avec deux majorations pour les géométries plus tenaces. Les coefficients 2 et 3 n'ont pas d'unité et ne sortent d'aucun modèle mécanique. À traiter comme un feu tricolore, jamais comme un critère de dimensionnement.
Deux questions reviennent régulièrement au sujet des matériaux déposés : où intervient leur capacité calorifique, et que fait-on des résistances thermiques d'interface. Réponses en clair, avec les hypothèses effectivement retenues dans les trois outils thermiques.
C'est le régime qui décide. En régime permanent, la capacité calorifique n'a aucun effet : une fois le champ de température stabilisé, seule la conductivité fixe la répartition de la chaleur. En régime transitoire, au contraire, c'est elle qui, avec la conductivité, fixe la diffusivité thermique α = k/(ρ·c_p) et donc la vitesse d'évolution de la température.
Les trois outils se répartissent donc ainsi :
| Outil | Régime | ρ·c_p du film |
|---|---|---|
| Simulateur d'anisotropie thermique | permanent | n'intervient pas |
| Balayage paramétrique | permanent | n'intervient pas |
| Comparateur de diffusion thermique | transitoire | intervient explicitement |
Pour cet outil-ci, ce n'est donc ni un oubli, ni une intégration implicite : l'équation résolue à l'étape 4 est ∇·(k∇T) = 0, où ρ·c_p n'apparaît pas — et le simulateur ne porte d'ailleurs pas ces champs dans ses paramètres. Le catalogue matériaux commun aux trois outils renseigne bien ρ et c_p pour chaque matériau, mais seul le comparateur de diffusion les consomme.
Dans le comparateur de diffusion, la capacité entre par la diffusivité effective du système film + substrat :
Deux hypothèses explicites y sont faites sur les couches déposées. D'abord, la fraction ablatée du film ne stocke rien : le milieu de fente (air, vide, aérogel) est considéré sans capacité utile à ces échelles, d'où le facteur (1 − ψ). Ensuite, ρ et c_p sont pris à 300 K et constants — aucune dépendance en température n'est modélisée.
Les valeurs proviennent du catalogue serveur partagé :
| Matériau | ρ (kg/m³) | c_p (J/kg·K) | Statut de la donnée |
|---|---|---|---|
| Ginestium™ | 2 100 | 710 | substitut — carbone pyrolytique synthétique orienté à 300 K ; pas une mesure Ginestium™ |
| Graphène CVD few-layer | 2 200 | 700 | littérature |
| Diamant CVD nanocristallin | 3 400 | 536 | mesuré sur film NCD à 302 K |
| Cuivre (témoin) | 8 960 | 385 | littérature, cuivre massif |
| Substrat | saisi par l'utilisateur (défauts 3 900 / 880) | — | |
Un matériau « personnalisé » n'a ni ρ ni c_p : dans ce cas le comparateur renvoie explicitement « pas de diffusivité » plutôt que de deviner une valeur.
La pondération par les épaisseurs écrase la contribution du film. Avec les réglages par défaut du comparateur — Ginestium™ 400 nm sur 10 µm de substrat — le film ne pèse que ≈ 1,6 % de la capacité du système ; portez le substrat à 2 mm et sa part tombe sous 0,01 %. La diffusivité effective est donc, en pratique, « la conductivité latérale du système divisée par la capacité du substrat ».
Conséquence directe : le front de diffusion est très peu sensible au c_p du film. Non pas parce qu'on le néglige — il est bien dans la formule — mais parce que la géométrie de l'empilement lui laisse peu de poids. Le film pèse par son k·t, au numérateur, bien plus que par sa capacité au dénominateur.
Elles ne sont pas ignorées, mais leur prise en compte est partielle et hors-plan uniquement.
Ce qui existe. Le paramètre R_th interface de ce simulateur (0 à 100 ×10⁻⁶ m²·K/W, 10 par défaut ; même curseur dans le comparateur de diffusion) entre en série dans la résistance hors-plan — c'est l'équation (10) du §09, rappelée ici :
Il s'agit d'une seule interface globale film / substrat. Le modèle géométrique ne comporte que deux couches : il n'y a ni empilement multicouche avec une résistance par interface, ni résistance latérale entre grains ou entre motifs.
Un ordre de grandeur utile : pour le graphène (k_z = 6 W/m·K sur 1,675 nm), le terme de film vaut ≈ 3×10⁻¹⁰ m²·K/W, cinq ordres de grandeur sous la résistance d'interface typique. Dans R_z, ce sont toujours l'interface et le substrat qui dominent, jamais le film — c'est vrai aussi du Ginestium™, dont le k_z élevé (1 000 W/m·K) rend le terme de film négligeable.
Ce qui est négligé, et c'est une hypothèse assumée. R_th,if n'entre ni dans l'homogénéisation k_xx/k_yy de l'étape 4, ni dans K_xx/K_yy de l'étape 7, ni dans la diffusivité α_eff du front de diffusion du comparateur, qui l'affiche d'ailleurs explicitement sous le curseur concerné. Le balayage paramétrique, lui, ne porte aucune résistance d'interface : c'est un outil purement dans le plan du film.
Pour un transport latéral, la résistance d'interface n'est pas une résistance en série sur le chemin de la chaleur : c'est un couplage entre deux chemins parallèles, le film et le substrat. La négliger revient exactement à supposer un couplage parfait selon z — film et substrat à la même température locale à chaque abscisse, ce qu'on appelle l'hypothèse de plaque équivalente. Elle est déjà implicite dans l'équation (9), qui attribue une seule conductivité au couple film + substrat.
Le régime où elle casse est identifiable : temps courts, ou film très conducteur sur substrat peu conducteur. Le film y prendrait de l'avance sur le substrat au lieu de le tirer avec lui — un effet que le modèle ne peut pas voir. La sensibilité à R_th interface se teste ici même, mais elle ne se lira que sur la R_th hors-plan.
Pas besoin de manier ces unités pour utiliser l'outil : les curseurs et champs numériques font l'essentiel du travail. Ce glossaire sert juste à décoder les résultats affichés.