Page 112 - 140

рація домішки у кінці n -ого елемента приймається рівною
e
нулю
K  . (2.139)
0
n 

Остання умова вимагає наступної модифікації матриці A і
вектора C
0
a  0 , a  0 , c  . (2.140)
( n n ) 1 ( n n ) 1 n

Для подальшого знаходження вектора вузлових невідомих
у певний момент часу необхідно розв'язати систему n алгебра-
їчних рівнянь з n невідомими. Для цього застосовуємо най-
більш поширений метод – метод виключення (метод Гауса).
Нами розроблена оригінальна програма, яка дає можли-
вість одержати несиметричну криву розподілу концентрацій
нафтопродуктів по довжині зони суміші для будь-якого моме-
нту часу їх змішування. Програма передбачає виконання на-
ступних операцій.
Обчислюється ефективний коефіцієнт дифузії. Залежно
від часу змішування знаходиться довжина "голови" суміші
L і визначається кількість скінчених елементів n . Обчис-
гс
e
люються значення коефіцієнтів a 1 a , 2 b , 1 b , 2 для формування
матриць A і B. Здійснюється автоматичне формування мат-
риць A і B, а також вектора-стовпчика початкового розподі-
лу концентрацій K . Виконується множення матриці B на
o
вектор K , одержується вектор-стовпчик C , здійснюється
o
модифікація вектора C і матриці A для врахування гранич-
них умов

K  1  5 , 0 , K n   0 . (2.141)

Розв'язується система алгебраїчних рівнянь методом Гау-
са, в результаті знаходиться вектор вузлових невідомих K .
Далі розрахунки багатократно повторюються залежно від ча-
су змішування рідин і вибраного приросту часу  .
- 106 -