Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,441

РАЗРАБОТКА МЕТОДОВ МОНТЕ-КАРЛО ДЛЯ РЕШЕНИЯ ЗАДАЧ АЭРОТЕРМОДИНАМИКИ ВОЗВРАЩАЕМЫХ КОСМИЧЕСКИХ АППАРАТОВ

Жаров В.А. 1, 2 Зея Мьо Мьинт 2 Поляков М.С. 2 Хлопков А.Ю. 2 Чжо Зин 2
1 ФГУП ЦАГИ «Центральный аэрогидродинамический институт»
2 ФГОУ ВПО «Московский физико-технический институт (государственный университет)»
Большую часть срока службы воздушно-космический аппарат находится на большой высоте, при свободномолекулярных условиях и экспериментальное исследование его движения довольно проблематично. Значение динамики разреженных газов обусловлено практической важностью решения широкого круга задач, связанных с современным этапом освоения космоса, развитием вакуумной, компьютерной и нано технологий. Первое применение статистических методов связывалось с непосредственным моделированием течений разреженных газов и методы прямого статистического моделирования Монте-Карло оказались здесь наиболее эффективными. Целью данной работы является разработка методов вычисления аэродинамических характеристик воздушно-космических аппаратов в высокоскоростном потоке разреженного газа. Представлены результаты расчетов аэродинамических характеристик спускаемых аппаратов. Полученные результаты могут быть использованы при проектировании перспективных космических аппаратов.
метод Монте-Карло
динамика разреженного газа
воздушно-космический аппарат
высотная гиперзвуковая аэродинамика
1. Баранцев Р.Г. Взаимодействие разреженных газов с обтекаемыми поверхностями. – М.: Наука, 1975.
2. Белоцерковский О.М., Хлопков Ю.И. Методы Монте-Карло в механике жидкости и газа. – М.: Азбука, 2008.
3. Воронич И.В., Зея Мьо Мьинт Влияние особенностей взаимодействия газа с поверхностью на аэродинамические характеристики космического аппарата // Вестник МАИ. – 2010. – Т. 17, № 3. – С. 59–67.
4. Зея Мьо Мьинт, Хлопков А.Ю., Чжо З. Основные подходы к построению методов Монте-Карло в вычислительной аэродинамике // Труды МАИ. –2011. – № 42. – 17 с.
5. Коган М.Н. Динамика разреженного газа. – М.: Наука, 1967.
6. Хлопков Ю.И. Статистическое моделирование в вычислительной аэродинамике. – М.:МФТИ, 2006.
7. Хлопков Ю.И., Зея М.М., Хлопков А.Ю., Чжо З. Методы Монте-Карло для определения аэротермодинамических характеристик гиперзвуковых воздушно космических систем // Materials digest of LI International Research and Practice Conference «Physical, Mathematical and Chemical Sciences: Theoretical Trends and Applied Studies». – London: IASHE. – 2013. – P. 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, № 2. – P. 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, № 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, № 3. – P. 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. – P. 1427–1433.
13. Lord R.G. Some Further Extensions of the Cercignani-Lampis Gas-Surface Interaction Model // Phys. Fluids. – 1995. – Vol. 7, № 5. – P. 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, № 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. – P. 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. – P. 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.


Библиографическая ссылка

Жаров В.А., Жаров В.А., Зея Мьо Мьинт, Поляков М.С., Хлопков А.Ю., Чжо Зин РАЗРАБОТКА МЕТОДОВ МОНТЕ-КАРЛО ДЛЯ РЕШЕНИЯ ЗАДАЧ АЭРОТЕРМОДИНАМИКИ ВОЗВРАЩАЕМЫХ КОСМИЧЕСКИХ АППАРАТОВ // Фундаментальные исследования. – 2013. – № 11-9. – С. 1819-1823;
URL: http://fundamental-research.ru/ru/article/view?id=33463 (дата обращения: 22.09.2020).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074