PL 
EnglishDomyślny
Українська💎 Dyfuzja w krysztale
Symulacja dyfuzji atomów w sieci krystalicznej mechanizmem wakancyjnym z drugim prawem Ficka i dyfuzyjnością Arrheniusa.
O dyfuzji w krysztale
Ta symulacja modeluje, jak atomy rozpuszczonej substancji rozprzestrzeniają się w ciele krystalicznym poprzez mechanizm wakancyjny, w którym atom przeskakuje do sąsiedniego pustego miejsca sieciowego. Cienkie źródło atomów jest umieszczone w centrum, a jego profil stężenia C(x,t) jest przesuwany w czasie poprzez numeryczne rozwiązywanie drugiego prawa Ficka, ∂C/∂t = D·∂²C/∂x², za pomocą jawnego schematu różnic skończonych na jednowymiarowej siatce 120 komórek.
Współczynnik dyfuzji jest ustalany przez prawo Arrheniusa D = D₀·exp(−Eₐ/k_BT). Cztery suwaki dostosowują temperaturę T (300–1800 K), energię aktywacji Eₐ (0,3–3,0 eV), współczynnik przedwykładniczy D₀ oraz skalę kroku czasowego. Przerywana krzywa Gaussa pokazuje rozwiązanie analityczne. Taka dyfuzja rządzi rzeczywistymi procesami, takimi jak nawęglanie stali, domieszkowanie półprzewodników i pełzanie w stopach turbinowych.
Najczęściej zadawane pytania
Co pokazuje ta symulacja?
Pokazuje atomy rozpuszczonej substancji dyfundujące na zewnątrz z cienkiego, centralnego źródła przez sieć krystaliczną. Górny panel koloruje atomy według lokalnego stężenia, a dolny panel przedstawia zmieniający się w czasie profil stężenia C(x,t) w funkcji położenia w nanometrach, wraz z analitycznym przewidywaniem gaussowskim.
Czym jest mechanizm wakancyjny?
W większości ciał krystalicznych atom porusza się, przeskakując do sąsiedniego wolnego miejsca w sieci, zamiast przeciskać się między atomami. Każdy skok musi pokonać barierę energii migracji, więc dyfuzja silnie zależy zarówno od liczby dostępnych wakancji, jak i od dostępnej energii cieplnej.
Jakie równanie rządzi profilem?
Profil podlega drugiemu prawu Ficka, ∂C/∂t = D·∂²C/∂x². Symulacja rozwiązuje je za pomocą jawnego schematu różnic skończonych, automatycznie dzieląc każdy krok tak, by współczynnik stabilności D·dt/dx² pozostawał na poziomie 0,4 lub niższym.
Jak obliczany jest współczynnik dyfuzji?
Dyfuzyjność wykorzystuje zależność Arrheniusa D = D₀·exp(−Eₐ/k_BT), gdzie stała Boltzmanna k_B = 8,617×10⁻⁵ eV/K. Podwyższenie temperatury lub obniżenie energii aktywacji gwałtownie zwiększa D, ponieważ znajduje się ona wewnątrz wykładnika.
Co kontrolują cztery suwaki?
Temperatura T (300–1800 K) oraz energia aktywacji Eₐ (0,3–3,0 eV) wspólnie ustalają dyfuzyjność poprzez prawo Arrheniusa. Współczynnik przedwykładniczy D₀ przeskalowuje ogólną szybkość przeskoków, a skala kroku czasowego przyspiesza lub spowalnia postęp symulowanego czasu.
Dlaczego profil ma kształt Gaussa?
Źródło punktowe lub cienkowarstwowe dyfundujące zgodnie z drugim prawem Ficka ewoluuje w krzywą Gaussa, której szerokość rośnie z czasem. Rozwiązanie analityczne to C(x,t) ∝ exp(−x²/4Dt), więc rozrzut jest w całości ustalany przez iloczyn Dt. Przerywana zielona krzywa nakłada się na ten dokładny wynik.
Co oznacza √(Dt)?
Wielkość √(Dt) to długość dyfuzji, miara tego, jak daleko typowo przemieściły się atomy. Połowa szerokości krzywej Gaussa skaluje się z nią, więc podwojenie czasu zwiększa rozrzut jedynie około 1,4-krotnie, a nie dwukrotnie.
Czy symulacja jest fizycznie dokładna?
Dyfuzyjność Arrheniusa i drugie prawo Ficka są zgodne z podręcznikami, a solver różnic skończonych ściśle odpowiada analitycznej krzywej Gaussa. Jest to jednak uproszczony, jednowymiarowy model kontinuum: pomija dyskretne skoki atomowe, klastrowanie defektów, granice ziaren oraz zależność dyfuzyjności od stężenia.
Dlaczego ogrzewanie tak dramatycznie przyspiesza dyfuzję?
Ponieważ temperatura występuje wewnątrz wykładnika prawa Arrheniusa, umiarkowany wzrost T powoduje duży wzrost D. Więcej atomów zyskuje wystarczającą energię cieplną, by pokonać barierę migracji, więc kilkaset kelwinów może zmienić dyfuzyjność o rzędy wielkości.
Gdzie dyfuzja w krysztale ma znaczenie w przemyśle?
Leży u podstaw wielu procesów materiałowych: nawęglanie i azotowanie utwardzające powierzchnie stali, domieszkowanie płytek krzemowych w celu wytworzenia tranzystorów, homogenizacja odlewanych stopów, spiekanie ceramiki oraz powolna deformacja pełzania, która ogranicza żywotność łopatek turbin silników odrzutowych w wysokiej temperaturze.