Исследование течений в вязком ударном слое при помощи схем высокого порядка аппроксимации

Тимченко, Сергей Викторович

Исследование течений в вязком ударном слое при помощи схем высокого порядка аппроксимации тема автореферата и диссертации по механике, 01.02.05 ВАК РФ

Тимченко, Сергей Викторович АВТОР

доктора физико-математических наук УЧЕНАЯ СТЕПЕНЬ

Ижевск МЕСТО ЗАЩИТЫ

1999 ГОД ЗАЩИТЫ

01.02.05 КОД ВАК РФ

Диссертация по механике на тему «Исследование течений в вязком ударном слое при помощи схем высокого порядка аппроксимации»

Автореферат диссертации на тему "Исследование течений в вязком ударном слое при помощи схем высокого порядка аппроксимации"

На правах рукописи

у;.

Тимченко Сергей Викторович

ИССЛЕд. 1 ШКЕ ТЕЧЕНИЙ В ВЯЗКОМ УДАРНОМ СЛОЕ ПРИ ПОМОЩИ СХЕМ ВЫСОКОГО ПОРЯДКА АППРОКСИМАЦИИ

01.02.0; - Механика жидкости, газа и плазмы

Автореферат

диссертации на соискание ученой степени доктора физико-математических наук

Ижевск - 1999

Работа выполнена в НИИ Прикладной математики и механики Томского государственного университета.

Научный консультант: доктор физико-математических наук

профессор С.В.Пейгин

Официальные оппоненты: доктор физико-математических наук

профессор А.В.Алиев, доктор физико-математических наук профессор И.М.Васенин, доктор технических наук профессор В.Н.Голубкин

Ведущая организация - Институт математического моделирования РАН, г. Москва

Защита диссертации состоится " 28" ¿1 ¿Щ г. в

'У? " часов на заседании диссертационного совета Д200.70.01 при Институте прикладной механики УрО РАН по адресу 426000. г. Ижевск, ул. М.Горького, 222

С диссертацией можно ознакомиться в библиотеке Института прикладной механики УрО РАН, г. Ижевск

■ Автореферат разослан _" дбКаЬРЯ__ 1999 г

Ученый секретарь диссертационного совета кан.ф.-м.н.

С.П.Копысов

2Го. 33V. во Г-^'1,

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

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

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

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

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

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

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

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

числе для сильно затупленных и выпукло-погнутых тел.

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

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

Автор защищает следующие основные результаты:

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

2. В сопряженной постановке решена задача о течении многокомпо-

нентного химически реагирующего газа в окрестности критической точки тела, движущейся вдоль траектории входа в атмосферу Земли перспективного аппарата NASA Х38.

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

4. Предложены модификации схемы Петухова для решения эллиптических уравнений.

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

6.Предложены новые подходы к созданию параллельных неявных трехмерных численных алгоритмов с высоким уровнем эффективности.

Практическая ценность и реализация работы. Данная диссертация была, выполнена при поддержке грантов РФФИ 93-013-17957, 95-01-00832, 98-01-00298. Разработанные в диссертации методы и программы внедрены в аэрокосмическом концерне Dassault Aviation (Франция), Российско-французском центре по прикладной математике и механике им. Ляпунова при МГУ и в национальном институте информатики и автоматики INRIA (Франция).

Апробация работы. Основные результаты диссертации докладывались на конференциях молодых ученых ИТПМ (г.Новосибирск, 1987,1988 гг.), на IV, У Всесоюзных конференциях по современным проблемам аэродинамики (г. Жданов, 1987г., п.Рыбачий, 1989 г.), на Школе молодых ученых по численным методам механики сплошной среды (г.Абакан, 1989 г.), на 7 Всесоюзной школе-семинаре "Современные проблемы механики жидкости и газа" (г.Иркутск, 1990 г.), на III Всесоюзной школе-семинаре по макрокинетике и химической и магнитной газодинамике Красноярск. 1990 г.), Школе по математическому моделированию (Абрау-Дюрсо, 1993 г.), конференциях "Применение математического моделирования для решения задач в науке и технике". (г.Ижевск. 1996, 1998гг.). Международной конференции по внутрикамерным процессам и горению ICOC-96 (г.С.-Петербург, 1996 г.), международных конференциях Parallel Computational Fluid Dynamics (г.Париж, Франция, 1993г., г.Киото, Япония, 1994 г., Г.Манчестер, Великобритания, 1997 г., г.Вильямсберг, США, 1999 г.), конференции "Всесибирские чтения по математике и механике" (г. Томск, 1997 г.), конференции "Математические модели и методы исследования" (г.Красноярск, 1997 г.), международной конференции по генетическим алгоритмам EVROGEN97, (г.Триест, Италия, 1997 г.. ), конференции "Фундаментальные и прикладные проблемы современной

механики"' (г. Томск, 1998 г.), Четвертом европейском конгрессе по вычислительной гидродинамике ЕССОМА8'98 (г.Афины, Греция, 1998 г.)

Объем работы: Диссертация состоит из введения, семи глав и заг ключения. Работа содержит 232 страницы машинописного текста, 97 рисунков, 24 таблицы и список литературы из 221 наименования. Общий объем 285 страниц.

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

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

В первой главе диссертации рассмотрена задача о течении однородного сжимаемого газа в окрестности критической точки при обтекании тела неравномерными потоками газа. В 1.1 приводится математическая постановка задачи, в которой исходные уравнения пространственного тонкого вязкого ударного слоя (ТВУС) и соответствующие граничные условия на поверхности тела и ударной волне записываются следующим образом (х1 и х2 - координаты на поверхности, ж3 - координата отсчитываемая по нормали к поверхности):

щ/ж^2^.;/2) = о, (1)

р( Би" + А><У) = - ^а^аЩ^р + Д3 (£ А»«в).

2 ОзР = -(1 + е)РА10и°и0

рВТ = ^-Рр + Г) + -£во„£>зиаА»гЛ

1+е \Ка ) А

Р = рТ, Ц-Т", Ва)3 = а0/ А = д/дх', Б* = и"«^? А,, П = В* + гг3£»3.

А" = £йе, йе = Т0 = - 1 е 7 1

7 + 1

х3^!1,*2): р{чг-0*х1)=р00ч10, 2Р = (1 + е)р00(и^с)2, (2)

- О = чК-^ги",

■ РооиЦТ - 0.5(и^)2 - Ва/>(«° - ОК - «£,)) = МгАТ^зГ,

-3-0: и° = «°(®1)а:2), Г = Гм(г\х2), (3)

.т

„3\ _ /-/„1 „2\

(/ш^СО*1,^).

Здесь все линейные размеры отнесены к характерному линейному размеру Я, а нормальная координата х3 к ей; У^и", еУ^и3 - физические компоненты вектора скорости в направлениях ха, р°00е~1Т0р/Т£0, е~1р°уср, Т0Т, о; - соответственно давление, плотность, температура, коэффициент вязкости и число Прандтля.

Рассматриваются следующие два основных типа течения.

1.Обтекание тела потоком типа дальнего следа. Параметры набегающего потока в этом случае могут быть определены по формулам

1-ехр(-6г2) 1

И= 1 — а ' * = ^=1 + с(1-У?) (4)

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

-1 ¡V оо

Для удобства численного решения и аналитического исследования задачи (1), (2), (3) осуществляется переход к новым зависимым и независимым переменным типа А.А.Дородницына:

1 ^ в

£° = ха, £=д $ Р\?а<1х\ Д = | ¿ =

рра

дфа

2Т

ТХуг'

1 др

р1 Р*

ру}л/а -

,/а^ Эх'*

Показано, что при примененном обезразмеривании влияние неравномерности на течение в окрестности критической точки проявляется только через комплекс Л = 2а6(1 + с)/(1 — а), входящий в граничное условие для градиента давления и являющийся параметром подобия для таких течений. Во втором случае влияние неравномерности также проявляется через единственный параметр <1.

В 1.2 рассмотрено асимптотическое решение поставленной выше задачи при больших числах Рейнольдса и сильном вдуве.

В слое вдува получено аналитическое решение задачи:

и1 = и2 = р2

1 +С&

1 - сл

Сг|*2 - 1|1УМ/2

1 4 — /

г, А, •/

усИ

И-II

с, =

* и + 1|'

1 - сл

< * < 1,

С2

А, - За 1

ЯпГ'а

где 5а,га,та определяются типом и интенсивностью неравномерности набегающего потока.

Анализ показывает, что отвечающее физическому условию конечности толщины слоя вдува решение существует только при -и* > —,За.

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

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

I [да -яв Г v2 2 gi а в QeJ >

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

«¿, = 14 = 1, G = —1, /JQ =

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

.of, к(к-Ъ) к\пк )

При Л < Л* реализуется безотрывное невязкое обтекание затупленного тела, движущегося в следе за другим телом. Из (8) следует важный вывод: при одной и той же неравномерности набегающего потока ответ па вопрос: будет ли иметь место безотрывное обтекание, либо у лобовой поверхности произойдет отрыв потока, зависит от геометрии тела, характеризуемой в окрестности критической точки тела параметром к. В частности, при к — 1 (осесимметричное тело) Л* максимально и равно 1,33, а при fc = 0 (плоское тело) Л*. = 0, то есть отрыв в этом случае происходит даже при малой интенсивности следа.

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

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

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

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

д(11е,к, с!) = У[ГТа^Ч°(Т1е*,к*), (9)

г, . Ке <1 ,, 1 + ка

Ие* = --к* = ---

1 + <Г 1+с1

где q° - тепловой поток в критическую точку тела, обтекаемого равномерным потоком газа.

В 1.4 исследовано совместное влияние неравномерности набегающего потока и нестационарности задачи на течение а окрестности линии торможения при прохождении тела через температурную неоднородность. При этом распределение температуры набегающем потоке в декартовой системе координат (у1, у2, у3), начало которой лежит в плоскости симметрии слоя, а ось Оуъ направлена по нормали к этой поверхности, описывается формулой:

Рисунок 1. Профили и1 при к = 1, #та = 0.15.

Рисунок 2. Профили и1, и2 и в около проницаемой поверхности при & 0.1,X = 0.25, Ые = 103, (I = 106, 3, 0.5 - линии 1, 2, 3, в = -0.2.

fj-lo

Too= 1-аеМ-(У3/Ь)2) (10)

Показано, что тепловой поток в критическую точку, отнесенный к своему значению при равномерном обтекании, является величиной консервативной и слабо зависит от числа Рейнольдса, температуры поверхности и геометрии тела (рис.3). Проведена оценка возможности решения лашюй^адачи-в-квазистационарнойтюстановкёТЛкжазано, что неучет нестационарных членов уравнений приводит как к количественным, так и к качественным ошибкам при расчете параметров течения в ударном слое. Из полученных численных решений следует, что зависимость решения от толщины температурной неоднородности можно значительно ослабить, если в качестве переменной использовать не т, а величину т* = (т-т0)/Ъ*/Ь + т0.

В 1.5 предложен эффективный численный метод для решения сопряженной задачи о течении многокомпонентного химически реагирующего газа в окрестности критической точки тела движущейся вдоль траектории входа в атмосферу Земли и решена задача для траектории перспективного аппарата NASA Х38.

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

Re—too

анализазадачи при больших числах Рейнольдса. Наличие закрутки прак-

тически не влияет на значение Л* при Ле > 103 и может значительно увеличивать область безотрывного обтекания при малых значениях числа Рейнольдса. Анализ полученных численных решений показывает, что наличие неравномерности в набегающем потоке газа может приводить, по сравнению с равномерным обтеканием, к качественному изменению характера течения в ударном слое на боковой поверхности. В частности, возможны режимы, при которых характерное местное число Рейнольдса, вычисленное по параметрам за ударной волной, значительно меняется по обводу тела с ростом продольной координаты вследствие чего режимы течения с ударном слое на боковой поверхности и в окрестности критической точки могут существенно различаться: если при равномерном обтекании режим течения не меняется по всей его длине, то при наличии неравномерности в набегающем потоке в окрестности критической точки реализуется режим вязкого ударкого слоя, в то время как на боковой поверхности при достаточно больших £ рке ясно вырисовывается зона невязкого течения, примыкающая к ударной волне, и вязкий пограничный слой вблизи поверхности тела.

В 2.3 рассмотрено неравномерное обтекание удлиненных осесимме-тричных тел. Показано, что при течении на боковой поверхности наряду с параметром Л влияние неравномерности на течение зависит также и от характерного радиуса следа г* ~ б-1/2. Вид этого влияния в сильной степени определяется соотношениями между Л, Яе и г*. Показано, что имеются целые классы тел, для которых условие А < Л* не является достаточным условием безотрывного обтекания, поскольку для них градиент давления на поверхности может стать положительныл не в критической точке, а при £ > 0. Это связано с тем, что дРы/д£, зависит как от градиента полного давления в набегающем потоке, так и от к - кривизны ударной волны. Поскольку при £ 1 для некоторых тел, например для гиперболоидов, к —> 0, то определяющим может стать влияние градиента полного давления, который, как рке отмечалось выше, для следа положителен. Таким образом, возникновение отрыва потока на боковой поверхности тела определяется соотношением порядков, с которым стремятся к нулю продольная кривизна поверхности тела и градиент полного давления в набегающем потоке. В частности, для фиксированного Л при малых радиусах следа (большие Ь) быстрее стремится к нулю неравномерность потока, и в этих случаях - монотонно убывающая

функция. Однако с увеличением г* ситуация меняется - распределение Рш({) теряет монотонный характер и на боковой поверхности возникает характерная "ложка" давления (рис.4). При этом величина Рш с хорошей

точностью может быть подсчитана по следующей формуле:

Pw = (pcoV? sin а — — JrpooVi cos'2 adrj (11)

Влияние неравномерности на течение при обтекании тела потоком из сверхзвукового источника носит качественно иной характер по сравнению с обтеканием потоком типа следа. В частности, в этом случае "ложка" давления в распределении Р„, не возникает. Р„,(£) - всегда монотонно убывающая функция, а отрыв потока может произойти лишь при малых d. вследствие появления на поверхности тела точки с нулевым давлением. В окрестности критической линии уменьшение d, ведет к уменьшению толщины ударного слоя и увеличению Сг По мере удаления от критической точки характер зависимости С, от d меняется - уменьшение d при достаточно больших f ведет к понижению абсолютных значений теплового потока. При этом важно отметить, что даже при достаточно малой неравномерности (большие d ) это изменение при расчете обтекания удлиненных тел необходимо учитывать. Например, если при d — 100 для £ < 5 отличие Cq от равномерного обтекания не превышает 0.5 %, то при £ = 50 данное отличие составляет уже 35 %.

Последний параграф данной главы посвящен численным и аналитическим решения уравнений вязкого ударного слоя в закрученном потоке газа. При малых числах Рейнольдса методом последовательных приближений получено приближенное аналитическое решение задачи. Расчеты показали, что наличие закрутки газа в набегающем потоке оказывает влияние не только на количественный, но и на качественный характер течения в вязком ударном слое. По мере увеличения Í) профиль продольной скорости и теряет свой монотонный характер и приобретает характерный максимум вблизи ударной волны. Наличие такого максимума объясняется сильным влиянием на течение газа в этой области ударного слоя центробежных сил, величина которых в данном месте становится значительной. При фиксированном параметре П значение максимума растет с увеличением числа Рейнольдса и стремится к некоторому. конечному пределу. Такое поведение решения хорошо согласуется с асимптотическим анализом течения при Re —» оо в рамках традиционной модели "невязкий ударный слой - пограничник слой".

Рисунок 3. Зависимости коэффициента теплообмена Сч(т) от времени. к = 0.1, вы = 0.25, т0 = 1.5, Ь = 0.45.

Рисунок 4. Распределение давления вдоль поверхности гиперболоида вращения. 6^-45°, Р.е = 10 (сплошные линии), 11е = 10° (штриховые), х х х - формула (11), Л = 0 (линия 5), Л = 2/3, Ь =1, 1.2. 1.5, 5 (линии 1-4).

кретизируется система координат на поверхности тела. Пусть {t1} - декартова система координат, в которой уравнение поверхности обтекаемого тела имеет вид:

JP(Í1,Í2,Í3)= О,

a f(sin a sin /3, sin а cos /3, cos а) - единичный вектор, совпадающий по своему направлению с вектором скорости набегающего потока (а -угол атаки, ,0 - угол скольжения). Обозначим радиус - вектор в этой системе через £ Радиус-вектор tp определяет критическую точку (точку, в котороÍÍ вектор Г перпендикулярен к поверхности). От системы координат {?'} перейдем к декартовой системе координат {г1} по формуле:

£- Í0 = B(A(i)) (12)

где В - диагональная матрица, определяющая преобразование сжатия (растяжения), А - матрица, задающая поворот осей координат таким образом, чтобы ось г2 была направлена по нормали к поверхности F внутрь тела, а оси г1 и z2 лежали в касательной плоскости к поверхности F в критической точке. Далее выберем точку С на оси г и введем сферическую систему координат г,хх,х<2\

z1 = rsin(mo:1)cosx2, г2 = г sin(ms:1) sin х2, (13)

г3 - = cos(mxí),

где т = т(х2) является известной функцией.

Для замыкания исходной системы и определения продольных градиентов давления Ра к уровнению импульсов в проекции на нормаль

применяются операторы (^1)_1а/ЗС1 и (S1)-^/^2-

В силу выбора системы координат на поверхности численное решение двухточечной начально-краевой задачи имеет ряд особенностей. Прежде всего, поскольку меняется в диапазоне от 0 до 2-к, то все геометрические характеристики, коэффициенты в уравнениях и граничных условиях, а также искомые величины являются периодическими функциями по координате £2. Кроме того, из-за вырожденности системы £2) при = 0 для получения решения уравнений пространственного тонкого вязкого ударного слоя в окрестности критической точки использовалась невырожденная криволинейная система координат {.(/'}: ее начало находилось в критической точке, ось Оу3 совпадала с нормалью к телу, а оси Оу1 и О у'2 образовывались пересечением касательной плоскости к поверхности тела в точке О с главными нормальными сечениями тела в данной точке. Получающаяся при этом система обыкновенных дифференциальных уравнений аналогична системе, решение которой было получено в первой главе диссертации. Далее решение и".

9а в критической точке в системе координат {у!} пересчитывалось в систему координат ().

Для численного решения системы осуществлялся переход от дифференциальных уравнений к разностным, производные по аппроксимировались - "разностями назад", производные по £2 при помощи раз-■ ностей против потока:

9/ _ /;+1 ~ /ы 0 г!3-2 - ЗД_1 + 3/, - /;+1, 2

а?" где2 + зд^ 1 р > ]

2Де2 +»••> ЗДС2 (при и < и)

в направлении С используется неявная разностная схема четвертого порядка аппроксимации. Таким образом, общий порядок аппроксимации используемой разностной схемы 0(А£г) + 0(Д£2)3 + 0(А()4.

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

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

дП-

^^СЭз+^Ил^+т^)^ 0" = 1,2) (14)

4 4 4 dui 4 du¡

Щ = E rj¡Ui, Qj = E q}iUi, Fj ~ E = £ "j.'T,

¿=0 ¿=D »=2 ¿=2 OÍJ

Щ = 1, Ui = i¿2i u2 — u3, u3 = =

притледующих-граничных-условиях-

¿Aj;u,=0, (j = 1,2,3-C = 0, j = 4-C = Ce) (15) ¡=0

Далее путем некоторой разностной аппроксимации дифференциальных операторов Fj и JV,- эта система уравнений в частных производных сводится к системе обыкновенных дифференциальных уравнений относительно координаты С Щ — Qj с последующим поэтапным решением набора соответствующих одномерных краевых задач. При этом значения щ в соседних точках, необходимые для аппроксимации частных производных в направлениях £ и т/ , либо лежат в области, в которой решение уже известно, либо берутся с предыдущей итерации.

В 4.2 излагается метод решения этой системы уравнений. Для построения двухточечной конечно-разностной схемы, аппроксимирующей систему ОДУ с шестым порядком точности, записываются четыре соот-

(1) (0)

ношения, связывающие искомые величины и] и и] в соседних узлах.

¿bM^bSM1', 7 = 1,2,3,4 (16)

Для этого используются следующие формулы /(1) _ /(0) = ДС(7/|(0) +32/<(1/4) + 12/'(1/2) + 32/(3/4) +7/'(1)) + 0(ДС7)

в которых / и /' заменяются при помощи исходной системы ОДУ, а значения в нецелых точках заменяются при помощи формул

„И _ 3 г+ (,(?) 15ó>f , Ю5 ffüf 105 ¿>f ,

W _ (аЛС)2 м , 4- 35 S° U"3] + Of

'з - 96 + 8 d„«a - 32 а 3 +16^лГ + 0(Л(:)-

= + ^^ - ^rf + + о(лс8).

Здесь введены обозначения

СЮ = С + 0Д<, f[0) = f(C(0)) /С) = + = /(/7+а) -

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

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

Пятая глава диссертации посвящена применению методов высокого порядка аппроксимации к решению эллиптических уравнений. В 5.1 предложен вариант метода установления для решения уравнений вязкого ударного слоя. В криволинейной ортогональной системе координат (.i\ у) нестационарные уравнения ВУС имеют безразмерный вид:

f №г) + pur) + ~(pvHir) = 0 (17)

иьк\ __ 1 ЭР 1 д (Н'1гц(ди ик

и2к\ дР

1 д /Н1ГудТ\ _ «к\2

Я Н\г ду \RePr ду ) Re \£>у Н\)

Р = 1^1рТ

П 9 . " 9 , 5 Т™ ЕГ 1 ,

= Ы + + , "1 = 1 + ук, >~ = гш 4- у соб а

3/ Н\ ах ду

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

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

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

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

Эти методы сочетают в себе естественный отбор среди строчных структур с упорядоченным, (хотя в чем-то и случайным) обменом информацией. Будучи вероятностными, генетические алгоритмы тем не менее не I являются просто еще одним вариантом случайного поиска, поскольку ; при отборе новых точек с ожидаемыми более хорошими возможностями, ] они эффективно используют предыдущую информацию. Данные алгоритмы носят итерационный характер и имеют дело с обработкой популяций индивидуумов Р(£) = {«1,..., х^} для итерации Ь (поколение 4). Каждый индивидуум представляет собой потенциальное решение задачи (испытание) и является в данном случае строкой (хромосомой) С, ; состоящей из Н.ем~2 элементов (генов), представляющих собой ординаты ; контрольных точек. Каждое решение х\ оценивается и определяется ме- ! ра его "пригодности". Затем формируется новая популяция (итерация или поколение Ь 4-1).

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

На следующем шаге отобранные таким образом индивидуумы подергаются преобразованиям с помощью "генетических операторов": му-ации и скрещивания. Оператор мутации т, создает нового индивидуума утем относительно малого изменения в одном индивидууме (гп,- : й

S), а оператор скрещивания cj - осуществляет более сильные трансформации и создает нового индивидуума путем комбинирования частей из нескольких (двух или больше) индивидуумов (с3- : S х ... х S S).

Сначала все особи из промежуточной популяции Р' попарно либо подвергаются скрещиванию (с вероятностью Рс), либо остаются без изменения. Наиболее простым является одноточечное скрещивание - каждая выбранная таким образом пара строк скрещивается следующим образом: случайным образом выбирается положение точки сечения (целое число к~ъ~п ромежуткс от-! _и l — l ^ где дл и на строки) ^Затсм,мух£м_

обмена всеми элементами между позициями k + 1 и I включительно, рождаются две новых строки. Пусть А = {у\,у1,уъ,у\) И А' = (у\, у'2,у'3, у'4) являются родителями, выбранными в процессе селекции. Тогда, (считая, что случайно выбранная точка сечения находится после первого гена), они производят двух детей В = (уиу'^уз, Уа) и В' = (у'ьу2,2/3,2/4)-После этого дети замещают родителей в промежуточной популяции Р'.

На следующем шаге все особи из промежуточной популяции Р' с некоторой вероятностью Рт подвергались мутации. Использовалась неоднородная мутация Михалевича. Если ген у,- подвергается мутации, то его новое измененное значение у\ выбирается внутри интервала [Min,-, Maxi\ следующим образом:

у\ = У1 + 8{М{-Уг){ (18)

где s - случайное число из интервала [0,1], д - номер поколения, G -максимальное число итераций, Ъ - параметр уточнения, М,- случайным образом выбирается из Мгпг и Max,, Mini и Махг - нижняя и верхняя границы возможного изменения значения переменной у,- Такая адаптивная мутация позволяла соблюдать в процессе реализации ГА (эволюции) необходимый баланс между двумя разномасштабными изменениями (мутациями) генов, так как на первоначальных шагах алгоритма в основном преобладали крупномасштабные изменения (обеспечивающие широкую область поиска), в то время как на заключительном этапе (за счет уменьшения масштаба мутаций) происходило уточнение решения.

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

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

В 6.2 первый из предложенных методов применен для решения задачи оптимизации затупления лобовой поверхности осесимметричного тела по конвективному тепловому потоку, полученному в рамках модели ламинарного пограничного слоя. Оптимальное решение, полученное на основе ГА сравнивается с решением, полученным с использованием традиционного подхода с применением модифицированного метода Давидона для поиска экстремума нелинейной функции многих переменных. Показано, что ГА позволяет достигнуть более лучших результатов: интегральный тепловой поток для решения с использованием ГА меньше, чем соответствующее значение для градиентого метода на 3-5%. На основании систематических численных расчетов сделан вывод, что предложенный метод решения вариационной задачи является устойчивым, робастым, имеет высокий уровень вычислительной эффективности и позволяет получать решение в широком диапазоне изменения параметров задачи.

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

<г(2,Т:2) если Т.^ <

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

Рисунок 5. Изолинии Сг 11е = 103, а = (3 =45°, Л = 10, Ь = 40.

1.50 X 100 000

Рисунок 6. Оптимальные формы затупления, полученные в рамках модели ПС (линии а) и ВУС (линия б), а также соответствующие им распределения теплового потока вдоль поверхности (линии виг) при Т*=1470К; 1500К; 1550К (линии 1,2,3)

оптимизационной задачи в рамках вязкого ударного слоя. Например, оптимальная форма, полученная в рамках пограничного слоя при Т* = 1478К с графической точностью совпадает с формой, полученной при помощи вязкого ударного слоя при Т* = 1500К.

В 6.4 рассматривается задача оптимизации параметров траектории входа в атмосферу Земли по величине суммарного конвективного теплового потока в точке торможения затупленного тела. В качестве исходной математической модели для расчета теплового потока используются уравнения тонкого вязкого ударного слоя с учетом неравновесного характера протекания химических реакций и многокомпонентной диффузии и сопряженного теплообмена на поверхности тела. Оптимальное решение находится при помощи ВГА-2.

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

Во всех случаях вычислительная область П разбивалась на N непересекающихся подобластей Qj (j = 1,..., N). Каждая подобласть дискре-тизируется на Л'у сеточных точек в направлении х2.

п = U л,-,

3=1

Оi - {(г\ж2,а:3)|0 < х1 < х\-,х) < х2 < г2+1;0 < г3 < х3},

i=i

где L - общее число сеточных точек по координате хг.

В большинстве случаев для MIMD-компыотеров равномерная загрузка (uniform, loadbalancing) процессоров и максимум компьютерной производительности достигается в случае, когда Kj = const и все подобласти

имеют одинаковый размер. Это является верным и для вычислительных систем с сетевой структурой в случае рабочих станций (компьютерных кластеров) с одинаковой производительностью. Для кластеров различной производительности равномерность загрузки в этом случае не достигается и для ее улучшения использовался специальный метод декомпозиций области, при котором число точек сетки К, размещаемых на каждом кластере было пропорционально его производительности "УУБ,

-ЛУЭ^-

Е WSj

»=1

где WS;- - машинозависимая константа, которая измеряется при помощи широкоизвестного эталонного теста Whetstone. Показано, что даже для случая, когда компьютеры в сети имеют существенно различную производительность эффективность распараллеливания достигает 90 %.

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

Основные результаты работы, выводы и рекомендации

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

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

3. В сопряженной постановке решена задача о течении многокомпонентного химически реагирующего газа в окрестности критической точки тела, движущейся вдоль траектории входа в атмосферу Земли перспективного аппарата NASA Х38.

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

скими методами исследовано течение в ударном слое около тела, обтекаемого закрученным потоком газа. Получено приближенное аналитическое решение задачи методом последовательных приближений, даю7 щее хорошую точность уже в первом приближении для Re < 50. Показано, что влияние закрутки набегающего потока приводит к смещению максимума теплового потока из критической точки на боковую поверхность тела и более раннему отрыву ударного слоя по сравнению со случаем равномерного.обтекания.____

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

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

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

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

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

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

сопряженного теплообмена.

10. На основе принципа геометрического параллелизма предложены подходы к созданию численных алгоритмов с высоким уровнем эффективности для решения трехмерных уравнений тонкого вязкого ударного слоя и турбулентного пограничного слоя на основе неявной конечно-разностной схемы высокого порядка точности на сети компьютеров с различной производительностью либо на многопроцессорной ЭВМ с распределенной памятью. Предложен вариант метода встречной прогонки тгего-адаптация-для-параллельной-вычислительной-системы-на-примерс— решения уравнений полного вязкого ударного слоя. Приведен асинхронный параллельный генетический алгоритм для решения задач оптимизации. Показано, что во всех случаях эффективность распараллеливания достигает 90% и более.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ ОПУБЛИКОВАНЫ В РАБОТАХ

1.Пейгин C.B., Тимченко C.B. Гиперзвуковой пространственный вязкий ударный слой в неравномерном потоке газа в окрестности критической точки// Изв. АН СССР. МЖГ. 1987, №6., С.136-145.

2.Пейгин C.B., Тимченко C.B. О численных и аналитических решениях уравнений вязкого ударного слоя в закрученном потоке газа// В сб.:Аэродинамика, Томск. 1987. С.71-79

3.Пейгин C.B., Тимченко C.B. Влияние неравномерности потока на трение и теплообмен в пространственном вязком ударном слое в окрестности критической точки// В сб.: Аэродинамика нестационарных процессов, Томск. 1988. С.83-94

4.Пейгин C.B., Тимченко C.B. Тонкий вязкий ударный слой около вращающегося затупленного тела, расположенного в дальнем следе// Ученые записки ДАГИ. 1989, Т.20, №2., С.83-94.

б.ПейгинС.В., Тимченко C.B. Нестационарный пространственный вязкий ударный слой в неравномерном потоке газа// Численные методы механики сплошной среды. Красноярск. 1989. 4.1. С.93-94.

6.Пейгнн C.B., Тимченко C.B. Обтекание удлиненных осесимметрич-ных тел сверхзвуковым неравномерным потоком газа// ПМТФ, 1989, №5., С.60-65.

7.Бородин А.И., Пейгин C.B., Тимченко C.B. Пространственный тонкий вязкий ударный слой в неравномерном потоке газа при отсутствии в течении плоскостей симметрии// Математическое моделирование. 1989, Т.1, №11.. С.51-57.

8.Пейгин C.B., Тимченко C.B. Нестационарный пространственный вязкий ударный слой в неравномерном потоке газа в окрестности критической точки// Изв. АН СССР. МЖГ. 1990, №3., С.175-180.

Э.Бородин А.И., Пейгин С.В., Тимченко С.В. Пространственное обтекание затупленных тел неравномерным гиперзвуковым потоком вязкого газа// Макроскопическая кинетика, химическая и магнитная газодинамика. Томск:Изд. Томск, унив., 1091, 4.1. С.126.

10.Пейгин С.В., Тимченко С.В. Обтекание удлиненных осесимметрич-ных тел, раположенных в дальнем следе// В сб.:Газовая динамика, Томск. 1991. С.68-74

Н.Бородин А.И., Пейгин С.В., Тимченко С.В. Пространственное обтекание затупленных тел неравномерным гиперзвуковым потоком вязкого газа// Теплофизика высоких температур. 1992, Т.ЗО, №1.

12.Peigin S.Y., Timchenko S.V. The Investigation of 3D Viscous Gas Flow over Complex Geometries Using Multiprocessor Transputer System. In: Parallel Computational Fluid Dynamics. New Trends and Advances. J.Periaux et al (ed). Proc. Parallel CFD'93 Conference. Elsevier. Amsterdam, 1995. P.13-20.

13.Пейгин C.B., Тимченко C.B., Четверушкин Б.Н., Чурбанова Н.В. Моделирование течений в пространственном вязком ударном слое на многопроцессорных вычислительных системах с распределенной памятью// Математическое моделирование. 1993, Т.5, №7.

14.Peigin S.Y., Timchenko S.V., Kazakov V.Yu. Application Multiprocessors Transputer System for Numerical Solution of 3D Nonequilibrium Hypersonic Viscous Gas Flow over Complex Bodies. In: Parallel Computational Fluid Dynamics. New Algorithms and Applications. Proceedings of Parallel CFD94 Conference. Elsevier Science B.V., 1995, Amsterdam. P.367-374.

15.Тимченко С.В. Расчет пространственного турбулентного пограничного слоя на сети ЭВМ различной производительности// Применение математического моделирования для решения задач в науке и технике. Ижевск. 1996. С.274-283

16.Тимченко С.В. Моделирование пространственных задач аэродинамики на многопроцессорных и многомашинных вычислительных комплексах/ / Всесибирские чтения по математике и механике. Томск:Изд. Томск, унив., 1997. Т.2., С.168.

17.Тимченко С.В. Применение схем высокого порядка для моделирования течений в ударном слое на многопроцессорных вычислительных комплексах// Математические модели и методы их исследования. Красноярск. 1997. С.175-176.

18.Timchenko S,V. The Effective Parallel Algorithm for Solution of Parabolic Partial Differential Equation System. In: Parallel Computational Fluid Dynamics. Resent Developments and Advances Using Parallel Computers. Proceedings of Parallel CFD'97 Conference. Elsevier Science B.V., 1998, Amsterdam. P.375-382

19.Peigin S.V., Timchenko S.V. Application of the Networked Computers

for Numerical Investigation of 3D Turbulent Boundary Layer over Complex Bodies. In: Parallel Computational Fluid Dynamics. Resent Developments and Advances Using Parallel Computers Proceedings of Parallel CFD'97 Conference. Elsevier Science B.V., 1998, Amsterdam. P.185-194

20.Пейгин C.B., Тимченко C.B. Метод высокого порядка аппроксимации для решения краевых задач параболического типа. Часть 1.//Математическое моделирование. 1998. Т.10, №4, С. 70-82.

-21:Пейгш1-С.-В.-—Peiiaux-J^,-Тимченко С.В^Лрименени

алгоритмов для оптимизации формы тела по тепловому потоку. //Математическое моделирование. 1998. Т.10, №9, С. 111-122.

22.Бородин А.И., Пейгин С.В., Тимченко С.В. Расчеты трехмерного пограничного слоя методом высокого порядка аппроксимации. //Математическое моделирование. 1998. Т.10, №10, С.79-87

23.Пейгин С.В., Тимченко С.В. Расчет пространственного турбулентного пограничного слоя на сети компьютеров различной производительности// ЖВМиМФ, 1998, т.38, №3, С. 510 - 519

24.Казаков В.Ю., Пейгин С.В., Тимченко С.В. Решение уравнений неравновесного вязкого ударного слоя для затупленных тел с каталитической поверхностью// ЖВМиМФ, 1998, т.38, №5, С. 860 - 869

25.Пейгин С.В., Тимченко С.В. Численное исследование обтекания осесимметричных затупленных тел сверхзвуковым потоком вязкого газа// Фундаментальные и прикладные проблемы современной механики. Томск:Изд. Томск, унив., 1998. С.240-241.

26.Mantel В., Peigin S., Periaux J., Timchenko S. A Heat Flux Optimization Using Genetic Algorithms// Computitional Fluid Dynamics'98. Proc. of the Forht Eoropean CFD Conference. J.Wiley&Sons. 1998. V.l, p.365-370

27.Иванов В.А., Пейгин С.В., Тимченко С.В. Численное исследование обтекания осесимметричных затупленных тел сверхзвуковым потоком вязкого газа// Изв. АН РАН МЖГ. 1999 г. №1. С.107-117

28.Borodin A., Mantel В., Peigin S., Periaux J., Timchenko S. Application of the Genetic Algorithm for Heat Flux Optimization Problem//MathematicaI survey for industry. 1999.

29.Казаков В.Ю. , Пейгин C.B., Тимченко С.В. Оптимизация траектории входа в атмосферу земли по интегральному тепловому пото-ку//ПМТФ. 1999. - в печати.

30.Peigin S., Periaux J., Timchenko S. Asynchrone Parallel Genetic Algorithm for Heat Flux Optimization Problem// In: Parallel Computational Fluid Dynamics. Development and Applications of Parallel Technology. Elsevier Science B.V., 1999. Amsterdam. P.377-384.

Содержание диссертации автор исследовательской работы: доктора физико-математических наук, Тимченко, Сергей Викторович

ВВЕДЕНИЕ.

1 ИССЛЕДОВАНИЕ ТЕЧЕНИЙ В ВЯЗКОМ УДАРНОМ СЛОЕ В ОКРЕСТНОСТИ ПРОСТРАНСТВЕННОЙ КРИТИЧЕСКОЙ ТОЧКИ.

1.1. Математическая постановка задачи

1.2. Асимптотическое решение уравнений пространственного вязкого ударного слоя в неравномерном потоке при больших числах Рейнольдса и сильном вдуве.

1.3. Численное решение задачи.

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

1.5. Химически неравновесный вязкий ударный слой на каталитической поверхности с учетом сопряженного теплообмена

2 ТЕЧЕНИЕ ОКОЛО ВРАЩАЮЩИХСЯ ОСЕСИММЕТРИЧ-НЫХ ТЕЛ ПРИ ИХ НЕРАВНОМЕРНОМ ОБТЕКАНИИ

2.1. Система уравнений и граничные условия

2.2. Обтекание вращающейся сферы неравномерным потоком газа.

2.3. Неравномерное обтекание удлиненных осесимметричных тел

2.4. Численные и аналитические решения уравнений вязкого ударного слоя в закрученном потоке газа

3 ОБТЕКАНИЕ ЗАТУПЛЕННЫХ ТЕЛ НЕРАВНОМЕРНЫМИ ГИПЕРЗВУКОВЫМИ ПОТОКАМИ ГАЗА ПОД УГЛАМИ АТАКИ И СКОЛЬЖЕНИЯ

3.1. Постановка задачи и метод решения.

3.2. Решение задачи при обтекании тела потоком типа дальнего следа

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

4 МЕТОД ШЕСТОГО ПОРЯДКА АППРОКСИМАЦИИ ДЛЯ РЕШЕНИЯ КРАЕВЫХ ЗАДАЧ ПАРАБОЛИЧЕСКОГО ТИПА

4.1. Математическая постановка задачи.

4.2. Метод решения.

4.3. Пример применения метода: уравнение Фокнера-Скэн.

4.4. Решение уравнений трехмерного ламинарного пограничного слоя на проницаемой поверхности.

5 ПРИМЕНЕНИЕ МЕТОДОВ ВЫСОКОГО ПОРЯДКА ТОЧНОСТИ ДЛЯ РЕШЕНИЯ ЭЛЛИПТИЧЕСКИХ УРАВНЕНИЙ

5.1. Метод установления для решения уравнений вязкого ударного слоя.

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

5.3. Метод переменных направлений для схемы высокого порядка аппроксимации.

6 ПРИМЕНЕНИЕ ГЕНЕТИЧЕСКИХ АЛГОРИТМОВ К РЕШЕНИЮ ЗАДАЧ ОПТИМИЗАЦИИ В ГИПЕРЗВУКОВОЙ АЭРОДИНАМИКЕ

6.1. Метод решения.

6.2. Оптимизация формы тела по тепловому потоку.

6.3. Оптимизация затупления тела в рамках различных газодинамических моделей.

6.4. Оптимизация траектории входа в атмосферу земли по интегральному тепловому потоку.

7 ПРИМЕНЕНИЕ МНОГОПРОЦЕССОРНЫХ ВЫЧИСЛИТЕЛЬНЫХ КОМПЛЕКСОВ К РЕШЕНИЮ ЗАДАЧ ГИПЕРЗВУКОВОЙ АЭРОДИНАМИКИ.

7.1. Параллельный алгоритм для решения пространственных уравнений ТВУС.

7.2. Моделирование течений в пространственном химически неравновесном вязком ударном слое на многопроцессорных вычислительных системах с распределенной памятью

7.3. Расчет пространственного турбулентного пограничного слоя на сети компьютеров различной производительности

7.4. Метод встречной прогонки и его параллельная реализация

7.5. Асинхронный параллельный генетический алгоритм

Введение диссертация по механике, на тему "Исследование течений в вязком ударном слое при помощи схем высокого порядка аппроксимации"

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

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

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

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

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

Как отмечалось в [6]-[7], при движении тела на верхнем участке траектории характерные числа Рейнольдса относительно невелики, и классическая модель разбиения возмущенной области на зону невязкого течения и пограничный слой становится неприменимой, поскольку эффекты молекулярного переноса оказываются существенными во всей области течения между ударной волной и поверхностью тела,

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

Уравнения ПВУС, предложенные в работе [9], являются, как показывает анализ [10]-[13], эллиптическими в дозвуковых областях течения, и поэтому для их решения невозможно прямое применение маршевых методов. В связи с этим для их решения требуются значительные машинные ресурсы, особенно при решении трехмерных задач. В то же время в целом ряде практически важных случаев обтекания тела при больших сверхзвуковых скоростях (М^ 1) можно использовать более простую газодинамическую модель течения - ТВУС. Эта модель была предложена в [29] на основании анализа плоского и осеоимметричного гиперзвуковых обтеканий тел с непроницаемой поверхностью в предположении Мое —)• оо, Reoo —> оо, г == (7 — 1)/(7 + 1) —>• 0, sRe ~ 0(1). Уравнения, описывающие течение в ТВУС, содержат все члены уравнений пограничного слоя Прандтля и гиперзвукового невязкого ударного слоя ([2]) и являются, с математической точки зрения, уравнениями параболического типа. Трехмерные уравнения ТВУС были получены Э.А.Гершбейном [30], уравнения нестационарного трехмерного тонкого вязкого ударного слоя получены в [31]. Как следует из проведенного в обзорах [6], [7] анализа многочисленных сравнений решений уравнений ТВУС с данными эксперимента и с расчетами, выполненными в более полной математической постановке, теория ТВУС дает хорошую точность получаемых решений всякий раз, когда вязкий ударный слой является тонким.

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

Первые приближенные аналитические решения задачи обтекания сферы в рамках модели ТВУС получены в [29], [33],[34], численное решение для различных режимов обтекания гиперболоидов получено в [32], течение в ТВУС с учетом вдува газа с поверхности исследовано в [35] ,[36].

Первые решения уравнений трехмерного тонкого вязкого ударного слоя были получены для течения в окрестности точки торможения. Эта точка является особой точкой системы уравнений ТВУС, причем, вследствие параболичности этих уравнений, решение в окрестности критической точки может быть получено независимо от решений в остальных областях, а система исходных уравнений в критической точке сводится к системе обыкновенных дифференциальных уравнений, зависящих только от одного геометрического параметра - отношения радиусов главных кривизн в критической точке. Решение уравнений ТВУС в окрестности критической точки тела с проницаемой поверхностью рассмотрено в [38],[39], эта же задача, но без учета масоообмена на поверхности была рассмотрена в [37]. Обтекание затупленных тел с проницаемой поверхностью под углами атаки и скольжения исследуется в [40]-[42]. При этом используется эффективный численный метод четвертого порядка точности по нормальной к поверхности тела координате, позволяющий проводить расчеты при отсутствии в течении плоскостей симметрии.

В первых трех главах диссертации в рамках модели трехмерного тонкого вязкого ударного слоя исследуется обтекание затупленных тел неравномерными потоками газа. Возникающая при этом картина течения имеет характерные особенности, которые делают необходимым детальный учет эффектов неравномерности в набегающем на тело потоке газа. Характерным примером неравномерного обтекания является обтекание моделей летательных аппаратов сильно недорасширенными струями в экспериментальных установках с целью получения достаточно высоких тепловых потоков, которые были бы близки к натурным [43], [44]. В ряде случаев, например, при использовании плазменных сверхзвуковых источников и отсутствии сопел с жесткими стенками, хорошей математической моделью таких течений является течение от сверхзвукового сферического источника. В работе [45] показано, что при истечении сверхзвуковой струи в вакуум на больших растояниях от среза сопла течение также близко к течению от источника и приведены формулы для определения параметров на источнике через характеристики реального сопла. Первые работы, в которых было рассмотрено обтекание затупленных тел неравномерными набегающими потоками газа относятся к середине 60-х годов: в работах [46], [47] в рамках уравнений Эйлера рассматривалось обтекание поверхности сильно недорасширенными струями, в работах [48], [49] исследовалось невязкое обтекание сферы расширяющимся потоком газа. При этом течение перед скачком уплотнения предполагалось течением от сферического источника [48], либо расчитывалось методом характеристик [49]; в ударном слое решение находилось методом интегральных соотношений Дородницына-Белоцерковского. В работе [50] в рамках гиперзвукового невязкого ударного слоя получено приближенное аналитическое решение для течения в окрестности критической точки тупого тела при его обтекании расходящимся потоком. При этом было показано, что влияние неравномерности на решение можно сильно ослабить, если при анализе задачи ввести эффективный радиус кривизны, зависящий от степени неравномерности потока. Задача об ударе неравномерного сверхзвукового потока невязкого газа о плоскую преграду рассматривалось в [51].

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

Обтекание тонких тел гиперзвуковым расширяющимся потоком невязкого газа рассмотрено в [53], [54]. Получена аналогия между обтеканием параллельным потоком и потоком от источника для случая, когда показатель адиабаты 7 = 2.

Задача обтекания затупленного тела потоком от источника с учетом эффектов молекулярного переноса в ударном слое в рамках уравнений Навье-Стокса и классического разделения возмущенной области на зону невязкого течения и пограничный слой рассмотрена в [55]. Эта же задача, но в рамках уравнений ПВУС была рассмотрена [56]-[58],[99]. При этом для численного интегрирования использозался метод, описанный в [21], [22]. Результаты проведенных расчетов показали, что усиление неравномерности (уменьшение растояния от тела до источника) приводит к уменьшению отхода ударной волны в критической точке, смещению звуковой линии к оси симметрии, распрямлению ударной волны и увеличению теплового потока на всей лобовой поверхности сферы. В [58] исследовались свойства подобия таких течений. Были получены универсальные распределения напряжения трения и теплового потока вдоль поверхности сферы в переменных работы [52] (в качестве независимой переменной использовалась rj = в/вв, где 9s - угловая координата звуковой точки. В рамках модели пограничного слоя данная задача рассматривалась в [59]-[62], причем решение в невязком ударном слое находилось с использованием метода Г.Г.Черного [2] в виде разложения искомых функций по параметру е, где £ - отношение плотностей на ударной волне, а коэффициенты трения и теплообмена на поверхности тела расчитывались по формулам [67], полученным с использованием теории локального подобия Лиза [68]. В [60] предложена приближенная формула для определения теплового потока через параметры неравномерности и тепловой поток в случае равномерного обтекания, которая, как показали сравнения с [56], имеет точность не хуже 10%. В работе [63] было рассмотрено обтекание сферического затупления расходящимся потоком газа при наличии сильного дозвукового вдува с поверхности в рамках невязкой постановки задачи. Теплообмен при взаимодействии недорас-ширенной струи с плоской преградой теоретически и экспериментально рассматривался в [64]- [66].

Другой интересный с практической точки зрения случай неравномерного обтекания возникает, если тело находится в следе за другим телом. Исследование течения в следе за телом, двяжущимся со сверхзвуковой скоростью, является классической задачей аэродинамики и нашло широкое отражение, например, в [4], [69]-[71]. Обтекание тела, находящегося в ближнем следе за другим телом, достаточно хорошо изучено экспериментальным образом (например, в работах [72], [73], [82]). С другой стороны, экспериментальные исследования моделей, летящих в дальнем следе за другим телом, весьма затруднены. Поэтому при исследовании таких течений на передний план выступают теоретические методы.

Ранее течение около поверхности сферы, находящейся в дальнем сверхзвуковом следе исследовалось в рамках уравнений Навье-Стокса в [55] и в рамках уравнений ПВУС в [56], [74]- [78], [99]. Результаты проведенных расчетов свидетельствуют о том, что с ростом интенсивности следа происходит увеличение отхода ударной волны в критической точке и размеров дозвуковой области, а при достаточно сильной неравномерности набегаюцего потока происходит качественная перестройка характера течения, связанная с образованием в ударном слое возвратно-циркуляционной зоны. В [74], [76] было рассмотрено обтекание сферы, летящей в следе за цилиндром, обтекание которого было расчитано в рамках уравнений Навье-Стокса в [5]. В [56], [74], [77] при расчете течения в ударном слое использовались также уравнения Эйлера.

Вопросы теплообмена и сопротивления тел, летящих в дальнем следе, рассмотрены в [79], [81]. При этом решение в невязком ударном слое находилось при помощи асимптотического метода [87], а параметры трения и теплообмена на поверхности расчитывались по формулам из работы [67]. В [79] также предложен асимптотический критерий, при выполнении которого осуществляется безотрывное обтекание. В ряде работ данная задача решалась в рамках модели тонкого вязкого ударного слоя (ТВУС). В [80] при помощи метода из работ [33]-[34] для малых значений числа Рейнольдса получено приближенное аналитическое решение, а в [85] асимптотическое решение в невязком ударном слое и слое вдуваемых газов. Численные решения задачи для течения около непроницаемой поверхности получены в [82], [83], для течения с учетом вдува с поверхности тела в [82], [84], около удлиненных гиперболоидов вращения в рамках турбулентного ТВУС - в [94].

В [88], [89] при помощи сеточно-характеристического метода расчитано течение около лобовой поверхности сферы, эллипсоида и плоского торца при их обтекании неравномерным потоком идеального газа, в [90] в предположении о наличии в ударном слое области невязкого течения и пограничного слоя расчитано обтекание затупленных конусов под углами атаки потоком типа дальнего следа.

При изучении движения тела через температурную неоднородность также необходимо учитывать эффекты неравномерности в набегающем потоке газа. Характерной особенностью такой задачи является её существенная нестационарность. Ранее некоторые нестационарные задачи течения в вязком ударном слое рассматривались в [31], [92] - в этих работах исследовались течения при нестационарном вдуве газа с поверхности тела. В ряде работ проводилось теоретическое исследование нестационарных эффектов, возникающих при сверхзвуковом движении тела в неоднородных средах. В [93]-[94] изучение велось в рамках модели идеального газа, в [96]-[99] для расчета обтекания сферы использовались уравнения вязкого ударного слоя, причем если в [96] решение ограничивалось окрестностью линии торможения, то в [97]-[99] оно было получено для всей лобовой поверхности. Приведенное в [97]-[99] сравнение с экспериментом [95] показывает удовлетворительную точность выбранной математической модели.

В работах [101]-[105] в качестве примера неравномерного набегающего потока рассматривалось закрученное течение газа. В [100] закрученное течение в гиперзвуковом ударном слое исследовалось в невязкой постановке. В рамках модели вязкого ударного слоя влияние внешней завихренности потока на теплообмен около непроницаемой поверхности тело в окрестности критической точки исследовалось в [101], а с учетом вдува с поверхности тела в [105]. В последней работе было получено также асимптотическое решение в невязком ударном слое, анализ которого позволил сделать вывод о возникновении в ударном слое при достаточно сильной неравномерности возвратно-циркуляционной зоны. В [102]-[103] закрученные течения изучались в рамках модели ламинарного пограничного слоя, в [104] -в рамках параболизованных уравнений Навье-Стокса.

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

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

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

В заключительном параграфе данной главы исследуется течение многокомпонентного химически реагирующего газа в окрестности критической точки тела, движущейся вдоль траектории входа в атмосферу Земли перспективного аппарата NASA Х38. Уравнения вязкого ударного слоя решаются совместно с уравнениями теплопроводности внутри покрытия и с учетом гетерогенных реакций на поверхности тела. Исследование вязкого ударного слоя в данной постановке на теплоизолированной поверхности рассматривалось в [116]-[118], в сопряженной постановке - в [119]-[121]. При этом во всех вышеперечисленных случаях рассматривалась траектория ДА предыдущего поколения (Space Shuttle) [122]. Траектория NASA Х38[123] обладает значительными особенностями, связанными с возможностью использования современных теплозащитных покрытий, позволяющих осуществлять торможение JIA на более низких высотах полета. Приводится температура поверхности в зависимости от высоты (времени) полета для различных каталитических активностей поверхности и различных теплофизических свойств покрытия. Показано, что необходимость учета влияния сопряженного теплообмена на температуру поверхности существенным образом зависит от характеристик покрытия ДА. Основные результаты этой главы опубликованы в работах [173], [175], [177], [180].

В данной работе получено численное решение в широком диапазоне определяющих параметров задачи, а в случае закрученного набегающего потока - для малых и умеренных чисел Речнольдса получено также приближенное аналитическое решение интегральным методом последовательных приближений. Исследовано совместное влияние врашения тела и неравномерности набегающего потока на характер течения в ударном слое, коэффициенты трения и теплообмена на поверхности тела. В этой же главе рассмотрено неравномерное обтекание удлиненных гиперболоидов вращения и обобщен невязкий критерий безотрывного обтекания из главы 1 на течение около боковой поверхности тела. Основные результаты этой главы опубликованы в работах [174], [176], [178], [182].

В третьей главе диссертации рассматривается задача о течении однородного газа в пространственном ТВУС около затупленных тел, обтекаемых неравномерными потоками типа дальнего следа и потока от сверхзвукового сферического источника под углами атаки и скольжения. При этом используется эффективный численный метод, обладающий четвертым порядком аппроксимации по поперечной координате и не требующий наличия в течении плоскостей симметрии. Проанализировано влияние геометрических характеристик течения, числа Рейнольдса и температуры поверхности на характер влияния параметров неравномерности набегающего потока на интегральные характеристики течения в ударном слое. Основные результаты этой главы опубликованы в работах [179], [181], [183].

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

В настоящее время имеется большое количество работ, посвященных использованию разностных схем применительно к расчету течений в пограничном слое. Достаточно подробный обзор таких работ изложен в [124],[125]. При этом наряду с традиционными схемами второго порядка аппроксимации по поперечной координате[126],[127], имеются работы, в которых предлагаются схемы четвертого порядка точности в указанном направлении [128]-[131].

Наиболее простой способ построения схем высокого порядка состоит в использовании многоточечных шаблонов. Однако, этот метод не на шел достаточно широкого применения как из-за сложности решения получающихся конечно-разностных уравнений, так и из-за проблем, возникающих при нестандартной записи схемы вблизи границы. Альтернативный способ построения схем высокого порядка аппроксимаций заключается в использовании эрмитовых или компактных аппроксимаций, главная идея которых заключается в построении разностных формул, включающих в себя как искомые функции, так и их производные. Некоторая аналогия этому подходу просматривается в теории интерполяции функций, когда при использовании полиномов Эрмита для увеличения точности аппроксимации в опорных точках кроме значений функций учитываются еще и значения их производных. Различные вопросы применения эрмитовых схем к решению разнообразных задач гидродинамики рассмотрены в обзоре [132]. Систематическое изложение теоретических основ компактных аппроксимаций с разностями против потока представлено в монографии [133].

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

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

Полученные при помощи данного метода решения модельного уравнения Фокнера-Скэн приведены в сравнении с результатами, полученными с помощью известных методов Келлера[127] и Петухова[128], обладающих соответственно вторым и четвертым порядками аппроксимации. В настоящей работе представлены также результаты применения этого метода к решению достаточно сложной и актуальной задачи - расчету течений в трехмерном ламинарном пограничном слое около затупленных тел с проницаемой поверхностью. Для широкого диапазона изменения определяющих параметров задачи проведено систематическое и детальное исследование влияния числа разбиений и типа используемой разностной сетки (равномерной либо неравномерной) на точность получаемого решения как в окрестности критической точки, так на боковой поверхности тела. Проанализирован характер сеточной сходимости на последовательности вложенных сеток и дана сравнительная оценка вычислительной эффективности полученного алгоритма и расчетов, полученных на основе методов Келлера[127] и Петухова[128]. Основные результаты этой главы опубликованы в работах [192], [194].

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

Метод глобальных итераций для решения двумерных задач впервые предложен в [9]. В работах [12]-[16] был предложен целый ряд усовершенствований метода: был изменен способ задания начальных данных, применены методы релаксации и сглаживания формы ударной волны. Все это позволило расширить класс тел, к которым мог быть применен метод глобальных итераций. Для решения трехмерных уравнений ПВУС этот метод впервые был применен в [17] при расчете обтекания сферически затупленных конусов под углами атаки.

В методе установления в качестве исходных используются нестационарные уравнения ПВУС, а искомое решение находится как предел при устремлении временной переменной к бесконечности. В [18], [19] этот метод используется для исследования течения совершенного газа в плоскости симметрии сферически затупленных конусов. При этом нестационарные уравнения решались явным методом Мак-Кормака [20]. В работах [21],[22] для решения нестационарных уравнений ПВУС в задаче обтекания сферы использовалась неявная по времени разностная схема с весами, имеющая второй порядок точности по пространственным координатам. Полученная таким образом система нелинейных разностных уравнений решается методом последовательных приближений, предложенным ранее в [23] для расчета течений невязкого газа. В [25] метод [21],[22] был обобщен на случай пространственного обтекания эллипсоида вращения под углом атаки. Обтекание тел более сложной формы при нулевом угле атаки рассматривалось в [26], а при ненулевом угле атаки в [27],[28].

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

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

В заключительном параграфе данной главы предложен метод переменных направлений для решения двумерных эллиптических уравнений, аппроксимируемых при помощи схем типа схемы Петухова и показано преимущество этого метода по сравнению с традиционным методом релаксации по строкам на примере уравнений Бюргерса. Здесь же приведен новый вывод конечно-разностных соотношений, полученных в оригинальной работе Петухова [128] из интегральных соображений, позволяющий выявить место этой схемы в ряду других компактных аппроксимаций. Основные результаты этой главы опубликованы в работах [189], [196], [199].

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

Предложенные в конце 60-ых годов, и получившие теоретическое обоснование в 1975 году [135] генетические алгоритмы (в силу своей универсальности и высокой эффективности) в дальнейшем получили чрезвычайно широкое распространение для решения задач поиска и оптимизации в самых разнообразных областях науки и техники [136], [138], [139]. ГА представляют собой алгоритмы поиска, основанные на механизмах природной селекции и генетики, и имеют дело с популяцией возможных решений задачи, которые подвергаются циклическому воздействию трех основных "генетических" операторов - селекции, скрещиванию и мутации. Будучи вероятностными, генетические алгоритмы тем не менее не являются просто еще одним вариантом случайного поиска, поскольку при отборе новых точек с ожидаемыми более хорошими возможностями, они эффективно используют предыдущую информацию.

ГА не являются единственно возможными методами нулевого порядка (т.е. не использующими градиентную информацию). Близкие идеи развиваются в работах, посвященных эволюционный методам и эволюционным стратегиям [137], а также управляемому случайному поиску [134] и некоторым другим методам, имеющим дело с наборами испытаний (возможных решений) [142].

В диссертации в рамках моделей пограничного слоя и вязкого ударного слоя исследуются формы затуплений, на которых реализуется минимум интегрального теплового потока к поверхности тела при заданных габаритных ограничениях и ограничении на максимум равновестной температуры поверхности. Для решения собственно задачи оптимизации затупления лобовой поверхности осесимметричного тела по конвективному тепловому потоку предлагается эффективный генетический алгоритм (ГА). На основании систематических численных расчетов сделан вывод, что предложенный метод решения вариационной задачи является устойчивым, робастым, имеет высокий уровень вычислительной эффективности и позволяет получать решение в широком диапазоне изменения параметров задачи. Проведено сравнение решений, получающихся при использовании различных газодинамических моделей. Оптимальное решение, полученное на основе ГА сравнивается с решением, полученным с использованием традиционного подхода с применением модифицированного метода Давидона для поиска экстремума нелинейной функции многих переменных. Показано, что ГА позволяет достигнуть более лучших результатов: интегральный тепловой поток для решения с использованием ГА меньше, чем соответствующее значение для градиентого метода на 3-5%.

Ранее задача об определении тела минимальных тепловых потоков была рассмотрена с использованием ряда упрощающих предположений. В [143] была поставлена задача об определении оптимального радиуса затупленного тела при заданном угле полураствора конуса и его объеме, в [144] для аналогичной задачи в расчет принималась излучаемая телом энергия. В работе [145] класс рассмотренных тел ограничивался формами с плоскими затуплениями, а в [146] был определен радиус сферического затупления и угол полураствора конуса минимальных тепловых потоков при заданных длине и миделе тела. В [147, 148] была поставлена и решена задача об определении формы тела с минимальным радиационным тепловым потоком к его поверхности. Приближенное решение аналогичной задачи в классе тонких тел для минимального полного (конвективного и радиационного) теплового потока было получено в [149, 150]. Обзор работ посвященных решению экстремальных задач при теплообмене тел, движущихся с гиперзвуковыми скоростями, а также постановка вариационной задачи для нахождения оптимальной формы пространственных тел, обладающих минимальным суммарным тепловым нагревом при движении вдоль траектории входа, приведен в [151].

В заключительном параграфе данной главы рассматривается задача оптимизации параметров траектории входа в атмосферу Земли по величине суммарного конвективного теплового потока в точке торможения затупленного тела. Как известно, при входе тела в плотные слои атмосферы Земли тепловые потоки к его поверхности настолько велики, что при создании космических аппаратов (КА) многоразового использования возникает сложная и комплексная задача их тепловой защиты. При этом решение задачи о снижении уровня тепловой нагрузки на КА подразумевает определение параметров траектории входа и выбор конструктивного решения К А (аэродинамические характеристики, тип системы теплозащиты и т.д.). Решение полной задачи оптимизации траектории по тепловой нагрузке на поверхность КА в общем виде включает в себя совместное решение задачи оптимального управления (механическая задача) для системы обыкновенных дифференциальных уравнений движения материальной точки в атмосфере [152, 153, 154, 155] в которую входят аэродинамические коэффиценты, требующие для своего определения расчетов пространственного обтекания всей поверхности КА (газодинамическая задача) при наличии физико-химических процессов (термодинамическая и химическая неравновесность, излучение и т.д.)

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

В настоящей работе приведены результаты решения задачи по определению траектории входа затупленного тела заданной формы в атмосферу Земли, вдоль которой суммарный тепловой поток в критическую точку тела имеет минимум, а равновесная температура поверхности тела не превышает наперед заданную величину, а также траектории, вдоль которой максимальная равновесная температура поверхности тела в критической точке минимальна. В качестве исходной газодинамической модели для определения теплового потока к поверхности тела используются уравнения тонкого (гиперзвукового) вязкого ударного (ТВУС) слоя с учетом химической неравновесности, многокомпонентного характера диффузии и сопряженного теплообмена на поверхности, обладающих хорошей точностью для условий входа в диапазоне высот от 100 до 50 км и скоростей от 8 до 2,5 км/с [7]. Показано, что решение этой задачи может быть сведено к решению последовательности задач без учета сопряженного теплообмена. Основные результаты этой главы опубликованы в работах [193], [198], [200] - [203].

Последняя глава диссертации посвящена построению некоторых параллельных алгоритмов гиперзвуковой аэродинамики. Интерес к пара-лельным методам вычислительной гидродинамики резко вырос за последнее десятилетие и связан с широким распространением сравнительно недорогих многопроцессорных систем [1], обладающих производительностью супер-компьютера и позволяющих уменьшить стоимость вычислений в несколько раз. С другой стороны, благодаря развитию компьютерных технологий, очень быстро растет производительность рабочих станций и персональных компьютеров. Некоторые эксперименты [157] показывают, что сетевой комплекс из 14 компьютеров IBM RS6000 оказывается быстрее чем суперкомпьютер Сгау-2. В настоящее время большинство научных учреждений обладает сетями достаточно производительных компьютеров, которые могут эффективно использоваться для решения вычислительных задач механики сплошной среды. Програмное обеспечение для этих сетей использует обычные принципы проектирования параллельных алгоритмов для MIMD-компьютеров. Несмотря на быстрый рост числа публикаций в этой области, необходимо отметить что существуют целый ряд разделов вычислительной гидродинамики, для которых пока не разработаны достаточно эффективные параллельные алгоритмы. Ряд последних достижений в этой области приведен в [162]-[165].

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

В разделе 7.4 предложен вариант метода встречной прогонки для схемы Петухова и на основе его параллельной реализации построен алгоритм для решения уравнений полного вязкого ударного слоя с эффективностью распараллеливания до 90 % и более. Ранее параллельный вариант метода, близкого к встречной прогонке для обычной трехдиа-гональной матрица приведен в [166], оригинальный подход к "распараллеливанию" прогонки представлен в [167].

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

Основные результаты этой главы опубликованы в работах [185]-[191], [203].

Основные результаты диссертации докладывались на конференциях молодых ученых ИТПМ (г.Новосибирск, 1987,1988 гг.), на IV, V Всесоюзных конференциях по современным проблемам аэродинамики (г. Жданов, 1987г., п.Рыбачий, 1989 г.), на Школе молодых ученых по численным методам механики сплошной среды (г.Абакан, 1989 г.), на 7 Всесоюзной школе-семинаре "Современные проблемы механики жидкости и газа" (г.Иркутск, 1990 г.), на III Всесоюзной школе-семинаре по макрокинетике и химической и магнитной газодинамике Красноярск, 1990 г.), Школе по математическому моделированию (Абрау-Дюрсо, 1993 г.), конференциях "Применение математического моделирования для решения задач в науке и технике". (г.Ижевск. 1996, 1998гг.), Международной конференции по внутрикамерным процессам и горению ICOC-96 (г.С.-Петербург, 1996 г.), международных конференциях Parallel Computational Fluid Dynamics (г.Париж, Франция, 1993г., г.Киото, Япония, 1994 г., г.Манчестер, Великобритания, 1997 г., г.Вильямсберг, США, 1999 г.), конференции "Всесибирские чтения по математике и механике" (г. Томск, 1997 г.), конференции "Математические модели и методы исследования" (г.Красноярск, 1997 г.), международной конференции по генетическим алгоритмам EVROGEN97, (г.Триест, Италия, 1997 г., ), конференции "Фундаментальные и прикладные проблемы современной механики" (г. Томск, 1998 г.), Четвертом европейском конгрессе по вычислительной гидродинамике ECCOMAS'98 (г.Афины, Греция, 1998 г.)

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

Заключение диссертации по теме "Механика жидкости, газа и плазмы"

ЗАКЛЮЧЕНИЕ

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

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

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

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

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

10. На основе принципа геометрического параллелизма предложены подходы к созданию численных алгоритмов с высоким уровнем эффективности для решения трехмерных уравнений тонкого вязкого ударного слоя и турбулентного пограничного слоя на основе неявной конечно-разностной схемы высокого порядка точности на сети компьютеров с различной производительностью либо на многопроцессорной ЭВМ с распределенной памятью. Предложен вариант метода встречной прогонки и его адаптация для параллельной вычислительной системы на примере решения уравнений полного вязкого ударного слоя. Приведен асинхронный параллельный генетический алгоритм для решения задач оптимизации. Показано, что во всех случаях эффективность распараллеливания достигает 90% и более.