розраховує значення концентрації в кожному вузлі розбиття області рішення на кінцеві елементи.
3.2 Використання ANSYS для отримання проміжних даних
Для визначення поля швидкостей в моделі скористаємося програмним продуктом ANSYS. Будемо вважати, що рух рідини не залежить від механічних домішок, тому можна вирішити задачу про рух рідини без домішок, використовуючи ANSYS. Потім, скориставшись отриманим розбиттям на кінцеві елементи і полем швидкостей, вирішимо завдання про розподіл механічних домішок. При цьому вхідними даними для ANSYS будуть геометричний опис області рішення, початкові і граничні умови. Вихідними параметрами будуть координати вузлових точок кінцевих елементів, вузлове опис елементів і значення вектора швидкості у вузлах. Вихідні параметри обробляються за допомогою програмного продукту, створеного спеціально для цієї мети, при цьому в програмному продукті відбувається рішення методом кінцевих елементів моделі розподілу механічних домішок. Вихідні параметри при передачі даних з ANSYS в програмний продукт мають вигляд, представлений в таблицях - координати вузлових точок, значення компонент вектора швидкості у вузлах, опис кінцевих елементів. br/>
Таблиця 3. Вид даних, що описують координати вузлів розбиття
Таблиця 4. Вид даних, що описують значення вектора швидкості у вузлах розбиття
Таблиця 5. Вид даних, що описують зв'язок номерів елементів з номерами вузлів розбиття
3.3 Обробка даних розрахунку за допомогою розробленого програмного продукту
При використанні ANSYS отримуємо в якості вихідних параметрів координати вузлових точок, опис кінцевих елементів вузловими точками і значення вектора швидкості в кожній вузловій точці. Для вирішення моделі був створений програмний продукт. У програмному продукті на основі даних, розрахованих в ANSYS, будується система лінійних рівнянь. br/>В
Рис. 5. Блок-схема ПЗ
Рішення проводиться найбільш оптимальним алгоритмом для обробки розріджених симетричних матриць. p align="justify"> Для зберігання матриць реалізований алгоритм їх компактного зберігання для мінімізації використовуваної програмою оперативної пам'яті. Оскільки матриця A (2.9) є розрідженій і симетричною, вона зберігається не повністю а тільки верхній трикутник і головна діагональ. Нульові елементи теж не зберігаються, що дозволяє істотно скоротити витрати пам'яті на зберігання даних при значних розмірах матриці A. Пошук потрібного елемента проводиться не шляхом циклічного перебору, а бінарним пошуком за номерами стовпців і рядків. p align="justify"> Вхідні дані:
файл із описом координат вузлових точок.
файл із описом кінцевих елементів вузловими точками.
