CC BY
84
^ВЕСТНИК
......рана ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
ХиуЮРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА
УДК: 532.59; 532.595 Г. А. Ахундова,
ассист., докторант
ИССЛЕДОВАНИЕ ПОГЛОЩЕНИЯ ИНТЕНСИВНОСТИ АКУСТИЧЕСКОЙ ВОЛНЫ В ОКЕАНЕ
THE ABSORPTION INTENSITY OF THE ACOUSTIC WAVES IN THE OCEAN
В статье исследуется влияние пузырьков газа на поглощение звуковой волны в жидкости. В отличие от ранее опубликованных работ, используется модель баротропной пузырьковой жидкости с учетом капиллярного эффекта. Описываются основные закономерности распространения и затухания звука в морской воде. Дается физическая интерпретация этих закономерностей. В работе исследуются два типа решения: первое - синусоидальное по координате и затухающее экспоненциально во времени возмущение, второе - синусоидальное по времени колебание, затухающее экспоненциально по длине возмущения, являющееся волновым аналогом вынужденных колебаний. Показано, что коэффициент поглощения, согласно классической теории затухания, возрастает с увеличением вязкости жидкой фазы. Результаты численных расчетов показали, что затухание звука в газожидкостной среде усугубляется тем, что поглощение наиболее сильно происходит именно на резонансной частоте.
In this article we investigate the influence of gas bubbles on the absorption of sound wave in the liquid. In contrast to published studies using barotropic model bubbly liquid with the capillary effects. Describes the basic laws ofpropagation and damping of sound in seawater. Physical interpretation of these laws are given. We study two types of solutions. The first decision - sinusoidal coordinate and decay exponentially in time perturbation. A second solution - sinusoidal oscillation time, decaying exponentially with the length of the disturbance is a wave analogue of forced oscillations. It is shown that the absorption coefficient according to the classical theory of damping increases with the viscosity of the liquid phase. The numerical results showed that the damping of sound in the gas-liquid medium is compounded by the fact that the absorption occurs is most strongly at the resonance frequency.
Ключевые слова: газ, пузырек, частота, звуковая волна, колебание, поглощение звука, плотность, поле давление, затухание колебаний, объемная концентрация.
Key words: gas, bubble, frequency, sound wave, oscillation, sound absorption, density, pressure field, fluctuation damping, volume concentration.
Введение
Изучение Мирового океана с каждым годом становится все более актуальной задачей как в научном, так и в техническом отношении. Каждый год возникают все новые проблемы и становится ясным, что многие, казалось бы, давно изученные вопросы, оказываются значительно более сложными, чем предполагалось, и все они еще далеки от своего решения [1] - [3].
В настоящее время широко используются дистанционные методы экспериментального исследования Мирового океана с помощью радиоволн. Последние практически совсем не проникают в морскую воду вследствие ее электропроводности и затухают на расстоянии, много меньшем
™ длины волны. Освещенность видимой части светового спектра падает в морской воде в тысячи раз о
Ц? на протяжении десятков или, в лучшем случае, одной-двух сотен метров и лишь звуковые волны могут распространяться в водной среде на сотни километров. Это несомненное преимущество акустических волн над другими видами излучения определило их широкое использование в приборах, предназначенных для работы в водах морей и океанов.
Распространение акустических колебаний в море также сопровождается уменьшением их интенсивности. Причины затухания звуковых волн в море не исчерпываются лишь вязкостью, теплопроводностью, релаксационными (молекулярными) процессами, т. е. поглощением. В такой неоднородной среде, как морская вода, важной причиной затухания является рассеяние волны
ВЕСТНИК«
ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА ШЦ^
МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА^
неоднородностями, рассредоточенными в толще воды. Особая роль здесь принадлежит таким рас-сеивателям, как газовые пузырьки.
Рассеяние звуковой волны неоднородностями среды зависит от формы и размеров неодно-родностей, от сжимаемости и плотности их вещества, и, соответственно, от волнового сопротивления. Если сжимаемость и плотность такие же, как у среды, то неоднородности не вызывают рассеяния, независимо от их размеров и формы. В противном случае, рассеяние, в той или иной мере, будет наблюдаться [4].
Газовый пузырек - это неоднородность, отличающаяся от воды как плотностью, так и сжимаемостью. Плотность жидких тел очень мало изменяется при изменении давления - внутренние силы стойко сопротивляются внешнему воздействию, стремящемуся уменьшить или увеличить объем, занимаемый данной массой вещества. Это сопротивление столь велико, что для решения многих задач жидкости можно считать несжимаемыми. Однако даже минимальные изменения плотности твердых и жидких сред, происходящие при изменении давления, в акустике играют первостепенную роль и ими нельзя пренебрегать.
Основная часть
Исследуем влияние пузырьков газа на поглощение звуковой волны в жидкости. При этом, в отличие от опубликованных работ, допустим, что состояния пузырьковой жидкости с учетом капиллярных эффектов описывается уравнением, предложенным в работах [5] и [6]. Запишем уравнения неразрывности, одномерного движения и состояния жидкости в виде:
#+^=0 (1) + = <2>
• £
р _ а 20( + £)р/р 0 _ £ •а 2оР/Р 0 . 5 = 2а <3>
Ро 1 _а юР/Р 0 V _а юР/Р о' КоР о'
Волны можно рассматривать как возмущения малой амплитуды, учет которых позволяет проводить линеаризацию уравнений (1) - (3):
др ди
Т+ р° ди = °' (4)
р° ^ Л = °' (5)
д1 дх
Р = С, С0 = . ^ , (6)
р V" 20Р0 3
где р, р, и - возмущения соответствующих характеристик акустической волной. Из уравнений (4) и (6) следует
1 др ди
—г — + р0 — = 0. (/)
С0 дt дх
Если это уравнение продифференцировать по а уравнение (5) - по х, то удастся исключить возмущение скорости и получить волновое уравнение
д2 Р _ д2 Р
ВО
Г
сс
X
ы
д?2 ° дх2
<0 £
В водной среде акустические колебания являются продольными. Рассмотрим наиболее простой, представляющий практический интерес случай плоских гармонических продольных волн в виде действительной части комплексного выражения:
p = p exp[i(k*x + Ш)] ; p0 = p0 + ipl,;
— 0 p
—0
<< 1; k* = k + ik**; ю, = ю + /ю,
(9)
где p - определитель амплитуды возмущении давления. Возможны два типа решения.
Тип I. Решение или волна первого типа, когда к* = k - действительное положительное число
(k > 0, к** = 0).
В этом случае будем иметь:
p = (p! + ip 1) exp [i (kx + Ш + i®„t)] = (p! + ip 1) exp (-m,,t) x |^cos (kx + roí) + i sin (kx + roí)];
Re {p} = exp (-w*J) p sin [ф + (kx + wt)] .
Здесь
— 0 p
V—02 —0 2 / -Г-0.
p* + p** ; ф = arctg(-p*/p,,).
Таким образом, решение первого типа представляет собой синусоидальное по координате и при ю** > 0 экспоненциально затухающее во времени возмущение, которое называют k-волной:
_/ \ h0 г / \ п I 2пГх + Уф(&)Л|
Р (k) = Р exp |_-Ю„ (k)t\ sin <|ф +-1—Ц£)-=4 ,
где уф(к) = ю(к)/к; Цк) = 2п/к; ф - начальная фаза.
Здесь уф( к) - фазовая скорость, или скорость перемещения фазы колебания; Цк) - длина волны; ю**(к) - декремент затухания колебаний во времени.
Другими словами, к-волны - это волны с однородной по длине, но изменяющейся во времени амплитудой, являющиеся аналогом свободных колебаний.
Тип II. Решения или волны второго типа, когда ю* = ю - действительное положительное число, т. е. ю > 0 и ю** = 0. В данном случае имеем:
Р = (Р* + ip**) ехР [ i (kx + + ik**x)] = (p* + ip**) exp (-k**x) x [cos (kx + roí) + i sin (kx + roí)];
Re jpj = exp (-k**x) p sin [ф + (kx + rot)].
Таким образом, решение второго типа представляет собой синусоидальное по времени колебание, возбуждаемое, например, каким-либо стационарным внешним источником монохроматических колебаний при x = 0, экспоненциально затухающей по длине амплитудой. Такие возмущения, являющиеся волновым аналогом вынужденных колебаний, называют ю-волнами:
p(ю) = p° (ю) exp(-£,, (co)x)sm 1ф + 2п[x + v0(a)t] 1 ; ^ (ю) = ю/^(ю); х(ю) = 2п /к(ю). v / I Х(ю) I
Здесь к** (ю) - декремент затухания колебаний по длине.
Другими словами, ю-волны - это волны со стационарными по времени, но изменяющимися по длине амплитудами. Случай к** > 0 соответствует режиму затухания амплитуд возмущения. Ограничимся только ю-волнами (ю > 0, ю** = 0).
Если воспользоваться выражением (9) в уравнении (8), то получим дисперсионное уравнение вида ю/к = C Поскольку поток количества движения за счет трения и других внешних сил в уравнении (2) не учитывался, то волновое число оказывается действительным положительным числом и фазовая скорость оказывается равной скорости звука: ю/к = C а декремент затухания равен нулю.
Рассмотрим иную постановку задачу с эффективной вязкостью в уравнении радиального движения для учета всех возможных диссипативных эффектов с использованием условия поли-тропии газа [7] - [10]:
3 Я _ Р 2 - Р - 2а / Я - 4ц е/Я , р2 =
Р? Р 20
ЯЯ+2Я _
0 р.
р о
20 У
(10)
Уравнение импульса среды, сохранения массы фаз и индивидуального пузырька можно представить в следующем виде:
ди ди др о о
р — + ри — + = 0' р » р°а1' р° = const' дt дх дх
Ф + Фц) = °' ^ + = 0' 1 (р2*з ) = 0.
Я* ^ Я/ ^ Л"2 >
(11)
Объемная концентрация газа а20 связана с количеством пузырьков в единице объема смеси равенством:
а = — пЯ п
(12)
п р где — = —.
п° Р°
Если выполнить линеаризацию уравнений (10) - (12), то в безразмерных переменных и параметрах запишем:
—о
В. = Я / ^о' р2 =р2/Р20' р = р / Ро' р 2 = Рг! Рго' и = и / ио'
Щ = R0/г, = Я02/V ^ V е/ = / р0о' 5 = 2а / R0p0' Ре = '
Рю V2
(13)
Рю = Ро + 2а / R0' т = г / Го' X = х / Л0.
Получим следующую систему уравнений относительно возмущений переменных:
5а,
а _ -а 10 дУ' а10 д^ _ Ре д Р ' а 2 " 3а 10« 20я , Р 2 2 '
— — —0
# + 3 § = о' 0 = Ре{(1 + 8)~Р 2 - р + й]
_ 4 -я
d т '
(14)
, _ = о• -2'
дт дт - т2
Исключив возмущения скорости, объемного газосодержания, давления и плотности газа, получим систему уравнений:
ё2Я п д2р d2Д „ г^Го ом -) ^Л
3а1Па,п^- = Ре •—^' —- = Ре -{Л [Я - 3£ (1 + 5 )- р}- 4 —
40^20 2
дт
дХ2 dт2
(15)
Будем искать решение системы (15) в виде действительной части комплексного выражения:
— —0 г , __ __ _ ч — 0 —0 —0
р = р ехр ^ (К*Х + Пт)]' р = р* + ip„'
—0 р
<< 1' к* = к + гк,,'
R = R0 ехр [г (К*X + Пт)]' R0 = R0 + ^ «'
— 0 Я
(16)
<< 1' ft = юR02/ V ^' к * = к ,Rn.
Получим систему линейных однородных алгебраических уравнений:
3а10а20П2 = Ре• р° • К*' Я° -П2 = Ре-{{ [5-3к(1 + 5)] + + 4Я°1 -П. (17)
Из условия существования нетривиального решения системы (17) вытекает следующая дисперсионная зависимость волнового числа К от частоты
во г
сс
X Ы
2
к
к _ За10а20 _. Q = + ^ ре.
(18)
^ + Q - 4/ ~ 4
Решив уравнение (18) с помощью выделения действительной и мнимой части решения, най-
— V .
дем выражения для фазовой скорости С ^ = С ^ ■ и декремента затухания колебаний по длине
"о
к„ = К„ / в виде:
т; О
с Г =-¡7 =
Qo
к
Оо =
+ Q)2 +16 ^/(0+0)476
к ** = ^0
3а10а 20
^ + Q
(19)
(20)
нате.
1(п + Q)2 +16 (П + Q) +16
Мнимая часть К** > 0 положительна. Поэтому амплитуды возмущения затухают по коорди-
На рис. 1 и 2 приведены зависимости коэффициента поглощения к„ = К**/ R0 и фазовой скорости Су = Су -\еГ / R0 от частоты ю = Ом^ / Я^. Равновесный радиус пузырьков принимался равным 10 мкм, а объемное газосодержание - 5 %.
Анализ выражения (20) показывает, что коэффициент поглощения, согласно классической теории затухания, возрастает с увеличением коэффициента вязкости и, что очень важно, частоты.
Как видно из рис. 1, характерным является диапазон частот, который иногда называют диапазоном непрозрачности из-за больших значений декремента затухания к**. Это означает, что в данном диапазоне волна будет сильнее поглощаться.
О.Е+ОО 2.Е+10 4.Е+10 6.Е+10 Частота, 1/с
Рис. 1. Зависимость декремента затухания от частоты
го £ и
Е? ей
.Е+10
Частота,
Рис. 2. Зависимость фазовой скорости от частоты
2
1
ВЕСТНИК^Г
ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА ШЧГ^.
МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА
В связи с тем, что в морской воде в течение длительного периода времени могут существовать во взвешенном состоянии пузырьки с миллимертовыми диаметрами, аномальное затухание в море могут испытывать волны, имеющие частоту, равную десяткам килогерц, что подтверждается также практическими наблюдениями.
Затухание звука на газовых пузырьках усугубляется еще и тем, что они не только рассеивают, но и поглощают часть энергии звуковой волны, и это поглощение наиболее сильно происходит именно на резонансной частоте.
Выводы
Отметим, что в приведенной постановке затухание колебаний пузырьков возможно только за счет вязкой диссипации и сжимаемости газожидкостной смеси, так как при выводе уравнения состояния пузырьковой жидкости считалось, что температура смеси постоянна, тепловая диссипация и сжимаемость несущей фазы не учитывались. В этом можно убедиться из анализа аналитического решения дисперсионного уравнения, представленого в виде формул для фазовой скорости и декремента затухания колебаний. Расчеты, выполненные на основе этого уравнения состояния, при изучении характера затухания звуковых волн в море показали, что интенсивность волны значительно убывает из-за наличия газа в жидкости, а демпфирование волны существенно возрастает с увеличением объемного газосодержания.
Влияние поверхностного натяжения на характеристики звуковой волны в жидкости с размерами пузырьков более 10 мкм несущественно и им можно пренебречь.
Список литературы
1. Безруков Ю. Ф. Океанология: в 2 ч. — Ч. I. Физические явления и процессы в океане / Ю. Ф. Безруков. — Симферополь: Таврический нац. ун-т им. В. И. Вернадского, 2006. — 159 с.
2. Доронин Ю. П. Физика океана / Ю. П. Доронин. — СПб., 2010. — 274 с.
3. Голодов Н. Ф. Методы прогнозирования температуры и солености моря для определения скорости звука в воде / Н. Ф. Голодов // Судовождение: сб. науч. тр. / Одесская национальная морская академия. — Вып. 15. — Одесса: ИздатИнформ, 2008. — С. 44-52.
4. Nigmatulin R. I. Mathimatical modelling of a single bubble and multibubble dynamics in a liquid / R. I. Nigmatulin, I. Sh. Akhatov // Proc. Int. Conf. on Multiphase Systems. — Ufa. — 2000. — Р. 294-301.
5. Нигматулин Р. И. Динамика импульсных волн в пузырьковых жидкостях. Сравнение теории с экспериментом / Р. И. Нигматулин, Д. А. Губайдуллин, А. А. Никифоров // Доклады Академии наук. — 2014. — Т. 456. — № 6. — С. 662.
6. Губайдуллин А. А. Введение в волновую динамику газожидкостных сред: учеб. пособие / А. А. Губайдуллин. — Тюмень: ТюмГНГУ, 2006. — 86 с.
7. Нагиев Ф. Б. Нелинейные колебания растворимых газовых пузырьков в жидкости / Ф. Б. Нагиев // Известия АН Аз. ССР. Серия: Физ.-техн. и матем. науки. — 1985. — № 1. — С. 136-140. Ы
8. Нагиев Ф. Б. Декременты затухания колебаний растворимых газовых пузырьков, радиально пуль- с сирующих в жидкости / Ф. Б. Нагиев // Известия АН Аз. ССР. Серия: физ.-техн. и матем. науки. — 1984. —
№ 4. — С. 125-130.
9. Хабеев Н. С. Динамика растворимых газовых пузырьков / Н. С. Хабеев, Ф. Б. Нагиев // Известия АН СССР. Механика жидкости и газа. — 1985. — № 6. — С. 52-59.
10. Nigmatulin R. Dynamics, heat-mass transfer of vapor-gas bubbles in a liquid. Printed in Great Britain / R. I. Nigmatulin, N. S. Khabeev, F. B. Nagiyev // Inter. J. Heat and Mass Transfer. — 1981. — Vol. 24. — № 6. — Р. 1033-1044.
