Коли рідкий метал або вода охолоджується нижче точки плавлення, кристалізація починається через нуклеацію — спонтанне утворення дрібних кристалічних зародків. Класична теорія нуклеації передбачає, що зародок мусить перевищити критичний радіус r* = 2γTeq/(ΔHf·ΔT), щоб стабільно рости, де γ — поверхнева енергія межі тверде-рідке, Teq — рівноважна точка плавлення, ΔHf — прихована теплота плавлення, а ΔT — переохолодження. Ця симуляція розміщує кристалічне зерно в центрі та розв'язує фазово-польове рівняння ∂φ/∂t = M[∇²φ + φ(1−φ)(φ−1/2+λu)] на сітці скінченних різниць, де φ — параметр порядку (0 = рідина, 1 = тверде тіло), M — рухливість межі, λ — константа теплового зв'язку, а u = T − Teq — безрозмірна температура.
Пов'язане рівняння теплової дифузії ∂u/∂t = D·∇²u + (1/2)∂φ/∂t враховує виділення прихованої теплоти в міру просування твердого фронту. Анізотропія кристала кодується через модуляцію поверхневої енергії функцією a(θ) = 1 − ε·cos(n·θ), де θ — кут нормалі до місцевої межі, ε — сила анізотропії, а n — кратність (n=6 для гексагонального льоду, n=4 для кубічних металів). Це породжує характерні дендритні гілки — витягнуті виступи вздовж енергетично переважних кристалографічних напрямків — що надають сніжинкам їхню шестикутну симетрію, а литим металам — галузисту мікроструктуру.
Часті запитання
Яка фазово-польова модель використовується в цій симуляції?
Симуляція розв'язує ∂φ/∂t = M[∇²φ + φ(1−φ)(φ−1/2+λu)], де φ — параметр порядку (0=рідина, 1=тверде тіло), M — рухливість межі, λ з'єднує теплове поле u=T−Teq з фазою. Подвійна яма φ(1−φ)(φ−1/2) стабілізує обидва стани. Рівняння ∂u/∂t = D·∇²u + (1/2)∂φ/∂t переносить приховану теплоту, що виділяється під час кристалізації.
Чому кристали ростуть у деревоподібному (дендритному) патерні?
Дендритний ріст виникає через морфологічну нестабільність: виступ на межі тверде-рідке росте швидше, оскільки ефективніше виділяє приховану теплоту в переохолоджений розплав. Конкуренція між поверхневим натягом та дифузією визначає характерний радіус кінчика дендриту та швидкість росту — розв'язок Іванцова.
Що таке переохолодження і чому воно впливає на форму кристала?
Переохолодження (ΔT) — це наскільки температура розплаву нижча за рівноважну Teq. Мале переохолодження дає компактні кристали; велике — швидкий дендритний ріст. Саме тому швидко загартовані метали мають тонку дендритну мікроструктуру, а повільно охолоджені — грубі рівноосні зерна.
Що контролює сила анізотропії ε?
Сила анізотропії ε визначає, наскільки енергія та рухливість межі змінюються з кристалографічною орієнтацією. Більша ε загострює переважні напрямки росту (6 рук сніжинки). При ε=0 межа ізотропна і кристал росте колом; зі збільшенням ε до 0.07 з'являються чіткі дендритні гілки з тонким бічним галуженням.
Як реалізована 6-кратна симетрія кристала?
Функція анізотропії a(θ) = 1 − ε·cos(n·θ) застосовується до кута градієнта межі θ = atan2(∂φ/∂y, ∂φ/∂x). При n=6 утворюються шість енергетично переважних напрямків через 60°, що відповідає симетрії льоду Ih. При n=4 — кубічна симетрія металів; n=3 — тригональна; n=8 — квазікристалічна.
Що таке нуклеація і як вона починає ріст кристала?
Нуклеація — утворення дрібного кристалічного зародка в розплаві. Зародок мусить перевищити критичний радіус r* = 2γTeq/(ΔHf·ΔT), щоб стати стабільним; менші зародки розчиняються назад. У симуляції маленьке тверде зерно розміщується в центрі для обходу бар'єру нуклеації та безпосереднього спостереження дендритного росту.
Що визначає теплову дифузивність у цій моделі?
Теплова дифузивність D визначає, наскільки швидко прихована теплота розсіюється від межі росту через ∂u/∂t = D·∇²u + (1/2)∂φ/∂t. Вища D швидше відводить теплоту, дозволяючи вищу швидкість кінчика. Нижча D створює більший тепловий шар, що уповільнює ріст і може спричинити розщеплення кінчика або «водоростеподібні» морфології.
Чому сніжинки шестикутно симетричні, але всі різні?
6-кратна симетрія сніжинок відображає кристалічну структуру льоду Ih, де молекули води утворюють водневі зв'язки з кутами 120°. Всі шість рук ростуть одночасно в однаковому середовищі, тому розвиваються однаково. Однак кожна сніжинка переміщується у повітрі з різною температурою і вологістю — тому кожна є унікальною в деталях, залишаючись симетричною.
Де на практиці застосовуються знання про дендритну кристалізацію?
Дендритна мікроструктура безпосередньо впливає на механічні властивості литих металів, сплавів для турбінних лопаток та деталей адитивного виробництва. Відстань між гілками дендриту (DAS) оптимізують через швидкість охолодження. Фазово-польові моделі використовуються промислово для проектування процесів лиття та прогнозування мікросегрегації в нікелевих суперсплавах.
Як відбувається конкуренція зерен у металах, що тверднуть?
Коли одночасно утворюються кілька зародків, кожен росте з дендритами вздовж власних кристалографічних осей. Де два дендритних поля зустрічаються, зерно, що росте швидше (з осями вздовж градієнта температури), блокує та поглинає сусіднє — «конкуренція зерен». Це визначає фінальну мікроструктуру: стовпчасті зерна в спрямовано затверділих сплавах, рівноосні в швидко литих матеріалах.