Недавно ставшие возможными наблюдения холодных и слабых белых карликов на телескопах GAIA и др. дают богатый материал для уточнения моделирования финальной стадии их эволюции, что в свою очередь может использоваться для изучения поведения вещества при экстремальных условиях. В настоящей работе было проведено моделирование тепловой эволюции белых карликов для разных составов атмосфер с использованием гидродинамического кода Блинникова С.И. и Дуниной-Барковской Н.В., позволяющего изучать эволюцию объекта до времен порядка возраста Вселенной. Новым в этом коде было введение уравнений состояния для неидеальной плотной плазмы, разработанных Потехиным А.Ю.. В моделировании использовалось приближение однокомпонентной плазмы, состоящей из углерода, либо кислорода. Ранее использовавшееся уравнение состояния Надёжина Д.К. с кулоновскими поправками И.Ибена для идеальной плазмы было рассчитано на горячие белые карлики и не учитывало эффектов, наступающих на поздних стадиях их эволюции, например, эффекта кристаллизации. Для расширения кода на область холодных и старых белых карликов были необходимы уравнения, включающие эффекты неидеальности в плотной плазме. Как и ожидалось, после наступления фазы кристаллизации, происходит небольшой выброс тепла, а затем скорость остывания увеличивается. Особенно этот эффект существенен для массивных белых карликов (масса ∼ 1.3 M⊙), для которых время остывания до Teff ∼ 3 · 10^3 K уменьшается почти в 10 раз.