AccueilChimie & MatériauxCroissance cristalline — Solidification dendritique

❄️ Croissance cristalline — Solidification dendritique

Simulation par champ de phase de la solidification cristalline. Observez des bras dendritiques croître avec une symétrie d'ordre 6, contrôlez la surfusion, la force d'anisotropie et la diffusivité thermique pour explorer la physique de la nucléation et de la croissance cristalline.

Chimie & Matériaux3DAvancé60 FPS
crystal-growth-3d ↗ Ouvrir en autonome

À propos de la croissance cristalline & de la solidification dendritique

Quand un métal liquide ou de l'eau refroidit sous son point de fusion, la solidification commence par nucléation — la formation spontanée de minuscules embryons cristallins. La théorie classique de la nucléation prédit qu'un noyau doit dépasser un rayon critique r* = 2γTeq/(ΔHf·ΔT) avant de pouvoir croître de manière stable, où γ est l'énergie de surface solide-liquide, Teq est le point de fusion d'équilibre, ΔHf est la chaleur latente de fusion, et ΔT est la surfusion. Ce simulateur place un cristal germe au centre et résout l'équation de champ de phase ∂φ/∂t = M[∇²φ + φ(1−φ)(φ−1/2+λu)] sur une grille en différences finies, où φ est le paramètre d'ordre du champ de phase (0 = liquide, 1 = solide), M est la mobilité de l'interface, λ est la constante de couplage thermique, et u = T − Teq est la température sans dimension.

L'équation de diffusion thermique couplée ∂u/∂t = D·∇²u + (1/2)∂φ/∂t capture la libération de chaleur latente à mesure que le front solide avance. L'anisotropie cristalline est encodée en modulant l'énergie de l'interface avec a(θ) = 1 − ε·cos(n·θ), où θ est l'angle local de la normale à l'interface, ε est la force d'anisotropie, et n est le nombre de plis (n=6 pour la glace hexagonale, n=4 pour les métaux cubiques). Cela produit les bras dendritiques caractéristiques — protubérances allongées le long des directions cristallographiques énergétiquement préférées — qui donnent aux flocons de neige leur symétrie d'ordre 6 et aux métaux coulés leur microstructure ramifiée.

Questions fréquemment posées

Qu'est-ce que le modèle de champ de phase utilisé dans cette simulation ?

La simulation résout ∂φ/∂t = M[∇²φ + φ(1−φ)(φ−1/2+λu)], où φ est le paramètre d'ordre (0=liquide, 1=solide), M est la mobilité de l'interface, λ couple le champ thermique u=T−Teq à la phase, et les termes à double puits φ(1−φ)(φ−1/2) stabilisent les deux états massifs. L'équation couplée ∂u/∂t = D·∇²u + (1/2)∂φ/∂t transporte la chaleur latente libérée pendant la solidification.

Pourquoi les cristaux poussent-ils selon un motif dendritique (en arbre) ?

La croissance dendritique naît d'une instabilité morphologique : une protubérance sur l'interface solide-liquide croît plus vite car elle libère la chaleur latente plus efficacement dans le liquide surfondu. La compétition entre la tension de surface (qui stabilise les interfaces plates) et la diffusion (qui les déstabilise) sélectionne un rayon de pointe de dendrite caractéristique et une vitesse de croissance connue sous le nom de solution d'Ivantsov.

Qu'est-ce que la surfusion et pourquoi affecte-t-elle la forme du cristal ?

La surfusion (ΔT) est de combien la température du liquide est en dessous du point de fusion d'équilibre Teq. Une faible surfusion produit des cristaux compacts à croissance lente ; une forte surfusion entraîne une croissance dendritique rapide avec des rayons de pointe fins. C'est pourquoi les métaux trempés rapidement ont des microstructures dendritiques fines tandis que les liquides refroidis lentement forment des grains équiaxes grossiers.

Que contrôle la force d'anisotropie ε ?

La force d'anisotropie ε contrôle à quel point l'énergie et la mobilité de l'interface varient avec l'orientation cristallographique. Un ε plus élevé accentue les directions de croissance préférées (six bras d'un flocon de neige) et supprime la croissance hors axe. À ε=0 l'interface est isotrope et le cristal croît en cercle ; augmenter ε vers 0,07 produit des bras dendritiques nets et bien définis avec un ramification latérale étroite.

Comment la symétrie cristalline d'ordre 6 est-elle mise en œuvre ?

La fonction d'anisotropie a(θ) = 1 − ε·cos(n·θ) est appliquée à l'angle du gradient de l'interface θ = atan2(∂φ/∂y, ∂φ/∂x). Fixer n=6 crée six directions énergétiquement préférées séparées de 60°, correspondant à la symétrie cristalline hexagonale (glace Ih). Choisir n=4 donne une symétrie cubique (d'ordre 4) correspondant à des métaux comme le fer ou l'aluminium ; n=3 donne trigonal, n=8 des motifs quasi-cristallins.

Qu'est-ce que la nucléation et comment démarre-t-elle la croissance cristalline ?

La nucléation est la formation d'un minuscule embryon cristallin dans le liquide. La théorie classique de la nucléation (CNT) montre qu'un noyau doit dépasser un rayon critique r* = 2γTeq/(ΔHf·ΔT) pour être stable ; les embryons plus petits se redissolvent dans le liquide. Dans cette simulation, une petite graine solide est placée au centre pour contourner la barrière de nucléation et observer directement la croissance dendritique anisotrope.

Qu'est-ce qui détermine la diffusivité thermique dans ce modèle ?

La diffusivité thermique D contrôle la vitesse à laquelle la chaleur latente diffuse loin de l'interface en croissance via ∂u/∂t = D·∇²u + (1/2)∂φ/∂t. Un D plus élevé évacue la chaleur plus vite, réduisant le halo thermique autour de la pointe de la dendrite et permettant une vitesse de pointe plus élevée. Un D plus faible crée une couche limite thermique plus grande qui ralentit la croissance et peut déclencher un fractionnement de pointe ou des morphologies en algue.

Pourquoi les flocons de neige ont-ils une symétrie d'ordre 6 mais sont-ils tous différents ?

La symétrie d'ordre 6 des flocons de neige reflète la structure cristalline hexagonale de la glace Ih où les molécules d'eau forment des anneaux liés par hydrogène avec des angles de liaison de 120°. Les six bras croissent simultanément dans le même environnement local, si bien qu'ils se développent de manière identique. Cependant, chaque bras trace une histoire de croissance unique alors que le flocon de neige culbute dans l'air avec une température et une humidité changeantes — rendant chaque flocon de neige unique dans le détail tout en étant symétrique en lui-même.

Quelles sont les applications réelles de la recherche sur la solidification dendritique ?

La microstructure dendritique affecte directement les propriétés mécaniques des métaux coulés, des alliages d'aubes de turbine, et des pièces imprimées en 3D. Contrôler l'espacement des bras de dendrite (DAS) via le taux de refroidissement optimise la résistance et la ductilité. Des modèles de champ de phase comme celui-ci sont utilisés industriellement pour concevoir des procédés de coulée, prédire la microségrégation dans les superalliages de nickel, et optimiser les paramètres de fabrication additive pour les composants aérospatiaux.

Comment se produit la compétition de grains dans les métaux en solidification ?

Quand plusieurs noyaux se forment simultanément, chacun fait croître des dendrites alignées avec ses propres axes cristallographiques. Là où deux champs de dendrites se rencontrent, le grain à croissance plus rapide (avec des axes préférés à peu près alignés avec le gradient de température) bloque et consomme le voisin à croissance plus lente. Cette compétition de grains détermine la microstructure finale — grains colonnaires dans les alliages solidifiés directionnellement, grains équiaxes dans les matériaux coulés rapidement.

⚙ Sous le capot

Un modèle de champ de phase de solidification cristalline dendritique avec surfusion et anisotropie réglables.

3DPhase-Field ModelDendritic GrowthCrystal SolidificationAnisotropy

3D · Three.js / WebGL renderer · 60 FPS target · runs fully client-side, no install

Qu'avez-vous trouvé ?

Ajouter les étapes de reproduction (facultatif)