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

005002886

Проценко Елена Анатольевна

МОДЕЛИРОВАНИЕ ТРАНСПОРТА НАНОСОВ В ПРИБРЕЖНОЙ ЗОНЕ ВОДОЕМА

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

АВТОРЕФЕРАТ

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

-1 ДЕК 2011

Таганрог-2011

005002886

Работа выполнена в Технологическом институте Южного федерального университета в г. Таганроге (ТТИ ЮФУ).

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

доктор физико-математических наук, профессор Сухинов Александр Иванович

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

доктор физико-математических наук, профессор Куповых Геннадий Владимирович

Защита состоится 29 декабря 2011 г. в 10.20 на заседании диссертационного совета Д 212.208.22 в Технологическом институте Южного федерального университета в г. Таганроге по адресу:

347928, Ростовская обл., г. Таганрог ГСП-17А, пер. Некрасовский 44, ауд. Д-406.

С диссертацией можно ознакомиться в зональной научной библиотеке ЮФУ по адресу: г. Ростов-на-Дону, ул. Пушкинская, 148.

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

347928, Ростовская обл., г. Таганрог ГСП-17А, пер. Некрасовский 44.

доктор физико-математических наук, профессор Чикин Алексей Львович

Ведущая организация: Северо-Кавказский государственный

технический университет

Автореферат разослан ноября 2011 г.

Ученый секретарь

диссертационного совета Д.212.208.22 доктор технических наук, профессор

А.Н. Целых

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

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

Проведенный анализ исследований по рассматриваемой проблеме показал, что:

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

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

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

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

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

1. Проведено исследование модели переноса многокомпонентной примеси в водной среде. Доказано достаточное условие единственности решения соответствующей начально-краевой задачи.

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

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

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

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

Научная новизна диссертационного исследования.

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

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

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

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

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

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

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

Программная реализация построенных численных алгоритмов базируется на языке «С++». Визуализация проводилась в среде разработки МАТНСАБ.

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

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

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

В рамках сформулированной в работе проблемы на защиту выносятся следующие результаты и положения:

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

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

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

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

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

прикладной математики в технике, медицине и экономике (Новочеркасск, 2 февраля 2007); Пятьдесят третья научно-теоретическая конференция ГОУВПО «ТГПИ». Секция математики. (Таганрог, 22 января 2010); Пятьдесят четвертая научно-теоретическая конференция ГОУВПО «ТГПИ». Секция математики. (Таганрог, 21 января 2011); Международная научно-практическая конференция: Молодежь и наука: реальность и будущее. (Невинномысск, НИЭУП, 3 марта 2010); Научный семинар кафедры ВМ ТТИ ЮФУ (Таганрог, 2011).

Структура и объем работы. Работа состоит из введения, 4 глав, заключения, списка литературы, содержащего 151 наименований. Работа содержит 33 рисунка, 2 таблицы. Полный объем диссертации составляет 196 страниц.

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

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

Глава I посвящена анализу моделей транспорта наносов и многокомпонентных примесей.

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

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

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

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

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

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

краевая задача имеет вид: (1-е)—+ = О, (2.1) О. = Ат1ц/Ч-, (2.2)

81 дх

где е - пористость грунта; Я - глубина дна, отсчитываемая от невозмущенного уровня водоема; <2 - расход наносов; у/в - параметр Шильдса для наклонного дна; Ли р - безразмерные постоянные (Л=!9,5; /?=3); йт - частота волны; / - время.

У* =0"«,±г*8т5/8тя,)/(д- д,)®*, гь =-0,6СС4|У4|, гУ„ = !/„((2.3)

где ть - касательное напряжение на дне; ты - критическое значение касательного напряжения, при котором начинается перемещение наносов; 5 - угол, составленный касательной к контуру дна в момент времени /; сра - угол естественного откоса грунта в воде; р,,р0 - плотности твердых частиц и воды; С - коэффициент сопротивления твердых частиц; 11ь - придонная скорость в области неразрушенных волн; с1- характеристика осадков, ИЬс - скорость на критической глубине Н; х„ - координата верхней границы наката; х0 - координата точки опрокидывания волн; п - постоянная (в данной работе п=0,33).

V /">,ч дН д( дНЛ дЬ „ Уравнение (2.1) представлено в виде:---а- + — = 0, (2.4)

3/ дх \ дх) дх где л = Ь = кть-, к = ^

Начальное условие для уравнения неразрывности: Н{х) = Бах при 1 = 1а. (2.6) Граничные условия определяются из физических соображений. В верхней границе наката, берег не подвергается деформациям: Н(х1!,1)=Н1Пхн =Н„/ЯЬ. (2.7) На границе «глубокой» воды глубина также не меняется: Н(хгл,1) = Н,,„ Нгл =;т1/2. ■ (2.8)

В разделе 2.4 на основе интегро-интерполяционного метода выполнено построение дискретной математической модели поверхностных волн от начальных возмущений по непрерывному аналогу. В результате построена однородная трехслойная разностная схема. Разделы 2.5 ~ 2.6 посвящены исследованию устойчивости и погрешности аппроксимации построенной разностной схемы.

В разделе 2.7 получена аппроксимация математической модели волновых гидродинамических процессов. Использована гидродинамическая модель, исходными уравнениями которой являются: одномерный аналог системы уравнений Навье-

Стокса: ((#+<?)*)' +(Н+£)ии'х (2.9)

* Р Р

и уравнение неразрывности для несжимаемой жидкости: £'+((//+£)«)' =0, (2.10) где £ - функция подъема уровня, и - скорость движения водной среды, ц - коэффициент турбулентного обмена, g - ускорение свободного падения, р - плотность жидкости, тр,ть -тангенциальное напряжение на поверхности и дне водоема.

Для аппроксимации уравнения (2.9) по временной переменной применены

схемы расщепления по физическим процессам.

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

Глава III посвящена дискретизации одномерной математической модели процессов перемещения наносов в прибрежной зоне и ее численной реализации.

Поставленная для уравнения параболического типа (2.4) - (2.5) начально-краевая задача решается приближенно конечно-разностным методом.

Введена пространственная сетка: /=1,2,..,^; +п-т, п=0Х-,Ц,

интервал (х

н'хгл) и временной интервал разбивались на равные части.

Для аппроксимации уравнения (2.4) использован метод баланса. В результате получена двухслойная неявная конечно-разностная схема:

_ И" _ ГГЛ+1 „„+1 , „ , „

= (3.1)

Доказана устойчивость конечно-разностной схемы (3.1) - (3.2). Исследование показало, что погрешность аппроксимации конечно-разностной схемы равна С\т+Иг).

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

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

На рисунках 2-3 представлена динамика изменения рельефа дна и функции возвышения уровня в случае наклонного дна. Изначально максимальная глубина водоема составляла 1 метр, минимальная - 50 см.

Рис. 2. Функция возвышения через время равное 1,5 и 2,5 периодам волны

Рис. 3. Изменение рельефа дна через время равное 1,5 и 2,5 периодам волны

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

Глава IV посвящена разработке двумерной модели транспорта наносов в мелководных водоемах.

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

приведена к виду: () = Атс1

- gradH

{(А-Р^ау где 0 ={&,£>,} - расход наносов, \q\~Q) ть

-gradH

(4-1)

япрв

■ касательное напряжение на дне; тЬс -

критическое значение касательного напряжения, при котором начинается перемещение наносов, А и /3 - безразмерные постоянные {А-Х9.5, /3=3), ш - частота волны, с1 - характеристика осадков, / - время; <р0 -угол естественного откоса грунта в воде. Уравнения процесса перемещения наносов приведены к виду:

(1 -е)~+ <№г(*гЛ= к-^-ёгас1Н д1 V "»Ро

с учетом ограничений на касательные напряжения на дне расчетной области:

д-1

к=-

Атлй

((я-Ро)^У

эш срй

-•¿гас1Н

-г»

(4.2)

(4.3)

где Я - глубина дна, отсчитываемая от невозмущенной поверхности водоема; е

- пористость грунта; х, у - горизонтальные декартовы координаты;

= - функция Хэвисайда.

[0, х<0

Уравнение (4.2) дополнено начальным условием: Н(х,у,0) = Н0(х,у). (4.4) Граничные условия определяются из физических соображений. В верхней границе наката, берег не подвергается деформациям: Н{х,у,0 = Нн. (4.5) На границе «глубокой воды» отсутствует поток, вызванный влиянием гравитационных сил: Н'0(х,у, 0 = 0.

Таким образом, построена непрерывная двумерная математическая модель формирования наносов в прибрежной зоне водоема (4.2) - (4.6).

В разделе 4.3. проведена дискретизация двумерной математической модели транспорта наносов. Использована равномерная прямоугольная расчетная сетка Ф = а>,хюххю,-.св,~{Г=пк„ 0<«<^,-1, 1, = -1)}., сох={х, =;м</ -и= (л =Л> 0 < у < А'; -1, ^ = к (л; -!)} . где п,и - индексы по временной координате и пространственным координатным направлениям Ох,Оу соответственно; - шаги; - количество

узлов; 1„1х,1у - длина расчетной области по временной координате и пространственным координатным направлениям Ох,Оу соответственно.

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

1.1-41

_ гЬс

тгп+а _ Нп+а

V» м-' 11 /г"

Яп+сг ип+сТ ^

'.1 ~-"'-и

V

51П<Э0

Яп+<т Ып*°

1г" ' -к"

К1,}+\П .2 ¡,1-1/2

Яп+сг ТТП+а А

и ~ пи-\

к

(4.7)

где сг — вес схемы

к" = 1+1/2,7

((Я -Ро)^)

МП.]

Получены следующие аппроксимации:

'¡+1/2,;

К

-I +-

2й„

I® 'и-и« К

(4.9) (4.10)

где 7,7 - единичные векторы, направленные вдоль осей Ох,Оу соответственно,

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

Разделы 4.6-4.7 посвящены исследованию устойчивости и погрешности аппроксимации разностной схемы. Устойчивость разностной схемы исследована на основе сеточного принципа максимума.

Построенные разностные схемы устойчивы при следующем ограничении;

. у\

¡\<ъ\п%{\-а)/

V \ х у;;

Оценка сверху для сеточной функции исходной задачи имеет вид:

(4.11)

(4.12)

1 ь Д1«0П 1 с

Исследован общий порядок погрешности аппроксимации математической

модели транспорта наносов: 0(Ъ, + И] + Щ}.

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

Для построения двумерной математической модели движения водной среды использована трехмерная гидростатическая модель, представленная в виде:

- системы уравнений Навье-Стокса:

г,\+ии%+ (4ЛЗ)

- уравнения неразрывности для несжимаемой жидкости: «1 + ^ + 4=0; (4.14)

- уравнения гидростатики:

Р = р&(г + £). (4.15)

Система уравнений (4.13)-(4.14) рассматривается при граничных условиях: на дне - условия непроницаемости и трения: ртц1п =^(0. РЛ^п =ТУА()> на поверхности задаются подъем уровня и ветровые напряжения: РЩ'„ = -г,.,(О, = -ту,р (/), м> = ; на боковых границах - условие свободного выхода: ип = где £ - функция подъема уровня, V = - вектор скорости движения водной

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

Для построения двумерной модели гидродинамики система уравнений (4.13) - (4.14) была проинтегрирована по глубине с учетом граничных условий. Для аппроксимации задачи по времени использован метод поправки к давлению. Для аппроксимации задачи по пространственным координатам использован ин-тегро-интерполяционный метод. Разработанная двумерная дискретная модель

гидродинамики мелководных водоемов представлена разностными схемами: - для расчета поля скорости на промежуточном временном шаге:

„+<г _ „ и+(7/2 _ п+аП

„"■нг/2 ,л+<т/2 .п-нт/2 »-нг/2

,,1+е-П „л+сг/2 пит/2_ п+а!2

. —"( (.1 /.. „„ \ Ц+1 / и_

-(и \и и

, и+ст/2 л-кг/2 2 —А,"*"11

"II ~ , -я \ 1,1

-{Н,А12,, + /4-1/2,/"- - ' + (Я, /+1/2 +#(Иу+1/2)д,

х

./МП _ »«г/2 \п1,Н12 ^ Ч/-1/2}Н-1,)АП -•• т

У

А2

У

- для расчета функции возвышения уровня:

й; р Д,

<:<н-1_ел / ет+1 _с»+1 гг"+1 Л

(#,„,+1/2 + £";+1/2) ,2 "{^1,1-41 + ^2

Ч 1 1

_ {Нц-112./ + + £-1/2,/)ЦМ/2,/

Л,

К

- для уточнения поля скорости:

п+<г_ п+о- сл+1 еп+1_еп+1

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

Следует отметить, что расчеты скорости без учета давления, функции возвышения уровня, транспорта наносов осуществлялись в два этапа:

- построение соответствующих сеточных уравнений;

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

В разделе 4.10 представлены результаты численных экспериментов.

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

Рис. 4. Начальный рельеф дна и положение свободной поверхности

На рисунке 5 показана геометрия дна и функция возвышения уровня, на рисунке 6 - изолинии глубины и скорость течения. Скорость течения жидкости максимальна в пиковой части области (на минимальных глубинах). Функция возвышения уровня принимает положительные значения в левой часта расчетной области (с наветренной стороны), и отрицательные в правой (Рис. 5).

Рис. 5. Геометрия дна и функция возвышения уровня через 20 мин

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

Рис. 6. Изолинии глубины и скорость течения через 20минут

Через 40 минут - присутствует волнообразная структура функции возвышения уровня не только в левой части, но и в области пиковых значений глубины. Через 60 минут - усиливаются колебания функции возвышения уровня в левой своей части и в области пиковых значений глубины. На дне расчетной области образовываются гряды. С наветренной стороны продолжает расширяться область с пиковыми значениями глубины (Рис. 7).

у. и

Рис. 7. Геометрия дна и функщ/я возвышения уровня через 60 минут

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

у, и

Рис. 8. Геометрия дна и функция возвышения уровня через 80минут

Рис. 9. Изолинии глубины и скорость течения через 80мин

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

у, «

Рис. 10. Геометрия дна и функция возвышения уровня через 100 минут

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

Таким образом, работоспособность метода подтверждена проведенным численным экспериментом.

В заключении изложены основные результаты и выводы.

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

Работа содержит следующие научные результаты:

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

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

3. Построена двумерная модель транспорта наносов в мелководных водоемах, включающая уравнения движения, неразрывности и перемещения наносов.

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

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

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

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

Благодарности. Автор выражает свою искреннюю благодарность научному руководителю - Сухинову Александру Ивановичу (руководителю ТТИ ЮФУ) за его за помощь и руководство.

Основное содержание диссертации изложено в следующих работах:

По теме исследования опубликовано 11 печатных работ, из них 5 статей в отечественных реферируемых журналах, входящих в список изданий, рекомендованный ВАК.

В изданиях из перечня ВАК:

1. Проценко Е.А. Модель и алгоритмы решения задачи о транспорте наносов // Известия ЮФУ. Технические науки. Тематический выпуск: Актуальные проблемы математического моделирования. - Таганрог: Изд-во ТТИ ЮФУ 2009 № 8 (97). С. 71-75.

2. Проценко Е.А., Сухинов A.A. Достаточное условие единственности решения задачи транспорта многокомпонентной взвеси // Естественные и технические науки, -М.: изд-во Спутник+, 2010. № 5 (49). С. 511-518.

3. Проценко Е.А. Двумерная конечно-разностная модель формирования наносов в прибрежной зоне водоема и ее программная реализация // Электронный научно-инновационный журнал Инженерный Вестник Дона, Северо-Кавказский научный центр высшей школы ЮФУ. - Ростов-на-Дону, 2010. №3. [Электронный ресурс]. Систем, требования: http://www.ivdon.ru/magazine/archive/n3y2010/224/

4. Сухинов А.И., Чистяков А.Е., Проценко Е.А. Построение дискретной математической модели транспорта наносов // Известия ЮФУ. Технические науки. Тематический выпуск: Актуальные проблемы математического моделирования. -Таганрог: Изд-во ТТИ ЮФУ, 2011. № 8 (121). С. 32-43.

5. Сухинов А.И., Чистяков А.Е., Проценко Е.А. Двумерная гидродинамиче-

екая модель, учитывающая динамическое перестроение геометрии дна мелководного водоема //Известия ЮФУ. Технические науки, Тематический выпуск: Актуальные проблемы математического моделирования. - Таганрог: Изд-во ТТИ ЮФУ, 2011. № 8 (121). С. 159-166.

В других изданиях:

6. Проценко Е.А., Сухинов А.И. Модифицированный попеременно-треугольный метод решения разностной краевой задачи Дирихле для уравнения эллиптического типа в прямоугольнике с линейной функцией источника // Методы и алгоритмы прикладной математик в технике, медицине и экономике: Материалы VII Международной научно-практической конференции. - Новочеркасск, 2 февраля 2007г.: В 2 ч. / Юж. Рос. гос. техн. ун-т (НПИ). - Новочеркасск: ЮРГТУ, 2007. -Ч. 1.С. 6-7.

7. Проценко Е.А. Исследование устойчивости разностной схемы при моделировании транспорта наносов на основе двумерной модели // Материалы III Международной научно-практической конференции: Молодежь и наука: реальность и будущее. В б томах. - Невинномысск: НИЭУП, 2010. Том V: Естественные и прикладные науки. С. 487-495.

8. Проценко Е.А. Пространственно-одномерная модель транспорта наносов в прибрежной зоне водоема // Актуальные проблемы социально-педагогических знаний: сборник научных статей. - Армавир: РИЦ АГПУ, 2010. С.82-95.

9. Проценко Е.А. Исследование аппроксимации и устойчивости разностной схемы для одномерной модели перемещения наносов в прибрежной зоне // Вестник Таганрогского государственного педагогического института. Физико-математические и естественные науки. Таганрог: Изд. отдел ГОУВПО «ТГГТИ», 2010. № 1.С. 218-224.

Ю.Проценко Е.А. Сухинов А. А. Исследование единственности решения задачи транспорта многокомпонентных взвесей с взаимными трансформациями фракций // В рамках Конгресса по интеллектуальным системам и информационным технологиям «AIS -IT 10». Научное изд. в 4-хтомах. -М.: Физматлит, 2010. Т.2. С. 147-156.

11. Проценко Е.А. Решение задачи о транспорте наносов в прибрежной зоне водоема на основе двумерной модели // Материалы международной научной конференции Славянского государственного педагогического университета. Украина, - Славянск, 2011,№ 11.-Т. 1.С. 109-118.

Личный вклад соискателя в работах, опубликованных в соавторстве: [2,10] - на основе разработанной модели переноса многокомпонентной примеси в водной среде получено достаточное условие единственности решения для задачи транспорта взвесей в случае выполнения условия несжимаемости жидкости. [4] -построение сеточных уравнений для задачи транспорта наносов, разработка алгоритма расчета коэффициентов сеточных уравнений. [5] - на основе гидрофизической модели транспорта наносов в прибрежной зоне построена дискретная модель, проведен анализ результатов численных экспериментов. [6] - применен модифицированный попеременно-треугольный метод для решения разностной краевой задачи Дирихле для уравнения эллиптического типа в прямоугольнике с линейной функцией источника.

Тип.ТТИ ЮФУ Заказ №34$тирЛ>0экз.

Введение.

ГЛАВА I. МОДЕЛИ ТРАНСПОРТА НАНОСОВ И ТРАНСПОРТА МНОГОКОМПОНЕНТНОЙ ВЗВЕСИ.

1.1. Введение в проблему формирования наносов в прибрежной зоне водоема.

1.2. Обзор существующих моделей транспорта наносов.

1.3. Модель транспорта многокомпонентных взвесей.36"

1.4. Исследование единственности решения задачи транспорта многокомпонентных взвесей.

1.5. Дискретная модель транспорта многокомпонентных взвесей.

1.6. Модель транспорта взвешенных материалов в водной среде.

Выводы по главе 1.

ГЛАВА П. ОДНОМЕРНАЯ МОДЕЛЬ ФОРМИРОВАНИЯ НАНОСОВ В ПРИБРЕЖНОЙ ЗОНЕ ВОДОЕМА.

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

2.2. Аналитическое решение волновой задачи.

2.3. Непрерывная математическая модель процессов перемещения наносов в прибрежной зоне.

2.4. Построение дискретной математической модели поверхностных волн.

2.5. Устойчивость разностной схемы.

2.6. Погрешность аппроксимации разностной схемы.

2.7. Построение дискретной модели волновых гидродинамических процессов.

2.8. Численная реализация одномерной математической модели гидродинамики.

Выводы по главе II.

ГЛАВА III. ДИСКРЕТНАЯ МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОЦЕССОВ ПЕРЕМЕЩЕНИЯ НАНОСОВ В ПРИБРЕЖНОЙ ЗОНЕ.

3.1. Дискретизация одномерной математической модели транспорта наносов.

3.2. Численная реализация математической модели.

3.3. Устойчивость разностной схемы.

3.4. Погрешность аппроксимации разностной схемы.

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

3.6. Результаты численных экспериментов.

Выводы по главе III.

ГЛАВА IV. ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ НАНОСОВ НА ОСНОВЕ ДВУМЕРНОЙ МОДЕЛИ.

4.1. Постановка непрерывной задачи.

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

4.3. Дискретизация двумерной математической модели транспорта наносов.

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

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

4.6. Устойчивость дискретной модели транспорта наносов.

4.7. Погрешность аппроксимации конечно-разностной схемы.

4.8. Двумерная математическая модель движения водной среды.

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

4.10. Результаты численных экспериментов.

Выводы по главе IV.

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

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

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

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

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

B.ГТ. Зенковичу, П. А. Каплину, O.K. Леонтьеву, В.И. Лымареву, В.В. Логинову, Г.А. Сафьянову, A.C. Судольскому и др. [31, 33-35, 52, 56, 69, 89].

Интенсивно развивающейся научной отраслью, достигшей определенных успехов, становится математическое моделирование водных экосистем. Проблемам математического моделирования гидродинамических процессов посвящены работы отечественных исследователей В. В. Вайтмана, К.В. Гришанина, И. Г. Кантаржи, A.B. Караушева, С.Ю. Кузнецова, И.О. Леонтьева, Г.И. Марчука, А.Л. Чикина и мн.др [19, 20, 26, 40, 43, 51, 62, 63, 103]. Представлены они и в ряде зарубежных публикаций [117, 121, 125, 128].

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

C.М. Анцыферова, Р.Д. Косьяна, Н.В. Пыхова [4, 32, 43, 129, 130]. В данных исследованиях выявлены основные особенности распределения осредненных значений концентрации взвешенных частиц.

Проведение натурных экспериментов в условиях реального моря сопряжено со значительными трудностями, поэтому количественная оценка распределения взвешенных наносов в прибрежной зоне нередко производится путем математического моделирования с использованием параметров, полученных в лабораторных условиях. Достаточно полный обзор подобных моделей приведен в работе С.М. Анцыферова, В.К. Дебольского, Т.М. Акивис [4].

При исследовании гидродинамических процессов прибрежной зоны водоемов были выявлены некоторые общие черты руслового и волнового потоков при воздействии движущейся воды на донные наносы. Поэтому в ряде случаев, возможен перенос научных достижений гидродинамики рек на область береговых процессов в водоемах и наоборот. Значительный вклад в теорию формирования руслового рельефа и транспорта наносов внесён Н.И. Алексеевским, М.А. Великановым, Н.Б. Барышниковым, КВ. Гришаниным, A.B. Караушевым, К.И. Россинским, Р. С. Чаловым и др. [1, 9, 26,39-40, 74, 100,101].

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

Процесс продуцирования турбулентной кинетической энергии, генерируемой деформированными и обрушающимися волнами, и ее влияние на взвешивание наносов, находятся в центре внимания отечественных и зарубежных исследований [5, 6, 13, 48, 49, 119-121, 124, 131, 138, 140]. Теоретическое обоснование теории турбулентности всесторонне изложено в работах Р.В. Озмидова, И. О. Хинце, A.M. Яглома и др. [68, 71, 99]. Численные методы исследования диффузии примеси в турбулентной среде развиты группой российских ученых под руководством Г.И. Марчука [65-67].

Широко известны исследования по изучению океанической турбулентности [67], моделированию циркуляции в океане [63]. Структурный анализ одного из важнейших направлений в исследовании турбулентных течений, связанного с конструированием моделей турбулентности, дан в работе И.А. Белова и С.А. Исаева [66], где рассмотрены методы теоретического и численного исследования турбулентных течений. Ретроспективный взгляд на алгебраический и статистический подходы моделирования турбулентного движения представлен в работах Ю.В.Лапина. Преимущественное внимание уделено вопросам применения моделей турбулентности в рамках сложившихся вычислительных технологий [48,49].

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

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

Проведенный анализ исследований по проблеме исследования показал, что

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

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

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

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

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

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

Объектом исследования являются процессы перемещения наносов на мелководье под воздействием поверхностных гравитационных волн.

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

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

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

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

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

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

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

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

Научная новизна диссертационного исследования.

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

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

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

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

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

Методы исследования. Основой теоретического исследования служат:

- основы комплексного берегового учения (исследования гидродинамических процессов береговой зоны [21, 30-34, 51, 89, 101], теория донного влечения наносов [4, 42, 46, 52, 78, 79]);

- методология математического моделирования и вычислительного эксперимента [58, 61, 82, 86];

- математические методы численного и аналитического анализа моделей (уравнения математической физики [18, 44, 45], численные методы [14, 60, 62, 85, 90, 97], теория разностных схем [22, 24, 27, 81, 84, 92]);

- методы теоретического и численного исследования турбулентности [5, 13, 70, 72, 91, 100].

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

Программная реализация построенных численных алгоритмов базируется на языке «С++». Визуализация и анализ решений, компьютерные эксперименты с моделью проводились в среде разработки МАТНСАО.

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

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

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

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

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

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

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

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

В рамках сформулированной в работе проблемы на защиту выносятся следующие результаты и положения:

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

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

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

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

Апробация работы. Результаты, полученные в рамках диссертационной работы, докладывались и обсуждались на конференциях и научных семинарах: VII Международная научно-практическая конференция: Методы и алгоритмы прикладной математики в технике, медицине и экономике (Новочеркасск, 2 февраля 2007); Пятьдесят третья научно-теоретическая конференция ГОУВПО «ТГПИ». Секция математики. (Таганрог, 22 января 2010); Пятьдесят четвертая научно-теоретическая конференция ГОУВПО «ТГПИ». Секция математики. (Таганрог, 21 января 2011); Международная научно-практическая конференция: Молодежь и наука: реальность и будущее. (Невинномысск, НИЭУП, 3 марта 2010); Научный семинар кафедры ВМ ТТИ ЮФУ (Таганрог, 2011).

Основное содержание диссертации изложено в следующих работах:

По теме исследования опубликовано 11 печатных работ, из них 5 статей в отечественных реферируемых журналах, входящих в список изданий, рекомендованный ВАК.

В изданиях из перечня ВАК:

1. Проценко Е.А. Модель и алгоритмы решения задачи о транспорте наносов // Известия ЮФУ. Технические науки. Тематический выпуск: Актуальные проблемы математического моделирования. - Таганрог: Изд-во ТТИ ЮФУ, 2009. № 8 (97). С. 71-75.

2. Проценко Е.А., Сухинов A.A. Достаточное условие единственности решения задачи транспорта многокомпонентной взвеси // Естественные и технические науки, -М.: изд-во Спутник+, 2010. № 5 (49). С. 511-518.

3. Проценко Е.А. Двумерная конечно-разностная модель формирования наносов в прибрежной зоне водоема и ее программная реализация // Электронный научно-инновационный журнал Инженерный Вестник Дона, Северо-Кавказский научный центр высшей школы ЮФУ. — Ростов-на-Дону, 2010. №3. [Электронный ресурс]. Систем, требования: http://www.ivdon.ru/magazine/archive/n3y2010/224/

4. Сухинов А.И., Чистяков А.Е., Проценко Е.А. Построение дискретной математической модели транспорта наносов // Известия ЮФУ. Технические науки. Тематический выпуск: Актуальные проблемы математического моделирования. - Таганрог: Изд-во ТТИ ЮФУ, 2011. № 8 (121). С. 32-43.

5. Сухинов А.И., Чистяков А.Е., Проценко Е.А. Двумерная гидродинамическая модель, учитывающая динамическое перестроение геометрии дна мелководного водоема // Известия ЮФУ. Технические науки. Тематический выпуск: Актуальные проблемы математического моделирования. - Таганрог: Изд-во ТТИ ЮФУ, 2011. № 8 (121). С. 159-166.

В других изданиях:

6. Проценко Е.А., Сухинов А.И. Модифицированный попеременно-треугольный метод решения разностной краевой задачи Дирихле для уравнения эллиптического типа в прямоугольнике с линейной функцией источника // Методы и алгоритмы прикладной математик в технике, медицине и экономике: Материалы УП Международной научно-практической конференции. - Новочеркасск, 2 февраля 2007г.: В 2 ч. / Юж. Рос. гос. техн. ун-т (НПИ). - Новочеркасск: ЮРГТУ, 2007. -Ч. 1. С. 6-7.

7. Проценко Е.А. Исследование устойчивости разностной схемы при1 моделировании транспорта наносов на основе двумерной модели // Материалы III Международной научно-практической конференции: Молодежь и наука: реальность и будущее. В 6 томах. - Невинномысск: НИЭУП, 2010. Том V: Естественные и прикладные науки. С. 487-495.

8. Проценко Е.А. Пространственно-одномерная модель транспорта наносов в прибрежной зоне водоема // Актуальные проблемы социально-педагогических знаний: сборник научных статей. - Армавир: РИЦ АГТТУ, 2010. С.82-95.

9. Проценко Е.А. Исследование аппроксимации и устойчивости разностной схемы для одномерной модели перемещения наносов в прибрежной зоне // Вестник Таганрогского государственного педагогического института. Физико-математические и естественные науки. Таганрог: Изд. отдел ГОУВПО «ТГПИ», 2010. № 1. С. 218-224.

10.Проценко Е.А. Сухинов А. А. Исследование единственности решения задачи транспорта многокомпонентных взвесей с взаимными трансформациями фракций // В рамках Конгресса по интеллектуальным системам и информационным технологиям «AIS -IT110». Научное изд. в 4-х томах. - М.: Физматлит, 2010. Т.2. С. 147-156.

11. Проценко Е.А. Решение задачи о транспорте наносов в прибрежной зоне водоема на основе двумерной модели // Материалы международной научной конференции Славянского государственного педагогического университета. Украина, - Славянск, 2011, № 11. - Т. 1. С. 109-118.

Личный вклад соискателя в работах, опубликованных в соавторстве: [2, 10] — на основе разработанной модели переноса многокомпонентной примеси в водной среде получено достаточное условие единственности решения для задачи транспорта взвесей в случае выполнения условия несжимаемости жидкости. [4] — построение сеточных уравнений для задачи транспорта наносов, разработка алгоритма расчета коэффициентов сеточных уравнений. [5] — на основе гидрофизической модели транспорта наносов в прибрежной зоне построена дискретная модель, проведен анализ результатов численных экспериментов. [6] - применен модифицированный попеременно-треугольный метод для решения разностной краевой задачи Дирихле для уравнения эллиптического типа в прямоугольнике с линейной функцией источника.

Структура и объем работы. Работа состоит из введения, 4 глав, заключения, списка литературы, содержащего 151 наименование. Работа содержит 33 рисунка, 2 таблиц. Полный объем диссертации составляет 196 страниц.