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

Актуальність теми. За своєю будовою земна кора - багатокомпонентна, містить тріщини та пустоти. Її складові відрізняються між собою щільністю, пористістю, мінералогічним складом, здатністю до деформування та ін.

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

У роботах Г.Н. Каменського, Н.К. Гіринського зазначається, що взаємодія горизонтальних водоносних горизонтів суттєво залежить від розмежовуючого їх слабкопроникливого прошарку. Його вплив враховано умовою перетоку.

У роботах І.В. Сергієнка, В.С. Дейнеки, В.В. Скопецького побудовані математичні моделі основних процесів, що протікають у багатокомпонентних грунтових середовищах, для випадків довільного розташування в просторі тонких включень (у вигляді крайових, початково-крайових задач для рівнянь у частинних похідних з умовами спряження та відповідних узагальнених задач).

У 70 - 80-х роках минулого століття були розроблені підходи до побудови ефективних різницевих та скінченно-елементних схем дискретизації основних класів задач математичної фізики з гладкими розв'язками. Значно менше було робіт стосовно побудови таких алгоритмів для задач з умовами спряження (з розривними розв'язками). При побудові розв'язків таких задач використову-валися різноманітні підходи, в тому числі підхід, що базується на продовженні рівняння нерозривності на лінію розриву за допомогою - функції Дірака. Для еліптичного рівняння з розривом розв'язку на прямій, що паралельна одній з осей координат, А.А. Самарський та В.Б. Андрєєв розробили та дослідили різницеву схему. І.М. Молчанов і Л.Д. Ніколенко побудували функціонал енергії для звичайного диференціального рівняння 2-го порядку з умовами спряження. Для з'ясування питання існування та єдиності узагальненого розв'язку використали клас неперервних функцій.

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

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

Зв'язок роботи з науковими програмами, планами, темами.

Дисертаційна робота виконувалася в рамках наукових тем та проектів Інституту кібернетики імені В.М.Глушкова НАН України: “Розробка та теоретичне обгрунтування моделей і високоточних обчислювальних алгоритмів дослідження проблеми екологічного стану грунтів і грунтових об'єктів” (№ 0197U015907); “Розробити математичну теорію, методи та програмні засоби розв'язування окремих класів задач прикладної математики” (№ 0198U005043); “Розробка та обгрунтування проблемно-орієнтованих методів та програмно-алгоритмічних засобів для розв'язання деяких класів задач прикладної математики на ПЕОМ” (№ 0199U001031); “Розробка математичних моделей та алгоритмів дослідження стаціонарних та динамічних процесів дифузії в неоднорідних середовищах” (№ 0102U006942); “Розроблення математичного забезпечення вирішення складних задач аналізу навколишнього середовища” (контракт № ОК_2001_3_ЦЕНТР_ІК від 14.09.2001 р.) (№ 0101U007894); “Розроблення математичного забезпечення аналізу процесів в грунтових середовищах на основі поєднання усномовного діалогу з високоточним моделюванням” (контракт № ОК_2002_5_ЦЕНТР_ІК від 5.07.2002 р.) (№ 0101U006938).

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

Для досягнення сформульованої мети в процесі досліджень поставлені та розв'язані наступні задачі:

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

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

– побудова високоточних обчислювальних алгоритмів методу скінченних елементів (МСЕ) дискретизації початково-крайових задач з розривними розв'язками;

– теоретичне обґрунтування розроблених алгоритмів та встановлення відповідних оцінок похибок наближених розв'язків;

– створення системи НАДРА автоматизованого розрахунку полів у багатокомпонентних середовищах;

– проведення експериментального дослідження розроблених обчислю-вальних алгоритмів та автоматизованої системи НАДРА.

Об'єкт дослідження. Процеси в багатокомпонентних грунтових середовищах.

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

Методи досліджень. Математичне моделювання та чисельні методи, зокрема, метод скінченних елементів та скінченно-різницеві схеми Кранка - Ніколсона. При побудові математичних моделей використані основні закони збереження. Класичні узагальнені задачі побудовано з використанням методу Гальоркіна та відповідних класів розривних функцій. Високоточні обчислювальні алгоритми розроблено на основі МСЕ та відповідних класів розривних функцій. Автоматизована система НАДРА створена з використанням засобів об'єктного програмування та реалізації розроблених високоточних обчислювальних алгоритмів.

Наукова новизна одержаних результатів. Отримані результати є новими та такими, що узагальнюють існуючі. Зокрема:

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

– побудовано класичні узагальнені задачі, визначені на класах розривних функцій;

– з використанням класів розривних функцій МСЕ побудовано високо-точні алгоритми знаходження наближених узагальнених розв'язків;

– доведено існування та єдиність наближених узагальнених розв'язків;

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

– розроблено різницеві схеми Кранка - Ніколсона для дискретизації відпо-відних задач Коші;

– отримано оцінки похибок наближених розв'язків за допомогою функцій МСЕ та схем Кранка - Ніколсона;

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

– створено автоматизовану систему НАДРА для математичного моделювання процесів у багатокомпонентних середовищах;

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

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

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

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

За створення математичних методів та програмного забезпечення для розв'язання деяких класів задач захисту навколишнього середовища у 2002 р. автору була присуджена премія Президента України для молодих вчених Національної академії наук України.

Особистий внесок здобувача. У роботах [1-3, 6, 10], що написані у співавторстві, автору дисертаційної роботи належать усі теоретичні та практичні результати, крім постановки задач та ідей теоретичних досліджень; у [4, 8, 9, 11] - розробка об'єктно-орієнтованої моделі автоматизованого комплексу дослідження процесів у багатокомпонентних середовищах, об'єктно-орієнтованих підсистем вводу інформації, задання фізичних параметрів, розв'язання задачі та виведення результатів у зручному для користувача вигляді, програмна реалізація комплексу в середовищі C++ Builder.

Апробація результатів дисертації. Основні результати дисертаційної роботи обговорювалися та доповідалися на конференціях: восьмій Міжнародній науковій конференції імені академіка М. Кравчука (Київ, 11-14 травня 2000 р.); Міжнародній конференції “Моделювання та оптимізація складних систем” (МОСС - 2001, Київ, 2001 р.); Міжнародній конференції “Обчислювальна та прикладна математика”, присвяченій 80-річчю академіка НАН України І.І. Ляшка (Київ, 9-10 вересня 2002 р.).

Автоматизована система НАДРА демонструвалася: на виставці ювілейної сесії загальних зборів НАН України, присвяченій 50-річчю з часу створення першої в континентальній Європі електронно-обчислювальної машини “МЭСМ” (Київ, 15, 16 січня 2002 р.); V виставці-ярмарку з міжнародною участю Екологія'2002 (Київ, 21-24 травня, система НАДРА відзначена Дипломом); Міждержавній виставці-ярмарку промислових технологій, наукових розробок, засобів виробництва, товарів та послуг країн СНД (Київ, 14-17 травня 2003 р.); VІ міжнародній виставці-ярмарку Екологія - 2003 (Київ, 20-23 травня, система НАДРА відзначена Дипломом).

Результати дисертації систематично доповідалися та обговорювалися на наукових семінарах відділу математичних методів і програмних засобів прикладної інформатики Інституту кібернетики імені В.М. Глушкова НАН України. У повному обсязі дисертаційна робота доповідалася та обговорю-валася у відділі математичних методів і програмних засобів прикладної інформатики Інституту кібернетики імені В.М. Глушкова НАН України та у відділі системного аналізу і керування Інституту космічних досліджень НАН України та НКА України.

Публікації. Основні результати дисертації опубліковано в 11 роботах, з яких 7 надруковано у наукових фахових виданнях, 2 - у збірниках доповідей міжнародних конференцій, 2 препринти (2 роботи опубліковано без співавторів).

Структура та обсяг роботи. Дисертація складається із вступу, чотирьох розділів, висновків, списку використаних джерел, що містить 137 найменувань, та додатків. Загальний обсяг дисертаційної роботи викладено на 180 сторінках друкованого тексту, обсяг основного тексту - 167 сторінок. Робота включає 30 рисунків, 11 таблиць.

ОСНОВНИЙ ЗМІСТ РОБОТИ

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

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

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

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

Лема 1. Якщо узагальнений розв'язок u(x,t) H існує, то він єдиний.

До розгляду введена лінійна множина функцій v(x) з базисом . За допомогою цього базису побудована підмножина функцій вигляду

Наближеним узагальненим розв'язком задачі (1)-(6) називається функція , яказадовольняє рівності (7), (8), де замість u(x,t) та v(x) використані відповідно функції U(x,t) та V(x).

Теорема 1. Нехай u(x,t) - узагальнений розв'язок початково-крайової задачі (1)-(6), а- її наближений узагальнений розв'язок із. Тоді існують такі додатні сталі , щомає місце нерівність

Теорема 2. Якщо класичний розв'язок задачі (1)-(6) досить гладкий на , то для наближеного узагальненого розв'язку справедлива оцінка

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

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

Теорема 3. Нехай u(x,t) - досить гладкий на класичний розв'язок початково-крайової задачі (1)-(6). Тоді існують додатні сталі , такі, що , для похибки наближеного розв'язку , отриманого з викорис-танням множини та різницевої схеми Кранка - Ніколсона, маємо

Теорема 4. Якщо u(x,t) - достатньо гладкий на , , класичний розв'язок початково-крайової задачі (1)-(6), то для похибки,, наближеного розв'язку , отриманого за допомогою функцій МСЕ з класута різницевої схеми Кранка - Ніколсона (14), (15), справедлива оцінка

З використанням функцій розв'язано модельний приклад. Відносна похибка отриманого розв'язку не перевищувала 5%.

У підрозділі 2.2 розглянуто задачу з розривним потоком. Початково-крайова задача визначена рівнянням (1), крайовими і початковою умовами (2), (3), (6) та умовами спряження

[u]=0, x=, t(0,T], (19)

x=, t(0,T]. (20)

Побудовано узагальнену задачу. Наведено теореми, аналогічні теоремам 1-3. При розв'язанні модельного прикладу з використанням функцій відносна похибка не перевищувала 7%.

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

Визначено рівняння (1) з початковою умовою (6), де

На кінцях відрізка [0,l] задано неоднорідні умови Діріхле, а в точці x= t(0,T] - неоднорідні умови спряження

[u]=(t), (21)

Побудовано узагальнену задачу у вигляді (7), (8), доведено лему про єдиність її розв'язку.

Наближений узагальнений розв'язок шукаємо у вигляді (9).

Встановлено оцінки (теореми, аналогічні теоремам 1, 2) для узагальненого розв'язку початково-крайової задачі та її наближеного узагальненого розв'язку.

Отримано відповідну задачу Коші вигляду (12), (13).

Для похибки , , встановлена оцінка, аналогічна (17), та оцінка вигляду (18).

Обмеження (21) можна розглядати як граничний випадок ( ) умови

Показано, що замість початково-крайової задачі (1), (6), (21), (22) (задачі 1) з головною умовою спряження (21) можна розв'язувати задачу (1), (6), (22), (23) (задачу ), де умова спряження (23) - природна. При цьому для розв'язку , знайденого МСЕ, задачі і класичного розв'язку u(x,t) задачі 1 маємо

Із використанням функцій розв'язано модельну задачу 1 та задачу . Значення відносної похибкидля задачі 1 не перевищувало 7.1%, а для задачі - 5.4 при (при).

У підрозділі 2.4 досліджено вісесиметричну задачу дифузії, яка в області

На кінцях відрізка задано змішані крайові умови, а в точці r= - неоднорідні умови спряження неідеального контакту

При t=0 задано початкову умову. Визначено класичний розв'язок. Узагальнена задача полягає у відшуканні функції , яка задовольняє рівності

Доведено єдиність узагальненого розв'язку. Аналогічно тому, як це зроблено для попередніх задач, введено множини , та визначено наближений узагальнений розв'язок. Встановлено оцінки, аналогічні (10), (11).

Отримано задачу Коші вигляду (12), (13), яка розв'язується з використанням різницевої схеми Кранка - Ніколсона (14), (15).

Для похибки, наближеного розв'язку, одержаного за допомогою множини і різницевої схеми Кранка - Ніколсона, отримано оцінки, аналогічні (17), (18).

З використанням функцій розв'язано модельну задачу. При цьому 5%.

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

У підрозділі 3.1 визначено відкриті, зв'язні строго ліпшицеві області , із; ,, - межа області ,, - порожня множина; , , , - бокова поверхня циліндра ,. Введено простори, деякі норми та напівнорми, необхідні для викладок матеріалу цього розділу.

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

Доведено лему про єдиність узагальненого розв'язку задачі (31)-(36).

Аналогічно до одновимірного випадку знаходиться наближений узагальнений розв'язок задачі (31)-(36) у вигляді

Лінійна множина функціймає базис, тобто

Зауваження 1. Якщо при кожному фіксованому на прямолінійних складових частини границі Г функція не є поліномом степеня не вище k, то, використовуючи класи функцій методу скінченних елементів, одержимо множину, що не є підмножиною.

Отримано задачу Коші

Теорема 5. Нехай u(x,t) - узагальнений розв'язок початково-крайової задачі (31)-(36), а U(x,t) - її наближений узагальнений розв'язок із . Тоді існують такі додатні сталі, що має місце нерівність

Якщо для знаходження наближеного узагальненого розв'язку використовується МСЕ, а саме множина функцій , які для всіх неперервні на та є повними поліномами степеня k змінних на кожному елементарному трикутнику скінченно-елементного розбиття областей, показана справедливість наступного твердження.

Теорема 6. Якщо класичний розв'язок u(x,t) початково-крайової задачі (31)-(36) неперервний разом зі своїми частинними похідними , ,(l=1,2) на кожній з областей і має неперервні частинні похідні , , , на (j=1,2), то для наближеного узагальненого розв'язкумає місце оцінка

Схема Кранка - Ніколсона має вигляд

Теорема 7. Нехай u(x,t) узагальнений розв'язок задачі (31)-(36), частинні похідні,, разом із - неперервні та обмежені на(j=1,2). Тоді існують додатні сталі , такі, що , для похибкинаближеного розв'язку , , одержаного з використанням множини та різницевої схеми (45), (46), справедлива нерівність

Теорема 8. Якщо класичний розв'язок u(x,t) початково-крайової задачі (31)-(36) досить гладкий на кожній з областей, , то для наближеного узагальненого розв'язку , отриманого за допомогою множинифункцій МСЕ та різницевої схеми Кранка - Ніколсона (45), (46), має місце оцінка

Розв'язана задача за допомогою МСЕ з використанням кусково-лінійних функцій. Відносна похибка чисельного розв'язку не перевищувала .

Підрозділ 3.3 присвячено задачі з неоднорідною головною умовою спряження. На області визначено рівняння (31), де , початкову умову (36), крайові умови (32), (33) та умови спряження

Побудовано узагальнену задачу.

Наближений узагальнений розв'язок U(x,t) має вигляд (39).

Зауваження 2. Якщо функції (або одна з них) при кожному фіксованому не є поліноміальними на відповідних прямолінійних складових відрізків , , то, використовуючи класи функцій МСЕ, одержуємо множину, що не є підмножиною. У дисертаційній роботі описано процес побудови функції вигляду

Отримано задачу Коші вигляду (41), (42).

Для узагальненого та наближеного узагальненого розв'язків одержано оцінки, аналогічні (43), (44).

Наближений розв'язок задачі Коші знаходиться за схемою (45), (46).

Для похибки також встановлено оцінки, аналогічні (47), (48).

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

У дисертаційній роботі показана можливість врахування головної умови (49) природною умовою з малим параметром . Розглянуто задачу, що задана рівностями (31)-(33), (34), (50) та обмеженням

У підрозділі 3.4 досліджено квазілінійну задачу з умовами власного зосередженого джерела. В області визначено рівняння

На границі задано умову Діріхле.

На умови спряження зосередженого власного джерела мають вигляд

[u]=0, t(0,T], (54)

=const>0, t(0,T]. (55)

Коефіцієнти задовольняють умову Ліпшиця, обмежені та гладкі на , .

При t=0 початкова умова має вигляд (36).

Узагальнений розв'язок задовольняє рівняння (37), (38), де

Введено до розгляду відповідну множину .

Наближений узагальнений розв'язок початково-крайової задачі (53)-(55), (36) задовольняє рівності (37), (38), де замість функцій та записуються відповідно та , і

Задача Коші має вигляд

де елементи матриці визначено у підрозділі 3.2, а елементи матриці та координати векторів , мають вигляд

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

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

Для наближеного розв'язку задачі Коші, знайденого з використанням класута різницевої схеми Кранка - Ніколсона, одержана оцінка, аналогічна встановленій в теоремі 8.

За допомогою запропонованих обчислювальних схем та класу функцій МСЕ розв'язана модельна задача. При цьому 7%.

У четвертому розділі розглянуто розроблений на основі поєднання математичного опису процесів у багатокомпонентних середовищах і об'єктно-орієнтованого програмування автоматизований комплекс НАДРА.

Комплекс НАДРА складається із системної компоненти, інтерфейсної складової та наступних підсистем: вводу геометричної інформації (ПВГІ); створення упорядкованих чотирикутних зон (ПЧЗ); скінченно-елементного розбиття (ПСЕР); вводу фізичних характеристик, характеристик зовнішнього впливу та тонких включень (ПВФП); формування систем алгебраїчних рівнянь методу скінченних елементів (ПФСАР); розв'язання систем алгебраїчних рівнянь методу скінченних елементів (ПРАЗ); виводу результату розв'язання системи алгебраїчних рівнянь методу скінченних елементів (ПВР).

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

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

У підрозділі 4.13 наведена математична модель усталеного руху рідини в насичено-ненасичених грунтових зонах із смітника до р. Дніпро. За допомогою комплексу НАДРА проведено розрахунки сформульованої задачі та наведено їх результати. Отримані результати використано при оцінці впливу очисних споруд на навколишнє середовище, що підтверджено відповідним актом про впровадження.

ВИСНОВКИ

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

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

2. Побудовані класичні узагальнені задачі, визначені на класах розривних функцій. Доведено єдиність узагальнених розв'язків таких задач.

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

4. Розроблена методика заміни головної неоднорідної умови спряження природною з малим параметром .

5. Встановлено оцінки похибок збуреного розв'язку та наближеного збуреного розв'язку.

6. Для дискретизації відповідних задач Коші розроблено різницеві схеми Кранка - Ніколсона. За допомогою них та функцій МСЕ отримано оцінки похибок наближених розв'язків.

7. Побудовано об'єктну модель, на основі якої розроблено підсистеми вводу геометричної інформації про розв'язувану задачу та інформації про взаємодію об'єкта з оточуючим середовищем. Створено автоматизовану систему НАДРА для моделювання процесів в багатокомпонентних середовищах.

8. З використанням розробленої системи проведено чисельні експерименти та розв'язана практична задача.

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

Дейнека В.С., Калинюк Н.А. Чисельна дискретизація параболічних рівнянь з умовами спряження та коефіцієнтами залежними від часової змінної //Волинський математичний вісник. - 2000. - № 6. - С. 48-54.

Калынюк Н.А., Дейнека В.С. Численная дискретизация цилиндрически-симметричной задачи диффузии с неоднородными условиями сопряжения //Компьютерная математика. - Киев: Ин-т кибернетики им. В.М. Глушкова НАН Украины, 2000. - С. 104-113.

Калынюк Н.А., Дейнека В.С. Численная дискретизация начально-краевой задачи для параболического уравнения с неоднородными главными условиями сопряжения //Компьютерная математика. - Киев: Ин-т кибернетики им. В.М. Глушкова НАН Украины, 2001. - С. 103-112.

Автоматизированная система НАДРА для исследования процессов в многокомпонентных средах / И.В. Сергиенко, В.С. Дейнека, Н.А. Калынюк, М.М. Радионов //Компьютерная математика. - Киев: Ин-т кибернетики им. В.М. Глушкова НАН Украины, 2002. - Вып. 1. - С. 3-11.

Калынюк Н.А. Подсистема ввода и обработки геометрической информации в комплексе НАДРА //Компьютерная математика. - Киев: Ин-т кибернетики им. В.М. Глушкова НАН Украины, 2002. - Вып. 1. - С. 12-23.

Калынюк Н.А., Дейнека В.С. Высокоточные вычислительные алгоритмы для параболического уравнения с условиями сопряжения и коэффициентами, зависящими от времени //Компьютерная математика. - Киев: Ин-т кибернетики им. В.М. Глушкова НАН Украины, 2002. - Вып. 2. - C. 31-40.

Калынюк Н.А. Численное решение нелинейного параболического уравнения с условиями сосредоточенного собственного источника //Компьютерная математика. - Киев: Ин-т кибернетики им. В.М. Глушкова НАН Украины, 2002. - Вып. 2. - C. 106-113.

Автоматизированный комплекс НАДРА для исследования процессов в многокомпонентных средах / И.В. Сергиенко, В.С. Дейнека, Н.А. Калынюк, М.М. Радионов : В 2 ч. - Киев, 2002. - Ч. 1. - 36 с. - (Препр. /НАН Украины. Ин-т кибернетики; 2002 - 3).

Автоматизированный комплекс НАДРА для исследования процессов в многокомпонентных средах / И.В. Сергиенко, В.С. Дейнека, Н.А. Калынюк, М.М. Радионов : В 2 ч. - Киев, 2002. - Ч. 2. - 36 с. - (Препр./ НАН Украины. Ин-т кибернетики; 2002 - 4).

Deyneka V.S., Kalinjuk N.A. The Higher Order Precision Schemes for Parabolic Equations with the Nonideal Contact Conjugation Conditions and Depending on Time Coefficients // VIII Міжнародна наукова конференція імені академіка М. Кравчука (Київ, травень, 2000). Матеріали конференції. - Київ, 2000. - С. 266.

Система НАДРА для дослідження процесів в багатокомпонентних середовищах / І.В. Сергієнко, В.С. Дейнека, Н.А. Калинюк, М.М. Радіонов //Збірка тез міжнар. конф. “Обчислювальна та прикладна математика”, присвяченої 80-річчю академіка НАН України І.І. Ляшка (Київ, 09.09-10.09. 2002). - Київ, 2002. - С. 90.

Размещено на Allbest.ru