УДК 533.6.011.8
М. А. Котов, В. В. Кузенов
ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ОБТЕКАНИЯ ПОВЕРХНОСТЕЙ ПЕРСПЕКТИВНЫХ ГИПЕРЗВУКОВЫХ ЛЕТАТЕЛЬНЫХ АППАРАТОВ
Кратко описан алгоритм численного решения уравнений Навье-Стокса, осредненных по Рейнольдсу в областях сложной пространственной формы. Приведены результаты расчета обтекания сверхзвуковым потоком воздуха гиперзвукового летательного аппарата Waverider. Изложены различные методы создания сложных моделей поверхностей перспективных гиперзвуковых летательных аппаратов X-43, Waverider, Х-51А.
E-mail: mikhail_kotov88@mail.ru; kuzenov@ipmnet.ru
Ключевые слова: газовая динамика, гиперзвуковые летательные аппараты, уравнения Навье-Стокса, численные методы, адаптивные декартовы сетки, геометрическое моделирование сложных поверхностей.
К одним из наиболее важных и сложных проблем аэротермодинамики аэрокосмических аппаратов можно отнести вопросы моделирования сложных (трехмерных) элементов течения у поверхности гиперзвуковых летательных аппаратов (ГЛА): отрывных течений и течений присоединения, ламинарно-турбулентных переходов, сложных ударно-волновых взаимодействий). Однако реальные физические эксперименты в этой области отличается высокой стоимостью и сопряжены со множеством технологических и технических трудностей. Поэтому большое значение для оптимизации конструкции и рабочих характеристик ГЛА имеет численное моделирование аэротермодинамических и теплофизических процессов. При рассмотрении сверх- и гиперзвуковых течений необходимо также, чтобы вычислительные коды, предназначенные для расчетов течений газа в рамках уравнений Эйлера, Навье-Стокса и Рейнольдса, включали в себя специализированные базы данных по термодинамическим, кинетическим и транспортным свойствам индивидуальных газов и газовых смесей.
На первом этапе из всей совокупности указанных задач интерес представляют лишь моделирование формы изучаемых объектов техники, их размеров и взаимное расположение, а не физических свойств этих тел. Построение таких поверхностей базируется на аналитической и дифференциальной геометрии, методах вычислительной математики. Усложняющей особенностью при построении геометрической модели ГЛА является необходимость учета многосвязности областей, составляющих поверхность тела, и сопряжений между элементами поверхности. Такой тщательный учет особенностей формы поверхности связан с решением последующей задачи построения расчетных сеток
высокого качества. Создание сложной модели поверхности таким способом представляет собой достаточно медленный, трудоемкий процесс и требует многократного редактирования модели поверхности.
На следующем этапе решения задачи оптимизации конструкции и рабочих характеристик ГЛА целесообразно создание вычислительных моделей, расчетных сеток, которые позволяют проводить расчетно-теоретические исследования аэротермодинамики перспективных ГЛА и космических аппаратов, предназначенных для входа в плотные слои атмосфер планет и возвращения на Землю [1, 2].
Постановка задачи. Цель настоящей работы — создание вычислительных аэротермодинамических моделей ГЛА, которые опираются на систему уравнений вязкого теплопроводного газа и могут быть применены для проведения численных экспериментов в широком диапазоне полетных условий. В них используются трехмерные уравнения Навье-Стокса, осредненные по Рейнольдсу, и они позволяют находить теплофизические характеристики потока в расчетной области сложной пространственной формы. Расчетная область включает в себя поле течения в невозмущенном потоке, за фронтом ударной волны, в следе за обтекаемым телом. Для решения рассматриваемой системы уравнений расчетная область представляется в виде регулярной декартовой (структурированной) сетки.
Одна из важных задач настоящей работы состоит также в геометрическом моделировании сложных поверхностей перспективных ГЛА типа Иурег-Х [4, 5], Waverider [6], Х-51А [7] в системе автоматизированного проектирования (САПР) SolidWorks [8].
Описание процесса геометрического моделирования тел сложной пространственной формы. Реальные ГЛА занимают некоторый конечный объем в пространстве, и для геометрического моделирования таких тел в математике используются объекты, называемые твердыми телами (или деталями). Основной метод моделирования геометрии тел использует форму простейших тел (прямоугольную призму, цилиндр, конус, шар, тор и др.), которые имеют одну оболочку. Далее путем удаления и/или добавления в определенных местах дополнительного объема (материала) получают тело требуемой формы. Существует другой вариант многомерного моделирования геометрии тел. Этот подход соответствует реальному технологическому процессу изготовления детали. В данном случае на передний план выходит не форма модели, а замысел проекта, т.е. последовательность выполняемых действий. Такой подход к моделированию использует САПР SolidWorks.
Система автоматизированного проектирования SolidWorks [8] предназначена для моделирования деталей и сборок в трехмерном пространстве с возможностью проведения различных видов экспресс-
Рис. 1. Общий вид геометрии поверхности модели экспериментального беспилотного ГЛА X-43
анализа, а также оформления конструкторской документации в соответствии с требованиями ЕСКД. Система использует графический интерфейс Microsoft Windows. В САПР SolidWorks имеются средства конвертации и сохранения в форматах, которые доступны для открытия и редактирования в других программах геометрического моделирования (Autodesk 3dsMax, AutoCAD и т.д.).
С использованием инструментов, входящих систему САПР SolidWorks, была построена геометрическая модель, приближенно соответствующая форме поверхности ГЛА X-43 [9, 10] (рис. 1). Эта форма поверхности ГЛА X-43 и поверхностная сетка на ней в дальнейшем была использована для проведения трехмерных термоаэрогазодинами-ческих расчетов течения сжимаемого газа для условий, при которых проводились полетные испытания ГЛА X-43 [11]. На рис.2 представлен вид окончательного варианта модели поверхности гиперзвукового аппарата типа Waverider [6]. На рис. 3 приведен вид предварительного варианта геометрической формы модели ГЛА Х-51А [12]. Отметим, что в дальнейшем планируется доработка формы модели ГЛА Х-51А и построение некоторых ее важных конструктивных элементов (ГПВРД, рули высоты и поворота), а также уточнение построенных элементов поверхности.
Рис.2. Изометрия поверхности модели экспериментального ГЛА Waverider
Рис. 3. Предварительный вариант поверхности модели экспериментального беспилотного ГЛА X-51А. Вид в изометрии
Математическая постановка рассматриваемой задачи сводится к системе трехмерных нестационарных осредненных по Рейнольдсу уравнений Навье-Стокса, которая выглядит следующим образом:
т + dXj(ри) = dt(ри) + dXj{рщи + ^) =0;
дд
(РЕ) + ^Г" (РЕиз + ujaä - Qj) = 0
(1)
д
дхА
где — время и декартовы координаты; р, щ и Е = Н/^ + и*/2 — осредненные среднемассовые плотность, компоненты скорости и полная энергия (Н = СрТ — энтальпия идеального газа (Ср — теплоемкость при постоянном давлении), 7 — отношение удельных теплоемкостей).
Тензор вязких напряжений , тепловой поток Оз и входящий в а^ тензор скоростей деформации определяются выражениями
aij = Sij Р
Sa —
2М Sij
1 д
3 дхк
ик S
1 _ t
ij I + °ij
Qj =
1
Re V Pr дх
Sij = 2
1 ( ди диa
+
дха дх
в которых тензор рейнольдсовых напряжений аТ и турбулентный поток теплоты О1- имеют следующий вид: аТ = 2^Т ( --^ ),
\ 3 дхк /
dt vt dh (7- 1) ,
Qt =--——. Здесь P = -ph — статическое давление; u и
j PrT dxj y
UT — коэффициенты молекулярной и турбулентной вязкостей; Pr и PrT — ламинарное и турбулентное числа Прандтля, Re — число Рей-нольдса. По повторяющимся индексам i,j,k везде в тексте предполагается суммирование. Система уравнений дополняется коэффициентом молекулярной вязкости и, кг/(м • с), вычисляемым по формуле Сазерленда
V = 1,458 • 10-6T1'5-1—-.
р ' 110,4 + T
Расчет турбулентных параметров течения может быть проведен на основе дифференциальной двухпараметрической q—ш модели Кокли [13]:
д д 1 д ^ (pq) + —. (pquj) —
dt
dx
Re dx,
д
(u + Prq ut )
д д 1 — (пш) + —— (ошп, ) — -1—7— dt ' dx j' Re dx
h = pq ( CfS 2 duk_
q 2 V ш 3 dxk
(u + Prw ut )
dq
dx jj
дш
=H
q
dx
= H ;
S — 2 ( Sij — 15ij d— — ш
ш ) ; H = p du
c'iCßS—2 шщг) — C-2
3 ij dxk
d-1 ; ut = ReC¡JJ'p(q-.
dx j ш
Пристеночная функция f имеет вид f = 1 — exp (—apqy •Re/^), где y — расстояние по нормали от твердой поверхности; параметр C1 = 0,045 + 0,405f, а остальные эмпирические константы принимают следующие значения: C2 = 0,92, Prq = 1, Ргш = 1,3, C^ = 0,09, а = 0,0065. В приведенных формулах параметр q отнесен к модулю скорости набегающего потока п\. Псевдозавихренность ш обезразме-ривается на величину uiL-1.
Для дальнейших построений систему уравнений (1) удобно представить в векторной форме в декартовой системе координат:
dQ + dF (Q) + dG (Q) + dH (Q) dt dx dy dz
= 0.
(2)
Вектор консервативных переменных (( и векторы потоков Г, СС, Н могут быть записаны в следующей форме:
/ Р \ ( ри \
Q =
pu
pV
pw pE
F =
puu + P — Œxx
puV — Oxy puW — O xz
\ (pE + P) u — uOxx — VOxy — WOxz + Qx j
G =
( pv
pvu - Cyx
pvv + P - Uyy pvw - GyZ
\ (рЕ + Р) и - поуХ - иоуу - та,уг + Яу )
/ рт \
рти - агх Н = ржи - агу
ртт + Р - о22 \ (рЕ + Р) т - иагх - иагу - тагг + Я )
Газовая среда, натекающая на ГЛА, рассматривается как совершенный газ с показателем адиабаты 7 = 1,4.
Для численного решения уравнений Навье-Стокса и дифференциальной двухпараметрической д-ш-модели вихревой вязкости Кокли использован явный интегроинтерполяционный метод, применение которого для уравнений (2) на равномерной декартовой сетке приводит к следующей конечно-разностной схеме:
Qn+l _Qn
Qi,j,k Qi,j,k
At
H n _ h n
+ Fi+l/2,j,k Fi-l/2,j,k +
+ Ax
(ТП _ (ТП тгп _ тгп
+ ^ъ,]+1/2,к ^ъ,]-1/2,к + Нъ,],к+1/2 Нъ,],к-1/2 = „ (з)
+ Ау + Аг = 0 (3)
где Ах, Ау и Аг — шаги сетки в направлениях осей х, у и г. Тройной целый нижний индекс (I, ], к) обозначает значения сеточных функций, отнесенные к центрам соответствующих трехмерных дискретных регулярных ячеек, а полуцелые нижние индексы относятся к значениям функций на границах ячеек. При численном решении гиперболической ("невязкой") части уравнения величины с полуцелыми индексами Яг±1/2)3)к, Яъ)]±1/2)к, Яг)3)к±1/2 находятся из решения соответствующих задач Римана [14] (для построения численной схемы следует решать задачу Римана на каждой из границ дискретных ячеек). Векторы потоков, которые входят в уравнение (2), определяются следующим образом:
Ръ±1/2,3,к = Р [Я1±1/2,3,к
Gn _ f~1
i,j±l/2,k = G
Q
i,j±l/2,k
tn _ tt
Hi,j,k±l/2 = h
Q i,j,k±l/2
Для нахождения потоковых величин к,
GHLL hhll
Gi,j±l/2,k' Hi,j,k±l/2
применяется ТУБ-монотонизированный вариант метода Хартена-Лакса-Ван Лира (ИЬЬ) (см., например, [14]). Функциональный вид
выражений для векторов потоков Р'г
n
i±l/2, j, к'
TT n
Gi ,j±l/2, k
H:
i,j,k±l/2
осно-
ван на обеспечении точного выполнения законов сохранения только на двух разрывах типа ударных волн (не рассматриваются разрывы типа контактных или тангенциальных):
трНЬЬ рп _ \ /VлПЬЬ ли N
Гг±1/2 ,з,к ГЯ,з,к = Ля \Чг±1/2,з, к ЧЯ,з,к) ,
грНЬЬ рп _ л (^лНЬЬ ^лп N
Гг±1/2,з,к ГЬ,з,к = ЛЬ \Чг±1/2,з,к ЧЬ,з,к) .
Отсюда находятся соотношения для потоковых величин ^±Ь/2 з к и ЯН±1/2 з к • Например, вдоль оси х эти потоки через грани расчетной ячейки вычисляются с помощью следующего выражения:
hll
i+1/2
Fl, 0 < Хь;
= < xrfl — xlfr + xlxr (qr — ql)
Л Л , Хь < 0 < Xr ;
Xr — xl
Fr, Xr < 0,
где Fl = F (ql), FR = F (qr). Реконструкция векторов консервативных переменных на гранях ячейки выполнялась согласно выражениям:
Ql = Qi + 0,5 min mod (Aq»-i, Адг) ;
Qr = Qi+i - 0,5 min mod (AQi+i, Aqi);
w ,ч f sign (a)min(Id , Ibl), ab > 0; min mod (a, b) = < n i ^ n
0, ab < 0,
где Aqi = qi — Qi-1, а функция minmod(^) покомпонентно применяется к паре векторных аргументов. Собственные значения Хь, Xr на вертикальных гранях ячейки вычислялись по вектору осредненных значений физических переменных U* = (p*,u*,v*,E*):
Хь = u* — c, XR = u* + c, c =
\
е- — iuül+iv*)! — 1) Y,
где
Рь + Ря * Рьиь + Ря ия
Р = -, и = -,
2 Рь + Ря
* РьУь + РяУя, „ * РьЕь + РяЕя V =-, Е =-.
Рь + Ря Рь + Ря
Вычисление потоковых величин СПНзЬ±ь1/2 к, НППзьк±1/2 по методу Хартена-Лакса-Ван Лира через грани расчетных ячеек вдоль осей г, у осуществлялось аналогично описанному способу.
При аппроксимации диффузионной составляющей векторов потоков СС, Н на грани расчетной ячейки применяется разностная схема типа центральных разностей второго порядка точности. В этом случае
производные вычислялись по формулам:
dU dx
_ Ui+l/2,j,k — Ui-l j = AX
dU_ dy
dU
dz
Ui+l/2,j+l/2,k + U j+l/2,k - Ui+l/2,j-l/2,k - Ui,j-l/2,k
jk = 4Ax
Ui+l/2,j,k+l/2 + U Ы = 4AZ ;
здесь и — вектор неконсервативных зависимых переменных решаемой задачи.
В настоящее время активно развивается направление вычислительной математики, связанное с использованием прямоугольных сеток при расчетах нестационарных аэрогидродинамических течений в областях сложной геометрической формы. В этом случае при дискретизации краевых условий используется представление искомой сеточной функции многомерным полиномом невысокой степени. Наиболее распространенным является использование метода наименьших квадратов, процедур билинейной и биквадратичной интерполяции.
В настоящей работе для решения системы трехмерных нестационарных осредненных по Рейнольдсу уравнений Навье-Стокса используются прямоугольные декартовы сетки (в этом случае отсутствует проблема, связанная с вырождением системы координат) и применен метод погруженной границы с фиктивными ячейками [15, 16] для аппроксимации краевых условий на криволинейных границах областей сложной пространственной формы. При этом предполагается, что восстановление значений газодинамических параметров в фиктивных узлах погруженной расчетной границы проводится в плоскостях ХУ (аналогично работам [15,16]), которые перпендикулярны оси симметрии ГЛА. Процесс определения решения в фиктивных узлах рассмотрим для гиперболической части уравнений Эйлера системы уравнений Навье-Стокса.
Для поиска стационарного решения уравнений Эйлера сформулируем краевые условия на границе расчетной области: считаем, что через левую границу (х = 0) вдувается сверхзвуковой невозмущенный поток, на правой (х = Ь) и верхней (х = К) расчетных границах ставится условие выхода сверхзвукового потока из расчетной области. На нижней границе (у = 0) задается условие симметрии. На поверхности
др дУт
обтекаемого тела задаются условия вида: — =0, Уп = 0, —— = 0,
дп дп
дР рУт2 , р - - -
=-(где К — локальный радиус кривизны криволинейной по-
дп К
верхности). Для аппроксимации краевых условий на криволинейной
границе используются три узла (xi,yi), (Х2,У2), (xG,yG) и одна точка на границе (xO,Уо), как показано на рис.4 [15-17]. На криволинейной границе твердого тела для уравнений Эйлера реализуется условие непротекания. В этом случае нормальная составляющая скорости un равна нулю:
uOO = uO cos 9 + vO sin 9 = 0,
где 9 — угол наклона нормали (см. рис. 4). Тангенциальная составляющая скорости равна uO = -uO sin 9 + vO cos 9. Значения величин uO и vO в точке на границе
(xO ,yO) связаны билинейными интерполяционными соотношениями с известными значениями в двух приграничных узлах (x1, y1), (x2, y2) и искомым значением в фиктивном узле (xG, yG):
i \ -i 1 xi yi
Рис. 4. Используемые при дискретизации краевых условий узлы расчетной сетки и вспомогательные точки
UO = ( 1 XO yo )
Vo = ( 1 Xo yo)
1 X2 У2
1 xG yG
Xi
X2
yi y2
i
1 XG yG
(1 Xo yo ) Пусть
1 Xi yi ^ -i ui
1 X2 y2 I U2
1 XG yG uG
( bi b2 bG ) = ( 1 Xo yo)
cos 0 +
1 Xi yi
1 X2 y2
1 XG yG
sin 9
= 0.
-i
тогда
b
uG cos 0 + vG sin 0 = —(ui cos 0 + vi sin 0) — — (u2 cos 0 + v2 sin 0).
Ь^"1 ' "х 7 Ьс Для однозначного определения искомых значений скоростей ис, vG в фиктивном узле (хс ,ус) необходимы дополнительные уравне-
дит
ния. Для замыкания системы используем краевое условие —— = 0
дп
на криволинейной границе в точке О. Аппроксимация этого краево-
го условия также осуществляется с помощью процедуры билинейной интерполяции:
uG sin 9 — vG cos 9 =
1
= ( 1 Xg yG)
Xi
У1
-1
1 X2 У2 0 cos 9 sin 9
ui sin 9 — vi cos 9 u2 sin 9 — v2 cos 9 (duT/dn)O
Пусть
( di d2 do) = ( 1 Xg yG )
1 xi yi 1 X2 У2 0 cos 9 sin 9
i
тогда uG sin 9—vG cos 9 = d1 (u1 sin 9 — v1 cos 9)+d2 (u2 sin 9 — v2 cos 9) + + do (duT/dn)O.
Введем вспомогательные обозначения:
RHS1 = — (ui cos 9 + v\ sin 9) — —2 (u2 cos 9 + v2 sin 9);
Og bG
RHS2 = d1 (u1 sin 9 — v1 cos 9) + d2 (u2 sin 9 — v2 cos 9). Получим два уравнения:
uG^1 cos 9 + vGt1 sin 9 = RHSi; ukG+1 sin 9 — vGt1 cos 9 = RHS2.
Откуда можно получить
u
k+l G
k+i G
= RH Si cos 9 + RHS2 sin 9; = RH Si sin 9 — RHS2 cos 9.
Величина ро определяется из краевого условия на границе др/дп = 0, которое аппроксимируется линейным соотношением
1 xi
yi
i
Pi
Pg = ( 1 Xg Ус) ( 1 X2 У2
0 cos 9 sin 9
р2
(др/дп)о
Значение Ео находим из краевого условия для давления на границе др/дп = р (ит)2 /К (где К — локальный радиус кривизны криволинейной поверхности):
р0 I 1( 2 , 2\
— + 2 (ио + ; 1
Eg =
PG (7 — 1) 1 xi
У1
Рс XG УС) • | 1 Х2 У2
0 cos 9 sin 9
Pi P2
(др/dn)o
Величину uTo в точке O на границе можно получить из уравнения
и,
o
1 xi yi \ i / ui sin 9 — vi cos 9 = ( 1 xo y o ) [ 1 x2 y2 I I u2 sin 9 — V2 cos 9 0 cos 9 sin 9 J \ (duT/dn)O
а величина ро может быть определена из уравнения
-i
1
Xi
yi
Pi
Po = ( 1 Xo yo) I 1 X2 y2
0 cos 0 sin 0
р2
(дР/дп)о
Для сохранения устойчивости необходимо предусмотреть механизм ограничения компонент вектора наклонов физических переменных
/ Рс \ ( Ро\
ис и = По
Ус ' 0 Уо
\ Ес ) \ Ео
С этой целью вводится дополнительная точка-образ I с координата-
= 2( хо
Ут ) V Уо
Sog = Ug — Uo, где UG =
ми (xj,yi), которые определяется выражением
XG Ус
. Далее на отрезке IO с помощью процедуры интерполяции
вычисляется вектор наклонов физических переменных по их значениям в узлах (Х1/У1), (х2,У2),. • • ,(х5,Уб):
Sro =
( Pi p2 p3 p4 p5
ui ui u3 u4 u5
Vi V2 V3 v4 v5
\ Ei e2 e3 E4 E5 )
х
х
11111
Xi X2 X3 X4 X5
yi y2 Уз y4 y5
-i
2 Jj i
2 x2
2
x3
X2 X2 45
v у2 у1 y¡ у4 У1 /
0
Xo — Xr
yo — yr
22 dj o ж j
\ yo — yJ J
2
Полученные наклоны Sog и Sio используются для восстановления окончательного вектора физических переменных UG в фиктивном узле (xG,yG): UG = UO + min mod (Sog,Sio ). Представленный метод погруженной границы с фиктивными ячейками, который используется
для расчета нестационарных, осредненных по Рейнольдсу уравнений Навье-Стокса на прямоугольных сетках, адаптируется к особенностям течения и границе обтекаемого тела генерацией и сгущением сетки по методу AMR [18-21]. В этом случае создается структурированная прямоугольная сетка, каждая расчетная ячейка которой представляется вершиной дерева и имеет заданный уровень иерархии.
Результаты численного моделирования. Расчетная область выполнена в виде прямоугольного параллелепипеда, внутри которого помещена модель ГЛА Waverider. При проведении численных расчетов использовалась адаптивная к границам ГЛА сетка, построенная в декартовой системе координат. Характерное число узлов регулярной расчетной сетки составляло 3 • 106. На левую границу параллелепипеда (вход в расчетную область) вдувался однородный поток воздуха, параметры которого соответствовали высоте H = 30 км от поверхности Земли. Расчеты проведены для углов атаки от 0° до 5°, числа Маха M = 6 (в набегающем потоке) и при следующих значения газодинамических параметров на входе: температура 227 K; давление 1200 Па; плотность 0,0184 кг/м3, скорость звука 302 м/с; динамическая вязкость 0,148 • 10-4 кг/(м-с). На рис. 5-8 (3-я и 4-я полосы обложки) приведены результаты расчетов, выполненных с помощью приведенной методики.
На рис. 5 и 7 показаны распределения чисел Маха для углов атаки 0° и 5°. Рис.6 и 8 соответствуют распределениям поля температуры при аналогичных углах атаки.
Приведенные в настоящей работе графические распределения могут быть описаны в виде следующей картины течения: вблизи носа и нижней поверхности ГЛА возникает головной скачок уплотнения и начинается формирование пограничного слоя, в области перехода от носовой поверхности к задней части ГЛА находится область волн разрежения (для углов атаки, больших 0), далее по потоку наблюдается "горло" следа (донная область течения) с примыкающими к нему хвостовыми скачками, замыкается течение ближним и дальним следом.
Из приведенных распределений следует, что экстремальные значения температуры и давления наблюдаются в носовой (T & 1300 K, p & 5500 Па) и донной областях течения (T & 1100 ... 900 K, p & 3 Па). В зависимости от угла атаки температура воздуха, обтекающего верхнюю поверхность ГЛА Waverider, достаточно сильно изменяется: для угла атаки а = 0 температура достигает значения T & 900 K; а для угла атаки 5 ° — значения T ~ 500 K. В обоих случаях наблюдается значительная по пространственным размерам область следа, в которой температура воздуха T & 1000 K. Здесь отметим, что для более правильного описания картины течения необходимо рассматривать газ, как многокомпонентную химически реагирующую смесь.
Иллюстрации к статье М.А. Котова, В.В. Кузенова
«ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ОБТЕКАНИЯ ПОВЕРХНОСТЕЙ ПЕРСПЕКТИВНЫХ ГИПЕРЗВУКОВЫХ ЛЕТАТЕЛЬНЫХ АППАРАТОВ»
Л
Рис.5. Распределение числа Маха вблизи поверхности экспериментального гиперзвукового летательного аппарата Waverider для угла атаки а = 0°
г л
Иллюстрации к статье М.А. Котова, В.В. Кузенова
«ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ОБТЕКАНИЯ ПОВЕРХНОСТЕЙ ПЕРСПЕКТИВНЫХ ГИПЕРЗВУКОВЫХ ЛЕТАТЕЛЬНЫХ АППАРАТОВ»
Рис.7. Распределение числа Маха вблизи поверхности экспериментального гиперзвукового летательного аппарата Waverider для угла атаки а = 5°
Рис. 8. Температурное поле вблизи поверхности экспериментального гиперзвукового летательного аппарата Waverider для угла атаки а = 5°
Выводы. Проведено численное моделирование задачи о натекании сверхзвукового вязкого потока газа на ГЛА типа Waverider. Определены аэротермодинамические параметры газа в окрестности носа и донной области ГЛА.
Изложены различные способы геометрического моделирования трехмерных объектов и методы создания сложных моделей поверхностей перспективных ГЛА в системе САПР SolidWorks. Приведены примеры построения модели поверхности ГЛА X-43 и примеры построения моделей поверхностей гиперзвуковых летательных аппаратов Waverider, Х-51А.
Работа проводилась в рамках программ фундаментальных исследований РАН и выполнена в Лаборатории радиационной газовой динамики Института проблем механики им. А.Ю. Ишлинского РАН (ИПМех РАН).
СПИСОК ЛИТЕРАТУРЫ
1. S u r z h i k o v S. T. 2D CFD/RGD Model of space vehicles //Proc. of the Int. Workshop on Radiation of High Temperature Gases in Atmospheric Entry. October 2003, Lisbon, Portugal, European Space Agency, SP-533. - 2003. - P. 95-102.
2. S u r z h i k o v S. T. Numerical Simulation of Heat Generated by Entry Space Vehicle // AIAA 2004-2379. - 2004. - 11 p.
3. Голованов Н. Н. Геометрическое моделирование. - М.: Физматлит, 2002. - 472 с.
4. R e u b u s h D. E., N g u y e n L. T., R a u s c h V. L., Review of X-43A Return to Flight Activities and Current Status // AIAA 2003-7085. - 2003. - 12 p.
5. E n g e l u n d W. C., Holland S. D., C o c k r e 11 C. E. et all. Propulsion System Airframe Integration Issues and Aerodynamic Database Development for the Hyper-X Flight Research Vehicle // ISOABE 99-7215. - 1999. - 12 p.
6. Charles E. Cockrell. Interpretation of Waverider performance data using computational fluid dynamics // AIAA Paper 93-2921. - 1993. - 11 p.
7. CD X-51A First Flight May 2010. V1.exb1. Pratt & Whitney A United Technologies Company.
8. SolidWorks. Практическое руководство. - М.: Бином, 2004.
9. Котов М. А., Кузенов В. В. Моделирование сложных поверхностей гиперзвуковых летательных аппаратов в САПР SolidWorks 2007 // Труды Всеросс. шк.-семинара "Аэрофизика и физическая механика классических и квантовых систем". - М., 2009. - С. 271-273.
10. К о т о в М. А., Кузенов В. В. Геометрическое моделирование сложных поверхностей гиперзвуковых летательных аппаратов в САПР // Сб. научн. трудов АФМ 2010.
11. Же л е з няко в а А. Л., Суржиков С. Т. Численное моделирование гиперзвукового обтекания модели летательного аппарата Х-43 / Препр. ИПМех им. А.Ю. Ишлинского РАН. - № 950. - 2010. - С. 20-25.
12. C a n a n W. J. Breathing new hope into hypersonics // AEROSPASE, November 2007. - 26 p.
13. Turbulence modeling methods for the compressible. AIAA Paper. - 1983. -No. 83-1693. Coacley T.J.
14. Куликовский А. Г., Погорелов Н. В., Семенов А. Ю. Математические вопросы численного решения гиперболических систем уравнений. -М.: Физматлит, 2001. - 607 с.
15. Мортиков Е. В. Применение метода погруженной границы для решения уравнений Навье-Стокса в областях сложной конфигурации // Вычислительные методы и программирование. - 2010. - Т. 11. - С. 32-42.
16. ВинниковВ.В., Ревизников Д. Л. Применение декартовых сеток для решения уравнений Навье-Стокса в областях с криволинейными границами // Математическое моделирование. - 2005. - Т. 17. - С. 15-30.
17. Винников В. В. Алгоритмическое описание метода погруженной границы с фиктивными ячейками // Вестник МАИ. - 2009. - Т. 16. № 6. - С. 193-198.
18. B e г g e г M., öliger J. Local adaptive mesh refinement for hyperbolic partial differential equations // J. Сотр. Phys. - 1984. - No. 53. - P. 484-512.
19. B e r g e r M. C o l e 11 a P. Local adaptive mesh refinement for shock hydrodynamics // J. Сотр. Phys. - 1989. - No. 82. - P. 64-84.
20. Супалов А. В. Разработка алгоритмов построения трехмерных локально-адаптивных сеток и их программная реализация. Дисс. . . . на соискание уч. степ. канд. физ.-мат. наук. ИВМ РАН. - М., 1993. - С. 1-100.
21. Bell J. et al. A three dimensional adaptive mesh refinement for hyperbolic conservatiin laws // SIAM. J. Sci. Statist. Comput. - 1994. - No. 15A. - P. 127-138.
Статья поступила в редакцию 26.02.2012
Михаил Алтаевич Котов окончил МИРЭА в 2011 г. Научный сотрудник лаборатории "Радиационная газовая динамика" ИПМех РАН. Автор шести научных работ в области аэромеханики и теплофизики.
M.A. Kotov graduated from the Moscow Institute of Radio Engineering, Electronics and Automatics in 2011. Researcher of Radiative Gas Dynamics Laboratory of the Institute of Applied Mechanics of the Russian Academy of Sciences. Author of six publications in the field of aeromechanics and thermal physics.
Виктор Витальевич Кузенов окончил в 1983 г. МГУ им. М.В. Ломоносова. Канд. техн. наук, старший научный сотрудник лаборатории "Радиационная газовая динамика" ИПМех РАН, доцент кафедры "Теплофизика" МГТУ им. Н.Э. Баумана. Автор более 110 научных работ в области теплофизики и радиационной газовой динамики.
V.V. Kuzenov graduated from the Lomonosov Moscow State University in 1963. Ph. D. (Eng.), senior researcher of Radiative Gas Dynamics Laboratory of the Institute of Applied Mechanics of the Russian Academy of Sciences, assoc. professor of "Thermal Physics" department of the Bauman Moscow State Technical University. Author of more than 110 publications in the field of thermal physics and radiative gas dynamics.