CC BY
97
НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. Н. Э. БАУМАНА
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 • 48211. Государственная регистрация №0421200025. ISSN 1994-0408
электронный научно-технический журнал
Динамика двухслойной жидкости, разделенной упругой
перегородкой с учетом сил поверхностного натяжения
# 11, ноябрь 2013
Б01: 10.7463/1113.0619258
Гончаров Д. А.
УДК 534-141
Россия, МГТУ им. Н.Э. Баумана zorghhh@gmail.com
Введение.
При работе маршевого двигателя разгонных блоков и иных космических аппаратов ввиду наличия длительных пассивных участков, во время которых реализуются условия микрогравитации [6], для исключения аварийных ситуаций необходимо обеспечивать бесперебойную подачу компонент в заборное устройство, для это цели широко применяются полупроницаемые капиллярные фазоразделители[4], обеспечивающие сплошность компонент при их подаче в заборное устройство. В этой связи, представляет интерес задача о выявлении тех режимов работы в условиях малой гравитации, при которых: а) давление на входе в заборное устройство является недостаточным для осуществления бесперебойной подачи компонент топлива; б) объём жидкости становится неустойчивым, когда нарушается устойчивость и односвязность всей механической системы (жидкость, совершающая малые движения). Настоящая работа служит выяснению тех условий работы, при которых система может становится неустойчивой
В работе [5] исследовалось движение идеальной, несжимаемой и нестратифицированной жидкости, совместно с упругим днищем. Схожая с [5] постановка задачи с иным подходом к разрешению дифференциального уравнения движения пластины реализовывались в работе [6]. Работа [8] посвящена исследованию движений стратифицированной жидкости совместно с упругим днищем, представляемым в виде пластины. В работе [11] схожая с нашей задача рассматривалась с применением операторных методов и исследованием свойств спектра. Тем не менее, в работе [11] не было определено критических значений собственных частот. В работе [9] рассмотрены различные технические реализации фазоразделяющих экранов. Работа [13] посвящена исследованию общих вопросов
разрешимости и устойчивости решений класса дифференциально-операторных уравнений. В работе [14] рассматривалась схожая задача без учета сил поверхностного натяжения.
В настоящей работе разработана методика расчёта собственных частот и форм колебаний механической системы, которая может рассматриваться в качестве приближенной для анализа динамики бака с фазоразделителем. В рамках рассматриваемой модели проницаемостью фазоразделителя пренебрегаем. Решена задача с условием на свободной поверхности жидкости при малой гравитации для двухслойной жидкости, разделенной упругой перегородкой, что соответствует возможным техническим реализациям расположения капиллярных фазоразделителей [4].
1. Постановка задачи
Свяжем цилиндрическую систему координат с центром фазоразделяющей мембраны.
Рис. 1. Модель сосуда с жидкостями, разделёнными мембраной
Тогда рассмотрим область 2 = и Е2, где
"1-
(1.1)
О < г < а, ( 0 < г < а,
О < г) < 271", В2: ; 0< 7] < 2п,
О < г < к^. < г < 0.
Обозначим через Г! невозмущенную свободную поверхность жидкости, а через Г2 — равновесную поверхность фазоразделяющей мембраны, —смоченные части
первой и второй жидкостями соответственно твердой стенки сосуда
(1.2)
П = Г, Л 51; у2 = Г2 Л 52
Тогда линеаризованные уравнения движения жидкости в размерном виде в осесимметричном случае примут вид
У2Ф(1) = 0в~.1,У2Ф(2) = 0в~2. (1.4)
Будем рассматривать задачу для угла смачивания в 90 градусов, в соответствии с [7, стр. 256] и [11], свободную невозмущенную поверхность полагаем плоской.
С граничными условиями:
дФ
(1]
дп
дФ
(1]
0 на 5г,
(2]
дФ
(2)
дп
дп
д п'2 дФ
0 на Б2,
(1)
д IV2 дФ
дш1
~дГ
на Г
1>
(1.5)
(1.6)
б1дФ' а1
(1)
в1 дП
— ——— ДГ№1 = 0 на Гх. а1 дп
Здесь Ф'"^ (где I = 1,2) — потенциалы скоростей первой и второй жидкостей соответственно; п — внешняя нормаль к границе области Е, заполненной жидкостью; — отклонение по нормали свободной возмущённой поверхности жидкости от невозмущённой Г1(; д12 — плотности жидкостей выше и ниже фазоразделяющей мембраны соответственно; — коэффициент поверхностного
натяжения; Д — дифференциальный оператор Лапласа; условие (1.7), записанное в
соответствии с [7, стр. 252] определяет динамическое условие на свободной поверхности, потенциал П внешних сил во время движения считаем постоянным. В случае перехода к безразмерным параметрам в соответствии с
1 1 _ / 1 _ ^ — ; '"разм* 2 ~~ ; -^разм > ^ — I Л J ^-разм* ® ~~
(да) Фразм,П = (^) празми Z = 1м
и полагая потенциал массовых сил П = Во ■ г, Во — число Бонда, характеризующее соотношение между внешними силами и силами поверхностного натяжения, Во = Я1в*12/&ь 3* ~ ускорение свободного падения, условие (1.7) в случае угла смачивания в 90 градусов в соответствии с [7, стр. 256]примет вид
Во
ЭФ
(1)
dz
— V2
ЭФ
(1)
dz
-ЛФ
(1)
= о.
z=h-.
где Я = со2 = — 12ыразм, о) — собственная частота колебаний системы. В дальнейшем будем применять условие (1.7) в виде (1.7)'.
Смещение мембраны удовлетворяет следующему уравнению движения в осесимметричной постановке [2], [3], [5]:
d2wz 1 dw2 ms 82w2
Эг;
г дт
ЭФ
№
01"
dt
z=О
(1.8)
г=0>
где — масса мембраны на единицу площади, а т — её натяжение. Мы полагаем, что мембрана находится под действием равномерного натяжения т, т.е. если провести линию по мембране в любом направлении, то сила взаимодействия между двумя частями, разделенными элементами линии, пропорциональна длине элемента и перпендикулярна его направлению; величина илы, действующей на элемент линии длины составляет величину т&е. Величину т считаем постоянной по площади мембраны.
Граничное условие для уравнения (1.8) имеют вид
^(г, 01г=а = О,
(1.9)
I
Запишем условие постоянства нижнего объёма w2rdrjdrdz = const
(1.10)
II. Осесимметричное движение жидкости
Рассмотрим применение метода Бубнова-Галёркина для нахождения собственных значений задачи (1.1) — (1.10).
Пусть состояние равновесия жидкостей устойчиво по первому приближению. Найдём собственные числа свободных осесимметричных колебаний системы из 2-х жидкостей, разделенных упругой мембраной.
Будем искать решение соответствующей задачи в виде
Ф® = £ (z) ехрОЧ.0-
(2.1)
п= 1
Применяя метод Фурье к задаче (1.4) — (1.10) в соответствии с [6] получим следующие выражения для потенциалов скоростей:
-i-
п=1
Jo (fon fl) w2sh( fon д) Bo + ffon\ V a ) ;fon q ;
Jo (.fo-n) ¡Во + (T)]Tsh(í К 0n a )- w2ch S0n~a)
Ch + 3h
eJ^X 2.2)
Xй»
Jo (fon a Jo (.fon)
ch
fon f . i ч
a
где /о —функция Бесселя 1-го рода порядка 0, fûn — те-й корень в уравнении (1.5), О, причём
= 1,8412, f02 = 5,3315, нт. д.; а* = -аразм — приведённый к безразмерному виду радиус рассматриваемого цилиндра.
Введём обобщённые координаты qm(t),m = 1,2,... и представим перемещения мемебраны в виде разложения в ряд по собственным формам колебаний мембраны без жидкости. В соответствии с [2] такое разложение примет вид
т= 1
где £оп— га_й корень функции Бесселя 1-го рода <01 = 2,40 48, (02 = 5,5201,...
Разложим /0 в РЯД по ортогональным функциям /0
г Уп 7 П 2^0-Л(О
Ц а Г}~ Ь °пт]о \аТ)'Впт~ - е V к у
(2.3)
порядка 0,
(2.4)
Удовлетворяя книметическим условиям на мембране (1.6) с учетом разложения, а также ортогональности функций Бесселя, а так же обозначая
г =
U
-J
d dr
'dh femi)
dr
U.
■/o (fo»i) ^^ = т /'(O C2-5)
Г- f r (io™ (fc 9 dr ■ zf (0),
ЭТЭТ1 Li Li-
(2.7)
dz
получим систему уравнении:
+ % Ci) + + S3Pi3 + - = o,
^1^21 ^22 + ) + S3+ 1 ■ ■ = 0,
(2.9)
+ + (.^ш^тт ^-т) + "" 0.
Составив определитель т-го порядка которой, получаем характеристическое уравнение. Важно отметить, что недиагональные члены определителя обратятся в ноль ввиду ортогональности функций Бесселя.
Ii
0 0
о
о о
tü2p*
= о
(2.10)
III. Результаты
Была определена зависимость частоты собственных колебаний механической системы от числа Бонда, характеризующего интенсивность гравитационных сил при различных значениях частоты колебаний свободной поверхности жидкости.
Далее, раскрывая полученный определитель для I-го тона, мы получим частотное уравнение относительно со. Ограничиваемся 5-ю членами ряда в уравнении относительно и>, полученном при раскрытии определителя (2.10). В полученном соотношении выражаем ш через число Бонда. Полученные зависимости при различных т, Macct мембраны: т% = 30 кг, плотности д2 = 70 размеры области
я = 1 .V, .'¡1 = l.v, .'!- = ! .'■■;.
Критическое значение частоты получаем из условия Во = 0. Из рис. 2 видим, что с уменьшением жесткости мембраны, снижается значение критической частоты.
..... 2т = 1500 Н/м, ш^р =2,26 с-1
... 3, т = 1800 Н/м, = 2,37 с-1 -- 4.т — 1900 Н/м, = 2,33 с"1
Рис. 2. Зависимость квадрата собственной частоты механической системы от числа
Бонда.
Список литературы
1. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. В 10 т. Т. 6. Гидродинамика. М.: Физматлит, 2006. 736 с.
2. Мартинсон Л.К., Малов Ю.И. Дифференциальные уравнения математической физики. М.: Изд-во МГТУ им. Н.Э. Баумана, 2006. 368 с.
3. Тихонов А.Н., Самарский А.А. Уравнения математической физики. М.: Наука, 1972. 736 с.
4. Багров В.В., Курпатенков А.В., Поляев В.М., Синцов А.Л., Сухоставец В.Ф. Капиллярные системы отбора жидкости из баков космических летательных аппаратов / под ред. В.М. Поляева. М.: УНПЦ «Энергомаш», 1997. 328 с.
5. Пожалостин А.А. Свободные колебания жидкости в жестком круговом цилиндрическом сосуде с упругим плоским дном // Известия высших учебных заведений. Сер. Авиационная техника. 1963. № 4. С. 25-32.
6. Петренко М.П. Собственные колебания жидкости со свободной поверхностью и упругого днища цилиндрической полости // Прикладная механика. 1969. Том 5, вып. 6. С. 44-50.
7. Бабский В.Г., Копаческий Н.Д., Мышкис А.Д., Слобожанин Л.А., Тюпцов А.Д. Гидромеханика невесомости / под ред. А.Д. Мышкиса. М.: Наука, 1976. 504 с.
8. Андронов А.В. Колебания идеальной стратифицированной жидкости в контейнере с упругим днищем // Вопросы волновых движений жидкости : сб. науч. тр. Краснодар: КубГУ,1987.
9. Сапожников В.Б., Меньшиков В.А., Партола И.С., Корольков А.В. Развитие идей профессора В.М. Поляева по применению пористо-сетчатых материалов для внутрибаковых устройств, обеспечивающих многократный запуск жидкостных ракетных двигателей // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. 2006. № 2. С. 78-88.
10. Михлин С.Г. Вариационные методы в математической физике. М.: Наука, 1969. 510 с.
11. Нго Зуй Кан. О движении идеальной жидкости, подверженной силам поверхностного натяжения, заполняющей сосуд с плоским днищем // Изв. АН СССР. Механика твердого тела. 1980. № 3. С. 143-154.
12. Бляшке В. Дифференциальная геометрия и геометрические основы теории относительности Эйнштейна: пер. с нем. М.: ОНТИ, 1935. 332 с.
13. Пашкова Ю.С. Малые колебания системы идеальных жидкостей, разделенных упругими перегородками / Симферопольский государственный университет. Симферополь, 1992. 28 с. Деп. в УкрИНТЭИ 02.10.92.
14. Гончаров Д.А. Осесимметричные колебания двухплотностной жидкости в цилиндрическом баке // Наука и образование. МГТУ им. Н.Э. Баумана Электрон. журн. 2012. № 4. Режим доступа: http://technomag.bmstu.ru/doc/362856.html (дата обращения 01.09.2013).
SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU
SCIENCE and EDUCATION
EL № FS77 - 48211. №0421200025. ISSN 1994-040S
electronic scientific and technical journal
Dynamics of two-layer liquid, divided by an elastic baffle considering surface tension forces (with corrections) # 11, November 2013 DOI: 10.7463/1113.0619258 Goncharov D.A.
Bauman Moscow State Technical University, 105005, Moscow, Russian Federation
zorghhh@gmail.com
Dynamics of an ideal incompressible fluid filling a cylindrical vessel divided by an elastic membrane with an allowance for surface tension forces was investigated in this work. Boundary conditions for surface tension at the wetting angle of 90 degrees were formulated; normal oscillations were determined. The Bubnov-Galerkin method was used to obtain an approximate frequency equation in an axially symmetric case (plane problem). Dependence of the natural frequency of oscillations on the Bond number which characterises the strength of gravity was analysed. A critical value of the natural oscillation frequency was also determined.
Publications with keywords: liquid, phase separating membrane, surface tension, contact angle Publications with words: liquid, phase separating membrane, surface tension, contact angle
References
1. Landau L.D., Lifshits E.M. Teoreticheskaya fizika. V10 t. T. 6. Gidrodinamika [Theoretical physics. In 10 vols. Vol. 6. Hydrodynamics]. Moscow, Fizmatlit, 2006. 736 p.
2. Martinson L.K., Malov Iu.I. Differentsial'nye uravneniia matematicheskoi fiziki [Differential equations of mathematical physics]. Moscow, Bauman MSTU Publ., 2006. 368 p.
3. Tikhonov A.N., Samarskiy A.A. Uravneniya matematicheskoy fiziki [Equations of mathematical physics]. Moscow, Nauka, 1977. 736 p.
4. Bagrov V.V., Kurpatenkov A.V., Poliaev V.M., et al. Kapilliarnye sistemy otbora zhidkosti iz bakov kosmicheskikh apparatov [The capillary system of selection of the liquid from the tanks of space vehicles ]. Moscow, UNPTs «Energomash» Publ., 1997. 328 p.
5. Pozhalostin A.A. Svobodnye kolebaniia zhidkosti v zhestkom krugovom tsilindricheskom sosude s uprugim ploskim dnom [Free oscillations of liquid in a rigid circular cylindrical vessel with an elastic flat-bottomed ]. Izvestiia vysshikh uchebnykh zavedenii. Ser. Aviatsionnaia tekhnika, 1963, no. 4. pp. 25-32.
6. Petrenko M.P. Sobstvennye kolebaniia zhidkosti so svobodnoi poverkhnost'iu i uprugogo dnishcha tsilindricheskoi polosti [Natural oscillations of a fluid with free surface and an elastic bottom of the cylindrical cavity]. Prikladnaia mekhanika, 1969, vol.5, no. 6, pp. 44-50.
7. Babskiy V.G., Kopacheskiy N.D., Myshkis A.D., Slobozhanin L.A., Tyuptsov A.D. Gidromekhanika nevesomosti [Hydromechanics of weightlessness]. Moscow, Nauka, 1976. 504 p.
8. Andronov A.V. Kolebaniia ideal'noi stratifitsirovannoi zhidkosti v konteinere s uprugim dnishchem [Oscillations of an ideal stratified fluid in a container with elastic bottom]. Voprosy volnovykh dvizhenii zhidkosti: sb. nauch. tr. [Questions of wave motions of fluid: collected papers]. Krasnodar, KubGU Publ., 1987.
9. Sapozhnikov V.B., Men'shikov V.A., Partola I.S., Korol'kov A.V. Razvitie idey professora V.M. Polyaeva po primeneniyu poristo-setchatykh materialov dlya vnutribakovykh ustroystv, obespechivayushchikh mnogokratnyy zapusk zhidkostnykh raketnykh dvigateley [Development of ideas of professor V.M. Polyaev on application of porous-meshed materials for internal tank devices providing repeated many times start-up of liquid propellant engines]. VestnikMGTU im. N.E. Baumana. Ser. Mashinostroenie [Herald of the Bauman MSTU. Ser. Mechanical Engineering], 2006, no. 2, pp. 78-88.
10. Mikhlin S.G. Variatsionnye metody v matematicheskoy fizike [Variational methods in mathematical physics]. Moscow, Nauka, 1969. 510 p.
11. Ngo Zuy Kan. O dvizhenii ideal'noy zhidkosti, podverzhennoy silam poverkhnostnogo natyazheniya, zapolnyayushchey sosud s ploskim dnishchem [On motion of ideal liquid which is exposed to the forces of surface tension and fills the vessel with flat bottom]. Izv. ANSSSR. Mekhanika tverdogo tela, 1980, no. 3, pp. 143-154.
12. Blaschke W. Differentsial'naya geometriya i geometricheskie osnovy teorii otnositel'nosti Eynshteyna [Differential geometry and geometric foundations of Einstein's theory of relativity]. Trans. from German. Moscow, ONTI, 1935. 332 p.
13. Pashkova Yu.S. Malye kolebaniya sistemy ideal'nykh zhidkostey, razdelennykh uprugimi peregorodkami [Small oscillations of system of ideal fluids separated by elastic walls]. Simferopol', 1992. 28 p. Dep. UkrlNTEI 02.10.92.
14. Goncharov D.A. Osesimmetrichnye kolebaniya dvukhplotnostnoy zhidkosti v tsilindricheskom bake [Axisymmetric oscillations of dual-density liquid in the cylindrical tank]. Nauka i obrazovanie MGTU im. N.E. Baumana [Science and Education of the Bauman MSTU], 2012, no. 4. Available at: http://technomag.bmstu.ru/doc/362856.html , accessed 01.09.2013.
