ПРОЦЕССЫ И АППАРАТЫ ХИМИЧЕСКОЙ ТЕХНОЛОГИИ
УДК 621.51-251
Аль Аззави Аус, М. Ю. Перухин, Т. Н. Мустафин
ПОСТРОЕНИЕ СЕТКИ РОТОРОВ ВИНТОВОГО КОМПРЕССОРА
Ключевые слова: винтовой компрессор, построение сетки.
В статье приведено описание построения сетки роторов винтового компрессора, а так же рассмотрены аспекты определения метрических коэффициентов для дальнейшего математического моделирования. Проанализированы методы решения систем линейных уравнений на ЭВМ и выбран наиболее оптимальный.
Keywords: screw compressor, building grids.
This article describes the Grid generation of screw compressor and also deals with aspects for the definitions of the metric coefficients for advanced mathematical modeling. The methods of solving systems of linear equations on the computer and choose the best
Одним из первых шагов в процессе численных решений дифференциальных уравнений с частными производными, описывающих физический процесс, является построение сетки. Сеткой называется процесс замены пространственной области с помощью системы узловых точек [1]. Физическую область необходимо покрыть сеткой, так чтобы обозначить все дискретные объемы или элементы, где могут быть применены законы сохранения. Хорошо построенная сетка значительно улучшает качество расчетов, и, наоборот, плохо построенная сетка является одной из основных причин плохого результата. Во многих приложениях, трудности с численным моделированием можно отнести к плохому качеству сетки. Существует два метода создания сетки: структурированный и неструктурированный подходы [2].
Для задачи построения сетки вала винтового компрессора был применен структурированный метод построения. Однако стоит заметить, что построение таким методом сетки в 3-Б связано с большими трудностями. Существует три основных математических метода построения сетки:
1. методкомплексной переменной;
2. алгебраический метод;
3. метод, связанный с решением дифференциальных уравнений частных производных.
Метод комплексной переменной имеет преимущество, так как используются аналитические или частично аналитические преобразования, в отличие от других методов, которые полностью численные. К сожалению, данный метод не может быть использован в трехмерном пространстве [2]. В отличие от метода комплексной переменной, алгебраические и дифференциальные методы могут быть использованы в сложных трехмерных (3-Б) задачах.
Алгебраическая сетка часто основывается на трансфинитной интерполяции. Она может представлять собой как процедуру многомерной интерполяции, так и логическую сумму одномерных интерполяций вдоль каждой вычислительной координаты. Метод широко используется для построения расчетных сеток в областях с гладкими границами, которые не сильно деформированы или в качестве на-
чального приближения для итерационного процесса построения эллиптической сетки. Наиболее подходящий метод численного построения сетки винта -трансфинитная интерполяция, в сочетании с процедурой адаптации и граничной ортогонализации.
После того как определены граничные условия для физического объекта, они могут быть пересчитаны и для координат сетки в расчетной (компьютерной) системе координат [1].
На рис.1 изображены исходные профили ведомого (а) и ведущего (б) валов.
г\ о
1 / (
) (
о О
б
Рис. 1 - Исходные профили ведомого (а) и ведущего (б) валов
Для построения положения сетки данных валов координаты профиля г1 и координаты окружности г2 были представлены в виде уравнения сплайн аппроксимации координат. В качестве аргумента использовалась длина кривой. В этом случае сетка строится постоянным шагом по координатам x и у и вращением по оси г.
Физические координаты можно написать в виде:
яЦ.ч. О = + ее И? *0 ё!П» *С
где в, С, ^-компьютерные координаты; /г-длина ротора; хс,ус - координаты двухмерной сетки, ограниченной профилем ротора, которые могут быть восстановлены из компьютерных координат по формуле:
а
[Оу 0*х Ох О^у Ох а*у Оу Лс1
' Д£ дп' ?КУ
>3ш
££,£Ги
'^у ф Ау
£у
[£*у Я йу1 ГЗЧ « 3у1
I
а* у
/ ах & & ал _ ох ^ ~ Дар"" ~ ар""ер/ " ^' ¡кг Гй'к
■Г
[а* к а/ а/
где <%, С;, Ид:,, а^ - функции перехода Эрмита, которые включают производные в конечных точках кривой для обеспечения ортогональности на границах; А Л (¡Е»Ф - коэффициенты (положительные числа меньше единицы), которые определяются путем проб и ошибок; X! - координата внутренней окружности.
Уравнение (1), после ряда преобразований, позволит определить необходимые метрические коэффициенты, которые, в свою очередь, позволят перевести исходные дифференциальные уравнения, записанные в физической системе координат, в новую компьютерную. Метрические коэффициенты имеют вид:
р - Л Ш\-г т Ш
[аЕу а/ ау1 Га1 у а/
гак ^ ^ аэт ц й щ 'щ, а? \ гау аЕк ас а£у ау аЕк а* а'у]
а/ га* ау &у асп 1ёГ аГ"а?" Ю-
аз щ щ
[ау а^я ах а1 у ах а£у ау аЕх1 а?'а^^а£" а^"~аГа^а?
оу ох а*у оу а'х ох а^у!
^[а^1' а^ 4 а^ ' а^ ер' ща^ ' а-^114
я? ^ " Д} а?" А» \
[Оу аЕх & ^ Ёс рву]
4 щ ' ач ' №}14
+
а/ ву Ау Ал
[ау а£я & а'у ах а£у ау а£к1 а? ' атр"4 "Щ' ~ а^' а^'~ ' а^|+
Я/ гЖх Ау Ау
1' ар4'
рак ау 5Г щ'
а*я ох а1 у ох а'у ау о*х]
а/ так ау ау ^"¡цщ 'ар'
где / - якобиан матрицы, который определяется по формуле:
1_
/— ^ П I Е XI Г"
г? ^ в^ г
По полученным формулам метрические коэффициенты определяются в каждой точке. Данный расчет связан с большим количеством одинаковых итераций. Такие расчеты при последовательном их решении могут занимать много времени. Решением этой проблемы может стать использование параллельного метода обработки данных. В этом случае задача делится на сотни и тысячи частей и решается параллельно. Для увеличения скорости расчетов также применяются ГПУ вычисления.
Для выбора наиболее оптимального подхода расчета коэффициентов на ЭВМ было использовано три метода расчета линейных алгебраических уравнений:
1) Линейный, в котором итерации происходят последовательно;
2) Параллельный, в котором несколько расчетов проходят параллельно;
3) CUDA - это параллельный метод с использованием ГПУ.
Программа написана на языке С++, и представляет собой Windows - оболочку. Изначально программа запрашивает все начальные параметры расчета и количество рассчитываемых точек. Затем производит расчет и выводит на экране 3D сетку.
Для определения оптимального метода расчета для данной определенной задачи были смоделированы все три метода и в дальнейшем определено время их работы [3-4]. На рисунке 2 представлены
100 200 300 400 Число задаваемых точек
500
-linear
-parallel
а
Cuda
8000
U 5000
s
1 4000
a —
e 2000
/
/
-» ■ -»■
зависимости времени от количества расчетных точек для всех методов для ведомого и ведущего валов. Из рисунка видно, что наиболее оптимальным по времени работы программы является классический параллельный метод.
При осуществлении расчета линейным методом задействуется только одно ядро многоядерного процессора, и программа занимает большее время, так как рассчитывается каждая точка последовательно. Параллельный метод вычисляет сразу несколько потоков в одно и то же время, такой расчет будет наоборот требовать не одно ядро, но займет меньше времени. Это наиболее эффективный метод, однако, в нашем случае он показал средний результат. Такой результат объясняется тем, что операция копирования и вставки информации в памяти и файлы занимает больше времени [3]. При более загруженном расчете, скорее всего он показал бы лучший результат.
Подытожив результаты можно сказать, что в процессе описанной в статье работы были смоделированы расчетные сетки роторов винтового компрессора, а так же определены все метрические коэффициенты для дальнейшего решения различных задач. Также к решению данной проблемы были применены различные методы решения систем линейных уравнений на ЭВМ и выбран наиболее оптимальный [3].
Литература
1. Three Dimensional Computational Fluid Dynamics and Solid Fluid Interaction Screw Compressors - Nikola Stosic, Ian Smith, Ahmed Kovacevic. Springer-Verlag Berlin Heidelberg New York.
2. Computational fluid mechanics and head transfer. John C. Tannehill, Dale A. Anderson, Richard H. Pletcher. Taylor & Francis, USA, l997.
3. Аль Аззави Аус, Перухин М.Ю., Мустафин Т.Н. Анализ вычислительных возможностей различных алгоритмов применительно к построению расчетных сеток. Вестник Казанского технологического университета. 20l4. Т. l7. №8. С. 306-307.
4. Мустафин Т.Н., Якупов Р.Р., Акшинская В.В., Хами-дуллин М.С., Хисамеев И.Г. Геометрический анализ зацепления роторов винтового компрессора. вестник Казанского технологического университета. 20l3. Т. l6. №l9. С. 273-277.
100 200 300 400 Число задаваемыхточек
500
-linear -♦—parallel б
-cuda
Рис. 2 - Графики зависимости времени работы программы от количества расчетных точек для ведомого (а) и ведущего (б) валов линейный метод, -0- параллельный метод, -■- метод CUDA)
© Аль Аззави Аус - магистрант каф. АССОИ КНИТУ, [email protected]; М. Ю. Перухин - доцент той же кафедры, [email protected]; Т. Н. Мустафин - к.т.н. каф. холодильной техники КНИТУ, [email protected].
© Al Azzawi Aws - the masters degree candidate KNRTU, [email protected]; M. Yu. Perukhin - associate professor, KNRTU, [email protected]; T. N. Mustafin - associate professor, KNRTU, [email protected].
l65