остійна величина. Більше складним реалізація граничної умови 3 роду, того что в цьом випадка Q rp є функцією різниці напорів у Розглянуто вузлі ї на границі. Для Завдання цієї граничної умови в загально віді у вузол уводиться додатковий фільтраційній Опір Ф гр, что розраховується залежних від конкретної ситуации по різніх методиках.
Если алгебраїчна Сума перерахованого вищє доданків балансового Рівняння дорівнює нулю, то зміна уровня у вузлі НЕ відбувається, тобто фільтраційній процес є стаціонарнім. У ІНШОМУ випадка доданок Q ємк характерізує величину ї ШВИДКІСТЬ Зміни рівнянь поверхні:
dXi * dYj t + dt t
Про ємк=m * * (Hij -Hij),
Dt
де: - крок за годиною; - гравітаційна або пружньою водовіддача кулі.
Запішемо балансова Рівняння для ij-того Вузли, вирази его СКЛАДОВІ через гідродінамічні параметри:
dX i * dY j t + dt t H il, j - H ij Hi + lj-Hij
m * * (Hij -Hij)=+ +
dt ФXi-l, j OXi ^ lj
Hij - 1-Hlj Hij + 1-Hij Hrp-Hij
+ ------------- + ------------ + ---------- +
OXij - 1 OYij + l ФГР
+ DХ (We, i + Wt, і + Wісп, i)
або после нескладних перетвореності:
dXi * dYi t + dt t Hi-l, j Ні + l, j Hi, j - 1
m=- * (Hi j - Hi j)=[------- - + -------- + -------- +
dt ФXi + l, j ФXi-l, j ФYi, j - 1
Hi, j + 1 Hrp 1
------- + ------- + dXi * (We, i + Wt, i + Wісп, i)] - Hi, j * (----------- +
ФYi, j + 1 ФГР ФXi-l, y
1. 1 січня l
+ + - + +)
ФXi + l, j ФXi, j - 1 ФYi, j + 1 ФГР
Позначівші в последнего виразі водоемінсть блоку через Ci, j доданок у квадратних дужках через CHi, j ї Последний доданок у круглих дужках через HXi, j, маємо:
+ dt t
С j * (Hy - Ну)=С * Ні, j - Hi, j.
Записати аналогічні Рівняння для всіх розрахункових вузлів сітки, одержимо систему рівнянь, для решение якої можна вікорістаті Різні обчислювальні схеми. ! Застосування класичності явної або неявної схем по різніх причини требует значний витрат машинного годині, тому доцільно вікорістаті набліжені ітератівні методи. Зміст їх Полягає в тому, что невідомі значення напору у вузлах обчислюють Із завданні точністю путем послідовніх набліжень. Використання ітераційніх методів віправдано такоже тім, что фільтраційній процес є нелінійнім у зоне Зміни УГВ. У сітковіх рівняннях це требует багаторазове підбору потужності потоку у верхніх блоках моделі, віходячі з вимог водного балансу. Ця задача вірішується Ефективний ітератівнім путем.
Чи не вдаючися у доповідних аналіз існуючіх чисельного методів решение систем алгебраїчніх рівнянь, відзначімо, что нами звертаючись РОЗРАХУНКОВА схема компонентної релаксації. У ее основу покладаючи ітераційна схема Зейделя, модіфікована введенням пріскорювального (релаксаційного) параметру.
У прийнятя нами Позначення на шкірному кроці ітерації Зейделя величина
t + dt, j обчіслюється в такий способ:
t + dt
- t + dt -Cij-C * Hij
Hi, j=- - для нестаціонарної задачі
Cij + H * Xi, j
C * Hij
Hi, j=для стаціонарної задачі
H * Xi, j
потім вікорістовується релаксаційній параметр:
= t + dt про t + dt =t + dt про t + dt
Ні, j=Hij + Топт. (Ні, j -Hi, j)
де:
Топт - релаксаційній параметр (1 е Т е 2);
про - індекс попередньої ітерації;
=- індекс поточної ітерації;
індекс проміжного значення напору, что обчіслює за схему Зейделя. Далі здійснюється Перехід до наступної ітерації. Обчислення трівають до виконан умови: dHi, je ерs2, де dHi, j - Різниця напору в іj-тому вузлі, отриманий на Попередній и поточній ітерації; eps2 - задана точність решение.
Для верхніх розрахункових блоків, у якіх перебуває УГВ, вместо схеми по компонентної релаксації вікорістається схема простої ітерації. Алгоритм розрахунку Виглядає в такий способ:
- Задається Початкове набліження УГВ....
