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

В отличие от классических подходов, основанных на уравнениях Навье - Стокса, которые, как правило, требуют итерационных процедур.

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. S. Chen and G. Doolen, Ann. Rev. Fluid Mech. 8, 2527 (1998).

2. S. Succi, The Lattice Boltzmann Equation for Fluid Dynamics and Beyond (Clarendon Press, Oxford, 2001).

3. Bhatnagar P.L., Gross E.P. and KrookM., A model for collision processes in gases. I. Small amplitude processes in charged and neutral one-component system, Phys. Rev., 94, 1954.

- P. 511-525.

4. D. Kandhai, A. Koponen, A. Hoekstra, M. Kataja, J. Timonen, andP.M.A. Sloot, Implementation aspects of 3D lattice-BGK: boundaries, accuracy, and a new fast relaxation method, J. Comput. Phys., 150, 1999. - P. 482-501.

5. D. d’Humieres, I. Ginzburg, M. Krafczyk, P. Lallemand, and L.-S. Luo, Multi-relaxation time lattice Boltzmann models in three dimensions, Phil. Trans. R. Soc. Lond. A 360, 2002.

- P. 437-451.

6. S. Chen, D. Martinez, andR. Mei, On boundary conditions in lattice Boltzmann method, Phys. Fluids 8, 1996. - P. 2527-2536.

7. Q. Zou and X. He, On pressure and velocity boundary conditions for the lattice Boltzmann BGK model, Phys. Fluids 9, 1997. - P. 1591-1598.

Сидоренко Борис Владимирович

Технологический институт федерального государственного образовательного учреждения высшего профессионального образования «Южный федеральный университет» в г. Таганроге.

E-mail: sidorenkobv@mail.ru.

347928, г. Таганрог, пер. Некрасовский, 44.

.: 8(8634)371-606.

Кафедра высшей математики; аспирант и ассистент.

Sidorenko Boris Vladimirovich

Taganrog Institute of Technology - Federal State-Owned Educational Establishment of Higher Vocational Education “Southern Federal University”.

E-mail: sidorenkobv@mail.ru.

44, Nekrasovskiy, Taganrog, 347928, Russia.

Phone: 8(8634)371-606.

The Department of Higher Mathematics; post-graduate student and assistant.

УДК 518.5.001.57

Т.В. Камышникова

ВЫВОД ДВУМЕРНОЙ МОДЕЛИ РАСПРОСТРАНЕНИЯ ЗАГРЯЗНЯЮЩИХ ПРИМЕСЕЙ В МЕЛКОВОДНОМ ВОДОЕМЕ

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

Модель «мелкой воды»; адвективно-диффузионное уравнение.

T.V. Kamyshnikova

RECEPTION OF TWO-DIMENSIONAL MODEL POLLUTING IMPURITY DISTRIBUTION IN SHALLOW WATER BASIN

Water basin hydro chemical characteristics ’ studying is actual problem because of variability as possible deviations ’ reflection in water basin functioning. The best method of this problem solving is mathematical modeling.

"Shallow water" model; the advection-diffusion equation.

Если на практике нет необходимости решать полные трехмерные уравнения,

, ,

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

Модель распространения загрязняющих примесей в мелководном водоеме включает в себя гидродинамическую задачу мелкой воды и задачу переноса при. -динамики мелкого водоема, которая была рассмотрена вне рамок данной работы, и результаты которой используются в данной работе при решении задачи распространения примеси в мелководном водоеме [3], [6]. Рассмотрим пространственнотрехмерную систему уравнений распространения загрязнений:

uX + v'y + Ч = 0 (!)

u'+(u 2)Х+(uv)y+(uw)Z =

= -Pp'x -v'x +u'yy +u'z'z ) + 2Q(vsin^- wsln#), (2)

vt+(uvX+(v 2)y+(vw)Z =

= -PPy -Py + + v”yy +v”z ) - 2Qu sin^ (3)

+(uw)X+(vw)'y +(w X =

= -—Pz -¥z +П{( + wZy + < ) + 2^ucosû

p px yy ’

(4)

Уравнение (1) - уравнение неразрывности, уравнения (2-4) - уравнения На-вье-Стокса для вязкой (в линейном приближении) несжимаемой (плотность р = const) жидкости во вращающейся с угловой скоростью

Q = Q^cosû- j + sinû-^ ),

где 1, ], к - единичные орты, р - полное гидростатическое давление; р - гравитационный потенциал; Т] - первый коэффициент вязкости в одном поле тяжести.

уФ = - g = - gk = const; Ро =Ро( л, y, Z, t); р = р0 + pg (д- z); yp=g(+д'у!-k), -h^z^д,

где д( X, y, t) - поднятие уровня свободной поверхности жидкости по отношению к невозмущенному состоянию; h(X, y, t) - высота столба жидкости под невозмущенной поверхностью.

Граничные условия:

du dv ,

Л^ = -тх, Л^ = -Ту на z(x,У>*Л dz dz

(5)

где Тх ,Т - составляющие касательного трения ветра. На дне и боковых границах условия прилипания:

и = V = w = 0. (6)

:

и (х, у, ¿,0) = и 0, v( х, у, г ,0) = V о,

w( х, у, ¿,0) = w 0, р( х, у, ¿,0) = р0, (7)

Интегрированием уравнений неразрывности и Навье - Стокса (1-4) по вертикальной координате г от — И до £ получены уравнения гидродинамической модели, учитывающей испарение жидкости и выпадение осадков ((w / р) - испаряющейся - выпадающий в виде осадков в единицу времени слой) [4]:

/ W

И, + и, + V. + — = 0,

(8)

и:+

Щи" + и" )-ки + И(f) + 2QVsin#

р\ xx yy ) И '•'sJx 5

ии

+

VU_

p

^wU TT , + C— + gHgx

pH

и

(9)

V'+

( UV1 f VV )

+

1 И ) x 1 И )

wV + C—

y РИ

+gHg'z

= n(V" + vyy)- к— + И (f) y + 2Q V Sin

\ xx yy ) JJ s >y ■>

p И

(10)

, -

верхности жидкости отнесены на счет обобщенной силы ветра о поверхность fs , а на дне - на счет обобщенной силы трения о дно

—Д- (ц + VJ.), и = |udz, V = |vdz, и = и(x, y, ¿, г) и V = v(х, у, г, г)

Н —И —И

- горизонтальные компоненты вектора скорости жидкости в точке (х, у, г) в

момент времени г . Н = И + £- полная глубина. Для описания поля некоторой субстанции обычно используется трехмерное адвективно-диффузионное уравнение [4].

— + V УБ — оБ = ^7 Б + — V— + /, (11)

дг дг

дБ

дг

где Б - концентрация субстанции; V = (и, V, w) - вектор скорости водного потока; Ц,У - коэффициенты горизонтальной и вертикальной турбулентной диффузии;

О Б - характеризует взаимодействие вещества с водой; / = / (х, у, г, г) - функция источников (стоков) веществ. Примесь предполагается пассивной, тх. задаваемое поле скорости (и, V, w) не зависит от Б. Рассмотрим краевую задачу для

уравнения (11) в области Q = {О(х, у, г), —И < г <£,г > 0} , где боковая поверхность О является объединением двух частей: О = О1 и О; О1 (х, у, г) -поверхность дна, О(х, у, г) - боковая цилиндрическая поверхность.

Запишем граничные условия по вертикали

V = 0, г = £( х, у, г), (12)

дп

V = 0, г = —И( х, у, г), (13)

дп

п - . :

= 0. (14)

дп

Начальное состояние характеризуется фоном

Б|г=0 = Б0. (15)

Остановимся на двумерном случае. В области О2 ={х, у е О, г > 0} поле концентрации примеси описывается уравнением [4]:

|<Б>+дхИ+дуН+ОБ)=м7ЧБ)+(/), <>6>

£ _ £

(Б) = |Б(х, у, г, г^г = БН, {/) = |/(х, у, г, г^г, которое получается ин-

— И — И

тегрированием по вертикали трехмерного уравнения (11) по вертикальной коорди-

дБ

нате г от —И до £ с условием---- = 0.

дп дО

Недивергентная форма этого уравнения:

дБ _дS _дS ё „-1*72^4 Г — + и — + V-—Г = рИ 1У2(БИ) + /,

Эt дх ду

где И = И + д. Иногда его используют в виде

дБ _д£ _д£ ^ тт-Ы^тг, ,

— + и — + V-—гБ = ^И У\8И) + /,

Эt Эх ду

(17)

(18)

где Б = Б, V = (и, V) = V. Действительно, интегрируя слагаемые уравнения (11) по вертикали почленно будем иметь:

І—Сї = — Ґ( БСї

і д 3/1 і

дґ[

. д 1 д-------аї =—

дх

- к і

{ дх дSv

Ь

I пБСі

<*)

д/

ї=д

\—Сї = —І \ vSdI к дУ ду І -{

д#( о, --^-(uS )

дх

ду ’

-дт ^)

д/

-- ^)

дх ’

ї=і

1=і

ї=-к

--(vS)

ду

ї=-к

grдSw д

-------Сї = —

-Т дУ дї

- к Ід2 S

Іш>Сї = - (wS)|ї=? - (wS)|ї=-к,

- к І 2

- Сї = — дх дх

V-к

\

дх

д 2 ^ д

Сї = —

ї ^Сї

J дУ

V - к (І

(і

V- к ^ У

д_#( дSЛ дх

- к дУ дУ V - к дУ

-І дУ V

К

дУ

V дх у ( ''

V дУ )

ї=і

ї=і

дк

дх

дк

дУ

ґдSл

V дх у

(дsЛ

V дУ)

?=- к

ї=-к

( дS Л . д ( і дsл Л (дS ' ґ дS Л

V— Сї =— І V—Сї = V - V

V дї) дї 1 -к дї ) V дУ ) ї=І V дї )

:=-к

ь ь ь

I oSdI = о | SdI ; | /Сї = /.

- И - И - И

С учетом граничных условий для компонентов скоростей задачи (1-4) и задачи (11-14) получим двумерное уравнение (18). Решение уравнения (18) для однородной несжимаемой жидкости ищется в области О с границей дО на участке дО1 - которой жидкость втекает, дП2- вытекает, дП0- непроницаемая часть границы [5], [6]. Имеют место граничные условия: на твердой границе дО0 ставится условие непротекания:

д8

= 0, (19)

дПо

дп

где П - нормаль к границе дО0. На участке дО1 стока в водоем задается значение

На участке выпуска воды, или на открытой границе дО 2 ставится условие III

:

где Vn - скорость течения по нормали П к дО2. Решение ищется при начальных

Входящие в уравнение (18) и в граничные условия (19-21) компоненты вектора скорости и форму свободной границы можно считать известными функциями [6].

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Гидрометеорология и гидрохимия морей СССР. - СПб.: Гидрометеоиздат, 1991. - Т.5.

2. Вольцингер Н.Е., Клеванный КА., Пеликовский ЕМ. Длиноволновая динамика прибрежной зоны. - J1.: Гидрометеоиздат, 1989. - 271 с.

3. Васильев B.C., Сухинов А.И. О моделировании ветровых течений в Таганрогском заливе с использованием квазиоптимальных сеток // Областная научно-техническая конференция, посвященная дню радио. - Ростов-на-Дону, 1993. - 26 с.

4. Васильев B.C., Целых А.Н. Принятие прогнозных решений в экологических задачах на основе методов численного моделирования. - Ростов-на-Дону: Изд-во Северо-

, 1999. - 48 с.

5. Колдоба А.В., Повещенко Ю.А., Самарская Е.А., Тишкин В.Ф. Методы математического моделирования окружающей среды. - М.: Наука, 2000. - 391 с.

6. . . -грязняющих примесей от автотранспорта в условиях городской застройки: дис. на соискание ученой степени канд тех наук. - Таганрог, 2003. - С. 182-191.

Камышникова Татьяна Владимировна

Технологический институт федерального государственного образовательного учреждения

высшего профессионального образования «Южный федеральный университет»

в г. Таганроге.

E-mail: kamyshtata@mail.ru.

347928, г. Таганрог, пер. Некрасовский, 44.

Тел.: 8(8634)371-606.

Кафедра высшей математики.

Kamishnikova Tatiana Vladimirovna

Taganrog Institute of Technology - Federal State-Owned Educational Establishment

of Higher Vocational Education “Southern Federal University”.

E-mail: kamyshtata@mail.ru.

44, Nekrasovskiy, Taganrog, 347928, Russia.

Phone: 8(8634)371-606.

The Department of Higher Mathematics.

(20)

данных:

(22)

- 236 c.