Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

THE SOLUTION OF THE TURBULENT DIFFUSION SEMI-EMPIRICAL EQUATION IN PROBLEMS OF POLLUTING IMPURITY TRANSFER BY GAUSS APPROACH

Duysebekova K.S. 1 Tayzhumanova Z.A. 1
1 Kazakh National University n.a. Al-Farabi
The analysis of the solution the semi-empirical solution of turbulent diffusion in problems of transfer polluting impurity is carried out by Gauss approach. It is supposed that there is a certain occupied industrial point in which is available one or several sources of the pollution which arrangement is determined by coordinates x, y, z. For calculation of impurity’s concentration average from a dot source used the solution semi-empirical equation by Gaussian function impurity’s distribution, received by a method of function of Green. Further the normalization of task’s key parameters is carried out, basic data are defined. By means of the received equation the model of an quantity’s assessment of the aerosol polluting substances arriving from a source with final and continuous duration of action is made. Thus, the considered settlement and analytical model of a technique of an concentration’s assessment of impurity is applicable to applied problems of expeditious control of the industrial region’s condition.
concentration of impurity
Gauss function
equation of turbulent diffusion
1. Bezuglaya E.U. Monitoring of a state pollution the atmosphere in the cities. Leningrad, Gitrometeoizdat, 1986. 200 p.
2. Makukhin V.L., Potemkin V.L. Modeling of transfer and transformation polluting impurity. Irkutsk, Russia 2012. Р. 286.
3. Naatz V.I., Naatz E.I. Mathematical models and numerical methods in problems of environmental monitoring the atmosphere. Moscow 2010. Р. 101–117.
4. Naatz V.I., Naatz E.I., Ryskalenko R.A. The parameterized models of the theory transfer in problems environmental monitoring the atmosphere and the principle of a minimax. the Stavropol state university, Stavropol 2009. Р. 132–133.
5. Koishybayev N.A. Assessment of influence physical and chemical processes on an ozone layer of Earth, and also change a condition an ozone layer. Republican state enterprise Kazakh research institute of ecology and climate. Almaty, 2006. Р. 133.

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

Пусть q(P, t) – функция, значение которой в момент времени t в точке P(x, y, z) совпадает со значениями мгновенной концентрации примеси, переносимой в атмосфере потоками воздуха. Предполагается, что функция q(P, t) непрерывно дифференцируема по x, y, z, t. Полуэмпирическое уравнение турбулентной диффузии записывают в виде [3]:

duys01.wmf (1)

где S – источник примеси, находящийся в точке P0 (x0, y0, z0) и производящий мгновенный выброс загрязняющих примесей в момент времени t0 в количестве Q0.

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

duys02.wmf (2)

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

а) момент времени, когда источник производит выброс загрязняющих веществ (ЗВ) – t0;

б) координаты источника P0 (x0, y0, z0) (м);

в) расстояние от источника до пункта наблюдения R (м) – расстояние от точки P0 – источника выбросов загрязняющих веществ до точки Р – пункта наблюдения.

г) Q0 (10–3 кг/с) – количество загрязняющих веществ, выброшенное источником в начальный момент времени t0;

д) турбулентность атмосферы в пограничном слое, характеризуемая коэффициентом турбулентной диффузии К = {Kx, Ky, Kz} (м2/с);

е) скорость ветра V = {Vx, Vy, Vz} (м/с).

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

duys04.wmf duys05.wmf (3)

Далее проводим нормирование переменных задачи. Положим, что

duys06.wmf (4)

где

duys07.wmf (5)

duys08.wmf, duys09.wmf – нормированные величины, принимающие значения в интервале (0,1). Пронормируем переменные и распределения задачи следующим образом: duys10.wmf, где duys11.wmf duys12.wmf duys13.wmf Аналогично получаем duys14.wmf, где duys15.wmf duys16.wmf duys17.wmf Коэффициенты турбулентной диффузии и скорости ветра пронормируем как обычно:

duys18.wmf duys19.wmf

duys20.wmf duys21.wmf (6)

где K* = max{Kx, Ky, Kz};

duys22.wmf duys23.wmf

duys24.wmf duys25.wmf (7)

где V* = max{ Vx, Vy, Vz }.

Далее нормируем переменные

duys26.wmf duys27.wmf duys28.wmf

duys29.wmf (8)

T = t* – t0 – временной интервал, в течение которого в пункт наблюдения будут поступать загрязнения.

С учетом вышеизложенного уравнение (5) можно записать таким образом:

duys31.wmf (9)

Рассмотрим случай, когда в точке P0 (x0, y0, z0) работает источник в течение конечного промежутка времени [ζ, ζ0 + T] и в точке наблюдения Р(x, y, z) происходит процесс накопления загрязняющих веществ. Предполагается, что в пункте наблюдения ведутся замеры (прием) концентрации поступивших примесей от источников в моменты времени t ∈ [t1, t2]. Функция S(P0, ζ) имеет смысл интенсивности источника, при этом функция S(P0, ζ) Δζ (10–3 кг/м3) определяет количество вещества, выброшенного в атмосферу за элементарный интервал времени Δζ в окрестности ζ. Согласно рассмотренной теории, импульс возмущения загрязняющих веществ S(P0, ζ) будет приниматься в точке Р в интервале [duys32.wmf]. В силу этого начальное возмущение S(P0, ζ) будет приниматься в точке Р во время duys33.wmf, а последнее – во время duys34.wmf. Объединение этих интервалов дает интервал duys35.wmf. Таким образом, становится видно, что интервал [ζ, ζ0 + T] соответствует интервалу duys36.wmf. Моменты времени ζ и duys37.wmf связаны простым преобразованием:

duys38.wmf (10)

и, следовательно, интервал [ζ, ζ0 + T] соответствует интервалу

ζ0 + t* – τ1 ≤ t ≤ ζ0 + T + t* + τ2.

В итоге имеем:

ζ0 + t* – τ1 ≤ t ≤ ζ0 + T + t* + τ2. (11)

Неравенство (11) определяет время возможного приема импульса возмущения загрязняющих веществ или, иными словами, время прохождения импульса возмущения загрязняющих веществ через точку Р – пункт наблюдения.

Теперь нужно построить интегральную форму Q(P0, P, t), определяющую количество загрязняющих веществ, накопившееся за время работы источника в точке Р. Для этого выберем на интервале [ζ, ζ0 + T] точку ζk и интервал Δζk следующим образом:

ζk = ζ0 + k∆ζ, (12)

где duys39.wmf и ∆ζk = ζk + 1 – ζk.

Введем в рассмотрение точку duys40.wmf и примем в точке P0 выброс:

duys41.wmf (13)

Тогда в точке Р реакция будет следующей: [3]

duys42.wmf (14)

В итоге можно записать:

duys43.wmf (15)

Это интегральная сумма, остается перейти к пределу при ∆ζ → 0 и сделать замену m(t) → ζ(t), получим:

duys44.wmf (16)

где ζ0 ≤ ζ < ζ0 + T,

ζ0 + t* – τ1 ≤ t ≤ ζ0 + T + t* + τ2

или t1 ≤ t ≤ t2, если принять

t1 = ζ0 + t* – τ1;

t2 = ζ0 + T + t* + τ2; ζ(t) = t – t*. (17)

Далее рассмотрим вариант, когда источник работает непрерывно, а в точке Р фиксируется количество загрязняющих веществ в течении определенного промежутка времени [t1, t2]. Необходимо определить

duys45.wmf (18)

Выражение (18) задает усредненное значение концентрации загрязняющих веществ, приходящееся на один кубический метр и накопившееся за время [t1, t2]. В данном случае интервал [t1, t2] задан и требуется определить соответствующий ему интервал [ζ1, ζ2]. Согласно (13) t1 = ζ0 + t* – τ1, t2 = ζ0 + T + t* + τ2. Отсюда находим ζ0 = t1 + τ1 – t* и ζ0 + T = t2 – t* – τ2. Обозначим ζ1 = ζ0 и ζ2 = ζ0 + T. В результате получаем ζ1 = t1 + τ1 – t* и ζ2 = t2 – t* – τ2. В итоге имеем следующие расчетные формулы:

duys48.wmf (19)

duys49.wmf (20)

ζ1 = t1 + τ1 – t*;

ζ2 = t2 – t* – τ2; t1 ≤ t ≤ t2. (21)

Проведем расчет концентрации примесей, поступающих в точку наблюдения от источника согласно методике (16), (18). Для этого необходимо, прежде всего, задать исходные данные и провести нормировку основных переменных и полей задачи. Затем необходимо вычислить прихода Р – максимума загрязняющих веществ, поступающих от источника. Затем задаем разбиение отрезка [t1, t2] и формируем массив {tl}, duys50.wmf duys51.wmf duys52.wmf [4]. После нормирования переменных duys53.wmf duys54.wmf duys55.wmf duys56.wmf duys57.wmf duys58.wmf dζ = (T – t0), duys59.wmf можно записать:

duys60.wmf (22)

Далее рассмотрим функцию источника. Допустим, что

duys61.wmf (23)

duys62.wmf

В итоге имеем [3]:

duys63.wmf (24)

duys64.wmf (25)

duys65.wmf S0 = Q0. (26)

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

Рецензенты:

Джомартова Ш.А., д.т.н., профессор, кафедра «Информационные системы», механико-математический факультет, Казахский национальный университет имени аль-Фараби, г. Алматы;

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

Работа поступила в редакцию 17.01.2014.