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