❄️ Crystal Growth — Dendritic Solidification
A phase-field model of dendritic crystal solidification with adjustable undercooling and anisotropy.
Podobne symulacje
O wzroście kryształów i krzepnięciu dendrytycznym
Gdy ciekły metal lub woda schładza się poniżej temperatury topnienia, krzepnięcie rozpoczyna się poprzez nukleację — spontaniczne powstawanie maleńkich zarodków krystalicznych. Klasyczna teoria nukleacji przewiduje, że zarodek musi przekroczyć promień krytyczny r* = 2γTeq/(ΔHf·ΔT), zanim będzie mógł stabilnie rosnąć, gdzie γ to energia powierzchniowa ciało stałe-ciecz, Teq to równowagowa temperatura topnienia, ΔHf to ciepło utajone topnienia, a ΔT to przechłodzenie. Ten symulator umieszcza zarodek krystaliczny w centrum i rozwiązuje równanie pola fazowego ∂φ/∂t = M[∇²φ + φ(1−φ)(φ−1/2+λu)] na siatce różnic skończonych, gdzie φ to parametr porządku pola fazowego (0 = ciecz, 1 = ciało stałe), M to ruchliwość granicy fazowej, λ to stała sprzężenia termicznego, a u = T − Teq to bezwymiarowa temperatura.
Sprzężone równanie dyfuzji cieplnej ∂u/∂t = D·∇²u + (1/2)∂φ/∂t uwzględnia uwalnianie ciepła utajonego w miarę postępu frontu ciała stałego. Anizotropia kryształu jest kodowana poprzez modulację energii granicy funkcją a(θ) = 1 − ε·cos(n·θ), gdzie θ to lokalny kąt normalnej granicy, ε to siła anizotropii, a n to liczba fałd (n=6 dla heksagonalnego lodu, n=4 dla metali kubicznych). Daje to charakterystyczne ramiona dendrytyczne — wydłużone wypustki wzdłuż energetycznie preferowanych kierunków krystalograficznych — które nadają płatkom śniegu ich sześciokrotną symetrię, a odlewom metalowym rozgałęzioną mikrostrukturę.
Najczęściej zadawane pytania
Czym jest model pola fazowego użyty w tej symulacji?
Symulacja rozwiązuje ∂φ/∂t = M[∇²φ + φ(1−φ)(φ−1/2+λu)], gdzie φ to parametr porządku (0=ciecz, 1=ciało stałe), M to ruchliwość granicy, λ sprzęga pole termiczne u=T−Teq z fazą, a wyrazy podwójnej studni φ(1−φ)(φ−1/2) stabilizują oba stany objętościowe. Sprzężone równanie ∂u/∂t = D·∇²u + (1/2)∂φ/∂t transportuje ciepło utajone uwalniane podczas krzepnięcia.
Dlaczego kryształy rosną we wzorze dendrytycznym (drzewiastym)?
Wzrost dendrytyczny wynika z niestabilności morfologicznej: wypustka na granicy ciało stałe-ciecz rośnie szybciej, ponieważ efektywniej uwalnia ciepło utajone do przechłodzonego stopu. Rywalizacja między napięciem powierzchniowym (stabilizującym płaskie granice) a dyfuzją (destabilizującą je) wybiera charakterystyczny promień czubka dendrytu i prędkość wzrostu, znane jako rozwiązanie Ivantsova.
Czym jest przechłodzenie i dlaczego wpływa na kształt kryształu?
Przechłodzenie (ΔT) to o ile temperatura stopu jest niższa od równowagowej temperatury topnienia Teq. Niskie przechłodzenie daje zwarte, wolno rosnące kryształy; wysokie przechłodzenie napędza szybki wzrost dendrytyczny z cienkimi promieniami czubków. Dlatego szybko hartowane metale mają drobną mikrostrukturę dendrytyczną, podczas gdy wolno chłodzone stopy tworzą grube ziarna równoosiowe.
Co kontroluje siła anizotropii ε?
Siła anizotropii ε kontroluje, jak bardzo energia granicy i ruchliwość zmieniają się z orientacją krystalograficzną. Wyższe ε wyostrza preferowane kierunki wzrostu (sześć ramion płatka śniegu) i tłumi wzrost poza osiami. Przy ε=0 granica jest izotropowa, a kryształ rośnie jako koło; zwiększanie ε w kierunku 0,07 daje ostre, dobrze zdefiniowane ramiona dendrytyczne z wąskim rozgałęzieniem bocznym.
Jak zaimplementowana jest 6-krotna symetria kryształu?
Funkcja anizotropii a(θ) = 1 − ε·cos(n·θ) jest stosowana do kąta gradientu granicy θ = atan2(∂φ/∂y, ∂φ/∂x). Ustawienie n=6 tworzy sześć energetycznie preferowanych kierunków oddalonych o 60°, odpowiadających heksagonalnej symetrii kryształu (lód Ih). Wybór n=4 daje symetrię kubiczną (4-krotną) odpowiadającą metalom takim jak żelazo czy aluminium; n=3 daje wzory trygonalne, n=8 quasi-krystaliczne.
Czym jest nukleacja i jak rozpoczyna wzrost kryształu?
Nukleacja to powstawanie maleńkiego zarodka krystalicznego w stopie. Klasyczna teoria nukleacji (CNT) pokazuje, że zarodek musi przekroczyć promień krytyczny r* = 2γTeq/(ΔHf·ΔT), aby być stabilny; mniejsze zarodki rozpuszczają się z powrotem w stopie. W tej symulacji mały zarodek ciała stałego jest umieszczony w centrum, aby ominąć barierę nukleacji i bezpośrednio obserwować anizotropowy wzrost dendrytyczny.
Co decyduje o dyfuzyjności cieplnej w tym modelu?
Dyfuzyjność cieplna D kontroluje, jak szybko ciepło utajone dyfunduje z dala od rosnącej granicy poprzez ∂u/∂t = D·∇²u + (1/2)∂φ/∂t. Wyższe D szybciej odprowadza ciepło, zmniejszając otoczkę termiczną wokół czubka dendrytu i pozwalając na wyższą prędkość czubka. Niższe D tworzy większą warstwę graniczną cieplną, która spowalnia wzrost i może wywołać rozszczepienie czubka lub morfologie „wodorostowe”.
Dlaczego płatki śniegu mają sześciokrotną symetrię, ale wszystkie są różne?
Sześciokrotna symetria płatków śniegu odzwierciedla heksagonalną strukturę kryształu lodu Ih, gdzie cząsteczki wody tworzą pierścienie połączone wiązaniami wodorowymi z kątami wiązań 120°. Wszystkie sześć ramion rośnie jednocześnie w tym samym lokalnym otoczeniu, więc rozwijają się identycznie. Jednak każde ramię kreśli unikalną historię wzrostu, gdy płatek śniegu koziołkuje w powietrzu przy zmieniającej się temperaturze i wilgotności — dzięki czemu każdy płatek śniegu jest unikalny w szczegółach, choć symetryczny sam w sobie.
Jakie są rzeczywiste zastosowania badań nad krzepnięciem dendrytycznym?
Mikrostruktura dendrytyczna bezpośrednio wpływa na właściwości mechaniczne odlewanych metali, stopów łopatek turbin i części drukowanych w 3D. Kontrolowanie odstępu ramion dendrytu (DAS) poprzez szybkość chłodzenia optymalizuje wytrzymałość i ciągliwość. Modele pola fazowego, takie jak ten, są używane przemysłowo do projektowania procesów odlewania, przewidywania mikrosegregacji w niklowych nadstopach oraz optymalizacji parametrów produkcji addytywnej dla komponentów lotniczych.
Jak zachodzi rywalizacja ziaren w krzepnących metalach?
Gdy wiele zarodków powstaje jednocześnie, każdy z nich rozwija dendryty ułożone wzdłuż własnych osi krystalograficznych. Tam, gdzie spotykają się dwa pola dendrytów, szybciej rosnące ziarno (z preferowanymi osiami mniej więcej zgodnymi z gradientem temperatury) blokuje i pochłania wolniej rosnącego sąsiada. Ta rywalizacja ziaren decyduje o ostatecznej mikrostrukturze — ziarna kolumnowe w stopach krzepnących kierunkowo, ziarna równoosiowe w szybko odlewanych materiałach.
Model pola fazowego dendrytycznego krzepnięcia kryształów z regulowanym przechłodzeniem i anizotropią.
3D · silnik renderujący Three.js / WebGL · cel 60 FPS · działa w całości po stronie klienta, bez instalacji