На основе известных дисперсионных кривых для свободных колебаний неоднородной жидкости определена ее неоднородность. В метриках пространств 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 (дата обращения: 14.12.2024).