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

ABOUT INFLUENCE OF ELASTICITY OF TUBE WALLS TO DYNAMICS OF NONLINEAR WAVES OF BUBBLY LIQUID

Bayazitova Y.R. 1
1 Sterlitamak branch of Bashkir State University
1004 KB
The mathematical model to study the evolution of nonlinear waves in a tube with elastic walls, filled with bubble liquid. Numerical calculations of the wave parameters for two models of the tube with elastic walls: in the first case it is assumed that Hooke´s law is valid for the entire pipe, in the second case considered radial inertia of the pipe walls. A comparative analysis of the results of numerical experiments. It is established that the processes typical for the bubbles in the propagation of waves in two-phase liquid affect the dynamics of waves at the bubble radius , with decreasing radius bubbles up dissipation on the bubbles and interfacial interaction negligible, radial inertia of the walls of the tube retains its impact on the dynamics of the wave. A similar trend is observed in the changes of the concentration of bubbles single radius: at low concentrations of bubbles on the evolution of pressure waves have a greater impact properties of elasticity of the walls of the tubes with increasing concentrations of bubble dynamics of waves is determined by the processes occurring in bubbly liquids. The speed of wave two-phase fluid less than the speed of sound, and it decreases with increasing concentration of bubbles. This trend has continued in the accounting for the radial inertia of the walls of the tube. It should be noted that the speed of waves in this case is close to the values of this parameter typical case no radial inertia of the elastic walls of the tube.
bubbly liquid
the elasticity
the wave of pressure
the speed
the amplitude
1. Nigmatulin R.I. Dinamika mnogofaznyh sred. M.: Nauka.T1, 2, 1987. 360 р.
2. Nigmatulin R.I. Projavlenie szhimaemosti nesushej fazy pri rasprostranenii volny v puzyr’kovoj srede / Nigmatulin R.I., Shagapov V.Sh., Vahitova N.K. // Dokl. RAN, 1989. T. 304, no. 5. pp. 1077–1088.
3. Bayazitova Ya.R. Features of the dynamics of the detonation waves in layered bubble systems/ Bayazitova Ya.R., Gimaltdinov I.K., Bayazitova A.R. // Tyumen State University Herald, 2013. no. 7. рр. 25–32.
4. Volobuev A.N., Techenie zhidkosti v trubkah s jelastichnymi stenkami // UFN, 1995. T. 165, no. 2. рр. 177–186, dostupno: http://mi.mathnet.ru/ufn1052.
5. Samarskij A.A., Popov Ju.P. Raznostnye shemy gazovoj dinamiki. M.: Nauka, 1975. 352.

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

Целью нашей работы является развитие теории волновых течений пузырьковых жидкостей в трубках с эластичными стенками.

Постановка задачи и система уравнений. Рассмотрим одномерные волновые возмущения, распространяющиеся в трубке малого радиуса, заполненной пузырьковой жидкостью.

Стенки трубки считаем эластичными. На рис. 1 представлено схематическое изображение системы, которое иллюстрирует трубку толщиной h, длиной L и радиусом a0 (L >> a0). Возмущения в системе возникают вследствие воздействия давлением по торцу трубки.

pic_54.wmf

Рис. 1. Схема задачи

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

bayazit001.wmf

bayazit002.wmf

bayazit003.wmf

bayazit004.wmf

bayazit005.wmf bayazit006.wmf

bayazit007.wmf bayazit008.wmf

здесь и далее ρi – плотность фазы; n – число пузырьков; S – площадь поперечного сечения трубки; υ – скорость; αi – объемное содержание фазы; а – радиус пузырьков; γ – показатель адиабаты для газа; q – тепловой поток; pg – давление газа.

Будем считать жидкость акустически сжимаемой, газ калорически совершенным:

bayazit009.wmf bayazit010.wmf

здесь bayazit011.wmf – начальное давление жидкости; R – газовая постоянная; Tg – температура газа.

При описании радиального движения в соответствии с уточнением, предложенным в [2], будем полагать, что w = wR + wA, где wR – определяется из уравнения Релея – Ламба, а wA – определяется из решения задачи о сферической разгрузке на сфере радиуса a в несущей жидкости в акустическом приближении:

bayazit012.wmf

bayazit013.wmf

где vl – вязкость жидкости; Cl – скорость звука в «чистой» жидкости.

Тепловой поток задаем следующим приближенным конечным соотношением [3]:

bayazit014.wmf

bayazit015.wmf

где bayazit016.wmf bayazit017.wmf

Температура газа в пузырьках меняется по закону

bayazit018.wmf.

Рассмотрим два случая: будем считать, что стенки трубки помимо эластичных свойств обладают массой, обусловливающей радиальную инерцию стенок при распространении волновых возмущений, и исследуем динамику нелинейных волн в эластичной трубке без инерционных свойств ее стенок. Закон Гука для эластичной трубки имеет вид [4]

bayazit019.wmf

В случае учета радиальной инерции стенок трубки считаем, что радиус трубки меняется согласно уравнению

bayazit020.wmf,

где E – модуль Юнга; h0 – толщина стенки трубки; b0 – начальное значение внутреннего радиуса трубки; ρw – плотность материала трубки.

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

bayazit021.wmf

bayazit022.wmf

bayazit023.wmf

bayazit024.wmf

bayazit025.wmf

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

bayazit026.wmf

bayazit027.wmf

bayazit028.wmf

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

bayazit029.wmf

bayazit030.wmf

Системы уравнений решаются численно по явной схеме с использованием метода конечных разностей [5].

Начальные и граничные условия. Нелинейные волны, возникающие в трубке, распространяются вдоль положительного направления оси z до противоположного торца трубки, ограниченного жесткой стенкой: z = L: υz = 0, 0 < a0 < R.

Начальные условия для остальных входных параметров соответствуют следующим:

t = 0: pl = p0, S = S0, pg = p0, wR = 0, a0 = 1 мм или 0,1 мм.

Возмущения в трубе возникают вследствие воздействия по торцу трубки давлением согласно условиям при t* = 1 мс и Δp = 1 атм:

bayazit031.wmf (1)

bayazit032.wmf (2)

Обсуждение результатов

pic_55.tif

Рис 2. Зависимость скорости нелинейных волн в пузырьковой жидкости от концентрации пузырьков

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

На рис. 3 представлены эпюры давления в жидкости в момент времени 1,5 мс, графики получены при различных радиусах пузырьков: на рис. 3, а эпюры построены при a0 = 1 мм, на рис. 3, б – при a0 = 0,1 мм приближаются к значению скорости для волн в трубке с жесткими стенками. Сплошная линия на рис. 3 соответствует случаю распространения нелинейных волн в трубке с поликарбонатными стенками с учетом радиальной инерции стенок трубки, пунктирная линия – случаю распространения волн в трубке с жесткими стенками.

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

pic_56.wmf аpic_57.wmf б

Рис. 3. Зависимость давления двухфазной жидкости от координаты в момент времени 1,5 мс при различных радиусах пузырьков: а – случай соответствует радиусу пузырьков 1 мм; б – случай – радиусу пузырьков 0,1 мм

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

На рис. 4 представлены графики изменения давления в жидкости от времени в фиксированной точке, находящейся в 500 см от торца трубки, по которому производилось первоначальное воздействие. Графики получены при различных значениях концентраций пузырьков: на рис. 4, а графики соответствуют α = 0,001, на рис. 4,б – α = 0,01. Пунктирная линия на рис. 4 соответствует случаю распространения нелинейных волн в трубке с поликарбонатными стенками с учетом радиальной инерции стенок трубки, сплошная линия – случаю распространения волн в трубке с эластичными стенками без учета инерции стенки. Из рис. 4 видно, что амплитуды волн давления при концентрации пузырьков α = 0,001 различаются на величину 0,1 атм в случаях распространения волн по трубке с эластичными поликарбонатными стенками с учетом и без учета радиальной инерции стенок, что обусловлено различными способами описания изменения поперечных размеров трубки.

pic_58.wmf аpic_59.wmf б

Рис. 4. Зависимость давления двухфазной жидкости от времени в точке, удаленной от нижнего торца трубы на расстояние z = 500 см, при различных объемных содержаниях пузырьков: а – случай соответствует концентрации пузырьков 0,001; б – случай – концентрации пузырьков 0,01

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

Заключение

По результатам исследования сделали следующие выводы: процессы, характерные для пузырьков при распространении волн двухфазной жидкости, оказывают влияние на динамику волн при радиусе пузырьков r = 1 мм, при уменьшении радиуса пузырьков до r = 0,1 мм диссипация на пузырьках и межфазное взаимодействие пренебрежимо малы, радиальная инерция стенок трубки сохраняет свое воздействие на динамику волны. Аналогичная тенденция наблюдается при изменении концентрации пузырьков единого радиуса: при малых концентрациях пузырьков на эволюцию волн давления оказывают большее влияние свойства эластичности стенок трубки, при увеличении концентрации пузырьков динамика волн определяется процессами, протекающими в пузырьковой жидкости. Скорость волн двухфазной жидкости меньше скорости звука, и она убывает с увеличением концентрации пузырьков. Эта тенденция сохраняется и при учете радиальной инерции стенок трубки. Следует отметить, что скорость волн в этом случае близка к значениям данного параметра, характерного для случая отсутствия радиальной инерции у эластичных стенок трубки.

Рецензенты:

Гималтдинов И.К., д.ф.-м.н., профессор, зав. кафедрой «Прикладная математика и механика» Стерлитамакского филиала БашГУ, г. Стерлитамак;

Михайлов П.Н., д.ф.-м.н., профессор, зав. кафедрой «Алгебра, геометрия и методика преподавания математики» Стерлитамакского филиала БашГУ, г. Стерлитамак.

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