Научная статья на тему 'Исследование процесса теплоотдачи в сетчатой матрице роторного теплообменника'

Исследование процесса теплоотдачи в сетчатой матрице роторного теплообменника Текст научной статьи по специальности «Механика и машиностроение»

CC BY
248
47
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
КОМПАКТНЫЕ ТЕПЛООБМЕННИКИ / МОДЕЛИРОВАНИЕ ТЕПЛОГИДРАВЛИЧЕСКИХ ПРОЦЕССОВ / КОЭФФИЦИЕНТ ТЕПЛООТДАЧИ ПОРИСТЫХ МАТРИЦ / COMPACT HEAT EXCHANGERS / HYDRAULIC PROCESS MODELING / HEAT TRANSFER COEFFICIENT OF POROUS MATRICES

Аннотация научной статьи по механике и машиностроению, автор научной работы — Алексеев Р. А., Костюков А. В.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Investigation of the heat transfer process in a mesh matrix of a rotary heat exchanger

The paper presents the results of computational modeling of thermal processes in mesh matrixes of rotary frame heat exchangers. The authors received the refined dependence of the mesh matrix heat transfer of the coolant flow characteristics.

Текст научной работы на тему «Исследование процесса теплоотдачи в сетчатой матрице роторного теплообменника»

• повышение грузоподъёмности автопоезда;

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

• возможность объединения прицепных звеньев для транспортировки длинномерных и крупногабаритных тяжёлых грузов;

• снижение разрушающего воздействия на почву до 40%.

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

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

Литература

1. Шухман С.Б., Соловьев В.П., Прочко Е.И. Гидрообъёмные передачи - перспектива для полноприводных АТС. «Автомобильная промышленность», 1997, № 6. - с. 21-23.

2. Пирковский Ю.В., Шухман С.Б. Теория движения полноприводного автомобиля (прикладные вопросы оптимизации конструкции шасси). - М.: ЮНИТИ-ДАНА: Элит-2000, 2001.-230 е., ил.

3. Коркин С.Н., Курмаев Р.Х., Петров С.Е. Экологические аспекты движения тяжёлых специальных автопоездов в условиях бездорожья // Материалы Международной научно-технической конференции Ассоциации автомобильных инженеров (ААИ) «Автомобиле- и тракторостроение в России: приоритеты развития и подготовка кадров», посвященной 145-летию МГТУ «МАМИ». - М.: МГТУ «МАМИ», 2010 г. - Книга 1. - с. 179 - 185.

4. Лепёшкин A.B., Курмаев Р.Х., Бочаров С.Е., Куру Д.С. Приводы активизированных прицепов автопоездов // Материалы Международной научно-технической конференции Ассоциации автомобильных инженеров (ААИ) «Автомобиле- и тракторостроение в России: приоритеты развития и подготовка кадров», посвященной 145-летию МГТУ «МАМИ». -М. : МГТУ «МАМИ», 2010 г. - Книга 1. - С. 227 - 231.

5. Лепешкин A.B. Математическая модель установившегося движения автопоезда с активизированным прицепным звеном, позволяющая оценить потери в трансмиссии. Известия МГТУ «МАМИ». Научный рецензируемый журнал. - М., МГТУ «МАМИ», №2 (12), 2011. -с. 27-41.

Исследование процесса теплоотдачи в сетчатой матрице роторного

теплообменника

доц. к.т.н. Костюков A.B., Алексеев P.A.

Университет машиностроения ronopolk88(a>,mail. ru, kostukovl23(a),yandex.ru

Аннотация. Приводятся результаты расчетных исследований тепловых процессов в сетчатых матрицах роторных каркасных теплообменников. Получена уточненная зависимость теплоотдачи сетчатой матрицы от характеристик потока теплоносителя.

Ключевые слова: компактные теплообменники, моделирование теплогид-равлических процессов, коэффициент теплоотдачи пористых матриц.

Одним из путей повышения эффективности применяемых в настоящее время в распре-

Раздел 1. Наземные транспортные средства, энергетические установки и двигатели, деленной энергетике микротурбин является установка в них теплообменников со сверхвысокой степенью регенерации (95-97%). По этому направлению идет американская компания Wilson, разрабатывающая микротурбину мощностью 300 кВт с электрическим КПД 50% [1]. Получение степени регенерации на уровне 95-97% при приемлемых габаритах возможно только в компактных роторных теплообменниках. В таких теплообменниках, как правило, используется пористая теплопередающая матрица. В частности, в роторном теплообменнике транспортного газотурбинного двигателя (ГТД) Горьковского автозавода (ГАЗ) применяется пористая матрица, образованная намоткой металлических сеток [2]. Имеющаяся информация по теплопередаче в таких сетках весьма невелика. В частности, в Николаевском кораблестроительном институте им. С.О. Макарова проводилось исследование теплогидравлических характеристик пористой сетчатой матрицы, приметаемой в теплообменниках ГТД ГАЗ [3,4,5,6,7]. Однако работы проводились в очень узком диапазоне температур матрицы и теплоносителей (в эксперименте сетчатая матрица обдувалась попеременно воздухом с температурами 50°С и 30°С соответственно, в этом же диапазоне металась температура сетчатой матрицы) и в полученной эмпирической зависимости (см. ниже) отсутствуют элементы, учитывающие температуру сетчатой матрицы.

Jk= 0.11-Rex'046 , (1)

n 0.00375 где: Rex - модифицированное число Рейнольдса: Rex = Re--—^-;

J, - фактора Колборна: J, = St- Рг2 3 =-——;

" 1 ^ -Re-Pr

ОС - коэффициент теплоотдачи матрицы;

dv - гидравлический диаметр матрицы;

X - коэффициент теплопроводности;

St - число Стэнтона;

Рг - число Прандтля;

Re - число Рейнольдса;

Р - пористость матрицы.

Следствием этого явилось значительное, на 3-4% (абсолютных), несовпадение экспериментально полученной на режиме с температурой матрицы, лежащей в диапазоне 230-600иС, и расчетной (с использованием зависимости (1)) степени регенерации теплообменника ГТД ГАЗ [2].

Уточнению эмпирической зависимости (1) теплоотдачи для пористой сетчатой матрицы посвящена данная статья.

Для выполнения поставленной задачи проводилось математическое моделирование теплогидравлических процессов в сетчатой матрице. В качестве объекта исследования была выбрана сетчатая матрица, приметаемая в роторных каркасных теплообменниках ГТД ГАЗ [2]. Основные геометрические параметры сетчатой матрицы: диаметр проволоки d = 0.13 мм, "Ячейка в свету" S = 0.2mm, пористость Р = 0.69, компактность

К = 9508 м2/м3, гидравлический диаметр dr = 0.291 мм

Схематичное изображение слоя сетки показано на рисунке 1.

Расчетная модель представляла собой фрагмент теплопередающей матрицы, состоящий из десяти слоев сеток (рисунок 2). Математическое моделирование базировалось на решении системы уравнений, включающей уравнения Навье-Стокса, энергии, неразрывности и состояния. Физические свойства стали и воздуха брались переменными в зависимости от температуры. В качестве материала теплопередающей матрицы элемента выбрана сталь 12X18Н9Т.

Зависимости её теплофизических свойств от температуры приведены в таблице 1.

О.ООГСб О.СХЮ75

Рисунок 2 - Исследуемый участок сетчатой матрицы

Таблица 1

Зависимость коэффициента теплопроводности, удельной теплоёмкости и плотности

стали 12Х18Н9Т от температуры

Температура, К Теплопроводность, Вт/м2 ■К Теплоемкость, Дж/кг-К

273 15 457

373 16 469

573 20 498

873 25 528

973 26 541

Теплофизические свойства воздуха были аппроксимированы в виде полиномов в зависимости от температуры.

Динамическая вязкость:

/1 = (5.0697 + 0.005006 • Т - 0.00001262 • Т2 )/\000000 Па ■ сек Теплоемкость:

С= 933.91 + 0.19019 • Т + 0.000016667 • Т2 Дж/кг ■ К Теплопроводность:

Л = (0.233 + 0.00891 • Т + 0.0000022857 • Г2)/1005да/л/2 • К Численный эксперимент проводился в области чисел Рейнольдса Яе = 10-70 (ламинарный режим течения). Продувка воздухом матрицы проводилась до момента полного ее нагрева или охлаждения. Температуры воздуха Тв и матрицы '/',, брались приближенными к их реальным значениям в роторном теплообменнике ГАЗ. В случае нагрева: Тв = 965К,Тм = 300^

В случае охлаждения: Тв = 500К,Тм = 965К

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

1. Для всякого выбранного расчетного момента времени I были получены средние интегральные величины температур матрицы (Тм) и воздуха в ней ( Тг), а также теплофизиче-

ских свойств теплоносителей.

2. Вычислялся коэффициент теплоотдачи и критерии подобия Прандтля:

АН-в „ |ы-С

а

Рг

Рт,.-\Тв-Тм\ X

где: АН - разности энтальпий воздуха на входе и на выходе; О - расход воздуха;

¥т 0 - полная площадь теплообмена в данном объеме сетчатой матрицы; |Ы - динамическая вязкость; С - коэффициент теплоемкости;

X - коэффициент теплопроводности. 3. Находился фактор Колборна:

Л

а •

X • Re- Рг°33

4. Вычислялся модифицированный критерий Рейнольдса:

0.00375

Rev = Re

Р3

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

Те

Д

возд вх ' ^ ® охлаждения

г„

"нагрева ^ мапгр / ^ возд_вх ' ^ ~охлаждения ^ возд_вх /^матр '

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

Jt_Te_0.502S2_lRsX.i „ -Jt_T5_S.65S4:_:RsX' Jk_Te_0. S3041_iReX/0:i3 Jk_Ts_C.954j4_:RäX' 02 Jt_Te_0.99123_|Rsx}

Ж

1

Jt(R=x) := O-U-fRex)

щ. |

1г_Те_0.53374_1Еех) ' k_Te_0 <50067JRix) 3 ; k_Ts_iiSS:5S_:RiX

ъъьтпЫ*^ t_It_O S«fflJR5x) 03

k Te_0 97367 JRsX

'¿/А

0 U- . — ReX.' 0.46 1

1 t |

0.1 0.: 0.3 0.4 0.5 0.6 0.7 0.S 0.9

R=x

0 Ol 02 5 3 04 0 5 06 ST 0.S 0 5 l R=X

a)

6)

Рисунок 3 - Полученные зависимости Jk = /(Rex), соответствующие различным значениям температурного фактора в случае: а) нагрева матрицы; б) охлаждения

матрицы

Тройки значений (Rex, Те, Jk) аппроксимировались следующим образом: сначала для каждого выбранного Te=const была проведена аппроксимация пар значений (Rex,l/jt) функцией:

Jk=A-ReJ.

Далее находились аппроксимационные зависимости полученных коэффициентов А и В от температурного фактора Те в виде:

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

А = М1- Те+ Л^ В = М2 • Те + Л^

В итоге была получены зависимости:

^ = (0.029 • Те+ 0.06) • Ке/00б31-Те-°-514,

Лохл = (0.019-Те+0.064)-Ке/(0143Те+°32). (2)

Поверка полученных зависимостей проводилась путем расчета в конечно-элементном комплексе теплогидравлических процессов в пористом теле с заданием процессов теплоотдачи зависимостями (2), а также формулой (1). Используемая расчетная методика приведена в работе [8].

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

*** Расчет по зависимостям (2) Численный эксперимент Расчет по зависимости (1)

б)

Рисунок 4 - Графики изменения во времени температур воздуха на выходе из расчетной модели и средних температуры матрицы в случае: а) нагрева матрицы;

б) охлаждения матрицы

Сравнение этих значений с результатами численного эксперимента по продувке фрагмента сетчатой матрицы (рисунок 2) показано в виде графиков на рисунке 4 а, б.

Как видно, расчет с применением формул (2) дает более близкое совпадение с численным экспериментом.

Выводы

1. Проведено расчетное исследование тепловых процессов в сетчатой матрице роторного каркасного теплообменника.

2. Для случаев нагрева и охлаждения получены зависимости фактора Колборна (связанного с коэффициентом теплоотдачи матрицы) от параметров потока теплоносителя с учетом температурного фактора.

3. Результаты поверочного расчета с применением полученных зависимостей показали хорошее совпадение с численным экспериментом.

Литература

1. Dr. David Gordon Wilson Professor Emeritus, MIT Chief Scientist, Wilson TurboPower Inc. «The high efficiency of our multistage ceramic-bladed turbine design can produce an attractive economic return even in the absence of CHP".

2. Плотников Д.А. Разработка и исследование дисковых секционных регенераторов автотракторных ГТД. Диссертация на соискание ученой степени кандидата технических наук, 1981 г.

3. Лебедь Н. Исследование гидравлических и теплообменных характеристик теплопередаю-щих поверхностей для вращающегося регенератора секционного типа. Научно-технический отчет, 1971 г.

4. Исследование роторного теплообменника малоразмерного регенеративного газотурбинного двигателя. Общероссийский н/т журнал "Полет" №1, 2005г., с.38-42. A.B. Костюков и др.

5. Расчет распределения температур в каркасе дискового роторного теплообменника при нестационарных граничных условиях. Общероссийский н/т журнал "Полет" №10, 2005г., с.14-18. A.B. Костюков и др.

6. Высокотемпературный роторный теплообменник для газотурбинных и героторных двигателей. Известия МГТУ "МАМИ" №2 (12), 2011г., с.23-26. A.B. Костюков и др.

7. Исследование системы охлаждения каркаса роторного теплообменника. Известия МГТУ "МАМИ" №2 (4), 2007г., с.61-65. A.B. Костюков и др.

8. Алексеев P.A., Костюков A.B. Повышение эффективности роторного теплообменника малоразмерного газотурбинного двигателя. Известия МГТУ "МАМИ" №1 (13), 2012 г., с.52-58.

9. Справочник по теплофизическим свойствам газов и жидкостей. Н.Б. Варгафтик, 1972г.

Исследование режимов горения энергоаккумулирующих веществ

к.т.н. проф. Костюков A.B., доц. Кузнецов В.В., доц. Ащеульников Е.К., Родькин К.П.

Университет машиностроения kostukovl23(ci),yandex. ru Аннотация. В статье рассматриваются модель горения порошкообразных энергоаккумулирующих веществ (ЭАВ) в цилиндрической камере сгорания. Такая камера может рассматриваться как форкамера водородного реактора для транспортных и стационарных энергоустановок на базе ГТД. Составлена программа расчета горения порошкообразных ЭАВ в потоке перегретого водяного пара. Получено распределение температуры по оси камеры при различных значениях коэффициента избытка окислителя - водяного пара. Определены размеры глубины фронта пламени при заданных условиях работы камеры сгорания на энергоаккумулирующих веществах.

Ключевые слова: энергоаккумулирующие вещества, модель горения в цилиндрической камере сгорания

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

i Надоели баннеры? Вы всегда можете отключить рекламу.