Роль физико-химических процессов при численном моделировании задач входа спускаемых космических аппаратов в атмосферу Земли и других планет.

Ответственный за направление д.ф.-м.н. Сахаров В.И.

Несмотря на значительное количество полетов высокоскоростных летательных аппаратов в атмосфере Земли, интерес к задачам газовой динамики, возникающим при их разработке, не ослабевает. Это обусловлено как фундаментальными газодинамическими аспектами рассматриваемой проблемы, так и практическими приложениями, связанными с необходимостью расчетов параметров аэродинамики и теплообмена при проектировании перспективных летательных аппаратов с целью сокращения количества дорогостоящих экспериментов в аэродинамических установках.

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

Молекулярное строение газа не рассматривается в классической газовой динамике, а величины, связанные со свойствами молекул (теплоемкости, коэффициенты переноса и др.), входят в уравнения газовой динамики в виде различных коэффициентов, для определения которых можно не привлекать теорию молекулярного строения. Высокие температуры исключают возможность полного описания задач аэродинамики в рамках адиабатического течения совершенного газа с постоянными теплоемкостями. Если в газе происходят физико-химические превращения: возбуждение колебательных и вращательных степеней свободы молекул, диссоциация (рекомбинация) молекул, ионизация (нейтрализация), возбуждение электронных уровней атомов и молекул, излучение и поглощение лучистой энергии, необходимо использовать представления о свойствах частиц, из которых состоит газ, рассматривать их внутреннее строение и взаимодействие. Все эти процессы в большей степени влияют на теплопередачу к обтекаемому таким высокотемпературным газом телу и меньше на сопротивление.

Характерной особенностью течений высокотемпературного многоатомного газа являются процессы диссоциации (рекомбинации) и ионизации (дезактивации). Выделение механизма реакций, играющую основную роль в рассматриваемом диапазоне температур важно для правильного расчета параметров в газе и сокращения объема необходимых вычислений.

Важным обстоятельством для численного моделирования задач входа космических тел в атмосферы Земли и других планет является выбор кинетических моделей процессов, протекающих в ударно нагретом газе около тела. Так, анализ кинетических процессов за фронтом головной ударной волны в воздухе при скоростях движения 6–11 км/с, и рассмотрение, наряду с процессами диссоциации, и процессов ионизации позволил правильно оценить длину релаксационной зоны ионизации за фронтом волны (работы С.А. Лосева, В.А. Полянского в НИИ Механики МГУ).

Существенное различие в константах скоростей химических реакций и реакций ионизации в разных моделях может оказывать влияние на макроскопические параметры в расчетах течений высокотемпературного газа. В настоящее время существуют многочисленные кинетические модели реакций в воздухе, в смесях, содержащих CO2, и в других газах.

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

На важность учета неравновесного возбуждения электронных состояний атомарных компонентов при рассмотрении задач обтекания тел со скоростями > 10 км/с указано в работах В.А.Горелова и др. (ЦАГИ), где предложена упрощенная ступенчатая модель ионизации атомов N и O, учитывающая резонансный обмен между электронными состояниями уровней этих компонентов. Применение модели в расчетах приводит к лучшему совпадению с экспериментом температуры электронов за фронтом волны.

При разработке программных комплексов, предназначенных для проведения расчетов гиперзвуковых течений газа в рамках уравнений Эйлера, Навье?Стокса и Рейнольдса, необходимо опираться на специализированные базы данных по термодинамическим, кинетическим и транспортным свойствам индивидуальных газов и газовых смесей. Используя модели газовой среды различной сложности от совершенного газа до многокомпонентной термически и химически неравновесной газовой среды, можно создавать для решения конкретных задач оптимальные рабочие программы, сочетающие достаточно полное описание рассматриваемого течения с минимальными затратами ресурсов вычислительных средств. Такие программные комплексы могут быть сориентированы на решение ряда актуальных задач внешней и внутренней аэротермодинамики и открыты для разработчиков в плане внедрения в него математического описания и программной реализации процессов, учет которых необходим для решения конкретных задач.