Аннотация:Рассмотрены проблемы численного моделирования термодинамически согласованной моделидинамики параметра порядка для фазовых превращений первого рода в рамках фазово-полевого подхода. Модель основана на расширенной необратимой термодинамике, использующей для описания состояния системы потенциалы Гиббса реальных систем. Численные алгоритмы для моделирования уравнений фазового поля как для локально-равновесных, так и локально-неравновесных моделей релаксационной динамики получены из условия гра диентной устойчивости. Выполненные оценки позволяют построить градиентно-устойчивый алгоритм для моделирования уравнений фазового поля в реальных термодинамических системах. В работе проведён анализ градиентно-устойчивых численных схем на устойчивость по малым отклонениям для уравнений фазового поля в Фурье-приближении. Выполненные численные расчеты для раствора Si-As демонстрируют устойчивость полученного алгоритма как при увеличении шага по времени на несколько порядков, так и при росте движущих сил из-за сильного переохлаждения системы.