Multiphysics simulations of microstructure influence on hysteresis and eddy current losses of electrical steel

Cet article utilise des simulations micromagnétiques et l'homogénéisation computationnelle sur des microstructures d'acier Fe-Si numérisées, imprimées par jet de liant, pour démontrer que l'optimisation de la taille des grains et de l'épaisseur de la phase aux joints de grains peut réduire efficacement les pertes par hystérésis et les pertes par courants de Foucault dans l'acier électrique.

L'essentiel

Imaginez que vous essayez de construire un moteur électrique super efficace. Le cœur de ce moteur est un métal spécial appelé acier électrique. Lorsque l'électricité traverse le moteur, ce métal agit comme un contrôleur de trafic pour les champs magnétiques. Cependant, tout comme une autoroute très fréquentée, le métal n'est pas parfait. À mesure que les champs magnétiques s'inversent, le métal « se fatigue » et perd de l'énergie sous forme de chaleur. C'est ce qu'on appelle la perte d'énergie, et cela rend votre moteur moins efficace.

Pendant longtemps, les scientifiques ont essayé d'améliorer ce métal en modifiant sa recette chimique. Mais récemment, une nouvelle façon de fabriquer du métal appelée Fabrication Additive (en gros, l'impression 3D de métal) a ouvert une nouvelle porte. Cet article explore ce qui se passe à l'intérieur de ce métal imprimé en 3D au niveau microscopique et comment réduire sa perte d'énergie.

Voici la décomposition de leur découverte, en utilisant des analogies de la vie quotidienne :

1. Les deux ennemis : l'hystérésis et les courants de Foucault

Pour comprendre le problème, imaginez que le métal est une foule géante de petits aimants (appelés domaines magnétiques) à l'intérieur d'un stade.

• Perte par hystérésis (La foule « collante ») : Imaginez que la foule essaie de se retourner pour faire face à une nouvelle direction. Certains sont têtus et restent accrochés à leur ancienne direction, ce qui rend le retournement de tout le groupe difficile. Vous devez pousser très fort (utiliser de l'énergie) pour les faire basculer. Cette « collante » est l'hystérésis. L'article a découvert que la « colle » entre les grains de métal (les joints de grains) agit comme un piège collant. Si les grains sont trop gros, la foule reste bloquée dans des endroits spécifiques, ce qui rend le retournement plus difficile.
• Perte par courants de Foucault (La foule en « court-circuit ») : Maintenant, imaginez que la foule court également sur la piste du stade. Si la piste est une boucle lisse et ouverte, ils peuvent courir vite et facilement. Mais s'il y a des murs ou des barrières, ils doivent courir en cercles ou heurter des obstacles, créant de la friction (chaleur). Dans le métal, ces chemins de course sont les courants électriques. Si le métal est une seule pièce géante et lisse, les courants circulent librement et créent beaucoup de chaleur. Si vous érigez des murs (isolants) entre les grains, les courants sont bloqués et ne peuvent pas courir aussi loin, ce qui réduit la chaleur.

2. L'expérience : Construire des jumeaux numériques

Les chercheurs n'ont pas seulement deviné ; ils ont construit des jumeaux numériques du métal.

• Ils ont pris de vrais échantillons de métal imprimés en 3D (certains avec du Bore, d'autres sans) et ont pris des photos haute puissance (images MEB).
• Ils ont ensuite créé deux types de modèles informatiques :
  1. Le modèle « Idéal » : Ils ont construit des grains parfaits générés par ordinateur, comme un puzzle de mosaïque.
  2. Le modèle « Réel » : Ils ont scanné les photos réelles du métal et les ont transformées en une carte numérique.

Ils ont utilisé ces cartes pour simuler le comportement de la foule magnétique et le mouvement des coureurs électriques.

3. Les grandes découvertes

En lançant des milliers de simulations, ils ont découvert des règles surprenantes pour ajuster le métal :

La taille de grain « Goldilocks » (Ni trop grande, ni trop petite)

• La découverte : Ils ont trouvé que si les grains (les « carreaux » individuels de la mosaïque) font environ 120 micromètres de large, la « collante » (perte par hystérésis) est à son niveau le plus bas.
• Le bémol : Cependant, rendre les grains plus gros fait que les « coureurs » (courants de Foucault) courent plus vite et perdent plus d'énergie.
• L'analogie : Pensez à une piste de danse. Si les carreaux du sol sont trop petits, les danseurs (aimants) trébuchent constamment sur les bords. Si les carreaux sont énormes, les danseurs peuvent tournoyer librement, mais la musique (l'électricité) voyage trop vite et provoque un désordre. Il faut un carreau de taille moyenne pour que la danse soit fluide sans que la musique ne s'emballe.

La stratégie du « Mur Épais »

• La découverte : L'espace entre les grains est rempli d'un matériau spécial (une phase de joint de grain). Les chercheurs ont découvert que rendre ce « mur » plus épais est une situation gagnant-gagnant.
• L'analogie : Imaginez que les grains sont des maisons et que la limite est la clôture.
  • Pour l'hystérésis : Une clôture plus épaisse agit comme une meilleure zone tampon, aidant la « foule » magnétique à changer de direction plus facilement sans rester coincée.
  • Pour les courants de Foucault : Une clôture plus épaisse est une meilleure barrière. Elle empêche les « coureurs » électriques de sauter d'une maison à l'autre. Si la clôture est épaisse et résistante, les coureurs restent coincés dans leurs propres maisons et ne peuvent pas créer une grande boucle génératrice de chaleur.
• Résultat : Des joints plus épais réduisent les deux types de pertes d'énergie.

4. Pourquoi cela importe (selon l'article)

L'article conclut qu'en optimisant simplement la microstructure — spécifiquement en contrôlant la taille des grains et en rendant les limites entre eux plus épaisses — nous pouvons réduire considérablement l'énergie gaspillée dans ces noyaux magnétiques.

Ils ont prouvé qu'il n'est pas nécessaire d'inventer une nouvelle formule chimique ; il suffit d'organiser les atomes existants selon un schéma plus intelligent. Leurs modèles informatiques ont montré que la stratégie du « mur épais » aide le matériau magnétique à changer de direction plus facilement (moins de collante) tout en bloquant simultanément les courants électriques qui causent la chaleur (moins de courts-circuits).

En résumé : Les chercheurs ont utilisé des simulations informatiques pour montrer que l'acier électrique imprimé en 3D fonctionne mieux lorsque les « grains » sont d'une taille moyenne spécifique et que les « clôtures » entre eux sont épaisses. Cette disposition rend le métal moins « collant » pour les aimants et meilleur pour bloquer la chaleur électrique, conduissant à des machines plus efficaces.

Résumé technique

Résumé technique : Influence microstructurale sur l'hystérésis et les pertes par courants de Foucault des aciers électriques fabriqués par fabrication additive

1. Énoncé du problème

L'efficacité des machines électriques, qui dépendent fortement des matériaux magnétiques doux comme les alliages Fe-Si, est fondamentalement limitée par les pertes d'énergie dans le noyau magnétique. Ces pertes se manifestent principalement sous forme de pertes par hystérésis (dominantes à basse fréquence) et de pertes par courants de Foucault. Alors que les aciers électriques conventionnels sont transformés par laminage, une teneur élevée en silicium (>4 % en poids) rend le matériau fragile, ce qui complique la fabrication traditionnelle. La fabrication additive (FA), spécifiquement l'impression par jet de liant (Binder Jet Printing - BJP), offre une nouvelle voie de transformation capable de contourner les limitations du laminage et a été observée comme améliorant les propriétés magnétiques, notamment par la réduction de la coercivité et des pertes, particulièrement lors de l'ajout de Bore. Cependant, les mécanismes sous-jacents reliant les caractéristiques microstructurales spécifiques des aciers Fe-Si-B produits par FA (ex. : taille de grain, morphologie de la phase aux joints de grains) à leurs performances magnétiques et électriques restent obscurs. Ce travail vise à élucider ces mécanismes via une approche de simulation multiphysique.

2. Méthodologie

Les auteurs emploient un cadre multiphysique combinant des simulations micromagnétiques et l'homogénéisation computationnelle pour analyser des microstructures numériques.

2.1 Reconstruction et numérisation de la microstructure

Deux flux de travail distincts ont été développés pour générer des domaines de simulation :

1. Reconstruction synthétique : Basée sur des descripteurs statistiques (taille de grain moyenne $d_G$, épaisseur du joint de grain $\ell_{GB}$) dérivés de la caractérisation expérimentale. Cette méthode utilise l'échantillonnage par disque de Poisson pour la génération de germes, la tessellation de Voronoi pour la formation des grains, et l'intégration d'une phase de joint de grain (GB) définie.
2. Numérisation directe : Des images SEM expérimentales de deux échantillons spécifiques (Échantillon A : 3 % en poids de Si, 0 % de B ; Échantillon B : 5 % en poids de Si, 0,25 % de B) ont été traitées. Cela a impliqué un lissage gaussien, une détection de frontière de phase par seuillage, une binarisation et un étiquetage par composantes connexes pour convertir les images en domaines de simulation étiquetés.

2.2 Simulation de la perte par hystérésis (Landau-Lifshitz)

Le comportement d'hystérésis a été simulé à l'aide de la théorie de Landau-Lifshitz (LL) implémentée dans le package mumax3 via la méthode des différences finies (FDM).

• Équations directrices : L'évolution de l'aimantation $\mathbf{m}(\mathbf{r})$ est régie par l'équation LL, minimisant une fonctionnelle d'énergie libre comprenant les termes d'échange, d'anisotropie, de magnétostatique et de Zeeman.
• Propriétés des matériaux : Les paramètres magnétiques ($K_u$, $M_s$) ont été définis comme des fonctions de la teneur en Silicium ($X_{Si}$). La phase GB a été modélisée avec des propriétés magnétiques considérablement réduites ($K_u$ et $M_s$ mis à l'échelle par 0,001) pour approximer la structure lamellaire mixte Fe-Si/Fe$_2$B.
• Pontage d'échelle : Pour assurer la validité numérique au sein du domaine mésoscopique (échelle micrométrique) tout en maintenant la pertinence physique, l'épaisseur du mur de domaine de Bloch a été artificiellement fixée à $\sim 4$ $\mu$m. Cela a nécessité une constante de gradient $\kappa_m$ environ 5 000 fois plus grande que les valeurs physiques, créant ainsi un « matériau modèle » qui capture les tendances de la dynamique des parois de domaines (nucléation, piégeage, migration) par rapport aux caractéristiques microstructurales.
• Calcul des pertes : La perte par hystérésis ($P_H$) a été dérivée de l'aire des boucles $B-H$ simulées sur plusieurs cycles de désaimantation.

2.3 Simulation des pertes par courants de Foucault (Homogénéisation magnéto-quasi-statique)

Les pertes par courants de Foucault ont été calculées à l'aide d'une approche magnéto-quasi-statique (MQS) basée sur les équations de Maxwell, implémentée via la méthode des éléments finis (FEM) dans le code MFM-FM (framework MOOSE).

• Homogénéisation : Le problème a été réduit à un problème de conduction stationnaire pour déterminer le tenseur de conductivité effective ($\sigma_{eff}$) de la microstructure. Cela a impliqué la résolution des potentiels électriques sous des conditions limites linéaires pour satisfaire la condition de Hill.
• Modélisation de la conductivité : La conductivité de l'intérieur des grains a été dérivée de données expérimentales pour les alliages Fe-Si. La conductivité de la phase GB a été supposée être celle de l'Fe$_2$B massif ($0,018$S m$^{-1}$), nettement inférieure à celle de l'intérieur des grains.
• Calcul des pertes : La perte par courants de Foucault ($P_E$) a été calculée en utilisant la conductivité effective et la fréquence AC appliquée, en considérant l'effet de profondeur de peau. À basse fréquence (jusqu'à 100 Hz), la perte a été approchée par une formulation de limite de basse fréquence.

3. Résultats clés

3.1 Comportement d'hystérésis et pertes

• Rôle de la phase GB : L'introduction d'une phase GB a considérablement modifié le comportement d'hystérésis. Elle a augmenté le champ coercitif ($H_c$) pour les deux échantillons (par exemple, de ~6,8 à ~15,2 kA/m pour l'échantillon A) et a réduit la rémanence ($B_r$). La phase GB agit comme un site de piégeage, restreignant les parois de domaines à des grains uniques et favorisant la nucléation aux frontières.
• Effets de la taille de grain ($d_G$) :
  • Coercivité : Pour 3 % en poids de Si, $H_c$ a montré une légère tendance croissante avec la taille de grain. Pour 5 % en poids de Si, $H_c$ a diminué lorsque la taille de grain augmentait.
  • Rémanence : $B_r$ a augmenté de manière constante avec de plus grandes tailles de grains en raison de la réduction de la fraction volumique de la phase GB à faible aimantation.
  • Perte par hystérésis ($W_H$) : Les pertes ont généralement diminué avec l'augmentation de la taille de grain, atteignant un minimum pour une taille de grain moyenne d'environ 120 $\mu$m.
• Effets de l'épaisseur de la phase GB ($\ell_{GB}$) : L'augmentation de l'épaisseur de la phase GB a conduit à une réduction de la coercivité et de la perte par hystérésis. Ceci est attribué à l'augmentation du volume de la phase GB résistive qui modifie le paysage énergétique magnétique.

3.2 Conductivité et pertes par courants de Foucault

• Conductivité effective ($\sigma_{eff}$) : La phase GB a considérablement réduit la conductivité effective globale du matériau. $\sigma_{eff}$ a diminué à mesure que l'épaisseur de la phase GB augmentait et a augmenté à mesure que la taille de grain augmentait (réduisant la fraction relative de la phase GB résistive).
• Flux de courant : Dans les microstructures présentant des phases GB résistives, les courants électriques étaient fortement déviés, circulant principalement à travers l'intérieur des grains et contournant la phase GB. Cela a créé des « goulots d'étranglement » et des chemins préférentiels.
• Perte par courants de Foucault ($W_E$) :
  • Taille de grain : Les pertes par courants de Foucault ont augmenté avec l'augmentation de la taille de grain, car les grains plus grands offrent moins de barrières résistives (GB) au flux de courant.
  • Épaisseur de la phase GB : L'augmentation de l'épaisseur de la phase GB a réduit les pertes par courants de Foucault en abaissant la conductivité effective du matériau massif.
  • Teneur en Silicium : Les échantillons avec une teneur plus élevée en Silicium (5 % en poids) ont présenté des pertes plus faibles que les échantillons à 3 % en poids, principalement en raison de la faible conductivité intrinsèque des alliages Fe-Si à haute teneur en Si.

4. Signification et affirmations

L'article affirme fournir une approche multiphysique validée pour comprendre les relations microstructure-propriété dans les aciers électriques fabriqués par ajout. Les principales contributions et affirmations incluent :

• Aperçu mécanistique : L'étude démontre que la phase GB n'est pas simplement une frontière passive, mais un composant actif qui influence le piégeage des parois de domaines, la nucléation et les chemins de flux de courant.
• Potentiel d'optimisation : Les résultats suggèrent que l'optimisation de la microstructure peut réduire les pertes dans les noyaux magnétiques. Plus précisément, les auteurs identifient un compromis : alors que des grains plus gros minimisent la perte par hystérésis, ils augmentent la perte par courants de Foucault. Inversement, augmenter l'épaisseur de la phase GB aide à réduire à la fois la perte par hystérésis et la perte par courants de Foucault, bien que cela puisse impacter l'induction de saturation.
• Validation méthodologique : La forte concordance entre les reconstructions synthétiques et les échantillons numérisés par SEM suggère que les modèles purement synthétiques, basés sur des descripteurs statistiques, sont suffisants pour obtenir des informations significatives sur le comportement des matériaux magnétiques doux, réduant potentiellement le besoin d'une caractérisation expérimentale extensive lors des premières étapes de conception.
• Directives de conception : Le travail propose des directives de conception spécifiques, telles que viser une taille de grain moyenne de ~120 $\mu$m pour une perte par hystérésis minimale et utiliser des phases GB plus épaisses pour atténuer les deux mécanismes de perte, à condition que la réduction de l'induction de saturation soit acceptable.

Les auteurs maintiennent un ton modeste, reconnaissant que le matériau modèle utilise des paramètres mis à l'échelle pour la stabilité numérique et que des travaux futurs sont nécessaires pour capturer plus précisément des caractéristiques complexes comme les pores isolés et les phases GB insulaires, ainsi que pour affiner l'homogénéisation des propriétés des phases GB elles-mêmes.