Предложена численная модель обтекания тел плазмой в рамках двухжидкостной модели магнитной гидродинамики. Разработаны метод учета граничных условий посредством перегрузки оператора индексации, численная модель течения бесстолкновительной плазмы. Оценен эффект нарушения непрерывности течения газа в задаче гиперзвукового обтекания тела.