РАСЧЕТ ТЕРМОДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК СИСТЕМЫ Fе – P МЕТОДОМ МОЛЕКУЛЯРНОЙ ДИНАМИКИ

Проблема дефосфоризации железоуглеродистых сплавов является актуальной для металлургической отрасли, так как высокая концентрация фосфора способствует появлению целого ряда крайне нежелательных явлений. Решению этой задачи посвящено множество экспериментальных работ, но справиться с ней окончательно до сих пор не удалось. Любые натурные эксперименты, направленные на исследование процесса удаления фосфора, требуют немалых материальных и временных затрат, но при этом не гарантируют получение нужного результата. Поэтому для поиска новых подходов к разрешению обозначенной проблемы гораздо рациональнее использование методов численного моделирования с привлечением вычислительных возможностей современных компьютеров. В настоящее время компьютерные эксперименты являются таким же признанным методом исследования, как теоретическое исследование и реальный эксперимент. Для изучения поведения атомов фосфора в железе с помощью численного эксперимента необходимо построить расчетную модель и апробировать ее путем вычисления различных характеристик, значения которых заранее известны. В качестве метода компьютерного моделирования в данной работе был выбран метод молекулярной динамики. С его помощью можно проводить эксперименты с заданными скоростями атомов и описывать динамику исследуемых процессов. Для описания межчастичного взаимодействия использовался потенциал, рассчитанный в рамках метода погруженного атома. Исследование проводилось на расчетной ячейке, имитирующей кристалл α-железа с атомами замещения фосфора. Построенная модель продемонстрировала удовлетворительные результаты при вычислении известных характеристик моделируемой системы. Были установлены зависимости изменения таких характеристик, как температурный коэффициент линейного расширения, температура плавления, скрытая теплота плавления и теплоемкость от концентрации атомов замещения фосфора, а также в отдельных случаях от величины приложенного внешнего давления. Расчеты показали, что, например, концентрация фосфора 0,5 % приводит к росту среднего теплового коэффициента линейного расширения на 9 %, снижению температуры и скрытой теплоты плавления на 5 %, а теплоемкости на 7 %. 

1. Дауд А.Д., Семин А.Е., Котельников Г.И., Щукина Л.Е. Дефосфоризация хромистых расплавов с использованием оксидов редкоземельных металлов // Изв. вуз. Черная металлургия. 2017. Т. 60. № 1. С. 54 – 59.

2. Робей Р., Уайтхед М. Внедоменная обработка чугуна с учетом конкретных производственных условий // МРТ. Металлургическое производство и технология металлургических процессов. 2014. № 1. С. 16 – 24.

3. Георгадзе А.Г., Гернер В.И., Елашвили М.И. и др. Условия дефосфоризаци жидкого металла в разливочном ковше // Литье и металлургия. 2012. № 3(67). С. 117 – 119.

4. Starostenkov M.D., Potekaev A.I., Grinkevich L.S. etc. Dynamics of edge dislocations in a low-stability fcc-system irradiated by highenergy particles // Russian Physics Journal. 2017. Vol. 59. No. 9. P. 1446 – 1453.

5. Markidonov A.V., Starostenkov M.D., Smirnova M.V. Selfdiffusion process in an FCC crystal caused by the passage of a shock wave // Russian Physics Journal. 2015. Vol. 58. No. 6. P. 828 – 832.

6. Markidonov A.V., Starostenkov M.D., Poletaev G.M. Transformation of nanopores in gold under conditions of thermoactivation and the effects of acoustic and shock waves // Bulletin of the Russian Academy of Sciences: Physics. 2015. Vol. 79. No. 9. P. 1089 – 1092.

7. Markidonov A.V., Starostenkov M.D., Tabakov P.Y. Splitting vacancy voids in the grain boundary region by a post-cascade shock wave // Materials Physics and Mechanics. 2013. Vol. 18. No. 2. P. 148 – 155.

8. Andersen H.C. Molecular dynamics simulations at constant pressure and/or temperature // The Journal of Chemical Physics. 1980. Vol. 72. No. 4. P. 2384 – 2393.

9. Daw M.S., Baskes M.I. Embedded-atom method: Derivation and application to impurities and other defects in metals // Physical Review B. 1984. Vol. 29. No. 12. P. 6443 – 6453.

10. Foiles S.M., Baskes M.I., Daw M.S. Embedded-atom-method functions for the FCC metals Cu, Ag, Au, Ni, Pd, Pt, and their alloys // Phys Rev B. 1986. Vol. 33. No. 12. P. 7983 – 7991.

11. Biersack J.P., Ziegler J.F. Refined universal potentials in atomic collisions // Nuclear Instruments and Methods. 1982. Vol. 194. P. 93 – 100.

12. Ackland G.J., Mendelev M.I., Srolovitz D.J. etc. Development of an interatomic potential for phosphorus impurities in α-iron // Journal of Physics: Condensed Matter. 2004. Vol. 16. P. 2629 – 2642.

13. Ko W.-S., Kim N.-J., Lee B.-J. Atomistic modeling of an impurity element and a metal–impurity system: pure P and Fe–P system // Journal of Physics: Condensed Matter. 2012. Vol. 14. P. 225002 – 225016.

14. Григорьев И.С., Мейлихов Е.З. Физические величины: Справочник. – М.: Энергоатомиздат, 1991. – 1232 с.

15. Belonoshko A.B., Skorodumova N.V., Rosengren A., Johansson B. Melting and critical superheating // Physical Review B. 2006. Vol. 73. P. 0122011 – 0122013.

16. Mazhukin V.I., Shapranov A.V., Perezhigin V.E. etc. Kinetic melting and crystallization stages of strongly superheated and supercooled metals // Mathematical Models and Computer Simulations. 2017. Vol. 9. No. 4. P. 448 – 456.

17. Yang H., Lu Y., Chen M., Guo Z. A molecular dynamics study on melting point and specific heat of Ni3Al alloy // Science in China. Series G: Physics, Mechanics & Astronomy. 2007. Vol. 50. No. 4. P. 407 – 413.

18. Таблица физических величин: Cправочник / Под ред. И.К. Кикоина. – М.: Атомиздат, 1976. – 1008 с.

19. Зайцев И.Д., Зозуля А.Ф., Асеев Г.Г. Машинный расчет физико-химических параметров неорганических веществ. – М.: Химия, 1983. – 256 с.

20. Зиновьев В.Е. Теплофизические свойства металлов при высоких температурах: Справочник. – М.: Металлургия, 1989. – 384 с.