Пространственные задачи сверхзвукового обтекания тел потоком вязкого газа

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

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

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

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

Исторически первой такой упрощенной моделью является модель пограничного слоя (НС) (предложенная Л. Прандтлем [1] еще в начале XX века), в которой все течение вязкого газа около поверхности обтекаемого тела при больших числах Рейнольдса разбивается на две части: на область пренебрежимо малого влияния вязкости (невязкое течение) и область пограничного слоя. Решение в невязкой области находится из уравнений газовой динамики (уравнения Эйлера), а в области вязкого течения — из уравнений Прандтля, представляющих собой предельную асимптотическую форму уравнений НС при числе Рейнольдса, стремящемся к бесконечности [2]. Такое приближение соответствует условиям задачи входа тел в атмосферу на низких участках траектории.

Классическая теория ламинарного пограничного слоя хорошо изложена в монографиях [3, 4], там же представлены некоторые результаты исследований двумерных задач теории ламинарного ПС. Обзор работ, относящихся к пространственным задачам ПС, дан в [5, 6]. Объективные трудности при интегрировании трехмерных уравнений ПС обусловили интерес к вырожденным пространственным течениям газа, когда все три компоненты скорости отличны от нуля, но зависят от двух, а иногда лишь от одной, переменных. К таким задачам относятся задачи о течении на линии торможения в окрестности критической точки двоякой кривизны [12, 13, 14, 15], обтекание цилиндра под углами атаки [17, 16, 18, 9], задача о течении в окрестности плоскости симметрии затупленных тел [19, 22, 25, 28, 43, 46, 45] и т. д.

Первые работы, в которых течение в пространственных ПС было рассмотрено в точной постановке без каких-либо упроп], аюш, их предположений, относятся к середине 60-х годов прошлого столетия: в работе [20 было рассмотрено обтекание затупленного по сфере конуса под углом атаки, в [21] проведено численное исследование ПС около эллипсоидов враш, ения при ненулевом угле атаки. Вопросы, связанные с расчетом ПС, возникаюш, его при обтекании вязким сверхзвуковым потоком газа под углом атаки сферического сектора, рассматривались в [23, 24. Численные и асимптотические решения трехмерных уравнений ламинарного ПС около затупленных тел с двумя плоскостями симметрии, обтекаемых под нулевым углом атаки при наличии вдува или отсоса получены в [48, 47]. Расчеты течения сжимаемого газа в ПС около длинных биэллиптических тел проведены в [29, 30]- обтекание под углом атаки и скольжения сегментально-конических тел и трехосных эллипсоидов с непроницаемой поверхностью — в [26]. Тела более сложной формы, моде-лируюш, ие сверхзвуковой аппарат рассмотрены в монографии [31 .

Все используемые методы при решении трехмерных по супдеству задач ПС можно разделить на две группы: приближенные аналитические или численные методы, используюпдие некоторые упрош, ающие предположения, и точные численные методы без какого-либо упрощения исходных уравнений пространственного ПС. К первой группе относятся метод осесимметричной аналогии [7], локально-автомодельное приближение [8], метод интегральных соотношений [9], многопараметрический метод Л. Г. Лойцянского [10], интегральный метод последовательных приближений [11]. Многочисленные исследования пространственных тече.

НИИ этими методами позволили прояснить качественную структуру течения, а в некоторых случаях и дать рекомендации для инженерных приложений. Обзор численных методов, относящихся ко второй группе, дан в [27, 32-.

Описываемый в первой главе настоящей работы метод решения трехмерных уравнений ПС относится к точным численным методам. Он позволяет построить полную картину течения в пространственном слое около проницаемой поверхности гладкого затупленного тела, обтекаемого гиперзвуковым потоком теплопроводного сжимаемого вязкого газа под различными углами атаки и скольжения, т. е. когда в поставленной задаче отсутствует какая-либо симметрия. Метод имеет повышенный порядок аппроксимации по координате, отсчитываемой поперек пограничного слоя. На основе этого метода проведено систематическое численное исследование обтекания под углами атаки и скольжения различных затупленных тел (трехосные эллипсоиды, эллиптические параболоиды и двуполостные гиперболоиды) с проницаемой (в том числе и с частично проницаемой) поверхностью. Проводится оценка влияния на структуру течения в пограничном слое, распределения компонент напряжения трения и теплового потока формы тела и температуры его поверхности, параметра вдува газа с поверхности, углов атаки и скольжения. Показано, что в зависимости от определяющих параметров задачи распределения теплового потока вдоль поверхности может носить качественно различный характер. Результаты, полученные с помощью этого метода и частично изложенные в главе I, опубликованы в работах 127, 128, 129, 130, 132, 146, 152, 151 .

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

Во второй главе данной работы рассматривается одна из самых известных моделей вязкого слоя — модель гиперзвукового или тонкого вязкого ударного слоя (ТВУС), предложенная И.К.СИе^'ом [100], в которой все течение разбивается на два слоя: область между поверхностью тела и задним фронтом ударной волны и непосредственно сама ударная волна, поперечные размеры которых при больших числах Маха считаются малыми. Уравнения ТВУС содержат все члены уравнений пограничного слоя и невязкого ударного слоя в гиперзвуковом приближении. Уравнения, описывающие область перехода через ударную волну, интегрируются и полученные соотношения (обобщенные условия Рэнкина-Гюгонио) используются, а качестве граничных условий на внешней границе вязкого ударного слоя. Многочисленные сравнения с экспериментом и расчетами в более полной постановке, анализ которых можно найти в обзорах 67, 6], свидетельствуют о удовлетворительной точности модели ТВУС при определении тепловых потоков и сопротивления трения там, где ударный слой является тонким.

Первые решения, полученные численными методами в рамках модели ТВУС, были получены для плоских и осесимметричных течений (подробный обзор двумерных решений уравнений ТВУС — в [67]). В дальнейшем появились работы, посвященные различным частным случаям пространственных течений: например, в [54, 33, 35] получено решение задачи в окрестности критической точки двоякой кривизны, в [34] - в окрестности плоскости симметрии эллиптического гиперболоида. Используя приближенные подходы в определении продольных составляюш-их градиента давления, которые в отличие от ПС заранее неизвестны, были проведены исследования при обтекании под нулевым углом атаки эллиптических параболоидов и гиперболоидов в [36], а в [107] - под углом атаки эллиптического параболоида. До публикации первых авторских статей в литературе не было представлено работ, посвяш-енных обш, ему случаю обтекания затупленных тел под углами атаки и скольжения.

Уравнения ТВУС являются параболическими уравнениями, также как и в ПС, поэтому наработанная методика расчета трехмерных уравнений ПС (первая глава) была применена и в данном случае. Излагаемый в работе строгий подход к расчету уравнений пространственного ТВУС позволил найти решение задачи в точной постановке. Для исследования используется цилиндрическая система координат с началом в критической точке на поверхности, что позволило выписать не содер-жап], ие особенности уравнения для продольной и окружной составляю-ш, их градиент давления и тем самым замкнуть задачу. Большое количество расчетов, проведенных в широком диапазоне изменения определя-юш-их параметров задачи, позволили сделать вывод о том, что данный численный метод решения уравнений пространственного ТВУС является эффективным и экономичным. Проведенные на основе этого метода исследования выявили сильную зависимость структуры течения в вязком слое от формы обтекаемого тела и величины углов атаки и скольжения, и как следствие, тепловых нагрузок на обтекаемую поверхность. Основные результаты, полученные при исследовании трехмерного тонкого вязкого ударного слоя в однородном газе около тел, обтекаемых под углами атаки и скольжения, опубликованы в работах 131, 133, 132, 136, 135, 139, 141, 140, 146 .

Однако модель ТВУС обладает супдественными недостатками. Так применение упрощенного уравнения импульсов в нормальном к поверхности направлении приводит к появлению на поверхности выпуклого тела точки с нулевым давлением (точки отрыва), за которую решение задачи не может быть продолжено. Кроме того, расчеты по модели ТВУС дают низкую точность вниз по потоку, где толщина ударного слоя становилась не малой. Частично эти ограничения были преодолены методом глобальных итераций по форме ударной волны (см. последний параграф главы 2). Ранее в работах [72, 73] были предприняты попытки сдвинуть точку отрыва ближе к миделю тела на основе метода двух приближений.

С этой же целью в третьей главе данной работы предлагается новая математическая модель течения — модель параболизованного вязкого ударного слоя (ПВУС). Основная идея предлагаемой модели состоит в том, чтобы, с одной стороны, использовать более точное (по сравнению с ТВУС) уравнение движения в проекции на нормаль к поверхности, служащее для определения давления, и тем самым существенно увеличить размер расчетной области по маршевой координате, а с другой стороны — сохранить достоинство модели ТВУС, состоящее в возможности прямого применения маршевых методов к решению исходных уравнений. На основе систематического исследования было определено, что область применимости модели ПВУС в асимптотическом смысле не отличается от области применимости модели ТВУС, но в рамках новой модели ПВУС существенно увеличивается расчетная область по маршевой координате (но сравнению с моделью ТВУС), что позволяет применять данную модель для расчетов практически на всей лобовой поверхности затупленных тел, вплоть до их миделя. Результаты, полученные при исследовании трехмерного ПВУС в однородном газе около тел, обтекаемых под углами атаки и скольжения, опубликованы в работах [142, 143, 144, 146, 150, 154 .

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

Ранее задача о течении многокомпонентного химически реагирующего газа в ТВУС в окрестности критической точки была решена в [56, 82. В ряде работ были рассмотрены частные случаи течения в вязком ударном слое, когда все три компоненты вектора скорости отличны от нуля, однако решение в силу той или иной симметрии задачи зависит лишь от двух переменных. К таким задачам, в частности, относятся задачи о течении в плоскости симметрии затупленного тела, об обтекании под углом атаки стреловидных крыльев бесконечного размаха и под нулевым углом атаки осесимметричных вращающихся тел (см. обзор [67]). Большое внимание уделено течениям с физико-химическими процеесса-ми в монографии [112]. Полностью данному вопросу посвящена работа 113.

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

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

Результаты, полученные при исследовании неравновесных течений в рамках моделей ТВУС и ПБУС около тел, обтекаемых под углами атаки и скольжения, опубликованы в работах [134, 132, 137, 138, 145, 148, 149, 158, 159 .

Стремление устранить перечисленные выше недостатки модели ТВУС привело к появлению модели полного вязкого ударного слоя (ВУС) [75. Уравнения ВУС включают в себя все члены уравнений Эйлера и пограничного слоя. В отличие от модели ТВУС в уравнениях ВУС члены с кривизной учитываются, уравнение импульсов по нормали берется полным (за исключением вязких членов), а форма ударной волны — не эквидистантна поверхности тела, а определяется в процессе решения. Отсюда следует расширение границ применимости модели ВУС: возможность расчета обтекания тел при небольших сверхзвуковых скоростях, а также в тех случаях, когда ударный слой не является тонким. Отсутствие в системе ВУС вторых производных по продольным координатам причисляет ее к системам эволюционного типа, что в принципе предполагает применение к ее решению маршевых методов. Однако, как показывает анализ [114, 61, 65], поскольку уравнения ВУС включают в себя все члены уравнений Эйлера, то невязкий оператор является гиперболическим в сверхзвуковых областях и эллиптическим в дозвуковых (а это пристеночная область и окрестность линии торможения). Для подавления возмущений против потока необходим механизм регуляризации решения. В настоящее время для этих целей широкое распространение получили метод установления и метод глобальных итераций (ГИ). в качестве такого в настоящей работе (глава пятая) используется метод установления, в котором в качестве исходных уравнений используются нестационарные уравнения ВУС, а необходимое решение получается как предел установления по времени.

В методе ГИ используется итерационный процесс, в котором в качестве глобальных итерируемых функций выступают неизвестная заранее форма ударной волны и продольные составляющие градиента давления, передающие возмущение вверх по потоку в дозвуковых областях. Метод ГИ впервые был использован при гиперзвуковом осесимметричном обтекании лобовой поверхности гиперболоида вращения в работе [75]. Для первого приближения использовалась модель тонкого (гиперзвукового) вязкого слоя. Проведенные расчеты показали хорошую сходимость метода, однако, для тел другой формы (с «толстым» ударным слоем) описанный в работе итерационный процесс расходился. В связи с этим в ряде работ [60, 61, 64, 65] (см. также [67, 6] и цитируемую там литературу) впервые был развит сходящийся метод ГИ для двумерных задач с произвольным «толстым» ударным слоем. В работе [66] описывается итерационная процедура метода ГИ для трехмерного случая, но конкретные расчеты приведены для тел вращения. В настоящей работе во второй, третьей и четвертой главах предлагаются различные модификации метода ГИ, позволившие провести расчеты обтекания вязким однородным газом и многокомпонентной химически реагирующей смесью газов трехосных эллипсоидов под различными углами атаки и скольжения, т. е. в самом общем пространственном случае. Так, в параграфе 3.4 проведена оценка приближенного подхода к вычислению продольных составляющих градиента давления в ИВУ С. А в последнем параграфе четвертой главы исследуется влияние уточненной формы ударной волны и точных продольных составляющих градиента давления на результаты численных решений.

В пятой главе наряду с решением задачи обтекания тел, поверхность которых задана аналитической формулой, рассматривается интересный с практической точки зрения случай, когда эта поверхность задана таблично с некоторой точностью. Для вычисления производных вдоль этой поверхности, входящих в коэффициенты уравнений ВУС, применялась оригинальная методика сглаживания, основанная на кубических В-сплайнах. Особенности метода решения пространственных уравнений ВУС, результаты численных исследований в рамках этой модели опубликованы в работах [147, 153, 156 .

В задачах, рассмотренных в четвертой главе, предполагалось, что поверхность движущегося по траектории спуска тела — теплоизолированная. В реальности же происходит интенсивное тепловое взаимодействие нагретого газового потока с теплозащитным покрытием летательного аппарата. Теоретическое решение этой задачи в общем случае должно быть получено при совместном рассмотрении процессов в газовой и твердой фазах [121], т. е. должно основываться на решении дифференциальных уравнений внешней газодинамики в многокомпонентных реагирующих газовых смесях в сопряжении с уравнением теплопроводности в защитной оболочке движущегося тела [122]. В шестой главе настоящей работы решена задача спуска по планирующей траектории в атмосфере Земли гладкого затупленного тела, обладающего аэродинамическим качеством и теплопроводной поверхностью. Траектория спуска представлена не только высотой и скоростью, как функциями времени полета, но и изменяющимися во времени углами атаки и скольжения. Уравнения в газовой фазе, которыми являются пространственные уравнения ПВУС для многокомпонентной смеси газов, решаются совместно с трехмерным уравнением нестационарной теплопроводности в твердой фазе. Ранее учет сопряженного теплообмена для ряда задач теории вязкого ударного слоя проводился в окрестности критической точки в [123], для осесимметричных течений — в [124], при пространственном обтекании с фиксированным углом атаки (в течении присутствует плоскость симметрии) — в [110]. Во всех упомянутых работах: потоком тепла вдоль покрытия пренебрегал ось.

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

Результаты, полученные при исследовании неравновесных течений около тел с теплопроводной защитной оболочкой опубликованы в работах [157, 160-.

Данная диссертация была выполнена при поддержке грантов РФФИ 93−013−17 957, 95−01−832, 98−01−298. Часть разработанных в диссертации методов, алгоритмов и программ внедрена в национальном институте информатики и автоматики INRIA (Франция).

Основные результаты по теме диссертации докладывались на семинаре в ЛФТИ АН СССР под руководством профессора Ю. П. Головачева, на семинаре в ТГУ под руководством профессора И. М. Васенина, на семинаре НИИ ПММ под руководством профессора И. Б. Богоряда.

Результаты исследований по теме диссертации докладывались на VI и VII Всесоюзных школах-семинарах по механике реагирующих сред (г.Междуреченск, 1986 г., г. Красноярск, 1988 г., г. Кемерово, 1990 г.), на XI Всесоюзной школе-семинаре по численным методам механики вязкой жидкости (г.Свердловск, 1988 г.), на II Всесоюзной школе-семинаре по термогазодинамике и макрокинетике (г.Томск, 1988 г.), на XII Всесоюзной конференции по вычислительной математике и математической физике (г.Одесса, 1987 г.), на IV и V Всесоюзных конференциях «Современные проблемы аэрогидродинамики» (г.Жданов, 1987 г., п. Рыбачье, 1989 г.), на III Всесоюзной школе-семинаре по макроскопической кинетике,.

16 химической и магнитной газодинамике (г.Томск, 1991 г.), на конференции «Фундаментальные и прикладные проблемы современной механики» (г.Томск, 1991, 1998, 2000 гг.), на II Japan-Soviet Unionjoint symposium on computational fluid dynamics (г.Цукуба, 1990 г.), на международных конференциях Computational Fluid Dynamics (г.Штутгардт, Германия, 1994 г., Г. Париж, Франция, 1996 г.).

ЗАКЛЮЧЕНИЕ

.

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

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

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

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

231 ла носит сложный, существенно несимметричный характер.

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

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

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