CC BY
59
4. Diniz P.S.R. Adaptive filtering algorithms and practical implementation. Third edition. -Springer Science + Business Media, 2008. - 627 p.
5. Волков Ю.И., Джиган M.B., Джиган O.B. Виртуальные инструменты для исследования переходных процессов в линейных электрических цепях // Вопросы радиоэлектроники. Серия «Общетехническая». - 2011. - №1. - C. 167-177.
Статью рекомендовал к опубликованию д.т.н., доцент А.В. Гуреев
Джиган Мария Виктгоровна - Национальный исследовательский университет «МИЭТ»; e-mail: [email protected]; 124498, Москва, Зеленоград, проезд 4806, д. 5; тел.: +74997109541; сту-.
Джиган Ольга Виктгоровна - e-mail: [email protected]; студентка.
Djigan Maria Victorovna - National Research University of Electronic Technology; e-mail: [email protected]; Bld. 5, Pas. 4806, Zelenograd, Moscow, 124498, Russia; phone: +74997109541; student.
Djigan Olga Victorovna - e-mail: [email protected]; student.
УДК 519.63: 532.5
E.E. Дегтярева, A.E. Чистяков
МОДЕЛИРОВАНИЕ ТРАНСПОРТА НАНОСОВ ПО ДАННЫМ
ЭКСПЕРИМЕНТАЛЬНЫХ ИССЛЕДОВАНИЙ В АЗОВСКОМ МОРЕ
Представлены результаты экспериментальных исследований в Таганрогском заливе, проведенных в октябре 2011 года, которые легли в основу математической модели транспорта наносов в прибрежной зоне водоемов. Предложенная модель учитывает следующие факторы: транспорт наносов, находящихся в придонном слое; переход донного материала во взвешенное состояние и осаждение; транспорт взвешенных частиц. Описание гидродинамических процессов выполнено на основе системы уравнений Навье-Стокса. Для описания коэффициента турбулентного обмена использовалась подсеточная модель Смогорин-
.
зондов ADCP и SBE19 Plus.
Акустический Дотеровский профилограф; гидрофизика моря; модель транспорта .
E.E. Degtyareva, A.E. Chistyakov MODELING SEDIMENT TRANSPORT BASED ON EXPERIMENTAL STUDIES IN AZOV SEA
This paper presents the results of experimental studies in the Taganrog Bay in October 2011. By results of experimental studies formed the basis of a mathematical model of sediment transport in coastal waters. The proposed model considers the following factors: sediment transport in the near-bed layer; suspension, diffusion, settling; suspended sediment transport. Description of the hydrodynamic processes performed on the basis of the Navier-Stokes equations. A subgrid model Smogorinskogo was used to describe the turbulent exchange coefficient. Experiments were performed using two modern devices ADCP and SBE 19 Plus.
Acoustic Doppler Current Profiler; hydro physics sea; model of sediment transport.
Экологическая система Таганрогского залива и Азовского моря является уникальной. Состояние прибрежного рельефа дна во многом определяется движением наносов под действием волн и течений. Осуществленные без предварительного математического моделирования проекты по сооружению гидротехнических
, -су и заилению, также в результате непродуманного сооружения и углубления ка, .
1. Зоны аноксии и гипоксии в Азовском море. В ходе экспедиционных исследований в 2001 г. была обнаружена зона аноксии и наблюдалась массовая гибель ихтиофауны в восточной части Азовского моря. Площадь обнаруженной об-
1 000 . .
-| | —-1 - I 1
35 36 37 38 39
Рис. 1. Измеренные значения концентрации кислорода (мкМ). Контурная линия
показывает зону гипоксии
Механизм возникновения зоны был гипоксии интерпретирован руководителем экспедиционных исследований АЛ. Сухиновым [1]. Было высказано предпо-, -вых (замкнутых) образований в структуре потоков (течений), названных впослед-S- , -, .
2. Математическое моделирование гидрофизических процессов в Азов. -коточных моделей гидрофизических процессов в мелководных водоемах [2]. Данные модели описывают движение водной среды с учетом следующих факторов: ветровые течения и трение о дно, стоки рек, испарение, сила Кориолиса, турбу-, - , .
Исходными уравнениями гидродинамики мелководных водоемов являются [2]:
- уравнение движения (Навье-Стокса):
u' + uux + vuy + wuy = —— p'x + (цu'x )X + (Mu'y )'y + (vu' )' + 2Q.(v sine- w cose),
v' + uv'x + vv'y + wv'y = -—— p'x + (Mv'x )X + (Mv'y )'y + (vv[ )' - 2Ш sin в,
w' + uw'x + vw'y + ww'y =—- p' + (p.w'x) x + (juw'y )'y + (vw')' + 2Q.u cose; (1)
- -:
P + (—u) x + (pv)' + (pu)' = 0, (2)
где u ={u, v, w} - компоненты вектора скорости, p - превышение давления над гидростатическим давлением невозмущенной жидкости, р - плотность, Q - угло-
вая скорость вращения земли, 0 - угол между вектором угловой скорости и вертикалью, - горизонтальная и вертикальная составляющая коэффициента турбу-
лентного обмена.
Система уравнений (1)-(2) рассматривается при соответственных граничных и начальных условиях. Составляющие тангенциального напряжения для дна, с учетом введенных обозначений, могут быть записаны следующим образом:
=рср (И)и\й\, гу =рср (\й\)у\й\. (3)
Коэффициент вертикального турбулентного обмена, неоднородный по глубине:
v = C 2 А2
Ґ ---\2
гдил
dz
Ґ —\2 fdV^
dz
(4)
где и, V - осредненные по времени пульсации горизонтальных компонент скорости, А - характерный масштаб сетки, С - безразмерная эмпирическая константа.
Отличительной особенностью разработанных математических моделей явля-,
трех уравнений движения, а не на основе гидростатического приближения. Следует также отметить, что при описании рельефа дна и изменения геометрии водной поверхности за счет сгонно-нагонных явлений осуществлялся учет частичной «заполненности» ячеек [3]. Для решения сеточных уравнений используется адаптивный модифицированный попеременно-треугольный итерационный метод вариационного типа [4].
Рис. 2. Результаты численного моделирования
( . 2) -
чие вихревой структуры течения в восточной части Азовского моря, в данном районе вода богата органическими примесями, источниками которых являются реки Дон и Кубань.
3. Экспедиционные исследования 2011 г. С целью обновления баз данных многолетних наблюдений состояния водной среды, а также калибровки трехмерных математических моделей транспорта взвешенных частиц, сотрудниками ТТИ ЮФУ совместно с сотрудниками геолого-географического факультета ЮФУ в 2011 . - . -ту НИС т/х «Платов» находилось современное оборудование: акустический профилограф ADSP WHS600, гидрофизический зонд SBE 19plus с датчиком мутности Seapoint и датчиком растворенного кислорода, пробоотборная система, GPS.
Принцип работы ЛБСР основан на эффекте Доплера, т.е. на определении смещения частоты акустического эхо-сигнала, характеризующего проекцию скорости рассеивателей в измерительном объеме вдоль направления луча. Прибор производит измерение проекций скорости течений (рис. 3) на четыре пространственные оси для каждого элемента глубины, что в сочетании с показаниями датчиков угловой ориентации позволяет производить расчеты проекций измеренной скорости на оси географической системы координат.
Рис. 3. Профили компонентов вектора скорости движения водной среды
ЛБСР обеспечивают одновременное измерение профилей скорости течения и интенсивности обратно рассеянного акустического сигнала (рис. 4). Последний параметр зависит от количества и акустических свойств обратных рассеивателей в ,
взвеси. Совместная обработка данных по течениям и интенсивности эхо-сигнала открывает широкие возможности для исследования пространственно-временного распределения потоков взвеси в мелководных акваториях.
Рис. 4. Интенсивность отраженного акустического эхо-сигнала
В ходе экспериментальных исследований при помощи зонда ББЕ 19р1ш сделаны замеры: профиля мутности (рис. 5), растворенного кислорода (рис. 6), солености (рис. 7), температуры (рис. 8).
Рис. 5 Функция Рис. 6. Функция Рис. 7. Функция Рис. 8. Функция
распределения распределения распределения распределения
мутности растворенного солености температуры
кислорода
, , -тическая модель транспорта наносов в Азовском море.
4. Математическая модель транспорта наносов. Основными факторами , -
бинированное воздействие волн и течений. С помощью математических моделей можно оценить подвижность наносов; моделировать переформирование дна под действием волн и течений; моделировать перенос наносов в окрестности сооружений [6]. Для описания транспорта наносов будем использовать математическую модель, учитывающую следующие факторы: транспорт наносов, находящихся в ; -ние; транспорт взвешенных частиц.
Уравнения процесса транспорта наносов без перехода во взвешенное состояние [6] записываются в следующем виде:
дН —
(1 -є')—- + йіу( кть) = йїу дг
Ьс
к = -
Ашй
((Р -Ро)§й)
Р
Ьс
%тадН
8іи^0
Р-1
к
^тайН
- ,
(5)
V
Ьс
8ІП^0
^тайН
-ти.
где Н - глубина дна, отсчитываемая от невозмущенной поверхности водоема; в
пористость грунта; х, у - горизонтальные декартовы координаты; - касательное
напряжение на дне; ?Ъс - критическое значение касательного напряжения, при котором начинается перемещение наносов, А и р - безразмерные постоянные, Ш -частота волны, й - характеристика осадков, р - плотность, ¥с - поток взвеси.
Модель (5) описывает транспорт наносов без перехода их во взвешенное состояние. В [7] рассматривается математическая модель транспорта наносов с переходом их во взвешенное состояние
-+ ид- + у^ + (w-н’о)^ = ^| | + ^|
дс дс дс „ „ Эс д ( дс^ д ( дс | д ( дс ^
■ ■ ■ +—I Ч-э^\ +1(х,у,г,г), (6)
дt дх ду дг дх \ дх) ду ^ ду ^
где с(х,у,г,1) - концентрация частиц в точке (х,у,г) в момент времени г;
и = {и, V, м) - компоненты вектора скорости; м0 - скорость осаждения частиц
или их гидравлическая крупность; /л,ц - коэффициенты диффузии; /(х,у,&) -функция, описывающая интенсивность и распределение источников. Взмучивание и осаждение частиц описывают граничные условия.
Вертикальный поток взвешенных частиц определяется разностью между процессом осаждения Б и процессом взмучивания частиц Е.
дс
¥с = В - Е, В = м0с, Е = --Ц —, (7)
дг
где ¥с- поток взвеси, ц - вертикальный коэффициент диффузии.
. -
ния донных осадков важное место занимает построение общей физической карти, . моделирования является получение достаточного количества экспериментальных данных для калибровки и верификации модели.
В результате экспериментальных исследований была построена математическая модель транспорта наносов и проведена первичная верификация данных. Математическая модель транспорта наносов учитывает следующие факторы: движение водной среды, сложную геометрию дна и береговой линии, транспорт наносов, находящихся в придонном слое, переход донного материала во взвешенное состояние и осаждение, транспорт взвешенных частиц. Описаны процессы срыва, переноса и осаждения частиц.
Представленная модель может быть использована для прогнозирования переформирования донной поверхности Азовского моря.
В заключение мы хотим выразить глубокую признательность нашему постоянному наставнику профессору АЛ. Сухинову за многочисленные плодотворные обсуждения и внимание к работе.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Сух иное Л.И. Прецизионные модели гидродинамики и опыт применения в предсказании и реконструкции чрезвычайных ситуаций в Азовском море // Известия ТРТУ. - 2006.
- № 3 (58). - С. 228-235.
2. . ., . ., . .
// -
тическое моделирование. - 2011. - Т. 23, № 3. - С. 3-21.
3. Самарский АЛ. Теория разностных схем. - М.: Наука, 1989.
4. Чистяков А.Е. Теоретические оценки ускорения и эффективности параллельной реализации ПТМ скорейшего спуска // Известия ЮФУ. Технические науки. - 2010. - № 6 (107).
- С. 237-249.
5. . . . -
. - .: ,
2010. - . 284-286.
6. . ., . ., . . ,
// -
вестия ЮФУ. Технические науки. - 2011. - № 8 (121). - С. 159-167.
7. . .
// .
науки. - 2010. - № 6 (107). - С. 258-259.
Статью рекомендовал к опубликованию д.ф.-м.н., профессор АЛ. Жорник.
Дегтярева Екатерина Евгеньевна - Технологический институт федерального государственного автономного образовательного учреждения высшего профессионального образования «Южный федеральный университет» в г. Таганроге; e-mail: [email protected]; 347928, г. Таганрог, пер. Некрасовский, 44; тел.: +790450б9б9б; кафедра высшей математики; .
Чистяков Александр Евгеньевич - e-mail: [email protected]; тел.: +78б34371б0б; кафедра высшей математики; доцент.
Degtyareva Ekaterina Evgenievna - Taganrog Institute of Technology - Federal State-Owned Autonomy Educational Establishment of Higher Vocational Education “Southern Federal University”; e-mail: [email protected]; 44, Nekrasovsky, Taganrog, 347928, Russia; phone: +790450б9б9б; the department of higher mathematics; post-graduate student.
Chistyakov Alexander Evgenjevich - e-mail: [email protected]; phone: +78б34371б0б; the department of higher mathematics; associate professor.
УДК 532.5.031
Кандалфт Хекмат
ТРЕХМЕРНАЯ МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ДВИЖЕНИЯ ВОДНОЙ СРЕДЫ ПРИ НАЛИЧИИ НА ПОВЕРХНОСТИ ЛЕДЯНОЙ ПЛАСТИНЫ
Работа посвящена развитию математических подходов к изучению гидродинамических процессов в однородной жидкости, поверхностью активными напряжениями и начальными возмущающими факторами. Разработана трехмерная математическая модель движения водной среды, при наличии ледовой пластины на поверхности водоема. Основными уравнениями математической модели являются: система уравнений Навье-Стокса и уравнение неразрывности. Для построения численного алгоритма применен метод расщепления по физическим процессам. Для данной задачи получено численное решение, позволяющее определять поведение жидкости при данных условиях. На основе интегро-интерполяционного метода выполнена аппроксимация трехмерной математической модели гидродинамики. Был разработан комплекс программ, предназначенный для построения трехмерных полей скоростей движения водной среды.
Ледяная пластина; ячейка; гидродинамика; уравнений движения жидкости; уравнений неразрывности; заполненность.
Hekmat Kandalft
TREE-DIMENSIONAL MATHEMATICAL MODEL OF A MOTION OF THE AQUATIC ENVIRONMENT, THE SURFACE OF WHICH IS COVERED BY
ICE PLATE
This article is dedicated to the development of mathematical approaches of the study of hydrodynamic processes in an homogeneous _ fluid, the surface of the active stresses and initial disturbances. A mathematical model of the motion of w in the reservoir is partially covered by ice plate. For this task, the numerical solution determines the behavior of the fluid under these conditions. It creates a three-dimensional discrete fnite-volume model describing the hydrodynamic processes, taking into account the ice cover on the surface. Designed program was developed to build three-dimensional velocity fields of motion of the aquatic environment in the event of mathematical modeling of the reservoir area covered by an ice layer. The numerical model is threedimensional; the decision is based on the Navier-Stokes approximations. To construct a numerical algorithm, a method of splitting by physical processes was used.
The ice plate; hydrodynamics; the equations of fluid motion; the equations of continuity; the occupancy rate.
