На правах рукописи

ВАСИН Андрей Васильевич

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ОТРЫВНЫХ

ТЕЧЕНИЙ И ИХ ВОЗДЕЙСТВИЙ НА ГИДРОТЕХНИЧЕСКИЕ

СООРУЖЕНИЯ

Специальность 05.13.18 – «Математическое моделирование, численные

методы и комплексы программ»

Автореферат диссертации на соискание ученой степени доктора технических наук

САНКТ–ПЕТЕРБУРГ

2013

Работа выполнена в ФГБОУ ВПО «Государственный университет морского и речного флота имени адмирала С.О. Макарова»

Научный консультант:

доктор технических наук, профессор Голоскоков Дмитрий Петрович

Официальные оппоненты:

доктор физико-математических наук, профессор, ФГБОУ ВПО «СанктПетербургский государственный университет гражданской авиации», профессор кафедры прикладной математики Береславский Эдуард Наумович;

доктор технических наук, профессор, ФГБОУ ВПО «ГУМРФ имени адмирала С.О. Макарова», профессор кафедры гидротехнических сооружений, конструкций и гидравлики Гапеев Анатолий Михайлович;

доктор физико-математических наук, профессор, ФГБОУ ВПО «СанктПетербургский государственный университет», заведующий кафедрой математического анализа Широков Николай Алексеевич.

Ведущая организация: ФГБОУ ВПО «Санкт-Петербургский государственный морской технический университет».

Защита диссертации состоится «26» декабря 2013 г. в 1400 часов на заседании диссертационного совета Д 223.009.03 в ФГБОУ ВПО «ГУМРФ имени адмирала С.О. Макарова» по адресу: 198035, Санкт-Петербург, ул. Двинская, д.5/7, к. 235.

С диссертацией можно ознакомиться в библиотеке ФГБОУ ВПО «ГУМРФ имени адмирала С.О. Макарова».

Автореферат разослан «24» сентября 2013 г.

Ученый секретарь диссертационного совета, кандидат технических наук, доцент Е.Г. Барщевский

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность проблемы. Проблема разработки правильных и удобных методов расчета воздействия потоков жидкости на гидротехнические сооружения возникает постоянно. Особый интерес представляют водоводы с поворотами, поскольку, именно, каналы такой формы присутствуют на реальных гидротехнических сооружениях. В поворотных отделах возникают ввиду вязкости вихревые или даже турбулентные течения, где развитые методы потенциального течения идеальной жидкости не применимы. Исследование же вязкостных моделей жидкости с использованием уравнений Навье-Стокса в настоящее время далеко от завершения, и носит полуэкспериментальный характер с определением большого числа опытных коэффициентов и разбиением сложных областей на более простые. Инженеры – проектировщики поставлены в трудное положение. С одной стороны достаточно давно выведены уравнения Навье-Стокса, с другой стороны аналитическому решению поддаются исключительные случаи, входящие в «золотой фонд гидродинамики» [1]. В расчетах гидравлических характеристик при процессах наполнения и опорожнения шлюзов исходили из основного уравнения одномерного потока жидкости. При этом галерею делили на участки и приближенно считали скорость на каждом участке постоянной и равной средней скорости на участке. Главное требование для применения одномерного уравнения движения о выравнивании средних скоростей основывалось на лабораторных исследованиях распределения скоростей в потоке по поперечному сечению. Одномерный подход не мог дать кинематических картин движения жидкости, оставались неисследованными зоны отрыва исследовалось влияние распределения скоростей на элементы конструкций гидротехнических сооружений. По этой причине разработка приближенных моделей и вычислительных методов, альтернативных решению исходных трехмерных уравнений гидродинамики, является актуальной проблемой.

Заметим, что в настоящее время при наличии уравнений Навье-Стокса, расчет напряженно-деформированного состояния затвора производится в основном на гидростатическую нагрузку, что совершенно недопустимо.

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

Степень разработанности. Хорошее согласование с опытом в областях простой формы дает модель Лаврентьева-Шабата [2] склеивания вихревых и потенциальных течений. Известно математическое развитие данной задачи, но приложения математической модели Гольдштика ограничиваются модельными формами камер и галерей. Для применения данной модели к течению жидкости в реальных водопропускных сооружениях с поворотами и препятствиями, для получения результатов, полезных для проектирования и эксплуатации ГТС, требуется построение модификаций модели Лаврентьева-Шабата, которые отсутствуют для областей произвольной формы.

Кроме того большие трудности представляет разработка численных методов для решения задачи Гольдштика, поскольку в областях сложной формы нет других методов получения числовых результатов, отличных от численных методов. Многие исследователи отдают предпочтение методу конечных элементов, имеются коммерческие продукты ANSIS и Flow3d, разработчики которых утверждают, что это абсолютно автоматизированные системы на пути модель – алгоритм – программа. В данной работе применяется метод граничных элементов. Известны преимущества и недостатки обоих методов, не вступая в полемику, скажем, что метод доказательства существования решения задачи Гольдштика, представленный в данной работе, опирается на методы теории потенциала и теории граничных интегральных уравнений. Именно поэтому, мы отдаем в дальнейшем предпочтение МГЭ, поскольку не возникает проблем с точностью численного решения краевой задачи (проблемы подобного рода в работе обусловлены только величиной шага разбиения и возможностями вычислительной техники). Все требования точности вычислений выполнены и проверяются на тестовых примерах.

Последнее является необходимым условием функционирования программы расчета и применения ее для проектирования и повышения эксплуатационных характеристик действующих ГТС.

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

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

1. исследование математической модели Лаврентьева-Шабата об отрывных течениях идеальной жидкости в галереях с препятствиями и поворотами;

2. приложение разработанной модели к построению кинематической картины движения жидкости в обводных галереях шлюзовых камер;

3. доказательство существования решения и доказательство сходимости итерационного процесса к точному решению;

4. получение достоверных численных результатов о распределении скоростей в обходных галереях шлюзовых камер, что является теоретическим обоснованием для проектирования и конструирования обходных галерей шлюзовых камер;

5. исследование распределения гидродинамических нагрузок в областях галерей, непосредственно примыкающих к гидрозатворам при различных режимах работы последних;

6. расчет динамических характеристик НДС плоского гидрозатвора стоечного и ригельного типов и сравнение их с характеристиками в статичном режиме;

7. разработка и построение алгоритма и программной реализации численных результатов;

8. оценка точности разработанной математической модели и программы;

сравнение разработанного алгоритма на основе метода граничных интерполяции;

9. представление полученных результатов в форме пользовательского графического интерфейса.

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

Основой теоретических методов являются положения теории граничных интегральных уравнений, гидродинамической модели Лаврентьева-Шабата идеальной жидкости, методов Фурье и Навье решения уравнений в частных производных и метода итераций. Расчеты производились с использованием оригинальных программ, созданных в символьной среде Maple.

математической модели Лаврентьева-Шабата к движению идеальной жидкости в обходных галереях шлюзов. Доказано существование решения обобщенной краевой задачи Гольдштика о сопряжении вихревых течений идеальной жидкости. Метод доказательства дает алгоритм для численной реализации решения краевой задачи на основе метода итераций и метода граничных элементов. Проведен сравнительный анализ воздействия динамических и статических нагрузок на характеристики НДС затворов затопленных обходных галерей.

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

Конструкции такого рода эксплуатируются на шлюзах Нижнесвирского гидроузла. На основе полученного распределения скоростей вычислены динамические нагрузки на затворы в галереях. Выполнен расчет динамических характеристик НДС гидрозатворов при различных режимах работы последних (ранее подобные исследования проводились исключительно в статике).

Полученные расчеты применялись для моделирования явления вибрации, возникающего в процессе подъема затвора. Разработан алгоритм и программная реализация алгоритма; полученные результаты представлены в форме пользовательского графического интерфейса.

Основные новые научные результаты, полученные в работе и выносимые на защиту.

В рамках математической модели Лаврентьева-Шабата об отрывных течениях идеальной жидкости построена кинематическая картина распределения скоростей в каналах произвольной формы.

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

Предложен алгоритм численной реализации модели на основе метода граничных элементов. Создана и зарегистрирована в ФСИС программа на ЭВМ для расчета отрывных областей методом граничных элементов в Для модельных областей дана оценка точности разработанного метода.

динамические нагрузки на гидрозатворы затопленных обходных галерей.

рассчитываемыми по известным методикам и СНиПам.

Построена математическая модель напряженно-деформированного состояния затворов обходных галерей под действием динамических В рамках рассматриваемой модели течения дано объяснение явления кавитации. Получены конкретные результаты и даны методические рекомендации по устранению явления кавитации.

Представление полученных результатов в форме пользовательского графического интерфейса.

Таким образом, диссертация удовлетворяет следующим пунктам из паспорта специальности:

1. Разработка новых математических методов моделирования объектов и явлений.

3. Разработка, обоснование и тестирование эффективных вычислительных методов с применением современных компьютерных технологий.

4. Реализация эффективных численных методов и алгоритмов в виде комплексов проблемно-ориентированных программ для проведения вычислительного эксперимента.

5. Комплексные исследования научных и технических проблем с применением современной технологии математического моделирования и вычислительного эксперимента.

7. Разработка новых математических методов и алгоритмов интерпретации натурного эксперимента на основе его математической модели.

Достоверность и апробация результатов.

Внедрение. Результаты диссертационной работы нашли практическое применение в проектных работах по инновационным направлениям на ЗАО «Канонерский судоремонтный завод», в СПКТБ «Ленгидросталь», в ГУМРФ они используются при подготовке специалистов по направлению «Прикладная математика и информатика», разработана и зарегистрирована в Реестре программ для ЭВМ в Федеральной службе по интеллектуальной собственности программа.

Публикации. Результаты работы опубликованы в научно-технических изданиях. Всего опубликовано 24 статьи (из них 22 статьи в журналах, рекомендованных ВАК), одна монография и 9 тезисов докладов. Все результаты, входящие в текст диссертации выполнены автором лично, в совместных работах личный вклад автора составляет 70 %.

Апробация. Основные положения диссертационной работы представлялись на научных семинарах кафедры математики и кафедры прикладной математики ГУМРФ и докладывались на Всероссийской научнометодической конференции СПб, СПГУВК, 1994; на III международной конференции «Новое в охране труда и окружающей среды», Алмааты, 1998; на международной конференции, посвященной юбилею курорта «Сестрорецк», СПб, 1999; на XI международной научной конференции, посвященной 70летию профессора В.П.Дьяконова, Смоленск, 2010; на Юбилейной XIII СанктПетербургской международной конференции «Региональная информатика (РИСанкт-Петербург, 24-26 октября 2012 г; на Российской летней школе «Математическое моделирование фундаментальных объектов и явлений в системах компьютерной математики» (ММ СКМ2), 6-10 сентября 2010 г., Казань-Яльчик, на Герценовских чтениях – 2013 г. «Некоторые актуальные проблемы современной математики и математического образования» 15- апреля 2013г, СПб; на Международной конференции «Нелинейные уравнения и комплексный анализ», 18-22 марта 2013 г., Уфа; на IV межвузовской научнопрактической конференции "Современные тенденции и перспективы развития водного транспорта России» 15-16 мая 2013 г. в ГУМРФ имени адмирала С.О.

Макарова.

Структура и объем работы. Диссертация представлена в форме рукописи, состоящей из введения, пяти глав, заключения и 3 приложений.

Общий объем рукописи 319 стр., в том числе 105 рисунков, 17 таблиц и список использованных источников из 174 наименований.

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

.

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

Вторая глава диссертации посвящена теоретическим аспектам модели Лаврентьева-Шабата [2], которая описывает плоское течение жидкости при обтекании препятствия. Данная модель наиболее всего отвечает течению жидкости в водопропускных галереях шлюзовых камер. В данной модели идеальной жидкости движение делится на две зоны: зону основного потока и зону отрывного течения. Течение в зоне основного потока предполагается потенциальным, а области отрывных вихревых течений присутствуют в малых частях за некоторым естественным препятствием, будь-то выступ, или полуприподнятый затвор обходной галереи, или каверна на дне камеры.

Условие постоянства завихренности является естественным упрощением (принцип Прандтля-Батчелора [3]) для предельного течения вязкой жидкости, когда вязкость стремится к нулю. То, что линия раздела отрывного течения, которая в данном случае является линией тока, не известна, представляет основную трудность в задаче Лаврентьева-Шабата.

Во втором параграфе доказывается, что задача Лаврентьева-Шабата сводится к некоторой краевой задаче Римана для аналитических функций, поэтому возможно применение комплексного метода граничных элементов [7,8].

В третьем параграфе исследуется вариант задачи Лаврентьева-Шабата нелинейного (квазилинейного) уравнения эллиптического типа Математические модели движения жидкости в прямолинейных каналах с выступами основаны на решении задачи Гольдштика. Для галерей с поворотами модель Лаврентьева-Шабата недостаточно точно отражает опытные измерения [5], а именно такие галереи сконструированы на действующих судоходных шлюзах. В четвертом параграфе объясняется, почему требуется усовершенствование модели Лаврентьева-Шабата. Данное уточнение заключается в сопряжении двух вихревых областей. В отличие от исходной задачи область основного потока носит ярко выраженный вихревой характер. В области отрывного течения возникает дополнительная вихревая область с собственной завихренностью, не совпадающей с завихренностью основного течения. В пятом параграфе доказывается существование решения обобщенной задачи Гольдштика.

Теорема 2.1. Пусть заданы ограниченная область D на плоскости, положительная константа 1 и неотрицательная функция на границе области D. Тогда для достаточно больших величин 0 существует функция, удовлетворяющая решению обобщенной задачи Гольдштика с непустой областью D0 D, где 1.

В доказательстве используется метод [3,11] непрерывной аппроксимации задачи с разрывной нелинейностью линейными разрывными задачами. Данный метод доказательства существования решения одновременно дает метод расчета и построения отрывной области. Важное замечание представляет утверждение о непустоте области D0. Известно, что в задаче Гольдштика нет единственности решения (величина завихренности служит спектральным параметром), и имеются тривиальные решения, когда одна из областей вырождается. Данные положения подробно обсуждаются в разделе про вихри Рэнкина.

В третьей главе диссертации разрабатывается алгоритм построения оператора Грина задачи Гольдштика и обобщенной задачи Гольдштика в итерационному процессу из линейных разрывных задач Пуассона для уравнения Лапласа, каждая из которых решается методом граничных элементов. В первом параграфе разрабатывается алгоритм решения задачи Дирихле для индикатора подобласти G G0. По принципу суперпозиции для линейных уравнений частное решение неоднородной задачи Дирихле выбирается в виде потенциала площади – свертки с фундаментальным решением для оператора Лапласа:

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

где - выбранное направление нормали к поверхности и v – неизвестная плотность на. Плотность потенциала находим из граничного сингулярного интегрального уравнения Последняя формула не имеет места в угловых точках области, требуются корректировки данного уравнения в зависимости от величины угла. На самом деле можно при численной реализации оставить формулу для гладких граничных точек, но аккуратно следить за тем, чтобы узловые точки при интерполяции не совпадали с угловыми точками границы области.

Во втором параграфе выполняется дискретизация численного решения сингулярного уравнения с разбиением границы на граничные элементы j граничными точками которые называются узлами. Интегрируемые функции интерполируются ступенчатыми функциями. Поскольку такого рода интерполяция применяется под знаком интеграла, то незачем использовать более сложную и громоздкую интерполяцию кусочно-линейными функциями [7, 8], на точность квадратурных формул это не влияет. Записывая уравнение в каждом из узлов zi, получим систему линейных уравнений относительно использованием математического пакета Maple. После вычисления искомой плотности значение функции u в точке z x, y, z, достаточно далекой от границы, может быть вычислено квадратурами для потенциала двойного слоя.

Для точки z x, y, близкой к границе, сказывается сингулярность интеграла потенциала двойного слоя, поэтому можно воспользоваться геометрическим смыслом потенциала двойного слоя от индикатора дуги.

В третьем параграфе дается решение смешанной задачи Пуассона для уравнения Лапласа в произвольной области. Специфика задач гидродинамики в том, что на входе и выходе канала, как правило, значения функции тока не известны. При этом условие Неймана, заданное на входе-выходе является однородным и не улучшает число обусловленности задачи. Проследим основные отличия при решении данной краевой задачи от рассмотренной в пунктах 3.1-3.2. Формула Грина для гармонических функций дает При предельном переходе при z (с теми же замечаниями об угловых точках) и дискретизации получим систему линейных уравнений относительно неизвестных значений гармонической функции и ее нормальной производной в точках z j :

Разделяя известные и неизвестные значения гармонической функции и ее нормальной производной, и учитывая, что в каждом узле – середине граничного элемента нам известно только значение функции или ее использовании пакета LinearAlgebra.

В четвертом параграфе конструируется алгоритм решения задачи Гольдштика (1) в произвольной области. Зададим начальную область D0, в которой функция 1 удовлетворяет уравнению Пуассона, а в области D1 гармонична, кроме того функция 1 непрерывно дифференцируема, удовлетворяет граничным условиям, но на линии раздела областей D0 и D может не обращаться в нуль. Найдем линию нулевого уровня функции 1, которая ограничивает некоторую область D1. Далее итеративно определим последовательность функций где G(, z ) – функция Грина для уравнения Пуассона с заданными краевыми условиями для полной области D. Последовательность функции n сходится к непрерывно дифференцируемой функции, а последовательность областей Dn сходится к области D, для которых выполняются условия задачи Гольдштика. На линии раздела вихревого и потенциального течений функция должна обращаться в нуль, поэтому на этом шаге будем корректировать и границу, и величину завихренности. На рис. 1 показана схема исправления линии раздела.

В пятом параграфе выполнена оценка точности разработанного метода граничных элементов для решения задачи Гольдштика. Известных примеров отрывной области вихревых течений, которые можно построить аналитически, очень мало. Проблема – в сильном нелинейном переплетении в уравнении функций тока и завихренности. Для нашей задачи лучше всего подходит тест на основе вихря Рэнкина [1]. Данное течение строится аналитически и описывается как При фиксированной завихренности 4Сe / R2 может быть два различных решения, в случае равенства – одно. На рис. 2 указана зависимость величины завихренности от радиуса окружности, разделяющей области вихревого и потенциального течений в круге при C 1, R 1. Предлагаемая итерационная схема должна давать при определенных условиях оба решения для вихрей Рэнкина. Для проверки точности метода в качестве начальной области вихревого течения возьмем область в форме круга определенного радиуса 0 r 1. Вне круга течение предполагается потенциальным. Процесс итераций с поиском линии уровня задается следующей итерационной формулой типа где Dr – круг радиуса r, а G(w, z ) – функция Грина для единичного круга.

Функции r являются решением разрывного уравнения Пуассона с граничным условием r ( w) 1 при | w | 1. Но решением задачи Гольдштика эти функции будут только при r a1 или r a2. Если же радиус области интегрирования не совпадает ни с одним из указанных, то формула (4) задает итерационную схему для нахождения радиуса круга, представляющего новую область интегрирования.

На рис.3 построен график зависимости радиуса нулевой линии уровня от радиуса круга интегрирования (при заданной завихренности 20 ). Точки 2 и 3 как раз представляют решения для двух вихрей Рэнкина. Соответствующие им значения радиуса области, равные a1 и a2, являются неподвижными точками данного отображения. Точки 1 и 4 соответствуют двум тривиальным решениям, когда в круге наблюдается потенциальное течение с пустой областью завихренности (1) и вихревое течение во всем круге (4). Вид данной кривой объясняет, почему итерационная схема дает решение только для правой ветки вихрей Рэнкина (из рис.2). Дело в том, что здесь работает стандартная теорема о неподвижной точке, когда производная отображения не превосходит по абсолютной величине 1. Для того чтобы получить сходимость для отрывных течений малого радиуса (для левой ветви рис.2) требуется несколько иная схема метода граничных интегральных уравнений с корректировкой величины завихренности на каждом шаге итеративного процесса. В противном случае, если начальный радиус взять из интервала 0 r a1, то вихря Рэнкина не получим, а решение сойдется к тривиальному потенциальному с пустой вихревой областью. Такой эффект замечен Гольдштиком, где для поиска второй неподвижной точки предполагается исследовать определенный функционал энергии. Им замечено, что к результату приводят локальные экстремумы функционала, а абсолютные максимум и минимум дают лишь тривиальные течения. Все это усложняет поиск второго решения. На самом деле, хороший результат дает применение метода простой итерации для обратной функции.

Рис. 2. Зависимость завихренности от Рис. 3. Зависимость радиуса нулевого Интерполяционные задачи на основе формулы Келдыша-Седова [2] исследуются в шестом параграфе. В седьмом параграфе сравнивается точность двух методов граничных элементов: использующего теорию граничных интегральных уравнений и использующего комплексную интерполяцию.

В четвертой главе построены основные отрывные течения в камерах различной формы. В первом параграфе примером описанной техники является модель плоского течения идеальной жидкости в прямоугольной камере с выступом, которая послужит эталоном для моделирования течения жидкости в реальных галереях шлюзовых камер. На рис. 4 представлены области отрывных течений для различных профилей скорости на входе и выходе. В левой части для функции тока имеем y / 7 и в правой части y 2 / 49 на входе и выходе, соответствующие большим скоростям на верхней или на нижней стенке на рис. 4.

Во втором параграфе исследуется течение идеальной жидкости в камере с прямоугольной траншеей с плоским дном. Именно для камер подобной формы уделено большое внимание в работах [2, 3, 9]. Главный вычислительный эксперимент заключается в моделировании кинематической картины в траншее прямоугольного сечения 1:5 с закреплением концов кривой, представляющей линию, отделяющую отрывное вихревое течение. Метод простой итерации дает решение с рис. 5. Для получения решения на рис. 6 требуется итерационная схема с корректировкой завихренности на каждом шаге. В данном случае получается несвязная вихревая зона. Из двух решений с большей вихревой зоной (рис. 5) и с меньшей – несвязной (рис. 6) физический смысл имеет именно второе. Метод обратной итерации с корректировкой завихренности как раз позволяет отсечь первое постороннее решение.

В третьем параграфе исследуется обобщенная задача Гольдштика.

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

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

завихренностями. Наблюдается два вихревых течения (первое с завихренностью 1 1.1 в основной части камеры, второе с завихренностью 11 11.596 представляет отрывное течение). Процесс итераций осуществлялся с корректировкой завихренности, поэтому вихревая область корректно определилась. Для сравнения на рис. 8 изображена кинематическая картина течения в камере без выступа. Здесь отсутствует отрывная область и картина течения иллюстрирует решение обычного уравнения Пуассона с заданными граничными условиями.

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

Рис. 9. Посторонние решения задачи Гольдштика.

В шестом параграфе дано первое приложение для проектных работ, а именно построена кинематическая картина течения на выходных каналах конструкции типа Нижнесвирского гидроузла, имеющих гасители скоростей потока. При изучении предложенных гасителей были произведены лабораторные исследования [5] скоростей потока в выходных отверстиях галерей как без применения гасителей, так и с ними. Но математического обоснования опыта Городенского не было предложено. На рис. 10 изображены построенные области отрывных течений в камере с выступом и без выступа.

Форма областей максимально приближена к формам каналов из обходных галерей шлюзовых камер. Стенки каналов представляют собой дуги окружностей, но в отличие от модельного случая, когда канал имеет кольцевую форму, теперь дуги окружностей не являются концентрическими. В нашем числовом опыте рассмотрены два вида выступа: с вогнутой стороной со стороны течения и с противоположной стороны. Данные различия не принципиальны в достижимом эффекте. Оба вида выступов достаточно хорошо выполняют функцию по гашению скоростей. Кроме того оба данных выступа лучше чем выступ округлой формы.

Рис. 10. Кинематика в галереях с двумя видами выступа.

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

Рис. 11. Эпюры скоростей. Пунктиром и сплошной линией обозначены эпюры скоростей в галереях с двумя видами выступа, а штрихпунктирной линией – В пятой главе построенные кинематические картины распределения скоростей находят приложение к расчету гидродинамических нагрузок на водопропускные ворота в обводных галереях шлюзовых камер при процессах наполнения и опорожнения шлюзовых камер через затопленные обходные галереи (рис. 12) Рис. 12. Схема шлюза. 1 – Подводящий канал; 2 – Камера шлюза; 3 – Верхние ворота; 4 – Водопроводная галерея; 5 – Плоский затвор.

Вихревая область, формирующаяся за затвором (выделено окружностью на рис.12), влечет целесообразным использование модели Лаврентьева-Шабата.

В численном эксперименте для модели водопроводной галереи построена кинематическая картина движения потока (рис. 13). На нижней стенке функция тока н 0, на верхней стенке функция тока в 20, на входе/выходе функция тока меняется линейно 7.4y, средняя скорость на входе и на выходе равна 7.4 м/с.

Рис. 13. Построенная при моделировании картина течения.

Полученная картина течения и, как следствие, распределение скоростей позволяет определить величину гидродинамического давления на затвор.

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

Кроме того, схема потенциального течения также дает, ввиду симметрии, одинаковое распределение скоростей на разных сторонах затвора, а это влечет одинаковые давления на разных сторонах затвора на одном уровне. С другой стороны, опытные измерения показывают, что в процессе подъема затвор испытывает вибрацию. Ни гидростатические, ни потенциальные модели движения жидкости не дают правильного объяснения этого явления.

Применение модели Лаврентьева-Шабата дает не симметричное распределение скоростей, что ведет к не симметричным гидродинамическим нагрузкам на затвор. На рис. 14 изображен график результирующей силы давления на затвор для полученного течения. Сразу заметим, что данный график давлений существенно отличается от гидростатических нагрузок, которые в основном используются при расчетах напряженно-деформированного состояния затворов. Величины результирующих давлений уменьшаются по мере открытия затвора, но главным является наличие постоянного действующего момента сил давлений, испытываемого затвором.

Рис. 14. Результирующее давление на затвор при различных режимах высоты В работе рассчитано давление на затвор для различных значений расхода, (рис. 15). Увеличение расхода влечет увеличение входных скоростей, что, в свою очередь, влечет увеличение абсолютных скоростей в области затвора.

Рис. 15. Зависимость давлений от входных скоростей. 1 – график для 16 ;

2 – график для 18 ; 3 – график для 20. Давление отмечено в кПа.

Отметим, что найденные величины давлений достаточно велики для учета гидростатическими нагрузками на гидрозатвор.

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

Дифференциальное уравнение «технической» теории упругости для прогиба имеет вид где 1i, i 1,..., K1; 2 j, j 1,..., K2. – коэффициенты жесткости. Дельтафункции ( y yi ) и ( x x j ) нужны для склейки решений по ребрам на отдельных участках. В данной главе мы отступаем от основного метода граничных элементов и, ввиду простой формы области, решение ищем в виде однократных рядов Фурье по обеим переменным:

По принципу суперпозиции решение линейного уравнения представим в виде трёх слагаемых — решения w 0 ( x, y) уравнения гладкой, неподкрепленной пластины решения w1 ( x, y) уравнения с правой частью – линейной комбинацией функций по y решения w 2 ( x, y), подобного предыдущему уравнению с заменой y на x :

дифференциальное уравнение коэффициентами.

Снова по принципу суперпозиции, но для обыкновенного линейного уравнения, общее решение уравнения представим в виде суммы трёх слагаемых, содержащих линейные комбинации функций Хевисайда где Z mk ( x) — фундаментальная система решений однородного уравнения, Ck — произвольные постоянные интегрирования; и функции Теперь для получения общего решения необходимо выразить числовые коэффициенты w m ( x j ), что осуществляется подстановкой значений x x j.

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

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

Представленные результаты графиков прогибов, существенно различные для различных режимов работы затвора, как на рис. 17 (слева прогибы при закрытом затворе, справа – при открытом на 30%.). Величины разносторонних прогибов соизмеримы с прогибами, испытываемыми при гидростатической нагрузке, в результате чего, по мере открытия затвора приводят к вибрации.

Рис. 17. Прогибы затвора в статике и динамике.

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

В шестой главе диссертации представлены результаты о создании пользовательского графического интерфейса в математических пакетах Maple (первый параграф), MatLab (второй параграф). Интерфейс предназначен для удобного ввода исходных данных и получения в результате вычислений картины течения и предполагается использовать в учебных целях на факультете информационных технологий и гидротехническом факультете ГУМРФ при моделировании течений идеальной жидкости. На рис.19 представлен пример интерфейса в Maple.

В Приложении 1, 2, 3 представлены листинги программ для интерфейсов в средах Maple и MatLab и программы «Метод интегральных уравнений для расчета вихревых течений в водопропускных галереях шлюза». В программе реализованы: ввод данных, их обработка и анализ, сохранение результата.

Выходные формы соответствуют требованиям федеральных органов. Область программирования: символьная среда Операционная система:

Windows98/NT/2000/2003/XP/Vista/7.

Заключение. В данной работе метод граничных элементов применяется для численного моделирования отрывных течений в задачах гидравлики. В качестве основной модели подробно исследуется модель Лаврентьева-Шабата, которая наиболее удобна для исследования течений идеальной жидкости в камерах с выступами, с кавернами, в поворотных галереях гидротехнических сооружений. Даны качественные оценки и получены количественные результаты распределения скоростей в затопленных водопропускных галереях судоходных шлюзов. Последнее дало возможность исследовать динамические характеристики, испытываемые затворами гидротехнических сооружений.

Получены следующие результаты Исследована математическая модель Лаврентьева-Шабата в применении водопропускные галереи.

Разработана математическая модель и построена кинематическая картина отрывных течений в каналах произвольной формы: в кавернах, в областях за затворами, в обходных поворотных галереях.

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

Разработан алгоритм численной реализации модели на основе метода граничных элементов. Создана и зарегистрирована в ФСИС программа на ЭВМ для расчета отрывных областей методом граничных элементов.

Для модельных областей на основе вихрей Рэнкина дана оценка точности разработанного метода. Дано сравнение двух методов граничных элементов: метода на основе граничных интегральных уравнений и метода комплексной интерполяции.

гидродинамические нагрузки на затворы затопленных обходных галерей.

Построена математическая модель напряженно-деформированного состояния затворов обходных галерей под действием динамических нагрузок при различных режимах работы. Дано сравнение полученных результатов с традиционных исследований в статичном режиме.

В рамках рассматриваемой модели течения дано объяснение явления вибрации затвора и кавитации за затвором.

Создан пользовательский интерфейс для ввода исходных данных и получения в результате вычислений картины отрывных течений для информационных технологий и гидротехнического в среде Maple и Список использованных источников Мелешко В.В. Динамика вихревых структур./ В.В.Мелешко., М. Ю.

Константинов. – Киев: Наукова Думка. - 1993. -290с.

Лаврентьев М.А. Проблемы гидродинамики и их математические модели./ Лаврентьев М.А., Шабат Б.В. - М.: Наука. Гл. ред. физ.-мат. лит. – 1977. с.

Гольдштик М.А. Вихревые потоки. / Гольдштик М.А. - Новосибирск:

Наука, Сиб. Отделение. – 1981. - 368 с.

Голоскоков Д.П. Численно-аналитические методы расчета упругих тонкостенных конструкций нерегулярной структуры./ Голоскоков Д.П. – СПб.: Изд-во Кардакова. - 2006. -218с.

Городенский Н.Б. Гаситель вредных скоростей воды в выходных галереях судоходных шлюзов./ Городенский Н.Б., Эндер Г.В. - Авторское свидетельство на изобретение №36916. - 31.05.1954.

Михайлов Б.К. Деформированность и устойчивость пространственных пластинчатых систем с разрывными параметрами./ Михайлов Б.К., Кипиани Г.О. - С-Петербург: Стройиздат. – 1996. - 443 с.

Громадка Т. Комплексный метод граничных элементов в инженерных задачах./ Громадка Т., Лей Ч. - М.: Мир. - 1990. – 303 с.

Афанасьев, К. Е. КМГЭ для решения плоских задач гидродинамики и его реализация на параллельных компьютерах. / К. Е. Афанасьев, С. В.

Стуколов. - Кемерово: Изд-во КемГУ, 2001. - 206 с.

Потапов Д.К. О существовании луча собственных значений для уравнений эллиптического типа с разрывными нелинейностями в критическом случае / Потапов Д.К.// Вестн. СПбГУ. Сер. 10 (Прикладная математика.

Информатика. Процессы управления). 2004. Вып. 4. С. 125–132.

10. Потапов Д.К. О решениях задачи Гольдштика./ Потапов Д.К.//Сибирский журнал вычислительной математики. – 2012. – т.15. – вып.4.

Вайнштейн И.И. Решение двух дуальных задач о склейке вихревых и 11.

потенциальных течений вариационным методом М.А. Гольдштика // Журнал Сибирского федерального университета. Математика и физика, 2011, N4(3), 320–331.

Основные результаты Васина А.В. в научных журналах, входящих в перечень ВАК:

Васин А.В. Об одном необходимом условии для интерполяции функциями класса Липшица./Васин А.В. //Записки научных семинаров ПОМИ им. В.А.

Стеклова. – 2004. – т. 315.

2. Vasin A.V. A necessary condition of interpolation by functions of the Lipschitz classes. /Vasin A.V. //Journal of Mathematical Sciences. – 2006. – v.134.

Васин А.В. Граничная интерполяция в слабых классах Липшица. / Васин А.В. //Исследования по линейным операторам и теории функций.36:

Записки научных семинаров ПОМИ им. В.А. Стеклова.– 2008. – т.355.

4. Vasin A.V. Boundary interpolation of weak Lipschitz classes./Vasin A.V.

//Journal of Mathematical Sciences. –2009. –v.156.

Васин А.В. Расчет методом граничных элементов потоков жидкости в модельных областях./Васин А.В., Болбин С.Н. // Журнал университета водных коммуникаций. – 2011. – Вып.2(10).

Васин А.В. К вопросу о течении жидкости в водопроводных галереях коммуникаций. – 2011. – Вып.3(11).

Васин А.В. Нахождение линии раздела областей с потенциальным и вихревым течением./Васин А.В., Тимофеева О.А. //Журнал университета водных коммуникаций. – 2012. – Вып. 2(14).

Васин А.В. Математическое моделирование опыта Н.Б.Городенского о университета водных коммуникаций.–2012.–Вып.3(15).

Васин А.В. Моделирование зон вихревых течений в обводных галереях шлюзовых камер./Васин А.В., Тимофеева О.А. //Гидротехническое строительство. –2013, № 2.

10. Васин А.В. Определение линии раздела областей вихревых течений. /Васин А.В.//Вестник С-Петербургского Университета. Серия 10. (Прикладная математика. Информатика. Процессы управления). –2013. –№1.

11. Васин А.В. Приложение задачи Гольдштика к расчету гидродинамических нагрузок на гидрозатворы обходных галерей. /Васин А.В.//Речной транспорт (XXI век). –2013. –№2.

12. Васин А.В. Оценка точности метода граничных элементов в задачах математического моделирования движения потока по обводным галереям шлюзов. / Васин А.В. //Морские интеллектуальные технологии. – 2013. – 13. Васин А.В. Построение математической модели динамических нагрузок на затворы обходных галерей. / Васин А.В., Голоскоков Д.П., Тимофеева О.А.

// Журнал университета водных коммуникаций. – 2013. – Вып.3(19).

14. Васин А.В. Сравнение влияний динамических и гидростатических нагрузок на гидрозатворы обводных галерей. /Васин А.В., Голоскоков Д.П., Тимофеева О.А.// Журнал университета водных коммуникаций,. – 2013. – Вып.4(20).

Васин А.В. Проекторы на Lp-пространства полианалитических функций.

15.

/Васин А.В. //Записки научных семинаров ПОМИ им. В.А. Стеклова.–1991.

A.V.Vasin. Projections on Lp–Spaces of Polyanalytic Functions. /Vasin 16.

А.V.//Journal of Mathematical Sciences.–1994. –v.71.

17. Vasin A.V. Bounded projections onto spaces of area integrable polyanalytical functions. /Vasin A.V.//Indiana University Mathematical Journal.–1994. –v.43.

Васин А.В. Полианалитические формы на компактных римановых 18.

поверхностях. /Васин А.В. //Записки научных семинаров ПОМИ им.В.А.Стеклова.–1997. –т.247.

19. A.V.Vasin. Polyanalytic forms on compact Riemann surfaces./VasinA.V.// Journal of Mathematical Sciences. –2000. –v.101.

Васин А.В. Предельное множество фуксовой группы и лемма Дынькина.

20.

/Васин А.В. //Записки научных семинаров ПОМИ им.В.А.Стеклова. –2003.

A.V.Vasin. The Limit Set of a Fuchsian Group and Dyn’kin’s Lemma.

21.

/A.V.Vasin.//Journal Mathematical Sciences. –2005.–v.129.

Васин А.В. Локальная гладкость аналитической функции в сравнении с 22.

А.Н.//Алгебра и анализ. –2013. – т.25.– вып.3.

Монографии:

течений./Васин А.В. – СПб.:Изд-во ГУМРФ имени С.О. Макарова. – 2013.

Публикации в других изданиях:

Васин А.В. О симметричности ортогональности в нормированных 24.

инженерных расчетах. – СПб: Изд-во СПГУВК. – 1999.

Васин А.В. Автоматизация моделирования отрывных течений в обводных 25.

галереях шлюзовых камер./Васин А.В. //Вестник ИТАРК.–2012, –№1(3).

Публикации в материалах конференций:

26.

поверхностях.//Материалы Всерос.научн.метод.конф, СПб,СПГУВК,1994.

Богатушин В.Я. Моделирование сложных динамических объектов в виде 27.

систем линейных дифференциальных уравнений./ Богатушин В.Я., Васин А.В. //Материалы 3 межд. конф «Новое в охране труда и окружающей среды», Алмааты, 1998.

Богатушин В.Я. Моделирование сложных динамических объектов в виде 28.

систем линейных дифференциальных уравнений./ Богатушин В.Я., Васин А.В. // Материалы межд. конф. посвящ. юбилею курорта «Сестрорецк», Болбин С.Н. Задача течения жидкости по обводному каналу шлюза./ Васин 29.

А.В., Болбин С.Н. //Системы компьютерной математики и их приложения:

материалы XI межд. научн. конференции, посвященной 70-летию профессора В.П.Дьяконова.-Смоленск: Изд-во СмолГУ, 2010. Вып.11-342с.

Васин А.В. О полианалитических функциях на компактных римановых 30.

поверхностях. /Васин А.В. //Системы компьютерной математики и их приложения: материалы XI межд. научн. конференции, посвященной 70летию профессора В.П.Дьяконова.-Смоленск: Изд-во СмолГУ, 2010.

Вып.11.-342с.

Васин А.В. Обобщение задачи Гольдштика о вихревых течениях./ Васин 31.

А.В.//Региональная информатика (РИ-2012). Юбилейная XIII СанктПетербургская международная конференция. Санкт-Петербург, 24- октября 2012г.: Материалы конференции. СПОИСУ. – СПб, 2012.-420с.

Васин А. Математические модели прикладных NP-полных задач./Васин 32.

А.В, Нырков А.П., Соколов С.С.//Труды Российской летней школы «Математическое моделирование фундаментальных объектов и явлений в системах компьютерной математики» (ММ СКМ2), 6-10 сентября 2010 г., Казань-Яльчик. – Казань: Изд-во «Фолиантъ», 2010.-216с.

Васин А. В. Математическое моделирование динамических нагрузок на 33.

гидрозатворы. /Васин А.В.//Некоторые актуальные проблемы современной математики и математического образования. Герценовские чтения – 2013.

Материалы научной конференции, 15-20 апреля2013г. – СПб.: Изд. РГПУ им. А.И. Герцена, 2013. – 292с.

34. Kislyakov S.V. The loss of smoothness twice for outer function compared to its module./Kislyakov S.V., Medvedev A., Vasin A.V. //Международная конференция «Нелинейные уравнения и комплексный анализ», 18-22 марта 2013г., Уфа. – Уфа, 2013.