Капиллярная турбулентность на поверхности нормального и сверхтекучего гелия-4 тема автореферата и диссертации по физике, 01.04.07 ВАК РФ
Абдурахимов, Леонид Викторович
АВТОР
|
||||
кандидата физико-математических наук
УЧЕНАЯ СТЕПЕНЬ
|
||||
Черноголовка
МЕСТО ЗАЩИТЫ
|
||||
2010
ГОД ЗАЩИТЫ
|
|
01.04.07
КОД ВАК РФ
|
||
|
На правах рукописи
00461234?
Абдурахимов Леонид Викторович
КАПИЛЛЯРНАЯ ТУРБУЛЕНТНОСТЬ НА ПОВЕРХНОСТИ НОРМАЛЬНОГО И СВЕРХТЕКУЧЕГО ГЕЛИЯ-4
01.04.07 — физика конденсированного состояния
Автореферат диссертации на соискание учёной степени кандидата физико-математических наук
1 1 НОЯ 2010
Черноголовка — 2010
004612347
Работа выполнена в Учреждении Российской академии наук Институт физики твердого тела РАН.
Научный руководитель: Официальные оппоненты:
Ведущая организация:
доктор физико-математических наук Александр Алексеевич Левченко
доктор физико-математических наук, профессор
Валерий Степанович Цой
доктор физико-математических наук Александр Иванович Дьяченко
"Учреждение Российской академии наук Институт теплофизики им. С.С. Кутателадзе СО РАН
Защита состоится " ,9 » 2010 г. в ч. (Р& мин. на засе-
дании диссертационного совета Д 002.100.01 при Учреждении Российской академии наук Институт физики твердого тела РАН по адресу: 142432, Московская обл., г. Черноголовка, ул. Институтская 2.
С диссертацией можно ознакомиться в библиотеке Учреждения Российской академии наук ИФТТ РАН.
Автореферат разослан " £ " 2010 г.
Учёный секретарь диссертационного совета' доктор физ.-мат. наук В. Н. Зверев
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Объект исследования и актуальность темы. В данной работе представлены результаты экспериментов по капиллярной турбулентности на поверхности нормального и сверхтекучего гелия-4. Интерес к проблеме капиллярной турбулентности объясняется несколькими причинами. Во-первых, капиллярные волны интересны сами по себе. При условиях нулевой гравитации именно они являются основным типом поверхностных волн. Также капиллярные волны играют важную роль в динамике морской поверхности, так как вязкостная диссипация существенна именно на масштабах капиллярного диапазона (в отличии от слабозатухающих гравитационных волн). Во-вторых, капиллярные волны интересны как модельная система для проверки предсказаний теории волновой турбулентности, которая описывает широкий класс-слабо нелинейных систем не только в гидродинамике, но и в других областях физики (плазма, нелинейная оптика, физика твердого тела и т.д.).
Эксперименты по изучению капиллярной турбулентности на поверхности квантовых жидкостей (нормального и сверхтекучего гелия-4, жидкого водорода) представляют особый интерес. В работах по турбулентности в системе капиллярных волн на поверхности жидкого водорода М.Ю. Браж-никовым (1-6] была использована оригинальная экспериментальная методика, основанная на возбуждении волн на заряженной поверхности водорода переменным электрическим полем и их последующей регистрации по осцилляциям мощности лазерного луча, отраженного от колеблющейся поверхности. Данная методика позволяет легко контролировать амплитуду, частоту и спектральную характеристику возбуждающей силы, что является преимуществом по сравнению с другими экспериментальными методиками [7-9]. Эту методику удобно использовать для исследования капиллярной турбулентности на поверхности нормального и сверхтекучего гелия-4. Ранее турбулентность на поверхности жидкого гелия-4 экспериментально не изучалась. Жидкий гелий-4 при температурах ниже 2.17К становится сверхтекучим и обладает чрезвычайно малой вязкостью по сравнению с классическими жидкостями. Можно считать, что сверхтекучий гелий-4 -жидкость, наиболее близкая к модели идеальной жидкости, поэтому ис-
пользование поверхности Не-Н для изучения капиллярной турбулентности предоставляет принципиальные преимущества по сравнению с традиционными экспериментами, например с водой.
В цели диссертационной работы входило:
1. наблюдение и изучение турбулентных спектров капиллярных волн на поверхности нормального и сверхтекучего гелия-4
2. анализ полученных экспериментальных результатов на предмет выяснения влияния ограниченной геометрии экспериментальной ячейки на турбулентность
3. исследование влияния свойств жидкости (в первую очередь, вязкости) на характеристики спектров
4. изучение влияния типа накачки (монохроматическая, шумовая) на характеристики спектров капиллярной турбулентности на поверхности квантовых жидкостей.
Для достижения этих целей требовалось решение следующих задач:
1. создание и отладка экспериментальной методики для изучения капиллярной турбулентности на поверхности нормального и сверхтекучего гелия-4
■ 2. наблюдение турбулентных спектров капиллярных волн при разных значениях температуры жидкого гелия-4, при разных амплитудах и типах накачки
3. установление зависимости характеристик турбулентных спектров от свойств жидкости
4. анализ полученных экспериментальных результатов на соответствие опубликованным теоретическим работам.
Научная новизна. В данной работе впервые получены турбулентные спектры капиллярных волн на поверхности нормального и сверхтекучего гелия-4 в диапазоне температур 1.7-4.2 К. Наблюдено и объяснено новое явление - формирование локального максимума в области высоких
частот в турбулентном спектре капиллярных волн на поверхности Не-Н в случае гармонической накачки. Впервые экспериментально установлена зависимость частоты высокочастотного края турбулентного распределения от свойств жидкости (плотности, кинематической вязкости, коэффициента поверхностного натяжения). Впервые экспериментально установлено влияние спектральной ширины накачки на форму турбулентного спектра в случае капиллярной турбулентности на поверхности квантовых жидкостей.
Практическая ценность. В данной работе в рамках теории волновой турбулентности предложена качественная модель, описывающая условия реализации дискретного турбулентного режима. Также получена формула зависимости положения высокочастотного края инерционного интервала от амплитуды и частоты гармонической накачки и свойств жидкости. Результаты экспериментов согласуются с теоретическими оценками.
Основные положения, выносимые на защиту, состоят в следующем:
1. создана экспериментальная установка для исследования капиллярной турбулентности на поверхности нормального и сверхтекучего гелия-4
2. впервые наблюдено формирование локального максимума в турбулентном спектре капиллярных волн на поверхности сверхтекучего гелия-4 при монохроматической накачке
3. предложена качественная модель, описывающая условия реализации дискретного турбулентного режима
4. впервые исследовано влияние свойств жидкости на положение высокочастотного края инерционного интервала
5. впервые изучено влияние спектральной характеристики возбуждающей силы на степенную зависимость турбулентного распределения капиллярных волн.
Личный вклад автора. Автор непосредственно участвовал в постановке задач исследований, в выполнении экспериментов и в обсуждении полученных результатов. Диссертационная работа выполнена в лаборатории квантовых кристаллов ИФТТ РАН в период с 2004 г. по 2010 г.
Апробация работы. Результаты работы докладывались и обсуждались на конференциях:
1. "International Symposium on Quantum Fluids and Solids QFS2007"
2. Международная конференция "Low Temperature Physics - LT25"
3. Международная конференция "Solitons, Collapses and Turbulence:
Achievements, Developments and Perspectives" (Черноголовка, 2009)
4. 'International Symposium on Quantum Fluids and Solids QFS2009"
5. "XXXV Совещание по физике низких температур HT-35"
6. Международная конференция "Cryocrystals and Quantum Crystals
CC2010"
7. "International Symposium on Quantum Fluids and Solids QFS2010"
8. Сессии Совета РАН по нелинейной динамике.
Публикации. Основные результаты диссертационной работы опубликованы в 6 статьях, список которых приведен в конце автореферата. Общее количество публикаций по теме диссертации —10. Работы, вошедшие в диссертацию, были выполнены при поддержке РФФИ (гранты № 06-02-17253, 07-02-00728, 09-02-01146), Минобрнауки (НОЦ "Статистическая гидродинамика"), гранта по государственной поддержке ведущих научных школ НШ-4415.2010.2, и программ Президиума РАН ("Квантовая макрофизика", "Квантовая физика конденсированных сред", "Фундаментальные проблемы нелинейной динамики").
Объем и структура диссертации. Диссертационная работа состоит из введения, четырех глав и заключения. Общий объем диссертации -страниц текста, включая -3/ рисунков и список литературы из З® наименований.
Основное содержание работы.
Глава 1 содержит введение в предмет исследования, дан краткий обзор теоретических и экспериментальных результатов по исследованию капиллярной турбулентности. Обсуждаются свойства жидкого гелия-4 в свя-
зи с экспериментами по капиллярной турбулентности. Сформулированы основные задачи исследований.
В главе 2 описаны экспериментальная установка, методика возбуждения и регистрации капиллярных волн.
Раздел 2.1. Приведено описание экспериментальной установки и метода возбуждения волн на заряженной поверхности жидкого гелия-4.
Для исследования капиллярной турбулентности на поверхности жидкого гелия использовалась экспериментальная методика, которая была развита ранее в экспериментах по изучению капиллярной турбулентности на поверхности жидкого водорода [5]. Эксперименты проводились в цилиндрической оптической ячейке, установленной в вакуумной полости гелиевого криостата. Конструкция ячейки схематически показана на рис. 1.
42
6"
Рис. 1. Устройство экспериментальной ячейки. 1- текстолитовый брусок, 2- радиоактивная мишень, 3- медный контейнер, 4- верхняя обкладка конденсатора, 5- кварцевые окна, 6- медный холодопровод, 7- капилляр набора
Газообразный гелий поступал из баллона по капилляру в ячейку и конденсировался в медный стакан. Набор гелия производился до тех пор,
пока поверхность жидкости не достигала края стакана. Свободная поверхность жидкого гелия заряжалась с помощью расположенной на дне стакана радиоактивной мишени, излучающей ^-электроны. Для этого к медному стакану прикладывалось постоянное положительное напряжение 17 « 600 В относительно верхней пластины. При этом под поверхностью жидкости накапливались положительные заряды, которые появлялись в результате ионизации ^-электронами тонкого слоя жидкости вблизи дна. Таким образом, заряженная поверхность жидкого гелия и верхняя металлическая пластина образовывали плоский конденсатор.
Колебания заряженной поверхности жидкого гелия, представляющие собой цилиндрические стоячие волны, возбуждались переменным напряжением, прикладываемым к медному стакану в дополнение к постоянному, на одной из частот собственных колебаний поверхности. Амплитуды переменного напряжения были много меньше и и составляли 1-20 В.
Эксперименты проводились при фиксированной температуре в диапазоне 1.7-4.2 К.
В разделе 2.2 дано описание методики регистрации волн на поверхности жидкости. Схема измерений представлена на рис. 2.
Поверхностные волны регистрировались оптическим методом. Лазерный луч падал на поверхность жидкого гелия-4 под скользящим углом а —0.15 рад таким образом, чтобы он находился в диаметральной плоскости цилиндрического медного стакана, в котором был сконденсирован жидкий гелий. Луч отражался от колеблющейся поверхности и фокусиро-
верхняя обкладка
контейнер
Рис. 2. Схема измерений поверхностных колебаний
вался линзой на фотодетектор. Сигнал с фотодетектора, пропорциональный мощности отраженного лазерного луча, поступал на вход усилителя. Усиленный сигнал с фотодетектора P(t) оцифровывался платой высокоскоростного АЦП и записывался в память компьютера. С помощью алгоритма FFT (Быстрое Фурье Преобразование), вычислялся спектр мощности отраженного лазерного луча Я2.
В проведенных экспериментах длины волн были меньше характерного размера пятна лазерного луча на поверхности (режим "широкого" луча). В этом случае, как было показано в [5], мощность отраженного лазерного луча является интегральной характеристикой участка поверхности жидкого гелия, ограниченного пятном лазерного луча:
где Я,(а + <р(г, ¿)) - коэффициент отражения, <р(г, ¿) - угол отклонения поверхности от положения равновесия (<р « а), С\ и Сг - некоторые константы. Отсюда, спектр мощности отраженного лазерного луча Р2 связан с угловыми амплитудами волн следующим выражением [5]:
где ко, - волновое число волны с частотой и.
При сравнении экспериментальных результатов с предсказаниями теории волновой турбулентности принято использовать понятие Фурье-образа корреляционной функции отклонений поверхности от положения равновесия 1и:
где 1%, - Фурье-образ отклонений поверхности от положения равновесия г)(£) в выбранной точке поверхности. Так как угловые амплитуды <ри связаны с амплитудами ц^ соотношением ¡ри = то окончательно получаем:
P(t) ~ / R(a + v>(r, t))da ~ Сх + С2 / p(r, t)da,
Р2 ~ ф2 /к2
1 U1 гы/ ""Ы»
1> =< ll >
Таким образом, спектр мощности отраженного лазерного луча пропорционален спектру волн 1Ы.
Главы 3, 4, 5 посвящены результатам экспериментальных исследований капиллярной турбулентности на поверхности квантовых жидкостей и их обсуждению.
В главе 3 представлены результаты экспериментов на поверхности сверхтекучего гелия-4, в которых впервые было наблюдено накопление волновой энергии в области высоких частот (формирование локального максимума в турбулентном спектре), вызванное влиянием ограниченной геометрии экспериментальной ячейки.
На рис. 3 в двойном логарифмическом масштабе изображен спектр капиллярных волн на поверхности сверхтекучего гелия-4, полученный в эксперименте при интенсивной гармонической накачке на частоте 79.7 Гц.
Рис. 3. Распределение Р% на поверхности сверхтекучего гелия при Т=1.7К. Усреднение по 16-и файлам. Продолжительность записи одного файла 5 секунд. Поверхность накачивалась синусоидальной силой на частоте 79.7 Гц, амплитуда напряжения накачки 11В. Штриховая линия соответствует степенному закону ~ и>~3'7, предсказанному теорией для случая узкополосной накачки.
На рисунке отчетливо видна серия гармоник. Основная гармоника расположена на частоте накачки шр ~ 79.7 Гц. Остальные гармоники нахо-
дятся на частотах, кратных частоте накачки и)„ = п ■ мр. В рамках теории Захарова [10] этот факт объясняется тем, что закон дисперсии капиллярных волн и>2 = {а/р)к3 является распадным, и поэтому основной вклад во взаимодействие волн вносят трехволновые процессы - распад волны на две с сохранением суммарного волнового вектора кз = + кг и суммарной частоты о^з = и>\ + Ш2 (т.е. импульса и энергии), а также обратный ему процесс слияния двух волн в одну. В инерционном интервале амплитуды гармоник зависят от частоты по степенному закону Р2 ~ с индексом т рз 3.7, что согласуется с предсказаниями теории для случая узкополосной накачки ¡11]. На частотах порядка 4кГц амплитуды волн значительно снижаются из-за вязкостной диссипации и выходят на уровень аппаратного шума.
100 1000 ы
Частота, Гц
Рис. 4. Турбулентный спектр после снижения амплитуды накачки. Частота накачки 79.7Гц, амплитуда напряжения накачки 10 В. На частотах порядка 2.5 кГц наблюдается локальный максимум (выделено штриховой кривой). Штриховая прямая линия соответствует закону ~
Когда амплитуда накачки была уменьшена, форма турбулентного спектра изменилась (рис. 4). Высокочастотный край турбулентного каскада сдвинулся в сторону низких частот, что согласуется с предыдущими экспериментами на водороде [2]. Но кроме того, впервые было обнаружено формирование локального максимума в спектре волн на частотах поряд-
ка частоты высокочастотного края турбулентного каскада шъ, что можно интерпретировать как накопление волновой энергии (показано штриховой кривой на рис. 4).
Рис. 5. Турбулентный спектр в случае накачки на шр/2тг — 34 Гц. Формирование максимума в спектре отчетливо видно на основном рисунке (отмечено окружностью) и на вставке в увеличенном масштабе.
Формирование локального максимума наблюдалось в ряде экспериментов с гармонической накачкой. Например, когда поверхность накачивалась синусоидальной силой на 34 Гц, был получен спектр волн с отчетливым локальным максимумом (рис. Б). Общей чертой проведенных экспериментов было то, что локальный максимум формировался около высокочастотного края турбулентного каскада иь, хотя форма локального максимума и его точное положение зависели как от частоты накачки, так и от амплитуды.
В случае шумовой накачки, локальный максимум не наблюдался. Так, на рис. 6 представлено турбулентное распределение, полученное при накачке поверхности шумовым сигналом в диапазоне частот 60-130 Гц.
Формирование локального максимума можно объяснить влиянием ограниченной геометрии экспериментальной ячейки. Собственные частоты поверхностных колебаний в цилиндрической ячейке не эквидистантны. Так, расстояние между двумя ближайшими резонансными частотами опи-
100
1000
1
10'
Частота, Гц
100 1000 Частота, П*
Рис. в. Турбулентный сяектр капиллярных волн в случае шумовой накачки на 60-130 Гц.
сывается выражением:
где В - диаметр медного стакана, «тир коэффициент поверхностного натяжения и плотность жидкого гелия-4. Как упоминалось выше, в силу закона сохранения частоты трех-волновых процессов +си2 = и3> частоты нелинейных гармоник эквидистантны:
где и)Р - частота накачки. Таким образом, между частотами нелинейных гармоник и собственными частотами поверхностных колебаний в экспериментальной ячейке существует расстройка. Если бы резонансы ячейки не были уширены, то трех-волновые процессы нелинейного взаимодействия были бы запрещены, так как законы сохранения и)\ + = не могли бы быть выполнены на дискретном ^-пространстве резонансных частот в силу расстройки [12,13]. Однако резонансы уширены благодаря вязкостному затуханию и нелинейному взаимодействию:
где бшц - вязкое уширение, 6шп1 - нелинейное уширение, - суммарное уширение резонанса экспериментальной ячейки. Понятно, что расстройка
и>р, 2и)р, Ъи>р,..
ди> = би» + б^тй,
важна только в области частот, где относительное уширение Г = 6ш/Аш меньше единицы:
т.е. уширение резонанса меньше, чем расстояние между двумя соседними резонансными частотами (см. рис. 7).
ДСО П'Ыритр й- »1 |
Рис. 7. Когда уширение резонанса 5и> меньше расстояния между резонан-сами Аи> некоторые нелинейные гармоники расположены вне резонансов поверхностных колебаний в экспериментальной ячейке.
Вязкое уширение 5ии связано с характерным временем вязкого затухания ти и определяется следующим образом [14]:
бит,.
Г"1 = 4:Рк1
чг^'
где V - кинематическая вязкость жидкого гелия. Относительное вязкое уширение Дш изображено на рис. 8. Нелинейное уширение <5ш„г определяется характерным временем нелинейного взаимодействия гп1. которое в случае гармонической накачки зависит от частоты как ~ ш1^ [15]. Таким образом,
¿ш,
п1 '
'Ы
<Л)
-1/6
где б (А) - возрастающая функция от амплитуды накачки А. Предполагается, что 6и}п1 .те <5ш„ вблизи высокочастотного края инерционного интервала
и>ь (см. график относительного нелинейного уширения на рис. 8). Таким образом, полное относительное уширение описывается следующей формулой:
(Зтг/0)(а/Ру/Ы/з
Рис. 8. Результаты оценок относительного уширения 5ш/Да> при умеренной амплитуде гармонической накачки. Возрастающая штриховая кривая
- вязкое относительное уширение ¿ш„/Ды, убывающая штриховая кривая
- нелинейное относительное уширение ¿штг;/Дш, жирная кривая с минимумом - полное относительное уширение 6и>/Аи>.
Полное относительное уширение имеет минимум вблизи высокочастотного края инерционного интервала иъ, и при умеренной амплитуде накачки меньше единицы (Г < 1) в области частот около шь (см. рис. 8). В этой области частот реализуется дискретный режим капиллярной турбулентности: расстройка частот гармоник и резонансных частот становится существенной и определяет специфическую форму турбулентного распределения (формирование локального максимума). При больших амплитудах накачки полное относительное уширение больше единицы, Г > 1, для всех частот, и система становится квазинепрерывной, реализуется кинетический режим турбулентности.
В главе 4 представлены результаты изучения влияния свойств жидкого гелия-4 на положение высокочастотного края инерционного интервала.
Инерционный интервал на низких частотах ограничен областью накачки, а на высоких — диапазоном частот, где происходит изменение механизма распространения энергии волны: перенос энергии по каскаду, обусловленный нелинейным взаимодействием волн, заменяется трансформацией механической энергии волн в тепло из-за вязкого трения. Высокочастотный край инерционного интервала можно определить как частоту и)Ь, при которой характерное время нелинейного взаимодействия волн т^ сравнивается по порядку величины с временем вязкого затухания г„, т„1 ~ т„. В рамках этого определения можно показать, что в случае узкополосной накачки изменение частоты высокочастотного края турбулентного каскада изъ описывается степенной функцией амплитуды волны на частоте накачки А, частоты накачки коэффициента кинематической вязкости и, плотности жидкости р и коэффициента поверхностного натяжения а.
Обратное время нелинейного взаимодействия капиллярных волн прямо пропорционально квадрату коэффициента нелинейного взаимодействия V ~ а>3/2/у/а [10]:
1 /ты ~ №пк/шк ~
Время вязкого затухания капиллярных волн возрастает с уменьшением коэффициента вязкости V ]14]:
1 / Ту = 4ук ~ V (р/сг)2^3.
Полагая, что /Ы/Ыр А2(ш/мр)~21/6, получаем
шь ~ (1)
Амплитудная зависимость граничной частоты была подробно исследована в экспериментах на жидком водороде [2]. Было показано, что экспериментальные результаты шь{А), полученные на разных частотах накачки шр, хорошо описываются полученным выше выражением. Продолжение исследований на поверхности жидкого гелия позволяет изучить зависимость положения частоты высокочастотного края шь от величины коэффициента кинематической вязкости жидкости к и поверхностного натяжения а, так как они значительно изменяются с повышением температуры в интервале ог1.8Кдо4.2К [16].
На рис. 9 показана частотная зависимость распределения при возбуждении поверхности гармонической силой на резонансной частоте собственных колебаний поверхности жидкости шр/2тг = 25 Гц при температуре жидкого гелия Т'=1.8 К.
Частота, Гц
Рис. 9. Распределение на поверхности сверхтекучего гелия при Т=1.8К. Стрелка указывает положение высокочастотного края инерционного интервала, шъ/Ъ-п = 750 ± 100 Гц
Угловая амплитуда волны на частоте накачки составляла приблизительно 0.015 радиан, при этом амплитуда возбуждающей волны равнялась 0.007 мм. На рисунке отчетливо видно, как степенное распределение переходит в более крутую зависимость на высоких частотах, т.е. наблюдается высокочастотный край инерционного интервала. Граница высокочастотного края оценивалась из экспериментальных зависимостей как частота, при которой отчетливо наблюдается отклонение зависимости от степенного распределения. Положение высокочастотного края на рисунке отмечено стрелкой и составляет шь/2-к — 750 ± 100 Гц.
Повышение температуры от 1.8 К до 2.3 К привело к кардинальному изменению в распределении Р% по сравнению с распределением, представленным на рис. 9 (гелий перешел из сверхтекучего состояния в нормальное состояние при Т=2.17К). На рис. 10.показано распределение Р%, которое было получено при температуре 2.3 К. Амплитуда волны на частоте накачки составляла 0.006 мм. Хорошо видно, что на частоте около 310Гц
(±50 Гд) наблюдается особенность в распределении Р% (отмечена стрелкой на рисунке), которая может быть интерпретирована как край инерционного интервала.
Рис. 10. Распределение Р% на поверхности нормального гелия при Т=2.3 К. Стрелкой обозначен высокочастотный край инерционного интервала, 27г = 310 ± 50 Гц
На распределении виден "пьедестал", идущий от низких частот до частоты порядка 1 кГц. Формирование "пьедестала" можно объяснить кипением нормального гелия.
На рис. 11 представлены результаты экспериментов, выполненных при температуре 4.2 К. . Амплитуда волны на частоте накачки составляла 0.004 мм. Интенсивность кипения жидкости значительно уменьшилась. На рисунке видно, что высокочастотный край распределения сдвинулся в сторону низких частот [у>ъ!Ъ\ = 250 ± 50 Гц) по сравнению с распределением, полученным при Т=2.3К, "пьедестал" уменьшился по амплитуде, а его высокочастотный край также переместился в сторону низких частот.
Оценка отношений частот высокочастотного края инерционного интервала при температурах 1.8К и 2.3К по формуле (1) с учетом изменений амплитуды волны на частоте накачки, коэффициента поверхностного натяжения, плотности жидкости и кинематической вязкости дает величину ^1.8/^2.3 = 2. В тоже время экспериментальная величина отношения W1.8/W3 близка к 2.4 ± 0.5. При температурах 2.3К и 4.2К теоретическая
Рис. 11. Распределение на поверхности нормального гелия при Т=4.2 К. Стрелкой обозначен высокочастотный край инерционного интервала, Шь/2тг = 250 ± 50 Гц.
оценка отношения равна ^2.3/^4.2 = 1-3, а установленное из экспериментальных зависимостей отношение равно 1.2 ± 0.3. С учетом того, что зависимость и>ъ от амплитуды волны на частоте накачки была подробно исследована и подтверждена ранее в экспериментах с жидким водородом [2], можно констатировать, что выражение (1) удовлетворительно описывает экспериментально определенное изменение высокочастотного края с ростом амплитуды волны на частоте накачки, значения частоты накачки, плотности, коэффициента кинематической вязкости и коэффициента поверхностного натяжения.
Глава 5 содержит результаты экспериментального изучения влияния спектральной характеристики возбуждающей силы на свойства турбулентного спектра.
Внутри инерционного интервала корреляционная функция отклонений поверхности от равновесного плоского состояния в Фурье представлении (турбулентный каскад) описывается степенной функцией частоты 1и ~ ш~т. Величина индекса т зависит от спектральной характеристики возбуждающей силы. При возбуждении поверхности жидкости низкочастотным шумом в широкой полосе Дш, то есть когда ширина полосы
превосходит характерную частоту накачки изр, Аш > и>р, корреляционная функция 1Ш описывается степенной функцией ы~т с индексом т = 17/6. Численный расчёт (17] дал величину показателя т близкую к теоретическому значению. Компьютерные моделирования [11| показали, что по мере уменьшения ширины полосы накачки шумом Аи/ на турбулентном каскаде формируется ряд пиков равноудалённых друг от друга, а их ширина линейно зависит от частоты. При накачке в узкой полосе Аш < мр уменьшение высоты пиков с ростом частоты описывается степенной функцией частоты с показателем степени, абсолютная величина которого на единицу превосходит величину показателя для случая накачки шумом в широкой полосе, т.е. то = 23/6 « 3.8. Экспериментальные исследования на поверхности жидкого водорода [4] показали, что при переходе от возбуждения поверхности шумом в широкой полосе к накачке гармонической силой на одной резонансной частоте ячейки величина показателя степени уменьшается, что качественно согласуется с результатами работы [11]. Однако изменение частотной зависимости корреляционной функции при переходе от широкополосной к узкополосной накачке, как это было сделано в [11], ранее не исследовалось. Поэтому одна го целей проведенных исследований заключалась в подробном изучении эволюции турбулентного каскада при изменении ширины полосы Аи> от широкополосной шумовой накачки к узкополосному шумовому возбуждению. В экспериментах на поверхности сверхтекучего гелия-4 при широкополосной шумовой накачке не удалось получить спектр капиллярных волн, который был бы достаточно широким для надежного определения показателя степенной зависимости в инерционном интервале. Поэтому, эксперименты по изучению эволюции турбулентного каскада при изменении ширины полосы шумовой накачки были проведены на поверхности жидкого водорода. Конструкция экспериментальной ячейки и порядок работы подробно описаны в [5]. Эксперименты проводились при температуре жидкого водорода около 15.5 К.
На рис. 12 показано распределение (тёмно-серый график), когда ширина полосы электрического сигнала У{£) нахачки приложенного к медному стакану составляла 64 Гц, от 39 Гц до 103 Гц, то есть поверхность возбуждалась широкополосным шумом. В частотном интервале от 200 Гц до
Частота, Гц
Рис. 12. Турбулентный каскад на поверхности жидкого водорода при накачке случайной силой в полосе 39-103 Гц при двух уровнях накачки. Прямая линия — степенная зависимость от17/6.
8 кГц сформировался турбулентный каскад, частотную зависимость которого можно описать степенной функцией с показателем т — 2.8±0.1. Для сравнения на рисунке прямой линией показана функция, прямо пропорциональная ш~17/6. Светло-серый график на рис. 12 соответствует спектру при накачке в том же самом частотном интервале, но в два с половиной раза уменьшенной амплитудой.
Когда ширина полосы накачки шумом Аш была уменьшена относительно средней частоты накачки и).р, так что Аш и и>р/2, в спектре появились несколько широких пиков. Распределение при накачке шумом в интервале 57-89 Гц показано на рис. 13. Первый пик на частотах 57-89 Гц соответствует волнам возбуждаемым непосредственно внешней силой. Второй и третий, расположенные в диапазонах «120-170 Гц и «180-250 Гц, соответствуют волнам родившимся в результате нелинейного взаимодействия. Расстояние между центрами пиков приблизительно равно средней частоте накачки шр = 73 Гц. Свыше 250 Гц сформировано непрерывное турбулентное распределение с инерционным интервалом до ~ 6 кГц.
При дальнейшем уменьшении ширины полосы накачки пики на турбулентном каскаде становятся более выраженными, минимумы более глу-
102 103 Частота, Гц
Рис. 13, Спектр Р2 колебаний поверхности при накачке шумом в полосе 57-89 Гц.
и 10е ¿10*
Й 104 102
Рис. 14. Спектр Р2 при накачке случайной силой в полосе 71-75 Гц. Прямая линия соответствует зависимости а*-23/6.
Частота., Гц
бокими. На рис. 14 показано распределение при возбуждении поверхности шумом в полосе шириной 4 Гц, от 71 Гц до 75 Гц. Около 30 пиков отчётливо проявились на турбулентном каскаде. Расстояние между центрами пиками сохранилось и равняется « 73 Гц. Инерционный интервал простирается до « 15 кГц. Частотная зависимость положения максимумов пиков в пределах инерционного интервала близка к степенной функции ^-з.8±о.1 Отметим, что разница показателей степенной зависимости для случаев узкополосной и широкополосной накачки равна 1 ± 0.2.
Частота, Гц
Рис. 15. Зависимость ширины Д/ пиков в показанном на рис. 14 спектре от частоты. Сплошная линия — линейная зависимость 0.027w.
Ширина пиков Д/ с ростом частоты возрастает. Экспериментальная зависимость ширины пиков от частоты при накачке шумом в полосе 7175 Гц представлена на рис.15. Очевидно, что наблюдаемый рост описывается линейной функцией и наблюдается удовлетворительное согласие между экспериментальными данными и предсказанной теорией линейной зависимостью.
Основные результаты и выводы. В диссертационной работе выполнены экспериментальные исследования капиллярной турбулентности на заряженной поверхности нормального и сверхтекучего гелия-4 и жидкого водорода. Использование заряженной поверхности квантовых жидкостей для изучения турбулентности на поверхности жидкости предоставляет принци-
пиальные преимущества по сравнению с традиционными экспериментами. Примененные новые идеи и методики позволили детально изучить турбулентность в системе капиллярных волн и получить новую информацию о влиянии ограниченной геометрии экспериментальной ячейки на турбулентные процессы.
1. Впервые наблюдено накопление волновой энергии на высоких частотах в турбулентном спектре капиллярных волн на поверхности сверхтекучего гелия-4 при гармонической накачке. Предложено качественное объяснение этого явления, учитывающее дискретность собственных частот колебаний поверхности в цилиндрическом резонаторе.
2. Исследовано положение высокочастотного края турбулентного распределения в зависимости от величины коэффициента кинематической вязкости, коэффициента поверхностного натяжения и плотности жидкости, и показано, что экспериментальные результаты хорошо согласуются с теоретическими оценками, выполненными в рамках теории слабой волновой турбулентности.
3. Экспериментально наблюдено, что изменение спектральной характеристики возбуждающей накачки шумом (ширины полосы) приводит к качественной перестройке турбулентного каскада в системе капиллярных волн на поверхности жидкого водорода. При широкополосной накачке турбулентное распределение хорошо описывается монотонной убывающей степенной функцией частоты с показателем близким к т = —2.8 ±0.1. В случае узкополосной накачки турбулентный каскад представляет собой ряд эквидистантных пиков, высота которых описывается степенной функцией с показателем —3.8 ± 0.1. Полученные экспериментальные результаты хорошо согласуются с выводами теории волновой турбулентности для капиллярных волн.
Основное содержание диссертационной работы изложено в следующих публикациях:
[1*] G.V. Kolmakov, A. A. Levchenko, M.Yu. Brazhnikov, L.V. Abdurakhimov, P.V. E. McClintock, and L.P. Mezhov-Deglin, Capillary turbulence on the surface of quantum fluids, Progress in Low Temperature Physics, Vol. 163 on Quantum Turbulence, edited by Halperin & Tsubota, ISBN 978-0-08-054810-4 (2008)
[2*] L.V. Abdurakhimov, M.Yu. Brazhnikov, A.A. Levchenko and L.P. Mezhov-Deglin, Developed Capillary Turbulence on the Surface of Normal and Superfluid 4He, J. Low Temp. Phys. 150, 426-430 (2008)
[3*] JI. В. Абдурахимов, M. Ю. Бражников, А. А. Левченко, Распределение вероятпости осцилляций поверхности жидкого водорода в турбулентном режиме, Письма в ЖЭТФ 88, 21-26 (2008)
[4*] Л. В. Абдурахимов, М.Ю. Бражников, А. А. Левченко, Капиллярная турбулентность на поверхности нормального и сверхтекучего 4Не, Физика Низких Температур 35, 127-133 (2009)
[5*] Л. В. Абдурахимов, М. Ю. Бражников, А. А. Левченко, Эволюция турбулентного каскада на поверхности жидкого водорода при изменении спектральной характеристики возбуждающей силы, Письма в ЖЭТФ 89,139-142 (2009)
[6*J L.V. Abdurakhimov, M.Yu. Brazhnikov, I.A. Remizov, A.A. Levchenko, Observation of wave energy accumulation in the turbulent spectrum of capillary waves on the He-II surface under harmonic pumping, Письма в ЖЭТФ 91, 291 (2010)
Список цитированной литературы:
[1] М.Ю. Бражников, Г.В. Колмаков, А.А. Левченко, Л.П. Межов-Деглин, Письма в ЖЭТФ 73, 443 (2001)
[2] М.Ю. Бражников, Г.В. Колмаков, А.А. Левченко, Л.П. Межов-Деглин, Письма в ЖЭТФ 74, 660 (2001)
[3] М.Ю. Бражников, Г.В. Колмаков, А.А. Левченко, Л.П. Межов-Деглин, Физика Низких Температур 27, 1183 (2001)
[4| М.Ю. Бражников, Г.В. Колмаков, А.А. Левченко, ЖЭТФ 122, 521 (2002)
[5] М.Ю. Бражников, А.А. Левченко, Л.П. Межов-Деглин, ПТЭ 6, 31 • (2002)
[6] G.V. Kolmakov, A. A. Levchenko, M.Yu. Brazhnikov, L. P. Mezhov-Deglin, A. N. Silchenko, P.V.E. McClintock, Phys. Rev. Lett. 93, 074501 (2004)
[7] W.B. Wright, R. Budakian, S. J.. Putterman, Phys. Rev. Lett. 76, 4528 (1996)
[8] Mark-Tiele Westra, Patterns and Weak Turbulence in Surface Waves, Eindhoven (2001)
[9] E. Falcon, C. Laroche, S. Fauve, Phys. Rev. Lett. 98, 094503 (2007)
[10] V. Zakharov, V. L'vov, G. Falkovich, Kolmogorov Spectra of Turbulence, vol. 1, Springer, Berlin (1992)
[И] Г.Е. Фалькович, А. Б. Шафаренко, ЖЭТФ 94,172 (1988)
[12] E. A. Kartashova, Physica D 46, 43 (1990)
[13] E.A. Kartashova, Physica D 54,125 (1991)
[14] Л.Д. Ландау, E. M. Лифшиц, Теоретическая физика, т.6, Гидродинамика, Наука, Москва (1988)
[15] Г. В. Колмаков, Письма в ЖЭТФ 83, 64 (2006)
[16] R. J. Donnelly, С. F. Barenghi, J. Phys. Chem. Ref. Data 27,1217 (1998)
[17] A.N. Pushkarev, V.E. Zakharov, Phys. Rev. Lett. 76, 3320 (1996)
Сдано в печать 21.09.10. Формат 60x90 1/16 Объём 1,5 пл. Заказ 256. Тираж 80 Отпечатано в типографии ИПХФ РАН 142432, Московская обл., г. Черноголовка, пр-т академика Семёнова, 5 Тел. 8(49652)2-19-38
Общая характеристика работы.
Предисловие.
1 Введение
1.1 Спектр малых колебаний заряженной поверхности жидкости, собственные колебания поверхности в цилиндрической ячейке.
1.2 Волновая турбулентность.
2 Описание экспериментальной методики
2.1 Экспериментальная установка
2.2 Методика регистрации волн.
3 Формирование турбулентного каскада на поверхности жидкого гелия-4. Наблюдение локального максимума в спектре капиллярных волн при гармонической накачке поверхности
4 Влияние свойств жидкости на положение высокочастотного края инерционного интервала
5 Влияние спектральной характеристики возбуждающей силы на свойства турбулентного спектра
6 Заключения и выводы Литература
Общая характеристика работы
Эксперименты по изучению капиллярной турбулентности на поверхности нормального и сверхтекучего гелия-4 и жидкого водорода представляют особый интерес. В работах по турбулентности в системе капиллярных волн на поверхности жидкого водорода М.Ю. Бражниковым [1-6] была использована оригинальная экспериментальная-методика, основанная на возбуждении волн на заряженной поверхности водорода переменным электрическим полем и их последующей регистрации по осцилляциям мощности лазерного луча, отраженного от колеблющейся поверхности. Данная методика позволяет легко контролировать амплитуду, частоту и спектральные характеристики возбуждающей силы, что является преимуществом по сравнению с другими экспериментальными методиками [7-10]. Эту методику удобно использовать для исследования капиллярной турбулентности на поверхности нормального и сверхтекучего гелия-4. Ранее турбулентность на поверхности жидкого гелия-4 экспериментально не изучалась. Жидкий гелий-4 при температурах ниже 2.17 К становится сверхтекучим и обладает чрезвычайно малой вязкостью по сравнению с классическими жидкостями, поэтому использование поверхности Не-Н для изучения капиллярной турбулентности предоставляет принципиальные преимущества по сравнению с традиционными экспериментами, например с водой. В цели диссертационной работы входило:
1. наблюдение и изучение турбулентных спектров капиллярных волн на поверхности нормального и сверхтекучего гелия-4
2. исследование влияния свойств жидкости (коэффициента кинематической вязкости, коэффициента поверхностного натяжения, плотности жидкости) на положение высокочастотного края инерционного интервала
3. изучение влияния спектральной ширины накачки на характеристики спектров капиллярной турбулентности.
Для достижения этих целей требовалось решение следующих задач:
1. создание и отладка экспериментальной методики для изучения капиллярной турбулентности на поверхности нормального и сверхтекучего гелия-4
2. наблюдение турбулентных спектров капиллярных волн при разных значениях температуры жидкого гелия-4, при разных амплитудах и типах накачки
3. установление зависимости характеристик турбулентных спектров от свойств жидкости (коэффициента кинематической вязкости, коэффициента поверхностного натяжения, плотности жидкости).
Научная новизна. В данной работе впервые получены турбулентные спектры капиллярных волн на поверхности нормального и сверхтекучего гелия-4 в диапазоне температур 1.7-4.2 К. Наблюдено и объяснено новое явление - формирование локального максимума в области высоких частот в турбулентном спектре капиллярных волн на поверхности Не-Н в случае гармонической накачки. Впервые экспериментально установлена зависимость частоты высокочастотного края турбулентного распределения от свойств жидкости (плотности, кинематической вязкости, коэффициента поверхностного натяжения). Детально изучено влияние спектральной ширины накачки на форму турбулентного спектра в случае капиллярной турбулентности.
Практическая ценность.В данной работе в рамках теории волновой турбулентности предложена качественная модель, описывающая условия реализации дискретного турбулентного режима. Также получена формула зависимости положения высокочастотного края инерционного интервала от амплитуды и частоты гармонической накачки и свойств жидкости. Результаты экспериментов согласуются с теоретическими оценками.
Основные положения, выносимые на защиту, состоят в следующем:
1. создана экспериментальная установка для исследования капиллярной турбулентности на поверхности нормального и сверхтекучего гелия-4
2. впервые наблюдено формирование локального максимума в турбулентном спектре капиллярных волн на поверхности сверхтекучего гелия-4 при монохроматической накачке
3. предложена качественная модель, описывающая условия реализации дискретного турбулентного режима
4. впервые исследовано влияние свойств жидкости (коэффициента кинематической вязкости, коэффициента поверхностного натяжения, плотности жидкости) на положение высокочастотного края инерционного интервала
5. впервые детально изучено влияние спектральной ширины возбуждающей силы на степенную зависимость турбулентного распределения капиллярных волн.
Личный вклад автора. Автор непосредственно участвовал в постановке задач исследований, в выполнении экспериментов и в обсуждении полученных результатов. Диссертационная работа выполнена в лаборатории квантовых кристаллов ИФТТ РАН в период с 2004 г. по 2010 г.
Апробация работы. Результаты работы докладывались и обсуждались на конференциях:
1. "International Symposium on Quantum Fluids and Solids QFS2007"
2. Международная конференция "Low Temperature Physics - LT25"
3. Международная конференция "Solitons, Collapses and Turbulence: Achievements, Developments and Perspectives" (Черноголовка, 2009)
4. "International Symposium on Quantum Fluids and Solids QFS2009"
5. "XXXV Совещание по физике низких температур HT-35"
6. Международная конференция "Cryocrystals and Quantum Crystals CC-2010"
7. "International Symposium on Quantum Fluids and Solids QFS2010"
8. Сессии Совета РАН по нелинейной динамике.
Публикации. Основные результаты опубликованы в 6 статьях [11-16], список которых приведен в1 конце диссертации. Общее количество публикаций по теме диссертации -10. Работы, вошедшие в диссертацию, были выполнены при поддержке РФФИ' (гранты № 06-02-17253, 07-02-00728, 09-02-01146), Минобрнауки (НОЦ "Статистическая гидродинамика"), гранта по государственной поддержке ведущих научных школ НШ-4415.2010.2, и программ Президиума РАН ("Квантовая макрофизика", "Квантовая физика конденсированных сред", "Фундаментальные проблемы нелинейной динамики").
Объем и структура диссертации. Диссертационная работа состоит из введения, четырех глав и заключения. Общий объем диссертации - 70 страниц текста, включая 31 рисунок и список литературы из 30 наименований.
6. Заключения и выводы
В диссертационной работе выполнены экспериментальные исследования капиллярной турбулентности на заряженной поверхности нормального и сверхтекучего гелия-4 и жидкого водорода.
1. Создана экспериментальная установка для изучения капиллярной турбулентности на поверхности нормального и сверхтекучего гелия-4.
2. Создана и отлажена экспериментальная методика возбуждения волн на поверхности жидкости накачкой заданной спектральной характеристики (монохроматическая, узкополосная, широкополосная) и их регистрации.
3. Впервые наблюдены турбулентные спектры капиллярных волн на поверхности нормального и сверхтекучего гелия-4, при разных амплитудах и типах накачки.
4. Впервые наблюдено накопление волновой энергии на высоких частотах в турбулентном спектре капиллярных волн на поверхности сверхтекучего гелия-4 при гармонической накачке. Предложено качественное объяснение этого явления, учитывающее дискретность собственных частот колебаний поверхности в цилиндрическом резонаторе.
5. Впервые исследовано положение высокочастотного края турбулентного распределения в зависимости от величины коэффициента кинематической вязкости, коэффициента поверхностного натяжения и плотности жидкости, и показано, что экспериментальные результаты хорошо согласуются с теоретическими оценками, выполненными в рамках теории слабой волновой турбулентности.
6. Впервые экспериментально наблюдено, что изменение спектральной характеристики возбуждающей накачки шумом (ширины полосы) приводит к качественной перестройке турбулентного каскада в системе капиллярных волн на поверхности жидкого водорода. При широкополосной накачке турбулентное распределение хорошо описывается монотонной убывающей степенной функцией частоты с показателем близким кт = — 2.8 ±0.1. В случае узкополосной накачки турбулентный каскад представляет собой ряд эквидистантных пиков, высота которых описывается степенной функцией с показателем —3.8 ±0.1. Полученные экспериментальные результаты хорошо согласуются с выводами теории волновой турбулентности для капиллярных волн.
Использование заряженной поверхности квантовых жидкостей для изучения турбулентности на поверхности жидкости предоставляет принципиальные преимущества по сравнению с традиционными экспериментами. Применённые новые идеи и методики позволили детально изучить турбулентность в системе капиллярных волн и получить новую информацию о влиянии ограниченной геометрии экспериментальной ячейки на турбулентные процессы. В дальнейших экспериментах на поверхности жидкого гелия-4 планируется детально изучить влияние граничных условий на турбулентность, для чего будет создана экспериментальная ячейка с изменяемой геометрией границы поверхности жидкости.
1. М. Ю. Бражников, Г. В. Колмаков, A.A. Левченко, Л. П. Межов-Деглин, Письма в ЖЭТФ 73, 443 (2001)
2. М.Ю. Бражников, Г. В. Колмаков, A.A. Левченко, Л. П. Межов-Деглин, Письма в ЖЭТФ 74, 660 (2001)
3. М.Ю. Бражников, Г. В. Колмаков, A.A. Левченко, Л. П. Межов-Деглин, Физика Низких Температур 27, 1183 (2001)
4. М.Ю. Бражников, Г.В. Колмаков, A.A. Левченко, ЖЭТФ 122, 521 (2002)
5. М.Ю. Бражников, A.A. Левченко, Л.П. Межов-Деглин, ПТЭ 6, 31 (2002)
6. G. V. Kolmakov, A. A. Levchenko, М. Yu. Brazhnikov, L. Р. Mezhov-Deglin, A.N. Silchenko, P.V.E. McClintock, Phys. Rev. Lett. 93, 074501 (2004)
7. W.B. Wright, R. Budakian, S.J. Putterman, Phys. Rev. Lett. 76, 4528 (1996)
8. Mark-Tiele Westra, Patterns and Weak Turbulence in Surface Waves, Eindhoven (2001)
9. E. Falcon, С. Laroche, S. Fauve, Phys. Rev. Lett. 98, 094503 (2007)
10. С. Falcon, Е. Falcon, U. Bortolozzo, S. Fauve, Europhys. Lett. 86, 14002 (2009)
11. L.V. Abdurakhimov, M.Yu. Brazhnikov, A.A. Levchenko and L. P. Mezhov-Deglin, Developed Capillary Turbulence on the Surface of Normal and Superfluid 4He, J. Low Temp. Phys. 150, 426-430 (2008)
12. JI.B. Абдурахимов, М.Ю. Бражников, А. А. Левченко, Распределение вероятности осцилляций поверхности жидкого водорода в турбулентном режиме, Письма в ЖЭТФ 88, 21-26 (2008)
13. Л.В. Абдурахимов, М.Ю. Бражников, А.А. Левченко, Капиллярная турбулентность на поверхности нормального и сверхтекучего 4Не, Физика Низких Температур 35, 127-133 (2009)
14. Л.В. Абдурахимов, М.Ю. Бражников, А. А. Левченко, Эволюция турбулентного каскада на поверхности жидкого водорода при изменении спектральной характеристики возбуждающей силы, Письма в ЖЭТФ 89, 139-142 (2009)
15. L.V. Abdurakhimov, M.Yu. Brazhnikov, I.A. Remizov, A.A. Levchenko, Observation of wave energy accumulation in the turbulent spectrum ofcapillary waves on the He-II surface under harmonic pumping, Письма в ЖЭТФ 91, 291 (2010)
16. Д.М. Черникова, Физика низких температур 2, 1374 (1976)
17. Г.Н. Ватсон, Теория бесселевых функций, Иностранная литература, Москва (1949)
18. V. Zakharov, V. L'vov, G. Falkovich, Kolmogorov Spectra of Turbulence, vol. 1, Springer, Berlin (1992)
19. B.E. Захаров, H.H. Филоненко, Ж. Прикл. Мех. Техн. Физ. 8, 62 (1967)
20. Г. Е. Фалькович, А. Б. Шафаренко, ЖЭТФ 94, 172 (1988)
21. I.V. Ryzhenkova, G.E. Falkovich, Sov. Phys. JETP 71, 1085 (1990)
22. E. A. Kartashova, Physica D 46, 43 (1990)
23. A.N. Pushkarev, V.E. Zakharov, Physica D 135, 98 (2000)
24. E.A. Kartashova, Physica D 54, 125 (1991)
25. Л.Д. Ландау, E. M. Лифшиц, Теоретическая физика, т.6, Гидродинамика, Наука, Москва (1988)
26. Г. В. Колмаков, Письма в ЖЭТФ 83, 64 (2006)
27. R. J. Donnelly, C.F. Barenghi, J. Phys. Chem. Ref. Data 27, 1217 (1998)
28. В. Б. Шикин, Ю. П. Монарха, Двумерные заряженные системы в гелии, Наука,Москва (1989)
29. А. N. Pushkarev, V. Е. Zakharov, Phys. Rev. Lett. 76, 3320 (1996)