Page 119 - 140

i 1 i
R   R , (2.168)

i 1 i i 1 i
де R  A н ( K  ) K    R . (2.169)

i 1
Із рівняння (2.169) знаходимо вектор K  , далі обчис-
люємо уточнений вектор концентрацій

i 1 i i 1
K   K   K  . (2.170)

Як початкове значення вектора концентрацій використо-
вуємо розв'язок відповідного лінійного рівняння турбулентної
дифузії для того ж моменту часу змішування.
Розроблена нами програма дає змогу одержати несимет-
ричну криву розподілу концентрацій по довжині зони суміші з
врахуванням різниці в'язкостей транспортованих нафтопроду-
ктів для будь-якого моменту часу їх взаємного змішування у
трубопроводі. Програма передбачає виконання таких обчис-
лювальних операцій.
Задається крок зміни часу змішування  , знаходиться
довжина "голови" суміші і кількість скінчених елементів n . З
e
використанням вектора початкового розподілу концентрацій
по довжині зони суміші K знаходиться розв'язок лінійного
o
лін
рівняння турбулентної дифузії K  . Обчислюється значення
ефективного коефіцієнта дифузії за в'язкістю продукту А D .
a
Виконується автоматичне формування і запис квадратних ма-
триць A і B формату n  n  1.
e
н
н
Знаходиться розв'язок нелінійного рівняння методом
Ньютона-Рафсона. На початку апроксимації приймають

i лін
K  K  . (2.171)


- 113 -