Кристаллизация металлов: зарождение и рост зёрен
При охлаждении жидкого металла образование твёрдой фазы включает две разные задачи: появление устойчивых зародышей и последующее продвижение границы раздела. Температура ниже равновесной температуры плавления создаёт термодинамический стимул к затвердеванию, но не устраняет затраты на новую поверхность. Поэтому мелкий кластер ещё не является устойчивым кристаллом. Частота появления зародышей и скорость их роста совместно влияют на количество и размеры зёрен; подменять одну характеристику другой нельзя.
Классическая сферическая модель
Рассмотрим гомогенное зарождение в объёме чистого переохлаждённого расплава. Обозначим γ положительную удельную энергию границы твёрдое–жидкое, Дж/м², а Δg положительную величину объёмного выигрыша свободной энергии, Дж/м³. Для сферического зародыша радиуса r изменение энергии равно ΔG(r) = 4πr²γ − (4π/3)r³Δg. Первый член описывает поверхность, второй — объём. При малом r доминируют затраты на границу; при большом r возрастает роль объёмного выигрыша.
Условие dΔG/dr = 0 даёт критический радиус r* = 2γ/Δg. Вершина энергетического барьера равна ΔG* = 16πγ³/(3Δg²). Зародыш меньше r* имеет тенденцию растворяться, а больше r* — расти при сохранении условий. В самой вершине барьера ΔG ещё положительна. Радиус, при котором суммарное изменение энергии обращается в нуль, равен 3γ/Δg и потому больше критического радиуса. Это полезная проверка различия между устойчивостью дальнейшего роста и общим энергетическим выигрышем образования кластера.
Числовой пример с условными величинами
Зададим γ = 0,200 Дж/м², температуру равновесного плавления T_m = 1500 К и объёмную скрытую теплоту плавления L_v = 1,00·10⁹ Дж/м³. Все параметры модельные; они не приписаны определённому металлу. Для небольшого относительно T_m переохлаждения ΔT используем линейную оценку Δg ≈ L_v ΔT/T_m. При ΔT = 100 К получаем Δg = 6,6667·10⁷ Дж/м³, r* = 6,000 нм и ΔG* = 3,01593·10⁻¹⁷ Дж.
| Переохлаждение, К | Δg, Дж/м³ | r*, нм | ΔG*, Дж |
|---|---|---|---|
| 50 | 3,3333·10⁷ | 12,000 | 1,20637·10⁻¹⁶ |
| 100 | 6,6667·10⁷ | 6,000 | 3,01593·10⁻¹⁷ |
Таблицу можно прокручивать по горизонтали →
Удвоение переохлаждения в этой модели уменьшает r* вдвое, а барьер вчетверо. При r = 6 нм поверхностный вклад составляет 9,04779·10⁻¹⁷ Дж, объёмный — −6,03186·10⁻¹⁷ Дж; их сумма совпадает с барьером. При r = 9 нм суммарная энергия равна нулю. Проверка единиц обязательна: если теплоту плавления взяли в Дж/кг, сначала умножают её на плотность и получают Дж/м³.
Почему реальное литьё отличается
В производственном расплаве зарождение часто происходит на стенках формы и частицах. Смачивание подложки меняет геометрию зародыша и энергетический барьер, поэтому сферический расчёт не определяет необходимую концентрацию модификатора. При росте выделяется теплота; её отвод и перенос растворённых элементов меняют локальные условия перед границей. На охлаждающей кривой эти процессы могут проявляться температурной задержкой или рекалесценцией, но для количественного вывода необходимы калиброванный датчик и учёт его инерции.
Рассчитанный барьер не задаёт время кристаллизации. Для частоты зарождения нужны кинетический множитель, число потенциальных центров и температурная зависимость подвижности атомов. Классическая модель предполагает резкую границу, постоянную изотропную γ и объёмные свойства обеих фаз; на размерах нескольких атомов эти предпосылки становятся сомнительными. При заметном переохлаждении следует рассчитывать Δg по термодинамическим данным, а не продолжать линейную формулу. Для сплавов дополнительно учитывают состав жидкой и твёрдой фаз. Практический выбор режима охлаждения опирают на температурные измерения и металлографию, сопоставляя форму зёрен с реальными условиями теплоотвода.
Первичные источники и область использования
- Kiel University — Homogeneous Nucleation: баланс поверхностной и объёмной энергии, критический радиус и зависимость от переохлаждения.