Проведено детальное моделирование тепловой эволюции эпизодически\r\nаккрецирующих нейтронных звёзд для разных теоретических моделей ядра,\r\nвнутренней коры и внешних оболочек нейтронной звезды с учётом глубокого\r\nпрогрева коры при разных предположениях об истории её формирования.\r\nРассмотрена как средняя долговременная эволюция, так и сравнительно\r\nкратковременная тепловая эволюция нейтронной звезды в течение активной\r\nфазы аккреции и после её завершения. Тепловая эволюция впервые\r\nпромоделирована согласованно с изменениями состава коры, принимая во\r\nвнимание, что исходная неаккрецированная кора постепенно замещается\r\nпереработанным аккрецированным материалом. Для неаккрецированной части\r\nкоры проведено сравнение канонической модели вещества в основном\r\nсостоянии и новой модели, учитывающей конечную температуру отжига коры. \r\nПроведено сравнение результатов моделирования между собой и с\r\nрезультатами, полученными ранее в рамках упрощённых моделей, а также с\r\nнаблюдениями. Показано, что долговременная тепловая эволюция\r\nэпизодически аккрецирующих нейтронных звёзд немонотонна, а\r\nкратковременная эволюция зависит не только от массы нейтронной звезды,\r\nтемпа акреции на неё и микрофизических моделей вещества в её коре и\r\nядре, но также от возраста и сценария формирования коры. Поэтому\r\nдетальные самосогласованные расчёты тепловой эволюции эпизодически\r\nаккрецирующих нейтронных звёзд можно использовать для тестирования таких\r\nмоделей и сценариев. Работа поддержана грантом РНФ 14-12-00316.