Настоящий обзор посвящен теоретическим исследованиям пространственных задач обтекания гладких тел сверхи гиперзвуковым потоком вязкого, теплопроводного газа в широком диапазоне изменения чисел Рейнольдса от малых, для которых применима модель сплошной среды, и до умеренно больших, когда поток уже можно разделить на невязкое течение и пограничный слой около обтекаемого тела.
Пространственный вязкий ударный слой.
Уравнения пограничного слоя качественно и количественно правильно описывают течение в тонком слое около поверхности тела при больших числах Рейнольдса. Однако в верхних слоях атмосферы, когда характерные числа Рейнольдса становятся Re^ < 103, классическая модель разбиения возмущенной области течения на невязкое течение и пограничный слой становится неприменимой. Поэтому для исследования течения при умеренных числах Рейнольдса широкое распространение получила теория вязкого ударного слоя. Основное преимущество этой теории состоит в пригодности соответствующих уравнений во всей возмущенной области течения от ударной волны до поверхности тела. В зависимости от числа Маха набегающего потока различают модель гиперзвукового или тонкого вязкого ударного слоя (ТВУС) (при Мл «1) и модель полного вязкого ударного слоя (ПВУС).
Гиперзвуковой вязкий ударный слой.
Впервые модель ТВУС для расчета гиперзвукового обтекания плоских и осесимметричных тел при умеренно больших числах Рейнольдса предложил в 60-е годы Н. К. Cheng. В предположении у 1, МХ «I, Re0 = p^V^Rj-> К = ?» Re0 =0(1) вся область течения разбивалась на два подслоя: вязкий ударный слой, примыкающий к телу и зона перехода через скачок уплотнения. Поскольку уравнения, описывающие течение в области перехода через скачок, допускают однократное интегрирование, то полученные таким образом соотношения, называемые обобщенными соотношениями Рэнкина-Гкмонио, можно использовать в качестве граничных условий на внешней границе вязкого ударного слоя. Поэтому задачи о течении в собственно ударном слое и в структуре ударной волны можно решать отдельно друг от друга. В отличие от классических условий Рэнкина-Гюгонио обобщенные условия учитывают влияние эффектов молекулярного переноса в области непосредственно за головным скачком уплотнения и тем самым влияют на течение во всем ударном слое.
Полный вязкий ударный слой.
Полные уравнения вязкого ударного слоя являются естественным обобщением модели ТВУС. Они содержат все члены, вносящие вклад во второе приближение пограничного слоя, включающее эффекты вытеснения пограничного слоя, вихревого взаимодействия, влияние продольной и поперечной кривизн, эффект скольжения и скачка температуры на обтекаемой поверхности, и все члены уравнений Эйлера. Уравнения рассматриваются во всем ударном слое от тела до ударной волны, которая является внешней границей области течения. Впервые такой подход для двумерных течений совершенного газа предложил, А И Толстых (см. О М Белоцерковский и ДР- [6]).
В качестве граничных условий к системе уравнений ПВУС на ударной волне используются обобщенные условия Рэнкина-Гюгонио, а на поверхности тела для тангенциальных компонент вектора скорости ставятся условия прилипания, для температуры — условия, соответствующие охлажденной или теплоизолированной поверхности. В случае малых и умеренных чисел Рейнольдса граничные условия видоизменяются с учетом эффектов скольжения и скачка температуры.
Сопоставление полных уравнений Навье-Стокса с системой уравнений ПВУС показывает, что в последних опущены члены, имеющие порядок.
O (Re0)-1 и выше. При К =.
Важным достоинством модели ПВУС является то, что в ней можно не предполагать, что ударный слой тонок, т. е. не требовать, чтобы параметр е был мал. Отсюда следует возможность применения рассматриваемой модели для расчета обтекания тел при небольших сверхзвуковых скоростях набегающего потока, включая режимы, для которых ударный слой не является тонким на всей наветренной части обтекаемого тела.
В большинстве задач обтекания в потоке имеются области как дозвуковых, так и сверхзвуковых скоростей. Например, при обтекании со сверхзвуковой скоростью тел имеющих затупленную носовую часть Это усложняет решение задачи..
Одним из первых методов решения уравнений ПВУС, позволивший преодолеть трудности, связанные с дозвуковым течением газа в некоторых областях ударного слоя, был метод глобальных итераций по всему полю течения, предложенный для двумерных задач R. Т. Davis’oM [53]. В рамках этого метода форма ударной волны и эллиптические члены в уравнениях, учитывающие распространение вверх по потоку, считаются известными из предыдущей итерации, причем для расчета первого приближения используется модель ТВУС. Однако попытки использовать метод глобальных итераций в этом виде для сферически затупленных конусов с малыми углами полураствора встретили ряд затруднений, т.к. модель ТВУС дает отрыв на подветренной стороне тела..
Для преодоления этих трудностей было предложено изменить способ задания начальных данных, а также применить метод релаксации при определении формы ударной волны. Кроме того, было введено сглаживание угла наклона ударной волны, а в окрестности разрыва кривизны поверхности тела использована специальная разностная запись продольных производных, учитывающая условие непрерывности физических составляющих вектора вдоль контура тела. Для получения устойчивого решения при расчете течения около тел, где ударный слой не является тонким, важным было также предложение совместного решения системы уравнений первого порядка (уравнение неразрывности и уравнение импульсов в проекции на нормаль) путем векторной прогонки..
Первые работы по применению метода глобальных итераций к решению трехмерных уравнений ПВУС выполнили A. L. Murray и С. Н. Lewis [63]. Был проведен расчет обтекания однородным газом сферически затупленных конусов с малыми углами полураствора под углами атаки а= 04−38°. В этих работах численное решение начиналось на сферическом затуплении в осесимметричной постановке Результаты, полученные на границе осесимметричной области, использовались в дальнейшем как начальные данные для решения трехмерной задачи, начинающейся от плоскости симметрии на наветренной стороне тела При этом необходимое для замыкания задачи значение градиента давления в окружном направлении аппроксимировалось по расчетным данным из таблиц невязкого сверхзвукового обтекания этого же тела. По завершению цикла по окружной координате делался шаг по маршевой координате вниз по потоку. В каждой расчетной области уравнения решались последовательно в следующем порядке, уравнение импульсов в проекции на окружную координату, уравнение энергии, уравнение импульсов в проекции на маршевую координату, уравнение неразрывности для определения отхода ударной волны, совместное решение уравнения неразрывности и уравнения импульсов в проекции на нормаль. На основании проведенных расчетов сделан вывод, что уравнения ПВУС позволяют с хорошей точностью получить полную картину обтекания при малых углах атаки. При больших углах атаки модель вязкого ударного слоя дает удовлетворительные результаты на наветренной стороне и непригодна для исследования течения на подветренной стороне тела..
В дальнейшем указанный метод был распространен на решение трехмерных уравнений ПВУС с учетом протекающих в потоке различных физико-химических процессов. S. Swaminathan, М. D. Кип и С. Н. Lewis [69] рассмотрели неравновесное течение диссоциированного воздуха около затупленных по сфере конусов, обтекаемых под небольшими (до 10°) углами атаки. Режимы полета соответствуют высотам 70.104, 83.82 км и числу Маха 27.18. В качестве химической модели применяется семикомпонентная гаювая смесь. Граничные условия на теле — полностью каталитическая поверхность Диффузионная модель ограничена бинарной диффузией. В работе [50] показано, что концентрация электронов слабо зависит от степени упрощения диффузионной модели Рассматривается влияние различных типов граничных условий (прилипания, скольжения) на концентрацию химических компонентов смеси. Сравнение показывает, что результаты расчетов с граничными условиями прилипания ближе к экспериментальным данным..
Другой метод, нашедший наряду с методом глобальных итераций достаточно широкое применение для решения системы уравнений ПВУС — это метод установления. В этом методе в качестве исходных используются нестационарные уравнения ПВУС, а искомое решение получается как предел установления по времени. В этом случае вначале установлением решалась задача в области затупления, а затем, также установлением по времени, находилось решение в узких перекрывающихся областях, последовательно перемещающихся вдоль тела. Возникающие при этом интенсивные осцилляции параметров поперек ударного слоя сглаживались методом четвертого порядка. Отмечено, что если для определения формы ударной волны и распределения давления по телу требуется порядка 103 шагов по времени, то для получения стационарного значения теплового потока требуется гораздо большее число шагов [24]..
Решение параболизованных уравнений Навье-Стокса.
Остановимся теперь на работах, использующих для решения рассматриваемых в обзоре задач, параболизованные уравнения Навье-Стокса (ПУНС). Несмотря на значительные вычислительные трудности, связанные с созданием устойчивых численных алгоритмов для этой модели, начиная с середины 70-х годов известно большое число работ (особенно зарубежных), в которых сверхзвуковое обтекание тел потоком вязкого газа изучается в рамках параболизованных уравнений Навье-Стокса. Сравнение проведенных в ее рамках расчетов с экспериментом позволяет считать, что для решения многих практически важных задач можно с успехом использовать такие параболизованные модели течения. Можно отметить большую эффективность модели параболизованных уравнений Навье-Стокса для решений струйных задач, течения в следе..
Система ПУНС получается из исходной системы уравнений Навье-Стокса пренебрежением эффектом молекулярного переноса в направлении основного течения газа в ударном слое (маршевом направлении) При этом возможны два случая. В первом из них данные эффекты учитываются лишь в направлении, перпендикулярном по отношению к поверхности тела (система уравнений ПУНС в приближении тонкого слоя), а во втором производится учет эффектов молекулярного переноса и в окружном направлении. В последнем случае система ПУНС позволяет изучать течения, в которых могут иметься возвратные в окружном направлении течения. Это становится особенно важным при расчете обтекания тел под большими углами атаки, когда, как показывает опыт, на подветренной стороне тела возникает протяженная зона возвратных поперечных течений, простирающаяся в окружном, направлении..
С использованием ПУНС в работе авторами J. D. Waskiewicz и С. Н Lewis [73] было рассмотрено обтекание сферически затупленных конусов с углом полураствора 6 = 1° и 10° при параметрах набегающего потока, а = 10°, 23°, 38°, Л^. =10.17, 18,22.7 и Re, = 104* 106. Вся исходная область течения разбивалась на две В первой из них, расположенной в окрестности затупления, где скорость газа в маршевом направлении была дозвуковой, решение находилось в рамках модели ПВУС методом, который предложили в работе [74] J. D. Waskiewicz, A. L. Murray и С. Н. Lewis Метод базируется на схеме с последовательным решением уравнений, с тем отличием, что дифференциальные уравнения первого порядка (уравнения неразрывности и движения в проекции на нормаль) решаются совместно Отмечено, что предложенная модификация существенно повышает устойчивость схемы и качество полученного решения В сверхзвуковой области решение продолжалось в рамках модели ПУНС методом S С Lubard’a и W. S. HelliweFa. Поскольку в процессе решения головная ударная волна выделялась, то для замыкания задачи и определения формы скачка использовалось уравнение неразрывности, в котором производные по нормали к скачку аппроксимировались односторонними разностями. В работе подробно исследован вопрос о влиянии различных способов определения окружного градиента давления на основные характеристики течения и размеры шага в маршевом направлении. На основании сравнения с экспериментом сделан вывод, что применение указанного комбинированного метода расчета позволяет удовлетворительно рассчитать структуру течения вязкого газа около сферически затупленных конусов под большими углами атаки при сверхи гиперзвуковых скоростях потока..
Как отмечают R. Т. Davis и S. G. Rubin, при расчете ПУНС точность получаемого решения в сильной степени зависит от правильного выбора маршевого направления Это становится особенно важным при больших углах атаки, когда линии тока достаточно сильно отклоняются от направлений, задаваемых при обычно используемых системах координат, связанных с телом. В связи с этим М. D. Kim, R. R Thareja и С Н Lewis [59] предложили выбирать в качестве маршевого направления линии тока внешнего невязкого течения на поверхности тела. Сравнение результатов расчета для умеренных углов атаки хорошо согласуется с решениями полученными ранее..
В работе [49] авторами [Bhutta, Lewis] разработан метод решения параболизованных уравнений Навье-Стокса, позволяющий рассчитывать трехмерные неравновесные гиперзвуковые течения при больших углах атаки около конусов с изломами образующей. Используется расщепление задачи о течении неравновесной газовой смеси на две: решение газодинамической задачи и расчет химических превращений В результате поочередного решения итерационным методом этих задач учитывается взаимное влияние газодинамических и химических процессов. Интегрирование уравнений производится маршевым методом при помощи схемы, построенной по типу «предиктор-корректор». Особое внимание уделяется адекватному учету сильного взаимодействия продольного течения с поперечным потоком в отрывной области на подветренной стороне тела. В качестве примера использования настоящего метода решения ПУНС проведены расчеты трехмерного обтекания затупленного по сфере конуса при числе Маха, равном 20, и угле атаки соответствующем 20 градусам Причем для оценки точности данного метода и его эффективности при больших углах атаки проведено сопоставление решений, полученных на трех различных сетках. Анализ результатов расчета трехмерного обтекания тела неравновесным гиперзвуковым потоком показывает, что при больших углах атаки диссоциация кислорода и азота, а также образование И (У протекают гораздо интенсивнее с подветренной стороны тела, чем с наветренной. Более сильная диссоциация газа по мнению авторов обусловлена тем, что на подветренной стороне тела формируется обширная область высокотемпературного газа. В этой области на подветренной стороне тела будет наблюдаться значительное увеличение концентрации свободных электронов. Сравнение результатов расчетов на сетках с различным числом узлов показывает, что измельчение сетки сопровождается вполне приемлемым (до 2-х раз) увеличением затрат машинного времени. Кроме того, сопоставление этих результатов позволяет заключить, что для достижения необходимой точности расчета распределений плотности теплового потока и коэффициента поверхностного трения, используемая в расчетах сетка должна иметь достаточное сгущение узлов в области отрыва поперечного потока Причем такое сгущение необходимо как по нормали к оси тела, так и в поперечном к потоку направлении. Разработанный метод решения трехмерных уравнений, алгоритм которого предусматривает нахождение решения на каждом маршевом шаге при помощи итерационной процедуры (т.е. методом установления). Метод позволяет рассчитывать пространственные течения с высокой точностью при относительно малых затратах машинного времени Это обсюятельство дает возможность вести интегрирование уравнений с большим шагом по маршевому направлению, что значительно сокращает суммарное время счета..
Наряду с анализом трехмерных течений с использованием точных моделей, основанных на уравнениях Навье-Стокса, достаточно широкое распространение получили различные приближенные подходы в расчете уравнений пространственного пограничного слоя, позволяющие в ряде случаев получить достаточно точное для инженерных приложений решение исходной задачи..
Метод осесимметричной аналогии..
В основе метода лежит предположение о том, что в системе координат, связанной с линиями тока внешнего течения (jc1 = const — линии тока невязкою течения на поверхности тела, х2 = const — их ортогональные траектории), скорость и2, поперечная линиям тока внешнего течения и ее производные малы по сравнению со скоростью вдоль линии тока Отметим, что поскольку поперечная скорость на поверхности тела и на внешней границе пограничного слоя равна нулю, предположение о малой интенсивности вторичного течения достаточно обосновано. Этот метод использует поле рассчитанного невязкого трехмерного течения, необходимое для вычисления поверхностных линий тока, метрических коэффициентов и граничных условий пограничного слоя..
В работе [55] рассмотрено обтекание сферически затупленного конуса с углом полураствора 15°. Приведено сравнение с экспериментальными данными. Режим обтекания соответствует, а =0°, 20°, Mv = 10.6, Ree = I.3M06 [1/м], радиус затупления равен 0.95 и 2.79 см. Отмечается хорошее согласование отношений тепловых потоков с экспериментальными данными как на наветренной, так и на подветренной сторонах тела..
Однако в литературе достаточно часто встречаются упоминания о хорошем согласовании приближенных подходов на наветренной стороне тела и о значительных погрешностях результатов, когда возникает отрыв потока..
Заключение.
Уравнения вязкого ударного слоя (ВУС) описывают течение во всем ударном слое с помощью единой композитной системы уравнений, одинаково пригодной как в невязкой, так и в вязкой областях течения Тем самым исключается проблема асимптотического сращивания решений уравнений пограничного слоя и уравнений Эйлера при не очень больших числах Рейнольдса (Re^ < 103). Эта система уравнений обычно решается методом установления в дозвуковых областях течения и маршевым методом в сверхзвуковых с той или иной аппроксимацией маршевого градиента давления в дозвуковом приповерхностном слое. В последнее время для решения двумерных уравнений ВУС развит метод глобальных итераций, позволяющий эффективно рассчитывать течение во всем ударном слое с наличием дои сверхзвуковых областей единообразным способом в широком диапазоне параметров без привлечения каких-либо других расчетных данных Развитие метода глобальных итераций на трехмерные задачи позволит значительно сократить затраты времени на решение и даст возможность практически осуществить серийные расчеты задач сверхзвукового обтекания в широком диапазоне чисел Рейнольдса и Маха.
Наиболее популярной моделью для решения подавляющего числа задач сверхи гиперзвукового обтекания, особенно у американских авторов, является модель ПУНС, для решения которой в настоящее время разработан ряд эффективных численных методов. Эта модель оказалась достаточно точной, сравнительно просто поддающейся эффективному численному решению и экономичной по затратам машинного времени по сравнению с решением и полных уравнений Навье-Стокса. В связи с чем она стала весьма полезной в инженерных приложениях.
Разработка численных методов остается в будущем в центре внимания математиков-вычислителей и специалистов по вычислительной аэродинамике. В процессе естественного отбора наибольшее распространение получат методы, дающие более точные результаты, обладающие высокой скоростью сходимости, легко программируемые и надежные, даже в ущерб их универсальности..
Необходимость разработки нового программного комплекса была обусловлена тем, что имеющиеся программы не позволяют надежно рассчитывать режимы течения с учетом физико-химических превращений в ударном слое, соответствующие высотам менее 40 км и числам Маха более 10 в набегающем потоке. Кроме того, применение известных программ для параметрических расчетов трехмерных течений в инженерной практике неэффективно в силу неоправданно больших временных затрат..
Заключение.
Применение моделей, основанных на уравнениях Навье-Стокса с учетом неравновесных физико-химических превращений в потоке, позволяет пересчитывать экспериментальные данные на реальные условия полёта гиперзвуковых летательных аппаратов, а также сократить количество дорогостоящих экспериментов по определению аэродинамических, тепловых и прочностных характеристик изделий..
Анализ известных численных методов решения задач сверхзвукового пространственного обтекания тел потоком воздуха позволяет сделать вывод, что наиболее удобными и универсальными являются конечно-разностные методы сеток. Среди них можно выделить методы «сквозного счета», позволяющие рассчитывать с приемлемой точностью поле течения около реальных летательных аппаратов без выделения неизвестных до решения задачи внутренних поверхностей разрыва газодинамических параметров.
На базе метода установления нестационарных уравнений газовой динамики и химической кинетики по времени с использованием принципа разделения задачи обтекания по физическим процессам разработан эффективный численный метод решения полных (параболизованных) уравнений Навье-Стокса с учетом неравновесных физико-химических процессов. Использование криволинейной системы координат, локально-близкой к линиям тока течения, позволяет адекватно учитывать поведение параметров в областях больших градиентов и расширить применимость метода для расчётов обтекания тел под большими углами атаки. При этом используемая консервативная форма записи уравнений газовой динамики позволяет наиболее точно удовлетворять разностным аналогам интегральных законов сохранения. Кроме того, для расчета обтекания удлиненных 1ел используется разделение всей области интегрирования на ряд взаимно перекрывающихся подобластей и последовательный расчет в каждой из них. Данное разделение области возможно в силу слабой передачи возмущений вверх по потоку при обтекании тел сверхзвуковым набегающим потоком вязкого газа..
Качественное и количественное согласование параметров неравновесных течений с экспериментальными и теоретическими данными из других источников свидетельствует о надежности используемого численного метода для определения параметров и структуры вязкого ударного слоя, в том числе тепловых потоков к поверхности тела, в широком диапазоне параметров набегающего потока Количественные отличия обусловлены, скорее всего, различными моделями химических реакций сравниваемых методов, включающими в себя константы скоростей химических реакций и граничные условия. При этом отмечается хорошее согласование концентраций, достигших сравнительно больших величин (атомарный кислород). Несколько большие различия наблюдаются в концентрациях компонентов, величины которых малы. Для заряженных частиц (электронов с~) можно отметить удовлетворительное, с точки зрения использования расчетных значений в прогнозировании затуханий радиосигнала в прибортовой плазме, согласование концентраций с экспериментальными данными.
Проведенное численное исследование пространственного неравновесного обтекания затупленного конуса под различными углами атаки выявило характерные особенности и закономерности трехмерных течений реального газа. Так, при изучении структуры ударного слоя для режимов с большими углами атаки было установлено, что максимальная степень ионизации газа с подветренной стороны тела до трех раз меньше, чем с наветренной. С другой стороны, ионизованный слой с подветренной стороны имеет гораздо большую (примерно в 13 раз) толщину, чем с наветренной. Отсюда следует, что ионизация и диссоциация газа происходит преимущественно с подветренной стороны тела. Эти же процессы оказывают решающее влияние на параметры газа, находящегося в следе за телом. Поперечное течение газа в отрывной области на подветренной стороне тела в значительной степени влияет на протекание химических реакций в этой зоне. Поэтому для детального изучения структуры течения и химического состава газа точность описания газодинамических и химических процессов в зоне отрывного течения имеет весьма важное значение.
Исследование прохождения радиосигнала в пространственном неравновесном потоке показало, что величина отхода ударной волны с сохранением интенсивности уровня заряженных частиц на подветренной стороне оказывает существенное воздействие на характер затухания по сравнению с наветренной стороной. Так, для больших углов атаки на подветренной стороне величина затухания постоянно увеличивается по длине конуса..
Анализ полученных результатов свидетельствует о возможности применения предложенного метода для расчета параметров пространственного химически неравновесного вязкого ударного слоя, включая характеристики прибортовой плазмы, а также определение аэродинамических и тепловых нагрузок летательных аппаратов, интервалов прекращения радиосвязи. Затраты машинного времени остаются приемлемыми и значительно меньшими, чем при использовании полностью неявных схем во всей рассматриваемой области без выделения головной ударной волны..
Основные условные обозначения.
ВУС — вязкий ударный слой, ДСК — декартовая система координат, КСК — криволинейная неортогональная система координат,.
НССК — нормированная сферическая системы координат,.
НСК — нормированная система координат,.
ПВУС — полный вязкий ударный слой-.
ПУНС — параболизованные уравнения Навье-Стокса,.
ССК — сферическая система координат-.
ТВУС — тонкий вязкий ударный слой-.
ЦСК — цилиндрическая система координат-.
А — произвольное целое число, задающее отношение наибольшего шага сетки к наименьшемус — местная скорость звукар. с = — - массовая концентрация / -ои компонентыР.
С — удельная теплоемкость смеси при постоянном давлении-.
С' - молярная теплоемкость / -ой компоненты при постоянном давлении,.
I) — скорость распространения волны по частицам газа-.
Div — дивергенция тензора-.
Dt — эффективный коэффициент диффузии-.
1 Р с =-— - внутренняя энергия единицы массы газау-р а~ - электронная компонентас — энергия колебательных степеней свободы- (уЛ S.
Ki 5 ?2 'ft r р= ] р I р у-р у — р р К у Р V2 2 к = 1.38−10″ 21 Дж/К° к к, К, L.
Л Sct т т, тх.
N. Р.
— полная внутренняя энергия единицы объема,.
— напряженности поля падающей, прошедшей и отраженной волн-.
— матрица фундаментального метрического тензора,.
— ковариантные базисные векторар е + — - энтальпия газаР.
— энтальпия / -ой компоненты- -энтальпия образования / -ой компонент-.
— полная энтальпия газа в однородном набегающем потоке-.
— постоянная Больцмана-.
— комплексное волновое число, характеризующее плазменную среду-.
— константы скоростей обратных реакций-.
— константы равновесия-.
— характерный линейный размер-.
— числа Льюиса-.
— молекулярный вес смеси- -молекулярный вес / -ой компоненты-.
— молекулярный вес набегающего потока-.
— число Маха-.
— определяет номер сеточного узла, с которого вводится неравномерная сетка-.
— давлениеfiC.
Pr =— = 0.72 — число ПрандтляX fdp q = Я.
— тепловой потокqJl{ - нормальная скорость скачкаR, 0,.
Raуниверсальная газовая постоянная-.
Re^. = ^ - число Рейнольдса-.
Rn — радиус затупления носка-.
R, Т — комплексные коэффициенты отражения и прохождения-.
St = - - , — число Стантона-.
S, j — компоненты тензора скоростей деформации,.
Sw — длина дуги на поверхности тела, отнесенная к радиусу затупления,.
Sct = — D, — числа Шмидта / -ой компонентыР t — время-.
Д?АФЛ — временной шаг по критерию Куранта-Фридрихса-Леви для линейных гиперболических уравнений в частных производных-.
Т — температураи, v, w — компоненты вектора скорости V-.
J, = -pDt grad у, — массовый поток / -ой компоненты вследствие диффузии-.
V — вектор скорости-.
V", Vt — проекции вектора скорости на нормаль и касательную плоскость к поверхности ударной волны-.
Wt — источник образования / -ой компонентыx, y, z — координаты в ДСК, с, т xt = —— - мольные концентрацииmt аугол атаки-.
Р — коэффициент прохождения радиосигналаР — угол поворота СК-.
8* - толщина ударного слоя-.
— характеристические потенциалы столкновенийу — показатель адиабаты, у, — эффективный показатель адиабаты-.
П, с, у, = = — - размерные концентрации компонентовтк т, rjt = -J—?- - концентрации компонентовmt.
Ф — функции для источниковых членов реакций-.
Л — коэффициент теплопроводности-.
Л — длина волны излученияi — коэффициент вязкости неравновесной смесид — коэффициенты вязкости чистых однокомпонентных) газовv — эффективная частота столкновений-.
0 — угол полураствора конусар — плотностьр = ррл — размерная плотностьpt — плотность / -ой компоненты, сг — коэффициент запаса (~0.9), сг — диаметры столкновений- - координаты в НССК, Qj2 2]* - интегралы столкновений, VL — символы Кристоффеля,.
Индексы верхние индексы.
Н — обозначены невязкие члены уравненийВ — вязкие члены уравнений, нижние индексы нижний индекс означает дифференцирование по соответствующему аргументуиндекс / = 1, 2, 3,4, 5, 6, 7 ассоциируется с компонентами О, N, NO, ()2, N2, NO*, е~ соответственно-.
О — соответствует параметрам в точке торможенияw — соответствует параметрам на стенкеоо — невозмущенный набегающий поток.