USER
Напиши подглаву 2.2.1. для реферата (без формул, только суть на 2-3 страницы А4): Министерство Науки и высшего образования и РФ
Федеральное государственное бюджетное учреждение высшего образования Тульский государственный университет
Кафедра РВ
Реферат
На тему: «Построение пространственной сетки, и её использование в ракетостроении »
Выполнили студенты группы 132331:
Ткачёва В.В.
Колпаков Ю.О.
Новиков И.А.
Проверил:
Ветров В.В.
Тула 2025
Содержание
Введение.................................................................................................................3
Теоретические основы построения сеток
1.1. Понятие пространственной сетки, ее виды........................................5
1.2. Основные методы построения сеток................................................6-7
1.3. Качество сеток и критерии их оценки............................8
Применение расчетных сеток в задачах ракетостроения
2.1. Подготовка данных
2.2. Пространственные сетки в реальных расчетах
2.2.1. Использование в аэротермодинамике..............................9-13
2.2.2. Использование при прочностных расчетах....................14-15
2.2.1. Использование для моделирования процессов горения в ракетных двигателях..................................................................16-17
Заключение.........................................................................................................18
Список литературы.............................................................................................19
Введение
Развитие ракетно-космической техники требует применения передовых методов численного моделирования для решения задач проектирования сложных конструкций и анализа их поведения в экстремальных условиях. Одним из ключевых инструментов в этом процессе является построение расчетных сеток, которые обеспечивают дискретизацию геометрических моделей и физических процессов.
Объект исследования: Численное моделирование процессов в летательных аппаратах.
Предмет исследования: Построение и использование пространственных сеток.
Цели исследования:
Проанализировать современные методы построения пространственных расчетных сеток, используемых в численном моделировании сложных геометрий ракетной техники.
Рассмотреть преимущества и ограничения различных видов сеток в контексте задач ракетостроения.
Изучить применение различных типов сеток для решения задач аэродинамического анализа, прочностного моделирования и теплового расчета ракетных и космических аппаратов.
Оценить влияние качества построенной сетки на точность численного моделирования.
Продемонстрировать примеры практического использования современных методов построения сеток в ключевых задачах ракетостроения.
//В настоящей работе рассматриваются актуальные подходы к генерации структурированных, неструктурированных и гибридных сеток, их преимущества и ограничения, а также их применение в задачах ракетостроения, включая аэродинамический анализ, прочностное моделирование и расчет тепловых нагрузок. Особое внимание уделяется современным методам, таким как триангуляция Делоне, метод продвигаемого фронта и адаптивные криволинейные сетки, которые повышают точность и эффективность моделирования. - хороший кусок текста, но пока не знаю, куда впихнуть. Наверное, пойдет в заключение//
Теоретические основы построения сеток
1.1. Понятие пространственной сетки, ее виды.
Пространственная сетка - это математическая модель, описывающая поверхность или весь объем некоторого объекта, представляющая из себя набор чаще всего одинаковых геометрических плоских или объемных фигур (например, треугольники или тетраэдры), каждая из которых будет наделена своими физическими свойствами, необходимыми для правильного описания физических процессов происходящих во всем теле.
Сетки делятся на несколько видов, каждый из которых имеет свои преимущества и области применения в ракетостроении:
Структурированные сетки: элементами являются кубы или параллелепипеды, которые располагаются упорядоченно. Они просты в генерации и эффективны для расчетов, но плохо подходят для сложных геометрий из-за возникающих больших погрешностей. Используются для моделирования внутренних компонентов ракеты, таких как топливные баки.
Неструктурированные сетки: элементами являются обычно тетраэдры, расположенные хаотично, что делает их гибкими для сложных форм корпуса ракеты. Однако их генерация сложнее, а расчеты будут занимать больше времени. Применяются, например, для аэротермодинамического моделирования.
Гибридные сетки: комбинируют структурированные и неструктурированные зоны, что оптимально для задач с разными масштабами, например, моделирования пограничного слоя и ударных волн вокруг ракеты. Но процесс создание таких сеток и программы для них сильно сложнее.
1.2. Основные методы построения сеток.
В настоящее время самыми распространенными методами для построения сеток являются методы, использующие треугольники и тетраэдры, а именно:
1. Ускоренная триангуляция Делоне - самый распространенный метод из всех, который имеет свои вариации (алгоритм с кэшированием поиска треугольников, алгоритм “Разделяй и властвуй” и др.).
Рисунок
Рис.1., триангуляция Делоне
Использует в основе разработку советского математика Б.Н.Делоне, заключающуюся в заполнение области треугольниками (в объеме - тетраэдрами), при которой для любого треугольника все точки из плоскости, за исключением точек, являющихся его вершинами, лежат вне окружности, описанной вокруг треугольника.
Метод довольно требовательный к производительности, но показывающий очень хорошие результаты построения.
Метод продвигаемого фронта - метод, использующийся сейчас в большинстве программ для построения сеток, приходящий на смену методу Делоне с целью повышения точности описания процессов при сокращении расчётного времени. Для метода так же разработаны различные вариации алгоритмом, отличающиеся качеством и временем построения сеток.
Различные методы на основе адаптивных криволинейных сеток - новые методы
Рисунок
Рис.2. Пример метод продвигаемого фронта
Суть метода заключается в построении равнобедренных треугольников с определённой площадью, начиная с ребра минимальной длины. Процесс повторяется до заполнения всей площади многоугольника треугольниками. Также обсуждаются корректировки для случаев, когда вершина треугольника близка к вершине многоугольника или пересекает его ребро. Для областей со сложным течением процесса предлагается измельчать треугольные сетки для более точного моделирования. Для эластичных границ предлагается отсекать увеличенную подобласть и строить для неё сетку отдельно, сохраняя конформность с остальной частью области.
Методы на основе адаптивных криволинейных сеток - различные новые подходы к построению сеток, использующие генерацию точек, создание треугольников и адаптацию сетки с использованием прямолинейных и криволинейных операций обмена рёбер. Но данные методы все еще проигрывают по соотношению качество-время с двумя предыдущими, но продолжают активно дорабатываться командами программистов и математиков.
1.3. Качество расчетных сеток и критерии их оценки.
Основная проблема при построении пространственных сеток заключается в том, что при дискретизации поверхностей трехмерного тела необходимо учитывать требования к качеству треугольников, точности и отсутствию пересечений между элементами сетки. Поверхностные триангуляционные сетки должны полностью описывать поверхность тела с заданной точностью; отклонение от геометрии свыше установленной нормы может негативно повлиять на достоверность расчетов.
При построении сетки следует учитывать искажение размеров, перекрытия, самопересечения, а также несогласованные разбиения границ. Необходимость точного представления геометрии оправдывает использование BREP (Boundary Representation) – аналитического представления корпуса, где каждая оболочка состоит из набора граней, которые, в свою очередь, описываются набором формул.
Применение расчетных сеток в задачах ракетостроения
2.1. Подготовка данных.
Подготовка исходных данных для построения расчетных пространственных сеток, предназначенных для разбиения сложных объектов на тетраэдрические элементы, является одним из ключевых и базовых этапов в процессе численного моделирования физических процессов. Этот этап определяет качество последующих расчетов и напрямую влияет на точность моделирования процессов, происходящих в ракетно-космических системах.
Процесс подготовки данных включает в себя несколько важных последовательных этапов, каждый из которых требует внимательного выполнения и контроля:
Во-первых, необходимо получить корректную геометрическую модель исследуемого объекта — ракеты или ее составной части. Геометрия может быть создана с нуля с помощью программ для проектирования или импортирована из сторонних источников в формате, поддерживаемом программным обеспечением для генерации расчетных сеток. Наиболее часто применяемые форматы включают STL, STEP и IGES. Выбор формата зависит от программного обеспечения и требований к точности модели.
Во-вторых, важно проверить геометрию на наличие ошибок, таких как пересечения поверхностей, незамкнутые контуры или слишком мелкие детали, которые могут затруднить процесс разбиения.
В-третьих, необходимо задать параметры сетки (размер элементов и др.). Выбор оптимальных параметров сетки требует понимания специфики моделируемых процессов, а также особенностей программного обеспечения.
2.2. Пространственные сетки в реальных расчетах.
2.2.1. Использование в аэротермодинамике. на основе статьи: Приведены результаты численного моделирования обтекания сверхзвуковым потоком идеального газа спускаемого космического аппарата "Fire II" при входе в атмосферу Земли. Моделирование проводилось на основе численного решения системы уравнений Навье-Стокса на неструктурированных сетках в трехмерной постановке с учетом сжимаемости газа.
Методами вычислительной аэродинамики исследованы особенности структуры полей течения и теплофизические процессы во всей возмущенной области вокруг спускаемого аппарата от головной ударной волны до дальнего следа, поведение местных аэродинамических и тепловых характеристик, а также структура сложных трехмерных вихревых течений у поверхности аппарата в условиях невозмущенного набегающего потока, соответствующих различным точкам возможных траекторий входа спускаемого аппарата в атмосферу Земли. Вычисления проведены в диапазоне изменения углов атаки а от 0 до 20°.
В рамках ряда моделей турбулентности исследовано влияние угла атаки а на поле течения, вихревые приповерхностные течения и местные аэродинамические и тепловые характеристики вдоль передней (носовой) и задней поверхностей. Результаты расчетов сопоставлены с данными летных испытаний, лабораторных экспериментов и приведенными в литературе данными расчетно-теоретического анализа.
Ключевые слова: численное моделирование, аэротермодинамика, газовая динамика, теплообмен, сверхзвуковые течения, спускаемые космические аппараты, уравнения Навье-Стокса, турбулентность, неструктурированные сетки.
При движении в атмосфере спускаемые космические аппараты (КА) подвергаются большим тепловым и динамическим нагрузкам, и при проектировании перспективных спукаемых аппаратов большое значение имеют расчетно-теоретические методы анализа полей течения и тепловых нагрузок, адекватность которых реальным условиям может быть установлена при всестороннем анализе возможных физических моделей и расчетных схем и сопоставлении их с данными натурных испытаний. Наиболее полная информация по характеристикам поля течения и тепловым нагрузкам получена в рамках программы полетов КА "Fire II" при входе в атмосферу Земли со скоростью,
превышающей вторую космическую. В ходе эксперимента проведены также измерения радиационного и суммарного тепловых потоков к передней и задней поверхностям.
В расчетном анализе условия невозмущенного набегающего потока соответствовали последнему участку траектории входа КА "Fire II" в атмосферу [1]. Вычисление тепловых потоков сильно усложнялось тем, что на заключительном этапе спуска течение около КА может переходить в турбулентное, наиболее интенсивное в области следа. Причиной значительного повышения теплового потока к поверхности аппарата может стать изменение режима течения. Отметим, что вычисление тепловых потоков для турбулентного режима течения более сложная задача, чем для ламинарного. Основная проблема заключается в выборе турбулентной модели, наиболее подходящей для исследуемых процессов обтекания. Обычно, при проектировании тепловой защиты неточность, связанная с неопределенностью в выборе турбулентной модели, компенсируется большим коэффициентом запаса. Одна из целей настоящей работы — это анализ различных турбулентных моделей для потока совершенного газа при небольших числах Рейнольдса и больших числах Маха, порядка M = 16, что соответствует условиям спуска. Для замыкания исходной системы уравнений применялись двухпараметрические к—ш [2, 3], BSL (baseline) к—ш [4, 5] и SST (shear stress transport) k—ш [5, 6] модели турбулентности, а также SSG (Sarkar, Speziale and Gatski) и LRR (Launder, Reece and Rodi) модели рейнольдсовых напряжений [7, 8] и проводились расчеты с использованием ламинарной модели. Данные летных испытаний [9-12], с которыми сравнивались результаты моделирования, включают в себя как термохимически неравновесные, так и равновесные режимы течения.
Наибольшим тепловым нагрузкам подвержена лобовая часть спускаемого аппарата; проведенные расчеты показывают, что для КА "Fire II" конвективный тепловой поток, приходящийся на заднюю поверхность, составляет всего 2... 3 % максимального теплового потока в окрестности передней критической точки. Несмотря на это, численное моделирование теплового состояния задней поверхности КА также необходимо при проектировании тепловой защиты.
Моделирование теплового состояния космического аппарата "Fire II" проводилось в ряде предыдущих исследований методами вычислительной аэродинамики [13, 14]. Авторы работы [15], используя расчетный код, основанный на уравнениях Навье-Стокса и учитывающий излучение, получили тепловые характеристики КА "Fire II" на участке траектории между 1634-й и 1651-й секундами спуска, что соответствует высотам от 77 до 37 км. В работе [16] проведено сравнение результатов расчетов по теплообмену для нескольких моделей турбулентности на задней и передней поверхностях летательного аппарата,
схожего по конфигурации с КА "Mars Pathfinder". Исследование КА такой же конфигурации проведено в работе [17], но с учетом перехода к турбулентному режиму течения и последующим использованием турбулентной модели.
В ряде работ изучаются переходное и турбулентное течения около КА "Fire II". В частности, в работе [18] был использован метод усреднения уравнений Навье-Стокса по Рейнольдсу (RANS, Reynolds-averaged Navier-Stokes) в совокупности с k—ш и Spalart-Allmaras турбулентными моделями для моделирования поля течения вокруг КА "FireII". В последующих работах [19-21] исследовались возможности и ограничения метода моделирования отсоединенных вихрей (DES, detached eddy simulation) при расчетах теплового состояния спускаемых КА. В работе [22] приведены результаты RANS моделирования по тепловым нагрузкам с учетом неравновесной химической кинетики, в результате которых получено хорошее соответствие с данными летных испытаний.
Угол входа КА "Fire II" в атмосферу составлял ~3°, скорость входа — 11327 м/с и превышала вторую космическую. В процессе движения в атмосфере угол атаки а возрастал, принимая значения 7,7°, 13° и достигая максимума 19,5° на заключительном этапе эксперимента. Увеличение угла атаки является причиной изменения физических характеристик КА "Fire II" в процессе спуска. Поэтому необходимы моделирование полей течения около КА и расчет местных характеристик при угле атаки, отличном от нуля. Результаты расчетов полей течения при а = 0 показывают, что течение в следе за аппаратом крайне асимметрично.
В большинстве предшествующих работ по вычислительной аэродинамике спускаемых КА основное внимание уделялось исследованию тепловых и динамических нагрузок на заднюю поверхность аппарата и аэротермодинамики следа за ним; кроме того, вычисления проводились в основном для а = 0. В настоящей статье приведены описания геометрической модели спускаемого аппарата и условий летного эксперимента, применяемых вычислительных моделей и методик построения расчетных сеток, представлены результаты трехмерного RANS моделирования поля течения для КА "Fire II" при различных углах атаки с использованием разных моделей турбулентности и проанализировано их влияние на структуру сложных трехмерных приповерхностных течений, местные тепловые и аэродинамические характеристики как задней, так и передней поверхностей КА "Fire II".
Геометрическая модель КА и условия набегающего потока. Цель программы полетов КА "Fire II" [10] состояла в исследовании нагрева летательного аппарата, имеющего характерную форму для спускаемых КА "Аполлон", входившего в земную атмосферу на скорости,
Рис. 1. Спускаемый КА "Fire II" с трехслойной тепловой защитой передней поверхности и держателем:
1... 3 — бериллиевые тепловые щиты; 4 — передняя поверхность аппарата; 5 — радиометрические окна; 6 — фенольно-асбестовые прослойки; 7, 15 — теплоизоляция; 8 — задняя поверхность (конус с углом 66°); 9 — радиометрическое окно пропускания волн оптического диапазона; 10, 14 — телеметрические антенны; 11 — антенна С -диапазона; 12 — механизм отделения спускаемого модуля; 13 — держатель; 16 — фенольно-асбестовая защита
превышающей вторую космическую скорость. Изучение поля течения около задней поверхности этого КА проводилось в условиях реального спуска в атмосфере Земли, что предпочтительнее проведения экспериментов в аэродинамической трубе (например, из-за отсутствия влияния хвостовой державки).
Космический аппарат "Fire II" представляет собой полностью автоматизированный летательный аппарат, состоящий из двигательной установки и спускаемого модуля. Компоновочная схема спускаемого модуля [10] показана на рис. 1. Модуль состоит из затупленной передней 4 и конической задней 8 поверхностей. Передняя часть модуля имеет многослойную структуру и включает три бериллиевых тепловых экрана 1, 2, 3, чередующихся с тремя фенольно-асбестовыми теплозащитными прослойками 6. Первые два слоя тепловой защиты отстреливаются в процессе спуска. Третий тепловой щит не отстреливается и определяет форму летательного аппарата на последнем участке траектории спуска, который исследуется в настоящей работе. Задняя поверхность представляет собой усеченный конус с углом 66°.
Задняя поверхность аппарата, выполненная в виде конуса (рис. 2, а), оснащена двенадцатью калориметрами, расположенными по окружности вдоль трех радиальных направлений, соответствующих углам
а
VR 0.609
Рис. 2. Приборное оснащение спускаемого КА "Fire II" (а) и внешний контур геометрической модели модуля, используемой для вычислений (б) (все размеры даны в сантиметрах)
р = 0; 120 и 240°. Значения температуры, зафиксированные датчиками на высоте 35 км, обрабатываются для получения теплового потока на поверхности и используются для сравнения с результатами, полученными методами вычислительной аэродинамики. Изменение угла атаки КА "Fire II" во время спуска определяли с помощью гироскопа и акселерометра, находящихся на борту спускаемого аппарата.
В работах [23-25] представлены результаты расчета изменения угла атаки во времени и основные показатели перемещения модуля, а также проанализированы колебания, возникающие при движении. Для определения углов атаки использовались данные, полученные с помощью гироскопа и акселерометра, а также результаты статических и аэродинамических измерений, выполненных в аэродинамической трубе. Углы атаки для некоторых периодов определялись методами компьютерного моделирования. Аэродинамические силы и моменты вы-
числялись на основе данных, полученных с помощью акселерометра и гироскопа.
Линия, очерчивающая внешний контур КА "Fire II", моделируемый в настоящей работе, приведена на рис. 2, б. Геометрическая модель модуля (см. рис. 2, б) имеет некоторые упрощения, аналогичные упрощениям, принятым в работе [1]; в частности, антенна 11 (см. рис. 1) в расчетах заменена плоской поверхностью.
В настоящей работе моделировались условия, задающие невозмущенный набегающий поток в точке конечного участка реальной траектории входа КА в атмосферу, которая соответствует одной из наименьших высот (35 км) и для которой были доступны измерения.
Условия набегающего потока, используемые при моделировании течения
около КА "Fire II"
Высота, км..........................................................................................................35
Число Маха........................................................................................................16
Температура набегающего потока, К........................................................237
Плотность набегающего потока, кг/м3 ......................................................0,0082
Скорость потока, км/с......................................................................................5
Температура поверхности модуля, К........................................................553,3
Угол атаки, град..................................................................................................0; 5; 10; 15; 20
Число Рейнольдса (х106)..............................................................................1,76
Вычислительные модели. Исходная система уравнений. Поле течения вокруг спускаемого КА моделировали, используя уравнения Навье-Стокса в трехмерной постановке, уравнение неразрывности и уравнение сохранения энергии.
Уравнения Навье-Стокса и уравнение неразрывности:
2 д i ди\ д дх \ дх) ду
dv ди дх ду
д_
дг
ды ди дх дг
; (1)
£ + div (pvV ) = -1 - 31 0" divV) +
^ д ( д^ д ду \ ду) дх
v и ^ дх ду J
д_
дг
ды дv ду дг
(2)
d.v (pwv )=- др - з I (, divv)+
2 д / ды\ д
дг у дг I дх
i ды ди
х г
д_ ду
ды дгиN ду дг/
(3)
дР + а1у (ру) = о. (4)
Уравнение сохранения энергии
дТ др
Рср-д£ + рс'рУ §гай Т = ^ (Л §гай Т) + + У gradр + /Ф, (5)
где диссипативная функция определяется как
Ф 2 I ди\2 2 (ду\2 2 (ди,\2 ( дь ди\2 дх) \ду ) \дг ) \дх ду)
дь дь\2 ( ди дь\2 2 (ди дь дь
\ду + дг) + \дг + дх) 3 \дх + ду + д^ ^
Для замыкания уравнений Навье-Стокса применялись двухпара-метрические модели турбулентности и модели рейнольдсовых напряжений.
В двухпараметрических моделях решались два раздельных уравнения переноса, определяющих турбулентную скорость и пространственный масштаб турбулентного движения (длину пути смешения). В таких моделях используется градиентная гипотеза Буссинеска, связывающая рейнольдсовые напряжения с усредненными градиентами скоростей и турбулентной вязкостью. В рамках настоящей гипотезы вводится предположение об изотропности турбулентной вязкости.
Турбулентная вязкость, в свою очередь, принимается пропорциональной величине пути перемешивания и масштабу турбулентной скорости. В рассматриваемых двухпараметрических моделях масштаб турбулентной скорости вычисляется через турбулентную кинетическую энергию к, которая находится из решения соответствующего уравнения переноса. Пространственный масштаб турбулентности определяется с использованием двух параметров потока, обычно это турбулентная кинетическая энергия к и скорость турбулентной диссипации е. Скорость диссипации турбулентной кинетической энергии определяется из уравнения переноса для е.
к—ш модель [2, 3]. Полуэмпирическая модель турбулентности выполняется на основе уравнений переноса турбулентной кинетической энергии к и специфической скорости диссипации (частоты турбулентных пульсаций) ш.
Одно из преимуществ к—ш модели — это применимость данной формулировки к моделированию пристеночных течений при больших продольных перепадах давления для низкорейнольдсовых вычислений. Модель не включает в себя сложную нелинейную демпфирующую функцию, учитывающую влияние вязкости на характеристики турбулентности в окрестности стенки, необходимую для к—е модели,
поэтому к—ш модель имеет менее высокую точность решения, но при этом метод не такой трудоемкий.
Модель предполагает, что турбулентная вязкость связана с турбулентной кинетической энергией и частотой турбулентных пульсаций (псевдозавихренностью) соотношением
к (7) / = р-. (7)
ш
Модель предполагает решение двух уравнений переноса: одно — для турбулентной кинетической энергии, другое — для частоты турбулентных пульсаций. Тензор турбулентных напряжений вычисляется из концепции изотропной турбулентной вязкости.
Уравнение для кинетической энергии турбулентности:
д (рк)
у • (pUk) = V
dt
Уравнение для частоты пульсаций:
р + р J Vk
Ok
Pk - ß'ркш. (8)
^ + V • (pUw) = V •
р + р ) Vw о.
<x-Pk - ßpw2. (9)
В этих уравнениях Pk представляет собой производство турбулентной кинетической энергии:
Pk = ptVU • (VU + VUT) - 2V • U(3p*V • U + pk); (10) константы модели имеют следующие значения:
в = 0,09; а = 5/9; в = 0,075; ак = 2; = 2. (11)
BSL (Baseline) k—ш модель [4, 5]. Главный недостаток k—ш модели — это сильная чувствительность к условиям набегающего потока. Вследствие зависимости от параметра частоты турбулентных пульсаций ш, принимаемого на входе в расчетную область, может быть получен значительный разброс в результатах моделирования. Этот нежелательный эффект был устранен благодаря применению гибридного метода. Метод, заключающийся в использовании k-ш модели в пристеночных областях и k—e модели [26] в областях, находящихся на достаточном удалении от стенки, был применен Ментером [4, 5]. Он заключается в некотором преобразовании k—e модели в k—ш формулировку и в последующем добавлении соответствующих уравнений. Таким образом, k—e модель дополняется стыковочной функцией 1 — F1. Данная функция принимает значение Fi = 1 вблизи поверхности и обращается в ноль за пределами пограничного слоя, т.е. на линии границы пограничного слоя и за его пределами k—e модель возвращается к первоначальной стандартной формулировке.
k-ш модель:
д (pk)
öt д (рш) т
V • (pUk) = V
i
Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
V • (риш) = V
р + ^ Vk Okl J
Pk - ß'ркш;
(12)
р + ) Vw
ш
ai kPk - ßipw2. (13)
Преобразованная к-е модель:
д (рк)
т
V • (pUk) = V
р + ^ Vk
Ok,2)
Pk - ß'pku;
(14)
д (рш) m
V • (риш) = V •
р + ) Vw Ош2
2р—1— УкУш + а2тРк - в2рш2. (15) Ош2 ш к
Теперь, если уравнения видоизмененной к—е модели дополнить функцией 1 — и добавить соответствующие уравнения к—ш модели, то получим ББЬ модель турбулентности:
д (рк)
т
д (рш)
т
V • (pUk) = V
V • (риш) = V •
р + Vk
Ok3 J
Pk - ß'pkw;
(16)
р + — )Уш +
(1 - ^)2р—1— УкУш + аз ШРк - взрш2. (17) Ош2 ш к
Коэффициенты новой модели — это линейная комбинация соответствующих коэффициентов моделей, лежащих в основе метода:
Фз = ЗД + (1 - *)Ф2. (18)
Следовательно, константы модели определяются следующим образом:
в' = 0,09; а1 = 5/9; в1 = 0,075; ак1 = 2; аш1 = 2;
a2 = 0,44; ß2 = 0,0828; Ok2 = 1; = 1/0,856.
(19)
Полученная модель сочетает в себе преимущества k—e (нечувствительность к параметрам свободного потока) и k—ш моделей (дает хорошие результаты в пограничном слое в непосредственной близости к стенке, особенно при наличии большого неблагоприятного градиента давления).
Shear Stress Transport (SST) k—ш модель [5, 6]. При применении гибридных моделей часть проблем удается избежать, однако, характеристики таких моделей существенно ухудшаются при наличии отрывных
течений. При использовании гибридной BSL модели турбулентности плохо предсказывается отрыв потока под действием неблагоприятного градиента давления. Модель дает завышенные значения турбулентной вязкости: k—£ и k—ш модели не учитывают перенос касательных напряжений.
Для устранения перечисленных недостатков была разработана модификация модели Ментера, так называемая SST (shear stress transport — перенос касательных напряжений) модель. SST модель переноса касательных напряжений является столь же экономичной, как и k-£ модель, но обеспечивает более высокое качество воспроизведения процессов для широкого ряда потоков и граничных условий. Она дает очень точное предсказание отрыва потока при неблагоприятных градиентах давления.
В этой модели точное поведение параметров переноса получено с помощью ограничения в формулировке турбулентной вязкости
Vt =--,-WWT, (20)
max(a1w, sF2)
где
Vt = pt/p; (21)
F2 — стыковочная функция, аналогичная функции F1 предыдущей модели.
Основу модели составляют те же два уравнения переноса, что и в случае BSL модели: уравнение переноса для частоты турбулентных пульсаций видоизмененной k-£ модели, дополненное стыковочной функцией 1—F1 (17), и уравнение переноса турбулентной кинетической энергии модели Уилкокса (16). Коэффициенты SST модели так же, как и в предыдущем случае, получаются как линейная комбинация коэффициентов k—ш и k—£ моделей (18), (19).
Существенное различие заключается в определении стыковочных функций (причина успеха метода), которые в данной модели представляют собой сложные функции безразмерной координаты стенки у и параметров турбулентного потока:
F1 = tanh(argl); (22)
( ( Vk 500v\ 4pk \
argi = min max ——, , —-- ; (23)
y ye шу у2ш J CDkwаШ2y2 J
CDkw = max ( 2p—1— VkVu, 1 • 10-10 ) ; (24)
V J
F2 = tanh(arg2); (25)
( 2Vk 500V .
arg2 = max —-, . (26)
\ в'шу у2ш
Модели рейнольдсовых напряжений. Двухпараметрические модели турбулентности, описанные ранее, предполагают справедливость градиентной гипотезы Буссинеска, линейно связывающей рейнольдсо-вые напряжения с усредненными градиентами скоростей и турбулентной вязкостью. В этом случае главные направления тензора скоростей деформации соответствуют главным направлениям тензора напряжений. Это разумно для весьма простого деформированного состояния, особенно, когда константы модели тщательно откалиброваны для рассматриваемых классов течений, но может быть не применимо при моделировании сложных полей деформаций, появляющихся при наличии закрутки потока или сложной геометрии. При таких условиях необходимо использовать более сложную связь между напряжением и скоростью деформации.
Модели переноса рейнольдсовых напряжений (Reynolds Stress Models) определяют турбулентные напряжения непосредственно по уравнениям переноса для каждого компонента тензора напряжений, что требует решения дополнительных шести взаимозависимых уравнений совместно с уравнением скорости турбулентной диссипации e, по которому определяется пространственный масштаб турбулентности.
Дифференциальное уравнение переноса рейнольдсовых напряжений имеет вид
дри & и ät
V • (pU & и & и) =
2k
= P + ф + V • р + -Csp— Vu & и - -8ре, (27)
3
2
3
где Р = -р(и ® и(Уи)т + (Уи)и ® и).
Уравнение переноса рейнольдсовых напряжений в тензорной записи можно представить как
дтц + = Pij + фгз + д
m
дхк
дхк
дт
V-
дхк
сг^ ) - е
-ij,
(28)
где Tij = pUiUj — тензор рейнольдсовых напряжений; eij = -5ijpe —
определение диссипации с использованием гипотезы Колмогорова о
2 к2 (дЩщ\
локальной изотропии; cij = -csp— —- — турбулентный перенос
3 e V dxk J
[Launder, Reece and Rodi].
Итак, получаем систему уравнений для модели рейнольдсовых напряжений:
Jt (pUUj) + дх (Uk рии) =
k д = Pij + фц +
дхъ
2 k2 \ ди и
р + ö Csp—
3 е / дх
H^j
2 Sijре (29)
2
k
i
Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
— дифференциальное уравнение переноса для тензора рейнольдсовых напряжений;
д(рс) д с д
—(ßUkс) = T(c£lP - ое2ре) +
, Pt р +--
О г
дс dxk
(30)
дЬ Охи к Охи
— уравнение для скорости диссипации турбулентной энергии;
фЦ = Фи! + ф'Ц2 (31)
— корреляционное соотношение для давления и деформации, где
Фт = —рс I Csla + C,
-',s2
aa — i a • aö t t ,
(32)
Фт2 = —Cr\Pa + Cr2pkS — Cr3pkS^a • a+
2
—с 3
2 3(
Cripk(aST + SaT — 2a • S^j + Crbpk (aWT + WaT) ; (33)
U ^ U 1 /п77 / Г,,ПТ\
a = —;--- о — тензор анизотропии; S = - ( VU + (VU) —
k
скорость деформации; W = 1 ^VU — (VU— завихренность.
Выше приведена общая форма записи для корреляции. В зависимости от значений соответствующих коэффициентов корреляция может быть линейной и квадратичной. По значениям коэффициентов корреляции различают три стандартных модели рейнольдсовых напряжений.
В моделях LRR-IP и LRR-QI ("IP" — Isotropisation of Production, "QI" — Quasi-Isotropic) [8] корреляция линейная (табл. 1).
SSG модель была разработана в работе [7]; при моделировании использовалась квадратичная корреляция (см. табл. 1).
Таблица 1
Значения констант для LRR и SSG моделей турбулентности
Модель Cs C£1 Cs2 Cs1 Cs2 Cr1 Cr2 Cr3 Cr4 Cr5
LRR-IP 0,22 1,45 1,9 1,8 0 0 0,8 0 0,6 0,6
LRR-QI 0,22 1,45 1,9 1,8 0 0 0,8 0 0,873 0,655
SSG 0,22 1,45 1,83 1,7 -1,05 0,9 0,8 0,65 0,625 0,2
Однако, как показывает практика, учет анизотропии касательных напряжений, т.е. использование полного тензора в моделях высокого порядка, для многих случаев не дает никаких преимуществ по сравнению с двухпараметрическими моделями, хотя требует значительно больших вычислительных ресурсов. К тому же модели рейнольдсовых
напряжений имеют худшую сходимость из-за сильной связи уравнений.
Далее в табл. 2 приведен краткий обзор используемых моделей с точки зрения вышеперечисленных основных преимуществ и недостатков.
Таблица 2
Краткий обзор примененных моделей турбулентности
Модель Преимущества Недостатки
к—ш 1. Применимость данной формулировки к моделированию пристеночных течений при больших продольных перепадах давления 2. Модель не включает в себя сложную нелинейную демпфирующую функцию, учитывающую влияние вязкости на характеристики турбулентности в окрестности стенки, необходимую для к—e модели, следовательно, метод не такой трудоемкий 1. Сильная чувствительность к условиям набегающего потока (в зависимости от параметра ш, принимаемого на входе в расчетную область, может быть получен значительный разброс в результатах моделирования) 2. Низкая точность решения в области свободного течения
BSL 1. Нечувствительность к параметрам свободного потока (в частности, к ш) 2. Сочетает преимущества к—ш и к— e моделей применительно к пристеночным течениям и области свободного течения соответственно 1. Модель плохо предсказывает отрыв потока под действием неблагоприятного градиента давления. 2. Модель дает завышенные значения турбулентной вязкости
SST 1. Обладает всеми преимуществами BSL модели 2. Модель учитывает перенос касательных напряжений 3. Дает очень точное предсказание отрыва потока при неблагоприятных градиентах давления Не изучены
LRR 1. В рамках моделей решается дифференциальное уравнение переноса для каждого компонента тензора рейнольдсовых напряжений 2. Данные модели дают самые лучшие результаты в случае сложных течений, при наличии закрутки потока или сложной геометрии 1. Как показывает практика, учет анизотропии касательных напряжений в моделях высокого порядка часто не дает никаких преимуществ по сравнению с двухпараметрическими моделями 2. Требуется значительно больше вычислительных ресурсов 3. Модели рейнольдсовых напряжений имеют худшую сходимость из-за сильной связи уравнений
SSG
Построение расчетных сеток. При решении задач вычислительной гидродинамики процесс построения сетки обычно наиболее трудоемкий, требующий больших временных затрат. Точность и эффек-
тивность численного исследования существенно зависят от используемой расчетной сетки.
В настоящей работе моделирование трехмерного поля течения около КА "Fire II" проводилось с использованием неструктурированных тетраэдральных сеток.
На рис. 3 приведены различные фрагменты неструктурированной расчетной сетки, применявшейся для моделирования внешнего обтекания модели спускаемого КА "Fire II". Из рисунка хорошо видна степень сгущения узлов сетки к обтекаемой поверхности. Размер ячеек, непосредственно прилегающих к поверхностям КА "Fire II" составляет 10-4 м, что является достаточным для определения местных аэродинамических и тепловых характеристик передней и задней поверхностей КА. Размерность сетки 1188329 тетраэдральных элементов.
Неструктурированные сетки имеют бесспорное преимущество в том, что:
1)для большинства задач, представляющих практический интерес, форма расчетной области очень сложна и построение структурированной расчетной сетки в этом случае становится чрезвычайно трудоемким даже при использовании многооболочного подхода. Таким образом, неструктурированные сетки могут быть приспособлены к любой
Рис. 3. Фрагменты трехмерной расчетной сетки (справа показана структура сетки вблизи обтекаемого модуля и на его передней и задней поверхностях)
желаемой геометрии (сетка имеет большую гибкость применительно к сложным геометриям);
процесс построения неструктурированных сеток может быть автоматизирован в высокой степени, что может существенно сократить временные затраты на генерацию сеток. Построение качественной неструктурированной сетки занимает максимум несколько часов. Для создания структурированной сетки иногда требуется несколько недель;
неструктурированные сетки проще генерировать.
Неструктурированный метод имеет и ряд недостатков:
точность данного метода невысока, так как тетраэдр не имеет параллельных граней. Для получения требуемой точности результатов необходимо гораздо больше расчетных элементов, чем в случае применения структурированной сетки;
сетка должна быть в высокой степени изотропна. Соседние элементы могут лишь очень незначительно отличаться по форме и размерам. Это также приводит к существенному увеличению числа необходимых расчетных элементов;
неструктурированная сетка зависит от качества самой геометрической модели; даже незначительная ошибка в задании геометрии может сделать построение таких сеток невозможным;
некоторые трудности возникают при интерполяции расчетных данных на неструктурированных сетках для наглядного представления результатов расчетов в требуемых плоскостях и на поверхностях;
решатель, использующий неструктурированные сетки, требует больше памяти и расчетного времени. Численные алгоритмы, основанные на неструктурированной топологии сетки, более трудоемкие и сложные. Это связано с увеличением числа ячеек и граней по сравнению с шестигранными сетками. Тетраэдральная сетка из N узлов имеет около 6N ячеек и 12N граней, в то время как шестигранная состоит из порядка N ячеек и 3N граней;
для записи и хранения неструктурированных сеток требуется больше информации (дополнительная память требуется для хранения информации о связях ячеек сетки).
Производительность современных компьютеров и развитие параллельных вычислений позволяют решать все перечисленные проблемы, возникающие при использовании неструктурированного метода.
Для генерации сетки использовался метод распространяющегося фронта (advancing front method).
Алгоритм включает следующие основные этапы:
построение двумерной триангуляции тем же методом на границах расчетной области и поверхности КА;
сетка на границах области определяется как начальный фронт;
для каждой элементарной площадки фронта определяется оптимальное положение нового узла (из условия непересечения и с учетом расположения всех уже существующих соседних узлов);
новый узел и элементарная площадка фронта определяют новый тетраэдрический элемент. По мере образования ячеек сетки новые элементарные площадки добавляются к фронту ("использованные" для построения тетраэдров площадки фронту уже не принадлежат). Таким образом, активный фронт распространяется внутрь расчетной области. Тетраэдры строятся от граней, образующих фронт, и наращиваются внутрь расчетной области.
Результаты расчетов. На рис. 4, 5 (см. 4-ю полосы обложки) приведены результаты численного моделирования в виде полей чисел Маха при различных углах атаки. В каждом случае данные представлены в центральном осевом сечении, в плоскости наклона. Наиболее яркие особенности структуры течения отмечены на рис. 4 (результаты соответствуют нулевому углу атаки а).
Структура течения включает в себя головную ударную волну, отстоящую от передней части летательного аппарата; обширную зону формирования возвратных течений, прилегающую к задней поверхности КА "Fire II" сразу за точкой отрыва пограничного слоя; зону рециркуляции за КА. Внешний невязкий поток над рециркуляционной зоной разворачивается в направлении от оси, образуя ударную волну повторного сжатия, отходящую от горловины потока. Поток, набегающий на лобовую поверхность спускаемого КА, испытывает сильное сжатие и последующее расширение при стекании с кромки переднего теплового щита, образуя область расширения потока.
Сама форма головной ударной волны существенно изменяется при вариации угла атаки (см. рис. 5, 4-я полоса обложки). Ударная волна симметрична относительно оси при нулевом угле атаки. При возрастании угла атаки головная ударная волна все больше прилегает к наветренной стороне (результат сильного сжатия потока) и отходит от подветренной стороны передней части КА. Застойная область вблизи носовой части капсулы также постепенно смещается вниз с ростом угла атаки.
Главный эффект от вариации угла атаки выражается в изменении формы и размера рециркуляционной области за капсулой и изменении структуры течения (см. рис. 5, 4-я полоса обложки). При а = 0 отрыв пограничного слоя происходит сразу после затекания потока на усеченный конус задней поверхности КА, при этом образуется осе-симметричная рециркуляционная зона, прилегающая к задней поверхности. Если угол набегания потока не равен нулю, то зона возвратных течений за КА "Fire II" наклонена и меньше в размерах. При а = 0° набегающий поток плавно обтекает наветренную сторону усеченного
конуса задней поверхности без образования отрывных течений, хотя на подветренной стороне отрывная зона сохраняется и не изменяет расположения. При а = 5° след за обтекаемым телом меньше, чем в случае нулевого угла атаки, главным образом из-за позднего отрыва потока на наветренной поверхности. При дальнейшем возрастании угла атаки эффект, связанный с более поздним отрывом потока, более чем компенсируется благодаря возрастанию наклона пограничного слоя (по отношению к усеченному конусу аппарата). Поэтому при значительных углах падения (а > 10°) ширина следа возрастает по сравнению с нулевым углом атаки. Похожие структуры течений для различных углов атаки были получены в работе [1]. Отличие заключается в том, что по результатам настоящей работы отрыв течения на подветренной поверхности конуса КА "Fire II" происходит несколько позже. Также некоторые расхождения наблюдаются в размерах рециркуляционной зоны.
Значительные изменения в структуре потока при вариации угла атаки также можно проследить по рис. 6 (см. 3-ю полосу обложки). Здесь для нескольких углов атаки структура течения представлена в виде линий тока, нанесенных на поле чисел Маха.
Как видно из рис. 6, при нулевом угле атаки ближний след состоит из осесимметричного тороидального вихря, а также видны две его отдельные петли в центральном осевом сечении. При а = 5° петля вихря наблюдается только на подветренной стороне и ее размер постепенно уменьшается при возрастании угла атаки. Вихревая структура формируется вплоть до значений углов атаки порядка а = 15°. Центр вихревого образования с ростом угла атаки перемещается вверх по потоку вдоль подветренной конусной поверхности.
На рис. 7 с помощью линий тока (линий тензора касательных напряжений) показана структура трехмерных вихревых течений у задней поверхности КА "Fire II" при разных углах атаки.
Расходящиеся по поверхности (не пересекающиеся линии тока) соответствуют безотрывному обтеканию задней поверхности КА. Схождение линий в одной точке соответствует отрыву потока.
При а = 0 вихревая структура осесимметрична и покрывает всю заднюю поверхность летательного аппарата и плоское основание усеченного конуса. Даже при маленьком угле атаки (а = 5°) размеры вихревой зоны значительно уменьшаются, при этом поток имеет безотрывной характер течения уже только на наветренной стороне. Симметрия относительно центрального осевого сечения (плоскости наклона) сохраняется при любых углах атаки.
Местные тепловые характеристики. Далее исследуется влияние различных моделей турбулентности на местные тепловые характеристики (плотность конвективных тепловых потоков) вдоль передней (носовой) и задней поверхностей (рис. 8, 9).
Иллюстрации к статье А.Л. Железняковой, С.Т. Суржикова «Численное моделирование поля течения при входе в атмосферу Земли спускаемого аппарата с аэродинамическим качеством»
Рис. 4. Структура поля течения при нулевом угле атаки
Рис. 5. Распределение числа Маха в поле течения КА «Fire II» в плоскости наклона при разных углах атаки
Иллюстрации к статье А.Л. Железняковой, С.Т. Суржикова «Численное моделирование поля течения при входе в атмосферу Земли спускаемого аппарата с аэродинамическим качеством»
Рис. 6. Структура течений около КА «Fire II» в плоскости наклона при разных углах атаки (линии тока нанесены на поле чисел Маха)
V
J
Рис. 7. Линии тензора касательных напряжений на задней поверхности модуля при разных углах атаки
Рис. 8. Распределения плотности теплового потока вдоль всей поверхности КА "Fire II" от передней критической точки для различных моделей турбулентности
Как и ожидалось, при нулевом угле атаки распределение местных характеристик осесимметрично. Плотность теплового потока достигает максимальных значений в области передней критической точки.
Плотность теплового потока на передней и задней поверхностях летательного аппарата изменяется от 5 Вт/см2 в зоне отрыва потока (соответствует s/D = 0,9, s — расстояние вдоль поверхности аппарата от передней критической точки; D — максимальный диаметр модуля) до 250Вт/см2 (максимальное значение зависит от используемой модели турбулентности) в области передней критической точки (см. рис. 8).
35
30 25
а 20
о CQ
15
10
0
7 я - а = 0 8' - avg а • ~(р = 0 * ~(р = 120( а - <р -240° 1 - SST 2-к-о) 3 - BSL 4- SSG 5 - LRR 6- Lam 3 щ Ш щ
1 U 7. * Я
ч ■ ш 1 8.
i
Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
/VI1 ? * ■ J д * / * ■ г revisits üVirfu»1 J
■ ■ ■ ■ 1 в ■ ■ glSilМШШШ9* 5 4
0,6 0,8 1,0 1,2 1,4 5/ В
Рис. 9. Распределения плотности теплового потока на задней поверхности КА для различных моделей турбулентности. На графическую область также нанесены данные летного эксперимента для калориметров, расположенных под углами ^ = 0; 120 и 240°, и результаты численного моделирования [1]
Из рис. 8, 9 видно, что значения местных тепловых характеристик на задней поверхности спускаемого КА, полученные с помощью различных моделей турбулентности, отличаются незначительно, в то время как распределения плотности тепловых потоков на передней поверхности очень различны.
Как отмечалось ранее, задняя поверхность спускаемого КА была оснащена двенадцатью калориметрами, расположенными по окружности вдоль трех радиальных направлений, соответствующих углам = 0; 120 и 240° (см. рис.2). Данные летного эксперимента [10], зафиксированные этими измерительными приборами, также помещены в графическую область распределения плотности теплового потока в зависимости от безразмерного параметра s/D (см. рис.9). Значения s/D = 0,64 и 1,43 соответствуют началу и концу задней конической поверхности КА "Fire II" соответственно.
Обратим внимание на значения плотности теплового потока, взятые из летного эксперимента [10] (см. рис.9). Показания датчиков получены в точке траектории спускаемого модуля на высоте 35 км. Угол атаки в данной точке траектории достигал максимального значения порядка 15°. Исходя из структуры течения при ненулевом угле атаки (отрывной характер течения на подветренной стороне и плавное обтекание на наветренной, см. рис. 5,6, 4-я и 3-я полосы обложки) можно было предположить, что показания датчиков на наветренной и подветренной сторонах капсулы должны значительно отличаться в зависимости от угла расположения. Например, датчик, расположенный в наветренной зоне плотного прилегания потока, должен регистрировать большие тепловые нагрузки, чем датчик, расположенный в зоне отрывного течения. Однако данные, полученные с помощью летного эксперимента (см. рис. 9), нанесенные на общий график черными треугольниками, кружками и ромбами, не показывают значительного отличия в зависимости от расположения по углу. Это может быть вызвано тем, что при входе в плотные слои атмосферы спускаемый модуль вращается с угловой скоростью порядка 160 об/мин. Возможно, с помощью датчиков были получены значения тепловых нагрузок, усредненные по углу. В связи с этим для сравнения численных результатов с экспериментальными данными [10] при ненулевом угле атаки были получены усредненные по углу значения плотности теплового потока для каждого рассматриваемого угла атаки. Если учесть, что в процессе спуска, кроме вращения, наблюдаются значительные колебания угла атаки, то для сравнения c экспериментальными данными [10] распределения плотности тепловых потоков были также усреднены по углам атаки. Получившиеся распределения для различных моделей турбулентности приведены на рис. 8, 9. Также на общую диаграмму
(см. рис. 9) нанесены усредненное по углам атаки от 0 до 20° распределения плотности теплового потока и распределение теплового потока для нулевого угла атаки [1].
При нулевом угле атаки в зоне отрыва потока (s/D = 0,9) наблюдаются минимальные значения плотности теплового потока. Результаты расчетов, приведенные в работе [1], указывают на минимальное значение тепловых потоков при s/D = 0,7 и, как следствие, на более ранний отрыв потока на подветренной стороне конуса КА "Fire II". При движении вдоль задней поверхности тепловой поток возрастает незначительно. Распределение, полученное в работе [1], при нулевом угле атаки дает более значительное возрастание плотности теплового потока при движении вдоль конуса задней поверхности КА "Fire II". Усреднение значений плотности теплового потока по углам атаки приводит к снижению этой характеристики приблизительно на 40% для задней поверхности модуля.
Из рис. 9 видно, что полученные результаты хорошо согласуются с экспериментальными данными.
Заключение. В работе приведены результаты численного моделирования обтекания спускаемого КА "Fire II" с использованием уравнений Навье-Стокса в трехмерной постановке при измении угла атаки от 0 до 20° и применении различных моделей турбулентности. Результаты расчетов показали, что структура течения значительно изменяется даже при маленьком угле атаки. Пограничный слой на наветренной стороне плотно прилегает к поверхности летательного аппарата. На подветренной стороне, напротив, при натекании на заднюю коническую поверхность КА "Fire II" происходит отрыв потока с образованием рециркуляционной зоны. Две противоположные тенденции наблюдаются для плотности теплового потока на задней конической поверхности капсулы. Плотность теплового потока возрастает в наветренной зоне плотного примыкания пограничного слоя с ростом угла атаки, в то время как в рециркуляционной зоне отрыва на подветренной стороне усеченного конуса с ростом угла атаки наблюдается снижение тепловых нагрузок. Это приводит к тому, что усредненные по углу (по всей поверхности аппарата) значения плотности теплового потока практически не изменяются при вариации угла атаки. Усреднение по углам атаки приводит к более равномерному распределению плотности теплового потока на задней поверхности КА "Fire II" и, следовательно, к несколько меньшим тепловым нагрузкам на задней поверхности КА по сравнению с результатами моделирования при а = 0.