ON MODELLING OF PERIODIC SURFACE WAVES OF AN IDEAL HEAVY INCOMPRESSIBLE FLUID

Research article
DOI:
https://doi.org/10.23670/IRJ.2023.138.199
Issue: № 12 (138), 2023
Suggested:
24.07.2023
Accepted:
22.11.2023
Published:
18.12.2023
371
19
XML
PDF

Abstract

An algorithm for modelling periodic surface waves of an ideal heavy incompressible fluid is presented, which, unlike most of the known ones, does not use hypotheses about the form of the potential function specifying the laminar flow. The analysis shows that the studied laminar flow is not a potential flow. The algorithm is based on the study of cyclic motion of a fluid particle within the laws of classical mechanics. The equations of motion using functions of velocity and pressure gradient, periodic in time, are obtained. It is shown that the implementation of these equations is equivalent to the execution of the Euler equation and the continuity equation. The model admits the existence of large amplitude waves embedded in the flow of periodic waves, the so-called rogue waves.

1. Введение

Моделирование периодических поверхностных волн идеальной тяжелой несжимаемой жидкости, которыми могут быть морские волны на глубокой воде, представляет трудность, так как нужно обеспечить одновременное удовлетворение нескольких условий. А именно, процесс должен быть периодическим, должны удовлетворяться уравнение Эйлера, являющееся следствием законов классической механики, и уравнение непрерывности, которое соответствует равенству притока и оттока несжимаемой жидкости в каждый момент времени для каждой элементарной области потока. Большинство известных моделей основано на постулировании предположений, например, что движение частицы жидкости является комбинацией кругового движения и поступательного, или о том, что поток является потенциальным, то есть, скорость в точках потока равна производной некоторой функции

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

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

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

2. Постановка задачи

Для элемента единичного объема идеальной несжимаемой жидкости плотности ρ (или, что то же самое, элемента идеальной несжимаемой жидкости массы ρ) ламинарного потока справедливо равенство

:

img
(1)

где img – вектор скорости элемента единичного объема жидкости плотности ρ, img – время, img – градиент давления img.

При действии силы тяжести возникает дополнительное слагаемое img, модуль вектора img – ускорение свободного падения img, то есть img:

img
(2)
Локальные системы координат

Рисунок 1 - Локальные системы координат

Так как в (2) входит производная скорости, то постоянная ее составляющая не влияет на равенство. Для периодической волны такой составляющей является перемещение контура волны в целом, то есть фазовая скорость. Рассмотрим траекторию движения частицы жидкости единичного объема (или, что то же самое, – массы img) без учета фазовой скорости (рис. 1), для краткости ниже будем называть ее «траекторией». Система координат наблюдателя img (рис. 1) перемещается с постоянной скоростью, равной фазовой скорости волны. 

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

img
(3)

Где угол img – угол с осью img от противоположного направления вертикально вниз направленной силы тяжести, измеряемый по часовой стрелке в точках отрицательной кривизны и против часовой в точках положительной кривизны. Угол img – угол вектора градиента давления img с осью img. img – знакопеременный радиус кривизны в рассматриваемой точке, того же знака, что знак кривизны

.

В левой части второго равенства (3) проекцию ускорения на ось img, перпендикулярную траектории, мы приравняли центростремительному ускорению, соответствующему радиусу кривизны img 

, равному img. То есть,

img

где img – единичные векторы вдоль осей img локальной системы координат (рис. 1) соответственно.

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

Так как волны периодические, и каждой точке траектории в рамках периода соответствует единственное значение img, img, то входящие в (3) величины – периодические функции img, причем, в диапазон img «укладывается» один период, то есть период меньше или равен img. Так как вторые слагаемые (3) имеют период img, то остальные слагаемые также имеют период img. То есть, имеют представление

:

img
(4)

где img – константы.

Из (3) очевидно, что в представлении (4) вторых слагаемых правых частей (3) параметр img, то есть для них оно совпадает с (5).

img
(5)

где img – константы. Из постулированного нами вида геометрии волны (рис. 1) (то есть вида граничной поверхности ламинарного потока) следует, что периодическая с периодом img функция img имеет в представлении (4) параметр img (5). Поэтому центростремительное ускорение img, равное левой части второго равенства (3), является периодической с периодом img функцией, если скорость img равна (6).

img
(6)

где img – константа img, а функция img имеет представление (5).

Действительно, для того, чтобы центростремительное ускорение img было конечно, необходимо, чтобы функция img была равна нулю в тех же точках img, что функция img:

img
(7)

Выполняя требование (7), почленно сравнивая представление (5) функции img с представлением (4) функции img при img, img, img, получим (6). В соответствии с изложенным ниже случай img рассмотрен как предельный от случая img.

Из двух предыдущих абзацев следует, что первое слагаемое правой части второго уравнения (3) также должно иметь представление (5). Это возможно, если

img

img
(8)

или

img
(9)

или

img
(10)

Рассмотрим правую часть первого равенства (3) варианта (8):

img
(11)

Из (11), (6) следует, что вариант (8) для выполнения первого равенства (3) невозможен: производная по img от img вида (5) не может иметь постоянного слагаемого правой части (11).

При вариантах (9), (10) правая часть первого равенства (3) – сумма тригонометрических функций img с постоянными множителями. Следовательно, левая часть первого равенства (3) также сумма тригонометрических функций img с постоянными множителями. Следовательно, так как img вида (5), то img – константа, то есть

img
(12)

где img – константа, img.

Если справедливо (9), то 

img
(13)

где img – константы, такие, что img, img, то есть

img
(14)

или

img
(15)

В вариантах (13), (14) или (15), давление img на оси img по глубине не меняется. Получаем противоречие с физическим смыслом. Поэтому для поверхностной траектории физически единственно возможным является вариант (10) когда уравнения (3) преобразуются к виду

img
(16)

где

img
(17)

img – константы,

img
(18)

Кроме равенств (17), (18) касательная скорость img и радиус img связаны величиной угловой скорости, то есть

img
(19)

где img – угловая скорость. С учетом принятых знаков img и img из (6), (19) получим

img
(20)

Так как img – константа, от img не зависит, из (19), (20) следует

img
(21)

Из первого равенства (16), (17), (21) следует

img
(22)

Таким образом

img
(23)
img
(24)

Из второго равенства (16), (23), (24) следует

img
(25)

то есть градиент давления направлен по нормали к траектории и является константой (когда траектория совпадает с поверхностью волны, давление img на ней равно атмосферному принятому постоянным).

Полученные представления справедливы не только при непрерывном изменении аргумента img в диапазоне: img, но постулированная нами геометрия (рис. 1) допускает наличие разрывов аргумента img в точках нулевого радиуса кривизны img. В этом варианте в точке img, img, функция img имеет скачек значений: предельное значение при стремлении img слева, img, равно img, а при стремлении img справа равно img. На длине половины волны каждая из функций (23), (24), (25) «распадается» на две функции с областями аргумента: img, img, заданных в местных системах координат (рис. 1), в которых при img и img значения радиуса кривизны img равны нулю (24), что соответствует равным по модулю и противоположным по знаку значениям параметра img. То есть, эти две функции «сшиты» в точках равного нулю радиуса img, где задаваемые ими участки траектории имеют совпадающую касательную и косую симметрию относительно точки img. При этом градиент давления (25) остается постоянным на всей «сшитой» траектории, так как изменение его знака происходит в местной системе координат. С помощью косой симметрии относительно точки img определяется половина волны, в свою очередь с помощью симметрии относительно точки img определяется целая волна (поэтому случай img в (7) соответствует пределу img).

Так как 

img
(26)

где img – период, то, задавая img и img, мы определяем траекторию. Действительно, зная img, определяем img из (26) и параметр img из (27).

img
(27)

Эти параметры определяют четверть траектории и, следовательно, всю траекторию через описанную выше симметрию.

Координаты точек поверхности волны img (рис. 1), соответствующие четверти волны: img, задаются выражениями:

img
(28)
img
(29)

imgimg

img

Высота траектории img равна 

img
(30)

длина траектории img равна 

img
(31)

Из (30), (31) следует, что при заданном значении img наибольшая амплитуда при img. В таблице 1 приведены значения максимальных длины и высоты волны в зависимости от периода img, что соответствует img img. Так как амплитуда траекторий одной длины волны img с глубиной уменьшается, то соответствующее значение img, img, уменьшается. Следовательно, заданные значения img определяют не только поверхностную, но и все последующие траектории по глубине ламинарного потока.

Действительно, длина волны img и значение img поверхностной траектории определяются (26), (27), (31), уменьшая значение img от максимального на поверхности до нуля, можно определить при постоянной img значение img каждой последующей траектории из (31). Так как при img значение скобок (31) стремится к нулю, то img, модули радиусов (24) стремятся к бесконечности.
Определение расстояния h между траекториями

Рисунок 2 - Определение расстояния h между траекториями

Определим положение каждой траектории по глубине. Рассмотрим сегмент img (рис. 2) с осью симметрии img на которой img (рис. 1). Ламинарный поток несжимаемой жидкости, скорость которого связана с радиусом пропорцией: img, (23), (24), может существовать вблизи img только в рамках сегмента img (рис. 2).

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

img

пропорционален радиусу дуги с центром img (рис. 2), и эти пути не пересекаются (21). Из этого условия определяется расстояние img.

Таблица 1 - Максимальные длина и высота волны в зависимости от периода T

Первые 16 траекторий волны периода 8 сек и длиной L = 49,94 м с параметрами – таблица 2

Рисунок 3 - Первые 16 траекторий волны периода 8 сек и длиной L = 49,94 м с параметрами – таблица 2

Таблица 2 - Параметры траекторий – рисунок 3

В предположении потенциальности потока проекции img скорости на оси img ортогональной системы координат, связанной с Землей, где ось img параллельна оси img, а ось img параллельна оси img (рис. 1), должны быть производными от некоторой гармонической функции img (32)
.
img
(32)

В точках оси img проекция img равна нулю: img, так как эта ось перпендикулярна траекториям (рис. 1) и фазовой скорости. Следовательно, в этих точках (32):

img
(33)

где img – константа, значения скорости постоянны и равны img, то есть img, где img – фазовая скорость равная константе. Очевидно, что построенное решение условиям (32), (33) не удовлетворяет, поэтому моделируемый поток не является потенциальным.

Для фазовой скорости img должно выполняться условие

img
(34)

где img – максимальное по модулю значение орбитальной скорости img (23) в точках оси img (рис.1), равное значению в самой высокой точке волны.

Покажем. что предложенная модель удовлетворяет уравнению Эйлера

и уравнению непрерывности
. Уравнение Эйлера и уравнение непрерывности сформулированы в глобальных, связанных с Землей координатах, относительно которых движется волна. Уравнение Эйлера соответствует выполнению законов классической механики, уравнение непрерывности – уравнение равенства притока и оттока несжимаемой жидкости в рассматриваемый элементарный объем.

Уравнения (16) – уравнения выполнения законов классической механики в местной (движущейся вместе с единичным объемом) системе координат. То есть, уравнения (16) заменяют собой (обеспечивают выполнение) уравнение Эйлера, так как уравнение Эйлера – следствие законов классической механики в глобальных, связанных с Землей координатах, а (16) – выполнение тех же законов в локальных координатах, перемещающихся вместе с частицей жидкости единичного объема.

При условии (34) кривую, четверть которой задана координатами (28), (29), а половина и целая описанной выше симметрией, можно рассматривать как траекторию перемещения частицы жидкости единичного объема, массы ρ в рамках одного цикла периодической волны. В то же время ее можно рассматривать как геометрическое место точек частиц жидкости единичного объема массы ρ в заданный момент времени img. При этом в каждой точке кривой в каждый момент времени приток жидкости равен оттоку.

Так как аналогичную кривую с той же длиной волны img можно провести через любую точку моделируемого потока, то бесконечный набор таких кривых, связанных соотношением радиусов img (рис. 2), задает геометрическое место точек частиц жидкости в заданный момент времени img одного цикла волны в целом. Где в каждой точке в каждый момент времени приток жидкости равен оттоку, следовательно, выполняется уравнение непрерывности

.

В качестве примера на (рис. 3) приведены первые 16 траекторий волны с периодом 8 сек и длиной L = 49,94 м на половине волны.

Значения таблицы 1 получены из условия img. Из того, что наблюдаемых высот волн больше значений третьего столбца таблицы 1 не зафиксировано, следует, что для волн большего периода и длины физически реализуется вариант поверхностной траектории при img.

3. О «волнах-убийцах»

Есть наблюдения одиночных волн большой амплитуды, «встроенных» в периодический волновой процесс «обычной» амплитуды, о так называемых «волнах-убийцах»

 (аномальные волны могут объединяться в серии, тогда «волн-убийц» будет несколько). Они имеют большое практическое значение, так как, по мнению ряда авторов, являются основной причиной гибели кораблей. Эти «волны-убийцы» могут быть в виде одиночного гребня или одиночной впадины. Покажем, что предлагаемая модель допускает существование таких гребней или впадин.

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

Действительно, непрерывность функции img соответствует непрерывности радиуса кривизны траектории в точке img, img. Этого достаточно для «склейки» гребня одной волны с бесконечно дифференцируемой траекторией со впадиной другой волны с бесконечно дифференцируемой траекторией, при этом уравнения (16) выполняются. Сила тяжести – бесконечно дифференцируемая функция времени, поэтому левая часть первого уравнения (16) также должна быть бесконечно дифференцируемой, как и правая. Однако, не учтенные в модели физические факторы, в частности воздействие ветра или разницы температур на разных высотах, могут привести к однократному или постепенному росту модуля разности величины потенциальной энергии в точках img и img, или модуля разности величины потенциальной энергии в точках img и img. Так как градиент давления перпендикулярен траектории, работу в рамках модели совершает только сила тяжести, следовательно, этот модуль будет равен произведению ускорения свободного падения img на разницу высот в точках img и img, или в точках img и img. То есть, приток энергии будет трансформироваться в увеличение высоты гребня или впадины. Этот аномальный гребень или впадина будет продолжать существовать «встроенным» в поток периодических волн «обычной» амплитуды, до тех пор пока бесконечная дифференцируемость по времени слагаемых уравнений (16) не будет восстановлена.

4. Заключение

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

Article metrics

Views:371
Downloads:19
Views
Total:
Views:371