УДК 378.4
А. А. Зайцев, С. Н. Иванов, А. В. Юдина ВЗАИМОДЕЙСТВИЕ АНТИЦИКЛОНА С БАРЬЕРОМ
Рассмотрено взаимодействие двумерного антициклона с барьером. Получено и решено его уравнение движения. Установлено, что при приближении к краю барьера антициклон ускоряется, создавая интенсивное течение. Изучена картина линий тока.
The article considers the interaction of a two-dimensional anticyclone with a barrier and generates and solves its motion equation. The authors claim that an anticyclone accelerates at the approach of the edge of a barrier forming an intensive stream. The article also offers the streamline pattern analysis.
Ключевые слова: двумерный вихрь, антициклон, взаимодействие, траектория, скорость течения, линии тока.
Keywords: two-dimensional vortex, anticyclone, interaction, trajectory, flow velocity, streamline.
Введение
Вопрос о взаимодействии вихрей с морскими берегами представляет несомненный практический интерес, однако детально разработанная теория этого взаимодействия отсутствует. В классической монографии [1] имеется следующее замечание. Взаимодействие циклона с антициклоном, модули интенсивностей которых одинаковы (вихревой дублет), ведет к поступательному движению этой пары с постоянной скоростью, обратно пропорциональной расстоянию между вихрями. При этом средняя линия между ними оказывается линией тока. Поскольку береговая линия всегда является линией тока, то отсюда можно сделать такой вывод. Течение, создаваемое одиночным вихрем вблизи берега, отражается от береговой линии. Отраженное течение можно рассматривать как течение, создаваемое мнимым вихрем противоположной интенсивности, который расположен в симметричной точке. Вихрь попадает в отраженное течение и начинает двигаться параллельно берегу с постоянной скоростью, обратно пропорциональной расстоянию до берега. Для берегов, имеющих сложную конфигурацию, задачи взаимодействия решить значительно трудней. В этой статье мы решаем ее для случая взаимодействия двумерного антициклона с барьером.
Постановка задачи. Уравнения движения
Вестник Российского государственного университета им. И. Канта. 2009. Вып. 4. С. 17—20.
Здесь рассматривается двумерная задача. Используются комплексные переменные: z и Zo обозначают соответственно комплексные координаты точки наблюдения и антициклона (случай циклона интереса не представляет, поскольку он уходит от кромки барьера), интенсивность которого равна -к, к>0. Допустим, что барьером служит луч 0<Rez<да, Imz=0 (рис. 1), а для начальной координаты антициклона выполнено условие Imz0>0. Требуется найти уравнение движения вихря и определить его траекторию.
Поле скоростей течения, создаваемого вихрем, определяется с помощью конформного отображения £= Г на верхнюю полуплоскость 1т^>0 и имеет вид [2; 3]
Рис. 1. Антициклон вблизи барьера
V=-
4П
1 1
ч г-^|z0z г + ^[z0z у
(1)
Вычитая собственное течение вихря, получаем поле скоростей того те чения, которое воздействует на вихрь:
( \
1 1
Vr=-
к
4П
(2)
Скорость вихря должна совпадать со скоростью V в точке z=z0. Переходя в формуле (2) к пределу при z—^о, получаем
к
4п
V 2г0 - | |
(3)
Получено уравнение движения антициклона.
Траектория антициклона
Уравнение траектории антициклона получается интегрированием уравнения (3). С этой целью полагаем z0=rexp(iф) и делаем подстановку в уравнение (3). После отделения действительной и мнимой частей получаем
к Г -- к (4)
Г =-
4пг2
Система (4) имеет первый интеграл
• V л гет— =а,
(5)
который и будет уравнением траектории. Из соотношения (5) следует, что траектория вихря симметрична относительно оси абсцисс.
С помощью равенства (5) несложно найти зависимость г от t. Действительно, с учетом (5) первое уравнение системы (4) принимает вид
1
1
. кл]г2 - ё2
г=—-г-,— ^
Его интегрирование дает нужную зависимость г=л/к 2t2 + 64п2ё4 /8пй (мы считаем, что вихрь находился в окрестности выступа в момент £=0). Она показывает, что с ростом £ от -ж до 0 расстояние до выступа убывает и в момент пересечения оси абсцисс достигает минимального значения й. Можно показать, что в том случае, когда вихрь находится вдали от выступа, расстояние до барьера равно 2й.
Формулы (4) и (6) позволяют найти зависимость модуля скорости
вихря от расстояния; она следующая: | V |=кл/г2 + 3ё4 /8пйг. Это равенство показывает, что вдали от кромки барьера скорость вихря минимальная и равна | V | тп=к/8лй, а вблизи максимальная и равна | V | тах=2| V | тап=к/4пй.
Поле течений и линии тока
Выражение для скорости течения, создаваемого движущимся вихрем, можно получить, если использовать формулу (1), равенство z0=rexp(iф) и зависимость г и ф от £, в том числе зависимость (6). Наибольший интерес представляет оценка модуля скорости течения вблизи кромки барьера.
Полагая в формуле (1) R=|z|<<r, получаем | V |=Ып(ф/2)/2л-\/гЯ , то есть модуль скорости течения растет обратно пропорционально корню квадратному из расстояния до кромки барьера. Это означает, что вблизи барьера скорость течения очень велика.
Представляет интерес картина линий тока. Ограничимся случаем, когда вихрь проходит в окрестности края барьера. Используя формулу (1), получаем следующее выражение для функции тока:
к
т=-------1п
2п
V? - 1*[ё
V? + ^л[ё
Вдоль линий тока эта величина является константой; следовательно, аргумент логарифма также будет константой. Переходя к полярной системе координат, получаем уравнение линий тока в виде
г - 2у[ёгзш((р/2) + ё _
г + 2^/dГ" зш((р/2) + ё
=а, 0<я<1.
Анализ этого уравнения показывает: 1) через каждую точку водного пространства проходит единственная линия тока, 2) все линии тока являются гладкими замкнутыми линиями, симметричными относительно оси абсцисс, 3) для каждого значения угла ф0 в промежутке 0<ф0<п линия тока касается границ сектора ф0<ф<2л-ф0 в двух симметричных точках; в случае ф0>2п/3 и только в этом случае линия тока будет выпуклым овалом.
Заключение
В работе установлено, что взаимодействие вихря с барьером приводит к его движению, которое ускоряется по мере приближения к краю барьера. Расстояние до барьера сокращается. При этом возникает интенсивное течение, модуль скорости которого обратно пропорционален корню квадратному из расстояния до края барьера, что может привести к неблагоприятным последствиям, хорошо известным на практике. Изучена картина линий тока.
Для решения поставленной задачи использованы методы теории функций комплексной переменной и теории динамических систем. Их сочетание должно позволить получить решение более сложных задач о взаимодействии вихрей.
Список литературы
1. Милн-Томсон Л. М. Теоретическая гидродинамика. М., 1964.
2. Лаврентьев М. А., Шабат Б. В. Методы теории функции комплексной переменной. М., 1973.
3. Лаврентьев М. А., Шабат Б. В. Проблемы гидродинамики и их математические модели. М., 1973.
Об авторах
А. А. Зайцев — канд. физ.-мат. наук, ст. науч. сотр., РГУ им. И. Канта.
С. Н. Иванов — канд. физ.-мат. наук, зав. лаб. морской метеорологии АО ИОРАН.
А. В. Юдина — студ., РГУ им. И. Канта.
Authors
A. Zaytsev — Dr., IKSUR.
S. Ivanov — Dr., Oceanology Institute of the Russian Academy of Sciences.
A. Yudina — student, IKSUR.