CC BY
24
ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ
УДК 532.593
Длинные волны в прямоугольном стенде, вызванные пластиной В. В. Трепачёв
(Донской государственный технический университет)
Проведено исследование формы поверхности тяжёлой жидкости в прямоугольном волновом стенде, покрытой длинными волнами над ровным дном. Волны вызваны горизонтальными гармоническими колебаниями волнопродуктора щитового типа. Получено точное решение краевой задачи с учётом диссипации энергии в жидкости. Результаты применимы для расчёта и лабораторного моделирования поверхностных и акустических волн.
Ключевые слова: тяжёлая несжимаемая жидкость, волнопродуктор, гармонические волны, длинные поверхностные волны, волновой стенд, диссипация энергии.
Введение. При рассмотрении теории длинных волн на поверхности тяжёлой жидкости можно в первом приближении пренебречь вертикальным ускорением частиц жидкости, это приводит к тому, что давление в жидкости практически равно статическому давлению, а скорость движения её частиц является функцией, которая не зависит от вертикальной координаты [1], [2]. Постановка задачи о волнах гармонического вида на мелкой воде в прямоугольном бассейне здесь сводится к рассмотрению краевой задачи в прямоугольнике для уравнения Гельмгольца [3], используемого в акустике. Траектории колебательного движения частиц жидкости представляют собой фактически отрезки прямых линий, параллельных дну бассейна. Наличие дна в таких условиях движения частиц вызывает дополнительное рассеяние энергии, учитываемое нами в рамках теории Рэлея
[1]—[5]. Данная работа дополняет результаты [1]—[5] и других исследований, т. к. в ней учитываются не только эффекты конечности всех размеров волнового стенда, но также эффекты диссипации энергии волн, наблюдаемые в ряде натурных экспериментов для волновых полей в океане [2].
Постановка задачи. Рассмотрим гармонические волновые колебания несжимаемой, тяжёлой жидкости под ровным дном в прямоугольном бассейне, имеющем длину I, ширину 2Ь и глубину
h. Опустим множитель е1ш1:, где ш — частота, t — время, i — мнимая единица. Постановка краевой задачи в линейном приближении теории длинных волн имеет вид
где ф(х, z) — потенциал скорости, к — волновое число, ^ — коэффициент трения о дно, д — ускорение свободного падения, с — скорость длинных волн, х — продольная координата, z — поперечная координата, и(г) — амплитуда колебаний пластины волнопродуктора на плоскости передней стенки х = 0.
д2ф/дх2 +д2ф/дг2 + к2ф = 0; к = ^(ш^^фшУдЛ ; ^ > 0; дф/дz = 0, г = ±Ь , 0 < х < I; ке к > 0; 1т к < 0; дф/дх = 0, |г| < Ь, х = I; I > 0, Ь > 0; дф/дх = и (г), |г| < Ь, х = 0; с =^[дїї,
(1)
(2)
(3)
(4)
Построение решения. Решение краевой задачи (1)—(4) отыскиваем в виде функционального ряда с неизвестными коэффициентами Am, но который удовлетворяет краевым условиям (2), (3) и имеет вид
ф (x, z) = 2 Am (e_2X-l/ • eiXmX + e-Xmx) cos [Ym (z - b)], (5)
m=0
где y m = mn/2b, а волновые числа xm продольных волн в (5) равны
Xm =>/k2 - Ym = ^(“7с2 - m2n/(2b)2) - /jw/с2 , c =Jgh . (6)
Подкоренное выражение в (6) имеет разрез на отрицательной части действительной оси.
Решение в виде (5) удовлетворяет также и волновому уравнению (1). Коэффициенты A находятся из краевого условия (4), которые представим соотношением
u (z) = ' 2 XmAm (e"2Xm"' - 1) C0s [Ym (Z - b)] , Ym = mn/2b . (7)
m=0
Система тригонометрических функций cos [Ym (z - b)] попарно ортогональна на отрезке |z| < b,
поэтому коэффициенты Am в случае m > 1 равны выражениям
b
Am = ('/(1 - ^ ) Xmb)\ U(Z) C0S Ym (Z - b)dz , Re (Xm )> ^ Im (Xm )< °. (8)
-b
Подставим (8) в (5), получим точное решение краевой задачи о длинных волнах в прямоугольном бассейне (1)—(4).
Форма свободной поверхности п (x, z) выражается через потенциал скорости
П (x, z) = ((-/“ + j)/g^(x,z). (9)
Коэффициенты Am имеют смысл амплитудных множителей поверхностных волн. Коэффициент A описывает амплитудный множитель «поршневой» волны. «Поршневая» волна имеет постоянное значение амплитуды вдоль направления оси z
b
Ao =(//(1 - e-2k) 2kb)\ u(z)dz . (10)
-b
Если u(-z) = -u(z), т. е. в случае нечётной по z амплитуды колебаний вибратора, то всегда амплитуда колебаний «поршневой» волны (10) имеет значение ноль. Можно подобрать вид чёт-
и I — U / \ U у— U
ной функции для амплитуды колебаний u (z), при которой амплитуда колебаний «поршневой»
волны равна нулю. Амплитуда колебаний «поршневой» волны отлична от нуля, если мгновенный расход жидкости q, вытесняемый вибратором в горизонтальном направлении, отличен от нуля
b
q = hj u(z)dz ф 0. (11)
-b
«Поршневая» волна распространяется по кратчайшему пути, соединяющему плоскость вибратора x = 0 и плоскость задней стенки бассейна x = l. Согласно (5) и ограничениям (8), построенное решение допускает предельный переход к бесконечной длине стенда l. В случае бесконечной длины стенда l в решении (5) присутствуют только волны, бегущие в одном направлении.
Бегущие волны в направлении от волнопродуктора к задней стенке стенда называют па-
дающими. Падающая волна определяется вторым слагаемым в круглой скобке формулы (5). В случае конечной длины волнового стенда l возникает отражённая волна от задней стенки бассейна, занимающей положение x = l. Амплитудный коэффициент отражённой волны описывает-
б
ся первым слагаемым в круглой скобке формулы (5). Он содержит параметр l. Предельный переход к бесконечной длине стенда l возможен при коэффициенте трения о дно j > 0.
Анализ волнового поля в стенде. Процесс возбуждения бегущих волн в волновом стенде возможен благодаря диссипации энергии. Рэлей предложил способ избавиться от паразитных свободных волн введением диссипативных сил, которые пропорциональны скоростям частиц жидкости [1]. В ряде работ, выполненных Вапняр Д. У., Михайловым Э. Н., Шапиро Н. Б., горизонтальный обмен количеством движения учитывается путём введения диссипативной массовой силы, которая пропорциональна горизонтальной скорости течения. Коэффициент трения j позволяет
[2] рассчитывать горизонтальный турбулентный обмен течений и внутренних волн в экваториальной зоне океана.
В монографии [2] приведены натурные оценки коэффициента трения j в пределах от
10-7 сек-1 до 10-5 сек-1. Известно, что в мелководных районах моря отмечена связь коэффициента
трения и угловой частоты колебаний. В лабораторных испытаниях волнопродукторов поверхност-
ных и внутренних волн, проводимых с участием автора данной статьи, отмечена высокая эффективность применения такой модели с использованием частотной зависимости коэффициента трения от масштаба волн. Отражённая волна имеет больший пробег, чем падающая волна, поэтому в волновом стенде достаточно большой длины создаются условия её подавления.
Пусть скорость движения пластины вибратора u (z) чётная функция следующего вида
u (z) = u cos yz , Y1 = 2b, (12)
где u(+b) = 0.
Используя (8), (10), находим коэффициенты разложения потенциала скорости (5):
uek u /X2ml Y
Am > ^A=-^ <13)
В отсутствии диссипации, т. е. в случае j = 0
k =^Т= С ' Xm =
4gh с \
“2 m2n2
“^2 , и = 0, (14)
с2 (2Ь)2
где к — действительное положительное число, называемое волновым числом «поршневой волны». «Поршневая волна» распространяется вдоль направления оси х. Появление волн, распространяющихся под острым углом к оси z , описывает условие излучения таких волн
. ш тп .. _
к = — >—, и = 0, т = 1,2,...,М , (15)
с 2Ь
где М — число волн, полученных под углом к оси х . Угол излучения ат равен
ат = агС5'п2Ьк' т = 1,2,...,М, и = 0. (16)
Чем больше ширина стенда 2Ь , тем большее количество волн М сможет излучаться под углом к оси х . Чем больше значение номера волны т, тем больше её угол наклона ат в (16). В случае к = тп/2Ь угол наклона становится наибольшим и принимает значение, равное п/2. Упругие волны, распространяющиеся в плоском слое, изучались в [3] применительно к проблеме использования слоя в качестве волновода акустических волн.
При отсутствии диссипации энергии ^ = 0 возможны резонансные явления, связанные с конечностью длины стенда I. Тригонометрическая функция, находящаяся в знаменателе коэффициента разложения потенциала скорости (13), обращается в ноль при выполнении условий
При частотах колебаний, определённых в (17) и (18), возникают резонансы с бесконечной амплитудой колебаний. В постановке задачи присутствует диссипация энергии ^ ф 0, которая не допускает возникновение резонансов с бесконечной амплитудой колебаний для распространяющихся волн, что дополняет известные результаты из теории плоских волн в волновом стенде [1].
Рассмотрим условие, которое противоположно по смыслу условию (15)
При выполнении (19) волны, распространяющиеся под углом к оси х, не излучаются. Соответственно резонансы типа (18) невозможны. При выполнении (19) колебания формы поверхности жидкости состоят из одной «поршневой волны» и счётного набора нераспространяющихся волн
[3]—[5]. На рисунке 1 изображена поршневая мода колебаний на оси волнового канала с учётом диссипации энергии.
Под нераспространяющейся волной понимается волна, которая и в отсутствии диссипации энергии имеет амплитуду, убывающую по экспоненциальному закону. Нераспространяющиеся волны локализуются в некоторой окрестности волнопродуктора. При этом они быстро затухают в направлении движения от волнопродуктора.
к, = — = іп /1, і = 1,2,..., ц = 0; ' с
(17)
— т2п2
~5-----------= І /1, І = 0,1,2,..., т = 1,2,...,М , ц = 0.
с 2Ь
(18)
(19)
-0,5
0,5
-1
0
1
1)
Рис. 1. Структура волнового поршневого поля на оси прямоугольного волнового стенда
Нераспространяющиеся волны могут также локализоваться и в окрестности отражающей стенки, что видно из формулы (5).
Заключение.
1. В работе показана возможность моделирования явления горизонтального турбулентного обмена для поля из падающих и отражённых поверхностных, длинных волн в прямоугольном стенде с произвольным соотношением сторон.
2. Показано, что волновое поле, состоящее практически только из бегущих волн, возникает в волновом стенде конечной длины благодаря влиянию диссипации энергии.
3. Получены новые условия отсутствия возбуждения поршневых, длинных волн горизонтальными гармоническими колебаниями щитового волнопродуктора.
Библиографический список
1. Сретенский, Л. Н. Теория волновых движений жидкости / Л. Н. Сретенский. — Москва: Наука, 1977. — 816 с.
2. Вапняр, Д. У. Планетарные волны и течения в экваториальной зоне океана / Д. У. Вапняр. — Киев: Наукова думка, 1976. — 222 с.
3. Лепендина, Л. Ф. Акустика / Л. Ф. Лепендина. — Москва: Высшая школа, 1978. — 448 с.
4. Трепачёв, В. В. Дифракция гармонической волны на вертикальном щите в канале / В. В. Трепачёв // Колебания и волны в жидкости и газе. — Горький: Изд-во ГПИ, 1990. — С. 89—95.
5. Трепачёв, В. В. Дифракция волн на пластине в узком канале / В. В. Трепачёв // Изв. вузов Сев.-Кавказ. региона. Естеств. науки. — 2002. — № 2. — С. 63—65.
Материал поступил в редакцию 25.05.2012.
References
1. Sretenskij, L. N. Teoriya volnovy'x dvizhenij zhidkosti / L. N. Sretenskij. — Moskva: Nauka, 1977. — 816 s. — In Russian.
2. Vapnyar, D. U. Planetarny'e volny' i techeniya v e'kvatorial'noj zone okeana / D. U. Vap-nyar. — Kiev: Naukova dumka, 1976. — 222 s. — In Russian.
3. Lependina, L. F. Akustika / L. F. Lependina. — Moskva: Vy'sshaya shkola, 1978. — 448 s. — In Russian.
4. Trepachyov, V. V. Difrakciya garmonicheskoj volny' na vertikal'nom shhite v kanale / V. V. Trepachyov // Kolebaniya i volny' v zhidkosti i gaze. — Gor'kij: Izd-vo GPI, 1990. — S. 89—95. — In Russian.
5. Trepachyov, V. V. Difrakciya voln na plastine v uzkom kanale / V. V. Trepachyov // Izv. vuzov Sev.-Kavkaz. regiona. Estestv. nauki. — 2002. — № 2. — S. 63—65. — In Russian.
LONG WAVES IN RECTANGULAR TANK GENERATED BY PLANE
V. V. Trepachev
(Don State Technical University)
The surface form of heavy liquid in the rectangular wave tank covered with long waves over the flat bottom is investigated. The waves are generated by the horizontal harmonic oscillations of the screen-formed wavemaker. The accurate solution to the boundary problem with regard to the energy dissipation in liquid is obtained. The results are applicable for the calculation and laboratory simulation of the surface and acoustic waves.
Keywords: heavy incompressible liquid, wavemaker, harmonic waves, long surface waves, wave tank, energy dissipation.
