Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

DEVELOPMENT OF MONTE-CARLO METHODS FOR SOLVING AERODYNAMICS PROBLEMS OF REENTRY SPACE VEHICLES

Zharov V.A. 1, 2 Zay Yar Myo Myint 2 Polyakov M.S. 2 Khlopkov A.Y. 2 Kyaw Zin 2
1 Central Aerohydrodynamic Institute
2 Moscow Institute of Physics and Technology (state university)
The most service time of aerospace vehicles is at high-altitude, in free molecular flow conditions and the experimental investigation for its moving in this condition is quite problematic. Significance of rarefied gas dynamic is explained in practical importance to solve a wide range of the problems connected with the present stage of space exploration, development of vacuum and nano technologies. The first application of statistical methods connected with direct simulation of rarefied gas flows and the direct simulation of Monte-Carlo (DSMC) method was appeared the most effective. The purpose of this work is to develop aerodynamic characteristics calculation method for space vehicles in high-speed rarefied gas flow. Calculation results of aerodynamic characteristics for reentry vehicle are presented. The results may be useful for perspective space vehicle project.
Monte-Carlo method
rarefied gas dynamics
aerospace vehicle
high-altitude hypersonic aerodynamics
1. Barantsev R.G. Vzaimodeistvie razrezhennykh gazov s obtekaemymi poverkhnostyami. Moscow, Science, 1975.
2. Belotserkovskii O.M., Khlopkov Yu.I. Metody Monte-Karlo v mekhanike zhidkosti i gaza. Moscow, Azbuka, 2008.
3. Voronich I.V., Zay Yar Myo Myint Vliyanie osobennostei vzaimodeistviya gaza s poverkhnost’yu na aerodinamicheskie kharakteristiki kosmicheskogo apparata // Vestnik MAI. 2010. Vol. 17, nо. 3. рр. 59–67.
4. Zay Yar Myo Myint, Khlopkov A.Yu., Kyaw Zin Osnovnye podkhody k postroeniyu metodov Monte-Karlo v vychislitel’noi aerodinamike // Trudy MAI. 2011, nо. 42. 17 p.
5. Kogan M.N. Dinamika razrezhennogo gaza. Moscow, Science, 1967.
6. Khlopkov Yu.I. Statisticheskoe modelirovanie v vychislitel’noi aerodinamike. Moscow, MIPT, 2006.
7. Khlopkov Yu.I., Zay Yar Myo Myint, Khlopkov A.Yu., Kyaw Zin Metody Monte Carlo dlya opredeleniya aerotermodinamicheskikh kharakteristik giperzvukovykh vozdushno kosmicheskikh sistem // Materials digest of LI International Research and Practice Conference «Physical, Mathematical and Chemical Sciences: Theoretical Trends and Applied Studies». London: IASHE. 2013. рр. 41–44.
8. Bird G.A. Molecular Gas Dynamics and the Direct Simulation of Gas Flows. Oxford: Clarendon Press, 1994.
9. Cercignani C., Lampis M. Kinetic Models for Gas-Surface Interactions // Transport Theory and Statistical Physics. 1971. Vol. 1, no. 2. рр. 101–114.
10. Haviland J.K., Lavin M.L. Application of the Monte Carlo Method to Heat Transfer in a Rarefied Gas // Phys. Fluids. 1962. Vol. 5, no. 11.
11. Ketsdever A.D., Muntz. E.P. Gas-Surface Interaction Model Influence on Predicted Performance of Microelectromechanical System Resistojet // Journal of Thermophysics and Heat Transfer. 2001. Vol. 15, no. 3. pp. 302–307.
12. Lord R.G. Application of the Cercignani-Lampis Scattering Kernel to Direct Simulation Monte Carlo Calculations // Proc. of 17th Int. Symp. on Rarefied Gas Dynamics. 1991. pp. 1427–1433.
13. Lord R.G. Some Further Extensions of the Cercignani-Lampis Gas-Surface Interaction Model // Phys. Fluids. 1995. Vol. 7, no. 5. pp. 1159–1161.
14. Padilla J.F. Assessment of Gas-Surface Interaction Models for Computation of Rarefied Hypersonic Flows // Ph.D. Dissertation. University of Michigan, 2008.
15. Padilla J.F., Boyd I.D. Assessment of Gas-Surface Interaction Models in DSMC Analysis of Rarefied Hypersonic Flow // AIAA Paper 2007–3891. 2007.
16. Santos W.F.N. Gas-Surface Interaction Effect on Round Leading Edge Aerothermodynamics // Brazilian Journal of Physics. 2007. Vol. 37, no. 2A.
17. Utah S. and Arai H. Monte Carlo Simulation of Reentry Flows Based Upon a Three-Temperature Model // Proc. of 23rd Int. Symp. on Space Technology and Science. 2002. Vol. 1. pp. 1209–1214.
18. Wadsworth D.C., Van Glider D.B., Dogra V.K. Gas-Surface Interaction Model Evaluation for DSMC Applications // Proc. of 23rd Int. Symp. on Rarefied Gas Dynamics. 2003. pp. 965–972.

Многоразовая воздушно-космическая система предназначена для решения широкого круга задач в космосе, в том числе:

  • Выведение на околоземную орбиту и возврат с орбиты людей и полезных грузов;
  • Транспортно-техническое обеспечение космических объектов различного назначения;
  • Проведение аварийно-спасательных работ на орбите;
  • Проведение научно-технических и технологических экспериментов в космосе;
  • Осуществление международного контроля за космическим пространством.

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

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

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

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

Целью данной работы является определение метода вычисления АДХ ВКА в высокоскоростном потоке разреженного газа. В работе представлены результаты расчетов аэродинамических характеристик перспективных воздушно-космических аппаратов, полученные методом прямого статистического моделирования (ПСМ).

Метод и алгоритм моделирования АДХ ВКА

Развитие численных методов в динамике разреженных газов связано в первую очередь с использованием метода прямого статистического моделирования процессов, описываемых кинетическим уравнением Больцмана [2, 6, 7],

garov01.wmf

Здесь f = f(t, x, y, z, xx, xy, xz) – функция распределения молекул по времени, координатам и скоростям. f, f1, f′, f1′ – функции распределения, соответствующие скоростям пары частиц до и после столкновения, garov02.wmf – относительная скорость, b – прицельное расстояние, e – азимутальный угол в плоскости, перпендикулярной плоскости столкновения. В практической реализации методы прямого статистического моделирования, основанные на подходах Бёрда [8] (моделирование динамики ансамбля молекул) и Хэвиленда [10] (моделирование индивидуальных траекторий молекул), оказались наиболее эффективными и их модификации с переменным успехом осуществляли победное шествие по вычислительной аэродинамике. Рассмотрим моделирование свободномолекулярного течения методом Монте-Карло. Алгоритм метода Монте-Карло выглядит следующим образом:

  1. Ввод данных.
  2. Выбор границы, которой стартует частица.
  3. Вычисление координат точки влета частицы в область.
  4. Вычисление скорости частицы.
  5. Вычисление координаты точки пересечения траектории частицы с поверхностью тела.
  6. Вычисление импульса и энергии, приносимых частицей.
  7. Вычисление скорости отраженной частицы.
  8. Вычисление импульса и энергии отраженной частицы.
  9. Выполнение пп. 4–7 до покидания молекулой расчетной области.
  10. Осреднение данных.

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

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

Роль законов взаимодействия молекул с поверхностью проявляется тем сильнее, чем более газ разрежен [5]. Граничным условиями для уравнения Больцмана являются условия, ввязывающие функцию распределения падающих и отраженных молекул [1].

В модели Максвелла плотность распределения отраженных молекул имеет вид

garov03.wmf

и ядро рассеяния [1, 2] имеет вид

garov04.wmf; garov05.wmf

Здесь полагается, что доля (1 – st) молекул отражается зеркально, а остальная часть st молекул – диффузно, параметр 0 ≤ st + 1 определяет коэффициент аккомодации касательной компоненты импульса st = (Pti – Ptr)/Pti.

Компоненты вектора скорости при диффузном отражении моделируются в локальной сферической системе координат, ось которой направлена вдоль вектора внешней нормали к поверхности, с помощью выражений [8]

garov06.wmf garov07.wmf

garov08.wmf

где a1, a2, a3, a4 – независимые случайные числа, равномерно распределенные в интервале (0, 1), q и j – полярный и азимутальный углы.

К. Черчиньяни и М. Лампис предложили феноменологическую модель (CL), которая также удовлетворяет принципу взаимности и является усовершенствованием максвелловской модели [9]. Модель основана на введении двух параметров, которые представляют собой коэффициент аккомодации sn = sEn по кинетической энергии, связанной с нормальной компонентой скорости, и коэффициент аккомодации касательной компоненты импульса st. Модель CL хорошо соответствует результатам лабораторных исследований с высокоскоростными молекулярными пучками [2]. Хотя сравнение ограничено лабораторными условиями, модель CL является теоретически обоснованной и относительно простой. Позднее появились модификации ядра рассеяния модели CL [12], однако они дают незначительное улучшение при сравнении с лабораторными экспериментами. В общем случае модель взаимодействия имеет несколько параметров произвольного физического смысла, которые позволяют добиться разумного согласия с результатами лабораторных исследований в некотором диапазоне условий. Универсальная модель должна использовать ядро рассеяния, полученное на основе физического эксперимента в широком диапазоне чисел Кнудсена и скоростей потока [2].

Использованное преобразование расширяет CL модель для учета обмена вращательной энергией между газом и поверхностью [12]. Модель в таком виде называется моделью Черчиньяни‒Лампис‒Лорда (CLL). Потом были предложены модификации модели [13] для учета обмена колебательной энергией и расширения диапазона состояний рассеянных молекул. Модель CLL в настоящее время получила широкое признание в работах [3, 11, 14–18].

Результаты расчетов

Рассмотрим результаты расчета коэффициентов силы сопротивления и подъемной для возвращаемого аппарата в свободномолекулярном потоке разреженного газа. Используются различные модели взаимодействия молекул с поверхностью (Максвелла и Черчиньяни‒Лампис‒Лорда, CLL). Представлены результаты расчета различным моделями взаимодействия газа с поверхностью (Максвелла и CLL) методом Монте-Карло. Значения параметров: температурный фактор tw = Tw/T∞ = 0,1; скоростное отношение s = 20; коэффициенты аккомодации тангенциального импульса и нормальной энергии στ, σn = 0,5. Расчет проводился с использованием 5·106 частиц.

pic_56.tif

Рис. 1. Геометрическое представление варианта возвращаемого аппарата

На рис. 2 представлены зависимости коэффициентов силы сопротивления Cx, подъемной силы Cy от угла атаки a от 0 до 90° для возвращаемого аппарата (рис. 1).

При уменьшении коэффициент аккомодации στ от 1 до 0.5, величина Cx увеличивается при 0° < α < 55° и значение достигает до 2,72 при α = 0°. Коэффициент Cy снижает в несколько раз по модулю при уменьшении στ от 1 до 0.5 при 0° < α < 80°. Можно объяснить, что при нулевой аккомодации все молекулы отражаются зеркально, и полной аккомодации отражаются диффузно. Зеркальные отраженные молекулы передают поверхности больший импульс, чем диффузно рассеянные от холодной стенки молекулы.

pic_57.tif pic_58.tif

Рис. 2. Зависимости Cx(α) и Cy(α) при στ = 0,5 и 1 (tw = 0,1)

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

Заключение

Представлены результаты расчетов аэродинамических сил сопротивления Cx, подъемной Cy возвращаемого аппарата методом Монте-Карло при различных значениях коэффициентов аккомодации с использованием различных моделей взаимодействия молекул с поверхностью. Исследовано влияние на АДХ особенностей модели взаимодействия молекул с поверхностью. Разработанные программные системы позволяют оперативно получать АДХ разрабатываемых и эксплуатируемых ВКА на орбите и на начальном участке траектории спуска и могут быть использованы при проектировании перспективных космических аппаратов.

Работа выполнена при поддержке РФФИ (Грант № 14-07-00564-а).

Рецензенты:

Боголепов В.В., д.ф.-м.н., действительный член РАЕН, главный научный сотрудник, ФГУП ЦАГИ, профессор МФТИ, г. Жуковский.

Липатов И.И., д.ф.-м.н., чл.-корр. РАН, действительный член РАЕН, начальник отдела ФГУП ЦАГИ, профессор МФТИ, г. Жуковский.

Работа поступила в редакцию 17.01.2014.