CC BY
52
Раздел II. Математическое моделирование физических
процессов
УДК 532.5.031
А.И. Сухинов, А.Е. Чистяков, Е.Е. Дегтярева
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ТРАНСПОРТА ДОННЫХ
ОТЛОЖЕНИЙ С УЧЕТОМ ГИДРОДИНАМИЧЕСКИХ ПРОЦЕССОВ
Рассматривается трехмерная математическая модель транспорта наносов. Предложенная модель учитывает следующие факторы: транспорт взвешенных частиц, переход донного материала во взвешенное состояние и осаждение. Для описания транспорта взвешенных частиц применяется уравнение диффузии-конвекции-реакции. В качестве вертикальной переменной используется безразмерная переменная а€[0,1]. При численной реализации модели был использован метод расщепления по пространственным координатам и . -грамм, предназначенный для расчетов полей скорости движения водной среды и транс.
Транспорт веществ; гидродинамика; уравнение диффузии-конвекции.
A.I. Sukhinov, A.E. Chistyakov, E.E. Degtyreva
MATHEMATICAL MODELLING BOTTOM SEDIMENTS TRANSPORT
IN VIEW OF HYDRODYNAMIC PROCESSES
This paper proposed the three-dimensional mathematical model of sediment transport. The model takes into account the following factors: transport of suspended particles, suspension and deposition. To describe the transport of suspended particles diffusion the equation of convection — reaction is applied. As the vertical variable is used dimensionless variable a € [0,1]. In the numerical implementation of the model we used the method of splitting on the spatial coordinates and the physical processes. The developed algorithms was built a program designed to calculate velocity fields and the aquatic environment of the impurity transport.
Sediment transport; hydrodynamics; diffusion-convection equation.
Введение. Перемещение взвесей и наносов под воздействием течений и волн влияет на состояние берегов и подходных каналов, поэтому научные знания о гидродинамических процессах очень важны в условиях расширения хозяйственной деятельности. Отсутствие учета расхода наносов может привести к авариям, переформированию донной поверхности и уменьшению пропускной способности подходных каналов [1]. Значительная часть наносов перемещается во взвешенном состоянии, поэтому изучение транспорта наносов невозможно без учета транспорта взвесей. Гидролито динамические процессы береговой зоны изучены в работах С.М. Анцыферова, Р.Д. Косьяна, Н.В. Пыхова [2-5] на основе многочисленных натурных экспериментов. В данных исследованиях выявлены основные особенности распределения осредненных значений концентрации взвешенных частиц.
Основным фактором взвешивания, перераспределения и транспорта донного материала является комбинированное воздействие волн и течений. Математическое моделирование помогает решать задачи прогнозирования образования наносов в прибрежной зоне моря. С помощью таких моделей можно оценить подвиж-
ность наносов; моделировать переформирование дна под действием волн и течений; моделировать перенос наносов в окрестности сооружений, находящихся на берегу или на некотором расстоянии от него [5].
Когда частица находится во взвешенном состоянии в воде, ее веса недостаточно для осаждения. Поэтому осаждаться могут только хлопья, образованные из частиц или хлопьев и частиц (рис. 1).
Рис. 1. Графическое представление механизмов горизонтального и вертикального
движения взвешенных частиц
Описание модели. Для описания транспорта вз вешенных частиц введем систему координат. В модели будем использовать декартовую систему координат в горизонтальной плоскости и а - координатную систему в вертикальном направле-
£ - Л
нии (рис. 2): а =-----, ха = х, уа = у , їа = Ї, здесь х, у, £, а = а = 1 на сво-
Н + 7
бодной поверхности, а = Ь = —1 на дне; Н = Н + 7 - общая глубина до свободной поверхности, Н=Н(х,у) - глубина водного объекта, 7 = 7(ґ, х, у) - возвыше-
( ).
Для описания транспорта взвешенных частиц будем предполагать, что расчетная область ^ вписана в параллелепипед.
П
і о=1=а
о=-1=Ь
Рис. 2. о - координатная система
Для описания транспорта взвешенных частиц использовано уравнение диф-
фузии-конвекции-реакции, которое может быть записано в следующем виде:
дС дС дС а - Ь, * дС
----ъи----ЪV-----\----(ю-ю )-----=
дг дх ду Н v да
д (п дС} д ( „ дС' ' '' г (1)
л (а -ЬЛ2 д к ~]
) Ч н ) да да
дх^ Н дх) ду^ Н ду где С - концентрация осадка [г/л или кг/лг3]; V = {и, V, ы] - составляющие поля вектора скорости [м/с]; - гидравлическая крупность или скорость осаждения по
о-координате в вертикальном направлении [м/с]; Н - глубина [м]; Бк, - гори-
зонтальный и вертикальный коэффициенты диффузии [л*2/с]; х, у - координаты в горизонтальном направлении; а- координата в вертикальном направлении; I - временная переменная [с].
Уравнение (1) рассматривается при граничных и начальных условиях. Начальная скорость устанавливается равной нулю, начальный уровень подъема воды, температуры, солености равен средним значениям подъема воды, температуры, солености исследуемого объекта;
- , :
Н
■ вблизи поверхности дна: а - Ь
Н
в, ^ + о£к = о,
да
’Л/-'
В -----------= Е - В + ш Ск,
V '■Ч 5 к ’
да
■ на открытых границах:
дС „ дС Л -----+ К----------= 0,
ді дп
(2)
(3)
(4)
где Е - поток эрозии [кг/л* /с]; Б - интенсивность осаждения осадка [кг/л* /с]; Ск - массовая концентрация взвеси [масса/объем]; п - единичный вектор нормали к открытой границе Гореп;
В =
0
т > т
Ь кр
(О С,
ґ \ Ґ \
1 тЬ • Е = < т
Т <т ’ Ь кр м -1
К Тт у тЕ
0
т<т
крЕ
т>т
Е
где т - касательное напряжение сдвига у дна [#/лг2 ]; - критическое касательное
напряжение у дна, которое оценивается на основе лабораторных испытаний как величина, со значениями от 0,05 до 0,15 [#/л*2]; ткрЕ - критическое касательное напряжение
эрозии, определенное как те = 0,015[рЬ -1000)°'73; рь - объемная плотность вещества в диапазоне плотности 1000-2000 ^/сг/л/3]; М - экспериментально-
наблюдаемая, постоянная эрозии, со значениями между 0,0002 и 0,002 Цк -Н-1с-1 ] . Коэффициент вертикального турбулентного обмена, неоднородный по глуби, :
у = С2 А21
( диЛ 2
+
го ^ і
(5)
где и, V - осредненные по времени пульсации горизонтальных компонент скорости, А - характерный масштаб сетки, С - безразмерная эмпирическая константа, значение которой обычно определяется на основе расчета процесса затухания однородной изотропной турбулентности.
Одним из преимуществ данной модели является усовершенствованная параметризация коэффициента турбулентного обмена по вертикали, подобранная исходя из сравнений данных экспериментальных измерений. Следует также отме-,
счет сгонно-нагонных явлений осуществлялся учет частичной «заполненности» ячеек [13-16].
Построение дискретной модели. Для дискретизации поставленной задачи введем равномерную прямоугольную сетку: са = юх хау хаа ха,
а : \х. = .к ,I = 0,N -1,к (И -1) = I ], а :{у. = ¡к , . = 0,N -1,к ( -1) = I ],
х [ г х’ ’ х ’х\х / ху у . J у* J ’у ’у\у / У_|’
а : |а, = кк ,к = 0,N -1,к (И -1) = I ], а,:\г = пт,п = 0И~,гИ, = I,],
а у к а ’ ’а ’а \ а / а \ , [ п ’ ’ , ’ , , у
где т - шаг по времени; кх, ку ка - шаги по пространству; - количество
временных слоев; Т - верхняя граница по времени; Их, Иу N а - количество узлов по пространству; 1х, 1у 1а - размеры параллелепипеда по координатным направлениям. Для аппроксимации будем использовать схемы расщепления. Согласно этой схеме на первом этапе осуществляется решение системы трехдиагональных алгебраических уравнений методом прогонки в одном из направлений, в результате чего находятся значения искомой функции на промежуточном (п + 1/2) -м временном слое. На втором этапе находится искомое решение на верхнем (п+1)-м вре-. , решение не зависит от величины временного шага.
Результаты численного эксперимента. Из результатов проведенного численного эксперимента (рис. 3) видна картина поля течения водной среды при об-
( ).
Рис. 3. Поле распределения давления от источника в момент времени г без учета концентрации и турбулентного обмена
Ниже приведены результаты численного моделирования транспорта веществ при обтекании водной средой волнореза. В данном эксперименте источник веществ задается точечной функцией и поле загрязняющих веществ задавалось ну.
Рис. 4. Поле конценртации от источника в момент времени I
Из рис. 4 видно, как изменится поле концентрации веществ через заданное время t. Функция концентрации веществ показана палитрой.
Выводы. В работе построена математическая модель транспорта донного , : взвесей со дна, конвективно-диффузионный перенос, а также осаждение донного материала. Для описания движения водной среды использованы модели гидродинамики, описываемые уравнением Навье^токса [8-11]. Следует также отметить, что при построении дискретной модели использовались разностные схемы [6], а при решении сеточных уравнений - эффективный метод решений сточных уравнений (адаптивный МПТМ) [12]. К достоинствам разработанной модели также можно отнести высокую точность, которая достигается за счет учета степени «заполненности» контрольных ячеек.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Сух иное А.И. Прецизионные модели гидродинамики и опыт применения в предсказании и реконструкции чрезвычайных ситуаций в Азовском море // Известия ТРТУ. — 2006.
— № 3 (58). — С. 228-235.
2. Анцыферов С.М., Дебольский В.К., Акивис Т.М. О влиянии состояния дна на формирование потока несвязных наносов // Водные ресурсы. — М.: Наука, 2004. — Т. 31, № 6.
— С. 675-682.
3. / . . . , . . ,
Н.В. Пыхова. — М.: Научный мир, 2003. — 320 с.
4. . ., . . -
// . .
— 2012. — № 2 (127). — С 112-118.
5. . ., . . -
//
программирование: Новые вычислительные технологии. — 2012. — Т. 13. — С. 290-297.
6. Самарский АЛ. Теория разностных схем. — М.: Наука, 1989.
7. Заковоротнова Е.Е. Выбор модели транспорта наносов в мелководном водоеме // Прикладные технологии гидроакустики и гидрофизики. - СПб.: Российская академия наук, 2010. - С. 284-286.
8. Алексеенко ЕМ., Сидоренко Б.В., Колгунова ОМ., Чистяков А.Е. Сравнительный анализ классических и неклассичнских моделей гидродинамики водоемов с турбулентным обменом // Известия ЮФУ. Технические науки. - 2009. - № 8 (97). - С. 6-18.
9. Лапин ДМ., Черчаго АЛ., Чистяков А.Е. Совместные экспедиционные исследования гидрофизических параметров Азовского моря на многоцелевой яхте «Буревестник» и НИС т/х «Платов» // Известия ЮФУ. Технические науки. - 2009. - № 8 (97). - С. 82-89.
10. . .
// . . -2009. - 8 (97).
- С. 75-82.
11. . ., . ., . . -
//
ЮФУ. Технические науки. - 2008. - № 10 (87). - С. 171-175.
12. . . -
// . . - 2010. - 6
(107). - С. 237-249.
13. Чистяков А.Е. Об аппроксимации граничных условий трехмерной модели движения
// . . - 2010. - 6 (107). - . 66-77.
14. . ., . ., . .
// .
Технические науки. - 2011. - № 8 (121). - С. 22-32.
15. . ., . ., . . -
тической модели транспорта наносов // Известия ЮФУ. Технические науки. - 2011.
- № 8 (121). - С. 32-44.
16. Сухинов А.И., Чистяков А.Е., Процент ЕЛ. Двумерная гидродинамическая модель,
// -
вестия ЮФУ. Технические науки. - 2011. - № 8 (121). - С. 159-167.
Статью рекомендовал к опубликованию д.ф.-м.н. Г.В. Куповых.
Сухинов Александр Иванович - Технологический институт федерального государственного автономного образовательного учреждения высшего профессионального образования «Южный федеральный университет» в г. Таганроге; e-mail: sukhinov@gmail.ru; 347928, г. Таганрог, пер. Некрасовский, 44, ГСП 17А; тел.: 88634310599; руководитель ТТИ ЮФУ.
Чистяков Александр Евгеньевич - e-mail: cheese_05@mail.ru; тел.: 88634371606; кафедра высшей математики; к.ф.-м.н.; доцент.
Дегтярева Екатерина Евгеньевна - e-mail: katerina.degtyareva@gmail.com; тел.: 89045069696; кафедра высшей математики; аспирантка.
Sukhinov Alexander Ivanovich - Taganrog Institute of Technology - Federal State-Owned Autonomy Educational Establishment of Higher Vocational Education “Southern Federal University”; e-mail: sukhinov@gmail.ru; GSP 17A, 44, Nekrasovsky, Taganrog, 347928, Russia; phone: 88634310599; head the TIT SFedU.
Chistyakov Alexander Evgenjevich - e-mail: cheese_05@mail.ru; phone: +78634371606; the department of higher mathematics; cand. of phis.-math. sc. associate professor.
Degtyareva Ekaterina Evgenevna - e-mail: katerina.degtyareva@gmail.com; phone: 89045069696; the department of higher mathematics; postgraduate student.
