Математическое моделирование волнового движения воды в узком глубоком непризматическом водохранилище с учетом упругости воды
Асп. Музаев Н.И. Кафедра математики. Северо-Кавказский горно-металлургический институт (государственный технологический университет Составлена математическая модель волнового движения воды в узком глубоком непризматическом водохранилище с учетом упругости воды. Модель представляет начально-краевую задачу математической физики для потенциала средней по ширине векторной скорости. В основном дифференциальном уравнении начально-краевой задачи в качестве переменных коэффициентов содержится ширина водохранилища, зависящая от продольной и вертикальной координат. Составленная математическая модель позволяет решить широкий класс прикладных задач, связанных с теорией колебаний и волн в узких глубоких непризматических водохранилищах. Предположим, что в прямоугольной системе координат xoyz часть пространства, ограниченная условиями 0 £ x £ l, – 1/2 B(x, z) £ y £ 1/2 B(x, z), –H £ z £ 0, представляет узкое глубокое непризматическое водохранилище. Ось oz направлена вертикально вверх, ось ox направлена в продольном, а ось oy – в поперечном направлении водохранилища. L – длина, B(x,z) – ширина, H – глубина водохранилища. Как правило, в горных условиях водохранилища строятся в узких глубоких каньонах ущелий рек. В связи с этим в дальнейшем будем считать, что ширина водохранилища B(x, z) намного меньше, чем ее длина. Кроме этого будем считать, что градиенты в поперечном направлении поля скоростей и гидродинамического давления намного меньше, чем градиенты в продольном и вертикальном направлении водохранилища. Ширина схематизированного водохранилища зависит от продольной и вертикальной координат B = B(x, z), т.е. рассматривается водохранилище с непризматической конфигурацией как в продольном, так и в вертикальном направлении. Для таких водохранилищ решение пространственной задачи волнового движения воды связано с большими математическими трудностями и в мире никем не решена. В связи с этим трехмерные дифференциальные уравнения гидродинамики интегрально усредняют по площади живого сечения воды, в результате получают одномерные дифференциальные уравнения движения воды в естественных водоемах. В связи с тем, что водохранилища в горных местностях являются глубокими и узкими, то, в отличие от теоретической гидравлики, трехмерные уравнения гидродинамики мы интегрально усредняем только по поперечной координате y, а вертикальную координату оставляем без изменений. В гидродинамике волнового движения жидкости дифференциальные уравнения используют в «отфильтрованном» виде, т.е. пренебрегают нелинейные члены как малые величины по сравнению с линейными членами. В проекциях на оси x, y и z эта система в «отфильтрованном» виде запишется так [1-3]: ; ; ; (1) , где приняты следующие обозначения: Vx , Vy и Vz – скорости в продольном, поперечном и вертикальном направлениях соответственно, зависящие от всех пространственных координат и времени t ; r – плотность; P – гидродинамическое давление; a – скорость звука в воде. Усредним интегрально систему дифференциальных уравнений (1) по поперечной координате y. ; ; . (2) . Обратимся к известной формуле дифференцирования под знаком интеграла: . (3) Интегралы, входящие в выражения (2), преобразуются так: ; . (4) В результате такого усреднения система~(2) запишется следующим образом: ; (5) ; (6) , (7) где приняты обозначения: , , . (8) Величины Ux , Uz и P представляют собой средние значения по ширине водохранилища соответственно Ux , Uz и P; q(x,z,t) – интенсивность боковой приточности, определяющаяся выражением: (9) Систему (5,6) в векторной форме можно записать так: , (10) , (11) где . Считая, что движение воды безвихревое, т.е. rot = 0, и вводя потенциал средней по ширине скорости , (12) из выражения (10) получаем интеграл Коши в линейном приближении: . (13) Компоненты средней скорости через потенциал скорости F(x, z, t) выражаются так: , . (14) В связи с тем, что потенциал скорости волнового движения жидкости определяется с точностью до произвольной функции, зависящей только от времени t, произвольную функцию f(t) можно считать тождественно равной нулю. На свободной волновой поверхности должно быть задано гидродинамическое давление . При отсутствии внешнего давления . Обозначив уравнения волновой поверхности через z = h(x, t), выражение (13) запишется так: . (15) Линеаризуя выражение (15), получаем: . (16) В линейном приближении очевидно равенство: . (17) Дифференцируя выражение (16) по t и подставляя в него (17), получаем: . (18) Из выражения (13) при f(t) = 0 для давления получается следующая его зависимость от потенциала скорости: . (19) Подставив выражения (14) и (19) в (11), получим следующее дифференциальное уравнение для потенциала скорости: . (20) Как известно, в классической теории двумерного волнового движения упругой жидкости, для потенциала скорости имеется следующее уравнение [1,3]: . (21) Сравнивая уравнения (20) и (21), легко заметить, что в полученном в данной работе уравнении дополнительно содержатся три слагаемых. Последние две слагаемые в левой части учитывают непризматическое очертание водохранилища как в плане, так и по глубине. Величина q(x, z, t) представляет интенсивность вытеснения воды обвально-оползневым массивом либо интенсивность вторжения селелавинообразного потока в водохранилище. Отметим, что в статье [4] получено дифференциальное уравнение для потенциала волнового движения несжимаемой жидкости в непризматическом водохранилище. В данной работе теория представляется более общей в связи с тем, что в ней учтена упругость воды, т.е. первое слагаемое уравнения(20). Список литературы 1. Ламб Г. Гидродинамика. М.: Гостехиздат, 1947. 2. Стокер Дж. Дж. Волны на воде. М.: Изд-во иностранной литературы, 1959. 3. Сретенский Л. Н. Теория волновых движений жидкости. М.: Наука, 1977. 4. Музаев И. Д., Созанов В. Г. К теории поверхностных гравитационных волн Коши – Пуассона в узких глубоких непризматических водоемах// Изв. вузов. Сев.-Кав. регион. Сер. ест. науки. Ростов-на Дону. 1995. № 3. Для подготовки данной работы были использованы материалы с сайта http://www.skgtu.ru/
Популярное: Почему человек чувствует себя несчастным?: Для начала определим, что такое несчастье. Несчастьем мы будем считать психологическое состояние... Как построить свою речь (словесное оформление):
При подготовке публичного выступления перед оратором возникает вопрос, как лучше словесно оформить свою... Как распознать напряжение: Говоря о мышечном напряжении, мы в первую очередь имеем в виду мускулы, прикрепленные к костям ... ©2015-2024 megaobuchalka.ru Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав. (249)
|
Почему 1285321 студент выбрали МегаОбучалку... Система поиска информации Мобильная версия сайта Удобная навигация Нет шокирующей рекламы |