РАВНОВЕСИЕ И КОЛЕБАНИЯ СВОБОДНОЙ ПОВЕРХНОСТИ ЖИДКОГО ТОПЛИВА В КОАКСИАЛЬНО-ЦИЛИНДРИЧЕСКИХ СОСУДАХВ УСЛОВИЯХ МИКРОГРАВИТАЦИИ

УДК 532.65/534-141

БОТ: 10.18698/2308-6033-2021-8-2099

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

© Юй Чжаокай, А.Н. Темнов МГТУ им. Н.Э. Баумана, Москва, 105005, Россия

Представлено решение задачи о равновесии и малых колебаниях идеальной жидкости в условиях микрогравитации. Дана количественная оценка влияния таких параметров, как угол смачивания, число Бонда, соотношение радиусов внутренней и внешней стенки сосуда, глубина жидкости в нем. Для сосудов в форме коаксиального цилиндра получены выражения потенциала скоростей жидкости и поля смещения свободной поверхности в виде ряда Бесселя. Дополнительно к аналитическим и экспериментальным данным, приведенным в литературе, доказана достоверность разработанного численного алгоритма и сделан вывод о том, что для Во > 5 и 5 > 0,2 при неизменном физическом состоянии смачиваемой поверхности форма свободной поверхности приближается к плоской и угол смачивания мало влияет на собственные частоты колебаний свободной поверхности жидкости. Полученные результаты будут полезны при решении задач по определению гидродинамических характеристик движения жидкого топлива в условиях космического пространства.

Ключевые слова: микрогравитация, сила поверхностного натяжения, капиллярная жидкость, коаксиальный цилиндр, метод Рунге — Кутты, ряд Бесселя

Введение. В баках современных космических аппаратов содержится значительная часть жидкого топлива. При проектировании подобных транспортных средств инженерам и конструкторам полезно знать конфигурацию и поведение поверхности раздела жидкости — газа в условиях микрогравитации, для того чтобы более точно описать динамические процессы в космических аппаратах при выполнении различных маневров. Следует заметить, что плескание жидкости может приводить к изменению положения центра масс, а значит, повлиять на устойчивость передвижения космических аппаратов. В российских и зарубежных монографиях [1-3] обобщены результаты исследования статики и динамики жидкости в условиях, близких к невесомости. Анализ приведенных в статьях [4-6] методов показывает, что в современной гидромеханике невесомости недостаточно полно разработаны методы численного моделирования поведения жидкости в двухсвязных сосудах. В статьях [7-10] представлены приближенные методы вычисления частот и форм колебаний жидкости в сосудах, имеющих форму кругового цилиндра и сферы, кото-

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

Цель настоящей работы — представить исследование равновесия и колебаний жидкости с учетом силы поверхностного натяжения в сосудах, имеющих форму коаксиального цилиндра. В статьях [11, 12] было рассмотрено решение задачи о равновесии жидкости в коаксиальном цилиндре в условиях невесомости и микрогравитации, но без учета влияния внутренней стенки сосуда, а в [13-16] приведено исследование колебаний капиллярной жидкости в круговом цилиндре на основе ряда Бесселя.

Определение форм равновесной свободной поверхности. В условиях микрогравитации = (10-6 -10-4)g0) поведение жидкого топлива определяют силы поверхностного натяжения, иначе — межмолекулярные силы на границе двух фаз. Вектор ускорения g действует параллельно продольной оси симметрии сосуда. Введем цилиндрическую систему координат Ог 0г, в которой будет представлено поперечное сечение сосуда (рис. 1). Используем длину дуги ^ в качестве переменной для описания формы свободной поверхности.

ё о Ч <хГ С -У г

1 111111111 ,¡,¡,¡,¡,¡,1,1,1,1 п — — - о I—Г—I о

Рис. 1. Основные обозначения параметров жидкости и система координат Ог0х:

2тт — наименьшая высота свободной поверхности; гтах — наибольшая высота свободной

поверхности; 20(г)— равновесная свободная поверхность; /(г,0,/) — отклонение возмущенной свободной поверхности от равновесной; у — линия смачивания; ^ — область,

которую занимает жидкость; г1 — радиус внутренней стенки сосуда; Г2 — радиус внешней

стенки сосуда; И0 — глубина жидкости около внутренней стенки сосуда

Для того чтобы приступить к исследованию колебаний капиллярной жидкости, необходимо предварительно решить задачу о форме равновесной свободной поверхности. Вывод условий равновесия гидромеханической системы газ — жидкость — твердая стенка и алгоритм решения задачи подробно рассмотрены в [17].

Запишем уравнение невозмущенной свободной поверхности жидкости в безразмерной форме:

2 Н р = Во + С,

где Нр = (к1 + к2) /2 — средняя кривизна равновесной свободной поверхности (к1 и к2 — главные кривизны поверхности); Во = р£г22 /о — число Бонда, характеризующее соотношение массовой силы и силы поверхностного натяжения ( р — плотность жидкости, о — коэффициент поверхностного натяжения раздела жидкости — газа); С — константа, определяемая физическими параметрами задачи и формой сосуда, которая зависит от начала системы координат.

В цилиндрической системе координат От§г средняя кривизна

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

( \

= — +

т

2H р = ЧТ

r dr

+ z2r J

где ( )т — первая производная функции по координате т; ( )* и ( )** — первая и вторая производные функции по длине дуги * .

Используя уравнения (1)-(2), получим систему дифференциальных уравнений, описывающих равновесную свободную поверхность:

При произвольном значении числа Bo система уравнений (3) не имеет аналитического решения и решается методом Рунге — Кутты при условиях:

r (0 ) = 5; z (0 )= 0; u (0 ) = sin (a0); v (0) = - cos (a0);

r

где 5 = r¡ / r2; a0 — угол смачивания; s0 — общая длина дуги.

При интегрировании системы уравнений (3) с граничными условиями (4) необходимо определять константу C, которая при Bo = 0 имеет аналитическую формулу [1] C = 2 cos a0 / (1 - 5), а при других

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

г

0,6 0,4 0,2 0 -0,2 -0,4 -0,6 -0,8

0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 г

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

-а =10° -а = 30° / -а = 60° / -а = 90° /

1 1 1

1111

0,6 0,4 0,2 0 -0,2 -0,4 -0,6 -0,8

0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 г

Рис. 3. Формы равновесной свободной поверхности при различных значениях числа Бонда и угле смачивания ао = 10°

-Во = 1 - Во = 5 > - Во = 10 /

| | 1 ____> 1111

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

уравнению Лапласа

АФ = 0

4 д2 д2 1 д 1 д2

где А = —- +--- +----+ —--2--оператор Лапласа в цилиндриче-

дг2 дг г дг г д02

ской системе координат; Ф — потенциал скоростей жидкости;

условию непротекания на смоченных поверхностях

дФ дФ

£- = 0

дг дг

Предположим, что возмущенная свободная поверхность ¥ будет представлена в виде

¥

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

д¥

^ = д¥ + уф .VF = 0,

Жг дг

^ д _ д _ 1 д _

где V = — е7^ +--ег +---е0 — оператор 1 амильтона.

дг дг г д0

Из уравнений (7)-(8) следует, что линеаризованное кинематическое условие на свободной поверхности имеет вид

д± + . -дФ = 0.

дг дг Жг дг

Используя интеграл Коши — Лагранжа

дФ Р р0

- + - + gz = —,

дг р р

где Р = Рст + Рд — полное давление жидкости на возмущенной свободной поверхности (Рст — статическое давление, Рд — динамическое давление); Р0 — статическое давление жидкости при г = 0, получим динамическое условие на свободной поверхности.

Из уравнения Эйлера для покоящейся жидкости определим давление на равновесной свободной поверхности:

Р0 - Рст = Pgz0 (г).

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

Рг - Р = 2аИв

где Рг — давление газа; Ив — кривизна, определяемая уравнением

2И в =

1Ж

г Жг

гг0г

1 д

+--

г дг

г!г

/ 2 \3/2

_1 А

г2 д0

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

Рг - Рст = О

1Ж

г Жг

гг0г

/ ,,N.1/2

а с помощью (12)-(14) получаем выражение для динамического давления:

Рд =-О

д

дг

г!г

I о\3/2

_1 д г2 д0

Из уравнений (7), (10)—(11) и (15) выведем динамическое условие на свободной поверхности с учетом влияния силы поверхностного натяжения:

дФ

1 д

— + &-Ч-^г

дг

г дг

г!г

I о\3/2

_1 А.

г2 д0

> = 0.

Введем характерный размер сосуда г2 , а также характерные значения для времени г* и потенциала Ф*:

г* =

д/рг23 / о; Ф* = г22 / г*.

Подставив величины х = г2 • х, г = г* • г, Ф = Ф* • Ф в уравнения (5)—(6), (9) и (16), получим безразмерную формулировку задачи, в которой опускаем знак « - » над буквами:

ЛФ = 0(Q) , ^ = 0 (Г! = 5, r2 = 1), — = 0 (Г! = 5, r2 = 1); dr dr

дФ nf [í^ df дФ dz0 дФ , ч

__ = 0(z = -H), - + -• ----0 (z = .-0(r)),

5Ф

— + B0 f- dt

1 д r dr

rfr

i ^3/2

_1 A

r2 de

i ^1/2

= 0 ( z = Zo(r)),

где H = h0 / r2.

Определение собственных частот и форм колебаний. Ограничимся в настоящей статье решениями, отвечающими функции cos e. Решение поставленной задачи для функций Ф(гД z, t) и f(r, e, t), удовлетворяющее уравнению Лапласа и условию непротекания, можно представить в виде ряда Бесселя:

ф(г ,e, z, t) = ZpnA(z nr)cosh (; + H )] cosecos(Wt);

n=1 cosh

f (r, e, t) = Y^nA(Z nr) cos e sin(Wt),

n=1

гда A (Cnr) = J1 (Cnr)- J1iZn) Y1 (znr), J1, Y1 — функции Бесселя

Y1

первого и второго рода первого порядка; числа Zn для различных соотношений радиусов 5 определены из трансцендентного уравнения,

получаемого из условия непротекания [18]; W = ю(pr23 / а) .

Подставив потенциал скоростей и поле смещений в кинематическое и в динамическое условие на свободной поверхности, умножим полученные выражения на функцию rA (Zmr) и проинтегрируем в

диапазоне [5, 1 ]. В результате получим:

Wqm Z nPn;

n=1

B° qm + ^,PmnZ nqn = W ^Y mnpn,

n=1 n=1

где

JVa (Z m,)( A (Z )slnh ^ Я)]

cosh[Zn (z0(r) + H)]} ^

- z0rA (Znr )-^-;-~> dr

cosh

jV Z ) dr );

J/A ( Z mr У

Pmn

cosh

A (Z nr) + 3A(Z nr)z0rz0rr , A (Z ) z0r (l + z0r )3/2 Zn (1 + z02r )5/2 Zlr2 (1 + z02r )3/2

dr

JVa2 ( Zmr ) dr

J/A

Jmn

A

cosh [Zn

dr

J Va2 ( Zmr ) dr

( z0r и г0гг — первая и вторая производные функции z0(г) по г ). Исключим q и получим задачу на собственные значения в виде

Z

n=1

Во а mnZ n +

аkn Zn - W2Ym V k=1 J

pn =

Запишем полученную задачу в матричной форме

(K - W2M ) P = 0,

Л >

где K = В° amnZn + SPmk akn Z3 n

V k=1 J

,M [Чтп ]mxn ' ^ {pn }nx1 .

Допустив m = и = N, можно получить собственные значения

W/2(/ = 1,2,—, N) и собственные векторы {рп ^х1. В дальнейшем

определим размерные собственные частоты ю и поле смещения / жидкости.

Если ао = 90°, то равновесная свободная поверхность плоская:

z0 (r) = z0 r = z0 rr = 0 [tanh (H ), n = m;

R = = J I' n = m;

Pmn = Ymn = I „ n Ф m; 10, n Ф m,

amn

а матрица К — диагональная, М — единичная матрица, и собственное значение = (Бо^ п + С33) tanh (Н) согласуются с аналитическими результатами для цилиндрических сосудов в [ 1, 2].

Обсуждение результатов. В программе МЛТЬЛБ получены формы равновесия и собственные частоты колебаний капиллярной жидкости для различных практически важных ситуаций. Все формы свободной поверхности жидкости показаны для правой части поперечного сечения сосуда. На рис. 2 и 3 видно, что при увеличении угла смачивания или числа Бонда равновесная свободная поверхность приближается к плоской. Следует отметить, что наличие внутренней стенки сосуда уменьшает искривление формы свободной поверхности. Собственные частоты при полной невесомости в зависимости от углов смачивания показаны на рис. 4. Красными звездочками подчеркнуты аналитические значения для угла смачивания ао = 90°. Видно, что численные значения очень хорошо согласуются с аналитическими и доказывается правильность разработанного алгоритма. Черной линией со звездочками показаны экспериментальные результаты при угле смачивания, близком к нулю [20].

личных значениях угла смачивания а0 и условиях Бо = 0, Н = 2

Численные результаты для угла смачивания ао = 10° близки к экспериментальным, и можно сделать вывод о том, что при уменьшении угла смачивания форма свободной поверхности получается более искривленной, а собственные частоты уменьшаются. Это свидетельствует о том, что собственные частоты обратно пропорциональны площади свободной поверхности [18]. Согласно представленному на рис. 5, для числа Бонда Бо > 5 и

5 > 0,2 угол смачивания мало влияет на собственные частоты колебания, а в ракетно-космической технике данные условия часто выполняются, что позволяет оценить собственные частоты с помощью аналитической формулы. При увеличении глубины жидкости собственные частоты возрастают и приближаются к предельным значениям (рис. 6), так как на большой глубине можно исключить влияние дна сосуда на собственные частоты. Собственная форма колебания жидкости для первой частоты при полной невесомости показана на рис. 7.

Рис. 5. Квадраты собственных частот Ж2 при разных значениях числа Бонда Во и глубине жидкости Н = 2

Ж2

0 0,2 0,4 0,6 0,8 1,0 1,2 1,4 1,6 1,8 н

Рис. 6. Квадраты собственных частот Ж2 в зависимости от глубины жидкости Н при разных значениях числа Бонда и 5 = 0,5

z

Рис. 7. Собственная форма колебания жидкости при Bo = 0, а0 = 10°, H = 2

Заключение. В настоящей работе численно определены формы равновесия, собственные частоты и формы колебаний свободной поверхности жидкости в условиях микрогравитации при любом соотношении радиусов стенок сосуда, имеющего форму коаксиального цилиндра. Для случая, когда допустимо отклонение свободной поверхности вдоль оси Oz, подробно приведен вывод кинематического и динамического условий на возмущенной свободной поверхности с учетом влияния силы поверхностного натяжения. Из результатов исследований следует, что угол смачивания а0 влияет на форму равновесной свободной поверхности и, следовательно, оказывает влияние на динамические характеристики колебания жидкости. При Bo > 5 и 5 > 0,2 форма равновесия свободной поверхности становится близка к плоской, и можно проводить оценку собственных частот колебаний по аналитическими формулами, получаемыми при ао = 90.

ЛИТЕРАТУРА

[1] Мышкис А.Д., Бабский В.Г., Жуков М.Ю., Копачевский Н.Д., Слобожанин Л.А., Тюпцов А.Д. Методы решения задачи гидромеханики для условий невесомости. Киев, Наукова Думка, 1992, 592 с.

[2] Dodge F.T. The new "Dynamic behavior of liquids in moving containers". NASA SP-106, 2000, 202 p.

[3] Ibrahim R.A. Liquid Sloshing Dynamics: Theory and Applications. 2nd ed. Cambridge, Cambridge Univ. Press, 2005, 948 p.

[4] Полевиков В.К. О методах численного моделирования равновесных капиллярных поверхностей. Дифференциальные уравнения, 1999, т. 35, № 7, с. 975-981.

[5] Cheng X., Wang Z. The equation and the numerical analysis of static fluid surface in revolving symmetrical tank under low gravity. Chinese Journal of Computational Physics, 2000, vol. 17 (3), pp. 273-279.

[6] Yang D., Yue B., Zhu L., Song X. Solving shapes of hydrostatic surface in rectangular and revolving symmetrical tanks under microgravity using shooting method. Chinese Journal of Space Science, 2012, vol. 32 (1), pp. 85-91.

[7] Chu W. Low-Gravity Fuel Sloshing in an Arbitrary Axisymmetric Rigid Tank. ASME. J. Appl. Mech., 1970, vol. 37 (3), pp. 828-837. https://doi.org/10.1115/L3408616

[8] Wang Z., Deng Z. Sloshing of Liquid in Spherical Tank at Low-gravity Environments. Chinese Journal of Space Science, 1985, vol. 5 (4), pp. 294-302.

[9] Wang Z., Deng Z. On the Sloshing of Liquid in a Partially Filled Rectangular Tank under Low-gravity Condition. Journal of Tsinghua University, 1986, vol. 26 (3), pp. 1-9.

[10] Utsumi M. Low-gravity propellant slosh analysis using spherical coordinates. Journal of Fluids and Structures, 1998, vol. 12 (1), pp. 57-83. https://doi.org/10.1006/jfls.1997.0125

[11] Siekmann J., Scheideler W., Tietze P. Static meniscus configurations in propellant tanks under reduced gravity. Computer Methods in Applied Mechanics and Engineering, 1981, vol. 28 (1), pp. 103-116.

[12] Zhou D. The Analysis of The Interface Shape of Rotating Liquid Annular Container at Zero-gravity Environment. Structure & Environment Engineering, 1994, no. 2, pp. 40-44.

[13] Dodge F.T., Garza L.R. Experimental and Theoretical Studies of Liquid Sloshing at Simulated Low Gravity. ASME. J. Appl. Mech., 1967, vol. 34 (3), pp. 555-562. https://doi.org/10.1115/1.3607743

[14] Wu W., Yue B. An Analytical Method for Studying the Sloshing Properties of Liquid in Cylindrical Tank under Low Gravity Environment. Journal of Astronautics, 2014, vol. 35 (4), pp. 397-403.

[15] Yue B., Wu W., Yan Y. Modeling and Coupling Dynamics of the Spacecraft with Multiple Propellant Tanks. AIAA J., 2016, vol. 54 (11), pp. 3608-3618. https://doi.org/10.2514/1J055110

[16] Yue B., Yu J., Wu W. Rigid and liquid coupling dynamics and hybrid control of spacecraft with multiple propellant tanks. Chinese Journal of Theoretical and Applied Mechanics, 2017, vol. 49 (2), pp. 390-396. URL: http://lxxb.cstam.org.cn/CN/10.6052/0459-1879-16-342 (дата обращения 25.04.2021).

[17] Юй Чжаокай, Темнов А.Н. Исследование равновесной свободной поверхности капиллярной жидкости в тороидальном сосуде. Инженерный журнал: наука и инновации, 2021, вып. 3. http://dx.doi.org/10.18698/2308-6033-2021-3-2060

[18] Моисеев Н.Н., Петров А.А. Численные методы расчета собственных частот колебаний ограниченного объема жидкости. Москва, Вычислительный центр АН СССР, 1966, 272 с.

[19] Labus T. Natural frequency of liquids in annular cylinders under low gravitational conditions. NASA TN D-5412, 1969, 28 p.

Статья поступила в редакцию 06.07.2021

Ссылку на эту статью просим оформлять следующим образом:

Юй Чжаокай, Темнов А.Н. Равновесие и колебания свободной поверхности жидкого топлива в коаксиально-цилиндрических сосудах в условиях микрогравитации. Инженерный журнал: наука и инновации, 2021, вып. 8.

http://dx.doi.org/10.18698/2308-6033-2021-8-2099

Юй Чжаокай — аспирант кафедры «Космические аппараты и ракеты-носители» МГТУ им. Н.Э. Баумана. e-mail: yuzhaokai933@mail.ru

Темнов Александр Николаевич — канд. физ.-мат. наук, доцент кафедры «Космические аппараты и ракеты-носители» МГТУ им. Н.Э. Баумана. e-mail: antt45@mail.ru

Equilibrium and oscillations of liquid fuel free surface in coaxial-cylindrical vessels under microgravity conditions

© Yu Zhaokai, A.N. Temnov Bauman Moscow State Technical University, Moscow, 105005, Russia

In the absence of significant mass forces, the behavior of liquid fuel under microgravity conditions is determined by surface tension forces, which are intermolecular forces at the interface of two phases. The paper posed and solved the problem of equilibrium and small oscillations of an ideal liquid under microgravity conditions, and also quantified the influence of various parameters: the contact angle a0, the Bond number, the ratio of the radii of the inner and outer walls of the vessel and the depth of the liquid. For the coaxial-cylindrical vessels, there were obtained expressions in the form of a Bessel series for the potential of the fluid velocities and the free surface displacement field. The study relies on the analytical and experimental data available in the literature and proves the reliability of the developed numerical algorithm. Findings of research show that for Bo > 5 and 5 > 0,2 r, with the physical state of the wetted surface being unchanged, the shape of the free surface tends to be flat and the contact angle has little effect on the intrinsic vibration frequency of the free surface of the liquid. The results obtained can be used to solve problems of determining the hydrodynamic characteristics of the movement of liquid fuel in outer space.

Keywords: microgravity, surface tension force, capillary fluid, coaxial cylinder, Runge — Kutta method, Bessel series

REFERENCES

[1] Myshkis A.D., Babskiy V.G., Zhukov M.Yu., Kopachevskiy N.D., Slobozha-nin L.A., Tyuptsov A.D. Metody resheniya zadachi gidromekhaniki dlya usloviy nevesomosti [Methods for solving the problems of hydromechanics for weightlessness conditions]. Kiev, Naukova Dumka Publ., 1992, 592 p.

[2] Dodge F.T. The new "Dynamic behavior of liquids in moving containers". NASA SP-106, 2000, 202 p.

[3] Ibrahim R.A. Liquid Sloshing Dynamics: Theory and Applications. 2nd ed., Cambridge, Cambridge Univ. Press, 2005, 948 p.

[4] Polevikov V.K. Differentsialnye uravneniya — Differential Equations, 1999, vol. 35, no. 7, pp. 975-981.

[5] Cheng X., Wang Z. The equation and the numerical analysis of static fluid surface in revolving symmetrical tank under low gravity. Chinese Journal of Computational Physics, 2000, vol. 17 (3), pp. 273-279.

[6] Yang D., Yue B., Zhu L., Song X. Solving shapes of hydrostatic surface in rectangular and revolving symmetrical tanks under microgravity using shooting method. Chinese Journal of Space Science, 2012, vol. 32 (1), pp. 85-91.

[7] Chu W. Low-Gravity Fuel Sloshing in an Arbitrary Axisymmetric Rigid Tank. ASME. J. Appl. Mech, 1970, vol. 37 (3), pp. 828-837. https://doi.org/10.1115/L3408616

[8] Wang Z., Deng Z. Sloshing of Liquid in Spherical Tank at Low-gravity Environments. Chinese Journal of Space Science, 1985, vol. 5 (4), pp. 294-302.

[9] Wang Z., Deng Z. On the sloshing of liquid in a partially filled rectangular tank under low-gravity condition. Journal of Tsinghua University, 1986, vol. 26 (3), pp. 1-9.

[10] Utsumi M. Low-gravity propellant slosh analysis using spherical coordinates. Journal of Fluids and Structures, 1998, vol. 12 (1), pp. 57-83. https://doi.org/10.1006/jfls.1997.0125

[11] Siekmann J., Scheideler W., Tietze P. Static meniscus configurations in propellant tanks under reduced gravity. Computer Methods in Applied Mechanics and Engineering, 1981, vol. 28 (1), pp. 103-116.

[12] Zhou D. The analysis of the interface shape of rotating liquid annular container at zero-gravity environment. Structure & Environment Engineering, 1994, no. 2, pp. 40-44.

[13] Dodge F.T., Garza L.R. Experimental and Theoretical Studies of Liquid Sloshing at Simulated Low Gravity. ASME. J. Appl. Mech., 1967, vol. 34 (3), pp. 555-562. https://doi.org/10.1115/1.3607743

[14] Wu W., Yue B. An analytical method for studying the sloshing properties of liquid in cylindrical tank under low gravity environment. Journal of Astronautics, 2014, vol. 35 (4), pp. 397-403.

[15] Yue B., Wu W., Yan Y. Modeling and coupling dynamics of the spacecraft with multiple propellant tanks. AIAA J., 2016, vol. 54 (11), pp. 3608-3618. https://doi.org/10.2514/1J055110

[16] Yue B., Yu J., Wu W. Rigid and liquid coupling dynamics and hybrid control of spacecraft with multiple propellant tanks. Chinese Journal of Theoretical and Applied Mechanics, 2017, vol. 49 (2), pp. 390-396. Available at: http://lxxb.cstam.org.cn/CN/10.6052/0459-1879-16-342 (accessed April 25, 2021).

[17] Yu Zhaokai, Temnov A.N. Inzhenerny zhurnal: nauka i innovatsii — Engineering Journal: Science and Innovation, 2021, iss. 3. http://dx.doi.org/10.18698/2308-6033-2021-3-2060

[18] Moiseev N.N., Petrov A. A. Chislennye metody rascheta sobstvennykh chastot kolebaniy ogranichennogo obema zhidkosti [Numerical methods for calculating the natural frequencies of vibrations of a limited volume of liquid]. Moscow, Computing Center of USSR Academy of Sciences Publ., 1966, 272 p.

[19] Labus T. Natural frequency of liquids in annular cylinders under low gravitational conditions. NASA TN D-5412, 1969, 28 p.

Yu Zhaokai, post-graduate student, Department of Spacecraft and Launch Vehicles, Bauman Moscow State Technical University. e-mail: yuzhaokai933@mail.ru

Temnov A.N., Cand. Sc. (Phys.-Math.), Assoc. Professor, Department of Spacecraft and Launch Vehicles, Bauman Moscow State Technical University. e-mail: antt45@mail.ru