На основе известных дисперсионных кривых для свободных колебаний неоднородной жидкости определена ее неоднородность. В метриках пространств C, L1, L2 найдена погрешность отыскиваемой неоднородности в зависимости от точности входной информации (точности измеряемых длин внутренних волн в неоднородной жидкости и точности задания их частот).
В океанологическом приближении рассматривается задача о свободных колебаниях неоднородной жидкости [1]:
(1)
где - вектор скорости в декартовой системе координат
, связанной с поверхностью Земли, p - гидродинамическое давление, ρ - плотность жидкости,
- ускорение свободного падения,
- сила Кориолиса,
- вектор угловой скорости вращения Земли, k - единичный орт, направленный по оси z (противоположно силе тяжести).
Граничные условия запишутся следующим образом:
условия непротекания на дне имеют вид
,
кинематическое граничное условие:
,
здесь - возвышение свободной поверхности,
динамическое граничное условие записывается в следующем виде
.
В линейном приближении решение задачи (1) ищется в виде бегущих волн. В результате этого рассматриваемая задача сводится к следующей задаче для амплитудной функции вертикальной скорости W:
(2)
В приближении Буссинеска с граничным условием типа «твердой крышки» (фильтрация поверхностных волн) в безразмерных величинах, после замены переменных
,
,
,
,
,
,
,
задача сведена к следующей (опущено обозначение безразмерности «~», заменено на z):
(3)
Будем считать, что неоднородность жидкости задана по линейному закону и значения ω таковы, что существует точка поворота. То есть, пусть на отрезке [0,1] в одной точке z0 функция
меняет знак, так что
,
,
,
.
Точка z0 называется точкой поворота, т.е. точка, в которой обращается в ноль переменный коэффициент в дифференциальном уравнении второго порядка системы (3).
Решение краевой задачи строим в явном виде:
,
;
,
;
,
,
.
Удовлетворяя граничным условиям, получаем частотное уравнение
.
Считая , заменим в этом уравнении функции Бесселя их асимптотическим представлением [2]
, при
,
и учтем, что при больших z справедливо .
Частотное уравнение примет вид
.
Отсюда выводим
.
Так как , где
, получим
,
Для линейного профиля скоростей интеграл вычисляется в явном виде, на основе которого при заданных значениях ,
, построены дисперсионные кривые.
Теперь ставится задача: считая, что известны значения ω и k с дисперсионных кривых, определить параметры неоднородности жидкости ,
.
Для их определения выписываем уравнение для двух пар (ω, k)
,
,
что приводит к нелинейной системе для определения параметров ,
.
Нелинейность системы обуславливает неоднозначность определения параметров ,
. Для однозначности определения этих параметров выписывается «избыточное» число уравнений для «избыточных» пар чисел ω и k, и параметры
,
определяются как минимум функционала
.
Изучается точность восстановленной функции неоднородности жидкости:
в пространстве C
;
в пространстве L1
в пространстве L2
Здесь - восстановленная функция стратификации по найденным значениям параметров
,
, а функция
- заданная стратификация жидкости.
При этом исследуется вопрос влияния на единственность и точность восстановления неоднородности жидкости тот факт, берутся ли пары чисел с одной дисперсионной кривой, или с разных. Также исследуется влияние близости между собой точек, взятых с дисперсионных кривых, на точность восстановления профиля стратификации.
СПИСОК ЛИТЕРАТУРЫ:
- Миропольский Ю.З. Динамика внутренних гравитационных волн в океане. - Л.: Гидрометеоиздат, 1981. - 301 с.
- Градштейн И.С., Рыжик И.М. Таблицы интегралов, сумм, рядов и произведений. - М.: Наука, 1971. - 1108 с.
Библиографическая ссылка
Козьменко Ю.Г., Потетюнко Э.Н. ОБРАТНАЯ СПЕКТРАЛЬНАЯ ЗАДАЧА ОПРЕДЕЛЕНИЯ ЛИНЕЙНОГО ПРОФИЛЯ СТРАТИФИКАЦИИ ЖИДКОСТИ // Успехи современного естествознания. 2010. № 5. С. 105-108;URL: https://natural-sciences.ru/ru/article/view?id=8142 (дата обращения: 02.04.2025).