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

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

СИСТЕМЫ МОМЕНТНЫХ УРАВНЕНИЙ И СЛЕДУЮЩИЕ ИЗ НИХ МОДЕЛИ НЕРАВНОВЕСНЫХ ТЕЧЕНИЙ

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

2 2 И10Л 2015

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

005570854

Москва-2015

005570854

Работа выполнена в Московском авиационном институте (Национальном исследовательском университете)

Научный консультант

действительный член РАН, доктор технических наук, профессор Ю.А.Рыжов Официальные оппоненты:

Галкин Владлен Сергеевич, доктор физико-математических наук, Центральный аэрогидродинамический институт им. Н.Е.Жуковского, главный научный сотрудник.

Козлов Виктор Владимирович, доктор физико-математических наук, профессор, Институт теоретической и прикладной механики им. С.А.Христиановича СО РАН, зав. лабораторией "Аэрофизических исследований дозвуковых течений".

Кузнецов Михаил Михайлович, доктор физико-математических наук, профессор, Московский государственный областной университет, профессор кафедры "Теоретической физики".

Ведущая организация: Институт теплофизики им. С.С.Кутателадзе СО РАН.

Защита состоится 30 октября 2015г. в 10 часов на заседании диссертационного совета Д 212.125.14 по защите докторских и кандидатских диссертаций при Московском авиационном институте, г. Москва, Волоколамское шоссе, д. 4

С диссертацией можно ознакомиться в библиотеке МАИ

Автореферат разослан «

Ученый секретарь

диссертационного совета Д 212.125.14 кандидат физико-математических наук

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

Актуальность темы.

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

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

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

Объект исследования - течения высокой степени динамической неравновесности.

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

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

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

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

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

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

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

Два упомянутых метода, по существу, альтернативны. Метод Чепмена-Энскога, включая его обобщенные версии, принципиально ограничен значениями Кп < 1. Система дифференциальных уравнений содержит только уравнения сохранения, т.е. неполную систему моментных уравнений второго порядка. Численная реализация получаемых моделей существенно затруднена

в виду высокой нелинейности дифференциальных уравнений, начиная со второго (барнеттова) приближения.

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

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

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

Приведенные доводы определили направление и содержание настоящей работы.

Цель работы - развитие метода моментных уравнений и расширение его на течения многоатомных газов. Основные задачи работы.

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

- Распространение метода на течения многоатомных газов.

- Разработка метода снижения коротковолновой неустойчивости системы моментных уравнений.

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

- Исследование свойств моделей первого и второго приближения, следующих из системы моментных уравнений.

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

- Разработка инженерных моделей неравновесных течений многоатомного однокомпонентного газа.

- Численное тестирование разработанных моделей. Методология и методы исследования.

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

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

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

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

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

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

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

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

На защиту выносятся:

- метод построения системы моментных уравнений для функции распределения общего вида;

- методы снижения коротковолновой неустойчивости моментной системы;

- модель граничных условий на твердой поверхности, не имеющая ограничений по числам Кнудсена;

- внепорядковая и гибридная инженерные модели. На защиту, также выносятся положения:

- об основной причине коротковолновой неустойчивости моментных систем;

- о физической неадекватности температур поступательных и внутренних степеней свободы, определенных 5-моментными системами уравнений. Научная и практическая значимость работы.

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

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

Достоверность результатов исследования подтверждена сравнением полученных расчетных данных с данными экспериментальных исследований разных авторов.

Апробация и внедрение результатов.

Результаты работы были представлены на следующих семинарах:

- видеосеминар по аэромеханике ЦАГИ - ИТПМ СО РАН - СПбГТУ -НИИМ МГУ, 4 февраля 2014г.;

- семинар Сектора кинетической теории газов Вычислительного центра им. А.А.Дородницына, 3 марта 2015г.;

- международный авиационно-космический научно-гуманитарный семинар им. С.М. Белоцерковского,19 марта 2015г..

Результаты работы включены в курс лекций по дисциплине "Динамика неравновесных сред", читаемый на кафедре "Аэродинамика ЛА" Московского авиационного института. Публикации.

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

Из указанных публикаций 10 статей и монография опубликованы лично автором.

Структура и объем работы.

Диссертация состоит из введения, семи разделов, заключения, списка сокращений и условных обозначений и списка использованных источников. Общий объем составляет 252 страницы, включая 51 рисунок. Список использованных источников содержит 106 наименований.

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

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

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

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

Рассматриваются центрированные моменты вида:

=т0$с1с^к.../с1сс1е и М?^ = \ес1с;ск.../с[с<1£. В этих выражениях:

та, е - масса и энергия внутренних степеней свободы молекулы; С1 -проекция тепловой скорости молекулы; / =

одночастичная функция распределения; N - порядок момента. Основные моменты в формальном и газодинамическом обозначениях: МР*РЦ; = Еп; Л/£>-2%; Общее моментное уравнение порядка N >2:

---------. . д

Э* дха г]к-' дха р дха дха

(1)

w(N)+

Здесь = --и.и.\У^:2)ГсЩс1е-...

- быстрота наработки момента в результате обратных столкновений; -

момент, вычисленный по распределению столкнувшихся молекул; г -среднее время свободного пробега молекулы; 1/г+ - средняя частота обратных столкновений. В уравнениях моментов, содержащих энергию внутренних степеней свободы, символ М'~ должен быть заменен символом

Повторяющийся в одночлене греческий индекс обозначает операцию свертки тензора. Подстрочный символ "звездочка" обозначает операцию симметрирования (компактного симметрирования) с исключением повторений тождественных по величине одночленов.

Уравнения моментов второго порядка:

1р+МиР)+Р Ь. Рц

д/ и дха дха дха

(2)

+ Э (и р т —¿п+-—+ ----1--^

с)/ оха дха т т т

Здесь 2 - быстрота передачи энергии от поступательных степеней свободы к внутренним.

Тензор напряжений может быть разложен на сферический и бездивергентный тензоры следующим образом:

Ри=р^+6ирт (4)

Здесь р™ - неравновесные напряжения; рт=Раа/3 - "механическое" давление.

С учетом (4) правая часть уравнения (2) может быть представлена как:

Ру =Р^ + 5иРт+ + Т+ т Т+ г (5)

Коэффициентом П обозначено отношение неравновесных напряжений, создаваемых молекулами после и до столкновений, т.е. П = /р™ .

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

Выражение г/(1-г/г+ п), фигурирующее в (5), имеет смысл времени

релаксации напряжений:

(6)

(

! т Рц

Уу у

В отношении релаксационного процесса принято еще одно положение: в процессе релаксации энергия теплового движения молекул не переходит в энергию их группового движения. На основании этого положения члены уравнений (3) и (5), отвечающие за энергообмен между поступательными и внутренними степенями свободы, обозначены единым символом Е: _ 1 [ г „Л _ 3 1Г т г „

Условие энергообмена в терминах температуры:

Здесь сп = Су - 3/2 /? - теплоемкость внутренних степеней свободы; су -изохорная теплоемкость; Я - удельная газовая постоянная. Поступательная (Г,) и внутренняя (Тп) температуры связаны с термодинамической температурой (7) соотношением:

сутЛв.Т,+саТа. (8)

Из (7) в совокупности с выражениями 5 следует:

2=3рцса 9 =Ъ (9)

2 Су Тд 2 2 Тд где 9 = Т,-Тп\ тв - время релаксации разности температур, определенное аналогично времени релаксации напряжений (6); у - показатель адиабаты.

Время релаксации те представлено в виде тд=кг . В настоящей

работе величина А рассматривается как свободный параметр модели или оценивается в соответствии с известными рекомендациями.

Уравнения моментов третьего порядка (рцк и щ следуют из общего

моментного уравнения (1). Вид релаксационных членов этих уравнений и времена релаксации т(р и та получены аналогично (5), (6).

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

м$ = ++ГЛ )- К2Р ; (10)

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

(12)

М$5) = спрТпУ/,;к+ЯрТ„0.. (13)

М$тп=Я2рТ.,Т„Т„. (14)

Здесь использованы переменные Ту = Ру/(рЯ), у/..к=2Ф^/Р и = /р,

отражающие структуру полученных выражений лучше, чем системные переменные Р^, /рик и щ.

Анализ выражений (10) + (15) с учетом основных свойств центрированных моментов позволил принять допущение: моменты высших порядков могут быть адекватно аппроксимированы линейной комбинацией моментов низших порядков. Под моментами низших порядков понимаются моменты, определенные системой, требующей замыкания.

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

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

Для определения времен релаксации тр, гр, тш принято допущение:

время релаксации неравновесной величины не зависит от степени неравновесности течения. Это допущение позволило выразить указанные величины посредством закона трения Стокса и закона Фурье. В результате получены следующие значения: тр = //(Г,)/рт ; тр = 3/2 Тр; та = тр. В последнем выражении использована аппроксимация Эйкена Рг = 4у/(9у- 5).

На основании принятых положений и допущений сформулирован метод построения системы моментных уравнений порядка для функции

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

1) Записываются моментные уравнения со второго до УУ-го порядка с использованием зависимости (1).

2) Правые части уравнений приводятся к виду релаксационного члена. Время релаксации неравновесной величины типа (рЦ,<Рук,Щ, •••) выражаются

поступательными и внутренними степенями свободы записывается с учетом

как

Член, описывающий энергообмен между

3) Время релаксации моментов второго и третьего порядка определяются с использованием закона трения Стокса и закона Фурье. Время релаксации моментов порядка N>2 принимаются в качестве свободных параметров модели.

4) Для моментов (// + 1)-го порядка записываются зависимости вида (10) ^ (13).

Описанным методом построена система моментных уравнений третьего порядка, называемая в дальнейшем 24-моментной системой или системой М24. Одна из форм записи, удобная для много процессорной численной реализации:

Об)

а! оха

*&. + А.{Ригв + Р1в) = 0; (17)

<я оха

д Р + Э и, Р ч, Эй)а_3(у-1) в . ,

Ч --^--д«~ГкрТГр' (19)

Э Э / \ Эи. Л Э7». 2 (¡5/Л

Э/ 7 Эдга и дха 2 дха 3

Э Э / ч ди, _ Э — ¿У,- + —К«;)+ «« ~ + За о? оха ОХд оха

£0.а (21)

Р У тр

Свойства 24-моментной системы в отношении коротковолновой неустойчивости такие же, как и у 20-моментной системы Трэда. На Рисунке 1 показаны профили скорости и температуры в плоской ударной волне одноатомного газа.

2.4 -

1.6 -

1.2 -

2 -

20

24

28

32

Рисунок 1. Профили скорости и температуры в ударной волне одноатомного газа. Система М24.

При М>1.8 возникает физически неадекватное искажение профиля (субскачок). На профилях ударной волны в многоатомных газах субскачок проявляется аналогичным образом.

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

Показано, что локальные выражения замыкающих моментов (10), (11) не удовлетворяют их моментным уравнениям (балансовым соотношениям).

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

М$ = КТ,.Я, + т

мп(а) -г т Р + жп(4) му - сп]п'и + ти

У

(22)

(23)

В этих выражениях ) и /и? ^ выполняют функции согласующих добавок.

Используемая методика выделения согласующих добавок аналогична методике "регуляризации" 13-моментной системы Грэда.

Подстановка (22) и (23) в соответствующие моментные уравнения позволяет получить уравнения согласующих добавок:

т1]Ы

+ и„

тцИ

дхп

тЩа Зы. Р Ъха

+ ЩТ.

ЭдгЛ

дТ.<

дх„

--Тг

р» х 5-3/ Кв

--Ь О*+ *--

2 Лт

тцк1

рТ4

(24)

4}

+ и

а ^

Эх,,

/я.

П(4)

Эы«

-

зть

Эх„

+ Л

Эх„ Эх„ ) 2 2 г„ /г

+ ^ 5-Зу ДА + а 2

рТ

__ти

£1(4)

.(25)

РГП4

В правых частях уравнений фигурируют времена релаксации моментов четвертого порядка г4 и ТП4, определяемые в соответствии с описанным выше методом построения системы моментных уравнений.

Система М24, дополненная зависимостями (24) и (25), образует систему 45-ти скалярных уравнений, условно называемую системой М45. Область течений, свободная от коротковолновой неустойчивости, у системы М45 несколько шире, чем у системы М24. На Рисунке 2 показаны профили скорости и температуры в ударной волне одноатомного газа, рассчитанные с помощью системы М45.

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

Ослабление влияние локальных выражений моментов пятого порядка на решения системы и сокращение ее до 24-х уравнений, будет достигнуто,

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

3.5-,

О 4 8 12

Рисунок 2. Профили скорости и температуры в ударной волне одноатомного газа. Система М45; Т/Тр= 0.5, т4/тр= 1.

Если рассматривать время релаксации как малую величину и принять Ту ~ 8уТ(, то выделение старших членов (24), (25) приводит к выражениям:

ук! х р Р

2Л

Эр... _1_Эр_

Эх, *** р Эх,

ВТ

а /

Г р

\ Р

Р„ _ 5 — ЗуЯр9 -

г

v р

2 А г

р /

+ 5-ЪуЯрв + " 2 Ит

г

v р

р /

; (26)

у т р

р

эг зс г Й61

2спе>.. —+—---—

" У" Эх 2 с т ¡г

а V р

Ри

Ярв

г ч с кт

V р

г г

+ К

Э®, 1 Э р

■СО,

\ V

Эх, * р Эх,

дТ + 8..0) — V адх

а )

Уравнения моментов третьего порядка с учетом (22), (23) примут вид:

^ й. \Ът% _ 2 <Рцк .

(28)

(29)

Остальные уравнения системы М24 остаются неизменными. Полученная система уравнений в дальнейшем называется системой М24с. Ее коротковолновую неустойчивость демонстрирует Рисунок 3.

Рисунок 3 Профили скорости в ударной волне одноатомного газа. Система М24с; т/Тр =0.612; т4/гр = 1.

При М„ > 3 профиль начинает «делиться» на два участка, что не соответствует природе явления. Субскачок в характерном для него виде не возникает даже при гиперзвуковых числах Маха.

Дальнейшее ослабление влияние локальных выражений моментов пятого порядка на решения системы будет достигнуто, если приближение Ту ~ буТ, дополнить приближением фук = [8^+ 5Л(р^ + буфк)[ 5. Уравнения (26) и (27) примут следующий вид:

Р

Эх.

Эх,

г

--т.

Эх„

рТ. 1 § 5-3у Яв рт_ " 2 Ит

1 5-Зг ( дТ . ЭГ 5 * ЭхГ + " Эх^

р

\

; (30)

р) )

„дш. г дТ + Т— + 5,/Ва — +

Эх, охп

35-3у г 6> + 2 2 трИ

г

■ (31)

р» I 85~ЪгКв " 2

Система М24с, в которой уравнения (28), (29) записаны с учетом (30), (31), в дальнейшем обозначается как М24сс. Профили скорости, рассчитанные с помощью системы М24сс представлены на Рисунке 4.

и

1 Х/А.О0

О 4 8 12 16 20

Рисунок 4. Профили скорости в ударной волне одноатомного газа. Система М24сс; г^ р =0.612; г4/тр=1.

Согласующие добавки (26), (27) приводят к появлению вторых производных <р и со. в их моментных уравнениях (28), (29). Известным

свойством вторых производных является сглаживание решения. Кривые Рисунка 3 сглажены вторыми производными моментов <р и со..

Согласующие добавки в виде (30), (31) не содержат момента <р, , и будучи

ук

подставленными в (2В) не повышают порядка дифференциального уравнения. Таким образом, сглаженный вид профилей Рисунка 4 следует отнести к достоинствам метода согласования.

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

Аппроксимирующая функция распределения системы М24 не содержит разложения в ряд, тем не менее, на профиле ударной волны, рассчитанной с помощью системы М24, при М„ =1.8 возникает субскачок, см. Рисунок 1. Характер коротковолновой неустойчивости системы М24 полностью соответствует характеру 20-моментной системы Трэда.

Согласование балансового и локального выражений замыкающих моментов позволяет исключить субскачок, см. Рисунок 2. Ослабление влияния локального выражения моментов пятого порядка, присутствующих в зависимостях замыкающего момента, в еще большей степени снижают коротковолновую неустойчивость, см. Рисунок 3 и Рисунок 4. Материалы Раздела 2 приводят к следующим выводам:

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

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

В третьем разделе получены первое и второе приближения системы

М24.

Линейные комбинации уравнений (16), (18), (19) позволяют записать уравнения температур:

Э Т д Т , ди„ 1 да„

а! дха охр сур оха

+ н + д"а I 2 д(Ра= 5~3у в ;

ЭI " дха 3 Ъхп ЗЛрдха 2 Ьт '

Э7Ь Э/

1 да)а 3(7-1) в

дТа

дха спр дха

Ит„

Два последних уравнения дают уравнение разности температур:

Ъв

Ъв

+ -Т,

ди„ 1 Э

3/ а дха 3 дхп р дха

2 1 —<ра--соа

ЗЛ'

са

в

~1гтп

(32)

(33)

(34)

(35)

Линейная комбинация уравнения (19) с его сверткой приводит к уравнению:

^ Р\! + -V. \»агц Г "а

Э и.

дх„ у и ' " дх„ '' 3 дх„ Эх,

Первые приближения (в смысле процедуры Чепмена-Энскога) системы М24 могут быть получены в предположении того, что время релаксации напряжений тр является малой величиной, по отношению к величине

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

2к V диа. 3 Яр дха '

фЬ = <Рша =

^Эх, у 3 дха

1Ч = ГЛ[.] =_Цр„К.

4 ^ Эх,- '

15

с

Рг 4

М

Э7Ь

ЭХ;

(37)

(38)

(39)

(40)

В этих выражениях надстрочные квадратные скобки обозначают порядок приближения. Коэффициент вязкости определен как ц — тррт.

В механике сплошной среды тензор напряжений разложен как: Ру = рд + ^р, (41)

где ру - термодинамические неравновесные напряжения; р = рКТ -термодинамическое давление. С учетом (4), (8) и (41)

= (42)

В первом приближении:

(43,

Это выражение содержит коэффициент объемной вязкости р' в явном виде: > 5-3у,

(44)

Первое приближение теплового потока:

(45)

Рг Э^-

Полученные уравнения представляют собой строгие (не содержащие

дополнительных допущений) первые приближения неравновесных величин,

следующие из системы М24. Эти уравнения в совокупности с уравнениями

сохранения массы и импульса позволяют построить две неэквивалентные

системы уравнений первого приближения.

Система, содержащая 5 моментных уравнений и называемая в

дальнейшем 5-моментной системой или 5-моментной моделью, имеет вид:

Ър | Ъриа _ д Эг дха

Эм; Э и, 1 дРа „ ,, ^

+ + = (46)

о/ оха р дха

дт

а т

диг

+ т~ + \r-Waa +

Э/ ' "" Эх„ ' ч/ ар дхр ' с,р дха

"Эх. ~у 3^* 2 "Т'а*/ Рг Эх,-

Температуры Т, и Тп в приближении, соответствующем приближению системы (46), могут быть определены следующим образом:

+ = (47)

'2 3 Яр дха

= г + (48)

2 Яр оха

На Рисунке 5 показаны профили температур в ударной волне

двухатомного газа. Значение параметра к составляло к = 2.5 (5 = 1 - модель

максвелловских молекул).

Рисунок 5. Профили температур в плоской ударной волне двухатомного газа. 5-моментная система; М„ = 5; 5 = 1; к = 2.5.

На профиле температуры внутренних степеней свободы имеет место физически неадекватная область. Этот модельный эффект наблюдается при всех значениях к> 1 и подробно рассмотрен в диссертационной работе. В

частности показано, что первого приближения разности температур в не достаточно для определения температур Т\ и Тп, в том же приближении. Это

положение распространено на модели высших приближений.

Термодинамическая температура определена качественно верно. Если

расчет поля течения ограничивается только этой температурой, то 5-

моментная модель может рассматриваться как модель, уточняющая

классическую модель НСФ в отношении коэффициента объемной вязкости.

Система, содержащая 6 моментных уравнений (двухтемпературная

модель), также является строгим первым приближением системы М24.

Эр | Эриа _ 0

Э/ Ъха

д и, Эм; 1 Э Р,а Л Э/ аха р дха

+ и + ЭЦ/? I 2 д(Ра= Э/ * дха 3 ар дха 3рЯ дха 2/г тр

ЭГп , ЭГП ( 1 Ъсоа _Ъ{Г-\)Т,-Та Э/ " дха спр дха 2А Тр

(49)

Рц ~ ~М

Г ~ди, ^ 2 ди ^ У дх, 1 3 дха у

15 дГ дТп

Здесь Р. = р™ + б..рт, р=р[т5). Коэффициент объемной вязкости, а

следовательно и первое приближение разности температур, не используется.

На Рисунке 6 совмещены профили температур 5-моментной и двухтемпературной моделей. Параметры ж и А моделей подобраны таким образом, чтобы наклон профиля плотности соответствовал экспериментальным данным Алсмейера (см. список литературы диссертации, [92]).

Несмотря на завышенную по сравнению с предыдущими тестами величину параметра А = 3.1, двухтемпературная модель обеспечивает качественно верное решение и дает более широкую, чем 5-моментная модель, область возмущений.

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

- ✓ _ ___

-

- ТЛЧ/ // К / '/ / /

-

-

• ^¿^!

- N / N /

- х/Я, 1 1 1 1 1 1

23 32 36 40

Рисунок 6. Профили температур в плоской ударной волне двухатомного газа; М„ = 5. Сплошные линии - двухтемпературная модель, 5 = 1, А = 3.1. Пунктирные линии - 5-моментная модель; 5 = 1, А = 4.

Рассмотрено второе (в смысле процедуры Чепмена-Энскога) приближение системы М24. В качестве малой величины принято среднее время свободного пробега молекулы Т. Строгое второе приближение системы соответствует модели Барнетта с коэффициентом объемной вязкости в виде (44).

Для получения смешанных производных групповой скорости и температуры процедура Чепмена-Энскога предусматривает

дифференцирование соответствующих уравнений модели Эйлера по координате геометрического пространства. Одной из особенностей уравнений Барнетта является наличие в них посторонних решений. Можно

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

Во избежание дополнительного дифференцирования уравнение неравновесной величины, например (35), может быть перегруппировано и представлено в следующем виде:

Л (д0 дв \

-А г,

3, V V

+

„/л

РаР Яр

диа | 1 Э дхл р дха

ЗЛ

Ра —

_1_

сп

-со„

(50)

Уравнения остальных неравновесных величин преобразуются аналогично.

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

3 Яр \Рвф ЭХр ' дха у

I сп

(51)

Здесь в квадратных скобках указан номер приближения, а символ с обратной стрелкой представляет собой первое приближение в более ранний момент времени. Для в справедливо уравнение:

= (52)

в +

ГВ9 д§ кд( дха;

р " 3""' дх„

Полная система уравнений модели второго приближения (релаксационной модели):

+ = о;

Э/ дха

61 дха р дха

ВТ дТ , _ Рав

Эг<

Эх» с^р дха

Э0 Э<9

1

д( "а Эх„ г„

ч [Эху- Эх, у 3 Эха

• г;'7'

+ Ру

= 1 (\5 ЪТ [ ' Э/ " Эха V 4 Эх,- '

Эй>,- Эй), 1 Г ЭГ _ Э/ дха тю I Эх,-

В этой системе:

-+

8 (у-1)Рг" спЯ'"2 4 + ])Рг'

а, =

63

!ч = —1-%—1--аь — — с л;

3 8 (/-1)Рг спЯ 4 2 л

( д2т, | ац эт эт 5 сх э Г эг^

дx¡дxj с1Т Эх,- дXj у 3 Эха ^ дха J

+

Г

ди. Эг/а + Эма Эм. ^ 0 Эм,- Э^

+ 2:

Эх. Эха Эх, Эха Эха Эха

/ /

^Эх^ Эха

Эи^ дхг.

+ --С

'ДЭхл

_ _ (з д (р

г ди,- Эи„ Л | + а

Э и„ э т

г, —=---1- а2 —— + —-

дха Эх,- ^ Эх,- дха

ди

дха Эх,-Эм,-

дт

Эх„

Эх,- [ р г

+

Эы„

У у

3 Э2м,- Э 2и ' +а4

у 2 дхадха

ЭхаЭх,-

\ >

Т

У У

Основные свойства полученной системы уравнений исследованы на примере решения задачи о профиле ударной волны. На Рисунке 7

представлены профили плотности, скорости и температуры при М„ = 5. Значения р, их, Т преобразованы к единичному отрезку.

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

0.8 0.6 0.4 0.2 0

5 10 15 20 25 30 °°

Рисунок 7. Профили плотности, скорости и температуры в плоской ударной волне одноатомного газа. Релаксационная модель, М„ =5, 5 = 1

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

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

В геометрическом пространстве использована декартова прямоугольная система координат Ох¡х2х3, связанная с обтекаемой поверхностью. Ось Ох1

принята в качестве внешней нормали поверхности. Скоростная система координат Ос]с2с2 связана с центром массы жидкой частицы.

Принят диффузный закон отражения молекул от поверхности с полной аккомодацией энергии. Функция распределения отраженных от поверхности молекул:

(2жКТ„У2

Здесь Тк - температура поверхности; с,- = с; + и{. Надстрочным символом "тильда" обозначены величины, относящихся к отраженным молекулам.

Падающие на поверхность молекулы условно разбиты на группы медленных и быстрых молекул. Медленные молекулы имеют малые по модулю скорости с,. Составляющие скорости с2 и с3 не ограничены.

Медленные молекулы испытывают последнее столкновение вблизи поверхности. Составляющие температуры Т22 и этих молекул приняты равными температуре газа Т, в граничной точке, а их групповая скорость -равной скорости газа и в той же точке. Составляющая температуры Тц обозначена Г., концентрация - й0. Функция распределения медленных молекул:

/о --(55)

(2 тгКу^Т.^Е 2КТ<

Надстрочный символ "крышка" относится к падающим молекулам.

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

2 п с2

/ =-Т/ ехр(-—с;=С;~Щ, й,=0. (56)

(2 лЯТУ2 2К1

Вторая функция "вырезает" (знак "минус") из первой медленные молекулы:

2 пп

-ехр(—-

С2 + ¿3 )

(57)

Для того чтобы все молекулы, имеющие скорость с1 = = 0, были исключены из распределения /, необходимо выполнение следующего условия:

/_(с1=0) = -/(с1=0), (58)

откуда следует: Т.

(59)

Параметр а, значение которого находится в пределах 0 < а < 1, принят в качестве свободного параметра модели. Он позволяет исключить и0 и Г.:

й0 = ап, Т, = а Т.

Функция распределения/в граничной точке представлена как:

/ =

/

при с, >0;

(60)

(61)

[/+/.+/» прис,<0. На Рисунке 8 моделирующие функции представлены схематически линиями уровня в пространстве тепловых скоростей. С1

Рисунок 8. Схема расположения моделирующих функций распределения в пространстве тепловых скоростей

Функция распределения в граничной точке содержит восемь неизвестных величин: п, п, и2, «3, й2, й3, Тп Т (сг - свободный параметр). Для их определения могут быть записаны пять скалярных уравнений сохранения и уравнения для рх2, р13, (¡\.

Система восьми уравнений может быть получена интегрированием функции распределения / с соответствующими молекулярными скоростями. Решение этой системы:

(я 1 + г? , ...

и,=-,—-х-Ри, 1 = 2,3,

У 2 г?+(1 + ст)г?2

(62)

(3 - 2а)(г92 -1) + г?(1 + )и'2

Д Г =--. (63)

3г? + 3-2ст

Здесь:

..' _ иг . \Т'- ■ п - Р]1 1 = 23-

я =

в

-; г9 = „|— =2Л ^(-агссоз! - ,

2ЯГ г ^ ^

7Г

1 г-^+о.би^-^АГ-^^ д =_^_

V 2-а2+0.5^Н'2 ' 2-<т2+0.5^и'2'

1-сг 1-СГ

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

Оценка параметра сг0 проведена по модели граничных условий с 13-моментной функцией распределения падающих молекул. Сравнение результатов расчетов, полученных в условиях близких к равновесию по двум моделям, дает значение а0 = 0.47.

Параметр а представлен в следующем виде:

°-47 V ' Л (ЛАЛ

о —-г; Кп = —. (64)

1 + /ЛСп I

Здесь Я - средняя длина свободного пробега; Ь* - характерный размер возмущенной области течения, вычисленный по нормали к поверхности; /3 — свободный параметр модели, значение которого подобрано в Разделе 6.

В пятом разделе построены инженерные модели неравновесных течений.

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

Р = грРт.

В (36) член, отвечающий за наработку напряжений р™, имеет вид

Вязкие свойства модели будут улучшены, если коэффициент вязкости (тензорный коэффициент вязкости) определить как р^ = Тр Ру.

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

Ми = тр(биРт+ХР^). (66)

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

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

ги =к,гм + {\-к,)г[2-\ (67)

Здесь г", гм, г[2'] - неравновесная величина в гибридном, моментном и внепорядковом (неполном втором) приближениях. Весовой коэффициент определен как: к1 = , где т] - параметр модели (77 = 2 и- 3), IV- степень

неравновесности течения:

У_1

К^ф-г1-?-• (68)

л/6"2 Г Р

На Рисунке 9 показаны профили скорости в ударной волне двухатомного газа при Мм = 5. Свободные параметры всех моделей выбирались из соображений совпадения наклона профиля плотности с экспериментальными данными (см. список литературы диссертации, [92]).

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

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

Рисунок 9. Профили скорости в ударной волне двухатомного газа, М^ = 5. Пунктирная линия - кинетическая модель; 1 - гибридная модель (5.57), у = 0.73, /г = Л(з), /7 = 2, ^ = 0.5; 2 - внепорядковая двухтемпературная модель, 5 = 0.98, /г = /2(5), ^ = 1; 3 - двухтемпературная модель, 5 = 1, /г =3.1; 4 - 5-моментная модель, 5 = 0.82, /г = 7.3.

И'1

Ли,-=-1.431 _-0.224^-, / = 2,3; (69)

W^7» рт

срЫ

Л Г. = -0.735 ^ (70)

Здесь координата х,, фигурирующая в напряжении р"'^ и тепловом потоке ^.рассматривается как внешняя нормаль к обтекаемой поверхности.

Сравнение расчетных и экспериментальных данных позволяют уточнить постоянные множители уравнений (69) и (70) применительно к течениям многоатомных газов. С этой целью рассмотрены граничные условия в виде:

«[i] [i]

'-2-3; <71)

Коэффициенты £иТ и рассматриваются как свободные

параметры.

В задаче о теплопередача в плоском слое неподвижного газа рассмотрена зависимость отношения теплового потока к его

свободномолекулярному пределу цт от числа Кнудсена. Результаты расчетов по различным моделям представлены на Рисунке 10.

Наилучшую согласованность с экспериментальными данными дает феноменологическая модель граничных условий и условия скольжения (72) при = 0.55.

Рисунок 10. Тепловой поток в плоском слое газа. Газ - воздух, Тк1/Т„2 = 1.2. О - экспериментальные данные. 1 - модель НСФ с граничными условиями прилипания; 2, 3, 4 - М24, релаксационная, 5-моментная и двухтемпературная модели; 2 - феноменологические граничные условия; 3 -граничные условия (70); 4 - граничные условия (72), Сг ~ 0-55.

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

ел,

1

0.8 0.6 0.4 0.2

0

0.01 0.1 1 ю Кп

Рисунок 11. Зависимость коэффициента трения от числа Кнудсена в плоском течении Куэтта. Газ - воздух. М = 1.4. О - экспериментальные данные; 1 - модель НСФ с граничными условиями прилипания; 2, 3, 4 -модель М24; 2 - феноменологические граничные условия; 3 - граничные условия (69), (70); 4 - граничные условия (71), (72), = 1.21, - 0.55; 5 -5-моментная модель, граничные условия (71), (72), £"„ = 1.21, = 0.55.

Граничные условия скольжения при £,=1.21, = 0.55 и Кп<10 обеспечивают удовлетворительные значения коэффициента трения. При больших значениях числа Кнудсена эти граничные условия дают коэффициент трения, превышающий свободномолекулярный предел. Расчеты показывают, что при М = 4 превышение свободномолекулярного предела

наблюдается уже при Кп ~ 1. Феноменологические граничные условия

лишены этого недостатка. Применительно к течению Куэтта при Кп —»°°

феноменологические граничные условия соответствуют точной свободномолекулярной функции распределения.

В задаче о структуре ударной волны анализировалась обратная ширина профиля (максимальный наклон профиля) и форма профиля плотности. На Рисунке 12 приведены экспериментальные (см. список литературы диссертации, [91; 92; 97]) и расчетные данные обратной ширины профиля плотности в ударной волне аргона.

и/5

Рисунок 12. Обратная ширина ударной волны в аргоне. Знак + -экспериментальные данные; 1 - модель первого приближения (НСФ) 5 = 1; 2 - релаксационная модель М„ < 5, 5 = 0.9; 3 - внепорядковая модель, 5 = 0.9, %=\\ 4 - внепорядковая модель, 5 = 0.82, ,£ = 2.3; 5 - гибридная модель 5 = 0.8.

Кривая 1 соответствует 5-моментной модели, совпадающей в данном случае с моделью НСФ. В сверхзвуковой области чисел Маха даже при 5 = 1 профиль слишком крутой, в гиперзвуковой области — слишком пологий. Удовлетворительный наклон профиля дает релаксационная модель (кривая 2).

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

На Рисунке 13 расчетные значения обратной ширины профиля плотности в азоте сопоставлены с экспериментальными данными (см. список литературы диссертации, [92; 98]). Кривая 1, соответствующая модели НСФ без коэффициента объемной вязкости, проходит далеко от области экспериментальных точек. Если рассматривать к как свободный параметр модели, то при 5 = 0.82 и к = 13 5-моментная модель (кривая 3) дает удовлетворительный наклон профиля.

0 6 _ /б

0.4 -

0.2-

2

4

6

8

10

Рисунок 13. Обратная ширина ударной волны в азоте. Знак + -экспериментальные данные [92; 98]; 1 - НСФ 5 = 1; 2 - 5-моментная модель 5 = 1, к =2.5; 3 - 5-моментная модель 5 = 0.82, Л = 7.3; 4 -двухтемпературная внепорядковая модель 5 = 0.81, й = А(5), ; ,£ = 2.3 5 -гибридная модель 5 = 0.8, А = А(5); 6 - релаксационная модель М00<6, 5 = 0.83, Л = /г(5)

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

Сравнение экспериментальных и расчетных профилей плотности приведены в диссертации.

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

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

Для анализа результатов расчета использованы экспериментальные данные (см. список литературы диссертации, [99]), полученные для обтекания тонкого клина. Расчеты проводились для двухатомного газа при = 10.15 и температуре поверхности пластины ТК=2.\6Т^. Начало прямоугольной

системы координат совмещено с носиком пластины. Ось Ох направлена вдоль пластины, ось Оу является нормалью к поверхности.

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

Руу/Р«,

Рисунок 14. Распределение нормального напряжения при у = О, М„ =10.15, =2.16. О - экспериментальные данные; 1 - 5-моментная модель; 2 - внепорядковая модель; 3 - гибридная модель; 4 - нормальное напряжение на поверхности пластины при свободномолекулярном обтекании.

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

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

распределение продольной составляющей скорости их в этой области. Вблизи носика локальные числа Маха составляют: 7.58 - 5-моментная модель; 3.21 - внепорядковая модель; 2.29 - гибридная модель.

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

- ч

-

- \\

- 3^ \ \

■ \

- 11111111

-4 -3 -2 -1 О

Рисунок 15. Распределение скорости их при у = 0. Область торможения потока перед носиком пластины. 1 - 5-моментная модель; 2 - внепорядковая модель; 3 - гибридная модель.

ЗАКЛЮЧЕНИЕ

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

Указанным методом построена замкнутая система моментных уравнений третьего порядка (24-моментная система). Тесты показали, что 24-моментная система содержит основной "дефект" моментных систем -коротковолновую неустойчивость. В остальном эта система дает результаты, удовлетворительно согласующиеся с экспериментальными данными.

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

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

3. Показано, что первое приближение 24-моментной системы уравнений приводит к двум неэквивалентным моделям течения многоатомного газа: к модели Навье-Стокса-Фурье с коэффициентом объемной вязкости (5-моментная модель) и к двухтемпературной модели.

Показано, что в условиях высокой неравновесности течения 5-моментная модель определяет поступательную и внутреннюю температуры

качественно неверно. Качественно верное определение температур в первом приближении дает только двухтемпературная модель.

Установлена причина неэквивалентности моделей. Показано, что первого приближения разности температур Т, - Та не достаточно для определения самих температур в том же приближении.

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

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

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

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

7. Разработана гибридная модель неравновесного течения. Модель лишена коротковолновой неустойчивости и позволяет получить физически адекватное поле течения в высоко неравновесной области.

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

СПИСОК ПУБЛИКАЦИЙ ПО ТЕМЕ ДИССЕРТАЦИИ

Статьи в журналах из перечня ВАК:

1. Никитченко, Ю.А. Модель течения газа вблизи поверхности раздела сред / Ю.А.Никитченко // Математическое моделирование. - 2003. - Т. 15. -№8.-С. 88-98.

2. Никитченко, Ю.А. Сравнение различных моделей течения газа в широком интервале чисел Кнудсена / Ю.А.Никитченко // Математическое моделирование. - 2004. - Т.16. - № 8. - С. 77-93.

3. Никитченко, Ю.А. Система моментных уравнений многоатомных газов / Ю.А.Никитченко // Полет. - 2010. - № 11. - С. 43-51.

4. Никитченко, Ю.А. Применение модели Навье-Стокса-Фурье к расчету гиперзвукового обтекания тонкой пластины / Ю.А.Никитченко // Вестник МАИ. - 2011. - Т. 18. - № 3. - С. 21-28.

5. Никитченко, Ю.А. Снижение коротковолновой неустойчивости системы моментных уравнений за счет ее расширения / Ю.А.Никитченко // Ученые записки ЦАГИ. - 2015. - Т. 46. - № 1. -С. 72-84.

6. Никитченко, Ю.А. Оценка эффективности методов распараллеливания вычислений в задачах динамики разреженных газов / Ю.А.Никитченко, Ю.А.Рыжов, С.А.Попов, О.Ю.Власов // Научный вестник МГТУ ГА.-2015.-№211.-С. 85-92.

7. Никитченко, Ю.А. Моментные модели для течений с большим числом Маха / Ю.А.Никитченко // Вестник МАИ. - 2014. - Т.21. - № 4. - С. 39-48.

8. Никитченко, Ю.А., Модели первого и второго приближений для течений многоатомных газов / Ю.А.Никитченко // Вестник МАИ. - 2012. - Т. 19,-№2.-С. 11-17.

9. Никитченко, Ю.А. Модели первого приближения для неравновесных течений многоатомных газов / Ю.А.Никитченко // Электронный журнал Труды МАИ. - 2014, - № 77.

10. Никитченко, Ю.А. Феноменологическая модель граничных условий на твердой поверхности / Ю.А.Никитченко // Вестник МАИ. - 2012. - Т. 19. -№3. - С. 5-14.

11. Никитченко, Ю.А. Инженерная модель неравновесного течения / Ю.А.Никитченко // Изв. Вузов. Авиационная техника. - 2014. - № 3. - С. 3740.

12. Рыжов, Ю.А. Гибридная модель гиперзвукового течения / Ю.А.Рыжов, Ю.А.Никитченко, С.А.Попов // Изв. Вузов. Авиационная техника. - 2015. - № 1. - С. 7-11.

13. Рыжов, Ю.А. Численное исследование гиперзвукового обтекания острой кромки на основе модели Навье-Стокса-Фурье / Ю.А.Рыжов, Ю.А.Никитченко, И.В.Парамонов // Электронный журнал Труды МАИ. -2012.-№55.-9 с.

Монография

Никитченко, Ю.А. Модели неравновесных течений: монография / Ю.А.Никитченко. - М.: Изд-во МАИ, 2013. - 160 с.

Учебное пособие

Свирщевский, С.Б. Кинетические методы в аэрогазодинамике: учебное пособие / С.Б.Свирщевский, Ю.А.Никитченко. - М.: Изд-во МАИ, 2001. - 63 с.

Множительный центр МАИ (НИУ) Заказ от. о? 201-5" г. Тираж (00 экз.