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

доктора технических наук
Беляев, Николай Николаевич
город
Киев
год
1996
специальность ВАК РФ
05.13.02
Автореферат по информатике, вычислительной технике и управлению на тему «Математическое моделирование распространения загрязняющих веществ в водной и воздушной средах»

Автореферат диссертации по теме "Математическое моделирование распространения загрязняющих веществ в водной и воздушной средах"

КИЇВСЬКИЙ УНІВЕРСИТЕТ імені ТАРАСА ШЕВЧЕНКА

МАТЕМАТИЧНЕ МОДЕЛЮВАННЯ ПОШИРЕННЯ ЗАБРУДНЮЮЧИХ РЕЧОВИН У ВОДНОМУ ТА ПОВІТРЯНОМУ СЕРЕДОВИЩАХ

05.13.02.— Математичне моделювання в наукових дослідженнях

На правах рукопису УДК 519.63:502.5(203):628.19

БІЛЯЄВ Микола Миколайович

АВТОРЕФЕРАТ

дисертації на здобуття наукового ступеня доктора технічних наук

Київ—1996

Робота виконана в Дніпропетровському державному технічному університеті залізничного транспорту .

Науюэвий консультант: доктор фізимо-иатематичних наук, професор ХРУЩ Віктор Кузьмич Офіційні опоненти; член-кореспондеит НАН України,

доктор фізико-математичних наук, професор СКОПЕЦЬКИЙ Василь Васильович; доктор технічних наук, професор ЛАВРИК Володимир Іванович; дійсний член Академії технологічних наук України,

. доктор фізшю-матеиатичних наук,

професор СЕЛЕЗОВ Ігор Тімофійович

Провідна організація - Національний технічний університет України

“Київський політехнічний інститут”.

Захист відбудеться “ 3 0 ” 1996 р. о год. на

засіданні спеціалізованої вченої ради Д 01.01,20 при Київському університеті імені Тараса Шевченка за адресою: 252127, Київ-127, проспект Академіка Глушкова, 6, факультет кібернетики, вуд. ^ О )

З дисертацією мояаи ознайомишся в науковій бібліотеці Киїасьюзкі уніперентету імені Тараса Шевченка (вул. Володимирська, 58).

Автореферат розісланий " ё?^(Я&И>Я 1996 р.

Вчений секретар спеціалізованої вченої ради кандидат фізимо-маїематичних наук,

доіккт

ішько ТІМ.

ЗАГАЛЬНА ХАРАКТЕРИСТИКА РОБОТИ

Актуальність теми. Прогноз динаміки змін екологічного стану навколишнього середовища можливий лише на основі методу математичного моделювання. Підвищення якості прогнозної інформації вимагає врахування при розрахунках поширення забруднюючих речовин (ЗР) багатьох факторів, що робить задачу прогнозування далеко не тривіальною, існуючі нормативні методики прогнозування якості водного та повітряного середовищ не дозволяють врахувати форму русел, джерел забруднення, рельєф, різноманітні форми викидів та ін, а тому мають суттєву обмеженість і низьку якість у застосуванні до практичних ситуацій. В даній роботі розв'язано важливу наукову проблему розробки методів моделювання, які адекватно відтворюють найбільш важливі фізичні процеси гідродинаміки та переносу ЗР у водоймах, підземних водах, атмосфері та повітряному просторі виробничих приміщень. Розроблено комп'ютерні методи, які дозволяють при обмежених ресурсах ПЕОМ • на новому якісному рівні вирішувати проблемні задачі антропогенного забруднення водного та повітряного середовищ (ВПС).

Мета роботи. Розробка нових ефективних методів розв'язування широкого класу практичних задач, які стосуються проблем охорони навколишнього середовища; створення регуляторш« комп'ютерних моделей і комплексів програм для розрахунку на ІІЕОМ процесів поширення ЗР у ВПС, що дозволяють на. новому рівні прогнозувати якість вод, атмосфери і повітряного простору всередині виробничих приміщень. Ця,мета передбачає розв'язаш(^мступних_задач:

• побудова ефективних неявних різницевих алгоритмів з корекцією величини схемної похибки для інтегрування рівнянь гідродинаміки водних та повітряних потоків, переносу забруднень у

водному та повітряному середовищах і розробка на їх основі спеціалізованого програмного забезпечення обчислювального експерименту - комплексів програм (КП) для прогнозу якості вод і повітря;

• розв'язування на базі розроблених методів прикладних задач про забруднення проточних водойм складної геометричної форми, підземних вод через скид шахтних та промислових стічних вод;

• розв'язування праісгичіїих задач про розсіювання ЗР в атмосфері від різних джерел з врахуванням типу техногенного вибросу, рельєфу місцевості, форми джерела, присутності будівель, метеумов;

е дослідження закономірностей поширення ЗР в робочих зонах приміщень і підземних виробках у випадку стаціонарних, циклічних

■ ' ' ' • ’ ' • та залпових викидів;

« теоретичне обгруїпування, підтвердження достовірності і точності розрахунків за розробленими методами.

Методика дослідження грунтується на чисельному експерименті з математичними моделями різного рівня, котрі описують перенесения ЗР у повітряному і водному середовищах. При розв'язуванні рівнянь на ПЕОМ використовуються нові різницеві алгоритми чисельного інтегрування з корекцією величини схемної похибки у процесі обчислювального експерименту.

Наукова новизна. Внаслідок проведених в роботі досліджень одержано такі нові результати:

• розроблено комплексний підхід до математичного описання різноманітних фізичних процесів поширення забруднюючих речовин у водному та повітряному середовищах, який грунтується на балі дискретних балансових співвідношень для контрольного елемента,

• розроблено нові ефективні дискретні алгоритми для чисель-

ного моделювання одіїр- і двовимірних гідродинамічних течій у відкритих руслах. Вони дозволяють розв'язувати багатофакгорні задачі про забруднення водойм з детальним урахуванням особливостей русла: геометричної форми берегів, наявності приток, островів, водозаборів, регулювання стоку та ін.;

• розроблено ефективні неявні алгоритми для розрахунку поширення забруднюючих речовин у повітрі. На основі алгоритмів створено комплекс комп'ютерних методів прогнозування якості атмосфери і повітряного середовища у робочих зонах виробничих приміщень;

• розроблено ефективний алгоритм розрахунку процесів гідродинаміки та поширення забруднень в атмосфері з урахуванням рельефу місцевості, метеофакторів, форми вогнищ забруднення, наявності забудови, типу викидів;

• запропоновано нові методики побудови нелінійних квазі-

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

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

• розв'язано ряд важливих прикладних задач для Придні-. провською регіону стосовно дослідження перенесення забруднень у

воді та повітрі. '

На захист виноситься:

• чисельні моделі і алгоритми для розрахунку

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

о КП "Прогноз-І", "Прогноз-2", "Прогноз-3", "Аспірація" для прогнозування якості атмосфери;

в КП "Акваторія", "Водойма", "Водойма-2", "Виїмка", "Водовідстійник" для прогнозування якості поверхневих і підземних вод; ■

, • результати багатоваріантних розрахунків поширення ЗР у по-

верхневих і підзгмшгх водах, атмосфері та повітряному середовищі виробничих приміщень;

в обгрунтування ефективності, економічності розроблених програмних засобів.

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

Практична цінність. Розроблено чисельні моделі і алгоритми реалізовані у вигляді КП, ям на новому якісному і кількісному рівні дозволяють оперативно прогнозувати антропогенне навантаження на навколишнє середовище від промислових підприємств чи інших об'єктів з метою вироблення науково обгрунтованих нриродоохорон-

них заходів. Запропоновані нові комп'ютерні технології є економічними, універсальними і мають широкий "робочий" діапазон, зручні для практичного використання, дозволяють замінити в багатьох випадках дороге фізичне моделювання. Розроблена у роботі теоретична база і методологія побудови алгоритмів може бути використана для створення нових ефективних методів розрахунку розсіювання ЗР у різних природних середовищах.

Реалізшфі результатів роботи. Розроблені комплекси програй використовуються в регіональному проектно-дослідницькому інституті “Дніпрогіпроводгосп” при проведенні розрахунків в рамках програми розробки комплексу заходів з метою запобігання шкідливому впливу скиду шахтних вод ВО “Павлоградвугілля” на природне середовище. Розроблені методи та комплекси програм використовуються для розв'язання екологічні« задач на Павлоградському механічному заводі, в Дніпропетровському НДІ великогабаритних шин, Українському науковому центрі технічної екології (Донецьк), в КБ “Південне” (Дніпропетровськ), в проектному інституті “Дніпрогіпрошахт” (Дніпропетровськ), ДержНДІ ХП (Шосгка), ВО “Десна” (Шосгка). Впровадження результатів дисертації в практику підтверджено актами впровадження. Результати роботи використовувались при виконанні ряду господарчодоговірннх тем. Робота виконана в межах науково-технічної проблеми ДКНТ і Держплану СРСР 085.04 "Створити і впровадити ефективні методи і засоби контролю забруднення навколишнього середовища" та плану НДР ДДУ "Чисельне моделювання течій рідшш і гачу на базі змінно-трикупшх різницевих схем і розробка комп'ютерних моделей забруднення повітряного середовища та водного басеґніу Дніпра" (№ Держ.

per 00194V010128-22.03.94). Матеріали робота увійшли в навчальний посібник "Чисельний розрахунок поширення аерозольних забруднень", використовуються в навчальному процесі Дніпропетровського державного університету, Дніпропетровського державного технічного університету залізничного транспорту, Донецького державного університету.

Апробація роботи. Результати дисертаційної роботи були представлені на 20-ти Всесоюзних, Міжнародних конференціях і семінарах, у тому числі: науковому семінарі "Чисельні методи механіки суцільного середовища" під кер. проф. В.М.КЬвені, Інститут теоретичної та прикладної механіки СВ AH СРСР (Новосибірськ, 1986); Першій і Другій Всесоюзних конференціях •"Актуальні питання охорони навколишнього середовища від антропогенного вплину" (Севастополь, 1989,1S30 ); Всесоюзній конференції "Математичне

моделювання і експериментальні методи дослідження фізико-хімічних процесів у суцільних середовищах" (п. Рибачьє, 1989 );

Всесоюзному семінарі "Проблеми фізико-хімічних взаємодій в механіці суцільних середовищ" (Ужгород, 1989); Міжнародній конференції "Теорія наближення та задачі обчислювальної математики"

(Дніпропетровськ, 1993); Першому міжнародному симпозіумі

"Ecosystem Health & Medicine" (Онтаріо, 1994); 6-th International Conference on Numerical Methods" (Будапешт, 1994); Першому міжнародному конгресі ІАС'94 " Intern al'onaJ Aerospace Congress" (Москва, 1994); 12-му Всесвітньому конгресі "Agricultural

Engineering" (Мілан, 1994); 45-th International Astronautical Congress (ІЄрусалим, 1994), 1-й Міжнародній конференції "Чисельні методи в гідравліці і гідродинаміці" (Донецьк, 1994): Міжнародній научно-

технічній конференції) "Екологія хімічних виробництв" (Сє-вєродонецьк, 1994); International Rubber Conference (Москва, 1994), International Congress "Water: Ecology and Technology" (Москва, 1994) Fifth International Congress of Fluid Mechanics (Каїр, 1995), науковому семінарі Інституту гідромеханіки HAH України (1995), науковому семінарі Інституту природокористування та екології НАН України (Дніпропетровськ, 1995).

Комплекси прикладних програм комп’ютерного моделювання переносу ЗР у повітрі демонструвались на виставці-ярмарці завершених нвукових робіт вузів України (Київ, 1993).

Публікації: за темою диссертації опубліковано 67 робіт.

Обсяг і структура роботи. Дисертація складеться з вступу, шести глав, висновку, списку літератури і додатку. Робота містить 345 сторінок машинописного тексту, 168 рисунків, 6 таблиць. Список літератури включає 306 назв.

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

Згяіст роботи.

У вступі визначена актуальність і мета роботи, викладасться короткий зміст диссертації і обговерено ста» досліджуваної проблеми.

У першій паві розглянуто математичні моделі розсіювання ЗР у атмосфері і повітрі виробничих приміщень. Для опису розсіювання ЗР у атмосфері п локальному масштабі викпристоиуегься нЕ-сгупне

рівняння відносно концентрації <р забруднюючої речовини

^ + Ці + 012 + ^ + оф = с)Іу(цУф) + £Чі(‘)3 (г - г,(і>), (1)

оі ох оу аг м

де и, V, іл' - компонента вектора швидкості повітряного потоку; -швидкість гравітаційного осідання частинок; сг - коефіцієнт хімічного розкладення забруднення; р.=(цх, уіу, - коефіцієнт турбулентної дифузії; чі - інтенсивність точкових джерел забруднення; П=(хі, у;, г,) -місцеположення точкових джерел забруднення; б(г-гі) - дельта-функція Дірака.

На основі рівняння (1) та розробленого алгоритму створено КП "Прогноз-1".

Для математичного моделювання далекого переносу забруднень у атмосфері використовується рівняння .

1^ + ^ + ^ + „ = (1І''(,і7^ + Е9,(08 X у-у ^

01 ОХ Оу М .

, Н

де ф(х, у) = — (<р(х,у,г)бг - усереднене по висоті Н значення конН ь

центрації ЗР у одиниці об'єму повітря. Коефіцієнт а в моделі (2) визначає наступні процеси: хімічне перетворення ЗР; вимивання опадами та гравітаційне осадження частинок ЗР; турбулентний обмін у приграничній зоні біля підстилаючої поверхні.

На основі математичної моделі (2) та розробленого алгоритму створено КП "Прогноз-2".

Для визначення полів концентрації забруднень всередині виробничих приміщень та робочих зон, розрахунку оптимального об'єму видаленого місцевим відсисшшям повітря для вибраної схеми установки відбору з урахуванням взаємною розташування техіюло-

гічного устаткування, використовується рівняння дф дікр Зуф д(\у — н>,)ф

ді * дх + ду дг ^ ^

= <1іу(|і Уф) + 2ч,0)8(г-г,(0) - " гі)

1-І £4 « і >

яке доповнюється формулою розрахунку швидкості повітря за рахунок роботи місцевих відсосів

де 4$) - інтенсивність точкового джерела вихиду ЗР; - об'ємна витрата точкового просторового стоку (місцевого відсосу); ^=(х^ yj, гр - координати відсосу. На базі моделі (3) створено КП "Аспірація".

У першій главі обговорюються моделі і підходи до визначення метеорологічних і фізичних параметрів, залучених у рівняння (1)-(3).

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

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

режешш монотонних профілів розподілу шуканих параметрів. Проведено аналіз відомих алгоритмів інтегрування за часом. Розглянуто нову організацію алгоритмів інтегрування рівнянь (1)-(3) за часовою координатою. Виділено неявну змінно-трикутну схему (ЗТС), яка забезпечує ефективний розрахунок, і показано методи підвищення порядку апроксимації за рахунок алгоритмів корекції потоків і використання антадифузійного кроку.

У п. 2.5 побудована базова чотирикрокова ЗТС для чисельного інтегрування двовимірного рівняння переносу забруднюючої речовини. В опереторній формі різницеве рівняння для моделі (2) на кожному кроці розщеплення моше бути записане у вигляді:

■ _ (4)

= \ «<Рк + М^срс + М^Фк + М^Фе) + ІД б, ^

де к = п+1/4; п+1/2; п+3/4; п+1 - верхній часовий шар; с = п; п+1/4; п+1/2; п+3/4 - нижній часовий шар;

ц - збурений коефіцієнт дифузії., ри = ц, -^1 + •

Різницеві оператори Ь* апроксимують конвекіивні похідні з врахуванням напрямку переносу ЗР

, . + ,.п+1 + п+1

Т+п+1_ 5и ф «і+иФц “іі'Рі-іо

----------------дїї-------— (5)

Швидкість переносу ЗР подасться у вигляді суми знакосталих

і?

+ и + |и| - и-|и

величин и =—=—гг—

¿л | А

Для апроксимації других похідних використовуються два часових шари (за методом В.К Сгульєва). На верхньому часовому шарі апрок-симусться та частина похідної, яка "узгоджується" з напрямом розділеної різниці при заміні конвективно! похідної, в опералгорній формі ця апроксимація затісується так.

На кожному кроці розщеплення враховується лише один напрямок переносу ЗР за координатними напрямками.

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

Різноманітна фізична постановка задач розглянутого класу (поширення ЗР від постійно діючих джерел викидів, циклічно діючих і тій.) припела до необхідності виділення у окремий розрахунковий крок джерельних доданків, що дозволило зробити алгоритм розрахунку більш гнучким і універсальним. Даний варіант ЗТС для рівняння (2) побудовано у п. 2.6. '

Чисельний розв'язок тривимірних задач переносу забруднень з допомогою ЗТС будується аналогічно описаній яиіце схемі; при ял-

роксимації конвекгавних похідних у розробленій схемі використовується проміжний шар за типом схеми Кранка-Нікольсона. В опера-торній формі різницеві рівняння мають вигляд:

де 0^1 - параметр, що визначає проміжний шар при обчисленні конвективних похідних. На кожному кроці розщеплення знак "±" вибирається у відповідності з напрямком конвективного переносу ЗР.

Для створення гнучкого алгоритму розв'язування задач з різною фізичною постановкою виділено окремий п'ятий крок для розрахунку змін концентрації ЗР під діянням джерел і стоків у розрахунковій області. Розрахункове співвідношення для цього кроку алгоритму одержано за методом балансу шляхом інтегрування рівняння з дельта-функцією на різницевій комірці

У дискретному витяді дельта-функція Дірака "розмазується" по об'єму різницевої комірки зі збереженням сумарної кількості викинутого (для джерела) чи відсмоктаного (для стоку) забруднення.

При значенні 5=1 схема мас перший порядок точності, а при 4^1/2 - другий; послідовність розв'язування різницевих рівнянь, незважаючи на їх некомутативність, забезпечує другий порядок точності за часом, що відповідає схемі, запропонованій ПІ.Марчуком При

і (Ъ с N 1

-----------+| 14 + ЬУ + Мі)|фцк 5.+ Фіікі1 ~ + 4---2

=|Ді (міс +м^у+м!) фЦк+|м++.м+ + МІ) <рЦк ^

(їцк+фцк)

односторонній апроксимації похідних забезпечується монотонність профіля розв'язку, однак'має місце велика схемна в'язкість. Розроблені у п. 2.4 методи корекції потоків дозволяють зменшити величину схемної в'язкості при односторонній апроксимації конвекгивних похідних. Методика розрахунку корекції потоків виділена у окремий модуль розроблених програмних засобів. Необхідно відзначити, що для розглянутого класу задач існує нестача у вхідної інформації (задания коефіцієнтів а, ц і т.д.), що обумовлено низкою причин і, зокрема, слабкою системою моніторингу навколишнього середовища. У цьому випадку при проведенні обчислювального експерименту, а також при виконанні "пілогших" розрахунків мас сенс вести розрахунок за схемами 1-го порядку точності без залучення програмного модуля корекції потоків. Як показали гєсгоеі розрахунки і порівняння з експериментальними і розрахунковими даними інших авторів, розрахунок за схемами 1-го порядку забезпечує необхідну для інженерних розрахунків точність. При наявності повної надійної бази вхідних даних для моделі з метою детального описання розподілу домішок має сенс виконувати обчислювальний експеримент з обліком корекції потоків.

У п. 2.8, 2.9 проведено методичні дослідження побудованих різницевих алгоритмів. Методом Неймана і методами теорії диференціальних наближень досліджено питання апроксимації та стійкості алгоритмів. Проведено порівняння результатів розрахунків за розробленими алгоритмами з аналітичними розв'язками рівнянь конвективного і дифузійного переносу ЗР, з розрахунками В.М.Ковєиі, експериментальними даними Р.В Озмідопа (рис. 1,а). Проведено методичні досліди чисельних алгоритмів, які дотопили зробити висновок про високу точність, надійність і економічність розроблених ЗТС для чи-

ссльного розв'язування задач про поширення ЗР у повітряному і водному середовищах.

У третій главі розглянуто (п. 3.1) загальні вимога, що ставляться до вентиляційних мереж промислових цехів, і підходи до теоретичного описання процесів вентиляції.

Побудовано комп'ютерні моделі конкретних задач і розглянуто техніку моделювання розсіювання ЗР у виробничих приміщеннях з допомогою КП "Аспірація". ДослЬшго розподіл забруднення в робочій зоні і визначено ефективність роботи місцевого відсосу при різному розташуванні його по відношенню до джерела утворення шкідливої речовини (п. 3.2). Розв'язок задачі приведено для одного і двох підсосів, розташованих різним чином по відношенню до джерела забруднення.

Ряд джерел викиду ЗР на виробництві діє циклічно на протязі технологічного процесу. У п. 3.3 розглянуто моделювання процесу видалення шкідливих домішок з робочої зони циклічно діючого устаткування (дозаторів). Методом обчислювального експерименту досліджено вплив розташування місцевих підсосів на якість повітря на робочих місцях; досліджено ефективність роботи відсосів; про-слідковано зміни інтенсивності забруднення повітря на робочих ' місцях з плином часу. У п. 3.4 проведено розрахунок поширення забруднень від системи розподілених у просторі приміщення джеред-забруднювачів. Илюструегься можливість застосування КП для проптозу якості повітряного середовища у приміщенні на сталі проектування місцевої витяжної мережі для джерел із змінною інтенсивністю викиду, розташованих "складним" чином у приміщенні. Зроблено аналіз якості повітря і досліджено ефективність роботи відсосів

У п. 3.5 розглянутої моделювання поширення ЗР від рухомого у робочому приміщенні джерела-забрудшовача зі сталою та змінною в часі інтенсивністю викиду. Виконано прогноз формування зони забруднення повітряного середовища і проведена оцінка ефективності роботи витяжних підсосів, розташованих по-різному в приміщенні.

Відомо, що забруднення повітря у робочих приміщеннях молсс бути визване залповим викидом від технологічного устаткування. У п. 3.6 показана можливість застосування КП "Аспірація" для розв'язання задач такого класу. Методом обчислювального експерименту "реалізується" "сценарій" технологічного викиду, визначається ефек-лгвність роботи відсосів у залежності від їх положення відносно джерел викиду забруднювачів, досліджується динаміка змін інтенсивності забруднення повітря в приміщенні з плином часу (рис. 4,а).

В четвертій гладі розглянуто математичні моделі конкретних задач і техніка дослідження процесів розсіюваній забруднень у атмосфері з допомогою КП "Прогноз-1", "Прошоз-2", "Пропшз-3".

Моделювання процесу далекого переносу ЗР у атмосфері (п. 4.1) здійснено в задачі про поширення антропогенних викидів міст і промислових центрів одного з найбільших регіонів України - Прид-ніпров'я-Донбас. Як джерела забруднення вибрані промислові населені пункти Придніпров'я і Донбасу з населенням понад ЗО тис. чоловік. Подаються результати забруднення атмосфери для різних напрямків вітру.

Однією з важливі« прикладних задач є прогноз забруднення повітряного басейну при аваріях чи інших надзвичайних ситуаціях, пов'язаних зі значним викидом токсичних речовин у атмосферу. Виразною особливістю таких процесів є їх нестацюиарність. Можли-

вість застосування КП "Прошоз-2" дня розв'язування задач даного класу проілюстрована у п. 4.2 і 4.3. В п. 4.2 розглянуто поширення шкідливих речовин у атмосфері при "залповому" викиді у випадку імовірної аварії на великому газосховищі. Досліджено динаміку поширення хмари ЗР при змінних напрямках переміщень повітряних мас. У я. 4.3 представлено результати моделювання прогнозу формування зони забрудненій повітря при значній аварії на ряді промислових об’єктів за допомогою КП "Прогноз-2". Розглядається статистичний "сценарій" забруднення повітряного’середовища при нестаціонарному процесі викиду забруднювачів.

У п. 4.4-4.6 показано застосування КП "Прогноз-Г для розв'язування задач у різних фізичних постановках для прогнозу забруднення атмосфери і підстилаючої поверхні у випадку стаціонарного і нестаціонарного викиду ЗР в повітря з урахуванням просторового переносу домішок. Додатково до математичного обгрунтування разроблених різницевих схем (глава П) у п. 4.4 зроблено порівняння результатів розрахунку за розробленими алгоритмами з даними експерименту і розрахунковими методиками М.Е. Берлянда. У даному параграфі представлено результати розрахунку просторового розсіювання домішки в атмосфері від стаціонарного джерела забруднення, розташованого на висоті Н. У п. 4.5 досліджено нестаціонарний просторовий перенос пилового забруднювача у атмосфері після вибуху. '

В п. 4.6 розроблено алгоритм розрахунку розсіювання ЗР в умовах забудови. Для опису гідродинаміки течії між будинками використовується модель вихрової течії нев’язкої рідини

дм ^¿Зсои^дсоу

— 4----Н---г

Й дх бу

(7)

а процес розсіювання ЗР описується рівнянням турбулентної дифузії домішки. Для чисельного інтегрування рівняння переносу зашгхре-ності використовується ЗТС розщеплення виду (4) без дифузійних, складових. Для розв'язування рівняння Пуассона використовується ідея встановлення розв'язку за фіктивною часовою координатою, а алгоритм розв'язування розглянуто в п. 5.10. При розрахунках припускається, що відрив потоку в кутових точках споруд приводить до утворення вихорів, інтенсивність яких розраховується. Розрахункова область формується на прямокутній різницевій сітці, де використовуються маркери, що модулюють геометричну форму будівель. Таким чкнон, форма біудівель подається як прямокутні багатокутники різної геометричної форми. На базі даної моделі створено КП "Прогноз-3", який може бути використаний також для розрахунку розсіювання ЗР з урахуванням обтікання рельєфу та провітрювання вкробничігх приміщень. Представлено результати розрахунку розсіювання ЗР в умовах забудови на промманданчику аглофабрики. Досліджено вплив типу викиду (низький, високій, вентиляційний), геометричної форми будівель на формування, інтенсивність зони забруднення атмосфери та процес аерації проммайдалчкка.

В п. 4.7 розв'язано задачу про забруднення повітряного басейну від задимленого хвостосховища. Досліджено вплив рельєфу місцевості* метеофаі ¡торів і місцеположення хвостосховища на фор-мувшшя, інтенсивність та розміри зони забруднення. 'Наведено результати дослідження забруднення атмосфери від хвостосховища, рочтанюишюіх) в ярузі, з урахуванням геометричної форми вогнища

забруднення та форми яруги.

В п. 4.8. наведено розв'язування задачі про провітрювання підземної камероподібної виробки складної форми при роботі в ній двох дизельних агрегатів. Визначено поле концентрації домішок у виробці, форма та розміри утвореної зони забруднення. Розв'язано задачу про провітрювання виробки після залпового викиду ЗР Проведено порівняння результатів розрахунку з експериментальними та чисельними даними Г.В.Калабіна, А. А. Бакланова, П.В. Амосова.

В п. 4.9 досліджено динаміку забруднення атмосфери при рухові готової хмари в кар'єрі. Дослідження еволюції форми зони забруднення, Її інтенсивності та розмірів виконано з урахуванням ме-теофакгорів, геометричної форми кар'єру та форми утвореної вибухом хмари (рис. 2). Розрахунок проведено для двох різних за інтенсивністю вибухів. Проведено порівняння розрахунків з даними Б.В.Плюхіна, А.А.Фалейчик (СВ АН СРСР).

У-п'ятій главі побудовано різницеві алгоритми і КП розрахунку гідродинаміки течій і транспорту забруднень у акваторії рік та підземних водах; вироблена постановка і техніка моделювання широкого кола задач прогнозу якості вод на основі розроблених КП.

. Для описання гідродинаміки течії і транспортування забруднень у проточних водоймах у роботі розроблено три моделі. Перша модель орієнтована на розв'язування задач прогнозу якості води у акваторіях рік великої протяжності. Дана модель (п. 5.1, п. 5.2) грутусться на одновимірних нестаціонарних рівняннях гідравліки, у яких шуканими величинами с об'ємна витрата води С> і площа асивого перерізу 8

б

де и = ;? =

0

а. - нахил дна русла; додадок q враховує притік підземних, дощових чи тплих вод а також притік рік у русло; Ь» - характерна ширино русла .

Дані про нестаціонарну гідродинаміку течії використовуються в задачі про поширенім забруднень, котрі попадають у русло при технологічних викидах забруднених вод від промиспопЕїя підприємств, розташованих у басейні даної ріки. Дм розрахукху ко'щеятрації забруднення в руслі використовується рімюша '

^+^+СТф8=^0+^л+&Ф-*)+9,<г,б(х-5!?) г (9)

де 9 - коннк.траціх забруднення; о.= ох +а3 - коефіціпгг, що вра-

хсзус хімічішй чи біологічний розпад звбрзеткп. і ззвгасегиш юяцея-' традії забруднення у потоці при осадженій еешких зебр>днень ка дно русла; ед - концентрація забруднень при стоку дощових вод; а ¿(і)

кільхість дощових вод, ясі попадають у русло; ® - кількість забруд-

• / ' ' • •

некь, які попадактп. у русло в точці хі зосередженого вкккду забруднених вод; <р? - концентрація зебруцягозачіз у притоці; Ор - витргта води у пр:ттоці, яка моделюється зосередженим дкерелом з геометричною координатою х,. ■

На базі диференціальної моделі поширкгая забруднень по доз-•¿сін.і русла (9) і рівклкь (8) створено КП "Акваторія".

У п. 5.3 побудовано різницеві алгоритми розрахунку однозн-мірних нестаціонарних течій у відкрігпгх руслах. Описано два алгоритми розрахунку - хвішй і неявний. Оскільки вихідна система

рівнянь (8) с нелінійною, то для побудови неявного алгоритму застосовується лінеаризація рівнянь відносно приросту розв'язку на малому часовому кроці, виділяються матриці конвективного переносу, які представляються у вигляді суми знакосталих матриць, що дозволяє адекватно описати напрям процесу переносу збурень. Схема розщеплення для організації неявного алгоритму виглядає так

Дня забезпечення другого порядку точності за часом використовується двоциклічне розщеплення.

Для чисельного інтегрування даної системи використовується неявна схема "проти потоку". При розрахунку гідродинаміки течії, коли необхідно враховувати дрібномасштабні зміни розв'язку за часом, загосовується явний алгоритм, записаний у вигляді

ади2

ді

ми3 д_

61 + дх

(л'-ди3) =0

де А = А+ + А" - знахостапі матриці: А = , В =

деА,-С-В); А2 ■= -¿•(а*): А3 = -¿(л-).

иП4і_иП (Аі+1/2)ПиГ-(<1/2)Пи"_1

1 м и+ щ ;

+

АХ;

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

Алгоритм розрахунку течії в основному руслі ріки складається з послідовного розрахунку течії на кожній з ділянок розбиття основного русла Коли необхідно обрахувати течію і по притоках, то кожна притока розглядається як нова дшпгка

Обчислене значення параметрів в і (} с базою для визначення розподілу ЗР по руслу ріки. Чисельне інтегрування рівняння транспортування забруднювачів у ріці аналогічне до підходу, «кий вимори-

і ' ' • ' -

стало у другій главі для побудови алгоритму розрахунку переносу домішки у повітряному середовищі. Неявна схема розв'язування рівняння (9) має вигляд

ах ■

Математичні моделі розрахунку гідродинаміки течій і транспортування забруднень у системі русел на базі моделей (8), (9)

п+1

Ьі- + б8(рр+1 =

Аі

Ах

орієнтовані на розв’язання задач прогнозу якості води у р. Дніпро. Для обгрунтування точності розробленої моделі вибрано задачу про моделювання проходження хвилі в горизонтальному каналі при зруйнуванні гребні. Результати розгляду задачі за розробленим алгоритмом порівнюються з розрахунком С. К Годунова і аполітичним розв'язком, зроблено висновок про достатню для інженерної практики точність розрахунків за розробленою моделлю.

У п. 5.5 КП "Акваторія" використана для розв'язування ряду прогнозних задач про перенесення забруднень по руслу ріки Дніпро для різноманітних фізичних ситуацій (залповий викид ЗР, викид від стаціонарних джерел тощо) (рис. 4,г). •

КП "Акваторія" засновано на однозимірній моделі гідродинаміки течії і транспортування домішок у акваторії рік великої прапажосгі з регулюванням стоку на греблях. На практиці часто необхідно здійснити прогноз забруднення поверхневих вод на деякій ділянці водотоку поблизу від джерел викиду ЗР у водойму. В цьому'випадку для адекватного прогнозу якості води в руслі необхідно прк розрахунку гідродинаміки течії і переносу домішки враховували геометрію русла ріки в плані, наявність приток (стоків), островів в руслі та ін. З цією метою у п. 5.6, 5.7 розглянуто математичну модель гідродинаміки течії і переносу ЗР у водоймах, де використовуються рівняння Нааье-Стокса і допущення теорії "мілкої вода":

Для описання процесу транспортування ЗР у водоймі править

і ■ . ’ .

рівняння ;

У п, 5.8, 5.9 побудовано різницевий алгоритм розрахунку гідродинаміки течій у водоймах. Реалізація алгоритму здійснюється на прямокутній рознесеній сітці: функція току визначається у вузлах різницевої сітки, згзіїхрешсть - у центрах, компонента швидкості и -на вертикальних, а компонента швидкості V - на горизонтальних межах різницевих осередків.

Різницева схема розщеплення, що має другий порядок точності за часом, для рівняння переносу зааихреності має вигляд

Для розв'язування рівняння Пуассона застосовується метод встановлення розв'язку за фіктивним часом

(12)

Схема розщеплення дня даного рівняння маг вигляд

За розраховшшші значеннями функції току ц/ у вузлах різницевої сітки визначаються компоненти швидкості на вертикальних і горизонтальних гранях різницевих осередків. Гранична умова типу Неймана задовольняється в ході ітерацій за рахунок визначення завихре-носгі у приграничшїх осередках. Знайдені значення компонент вектора швидкості води використовуються для розрахунку транспортування ЗР у водоймі. У п. 5.10 побудована неявна різницева схема розщеплення рівняння (12), що маг, вигляд:

Чисельний алгоритм розрахунку гідродинаміки гсчії і переносу забруднень у водоймах реалізовано у вигляді КП "Водойма". Для моделювання гідродинаміки і тралсіюргу ЗР у водоймах зі складною

геометрісю берегів, при наявності приток (рис. 4,в), островів в розробленому алгоритмі використовуються маркери, які дозволяють виділити різницеві осередки, що стосуються русла водойми, приток і суші.

Для розрахунку гідродинаміки течі! та транспорту ЗР в руслах, що мають складну геометрію границь при великих числах Не, побудовано КП "Водойиа-2", який допомагає реалізувати на ПЕОМ модель вихрової відривної течії нев'язкої рідіти. (7) (п. 5.10.1)

У п. 5.11 проведено методичні дослідження розробленого чисельного алгоритму для розрахунку гідродинаміки течії у руслах. Методичне дослідження проведено на прикладі задачі про сбтігаЗіня зворотнього уступу плоским потоком нестисливої в'язко! рідіппі, течії в порожшні прямокутного перерізу з верхньою рухомою границею та втікання струменя у прямокутну камеру.. Результати розрахунку за розробленим алгоритмом порівнюються з широко відомими у літературі експериментальними даними і розрахунками інших азторіз (рис. 1,6). ■ _

У п. 5.12-5.14 КП "Водойма", "ВодоЙма-2" застосовувався для

і ■ ' ' ' дослідження в різшіх фізичних постановках задач забруднення води у

водоймах.

У п. 5.12 наводяться результати розрахунку поширення ЗР у акваторії ріки при викиді забруднювачів з острова. Моделюється ступінь забруднення акваторії водойми у випадку зміни місця викиду ЗР на острові. Обраховувалась динаміка змін інтенсивності, розмірів і форми зони антропогенного забруднення води (рис. 4,5).

У п. 5.13 подано результати математичного моделювання гідродинаміки течії і поширення ЗР у водоймі при викиді домішки з

півострова у акваторію. Прогноз стану водойми і еволюції зони забруднення акваторії виконано для двох різних розмірів півострова і положень місця викиду забруднювачів з півострова

У п. 5.14 представлено результати розрахунку гідродинаміки течії і транспорту ЗР у водоймі зі складною геометрією берегів, при наявності притоків. Зроблено прогноз розмірів, форми та інтенсивності зони забруднення основного русла водойми. Розглядалося розв'язання двох прогнозних задач. Перша задача моделювала ситуацію попадання у основне русло водойми забруднювачів з водою з притоки, причому ступінь забруднення води у притоці значно перевершує величину фонового забруднення води у головному потоці. Постановка другої задачі моделювала ситуацію "розбавлення" забрудненої воді} у головному руслі водою з притоки. Результати розрахунків знайдено за допомогою КП "Водойма" та "Водойма-2" (рис. 3,4,в).

Відомо, що на островах, розташованих у руслах рік, можуть знаходитись промислові підприємства, на території яких нерідко збере-гаються усілякі віходи, сировина і таке інше. Під впливом атмосферних опадів і при відсутності на поверхні землі захисного екрану відбувається надходження забруднювачів у грунт, фільтрація їх і виніс у акваторію ріки. У п. 5.15 побудована математична модель і алгоритм розв'язування даної задачі, яка названа спряженою. Покладаємо, що забруднення у грунті острова поширюються за рахунок дифузії. На границі поділу (береговій лінії) забруднення з острова потрапляє в акваторію за рахунок дифузії, а у рідині перенесення забруднень відбувається двома фізичними процесами: дифузією і конвекціпо У п. 5.15 даються результат проінозу розмірів, форми та інтенсивності зони забруднені« акваторії для розглянуті о механізму забруднення

водойми.

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

дойми. В моделі допускається, що всередині виїмки процес поощрения забруднювачів обумовлений дифузією, а на решті області течії -конвекцією і дифузією. Вихідним рівнянням для описання цього процесу поширення ЗР у водному середовищі є рівняння виду (1) без джерельних складових і величини се, . Для концентрації ЗР т дні виїмки ставляться граничні умови даох видів: і?“фр де (р9 - звдше

• „ би . .

значения концентрації ЗР; чи ■— = я0, де до - відоме зяаяешм потоку ЗР. На основі розглянутої моделі створено КП "Виїмка". Чиспоягіі розв'язок задачі спільно з розгаянутиші крайовими умовами • здійснюється за допомогою неявної ЗТС, побудованої у глазі її. У п. 5.16 проведено параметричні 'дослідження вшиту швидкості і напрямку водного потоку і розмірів шгїмки на ступінь забруднення водного середовища. '

В п. 5.17 розглянуто один з аспектів проблеми моделювання забруднення підземних вод: міграція ЗР у підземних водах від водовідстійників. У комп'ютерній моделі розрахунку гідродинаміки течії і переносу ЗР у підземних водах використовуються усереднені по -товщині водоносного шару Ь рівняння: 1

ти розв’язування задачі про такий механізм забруднення ехвяхорй по-

зо

1

де ,ф(х,у) = |ф(х,у,г)сІ2 - концентрація забруднювачів, 2і

1

Н = -|Н(х,у,г)аг, - п'єзометричний напір, г{, г2 - нижня та верхня

позначки водоносного шару.

В побудованому алгоритмі розрахунку гідродинаміки фільтраційної течії і транспорту ЗР у підземних водах, котрий реалізовано у КП "Водовідстійник”, використовуються маркери, які дозволяють виділити різницеві осередки, що належать тій частині розрахункової області, де розташована ріка чи якась інша водойма, джерела ЗР і осередки, що відповідають області фільтраційної течії. Використання маркерів дозволило створити універсальний алгоритм розрахунку, що дає можливість чисельного моделювання для джерел забруднення і водойм, які мають довільну" форму в плані і розташовані різним чином по відношенню один до одного в розрахунковій області.

Для розв’язування рівняння (16) застосовується ЗТС з викори- . станням ідеі встановлення розв'язку за фіктивним часом:

Для числового інтегрування рівняння (17) використовується ЗТС, побудована у главі II.

В п.5.17 КІТ 'водовідстійник" використано для розв'язування прогнозної задачі про забруднення підземних вод від одного і двох водовідстійників, розташованих поблизу водойми.

В шостій піаді розроблені методи використано для розв'язування проблемних задач з природоохоронної тематики для Придніпровського регіону.

В п. 6.1 досліджено динаміку мінералізації ділянки р. Самара і підземіпгх вод навколо двох відстійників шахтних вод на б. Косьміїша та б. Свідовок ВО “Павлоградпугілля”. Внаслідок проведеного розрахунку досліджено еволюцію, форму та розміри зони забруднення підземних вод в регіоні; визначено час, коли засалені підземні вода досягнуть русла р. Самара і розпочнеться процес П мінералізації; одержано оцінку кількості домішки, яка потрапляє в р. Самара завдяки наявності поблизу водовідстійників.

В п. 6.2 розв'язано задачу про забруднення ділянки р. Самара при раптовому викиді шахтних вод з відстійника на б. Свидовок. Досліджено динаміку руху забруднювача вздовж водотоку та зміну інтенсивності забруднення ріки з плином ча<у; визначено час, коли забруднена вода по руслу досягне найближчого населеного пункту -села Вербки. Одержані тут результати розрахунків використано в регіональному проектно-дослідницькому інституті “Дніпроводгосп” для розробки заходів щодо зниження техногенного навантаження від ПО “Панлоградвугілля” на природне середовище.

В п. 6.3 виконано розрахунок зміни рівня мінералізації р. Дніпро на ділянці впадіння в неї р. Самара при прогресуючому рівні засолення р. Самара. Одержано прогнозну оцінку форми, рівня та інтенсивності зони засолення р. Дніпро на проміжку водотоку біля м. Дніпропетровськ. Виявлено підзони з високим рівнем домішки, розташовані поблизу водозаборів Придніпровського району.

В п. 6.4 розв'язано задачі про поширення шкідливих викидів від п' яти відвалів шахт ВО “Павлоградвугілля”. Шляхом обчислювального експерименту при різних напрямах вітру визначено підзони забруднення, виникнення яких обумовлено взаємовпливом викидів від відвалів; одержано прогнозну оцінку рівня забруднення атмосфери над екологічно значимими об'єктами; населеннями пунктами, сільгоспугіддями, лісопарком та ін,

В п. 6.5 розв'язано задачу про забруднення повітряного середовища від відвалу шахти “Героїв космосу” (ВО “Павлоградвугілля”) під дією метефакгорів та з урахуванням рельєфу місцевості, а також геометричної форми відвалу. Одержано оцінку рівня забруднення повітряного середовища за різних метеоумов над с.-г. угіддями, садом та селом Вербки, що розташовані поблизу відвалу. Одержані результати та розроблені методи поширення ЗР у повітряному середовищі використовуються в практиці проектного інституту “Дніпрогіпроцшхт” при аналізі антропогенного забруднення від відвалів та технологічних комплексів шахт Західного Донбасу.

В п. 6.6 представлено розв'язування задачі про поширеній домішок в робочій зоні шприц-машини (дослідний завод НДІ великогабаритних шин) Методом обчислювального експерименту визначено ефективність роботи аспіраційної системи по видаленню ЗР у ва-

падку вибраної схеми розташування місцевих відсосів.

В п. 6.7 розв'язано задачу про поширення пилових викидів в робочій зоні дозаторіп будівельних матеріалів (Дніпропетровських комбінат ТІромбудмотеріали”). Досліджено ефективність роботи аспіраційної мережі у випадку циклічного викиду ЗР і визначено необхідний час роботи відсосів для видалення ЗР з робочої зони.

В додатку приводяться документи про практичне використання результатів дисертаційної роботи.

Основні наукові результати.

Одержані в' дисертаційній роботі результати спрямовані на розв’язанім важливої наукової проблеми - розробки ефехпшшх'Ма-тематичних методів прогнозування забруднення водного та повітряного середовищ при інтенсивних техногенних навантаженнях.

Основний підсумок роботи. ’

1. Розроблено комплексний підхід до математичного описання різноманітних фізігших процесів поширення забруднюючих речовин у водному та повітряному середовищах, який грунтується на базі дискретних балансових співвідношень для коїгтрольного елемента На базі розробленого підхода побудовано ефективні неявні різницеві алгоритми для розрахунку гідродинаміки та переносу ЗР у повітряному та водному середовищах, які відповідають сучасним напрямкам розвитку чисельних методів типу ТУО. Алгоритми мають логічну простоту, економічні, зручні для програмування і вимагають невеликих витрат машинного часу при реалізації на ПЕОМ. Нойий якісний рівень чисельного моделювання грунтується на розробці алгоритмів корекції потоків, які забезпечують підвищений порядок точності, мінімальні дифузійні похибки при збереженні монотонності розподілу

основних параметрів.

2. Розроблено нові ефективні комп'ютерні методи розрахунку гідроміграційних процесів у протічних водоймах:

* ' •

• неявний метод, що грунтується на чисельному інтегруванні

рівнянь Навьє-Стокса (планова задача), ■

• неявний метод, заснований на чисельному інтегруванні рівнянь вихрових відривних течій ненязкої рідини (планова задача);

• неявний і явний методи, що базуються на чисельному інтегруванні рівнянь річкової гідравліки.

ГГерші два методи дозволяють розрахувати гідродинаміку течії, концентрацію ЗР у водоймі, визначити форму і розміри зони техногенного забруднення води з урахуванням геометричної форми русла, наявності приток, водозаборів та островів зі складною геометричною формою, форми вогнища забруднення, типу і характеру викиду. Розроблені методи дозволили виконати комплекс досліджень, пов'язаний з визначенням впливу даних факторів на інтенсивність забруднення водойм. Методи розрахунку, що базуються на рівняннях річкової гідравліки, дають можливість розрахувати гідродинаміку та концентрацію ЗР у водотоках з урахуванням типу викиду, приток, регулювання стоку та ін. На базі даних методів досліджено динаміку забруднення р. Самари при скиді шахтових вод з відстійника на б. Свидовок.

3. Розроблено ефективний алгоритм розрахунку і розв'язано просторову задачу прогнозування забруднення акваторії водойми при поширенні шкідливих речовин з могильника відходів на дні водойми.

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

плані довільну форму, який було застосовано для розв'язування прогнозної задачі засолення підземних вод і р. Самари поблизу водовідстійників ВО ‘ТТавлоградвугшія’’.

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

6. Розроблено ефективний комп'ютерний метод прогнозування якості повітряного середовища всередині робочих зон. Метод дає можливість визначити концентрацію ЗР п робочих зонах, потрібні об’ смя повітря і режими роботи підсосів для забезпечення необхідних санітарно-гігієнічних норм; аналізувати на стадії проекіувеикя аспірзційної мережі ефективність її роботи дня вкбреяої схеш розміщення місцевих відсосів. Метод було використаю для розв’язання широкого кола задач про забруднення робочих зон на підприємствах шинної та хімічної промисловості.

7. Розроблено комплекс методів прогнозування якості

атмосфери: • ‘

• побудовано ефективний комп'ютерний метод розрахунку

гідродинаміки повітряного пстоху та поширення ЗР в атмосфері з урахуванням рельєфу місцевості, забудови, форми джерел забруднення, типу викиду. Метод дозволяє також розв'язувати задачі про провітрювання виробничих приміщень та виробок; .,

• побудовано ефективні комп'ютерні методи розрахунку поширення ЗР в атмосфері в локальному масштабі, а такоік у випадку далекого переносу ЗР.

Розроблені методи прогнозування якості атмосфери було використано для розп'язуноння прикладних задач про забруднення попітря

від джерел гірничопромислового комплексу.

8. Розроблено спеціалізоване програмне забезпечення - КП “Акваторія”, “Водойма”, “Водойма-2”, “Виїмка”, “Водовідстійник” “Прогноз-1”, “Прогноз-2”, “Протоз-3”, “Аспірація" для аналізу і прогнозу якості водного та повітряного середовищ.

9. Розв'язано ряд важливих прикладних задач про поширеній ЗР у повітряному та водному середовищах, що проілюструвало ефективність та універсальність розроблених КП. Результати роботи впроваджено в практиці низки установ та вузів України.

Основні положення і результати'диссертаційної роботи надруковані у 67 роботах, у тому числі

1. Беляев H.H., ХрущВ.К. Численный расчет распространения аэрозольных загрязнений. Учебное пособие. Днепропетровск: ДГУ, 1990, - 80с.

2. Беляев H.H., Хрущ В.К. Расчет принудительной вентиляции карьеров от вредных примесей // Математические методы тепломассопереноса. - Днепропетровск: ДГУ, 1985. - С. 60-63.

3: Беляев Н.Н. Численні, й расчет распространения аэрозоля от нескольких источников загрязнения // Численные решения задач механики жидкости и газа. Днепропетровск: ДГУ, 1988. - С. 73-80.

4. Беляев H.H., ХрущВ.К. Численное решение распространения

аэрозоля при точечном выбросе // Гидроаэромеханика и теория упругости. - Днепропетровск: ДГУ, 1988. - С. 11-15. .

5. Беляев H.H., ХрущВ.К. Расчет распространения выброса аэрозоля в атмосфере // Математическое моделирование процессов тепломассопереноса. - Днепропетровск: ДГУ, 1988. - С. 120-123.

6. Беляев H.H., ХрущВ.К. Численное моделирование удаления примесей из рабочей зоны технологического оборудования // Расчет течений жидкости и газа Днепропетровск: ДГУ, 1989. - С. 145-163.

7. Беляев H.H. Расчет переноса спор растений в атмосфере // Гидроаэромеханика, и теория упругости. - Днепропетровск. ДГУ, 1989. -С. 15-18.

Ü. Беляев H.H., ХрущВ.К. Расчет вентиляции помещения от вредных примесей // Гидроаэромеханика и теория упругости. -Днепропетровск: ДІ’У, 1990. - С. 23-26

9. Беляев H.H., Хрущ В.К. Численное моделирование вентиляции помещения от аэрозольных загрязнителей. // Моделирование и методы расчета процессов тепломассопереноса. - Днепропетровск: ДГУ, 1990. - С. 155-163

10. Беляев H.H. Расчет сноса химикатов ветром при опрыскивании растений // Численно-аналитическое исследование процессов теплообмена. - Днепропетровск: ДГУ, 1990. - С. 41-45.

11. Беляев H.H. Расчет процесса распространения примеси в помещении // Численные методы и математическое моделирование тепломассопереноса - Днепропетровск: ДГУ, 1991. - С. 142-148.

12. Беляев H.H., Побединская Л.В. К расчету распространения

аэрозольных примесей от движущегося источника // Численные методы и математическое моделирование тепломассопереноса. -Днепропетровск: ДГУ, 1991. - С. 95-101. ,

13. Беляев H.H., Побединская Л.В. К расчету распрссхрапекиз

примеси в помещении от движущегося истоштка-зсгрязиителя // Математическое моделирование в механике жидкости и газа. -Днепропетровск: ДГУ, 1992. - С. 103-108. ■

14. Беляев H.H. Математическое моделирование процессов распространения ядохимикатов при обработке сельскохозяйственных объектов и ' его компьютерная реализация // Вестник Дмепропетропского университета Биология и екологая. -Днтропеггровськ, ДДУ, 1993. - С. 28.

15. Беляев H.H. Компьютерное моделирование распространения,

вредных веществ в paßo4ifx зонах технологичесюи агрегатов // Иза. вузов и высших энергетических объединений СНГ. Сер. Энергетика,

1994.-№ 5-6. - С. 70-75. . <

16. Антонова В.А., Беляев H.H., Хрущ В.К. Моделирование уноса загрязнений с углублений на дне водоема // Изв. зузоа и высших энергетических объединений СНГ. Сер. Энергетика, 1994. -Xs 7-8. -С. 104-107.

17. Беляев H.H., Хрущ В К. Математическая модель

нестационарного пространственного переноса загрязнений' в атмосфере // Изв. вузов и высших энергетических объединений СНГ. Сер. Энергетика 1994. -№ 1-2. - С. 134-141. ■ '

18. Беляев H.H. и др. Прогнозирование экологического ущерба при ведении боевых действий неядерными средствами. И Космическая техника. Ракетное вооружение. Научно-технический сборник, вып 3. - Днепропетровск, 1994. - С. 37-49.

19. Беляев H.H. и др. Математическое моделирование гидродинамики ^ счсний и переноса загрязнений п водном бассейне

Днепра II Космическая техника Ракетное вооружение. Научнотехнический сборник, вып. 3. - Днепропетровск, 1994. - С. 50-58.

20. Беляев H.H., Малевич Ю. А., Хрущ В.К. Численное моделирование гидродинамики течения и переноса загрязняющих веществ в мелких водоемах. // Изв. вузов и высших энергетических объединений СНГ. Сер. Энергетика, 1994. - № 11-12. - С. 81-85.

21. Беляев H.H., Хрущ В.К. Численное моделирование удаления примеси местными отсосами II Гидромеханика Киев: Наукова думка,

1994. - Вып. 68. - С. 26-30.

22. Беляев H.H., Малевич Ю.А., Хрущ В.К Моделирование загрязнения подземных вод в окрестности водоотстойника II Изв. вузов и высших энергетических объединений СНГ. Сер. Энергетика,

1995. - № 1-2. - С. 82-85.

23. Беляев М.М. та ін. Комп'ютерне моделювання транспорту

забруднюючих речовин в акваторії рік І/ Хімічна промисловість України. - Київ. - . - 1995. - ЛУЇ. - С. 40-42.

24. Беляев H.H. и др. Компьютерные методы прогноза качества

поверхностных и подземных вод. // Хімічна промисловість України. -KitJ'b. - . . - 1995. - № 3. - C. 77-80.

25. Беляев H.H. Компьютерно-вычислительные системы прогноза загрязнения воздуха в промышленных регионах Украины. // Биологические и технические системы регулирования. -Днепропетровск: ДГУ. - 1995 - С. 6-9.

26. Беляев H.H., Хрущ В.К Компьютерное моделирование рассеивания загрязняющих веществ в атмосфере, водоемах и подземных водах. // Сучасні технології у важкому та транспортному машинобудуванні. Від ділення важкого та транспортного машинобудування Академії інженерних наук України. - Дніпропетровськ. -

1995. -С. 95-104.

' 27. Беляев H.H.и др. Прогноз качества атмосферы при аварий-

ных выбросах II Новые разработки в машиностроительном производстве. - Академия инженерных наук Украины. - Днепропетровск,

1996. г С. 84-88.

28. Беляев Н.Н.и др. Математическое моделирование загрязнения приземного слоя от отвалов // Новые разработки в машиностроительном производстве. - Академия инженерных наук Украины. -Днепропетровск, 1996. - С,81-83.

29. Belyaev N N. et ai. The mathematical simulation of the pollution of the ground waters. // International Congress "Water: Hcoloyy and Technology". Moscow, September 6-9, 1994. - Vol. 1. -P. 90.

30. Belyaev N N. et al. The prediction of the River Dniepr pollution by the method of the mathematical modelling. II (see in previous book). -P. 91.

31. Belyaev N N. et al. The use of the space investigations алгі the results of the mathematical modelling for creation of the computer system for transboundary pollutant transfer calculation. // IAC'94 International Aerospace Congress. Theory, Applications, Technologies. August 15-19, 1994. Moscow, Russia - P. 79.

32. Pobedinskaya L.V., Belyaev N.N. The application of the computer experiment for the calculation of the pollutant dispersion in the industrial areas of the tire plants. // IRC'94, Environment Protection. - Volume 3. -Moscow. - 1994. - P. 652-654.

33. Belyaev N.N., Khnitsch V.K. The joint use of the space probe results and the results of the computer modelling for creation the computer system for calculation of long-range pollutant transfer. // AIAA Paper,

1994. - № 1AA-94-1AA. 33.650. - Washington, USA. - 3P.

34. Belyaev N.N., Khrutsch V.K. Numerical modelling of the flows and conjugate heal transfer under microgravitation in closed areas. // AIAA Paper, 1994. - }h IAF-94-J.5.260. - Washington, USA. - 4P.

35. Belyaev N.N., Khnitsch V.K Numerical modelling of the gydro-dynamics of flows in the case of the regulation of the run off in river beds.

И 1CFM5/95, Proc. 5th Int. Conf. Fluid Mechanics, Cairo, Egypt, 2-5 Jan.

1995. - Cairo, Cairo University Publication. - 1995. - Vol. 1. - P.339-347.

36. Belyaev N.N., Khrutsch V.K. Implicit monotonic difference-schemes of splitting of the second oder of accuracy of solving the Euler and Navier-Stokes equation // 1CFM5/95, Proc. 5th Int. Conf. Fluid Mechanics, Cairo, Egypt, 2-5 Jan. 1995. - Cairo, Cairo University Publication.

- 1995. -Vol. 1.- P. 369-377.

37. Belyaev N.N., Khnitsch V.K. Computer system for analysis and

prediction of air pollution in plant shops. // Доповіді HAH України. -Київ - 1995. -№ 2. -C 130-131. ,

Belyaev N N. Mathematical simulation of the pollutant transfer in the aqueous medium and atmosphere. The manuscript is represented for the Doctoral Degree in the specialty 05.13.02 - mathematical simulation in science investigations Taras Shevchenko Kiev University, Kiev, 199<?. 67 scientific publications are defended which consist of the unity of theoretical investigations on the mathematical simulation of the pollution processes

in water bodies, ground waters, atmosphere, air of the industrial rooms. The new difference schemes to calculate both the hydrodynamics of flows in water bodies and pollutant transfer in the atmosphere, water bodies, ground waters were developed. The software capable of predicting pollution transfer under varying factors was developed. The wide range of the problems of the environment protection solved with the software developed is presented.

Беляев H.H. Математическое моделирование распространения загрязняющих веществ в водной и воздушной средах. Диссертация на соискание ученой степени доктора технических наук . по специальности 05.13.02. - математическое моделирование в научных исследованиях. Киевский университет имени Тараса Шевченко, Киев,

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

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

и

ср-ІО'9 г/м3

10

У, м

• -V

( /

і / И

<! у

Не

РИС. І,а Розподіл їошкїгграшї ЗР: - - розрахуйся тл схемою 1-го порядку точності; — - ротрпхуігог з впхорн-сташіям янтидпфузіПнопо кроку; о-амалгпічшій розв язок; О - експеримент (Р.В.Озмідов).

0 200 400

Рис. 1,6. Залежність від числа Не геличшм області зворотної течії за

уступом: ----- - розрахупої за

розробленим алгоритмом; X -розрахунок (І.О.Бслов та іл.); © -експеріїмеїгг (В.ЛгтаІу)

Рис. 2. Раїполш канцеїпрацції ЧРпісля вибуху п кар'єрі

Ряс. 3. Рочподіл конпентршіії ЗР при внесенні домішки з двома притоками.

(1=0,40); б - при викиду ЗР з острова (12,8); в - внесення ЗР з притокою; г - положення "плями" ЗР в руслі річки для двох моментів ■' • часу. ' ■ . ;