Численное и физическое моделирование течения жидкости в вихревом теплогенераторе

МАШИНОСТРОЕНИЕ

УДК 621.1

Ю. М. Ахметов, Р. Р. Калимуллин, В. А. Целищев

ЧИСЛЕННОЕ И ФИЗИЧЕСКОЕ МОДЕЛИРОВАНИЕ ТЕЧЕНИЯ ЖИДКОСТИ В ВИХРЕВОМ ТЕПЛОГЕНЕРАТОРЕ

Рассматривается актуальная проблема численного моделирования многофазного вихревого течения жидкости. Представлены обобщенные результаты работ по моделированию процессов в вихревых теплогенераторах. Приведены результаты численного моделирования внутренних гидравлических и термодинамических процессов. Численное моделирование; вихревой теплогенератор; энергетика; кавитация

ВВЕДЕНИЕ

Основным направлением развития энергетики страны является создание новых установок эффективного преобразования энергии в энергию тепла с малыми потерями и экологически безопасным влиянием на окружающую среду. Одним из многих существующих методов преобразования энергии является вихревой метод.

Обращаясь к трудам различных авторов, изучавших вихревой эффект, можно выявить множество разногласий как в теоретических, так и в экспериментальных исследованиях. На данный момент среди ученых нет единого мнения

о природе возникновения температурного повышения в вихревом теплогенераторе. Сложность изучения данного явления связана с видом движения жидкости в теплогенераторе, так как закрученный поток относится к группе пространственных течений и в поле массовых и центробежных сил происходит непрерывное изменение структуры потока. Такое положение вещей является причиной усложнения механизма протекающих в закрученном потоке процессов и трудностей выявления закономерностей, управляющих этими процессами.

Полученные результаты НИР будут применены для разработки новых нетрадиционных источников энергии, в том числе более эффективных видов вихревых теплогенераторов.

1. СОСТОЯНИЕ ВОПРОСА

В настоящее время все большее значение приобретает использование вихревых теплогенераторов. На сегодняшний день уже существует множество видов теплогенераторов, работающих на разных принципах подогрева [1].

С начала 80-х годов прошлого столетия одним из направлений кафедры прикладной гид-

ромеханики является исследование гидродинамических и тепловых процессов нестационарного течения несжимаемых жидкостей с целью разработки высокоэффективных принципов использования преобразования энергии. В рамках реализации планов исследований учебного научного инновационного центра «Гидропневмоавтоматика» и научно-образовательного центра «Технологии использования новых и возобновляемых источников энергии» создан стенд «Гидродинамическое моделирование высокоскоростного многофазного течения жидкости» (рис. 1). Основные принципиальные решения были отработаны на разработанном прототипе вихревого теплогенератора. Стенд предназначен для проведения фундаментальных и прикладных исследований, научно-технических и учебно-лабораторных работ по изучению нестационарных процессов течения высоконапорных гидродинамических течений несжимаемой жидкости однофазной и многофазной структуры, а также вихревых и кавитационных испытаний центробежных насосов.

Рис.1. Экспериментальный стенд

Контактная информация: (347) 273-09-44

2. ПОСТАНОВКА ЗАДАЧИ

Настоящая работа посвящена:

• проведению численного моделирования течения жидкости в вихревом теплогенераторе на базе математического пакета «COSMOS FloWorks 2009»;

• проведению экспериментального исследования вихревого теплогенератора и обработке и анализу полученных результатов;

• разработке имитационной математической модели на основании физических и численных экспериментов вихревого теплогенератора.

Первым шагом в разработке численной модели является такая постановка научной задачи, которая может эффективно ответить на вопросы, стоящие перед моделированием. Основными интересующими аспектами является описание течения в проточной части вихревого теплогенератора, выявление особенностей течения и определение возможности увеличения величины влияния нагрева. Для этого необходимо детально рассмотреть картину течения жидкости как во всем объекте, так и в отдельных частях, в связи с чем наиболее значимыми из представленных частей теплогенератора для достижения поставленных целей можно выделить:

• конфузор (увеличивающий скорость потока жидкости);

• улитка (позволяющая закручивать поток жидкости);

• устройство торможения потока (позволяющее уменьшить величину закрутки потока практически до нуля и увеличить время контакта жидкости с внутренними стенками винтовой головки);

• устройство возврата потока (возвращающее поток жидкости).

Для разработки твердотельной модели использовался пакет «SolidWorks 2009», твердотельная модель представляет собой сборку, аналогичную реальному объекту, выполненную в соответствии с особенностями проточной части и конструктивным исполнением.

3. ФОРМИРОВАНИЕ МАТЕМАТИЧЕСКОЙ МОДЕЛИ ЧИСЛЕННОГО МОДЕЛИРОВАНИЯ

Для анализа течения жидкости вихревой теплогенератор условно разделен на 5 участков (рис. 2):

• насос,

• конфузор,

улитка (завихритель),

участок трубы от улитки до крестовины,

крестовина.

Рис. 2. Т вердотельная модель вихревой трубы

Обобщенная имитационная математическая модель вихревого теплогенератора в первом приближении основана на приведенных ниже уравнениях.

Уравнение Навье-Стокса для вязкой несжимаемой жидкости [2]: дУ. дУ.

г- + V —L = F -1-Эр + V -V2 • У.

dt к Эх, 1

dV,

Эх

- = 0,

, р Эх.

к v і

(i, к = 1,2,3).

к

Уравнение неразрывности [2]: V • V = 0. Уравнение энергии [2]:

Р-

f v2 ^

U + —

v 2/

= pF • V + div

Коэффициент турбулентности потока ^t [2]:

m t = Qp

к2

где к - турбулентная кинетическая энергия; є -турбулентная энергия диссипации; Сц - константа, определяющая турбулентную вязкость.

Для нахождения к и є требуется решить два дополнительных уравнения - уравнение турбулентной кинетической энергии и уравнение скорости диссипации. Вид данных уравнений зависит от выбора модели турбулентности. Для стандартной к - є модели:

уравнение для массовой плотности турбулентной энергии к [2]:

э(Рк) Эt

+ V(p Vk ) = V

ff mt ^

m+—

VV Sк J

\

Vk

+ m tG -pe;

уравнение для скорости диссипации турбулентной энергии [2]:

e

д

дt = У

+ У(рУе) =

(( тл

т+— а

1 е е2

Уе + т о - с2р—,

'е у

где

О = D

k

V

ЭХт

k

УУ + рк-

V

6;

1 У

ЭУ , эу

Эх-,

+ -

х2 Эх1

Принимаются следующие значения: як = 1; е = 1,3; Сц = 0,09; С1 = 1,44; С2 = 1,92;

Удельное количество энергии [3, 4]:

• ср -ЭТ-, к = у1, 1г, кг ,

Чк = -

^ + ¿1 ^ ЧРГ а с у

Р Эхк

где ас = 0,9; Рг - число Прандтля; ср - удельная теплоемкость при постоянном давлении; Т -температура текучей среды.

Нагрев жидкости при прохождении через насос за один цикл определяется уравнением

[5]:

_у

с „

ДТ = Т - Т = — (Р - Р )—П

У1Х УХ V У1Х УХ /

п

+

+ £

У1Х 2ухГ ч-4 У1Х ух /’

где ТУ1Х и ТУХ - соответственно температура на входе и выходе насоса; V - удельный объем жидкости:

1

П=-; р

Ср - удельная теплоемкость воды; РУ1Х и РУХ -давление на входе и выходе насоса соответственно; Di - коэффициент Джоуля-Томсона.

Учет нагрева жидкости в результате трения и кавитационного эффекта в первом приближении можно описать уравнениями гидравлических потерь давления и числом кавитации.

Уравнение зависимости изменения температуры от потерь давления в местных сопротивлениях и по длине тракта генератора:

Др

Д = К •% •---------,

С • р р

где Ар - перепад давления на улитке; Ср - теплоемкость; р - плотность; К - эмпирический коэффициент.

Уравнение зависимости безразмерного числа кавитации от скорости течения жидкости [6]:

2 • (Р - Рн.п.) плл

% =----------2-----,

р • и2

где р - давление на участке; и - скорость течения жидкости на участке; р - плотность; рип -давление насыщенных паров.

При достижении числа кавитации предельно допустимого (критического) значения в рассматриваемом местном сопротивлении начинается кавитация. В первом приближении для местных сопротивлений, вызванных изменением сечения потока, можно предложить зависимость для критического числа кавитации [10]:

%„ =Х + 2 V?.

Зная критическое число кавитации %кр для рассматриваемого местного сопротивления, можно определить предельную допустимую скорость перед сопротивлением по формуле [6]:

и <

^пр -

2

. ре кр

Обобщенная имитационная математическая модель вихревого движения и нагрева в теплогенераторе состоит из уравнений Навье-Стокса для вязкой несжимаемой жидкости в случае изотермического движения, а также ряда уравнений, описывающих нагрев жидкости при условии эффекта кавитации.

Эмпирическое выражение изменения температуры на элементе вихревой трубы как функция от перепада давления:

ДР1

Ср р

1 = у1, &, кг.

Суммарный нагрев теплогенератора:

Д1х =

ДРх

Срр

-х

х(К % к + К у1 % у1 + К1г % 1г + Ккг% кг ) .

Суммарные потери в теплогенераторе

ДР£=ДРк +ДРу + ДР +ДРкг,

где АРк - перепад давления в конфузоре; АРу -перепад давления в улитке; АР1г - перепад давления на участке трубы от улитки до крестовины; АРкг - перепад давления в крестовине.

Для численного моделирования в качестве объекта выбрана цилиндрическая вихревая труба. Основные параметры моделируемой вихревой трубы полностью соответствуют имеющейся экспериментальной вихревой установке, на которой проводились натурные исследования. Расчетная область определяется координатами

с

р

x, y, z декартовой системы, задающимися путем твердотельного моделирования проточной части в пакете «SolidWorks». Геометрия твердотельной модели соответствует геометрии трубы экспериментального стенда.

Расчетная область (Computational Domain) -прямоугольный параллелепипед, ограничивающий область, в которой производится расчет (рис. 3).

Для создания сетки расчетной области использовался программный продукт «COSMOS FloWorks», предназначенный для создания поверхностных и внутренних сеток, а также решения представленной модели.

Рис. 3. Область расчетной сетки

Количество объемных элементов в расчетной сетке составило 880 000 элементов.

Численные расчеты производились в компьютерном пакете по вычислительной гидродинамике «COSMOSFloWorks». На твердые поверхности ставилось граничное условие Wall -твердая непроницаемая адиабатическая стенка. На входной поверхности задавалось граничное условие Inlet Volume Flow - входной объемный расход, с заданием температуры потока и параметров турбулентности. На выходных поверхностях задавались граничные условия по давлению, а по температуре и параметрам турбулентности устанавливался нулевой градиент.

1. Начальные условия:

• жидкость несжимаема;

• давление на входе в теплогенератор постоянное Рвх=6,5^105 Па;

• температура на входе в теплогенератор Твх=20 °С.

2. Граничные условия:

• модель турбулентности к-е, где к = = 1 Дж/кг, е =1 Вт/кг;

• максимальный объемный расход насоса Gmax=50 м3/час;

• давление на выходе из теплогенератора Рвых=3>105 Па.

4. АНАЛИЗ И ТРАКТОВКА РЕЗУЛЬТАТОВ

В результате проведения численного моделирования получен ряд характеристик, которые представлены в виде цветовых полей распределения параметров, в форме векторных полей, в форме изолиний и изоповерхностей, а также в виде графиков. Численное моделирование не отражает реальную физику происходящих процессов в вихревом теплогенераторе, однако дает качественную картину распределений параметров как статических, так и динамических, а также дает возможность визуализировать картины течений линий тока жидкости по внутренней поверхности теплогенератора [7, 8].

Расчет проводился на один цикл проходного потока, создаваемого напором насоса. Максимальный объемный расход насоса Qmax = = 50 м3/час, объем жидкости У = 0,02 м3. Таким образом, время прохода жидкости через теплогенератор за один цикл равно

У 0 02

1 =----= —— = 0,0004 часа .

0тах 50

Выбрана наиболее длинная линия тока, проходящая по всему теплогенератору без байпасной линии (рис. 4), то есть на графиках, приведенных далее, представлены характеристики, построенные по тракту вихревого теплогенератора.

Соответствие длины линии тока по пройденным участкам вихревого теплогенератора следующее:

• конфузор - 0,272 м;

• улитка - 0,100 м;

• вихревая камера до крестовины -1,470 м;

• крестовина - 0,100 м;

• участок трубы после крестовины -

0,155 м.

При прохождении жидкости в теплогенераторе за один цикл разность температур между входной точкой и точкой на выходе из теплогенератора составила АТ = 0,1 °С (рис. 5, 6).

і

Рис. 5. Температура жидкости в улитке

Температура

20,13 20,12 20,11 20,1 20,09 20,08 20,07 20,06 20,05 20,04

0 0, 1 1, 2 Длина линии тока

Рис. 6. Распределение температуры жидкости в теплогенераторе

При моделировании циклического течения жидкости за время равное 16 секундам, при этом выходные данные по температуре были входными данными следующего цикла, наблюдалась картина линейного повышения температуры.

Рис. 7. Статическое давление жидкости в улитке

Распределение статического давления в теплогенераторе показывает неравномерное поле (рис. 7-9). В центре основной трубы давление значительно ниже, чем на периферии. Таким образом, существует некая граница двух вихревых токов, которые непосредственно взаимодействуют. Дальнейшее движение жидкости по

У

гТ

^ к

тракту показывает выравнивание статического давления.

Статическое давление

Длина линии тока

Рис. 8. Распределение статического давления по тракту теплогенератора

На рис. 10 представлен график полной скорости течения жидкости в вихревой трубе. Видно, что скорость течения в периферийной области превышает скорость центрального течения в 2 раза.

Скорость

Длина линии тока

Рис. 10. График распределения скорости жидкости по тракту теплогенератора

Некоторые авторы, например Ю. С. Потапов, утверждают, что при попадании жидкости в вихревую камеру (улитку) часть потока, нагреваясь, проходит через байпасную линию и встречается с основным потоком на участке

после крестовины [9]. Однако при анализе линий тока (рис. 11, 12) жидкости в вихревом теплогенераторе, которые строились на основе численного моделирования, выяснилось, что весь поток жидкости проходит через улитку. Когда основной поток достигает крестовины, часть подогретой жидкости проходит через байпасную линию, тем самым подогревая следующую «порцию» жидкости. Этот эффект подобен положительной обратной связи, но, к сожалению, на сегодняшний день полностью не изучен.

Рис. 11. Направление линий тока жидкости в вихревом теплогенераторе

Рис. 12. Направление линий тока жидкости в улитке

Результаты решения математической модели в пакете «CosmosFloWorks» показали характер внутреннего течения рабочей жидкости в вихревом теплогенераторе. Однако для подтверждения полученных результатов требуются экспериментальные данные натурных экспериментов.

Наиболее интересным участком является крестовина вихревого теплогенератора.

5 0,1

,,Д41У а

ГГ 1 И 1

Л11

/ * 10 2 II 31 0 4< 10 5 О 0

Время {с|

Рис. 13. График разности температуры на участке «вход в крестовину - выход из крестовины»

Из графика (рис. 13) видно, что скорость нагрева жидкости на крестовине увеличивается с ростом температуры. Так как крестовина является существенным местным сопротивлением, способствующим резкому увеличению скорости потока рабочей жидкости, можно сделать предположение о наличие кавитационного эффекта на данном участке. При аппроксимации аналитического графика получено эмпирическое уравнение интенсивности (скорости) нагрева жидкости на крестовине в вихревом теплогенераторе

АТ = 0,087 • Ьп

где АТ - разность температур на участке крестовины, ^ - время.

Учитывая специфику решения математической модели, сравнение данных натурного эксперимента с результатами численного моделирования производилось для интегральных параметров численного решения с целевыми функциями.

На графике (рис. 14) сопоставлены результаты натурного и численного экспериментов за время равное 16 секундам на участке «выход из крестовины». Так как математический пакет «CosmosFloWorks» не предусматривает моделирование циклических процессов, при которых меняются граничные условия в течение времени, моделирование проводилось решением алгебраических уравнений имитационной математической модели.

Расхождение температурных данных натурного эксперимента с результатами численного моделирования объясняется наличием интенсивного теплоотвода в окружающую среду и эффекта кавитации в вихревом теплогенераторе при проведении натурного эксперимента. Численное исследование не учитывает теплообмен с окружающей средой и повышение температуры от эффекта кавитации, которое стано-

вится интенсивнее с увеличением количества прохождения циклов рабочей жидкости через вихревой теплогенератор.

Следовательно, на данном этапе исследований является целесообразным введение поправочного коэффициента, зависимого от текущей температуры для каждого участка теплогенератора.

Уравнение зависимости поправочного коэффициента от температуры при аппроксимации на участке «выход из крестовины»:

Ккг = 0,37 • ґ0’33,

где Ккг - поправочный коэффициент; ґ - температура.

22,5---------------------------------------------

22___________________________________________—

о 21,5 - ------

о.

<е

0 2 4 6 8 Ю 12 14 16

Время

ч- при численном моделировании при натуральном эксперименте

Рис. 15. График поправочного коэффициента на участке «выход из крестовины»

-J

м

/

/  .1 г

ж .*

г- Ґ-

О 2 4 6 8 10 12 14 16

Время

при численном моделировании при натуральном эксперименте

Рис. 14. График сопоставления температур на участке «выход из крестовины»

Таким образом, в результате подстановки поправочного коэффициента в имитационную математическую модель температурная интегральная характеристика, полученная при численном моделировании, имеет расхождение

с натурным экспериментом в пределах ±0,3% (рис. 15).

Сравнение результатов расчета с экспериментами подтверждает правильность теоретического подхода к определению статических характеристик.

Анализ адекватности математической модели показывает достаточно точную сходимость результатов ее решения с натурным экспериментом, и модель рекомендуется для проведения численных исследований вихревого эффекта и проектирования вихревых теплогенераторов.

ВЫВОДЫ

1. Разработана и решена в пакете Cosmos FloWorks система уравнений имитационной математической модели процесса в вихревом теплогенераторе в трехмерной постановке, с использованием стандартной ¿-s-модели турбулентности и уравнения нагрева жидкости. Результаты вычисления показали наличие поля температур, подтверждающее возможность теплообмена и возникновения кавитационных эффектов как условия нагрева жидкости.

2. Проведены экспериментальные исследования процессов высоконапорных вихревых течений в вихревом теплогенераторе на натурном стенде. Время нагрева 20 литров жидкости до 100 °С составило 630 с. В процессе проведения натурного эксперимента было выявлено выделение конкурентной газовой фазы, предположительно из-за эффекта кавитации.

По результатам экспериментальных и численных исследований была проведена идентификация имитационной математической модели. При сопоставлении температурной характеристики был введен эмпирический поправочный коэффициент для каждого участка вихревого теплогенератора. Таким образом, с учетом поправочного коэффициента расхождение температурной составляющей составило ±0,3%.

Исследования выполнены в соответствии с планом работ по направлению «Новые и возобновляемые источники энергии» федеральной целевой программы «Научные и научнопедагогические кадры инновационной России» на 2009-2013 годы.

СПИСОК ЛИТЕРАТУРЫ

1. Меркулов А. П. Вихревой эффект и его применение в технике. М.: Машиностроение; Самара: Оптима, 1997. 292 с.

2. Лойцянский Л. Г. Механика жидкости и газа: Учеб. для вузов. М.: Дрофа, 2003. 840 с.

3. Абрамович Г. Н., Степанов Г. Ю. Гидродинамика закрученного потока в круглой трубе с внезапным увеличением поперечного сечения и при истечении через насадок Борда // Механика жидкости и газа.1994. № 3. С. 51-66.

4. Кныш Ю. А. Физическая модель явления энергопереноса в вихревой трубе // Вихревой эффект и его применение в технике. Куйбышев: КуАИ. 1988. С. 71-74.

5. Ким Д. П., Рахматуллин Ш. И. Степень влияния дифференциального напора, подачи и адиабатического сжатия нефти на ее нагрев в центробежном насосе // Нефтегазовое дело. 2005.

6. Идельчик И. Е. Справочник по гидравлическим сопротивлениям. М.: Машиностроение, 1975. 559 с.

7. Полянский А. Ф., Скурин Л. И. Моделирование течений жидкости и газа в вихревой трубе и струе // Математическое моделирование. 2001. Т. 13, № 7. С. 116-120.

8. Калимуллин Р. Р., Ахметов Ю. М. Моделирование течения жидкости в вхревом теплогенераторе // Динамика машин и рабочих процессов: Тр. Все-российск. науч.-техн. конф. Уфа: УГАТУ, 2009. С. 8-10.

9. Фоминский Л. П. Теплогенератор Потапова -работающий реактор холодного ядерного синтеза // Электрик. 2001. № 1. С. 19-21.

10. Червяков В. М., Юдаев В. Ф. Гидродинамические и кавитационные явления в роторных аппаратах: монография. М.: Машиностроение-1, 2007. 128 с.

ОБ АВТОРАХ

Ахметов Юрий Мавлютович,

доц каф. прикл. гидромех., зам. ген. дир. НИИТ. Дипл. инж.-мех. по авиац. двиг. (УАИ, 1959). Канд. техн. наук по тепл. двиг. (МАИ, 1978). Иссл. в обл. газо-гидр. течений и систем упр. энерг. установок.

Калимуллин Радик Рифкатович, асп. той же каф. (УГАТУ, 2009). Дипл. инж. техн. и технол. по энергомашиностр. (УГАТУ, 2008). Иссл. в обл. высоконапор. многофаз. течения жидкости.

Целищев Владимир Александрович, зав., проф. той же каф. Дипл. инж.-мех. по гидравл. машинам (УГАТУ,1982). Д-р техн. наук по тепловым двигателям (УГАТУ, 2000). Иссл. в обл. автоматики энергетических установок.