КАРТОЧКА ПРОЕКТА ФУНДАМЕНТАЛЬНЫХ И ПОИСКОВЫХ НАУЧНЫХ ИССЛЕДОВАНИЙ,
ПОДДЕРЖАННОГО РОССИЙСКИМ НАУЧНЫМ ФОНДОМ

Информация подготовлена на основании данных из Информационно-аналитической системы РНФ, содержательная часть представлена в авторской редакции. Все права принадлежат авторам, использование или перепечатка материалов допустима только с предварительного согласия авторов.

ОБЩИЕ СВЕДЕНИЯ

Номер 21-79-10012

НазваниеФазовые превращения в метастабильной жидкости с приложениями к кристаллизации бинарных материалов

РуководительТоропова Любовь Валерьевна, Кандидат физико-математических наук

Организация финансирования, регион федеральное государственное автономное образовательное учреждение высшего образования "Уральский федеральный университет имени первого Президента России Б.Н. Ельцина", Свердловская обл

Период выполнения при поддержке РНФ

07.2021 - 06.2024

Конкурс№61 - Конкурс 2021 года «Проведение исследований научными группами под руководством молодых ученых» Президентской программы исследовательских проектов, реализуемых ведущими учеными, в том числе молодыми учеными.

Область знания, основной код классификатора 09 - Инженерные науки, 09-203 - Фазовые равновесия и превращения

Ключевые словаФазовые переходы, тепломассоперенос, нуклеация, рост кристаллов, дендриты, двухфазная зона, математическое моделирование, компьютерные симуляции

Код ГРНТИ27.35.45

ИНФОРМАЦИЯ ИЗ ЗАЯВКИ

Аннотация
Проект направлен на развитие теории фазовых превращений – теоретическое описание совместного протекания объёмной и направленной кристаллизации, а также создание программных систем для проведения компьютерных симуляций. Конечной целью комплексного математического моделирования является выбор оптимальных режимов реализации и управление процессами кристаллизации для получения материалов с заданными свойствами и микроструктурой. Решение проблемы будет основано на объединении моделей нескольких принципиально разных этапов кристаллизации: образование жизнеспособных зародышей, их последующий рост и агрегирование, их дальнейшее укрупнение, а также рост твёрдой фазы с образованием ветвящихся дендритных кристаллов, ответственных за структуру и свойства материалов. Исследуемые модели станут единой теоретической базой программных комплексов, разработанных в рамках проекта для симуляции кристаллизации вещества при наличии совместного (объёмного и направленного) фазового превращения в метастабильных жидкостях. Актуальность исследований обусловлена широким применением теории объёмной и направленной кристаллизации в металлургии, химической промышленности, биотехнологической сфере и различных областях прикладной физики. Так, например, динамика метастабильной области двухфазного состояния вещества определяет микроструктуру формирующейся в результате кристаллизации твёрдой фазы материала (например, образование примесных ячеек, слоистой ликвации примеси, переходы между монокристаллической и поликристаллической структурами). Математические задачи, поставленные в настоящем проекте, описывают реальные процессы фазовых переходов, которые осложнены конвективными течениями расплава, процессами зарождения и роста кристаллов в метастабильной области фазового превращения. Математическое моделирование кристаллизации позволяет существенным образом оптимизировать изучаемые процессы, управлять ими и получать материалы с заданными расчётными характеристиками. Разрабатываемые программные комплексы призваны объединить и систематизировать современные математические модели по отдельным направлениям; разработать на их основе средство для проведения комплексных (многомасштабных) компьютерных симуляций с применением современных алгоритмов. Это позволит верифицировать накопленный задел и вновь разработанные модели по отдельности, а также оценить синергетический эффект от их совместного использования. Развитые в рамках проекта математические модели, а также выведенные на их основе аналитические зависимости и проведённые компьютерные симуляции будут применены для описания реальных экспериментальных данных о росте кристаллов в пересыщенных растворах и переохлаждённых расплавах. Так, например, найденный профиль доли твёрдой фазы в двухфазной области позволит вычислить междендритное расстояние и определить пористость кристаллизующихся материалов. Учёт конвективного тепломассопереноса в жидкой фазе, а также совместного протекания объёмного и направленного фазового превращения является новизной проекта. Научная идея состоит в объединении математических моделей направленной и объёмной кристаллизации, а также модели конвективного тепломассопереноса с теорией кристаллизации в области фазового перехода, которая заполнена растущими элементами твёрдого материала. Постановка заявленной проблемы заключается в формулировке замкнутой системы интегро-дифференциальных уравнений, граничных и начальных условий для описания задач проекта. Решение сформулированных задач будет основано на многолетнем опыте работы коллектива авторов проекта в данной области, его научном заделе и полученных ранее результатах. Решение задач проекта будет осуществляться с помощью ранее развитых методов и подходов для решения аналогичных моделей структурно-фазового превращения с метастабильной областью, классических аналитических методах решения задач с подвижными границами (методы малого параметра, дифференциальных рядов, интегральных преобразований, седловой точки, теории устойчивости), а также численного моделирования. Программные системы будут реализованы на базе свободного программного обеспечения (OpenFOAM) в качестве базы для создания гибкого (в отношении возможностей для модификации моделей) программного комплекса, позволяющего рассчитывать динамику фазовых переходов для различных конвективных течений (естественная и вынужденная конвекция).

Ожидаемые результаты
В рамках проекта будет развита теория совместного протекания направленного и объёмного фазового перехода в метастабильной движущейся двухфазной области с учётом зарождения и роста в ней кристаллов, эволюции дендритоподобных структур, наличия вынужденной и естественной конвекции. Развиваемая теория позволит определить законы движения границ фазового перехода, установить функцию распределения кристаллов по размерам и закон снятия переохлаждения в двухфазной области, найти температурное поле и поле концентрации примеси, определить долю твёрдой фазы, ответственную за междендритное расстояние и пористость материала. Развиваемая теория и проведённое на её основе компьютерное моделирование будут сопоставлены с экспериментальными данными по кристаллизации реальных расплавов из переохлаждённой жидкости (например, Al-Li, Ti-Al, Fe-Ni). Эти сплавы перспективны для использования в различных отраслях промышленности. Например, литий придает уникальные свойства алюминиевым сплавам - повышение прочности и модуля упругости при снижении плотности. Его главное преимущество в том, что многие алюминий-литиевые сплавы можно сваривать и перейти в ближайшем будущем к сварным соединениям в авиационной промышленности. Кроме того, данные сплавы обладают наилучшим комплексом механических, эксплуатационных и коррозионных характеристик, позволяющих конкурировать с традиционными алюминиевыми сплавами и полимерными композиционными материалами. Поэтому определение их свойств и микроструктуры с помощью теоретического и компьютерного моделирования существенно снизит экономические издержки. Эффектом в социальной сфере станет создание новых рабочих мест в металлургической промышленности для производства алюминий-литиевых сплавов. Настоящее исследование является оригинальным, новым и соответствует мировому уровню, что подтверждается предыдущими публикациями коллектива в высокорейтинговых журналах по данному научному направлению. Развитие теории совместного объёмного и направленного фазового превращения позволит определить микроструктуру твёрдой фазы в чистых (однокомпонентных) и бинарных расплавах в условиях конвективного тепломассопереноса вблизи межфазной границы твёрдая фаза – жидкость. Будут сформулированы новые математические модели фазового превращения в подвижной метастабильной области при наличии конвективного тепломассопереноса, найдены квазистационарные и нестационарные решения, описывающие температурное и концентрационное поля, будет найдена скорость кристаллизации, выполнен линейный анализ морфологической устойчивости при конвективном тепломассопереносе, выведены уравнения для возмущений температурного и концентрационного полей, определено дисперсионное соотношение и кривая нейтральной устойчивости. Такой анализ позволит найти области теплофизических параметров, ответственные за устойчивую (неустойчивую) кристаллизацию, формирование автоколебательного режима и неустойчивости жёсткого типа. Будут рассчитаны характерные размеры ячеистой микроструктуры и слоистой ликвации примеси в затвердевшем материале. В рамках модели будет проанализировано влияние интенсивности конвективного течения расплава на морфологическую устойчивость и выполнено сравнение с классической теорией устойчивости Маллинза и Секерки. В рамках проекта также будет сформулирована новая математическая модель процесса затвердевания с областью фазового превращения и граничными условиями конвективного тепломассопереноса на межфазных поверхностях. Будут построены аналитические и численные решения разработанной модели и определено влияние скорости трения на межфазной границе с расплавом на температурное и концентрационное распределения, долю твёрдой фазы в области фазового перехода, законы движения её границ и микроструктуру закристаллизовавшегося материала. Будет разработан программный комплекс для трёхмерных симуляций одновременной реализации объёмной и направленной кристаллизации, основанный на методе энтальпии-пористости для описания структурно-фазовых переходов. Программный комплекс будет учитывать наличие вынужденной и естественной конвекции расплава, движение области фазового превращения, нестационарное зарождение и рост в ней кристаллов. С помощью разработанного программного комплекса будут проведены серии практико-ориентированных симуляций кристаллизации конкретных материалов (например, Al-Li, Ti-Al, Fe-Ni, органические и неорганические материалы с фазовыми переходами – PCM или phase change materials). Результаты симуляций позволят верифицировать развитые теоретические модели и объяснить новые эмпирические данные по кристаллизации сильно переохлаждённых капель расплава на установках электромагнитной левитации.

ОТЧЁТНЫЕ МАТЕРИАЛЫ

Аннотация результатов, полученных в 2021 году
Теория. Сформулирована самосогласованная модель тепломассопереноса, описывающая процесс совместного протекания объёмного и направленного затвердевания бинарных расплавов с целью оптимизации и управления процессами кристаллизации и получения изделий с заданными свойствами и микроструктурой. Модель учитывает направленный фазовый переход, протекающий под действием температурных градиентов в жидкой и твёрдой фазах материала, которые обеспечивают движение области фазового превращения в заданном пространственном направлении, а также объёмное фазовое превращение, при котором происходит зарождение и рост кристаллитов. Модель является интегро-дифференциальной системой уравнений тепломассопереноса с подвижными границами фазового превращения: твёрдая фаза – двухфазная область и двухфазная область – жидкая фаза. Для решения сформулированной модели были разработаны приближённые аналитические методы. Используя формулу Шейла для концентрации примеси бинарной системы была определена постоянная скорость кристаллизации как функция температурных градиентов, также получено нелинейное дифференциальное уравнение для определения доли твёрдой фазы в двухфазной области как функции её переохлаждения. Из кинетического уравнения для функции распределения и граничных условий к нему было выведено второе уравнение, связывающее долю твёрдой фазы, функцию распределения кристаллов по размерам и переохлаждение. Здесь было рассмотрено два процесса роста кристаллов: с кинетическим уравнением первого и второго порядков по пространственной переменной. Уравнение первого порядка пренебрегает флуктуациями в скоростях роста кристаллов, а уравнение второго порядка учитывает такие флуктуации (отметим, что флуктуации становятся особенно важными на начальных этапах роста зарождающихся кристаллов). Для построения решений этих двух моделей было развито два аналитических подхода. Подстановка этого решения в первое выведенное уравнение позволило получить нелинейное дифференциальное уравнение, определяющее переохлаждение в двухфазной области в зависимости от её координаты. Из этого уравнения было определено переохлаждение, а его подстановка в другие уравнения позволила получить явные аналитические решения для функции распределения кристаллов по размерам и динамики снятия переохлаждения/пересыщения в области фазового перехода, законов её движения, температурного и концентрационного распределений в твёрдой, жидкой фазах и метастабильной зоне, а также доли твёрдой фазы в ней. Моделирование. Выполнено математическое описание модели энтальпии-пористости, а также анализ необходимых ограничений и допущений в модели, связанных с рассмотрением ламинарных течений и ньютоновской модели жидкости. Постановка вычислительной задачи выполнялась в терминах метода конечных объемов, также проводилась дискретизация вычислительной области и уравнений гидродинамики. Построен вычислительный алгоритм OpenFOAM для анализа физического состояния системы с учетом начальных и граничных условий в случае кондуктивного и конвективного теплопереноса. Выполнено компьютерное моделирование процессов кристаллизации с двухфазной областью, построены профили температуры и концентрации примеси в двухфазном слое, доля твёрдой фазы в нём. Эксперимент. Развитая теория кристаллизации с двухфазной областью сопоставлена с экспериментальными данными. Сопоставления были сделаны по микроструктуре и морфологии кристаллов кремния и кремниево-алюминиевого сплава, показавшие дендритный механизм роста твёрдой фазы, определён структурно-фазовый состав затвердевающих материалов в условиях объёмной и направленной кристаллизации. Конференции. Результаты проекта были доложены на 4 научных конференциях, сделано 7 докладов на профильных конференциях: Национальный Суперкомпьютерный Форум (НСКФ-2021), 30 ноября-3 декабря 2021 г., г. Переславль-Залесский, Россия; UK-Russia Workshop «Phase Transition in External Fields», 14-16 февраля 2022 г., Лондон, Великобритания; IX Международная конференция «Кристаллизация: компьютерные модели, эксперимент, технологии» (КРИС-2022), 6-9 апреля 2022 г., г. Ижевск, Россия; IX Международная молодежная научная конференция Физика. Технологии. Инновации ФТИ-2022, 16-20 мая 2022 г., г. Екатеринбург, Россия). СМИ. О результатах проекта подготовлено 5 сообщений в СМИ: 1) Портал «Научная Россия» https://scientificrussia.ru/articles/fiziki-urfu-pridumali-kak-uskorit-sozdanie-splavov-s-zadannymi-svojstvami 2) Портал «Российского научного фонда» https://rscf.ru/news/presidential-program/matematicheskaya-model-uskorit/ 3) Портал «Яндекс Дзен» https://zen.yandex.ru/media/uralscience/matematicheskaia-model-uskorit-sozdanie-splavov-s-zadannymi-svoistvami-62320f9cb519495abe123f76?& 4) Портал Уральского федерального университета https://urfu.ru/ru/news/40382/ 5) Портал Дирекции информационных технологий Уральского федерального университета https://dit.urfu.ru/ru/novosti/37478/

Публикации

1. Торопова Л.В., Александров Д.В., Реттенмайр М., Лиу Д. Microstructure and morphology of Si crystals grown in pure Si and Al-Si melts Journal of Physics: Condensed Matter, том 34, выпуск 9, номер статьи 094002 (год публикации - 2022) https://doi.org/10.1088/1361-648X/ac3792

2. Торопова Л.В., Асеев Д.Л., Осипов С.И., Иванов А.А. Mathematical modeling of bulk and directional crystallization with the moving phase transition layer Mathematical Methods in the Applied Sciences, том 45, стр. 1-11 (год публикации - 2021) https://doi.org/10.1002/mma.7864

3. Маковеева Е.В., Торопова Л.В., Александров Д.В. NucleationGibbsThomsonEffect_OpenFOAM -, 2021669295 (год публикации - )

4. Торопова Л.В., Осипов С.И., Стародумов И.О., Галенко П.К. PhaseTransitionSimulation_OpenFOAM -, 2022617806 (год публикации - )

5. - Наука и IT-технологии Официальный портал Дирекции информационных технологий Уральского федерального университета (https://dit.urfu.ru/ru/), - (год публикации - )

6. - Физики УрФУ придумали, как ускорить создание сплавов с заданными свойствами Портал «Научная Россия» (https://scientificrussia.ru/)), - (год публикации - )

7. - Математическая модель ускорит создание сплавов с заданными свойствами Портал «Российского научного фонда» (https://rscf.ru/), - (год публикации - )

8. - Математическая модель ускорит создание сплавов с заданными свойствами Портал «Яндекс Дзен» (https://zen.yandex.ru/), - (год публикации - )

9. - Математическая модель ускорит создание сплавов с заданными свойствами Портал Уральского федерального университета (https://urfu.ru/ru/), - (год публикации - )

Аннотация результатов, полученных в 2022 году
Развита теория направленной и объёмной кристаллизации, выполнено компьютерное моделирование и сделано сопоставление полученных результатов с экспериментом по затвердеванию бинарных расплавов в условиях конвективного тепломассопереноса в жидкой фазе. Теория. Построена линейная теория морфологической устойчивости плоского фронта кристаллизации с учётом конвективных движений в жидкости. Рассмотрены случаи медленной и интенсивной конвекции, для обеих моделей выведены квазистационарные распределения температуры и концентрации примеси, определена скорость движения фронта кристаллизации при наличии конвекции, сформулирована система динамических уравнений для морфологических возмущений. Получены дисперсионные соотношения, определяющие частоту возмущений как как функцию длины волны и других параметров процесса в случаях медленной и интенсивной конвекции. Построена теория направленной кристаллизации бинарного расплава с двухфазной областью в присутствии слабого конвективного течения жидкости. Построено полное аналитическое решение нелинейных уравнений двухфазной области в параметрической форме, где доля твёрдой фазы представляет собой параметр решения. Аналитически определены распределения температуры и концентрации примеси, проницаемость двухфазной области и среднее междендритное расстояние в ней, а также скорость затвердевания и толщина двухфазного слоя. Построена теория направленной кристаллизации бинарного расплава с двухфазной областью в присутствии интенсивного конвективного течения жидкости. Разработана нелинейная модель тепло- и массопереноса, описывающая сценарий квазистационарной кристаллизации с двухфазной областью. Построено точное аналитическое решение нелинейных уравнений тепломассопереноса в параметрической форме, где температура является параметром. Также аналитически найдены толщина и распределение температуры в жидкой фазе. Разработанная теория хорошо описывает эксперименты по кристаллизации раствора NaCl-H2O. Развит аналитический метод решения интегро-дифференциальной системы кинетического и балансового уравнений, описывающих эволюцию ансамбля кристаллов на промежуточной стадии процесса объёмной кристаллизации. Функция распределения частиц по размерам, а также динамика уменьшения степени метастабильности жидкости определены аналитически с помощью методов седловой точки и преобразования Лапласа. Теория позволяет получить функцию распределения кристаллов по размерам которая устанавливается в переохлажденной (пересыщенной) жидкости в начале заключительной стадии фазового превращения, когда определяющую роль играют коалесценция, коагуляция и фрагментация кристаллов. Исследовались нестационарные режимы кристаллизации, в которых происходит зарождение двухфазной области перед плоским фронтом. Построена приближённая аналитическая теория, определяющая время зарождения двухфазной области. Теория хорошо согласуется с результатами численного моделирования. Далее изучалась динамика процесса, когда перед фронтом кристаллизации располагается двухфазная область, где идёт зарождение и рост частиц твёрдой фазы. Разработана теория движения границы раздела твёрдой и жидкой фаз в переохлажденном расплаве при наличии зарождения и роста кристаллов. Сформулирована и решена интегродифференциальная система уравнений тепло- и массопереноса в двухфазном и жидком слоях, разделённых движущейся межфазной границей. Моделирование и эксперимент. Выведена функция формы вершины дендритного кристалла и проведено сопоставление теории с экспериментальными данными и численным моделированием. Выполнено сравнение с численным моделированием для тяжёлой воды и чистого никеля. Сопоставление с никелем выполнено с помощью компьютерного моделирования энтальпийным методом. Также выполнено сравнение теории с экспериментальными данными для сукцинонитрила. Показано, что теоретически найденная форма вершины дендрита с учётом конвекции и примеси находится в хорошем согласии с экспериментами и расчётами. На основе развитой теории было выполнено компьютерное моделирование обтекания вершины дендритного кристалла потоком вязкой жидкости. Показано, что механизмы тепло- и массопереноса в наклонных течениях жидкости могут быть существенно различными. Конференции. Результаты проекта были доложены на 3 научных конференциях, сделано 3 доклада на профильных конференциях: XXXI Всероссийская школа-конференция "Математическое моделирование в естественных науках", 5-8 октября 2022, г. Пермь, Россия. доклад "Зарождение и рост кристаллов на промежуточной стадии фазового перехода"; Национальный Суперкомпьютерный Форум (НСКФ-2022), 29 ноября-02 декабря 2022, г. Переславль-Залесский, Россия, доклад "Математическое и компьютерное моделирование динамики криволинейного фронта рекалесценции"; Международная молодежная научная конференция Физика. Технологии. Инновации; 15-19 мая 2023, г. Екатеринбург, Россия. Доклад " the role of fluid velocity on the shape of dendritic tips". СМИ. О результатах проекта подготовлено 8 сообщений в СМИ: 1. Портал Министерства науки и высшего образования Российской Федерации (https://www.minobrnauki.gov.ru/) - Российские физики объяснили причины аномального поведения расплавов https://www.minobrnauki.gov.ru/press-center/news/nauka/53658/ 2. Портал Уральского федерального университета (https://urfu.ru/ru/) - Физики выявили новые свойства затвердевания вещества - https://urfu.ru/ru/news/42890/ 3. Портал Поиск (https://poisknews.ru/) - Учёные придумали, как управлять свойствами сплава железа и никеля – Scientific Reports https://poisknews.ru/themes/physics/fiziki-pridumali-kak-upravlyat-svojstvami-splava-zheleza-i-nikelya-scientific-reports/?utm_source=yxnews&utm_medium=desktop 4. Портал «Российского научного фонда» (https://rscf.ru/) - Физики придумали, как управлять свойствами сплава железа и никеля https://rscf.ru/news/release/fiziki-pridumali-kak-upravlyat-svoystvami-splava-zheleza-i-nikelya-/ 5. Портал Научная Россия (https://scientificrussia.ru/) - Физики урфу рассчитали, как управлять свойствами сплава железа и никеля https://scientificrussia.ru/articles/fiziki-urfu-rasscitali-kak-upravlat-svojstvami-splava-zeleza-i-nikela 6. Портал Уральского федерального университета (https://urfu.ru/ru/) Физики предложили теорию затвердевания сплавов никеля и железа https://urfu.ru/ru/news/44155/ 7. Nano News Net - on-line издание, посвященное вопросам наноиндустрии (https://www.nanonewsnet.ru/) - Физики УрФУ рассчитали, как управлять свойствами сплава железа и никеля https://www.nanonewsnet.ru/news/2022/fiziki-urfu-rasschitali-kak-upravlyat-svoistvami-splava-zheleza-nikelya 8. ИНВУР – Инновационный портал Уральского федерального округа (http://www.invur.ru/) - Физики УрФУ рассчитали, как управлять свойствами сплава железа и никеля http://www.invur.ru/index.php?page=news&id=132061

Публикации

1. Александров Д.В., Торопова Л.В. The role of incoming flow on crystallization of undercooled liquids with a two-phase layer Scientific Reports, 12, 17857 (год публикации - 2022) https://doi.org/10.1038/s41598-022-22786-w

2. Е.В. Павлюк, Д.В. Александров, Н.В. Кропотин, Л.В. Торопова, И.О. Стародумов Моделирование процесса направленного затвердевания/плавления методом энтальпии-пористости Расплавы, 2023, № 3, с. 250-266 (год публикации - 2023)

3. Торопова Л.В., Александров Д.В. Dynamical law of the phase interface motion in the presence of crystals nucleation Scientific Reports, 12, 10997 (год публикации - 2022) https://doi.org/10.1038/s41598-022-15137-2

4. Торопова Л.В., Маковеева Е.В., Осипов С.И., Малыгин А.П., Йанг Й., Александров Д.В. Nucleation and Growth of an Ensemble of Crystals during the Intermediate Stage of a Phase Transition in Metastable Liquids Crystals, 12 (7), 895 (год публикации - 2022) https://doi.org/10.3390/cryst12070895

5. Галенко П.К., Кропотин Н.В., Стародумов И.О., Торопова Л.В. ImpurityDistribution_OpenFOAM -, 2022669416 (год публикации - )

6. Павлюк Е.В., Стародумов И.О., Галенко П.К., Кропотин Н.В., Торопова Л.В. PostPr_OF -, 2023619165 (год публикации - )

7. - Физики УрФУ рассчитали, как управлять свойствами сплава железа и никеля Nano News Net - on-line издание, посвященное вопросам наноиндустрии (https://www.nanonewsnet.ru/), - (год публикации - )

8. - Физики УрФУ рассчитали, как управлять свойствами сплава железа и никеля ИНВУР – Инновационный портал Уральского федерального округа (http://www.invur.ru/), - (год публикации - )

9. - Российские физики объяснили причины аномального поведения расплавов Портал Министерства науки и высшего образования Российской Федерации (https://www.minobrnauki.gov.ru/), - (год публикации - )

10. - Физики выявили новые свойства затвердевания вещества Портал Уральского федерального университета (https://urfu.ru/ru/), - (год публикации - )

11. - Физики придумали, как управлять свойствами сплава железа и никеля Портал «Российского научного фонда» (https://rscf.ru/), - (год публикации - )

12. - Физики предложили теорию затвердевания сплавов никеля и железа Портал Уральского федерального университета (https://urfu.ru/ru/), - (год публикации - )

13. - Физики УрФУ рассчитали, как управлять свойствами сплава железа и никеля Портал Научная Россия (https://scientificrussia.ru/), - (год публикации - )

14. - Учёные придумали, как управлять свойствами сплава железа и никеля – Scientific Reports Портал Поиск (https://poisknews.ru/), - (год публикации - )