автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Комплексная система моделирования гидрофизических характеристик замкнутых соленых стратифицированных озер

кандидата физико-математических наук
Якубайлик, Татьяна Валерьевна
город
Красноярск
год
2014
специальность ВАК РФ
05.13.18
цена
450 рублей
Диссертация по информатике, вычислительной технике и управлению на тему «Комплексная система моделирования гидрофизических характеристик замкнутых соленых стратифицированных озер»

Автореферат диссертации по теме "Комплексная система моделирования гидрофизических характеристик замкнутых соленых стратифицированных озер"

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

Якубайлик Татьяна Валерьевна

КОМПЛЕКСНАЯ СИСТЕМА МОДЕЛИРОВАНИЯ ГИДРОФИЗИЧЕСКИХ ХАРАКТЕРИСТИК ЗАМКНУТЫХ СОЛЕНЫХ СТРАТИФИЦИРОВАННЫХ ОЗЕР (НА ПРИМЕРЕ ОЗ. ШИРА)

05.13.18 - Математическое моделирование, численные методы и комплексы программ

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

- 3 ИЮЛ 2014

Новосибирск — 2014 005550297

005550297

Работа выполнена в Федеральном государственном бюджетном учреждении науки Институте вычислительного моделирования Сибирского отделения Российской академии наук, г. Красноярск

Научный руководитель: кандидат физико-математических наук, доцент

Компаниец Лидия Алексеевна

Официальные оппоненты: Воеводин Анатолий Федорович,

доктор физико-математических наук, профессор, ИГиЛ СО РАН, г. Новосибирск, главный научный сотрудник

Платов Геннадий Алексеевич, доктор физико-математических наук, ИВМиМГ СО РАН, г. Новосибирск, ведущий научный сотрудник

Ведущая организация: Федеральное государственное автономное образовательное учреждение высшего профессионального образования «Южный федеральный университет», г. Ростов-на-Дону

Защита состоится 16 сентября 2014 года в 11:30 на заседании диссертационного совета ДМ 003.046.01 на базе Федерального государственного бюджетного учреждения науки Института вычислительных технологий Сибирского отделения Российской академии наук по адресу: 630090, г. Новосибирск, пр. Академика Лаврентьева, 6, конференц-зал ИВТ СО РАН.

Материалы по защите диссертации размещены на официальном сайте ИВТ СО РАН: Ь11р:/Ау\у\у.{с1.п5с.ги/я1Чег)аее.пЬг)?Раее1Р=17

С диссертацией можно ознакомиться в библиотеке ИВТ СО РАН.

Автореферат разослан «11» июня 2014 г.

Учёный секретарь диссертационного совета, к. ф.-м. н.

Лебедев А. С.

Общая характеристика работы

Актуальность. Соленые озера являются важными природными объектами, в которых заключены разнообразные запасы минеральных и биологических ресурсов. Они представляют собой большую экономическую ценность, так как являются одним из основных мест рекреационного туризма, медицинским ресурсом по запасам лечебной грязи и соли. Мировым сообществом признана необходимость «изучения, рационального использования и сохранения природы соленых озер» (Манифест, принятый участниками V Международного симпозиума по изучению соленых озер, Боливия, 1991 г.). Современные соленые озера имеют научное значение как природные седи-ментологические лаборатории, которые не только хранят записи обо всех па-леоклиматических изменениях, но и являются природной площадкой для изучения солеобразующих процессов и изменений окружающей среды.

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

Логическим продолжением может служить объединение нескольких направлений для исследований конкретного водоема, что и обосновывает актуальность данной работы. В диссертационной работе предлагается комплексный подход к исследованию неглубоких замкнутых соленых стратифицированных озер, разработанный и опробованный на изучении гидрофизических характеристик озера Шира. Идея исследования заключается в применении технологии вычислительного и натурного экспериментов для построения моделей, описывающих гидрофизические характеристики изучаемого уникального водного объекта. Методология вычислительного эксперимента разработана в трудах академиков А. А. Самарского, Н. Н. Яненко, Ю. И. Шокина.

Цель работы состоит в создании на основе технологии вычислительного эксперимента комплексной системы моделирования соленых стратифицированных озер, ориентированной на изучение их гидрофизических характеристик и применении ее для изучения ветровых течений в реальном водоеме (озере Шира, расположенном в республике Хакасия).

В соответствии с этой целью в работе поставлены и решаются следующие задачи:

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

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

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

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

5. Дать оценку гидрофизического режима озера Шира в летний период, используя выявленные закономерности.

Объектом исследований являются физические процессы, специфические для замкнутых стратифицированных соленых водоемов.

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

На защиту выносятся следующие результаты, соответствующие паспорту специальности 05.13.18 - «Математическое моделирование, численные методы и комплексы программ» по физико-математическим наукам:

1. Новые аналитические решения для трехмерной модели движения жидкости с учетом бокового обмена в водоемах простейшей геометрии (п. 2 паспорта специальности — Развитие качественных и приближенных аналитических методов исследования математических моделей).

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

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

специальности - Разработка систем компьютерного и имитационного моделирования).

4. Результаты комплексного исследования гидрофизических характеристик оз. Шира с применением математического моделирования, численного экспермента и натурных наблюдений {п. 5 паспорта специальности - Комплексные исследования научных и технических проблем с применением современной технологии математического моделирования и вычислительного эксперимента).

Научная новизна выносимых на защиту результатов заключается в том, что:

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

2. В рамках технологии вычислительного эксперимента построена оригинальная комплексная система для трехмерного моделирования течений в соленых замкнутых стратифицированных озерах на основе свободно распространяемого пакета ОЕТМ.

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

4. На основании полного анализа гидрофизических характеристик определен характер течений в озере Шира в летний период, что дает новые знания о процессах перемешивания водной среды.

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

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

Обоснованность и достоверность полученных результатов подтверждается строгостью математических выкладок; количественным и каче-

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

Использование результатов диссертации. Результаты, полученные в диссертации, были использованы при выполнении гранта РФФИ № 07 - 01 -00153 — а, № 13 — 05 — 00853 и междисциплинарных интеграционных проектов СО РАН № 95 - 2009, 56 - 2012.

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

Представление результатов диссертации. Результаты работы докладывались и обсуждались на следующих конференциях: Международной конференции «Вычислительные и информационные технологии в науке, технике и образовании» (Усть - Каменогорск, Казахстан, 2003); Международной конференции «Вычислительные и информационные технологии в науке, технике и образовании» (Казахстан 2004); IX Всероссийской конференции «Современные методы математического моделирования природных и антропогенных катастроф» (Барнаул, 2007); международной конференции «Современные проблемы математического моделирования и вычислительных технологий -2008» (Красноярск, 2008); Всероссийской конференции «Новые математические модели механики сплошных сред: построение и изучение», приуроченной к 90-летию академика JI. В. Овсянникова (Новосибирск, 2009); Международной конференции «Mathematical and Informational Technologies» MIT-2009, 2011 (Копаоник, Сербия, 2009; Врнячка Баня, Сербия, 2011); X, XI Всероссийской конференции «Прикладные технологии гидроакустики и гидрофизики» ГА-2010, ГА-2012 (Санкт - Петербург, 2010, 2012); XV, XVI, XVII Международной научно-практической конференции «Решетневские чтения» (Красноярск, 2011, 2012, 2013); VII, VIII Всесибирском конгрессе женщин - математиков (Красноярск, 2012, 2014), Международной конференции «Информационно- вычислительные технологии и математическое моделирование (ИВТ&ММ)» (Кемерово, 2013), Международной научной конференции «Методы создания, исследования и идентификации математических моделей» (Новосибирск, 2013).

Публикации. По теме диссертации опубликовано 24 работы, в том числе (в скобках в числителе указан общий объем этого типа публикаций в печатных листах, в знаменателе — объем, принадлежащий лично автору): 2 монографии (16,25/4,06), 8 статей в периодических изданиях, рекомендованных ВАК для представления основных научных результатов диссертаций на соискание ученой степени доктора или кандидата наук (5,0/1,5), 2 - в рецензируемых журналах (1,18/0,52), 8 работ в трудах международных и всероссийских конференций (2,13/1,01), 1 препринт (3/0,33), 9 из них опубликова-

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

Личный вклад автора. Созданная система имитационного моделирования трехмерных течений в соленых замкнутых стратифицированных озерах является итогом исследований, проводимых в ИВМ СО РАН (г. Красноярск), ИБФ СО РАН (г. Красноярск), АО ИО РАН им. Ширшова (г. Калининград).

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

В основе системы вычислительного моделирования лежит пакет GETM, разработанный тремя авторами: Hans Burchard, Karsten Bolding and Lars Umlauf и распространяемый под GNU лицензией. Личный вклад автора заключается в калибровке численной модели с помощью данных натурного эксперимента; в проведении компьютерного моделирования ветровых течений как для водоемов простейшей формы, так и для водоема с батиметрией озера Шира; в интерпретации результатов численного моделирования; в создании набора программ, который позволяет формировать входные файлы, а также извлекать и обрабатывать выходные данные (написанных автором в соавторстве с О. Э. Якубайликом). Все результаты неоднократно обсуждались с Л. А. Компаниец.

Измерения скоростей течений в озере Шира в летний период с использованием допплеровских профилографов проводились как автором лично, так и совместно с О. С. Володько, В. И. Барановым. Измерения других гидрофизических характеристик озера проводились совместно с А. П. Толомеевым. Личный вклад автора состоит в проведении измерений и обработке результатов измерений трехмерного вектора скорости течений, а также обработке результатов измерений температуры и солености.

Структура диссертации. Диссертация состоит из введения, трех глав, заключения и приложения. Работа состоит из 122 страниц основного текста, включая 72 рисунка. Список использованных источников включает 101 наименование.

Содержание работы

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

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

В параграфе 1.1 проведен обзор моделей для определения ветрового течения в замкнутых водоемах. Обоснован выбор одной из них, взятой в качестве базовой для моделирования ветрового течения в различных водоемах (морях, водохранилищах, озерах).1

Это трехмерная математическая модель, описывающая движение жидкости в водоеме в приближениях Буссинеска и гидростатики

ди ди ди ди дt дх ду дг

д I ь,ди\ д 1 „ ди\ д / ди\ а д Г

г

ду дv ду ду дЬ дх ду дг

д (д / „дгл д г „ду\ а д Г

г

ди ду дю дх ду дг

дТ дТ дТ дТ д ( дТ\ д ( „дТ\ д ( тдТ\ дБ дБ дБ дБ д ( д ( д / сЭ5ч

Здесь и, V и и» являются компонентами скорости относительно направлений х, у и г соответственно. Вертикальная координата г находится в диапазоне от дна —Н(х,у) до поверхности ^(Г, х,у), t - время, К^ - коэффициент вертикального турбулентного обмена, I - параметр Кориолиса, д - ускорение свободного падения; р - плотность и р0 - средняя плотность. Горизонтальное перемешивание параметризуется членами, содержащими коэффициент горизонтального турбулентного обмена К", Т - температура, 5 - соленость, К?, К^, — коэффициенты вертикальной и горизонтальной турбулентной диффузии.

(1)

1 Марчук, Г. И. Математическое моделирование циркуляции океана / Г. И. Марчук, А. С. Саркисян. - М.:Наука, 1988.

Система уравнений дополнена граничными условиями. На свободной поверхности при г = х,у) ставятся следующие условия:

пКм^_тХ пКм^__ у Э£ 3£ Э£ дт_ дз_

• РоК* ' ™ ~ аг + идх + У; ду' Эг_и' дг

Поверхностные напряжения т* и ту вычисляются как функции скорости и направления ветра.

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

В общем случае система замыкается при помощи уравнения состояния:

Р = Р(Т,Б).

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

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

Решение в комплексной форме для прямоугольного бассейна постоянной глубины имеет вид

^ = £™=1£т=1(с1 п.те^г + с2пте~^г + сЗптс1п,т) зт™х зт^у, (2)

ап,т ^гм

Здесь = и + аиЬ - длина и ширина бассейна, ( - мнимая единица, с1п т, с2п т, сЗп т, <1п т - числовые коэффициенты, вычисленные из системы уравнений, полученной подстановкой (2) в уравнения движения, записанные в комплексном виде, и в уравнение для свободной поверхности с учетом граничных условий. Примеры полученного при таких условиях решения показаны на рисунке 1.

Приведены полученные аналогичным образом решения для прямо-

угольного бассейна с ровным дном, движение в котором возбуждается ветром с напряжением трения, задаваемым формулами: тх/р0 = —у, ту/р0 = х.

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

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

Рисунок 1 - Годограф скорости а) в центре бассейна и б) вблизи берега. Непрерывной линией показано решение по модели Экмана, ромбиками - по модели с учетом бокового обмена

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

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

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

элементные), реализации численного алгоритма. Отдельные варианты включают в себя возможность исследования биологических процессов. Большинство из них находится в свободном доступе в сети Internet, запрещено для коммерческого использования, но открыто для исследовательских целей. Программу GETM (General Estuarine Transport Model) попытались сделать так, чтобы она была наиболее общей, то есть позволяла рассчитывать большое разнообразие течений в водных объектах, включая осушение и затопление. Программа GOTM (General Ocean Turbulence Model) является составной частью GETM и предназначена для расчета вертикального перемешивания.

Она использовалась в работах европейских ученых для расчетов стратификации во фьордах Норвегии, для расчета переноса примесей в северозападной части Северного моря и т. д., в работах российских исследователей для расчета водоемов-охладителей Шатурской ГРЭС.

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

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

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

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

В параграфе 3.1 разработана структура программно-технологического

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

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

Комплексная система моделирования гидрофизических процессов в соленых стратифицированных замкнутых озерах

Аналитические Вычислительный эксперимент Натурные наблюдения

Рисунок 2 - Схема комплексной системы моделирования

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

—Т 02.07.16м

• Т 03.07.22 м Т 04.07.24 м

—Т 04 07 13 м —Т 0607.4 м

- Т 06.07.22 м —Т 07.07.5 м —S 02.07.16 м

S 03.07.22 м

• Ь 04.07. 24 и

• S 04.07. 13 м

- S 06.07.4м S 06.07.22 м S 07.07.5 и

Рисунок 4 - Типичное распределение температуры и солености (слева), плотности воды (справа), в зависимости от глубины, в летний период на оз. Шира

В параграфе 3.2 кратко сформулированы методики проведения натурных экспериментов, которые выполнены на озере Шира в 2009 - 2013 гг. Представлены результаты обработки данных натурных наблюдений, в том числе авторские.

Автором построена цифровая модель рельефа, которая приведена в виде карты изолиний (рисунок 3, слева). Обработаны данные измерений температуры и кондуктивности. Пример результатов обработки представлен на рисунке 4, на котором слева изображен график распределения температуры и солености в различных точках озера, типичного в летний период. Разным цветом показано измерение температуры и солености в точках с различной глубиной в период со второго по седьмое июля 2011 г.

Рисунок 3 - Цифровая модель рельефа оз. Шира (слева), картина установки датчиков на озере в июне - июле 2013 г. (справа)

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

Картина колебания изотерм на рисунке 5 (типичная для течения в присутствии внутренних волн") подтверждает возможность их возникновения в озере Шира при определенных метеорологических условиях.

Температура. °С. Соленость, Пп.

2 Hutter, К. Physics of lakes / К. Hutter, Y. Wang, I.P. Chubarenko. - Springer-Verlag Berlin Heidelberg, 2011.-V 1.

Лиги [9.06.13 20.06.13

Промя 05:32 23>17

Рисунок 5 — Данные с термисторных распределенных датчиков в точках Т1 (слева) и ТЗ (справа), изображенных на рисунке 3

В параграфе 3.3 выполнен вычислительный эксперимент на основе комплексной системы для конкретного уникального водоема (озера Шира). Приведены результаты численных расчетов при постоянном ветре и проведено сравнение с натурными данными. Начальное распределение температуры и солености соответствовало реальному распределению этих величин в оз. Шира. При постоянном южном ветре 5 м/с в расчетах были получены следующие результаты.

Показано, что на расчетный период времени, даже при наличии горизонтальных скоростей на дне, нижний и верхний слои не перемешиваются, что иллюстрируют графики температуры на рисунке 6. Как и в натурных наблюдениях, в численных расчетах горизонтальные скорости два раза меняют знак в распределении по глубине (рисунки 7 и 9, рисунки 8 и 10).

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

Рисунок 6 - Характерные профили распределения температуры по глубине в центральной точке бассейна в начале и в конце расчета

Дпта 19.06.13 Нремя 05:32

Рисунок 7 - Профили скоростей в точке глубиной 16 м, полученные в результате расчетов

Рисунок 8 - Профили скоростей в точке глубиной 5,8 м. полученные в результате расчетов

-»0 -20 20 40 (Л -100 .so 50 100

и, ММ/С .___ V, «м/с

г>

/у -•-1200

и

У

-Я, м J --M-J •Н. м --Н—

Рисунок 9 - Измерения 03.07.10 г. Точка имеет широту 54°30.392', долготу 90°09.249', глубина16 м. ADCP 1200 установлен излучателями вниз. ADCP 600 установлен излучателями вверх

Рисунок 10 - Измерения 04.07.10 г. Точка имеет широту 54°30.295\ долготу 90°15.195', глубина 5.5 м. АЭСР 1200 установлен излучателями вниз, АЭСР 600 установлен излучателями вверх

Пример обработки результатов одного из вычислительных экспериментов приведен на рисунке 11, когда расчет начинался при стационарном ветре 6 м/с в течение 6 часов, затем 1 час при ветре 11 м/с, далее расчет

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

Рисунок 11 - Изотермы, полученные в расчетной точке, ближайшей к Т1 (слева) и к ТЗ (справа) (рисунок 5) в бассейне с батиметрией оз. Шира при переменном ветре

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

В заключении перечислены основные результаты:

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

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

3. Комплексное исследование гидрофизических характеристик озера Шира с применением средств математического моделирования и технологии вычислительного эксперимента позволило: обосновать возможность появления внутренних волн в озере Шира в летний период; выявить существенно трехмерный характер течения в озере с образованием при определенных метеорологических условиях горизонтальных вихревых структур.

В приложении описаны приборы и методики проведения натурных измерений, результаты обработки полученных данных. Изложена технология уточнения батиметрии озера при помощи эхолота и допплеровских профило-графов в режиме «bottom track». Представлены данные трехмерного вектора скорости течения, полученные с допплеровских профилографов по поперечным разрезам озера Шира. Построены графики точечных измерений темпе-

ратуры и кондуктивности, выполненные при помощи CTD (conductivity, temperature, depth) зондов. Отображены результаты измерений CTD зондом IDRONAUT по галсам с системой промерных станций на расстоянии 500 м для получения информации о температуре и солености по поперечным разрезам озера. Приведены результаты обработки измерений температурных датчиков.

Публикации автора по теме диссертации

Монографии:

1. Модели экмановского типа в задачах гидродинамики / Л. А. Компаниец, Т. В. Якубайлик, Л. В. Гаврилова, К. Ю. Гуревич - Новосибирск: Наука. - 2007.

2. Аналитические решения для задач стационарного ветрового движения жидкости / Л. А. Компаниец, Т. В. Якубайлик, Л. В. Гаврилова,

0. С. Володько - Красноярск: Изд-во СФУ. - 2012.

В рецензируемых журналах, рекомендуемых ВАК:

1. Компаниец, Л. А. Аналитическое решение одной модели ветрового движения жидкости / Л. А. Компаниец, Т. В. Якубайлик // Вычислительные технологии,- 2003. - Т.8. - № 5. - С. 78-83.

2. Белолипецкий, В. М. Математические модели и численные алгоритмы для исследования задач гидрофизики / В. М. Белолипецкий, С. Н. Генова, Л. А. Компаниец, Т. В. Якубайлик // Вычислительные технологии,- 2004-Т.9. - (Специальный выпуск, посвященный 30-летию Института вычислительного моделирования Сибирского отделения РАН) - С. 22 - 28.

3. Компаниец, Л. А. Сравнительный анализ двух моделей двухслойной жидкости в приближении Экмана / Л. А. Компаниец, Т. В. Якубайлик, К. Ю. Гуревич, Л. В. Гаврилова, Е. А. Кирилюк // Журнал Сибирского федерального университета.-2008.-Т. 1. -№2. - С. 197-209.

4. Компаниец, Л. А. Аналитическое решение одной модели ветрового течения вязкой жидкости / Л. А. Компаниец, О. С. Питальская, Т. В. Якубайлик // Вычислительные технологии. - 2009. - Т. 14. - № 4. - С. 46 -57.

5. Компаниец, Л. А. О модели ветрового движения двухслойной вязкой жидкости / Л. А. Компаниец, Л. В. Гаврилова, Т. В. Якубайлик // Компьютерные исследования и моделирование. - 2009. - Т. 1. - №4. - С. 381 -390.

6. Компаниец, Л. А. Аналитическое решение одной модели ветрового движения вязкой жидкости (трехмерный случай) / Л. А. Компаниец, Т. В. Якубайлик, О. С. Питальская // Вычислительные технологии - 2011. -Т.16. -№ 3. - С. 50-63.

7. Компаниец, Л. А. Анализ характеристик озера Шира на основе натурных наблюдений / Л. А. Компаниец, Т. В. Якубайлик, О. С. Володько // Вестник Бурятского государственного университета: математика, информатика. -Улан-Удэ: Изд-во БГУ. - 2012. - выпуск 9. - С. 167 - 176.

8. Баранов, В. И. Пространственно-временная изменчивость основных характеристик озера Шира в сезоне 2911-2012 гг. / В.И.Баранов,

H. Н. Голенко, Л. А. Компанией, В. Т. Пака, Т. В. Якубайлик // Вестник Бурятского государственного университета. Математика. Информатика. - Т.9. -2013.-С. 148-156.

В рецензируемых журналах:

I. Гаврилова, Л. В. Аналитическое решение одной модели движения однородной жидкости в мелком водоеме (2D и 3D случай) / Л. В. Гаврилова, Л. А. Компанией, Т. В. Якубайлик // Вычислительные технологии. - 2003. -Т. 8. - Региональный вестник Востока. - 2003. - № 3. - (Совм. выпуск по материалам Международной конференции «Вычислительные и информационные технологии в науке, технике и образовании». Ч. 4.). -С. 125-131.

2. Компанией, Л. А. Аналитическое решение одной модели ветрового движения двухслойной жидкости./ Л. А. Компанией, Т. В. Якубайлик // Вычислительные технологии (совместный выпуск с Вестник КазНУ, №3 (42), 2004г.). - Алматы. - Новосибирск. - 2004. - Т.9. - Часть II. - С. 372 - 383.

В трудах конференций:

1. Yakubailik, Т. V. On some analtical solutions for the model of the wind-induced motion of the viscous incompressible liquid (the case of three dimensions) / Т. V. Yakubailik, L. A. Kompaniets // Mathematical and Informational Technologies - MIT 2009 (August, 27- 31, 2009, Kopaonik, Serbia; August, 31 - September, 5, 2009, Budva, Montenegro) // Conference Information, Univerzitet u PriSti-ni, Prirodno-matematicki fakultet,KosovskaMitrovica.-2009.- P. 113-114.

2. Компанией, Л. А. Опыт использования современных измерительных приборов для определения гидродинамических режимов водоема / Л. А. Компаниец, Т. В. Якубайлик, О. С. Питальская // Труды X Всероссийской конференции «Прикладные технологии гидроакустики и гидрофизики» ГА-2010. - СПб: Наука. - 2010. - С. 310 - 312.

3. Компаниец, Л. А. Численное моделирование гидрофизических процессов в озере Шира в летний период / Л. А. Компаниец, Т. В. Якубайлик // Zbornik radova Konferencije MIT 2011. - Beograd: Alfa univerzitet; Kosovska Mitrovica: Druätvo matematicara Kosova i Metohije, Serbia; Novosibirsk: Institute of Computational Technologies, Siberian Branch of the Russian Academy of Sciences, 2012 (Kraljevo: Graficolor). - 404 str. - ISBN 978-86-83237-90-6 (AU). -P. 198-202.

4. Володько, О. С. Исследование гидротермических режимов непроточного стратифицированного водоема с использованием допплеровских про-филографов / О. С. Володько, Л. А. Компаниец, Т. В. Якубайлик // Труды XI Всероссийской конференции «Прикладные технологии гидроакустики и гидрофизики» ГА-2012 - СПб: Наука. - 2012. - С. 271 - 274.

5. Якубайлик, Т. В. Теоретическое обоснование возможности появления внутренних волн в оз. Шира / Т. В. Якубайлик, Ю. Б. Гульденбальк, Л. А. Компаниец // Избранные труды международной конференции по изме-

рениям, моделированию и информационным системам для изучения окружающей среды: ЕЫУИЮМга. - 24 июня-2 июля. - Иркутск. - 2012. - С. 104107.

6. Гульденбальк, Ю. Б. Анализ натурных данных для одного из сценариев течения в озере Шира в летний период / Ю. Б. Гульденбальк, Л. А. Компанией, Т. В. Якубайлик // Актуальные проблемы науки и образования: прошлое, настоящее, будущее: Сборник научных трудов по материалам Международной заочной научно-практической конференции 29 марта 2012 г. - Тамбов. - 2012. - С 59 - 64.

7. Якубайлик, Т. В. Калибровка численной модели для расчета гидрофизических характеристик озера Шира с использованием натурных данных // Материалы конференции «VII всесибирский конгресс женщин - математиков». - Красноярск: СибГТУ. - 2012. - С. 245 - 249.

8. Якубайлик, Т. В. Численное моделирование поведения внутренних волн в прибрежной зоне неглубокого стратифицированного водоема / Т. В. Якубайлик, Л. А. Компанией И Теоретические и прикладные вопросы науки и образования. Сборник научных трудов по материалам Международной научно-практической конференции 31 августа 2013 г. - Тамбов: Изд-во ТРОО.-2013.-С. 54-58.

Местные издания:

1. Компанией, Л. А. Теоретические и экспериментальные исследования скоростного режима озера Шира в летний период / Л. В. Гаврилова, Л. А. Компаниец, Т. В. Якубайлик // Препринт ИВМ СО РАН. - Красноярск. -2012.-№12.-3.

Свидетельства о государственной регистрации программ:

1. Якубайлик О. Э., Якубайлик Т. В. Программа для систематизации полей температуры по результатам численного моделирования замкнутых соленых стратифицированных озер // Федеральная служба по интеллектуальной собственности (Роспатент). Свидетельство о государственной регистрации программы для ЭВМ № 2014610077 от 9 января 2014 г.

2. Якубайлик О. Э., Якубайлик Т. В. Программа формирования данных батиметрии для численного моделирования течений в соленых замкнутых стратифицированных озерах // Федеральная служба по интеллектуальной собственности (Роспатент). Свидетельство о государственной регистрации программы для ЭВМ № 2014610113 от 9 января 2014 г.

3. Якубайлик О. Э., Якубайлик Т. В. Программный комплекс обработки данных результатов численного моделирования замкнутых соленых стратифицированных озер // Федеральная служба по интеллектуальной собственности (Роспатент). Свидетельство о государственной регистрации программы для ЭВМ № 2014610060 от 9 января 2014 г.

Подписано в печать 06.06.2014. Формат 60*84/16 Усл. псч. л. 1. Тираж 100 экз. Отпечатано на ризографе ИВМ СО РАН 660036, Красноярск, Академгородок, 50/44

Текст работы Якубайлик, Татьяна Валерьевна, диссертация по теме Математическое моделирование, численные методы и комплексы программ

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ УЧРЕЖДЕНИЕ НАУКИ ИНСТИТУТ ВЫЧИСЛИТЕЛЬНОГО МОДЕЛИРОВАНИЯ СИБИРСКОГО ОТДЕЛЕНИЯ РОССИЙСКОЙ АКАДЕМИИ НАУК

На правах рукописи 0420146061Ь 0 г ^

Якубайлик Татьяна Валерьевна

Комплексная система моделирования гидрофизических характеристик замкнутых соленых стратифицированных озер (на примере оз. Шира)

05.13.18 — Математическое моделирование, численные методы и комплексы

программ

Диссертация на соискание ученой степени кандидата физико-математических наук

Научный руководитель

к.ф.-м.н., доцент Компаниец Лидия Алексеевна

Красноярск — 2014

Оглавление

Введение..............................................................4

Глава 1. Математические модели ветровых течений в стратифицированных водоемах. Аналитические решения упрощенных моделей.................................................... 20

1.1 Основные уравнения ветровых течений жидкости в замкнутых стратифицированных водоемах........................................ 20

1.2 Аналитические решения трехмерных уравнений стационарного ветрового течения ..................................................... 25

1.2.1 Однородная жидкость ............................................. 25

1.2.2 Двухслойная жидкость ............................................ 38

Глава 2. Численная модель ветровых течений

в стратифицированных водоемах ..................................... 47

2.1 Обзор численных моделей для расчета ветровых течений

в стратифицированных водоемах ...................................... 47

2.2 Комплекс программ СЕТМ (СОТМ) .................................. 48

2.2.1 Алгоритм для расчета уравнений движения жидкости ........... 49

2.2.2 Алгоритм для расчета полей температуры и солености ........... 54

2.2.3 Определение коэффициентов горизонтального турбулентного обмена и горизонтальной турбулентной диффузии.

Уравнение состояния .............................................. 55

2.2.4 Формат ввода и вывода данных ................................... 55

2.3 Адаптация и верификация численного алгоритма для расчета течений

в неглубоких замкнутых стратифицированных водоемах ............. 57

2.3.1 Выбор сеток на основе тестовых расчетов ........................ 57

2.3.2 Выбор разностных аппроксимаций членов переноса .............. 58

2.3.3 Выбор уравнения состояния ....................................... 59

2.3.4 Тестовый расчет: прямоугольный бассейн с ровным дном ........ 61

2.3.5 Тестовый расчет: цилиндрический бассейн с ровным дном ....... 67

2.3.6 Тестовый расчет: бассейн с ровным дном и береговой линией реального озера.................................................... 71

2.3.7 Тестовый расчет: параболический бассейн ........................ 75

Глава 3. Комплексная система моделирования гидрофизических процессов в соленых стратифицированных замкнутых озерах. Применение к моделированию гидрофизического режима озера Шира............................................................... 81

3.1 Структура комплексной системы ...................................... 81

3.2 Общее описание натурного эксперимента..............................84

3.3 Тестовый расчет: водоем с батиметрией озера Шира, однородная жидкость............................................................... 87

3.4 Численные расчеты для озера Шира при постоянном ветре. Сравнение с натурным экспериментом ................................ 90

3.5 Численные расчеты для озера Шира при переменном ветре. Сравнение с натурным экспериментом ................................ 94

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

Список литературы ............................................... 98

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

1 Технология измерения батиметрии ................................... 111

2 Технология измерения скорости течения ............................. 114

3 Технология измерения температуры ................................ 118

Введение

Актуальность. Соленые озера являются важными природными объектами, в которых заключены разнообразные запасы минеральных и биологических ресурсов. Они представляют собой большую экономическую ценность, так как являются одним из основных мест рекреационного туризма, медицинским ресурсом по запасам лечебной грязи и соли. Мировым сообществом признана необходимость «изучения, рационального использования и сохранения природы соленых озер» [24]. Современные соленые озера имеют научное значение как природные седиментологические лаборатории, в которых не только записаны все палеоклиматические изменения, но они являются и природной площадкой для изучения солеобразующих процессов и изменений окружающей среды.

Экологическое состояние водных объектов зависит от большого числа разнообразных факторов и процессов: гидрофизических, гидробиологических, гидрохимических, метеорологических и антропогенных. Гидрофизические процессы в значительной мере формируют среду обитания гидробионтов, определяют перенос и седиментацию веществ, интенсивность загрязнения и самоочищения водоёмов

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

Получили свое развитие и дистанционные аэрокосмические методы исследования озер. Использование дистанционных материалов для контроля со-

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

Теоретические исследования в гидрофизике и, в частности, в лимнологии (в построении и исследовании математических моделей для однородной и неоднородной жидкости) представлены в целом ряде монографий российских и зарубежных ученых: [19, 40, 41, 42, 44, 45, 49, 70]. В конце XX, начале XXI века были созданы первые модели соленых озер.

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

Численные модели для расчета течений и температурного режима в озерах, в том числе соленых, начали создаваться в 60-х годах XX века. Так, в [26] предложена трехмерная численная модель для Каспийского моря. В [10, 11, 12] рассматривались малоразмерные численные модели для озера Шира. В [43] предложена численная реализация модели Экмана. В 70-е годы были разработаны первые математические модели экосистем для озер Байкал, Севан, Ладожского и Чудского, американских Великих озер. Они изложены, например, в работах [1, 3, 6, 15, 70].

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

дерландах, лаборатория Кориолис во Франции и т. д.). В последние сорок — пятьдесят лет начали уделять значительное внимание лабораторным исследованиям нелинейных внутренних волн [14].

Измерения скорости течения проводились, начиная с прошлого века, с использованием простейших ультразвуковых измерителей скорости. Результаты таких измерений отражены, например, в монографии [47]. С развитием измерительной техники появилась возможность проводить достаточно точные натурные наблюдения гидрофизических характеристик не только крупных водных объектов (таких как моря, заливы, крупные соленые и пресные озера), но и более мелких водоемов (таких как внутренние соленые и пресные озера диаметром до 10 км и небольшой глубины), в которых измеряемые величины (в частности, скорости течений) могут быть меньше в десятки раз. Сейчас во всем мире широко применяются допплеровские профилографы для измерения трехмерного вектора скорости течения по всей глубине водного объекта. В нашей стране их использовали, в частности, для исследования течений в Балтийском, Черном, Печерском морях.

Исследования гидрофизики озер в настоящее время проводятся по нескольким направлениям: нахождение аналитических решений упрощенных моделей, численное моделирование, лабораторные и натурные эксперименты. Причем, большинство исследований относится к пресным озерам (Кинерет [67], Боденское [70], Ладожское [1, 3, 6], Онежское [23], Чудское [15], Севан [1] Байкал [50]) или к соленым крупным озерам (Каспийское [26], Азовское [52], Мертвое моря [69]). В то же время, исследованию неглубоких замкнутых стратифицированных озер уделялось недостаточно внимания.

Научные исследования на оз. Шира начали проводить в 50-е годы XX века. Экспедиция МГУ впервые провела измерения некоторых биологических, химических и гидрофизических характеристик озера. Это было отправной точкой для перехода к широким исследованиям в озере. В последующие годы была уточнена формула для определения солености воды. Была построена одномерная модель для распространения химических, биологических элементов и их

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

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

Цель работы состоит в создании на основе технологии вычислительного эксперимента комплексной системы моделирования соленых стратифицированных озер, ориентированной на изучение их гидрофизических характеристик и применении ее для изучения ветровых течений в реальном водоеме (озере Шира, расположенном в республике Хакасия).

В соответствии с это целью в работе поставлены и решаются следующие задачи:

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

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

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

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

5. Дать оценку гидрофизического режима озера Шира в летний период, используя выявленные закономерности.

Объектом исследований являются физические процессы, специфические для замкнутых стратифицированных соленых водоемов.

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

На защиту выносятся следующие результаты, соответствующие паспорту специальности 05.13.18 — «Математическое моделирование, численные методы и комплексы программ» по физико-математическим наукам:

1. Новые аналитические решения для трехмерной модели движения жидкости с учетом бокового обмена в водоемах простейшей геометрии (п. 2 паспорта специальности — Развитие качественных и приближенных аналитических методов исследования математических моделей).

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

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

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

4. Результаты комплексного исследования гидрофизических характеристик оз. Шира с применением математического моделирования, численного эксперимента и натурных наблюдений (п. 5 паспорта специальности — Комплексные исследования научных и технических проблем с применением современной технологии математического моделирования и вычислительного эксперимента).

Научная новизна выносимых на защиту результатов заключается в том, что:

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

2. В рамках технологии вычислительного эксперимента построена оригинальная комплексная система для трехмерного моделирования течений в соленых замкнутых стратифицированных озерах на основе свободно распространяемого пакета СЕТМ.

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

4. На основании полного анализа гидрофизических характеристик определен характер течений в оз. Шира в летний период, что дает новые знания о процессах перемешивания водной среды.

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

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

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

Использование результатов диссертации. Результаты, полученные в диссертации, были использованы при выполнении гранта РФФИ № 07—01— 00153-а, № 13-05-00853 и междисциплинарных интеграционных прое