Процессы установления нелинейных колебаний  газового пузырька в жидкости

Процессы установления нелинейных колебаний газового пузырька в жидкости

Соседко Е.В.

Тихоокеанский океанологический институт им. В.И. Ильичева ДВО РАН

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

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

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

возбуждение субгармонической составляющей носит пороговый характер, то величина и вариации порогового значения амплитуды давления рт представляют значительный интерес, поскольку субгармонический сигнал возникает только в определенной области значений определяющих параметров рт и а>р - выше порога. Однако, начиная с работы Непайерса [1], отмечается наличие небольшой субгармонической составляющей и ниже порога [2-4], так что нарастание этой спектральной компоненты не имеет, строго говоря, порогового характера.

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

колебаний близким к удвоенному периоду внешнего поля (неустойчивость удвоения периода). В этом случае обращается в нуль один из двух показателей Ляпунова характеризующий линейную устойчивость указанных колебаний. При включении внешнего поля, передний фронт импульса, проходя через пузырек, помимо вынужденных, возбуждает и собственные колебания с периодом 2 Т. В окрестности порога затухание собственных колебаний, определяемое показателем линейной устойчивости, благодаря параметрической перекачке энергии в эту компоненту, очень мало. Оно обращается в нуль для порогового значения амплитуды накачки. В течение продолжительного времени пузырек продолжает совершать колебания на собственной частоте, не совпадающей с частотой внешнего воздействия. По этой причине время переходных процессов может превышать длительность импульса накачки, а соответствующая компонента в спектре излучения интерпретироваться как появление субгармонической компоненты ниже порога.

Для рассмотрения обсуждаемого эффекта использовано уравнение Рэлея-Плессета (РП), описывающего пульсации пузырька в поле давления Рр = ртсоз(а1^+а).

ЯЯ + 3 Я2 + р 2 Ро

1 -

' Яо ^ Я

V

+ 2^Яо Я = -Рр,

Уо

здесь Яо , Я - равновесный и текущий радиусы пузырька, Ро и р0 - равновесные значения давления и плотности жидкости, у - показатель политропы, Ь - затухание, эффективно учитывающее диссипативные процессы вязкости и теплопроводности, а также радиационные потери.

Для решения уравнения РП часто используют асимптотические методы, наиболее известным из которых является метод Боголюбова. Асимптотическое разложение решения в окрестности субгармонического резонанса \тр -2О0\ << а>р имеет следующий вид [6]:

(Я - Яо)/ Я о = а ^(/П о X + /З) + рр-^*). +

р0Яо

+еи1

здесь О0 = (3уР0/р0Я02)1/2 - собственная частота пузырька; е - безразмерный малый параметр, вводимый для обозначения порядка нелинейных членов, и1(а,З,Х), и2(а,З,Х) -члены разложения высших порядков.

Медленно меняющиеся амплитуда а и фаза и колебаний определяются из системы «укороченных» уравнений, следующей из требования отсутствия секулярных членов в

разложении. Учет в (1) нелинейных членов до второго порядка включительно приводит к следующему уравнению для 2=аехр(г))ехр[-г(ар/2- &о)1]:

2 = -82-/Л0.2-гд^т.exp(a)2*, Л0 = сор /2-О

Pks

0'

где рк5 = 480 0р0 Я^У 1 - порог субгармонического резонанса. Решение системы уравнений с постоянными коэффициентами (3) имеет вид

2 ()=)2

&

' . 2ЛО л 1 - /

V

Л -Л2

12 *

8рп

рк5

+

+ exp

1

*0)

/1 . 2Л0 л 1 + /

Л -Л2

+ /2*

дрг,

рк5

где Л 2 = 8

1 +

р1 )-(ЛО / 8 д

показатели линейной устой-

чивости. Вблизи порога рт=ръ-Ар, Лр>0, Ар<<ръ,

2 2 ръ=р1

1 + (Л0/8)2

показатель Л1 « 8 (Лрр^р2) мал, поэтому первое слагаемое в формуле (4) будет доминировать на временах г = Л- 1 ~ 8 (Лрръ /р^ ) I >> и описывать

субгармоническую компоненту излучения пузырька.

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

а (г) ^^г + ))) = а(0) ( / рь) exp \_-8t {Лрръ)р1)_ х

хЯс {exp ((р / 2)г + агсг^(2tga))] х

х0,51 exp г

Затягивание собственных колебаний пузырька, возбужденных передним фронтом импульса накачки, позволяет объяснить экспериментально наблюдаемую генерацию субгармонического сигнала ниже порога. В эксперименте [2] зондирование области, в которой с помощью электролиза создавались пузырьки с размерами (5-10)* 10-6 м, производилось ультразвуковыми импульсами с частотой заполнения 1,2 Мгц и

длительностью 10 мкс. Проводились измерения уровня сигнала на частоте 0,6 Мгц, которые показали на отсутствие явно выраженного порога генерации в диапазоне амплитуд зондирующих импульсов (0,25-1)*105 Па.

Из формул (3)-(5) можно получить спектральную плотность давления, излученного пузырьком на частоте субгармоники при воздействии частотнозаполненого импульса длительностью Т/тр:

Б

1 -^^^\-Ттр (д-Лу))) + 2(1 - ^(АПТШр))expГ-Т^

((8-\)/8) + (ДП/д)2

Коэффициент затухания пузырька д, резонансного на частоте 0,6 МГц, согласно

[106] составляет Ь ~ 4*105 (1/с). Порог субгармонического резонанса при этом будет равен ~ 0,8* 105 Па. Поскольку значение концентрации пузырьков не определялось в ходе обсуждаемых экспериментов, эта величина являлась подгоночным параметром. Качественное согласие с экспериментом являлось определенно аргументом в пользу предлагаемого объяснения, количественное же сопоставление требовало последовательного суммирования вклада отдельных пузырьков и численного решения уравнения РП, поскольку значительная величина порога ограничивает точность используемого асимптотического разложения.

Случай переходных процессов вблизи субгармонического резонанса был проанализирован в [5,6] для небольшой амплитуды внешней силы. Особенностью настоящей работы является изучение численными методами нелинейных колебаний пузырьков в широком интервале определяющих параметров - амплитуды накачки и радиусов пузырьков.

Выбор значений определяющих параметров диктуется следующими обстоятельствами: анализируется окрестность субгармонического резонанса для следующих значений частоты накачки п = 1.2 Мц и амплитуды накачки ^ = pm/pks, изменяющейся в пределах от 0.5* 105 до 2*105 Па (ръ ~ 0.8* 105 Па для ?=1.4 и добротности Q=П0/д=10). Длительность временной эволюции траекторий составляет 250 периодов внешнего поля. Причем после 150 периодов была отключена внешняя сила. Размеры газовых пузырьков (от 2.66 до 7.98 мкм) задавались в зависимости от частоты и отвечали дорезонансному, резонансному и зарезонансному возбуждению на субгармонике. Результат решения представлен на рис. 1 в виде зависимости Я от X', где Х'=трХ.

Рис. 1. Временная эволюция пульсаций стенки пузырька при следующих значениях определяющих параметров: добротность Q = 10, частотa п = 12, амплитуда накачки ии = [0.5; 1.25; 2.0]. Радиусы пузырьков Я = [2.66; 5.32; 7.98]. По оси х отложено: С = юр(, по оси у отложено: ии = р„/рки , по оси г отложено: (Я-Яо)/Яо

В результате изменения амплитуды накачки и радиусов пузырьков значительные

изменения претерпевает и форма акустического давления, переизлученного пузырьком: • • •

Р(()г = р(Я(()2 Я(() + 2Я(()Я(()2). На рис. 2 представлены 9 спектров плотности мощности переизлученного пузырьком звука при различных значениях определяющих параметров - амплитуды накачки и радиусов пузырьков. Исследования проводились при различных величинах добротности Q, в настоящей работе Q была положена равной 10. При малой добротности ^ < 5) характер переходных процессов имеет похожий вид, существует лишь единственное отличие - колебания устанавливаются быстрее.

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

На рис. 3 показана зависимость отклика сигнала на субгармонике от длительности внешней силы для различных значений радиусов пузырьков.

Рис. 2. Спектр плотности мощности переизлученного пузырьком звука при следующих значениях определяющих параметров: добротность Q = 10, частота п = 1.2, амплитуда накачки ^ = [0.5; 1.25; 2.0] - изменяется слева направо. Радиусы пузырьков Я0 = [2.66; 5.32; 7.98] - изменяются сверху вниз.

5С/-

т---1-1-1-1-г

т

ЦП,"

Рис. 3. Зависимость отклика сигнала на субгармонике от длительности внешней силы для различных значений радиусов пузырьков. Добротность Q = 10, частотa накачки rjs = 12, амплитуда накачки ss = 1.25. График a -радиус пузырьков R0 = 5.32. График b - радиус пузырьков R0 = 2.66. График c - радиус пузырьков R^ = 7.98.

По оси абсцисс отложена Timp - длительность внешней силы, по оси ординат - SGr - значение отклика сигнала на субгармонике.

5Сг

Я*-_ а

0,00-

0,20-

0,15-

0,10-

0,05-

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

Рис. 4. Зависимость отклика сигнала на субгармонике от амплитуды накачки для различных значений радиусов пузырьков. Добротность Q = 10, частотa накачки п = 12, амплитуда накачки ^ = изменялась от 0.5 до 2. График а - экспериментальные данные из статьи [2]. График ЬI радиус пузырьков Я0 = 5.32. График Ь2 - радиус пузырьков И0 = 2.66. График Ь3 -радиус пузырьков Я = 7.98. По оси абсцисс отложена амплитуда, по оси ординат - БОг - значение отклика сигнала на субгармонике.

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

В заключение была проведена оценка зависимости отклика сигнала на субгармонике (внешняя частота п = 1.2) от амплитуды накачки для различных значений радиусов пузырьков. Численный эксперимент подтверждает результаты экспериментальных данных [2], что продемонстрировано на рис. 4.

При амплитудах акустического сигнала, меньших э*=рт/рь эта зависимость близка к степенной или экспоненциальной, при 5>5* - наблюдается насыщение. Таким образом,

представляет собой пороговое значение для развития субгармоники. Исходя из этого,

*

можно сделать три вывода. Первое, субгармоника возникает при рт<5 Второе, субгармоника резко усиливается при рт~$рк$. Третье, субгармоника имеет насыщение при

Pm>s*Pks , т.е. вне зависимости от амплитуды поля уровень субгармоники не увеличивается.

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

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

Автор выражает благодарность Максимову А. О. за плодотворное обсуждение.

Работа выполнена при поддержке грантов РФФИ 04-02-16412 и ДВО РАН 04-3-А-

07-038.

СПИСОК ЛИТЕРАТУРЫ:

[1] Neppiras E.A. Subharmonic and other low-frequency emission from bubbles in sound-irradiated liquids // J. Acoust. Soc. Am. 1969. V. 46. N. 3. P. 587-601.

[2] Мансфельд А.Д., Рейман А.М. Особенности обнаружения газовых пузырьков в неоднородных нелинейных средах // В кн.: Ультразвуковая диагностика. Горький: ИПФ АН. 1983. С. 151-161.

[3] Кирилов А.Г., Мансфельд А.Д., Рейман А.М., Чичагов П.К. // В кн. Проблемы нелинейной акустики. 2 часть. Новосибирск: Ин-т Гидродинамики. 1987. C. 32-34.

[4] Losberg O., Hovem J.M., Aksum B. Experimental observation of subharmonic oscillations in Infoson bubbles // J. Acoust. Soc. Am. 1996. V. 99. N. 3. P. 1366-1369.

[5] Максимов А.О., Соседко Е.В. Переходные процессы в окрестности порога субгармонического резонанса // Письма в ЖТФ. 1999. Т. 25. No 17, С. 1-6.

[6] Соседко Е.В. Особенности нелинейный резонансов и их проявления в акустике микронеодногродных сред // Диссертация на соискание ученой степени к.ф-м.н. ТОИ ДВО РАН, Владивосток: 2003. 132c.

[7] Gong X.F. , Lu R.R. , Yu J.F. , and Shi W.T. Measurement of Subharmonics in Liquid Containing Encapsulated Microbubbles / Proceedings 16th International Symposium on nonlinear acoustics, 2002.