ння вихідної інформації, що характеризує склад і властивості компонентів суміші.
У блоці 2 розраховуються індивідуальні коефіцієнти Fi (F (I)) по рівнянню (1.22), bi (В (I)) по (1.20), а також сумарний коефіцієнт bm (ВМ) по (1.19).
У блоках 4 і 5, які структурно входять в цикли по I і по J, розраховуються знаменник і чисельник рівняння (1.21), а в блоці 8 завершується його рішення визначенням коефіцієнта Fm (FM), що враховує параметри бінарних взаємодій компонентів суміші.
Блоки 3, 6 і 7 організовують розрахункові цикли по I і J по числу компонентів суміші від 1 до N.
Висновок і передачу розрахованих величин у відповідні підпрограми і процедури здійснює блок 9.
На підставі наведеної блок-схеми розроблена і реалізована підпрограма MIXTURE.
Для розрахунку тиску суміші [див. рівняння (1.1)] за відомою щільності і температурі може бути використана блок-схема (рис. 1.35).
Рис. 1.35. Блок-схема алгоритму розрахунку тиску
У блоці 1 вводяться значення щільності, температури, універсальна газова постійна, константи рівняння Редліха-Квонг і коефіцієнти, попередньо розраховані в підпрограмі MIXTURE.
У блоці 2 по рівнянню (1.18) визначається значення коефіцієнта стисливості суміші, а в блоці 3 вирішується рівняння стану реального речовини щодо тиску.
Обчислення щільності з використанням рівняння (1.18) проводиться відповідно до блок-схемою (рис. 1.36).
Рис. 1.36. Блок-схема алгоритму розрахунку щільності суміші
Рівняння Редліха-Квонг може бути вирішено в явному вигляді тільки щодо коефіцієнта стисливості і тиску. Тому для визначення щільності застосовується один з ітераційних методів, зокрема метод Ньютона, яким вирішується рівняння
. (1.23)
Для локалізації коренів рівняння стану приймаємо первинні значення щільності пари і рідини, які ідентифікуються цілочисельний змінної М: для рідини М=0, для пара М=1 (див. рис. 1.36, блок 2).
У блоках 4 і 6 забезпечується реалізація методу Ньютона, для чого при заданих можливих значеннях щільності при постійній температурі за допомогою рівняння (1.23) знаходимо приватну похідну, уточнене значення щільності і ступінь збіжності DELTA рішення рівняння (1.23).
Після порівняння (у блоці 7) DELTA із заданою точністю рішення EPS приймаємо рішення або про звернення до блок 3 для уточнення вихідного значення щільності і повторення розрахунків, або, при досягненні необхідної збіжності, про передачу інформації в блок 8 виведення інформації, на ніж рішення задачі закінчується.
Література
1.Архаров AM, Марфенін І.В., Микулин Є.І. Кріогенні системи. Т. 1. Основи теорії та розрахунку.- М .: Машинобудування, 2006. - 576 с.
. Теплофізичні якості неону, аргону, криптону і ксенону/Под ред. В.А. Рабиновича.- М .: Изд-во стандартів, 1976. - 636 с.
. Термодинамічні властивості азоту В.В.. Сичов, А.А. Вассерман, А.Д. Козлов та ін. - М .: Изд-во стандартів, 2007. - 352 с.
. Термодинамічні властивості повітря В.В.. Сичов, А.А. Вассерман, А.Д. Козлов та ін. - М .: Изд-во стандартів, +2008.- 346 с.
. Термодинамічні властивості гелію - 4 В.В.. Сичов, А.А. Вассерман, А.Д. Козлов та ін. - М .: Изд-во стандартів, +2004.- 307 с.
. Термодинамічні властивості кисню В.В.. Сичов, А.А. Вассерман, А.Д. Козлов та ін. - М .: Изд-во стандартів, +2001.- 312 с.
. Теплофізичні якості кріопродуктов Л.А.. Акулов, Є.І. Борзенко, В.Н. Новотельная, А.В. Зайцев.- СПб .: Політехніка, 2001. - 243 с.
. Довідник по фізико-технічних основах кріогеникі /м.п. Малков, І.Б. Данилов, Л.Г. Зельдович, А.Б. Фрадков; Под ред. М.П. Малкова.- 3-е изд., Перераб. і доп.- М .: Вища школа, 2005. - 432 с.
