Диффузия управляет эволюцией микроструктуры практически в каждом термомеханическом процессе практического значения: цементации, гомогенизации, термической обработке, сварке, деградации покрытий и ликвации при затвердевании. Физика принципиально хорошо изучена — законы Фика, градиенты химического потенциала, подвижность. Сложность состоит в применении этой физики к реальным многокомпонентным промышленным сплавам с шестью, восемью или десятью элементами, где диффузия каждого элемента зависит от концентраций всех остальных, а межфазные границы движутся по мере образования и растворения фаз в процессе.
DICTRA (DIffusion Controlled TRAnsformation — Превращения, Контролируемые Диффузией) — ответ Thermo-Calc на этот вызов. Это эталонное программное обеспечение для симуляции многокомпонентных диффузионно-контролируемых превращений — инструмент, к которому материаловеды обращаются, когда вопрос стоит не просто «какие фазы стабильны?», а «как быстро они образуются и какие профили концентраций возникают?»
Данное руководство объясняет, что такое DICTRA, как он работает, какие задачи решает и какие базы данных мобильности обеспечивают его функционирование.
Что такое DICTRA?
DICTRA — дополнительный модуль к Thermo-Calc, симулирующий диффузионно-контролируемые реакции в многокомпонентных сплавах. Требует:
- Лицензии на Thermo-Calc (DICTRA — платный дополнительный модуль)
- Хотя бы одной термодинамической базы данных (например, TCFE для сталей, TCNI для никелевых сплавов)
- Хотя бы одной базы данных мобильности (например, MOBFE для сталей, MOBNI для никелевых сплавов) — кинетических данных, необходимых DICTRA сверх термодинамики
Термодинамическая база данных сообщает DICTRA, какие фазы стабильны и каков их состав при локальном равновесии. База данных мобильности сообщает, насколько быстро атомы перемещаются через каждую фазу. Комбинация — термодинамика CALPHAD плюс кинетика, оцениваемая методом CALPHAD, — позволяет DICTRA симулировать реальные промышленные сплавы без подгоночных параметров.
Ключевой принцип: DICTRA основан на численном решении многокомпонентных уравнений диффузии, совмещённых с термодинамикой CALPHAD на подвижных межфазных границах. Предполагается локальное равновесие на межфазных границах — каждая граница смещается так, чтобы химические потенциалы всех компонентов были равны по обе стороны на каждом шаге по времени. Это физически обоснованное допущение для большинства задач диффузии в твёрдом теле при повышенных температурах.
Геометрическое ограничение: симуляции DICTRA одномерны. Поддерживаются три геометрии:
- Плоская — бесконечная пластина (цементация плоского образца, профиль концентраций в диффузионной паре)
- Цилиндрическая — бесконечный цилиндр (диффузия через стенку трубы, растворение стержневидного выделения)
- Сферическая — сфера (растворение сферической частицы, рост сферического выделения)
Это одномерное ограничение достаточно для обширного класса практически важных задач, а три геометрических варианта расширяют охват на многие геометрии с круговым или сферическим поперечным сечением.
Физическая основа: многокомпонентная диффузия
В бинарном сплаве A-B диффузия A относительно B описывается одним коэффициентом диффузии. В многокомпонентном сплаве из N компонентов диффузия каждого компонента зависит от градиентов химического потенциала всех остальных — формируется матрица коэффициентов диффузии размером N×N.
Внедиагональные члены этой матрицы — перекрёстные коэффициенты диффузии — физически означают, что градиент концентрации хрома может приводить диффузию никеля даже при отсутствии градиента самого никеля. В классическом опыте Даркена (знаменитый эксперимент 1949 года с железными парами с разным содержанием кремния) углерод диффундировал вверх по своему собственному градиенту концентраций — из образца с меньшим содержанием углерода в образец с бо́льшим — поскольку кремний обратил знак градиента химического потенциала углерода. Эффекты перекрёстной диффузии реальны, количественно важны в многокомпонентных сплавах, и DICTRA учитывает их в полной мере.
DICTRA формулирует задачу в системе отсчёта, связанной с решёткой (lattice-fixed frame of reference), что корректно описывает потоки вакансий, дрейф Киркендалла и движение межфазных границ из-за разницы молярных объёмов фаз.
Матрица коэффициентов диффузии в DICTRA не хранится напрямую. Вместо этого в базе данных мобильности хранятся параметры мобильности для каждого элемента в каждой фазе. Матрица коэффициентов диффузии вычисляется из этих мобильностей в сочетании с термодинамическим фактором (второй производной энергии Гиббса по составу) — базы данных CALPHAD обеспечивают этот термодинамический фактор. Именно это сопряжение принципиально отличает DICTRA от чисто кинетических инструментов симуляции: термодинамика и кинетика обрабатываются согласованно в рамках метода CALPHAD.
Что симулирует DICTRA: категории применений
1. Химико-термическая обработка поверхности — цементация, азотирование, нитроцементация
Цементация — наиболее исторически важное применение DICTRA. Процесс подразумевает воздействие на сталь углеродсодержащей атмосферой при повышенной температуре, вследствие чего углерод диффундирует в поверхность, повышая твёрдость. То, что выглядит просто в бинарном приближении, становится сложным в реальном сплаве:
- Сталь содержит Cr, Mo, Mn, Ni, Si в дополнение к Fe и C
- Коэффициент диффузии углерода сильно зависит от концентраций легирующих элементов
- Карбиды (M₂₃C₆, M₇C₃, цементит) образуются и растворяются по мере эволюции углеродных профилей
- Положения границ карбид/аустенит движутся в ходе процесса
DICTRA обрабатывает всё это без подгоночных параметров. Входные данные: состав стали, атмосфера (активность углерода), профиль температуры и время. Выходные данные:
- Профиль концентрации углерода по глубине в зависимости от времени
- Профиль обеднения по хрому (Cr расходуется на образование карбидов)
- Профиль фазовых долей (сколько карбида на каждой глубине и какого типа)
- Положение межфазных границ как функция времени
Промышленная значимость: Симуляции цементации в DICTRA верифицированы по экспериментальным измерениям для широкого ряда нержавеющих и инструментальных сталей и высокотемпературных сплавов. Согласие, как правило, отличное — поскольку модель не содержит подгоночных параметров сверх термодинамических и кинетических баз данных. Это позволяет использовать DICTRA предиктивно — проектировать режимы цементации или составы сплавов для заданной глубины цементованного слоя и твёрдости без экспериментальных итераций.
Азотирование и нитроцементация следуют той же вычислительной схеме. Химия отличается (азотная атмосфера, образование нитридов железа и легирующих элементов вместо карбидов), но DICTRA обрабатывает оба процесса в одном фреймворке при наличии соответствующих данных в базе данных мобильности.
2. Гомогенизация литой микроструктуры
Затвердевание порождает микроликвацию — области, обогащённые поздно затвердевающими элементами (межосные области), и области, обеднённые ими (оси дендритов). Гомогенизирующий отжиг при повышенной температуре позволяет диффузии в твёрдом теле уменьшить эти композиционные градиенты. DICTRA симулирует этот процесс.
Ключевые вопросы, на которые отвечает DICTRA:
- Как долго должна длиться изотермическая выдержка при заданной температуре, чтобы снизить неоднородность состава до допустимых норм?
- Какая температура оптимизирует кинетику гомогенизации без риска оплавления?
- Как эволюционирует распределение растворённых элементов при медленном охлаждении после гомогенизирующего отжига?
Настройка: Определяется единичная ячейка, представляющая литую микроструктуру, — как правило, половина межосного расстояния в плоской геометрии, с начальным профилем состава из расчёта затвердевания по Шейлю (который TC-Python или графический режим Thermo-Calc могут предоставить как начальное условие для DICTRA). Затем DICTRA симулирует диффузию в твёрдом теле из этого начального состояния через цикл термической обработки.
Области применения: Слитки жаропрочных никелевых сплавов (гомогенизация перед ковкой), заготовки из алюминиевых сплавов (гомогенизация перед прессованием), слябы из дуплексных нержавеющих сталей (устранение δ-феррита затвердевания перед горячей прокаткой).
3. Диффузионно-контролируемые фазовые превращения
Образование и распад аустенита в сталях: Превращение аустенит/феррит в стали контролируется диффузией, когда происходит в межкритическом температурном интервале, где углерод и замещающие элементы должны перераспределяться между двумя фазами по мере движения границы. DICTRA симулирует:
- Рост аустенита в феррите при межкритическом отжиге (производство двухфазных и TRIP-сталей)
- Рост феррита в аустените при медленном охлаждении
- Профиль углерода в обеих фазах по мере движения границы
- Влияние легирующих элементов (Mn, Si, Cr, Mo) на кинетику превращения
Это задача с подвижной границей — положение границы аустенит/феррит меняется по мере перераспределения углерода, и DICTRA самосогласованно отслеживает это движение.
Аллотропные превращения в титане: Превращение α (ГПУ) / β (ОЦК) в Ti-сплавах обрабатывается аналогично — DICTRA симулирует диффузию β-стабилизирующих элементов (Al, V, Mo, Cr) по мере движения границы α/β при термической обработке, обеспечивая профили концентраций и кинетику превращения.
Рост перлита: DICTRA может симулировать рост перлита из аустенита — совместное эвтектоидное превращение, при котором углерод диффундирует перед продвигающимся фронтом перлита, питая одновременно ламели феррита и цементита.
4. Взаимодиффузия в системе покрытие/подложка
Высокотемпературные покрытия на металлических деталях — жаростойкие покрытия турбинных лопаток, диффузионные барьеры, оцинкованные стали — деградируют со временем по мере взаимодиффузии элементов между покрытием и подложкой. DICTRA количественно описывает эту деградацию.
Покрытие MCrAlX / подложка из никелевого суперсплава: Классический случай: покрытие MCrAlY (M = Ni, Co или NiCo) на монокристаллическом жаропрочном никелевом суперсплаве. Элементы диффундируют в обоих направлениях:
- Al диффундирует из покрытия в подложку (Al, необходимый для образования защитной Al₂O₃, расходуется)
- Тугоплавкие элементы (Re, W, Ta) диффундируют из суперсплава в покрытие, потенциально образуя хрупкие ТКУ-фазы (топологически плотноупакованные)
DICTRA предсказывает:
- Профиль концентрации Al в покрытии и подложке в зависимости от времени при рабочей температуре
- Скорость обеднения β-фазы (NiAl) в покрытии (переход β → γ истощает резерв Al для окислительной защиты)
- Состав и фазовый состав зоны взаимодиффузии
- Оценку ресурса — момент полного исчерпания β-фазы
Оцинкованная сталь: DICTRA симулирует взаимодиффузию Fe-Zn при отжиге оцинкованной стали — образование и рост интерметаллических фаз в диффузионной зоне Fe-Zn, — позволяя выбрать температуру и время процесса для достижения целевой микроструктуры покрытия.
Градиентное спекание твёрдых сплавов (ВК): DICTRA симулирует намеренное создание композиционных градиентов (обогащённых связкой поверхностных зон) в твёрдых сплавах WC-Co при спекании — проектирование режущих инструментов с «двойными свойствами»: вязкий поверхностный слой и твёрдый сердечник.
5. Термическая обработка после сварки (ТТОС)
При сварке разнородных материалов — аустенитная нержавеющая с ферритной сталью, жаропрочный никелевый суперсплав со сталью, разнородные варианты суперсплавов — в зоне соединения возникают крутые градиенты состава. Во время ТТОС элементы диффундируют через эти градиенты, потенциально образуя вредные фазы.
Что рассчитывает DICTRA:
- Образуются ли нежелательные фазы (хрупкие интерметаллиды, карбиды) в зоне соединения при ТТОС
- Как долго должна длиться ТТОС, чтобы перераспределить состав и снизить риск
- Миграция углерода — эффект Даркена приводит к диффузии углерода с ферритной стороны на аустенитную в сварных стыках «сталь-сталь», что вызывает обезуглероживание с одной стороны и науглероживание с другой
6. Рост и растворение частиц
DICTRA симулирует рост или растворение отдельных частиц выделений — как правило, используется сферическая геометрия.
Применения:
- Растворение карбидов при аустенитизации стали — как долго при данной температуре растворяются все карбиды?
- Растворение γ’ в никелевых суперсплавах при термической обработке
- Рост покрытий или интерметаллических слоёв на границах раздела
Критерии остановки (добавлены в Thermo-Calc 2024a, консольный режим): Вместо фиксированного времени симуляции пользователи могут задать условие остановки — например, прекратить симуляцию цементации, когда содержание углерода на глубине 0,5 мм достигнет целевого значения. Это делает DICTRA непосредственно применимым для оптимизации процессов: найти время, необходимое для достижения целевого профиля, а не симулировать до фиксированного момента и затем анализировать результат.
Базы данных мобильности — кинетический каркас
DICTRA требует базы данных мобильности в дополнение к термодинамической. База данных мобильности содержит параметры атомной мобильности — кинетический эквивалент термодинамических параметров энергии Гиббса, оцениваемый методом CALPHAD.
Структура данных о мобильности
Атомные мобильности оцениваются поэлементно для каждой фазы, в которой они актуальны. Для каждого элемента i в фазе φ:
M_i^φ = M_i^φ(T, x) — функция температуры и состава
Мобильность выражается в форме Аррениуса: M_i = A × exp(−Q/RT), где Q — энергия активации диффузии, A — предэкспоненциальный множитель. Оба параметра могут зависеть от состава в многокомпонентных системах — зависимость описывается полиномами Редлих-Кистера по аналогии с термодинамическими параметрами взаимодействия.
Из атомных мобильностей в каждой точке симуляции вычисляется матрица коэффициентов диффузии D как комбинация мобильностей и термодинамических факторов:
D_ij = Σ_k (M_k × (δ_ik − x_i) × ∂μ_k/∂x_j)
где μ_k — химические потенциалы из термодинамической базы данных. Именно это сопряжение обеспечивает физическую самосогласованность DICTRA: база данных мобильности предоставляет атомные транспортные свойства, а термодинамическая база данных — движущие силы.
Доступные базы данных мобильности
| База данных | Система сплавов | Парная термодинамическая БД |
|---|---|---|
| MOBFE9 | Сталь и Fe-сплавы | TCFE15 |
| MOBNI (MOBNI6) | Жаропрочные Ni-сплавы | TCNI |
| MOBAL | Al-сплавы | TCAL |
| MOBTI6 (2026a) | Ti и TiAl-сплавы | TCTI7 |
| MOBHEA4 (2026a) | Высокоэнтропийные сплавы | TCHEA |
| MOBMG | Mg-сплавы | TCMG |
| MOBCU | Cu-сплавы | Различные |
Каждая база данных мобильности разработана в паре с конкретной термодинамической базой. MOBFE9 разрабатывается совместно с TCFE15, MOBNI — совместно с TCNI. Использование несовместимых версий баз данных может привести к некорректным результатам.
Охват базы данных и ограничения
База данных мобильности указывает, для каких фаз имеются данные о диффузии. Фазы, не перечисленные в базе данных мобильности, трактуются как «диффузия отсутствует» (diffusion NONE) — DICTRA учитывает их в расчёте термодинамического равновесия, но считает неподвижными (без диффузии внутри них). Это физически обоснованно для многих применений — карбиды в сталях, например, как правило, имеют очень медленную диффузию, и изменения их состава определяются потоком на границах, а не внутренней диффузией.
Чтобы диффузия внутри фазы симулировалась точно, фаза должна быть явно внесена в базу данных мобильности с оценёнными параметрами. Система на основе Fe (MOBFE) охватывает BCC_A2 (феррит), FCC_A1 (аустенит), ЦЕМЕНТИТ, FE₄N, HCP_A3 и ЖИДКОСТЬ. Система на основе Ni (MOBNI) охватывает FCC_A1, L1₂ (γ’), BCC_A2, B2 и ЖИДКОСТЬ.
Настройка симуляции: ключевые понятия
Области и ячейки
Симуляция DICTRA организована вокруг областей (regions) — пространственных доменов с заданной геометрией и составом. Диффузионная пара включает две области (два сплава). Симуляция цементации — одну область (образец стали). Каждая область может содержать несколько фаз в локальном равновесии.
Ячейка представляет одномерный пространственный охват симуляции с заданными граничными условиями и начальным профилем состава.
Граничные условия
Граничные условия DICTRA задают поведение на поверхностях расчётной области:
Фиксированный состав: Состав на поверхности зафиксирован на заданном значении. Используется для цементации (активность углерода на поверхности фиксируется составом газовой атмосферы).
Нулевой поток: Ни один элемент не пересекает границу. Используется для плоскостей симметрии (ось дендрита при гомогенизации) или запечатанных поверхностей.
Фиксированный поток: Заданный поток элемента поступает через поверхность или уходит через неё. Используется, когда скорость переноса в газовой фазе ограничена.
Граничное условие по активности: Состав поверхности поддерживается при заданном химическом потенциале (активности). Связано с фиксированным составом, но задаётся термодинамически.
Профили температуры
Температура в DICTRA может задаваться как константа или как функция времени, используя любую комбинацию:
- Линейные нагревы/охлаждения (управляемый нагрев или охлаждение)
- Изотермические выдержки
- Произвольные кривые термического цикла (для реалистичных промышленных режимов)
Это позволяет симулировать полные промышленные циклы термической обработки: нагрев → выдержка → охлаждение → старение — в одном запуске DICTRA.
Измельчение расчётной сетки
DICTRA использует метод конечных объёмов с настраиваемым шагом сетки. Узлы сетки могут располагаться равномерно или геометрически (гуще у межфазных границ, где крутые концентрационные градиенты). Правильное измельчение сетки важно для точности, особенно вблизи подвижных межфазных границ.
Модели DICTRA для различных физических ситуаций
Модель чёткой границы (подвижная граница)
Модель DICTRA по умолчанию. Фазы разделены чёткими, хорошо определёнными границами. Положение границы смещается по мере роста или уменьшения фаз. На границе предполагается локальное термодинамическое равновесие — составы по обе стороны границы определяются термодинамической базой данных на каждом шаге по времени. Модель физически уместна, когда:
- Границы движутся по кинетике, контролируемой диффузией
- Межфазные границы атомно резкие относительно масштаба диффузии
Применения: Превращение аустенит/феррит в стали, цементация с образованием карбидов, взаимодиффузия покрытие/подложка, рост и растворение частиц.
Модель гомогенизации
Для многофазных областей, где фазы перемешаны в масштабе меньше разрешения симуляции (например, двухфазная область α+β при гомогенизации), модель подвижной границы не может точно представить геометрию. Модель гомогенизации трактует многофазную смесь как единую эффективную среду:
- Локальное равновесие выполняется в каждом элементарном объёме (фазовые доли и составы из CALPHAD)
- Диффузия усредняется по смеси с использованием локально усреднённых кинетических свойств
- Вычислительно преобразует многофазную задачу в задачу эффективной однофазной диффузии
Модель уместна, когда масштаб двухфазной смеси (например, межосное расстояние при гомогенизации, межламельное расстояние в перлите) значительно меньше диффузионных расстояний, представляющих интерес.
Применения: Гомогенизация литых микроструктур (двухфазная структура в масштабе дендрита), диффузия через кермет-покрытия, градиентное спекание твёрдых сплавов.
Интеграция с TC-PRISMA и TC-Python
DICTRA и TC-PRISMA совместно
DICTRA и TC-PRISMA (модуль осаждения) описывают различные масштабы кинетического поведения:
- DICTRA: дальнодействующая диффузия — профили концентраций на расстояниях от микрон до миллиметров, подвижные межфазные границы, профили поверхностной обработки
- TC-PRISMA: короткодействующая кинетика — зародышеобразование, рост и огрубление выделений внутри однофазной матрицы
Оба могут использоваться последовательно: DICTRA рассчитывает профиль концентраций после термической обработки, а TC-PRISMA рассчитывает, как зарождаются и растут выделения в этом профиле. Например, при моделировании цементованного слоя стали — DICTRA предоставляет профили углерода и хрома, а TC-PRISMA моделирует кинетику выделения карбидов на каждой глубине.
DICTRA в TC-Python
Все расчёты DICTRA доступны через TC-Python, обеспечивая автоматизированную диффузионную симуляцию в исследовательских и оптимизационных рабочих процессах:
- Перебор режимов цементации для ряда составов сталей и температур
- Оптимизация циклов гомогенизации нового сплава перебором комбинаций температура-время
- Включение результатов симуляций DICTRA в обучающие наборы для машинного обучения
- Стыковка выходных данных DICTRA с конечно-элементными моделями
Валидация: почему DICTRA можно доверять
Доверие к DICTRA основывается на той же базе CALPHAD, что и термодинамика Thermo-Calc:
Методология разработки баз данных: Параметры атомных мобильностей в MOBFE, MOBNI и других базах данных оцениваются по экспериментальным измерениям — коэффициентам трассерной диффузии, коэффициентам взаимодиффузии из опытов с диффузионными парами и измерениям профилей концентраций. Процесс оценки — та же методология CALPHAD, применённая к кинетике: оптимизация параметров для воспроизведения известных экспериментальных данных с последующей экстраполяцией на многокомпонентные системы.
Отсутствие подгоночных параметров в применениях: После оценки баз данных симуляция цементации или гомогенизации в DICTRA не использует свободных параметров — базы данных обеспечивают всю физику. Это принципиально отличает DICTRA от инструментов симуляции, требующих параметров подгонки, регулируемых для соответствия каждому новому экспериментальному набору данных.
Обширная литературная валидация: DICTRA верифицирован в сотнях опубликованных исследований по экспериментальным измерениям:
- Профилей углерода, азота и хрома в цементованных сталях и высокотемпературных сплавах
- Профилей концентраций в диффузионных парах (Fe-Cr-Ni, Ni-Al-Cr, Ti-Al-V и многих других)
- Взаимодиффузии в системах MCrAlY / суперсплав
- Кинетики гомогенизации в никелевых суперсплавах и алюминиевых сплавах
Согласие между симуляциями DICTRA и экспериментом, как правило, укладывается в погрешность измерений при использовании соответствующих баз данных — весомая валидация сопряжённого термодинамического/кинетического подхода CALPHAD.
Ограничения
Понимание того, где DICTRA неуместен:
Только одно измерение: DICTRA не способен симулировать подлинно двух- или трёхмерные диффузионные поля — например, диффузию вокруг угла, через неоднородную геометрию шва или вокруг частицы нерегулярной формы. Для трёхмерных диффузионных задач требуются фазово-полевые методы (MICRESS, OpenPhase, MOOSE) или МКЭ-коды с термодинамикой, связанной с DICTRA.
Допущение о локальном равновесии: DICTRA предполагает термодинамическое локальное равновесие на межфазных границах. Это допущение нарушается при:
- Очень высоких скоростях охлаждения (мартенситные превращения)
- Массивных превращениях, где миграция границ очень быстрая
- Системах со значительным захватом растворённых атомов при высоких скоростях роста
Для таких случаев более уместны модели пара-равновесия или квази-пара-равновесия.
Покрытие баз данных: DICTRA может симулировать только фазы и системы сплавов, для которых оценены данные о мобильности. Фаза с пометкой «диффузия NONE» не имеет кинетического описания — DICTRA учитывает её равновесное присутствие, но не может симулировать диффузию внутри неё.
Отсутствие диффузии по границам зёрен: DICTRA моделирует решёточную диффузию. Диффузия по границам зёрен — как правило, на 2–5 порядков быстрее решёточной — в стандартную модель не включена. Для процессов, где диффузия по границам зёрен важна (например, низкотемпературная цементация нержавеющих сталей), требуются дополнительные подходы к моделированию.
Итог
DICTRA — инструмент вычислительной металлургии для количественной симуляции многокомпонентных диффузионно-контролируемых превращений: от химико-термической обработки при цементации до взаимодиффузии покрытий до гомогенизации микроликвации затвердевания. Его сила — в строгом сопряжении термодинамики CALPHAD (базы данных энергий Гиббса) с кинетикой, оцениваемой методом CALPHAD (базы данных мобильности), что обеспечивает симуляции реальных промышленных сплавов без подгоночных параметров.
Для любой инженерной задачи материаловедения, где вопрос стоит не только «какие фазы стабильны?», но и «как эволюционируют профили концентраций и как быстро движутся межфазные границы?», DICTRA является стандартным инструментом. Его более чем 30-летняя валидация в первичной литературе в сочетании с постоянно обновляемыми термодинамическими и кинетическими базами данных Thermo-Calc делает его наиболее доверенным пакетом для симуляции диффузии в вычислительном материаловедении.
По вопросам лицензирования DICTRA и Thermo-Calc обращайтесь в Telegram: t.me/DoCrackMe
Смотрите также: Thermo-Calc 2026a — Что нового: полный обзор | TC-Python — Автоматизация расчётов Thermo-Calc на Python | Thermo-Calc vs FactSage vs Pandat — Сравнение CALPHAD-программ
