💎 Дифузія в Кристалі
Симулюйте дифузію атомів через кристалічну гратку за вакансійним механізмом.
Про дифузію в кристалі
Ця симуляція моделює, як атоми розчиненої речовини поширюються крізь кристалічне тверде тіло вакансійним механізмом, коли атом перестрибує в сусідній порожній вузол ґратки. Тонке джерело атомів розміщене в центрі, і його профіль концентрації C(x,t) розвивається в часі шляхом чисельного розв'язання другого закону Фіка, ∂C/∂t = D·∂²C/∂x², явною схемою скінченних різниць на одновимірній сітці зі 120 комірок.
Коефіцієнт дифузії задається законом Арреніуса D = D₀·exp(−Eₐ/kₕT). Чотири повзунки регулюють температуру T (300–1800 К), енергію активації Eₐ (0,3–3,0 еВ), передекспоненційний множник D₀ та масштаб кроку часу. Пунктирна крива Гаусса показує аналітичний розв'язок. Така дифузія керує реальними процесами: цементацією сталі, легуванням напівпровідників та повзучістю в сплавах турбін.
Часті запитання
Що показує ця симуляція?
Вона показує, як атоми розчиненої речовини дифундують назовні з тонкого центрального джерела крізь кристалічну ґратку. Верхня панель забарвлює атоми за локальною концентрацією, тоді як нижня панель показує еволюцію профілю концентрації C(x,t) залежно від положення в нанометрах, поряд з аналітичним прогнозом Гаусса.
Що таке вакансійний механізм?
У більшості кристалічних твердих тіл атом рухається, перестрибуючи в сусідній порожній вузол ґратки, а не протискаючись між атомами. Кожен стрибок має подолати енергетичний бар'єр міграції, тому дифузія сильно залежить як від наявності вакансій, так і від доступної теплової енергії.
Яке рівняння керує профілем?
Профіль підпорядковується другому закону Фіка, ∂C/∂t = D·∂²C/∂x². Симуляція розв'язує його явним оновленням скінченних різниць, автоматично поділяючи кожен крок так, щоб коефіцієнт стійкості D·dt/dx² залишався не вище 0,4.
Як обчислюється коефіцієнт дифузії?
Дифузійність використовує співвідношення Арреніуса D = D₀·exp(−Eₐ/kₕT), зі сталою Больцмана kₕ = 8,617×10⁻⁵ еВ/К. Підвищення температури або зниження енергії активації різко збільшує D, оскільки вона перебуває всередині експоненти.
Що контролюють чотири повзунки?
Температура T (300–1800 К) та енергія активації Eₐ (0,3–3,0 еВ) разом визначають дифузійність через закон Арреніуса. Передекспоненційний множник D₀ перемасштабовує загальну швидкість стрибків, а масштаб кроку часу прискорює чи сповільнює перебіг симульованого часу.
Чому профіль має форму Гаусса?
Точкове або тонкошарове джерело, що дифундує за другим законом Фіка, еволюціонує в гаусову криву, ширина якої зростає з часом. Аналітичний розв'язок C(x,t) ∝ exp(−x²/4Dt), тож поширення повністю визначається добутком Dt. Пунктирна зелена крива накладає цей точний результат.
Що означає √(Dt)?
Величина √(Dt) — це довжина дифузії, міра того, наскільки далеко зазвичай мігрували атоми. Напівширина гаусової кривої масштабується з нею, тож подвоєння часу збільшує поширення лише приблизно в 1,4 раза, а не вдвічі.
Чи фізично точна ця симуляція?
Дифузійність за Арреніусом та другий закон Фіка є підручниково точними, а розв'язувач скінченних різниць добре узгоджується з аналітичною гаусовою кривою. Однак це спрощена одновимірна континуальна модель: вона нехтує дискретними атомними стрибками, кластеризацією дефектів, межами зерен та концентраційно залежною дифузійністю.
Чому нагрівання так різко прискорює дифузію?
Оскільки температура входить в експоненту закону Арреніуса, помірне підвищення T дає велике зростання D. Більше атомів отримують достатньо теплової енергії, щоб подолати бар'єр міграції, тож кілька сотень кельвінів можуть змінити дифузійність на порядки величини.
Де дифузія в кристалі важлива в промисловості?
Вона лежить в основі багатьох матеріалознавчих процесів: цементації та азотування для зміцнення поверхні сталі, легування кремнієвих пластин для виготовлення транзисторів, гомогенізації литих сплавів, спікання кераміки та повільної повзучості, яка обмежує термін служби лопаток турбін реактивних двигунів за високої температури.